Estimating the mechanical anisotropy of the Iranian lithosphere using the wavelet coherence method

We calculated anisotropic wavelet coherence between Bouguer anomaly and topography in order to map the anisotropy of the effective elastic thickness of the Iranian lithosphere (Te). An orthotropic elastic plate model is used for inverting the anisotropic wavelet coherence to compute the mechanical anisotropy through the weak axis of the Te. Anisotropy of the Te-weak axis and the strength of the anisotropic parameter, namely the anisotropic coherence effect over the study area are estimated by restricting the rotated Morlet wavelet (fan wavelet) geometry over an azimuthal range of 90°. Large-scale Te variations have been shown to be associated with phenomena, such as mountain belts, subduction zones, craton boundaries, fault zones, and seismogenic regions. Although the correlation between the major tectonic features of the Iranian lithosphere and the distribution of the Te-weak axis is not general or precise, in some regions, such as the Central Iran Blocks, and the Alborz, Kopeh Dagh, Zagros, and Makran orogenic belts, the weak axis has a uniform or slowly varying pattern which changes over their boundaries. A perpendicular alignment between seismic anisotropy measurements in Iran and the Te-weak directions suggests a lithospheric origin for anisotropy. The correlation between averaged stress directions and the weak axis of the Te in Iran indicates that the present day stress field and the fossil strain are still related. Correlation between these factors suggests vertically coherent deformation of the lithosphere in Iran resulting from the multiply convergent orogenic processes. The complex mechanical anisotropy pattern of the Iranian lithosphere results from the interaction of many pre-existent structures which dominantly control the mechanical anisotropy of the lithosphere

Spatial variations in the effective elastic thickness of the lithosphere in Southeast Asia

As a proxy for long-term lithospheric strength, detailed information on lateral effective elastic thickness (Te) variations can aid in understanding the distribution pattern of surface deformation and its response to long-term forces. Here we present high-resolution maps of spatial variations of Te for the complex SE Asian region by analyzing the coherence of topography and Bouguer gravity anomaly data. We find that after considering the gravity deficit of less dense sediment, the recovered Te maps are more representative of the geology, particularly in elongated rift basins. The results show that the Te variation pattern in SE Asia, in general, agrees well with its tectonic provinces and major tectonic boundaries. The oceanic basins, the Indosinian suture zones between the Indochina and Sibumasu blocks, and the Makassar Strait are characterized by low Te, while moderate and high Te values are recovered in the Khorat plateau, West Burma, the Singapore Ridge, the Con Song Swell, Borneo, the northern Australian margin and the Molucca Sea. The Te pattern in the south Indonesian margin is complicated by the approach and collision of oceanic plateaus and seamounts with the fore-arc region. The heterogeneous strength features are consistent with the complex assemblage of different tectonic units, and significant deformation during Cenozoic tectonic events. In the Indochina Peninsula, the extruded displacement during the India-Eurasia collision might have been partitioned and absorbed by the combined mechanism of the extrusion and viscous tectonic models. As a result, the offshore displacements of the major strike-slip faults in the South China Sea are much smaller than originally assumed, thus having less effect on the development of the South China Sea than other mechanisms such as the slab pull of the proto-South China Sea. Since the displacement driven by the boundary tectonic forces has been greatly absorbed and decreased by subduction and deformation in the active margins and adjacent weak regions, the motion velocity of the interior regions is greatly lower than the boundary active margins, and they are largely free of seismicity and volcanism. Our results suggest that East Borneo might share a similar crustal basement, and represent a broad tectonic zone of the destroyed Meso-Tethys Ocean extending from West-Middle Java, through East Borneo to northern Borneo of the Sarawak and Sabah. The Indosinian zones between the Indochina and Sibumasu blocks might extend further southeastward across Billiton Island to offshore of southern Borneo, and the Singapore platform and SW Borneo might belong to the same block. The results also show that the internal load fraction F is high in the coastal area of South China, the northern margin of the South China Sea, and the coastal area of Indochina, which, in general, agrees with the distribution of a high-velocity lower crustal layer and Late Cenozoic basaltic rocks

Slip distribution and tectonic implication of the 1999 Chi‐Chi, Taiwan, Earthquake

We report on the fault complexity of the large (M_w = 7.6) Chi‐Chi earthquake obtained by inverting densely and well‐distributed static measurements consisting of 119 GPS and 23 doubly integrated strong motion records. We show that the slip of the Chi-Chi earthquake was concentrated on the surface of a ”wedge shaped” block. The inferred geometric complexity explains the difference between the strike of the fault plane determined by long period seismic data and surface break observations. When combined with other geophysical and geological observations, the result provides a unique snapshot of tectonic deformation taking place in the form of very large (>10m) displacements of a massive wedge‐shaped crustal block which may relate to the changeover from over‐thrusting to subducting motion between the Philippine Sea and the Eurasian plates

