Crustal Heterogeneity Research Articles

Crustal Structure of Southern California from Velocity and Attenuation Tomography

Recent studies of Southern California have employed velocity and attenuation tomography to elucidate two aspects of crustal structure. Low-frequency velocity tomography reveals large-scale heterogeneity that can be approximated by a discrete set (𝐾 ≤ 10) of tectonic regions, each characterized by an isostatically balanced lithospheric column reflecting the composition and tectonic history of the region. The boundaries between tectonic regions are typically localized structures expressed at the surface by major faults, topographic fronts, and geochemical transitions. The efficacy with which tomographic regionalization (TR) represents the subsurface geography of major tectonic units in Southern California indicates that the TR idealization works surprisingly well in this part of the Pacific‐North America plate boundary. High‐frequency attenuation tomography in Southern California supports the hypothesis that the decay of P and S waves propagating through the upper and middle crust is dominated by elastic scattering from small-scale heterogeneities of high fractal dimension, rather than by anelastic dissipation. The amplitude of the scattering varies systematically among the tectonic regions and correlates with the degree of recent faulting and deformation within a region. These results motivate speculation that small-scale crustal heterogeneities responsible for high-frequency attenuation are generated within the seismogenic zone through a dynamic interplay between distributed deformation and localized fracturing.

Seismic Evidence for Velocity Heterogeneity Along ∼40 Ma Old Oceanic Crustal Segment Formed at the Slow‐Spreading Equatorial Mid‐Atlantic Ridge From Full Waveform Inversion of Ocean Bottom Seismometer Data

AbstractIn slow spreading environments, oceanic crust is formed by a combination of magmatic and tectonic processes. Using full waveform inversion applied to active‐source ocean bottom seismometer data, we reveal the presence of a strong lateral variability in the 40–48 Ma old oceanic crust formed at the slow‐spreading Mid‐Atlantic Ridge in the equatorial Atlantic Ocean. Over a 120 km‐long section between the St Paul fracture zone (FZ) and the Romanche transform fault (TF), we observe four distinct 20–30 km long crustal segments. The segment affected by the St Paul FZ consists of three layers, an ∼2 km thick layer with a P‐wave velocity <6 km/s, a 1.5 km thick middle crust with a velocity of 6–6.5 km/s, and an underlying layer where velocity is ∼7 km/s, representing the lower crust. The segment associated with an abyssal hill morphology contains a high velocity of ∼7 km/s at 2–2.5 km below the basement, indicating the presence of primitive gabbro or serpenized peridotite. The segment associated with a low basement morphology seems to have 5.5–6.5 km/s velocity starting near the basement extending down to ∼4 km depth, indicating chemically distinct crust. The segment close to the Romanche TF, a velocity 4.5–5 km/s near the seafloor increasing to 7 km/s at 4 km depth indicates a magmatic origin. The four distinct crustal segments have a good correlation with the overlying seafloor morphology. These observed strong crustal heterogeneities could result from alternate tectonic and magmatic processes along the ridge axis, possibly modulated by thermal and/or chemical variations in the mantle during their formation along the ridge segment.

Evolution of the transtensional Barreirinhas pull-apart system in the Brazilian Equatorial margin and its correlation with the African conjugate counterpart

The Barreirinhas pull-apart system encompasses marginal basins in divergent and transform margin segments in the central sector of the Brazilian Equatorial Margin and its African conjugate counterpart. This ancient pull-apart system evolved through transtensional strike-slip motion within a highly heterogeneous crystalline basement affected by multiple rift phases. The geometry and development of pull-apart structural elements during the final rifting phase before continental breakup and the mechanisms and extent to which they were influenced by preexisting crustal heterogeneities are comprehensively addressed using an extensive database of potential field (magnetic and gravity) and 2D seismic reflection data. We also assess the lithospheric thermomechanical conditions and their influence on transtensional extension throughout Curie Point Depth, Heat Flow, and Moho depth, derived from potential field data and published seismological models. Plate reconstruction of Brazilian and African equatorial margins based on gravity patterns and comparison with sandbox analog models allow a 3D synoptic model to reveal the Barreirinhas pull-apart system evolution during the Equatorial Atlantic opening. During the rift phase I, the location of major grabens was controlled by favorably oriented Neoproterozoic shear zones, while the cooler, stronger, and thicker crust beneath cratonic areas formed the western barrier to strike-slip rift activity during rift phase II. This same geological domain anchored the onset of the pull-apart system in the last rift phase III, whose principal displacement zones developed along the extensive oceanic fracture zones linked by sigmoidal fault systems. Toward the end of the rift phase, a large asymmetric lozangle to a lazy-Z-shaped, pull-apart basin developed above low overlapping ∼90°, releasing stepover in oblique transtensional strike-slip motion.

