Crustal Model Research Articles

The Thermal Structure of Central Te Riu a Māui/Zealandia Continent as Determined From the Depth to Base of Magnetic Sources

AbstractTo investigate the thermal structure of central Te Riu a Māui/Zealandia continent we estimate the depth to base of magnetic sources (DBMS), using an updated compilation of magnetic data. We calculate thousands of power spectral density estimates and solve an analytic solution for the fractal distribution of magnetization using a Bayesian method constrained by a sediment thickness model. Our DBMS results are broadly similar to a global model yet show complex relationships with independent models of crust and elastic thicknesses, and depth to base of seismicity (D90). Across central Zealandia, DBMS occurs above, coincident with, and below the Moho reflecting that factors additional to temperature influence DBMS. Away from subduction zones crustal thickness likely modulates DBMS while plate margins show varying influence dependent on margin characteristics. In the Taupō Volcanic Zone we interpret DBMS to reflect Curie point temperature at an average depth of 9–11 km. Error analysis shows that the choice of prior influences the posterior standard deviation and can result in counter intuitive error versus DBMS relationships. We use the DBMS as a basal temperature isotherm in heat flow models to compare with heat flow derived from borehole temperature measurements. Using the standard Curie point temperature of 580°C we find that DBMS derived heat flow models over estimate heat flow in areas underlain by plutonic rocks. We propose this reflects a higher titanomagnetite content in these rocks, although detailed studies of the magnetic mineralogy of Zealandia's plutonic rocks are required.

A New Lithospheric Density and Magnetic Susceptibility Model of Iran, Starting From High‐Resolution Seismic Tomography

AbstractWe propose a new model for the crust and upper mantle in Iran by joint inversion of gravity and magnetic fields, constrained with a seismic tomography model. We then calculate shear modulus from the Vs velocities and densities. The crust and mantle tomography model is first converted to a density cube through empirical and petrological velocity‐density relations. The starting susceptibility is assigned to a two‐layer homogeneous model, above a heat flow‐derived Curie depth. Considering the uncertainties in the density‐velocity relations, and the starting layered susceptibility variation, we refine the model by a constrained inversion of the gravity and magnetic fields with a Bayesian approach, producing the final 3D density and susceptibility model. The area is tectonically active with high seismicity and active faulting which are regulated by the crustal density and rigidity variations. Higher rigidity matches lower seismicity and extended deserts and basins, suggesting the control of their development. The Neo‐Tethys suture, extending ∼1,500 km long, as well as the Paleo‐Tethys suture match crustal scale density variations, defining characteristics of the lighter Arabian plate and denser Eurasian crust. The South Caspian Basin is enigmatic, due to focusing on the seismicity along all its borders, but with relatively low average rigidity, which is contrary to what is observed for Iran, where the reduced rigidity correlates with higher seismicity. The 3D density model will be useful for numerical geodynamic models and obtaining geologic inferences from the crustal‐scale units.

Reference seismic crustal model of the Dinarides

Abstract. Continental collision zones are structurally one of the most heterogeneous areas intermixing various different units within a relatively small space. A good example of this is the Dinarides, a mountain chain situated in the central Mediterranean, where thick carbonates cover older crystalline basement units and remnants of subducted oceanic crust. This is further complicated by the highly variable crustal thickness ranging from 20 to almost 50 km. In terms of spatial extension, this area is relatively small but covers tectonically differentiated domains making, any seismic or geological analysis complex, with significant challenges in areas that lack seismic information on crustal structure. Presently there is no comprehensive 3D crustal model of the Dinarides (and surrounding areas). Using the compilations of previous studies and employing kriging interpolation, we created a vertically and laterally varying crustal model defined on a regular grid for the wider area of the Dinarides, also covering parts of Adriatic Sea and the SW part of the Pannonian Basin. The model is divided by three interfaces, Neogene deposit bottom, carbonate rock complex bottom and Moho discontinuity, with seismic velocities (P and S waves) and density defined at each grid point. To validate the newly derived model, we calculated travel times for an earthquake recorded on several seismic stations in the Dinarides area. The calculated travel times show significant improvement when compared to the simple 1D model used for routine earthquake location in Croatia. The model derived in this work represents the first step towards improving our knowledge of the crustal structure in the complex area of the Dinarides. We hope that the newly assembled model will be useful for all forthcoming studies (e.g., as a starting model for seismic tomography, as a model for earthquake simulations) which require knowledge of the crustal structure.

