scholarly journals Upper mantle structure of Mars from InSight seismic data

Science ◽  
2021 ◽  
Vol 373 (6553) ◽  
pp. 434-438 ◽  
Author(s):  
Amir Khan ◽  
Savas Ceylan ◽  
Martin van Driel ◽  
Domenico Giardini ◽  
Philippe Lognonné ◽  
...  

For 2 years, the InSight lander has been recording seismic data on Mars that are vital to constrain the structure and thermochemical state of the planet. We used observations of direct (P and S) and surface-reflected (PP, PPP, SS, and SSS) body-wave phases from eight low-frequency marsquakes to constrain the interior structure to a depth of 800 kilometers. We found a structure compatible with a low-velocity zone associated with a thermal lithosphere much thicker than on Earth that is possibly related to a weak S-wave shadow zone at teleseismic distances. By combining the seismic constraints with geodynamic models, we predict that, relative to the primitive mantle, the crust is more enriched in heat-producing elements by a factor of 13 to 20. This enrichment is greater than suggested by gamma-ray surface mapping and has a moderate-to-elevated surface heat flow.

2013 ◽  
Vol 31 (4) ◽  
pp. 619 ◽  
Author(s):  
Luiz Eduardo Soares Ferreira ◽  
Milton José Porsani ◽  
Michelângelo G. Da Silva ◽  
Giovani Lopes Vasconcelos

ABSTRACT. Seismic processing aims to provide an adequate image of the subsurface geology. During seismic processing, the filtering of signals considered noise is of utmost importance. Among these signals is the surface rolling noise, better known as ground-roll. Ground-roll occurs mainly in land seismic data, masking reflections, and this roll has the following main features: high amplitude, low frequency and low speed. The attenuation of this noise is generally performed through so-called conventional methods using 1-D or 2-D frequency filters in the fk domain. This study uses the empirical mode decomposition (EMD) method for ground-roll attenuation. The EMD method was implemented in the programming language FORTRAN 90 and applied in the time and frequency domains. The application of this method to the processing of land seismic line 204-RL-247 in Tacutu Basin resulted in stacked seismic sections that were of similar or sometimes better quality compared with those obtained using the fk and high-pass filtering methods.Keywords: seismic processing, empirical mode decomposition, seismic data filtering, ground-roll. RESUMO. O processamento sísmico tem como principal objetivo fornecer uma imagem adequada da geologia da subsuperfície. Nas etapas do processamento sísmico a filtragem de sinais considerados como ruídos é de fundamental importância. Dentre esses ruídos encontramos o ruído de rolamento superficial, mais conhecido como ground-roll . O ground-roll ocorre principalmente em dados sísmicos terrestres, mascarando as reflexões e possui como principais características: alta amplitude, baixa frequência e baixa velocidade. A atenuação desse ruído é geralmente realizada através de métodos de filtragem ditos convencionais, que utilizam filtros de frequência 1D ou filtro 2D no domínio fk. Este trabalho utiliza o método de Decomposição em Modos Empíricos (DME) para a atenuação do ground-roll. O método DME foi implementado em linguagem de programação FORTRAN 90, e foi aplicado no domínio do tempo e da frequência. Sua aplicação no processamento da linha sísmica terrestre 204-RL-247 da Bacia do Tacutu gerou como resultados, seções sísmicas empilhadas de qualidade semelhante e por vezes melhor, quando comparadas as obtidas com os métodos de filtragem fk e passa-alta.Palavras-chave: processamento sísmico, decomposição em modos empíricos, filtragem dados sísmicos, atenuação do ground-roll.


2021 ◽  
Author(s):  
Adam Cygal ◽  
Michał Stefaniuk ◽  
Anna Kret

AbstractThis article presents the results of an integrated interpretation of measurements made using Audio-Magnetotellurics and Seismic Reflection geophysical methods. The obtained results were used to build an integrated geophysical model of shallow subsurface cover consisting of Cenozoic deposits, which then formed the basis for a detailed lithological and tectonic interpretation of deeper Mesozoic sediments. Such shallow covers, consisting mainly of glacial Pleistocene deposits, are typical for central and northern Poland. This investigation concentrated on delineating the accurate geometry of Obrzycko Cenozoic graben structure filled with loose deposits, as it was of great importance to the acquisition, processing and interpretation of seismic data that was to reveal the tectonic structure of the Cretaceous and Jurassic sediments which underly the study area. Previously, some problems with estimation of seismic static corrections over similar grabens filled with more recent, low-velocity deposits were encountered. Therefore, a novel approach to estimating the exact thickness of such shallow cover consisting of low-velocity deposits was applied in the presented investigation. The study shows that some alternative geophysical data sets (such as magnetotellurics) can be used to significantly improve the imaging of geological structure in areas where seismic data are very distorted or too noisy to be used alone


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Gaochun Wang ◽  
Hans Thybo ◽  
Irina M. Artemieva

