| 1. |
- Abrehdary, Majid, et al.
(author)
-
The spherical terrain correction and its effect on the gravimetric-isostatic Moho determination
- 2016
-
In: Geophysical Journal International. - : Oxford University Press. - 0956-540X .- 1365-246X .- 1687-885X .- 1687-8868. ; 204:1, s. 262-273
-
Journal article (peer-reviewed)abstract
- In this study, the Moho depth is estimated based on the refined spherical Bouguer gravity disturbance and DTM2006 topographic data using the Vening Meinesz-Moritz gravimetric-isostatic hypothesis. In this context, we compute the refined spherical Bouguer gravity disturbances in a set of 1 degrees x 1 degrees blocks. The spherical terrain correction, a residual correction to each Bouguer shell, is computed using rock heights and ice sheet thicknesses from the DTM2006 and Earth2014 models. The study illustrates that the defined simple Bouguer gravity disturbance corrected for the density variations of the oceans, ice sheets and sediment basins and also the non-isostatic effects needs a significant terrain correction to become the refined Bouguer gravity disturbance, and that the isostatic gravity disturbance is significantly better defined by the latter disturbance plus a compensation attraction. Our study shows that despite the fact that the lateral variation of the crustal depth is rather smooth, the terrain affects the result most significantly in many areas. The global numerical results show that the estimated Moho depths by the simple and refined spherical Bouguer gravity disturbances and the seismic CRUST1.0 model agree to 5.6 and 2.7 km in RMS, respectively. Also, the mean value differences are 1.7 and 0.2 km, respectively. Two regional numerical studies show that the RMS differences between the Moho depths estimated based on the simple and refined spherical Bouguer gravity disturbance and that using CRUST1.0 model yield fits of 4.9 and 3.2 km in South America and yield 3.2 and 3.4 km in Fennoscandia, respectively.
|
|
| 2. |
- Abreu, Rafael, et al.
(author)
-
Micropolar modelling of rotational waves in seismology
- 2017
-
In: Geophysical Journal International. - : OXFORD UNIV PRESS. - 0956-540X .- 1365-246X. ; 210:2, s. 1021-1046
-
Journal article (peer-reviewed)abstract
- In this contribution we study elastic wave propagation via the introduction of the micropolar theory. As a generalization of a classical linear elastic medium, a micropolar medium allows each particle to have intrinsic rotational degrees of freedom (spin). We perform numerical experiments using the Pseudospectral method. We find analytical harmonic micropolar solutions for different problem configurations, which result in waveform differences between the classical linear elastic and micropolar media. In contrast to linear elastic media, wave propagation in micropolar media is dispersive. We study how the spin waveform depends on the micropolar elastic parameters and frequency content of the simulation. The micropolar effect on numerical seismograms has a direct implication on the phase, amplitude and arrival time. For frequencies lower than the cut-off frequency, the spin waveform has the same amplitude as the macrorotation field. For frequencies higher than the cut-off frequency, the amplitude of the spin waveform decreases with increasing frequency, so that then it is no longer comparable to the amplitude of macroscopic rotations. When both frequencies are equal there is no wave propagation. This work attempts to clarify the theory of micropolar media for its applications in seismology. We argue that micropolar theory should be further investigated for its potential uses in seismology to, for example, describe energy dissipation, seismograms recorded with rotational seismometers and rupture processes.
|
|
| 3. |
- Abril, Claudia, et al.
(author)
-
Relocating earthquakes with empirical traveltimes
- 2018
-
In: Geophysical Journal International. - : Oxford University Press. - 0956-540X .- 1365-246X. ; 214:3, s. 2098-2114
-
Journal article (peer-reviewed)abstract
- A strategy is proposed to incorporate effects of 3-D velocity variations on earthquake locationsusing empirical traveltimes (ETTs). Traveltime residuals are interpolated from those predictedusing a 1-D velocity reference, mapped on to the hypocentres of corresponding earthquakesfor each station in a network. First, station corrections are computed by averaging the residualsover a fixed scale. Then, summary-ray variograms are used to estimate uncertainty and that,in turn, is used to tune a local fit to neighbouring residuals to refine the corrections. Resulting3-D traveltime estimates are then used as a description of the forward problem in a nonlineargrid-search relocation. Data are weighted according to the estimated uncertainty. Data fromthe Icelandic Southern Iceland Lowlands (SIL) national seismic network are used to test thestrategy. ETTs are estimated forP- andS-waves at 65 stations in the SIL network, basedon four million arrival time readings of 300.000 events registered between 1990 and 2012.ETTs are strongly correlated for the two wave types. The spatial variations of the predictedcorrections are consistently comparable or somewhat less forS-waves thanP-waves. Thisfeature suggests variations of theVP/VSratio in the Icelandic crust. Error estimates are alsostrongly correlated for the two wave types and between nearby stations. Relocations aretested by comparison with explosions and small populations of well-located events withindenser subnetworks. Relocations result in modestly enhanced clustering of explosions andearthquakes and significantly improved depth estimates. Estimates of the random relocationerror are statistically better behaved than those of the SIL system. They are in general reduced,as is expected since 3-D heterogeneity is now partly taken into account.
|
|
| 4. |
- Andersson, Fredrik, et al.
(author)
-
Flawless diamond separation in simultaneous source acquisition byseismic apparition
- 2017
-
In: Geophysical Journal International. - : Oxford University Press (OUP). - 0956-540X .- 1365-246X. ; 209:3, s. 1793-1739
-
Journal article (peer-reviewed)abstract
- In this paper we prove that the recently introduced method of signal apparition optimally separates signals from interfering sources recorded during simultaneous source seismic data acquisition. By utilizing a periodic sequence of source signatures along one source line, that wavefield becomes separately partially visible in the spectral domain where it can be isolated from interfering signals, processed, and subtracted from the original recordings, thereby separating the wavefields from each other. Whereas other methods for simultaneous source separation can recover data in a triangle-shaped region in the spectral domain, signal apparition allows for the exact separation of data in a diamond-shaped region that is twice as large thereby enabling superior reconstruction of separated wavefields throughout the entire data bandwidth.
|
|
| 5. |
- Bazargan, Mohsen, et al.
(author)
-
Evolution of the statistical distribution of crystal orientations in time- and space-varying viscous flows
- 2019
-
In: Geophysical Journal International. - : Oxford University Press (OUP). - 0956-540X .- 1365-246X. ; 218:2, s. 773-786
-
Journal article (peer-reviewed)abstract
- Magmas and other viscously deforming fluids in the Earth frequently contain embedded crystals or other solid inclusions. These inclusions generally rotate about their own axis and, under certain conditions, align themselves in a direction dictated by the details of the flow. This rotational behaviour has been studied extensively for homogeneous flows. Here, we couple the crystal rotation dynamics with the fluid mechanical Navier-Stokes equations for the large-scale flow, thus allowing the analysis of crystal rotations in settings that are variable in both space and time. The solution is valid provided that the intercrystal spacing is sufficiently large to preclude interaction between crystals. Additionally, we derive an evolution equation for the probability density function (PDF) of crystal orientations based on the fundamental concept of conservation of generic properties in continuum mechanics. The resulting system of equations is extensively tested against previous analytical and numerical solutions. Given the focus on method validation, we limit the fluid mechanics to simple systems with analytical solutions for the velocity field. Even for the simple examples computed, all of which are characterized by fluid flow that is constant in time, the crystal orientation patterns are spatially complex and change in time. Pressure-driven flow in a channel results in coherent bands of crystal orientations with band thickness decreasing towards the channel walls. In corner flow constrained by two mutually perpendicular walls, the pattern of crystal orientations does not exhibit any significant similarity with the flow field. Given that there is no local one-to-one correspondence between the flow and the PDF pattern, a combined and larger-scale solution of the two systems is generally required. The simple flow examples shown demonstrate the viability of this new approach. Application to more complex flow geometries which may commonly occur in nature is deferred to future studies.
|
|
| 6. |
- Biete, Cristina, et al.
(author)
-
The influence of inherited continental margin structures on the stress and strain fields of the south-central Taiwan fold-and-thrust belt
- 2019
-
In: Geophysical Journal International. - : OXFORD UNIV PRESS. - 0956-540X .- 1365-246X. ; 219:1, s. 430-448
-
Journal article (peer-reviewed)abstract
- In this paper we test whether or not structural and morphological features inherited from the Eurasian continental margin are affecting the contemporary stress and strain fields in south-central Taiwan. Principal stress directions (sigma(1), sigma(2) and sigma(3)) are estimated from the inversion of clustered earthquake focal mechanisms and the direction of the maximum compressive horizontal stress (S-H) is calculated throughout the study area. From these data the most likely fault plane orientations and their kinematics are inferred. The results of the stress inversion are then discussed together with the directions of displacement, compressional strain rate and maximum shear strain rate derived from GPS data. These data show that there is a marked contrast in the direction of SH from north to south across the study area, with the direction of SH remaining roughly subparallel to the relative plate motion vector in the north, whereas in the south it rotates nearly 45 degrees counter-clockwise. The direction of the horizontal maximum compression strain rate (epsilon(H)) and associated maximum shear planes, together with the displacement field display an overall similar pattern between them, although undergoing a less marked rotation. We interpret the southward change in the SH, eH and the dextral maximum shear plane directions, together with that of the horizontal displacement field to be related to the reactivation of east-northeast striking faults inherited from the rifted Eurasian margin and to the shelf/slope break. Inherited faults in the basement are typically reactivated as strike-slip faults, whereas newly formed faults in the fold-and-thrust belt are commonly thrusts or oblique thrusts. Eastwards, the stress inversions and strain data show that the western flank of the Central Range is undergoing extension in the upper crust. SH in the Central Range is roughly parallel to the relative plate convergence vector, but in southwestern Taiwan it undergoes a marked counter-clockwise rotation westwards across the Chaochou fault. Farther north, however, there is no significant change across the Lishan fault. This north to south difference is likely due to different margin structures, although local topographic effects may also play a role.
|
|
| 7. |
|
|
| 8. |
- Chen, Chaojian, et al.
(author)
-
Exact solutions of the vertical gravitational anomaly for a polyhedral prism with vertical polynomial density contrast of arbitrary orders
- 2018
-
In: Geophysical Journal International. - : Oxford University Press. - 0956-540X .- 1365-246X. ; 214:3, s. 2115-2132
-
Journal article (peer-reviewed)abstract
- We present general closed-form solutions for the vertical gravitational anomaly caused by a polyhedral prism with mass density contrast varying with depth. Our equations are the first ones to implement a polynomial vertical mass density contrast of arbitrary order. Singularities in the gravity field which arise when the observation site is close to or in the anomalous polyhedral prism are removed in our analytic expressions. Therefore, the observation site can be located outside, on the faces of or inside the anomalous mass bodies. A simple prismatic body of anomalous density is adopted to test the accuracy of our newly developed closed-form solution. Cases of constant, linear, quadratic, cubic and quartic polynomial orders of mass density contrast are tested. For cases of constant, linear, quadratic and cubic polynomial orders, the relative errors between our results and other published exact solutions are less than 10(-11)%. For the case of quartic polynomial order, relative errors less than 10(-10)% are obtained between our solutions and those computed by a high-order Gaussian quadrature rule (512 x 512 x 512 = 134 217 728 quadrature points), where our new analytic solution needs significantly less computational time (0.0009 versus 31.106 s). These numerical experiments not only verified the accuracy of our new formula but also demonstrated their potential in computing exact gravity anomalies for complicated mass density distributions in the Earth.
|
|
| 9. |
- Cherevatova, M., et al.
(author)
-
A multi-resolution approach to electromagnetic modelling
- 2018
-
In: Geophysical Journal International. - : Oxford University Press. - 0956-540X .- 1365-246X. ; 214:1, s. 656-671
-
Journal article (peer-reviewed)abstract
- We present a multi-resolution approach for 3-D magnetotelluric forward modelling. Our approach is motivated by the fact that fine-grid resolution is typically required at shallow levels to adequately represent near surface inhomogeneities, topography and bathymetry, while a much coarser grid may be adequate at depth where the diffusively propagating electromagnetic fields are much smoother. With a conventional structured finite difference grid, the fine discretization required to adequately represent rapid variations near the surface is continued to all depths, resulting in higher computational costs. Increasing the computational efficiency of the forward modelling is especially important for solving regularized inversion problems. We implement a multi-resolution finite difference scheme that allows us to decrease the horizontal grid resolution with depth, as is done with vertical discretization. In our implementation, the multi-resolution grid is represented as a vertical stack of subgrids, with each subgrid being a standard Cartesian tensor product staggered grid. Thus, our approach is similar to the octree discretization previously used for electromagnetic modelling, but simpler in that we allow refinement only with depth. The major difficulty arose in deriving the forward modelling operators on interfaces between adjacent subgrids. We considered three ways of handling the interface layers and suggest a preferable one, which results in similar accuracy as the staggered grid solution, while retaining the symmetry of coefficient matrix. A comparison between multi-resolution and staggered solvers for various models shows that multi-resolution approach improves on computational efficiency without compromising the accuracy of the solution.
|
|
| 10. |
- Eshagh, Mehdi, 1977-
(author)
-
A theoretical discussion on Vening Meinesz-Moritz inverse problem of isostasy
- 2016
-
In: Geophysical Journal International. - : Oxford University Press (OUP). - 0956-540X .- 1365-246X. ; 207:3, s. 1420-1431
-
Journal article (peer-reviewed)abstract
- The Moho surface can be determined according to isostatic theories and the recent Vening Meinesz-Moritz (VMM) theory of isostasy has been successful for this purpose. In this paper, we will study this method from a theoretical prospective and try to find its connection to the Airy-Heiskanen (AH) and Vening Meinesz original theories. We develop Jeffrey’s inverse solution to isostasy according to the recent developments of the VMM method and compare both methods in similar situations. We will show that they are generalisations of the AH model in a global and continuous domain. In the VMM spherical harmonic solution for Moho depth, the mean Moho depth contributes only to the zero-degree term of the series, whilst in Jeffrey’s solution it contributes to all frequencies. We improve the VMM spherical harmonic series further so that the mean Moho can contribute to all frequencies of the solution. This modification makes the VMM global solution superior to the Jeffrey one, but in a global scale, the difference between both solutions is less than 3 km. Both solutions are asymptotically-convergent and we present two methods to obtain smooth solutions for Moho from them.
|
|