Bayesian inversion for modelling 3-D density structures in the eastern margin of Bayan Har block and its tectonic implications

SUMMARY Recently, the eastern margin of the Tibetan plateau has experienced three moderate earthquakes. Deep crustal structure is a critical factor for understanding the seismotectonic environment and deformations in this region. Although several multidisciplinary geophysical approaches have been used to develop crustal structure models for this area, substantial inconsistencies still persist in the results due to the intricate structural properties of crust. In this study, we introduce a novel approach to construct a 3-D crustal model with assimilated density structure (MADS) at a resolution of 0.2° × 0.2° × 5 km. This model is built by integrating data from various source, including multisource gravity anomalies and diverse available crustal structure reference models. We use a model-assimilated gravity inversion process with Bayesian parameter optimization. First, we combine the terrestrial gravity profiles data set with published global gravity field models to create a data set rich in reliable high-frequency anomaly information. Secondly, we incorporate three reference models derived from different seismological methods as prior constraints for our model. These models encompass seismic tomography, surface wave dispersion and receiver function data. We optimize the hyper-parameters of these constraints using the Bayesian criterion. The results demonstrate that the MADS not only captures significant changes in the crustal density but also discerns subtle variations in the upper and middle crust, thereby providing detailed insights into the morphologies of major faults. For instance, the central section of the Longmenshan fault is revealed as a high-angle deep thrust feature, while the frontal section of the Longmenshan fault appears as a low-angle mid-deep thrust feature, and the Xianshuihe fault exhibits a vertical deep subduction feature. Additionally, our findings indicate a correlation between the locations of moderate-to-large earthquakes in this region and the high density-gradient zones or asperities with high density within MADS. We believe that the insights into density characteristics offered by the new MADS model can shed light on the study of asperities associated with recent moderate earthquakes and enhance our understanding of deformation in this region.

Using Tidally‐Driven Elastic Strains to Infer Regional Variations in Crustal Thickness at Enceladus

AbstractConstraining the spatial variability of the thickness of the ice shell of Enceladus (i.e., the crust) is central to our understanding of the internal dynamics and evolution of this small Saturnian moon. In this study, we develop a new methodology to infer regional variations in crustal thickness using measurements of tidally‐driven elastic strain that could be made in the future. As proof of concept, we recover thickness variations from synthetic finite‐element crustal models subjected to diurnal eccentricity tides. We demonstrate recovery of crustal thickness to within ∼2 km of true values across the crust with ∼10% error in derived spherical harmonic coefficients at degrees l ≤ 12. Our computed uncertainty is significantly smaller than the inherent ∼10 km ambiguity associated with crustal thickness derived solely from gravity and topography measurements. Therefore, future measurements of elastic strain can provide a robust approach to probe crustal structure at Enceladus.

Late Cryogenian Circum-Rodinia Syn-Subduction Extension: Insights from Highly Fractionated S-Type and A-Type Granitoids in the Northern Tarim Craton

