Pub Date : 2025-10-24DOI: 10.1016/j.gsf.2025.102187
Jianping Li , Xing Ding , Huayong Chen
Partition coefficients for Cl between felsic melts and a supercritical aqueous fluid (∼4–16 wt.% NaCleq) were experimentally determined to better constrain Cl behavior during magmatic fluid exsolution in upper-crustal magma chambers. Experiments were conducted at 850 °C, 200 MPa, and oxygen fugacity near NNO + 0.5, using a range of melt and fluid compositions. At constant total chlorinity of 1 mol/kg H2O, values range from 11.3 to 21.1, negatively correlated with both the melt’s aluminum saturation index (ASI) and the HCl/total Cl ratio in the fluid. For a fixed melt composition (ASI = 1.02), values increase linearly from 18.7 to 60.1 as total chlorinity rises from 1 to 4 mol/kg H2O. Rayleigh fractionation modeling of fluid exsolution from upper-crustal magmas using these data indicates that during progressive crystallization, chlorinity of exsolved fluids rapidly decline before stabilizing at ∼1 mol/kg H2O (∼4 wt.% NaCleq), regardless of initial fluid chlorinity or H2O content in melt. This implies that the majority of exsolution fluids released from felsic magmas in the upper crust are of low salinity (∼1 mol/kg H2O). Copper transfer modeling further suggests that efficient metal extraction occurs in Cl- and H2O-rich magmas, particularly where early H2O saturation is achieved, thus favoring the formation of high-grade porphyry copper deposits.
{"title":"Fluid-silicate melt Cl partition and its implications on magmatic fluid exsolution and hydrothermal ore genesis","authors":"Jianping Li , Xing Ding , Huayong Chen","doi":"10.1016/j.gsf.2025.102187","DOIUrl":"10.1016/j.gsf.2025.102187","url":null,"abstract":"<div><div>Partition coefficients for Cl between felsic melts and a supercritical aqueous fluid (∼4–16 wt.% NaCl<sub>eq</sub>) were experimentally determined to better constrain Cl behavior during magmatic fluid exsolution in upper-crustal magma chambers. Experiments were conducted at 850 °C, 200 MPa, and oxygen fugacity near NNO + 0.5, using a range of melt and fluid compositions. At constant total chlorinity of 1 mol/kg H<sub>2</sub>O, <span><math><mrow><msubsup><mtext>D</mtext><mrow><mtext>Cl</mtext></mrow><mrow><mtext>fluid/melt</mtext></mrow></msubsup></mrow></math></span> values range from 11.3 to 21.1, negatively correlated with both the melt’s aluminum saturation index (ASI) and the HCl/total Cl ratio in the fluid. For a fixed melt composition (ASI = 1.02), <span><math><mrow><msubsup><mtext>D</mtext><mrow><mtext>Cl</mtext></mrow><mrow><mtext>fluid/melt</mtext></mrow></msubsup></mrow></math></span> values increase linearly from 18.7 to 60.1 as total chlorinity rises from 1 to 4 mol/kg H<sub>2</sub>O. Rayleigh fractionation modeling of fluid exsolution from upper-crustal magmas using these data indicates that during progressive crystallization, chlorinity of exsolved fluids rapidly decline before stabilizing at ∼1 mol/kg H<sub>2</sub>O (∼4 wt.% NaCl<sub>eq</sub>), regardless of initial fluid chlorinity or H<sub>2</sub>O content in melt. This implies that the majority of exsolution fluids released from felsic magmas in the upper crust are of low salinity (∼1 mol/kg H<sub>2</sub>O). Copper transfer modeling further suggests that efficient metal extraction occurs in Cl- and H<sub>2</sub>O-rich magmas, particularly where early H<sub>2</sub>O saturation is achieved, thus favoring the formation of high-grade porphyry copper deposits.</div></div>","PeriodicalId":12711,"journal":{"name":"Geoscience frontiers","volume":"17 1","pages":"Article 102187"},"PeriodicalIF":8.9,"publicationDate":"2025-10-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145475008","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-23DOI: 10.1016/j.gsf.2025.102172
Aline Costa do Nascimento , Davis Carvalho de Oliveira , Esa Heilimo , Marco Antonio Galarza , Eleilson Oliveira Gabriel , Martin Whitehouse , Matti Kurhila , Cláudio Nery Lamarão
Archean sanukitoids provide crucial insights into crust-mantle interactions during the early Earth’s geodynamic evolution. However, the role of crustal contamination in their genesis remains uncertain. Sanukitoids identified in the Sapucaia subdomain of the southern Carajás Province are represented by two plutons Água Limpa and Água Azul, collectively referred to as the Água Limpa sanukitoid suite. These plutons are compositionally similar to low-Ti sanukitoids (< 0.63 % TiO2) and their zircon isotopic data record a short period of magmatic activity around ca. 2.87 Ga. Sanukitoids zircons reveal εHf(t) values ranging from –3.31 to + 2.14, Hf and Nd model ages between 2.91 Ga and 3.28 Ga, whole-rock εNd(t) values from –1.64 to + 1.76, and δ18O values ranging from 5.0 ‰ to 7.6 ‰. The Pb isotopic compositions in K-feldspar (µ > 10) suggests a Mesoarchean mantle source affected by slight crustal contribution and/or contamination. Result of geochemical modelling indicates that the sanukitoids were formed by ∼ 15 % partial melting of mantle peridotite previously enriched by ∼ 30 % of slab-derived melts, with orthopyroxene, garnet, clinopyroxene, phlogopite, and magnetite as residual phases. The integration of our data with previously published results leads us to suggest that modern-style plate tectonics may have initiated along the northern Carajás Province during the Mesoarchean, while the Rio Maria domain to the south remained dominated by mantle plume-driven crustal growth and vertical tectonics.
{"title":"Crustal contamination and insights into diachronous changes toward modern tectonics: evidence from Hf-Nd-O-Pb isotope of sanukitoids from the Amazonian craton","authors":"Aline Costa do Nascimento , Davis Carvalho de Oliveira , Esa Heilimo , Marco Antonio Galarza , Eleilson Oliveira Gabriel , Martin Whitehouse , Matti Kurhila , Cláudio Nery Lamarão","doi":"10.1016/j.gsf.2025.102172","DOIUrl":"10.1016/j.gsf.2025.102172","url":null,"abstract":"<div><div>Archean sanukitoids provide crucial insights into crust-mantle interactions during the early Earth’s geodynamic evolution. However, the role of crustal contamination in their genesis remains uncertain. Sanukitoids identified in the Sapucaia subdomain of the southern Carajás Province are represented by two plutons Água Limpa and Água Azul, collectively referred to as the Água Limpa sanukitoid suite. These plutons are compositionally similar to low-Ti sanukitoids (< 0.63 % TiO<sub>2</sub>) and their zircon isotopic data record a short period of magmatic activity around ca. 2.87 Ga. Sanukitoids zircons reveal <em>ε</em><sub>Hf</sub>(<em>t</em>) values ranging from –3.31 to + 2.14, Hf and Nd model ages between 2.91 Ga and 3.28 Ga, whole-rock <em>ε</em><sub>Nd</sub>(<em>t</em>) values from –1.64 to + 1.76, and <em>δ</em><sup>18</sup>O values ranging from 5.0 ‰ to 7.6 ‰. The Pb isotopic compositions in K-feldspar (<em>µ</em> > 10) suggests a Mesoarchean mantle source affected by slight crustal contribution and/or contamination. Result of geochemical modelling indicates that the sanukitoids were formed by ∼ 15 % partial melting of mantle peridotite previously enriched by ∼ 30 % of slab-derived melts, with orthopyroxene, garnet, clinopyroxene, phlogopite, and magnetite as residual phases. The integration of our data with previously published results leads us to suggest that modern-style plate tectonics may have initiated along the northern Carajás Province during the Mesoarchean, while the Rio Maria domain to the south remained dominated by mantle plume-driven crustal growth and vertical tectonics.</div></div>","PeriodicalId":12711,"journal":{"name":"Geoscience frontiers","volume":"16 6","pages":"Article 102172"},"PeriodicalIF":8.9,"publicationDate":"2025-10-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145358042","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-23DOI: 10.1016/j.gsf.2025.102185
Qunjia Zhang , Jiacheng Mei , Lujun Lin , Lei Liu , Hanjie Wen , Le Wang
Carbonate clay-type lithium ore holds significant potential due to its extensive reserves, broad distribution, and relatively easy extraction. However, it presents significant beneficiation challenges due to its coexistence with karst-type bauxite, which often results in mixed, low-lithium ores. For the first time, spectral–texture features derived from hyperspectral imaging (HSI) are jointly modeled with a deep learning framework to explore the feasibility of pre-sorting carbonate clay-type lithium ores. Initially, the spectral responses reflecting mineral composition and the texture features characterizing structural differences were analyzed to evaluate the feasibility of using HSI for ore sorting. Furthermore, the influence of band selection, data standardization, and water absorption regions on pre-sorting performance was systematically investigated through comparative analysis of multiple dataset configurations. Two classification schemes, primary ore types classification and multi grade classification, were employed to assess sorting accuracy and identify key influencing factors. The proposed model, A2ST-OSNet, achieves excellent results in both ore localization and classification through staged data input with varying dimensions and a modular design. Results revealed that joint modeling of spectral and texture features enables efficient and accurate pre-sorting, whereas models relying solely on either spectral or texture features were insufficient, as discriminative information for ore sorting is jointly determined by mineral composition and structural characteristics. Moreover, refined feature extraction and fusion strategies, including spectral–texture feature selection and attention mechanisms, proved critical in enhancing classification performance. The proposed approach offers valuable technical support for ore beneficiation and tailings reutilization, contributing to sustainable resource utilization and providing an effective solution for the efficient recycling of carbonate clay-type lithium ores.
{"title":"Sorting carbonate clay-type lithium ores using a deep learning model with adaptive spectral-texture feature fusion","authors":"Qunjia Zhang , Jiacheng Mei , Lujun Lin , Lei Liu , Hanjie Wen , Le Wang","doi":"10.1016/j.gsf.2025.102185","DOIUrl":"10.1016/j.gsf.2025.102185","url":null,"abstract":"<div><div>Carbonate clay-type lithium ore holds significant potential due to its extensive reserves, broad distribution, and relatively easy extraction. However, it presents significant beneficiation challenges due to its coexistence with karst-type bauxite, which often results in mixed, low-lithium ores. For the first time, spectral–texture features derived from hyperspectral imaging (HSI) are jointly modeled with a deep learning framework to explore the feasibility of pre-sorting carbonate clay-type lithium ores. Initially, the spectral responses reflecting mineral composition and the texture features characterizing structural differences were analyzed to evaluate the feasibility of using HSI for ore sorting. Furthermore, the influence of band selection, data standardization, and water absorption regions on pre-sorting performance was systematically investigated through comparative analysis of multiple dataset configurations. Two classification schemes, primary ore types classification and multi grade classification, were employed to assess sorting accuracy and identify key influencing factors. The proposed model, A<sup>2</sup>ST-OSNet, achieves excellent results in both ore localization and classification through staged data input with varying dimensions and a modular design. Results revealed that joint modeling of spectral and texture features enables efficient and accurate pre-sorting, whereas models relying solely on either spectral or texture features were insufficient, as discriminative information for ore sorting is jointly determined by mineral composition and structural characteristics. Moreover, refined feature extraction and fusion strategies, including spectral–texture feature selection and attention mechanisms, proved critical in enhancing classification performance. The proposed approach offers valuable technical support for ore beneficiation and tailings reutilization, contributing to sustainable resource utilization and providing an effective solution for the efficient recycling of carbonate clay-type lithium ores.</div></div>","PeriodicalId":12711,"journal":{"name":"Geoscience frontiers","volume":"17 1","pages":"Article 102185"},"PeriodicalIF":8.9,"publicationDate":"2025-10-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145418362","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-22DOI: 10.1016/j.gsf.2025.102183
Jian Kang , Xie-Yan Song , Wei Xie , Yu-Feng Deng , Bo Wei , Stephen J. Barnes , Wen Chen , Wei-Dong Li
Experimental studies have demonstrated that olivine can be wetted by sulfide liquid under specific conditions. In this study, we investigate elongated sulfide-olivine patches in lherzolite from the Dahuangshan mafic–ultramafic complex (Central Asian Orogenic Belt, NW China), which display textural evidence indicative of this phenomenon. Two-dimensional (2D) microbeam X-ray fluorescence (micro-XRF) mapping and three-dimensional (3D) high-resolution X-ray computed tomography (HRXCT) scanning reveal that the sulfide-olivine patches are 1.0‒4.8 cm long, 0.3‒2.3 cm wide and 0.2‒2.0 cm thick, and have sharp boundaries in the lherzolite matrix. These patches consist entirely of olivine (Fo7–81.9) embedded within interstitial sulfides. The sulfides dominantly comprise pyrrhotite (>94 vol%), minor pentlandite (< 6 vol%) and trace chalcopyrite, and contain a total of 1.1‒1.9 wt.% (Ni + Cu + Co). Low olivine-olivine-sulfide dihedral angles (averaging 44.9°) indicate that olivine was wetted by the sulfide liquid in these patches. The crystal size distributions (CSDs) of the olivine grains in the sulfide-olivine patches are different from those of the olivine in the lherzolite matrix. These observations can be explained by the entrainment of pre-existing sulfide-olivine cumulates into flowing magma. The elongated shape and parallel distribution of the sulfide-olivine patches indicate that the magma flow was laminar. These findings support the hypothesis that sulfide liquid can wet olivine and under the right conditions can be transported and deposited as sulfide-olivine aggregates within magma conduits.
{"title":"Wetting of olivine by sulfide liquid: an example from the Central Asian Orogenic belt in NW China","authors":"Jian Kang , Xie-Yan Song , Wei Xie , Yu-Feng Deng , Bo Wei , Stephen J. Barnes , Wen Chen , Wei-Dong Li","doi":"10.1016/j.gsf.2025.102183","DOIUrl":"10.1016/j.gsf.2025.102183","url":null,"abstract":"<div><div>Experimental studies have demonstrated that olivine can be wetted by sulfide liquid under specific conditions. In this study, we investigate elongated sulfide-olivine patches in lherzolite from the Dahuangshan mafic–ultramafic complex (Central Asian Orogenic Belt, NW China), which display textural evidence indicative of this phenomenon. Two-dimensional (2D) microbeam X-ray fluorescence (micro-XRF) mapping and three-dimensional (3D) high-resolution X-ray computed tomography (HRXCT) scanning reveal that the sulfide-olivine patches are 1.0‒4.8 cm long, 0.3‒2.3 cm wide and 0.2‒2.0 cm thick, and have sharp boundaries in the lherzolite matrix. These patches consist entirely of olivine (Fo<sub>7</sub><sub>–</sub><sub>81.9</sub>) embedded within interstitial sulfides. The sulfides dominantly comprise pyrrhotite (>94 vol%), minor pentlandite (< 6 vol%) and trace chalcopyrite, and contain a total of 1.1‒1.9 wt.% (Ni + Cu + Co). Low olivine-olivine-sulfide dihedral angles (averaging 44.9°) indicate that olivine was wetted by the sulfide liquid in these patches. The crystal size distributions (CSDs) of the olivine grains in the sulfide-olivine patches are different from those of the olivine in the lherzolite matrix. These observations can be explained by the entrainment of pre-existing sulfide-olivine cumulates into flowing magma. The elongated shape and parallel distribution of the sulfide-olivine patches indicate that the magma flow was laminar. These findings support the hypothesis that sulfide liquid can wet olivine and under the right conditions can be transported and deposited as sulfide-olivine aggregates within magma conduits.</div></div>","PeriodicalId":12711,"journal":{"name":"Geoscience frontiers","volume":"17 1","pages":"Article 102183"},"PeriodicalIF":8.9,"publicationDate":"2025-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145418515","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-22DOI: 10.1016/j.gsf.2025.102186
Zishuo Li , Yanbin Yao , Xueguang Zhu
Fluid distribution and migration characteristics hold significant importance in evaluating the quality of reservoir rocks. Two-dimensional nuclear magnetic resonance (2D NMR) measurements have been widely applied to identify specific hydrocarbon contents within rocks. However, due to the complexity of the pore and fluid system, this technique was not fully exploited and was limited by existing data processing methods. In this study, a novel centroid method was developed to enhance the quantification of 2D NMR data for fluids in reservoir rocks. This method calculates the centroid of the 2D NMR map, which correlates with the average pore size derived from the imbibition process. To validate its effectiveness, the method was applied to analyze the T1 (longitudinal relaxation time) – T2 (transverse relaxation time) relationships taken during imbibition processes in three different reservoir rocks. Results demonstrate that the position of centroid can be used to analyze the dominant type of water in pores involved in the imbibition process. Besides, this method can also successfully assess state changes for different water types inside samples by tracking centroids’ movements and fluctuations in centroid T1/T2 ratios, as well as utilizing the 2D NMR map’s signal partitioning capacity. Compared to other approaches, it elucidates the distribution and migration characteristics of water in different types of pores and provides significant advantages in the quantitative processing and comparative analyses of 2D NMR data across various water-bearing conditions. Furthermore, it also demonstrates significant potential for investigating interactions and dynamics of multiphase fluids in unconventional reservoirs.
{"title":"Revealing fluid dynamics in unconventional reservoir rocks: A novel approach of T1-T2 analysis","authors":"Zishuo Li , Yanbin Yao , Xueguang Zhu","doi":"10.1016/j.gsf.2025.102186","DOIUrl":"10.1016/j.gsf.2025.102186","url":null,"abstract":"<div><div>Fluid distribution and migration characteristics hold significant importance in evaluating the quality of reservoir rocks. Two-dimensional nuclear magnetic resonance (2D NMR) measurements have been widely applied to identify specific hydrocarbon contents within rocks. However, due to the complexity of the pore and fluid system, this<!--> <!-->technique was not fully exploited and was limited by existing data processing methods. In this study, a novel centroid method was developed to enhance the quantification of 2D NMR data for fluids in reservoir rocks. This method calculates the centroid of the 2D NMR map, which correlates with the average pore size derived from the imbibition process. To validate its effectiveness, the method was applied to analyze the <em>T</em><sub>1</sub> (longitudinal relaxation time) – <em>T</em><sub>2</sub> (transverse relaxation time) relationships taken during imbibition processes in three different reservoir rocks. Results demonstrate that the position of centroid can be used to analyze the dominant type of water in pores involved in the imbibition process. Besides, this method can also successfully assess state changes for different water types inside samples by tracking centroids’ movements and fluctuations in centroid <em>T</em><sub>1</sub>/<em>T</em><sub>2</sub> ratios, as well as utilizing the 2D NMR map’s signal partitioning capacity. Compared to other approaches, it elucidates the distribution and migration characteristics of water in different types of pores and provides significant advantages in the quantitative processing and comparative analyses of 2D NMR data across various water-bearing conditions. Furthermore, it also demonstrates significant potential for investigating interactions and dynamics of multiphase fluids in unconventional reservoirs.</div></div>","PeriodicalId":12711,"journal":{"name":"Geoscience frontiers","volume":"17 1","pages":"Article 102186"},"PeriodicalIF":8.9,"publicationDate":"2025-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145418514","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Present models for the late Paleoproterozoic evolution of the Aravalli orogen (NW India) postulate the existence of a continental magmatic arc active between 1875 Ma and 1810 Ma, followed by a phase of post-collisional magmatism between 1730 Ma and 1700 Ma. However, the tectono-magmatic processes occurring between these two events remain cryptic. In this study, evidence for an intervening magmatic phase is revealed based on the investigation of granitoids exposed in the southern part of the Aravalli orogen. U-Pb zircon dating of these granitoids (granites to tonalites) yielded emplacement ages between 1770 Ma and 1760 Ma. Whole-rock geochemical data indicate a strongly peraluminous, S-type, high-K calc-alkaline character, with magnesian to ferroan signatures and a syn-collisional tectonic affinity. The REE patterns are predominantly highly fractionated, displaying depleted HREE profiles and moderate to weak negative Eu anomalies. The geochemical data further suggest derivation of the granitoids by partial melting of meta-greywackes at temperatures > 800 °C. Subchondritic εHf(t) values (−11.0 to −2.6) further indicate reworking of a heterogeneous crust. The results of this and previous studies collectively indicate that the Aravalli orogen evolved through three distinct late Paleoproterozoic tectono-magmatic phases: (1) subduction-related magmatism at 1875‒1810 Ma, (2) syn-collisional S-type plutonism at ca. 1770 Ma, and (3) post-collisional extension-related A-type magmatism at ca. 1720 Ma. Globally, Paleoproterozoic S-type granites were predominantly derived by anataxis of Archean crust. Additionally, the data suggest that the northern margin of proto-India collided with fragments of the Columbia supercontinent at ca. 1770 Ma, postdating Columbia’s maximum packing time (1900–1850 Ma).
{"title":"First evidence of ca. 1.77 Ga S-type magmatism, Aravalli orogen: Implications for the late Paleoproterozoic geodynamic evolution of NW India","authors":"Parampreet Kaur , Armin Zeh , Naveen Chaudhri , Prabhakar Dutta , Swati Sharma","doi":"10.1016/j.gsf.2025.102182","DOIUrl":"10.1016/j.gsf.2025.102182","url":null,"abstract":"<div><div>Present models for the late Paleoproterozoic evolution of the Aravalli orogen (NW India) postulate the existence of a continental magmatic arc active between 1875 Ma and 1810 Ma, followed by a phase of post-collisional magmatism between 1730 Ma and 1700 Ma. However, the tectono-magmatic processes occurring between these two events remain cryptic. In this study, evidence for an intervening magmatic phase is revealed based on the investigation of granitoids exposed in the southern part of the Aravalli orogen. U-Pb zircon dating of these granitoids (granites to tonalites) yielded emplacement ages between 1770 Ma and 1760 Ma. Whole-rock geochemical data indicate a strongly peraluminous, S-type, high-K calc-alkaline character, with magnesian to ferroan signatures and a <em>syn</em>-collisional tectonic affinity. The REE patterns are predominantly highly fractionated, displaying depleted HREE profiles and moderate to weak negative Eu anomalies. The geochemical data further suggest derivation of the granitoids by partial melting of <em>meta</em>-greywackes at temperatures > 800 °C. Subchondritic <em>ε</em><sub>Hf</sub>(<em>t</em>) values (−11.0 to −2.6) further indicate reworking of a heterogeneous crust. The results of this and previous studies collectively indicate that the Aravalli orogen evolved through three distinct late Paleoproterozoic tectono-magmatic phases: (1) subduction-related magmatism at 1875‒1810 Ma, (2) <em>syn</em>-collisional S-type plutonism at ca. 1770 Ma, and (3) post-collisional extension-related A-type magmatism at ca. 1720 Ma. Globally, Paleoproterozoic S-type granites were predominantly derived by anataxis of Archean crust. Additionally, the data suggest that the northern margin of proto-India collided with fragments of the Columbia supercontinent at ca. 1770 Ma, postdating Columbia’s maximum packing time (1900–1850 Ma).</div></div>","PeriodicalId":12711,"journal":{"name":"Geoscience frontiers","volume":"17 1","pages":"Article 102182"},"PeriodicalIF":8.9,"publicationDate":"2025-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145361149","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-13DOI: 10.1016/j.gsf.2025.102181
Hu Fu , Mona Alariqi
Critical minerals like copper, lithium, cobalt, nickel, and rare earth elements form the backbone of low-carbon technologies and are central to the success of global energy transitions. Their availability and the security of their supply chains determine the scalability of renewable energy systems, electric vehicles, battery storage, and hydrogen technologies. For member countries of the International Energy Agency (IEA), which plays a pivotal role in global energy markets, ensuring resilient access to these minerals is inseparable from the broader challenge of decoupling economic growth from carbon emissions. This study examines the dynamics of energy and carbon decomposition as mechanisms for decoupling economic growth from energy-related emissions across certain IEA countries between 1995 and 2022. The analysis employs a decomposition framework that incorporates value-added carbon intensity, value-added energy intensity, and CO2 transport and storage while accounting for the enabling role of critical mineral availability. The results reveal that improvements in energy decomposition significantly strengthen the decoupling of growth from emissions, whereas increases in carbon decomposition weaken it. Similarly, higher value-added energy intensity and carbon intensity are positively associated with decoupling, while expanded CO2 transport and storage capacity tend to reduce its effectiveness. Notably, integrating into the analysis considerations related to mineral supply demonstrates that stable and diversified access to critical resources magnifies the benefits of energy decomposition while mitigating the risks that are linked to carbon intensity. These findings underscore the dual importance of policy frameworks that advance energy efficiency and decomposition and strategies that secure critical mineral supply chains to ensure clean technologies’ scalability.
{"title":"Critical mineral supply chains and the economics of energy transition: A carbon decomposition perspective of growth and decoupling","authors":"Hu Fu , Mona Alariqi","doi":"10.1016/j.gsf.2025.102181","DOIUrl":"10.1016/j.gsf.2025.102181","url":null,"abstract":"<div><div>Critical minerals like copper, lithium, cobalt, nickel, and rare earth elements form the backbone of low-carbon technologies and are central to the success of global energy transitions. Their availability and the security of their supply chains determine the scalability of renewable energy systems, electric vehicles, battery storage, and hydrogen technologies. For member countries of the International Energy Agency (IEA), which plays a pivotal role in global energy markets, ensuring resilient access to these minerals is inseparable from the broader challenge of decoupling economic growth from carbon emissions. This study examines the dynamics of energy and carbon decomposition as mechanisms for decoupling economic growth from energy-related emissions across certain IEA countries between 1995 and 2022. The analysis employs a decomposition framework that incorporates value-added carbon intensity, value-added energy intensity, and CO<sub>2</sub> transport and storage while accounting for the enabling role of critical mineral availability. The results reveal that improvements in energy decomposition significantly strengthen the decoupling of growth from emissions, whereas increases in carbon decomposition weaken it. Similarly, higher value-added energy intensity and carbon intensity are positively associated with decoupling, while expanded CO<sub>2</sub> transport and storage capacity tend to reduce its effectiveness. Notably, integrating into the analysis considerations related to mineral supply demonstrates that stable and diversified access to critical resources magnifies the benefits of energy decomposition while mitigating the risks that are linked to carbon intensity. These findings underscore the dual importance of policy frameworks that advance energy efficiency and decomposition and strategies that secure critical mineral supply chains to ensure clean technologies’ scalability.</div></div>","PeriodicalId":12711,"journal":{"name":"Geoscience frontiers","volume":"17 1","pages":"Article 102181"},"PeriodicalIF":8.9,"publicationDate":"2025-10-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145361146","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-12DOI: 10.1016/j.gsf.2025.102180
Yuzhou Ge , Zhendong Luan , Lianfu Li , Shichuan Xi , Zengfeng Du , Xin Zhang
<div><div>Submarine hydrothermal activities release a large amount of greenhouse gases such as CO<sub>2</sub> and CH<sub>4</sub> into the ocean, influencing the global carbon cycle. Carrying out high-temperature and high-pressure hydrothermal experiments to simulate these geochemical processes is a prerequisite for clarifying the source of carbon-containing substances in the hydrothermal fluid. In situ monitoring of the change in carbon isotope composition of CO<sub>2</sub> is essential for high-temperature and high-pressure simulation experiments, but it is also a great technical challenge. In recent years, laser Raman spectroscopy has attracted wide attention as a supplementary means to mass spectrometry for measuring the <sup>13</sup>C/<sup>12</sup>C value of CO<sub>2</sub>. However, the existing research is limited to the Raman spectroscopy study of the carbon isotope composition of supercritical/liquid CO<sub>2</sub>, and there is little research on dissolved CO<sub>2</sub> in solution. In this study, we systematically studied the Raman spectral characteristics of dissolved <sup>13</sup>CO<sub>2</sub> and <sup>12</sup>CO<sub>2</sub> in the H<sub>2</sub>O ± <sup>13</sup>CO<sub>2</sub> ± <sup>12</sup>CO<sub>2</sub> system at 25–300 °C and 10–350 bar. The results show that the peak position of the Raman characteristic band of <sup>13</sup>CO<sub>2</sub> (aq) is 1367–1370 cm<sup>−1</sup>, which is 14–17 cm<sup>−1</sup> lower than that of <sup>12</sup>CO<sub>2</sub> (aq), and the full width at half maximum is 2–3 cm<sup>−1</sup> smaller than that of <sup>12</sup>CO<sub>2</sub>, which indicate the <sup>13</sup>CO<sub>2</sub> (aq) and <sup>12</sup>CO<sub>2</sub> (aq) can be identified by Raman spectroscopy. On this basis, we proposed the optimal likelihood curve fitting method (OLCF) for the first time to decompose the overlapping bands and accurately obtain the peak height ratio (H<sub>13</sub>/H<sub>12</sub>) of dissolved <sup>13</sup>CO<sub>2</sub> and <sup>12</sup>CO<sub>2</sub>. It has been shown that the G-factor ratios (G<sub>13</sub>/G<sub>12</sub>) are significantly affected by temperature and the relative content of <sup>13</sup>CO<sub>2</sub> and <sup>12</sup>CO<sub>2</sub>. Obtaining an appropriate G-factor ratio is a prerequisite for accurately determining the <sup>13</sup>C/<sup>12</sup>C of dissolved CO<sub>2</sub>. Based on the functional relationship between the H<sub>13</sub>/H<sub>12</sub> and the <sup>13</sup>C/<sup>12</sup>C, we established an empirical equation to quickly estimate the <sup>13</sup>C/<sup>12</sup>C value, thereby assisting in selecting an appropriate G-factor ratio. The calibrated G-factor can be well used to determine the <sup>13</sup>C/<sup>12</sup>C of dissolved CO<sub>2</sub> in the hydrothermal experiments with <sup>13</sup>C labelled. In-situ monitoring experiments show that the phase separation of the hydrothermal fluid hardly causes changes in the carbon isotope composition of dissolved CO<sub>2</sub>.
{"title":"In situ Raman spectroscopic measurement of the 13C/12C ratio of dissolved CO2 at high temperatures and pressures: Method and implications","authors":"Yuzhou Ge , Zhendong Luan , Lianfu Li , Shichuan Xi , Zengfeng Du , Xin Zhang","doi":"10.1016/j.gsf.2025.102180","DOIUrl":"10.1016/j.gsf.2025.102180","url":null,"abstract":"<div><div>Submarine hydrothermal activities release a large amount of greenhouse gases such as CO<sub>2</sub> and CH<sub>4</sub> into the ocean, influencing the global carbon cycle. Carrying out high-temperature and high-pressure hydrothermal experiments to simulate these geochemical processes is a prerequisite for clarifying the source of carbon-containing substances in the hydrothermal fluid. In situ monitoring of the change in carbon isotope composition of CO<sub>2</sub> is essential for high-temperature and high-pressure simulation experiments, but it is also a great technical challenge. In recent years, laser Raman spectroscopy has attracted wide attention as a supplementary means to mass spectrometry for measuring the <sup>13</sup>C/<sup>12</sup>C value of CO<sub>2</sub>. However, the existing research is limited to the Raman spectroscopy study of the carbon isotope composition of supercritical/liquid CO<sub>2</sub>, and there is little research on dissolved CO<sub>2</sub> in solution. In this study, we systematically studied the Raman spectral characteristics of dissolved <sup>13</sup>CO<sub>2</sub> and <sup>12</sup>CO<sub>2</sub> in the H<sub>2</sub>O ± <sup>13</sup>CO<sub>2</sub> ± <sup>12</sup>CO<sub>2</sub> system at 25–300 °C and 10–350 bar. The results show that the peak position of the Raman characteristic band of <sup>13</sup>CO<sub>2</sub> (aq) is 1367–1370 cm<sup>−1</sup>, which is 14–17 cm<sup>−1</sup> lower than that of <sup>12</sup>CO<sub>2</sub> (aq), and the full width at half maximum is 2–3 cm<sup>−1</sup> smaller than that of <sup>12</sup>CO<sub>2</sub>, which indicate the <sup>13</sup>CO<sub>2</sub> (aq) and <sup>12</sup>CO<sub>2</sub> (aq) can be identified by Raman spectroscopy. On this basis, we proposed the optimal likelihood curve fitting method (OLCF) for the first time to decompose the overlapping bands and accurately obtain the peak height ratio (H<sub>13</sub>/H<sub>12</sub>) of dissolved <sup>13</sup>CO<sub>2</sub> and <sup>12</sup>CO<sub>2</sub>. It has been shown that the G-factor ratios (G<sub>13</sub>/G<sub>12</sub>) are significantly affected by temperature and the relative content of <sup>13</sup>CO<sub>2</sub> and <sup>12</sup>CO<sub>2</sub>. Obtaining an appropriate G-factor ratio is a prerequisite for accurately determining the <sup>13</sup>C/<sup>12</sup>C of dissolved CO<sub>2</sub>. Based on the functional relationship between the H<sub>13</sub>/H<sub>12</sub> and the <sup>13</sup>C/<sup>12</sup>C, we established an empirical equation to quickly estimate the <sup>13</sup>C/<sup>12</sup>C value, thereby assisting in selecting an appropriate G-factor ratio. The calibrated G-factor can be well used to determine the <sup>13</sup>C/<sup>12</sup>C of dissolved CO<sub>2</sub> in the hydrothermal experiments with <sup>13</sup>C labelled. In-situ monitoring experiments show that the phase separation of the hydrothermal fluid hardly causes changes in the carbon isotope composition of dissolved CO<sub>2</sub>. ","PeriodicalId":12711,"journal":{"name":"Geoscience frontiers","volume":"17 1","pages":"Article 102180"},"PeriodicalIF":8.9,"publicationDate":"2025-10-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145361145","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-08DOI: 10.1016/j.gsf.2025.102179
Andrea Maffeis , Maria Luce Frezzotti , Rosario Esposito , Marco G. Malusà , Alessandro Aiuppa , Andrea Luca Rizzo , Simona Ferrando
This study investigates the petrological and metasomatic processes that lead to carbon enrichment in peridotites from Sal Island, Cape Verde. Geochemical and mineralogical analyses reveal a heterogeneous lithospheric mantle, consisting of harzburgites showing ultrarefractory compositions indicative of 20 %–40 % melting degrees, as well as fertile spinel lherzolites. Evidence of metasomatism is demonstrated by the formation of reaction coronae around dissolving orthopyroxene, consisting of olivine, clinopyroxene, spinel, and interstitial phonolitic glass, together with trachytic/phonolitic glass + carbonate (calcite, aragonite, and dolomite) microveins associated with CO2 fluid-rich melt inclusions (Type I and II) cutting through olivine and orthopyroxene. The widely differing proportions of silicate and carbonate components in inclusions likely reflect heterogeneous trapping of melt/fluid and degassing CO2. Thermobarometric data indicate equilibration temperatures from 950 to 1060 °C in harzburgites and up to 1200 °C for reaction coronas in harzburgites and lherzolites, with pressures reaching the aragonite stability field (∼2.2–3.5 GPa, or 66–106 km depth). These observations indicate the infiltration at the base of the lithosphere of a silicate-carbonate melt enriched in alkalies, Al, and volatiles (Cl, S, F, N, P). In microveins, the silicate glass composition (e.g., K and Ti content) is consistent with experimental partial melts derived from carbonated sediments with a minor addition of a carbonated eclogite. Enrichments in major and trace elements in clinopyroxene in harzburgites and lherzolites suggest at least two significant metasomatic events involving alkali-rich silicate-carbonate melts at the base of the lithosphere, and CO2-rich fluid, alkali-rich silicate melts in the deep lithosphere, close to pressure conditions of the carbonate ledge. The introduction of recycled carbon into the upper mantle beneath the Cape Verde archipelago likely occurred during the multiple subduction events that affected the region in the half a billion years leading to the Pangea assembly. Major mobilisation of crustal components, generation of carbonate-rich melts, and subsequent lithospheric metasomatism were triggered by the Oligocene thermal perturbation associated with the Cape Verde mantle plume.
{"title":"Carbon enrichment processes in the oceanic upper mantle preserved in peridotites from Sal Island (Cape Verde)","authors":"Andrea Maffeis , Maria Luce Frezzotti , Rosario Esposito , Marco G. Malusà , Alessandro Aiuppa , Andrea Luca Rizzo , Simona Ferrando","doi":"10.1016/j.gsf.2025.102179","DOIUrl":"10.1016/j.gsf.2025.102179","url":null,"abstract":"<div><div>This study investigates the petrological and metasomatic processes that lead to carbon enrichment in peridotites from Sal Island, Cape Verde. Geochemical and mineralogical analyses reveal a heterogeneous lithospheric mantle, consisting of harzburgites showing ultrarefractory compositions indicative of 20 %–40 % melting degrees, as well as fertile spinel lherzolites. Evidence of metasomatism is demonstrated by the formation of reaction coronae around dissolving orthopyroxene, consisting of olivine, clinopyroxene, spinel, and interstitial phonolitic glass, together with trachytic/phonolitic glass + carbonate (calcite, aragonite, and dolomite) microveins associated with CO<sub>2</sub> fluid-rich melt inclusions (Type I and II) cutting through olivine and orthopyroxene. The widely differing proportions of silicate and carbonate components in inclusions likely reflect heterogeneous trapping of melt/fluid and degassing CO<sub>2</sub>. Thermobarometric data indicate equilibration temperatures from 950 to 1060 °C in harzburgites and up to 1200 °C for reaction coronas in harzburgites and lherzolites, with pressures reaching the aragonite stability field (∼2.2–3.5 GPa, or 66–106 km depth). These observations indicate the infiltration at the base of the lithosphere of a silicate-carbonate melt enriched in alkalies, Al, and volatiles (Cl, S, F, N, P). In microveins, the silicate glass composition (e.g., K and Ti content) is consistent with experimental partial melts derived from carbonated sediments with a minor addition of a carbonated eclogite. Enrichments in major and trace elements in clinopyroxene in harzburgites and lherzolites suggest at least two significant metasomatic events involving alkali-rich silicate-carbonate melts at the base of the lithosphere, and CO<sub>2</sub>-rich fluid, alkali-rich silicate melts in the deep lithosphere, close to pressure conditions of the carbonate ledge. The introduction of recycled carbon into the upper mantle beneath the Cape Verde archipelago likely occurred during the multiple subduction events that affected the region in the half a billion years leading to the Pangea assembly. Major mobilisation of crustal components, generation of carbonate-rich melts, and subsequent lithospheric metasomatism were triggered by the Oligocene thermal perturbation associated with the Cape Verde mantle plume.</div></div>","PeriodicalId":12711,"journal":{"name":"Geoscience frontiers","volume":"17 1","pages":"Article 102179"},"PeriodicalIF":8.9,"publicationDate":"2025-10-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145323484","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-08DOI: 10.1016/j.gsf.2025.102173
Petar Glišović , Alexander Braun
<div><div>We investigate Earth’s evolution through thermally coupled core-mantle models spanning 4.5 billion years. These models employ a spherical pseudo-spectral approach to solve the conservation equations for mass, momentum, and energy within a compressible, self-gravitating mantle. The methodology incorporates time-dependent treatments for core-mantle coupling, dislocation and diffusion creep mechanisms, internal heating, and thermal conductivity. Using 3-D numerical simulations, we evaluate the sensitivity of mantle cooling, viscosity structure, and inner-core growth to variations in lithospheric viscosity, diffusion viscosity, mechanical surface boundary conditions, and initial core-mantle boundary and core liquidus temperatures. Results underscore the central role of lithospheric viscosity, particularly near an effective value of ∼10<sup>22</sup> <span><math><mrow><mtext>Pa</mtext><mo>·</mo><mtext>s</mtext></mrow></math></span>, in producing a mantle cooling pattern consistent with petrological constraints, characterized by net warming prior to ∼3 billion years ago (Ga) followed by long-term cooling, as predicted by low-Urey-ratio thermal evolution models. Notably, one model with lithospheric viscosity allowed to vary between 10<sup>18</sup> and 10<sup>24</sup> <span><math><mrow><mtext>Pa</mtext><mo>·</mo><mtext>s</mtext></mrow></math></span> exhibits nonlinear rheological feedbacks that trigger an early-stage thermal rebound. This behavior results from a relatively abrupt increase in lithospheric viscosity which redirects the mantle onto a sustained warming trajectory that departs from the expected monotonic cooling. This example also demonstrates how nonlinear parameter interactions can produce non-monotonic thermal evolution. However, lithospheric viscosity alone cannot fully account for present-day observations of heat flux, inner-core radius, and depth-dependent viscosity profiles. We find that varying the activation enthalpy ratio for grain-growth-controlled diffusion viscosity modifies the radial viscosity structure while leaving the overall cooling pattern intact. Furthermore, surface boundary conditions permitting viscous coupling between rigid surface plates and underlying mantle flow — specifically in our plate-like (PL) model — yield the most acceptable mantle cooling rates and dynamic evolution. This PL configuration also facilitates more realistic coupling between surface kinematics and internal convection, allowing plate velocities to emerge from the flow dynamics rather than being imposed. The PL model exhibits patterns which are similar to independently estimated present-day mantle viscosity profiles, including features such as the lithosphere-asthenosphere gradient and the viscosity jump at the 660 km discontinuity. The PL model also exhibits persistent large-scale lateral temperature anomalies, consistent with previous billion-year convection studies, and illustrates how plate-like surface coupling promotes the emergenc
{"title":"Modelling 4.5 billion years of Earth’s thermal evolution: Insights from core-mantle coupling, lithospheric viscosity, grain-size-dependent rheology, and surface boundary conditions","authors":"Petar Glišović , Alexander Braun","doi":"10.1016/j.gsf.2025.102173","DOIUrl":"10.1016/j.gsf.2025.102173","url":null,"abstract":"<div><div>We investigate Earth’s evolution through thermally coupled core-mantle models spanning 4.5 billion years. These models employ a spherical pseudo-spectral approach to solve the conservation equations for mass, momentum, and energy within a compressible, self-gravitating mantle. The methodology incorporates time-dependent treatments for core-mantle coupling, dislocation and diffusion creep mechanisms, internal heating, and thermal conductivity. Using 3-D numerical simulations, we evaluate the sensitivity of mantle cooling, viscosity structure, and inner-core growth to variations in lithospheric viscosity, diffusion viscosity, mechanical surface boundary conditions, and initial core-mantle boundary and core liquidus temperatures. Results underscore the central role of lithospheric viscosity, particularly near an effective value of ∼10<sup>22</sup> <span><math><mrow><mtext>Pa</mtext><mo>·</mo><mtext>s</mtext></mrow></math></span>, in producing a mantle cooling pattern consistent with petrological constraints, characterized by net warming prior to ∼3 billion years ago (Ga) followed by long-term cooling, as predicted by low-Urey-ratio thermal evolution models. Notably, one model with lithospheric viscosity allowed to vary between 10<sup>18</sup> and 10<sup>24</sup> <span><math><mrow><mtext>Pa</mtext><mo>·</mo><mtext>s</mtext></mrow></math></span> exhibits nonlinear rheological feedbacks that trigger an early-stage thermal rebound. This behavior results from a relatively abrupt increase in lithospheric viscosity which redirects the mantle onto a sustained warming trajectory that departs from the expected monotonic cooling. This example also demonstrates how nonlinear parameter interactions can produce non-monotonic thermal evolution. However, lithospheric viscosity alone cannot fully account for present-day observations of heat flux, inner-core radius, and depth-dependent viscosity profiles. We find that varying the activation enthalpy ratio for grain-growth-controlled diffusion viscosity modifies the radial viscosity structure while leaving the overall cooling pattern intact. Furthermore, surface boundary conditions permitting viscous coupling between rigid surface plates and underlying mantle flow — specifically in our plate-like (PL) model — yield the most acceptable mantle cooling rates and dynamic evolution. This PL configuration also facilitates more realistic coupling between surface kinematics and internal convection, allowing plate velocities to emerge from the flow dynamics rather than being imposed. The PL model exhibits patterns which are similar to independently estimated present-day mantle viscosity profiles, including features such as the lithosphere-asthenosphere gradient and the viscosity jump at the 660 km discontinuity. The PL model also exhibits persistent large-scale lateral temperature anomalies, consistent with previous billion-year convection studies, and illustrates how plate-like surface coupling promotes the emergenc","PeriodicalId":12711,"journal":{"name":"Geoscience frontiers","volume":"17 1","pages":"Article 102173"},"PeriodicalIF":8.9,"publicationDate":"2025-10-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145323483","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
扫码关注我们
求助内容：
应助结果提醒方式：
京公网安备 11010802042870号