Journal of Earth System Science最新文献

{"title":"Draft Earthquake Zone Map of India","authors":"S T G Raghukanth, Bhargavi Podili, K P Sreejaya, I D Gupta, A D Roshan, R Sinha, S Chopra, D Srinagesh, Alpa Sheth, R Goswami, H S Mandal, Ram Jivan Singh, J K Chaudhary, S Arun Kumar, C V R Murty","doi":"10.1007/s12040-024-02368-2","DOIUrl":"https://doi.org/10.1007/s12040-024-02368-2","url":null,"abstract":"<p>The paper describes the procedure employed for developing a new earthquake zone map of India as part of the seventh revision of the Indian Earthquake Standard IS 1893 (Part 1). This new zone map is based primarily on a probabilistic earthquake hazard analysis performed at a grid spacing of 0.1°×0.1° in longitudes and latitudes of the entire country. But, for grid locations with small probabilistic hazard estimates, a minimum level of hazard has been estimated deterministically for the most likely maximum magnitude of an earthquake on the nearest mapped fault. Based on the results, the Indian landmass is grouped into five zones, designated as ‘earthquake zones II, III, IV, V, and VI.’ The peak ground accelerations corresponding to a return period of 2475 yr in these zones are estimated as 0.15, 0.30, 0.45, 0.60, and 0.75<i>g</i>, which also include the site amplification effect. Common normalized response spectra are recommended for all five zones, one for each of the three different site soil conditions, as an interim measure.</p>","PeriodicalId":15609,"journal":{"name":"Journal of Earth System Science","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-08-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142195888","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Integrating physiographical and geophysical analyses for the remediation of a water-filled abandoned coal mining site in Chasnala Colliery, Jharkhand, India","authors":"Rajwardhan Kumar, Amit Bera, Saurabh Srivastava, Sanjit Kumar Pal","doi":"10.1007/s12040-024-02377-1","DOIUrl":"https://doi.org/10.1007/s12040-024-02377-1","url":null,"abstract":"<p>This study explores the comprehensive approach of utilizing physiographical and geophysical electrical resistivity tomography (ERT) investigation to evaluate and address the challenges associated with abandoned, unplanned water-filled galleries in Chasnala Colliery, located in Jharkhand, India. The integrated methodologies facilitate a thorough examination of subsurface conditions, encompassing factors such as geological stability, hydrological fluctuations, and environmental considerations. Utilizing physiographical analysis is of utmost importance in identifying locations with potential risks and developing appropriate site-specific reclamation procedures in the study area. The ERT analysis has successfully confirmed the findings of the physiographical study, revealing the presence of five distinct underground galleries, namely, GL1, GL2, GL3, GL4, and GL5, that are likely submerged in water. These galleries establish connections between the underground spaces and the groundwater, as indicated by their low resistivity values of ~50 Ωm or less. The resistivity measurements exhibit variations that can be attributed to fluctuations in the underground water content. The Wenner, Schlumberger, and dipole–dipole arrays have adeptly discerned the existence of water-filled underground galleries with commendable accuracy. However, the joint array configuration stands out as the pre-eminent choice among these standards due to its unparalleled technical robustness. The findings concurred with the notable correlation between water-filled galleries’ spatial arrangement and shallow groundwater level. The integration of physiographical and ERT data improves the precision of subsurface characterization, facilitating informed decision-making for efficient water management and site rehabilitation in the context of opencast mining.</p>","PeriodicalId":15609,"journal":{"name":"Journal of Earth System Science","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-08-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142195890","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Advancing earth science in geotechnical engineering: A data-driven soft computing technique for unconfined compressive strength prediction in soft soil","authors":"Ishwor Thapa, Sufyan Ghani","doi":"10.1007/s12040-024-02374-4","DOIUrl":"https://doi.org/10.1007/s12040-024-02374-4","url":null,"abstract":"<p>This study presents a pioneering approach that combines artificial intelligence and a nature-inspired optimization algorithm to predict soil unconfined compressive strength (UCS). The traditional laboratory-based method of UCS measurement, involving soil sample preparation, is time-consuming, labour-intensive, and prone to low accuracy. In this work, we propose a non-destructive soil UCS measurement technique utilizing robust AI-based models based on ensemble learning and hybrid learning techniques. Support vector machine (SVM) coupled with particle swarm optimization (PSO), extreme gradient boost (XGB), K-nearest neighbour (KNN), and nature-inspired optimization algorithm-based six hybrid ANFIS models, employing input features from experimental data, were adopted for UCS prediction. Model performance was assessed using standard metrics such as root mean square error, mean absolute error, variance account factor (VAF), expanded uncertainty (<i>U</i>₉₅), and coefficient of determination (<i>R</i>²) between predicted and actual unconfined compressive strength. The study employed 274 data points generated in our laboratory. Sensitivity analysis and Pearson correlation techniques were employed to select relevant elements as input features. Fine content, coarse content, liquid limit, plastic limit, plasticity index, and cohesion of soil were identified as the most effective configurations for accurate soil UCS predictions. XGB demonstrated the highest prediction efficiency in the training and testing phase, achieving an impressive <i>R</i>² of 99.2 and 96.8%, respectively. The results also emphasize the importance of the selected features. The experimental validation accuracy of 97% for the developed XGB model, whose data were not used during model calibration and verification, confirmed the generalization capability of the models. This study provides valuable insights for policymakers and industry stakeholders, facilitating optimized soil unconfined strength management practices.</p>","PeriodicalId":15609,"journal":{"name":"Journal of Earth System Science","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-08-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142195889","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Tectonomagmatic evolution of Khairagarh Group in Sitagota syncline, Dongargarh Supergroup, Bastar Craton: Insight into Paleoproterozoic crust mantle processes","authors":"Sunil Kumar Khare, Deepanker Asthana, A S Venkatesh","doi":"10.1007/s12040-024-02366-4","DOIUrl":"https://doi.org/10.1007/s12040-024-02366-4","url":null,"abstract":"<p>This contribution presents for the first time a digital elevation map and 1:50,000 scale geological map of Sitagota syncline, Khairagarh Group, which is spread in around 1000 km² area in the north Bastar Craton (survey of India toposheets 64 C/11 and C/15). We report for the first time, exposures of Algoma-type banded iron formation, intertrappean shale, and oxide and sulphide mineralization in Mangikhuta basalt. Mafic enclaves are reported in the Dongargarh granite. Geochemistry and petrogenetic study of Mangikhuta and Kotima volcanics of Khairagarh Group is presented. Although field investigation and digital elevation map reveal Khairagarh volcano-sedimentary sequence underwent more than one phase of orogeny, the ubiquitous presence of very low-grade metamorphic mineral assemblages in volcanic rocks indicates they did not undergo high P–T transformation and most of the alteration and metamorphism took place at near-surface conditions. Our tectonomagmatic model proposes the occurrence of a rift basin in the north Bastar Craton from 2.46 to 2.2 Ga, resulting in sedimentation and high-Mg basalt to basaltic-andesite magmatism. The genesis of Sitagota syncline is attributed to closure and deformation of this rift basin due to compressive forces, probably related to Paleoproterozoic Dongargarh Kotri mobile belt and Mesoproterozoic central Indian tectonic zone. Tectonomagmatic and geochronological similarity of Khairagarh Group to Lower Wyloo Group of Ashburton basin in Pilbara Craton and Hekpoort and Ongeluk basalt formations of Transvaal basin in Kaapvaal Craton indicates Bastar Craton was part of Vaalbara supercontinent in Paleoproterozoic times.</p>","PeriodicalId":15609,"journal":{"name":"Journal of Earth System Science","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-08-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142195982","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Regional landslide susceptibility zonation utilizing bivariate statistical techniques in the northwestern Himalayas, Jammu and Kashmir, India","authors":"Imran Khan, Harish Bahuguna, Ashutosh Kainthola","doi":"10.1007/s12040-024-02367-3","DOIUrl":"https://doi.org/10.1007/s12040-024-02367-3","url":null,"abstract":"","PeriodicalId":15609,"journal":{"name":"Journal of Earth System Science","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2024-08-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141920319","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Probabilistic seismic hazard analysis of the Coimbatore region, Tamil Nadu using a logic-tree approach","authors":"Manoharan Sambath, Sembulichampalayam Sennimalai Chandrasekaran, Sandeep Maithani, Ganapathy Pattukandan Ganapathy","doi":"10.1007/s12040-024-02356-6","DOIUrl":"https://doi.org/10.1007/s12040-024-02356-6","url":null,"abstract":"<p>The Coimbatore corporation area is comprised of very densely occupied residential and commercial buildings which are prone to future earthquakes. Probabilistic Seismic Hazard Analysis (PSHA) was carried out for the study region using the Classical Cornell approach and the logic-tree approach. A combination of 45 linear/fault sources and an areal source with a 500 km radius has been considered for the study. An updated earthquake catalogue has been compiled from various works of literature and authorized organizations. The collected earthquake catalogue of various magnitude scales has been homogenized into a uniform moment magnitude scale <span>(left({M}_{w}right))</span>. Fore-shocks and after-shocks have been removed from independent events using one of the declustering algorithms. The seismicity parameters have been evaluated using the Guttenberg–Richter recurrence law. A hybrid GMPE composed of three attenuation relationships was used to obtain the ground motion parameters for the study region. The contour maps of Peak Ground Acceleration (PGA) and Peak Spectral Acceleration (PSA) for the bed-rock condition have been presented in terms of 10 and 2% Probability of Exceedance (PoE) for the return period of 475 and 2475 yr, respectively. The Uniform Hazard Response Spectra (UHRS) for Coimbatore city has been compared with (IS 1893-I-(2016) Criteria for earthquake resistant design of structures. Part 1: General provisions and buildings; Bureau of Indian Standards). As a result of deaggregation, the predominant hazard has been found within a 100 km distance and no hazards have been observed from a long distance as a controlling scenario from the analysis.</p>","PeriodicalId":15609,"journal":{"name":"Journal of Earth System Science","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-08-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141938738","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Time series modelling of drought in a semi-arid region of south Gujarat, India","authors":"K A Jariwala, P G Agnihotri","doi":"10.1007/s12040-024-02357-5","DOIUrl":"https://doi.org/10.1007/s12040-024-02357-5","url":null,"abstract":"<h3 data-test=\"abstract-sub-heading\">Abstract</h3><p>Gujarat is situated within a climatic spectrum ranging from arid to semi-arid, characterised by periodic and recurrent drought phenomena coupled with an enduring and consistent challenge of water scarcity. Drought events in the Gujarat region place substantial stress on its water resource infrastructure, thereby impacting not only the hydrological aspects but also exerting notable repercussions on the intricate interplay between agriculture, economics, and the broader societal domain. Despite the implementation of numerous legislations and protocols aimed at mitigating the impact of drought and ameliorating the aftermath, the state of Gujarat experiences a recurring pattern of drought, resurfacing approximately every three years and inflicting substantial disruption upon the lives of its residents. This research focuses on analysing the current drought scenarios using the standard precipitation index (SPI) and meteorological drought modelling using autoregressive integrated moving average with exogenous variable (ARIMAX) for identifying future drought conditions, which is a statistical modelling approach. Analysis of reference locations was carried out to identify the best model and the same model used for drought modelling and forecasting in the remaining locations. A meteorological drought risk map and drought frequency map were prepared using forecasting results. Model validation was done by computing RMSE and <i>R</i>² of multiple locations in the entire region.</p><h3 data-test=\"abstract-sub-heading\">Research highlights</h3><ul>\u0000<li>\u0000<p>The study focuses on the recurring drought challenges faced by Gujarat due to its arid to semi-arid climate, leading to persistent water scarcity issues</p>\u0000</li>\u0000<li>\u0000<p>The research explores the intricate effects of droughts across hydrology, agriculture, economics, and society, highlighting their far-reaching consequences</p>\u0000</li>\u0000<li>\u0000<p>Despite mitigation efforts, a consistent drought pattern emerges approximately every three years, posing ongoing challenges to water resource management and societal stability</p>\u0000</li>\u0000<li>\u0000<p>The research employs a two-pronged strategy – utilising the standardized precipitation index (SPI) for assessing current drought conditions and leveraging the ARIMAX model for future predictions</p>\u0000</li>\u0000<li>\u0000<p>The ARIMAX model is carefully chosen and fine-tuned using reference locations to ensure accurate predictions across diverse geographical areas</p>\u0000</li>\u0000<li>\u0000<p>A significant outcome of the research is the development of drought risk maps, providing spatial insights into vulnerability levels and enabling targeted mitigation strategies</p>\u0000</li>\u0000<li>\u0000<p>The research rigorously validates the accuracy of its predictions using metrics such as root mean squared error (RMSE) and coefficient of determination (<i>R</i>²), enhancing the credibility of the findings</p>\u0000</li>\u0000<li>\u0000<p>The study reveals nuanced","PeriodicalId":15609,"journal":{"name":"Journal of Earth System Science","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-08-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141938737","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Magnetic anomalies of randomly magnetized finite-strike listric fault sources","authors":"M Deepak, B Ramamma, V Chakravarthi","doi":"10.1007/s12040-024-02358-4","DOIUrl":"https://doi.org/10.1007/s12040-024-02358-4","url":null,"abstract":"<p>A generalized forward modelling equation to calculate the magnetic anomalies of a randomly magnetized finite-strike listric fault source is derived using Poisson's relation. This new equation combines both analytic and numeric approaches to realize forward modelling of the anomalous source in any component. Polynomial functions are adopted to simulate the geometry of the curved fault plane between the displaced hanging wall and the footwall of the fault morphology. The utility of the derived equation is epitomized with a theoretical model of a limited-strike listric fault morphology by computing the anomaly in the vertical, horizontal and total field components. It is demonstrated that the magnitude of the anomalous field (in any component) does not remain the same but changes with the profile offset, albeit the anomalous source remains the same. The effect of structure dimensionality (2D <i>vs.</i> 2^1/2D) on the magnitude of the anomalous field is also discussed.</p>","PeriodicalId":15609,"journal":{"name":"Journal of Earth System Science","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-08-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141938739","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A systematic approach to tuning a neural network model and its application in estimating layer parameters from VES Schlumberger data","authors":"Abhirup Chaudhuri, S Venkateshwara Rao, Ankit Singh, M Pradeep Kumar, Debasis Atta","doi":"10.1007/s12040-024-02349-5","DOIUrl":"https://doi.org/10.1007/s12040-024-02349-5","url":null,"abstract":"<p>Interpretation of vertical electrical sounding (VES) data is inherently difficult due to its ambiguity and non-linear characteristics. Conventional least square-based methods rely on a well-defined <i>apriori</i> model, constrained by ground information (available borehole logs and field observations of exposed lithological sections) for meaningful interpretations. In this work, a back-propagation neural network-based model was developed to estimate layer resistivities and thickness from given apparent resistivities. The model was trained and validated on noise-infused synthetic datasets. Since the effectiveness and generalisation of any model depend on its hyperparameter settings, we investigated effective methods for estimating hyperparameters such as learning rate, momentum, model architecture and learning rate scheduling parameters. It is well known that the optimal values of hyperparameters are not entirely independent of each other. Thus, any change in one hyperparameter changes the optimal range of all other hyperparameters, and thus, tuning any hyperparameter individually is futile. This warrants a joint hyperparameter tuning along with network architecture, which was carried out using a modified version of meta-heuristic black hole algorithm. The modifications include randomly flipping one or more coordinates of the population stars (solutions) whose cost function was above a threshold value decided by a mutation rate parameter. This helped in boosting the exploration capability of the algorithm and prune trajectories with higher cost functions. It is demonstrated that with a finely tuned neural network model, reasonable resistivity model parameters which interpret the ground conditions fairly well could be obtained. The model was tested on resistivity-sounding data with associated borehole lithologs and was found to be giving reasonable results. The same model was used to estimate the overburden thickness consisting of topsoil and deposited silts in Bhadradri–Kothagudem district, India. The layer thickness was consistent with those seen in cutbanks in the area. The methods of optimal hyperparameter set estimation are not exclusive to this model and can be used to train models accomplishing other geophysical tasks.</p>","PeriodicalId":15609,"journal":{"name":"Journal of Earth System Science","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-08-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141938735","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Understanding the mudcrack patterns: How layer thickness, grain size, and secondary layers affect their formation – A laboratory study","authors":"Tushar Todkar, Puspendu Saha, Santanu Misra","doi":"10.1007/s12040-024-02337-9","DOIUrl":"https://doi.org/10.1007/s12040-024-02337-9","url":null,"abstract":"<h3 data-test=\"abstract-sub-heading\">Abstract</h3><p>This paper investigates the initiation and evolution of polygonal crack patterns in desiccating soil layers with varying thicknesses and grain sizes, both with and without a secondary saturated sand layer placed beneath the soil layer. Single-layer experiments involved soil samples within a specific range of grain sizes (0–100, 100–300, 300–500, and 500–850 μm), maintaining soil layer depths at 10, 20, 30, 40, 50, and 60 mm. In double-layer experiments, a saturated sand layer was introduced below the soil layer to check the role of relative layer thickness on crack patterns. Time-lapse photography captured surface crack development during desiccation, allowing for the measurement of geometric parameters like crack width, crack intensity factor (CIF), and intersection angle (CIA). The single-layer models indicate an increase in crack width and CIF with greater layer thicknesses, while CIA decreases with increased layer thickness. Additionally, experiments with finer grain sizes exhibit relatively wider cracks, along with higher CIF and CIA. In double-layer models with varying thicknesses of individual layers, crack growth is found to be independent of the upper soil layer's thickness. Instead, crack propagation is controlled by the lower sand layer, as the supply of water from the lower sand layer to the upper soil layer facilitates prolonged desiccation, resulting in larger CIF values. Using digital image processing and the box-counting method, we calculated the fractal dimensions (<i>D</i>) of the cracks were calculated. <i>D</i> demonstrates positive relationships with grain size in both single- and double-layer experiments, suggesting a self-similar evolution of crack patterns in the models with coarse-grained soil.</p><h3 data-test=\"abstract-sub-heading\">Research highlights</h3>\u0000<ul>\u0000<li>\u0000<p>The role of basal water-rich sand layer in the growth of mud-cracks is investigated.</p>\u0000</li>\u0000<li>\u0000<p>Layer thickness and grain size are two additional variables during the experiments.</p>\u0000</li>\u0000<li>\u0000<p>Single mud-layer experiments were used to compare the double-layer experiments.</p>\u0000</li>\u0000<li>\u0000<p>Basal wet-sand delays desiccation, resulting in a higher crack intensity factor.</p>\u0000</li>\u0000<li>\u0000<p>Finer grains produce denser cracks irrespective of the imposed variables.</p>\u0000</li>\u0000<li>\u0000<p>Low fractal dimensions in course-grained soils imply a self-similar crack pattern.</p>\u0000</li>\u0000</ul>","PeriodicalId":15609,"journal":{"name":"Journal of Earth System Science","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-08-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141869555","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}