Geology Matters for Antarctic Geothermal Heat

AbstractGeothermal heat plays a vital role in Antarctic ice sheet stability. The continental geothermal heat flow distribution depends on lithospheric composition and ongoing tectonism. Heat‐producing elements are unevenly enriched in the crust over deep time by various geological processes. The contribution of crustal heat production to geothermal heat flow is widely recognized; however, in Antarctica, crustal geology is largely hidden, and its complexity has frequently been excluded in thermal studies due to limited observations and oversimplified assumptions. Li and Aitken (2024), https://doi.org/10.1029/2023GL106201 take a significant step forward, focusing on Antarctic crustal radiogenic heat. Utilizing gravity inversion and rock composition data, they show that the crustal heterogeneity introduces considerable variability to heat flow. However, modeling crustal heat production proves challenging because it lacks distinct associations with geophysical observables and has a narrow spatial association. Robust quantification of geothermal heat production and heat flow must incorporate explicit aspects of geology.

Compositional Diversity of Early Mesozoic Granites in South Qinling: Derivation from Heterogenous Basement Rocks in the Orogenic Belt

Granitic rocks forming in the syn- to post-orogenic stages can trace the compositional and structural complexity of the crust beneath an orogenic belt. The Qinling orogenic belt undertook multiple stages of tectonics and magmatism, resulting in the multifaceted evolution and compositional diversity of the crust. In the present study, the Guangtoushan and Miba plutons in South Qinling were chosen to reveal the crustal heterogeneity in study area via isotopic geochemistry and zircon geochronology. The Guangtoushan pluton was emplaced between ~215 Ma and ~202 Ma and the Miba pluton formed at ~213 Ma, as constrained by zircon U-Pb isotopic dating. Granitic rocks of the Miba pluton are characterized by amphibole bearing and homogeneous composition, with relatively depleted Sr-Nd isotopic compositions (initial 87Sr/86Sr values of 0.7060 to 0.7084 and initial εNd values of −5.4 to −9.5) and high Pb isotopic values. The Guangtoushan pluton contains muscovite and complex inherited zircon grains and has variable Sr-Nd isotopic composition (initial 87Sr/86Sr values of 0.7050 to 0.7091 and initial εNd values of −4.5 to −12.9) and low Pb isotopic values. Felsic magmas of the Guangtoushan pluton should be derived mainly from meta-sedimentary rocks beneath South Qinling, while the Miba pluton originated primarily from partial melting of meta-igneous rocks. The compositional diversity recorded in the Early Mesozoic plutons was caused by the heterogeneous crust, and partial melting was induced by heating of the up-welling asthenosphere in a post-collision setting.

Three‐Dimensional Electrical Resistivity Structure Beneath a Strain Concentration Area in the Back‐Arc Side of the Northeastern Japan Arc

