| 1. |
- Abdollahi, Somayeh, et al.
(författare)
-
3D joint inversion of gravity data and Rayleigh wave group velocities to resolve shear-wave velocity and density structure in the Makran subduction zone, south-east Iran
- 2019
-
Ingår i: Journal of Asian Earth Sciences. - : PERGAMON-ELSEVIER SCIENCE LTD. - 1367-9120 .- 1878-5786. ; 173, s. 275-290
-
Tidskriftsartikel (refereegranskat)abstract
- In this study, we developed a method to invert jointly Rayleigh wave group velocities and gravity anomalies for velocity and density structure of the lithosphere. We applied the method to the Makran accretionary prism, SE Iran. The reason for using different data sets is that each of these data sets is sensitive to different parameters. Surface wave group velocities are sensitive mainly to shear wave velocity distribution in depth but do not well resolve density variations. Therefore, joint inversion with gravity data increases the resolution of density distribution. Our approach differs from others mainly in the model parameterization: Instead of subdividing the model into a large number of thin layers, we invert for the properties of only four layers: thickness, P- and S-wave velocities and densities and their vertical gradients in sediments, upper-crust, lower-crust and upper mantle. The method is applied first to synthetic models in order to demonstrate its usefulness. We then applied the method to real data to investigate the lithosphere structure beneath the Makran. The resulting model shows that Moho depth increases from Oman Sea (18-33 km) and Makran fore-arc (33-37 km) to the volcanic-arc (44-46 km). The crustal density is high in the Oman Sea as should be expected for the oceanic crust. We also find a high-velocity anomaly in the upper mantle under the Oman Sea corresponding to the subducting slab. The crust under the fore-arc, volcanic-arc and back-arc settings of Makran subduction zone is characterized by low-velocity zones.
|
|
| 2. |
- Abdollahi, Somayeh, et al.
(författare)
-
Crustal and upper mantle structures of Makran subduction zone, SE Iran by combined surface wave velocity analysis and gravity modeling
- 2018
-
Ingår i: Tectonophysics. - : ELSEVIER SCIENCE BV. - 0040-1951 .- 1879-3266. ; 747, s. 191-210
-
Tidskriftsartikel (refereegranskat)abstract
- The inversion of Rayleigh wave group velocity dispersion curves is challenging, because it is non-linear and multimodal. In this study, we develop and test a new Rayleigh wave dispersion curve inversion scheme using the Shuffled Complex Evolution (SCE) algorithm. Incorporating this optimization algorithm into the inverse procedure not only can effectively locate the promising areas in the solution space for a global minimum but also avoids its wandering near the global minimum in the final stage of search. In addition, our approach differs from others in the model parameterization: Instead of subdividing the model into a large number of thin layers, we invert for thickness, velocities and densities and their vertical gradients of four layers, sediments, upper-crust, lower-crust and upper mantle. The proposed inverse procedure is applied to non-linear inversion of fundamental mode Rayleigh wave group dispersion curves for shear and compressional wave velocities. At first, to determine the efficiency and stability of the SCE method, two noise-free and two noisy synthetic data sets are inverted. Then real data for Makran region in SE Iran are inverted to examine the usage and robustness of the proposed approach on real surface wave data. In a second step, we applied 3D Gravity Modeling based on surface wave analysis results to obtain the density structure and thickness of each layer. The reason for using both types of data sets, is that gravity anomaly has a bad vertical resolution and surface wave group velocities are good for placing layer limits at depth, but they are not very sensitive to densities. Therefore, using gravity data increases the overall resolution of density distribution. In a final step, we used again the SCE method to invert the fundamental mode Rayleigh wave group dispersion curves based on the gravity results. Gravity results like thicknesses and sediment densities have been used to constrain the limit of search space in the SCE method. Results show a high shear and compressional velocity under the Gulf of Oman which reduce to the North of the Makran region. The Moho depth of the Oman Gulf is about 18-28 km and it increases to 46-48 km under the Taftan-Bazman volcanic-arc. The density image shows an average crustal density with maximum values under the Gulf of Oman decreasing northward to the Makran.
|
|
| 3. |
- Amini, Samar, et al.
(författare)
-
Fault slip and identification of the second fault plane in the Varzeghan earthquake doublet
- 2018
-
Ingår i: Journal of Seismology. - : SPRINGER. - 1383-4649 .- 1573-157X. ; 22:4, s. 815-831
-
Tidskriftsartikel (refereegranskat)abstract
- An intraplate earthquake doublet, with 11-min delay between the events, devastated the city of Varzeghan in northwestern Iran on August 11, 2012. The first Mw 6.5 strike-slip earthquake, which occurred after more than 200 years of low seismicity, was followed by an Mw 6.4 oblique thrust event at an epicentral separation of about 6 km. While the first event can be associated with a distinct surface rupture, the absence of a surface fault trace and no clear aftershock signature makes it challenging to identify the fault plane of the second event. We use teleseismic body wave inversion to deduce the slip distribution in the first event. Using both P and SH waves stabilize the inversion and we further constrain the result with the surface rupture extent and the aftershock distribution. The obtained slip pattern shows two distinct slip patches with dissimilar slip directions where aftershocks avoid high-slip areas. Using the estimated slip for the first event, we calculate the induced Coulomb stress change on the nodal planes of the second event and find a preference for higher Coulomb stress on the N-S nodal plane. Assuming a simple slip model for the second event, we estimate the combined Coulomb stress changes from the two events on the focal planes of the largest aftershocks. We find that 90% of the aftershocks show increased Coulomb stress on one of their nodal planes when the N-S plane of the second event is assumed to be the correct fault plane.
|
|
| 4. |
- Amini, Samar
(författare)
-
Source analysis of multiplet earthquakes (two case studies in Iran)
- 2020
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Multiplet earthquakes are large earthquakes of similar magnitude which occur close in time in the same limited geographical area. They are not common but they considerably increase the potential hazard in the area in which they occur. This thesis studies source properties and triggering mechanisms of two sets of multiplet events in Iran, which both occurred in unexpected areas, but close to some major active fault systems. The first multiplet is an earthquake doublet (Mw 6.5 and Mw 6.4) which occurred in northwestern Iran and caused more than 300 fatalities and significant injuries. In paper I, a teleseismic body-waveform inversion was used to deduce the slip distribution pattern on the fault plane of the first mainshock. The estimated slip pattern was utilized to calculate the Coulomb stress changes on the second fault plane and on the following aftershocks. Based on this analysis, the ambiguity between the primary and auxiliary fault plane of the second mainshock could be resolved. The second set of events is a triplet (Mw 6.1 - 6.0) that occurred in eastern Iran, close to the Kerman province. In paper II, the rupture propagation patterns of the three mainshocks were analyzed using Empirical Green’s Function (EGF) deconvolution. Two different approaches were used: One, the analysis of the azimuthal variation of the apparent rupture duration based on the width of the observed relative source time functions, and the second, the analysis of along-strike rupture directivity by assessing azimuthal variations of the relative amplitude spectra. The second approach was also used to investigate the rupture directivity of the largest aftershocks in the sequence (Mw 5 - 5.5). A clear tendency for rupture propagation towards the northwest was observed for the sequence, which suggests that the regional stress field has a central role in controlling the rupture propagation direction. In paper III, the slip distribution patterns of the triplet earthquakes were analyzed using teleseismic body-waveform inversion, and the similarities and differences in the rupture processes of the three mainshocks were investigated. Using the Coulomb stress analyses, the stress interactions between the mainshocks were examined, leading to identification of the primary and auxiliary planes. Finally, we suggest that the hazard estimates in complex continental regions such as Iran need to consider the probability of multiplets, which might allow a reduction of the seismic risk associated to the occurrence of further large earthquakes subsequent to a devastating earthquake.
|
|
| 5. |
- Basir, Hadi Mahdavi, et al.
(författare)
-
Acoustic wave propagation simulation by reduced order modelling
- 2018
-
Ingår i: Exploration Geophysics. - : CSIRO PUBLISHING. - 0812-3985 .- 1834-7533. ; 49:3, s. 386-397
-
Tidskriftsartikel (refereegranskat)abstract
- Wave propagation simulation, as an essential part of many algorithms in seismic exploration, is associated with high computational cost. Reduced order modelling(ROM) is a known technique in many different applications that can reduce the computational cost of simulation by employing an approximation of the model parameters. ROM can be carried out using different algorithms. The method proposed in this work is based on using the most important mode shapes of the model, which can be computed by an efficient numerical method. The numerical accuracy and computational performance of the proposed method were investigated over a simple synthetic velocity model and a portion of the SEG/EAGE velocity model. Different boundary conditions were discussed, and among them the random boundary condition had higher performance for applications like reverse time migration (RTM). The capability of the proposed method for RTM was evaluated and confirmed by the synthetic velocity model of SEG/EAGE. The results showed that the proposed ROM method, compared with the conventional finite element method (FEM), can decrease the computational cost of wave propagation modelling for applications with many simulations like the reverse time migration. Depending on the number of simulations, the proposed method can increase the computational efficiency by several orders of magnitudes.
|
|
| 6. |
- Basir, Hadi Mahdavi, et al.
(författare)
-
Modified imaging condition for reverse time migration based on reduction of modelling time
- 2018
-
Ingår i: Exploration Geophysics. - : CSIRO PUBLISHING. - 0812-3985 .- 1834-7533. ; 49:4, s. 494-505
-
Tidskriftsartikel (refereegranskat)abstract
- Reverse time migration (RTM) is considered as a high-end imaging algorithm due to its ability to image geologically complex environments. However, this algorithm suffers from very high computational costs and low-frequency artefacts. The former drawback is the result of the intensive computations and huge memory allocation involved in RTM. Wave propagation modelling, as a kernel of RTM, demands intensive computations, and conventional imaging conditions are associated with huge memory allocation. In this paper, a modification of imaging condition is proposed that improves the efficiency of RTM as a reduction of computational cost, memory (RAM) allocation and low-frequency artefacts. The proposed imaging condition is similar to the conventional imaging condition but with the reduction of modelling time to near half the maximum time of recording. As the main idea of the proposed imaging condition, the impact of wave propagation modelling time is investigated on the quality of RTM and illumination of reflectors. The performance of the proposed method is considered using two synthetic models (SEG/EAGE and BP) and a real dataset from an Iranian oilfield in the south of Iran. Results showed that the new imaging condition can properly image the reflectors and enhance the efficiency of RTM. By using the proposed imaging condition, we achieved similar to 25% increase in CPU performance and 50% decrease in the memory allocation. Despite the improvement of the performance, results showed that the proposed imaging condition had no significant effect on the illumination.
|
|
| 7. |
- Basir, Hadi Mandavi, et al.
(författare)
-
Reverse time migration by Krylov subspace reduced order modeling
- 2018
-
Ingår i: Journal of Applied Geophysics. - : Elsevier BV. - 0926-9851 .- 1879-1859. ; 151, s. 298-308
-
Tidskriftsartikel (refereegranskat)abstract
- Imaging is a key step in seismic data processing. To date, a myriad of advanced pre-stack depth migration approaches have been developed; however, reverse time migration (RTM) is still considered as the high-end imaging algorithm. The main limitations associated with the performance cost of reverse time migration are the intensive computation of the forward and backward simulations, time consumption, and memory allocation related to imaging condition. Based on the reduced order modeling, we proposed an algorithm, which can be adapted to all the aforementioned factors. Our proposed method benefit from Krylov subspaces method to compute certain mode shapes of the velocity model computed by as an orthogonal base of reduced order modeling. Reverse time migration by reduced order modeling is helpful concerning the highly parallel computation and strongly reduces the memory requirement of reverse time migration. The synthetic model results showed that suggested method can decrease the computational costs of reverse time migration by several orders of magnitudes, compared with reverse time migration by finite element method.
|
|
| 8. |
- Eken, Tuna, 1978-, et al.
(författare)
-
S and P velocity heterogeneities within the upper mantle below the Baltic Shield
- 2008
-
Ingår i: Tectonophysics. - : Elsevier BV. - 0040-1951 .- 1879-3266. ; 462:1-4, s. 109-124
-
Tidskriftsartikel (refereegranskat)abstract
- Upper mantle structure beneath the Baltic (Fennoscandian) Shield is investigated using non-linear tomographic inversion of relative arrival-time residuals. 52 selected teleseismic earthquakes recorded by 45 broadband stations of the Swedish National Seismological Network (SNSN) provide 1532 good quality S-wave relative arrival times. SV and SH arrival-time residuals were initially analyzed independently, providing two separate models. These reveal several consistent major features, many of which are also consistent with P-wave results. Lateral velocity variations of ± 3–4% are observed to depths of at least 470 km. The correlation between the SH and SV models is investigated and shows a pattern of minor but significant differences down to around 150–200 km depth, below which the models are essentially similar. Direct cell by cell comparison of the model velocities reveals a similar pattern, with velocity differences between the models of up to 4%. Numerical tests show that differences in the two S-wave models can only be partially attributed to noise and limited resolution, and some features are attributed to the effect of large scale anisotropy. One of the significant and sharp discrepancies between the S models coincides with a presumed boundary between Archean and Proterozic domains, suggesting different anisotropic characteristics in the two regions.
|
|
| 9. |
- Enferadi, S., et al.
(författare)
-
Feasibility study of earthquake early warning in Tehran, Iran
- 2021
-
Ingår i: Journal of Seismology. - : Springer Nature. - 1383-4649 .- 1573-157X. ; 25:4, s. 1127-1140
-
Tidskriftsartikel (refereegranskat)abstract
- In this study, we examine the scientific feasibility of an Earthquake Early Warning System in Tehran, Iran, by the integration of the Tehran Disaster Mitigation and Management Organization (TDMMO) accelerometric network and the PRobabilistic and Evolutionary early warning SysTem (PRESTo). To evaluate the performance of the TDMMO-PRESTo system in providing the reliable estimations of earthquake parameters and the available lead-times for The Metropolis of Tehran, two different approaches were analyzed in this work. The first approach was assessed by applying the PRESTo algorithms on waveforms from 11 moderate instrumental earthquakes that occurred in the vicinity of Tehran during the period 2009-2020. Moreover, we conducted a simulation analysis using synthetic waveforms of 10 large historical earthquakes that occurred in the vicinity of Tehran. We demonstrated that the six worst-case earthquake scenarios can be considered for The Metropolis of Tehran, which are mostly related to the historical and instrumental events that occurred in the southern, eastern, and western parts of Tehran. Our results indicate that the TDMMO-PRESTo system could provide reliable and sufficient lead-times of about 1 to 15s and maximum lead-times of about 20s for civil protection purposes in The Metropolis of Tehran.
|
|
| 10. |
- Gregersen, S., et al.
(författare)
-
Uniqueness of modeling results from teleseismic P-Wave tomography in Project Tor
- 2010
-
Ingår i: Tectonophysics. - 0040-1951 .- 1879-3266. ; 481:1-4, s. 99-107
-
Tidskriftsartikel (refereegranskat)abstract
- Within Project Tor. which is about Teleseismic Tomography across the Tomquist Zone in Germany-Denmark-Sweden, we have confirmed very significant deep lithosphere differences And modeling is substantiated via completely independent methods. In 1996-1997 our 130 seismographs constituted the largest seismic antenna ever in Europe. The Tor area was chosen along a well studied crustal profile of an earlier project, and the modeling efforts were concentrated on the deep lithosphere and asthenosphere differences to depths around 300 km The Tor data have been subjected to P-wave travel time tomography. surface wave and receiver function analysis as well as anisotropy and scattering measurements An important goal of the project was to make several independent inversions of the tomography data. and compare the results in an attempt to evaluate uniqueness, resolution and accuracy of these inversions. The comparisons of this paper involve more diversity in methods than any previous comparison. The geological outcome is a substantiation of earlier statements that, "The transition is interpreted to be sharp and steep in two places It goes all through the lithosphere at the northern rim of the Tornquist Zone near the border between Sweden and Denmark, and here the lithosphere difference is large to depths more than 200 km. The other lithosphere difference. of smaller scale, is found near the southern edge of the Ringkobing-Fyn High near the border between Denmark and Germany Also this transition is sharp and steep. and goes all through the lithosphere to depths around 120 km. These two sharp transitions divide the Tor region into 3 different lithosphere structures distinguishable in P-wave travel time tomography. surface wave dispersion. P- and S-wave anisotropy and partly in P-wave scattering" The mentioned broad-scale features are judged to be unambiguously determined, with well-described resolution and accuracy Unfortunately a detail like the slope of the subcrustal lithosphere transition right under the Tronquist Zone cannot be constrained even if this is where the resolution is best. and the curiosity largest.
|
|
