| 1. |
- Eshagh, Mehdi, Professor, 1977-, et al.
(författare)
-
Comparison of gravimetric and mantle flow solutions for sub-lithopsheric stress modeling and their combination
- 2018
-
Ingår i: Geophysical Journal International. - : Oxford University Press (OUP). - 0956-540X .- 1365-246X. ; 213:2, s. 1013-1028
-
Tidskriftsartikel (refereegranskat)abstract
- Based on Hager and O’Connell’s solution to mantle flow equations, the stresses induced by mantle convection are determined using the density and viscosity structure in addition to topographic data and a plate velocity model. The solution to mantle flow equations requires the knowledge of mantle properties that are typically retrieved from seismic information. Large parts of the world are, however, not yet covered sufficiently by seismic surveys. An alternative method of modeling the stress field was introduced by Runcorn. He formulated a direct relation between the stress field and gravity data, while adopting several assumptions, particularly disregarding the toroidal mantle flow component and mantle viscosity variations. A possible way to overcome theoretical deficiencies of Runcorn’s theory as well as some practical limitations of applying Hager and O’Connell’s theory (in the absence of seismic data) is to combine these two methods. In this study, we apply a least-squares analysis to combine these two methods based on the gravity data inversion constraint on mantle flow equations. In particular, we use vertical gravity gradients from the Gravity field and steady state Ocean Circulation Explorer that are corrected for the gravitational contribution of crustal density heterogeneities prior to applying a localized gravity-gradient inversion. This gravitational contribution is estimated based on combining the Vening Meinesz-Moritz and flexural isostatic theories. Moreover, we treat the non-isostatic effect implicitly by applying a band-limited kernel of the integral equation during the inversion. In numerical studies of modeling, the stress field within the South American continental lithosphere we compare the results obtained after applying Runcorn and Hager and O’Connell’s methods as well as their combination. The results show that, according to Hager and O’Connell’s (mantle flow) solution, the maximum stress intensity is inferred under the northern Andes. Additional large stress anomalies are detected along the central and southern Andes, while stresses under most of old, stable cratonic formations aremuch less pronounced or absent. A prevailing stress-vector orientation realistically resembles a convergent mantle flow and downward currents under continental basins that separate Andean Orogeny from the Amazonian Shield and adjacent cratons. Runcorn’s (gravimetric) solution, on the other hand, reflects a tectonic response of the lithosphere to mantle flow, with the maximum stress intensity detected along the subduction zone between the Nazca and Altiplano plates and along the convergent tectonic margin between the Altiplano and South American plates. The results also reveal a very close agreement between the results obtained from the combined and Hager and O’Connell’s solutions. © The Author(s) 2018. Published by Oxford University Press on behalf of The Royal Astronomical Society.
|
|
| 2. |
- Eshagh, Mehdi, Professor, 1977-, et al.
(författare)
-
Lithospheric stress, strain and displacement changes from GRACE-FO time-variable gravity : case study for Sar-e-Pol Zahab Earthquake 2018
- 2020
-
Ingår i: Geophysical Journal International. - : Oxford University Press (OUP). - 0956-540X .- 1365-246X. ; 223:1, s. 379-397
-
Tidskriftsartikel (refereegranskat)abstract
- Temporal variations in the Earth’s gravity field can be used for monitoring of lithospheric deformations. The network of continuously operating gravity stations is required for this purpose but a global coverage by such network is currently extremely sparse. Temporal variations in long-wavelength part of the Earth’s gravity field have been, however, observed by two satellite missions, namely the Gravity Recovery And Climate Experiment (GRACE) and the GRACE Follow-On (GRACE-FO). These satellite gravity observations can be used to study long-wavelength deformations of the lithosphere. Consequently, considering the lithosphere as a spherical elastic shell and solving the partial differential equation of elasticity for it, the stress, strain and displacement inside the lithosphere can be estimated. The lower boundary of this shell is assumed to be stressed by mantle convection, which has a direct relation to the Earth’s gravity field according to Runcorn’s theory. Changes in gravity field lead to changes in the sublithospheric stress and the stress propagated throughout the lithosphere. In this study, we develop mathematical models in spherical coordinates for describing the stress propagation from the sublithosphere through the lithosphere. We then organize a system of observation equations for finding a special solution to the boundary-value problem of elasticity in the way that provides a stable solution. In contrast, models presented in previously published studies are ill-posed. Furthermore, we use constants of the solution determined from the boundary stresses to determine the strain and displacements leading to these stresses, while in previous studies only the stress has been considered according to rheological properties of the lithosphere. We demonstrate a practical applicability of this theoretical model to estimate the stress–strain redistribution caused by the Sar-e-Pol Zahab 2018 earthquake in Iran by using the GRACE-FO monthly solutions.
|
|
| 3. |
- Pitonak, Martin, et al.
(författare)
-
Spectral combination of spherical gravitational curvature boundary-value problems
- 2018
-
Ingår i: Geophysical Journal International. - : Oxford University Press (OUP). - 0956-540X .- 1365-246X. ; 214:2, s. 773-791
-
Tidskriftsartikel (refereegranskat)abstract
- Four solutions of the spherical gravitational curvature boundary-value problems can be exploited for the determination of the Earth’s gravitational potential. In this paper we discuss the combination of simulated satellite gravitational curvatures, that is, components of the third-order gravitational tensor, by merging these solutions using the spectral combination method. For this purpose, integral estimators of biased-and unbiased-types are derived. In numerical studies, we investigate the performance of the developed mathematical models for the gravitational field modelling in the area of Central Europe based on simulated satellite measurements. First, we verify the correctness of the integral estimators for the spectral downward continuation by a closed-loop test. Estimated errors of the combined solution are about eight orders smaller than those from the individual solutions. Second, we perform a numerical experiment by considering the Gaussian noise with the standard deviation of 6.5 x 10(-17) m(-1) s(-2) in the input data at the satellite altitude of 250 km above the mean Earth sphere. This value of standard deviation is equivalent to a signal-to-noise ratio of 10. Superior results with respect to the global geopotential model TIM-r5 (Brockmann et al. 2014) are obtained by the spectral downward continuation of the vertical-vertical-vertical component with the standard deviation of 2.104 m(2) s(-2), but the root mean square error is the largest and reaches 9.734 m(2) s(-2). Using the spectral combination of all gravitational curvatures the root mean square error is more than 400 times smaller but the standard deviation reaches 17.234 m(2) s(-2). The combination of more components decreases the root mean square error of the corresponding solutions while the standard deviations of the combined solutions do not improve as compared to the solution from the vertical-vertical-vertical component. The presented method represents a weight mean in the spectral domain that minimizes the root mean square error of the combined solutions and improves standard deviation of the solution based only on the least accurate components.
|
|