Late Cryogenian–Ediacaran magmatism represents the latest Precambrian tectonothermal event in the Tarim Craton. However, its geodynamic setting and geological significance are controversial. Here, we report the geochronology, whole-rock geochemistry, and Sr-Nd-Hf isotopic compositions of newly identified late Cryogenian A-type and highly fractionated S-type granites from two locations in the northern Tarim Craton. LA-ICP-MS zircon U-Pb analyses yield ages of 642 ± 7 Ma for a syenogranite and 643 ± 4.5 Ma for a mylonitized granite. The syenogranite is weakly peraluminous and shows an A-type granite affinity, as indicated by its high K2O + Na2O contents (8.35–8.64 wt.%), high field strength elements (Zr + Nb + Ce + Y = 435.8 − 463.4 × 10−6), Ga/Al ratios (2.79−2.83), and zircon saturation temperatures (829–844 °C). In contrast, the mylonitized granite contains Al-oversaturated minerals (e.g., garnet) and has high a differentiation index (DI = of 98.9–99.4), with lower zircon saturation temperatures (786–792 °C); the samples display high SiO2 contents (72.99–74.00 wt.%) and A/CNK values (1.16–1.17) and low Nb/Ta and Zr/Hf ratios and are enriched in Rb and depleted in Ba, Sr, which all point to a highly fractionated S-type granite affinity. The granites are characterized by elevated large-ion lithosphere elements (LILEs) and flat high-field-strength elements (HFSEs) patterns, with deep Nb and Ta troughs and pronounced negative Eu anomalies (Eu/Eu* = 0.17–0.38). They show apparently negative εNd(t) values (−10.1 to −9.8 and −6.8 to −7.9, respectively) and εHf(t) values (−9.66 to −1.77 and −33.5 to −1.3, respectively) with Paleoproterozoic crustal model ages, indicating that they were mainly generated by the partial melting of mature crustal materials with a minor contribution from a mantle-derived magmatic source. By integrating with previously published geological, sedimentological, and structural data, we suggest that the granites formed due to a high-temperature gradient in a syn-subduction extensional setting that was probably induced by northward slab rollback of the Paleo-Asian Oceanic lithosphere. Our new data highlight an upper-plate extension in the northern Tarim Craton that constitutes the northern periphery of the Rodinia supercontinent. The linear distribution of late Cryogenian magmatic rocks provides critical evidence for the orogen strike extension of the terminal suture between the Tarim Craton and southwestern Altaids.

The gravitational signature of the dynamics of oceanization in the Gulf of Aden

We perform a new gravity analysis in the Gulf of Aden with the aim of finding new constraints on the geodynamic evolution of the area. Our analysis is developed within the frame of the recent GO_CONS_GCF_2_TIM_R6 global gravity model solution (Brockmann et al., 2021) that reflects the Earth's static gravity field as observed by GOCE (Gravity field and steady-state Ocean Circulation Explorer). We analyzed the solution at different harmonic degrees, to account for different depths and dimensions of the sources. Terrain correction has been performed by means of a spherical tesseroidal methodology (Marotta and Barzaghi, 2017) and the obtained corrected Geodetic Residual Gravitation has been compared to the Model Residual Gravitation predicted by means of a 2D visco-plastic finite element thermo-mechanic model that simulates the dynamics of the Gulf of Aden, from rifting to oceanization. In order to perform the comparison between observed and predicted gravitational features, data have been extracted along five profiles crossing the Gulf of Aden at different sectors, from the south-east to the north west. Via the Model Residual Gravitation we reproduce its geodetic counterpart, obtaining a characteristic hat-shaped pattern, with a central plateau portraying the highest values, flanked by two regions where the residual gravitation diminishes, finally reaching the lowest values at the far continental domains. The maximum variations of the residual gravitation values, from ridge to continental domain, range between 280 and 420 mGal and the steepest decrease occurs over distances of 200 km for the thick crust model and of 100–150 km for the thin crust model. Finally, the comparison between Geodetic and Model Residual Gravitation allows to further constrain the dynamics of the Gulf of Aden oceanization occurring by a slow passive rift of a hot 150 km thick lithosphere characterized by an initial 40 km thick crust.

Crustal and tectonic structure of Northern Mozambique inferred by 2D gravity modeling

