P. Schmoldt, PhD - MTNet - DIAS
P. Schmoldt, PhD - MTNet - DIAS P. Schmoldt, PhD - MTNet - DIAS
8. Recovering a synthetic 3D subsurface model using lower-dimensional inversion schemes Distance from the centre of the mesh (km) 300 200 100 0 −100 −200 C−line 11−line 01−line 02−line 03−line 04−line 05−line 06−line 07−line 08−line 09−line 10−line D−line E−line F−line G−line H−line I−line J−line K−line L−line M−line C−line 11−line 00−line D−line E−line F−line 01−line 02−line 03−line 04−line 05−line 06−line 07−line 08−line 09−line 10−line N−line G−line H−line I−line J−line K−line L−line M−line −300 −300 −200 −100 0 100 200 300 Distance from the centre of the mesh (km) Fig. 8.4.: Location of magnetotelluric (MT) recording sites on top of the synthetic 3D subsurface model; North is located to the top of the figure. Displayed are 144 station arranged in a grid (inverted triangles), as well as 13 additional sites, which are a projection of the PICASSO Phase I stations. Locations of these 13 additional stations (pic001 – pic020, situated along the ’pic-line’ in this figure) are a translation of the PICASSO Phase I stations with the respective indices (cf. Chapter 9) onto the synthetic 3D model with pic001 and pic020 located in the North and South, respectively. Nomenclature for the other stations follows the grid system, i.e the name of a station is the combination of the two lines intersecting at the respective location, e.g. the centre-top station is named ‘C00’. Background colours indicate the relatively conductive (red) and resistive (yellow) regions of the crust; the location of the resistivity interface at mantle depth is indicated by the dashed line (cf. Fig. 8.3). Profiles 3D-crust and 3D-mantle (not displayed) coincide with the conductivity interfaces at mantle and crust depth, respectively. (in respect to Figures 8.4 and 8.5) to accommodate plotting of multiple maps side by side. Using Niblett-Bostick depth estimation (Sec. 6.3.1) yields that on the resistive side periods greater than 18 s (appr. 1.26 in log-scale) penetrate into the mantle, whereas on the conductive side penetration into the mantle is first achieved by periods of 72 s (appr. 1.86 in log-scale). In Figure 8.6 ‘XY data’ refers to the TE mode and ‘YX data’ refers to TM mode for the dataset adjusted to the crustal strike direction (N45W), whereas the ‘XY data’ refers to the TM mode and ‘YX data’ refers to TE mode for the dataset adjusted to the mantle strike direction (N45E). The similarity of crustal TE and mantle TM data (and vice versa) is due to the 90 degrees difference between the two strike directions, resulting in swapping of the two modes. Different colour scales are used to display values in maps of different modes and periods (rather than using uniform colour scale values) to highlight structures at the respective period. Note that, in order to enable plotting of all phase data in the first quadrant (0 - 90 degrees), 180 degrees are added to the respective YX phase 174 pic−line pic−line 00−line N−line
Distance from the centre of the mesh (km) 200 100 0 −100 −200 3D−crust 3D−crust−west 3D−crust 3D−crust−east 8.2. Generating synthetic 3D model data −200 −100 0 100 200 Distance from the centre of the mesh (km) Fig. 8.5.: Location of profiles on top of the synthetic 3D model (cf. Figs. 8.3 8.4); North is located to the top of the figure. Stations are associated with profiles of the same colour; for the profile labelled ‘3D-crust’ three different sets of stations are used: red (3Dcrust, using stations representing the PICASSO Phase I recording sites), blue (3D-crust-NS), and green (3D-crust-EW). Every profile contains at least one station on top of each of the four electric conductivity regions (two at crustal depth, two at mantle depth), to assure that data of every profile is affected by the oblique strike directions. Profiles are constructed orthogonal to the geoelectric strike direction of either crust or mantle and station locations are projected onto the profile. Background colours indicate the conductivity distribution at crustal depth (cf. Fig. 8.4). data. At shorter periods (0.01 s, top-left figures), responses are dominated by characteristics of crustal structures, therefore values for the two modes are alike and ρa = ρ; hence ρa = 50 Ωm (appr. 1.7 in log-scale) in the northeastern half and ρa = 200 Ωm (appr. 2.3 in log-scale) in the southwestern half. For longer periods (100 s, top-right figures) mantle structures start to add observable contributions to the response data. Values of ρa at 100 s are similar to values of ρa at 0.01 s (note the different colour scale), but φ differs significantly. In general, φ(100 s) < φ(0.01 s) ≈ 45 degrees owing to the more resistive nature of mantle regions in respect to the crust. At 100 s periods, skin depth (cf. Sec. 3.3) for stations on the resistive side of the crustal fault is appr. 70 km, whereas on the conductive side it is appr. 35 km. Thus, at 100 s periods all sites are sensitive to electric properties of the mantle. Phase anomalies at 100 s exhibit a point symmetry in regards to the centre of the station array, i.e. where crust and mantle interfaces intersect. Phase anomalies are a superposition of effects from crustal and mantle structures with 175
- Page 160 and 161: 6. Using magnetotellurics to gain i
- Page 162 and 163: 6. Using magnetotellurics to gain i
- Page 164 and 165: 6. Using magnetotellurics to gain i
- Page 168 and 169: Part II Geology of the study area I
- Page 170 and 171: 7. Geology of the Iberian Peninsula
- Page 172 and 173: 7. Geology of the Iberian Peninsula
- Page 174 and 175: 7. Geology of the Iberian Peninsula
- Page 176 and 177: 7. Geology of the Iberian Peninsula
- Page 178 and 179: 7. Geology of the Iberian Peninsula
- Page 180 and 181: 7. Geology of the Iberian Peninsula
- Page 182 and 183: 7. Geology of the Iberian Peninsula
- Page 184 and 185: 7. Geology of the Iberian Peninsula
- Page 186 and 187: 7. Geology of the Iberian Peninsula
- Page 188 and 189: 7. Geology of the Iberian Peninsula
- Page 190 and 191: 7. Geology of the Iberian Peninsula
- Page 192 and 193: 7. Geology of the Iberian Peninsula
- Page 194 and 195: 7. Geology of the Iberian Peninsula
- Page 196 and 197: 7. Geology of the Iberian Peninsula
- Page 198 and 199: 7. Geology of the Iberian Peninsula
- Page 200 and 201: 7. Geology of the Iberian Peninsula
- Page 202 and 203: 7. Geology of the Iberian Peninsula
- Page 205 and 206: Recovering a synthetic 3D subsurfac
- Page 207 and 208: direction direction Depth: 12 - 30
- Page 209: 8.2. Generating synthetic 3D model
- Page 213 and 214: 3D N45W 3D-crust TE Rho TE Phi Peri
- Page 215 and 216: 8.3. Inversion of 3D model data sch
- Page 217 and 218: Model variation RMS misfit Optimal
- Page 219 and 220: Profile: 3D-crust (TM-only) Depth (
- Page 221 and 222: Parameter Value 8.3. Inversion of 3
- Page 223 and 224: Depth (km) 10 -2 10 -1 10 0 10 1 10
- Page 225 and 226: Depth (km) 10 -2 10 -1 10 0 10 1 10
- Page 227 and 228: Step 1: Isotropic 2D inversion Step
- Page 229 and 230: 8.3. Inversion of 3D model data par
- Page 231 and 232: 8.4. Summary and conclusions bution
- Page 233 and 234: Regularisation order Smoothing para
- Page 235 and 236: S N 1% 0 Depth (km) 3% Depth (km) 1
- Page 237 and 238: 9.1. Profile location Data collecti
- Page 239 and 240: Location (degrees) Recording period
- Page 241 and 242: Geological region Stations Geologic
- Page 243 and 244: 9.4. Segregation of data acquired w
- Page 245 and 246: Phase (degrees) 135 90 45 0 Z xy -4
- Page 247 and 248: 0 km 10 km 30 km 100 km 300 km Dept
- Page 249 and 250: 0 km 10 km 30 km 100 km 300 km Dept
- Page 251 and 252: 0 km 10 km 30 km 100 km 300 km Dept
- Page 253 and 254: 0 km 10 km 30 km 100 km 300 km Dept
- Page 255 and 256: Pseudo-sections crustal strike dire
- Page 257 and 258: 9.8. Analysis of vertical magnetic
- Page 259: 9.8. Analysis of vertical magnetic
Distance from the centre of the mesh (km)<br />
200<br />
100<br />
0<br />
−100<br />
−200<br />
3D−crust<br />
3D−crust−west<br />
3D−crust<br />
3D−crust−east<br />
8.2. Generating synthetic 3D model data<br />
−200 −100 0 100 200<br />
Distance from the centre of the mesh (km)<br />
Fig. 8.5.: Location of profiles on top of the synthetic 3D model (cf. Figs. 8.3 8.4); North is located to the top of the figure. Stations<br />
are associated with profiles of the same colour; for the profile labelled ‘3D-crust’ three different sets of stations are used: red (3Dcrust,<br />
using stations representing the PICASSO Phase I recording sites), blue (3D-crust-NS), and green (3D-crust-EW). Every profile<br />
contains at least one station on top of each of the four electric conductivity regions (two at crustal depth, two at mantle depth), to<br />
assure that data of every profile is affected by the oblique strike directions. Profiles are constructed orthogonal to the geoelectric strike<br />
direction of either crust or mantle and station locations are projected onto the profile. Background colours indicate the conductivity<br />
distribution at crustal depth (cf. Fig. 8.4).<br />
data.<br />
At shorter periods (0.01 s, top-left figures), responses are dominated by characteristics<br />
of crustal structures, therefore values for the two modes are alike and ρa = ρ; hence<br />
ρa = 50 Ωm (appr. 1.7 in log-scale) in the northeastern half and ρa = 200 Ωm (appr.<br />
2.3 in log-scale) in the southwestern half. For longer periods (100 s, top-right figures)<br />
mantle structures start to add observable contributions to the response data. Values of<br />
ρa at 100 s are similar to values of ρa at 0.01 s (note the different colour scale), but φ<br />
differs significantly. In general, φ(100 s) < φ(0.01 s) ≈ 45 degrees owing to the more<br />
resistive nature of mantle regions in respect to the crust. At 100 s periods, skin depth (cf.<br />
Sec. 3.3) for stations on the resistive side of the crustal fault is appr. 70 km, whereas<br />
on the conductive side it is appr. 35 km. Thus, at 100 s periods all sites are sensitive to<br />
electric properties of the mantle. Phase anomalies at 100 s exhibit a point symmetry in<br />
regards to the centre of the station array, i.e. where crust and mantle interfaces intersect.<br />
Phase anomalies are a superposition of effects from crustal and mantle structures with<br />
175