Backus-Gilbert style inversions for mantle anisotropy using normal mode data 

Author(s):  
Federica Restelli ◽  
Paula Koelemeijer ◽  
Christophe Zaroli

<p>Seismic tomography is essential for imaging the Earth’s interior in order to better understand the dynamic processes at work. However, robust physical interpretation of tomographic images remain difficult as the inverse problem is under-determined, model amplitudes are biased and uncertainties are usually not quantified.</p><p>Commonly-used techniques, such as damped least-square inversions, break the non-uniqueness of the model solution by adding a subjective, ad hoc, regularization, which can lead to biased amplitudes and potential physical misinterpretations. The SOLA method (Zaroli, 2016; Zaroli et al., 2017), based on a Backus-Gilbert approach, removes the non-uniquess by averaging, rather than introducing a subjective regularization. The method explicitly constrains the amplitudes to be unbiased and the computation of the model resolution and uncertainty is inherent and efficient. Instead of aiming to minimize the data fit, the SOLA approach aims to minimize the size of the averaging volume <!-- Think it is clear enough without the extra sentence – as averaging is mentioned before -->and the associated uncertainties.</p><p>We aim to build a new tomographic model of the Earth’s mantle using the SOLA method. We focus our observations on normal mode data, the standing waves of the Earth observed after very large earthquakes, which are not affected by an uneven data distribution. As normal modes are sensitive to multiple seismic parameters, we treat the sensitivity to different parameters as so called “3D noise” within the SOLA framework. We are specifically interested in constraining seismic anisotropy, which provides more direct information on mantle flow.</p><p>Here, we report on some forward modelling results, fundamental to understanding normal mode sensitivity to seismic anisotropy at different depths and identifying which modes to focus on during inversions. We also show our initial work towards building a new tomography model, including the calculation of 3D noise and target kernels.</p><p> </p>

2020 ◽  
Author(s):  
Eugene Humphreys ◽  
Jorge Castellanos ◽  
Robert Clayton ◽  
Jonathan Perry-Houts ◽  
YoungHee Kim ◽  
...  

<p>Azimuthal anisotropy in the NW U.S. crust is derived using 3-17 s Rayleigh waves derived using ambient noise from about 300 broadband stations. Velocity is resolved between all station pairs in close proximity, and velocity as a function of azimuth is determined for each station. Azimuthal anisotropy orientations point strongly toward tomographically-imaged high-velocity structures in the underlying mantle, but show no relation to the underlying mantle anisotropy field. We suggest that the crustal anisotropy is decoupled from lateral tectonic forces and is created by upper mantle vertical loading, which in turn generates lateral pressure gradients that drive channelized flow in the ductile mid and lower crust. This idea is tested with geodynamic modeling. Using reasonable values for crustal viscosity and mantle buoyancy structure, we find that the local buoyancy sources within the upper mantle will drive the viscous crustal flow in a manner that reproduces well the imaged crustal anisotropy. We conclude that mantle vertical loading, acting independently from mantle flow, can actively control crustal deformation on a scale of several hundred kilometers.</p>


We pursue an abstract investigation of the theory of the infinitesimal free elasticgravitational oscillations of a fairly general rotating Earth model. By considering in some detail the transition to the non-rotating case, we are able to delineate certain of the intrinsic effects of rotation on the normal mode eigensolutions, and to show how profoundly rotation alters the fundamental mathematical and physical properties of these eigensolutions. In particular, we show that the displacement eigenfunctions of a rotating Earth model are not mutually orthogonal, and that the corresponding normal modes of oscillation cannot in general be represented by pure standing waves. We consider the excitation of the normal modes of oscillation of a rotating Earth model by a transient imposed body force distribution, and we show that the complex dynamical amplitude of each normal mode may, in many geophysical applications, be determined separately, in spite of the lack of orthogonality among the displacement eigenfunctions. The calculation of the associated static response after the decay of the normal modes of oscillation is, on the other hand, complicated considerably by the absence of orthogonality. We specifically examine the influence of rotation on the zero-frequency rigid body translational and rotational modes of any non-rotating Earth model, and show how to account for the corresponding rigid body modes of any rotating Earth model in excitation calculations.


A gross Earth datum is a single measurable number describing some property of the whole Earth, such as mass, moment of inertia, or the frequency of oscillation of some identified elastic-gravitational normal mode. We suppose that a finite set G of gross Earth data has been measured, that the measurements are inaccurate, and that the variance matrix of the errors of measurement can be estimated. We show that some such sets G of measurements determine the structure of the Earth within certain limits of error except for fine-scale detail. That is, from some setsG it is possible to compute localized averages of the Earth structure at various depths. These localized averages will be slightly in error, and their errors will be larger as their resolving lengths are shortened. We show how to determine whether a given set G of measured gross Earth data permits such a construction of localized averages, and, if so, how to find the shortest length scale over which G gives a local average structure at a particular depth if the variance of the error in computing that local average from G is to be less than a specified amount. We apply the general theory to the linear problem of finding the depth variation of a frequency-independent local elastic dissipation ( Q ) from the observed damping rates of a finite number of normal modes. We also apply the theory to the nonlinear problem of finding density against depth from the total mass, moment and normal-mode frequencies, in case the compressional and shear velocities are known.


