Pub Date : 2024-04-10DOI: 10.1016/j.pepi.2024.107194
Juntong Liu , Yinshuang Ai , Ying Chen , Jianshe Lei
The Wudalianchi-Erkeshan-Keluo (WEK) volcanic belt is a significant component of intraplate volcanism in Northeast China and is composed of the Wudalianchi, Erkeshan, and Keluo volcanic clusters. Using joint inversion of receiver functions and ambient noise, we construct a high-resolution 3-D S-wave velocity model of the WEK volcanic belt and its adjacent region, taking advantage of a deployed dense seismic array around this volcanic belt. There is a prominent low-velocity anomaly at 8–15 km depth beneath the Wudalianchi volcanic cluster, suggesting the presence of a crustal magma chamber. Low-velocity anomalies are also observed at 30–35 km depth beneath the Erkeshan volcanic cluster and 30–40 km depth beneath the Keluo volcanic cluster, resulting in discontinuous velocity structures at the Moho discontinuity. We further find a distinct low-velocity anomaly in the uppermost mantle beneath the WEK volcanic belt. Combined with previous geophysical and geochemistry studies, we propose a magma system scenario for the WEK volcanic belt. The upwelling molten material from the asthenosphere accumulated in the uppermost mantle and the magma chamber was formed, which provided the same uppermost mantle magma sources for the WEK volcanic belt.
{"title":"The crust-mantle velocity structure beneath the Wudalianchi-Erkeshan-Keluo volcanic belt by joint inversion of receiver functions and ambient noise","authors":"Juntong Liu , Yinshuang Ai , Ying Chen , Jianshe Lei","doi":"10.1016/j.pepi.2024.107194","DOIUrl":"https://doi.org/10.1016/j.pepi.2024.107194","url":null,"abstract":"<div><p>The Wudalianchi-Erkeshan-Keluo (WEK) volcanic belt is a significant component of intraplate volcanism in Northeast China and is composed of the Wudalianchi, Erkeshan, and Keluo volcanic clusters. Using joint inversion of receiver functions and ambient noise, we construct a high-resolution 3-D S-wave velocity model of the WEK volcanic belt and its adjacent region, taking advantage of a deployed dense seismic array around this volcanic belt. There is a prominent low-velocity anomaly at 8–15 km depth beneath the Wudalianchi volcanic cluster, suggesting the presence of a crustal magma chamber. Low-velocity anomalies are also observed at 30–35 km depth beneath the Erkeshan volcanic cluster and 30–40 km depth beneath the Keluo volcanic cluster, resulting in discontinuous velocity structures at the Moho discontinuity. We further find a distinct low-velocity anomaly in the uppermost mantle beneath the WEK volcanic belt. Combined with previous geophysical and geochemistry studies, we propose a magma system scenario for the WEK volcanic belt. The upwelling molten material from the asthenosphere accumulated in the uppermost mantle and the magma chamber was formed, which provided the same uppermost mantle magma sources for the WEK volcanic belt.</p></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":null,"pages":null},"PeriodicalIF":2.3,"publicationDate":"2024-04-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140618181","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-03-31DOI: 10.1016/j.pepi.2024.107185
Gabriel Negrucci Dragone, Mauricio de Souza Bologna
The origin of large igneous provinces (LIPs) is still an enigma but likely involves magma storage and pathways spread throughout the crust, requiring indirect methods for its study. Here, we present 3-D resistivity models derived from the inversion of broadband (∼0.0001–3000 s) magnetotelluric data with 9–13 km lateral spacing in the central Paraná Magmatic Province, an expressive Early Cretaceous LIP in South America. Our results map in greater detail the previously interpreted LIP magma conduit and support, in contrast with seismological models, significant magmatic underplating to explain the observed conductivity near the LIP central axis. The potential axial lava feeder appears as a pair of crustal conductors (5–15 km; >0.1 S/m) parallel to the region of maximum thickness of both pre-volcanic sedimentary rocks and erupted tholeiitic basalts along an extension of at least 800 km. We propose the high conductivity is due to graphite films of precipitated carbon during the ascension of carbon-bearing fluids released by crystallizing magmas underplated at the base of the crust. The association of high conductivity with underplating is supported by high Vp/Vs ratios close to the conductive lineament, by a lower-crustal zone of high P-wave velocities at the basin axis attributed to mafic intrusions, and by a residual gravity high interpreted as gabbros underplated/intruded in the lower crust. Moreover, the conductive lineament is spatially associated with intracrustal high densities inferred from geoid inversion and upper-crustal high P-wave velocities. Early CO2 release during crystallization of underplated magma before eruption could explain the time gap between the Weissert ocean anoxic event and the volcanism. Our study advances in the controversial topic of magmatic intrusive components in the Paraná LIP with implications for LIP generation and paleo-climate studies.
{"title":"Magmatic underplating, plumbing system, and carbon-enhanced electrical conductivity in the Paraná Magmatic Province","authors":"Gabriel Negrucci Dragone, Mauricio de Souza Bologna","doi":"10.1016/j.pepi.2024.107185","DOIUrl":"https://doi.org/10.1016/j.pepi.2024.107185","url":null,"abstract":"<div><p>The origin of large igneous provinces (LIPs) is still an enigma but likely involves magma storage and pathways spread throughout the crust, requiring indirect methods for its study. Here, we present 3-D resistivity models derived from the inversion of broadband (∼0.0001–3000 s) magnetotelluric data with 9–13 km lateral spacing in the central Paraná Magmatic Province, an expressive Early Cretaceous LIP in South America. Our results map in greater detail the previously interpreted LIP magma conduit and support, in contrast with seismological models, significant magmatic underplating to explain the observed conductivity near the LIP central axis. The potential axial lava feeder appears as a pair of crustal conductors (5–15 km; >0.1 S/m) parallel to the region of maximum thickness of both pre-volcanic sedimentary rocks and erupted tholeiitic basalts along an extension of at least 800 km. We propose the high conductivity is due to graphite films of precipitated carbon during the ascension of carbon-bearing fluids released by crystallizing magmas underplated at the base of the crust. The association of high conductivity with underplating is supported by high Vp/Vs ratios close to the conductive lineament, by a lower-crustal zone of high P-wave velocities at the basin axis attributed to mafic intrusions, and by a residual gravity high interpreted as gabbros underplated/intruded in the lower crust. Moreover, the conductive lineament is spatially associated with intracrustal high densities inferred from geoid inversion and upper-crustal high P-wave velocities. Early CO<sub>2</sub> release during crystallization of underplated magma before eruption could explain the time gap between the Weissert ocean anoxic event and the volcanism. Our study advances in the controversial topic of magmatic intrusive components in the Paraná LIP with implications for LIP generation and paleo-climate studies.</p></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":null,"pages":null},"PeriodicalIF":2.3,"publicationDate":"2024-03-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140350196","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-03-30DOI: 10.1016/j.pepi.2024.107186
Niraj Kumar, Om Prakash, V.M. Tiwari, A.P. Singh
The fold belts, which border the cratons, are the building blocks for understanding the origin and modification of older continents. The Aravalli-Delhi Fold Belt (ADFB) in northwestern India provides evidence of crustal block accretion due to continental collision, whereas Malani Igneous Suite (MIS) modified with widespread magmatism. Imaging of upper mantle structure and lithospheric modifications, if any, of such regions has been intricate owing to reworking by subsequent superimposed tectono-magmatic processes. We applied a 2D modelling approach to model the lithospheric architecture along a 1000-km long WNW-SSE geotransect across northwestern India, which has been deformed in the past. Our modelling technique combines terrestrial gravity anomaly, heat flow data, satellite-based geoid, and topographic datasets using the basic premise of thermal steady-state and local isostasy. The overall 38 to 40 km thick crustal geometry underneath the Marwar Block had the maximum lithological heterogeneity. The region surrounding the MIS is characterised by 8–10 km thick high-density (2.78 g/cm3) sills deposited in the upper crust down to 9 km depth and another 10–15 km thick high-density (3.05 g/cm3) mafic mantle material near the Moho. About 42 km thick crust, including an 8 to 10-km thick high-density (3.05 g/cm3) underplated layer at its bottom, characterises the high-relief ADFB. The Vindhyan region of Bundelkhand craton is defined by a ∼ 1 km thick trap, having Moho extending at a depth of ∼40 km. The lithospheric thickness varies substantially from ∼143–168 km underneath the Marwar block, which thins to ∼135 km under the ADFB and thickens gradually to ∼150–165 km beneath the Vindhyan region. Substantial crustal density differences in distinct crustal domains, when integrated with the thin lithosphere, reinforce the concept that tectono-magmatic processes might have modified the lithosphere in NW India.
{"title":"Large variations in the lithospheric thickness of northwestern India: Imprints of collisional and thermal reworking","authors":"Niraj Kumar, Om Prakash, V.M. Tiwari, A.P. Singh","doi":"10.1016/j.pepi.2024.107186","DOIUrl":"https://doi.org/10.1016/j.pepi.2024.107186","url":null,"abstract":"<div><p>The fold belts, which border the cratons, are the building blocks for understanding the origin and modification of older continents. The Aravalli-Delhi Fold Belt (ADFB) in northwestern India provides evidence of crustal block accretion due to continental collision, whereas Malani Igneous Suite (MIS) modified with widespread magmatism. Imaging of upper mantle structure and lithospheric modifications, if any, of such regions has been intricate owing to reworking by subsequent superimposed tectono-magmatic processes. We applied a 2D modelling approach to model the lithospheric architecture along a 1000-km long WNW-SSE geotransect across northwestern India, which has been deformed in the past. Our modelling technique combines terrestrial gravity anomaly, heat flow data, satellite-based geoid, and topographic datasets using the basic premise of thermal steady-state and local isostasy. The overall 38 to 40 km thick crustal geometry underneath the Marwar Block had the maximum lithological heterogeneity. The region surrounding the MIS is characterised by 8–10 km thick high-density (2.78 g/cm<sup>3</sup>) sills deposited in the upper crust down to 9 km depth and another 10–15 km thick high-density (3.05 g/cm<sup>3</sup>) mafic mantle material near the Moho. About 42 km thick crust, including an 8 to 10-km thick high-density (3.05 g/cm<sup>3</sup>) underplated layer at its bottom, characterises the high-relief ADFB. The Vindhyan region of Bundelkhand craton is defined by a ∼ 1 km thick trap, having Moho extending at a depth of ∼40 km. The lithospheric thickness varies substantially from ∼143–168 km underneath the Marwar block, which thins to ∼135 km under the ADFB and thickens gradually to ∼150–165 km beneath the Vindhyan region. Substantial crustal density differences in distinct crustal domains, when integrated with the thin lithosphere, reinforce the concept that tectono-magmatic processes might have modified the lithosphere in NW India.</p></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":null,"pages":null},"PeriodicalIF":2.3,"publicationDate":"2024-03-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140350195","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-03-27DOI: 10.1016/j.pepi.2024.107184
Qiang Fu , Maodu Yan , Mark J. Dekkers , Chong Guan , Liang Yu , Wanlong Xu , Zunbo Xu , Miaomiao Shen , Bingshuai Li
Remagnetization is a common yet notorious phenomenon that interferes with paleogeographic reconstruction. Classical paleomagnetic field tests are helpful in detecting remagnetization but their diagnostic power is limited: remagnetization may occur before folding, the tilting age may be ambiguous, or protracted remagnetization may yield dual polarities. Rock magnetic information can provide other constraints on our understanding of the origin of natural remanent magnetization (NRM). Here we focus on the rock magnetic properties of acknowledged remagnetized limestones and unremagnetized rocks of the Zaduo area in the Eastern Qiangtang Terrane, Tibetan Plateau (China). Chemical remanent magnetization is suggested as a more frequent mechanism than the thermoviscous resetting of the NRM. The secondary NRM resides in authigenic magnetite of stable single domain and superparamagnetic (SP) size which grew during post-depositional burial processes. Both high-field and low-field thermomagnetic runs reveal the alteration of existing iron sulfides to magnetite in the remagnetized limestones. NRM decay curves show that the maximum unblocking temperature of the remagnetized samples is significantly lower than that of the unremagnetized samples. Component analysis of acquisition curves of the isothermal remanent magnetization (IRM) reveals a hard component that represents SP magnetite in remagnetized limestones. This component is absent in unremagnetized rocks. End-member modelling reveals a convex curve in the coefficient of determination versus the number of end-members plot for the unremagnetized limestones, whereas the remagnetized rocks exhibit both near-linear and convex shapes. In addition, quantitative analysis of the hysteresis loop shape for different lithologies indicates its validity in detecting remagnetization. Furthermore, we show the differences in the hysteresis data distributions of the two rock types on the Day plot, the Néel diagram, the Borradaile diagram, and the Fabian diagram. Our research emphasizes that rock magnetic properties can serve as tools to diagnose remagnetization in magnetite-dominated rocks. We recommend a comprehensive rock magnetic study to discriminate remagnetization, involving the Day plot, Fabian diagram, thermal demagnetization curves, IRM component analysis and end member modelling.
再磁化是一种常见但却臭名昭著的现象,它干扰了古地理重建。经典的古磁场测试有助于检测再磁化现象,但其诊断能力有限:再磁化可能发生在褶皱之前,倾斜年代可能模糊不清,或者持久的再磁化可能产生双重极性。岩石磁性信息可以为我们了解天然剩磁(NRM)的起源提供其他限制。在此,我们重点研究了中国青藏高原东羌塘地层扎多地区公认的再磁化灰岩和未磁化岩石的岩石磁性。化学剩磁被认为是比热粘重置非磁化岩更常见的机制。二次非磁性磁化位于稳定单畴和超顺磁性(SP)尺寸的自生磁铁矿中,这种磁铁矿是在沉积后的埋藏过程中生长的。高磁场和低磁场热磁运行均显示,在再磁化的灰岩中,现有的硫化铁转变为磁铁矿。近红外衰减曲线显示,再磁化样品的最大解块温度明显低于未磁化样品。对等温剩磁(IRM)采集曲线的成分分析表明,在再磁化灰岩中存在代表 SP 磁铁矿的硬成分。在未磁化的岩石中则没有这种成分。端元建模显示,未磁化灰岩的确定系数与端元数量的关系曲线呈凸状,而再磁化岩石则呈现近线性和凸状两种形状。此外,对不同岩性的磁滞环形状进行的定量分析表明,磁滞环在检测再磁化方面是有效的。此外,我们还在 Day 图、Néel 图、Borradaile 图和 Fabian 图上显示了两类岩石磁滞数据分布的差异。我们的研究强调,岩石磁性可以作为诊断磁铁矿为主的岩石再磁化的工具。我们建议进行全面的岩石磁性研究来判别再磁化,包括戴伊图、费边图、热退磁曲线、IRM成分分析和末端构件建模。
{"title":"Remagnetization of magnetite-bearing rocks in the Eastern Qiangtang Terrane, Tibetan Plateau (China): Mechanism and diagnosis","authors":"Qiang Fu , Maodu Yan , Mark J. Dekkers , Chong Guan , Liang Yu , Wanlong Xu , Zunbo Xu , Miaomiao Shen , Bingshuai Li","doi":"10.1016/j.pepi.2024.107184","DOIUrl":"10.1016/j.pepi.2024.107184","url":null,"abstract":"<div><p>Remagnetization is a common yet notorious phenomenon that interferes with paleogeographic reconstruction. Classical paleomagnetic field tests are helpful in detecting remagnetization but their diagnostic power is limited: remagnetization may occur before folding, the tilting age may be ambiguous, or protracted remagnetization may yield dual polarities. Rock magnetic information can provide other constraints on our understanding of the origin of natural remanent magnetization (NRM). Here we focus on the rock magnetic properties of acknowledged remagnetized limestones and unremagnetized rocks of the Zaduo area in the Eastern Qiangtang Terrane, Tibetan Plateau (China). Chemical remanent magnetization is suggested as a more frequent mechanism than the thermoviscous resetting of the NRM. The secondary NRM resides in authigenic magnetite of stable single domain and superparamagnetic (SP) size which grew during post-depositional burial processes. Both high-field and low-field thermomagnetic runs reveal the alteration of existing iron sulfides to magnetite in the remagnetized limestones. NRM decay curves show that the maximum unblocking temperature of the remagnetized samples is significantly lower than that of the unremagnetized samples. Component analysis of acquisition curves of the isothermal remanent magnetization (IRM) reveals a hard component that represents SP magnetite in remagnetized limestones. This component is absent in unremagnetized rocks. End-member modelling reveals a convex curve in the coefficient of determination versus the number of end-members plot for the unremagnetized limestones, whereas the remagnetized rocks exhibit both near-linear and convex shapes. In addition, quantitative analysis of the hysteresis loop shape for different lithologies indicates its validity in detecting remagnetization. Furthermore, we show the differences in the hysteresis data distributions of the two rock types on the Day plot, the Néel diagram, the Borradaile diagram, and the Fabian diagram. Our research emphasizes that rock magnetic properties can serve as tools to diagnose remagnetization in magnetite-dominated rocks. We recommend a comprehensive rock magnetic study to discriminate remagnetization, involving the Day plot, Fabian diagram, thermal demagnetization curves, IRM component analysis and end member modelling.</p></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":null,"pages":null},"PeriodicalIF":2.3,"publicationDate":"2024-03-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140403961","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-03-22DOI: 10.1016/j.pepi.2024.107179
Zhu Zhang , Wen-Che Yu , Hsin-Hua Huang
Some seismic evidence suggests that the mantle transition zone (MTZ) may become highly hydrous and anisotropic, particularly in the vicinity of subduction zones. The two-dimensional path-integrated anisotropy from the upper mantle to the MTZ has been well established beneath the northwestern region of South America. However, explicit details of azimuthal anisotropy on the MTZ boundaries remains ambiguous. Therefore, we attempted to constrain the azimuthal anisotropy on the MTZ boundaries by implementing the P-to-S anisotropic receiver function analysis. We detected significant seismic evidence of azimuthal anisotropy on the 410-km discontinuity, but weak anisotropy on the 660-km discontinuity. The synthetic waveform modeling indicated the fast symmetry axis of anisotropy trends 50° from the north and plunges 40° downwards from horizontal with an anisotropy strength of 4.0% near 410 km depth. The direction of anisotropy suggests the mantle material moves downwards and towards the subducting Nazca slab near the depth of 410 km. The increased anisotropy strength around the 410 km suggests the hydrous wadsleyite may attribute to anisotropy in the upper MTZ. The lack of detectable seismic anisotropy near the depth of 660 km could be caused by the insufficient amount of aligned anisotropic minerals, even though the mantle material continues moving downwards.
{"title":"Constraining seismic anisotropy on the mantle transition zone boundaries beneath the subducting Nazca slab","authors":"Zhu Zhang , Wen-Che Yu , Hsin-Hua Huang","doi":"10.1016/j.pepi.2024.107179","DOIUrl":"10.1016/j.pepi.2024.107179","url":null,"abstract":"<div><p>Some seismic evidence suggests that the mantle transition zone (MTZ) may become highly hydrous and anisotropic, particularly in the vicinity of subduction zones. The two-dimensional path-integrated anisotropy from the upper mantle to the MTZ has been well established beneath the northwestern region of South America. However, explicit details of azimuthal anisotropy on the MTZ boundaries remains ambiguous. Therefore, we attempted to constrain the azimuthal anisotropy on the MTZ boundaries by implementing the P-to-S anisotropic receiver function analysis. We detected significant seismic evidence of azimuthal anisotropy on the 410-km discontinuity, but weak anisotropy on the 660-km discontinuity. The synthetic waveform modeling indicated the fast symmetry axis of anisotropy trends 50° from the north and plunges 40° downwards from horizontal with an anisotropy strength of 4.0% near 410 km depth. The direction of anisotropy suggests the mantle material moves downwards and towards the subducting Nazca slab near the depth of 410 km. The increased anisotropy strength around the 410 km suggests the hydrous wadsleyite may attribute to anisotropy in the upper MTZ. The lack of detectable seismic anisotropy near the depth of 660 km could be caused by the insufficient amount of aligned anisotropic minerals, even though the mantle material continues moving downwards.</p></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":null,"pages":null},"PeriodicalIF":2.3,"publicationDate":"2024-03-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140277834","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-03-19DOI: 10.1016/j.pepi.2024.107181
Seiji Tsuboi , Rhett Butler
We perform an adjoint inversion by using antipodal PKPab phases to estimate the heterogeneous Vp structure at the base of the lowermost mantle. We have carefully examined antipodal stations with high S/N ratios during the past 30 years and selected 20 source-receiver pairs with the epicentral distances >179.0 degree and the Mw <7.0. We have used the spectral-element method on the Earth Simulator of JAMSTEC to calculate the synthetic seismograms for a heterogeneous mantle Vp model with the accuracy of period about 8 s (110 km wavelength at the CMB). We have set up the time window to retrieve the PKPab phase from the vertical component of both observed and computed seismograms to calculate the adjoint source to obtain the sensitivity kernels of Vp in the mantle. The computed Vp sensitivity kernel for each event shows the characteristic annulus pattern in the lowermost mantle, which covers a large area of the CMB. The twenty PKPab earthquake-station pairs in this antipodal study contribute the equivalent of about 3140 measurements at the CMB—compared with 1871 previously studied—and provide new data. Therefore, although the number of individual source-receiver pairs is not large, the summed sensitivity kernels of the PKPab phase for Vp structure at the base of the mantle may be sufficient to model heterogeneity of the Vp structure at the CMB. We summed each event kernel to set up a sensitivity kernel of Vp in the lowermost mantle and iterated the inversion to estimate a heterogeneous structure in D″. Although we have iterations that dominantly affect both South America and the South Pacific, the summary final model shows features within South of Africa, South Pacific, SE Australia, and Central & South America. Keeping the Vs model fixed, we map the CMB Vp heterogeneity measured by the parameter Rs,p = dlnVs/dlnVp and find qualitative, proximal correspondence with the degree-2 pattern of Large Low Velocity Shear Provinces observed in shear tomographic models: in the Pacific and Africa Rs,p > 2.0, whereas the surrounding edges of the edges of the Pacific show Rs,p < 2.0.
{"title":"Adjoint inversion of antipodal PKPab waveforms for P wave velocity anomaly at the base of the lower mantle","authors":"Seiji Tsuboi , Rhett Butler","doi":"10.1016/j.pepi.2024.107181","DOIUrl":"10.1016/j.pepi.2024.107181","url":null,"abstract":"<div><p>We perform an adjoint inversion by using antipodal <em>PKPab</em> phases to estimate the heterogeneous <em>Vp</em> structure at the base of the lowermost mantle. We have carefully examined antipodal stations with high S/N ratios during the past 30 years and selected 20 source-receiver pairs with the epicentral distances >179.0 degree and the Mw <7.0. We have used the spectral-element method on the Earth Simulator of JAMSTEC to calculate the synthetic seismograms for a heterogeneous mantle <em>Vp</em> model with the accuracy of period about 8 s (110 km wavelength at the CMB). We have set up the time window to retrieve the <em>PKPab</em> phase from the vertical component of both observed and computed seismograms to calculate the adjoint source to obtain the sensitivity kernels of <em>Vp</em> in the mantle. The computed <em>Vp</em> sensitivity kernel for each event shows the characteristic annulus pattern in the lowermost mantle, which covers a large area of the CMB. The twenty <em>PKPab</em> earthquake-station pairs in this antipodal study contribute the equivalent of about 3140 measurements at the CMB—compared with 1871 previously studied—and provide new data. Therefore, although the number of individual source-receiver pairs is not large, the summed sensitivity kernels of the <em>PKPab</em> phase for <em>Vp</em> structure at the base of the mantle may be sufficient to model heterogeneity of the <em>Vp</em> structure at the CMB. We summed each event kernel to set up a sensitivity kernel of <em>Vp</em> in the lowermost mantle and iterated the inversion to estimate a heterogeneous structure in D″. Although we have iterations that dominantly affect both South America and the South Pacific, the summary final model shows features within South of Africa, South Pacific, SE Australia, and Central & South America. Keeping the Vs model fixed, we map the CMB <em>Vp</em> heterogeneity measured by the parameter <em>R</em><sub><em>s,p</em></sub> = <em>dlnVs/dlnV</em>p and find qualitative, proximal correspondence with the degree-2 pattern of Large Low Velocity Shear Provinces observed in shear tomographic models: in the Pacific and Africa <em>R</em><sub><em>s,p</em></sub> > 2.0, whereas the surrounding edges of the edges of the Pacific show <em>R</em><sub><em>s,p</em></sub> < 2.0.</p></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":null,"pages":null},"PeriodicalIF":2.3,"publicationDate":"2024-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0031920124000396/pdfft?md5=6364b7763d0a4d7b6f2c96d7852286ff&pid=1-s2.0-S0031920124000396-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140182069","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-03-19DOI: 10.1016/j.pepi.2024.107182
F. Alejandro Nava , Lenin Ávila-Barrientos
Seismic clusters in background seismicity have been associated with high stress levels and can be an important precursor to large earthquakes, but there is not a unanimous concept of cluster and most cluster identification methods are cumbersome and involve a priori assumptions. We propose a simple definition of seismic cluster and a straightforward method of identification involving a minimum of parameters that can be objectively determined in a data-driven way according to a principle of low random occurrence. As an illustration, definition and method were applied to the identification of cluster activity from October 1979 to March 2010 in northern Baja California, Mexico, between 118°W to 113°W and 30°N to 33°N, a tectonically complex seismic region with several fault systems. Twenty-one clusters were identified, of which 17 located around the places at the northeastern corner of the study area that would be ruptured on April 4, 2010 by the El Mayor-Cucapah Mw 7.1 earthquake, the largest recorded earthquake in Baja California, Mexico, and the four others occurred within 9 km from its epicenter. Clustering also became slightly more frequent as the time of the earthquake approached, so that if the clustering survey had been carried out before the whole northern Baja California area, the clustering might have identified the future epicentral region as a region of interest to be closely monitored (this earthquake featured foreshock activity starting some 15 days before the main event). Although the reliability of clusters as precursors to large earthquakes is still to be studied, it is certainly useful to have a reliable and simple method to identify and characterize them.
{"title":"Data-driven identification of earthquake clusters: Clusters before the 2010 El mayor-Cucapah earthquake MW 7.1, Baja California, Mexico","authors":"F. Alejandro Nava , Lenin Ávila-Barrientos","doi":"10.1016/j.pepi.2024.107182","DOIUrl":"10.1016/j.pepi.2024.107182","url":null,"abstract":"<div><p>Seismic clusters in background seismicity have been associated with high stress levels and can be an important precursor to large earthquakes, but there is not a unanimous concept of cluster and most cluster identification methods are cumbersome and involve a priori assumptions. We propose a simple definition of seismic cluster and a straightforward method of identification involving a minimum of parameters that can be objectively determined in a data-driven way according to a principle of low random occurrence. As an illustration, definition and method were applied to the identification of cluster activity from October 1979 to March 2010 in northern Baja California, Mexico, between 118°W to 113°W and 30°N to 33°N, a tectonically complex seismic region with several fault systems. Twenty-one clusters were identified, of which 17 located around the places at the northeastern corner of the study area that would be ruptured on April 4, 2010 by the El Mayor-Cucapah <em>M</em><sub><em>w</em></sub> 7.1 earthquake, the largest recorded earthquake in Baja California, Mexico, and the four others occurred within 9 km from its epicenter. Clustering also became slightly more frequent as the time of the earthquake approached, so that if the clustering survey had been carried out before the whole northern Baja California area, the clustering might have identified the future epicentral region as a region of interest to be closely monitored (this earthquake featured foreshock activity starting some 15 days before the main event). Although the reliability of clusters as precursors to large earthquakes is still to be studied, it is certainly useful to have a reliable and simple method to identify and characterize them.</p></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":null,"pages":null},"PeriodicalIF":2.3,"publicationDate":"2024-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140182033","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-03-19DOI: 10.1016/j.pepi.2024.107180
Meysam Mahmoodabadi , Farzam Yamini-Fard , Mohammad Tatar , Ahmad Rashidi
The Eastern Iranian Mountain Ranges (EIR) emerged as a consequence of the Late Cretaceous collision between the Afghan and Lut blocks. However, the response of the uppermost mantle to this collision remains enigmatic. Additionally, although petrological evidence suggests that post-collisional delamination is possible, it has not been conclusively identified in prior regional seismic imagery. This observation leads us to further explore this possibility using a dense seismic network. To gain insight into the geodynamic implications for eastern Iran and address knowledge gaps, we extensively investigated the seismic structure of the uppermost mantle beneath the EIR using a dense seismic network of 34 temporary stations, complemented by data from nine additional local permanent stations. By meticulously analyzing 6589 relative arrival time residuals from teleseismic records with favorable signal-to-noise ratios, we applied a non-linear tomography method to map P-wave velocity perturbations in a relative sense. Our tomographic images unveiled distinct instances of rapid high-velocity anomalies beneath low-velocity regions in the shallow mantle, suggesting the potential occurrence of lithospheric dripping, followed by subsequent asthenospheric upwelling. This observation offers a plausible explanation for the observed post-collisional magmatism over the Lut Block. Furthermore, to maintain the approximately 1.5-km positive residual topography across the EIR, beyond the influence of crustal properties, additional support from the hot and buoyant asthenosphere becomes crucial, particularly in the absence of a substantial lithospheric mantle.
{"title":"Post-collisional lithospheric delamination in eastern Iran, revealed by non-linear teleseismic tomography and residual topography","authors":"Meysam Mahmoodabadi , Farzam Yamini-Fard , Mohammad Tatar , Ahmad Rashidi","doi":"10.1016/j.pepi.2024.107180","DOIUrl":"10.1016/j.pepi.2024.107180","url":null,"abstract":"<div><p>The Eastern Iranian Mountain Ranges (EIR) emerged as a consequence of the Late Cretaceous collision between the Afghan and Lut blocks. However, the response of the uppermost mantle to this collision remains enigmatic. Additionally, although petrological evidence suggests that post-collisional delamination is possible, it has not been conclusively identified in prior regional seismic imagery. This observation leads us to further explore this possibility using a dense seismic network. To gain insight into the geodynamic implications for eastern Iran and address knowledge gaps, we extensively investigated the seismic structure of the uppermost mantle beneath the EIR using a dense seismic network of 34 temporary stations, complemented by data from nine additional local permanent stations. By meticulously analyzing 6589 relative arrival time residuals from teleseismic records with favorable signal-to-noise ratios, we applied a non-linear tomography method to map <em>P</em>-wave velocity perturbations in a relative sense. Our tomographic images unveiled distinct instances of rapid high-velocity anomalies beneath low-velocity regions in the shallow mantle, suggesting the potential occurrence of lithospheric dripping, followed by subsequent asthenospheric upwelling. This observation offers a plausible explanation for the observed post-collisional magmatism over the Lut Block. Furthermore, to maintain the approximately 1.5-km positive residual topography across the EIR, beyond the influence of crustal properties, additional support from the hot and buoyant asthenosphere becomes crucial, particularly in the absence of a substantial lithospheric mantle.</p></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":null,"pages":null},"PeriodicalIF":2.3,"publicationDate":"2024-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140181982","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The hyper-oblique indentation of the Indian plate beneath the Burmese sliver gives rise to the Indo-Burma Ranges (IBR) in the eastern part of the Indian subcontinent. This geological formation encompasses one of the enigmatic hotspots and includes the densely populated regions of India and Myanmar. Harmonic Decomposition (HD) of the receiver functions, derived from the Multi-Taper Correlation (MTC) technique, is used to model seismic anisotropy and morphological crustal deformation caused by subduction underneath IBR. We used teleseismic earthquake data from eleven broadband seismic stations installed within the IBR and its foredeep region. The findings indicate that the attitude of the fast symmetry axis or dipping direction of the interface is influenced by the trend of regional geological features and absolute plate motion, with the IBR exhibiting NNE-SSW and N-S directions and the Himalayan region showing NE-SW and E-W directions. Our results reveal that the coupling of the Indian plate with the Burmese and Eurasian plates induces lithospheric fabrics that align perpendicular to the coupling direction, resulting in anisotropy in the brittle upper crust. Directional analysis of the HD model for the interfaces at the middle or lower crust reveals the strike of the fast symmetry axis in the NNE-SSW direction, which suggests the alignment of minerals and partial melt in the direction of the major shear stress. The interface across the Moho reflects four-lobed periodicity, that is, 90o ambiguity in the strike direction of the fast symmetry axis, varying from the E-W to the N-S directions. The ambiguity indicates the possibility of the 2D-induced entrained mantle flow along the subducting Indian plate and the 3D toroidal flow parallel to the trend of the IBR.
{"title":"Deciphering the crustal anisotropy and mantle flow beneath the indo-Burma ranges from the harmonic decomposition of the receiver functions","authors":"Hitank Kasaundhan, Dhiraj Kumar Singh, Mohit Agrawal","doi":"10.1016/j.pepi.2024.107183","DOIUrl":"10.1016/j.pepi.2024.107183","url":null,"abstract":"<div><p>The hyper-oblique indentation of the Indian plate beneath the Burmese sliver gives rise to the Indo-Burma Ranges (IBR) in the eastern part of the Indian subcontinent. This geological formation encompasses one of the enigmatic hotspots and includes the densely populated regions of India and Myanmar. Harmonic Decomposition (HD) of the receiver functions, derived from the Multi-Taper Correlation (MTC) technique, is used to model seismic anisotropy and morphological crustal deformation caused by subduction underneath IBR. We used teleseismic earthquake data from eleven broadband seismic stations installed within the IBR and its foredeep region. The findings indicate that the attitude of the fast symmetry axis or dipping direction of the interface is influenced by the trend of regional geological features and absolute plate motion, with the IBR exhibiting NN<em>E</em>-SSW and N-S directions and the Himalayan region showing NE-SW and E-W directions. Our results reveal that the coupling of the Indian plate with the Burmese and Eurasian plates induces lithospheric fabrics that align perpendicular to the coupling direction, resulting in anisotropy in the brittle upper crust. Directional analysis of the HD model for the interfaces at the middle or lower crust reveals the strike of the fast symmetry axis in the NNE-SSW direction, which suggests the alignment of minerals and partial melt in the direction of the major shear stress. The interface across the Moho reflects four-lobed periodicity, that is, 90<sup>o</sup> ambiguity in the strike direction of the fast symmetry axis, varying from the <em>E</em>-W to the N-S directions. The ambiguity indicates the possibility of the 2D-induced entrained mantle flow along the subducting Indian plate and the 3D toroidal flow parallel to the trend of the IBR.</p></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":null,"pages":null},"PeriodicalIF":2.3,"publicationDate":"2024-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140182147","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-03-18DOI: 10.1016/j.pepi.2024.107178
Florence D.C. Ramirez , Kate Selway , Clinton P. Conrad , Valerie Maupin , Maxim Smirnov
Fennoscandia is continuously uplifting in response to past deglaciation, termed glacial isostatic adjustment or GIA, and its mantle viscosity is well constrained from ice sheet and sea level data. Here, we compare those GIA-constrained viscosities for the Fennoscandian upper mantle with geophysically-constrained viscosities. We construct the upper mantle viscosity structure of Fennoscandia by inferring temperature and water content from seismic and magnetotelluric (MT) data. Using a 1-D MT model for Fennoscandian cratons together with a global seismic model, we infer an upper mantle viscosity (below 250 km) of ∼1021±2 Pa·s, which encompasses the GIA-constrained viscosities of 1020 − 1021 Pa·s. The GIA viscosities are better matched if the Fennoscandian upper mantle is a wet harzburgite or a dry pyrolite, where pyrolite is ∼10 times more viscous than harzburgite. Using the average temperatures and water contents for harzburgitic upper mantle, the GIA viscosities require 1–4 mm grain sizes indicating a diffusion creep regime. In northwestern Fennoscandia, where a high-resolution 2-D resistivity model is available, greater inferred mantle water content implies viscosities that are 10–100 times lower than those for the Fennoscandian Craton. Our work suggests that the combination of seismic and MT observations can improve upper mantle viscosity estimates, especially for regions with laterally-varying viscosity structures or where GIA constraints are not available. Although our method represents an important step forward, viscosity uncertainty can be further reduced by incorporating additional constraints on rock composition, grain size and mantle stress, as well as more accurate geophysical data, into the viscosity calculation.
芬诺斯坎迪亚因过去的冰川消融而不断隆起,这被称为冰川等静力调整(GIA),其地幔粘度可从冰盖和海平面数据中得到很好的约束。在这里,我们将冰川等静力调整约束下的芬诺斯坎迪亚上地幔粘度与地球物理约束下的粘度进行了比较。我们通过地震和磁法(MT)数据推断温度和含水量,构建了芬诺斯坎迪亚的上地幔粘度结构。利用芬诺斯坎迪亚陨石坑的一维MT模型和全球地震模型,我们推断上地幔(250公里以下)的粘度为Pa-s,包含了GIA约束的Pa-s粘度。如果芬诺斯坎迪亚上地幔是湿润的哈兹堡垒岩或干燥的辉绿岩(辉绿岩的粘度是哈兹堡垒岩的约 10 倍),则 GIA 粘度更匹配。利用哈兹堡岩上地幔的平均温度和含水量,GIA 粘度要求 1-4 毫米的晶粒大小,这表明了一种扩散蠕变机制。在有高分辨率二维电阻率模型的芬诺斯坎迪亚西北部,推断的地幔含水量更大,这意味着粘度比芬诺斯坎迪亚克拉通的粘度低 10-100 倍。我们的工作表明,地震观测和地幔测量相结合可以改进上地幔粘度的估算,特别是对于粘度结构横向变化的地区或没有全球地幔影响约束的地区。尽管我们的方法是向前迈出的重要一步,但通过在粘度计算中加入更多有关岩石成分、粒度和地幔应力的约束条件,以及更精确的地球物理数据,可以进一步降低粘度的不确定性。
{"title":"Lateral and radial viscosity structure beneath Fennoscandia inferred from seismic and magnetotelluric observations","authors":"Florence D.C. Ramirez , Kate Selway , Clinton P. Conrad , Valerie Maupin , Maxim Smirnov","doi":"10.1016/j.pepi.2024.107178","DOIUrl":"10.1016/j.pepi.2024.107178","url":null,"abstract":"<div><p>Fennoscandia is continuously uplifting in response to past deglaciation, termed glacial isostatic adjustment or GIA, and its mantle viscosity is well constrained from ice sheet and sea level data. Here, we compare those GIA-constrained viscosities for the Fennoscandian upper mantle with geophysically-constrained viscosities. We construct the upper mantle viscosity structure of Fennoscandia by inferring temperature and water content from seismic and magnetotelluric (MT) data. Using a 1-D MT model for Fennoscandian cratons together with a global seismic model, we infer an upper mantle viscosity (below 250 km) of ∼10<sup>21±2</sup> Pa·s, which encompasses the GIA-constrained viscosities of 10<sup>20</sup> − 10<sup>21</sup> Pa·s. The GIA viscosities are better matched if the Fennoscandian upper mantle is a wet harzburgite or a dry pyrolite, where pyrolite is ∼10 times more viscous than harzburgite. Using the average temperatures and water contents for harzburgitic upper mantle, the GIA viscosities require 1–4 mm grain sizes indicating a diffusion creep regime. In northwestern Fennoscandia, where a high-resolution 2-D resistivity model is available, greater inferred mantle water content implies viscosities that are 10–100 times lower than those for the Fennoscandian Craton. Our work suggests that the combination of seismic and MT observations can improve upper mantle viscosity estimates, especially for regions with laterally-varying viscosity structures or where GIA constraints are not available. Although our method represents an important step forward, viscosity uncertainty can be further reduced by incorporating additional constraints on rock composition, grain size and mantle stress, as well as more accurate geophysical data, into the viscosity calculation.</p></div>","PeriodicalId":54614,"journal":{"name":"Physics of the Earth and Planetary Interiors","volume":null,"pages":null},"PeriodicalIF":2.3,"publicationDate":"2024-03-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0031920124000360/pdfft?md5=409afc5998d5f5cfe7b9883a6bc4c5e7&pid=1-s2.0-S0031920124000360-main.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140181980","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}