AbstractAll models of the magmatic and plate tectonic processes that create continental crust predict the presence of a mafic lower crust. Earlier proposed crustal doubling in Tibet and the Himalayas by underthrusting of the Indian plate requires the presence of a mafic layer with high seismic P-wave velocity (Vp > 7.0 km/s) above the Moho. Our new seismic data demonstrates that some of the thickest crust on Earth in the middle Lhasa Terrane has exceptionally low velocity (Vp < 6.7 km/s) throughout the whole 80 km thick crust. Observed deep crustal earthquakes throughout the crustal column and thick lithosphere from seismic tomography imply low temperature crust. Therefore, the whole crust must consist of felsic rocks as any mafic layer would have high velocity unless the temperature of the crust were high. Our results form basis for alternative models for the formation of extremely thick juvenile crust with predominantly felsic composition in continental collision zones.


2019 ◽  
Vol 16 (6) ◽  
pp. 1017-1031 ◽  
Author(s):  
Yong Hu ◽  
Liguo Han ◽  
Rushan Wu ◽  
Yongzhong Xu

Abstract Full Waveform Inversion (FWI) is based on the least squares algorithm to minimize the difference between the synthetic and observed data, which is a promising technique for high-resolution velocity inversion. However, the FWI method is characterized by strong model dependence, because the ultra-low-frequency components in the field seismic data are usually not available. In this work, to reduce the model dependence of the FWI method, we introduce a Weighted Local Correlation-phase based FWI method (WLCFWI), which emphasizes the correlation phase between the synthetic and observed data in the time-frequency domain. The local correlation-phase misfit function combines the advantages of phase and normalized correlation function, and has an enormous potential for reducing the model dependence and improving FWI results. Besides, in the correlation-phase misfit function, the amplitude information is treated as a weighting factor, which emphasizes the phase similarity between synthetic and observed data. Numerical examples and the analysis of the misfit function show that the WLCFWI method has a strong ability to reduce model dependence, even if the seismic data are devoid of low-frequency components and contain strong Gaussian noise.


Geophysics ◽  
2006 ◽  
Vol 71 (3) ◽  
pp. V79-V86 ◽  
Author(s):  
Hakan Karsli ◽  
Derman Dondurur ◽  
Günay Çifçi

Time-dependent amplitude and phase information of stacked seismic data are processed independently using complex trace analysis in order to facilitate interpretation by improving resolution and decreasing random noise. We represent seismic traces using their envelopes and instantaneous phases obtained by the Hilbert transform. The proposed method reduces the amplitudes of the low-frequency components of the envelope, while preserving the phase information. Several tests are performed in order to investigate the behavior of the present method for resolution improvement and noise suppression. Applications on both 1D and 2D synthetic data show that the method is capable of reducing the amplitudes and temporal widths of the side lobes of the input wavelets, and hence, the spectral bandwidth of the input seismic data is enhanced, resulting in an improvement in the signal-to-noise ratio. The bright-spot anomalies observed on the stacked sections become clearer because the output seismic traces have a simplified appearance allowing an easier data interpretation. We recommend applying this simple signal processing for signal enhancement prior to interpretation, especially for single channel and low-fold seismic data.


1983 ◽  
Vol 73 (6A) ◽  
pp. 1499-1511
Author(s):  
Paul Silver

Abstract A method is proposed for retrieving source-extent parameters from far-field body-wave data. At low frequency, the normalized P- or S-wave displacement amplitude spectrum can be approximated by |Ω^(r^,ω)| = 1 − τ2(r^)ω2/2 where r^ specifies a point on the focal sphere. For planar dislocation sources, τ2(r^) is linearly related to statistical measures of source dimension, source duration, and directivity. τ2(r^) can be measured as the curvature of |Ω^(r^,ω)| at ω = 0 or the variance of the pulse Ω^(r^,t). The quantity ωc=2τ−1(r^) is contrasted with the traditional corner frequency ω0, defined as the frequency at the intersection of the low- and high-frequency trends of |Ω^(r^,ω)|. For dislocation models without directivity, ωc(P) ≧ ωc(S) for any r^. A mean corner frequency defined by averaging τ2(r^) over the focal sphere, ω¯c=2&lt;τ2(r^)&gt;−1/2, satisfies ωc(P) &gt; ωc(S) for any dislocation source. This behavior is not shared by ω0. It is shown that ω0 is most sensitive to critical times in the rupture history of the source, whereas ωc is determined by the basic parameters of source extent. Evidence is presented that ωc is the corner frequency measured on actual seismograms. Thus, the commonly observed corner frequency shift (P-wave corner greater than the S-wave corner), now viewed as a shift in ωc is simply a result of spatial finiteness and is expected to be a property of any dislocation source. As a result, the shift cannot be used as a criterion for rejecting particular dislocation models.


Sign in / Sign up

Export Citation Format

Share Document