The northern region of Mozambique has a complex geological history, with an evolution that spans from the Precambrian Era to the Phanerozoic Era. In this work, we have integrated gravity and geothermal data to delineate the geotectonic evolution of the region, by estimating the thickness of the crust and the lithosphere through which was essential to generate a representative crustal model. It was necessary to complement the knowledge of structural geometry and tectonic evolution of the region. The data used in this study are the Bouguer and geoid anomalies, topography data, and radiogenic heat. These data were pre-processed, topography and geoid anomaly data were filtered by low-pass filter in the frequency and harmonic domains to remove undesirable effects associated with the sources. The data were used to estimate the thickness of the crust and lithosphere, as well as to determine the mean density distribution within the mantle. This was achieved by using a one-dimensional approach, considering the principle of local isostatic compensation, associated with equations governing the distribution of temperature in the crust. The Bouguer anomaly was used to generate a representative crustal 2D model of this region. The results showed that the crust is thinner in Nampula and Cabo Delgado provinces, with thickness ranging from 27 to 31 km, whereas in Niassa varies between 33 and 39 km. The analysis of lithospheric thickness indicates that the provinces of Nampula and Cabo Delgado present a thinning of the lithosphere, with values ranging from 150 to 165 km. Rather than Niassa province which exhibits a thicker lithosphere, ranging from 165 to 195 km. The obtained results underwent a comparative analysis with prior investigations, unveiling a noteworthy concurrence among these findings.

S-wave velocity structure and Vp/Vs ratios beneath Liyuexi Trough, Nansha Block, the South China Sea: Implications for crustal lithology and tectonic nature of the southern continental margin

By processing the vertical component data of OBS2019–2, a wide angle seismic line situated adjacent to the Liyuexi Trough in the Nansha Block of the South China Sea (SCS), we obtained a crustal model of P-wave velocity structure beneath the line and acquired information about the deep crustal structure in the Liyuexi Trough region. However, the lithology remains to be fully constrained. This paper presents the results of processing the horizontal component data of OBS2019–2 to obtain the S-wave velocity structure model and Vp/Vs ratios. The results show that the Vp/Vs ratios of sedimentary layer range from 1.90 to 4.60, while the corresponding S-wave velocities range from 0.4 to 1.8 km/s. As for the crust, the Vp/Vs ratios fall between 1.70 and 1.89, while the S-wave velocities range between 2.9 and 4.2 km/s. The lithological compositions of oceanic crust from top to bottom may consist of basalt, diabase, and gabbro. In contrast, the lithologies of continental crust transition from felsic to mafic downwards, including andesite, granodiorite, granite, greenschist, and anorthosite. The rock compositions of the Continent-Ocean Transition (COT) region are intermediate between the continental and oceanic crusts, exhibiting transitional features. The rock composition of the high-velocity layer (HVL) in the lower crust may consist of amphibolite, and the presence of magma intrusions and underplating in this region could be associated with the Zhongnan-Liyue Fault. Compared to the Dongsha Block, the lithologies of the Liyue Block and Zhongsha Block are more similar, indicating a higher possibility of conjugation between the Liyue Block and the Zhongsha Block.

Deformation of the Juan de Fuca Plate Beneath the Central Cascadia Continental Margin (44°-45°N) in Response to an Upper Plate Load

A 3D crustal model for the central Cascadia continental shelf and Coast Range between 44°N and 45°N shows that the crystalline crust of the forearc wedge beneath the coastline is characterized by a NW-trending, vertical slab of high-velocity rock interpreted to represent the dike complex that fed the Yachats Basalt, which was intruded into the forearc approximately 37 million years ago. A spatial correlation is observed between downward deflection of the crust of the subducting Juan de Fuca plate, inferred from inversion of PmP arrivals to image the Moho surface, and the high velocity (and consequently high density) anomaly underlying the Yachats Basalt. Apparent subsequent rebound of the subducting plate at greater depth suggests a primarily elastic response of the subducting plate to this load. Calculations for a range of plausible values for the magnitude of the load and the width and depth of the depression indicate that the effective elastic thickness of the subducted Juan de Fuca plate is < 6 km. Although our simple analytical models do not include partial support of the load of the slab by the adjacent upper plate crust or time dependence to account for the motion of the slab beneath the load, incorporation of those effects should decrease rather than increase the apparent strength of the subducted plate. We conclude that the subducted Juan de Fuca plate beneath the central Oregon margin is elastically thin and has the potential to store elastic strain energy before rupturing. Our model of a well-defined, focused and static upper plate load that locally deforms the subducted plate within the nominally seismogenic or transitional part of the Cascadia plate boundary may be unique in providing a relatively straightforward scenario for estimating the mechanical properties of the subducted Juan de Fuca plate. We extrapolate from these results to speculate that elastic deformation of the subducting plate may contribute to the low level of seismicity throughout much of the Cascadia forearc in the inter-seismic period between great earthquakes but note that our local results do not preclude faulting or elasto-plastic deformation of a thin and weak plate as it subducts. These results also suggest that the subducting plate should deform in response to larger scale variations in upper plate thickness and density.

