Arc-continent collisions are a fundamental plate tectonic process that control continental growth, orogen development and the distribution of natural hazards and natural resources. This process actively occurs in Taiwan where the Luzon Volcanic Arc collides obliquely with the passive margin of the Eurasian Plate. An important characteristic of the arc-continent collision in Taiwan is that oblique convergence has produced an orogen that decreases in age from north to south. Investigation of the temporal evolution of the Taiwan orogen is therefore made possible through studies at different latitudes on the island. The first long-period magnetotelluric (MT) measurements in Taiwan were recorded at 82 locations on three profiles across south, central and north Taiwan during 2006-2007. These MT data were collected as a component of the TAIGER (Taiwan Integrated Geodynamics Research) project and are analyzed and interpreted in this thesis. The TAIGER project was initiated to: 1) collect multi-technique high resolution geophysical data required to unambiguously distinguish between end-member tectonic models proposed for central Taiwan, and 2) to investigate the orogen evolution by comparison of data collected at latitudes ranging from south to north Taiwan. In this thesis, the central Taiwan TAIGER MT data are shown to be inconsistent with orogen-scale thin-skinned tectonic models and instead support predictions of lithospheric-scale deformation (i.e. thick-skinned tectonics) beneath the Central Ranges. Similarity between resistivity models of central and southern Taiwan indicate that thick-skinned tectonics is occurring in both locations, and is therefore a persistent mode of deformation that operates as the orogen develops. The resistivity model for northern Taiwan is shown to be consistent with dewatering of the subducting Philippine slab, and with deformation described by the subducting indenter tectonic model. A global context for the Taiwan results is provided by comparison to the South Island of New Zealand that shares many similarities with the tectonic setting in Taiwan. MT resistivity models support lithospheric deformation beneath the Southern Alps, and the occurrence of fluid-based conductors in the Marlborough region caused by dewatering of the Pacific slab. / Geophysics

No description available.

As part of the southern Ontario POLARIS project, this thesis uses magnetotelluric methods to investigate the lithospheric architecture of the Proterozoic Grenville Province and its margin with the Archean Superior Province. The first multi-dimensional crustal and lithospheric resistivity images for this region are presented. The resistivity structure of the Phanerozoic sedimentary rocks in the lower Great Lakes region was determined using 1-D methods. The responses are strongly affected by a 20-23 S conductive layer within the sedimentary rocks, interpreted to be associated with Upper Ordovician shale units. This layer excludes resolution of resistivity structure of underlying crust. The resistivity structure of the Precambrian crust and lithosphere was determined using 2-D methods. Different strike azimuths were determined for the crust, the upper lithospheric mantle and the deeper mantle layer. The crustal resistivity model for a profile from 50oN79oW to 43oN76oW images resistive Laurentian margin rocks dipping southeast to the base of the crust, bounded by the Grenville Front and the Central Metasedimentary Belt Boundary Zone. In a 2-D model of the mantle lithosphere for the same profile, a conductor at 70-150 km depth, located along-strike from the Mesozoic Kirkland Lake and Cobalt kimberlite fields, is interpreted to be due to mantle re-fertilization. Results from multiple MT profiles indicate conductive (<10 Ω.m) lithospheric mantle beneath the Central Metasedimentary Belt and show that the northwestern Grenville Province is characterized by large-scale, resistive lithosphere (>10,000 Ω.m) extending for about 300 km beneath the Grenville Province and 800 km along strike. Lithospheric thickness is interpreted to be 280 km; local decreases in this depth are attributed to refertilization of the lower mantle lithosphere by fluids associated with Cretaceous kimberlite magmatism. Anisotropic 2-D modeling reveals minimal electrical anisotropy (<10%) at mantle depths in contrast to the factor of 15 anisotropy determined in earlier 1-D studies. This result suggests that observed MT response anisotropy is caused by large-scale structures. Strike direction in the upper lithospheric mantle is interpreted to be related to the Archean fabric of the Superior craton and in the deeper, conductive, mantle it is interpreted to have been established in the Cretaceous.