AbstractIntraplate earthquakes occur more frequently in the Japanese islands than in other regions. Large intraplate earthquakes in the island arc preferentially occur in strain concentration zones detected by geological and geodetic studies. Crustal heterogeneity plays a crucial role in generating large intraplate earthquakes and strain concentrations. Thus, we elucidated the crustal heterogeneities beneath a strain concentration area on the back‐arc side of the northeastern Japan Arc based on electrical resistivity, which is sensitive to weak zones in the crust. By deploying wideband magnetotelluric surveys, we revealed the three‐dimensional electrical resistivity structure in the crust, suggesting the coexistence of two different types of strain‐concentration mechanisms in the strain‐concentration area. The shallow conductive layers and lower‐crustal conductors appear to act as low‐elastic‐modulus and low‐viscosity areas, respectively, and are responsible for the strain concentration. We found shallow and lower‐crustal conductors in the strain concentration zone revealed by geological studies. Those conductive areas are considered to act as mechanically weak zones and cause stress loading on the brittle pars of pre‐existing faults, resulting in large intraplate earthquakes. The resultant earthquakes presumably contribute to strain accumulation on a geological timescale. In addition, a spatial correlation between the epicenters of large intraplate earthquakes and edges of lower‐crustal conductors implies the contribution of the fluids in the lower crust to the generation of large earthquakes. We also identified vertical conductors ranging from the lower crust to Quaternary volcanoes, which may indicate trans‐crustal magmatic systems under these volcanoes.

Imaging of moho topography with conditional generative adversarial network from observed gravity anomalies

Accurate estimation of Moho topography plays a crucial role in understanding Earth’s structure, geodynamic processes, and resource exploration. This study presents a novel approach that utilizes conditional Generative Adversarial Networks (cGAN) to reveal Moho topography based on observed gravity anomalies. Synthetic training datasets of Moho topography were generated using the FFT filtering method due to the scarcity of true datasets. Spherical prism-based forward gravity modeling was employed to evaluate the resulting gravity anomalies. We compared the performance of our developed deep learning algorithm cGAN (conditional Generative Adversarial Networks) with a traditional inversion technique using various synthetic datasets, and a real case study in southern peninsular India, a geologically diverse region comprising ancient continental tectonic blocks. Bott’s inversion scheme was employed as a verification method for the Moho surface estimation using the presented deep learning model. Using spherical prism-based forward gravity modeling, observed gravity anomalies were corrected for multiple factors such as topography, bathymetry, sediments, crustal heterogeneities, and mantle heterogeneities. By removing these effects, we isolated the gravity contribution solely related to pure Moho undulation. The mean Moho depth and density contrast between the crust and mantle were derived from seismic constraints for improving estimation accuracy. The findings demonstrate the potential of the cGAN and spherical prism-based gravity modeling approach in accurately estimating the Moho topography, offering insights into Earth’s subsurface structures and enhancing our understanding of geodynamic processes and resource exploration efforts.

Estimation of Source and Spectral Decay Parameters for Local Earthquakes in Siang Region of Arunachal Himalaya and Its Implication to the Tectonics and Crustal Heterogeneity

Estimation of Source and Spectral Decay Parameters for Local Earthquakes in Siang Region of Arunachal Himalaya and Its Implication to the Tectonics and Crustal Heterogeneity

Detrital isotopic record of a retreating accretionary orogen: An example from the Patagonian Andes

Abstract U-Pb zircon geochronology and isotopic records have played an influential role in our understanding of convergent margin dynamics. Orogenic cyclicity models link tectonic regimes with magmatic isotopic signatures in advancing orogens, relating compressional regimes with evolved signatures and extension with juvenile signatures; however, such frameworks may not apply for retreating orogens, which commonly produce substantial crustal heterogeneities during backarc rifting and ocean spreading. We explore the Mesozoic to Cenozoic Patagonian Andes tectonic evolution, combining U-Pb zircon ages, bulk rock εNd, and new detrital zircon εHf from the retroarc basin to understand the associated magmatic arc evolution during retreat and advance of the margin. Our results reveal a protracted phase of isotopically juvenile magmatism between 150 and 80 Ma, which began during backarc extension and persisted long after the margin switched to a contractional regime. We propose that the prolonged juvenile isotopic trend started mainly due to trenchward migration of the arc during backarc extension (150–120 Ma) and persisted due to partial melting of underthrusted juvenile attenuated and oceanic crust during backarc basin closure (120–80 Ma). This interpretation implies that tectonic stress alone does not predict isotopic trends, and factors like assimilation or the composition of underthrusted crust are important controls on magmatic isotopic composition, especially in retreating and transitional orogens.