Active megadetachment beneath the western United States

Geodetic data, interpreted in light of seismic imaging, seismicity, xenolith studies, and the late Quaternary geologic history of the northern Great Basin, suggest that a subcontinental-scale extensional detachment is localized near the Moho. To first order, seismic yielding in the upper crust at any given latitude in this region occurs via an M7 earthquake every 100 years. Here we develop the hypothesis that since 1996, the region has undergone a cycle of strain accumulation and release similar to “slow slip events” observed on subduction megathrusts, but yielding occurred on a subhorizontal surface 5–10 times larger in the slip direction, and at temperatures >800°C. Net slip was variable, ranging from 5 to 10 mm over most of the region. Strain energy with moment magnitude equivalent to an M7 earthquake was released along this “megadetachment,” primarily between 2000.0 and 2005.5. Slip initiated in late 1998 to mid-1999 in northeastern Nevada and is best expressed in late 2003 during a magma injection event at Moho depth beneath the Sierra Nevada, accompanied by more rapid eastward relative displacement across the entire region. The event ended in the east at 2004.0 and in the remainder of the network at about 2005.5. Strain energy thus appears to have been transmitted from the Cordilleran interior toward the plate boundary, from high gravitational potential to low, via yielding on the megadetachment. The size and kinematic function of the proposed structure, in light of various proxies for lithospheric thickness, imply that the subcrustal lithosphere beneath Nevada is a strong, thin plate, even though it resides in a high heat flow tectonic regime. A strong lowermost crust and upper mantle is consistent with patterns of postseismic relaxation in the southern Great Basin, deformation microstructures and low water content in dunite xenoliths in young lavas in central Nevada, and high-temperature microstructures in analog surface exposures of deformed lower crust. Large-scale decoupling between crust and upper mantle is consistent with the broad distribution of strain in the upper crust versus the more localized distribution in the subcrustal lithosphere, as inferred by such proxies as low P wave velocity and mafic magmatism

Fractal and Morlet-wavelet analyses of M ≥ 6 earthquakes in the South-North Seismic Belt, China

The M ≥ 6 earthquakes occurred in the South-North Seismic Belt, Mainland China (longitudes from 98 - 107°E and latitudes from 21 - 41°N) during 1900 - 2016 are taken to measure the multifractal dimensionsspatial distribution and time sequence of events and the dominant periods. The multifractal dimensions, Dq, are measured from the log-log plots of Cq(r) versus r and Cq(t) versus t, where Cq(r) and Cq(t) are the generalized correlation integrals for the epicentral distribution and time sequence of events, respectively. r and t are the epicentral distance and inter-event time, respectively, at positive q. The log-log plot of Cq(r) versus r shows a linear por­tion when log(rl) ≤ log(r) ≤ log(ru). The rl and ru values are, respectively, 120 and 560 km for M ≥ 6 events, 100 and 560 km for M ≥ 6.5 events, and 63 and 560 km for M ≥ 7 events. The rl value decreases with the lower-bound magnitude. Dq monotonically decreases with increasing q. The Dq values are between 1.618 and 1.426 for M ≥ 6 events, between 1.562 and 1.108 for M ≥ 6.5 events, and between 1.365 and 0.841 for M ≥ 7 events. The log-log plot Cq(t) versus t show a linear distribution when log(tl) ≤ log(t) ≤ log(tu), where tl and tu are, respectively, 5 and 50.1 years for M ≥ 6 events, 5 and 50.1 years for M ≥ 6.5 events, and 16 and 63.1 years for M ≥ 7 event, thus sug­gesting that the time sequences of earthquake in the study region are multifractal. The Dq values are between 0.830 and 0.703 for M ≥ 6 events, between 0.835 and 0.820 for M ≥ 6.5 events, and between 0.786 and 0.685 for M ≥ 7 events. The Morlet wavelet technique is applied to analyze the dominant periods of temporal variations in num­bers of yearly earthquakes for the three magnitude ranges, i.e., M ≥ 6, M ≥ 6.5, and M ≥ 7. The resultant dominant period is 2.94 years for M ≥ 6 events and cannot be evaluated for M ≥ 6.5 and M ≥ 7 events

Recent tectonic reorganization of the Nubia-Eurasia convergent 2 boundary heading for the closure of the western Mediterranean

