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

8. Recovering a synthetic 3D subsurface model using lower-dimensional inversion schemes f (degrees) r A (log 10(Wm)) f (degrees) r A (log 10(Wm)) Period = 0.01 s Period = 100 s XY data YX data XY data YX data 2.3 1.7 46 44 Period = 1,000 s Period = 60,000 s XY data YX data XY data YX data 2.7 2.2 42 17 Crustal strike: XY data = TM mode; YX data = TE mode Mantle strike: XY data = TE mode; YX data = TM mode 2.3 1.7 46 44 2.9 2.1 33 28 Fig. 8.6.: Maps of apparent resistivity (ρa) and impedance phase (φ) of the 3D synthetic model (Fig. 8.3) at four periods representing different regimes; see text for details. Note that different colour scales are used to highlight structures of each plot. Plots are rotated anticlockwise by 45 degrees in respect to Figures 8.4 and 8.5 in order to accommodate plotting of multiple figures side-by-side; as a result North is located towards the top-left. the anomaly magnitude (i.e. the absolute difference of φ values from 45 degrees) being controlled by the resistivity difference between crust and mantle (highest in the northern quadrant: 50 to 1000 Ωm) and the induction depth (greater in the southwestern half). Response data for periods around 1000 s (bottom left figures) also exhibit point symmetry with the strongest phase anomalies (lowest φ values) located in the northern quadrant. In the southwestern half, phase values are closer to 45 degrees and apparent resistivity values are closer to values of the synthetic model mantle owing to the higher resistivity of local crustal structures and the resulting greater induction depth. Responses for the longest periods of the dataset (6∗10 4 s, bottom-right figures), are still affected by distortion of crustal structures. Phase data, for both modes, differ between the four quadrants (N, S, E, W), thus indicating effects from different induction depths and resistivity contrasts of the regions. ρxy data are similar to mantle values of the synthetic 176 2.5 1.6 39 28 2.9 2.7 43 39 2.5 1.6 37 30 3.1 2.4 42 40

3D N45W 3D-crust TE Rho TE Phi Period (s) Period (s) TM Rho TM Phi Period (s) Period (s) 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 syn020 syn019 syn017 syn015 syn013 syn011 syn009 syn007 syn006 syn005 syn004 syn003 syn002 syn001 0 50 100 150 Distance (Km) 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 Ohm.m 794 533 358 240 161 108 73 49 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 Ohm.m 1585 970 594 364 223 136 83 51 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 Degrees Degrees 90 80 70 60 50 40 30 20 10 90 80 70 60 50 40 30 20 10 8.2. Generating synthetic 3D model data 3D N45E 3D-mantle TE Rho Period (s) TE Phi Period (s) TM Rho Period (s) TM Phi Period (s) 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 syn020 syn019 syn017 syn015 syn013 syn011 syn009 syn007 syn006 0 50 100 Distance (Km) Fig. 8.7.: Pseudosections of the profiles 3D-crust (left-hand side) and 3D-mantle (right-hand side) with stations representing the MT recording sites of the PICASSO Phase I investigation; see Figures 8.4 and 8.5 for the location of stations and profiles. The northern end of the profile (NE for 3D-crust profile, and NW for 3D-mantle profile) is shown to the right. model (1000 Ωm in the northwestern half, and 500 Ωm in the southeastern half), however, ρyx values are significantly different. YX data at long periods exhibit significant distortion effects, originating from crustal structures and the oblique strike directions between the two depth regions. Among other, ρyx data of the southern quadrant are higher than for the northern quadrant resulting in an issue concerning, in particular, profiles using the pic-stations (pic001 - pic020, denoted by stars in Figure 8.4). The circumstance that YX responses of the southern mantle region are significantly more resistive than responses in the northern mantle region is reflected in TM mode data of the crustal strike direction and TE data of the mantle strike direction (see also Figure 8.7). Forward response data can be displayed using so-called pseudosections, i.e. gridded resistivity–period values beneath the profile, in which period is used as a depth proxy. The pseudosection for the pic-line (stations syn001 – syn020 in Figure 8.7) is displayed here because these stations form the most challenging combination by containing stations from each resistivity region of the synthetic 3D model, and, moreover, because the station locations are representations of the PICASSO Phase I recording sites (cf. Sec. 9.1). Due to the characteristics of the 3D model, two profiles are created with data arranged according to the geoelectric strike directions at crustal depth (N45W, left-hand side plot in Figure 8.7) and mantle depth (N45E, right-hand side plot in Figure 8.7); see Figure syn005 syn004 syn003 syn002 syn001 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 10 -2 10 -1 10 0 10 1 10 2 10 3 10 4 10 5 Ohm.m Degrees Ohm.m Degrees 1585 970 594 364 223 136 83 51 90 80 70 60 50 40 30 20 10 794 533 358 240 161 108 73 49 90 80 70 60 50 40 30 20 10 177