Fold localization at pre-existing normal faults: field observations and analogue modelling of the Achental structure, Northern Calcareous Alps, Austria

Abstract. Within the Northern Calcareous Alps (NCA) fold-and-thrust belt of the Eastern Alps, multiple pre-shortening deformation phases have contributed to the structural grain that controlled localization of deformation at later stages. In particular, Jurassic rifting and opening of the Alpine Tethys led to the formation of extensional basins at the northern margin of the Apulian plate. Subsequent Cretaceous shortening within the Northern Calcareous Alps produced the enigmatic Achental structure, which forms a sigmoidal transition zone between two E–W-striking major synclines. One of the major complexities of the Achental structure is that all structural elements are oblique to the Cretaceous direction of shortening. Its sigmoidal form was, therefore, proposed to be a result of forced folding at the boundaries of the Jurassic Achental basin. This study analyses the structural evolution of the Achental structure through integrating field observations with crustal-scale physical analogue modelling to elucidate the influence of pre-existing crustal heterogeneities on oblique basin inversion. From brittle–ductile models that include a weak basal décollement, we infer that oblique shortening of pre-existing extensional faults can lead to the localization of deformation at the pre-existing structure and predicts thrust and fold structures that are consistent with field observations. Consequently, the Achental low-angle thrust and sigmoidal fold train was able to localize at the former Jurassic basin margin, with a vergence opposite to the controlling normal fault, creating the characteristic sigmoidal morphology during a single phase of NW-directed shortening.

Deciphering the low-frequency seismic signals in the Weiyuan Shale gas field: implications for reservoir and structural heterogeneity

SUMMARY Hydraulic fracturing (HF) often stimulates the local earthquake productivity which provides a unique opportunity to characterize the crustal heterogeneities, reservoir properties and fluid injection effects. However, the velocity models acquired solely based on the arrival time records are often undermined due to the seismic network coverage and interpolation techniques. Instead, we adopt the waveform-based approach to apprehend; (1) structural heterogeneities, (2) reservoir distribution and (3) signatures of the injected fluid in the Weiyuan shale gas field. We categorize the waveforms into dominant high and low frequencies based on the qualitative inspection and frequency index analysis of the seismic waveforms. We first inspect the waveform to access the potential controlling mechanisms (source, site and path effects) at both single and multiple stations in different azimuthal orientations. As a result, we find the path effect as a dominant factor to influence the waveform characteristics, for example S-wave amplitude, and frequency. Subsequently, to localize the path effect, we conduct an in-depth examination of events within 10 km of each seismic station and classify the waveform records using their frequency indices. Notably, certain stations record a significant proportion of low-frequency waveforms (LFWs, up to 20 per cent), while others have limited occurrences (∼1 per cent) indicating suspected anomalous zones. Afterward, we identify two suspected anomalous zones based on LFWs intensity and ray tracing map. Both zones are in close proximity to fault zones and preserved reservoirs with no HF activities, where fault damage zones or the fluid-rich reservoir may contribute to our observed LFWs.

Crustal Heterogeneity of Antarctica Signals Spatially Variable Radiogenic Heat Production

AbstractGeothermal heat flow (GHF) is a key basal boundary condition for Antarctic ice‐sheet flow. Large‐scale variations are resolved by several recent models but knowledge of the smaller‐scale variations, crucial for ice sheet dynamics, is limited by unresolved variations in crustal radiogenic heat production. To define this at continent‐scale we use 3D gravity inversion constrained by seismic Moho estimates to identify variations in crustal composition and geometry beneath thick ice. Geochemically‐defined empirical relationships between density and heat production capture the global average trend and its variability, and allow to estimate from upper‐crust density spatial variations in radiogenic heat production. Significant variations are observed typically 1.2–1.6 μW/m3, and as high as 2 μW/m3 in West Antarctica. The contribution to GHF from these heat‐production variations is similarly variable, typically 16–24 mW/m2 and up to 60 mW/m2. The mapped variations are significant for correctly representing GHF in Antarctica.