: In the western Mediterranean area, after a long period (late Paleogene-Neogene) of Nubian northward subduction beneath Eurasia, subduction is almost ceased as well as convergence accommodation in the subduction zone. With the progression of Nubia-Eurasia convergence, a tectonic reorganization is therefore necessary to accommodate future contraction. Previously-published tectonic, seismological, geodetic, tomographic, and seismic reflection data (integrated by some new GPS velocity data) are reviewed to understand the reorganization of the convergent boundary in the western Mediterranean. Between northern Morocco, to the west, and northern Sicily, to the east, contractional deformation has shifted from the former subduction zone to the margins of the two backarc oceanic basins (Algerian-Liguro-Provençal and Tyrrhenian basins) and it is now active in the south-Tyrrhenian (northern Sicily), northern Liguro-Provençal, Algerian, and Alboran (partly) margins. Compression and basin inversion has propagated in a scissor-like manner from the Alboran (c. 8 Ma) to the Tyrrhenian (younger than c. 2 Ma) basins following a similar propagation of the subduction cessation and slab breakoff, i.e., older to the west and younger to the east. It follows that basin inversion is rather advanced in the Algerian margin, where a new southward subduction seems to be in its very infant stage, while it has still to properly start in the Tyrrhenian margin, where contraction has resumed at the rear of the fold-thrust belt and may soon invert the Marsili oceanic basin. GPS-derived strain rates higher in the Tyrrhenian margin than in the Algerian boundary suggest that this latter manner of contraction accommodation (contraction resumption at the rear of the orogenic wedge) is more efficient than subduction inception and basin inversion along newly-generated reverse faults (Algeria), but the differential strain rates may also be explained with the heterogeneous distribution of GPS stations. Part of the contractional deformation may have shifted toward the north in the Liguro-Provençal basin possibly because of its weak rheological properties compared with the area between Tunisia and Sardinia, where no oceanic crust occurs and seismic deformation is absent or limited compared with the adjacent strands of the Nubia-Eurasia boundary. The tectonic reorganization of the Nubia-Eurasia boundary in the study area is still strongly controlled by the inherited tectonic fabric and rheological attributes, which are both discontinuous and non-cylindrical along the boundary. These features prevent, at present, the development of long and continuous thrust faults. In an extreme and approximate synthesis, the evolution of the western Mediterranean is inferred as being similar to a Wilson Cycle in the following main steps: (1) northward Nubian subduction with Mediterranean backarc extension (since ~35 Ma); (2) progressive cessation, from west to east, of Nubian main subduction (since ~15 Ma); (3) progressive compression, from west to east, in the former backarc domain and consequent basin inversion (since ~8-10 Ma); (4) possible future subduction of former backarc basins

Tertiary-Quaternary subduction processes and related magmatism in the Alpine-Mediterranean region

During Tertiary to Quaternary times, convergence between Eurasia and Africa resulted in a variety of collisional orogens and different styles of subduction in the Alpine-Mediterranean region. Characteristic features of this area include arcuate orogenic belts and extensional basins, both of which can be explained by roll-back of subducted slabs and retreating subduction zones. After cessation of active subduction, slab detachment and post-collisional gravitational collapse of the overthickened lithosphere took place. This complex tectonic history was accompanied by the generation of a wide variety of magmas. Most of these magmas (e.g. low-K tholeiitic, calc-alkaline, shoshonitic and ultrapotassic types) have trace element and isotopic fingerprints that are commonly interpreted to reflect enrichment of their source regions by subduction-related fluids. Thus, they can be considered as ‘subduction-related’ magmas irrespective of their geodynamic relationships. Intraplate alkali basalts are also found in the region generally postdated the ‘subduction-related’ volcanism. These mantle-derived magmas have not been, or only slightly, influenced by subduction-related enrichment. This paper summarises the geodynamic setting of the Tertiary-Quaternary “subduction-related” magmatism in the different segments of the Alpine-Mediterranean region (Betic-Alboran-Rif province, Central Mediterranean, the Alps, Carpathian-Pannonian region, Dinarides and Hellenides, Aegean and Western Anatolia), and discusses the main characteristics and compositional variation of the magmatic rocks. Radiogenic and stable isotope data indicate the importance of continental crustal material in the genesis of these magmas. Interaction with crustal material probably occurred both in the upper mantle during subduction (‘source contamination’) and in the continental crust during ascent of mantle-derived magmas (either by mixing with crustal melts or by crustal contamination). The 87Sr/86Sr and 206Pb/204Pb isotope ratios indicate that an enriched mantle component, akin to the source of intraplate alkali mafic magmas along the Alpine foreland, played a key role in the petrogenesis of the ‘subduction-related’ magmas of the Alpine-Mediterranean region. This enriched mantle component could be related to mantle plumes or to long-term pollution (deflection of the central Atlantic plume and recycling of crustal material during subduction) of the shallow mantle beneath Europe since the late Mesozoic. In the first case, subduction processes could have had an influence in generating asthenospheric flow by deflecting nearby mantle plumes due to slab roll-back or slab break-off. In the second case, the variation in the chemical composition of the volcanic rocks in the Mediterranean region can be explained by “statistical sampling” of the strongly inhomogeneous mantle followed by variable degrees of crustal contamination