In this thesis, four studies using different geophysical electromagnetic methods are presented. In the first study dealing with airborne measurements, the noise response due to the rotation of the aircraft and the aircraft itself as a metallic conductive body on the Earth's electromagnetic response in very low frequency and low frequency band was investigated. The magnetic fields are independent of the aircraft in the VLF band and part of the LF band. But at higher frequencies (above 100 kHz), the signals are more influenced by the aircraft. The aircraft also generates its own noise frequencies which are mixed with the radio transmitter signals. The second and third studies are applications of radio-, controlled source-magnetotellurics and electrical resistivity tomography methods at a quick-clay landslide site in southwest Sweden. The data are processed and modeled in 2D and 3D, and the models are compared with high-resolution seismic and geotechnical data. The obtained results were further validated and refined by performing synthetic tests in the second study. The third study shows that the 3D models provide larger and more continuous volume of the quick clay structure than traditional 2D models. Both studies have shown that integrated application of geophysical methods for landslides is ideal. Quick clays often overlie the coarse-grained layers showing an increase of resistivity values in the models. In the fourth study, a new audio magnetotelluric data acquisition technique is developed and is named moving magnetotellurics (MMT). In this new technique, the magnetic sensors are placed on the ground and only 15 to 20 minutes data are acquired for each station, which usually is enough to cover the frequency range 30-300 Hz. The new technique is more efficient and convenient than the traditional magnetotelluric method, and test measurements have shown that it is an applicable method in shallow depth studies.

No description available.

In this work two ways of processing controlled-source magnetotelluric (MT) data were tried and compared against each other. The aim was to evaluate the differences between a multivariate-based processing method to a bivariate processing method. The software KMSProTF represents conventional processing using bivariate and robust statistics. MsDEMPCA utilizes a multi-variate Criss-Cross regression scheme to improve the condition of the data-matrix before robustly decomposing it into principal components. Data from the FENICS-19 survey in northern Scandinavia was processed to transfer functions (TF) using the respective method. The TFs were visually interpreted in KMSProTF. There were no significant differences found between the methods. In addition a calibration between instruments was carried out, which caused an exclusion of parts of the data-set.

The radio magnetotellucic (RMT) method originates from (Tikhonov, 1950) and (Cagniard,1953) who discovered the possibility to estimate resistivity in the Earth’s subsurface using naturally occurring telluric currents. This through measuring the surface impedance of electromagnetic waves. Later (Paal, 1965) makes use of electromagnetic plane-waves originating from radio transmitters that operates in the very low frequency range (10 − 30kHz). These plane waves artificially induce telluric currents which can then be used to estimate resistivity in shallow subsurface. It is hence the source of the electromagnetic wave inducing the telluric current that named the RMT method accordingly. It’s then shown by (Bastani, 2001) that the signal can be measured in an even broader band (10 − 250kHz) to gain a better vertical resolution. Finally, frequency of the electromagnetic wave is affecting the skin depth and one can therefore estimate the depth of resistivity layers in the subsurface. In this thesis a data processing algorithm has been developed, based on the work in (Bastani, 2001), that processes RMT field data. Impedance, resistivity among other quantities is calculated. One can then export this data and use it in existing modeling software. New processing parameters has been implemented and its effect on the data set investigated. The software is implemented in matlab, and tested with synthetic data and data measured at Blötberget, Sweden. However, work is still to be done due to field equipment malfunction during measurements at Blötberget. This rendering an, at least partly, faulty data set. Therefore one cannot completely exclude calibration issues until a clean data set is taken. / Den radio magnetotelluriska (RMT) metoden har sitt ursprung i (Tikhonov, 1950) och (Cagniard, 1953) arbete där strukturen hos geologisk resistivitet uppskattas med hjälp av naturligt förekommande elektriska jordströmmar. Detta genom att mäta den elektromagnetiska vågimpedansen vid markytan. Sedan använder (Paal, 1965) sig av elektromagnetiska planvågor från lågfrekventa radiosändare (10 − 30kHz), som artificiellt inducerar jordströmmar, för att uppskatta resistiviteten i mer ytligt underliggande geologisk struktur. Det är således källan till den inducerade jordströmmen som namngett RMT-metoden. Det är sedan påvisat av (Bastani,2001) att man kan använda sig av frekvenser i ett bredare spektrum (10 − 250kHz) för att få ökad vertikal upplösning. Slutligen påverkar frekvensen hos den elektromagnetiska vågen inträngningsdjupet och på så sätt tillåts det att uppskatta resistiviteten vid ett särskilt djup. I denna studie har en databehandlingsalgorithm utvecklats, baserat på arbete utfört av (Bastani, 2001), som bearbetar RMT fältdata. Impedans och resistivitet samt andra storheter beräknas för att sedan kunna exporteras och användas i existerande modelleringsprogramvara. Programvaran är skriven i matlab, och dess funktionalitet har prövats med syntetisk data samt data från Blötberget, Sverige. Även påverkan av olika processparametrar har undersökts. Detta har gjorts med framgång, men arbete återstår för att fullständigt verifiera korrekt implementering då utrustningen som användes för insamling av fältdata vid Blötberget felade. Därför kan kalibreringsproblem i programvaran inte helt uteslutas.