3D<br />

N45W<br />

3D-crust<br />

TE Rho<br />

TE Phi<br />

Period (s)<br />

Period (s)<br />

TM Rho<br />

TM Phi<br />

Period (s)<br />

Period (s)<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

<br />

syn020<br />

syn019<br />

syn017<br />

syn015<br />

syn013<br />

syn011<br />

syn009<br />

syn007<br />

syn006<br />

syn005<br />

syn004<br />

syn003<br />

syn002<br />

syn001<br />

0 50 100 150<br />

Distance (Km)<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

Ohm.m<br />

794<br />

533<br />

358<br />

240<br />

161<br />

108<br />

73<br />

49<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

Ohm.m<br />

1585<br />

970<br />

594<br />

364<br />

223<br />

136<br />

83<br />

51<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

Degrees<br />

Degrees<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

20<br />

10<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

20<br />

10<br />

8.2. Generating synthetic 3D model data<br />

3D<br />

N45E<br />

3D-mantle<br />

TE Rho<br />

Period (s)<br />

TE Phi<br />

Period (s)<br />

TM Rho<br />

Period (s)<br />

TM Phi<br />

Period (s)<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

<br />

syn020<br />

syn019<br />

syn017<br />

syn015<br />

syn013<br />

syn011<br />

syn009<br />

syn007<br />

syn006<br />

0 50 100<br />

Distance (Km)<br />

Fig. 8.7.: Pseudosections of the profiles 3D-crust (left-hand side) and 3D-mantle (right-hand side) with stations representing the MT<br />

recording sites of the PICASSO Phase I investigation; see Figures 8.4 and 8.5 for the location of stations and profiles. The northern<br />

end of the profile (NE for 3D-crust profile, and NW for 3D-mantle profile) is shown to the right.<br />

model (1000 Ωm in the northwestern half, and 500 Ωm in the southeastern half), however,<br />

ρyx values are significantly different. YX data at long periods exhibit significant distortion<br />

effects, originating from crustal structures and the oblique strike directions between the<br />

two depth regions. Among other, ρyx data of the southern quadrant are higher than for<br />

the northern quadrant resulting in an issue concerning, in particular, profiles using the<br />

pic-stations (pic001 - pic020, denoted by stars in Figure 8.4). The circumstance that YX<br />

responses of the southern mantle region are significantly more resistive than responses in<br />

the northern mantle region is reflected in TM mode data of the crustal strike direction and<br />

TE data of the mantle strike direction (see also Figure 8.7).<br />

Forward response data can be displayed using so-called pseudosections, i.e. gridded<br />

resistivity–period values beneath the profile, in which period is used as a depth proxy.<br />

The pseudosection for the pic-line (stations syn001 – syn020 in Figure 8.7) is displayed<br />

here because these stations form the most challenging combination by containing stations<br />

from each resistivity region of the synthetic 3D model, and, moreover, because the station<br />

locations are representations of the PICASSO Phase I recording sites (cf. Sec. 9.1).<br />

Due to the characteristics of the 3D model, two profiles are created with data arranged<br />

according to the geoelectric strike directions at crustal depth (N45W, left-hand side plot<br />

in Figure 8.7) and mantle depth (N45E, right-hand side plot in Figure 8.7); see Figure<br />

syn005<br />

syn004<br />

syn003<br />

syn002<br />

syn001<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

10 -2<br />

10 -1<br />

10 0<br />

10 1<br />

10 2<br />

10 3<br />

10 4<br />

10 5<br />

Ohm.m<br />

Degrees<br />

Ohm.m<br />

Degrees<br />

1585<br />

970<br />

594<br />

364<br />

223<br />

136<br />

83<br />

51<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

20<br />

10<br />

794<br />

533<br />

358<br />

240<br />

161<br />

108<br />

73<br />

49<br />

90<br />

80<br />

70<br />

60<br />

50<br />

40<br />

30<br />

20<br />

10<br />

177

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

Saved successfully!

Ooh no, something went wrong!