Martian Crustal Model From a Joint Inversion of Receiver Functions and Apparent Shear Wave Velocity

AbstractReceiver function (RF) based studies have revealed the layered Martian crustal structure beneath the InSight, which provides key information about the composition, rock alteration, deep magma process, and thermal history of Mars. However, limited by the frequency and trade‐offs between the thickness and velocity, inverting RFs alone is difficult to provide an accurate image of the upper crust. Here, by achieving stable measurements of the apparent shear wave velocities (Vs,app) at different frequencies for high‐quality marsquakes, we perform a joint inversion with the RFs. Our crustal models differ from previous models by having an extra shallow structure in the upper crust, either a layer with a sharp velocity jump interface at 4.0 ± 0.9 km depth or a gradient structure with shear wave velocity increasing gradually from 1.1 ± 0.1 km/s near the surface to 1.7 ± 0.3 km/s at a depth of 6.0 ± 2.4 km, which well fits both RFs and Vs,app's. Such an upper crustal velocity model is consistent with the depth‐decaying of porosity due to self‐compaction and indicates a highly porous upper crust with a surface porosity larger than 30%. Besides, a much higher porosity decay constant than that scaled from the Moon is derived, suggesting different mechanical properties of the lunar and Martian crust.

Crustal and uppermost mantle structure near the Gloria Fault, North Atlantic, from ocean bottom seismometer surface wave observations

SUMMARY In this work, we present both 1-D and 3-D shear wave velocity (Vs) models of the oceanic crust and uppermost mantle below the Deep OCean Test ARray area, located ∼ 70 km north of the central section of the Gloria Fault, in the eastern North-Atlantic Ocean. The velocity models are inferred from the dispersion of surface waves recorded on ocean bottom seismometers. Dispersion measurements are obtained from the analysis of ambient seismic noise at short periods (< 14 s) and teleseismic surface waves at long periods (> 14 s) using the two-station method. The 1-D Vs model is inferred from the joint inversion of Rayleigh wave phase velocities and Love wave group and phase velocities. The 3-D tomographic model is obtained by inversion of 2-D Love wave group velocity maps as a function of depth, further constrained by the average of Love wave phase velocities obtained from ambient noise (4–9 s) and the average Rayleigh and Love wave phase velocities calculated from teleseismic data (14–44 s). The 1-D Vs model shows a sediment layer with a low velocity of 1.05 km s−1, similar to previous studies in the region. Below the sediments, we find an oceanic crust with velocities ranging from 3.3 to 4.5 km s−1. The model reaches an unusually high velocity of 4.9 km s−1 in a 20 km thick layer at depths between 16 and 36 km. We interpret this fast velocity layer as indicative of the presence of harzburgite, a residue of enhanced melting that might have been formed by the proximity between the Mid-Atlantic Ridge and the Azores mantle plume. At greater depths the velocity decreases, forming a low-velocity zone that reaches a minimum at ∼ 70 km depth, which we interpret as the maximum depth for the lithosphere–asthenosphere boundary. The 3-D model shows a structure that is mostly horizontally layered, with Vs isocontours at 3.5–4.5 km s−1 highlighting oscillations of the crustal structure with wavelengths of ∼25–30 km. These oscillations may be due to changes in the rate of mantle upwelling and magma supply rate.

Ground-Motion Model for Significant Duration Constrained by Seismological Simulations