Crustal heterogeneity effects on coseismic deformation: numerical simulation of the 2008 MW 7.9 Wenchuan earthquake

Coseismic deformation of large earthquakes causes significant property damages and fatalities, which requires quantitative research of multiple disciplines such as geodesy, geological investigation, seismic tomography, and seismic dislocation theory. The finite element method accounts for material heterogeneity and geometric complexity, making it suitable for studying the coseismic deformation of large earthquakes. This paper develops a parallel elastic finite element program that utilizes split nodes and high-performance parallel computing technology on the FELAC software platform to study the coseismic deformation of large earthquakes. We verify the accuracy of the parallel elastic finite element program by comparing its results with the analytical solutions from seismic dislocation theory for four ideal earthquake cases. Finally, we established parallel elastic finite element models to study the coseismic deformation of the 2008 Wenchuan earthquake. The simulation results are consistent with the GPS and InSAR data. Coseismic surface deformation results are significantly influenced by medium regional heterogeneity with different layered structures besides the Longmenshan fault. The finite element program lay the foundation for the inversion of the coseismic fault rupture process based on the heterogeneous medium model and complex geometric model.

Lithotectonic architecture of the Paleoproterozoic intracratonic Kaladgi rift basin, southwestern India: Inference from a magnetotelluric study

Peninsular India is traversed by several major Paleoproterozoic rift basins, among which the Kaladgi rift basin is a prominent one. Here, we present results from Magnetotelluric (MT) studies in the eastern section of this rift basin. The data were acquired along a ∼120 km long profile passing through Belavanik-Bagalkot-Ukkali regions covering the Paleoproterozoic intracratonic Kaladgi rift basin, the Cretaceous Deccan volcanic province and the Archean gneisses/granites/greenstones of the Dharwar craton. The lithospheric electrical structure was obtained from joint inversion of TE- and TM-modes data using a 2-D nonlinear conjugate gradient algorithm. The results reveal that the Kaladgi sediments are ∼500–1000 m in thickness overlying a highly fractured crystalline basement. Proterozoic sediments are exposed beneath the Deccan basalt cover and in the shallow fractured zones. The conductive-resistive transitions in the model at the crustal depths correspond to major tectonic faults and deformed crustal rocks. The transition boundary at ∼50 km depth defines the electrical Moho. At the same time, a ∼25 km wide, steep conductive feature from ∼75 km to upper mantle lithosphere depths (∼180–200 km) is interpreted as the trace of the Chitradurga Suture Zone (CSZ). The observed crustal heterogeneity, conductivity variations, and resistive lithosphere are attributed to the geological history of the region, including Archean collisional events and subsequent reactivation associated with the Reunion hotspot. The crustal conductors are associated with mafic intrusions from the underplated basalts, and also the high conductivity may be due to the sulfide mineralization in the fault zones. The MT study provides valuable insights into the deep tectonic framework of the region.

Crustal structure of the Xisha block and implications for along-strike segmentation in the northern South China Sea continental margin

Continental blocks complicate the pre-rift lithospheric inheritance of passive continental margins, hence exert unique influences on rifting process. Along-strike rifting contrasts of the continental margins may reflect the pre-rift crustal heterogeneity and segmented rifting processes. A wide-angle seismic profile acquired on the Xisha Block, the northwestern margin of the South China Sea (SCS), provides the first constrains on the along-strike crustal structures of the continental block. Seismic imaging of the profile reveals a continental crust with relatively uniform thickness of about 22 km and the thickness ratio of the upper crust to lower crust at 0.8. Crustal stretching factor is approximately 1.5, and high-velocity anomalies are absent in the lower crust. Crustal velocities of the upper and lower crust are approximately 5.4–6.5 km/s and 6.6–6.9 km/s, respectively. The results indicate that the Xisha Block underwent slight and uniform extension during rifting, and the crust has not suffered from strong magmatic activities that formed extensive magmatic underplating. As a rigid continental block, the Xisha Block has retarded the rifting onset and shortened the process from lithospheric rifting to breakup in the northwestern SCS margin. Hypothetic models for the rheological heterogeneity of the inherited crust are proposed tentatively to interpret the contrasting distribution of the high velocities in the lower crust along the northern SCS margin.

