Using Empirical Mode Decomposition (EMD) for the processing of marine MT data

Magnetotelluric (MT) method determines a frequency dependent impedance tensor using the spectra of associated time-varying horizontal electric and magnetic fields measured at the Earth’s surface. In this abstract, we present a dynamic time series analysis method dealing the non-stationary MT data to infer the impedance tensor. Most current methods to determine the spectra use Fourier transform based procedure and, therefore, assume that the signals are stationary over the record length. We introduce a new method for dealing with non-stationarity of the MT time series based upon empirical mode decomposition (EMD) method, a dynamic time series analysis method. Using EMD complicated data sets can be decomposed into a finite and small number of "intrinsic mode functions" (IMFs), which are mono-component signals and allow the calculation of physical meaningful instantaneous frequencies. EMD has no bias due to non-stationary of geomagnetic time series, since the IMFs are based entirely on signal characteristics and not on any given set of base functions such as sines and cosines in the Fourier transform or wavelets in the Wavelet transform. We use the EMD method to decompose MT data into IMFs and calculate the instantaneous frequencies and spectra to determine the impedance tensor. The method is tested in synthetic and real marine MT data sets, the obtained estimate results are reliable compared to frequently-used BIRRP processing method. Furthermore, new method has the possibility of noise visualization and filtering, which is especially important in marine applications, where noise free time segments maybe short

Cyclic volcanism at convergent margins: linked to aarth orbital parameters or climate changes?

EGU2010-13373 The frequency of volcanic activity varies on a wide rangeof spatial and temporal scales, from <1 yr. periodicities in single volcanic systems to periodicities of 106 yrs. in global volcanism. The causes of these periodicities are poorly understood although the long-term global variations are likely linked to plate-tectonic processes. Here we present evidence for temporal changes in eruption frequencies at an intermediate time scale (104 yrs.) using the Pleistocene to recent records of widespread tephras of sub-Plinian to Plinian, and occasionally co-ignimbrite origin, along the Pacific Ring of Fire, which accounts for about half of the global length of 44,000 km of active subduction. Eruptions at arc volcanoes tend to be highly explosive and the well-preserved tephra records from the ocean floor can be assumed to be representative of how eruption frequencies varied with time. Volcanic activity along the Pacific Ring of Fire evolved through alternating phases of high and low frequency; although there is modulation by local and regional geologic conditions, these variations have a statistically significant periodicity of 43 ka that overlaps with the temporal variation in the obliquity of the Earth’s rotation axis, an orbital parameter that also exerts a strong control on global climate changes. This may suggest that the frequency of volcanic activity is controlled by effects of global climate changes. However, the strongest physical effects of climate change occur at 100 ka periods which are not seen in the volcanic record. We therefore propose that the frequency of volcanic activity is directly influenced by minute changes in the tidal forces induced by the varying obliquity resulting in long-period gravitational disturbances acting on the upper mantle

Adaption and GPU based parallelization of the code TEMDDD for the 3D modelling of CSEM data

The finite difference time domain code TEMDDD was modified for the 3D forward modeling of marine CSEM data. After changes in the code, which make it possible to create model geometries typically encountered in marine CSEM experiments, parts of the code have been parallelized using massive parallelization on graphic cards. Parts of the singular value decomposition, which is the most time consuming part of the code, have been successfully ported with massive speed-ups (8-12x faster) observed as compared to the standard code. The full parallelization of the code is still work in progress

Electrical resistivity image of the South Atlantic continental margin derived from onshore and offshore magnetotelluric data

We present a deep electrical resistivity image from the passive continental margin in Namibia. The approximately 700 km long magnetotelluric profile follows the Walvis Ridge offshore, continues onshore across the Kaoko Mobile Belt and reaches onto the Congo Craton. Two-dimensional inversion reveals moderately resistive material offshore, atypically low for oceanic lithosphere, reaching depths of 15–20 km. Such moderate resistivities are consistent with seismic P wave velocity models, which suggest up to 35 km thick crust. The Neoproterozoic rocks of the Kaoko Mobile Belt are resistive, but NNW-striking major shear-zones are imaged as subvertical, conductive structures in the upper and middle crust. Since the geophysical imprint of the shear zones is intact, opening of the South Atlantic in the Cretaceous did not alter the middle crust. The transition into the cratonic region coincides with a deepening of the high-resistive material to depths of more than 60 km

Controlled-source electromagnetic and seismic delineation of sub-seafloor fluid flow structures in a gas hydrate province, offshore Norway

Deep sea pockmarks underlain by chimney-like or pipe structures that contain methane hydrate are abundant along the Norwegian continental margin. In such hydrate provinces the interaction between hydrate formation and fluid flow has significance for benthic ecosystems and possibly climate change. The Nyegga region, situated on the western Norwegian continental slope, is characterized by an extensive pockmark field known to accommodate substantial methane gas hydrate deposits. The aim of this study is to detect and delineate both the gas hydrate and free gas reservoirs at one of Nyegga's pockmarks. In 2012, a marine controlled-source electromagnetic (CSEM) survey was performed at a pockmark in this region, where high-resolution three-dimensional seismic data were previously collected in 2006. Two-dimensional CSEM inversions were computed using the data acquired by ocean bottom electrical field receivers. Our results, derived from unconstrained and seismically constrained CSEM inversions, suggest the presence of two distinctive resistivity anomalies beneath the pockmark: a shallow vertical anomaly at the underlying pipe structure, likely due to gas hydrate accumulation, and a laterally extensive anomaly attributed to a free gas zone below the base of the gas hydrate stability zone. This work contributes to a robust characterization of gas hydrate deposits within sub-seafloor fluid flow pipe structures

High-resolution resistivity imaging of marine gas hydrate structures by combined inversion of CSEM towed and ocean-bottom receiver data

We present high-resolution resistivity imaging of gas hydrate pipe-like structures, as derived from marine controlled-source electromagnetic (CSEM) inversions that combine towed and ocean-bottom electric field receiver data, acquired from the Nyegga region, offshore Norway. 2.5-D CSEM inversions applied to the towed receiver data detected four new prominent vertical resistive features that are likely gas hydrate structures, located in proximity to a major gas hydrate pipe-like structure, known as the CNE03 pockmark. The resistivity model resulting from the CSEM data inversion resolved the CNE03 hydrate structure in high resolution, as inferred by comparison to seismically constrained inversions. Our results indicate that shallow gas hydrate vertical features can be delineated effectively by inverting both ocean-bottom and towed receiver CSEM data simultaneously. The approach applied here can be utilized to map and monitor seafloor mineralization, freshwater reservoirs, CO2 sequestration sites and near-surface geothermal systems

Reducing motion noise in marine magnetotelluric measurements by means of tilt records

The analysis of marine magnetotelluric data is often complicated by disturbing signals that are caused by small-scale periodic movements of the instrument. The motion-induced noise leads to a bias and/or severe scattering in the derived magnetotelluric transfer functions. Both the motion itself and its effects on the magnetic and telluric time-series are investigated in this study using an 80 d magnetotelluric data set that includes dynamic tilt records measured in the Pacific Ocean off Costa Rica. We apply a standard motion removal technique as well as a newly developed method to correct for motion-induced noise. The resulting magnetotelluric transfer functions are of significantly better quality than the uncorrected ones. Furthermore, the study of the properties of motion noise leads to conclusions about the optimal processing approach even in case of data sets where an explicit correction for that noise is not possible