ABSTRACT A duration ground-motion model for crustal earthquakes based on the normalized Arias intensity (IA) is developed. Two sets of seismological simulations are used to constrain the form and scaling of the duration model. Simulations using a 3D crustal model show that an additive model for the source, path, and site terms captures the physical behavior of duration better than a multiplicative model for the site term. Stochastic finite-fault simulations are used to constrain the saturation of the large-magnitude scaling at short distances. The duration model is developed in two parts: a duration model for the time interval between 5% and 75% of the normalized Arias intensity (D5−75) and a duration model for the ratio of the D5−X/D5−75 duration for X values from 10 to 95. Together, these two models provide a more complete description of the evolution of the seismic energy with time than a single duration metric. A new aspect of the statistical model for duration is the inclusion of a random effect for the path term in addition to random effects for the source and site terms. The source and site random effects are modeled as scale factors on the duration, whereas the path-term random effect is a scale factor on the distance slope. The distribution of the duration residuals has a skewness that is between the skewness of a lognormal distribution and the symmetry of a normal distribution. The final duration aleatory variability is modeled by a power-normal distribution with an exponent of 0.3, which accounts for the amplitude dependence of the aleatory variability of the duration with smaller aleatory variability for large-magnitude events and larger aleatory variability for small-magnitude events as compared to the variability from a lognormal distribution.

Zircon U–Pb–Hf isotope and geochemical constraints on the petrogenesis and tectonic setting of Mesoarchean granitoids from the Carajás province, Amazonian craton, Brazil

Based on U–Pb–Hf isotope and whole-rock geochemical data, this paper investigates the origin and tectonic setting of Mesoarchean granitoids from the Ourilândia do Norte area, Carajás province (northern Brazil) and compares them with other Archean units. The following units were selected: (i) tonalite-trondhjemite-granodiorite (TTG) from the Mogno and Rio Verde suites, represented by a 2.92 Ga tonalitic xenolith [εHf(2.92 Ga) = +2.0 to –0.2] and a 2.92 Ga porphyritic trondhjemite [εHf(2.92 Ga) = +2.3 to –3.5], respectively; (ii) the Ourilândia sanukitoid suite, composed of the 2.92 Ga Arraias granodiorite [εHf(2.92 Ga) = +1.9 to –4.4] and 2.88 Ga Ourilândia tonalite-granodiorite complex [εHf(2.88 Ga) = +3.2 to –1.0]; and (iii) the Boa Sorte granite (Canaã dos Carajás granitic suite), which provides four zircon populations dated at 3.04 Ga, 2.98 Ga, 2.93 Ga and 2.88 Ga. The 2.88 Ga population represents the granitic crystallization age, with εHf(2.88 Ga) values ranging from –0.8 to –4.1. The Paleoarchean geological record from the Carajás province is scarce and restricted to detrital zircons (<3.7 Ga), rare xenocrystals (<3.7 Ga) and crustal model ages (Hf-TDMC = 3.6–3.1 Ga), which indicate Paleoarchean crust extraction from the mantle. The zircon U–Pb–Hf isotope and whole-rock geochemistry data from the Mesoarchean granitoids support the hypothesis that the Carajás province experienced mantle enrichment and progressive crustal reworking during 3.1–2.8 Ga, similar to the Mesoarchean record from the Pilbara and Kaapvaal cratons. A proposed three-stage tectonic model explains the petrogenesis and Hf isotope signatures of the studied granitoids. Regardless of the extent of the Paleoarchean geological record, the available data allow us to postulate that the formation of Paleoarchean crust in the Carajás province involved a long-lived dome-and-keel setting (3.6–3.1 Ga), and that this crust was later recycled, allowing mantle enrichment during Mesoarchean low-angle subduction (3.07–2.92 Ga). Then, a short-lived collision defined by peak regional metamorphism (2.89–2.84 Ga) and associated with crustal thickening and slab breakoff gave rise to mantle- and crust-derived magmas at 2.88 Ga.

SKS Polarization Anomalies Due to the Coriolis Force