Electronics ◽  
2021 ◽  
Vol 10 (5) ◽  
pp. 577
Author(s):  
Luca Schirru ◽  
Tonino Pisanu ◽  
Angelo Podda

Space debris is a term for all human-made objects orbiting the Earth or reentering the atmosphere. The population of space debris is continuously growing and it represents a potential issue for active satellites and spacecraft. New collisions and fragmentation could exponentially increase the amount of debris and so the level of risk represented by these objects. The principal technique used for the debris monitoring, in the Low Earth Orbit (LEO) between 200 km and 2000 km of altitude, is based on radar systems. The BIRALET system represents one of the main Italian radars involved in resident space objects observations. It is a bi-static radar, which operates in the P-band at 410–415 MHz, that uses the Sardinia Radio Telescope as receiver. In this paper, a detailed description of the new ad hoc back-end developed for the BIRALET radar, with the aim to perform slant-range and Doppler shift measurements, is presented. The new system was successfully tested in several validation measurement campaigns, the results of which are reported and discussed.


1998 ◽  
Vol 06 (04) ◽  
pp. 435-452 ◽  
Author(s):  
Robert P. Gilbert ◽  
Zhongyan Lin ◽  
Klaus Hackl

Normal-mode expansions for Green's functions are derived for ocean–bottom systems. The bottom is modeled by Kirchhoff and Reissner–Mindlin plate theories for elastic and poroelastic materials. The resulting eigenvalue problems for the modal parameters are investigated. Normal modes are calculated by Hankel transformation of the underlying equations. Finally, the relation to the inverse problem is outlined.


Author(s):  
S. Y. Chen ◽  
M. S. Ju ◽  
Y. G. Tsuei

Abstract A frequency-domain technique to extract the normal mode from the measurement data for highly coupled structures is developed. The relation between the complex frequency response functions and the normal frequency response functions is derived. An algorithm is developed to calculate the normal modes from the complex frequency response functions. In this algorithm, only the magnitude and phase data at the undamped natural frequencies are utilized to extract the normal mode shapes. In addition, the developed technique is independent of the damping types. It is only dependent on the model of analysis. Two experimental examples are employed to illustrate the applicability of the technique. The effects due to different measurement locations are addressed. The results indicate that this technique can successfully extract the normal modes from the noisy frequency response functions of a highly coupled incomplete system.


Author(s):  
Enbo Fan ◽  
Yumei He ◽  
Yinshuang Ai ◽  
Stephen S. Gao ◽  
Kelly H. Liu ◽  
...  

2021 ◽  
Author(s):  
David Schlaphorst ◽  
Graça Silveira ◽  
João Mata ◽  
Frank Krüger ◽  
Torsten Dahm ◽  
...  

<p>The Madeira and Canary archipelagos, located in the eastern North Atlantic, are two of many examples of hotspot surface expressions, but a better understanding of the crust and upper mantle structure beneath these regions is needed to investigate their structure in more detail. With the study of seismic anisotropy, it is possible to assess the rheology and structure of asthenosphere and lithosphere that can reflect a combination of mantle and crustal contributions.</p><p>Here, as part of the SIGHT project (SeIsmic and Geochemical constraints on the Madeira HoTspot), we present the first detailed study of seismic anisotropy beneath both archipelagos, using data collected from over 60 local three-component seismic land stations. Basing our observations on both teleseismic SKS and local S splitting, we are able to distinguish between multiple layers of anisotropy. We observe significant changes in delay time and fast shear-wave orientation patterns on short length-scales on the order of tens of kilometres beneath the western Canary Islands and Madeira Island. In contrast, the eastern Canary Islands and Porto Santo the pattern is much more uniform. The detected delay time increase and more complex orientation patterns beneath the western Canary Islands and Madeira can be attributed to mantle flow disturbed and diverted on small-length scales by a strong vertical component. This is a clear indication of the existence of a plume at each of those archipelagos, nowadays exerting a strong influence on the western and younger islands. We therefore conclude that a plume-like feature beneath Madeira exists in a similar way to the Canary Island hotspot and that regional mantle flow models for the region should be reassessed.</p><p>This is a contribution to project SIGHT (Ref. PTDC/CTA-GEF/30264/2017). The authors would like to acknowledge the financial support FCT through project UIDB/50019/2020 – IDL.</p>


Sign in / Sign up

Export Citation Format

Share Document