C. Lyddon, Nguyen Chien, G. Vasilopoulos, Michael Ridgill, S. Moradian, Agnieska Olbert, Thomas Coulthard, A. Barkwith, Peter Robins
{"title":"Thresholds for estuarine compound flooding using a combined hydrodynamic–statistical modelling approach","authors":"C. Lyddon, Nguyen Chien, G. Vasilopoulos, Michael Ridgill, S. Moradian, Agnieska Olbert, Thomas Coulthard, A. Barkwith, Peter Robins","doi":"10.5194/nhess-24-973-2024","DOIUrl":"https://doi.org/10.5194/nhess-24-973-2024","url":null,"abstract":"Abstract. Estuarine compound flooding can happen when extreme sea level and river discharges occur concurrently, or in close succession, inundating low-lying coastal regions. Such events are hard to predict and amplify the hazard. Recent UK storms, including Storm Desmond (2015) and Ciara (2020), have highlighted the vulnerability of mountainous Atlantic-facing catchments to the impacts of compound flooding including risk to life and short- and long-term socio-economic damages. To improve prediction and early warning of compound flooding, combined sea and river thresholds need to be established. In this study, observational data and numerical modelling were used to reconstruct the historic flood record of an estuary particularly vulnerable to compound flooding (Conwy, North Wales). The record was used to develop a method for identifying combined sea level and river discharge thresholds for flooding using idealised simulations and joint-probability analyses. The results show how flooding extent responds to increasing total water level and river discharge, with notable amplification in flood extent due to the compounding drivers in some circumstances, and sensitivity (∼ 7 %) due to a 3 h time lag between the drivers. The influence of storm surge magnitude (as a component of total water level) on the flooding extent was only important for scenarios with minor flooding. There was variability as to when and where compound flooding occurred; it was most likely under moderate sea and river conditions (e.g. 60th–70th and 30th–50th percentiles) and only in the middle-estuary zone. For such cases, joint-probability analysis is important for establishing compound flood risk behaviour. Elsewhere in the estuary, either the sea state (lower estuary) or river flow (upper estuary) dominated the hazard, and single-value probability analysis is sufficient. These methods can be applied to estuaries worldwide to identify site-specific thresholds for flooding to support emergency response and long-term coastal management plans.\u0000","PeriodicalId":18922,"journal":{"name":"Natural Hazards and Earth System Sciences","volume":null,"pages":null},"PeriodicalIF":4.6,"publicationDate":"2024-03-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140222938","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Sophie Barthelemy, B. Bonan, J. Calvet, Gilles Grandjean, D. Moncoulon, Dorothée Kapsambelis, S. Bernardie
{"title":"A new approach for drought index adjustment to clay-shrinkage-induced subsidence over France: advantages of the interactive leaf area index","authors":"Sophie Barthelemy, B. Bonan, J. Calvet, Gilles Grandjean, D. Moncoulon, Dorothée Kapsambelis, S. Bernardie","doi":"10.5194/nhess-24-999-2024","DOIUrl":"https://doi.org/10.5194/nhess-24-999-2024","url":null,"abstract":"Abstract. Clay shrinkage, which consists of a reduction in the volume of clay soils during dry periods, can affect buildings and cause subsidence damage. In France, losses due to subsidence are estimated at more than EUR 16 billion for the period 1989–2021 (CCR, 2021) and are expected to increase under the effect of climate warming. This work aims to improve the current understanding of the conditions triggering subsidence by proposing an innovative drought index. We use a daily soil wetness index (SWI) to develop a new annual drought index that can be related to subsidence damage. The SWI is derived from simulations of soil moisture profiles from the interactions between soil–biosphere–atmosphere (ISBA) land surface model developed by Météo-France. The ability of the drought index to correlate with insurance claim data is assessed by calculating the Kendall rank correlation over 20 municipalities in France. The insurance data, aggregated by year and municipality, are provided by the Caisse Centrale de Réassurance (CCR). A total of 1200 configurations of the drought index are considered. They are generated by combining different calculation methods, ISBA simulation settings, soil model layers, and drought percentile thresholds. The analysis includes a comparison with the independent claim data of six additional municipalities and with a record of official “CatNat” (the French national natural disaster compensation scheme) decrees, useful for the analysis. The best results are obtained for drought magnitudes based on SWI values of the 0.8-to-1.0 m deep soil layer, an ISBA simulation with interactive leaf area index (LAI), and consideration of low drought SWI percentile thresholds. Comparison with claim data shows that drought magnitude is able to identify subsidence events while being spatially consistent. This drought magnitude index provides more insight into subsidence triggers while benefiting from advanced land surface modeling schemes (interactive LAI, multilayer soil). This work paves the way for more reliable damage estimates.\u0000","PeriodicalId":18922,"journal":{"name":"Natural Hazards and Earth System Sciences","volume":null,"pages":null},"PeriodicalIF":4.6,"publicationDate":"2024-03-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140223266","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
C. Scaini, A. Tamaro, Baurzhan Adilkhan, Satbek Sarzhanov, V. Ismailov, R. Umaraliev, Mustafo Safarov, V. Belikov, Japar Karayev, E. Fagà
{"title":"A new regionally consistent exposure database for Central Asia: population and residential buildings","authors":"C. Scaini, A. Tamaro, Baurzhan Adilkhan, Satbek Sarzhanov, V. Ismailov, R. Umaraliev, Mustafo Safarov, V. Belikov, Japar Karayev, E. Fagà","doi":"10.5194/nhess-24-929-2024","DOIUrl":"https://doi.org/10.5194/nhess-24-929-2024","url":null,"abstract":"Abstract. Central Asia is highly exposed to a broad range of hazardous phenomena including earthquakes, floods and landslides, which have cause substantial damage in the past. However, disaster risk reduction strategies are still under development in the area. We provide a regional-scale exposure database for population and residential buildings based on existing information from previous exposure development efforts at the regional and national scale. Such datasets are complemented with country-based data (e.g., building census, national statistics) collected by national representatives in each Central Asian country (Kazakhstan, the Kyrgyz Republic, Tajikistan, Turkmenistan, Uzbekistan). We also develop population and residential-building exposure layers for the year 2080, which support the definition of disaster risk reduction strategies in the region.\u0000","PeriodicalId":18922,"journal":{"name":"Natural Hazards and Earth System Sciences","volume":null,"pages":null},"PeriodicalIF":4.6,"publicationDate":"2024-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140225010","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
John Sykes, Håvard Toft, P. Haegeli, Grant Statham
{"title":"Automated Avalanche Terrain Exposure Scale (ATES) mapping – local validation and optimization in western Canada","authors":"John Sykes, Håvard Toft, P. Haegeli, Grant Statham","doi":"10.5194/nhess-24-947-2024","DOIUrl":"https://doi.org/10.5194/nhess-24-947-2024","url":null,"abstract":"Abstract. The Avalanche Terrain Exposure Scale (ATES) is a system for classifying mountainous terrain based on the degree of exposure to avalanche hazard. The intent of ATES is to improve backcountry recreationist's ability to make informed risk management decisions by simplifying their terrain analysis. Access to ATES has been largely limited to manually generated maps in high-use areas due to the cost and time to generate ATES maps. Automated ATES (AutoATES) is a chain of geospatial models which provides a path towards developing ATES maps on large spatial scales for relatively minimal cost compared to manual maps. This research validates and localizes AutoATES using two ATES benchmark maps which are based on independent ATES maps from three field experts. We compare the performance of AutoATES in two study areas with unique snow climate and terrain characteristics: Connaught Creek in Glacier National Park, British Columbia, Canada, and Bow Summit in Banff National Park, Alberta, Canada. Our results show that AutoATES aligns with the ATES benchmark maps in 74.5 % of the Connaught Creek study area and 84.4 % of the Bow Summit study area. This is comparable to independently developed manual ATES maps which on average align with the ATES benchmark maps in 76.1 % of Connaught Creek and 84.8 % of Bow Summit. We also compare a variety of DEM types (lidar, stereo photogrammetry, Canadian National Topographic Database) and resolutions (5–26 m) in Connaught Creek to investigate how input data type affects AutoATES performance. Overall, we find that DEM resolution and type are not strong indicators of accuracy for AutoATES, with a map accuracy of 74.5 % ± 1 % for all DEMs. This research demonstrates the efficacy of AutoATES compared to expert manual ATES mapping methods and provides a platform for large-scale development of ATES maps to assist backcountry recreationists in making more informed avalanche risk management decisions.\u0000","PeriodicalId":18922,"journal":{"name":"Natural Hazards and Earth System Sciences","volume":null,"pages":null},"PeriodicalIF":4.6,"publicationDate":"2024-03-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140227940","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Alan Demortier, M. Mandement, V. Pourret, O. Caumont
{"title":"Assimilation of surface pressure observations from personal weather stations in AROME-France","authors":"Alan Demortier, M. Mandement, V. Pourret, O. Caumont","doi":"10.5194/nhess-24-907-2024","DOIUrl":"https://doi.org/10.5194/nhess-24-907-2024","url":null,"abstract":"Abstract. Spatially dense surface pressure observations from personal weather stations (PWSs) are able to describe pressure patterns at the surface, such as those associated with convective events, in more detail than with standard weather stations (SWSs) only. In this study, the benefit of assimilating PWS observations with the 3DVar and the 3DEnVar data assimilation schemes of the AROME-France model is evaluated over a 1-month period and during a heavy precipitation event in the South of France. Observations of surface pressure from PWSs are bias-corrected, quality-controlled, and thinned with a spacing equal to the horizontal dimension of an AROME-France grid cell. Over France, almost half of the 55 187 available PWS observations are assimilated, which is 129 times more than the number of assimilated SWS observations. Despite the limited advantages found from their assimilation with the 3DVar assimilation scheme, the 3DEnVar assimilation scheme shows systematic improvement and reduces by 10.3 % the root-mean-square deviation in surface pressure between 1 h model forecasts and SWS observations over France. Significant improvement is observed over the first 9 h of the forecasts in mean sea level pressure. Finally, when PWS observations are assimilated with the 3DEnVar assimilation scheme, a surface pressure anomaly generated by a mesoscale convective system – observed by PWSs and not visible without them – is successfully assimilated. In that case, the forecasts of location and temporal evolution of the mesoscale convective system as well as rainfall are closer to the observations when PWS observations are assimilated.\u0000","PeriodicalId":18922,"journal":{"name":"Natural Hazards and Earth System Sciences","volume":null,"pages":null},"PeriodicalIF":4.6,"publicationDate":"2024-03-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140229045","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Simulation analysis of 3D stability of a landslide with a locking segment: a case study of the Tizicao landslide in Maoxian County, southwest China","authors":"Yuntao Zhou, Xiaoyan Zhao, Guangze Zhang, Bernd Wünnemann, Jiajia Zhang, M. Meng","doi":"10.5194/nhess-24-891-2024","DOIUrl":"https://doi.org/10.5194/nhess-24-891-2024","url":null,"abstract":"Abstract. Rock bridges, also known as locking masses in landslides, affect the three-dimensional (3D) stability and deformation patterns of landslides. However, it is always difficult to simulate rock bridges with continuous grid models in 3D landslides due to their discontinuous deformations. Tizicao landslide, located in Maoxian County, southwest China, is a typical landslide with a super-large rock mass volume of about 1388.2 × 104 m3 and a locking segment. To explore a better rock bridge model used to simulate 3D stability and deformations of the Tizicao landslide, this study introduced three rock bridge models into the FLAC3D program, including the intact rock mass model (IRMM), the Jennings model (JM), and the contact surface model with high strength parameters (CSM-HSP). The CSM-HSP model was eventually used in the FLAC3D program to obtain the 3D deformation characteristics of the landslide. In addition, the two-dimensional (2D) stability of the Tizicao landslide was analyzed using the GeoStudio program. The simulation results indicate that the Tizicao landslide is generally stable under current conditions owing to the existence of the locking segment in its southern front. This inference is consistent with the field deformation and monitoring data. It was found that the general stability and local deformations of the landslide are influenced by the locking segment according to the comparison between the 2D and 3D stability. There was a linear relationship between the locking ratio and the factor of safety (Fos), which applied to the 2D stability analysis of the landslides with a locking segment each, while there existed an approximate quadratic parabola suitable for the 3D stability of the landslides. Finally, this study analyzed the laws of the 3D Fos varying with the locking ratio and strength parameters of the locking masses and the sliding surface. Furthermore, it explored the advantages and disadvantages of the three rock bridge models in the simulation of the 3D stability of landslides with a locking segment.\u0000","PeriodicalId":18922,"journal":{"name":"Natural Hazards and Earth System Sciences","volume":null,"pages":null},"PeriodicalIF":4.6,"publicationDate":"2024-03-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140239719","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Flash flood detection via copula-based intensity–duration–frequency curves: evidence from Jamaica","authors":"Dino Collalti, Nekeisha Spencer, Eric Strobl","doi":"10.5194/nhess-24-873-2024","DOIUrl":"https://doi.org/10.5194/nhess-24-873-2024","url":null,"abstract":"Abstract. Extreme rainfall events frequently cause hazardous floods in many parts of the world. With growing human exposure to floods, studying conditions that trigger floods is imperative. Flash floods, in particular, require well-defined models for the timely warning of the population at risk. Intensity–duration–frequency (IDF) curves are a common way to characterize rainfall and flood events. Here, the copula method is employed to model the dependence between the intensity and duration of rainfall events flexibly and separately from their respective marginal distribution. Information about the localization of 93 flash floods in Jamaica was gathered and linked to remote-sensing rainfall data, and additional data on location-specific yearly maximum rainfall events were constructed. The estimated normal copula has Weibull and generalized extreme value (GEV) marginals for duration and intensity, respectively. Due to the two samples, it is possible to pin down above which line in the intensity duration space a rainfall event likely triggers a flash flood. The parametric IDF curve with an associated return period of 216 years is determined as the optimal threshold for flash flood event classification. This methodology delivers a flexible approach to generating rainfall IDF curves that can directly be used to assess flash flood risk.\u0000","PeriodicalId":18922,"journal":{"name":"Natural Hazards and Earth System Sciences","volume":null,"pages":null},"PeriodicalIF":4.6,"publicationDate":"2024-03-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140244689","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
T. Schmid, R. Portmann, Leonie Villiger, K. Schröer, D. Bresch
{"title":"An open-source radar-based hail damage model for buildings and cars","authors":"T. Schmid, R. Portmann, Leonie Villiger, K. Schröer, D. Bresch","doi":"10.5194/nhess-24-847-2024","DOIUrl":"https://doi.org/10.5194/nhess-24-847-2024","url":null,"abstract":"Abstract. Severe hailstorms result in substantial damage to buildings and vehicles, necessitating the quantification of associated risks. Here, we present a novel open-source hail damage model for buildings and cars based on single-polarization radar data and 250 000 geolocated hail damage reports in Switzerland from 2002 to 2021. To this end, we conduct a detailed evaluation of different radar-based hail intensity measures at 1 km resolution and find that the maximum expected severe hail size (MESHS) outperforms the other measures, despite a considerable false-alarm ratio. Asset-specific hail damage impact functions for buildings and cars are calibrated based on MESHS and incorporated into the open-source risk modelling platform CLIMADA. The model successfully estimates the correct order of magnitude for the number of damaged building in 91 %, their total cost in 77 %, the number of damaged vehicles in 74 %, and their total cost in 60 % of over 100 considered large hail events. We found considerable uncertainties in hail damage estimates, which are largely attributable to limitations of radar-based hail detection. Therefore, we explore the usage of crowdsourced hail reports and find substantially improved spatial representation of severe hail for individual events. By highlighting the potential and limitations of radar-based hail size estimates, particularly MESHS, and the utilization of an open-source risk modelling platform, this study represents a significant step towards addressing the gap in risk quantification associated with severe hail events in Switzerland.\u0000","PeriodicalId":18922,"journal":{"name":"Natural Hazards and Earth System Sciences","volume":null,"pages":null},"PeriodicalIF":4.6,"publicationDate":"2024-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140251764","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Ashok Dahal, H. Tanyaş, C. V. Westen, M. Meijde, P. Mai, Raphael Huser, L. Lombardo
{"title":"Space–time landslide hazard modeling via Ensemble Neural Networks","authors":"Ashok Dahal, H. Tanyaş, C. V. Westen, M. Meijde, P. Mai, Raphael Huser, L. Lombardo","doi":"10.5194/nhess-24-823-2024","DOIUrl":"https://doi.org/10.5194/nhess-24-823-2024","url":null,"abstract":"Abstract. Until now, a full numerical description of the spatio-temporal dynamics of a landslide could be achieved only via physically based models. The part of the geoscientific community in developing data-driven models has instead focused on predicting where landslides may occur via susceptibility models. Moreover, they have estimate when landslides may occur via models that belong to the early-warning system or to the rainfall-threshold classes. In this context, few published research works have explored a joint spatio-temporal model structure. Furthermore, the third element completing the hazard definition, i.e., the landslide size (i.e., areas or volumes), has hardly ever been modeled over space and time. However, technological advancements in data-driven models have reached a level of maturity that allows all three components to be modeled (Location, Frequency, and Size). This work takes this direction and proposes for the first time a solution to the assessment of landslide hazard in a given area by jointly modeling landslide occurrences and their associated areal density per mapping unit, in space and time. To achieve this, we used a spatio-temporal landslide database generated for the Nepalese region affected by the Gorkha earthquake. The model relies on a deep-learning architecture trained using an Ensemble Neural Network, where the landslide occurrences and densities are aggregated over a squared mapping unit of 1 km × 1 km and classified or regressed against a nested 30 m lattice. At the nested level, we have expressed predisposing and triggering factors. As for the temporal units, we have used an approximately 6 month resolution. The results are promising as our model performs satisfactorily both in the susceptibility (AUC = 0.93) and density prediction (Pearson r = 0.93) tasks over the entire spatio-temporal domain. This model takes a significant distance from the common landslide susceptibility modeling literature, proposing an integrated framework for hazard modeling in a data-driven context.\u0000","PeriodicalId":18922,"journal":{"name":"Natural Hazards and Earth System Sciences","volume":null,"pages":null},"PeriodicalIF":4.6,"publicationDate":"2024-03-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140258070","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Study on seismic risk assessment model of water supply systems in mainland China","authors":"Tianyang Yu, Banghua Lu, Hui Jiang, Zhi Liu","doi":"10.5194/nhess-24-803-2024","DOIUrl":"https://doi.org/10.5194/nhess-24-803-2024","url":null,"abstract":"Abstract. Using the peak ground acceleration (PGA) under four exceedance probabilities calculated via a probabilistic seismic hazard analysis method for China, the probability density function of PGA was obtained by fitting with the Cornell seismic hazard exceedance probability PGA function model. Combined with the seismic fragility function of the water supply system calculated on the basis of the empirical matrix of actual earthquake damage and the exposure of fixed assets, the expected seismic loss and expected loss rate models of the water supply system were established, and the classification standard with the expected seismic loss rate of the water supply system as the index was proposed. The seismic fragility of the water supply system was classified, and the exposure of the water supply system was analyzed. The total fixed assets in the Statistical Yearbook of Urban Water Supply, henceforth called Water Supply Yearbook, were taken as the exposure to earthquake in the region. The accuracy of the fragility model in this paper was verified through the actual earthquake damage losses in Deyang. Taking the water supply system of 720 cities in mainland China as an example, the distribution maps of expected seismic loss and expected loss rate were calculated and drawn. The expected loss rate model was verified by the key earthquake prevention areas in mainland China. The assessment model based on the expected loss and expected loss rate was taken as the seismic risk assessment model of water supply systems in mainland China.\u0000","PeriodicalId":18922,"journal":{"name":"Natural Hazards and Earth System Sciences","volume":null,"pages":null},"PeriodicalIF":4.6,"publicationDate":"2024-03-07","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140259601","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}