| 21. |
|
|
| 22. |
- Buntin, Laura M., et al.
(författare)
-
Improved accuracy of plane-wave electromagnetic modelling by application of the total and scattered field decomposition and perfectly matched layers
- 2023
-
Ingår i: Geophysical Journal International. - : Oxford University Press (OUP). - 0956-540X .- 1365-246X. ; 235:2, s. 1201-1217
-
Tidskriftsartikel (refereegranskat)abstract
- In 2-D magnetotelluric modelling, the standard application of Dirichlet boundary conditions (BC) may severely diminish the solution accuracy, because the unknown scattered part of the electromagnetic field is erroneously reflected at the domain boundary. Therefore, we adapt the total and scattered field decomposition (TSFD) to geophysical modelling, enabling the application of fully absorbing boundary methods, here perfectly matched layers (PML), to the scattered field. Our novel TSFD divides the modelling domain into two regions. In the total-field region containing the area of interest, the solution is computed for the total field. In the scattered-field region containing the boundaries, the solution is obtained for the scattered field, which is fully absorbed by PML at the boundaries. The plane-wave source is excited at the TSFD interface between both regions. Thus, boundary reflections are eradicated leading to superior solution accuracy, and boundaries can be placed closer to the receivers, shrinking the computational problem. Especially for challenging models with strong lateral changes, the solution accuracy of the TSFD is superior to that of the standard Dirichlet approach. Owing to the linearity of Maxwell's equations, the inaccuracy introduced to the electric and magnetic fields by using Dirichlet BC can be expected to partly cancel out in the magnetotelluric transfer functions, for example the impedance tensor. In this work, we quantify this cancellation effect. The inaccuracy is less than typical measurement errors in the vast majority of apparent resistivity and phase data, even, when the primary fields are strongly inaccurate.
|
|
| 23. |
- Burchardt, Steffi, et al.
(författare)
-
Sinking of anhydrite blocks within a Newtonian salt diapir : modelling the influence of block aspect ratio and salt stratification
- 2012
-
Ingår i: Geophysical Journal International. - 0956-540X .- 1365-246X. ; 188:3, s. 763-778
-
Tidskriftsartikel (refereegranskat)abstract
- 2-D Finite Differences models are used to analyse the strain produced by gravity-driven sinking of dense rectangular inclusions through homogeneous and vertically stratified Newtonian salt. We systematically modelled the descent of dense blocks of different sizes and initial orientations (aspect ratios) representing the Main Anhydrite fragments documented within, for example, the Gorleben salt diapir. Model results demonstrate that size of the blocks is a governing parameter which dictates the amount of strain produced within the block and in the surrounding host salt. Initial block orientation (aspect ratio), on the other hand, causes fundamental differences in block deformation, while the resulting structures produced in the salt are principally the same in all models with homogeneous salt, covering shear zones and folding of passive markers. In models with vertically stratified salt with different viscosities, block descent takes place along complex paths. This results from greater strain accommodation by the salt formation with the lowest viscosity and an asymmetrical distribution of initial vertical shear stresses around the block. Consequently, in these models, block strain is lower compared with the models with homogeneous salt (for the same viscosity as the high-viscosity salt), and sinking is accompanied by block rotation. The latter causes diapir-scale disturbance of the pre-sinking salt stratigraphy and complex sinking paths of the blocks. In particular, vertically oriented blocks sink into high-viscosity salt and drag with them some low-viscosity salt, while horizontal blocks sink in the low-viscosity salt. The resultant sinking velocities vary strongly depending on the sinking path of the block. Based on model results and observed structural configuration within the Gorleben salt diapir, we conclude that the internal complexity of a salt diapir governs its post-ascent deformation. Salt structure and its interaction with dense blocks should hence be considered in the assessment of the long-term stability of storage sites for hazardous waste.
|
|
| 24. |
- Chemia, Z., et al.
(författare)
-
Numerical modelling of rise and fall of a dense layer in salt diapirs
- 2008
-
Ingår i: Geophysical Journal International. - 0956-540X .- 1365-246X. ; 172:2, s. 798-816
-
Tidskriftsartikel (refereegranskat)abstract
- Numerical models are used to study the entrainment of a dense anhydrite layer by a diapir. The anhydrite layer is initially horizontally embedded within a viscous salt layer. The diapir is down-built by aggradation of non-Newtonian sediments (n = 4, constant temperature) placed on the top of the salt layer. Several parameters (sedimentation rate, salt viscosity, perturbation width and stratigraphic position of the anhydrite layer) are studied systematically to understand their role in governing the entrainment of the anhydrite layer. High sedimentation rates during the early stages of the diapir evolution bury the initial perturbation and, thus, no diapir forms. The anhydrite layer sinks within the buried salt layer. For the same sedimentation rate, increasing viscosity of the salt layer decreases the rise rate of the diapir and reduces the amount (volume) of the anhydrite layer transported into the diapir. Model results show that viscous salt is capable of carrying separate blocks of the anhydrite layer to relatively higher stratigraphic levels. Varying the width of the initial perturbation (in our calculations 400-800 m), from which a diapir triggers, shows that wider diapirs can more easily entrain an embedded anhydrite layer than the narrower diapirs. The anhydrite layer is entrained as long as rise rate of the diapir exceeds the descent rate of the denser anhydrite layer. We conclude that the four parameters mentioned above govern the ability of a salt diapir to entrain an embedded dense layer. However, the model results show that the entrained blocks inevitably sink back if the rise rate of the diapir is less than the rate of descent of the anhydrite layer or the diapir is permanently covered by a stiff overburden in case of high sedimentation rates.
|
|
| 25. |
- Chen, Chaojian, et al.
(författare)
-
Exact solutions of the vertical gravitational anomaly for a polyhedral prism with vertical polynomial density contrast of arbitrary orders
- 2018
-
Ingår i: Geophysical Journal International. - : Oxford University Press. - 0956-540X .- 1365-246X. ; 214:3, s. 2115-2132
-
Tidskriftsartikel (refereegranskat)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.
|
|
| 26. |
- Cherevatova, M., et al.
(författare)
-
A multi-resolution approach to electromagnetic modelling
- 2018
-
Ingår i: Geophysical Journal International. - : Oxford University Press. - 0956-540X .- 1365-246X. ; 214:1, s. 656-671
-
Tidskriftsartikel (refereegranskat)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.
|
|
| 27. |
- Eggertsson, Gunnar, et al.
(författare)
-
Earthquake or blast? : Classification of local-distance seismic events in Sweden using fully connected neural networks
- 2024
-
Ingår i: Geophysical Journal International. - : Oxford University Press. - 0956-540X .- 1365-246X. ; 236:3, s. 1728-1742
-
Tidskriftsartikel (refereegranskat)abstract
- Distinguishing between different types of seismic events is a task typically performed manually by expert analysts and can thus be both time and resource expensive. Analysts at the Swedish National Seismic Network (SNSN) use four different event types in the routine analysis: natural (tectonic) earthquakes, blasts (e.g. from mines, quarries and construction) and two different types of mining-induced events associated with large, underground mines. In order to aid manual event classification and to classify automatic event definitions, we have used fully connected neural networks to implement classification models which distinguish between the four event types. For each event, we bandpass filter the waveform data in 20 narrow-frequency bands before dividing each component into four non-overlapping time windows, corresponding to the P phase, P coda, S phase and S coda. In each window, we compute the root-mean-square amplitude and the resulting array of amplitudes is then used as the neural network inputs. We compare results achieved using a station-specific approach, where individual models are trained for each seismic station, to a regional approach where a single model is trained for the whole study area. An extension of the models, which distinguishes spurious phase associations from real seismic events in automatic event definitions, has also been implemented. When applying our models to evaluation data distinguishing between earthquakes and blasts, we achieve an accuracy of about 98 per cent for automatic events and 99 per cent for manually analysed events. In areas located close to large underground mines, where all four event types are observed, the corresponding accuracy is about 90 and 96 per cent, respectively. The accuracy when distinguishing spurious events from real seismic events is about 95 per cent. We find that the majority of erroneous classifications can be traced back to uncertainties in automatic phase picks and location estimates. The models are already in use at the SNSN, both for preliminary type predictions of automatic events and for reviewing manually analysed events.
|
|
| 28. |
- Eken, Tuna, 1978-, et al.
(författare)
-
Effects of seismic anisotropy on P-velocity tomography of the Baltic Shield
- 2012
-
Ingår i: Geophysical Journal International. - 0956-540X .- 1365-246X. ; 188:2, s. 600-612
-
Tidskriftsartikel (refereegranskat)abstract
- We investigate possible effects of neglecting seismic anisotropy on standard isotropic P-velocity tomographic images of the upper mantle beneath the Baltic shield. Isotropic inversions of teleseismic P- and S-wave traveltimes exhibit alternating high- and low-velocity heterogeneities down to depths of over 400 km. Differences in tomographic inversions of SV- and SH-wave traveltimes are distinct down to depths of about 200 km and are associated with anisotropy of the lithospheric mantle. Anisotropic structures of the upper mantle affect both the P and S traveltimes, shear-wave splitting as well as the P polarization directions. Joint inversion for isotropic and anisotropic velocity perturbations is not feasible due to the limited 3-D ray coverage of available data. Therefore, we correct the input traveltimes for anisotropic contributions derived from independent analyses and then perform standard isotropic inversions. These corrections are derived either directly from directional deviations of P-wave propagation or are calculated in anisotropic models retrieved by joint inversions of body-wave anisotropic parameters (P-residual spheres and shear-wave splitting). These anisotropic models are also used to fit backazimuth variations of P-wave polarization directions. General features of tomographic images calculated from the original and the anisotropy-corrected data are similar. Amplitudes of the velocity perturbations decrease below similar to 200 km depth, that is in the sub-lithospheric mantle. In general, large-scale anisotropy related to the fabrics of the continental mantle lithosphere can contaminate tomographic images in some parts of models and should not be ignored.
|
|
| 29. |
- Eken, Tuna, et al.
(författare)
-
Upper-mantle structure of the Baltic Shield below the Swedish National Seismological Network (SNSN) resolved by teleseismic tomography
- 2007
-
Ingår i: Geophysical Journal International. - 0956-540X .- 1365-246X. ; 169:2, s. 617-630
-
Tidskriftsartikel (refereegranskat)abstract
- Upper-mantle structure under the Baltic Shield is studied using non-linear high resolution teleseismic P-phase tomography. Observed relative arrival-time residuals from 52 teleseismic earthquakes recorded by the Swedish National Seismological Network (SNSN) are inverted to delineate the structure of the upper mantle. The network consists of 47 (currently working) three-component broad-band stations located in an area about 450 km wide and 1450 km long. In order to reduce complications due to possible significant three-dimensionality of Earth structure, events chosen for this study lay close to in-line with the long-axis of the array (±30°) . Results indicate P-wave velocity perturbations of ±3 per cent down to at least 470 km below the network. The size of the array allows inversion for structures even at greater depths, and lateral variations of velocity at depths of up to 680 km appear to be resolved. Below the central part of the array (60°–64° N), where ray coverage is best, the data reveals a large region of relatively low velocity at depths of over about 300 km. At depths less than about 250–300 km, the models include a number of features, including an apparent slab-like structure dipping gently towards the north.
|
|
| 30. |
- Elming, SA, et al.
(författare)
-
Results of magnetotelluric and gravimetric measurements in western Nicaragua, central America
- 1997
-
Ingår i: GEOPHYSICAL JOURNAL INTERNATIONAL. - : BLACKWELL SCIENCE LTD. - 0956-540X. ; 128:3, s. 647-658
-
Tidskriftsartikel (övrigt vetenskapligt/konstnärligt)abstract
- Magnetotelluric and gravity data have been collected within a ca. 170 km long traverse running from the Pacific coast of Nicaragua in the west to the Nicaraguan Highland in the east. This part of Nicaragua is characterized by sedimentary rocks of the Paci
|
|