There are very few software packages for model-building and visualization in electromagnetic geophysics, particularly when compared to other geophysical disciplines, such as seismology. The purpose of this thesis is to design, develop, and test a geophysical model-building interface that allows users to parameterize the 2D magnetotellurics problem. Through the evaluation of this interface, feedback was collected from a usability specialist and a group of geophysics graduate students to study the steps users take to work through the 2D forward-modeling problem, and to analyze usability errors encountered while working with the interface to gain a better understanding of how to build a more effective interface. Similar work has been conducted on interface design in other fields, such as medicine and consumer websites. Usability Engineering is the application of a systematic set of methods to the design and development of software with the goal of making the software more learnable, easy to use, and accessible. Two different Usability Engineering techniques — Heuristic Evaluation and Thinking Aloud Protocol — were involved in the evaluation of the interface designed in this study (FEM2DGUI). Heuristic Evaluation is a usability inspection method that employs a usability specialist to detect errors based on a known set of guidelines and personal experience. Thinking Aloud Protocol is a usability evaluation method where potential end-users are observed as they verbalize their every step as they work through specific scenarios with an interface. These Usability Engineering methods were combined in a effort to understand how the first prototype of FEM2DGUI could be refined to make it more usable and to understand how end-users work through the forward-modeling problem. The Usability Engineering methods employed in this project uncovered multiple usability errors that were corrected through a refinement of the interface. Discovery of these errors helped with refining the system to become more robust and usable, which is believed to aid users in more efficient model-building. Because geophysical model-building is inherently a difficult task, it is possible that other model-building graphical user interfaces could benefit from the application of Usability Engineering methods, such as those presented in this research.â / Master of Science

京都大学 / 0048 / 新制・課程博士 / 博士(理学) / 甲第20921号 / 理博第4373号 / 新制||理||1628(附属図書館) / 京都大学大学院理学研究科地球惑星科学専攻 / (主査)教授 大志万 直人, 准教授 吉村 令慧, 教授 中西 一郎 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DGAM

This thesis presents two new algorithms for doing constrained Magnetotelluric (MT) inversion based on an existing Occam 2D inversion program. The first algorithm includes borehole resistivity logs as prior information to constrain resistivity directly in the vicinity of boreholes. The second algorithm uses reflection seismic data as prior constraints to transfer structural information from seismic images to 2D resistivity models. These two algorithms are efficient (proved through tests of synthetic examples) and widely applicable. In this thesis, they have been successfully applied to the COSC (Collisional Orogeny in the Scandinavian Caledonides) MT data. The COSC project aims to study the mountain belt dynamics in central Sweden by drilling two 2.5 km deep boreholes. MT data were collected to locate the main décollement that separates the overlying Caledonian allochthons and the underlying Precambrian basement, as the main décollement is associated with very conductive Alum shale. The previous interpretation based on part of the COSC seismic profile (CSP) was that the main décollement was located along a reflection with depth of 4.5 km underneath Åre and ~3 km underneath Mörsil, in central Jämtland. The MT resistivity model reveals a very conductive layer in the central and western parts of the profile, the top of which coincides with the first seismic reflection. This means that the first conductive alum shale layer occurs at less than 1 km depth, supporting a new interpretation of the main décollement at shallower depth. In a re-interpretation of the CSP data based on the MT model, the main décollement occurs a few hundred metres below the top of the conductor and is coincident with a laterally continuous seismic reflection. Further, the overlying seismic reflections resemble imbricated alum shale of the Lower Allochthon. MT inversion using seismic constraints from CSP gives further support to the new interpretation. Moreover, MT investigations were conducted in the Alnö alkaline and carbonatite ring-intrusion complex in Sweden. 2D and 3D resistivity models inverted from MT data together with resistivity and porosity laboratory measurements delineate a fossil magma chamber as a resistive anomaly surrounded by electrically conductive up-doming and ring-shaped faults and fractures.