Gravity evidence for a heterogeneous crust of Mercury

We modeled gravity data to explore Mercury’s internal structure and show the presence of crustal heterogeneities in density. We first evaluated the lithospheric flexure occurring in the spherical harmonic degree range 5–80, according to the flexural isostatic response curve. We thus estimated a mean elastic lithosphere thickness of about 30 pm 10 km and modeled the crust-mantle interface, which varies from 19 to 42 km depth, according to a flexural compensation model. The isostatic gravity anomalies were then obtained as the residual field with respect to the contributions from topography and lithospheric flexure. Isostatic anomalies are mainly related to density variations in the crust: gravity highs mostly correspond to large-impact basins suggesting intra-crustal magmatic intrusions as the main origin of these anomalies. Isostatic gravity lows prevail, instead, above intercrater plains and may represent the signature of a heavily fractured crust.

Convection and segregation in heterogeneous orogenic crust with a VOF method – II: how to form migmatite domes

SUMMARY Migmatites and granitic-gneisses exhumed in Archean to Phanerozoic segments are former partially molten crustal roots, display typical domes structures ranging in size from kilometres to decakilometres, and are often interpreted as resulting from the development of diapiric or convective gravitational instabilities. In previous work (part I), we determined various regimes of gravity-driven segregation, by considering a thick continental crust heated from below and containing melt related heterogeneities. These heterogeneities, represented by inclusions of distinct densities and viscosities with respect to the ambient partially molten material, can be entrained into convection cells (in the ‘suspension’ and ‘layering’ regimes) and/or accumulate as clusters (in the ‘layering’ and ‘diapirism’ regimes). Here we further investigate the specific conditions that allow for the formation and preservation of domes resulting from diapirism at the top of convective cells. We show that both the cessation of basal heating and the freezing of the buoyant inclusions density favour their stacking and preservation at ca. 15 km depth, within about 10 Myr. The buoyant inclusions form domes, 5–20 km in size, that also record several convective cycles at velocities ranging from 0.5–4 cm yr−1. 3-D models demonstrate their radial geometrical nature. The influence of the size and concentration of the inclusions is also assessed, complementing the characteristics of crustal heterogeneity in driving its differentiation and the formation of migmatite domes.

Heterogeneous Strain Distribution in the Malawi (Nyasa) Rift, East Africa: Implications for Rifting in Magma‐Poor, Multi‐Segment Rift Systems

AbstractHalf‐graben basins bounded by border faults typify early‐stage continental rifts. Deciphering the role that intra‐rift faults play in rift basin development is challenging as patterns of early‐stage faulting are commonly overprinted by subsequent deformation; yet the characterization of these faults is crucial to understand the fundamental controls on their evolution, their contribution to rift opening, and to assess their seismic hazard. By integrating multiple offshore seismic reflection data sets with age‐dated drill core, late‐Quaternary and cumulative faulting patterns are characterized in the Central and South Basins of the Malawi (Nyasa) Rift, an active, early‐stage rift system. Almost all intra‐rift faults offset a late‐Quaternary lake lowstand surface, suggesting they are active and should be considered in hazard assessments. Fault throw profiles reveal sawtooth patterns indicating segmented slip histories. Observed extension on intra‐rift faults is approximately twice that predicted from hanging wall flexure of the border fault, suggesting that intra‐rift faults accommodate a proportion of the regional extension. Cumulative and late‐Quaternary throws on intra‐rift faults are correlated with throw measured on the border fault in the Central Basin, whereas an anticorrelation is observed in the South Basin. Viewed in a regional context, these differences do not relate solely to the proposed southward younging of the rift. Instead, it is inferred that the distribution of extension is also influenced by variations in lithospheric structure and crustal heterogeneities that are documented along the rift axis.