Numerical modelling and observations of nuclear-explosion coda wavefields

Frequency-dependent earthquake coda attenuation values are often reported; however such measurements usually depend on the types of the attenuation models employed. In this thesis, I use numerical modeling of Peaceful Nuclear Explosion (PNE) codas at far regional to teleseismic distances to compare two of such models, namely the conventional frequency-dependent attenuation with parameters (Q0, ¦Ç) defined by Qcoda(f) = Q0f¦Ç and frequency-independent effective attenuation (Qe) with geometrical attenuation (¦Ã). The results favour strongly the (¦Ã, Qe) model and illustrate the mechanisms leading to apparent Qcoda(f) dependencies. Tests for variations of the crustal velocity structures show that the values of ¦Ã are stable and related to lithospheric structural types, and the inverted Qe values can be systematically mapped into the true Swave attenuation factors within the crust. Modeling also shows that ¦Ã could increase in areas where relatively thin attenuating layers are present within the crust; such areas could likely be related to younger and active tectonics. By contrast, when interpreted by using the traditional (Q0,¦Ç) approach, the synthetic coda shows a strong and spurious frequency dependence with ¦Ç ¡Ö 0.5, which is also similar to many published observations. Observed Lg codas from two Peaceful Nuclear Explosions located in different areas in Russia show similar values of ¦Ã ¡Ö 0.75¡¤10-2 s-1, which are also remarkably close to the independent numerical predictions in this thesis. At the same time, coda Qe values vary strongly, from 850 in the East European Platform to 2500 within the Siberian Craton. This suggests that parameters ¦Ã and Qe could provide stable and transportable discriminants for differentiating between the lithospheric tectonic types and ages, and also for seismic coda regionalization in nuclear-test monitoring research

Tectonic evolution and current deformation of the NW Sicily Channel and the Lampedusa Plateau based on multi-resolution seismic profiles analysis

The thesis deals with the tectonic evolution and the active deformation of the NW part of the Sicily Channel and of the offshore area around the Lampione-Lampedusa islands based on the analysis of multi- and single-channel seismic reflection profiles calibrated with well-log data

Contributions to the Study of Lithospheric Deformation and Seismicity in Stable Continental Regions

Recently, the field of geophysics has seen increasing recognition of the unique character of deformation and seismicity in stable continental regions (SCRs). However several important questions remain understudied. What controls the locations of earthquakes in SCRs? How well do observations, in SCRs, of elastic strain accumulation and release correlate with each other? How well do they correlate with stresses and geological proxies for rheological variation? The ultimate goal of this study was to better understand stable continental regions like southern Africa, where large earthquakes occur despite not being near plate boundaries, for example the 2017 Mw 6.5 earthquake in Moiyabana, Botswana. One way of studying the stress and strain in stable continental regions is by understanding the surface deformation of the region. This deformation is easily studied using global navigation satellite system (GNSS) velocity data. One of the biggest difficulties when it comes to GNSS data is that it isn't collected on a regular grid, but rather as irregular data points that need to be interpolated. This research investigated multiple interpolation methods and recommended two methods that best replicate the original velocity field (using a well populated dataset from Southeast Asia). These interpolated GNSS data can then be used to determine deviatoric strain in a region, which can in turn be fed into numerical stress models. However, limited GNSS data exist across southern Africa, and therefore topographic data was used to calculate the gravitational potential energy, and in turn the body stress and deviatoric stress for the region. This study also investigated how this deviatoric stress (or deviatoric strain) can be more accurately calculated on a spherical rather than a flat surface, which is particularly important over large study areas. Across southern Africa, data show that deviatoric stress lined up with stress data within mobile belts. This suggests that in these weaker mobile belt crust (such as the Namaqua-Natal and Damara-Chobe belts), gravitational collapse is the dominant driver of deformation, which is in line with conclusions that have been made in previous literature. In other regions, deviatoric stress vectors and stress data do not coincide and therefore there are other forces at play. These observations are obviously restricted by limited data coverage; it remains an open question if areas that have increased deviatoric stress due to gravitational collapse, which are also aligned with the orientation of weak zones, will have elevated strain in the long term