ABSTRACT The Earth’s Coriolis force has been well-known to impact surface waves and normal modes, which is essential to accurately interpret these waves. However, the Coriolis force on body waves has been assumed to be negligible and mostly ignored. It has been previously shown that the Coriolis force impacts polarizations of shear waves, whereas the wavefronts remain unaffected. We expand on the potential influences of Earth’s Coriolis force on shear-wave polarization measurements by conducting 3D numerical simulations for elastic waves generated by earthquake and explosive sources in a radially symmetric, and 3D mantle and crustal models. The Coriolis force can produce polarization anomalies of mantle shear waves up to 7° and core phases, such as SKS and SKKS, up to 4°. Uncorrected shear-wave polarizations due to the Coriolis force can cause an additional source of error (5°–10° in fast direction, and 0.2–0.3 s delay time depending on the method and seismic phase), inaccurate interpretation of station misalignments, and imprecise estimates of the core–mantle boundary topography. We show how to correct for the Coriolis force on teleseismic shear waves using 1D ray tracing for well-isolated phases. We recommend the use of full waveform simulations to accurately account for earthquake sources parameters, poorly isolated phases that could include interfering phase arrivals within the measurement time window, and the effect of the Coriolis force on the polarizations of shear waves.

Crustal structure across the Central part of Narmada-Son Lineament, India based on the interpretation of aeromagnetic and gravity data: Geological implications

Narmada-Son Lineament (NSL) is a prominent geo-tectonic feature located in central India, characterised by complex geological structures and evolution history. In this study, we present a detailed joint interpretation of the aeromagnetic and gravity anomaly maps to delineate the crustal structure in order to address the thermo-tectonic activity below the NSL region. The NW-SE linear trend observed in the aeromagnetic anomaly map correlates well with the major tectonic elements of the NSL and the Vindhyan basin. The low-pass filtered regional gravity anomaly map and corresponding residual anomaly map revealed several crustal-scale structural features in the study area. The Curie depth map of the region prepared from the magnetic anomaly map shows the variation between 24-48 km. The corresponding heat flow map correlates well with the trends in surface heat flow data. The gravity-derived Moho depth map reveals that the values vary between 30-40 km in the NSL region with deeper Moho below the Vindhyan basin. Further, the crustal models obtained from the constrained potential field modelling for two long transects across the NSL in the present study area revealed i) the absence of crustal underplating in the eastern part of NSL below Hirapur-Mandla region, ii) the presence of mid-crustal horst like structure below NSL in both the transects, iii) undulating Moho with higher density crustal rocks within the NSL compared to the surrounding areas. The modelled crustal structure and the analysis of Curie depth as well as heat flow maps of the study area, reveal that the NSL region was affected by repeated tectonic activities since the Proterozoic times, and the present-day tectonics of this belt is greatly influenced by the reactivation of pre-existing faults along the NSL.

Ambient Noise Tomography Using a Nodal Seismic Array Reveals Evidence for Igneous Intrusion Contributed to Ore Deposits in South China

Abstract The causes of the large-scale and intense magmatism in South China and its coeval metallogenic processes in the Mesozoic have been debated, due partly to the lack of a high-resolution crustal model. In this study, we deploy ∼500 nodal seismometers around the Nanling metallogenic belt in South China, including one linear array across the metallogenic belt and one 2D array around the Fankou ore deposit: a deposit that hosts large amounts of lead and zinc in China. The derivation of the Moho and the intracrustal structures with receiver function imaging using the linear array was conducted in a separate study. Here, we focus mainly on the shallow shear wavespeed structure in the upper crust (from the surface to about 8 km depth) using ambient noise tomography with the 2D array data. We find low-shear wavespeed anomalies (VS&amp;lt;2.8 km/s) associated with a basin down to about 3.0 km and several high-shear wavespeed anomalies (VS&amp;gt;3.6 km/s) below 6.0 km adjacent to the Fankou ore deposit, which we interpret as the crystallized igneous intrusions in the Mesozoic. Gravity modeling shows that the high density anomaly derived using the wavespeed model is required to improve the data fitting. Considering the temperature derived from a typical geothermal gradient is insufficient for mineralization, we propose that the fluid and heat brought by the igneous intrusion might contribute to the mineral concentration of the Fankou ore deposit. This observation is consistent with geochemical and petrological evidence of diabase and quartz diorite veins in the study region. The experiment shows that ambient noise tomography using dense nodal arrays has the potential in resolving the favorable conditions for ore formation and may play a role in future ore exploration.

Regional crustal structure of Indonesia from receiver functions