Delineation of partial melts and crustal heterogeneities within the crust beneath Kumaon Himalaya, India from Lg wave attenuation

The crustal seismic attenuation or the Q structure is studied by using the Fourier spectra of Lg-wave along the Tanakpur- Dharchula- Dharma transect in the Kumaon Himalaya. The 1 Hz Lg Q (Q0) values are computed between different pairs of two stations and the observed values are later utilized to calculate the lateral variation in the Q0 values by following a back projection algorithm. This computation of Q0 values utilizes five regional distance earthquakes having moment magnitude (Mw) ≥ 4.0, which lie along the great circle path of the transect. Three of the five earthquakes occurred in the Tibetan plateau and the and the others occurred to the southwest on the Indian shield and are well recorded at all the 32 broadband seismographs operated between September 2018 and March 2022. The estimate Qo values range from 63 ± 2 and 203 ± 25, with the lowest value in the Lesser Himalaya and the highest across part of the Indo Gangetic Plain and Siwalik Himalaya. The Q0 model has low values ∼200 along the profile in the Indo Gangetic Plain and the Siwalik Himalaya, and are correlated with 2–5 km thick sedimentary layers below the Himalaya and the adjoining Indo-Gangetic Plain. We observe two distinctly different Q0 values to the northeast in the Lesser Himalaya tectonic unit. The region lying between the South Almora Thrust (SAT) and the Berinag Thrust (BT) shows extremely low Q0 values (∼60) but increases further north towards the Vaikrita Thrust (VT) to ∼200. The possible explanation for observing such huge variation of the Q0 values within a single tectonic unit may be the presence of fluid rich ramp structures, which introduces crustal heterogeneities and traps the aqueous fluids or partial melts lying within the crust. The Lg Q0 values decrease to the North and become ∼166 for station pairs in the Higher Himalaya and Tethys Himalaya tectonic units. The low Q0 values observed in this region may be correlated with low viscous partial melts in the form of Miocene leucogranite plutons, which resulted out of the Indo-Asian collision. The attenuation structure along the profile in the Kumaon Himalaya can be used to estimate ground motions of future earthquakes in the area and can contribute to seismic hazard assessment in the Himalaya and neighbouring regions.

Three-dimensional P-wave model of the shallow crustal structure, a complementary method for detecting a trapped hydrocarbon: a case study in the DehDasht region, SW Iran

Abstract Local P-wave tomography is an efficient method to study geologically complex areas where the seismic exploration methods are not ideal for unraveling the shallow crustal heterogeneity due to the great thickness of evaporitic deposits. Despite the complex geological features in the salt-rich DehDasht region, SW Iran, we used >11 000 micro-earthquake events, which have been recorded by a temporary seismic network (deployed between 18 October 2016 and 1 July 2017), to derive the three-dimensional velocity structure based on the first arrival time. We selected a subset of events (1571 micro-earthquakes) by various strict criteria for our processing, and then the 1D velocity model was calculated by the computer program VELEST. Afterward, the 3D initial model of the inversion procedure with 1.5-km horizontal and 1-km deep intervals was parametrized using the calculated 1D model. Finally, the observed data (first arrival P-wave traveltimes and events locations) was inverted with an optimum regularization parameter and iteration using the computer program SIMULPS14. Our tomographic results indicate the DehDasht Basin as a relatively low-velocity zone filled out dominantly by the Gachsaran Formation and surrounded by the high-velocity Asmari-Pabdeh-Sarvak Formations. The basin has a bowl shape that is elongated in the NW–SE direction or an oval on a horizontal view. The depth of the basin varies between 3 and 5 km and contains many folding-faulting systems, which lead to locally low-velocity patches. Moreover, some evaporate deposits, which are overlying the Gachsaran Formation, emerge as a thin low-velocity layer (e.g. Aghajari, etc.).