Characterizing the crustal structure of Indonesia is important to gain a better understanding of its geodynamic evolution and improve seismic hazard assessments in the area. However, a unified crustal model of the entire Indonesian region and its surroundings is lacking. We present new maps of crustal thickness and bulk Vp/Vs ratio in Indonesia and the surrounding area that are obtained using P-wave receiver functions at 36 seismic stations from several permanent regional networks. The measured crustal thickness varies from ∼24 km to ∼38 km. The thickest crust, ∼38 km, is beneath Flores Island, southern Maluku, and neighboring northernmost Australia, whereas the thinnest crust, ∼24 km, is found under eastern Malaysia. Thus, crustal thickness varies by ∼14 km (from ∼24 km to ∼38 km) despite the small changes in elevation at the measurement points. The Vp/Vs ratios are 1.79±0.11, with high values (>1.85) found along the Banda-Sunda arc-trench system. We attribute these high values to: (1) the presence of mafic island arc and oceanic crust and (2) partial melting within this volcanic region, which causes a larger decrease in S-wave velocities compared with P-wave velocities. The comparison of the seismic properties of Indonesian island arc crust, particularly the Vp/Vs ratio, with laboratory measurements and the petrology of the exhumed Talkeetna island arc, Alaska, allows us to infer the crustal composition of Indonesian island arc crust.

A High‐Resolution Shear Velocity Model of the Crust and Uppermost Mantle Beneath Westernmost Mediterranean Including Radial Anisotropy

AbstractUsing seismic data from 1,186 stations deployed across the westernmost Mediterranean, we construct a high‐resolution 3‐D radially anisotropic model from a joint inversion of receiver functions and Rayleigh and Love wave dispersions. The Rayleigh and Love data are extracted from both ambient noise interferograms and earthquake waveforms, and a new three‐station ambient noise interferometry method is used to further improve the data coverage. The obtained crustal thickness map reproduces independently the Moho depth maps compiled from previous data, showing crustal roots beneath the Pyrenean‐Cantabrian range and the Gibraltar Arc and thicker crust beneath the Variscan Iberian Massif than in the area affected by Alpine orogenesis. The Vsv crustal model outlines the Iberian Massif as a high shear wave velocity block and shows extremely low velocities in the Gulf of Cadiz and Gibraltar Arc. At mantle levels, a sharp boundary between the Aquitanean Basin and the Massif Central is imaged, with low Vsv beneath the Massif Central probably reflecting a remaining signature of the magma. To the south of Iberia, the geometry of the Alboran slab is captured by the model, while beneath the Atlas Mountains, widespread low Vsv and positive radial anisotropy is observed, favoring the edge‐driven convection model explaining the lithospheric thinning. The most relevant contribution of this work is mapping, for the first time at this scale, the radial anisotropy anomalies at crustal level, with higher values observed in the Alpine Iberia, dominated by present‐day extensional tectonics.

Moho depth estimation of northern of East African Rift System

Abstract The northern part of the East African Rift System is characterized by depleted Moho depth and thermally thinned lithosphere. This research aims to determine the Moho depth of the study area through non-linear gravity inversion and cross-validation with seismic Moho estimates. The study utilized gravity data to obtain the gravity anomaly of the Moho interface, a topographic grid for removing topographic effects, a crustal model to determine total sediment thickness and its gravitational effect, and seismic Moho depth for constraining the forward model and cross-validation. The estimated Moho depth of the study area ranges between 5 km (in the Indian Ocean) to 45 km (in the Ethiopian Highlands), with slight variation compared to seismic Moho relief. This is because the reference level, calculated for the thinner part of the study region, underestimates the entire area. Upwelling magma in the Eastern branches of the EARS may also incur slight variation in the estimated Moho depth; rifting, volcanism, melt intrusion, magmatic uplift, and tectonic setting all influence the Moho depth of the study area. Furthermore, reverberations affect most seismic Moho estimations in the region. The slight variation can be mitigated by improving the gravity network for accurate validation and precise heat flow measurement to correctly identify magmatic anomalies and density contrasts. Additionally, applying reverberation removal techniques in the study region could improve seismic Moho estimation.