{"title":"Attenuation of coda waves on frequency and lapse time in Izmir and Western Anatolia using local earthquakes","authors":"Berkay Kalkar, Elcin Gok","doi":"10.1007/s10950-025-10315-x","DOIUrl":"10.1007/s10950-025-10315-x","url":null,"abstract":"<div><p>This study investigated frequency dependent attenuation of seismic waves from earthquakes in Izmir and its surroundings. Earthquakes with a magnitude of M ≥ 2.3, recorded by the local seismic network in the region between 2008–2023, were used. The network of 26 accelerometers is predominantly deployed near the gulf to monitor tectonic features characterized by high seismic activity. The attenuation properties related to seismic wave propagation were evaluated by applying Coda Q analysis. Attenuation in coda waves is due to wave propagation between the inelastic medium and discontinuities or scattering of heterogeneities. The frequency dependence of the seismic quality factor for the coda wave was determined using Q<sub>0</sub> multiplied by f<sup>n</sup>. Various center frequencies (1,2,4,8,16) were selected to calculate the coda wave quality factor, and bandpass filtering was applied. The CODAQ algorithm in the SEISAN software (Havskov et al. 2020) was used in the study. It was observed that elevated Q values in the study area are associated with geologically stable zones, whereas lower Q<sub>0</sub> values reflect areas of heightened seismic activity. The complex tectonic regime of Izmir and its surrounding areas aligns with the low Q<sub>0</sub> values, suggesting significant heterogeneity in the crust. We found the result Q<sub>c</sub> = 61,9. f<sup>0,9</sup> with investigating the Q change examined using different lapse times (20 s,30 s and 40 s) and geometrical spreading (0.5 and 1) parameters.</p></div>","PeriodicalId":16994,"journal":{"name":"Journal of Seismology","volume":"29 4","pages":"963 - 982"},"PeriodicalIF":2.0,"publicationDate":"2025-08-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145236734","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}
Xueliang Duan, Yanqiong Liu, Bian Zhu, Liting Du, Long Zhang, Cheng Peng, Yonghao Li, Zhenghua Zhou
{"title":"Analysis of acceleration dominant azimuth distribution characteristics of Hualien MS7.3 earthquake","authors":"Xueliang Duan, Yanqiong Liu, Bian Zhu, Liting Du, Long Zhang, Cheng Peng, Yonghao Li, Zhenghua Zhou","doi":"10.1007/s10950-025-10305-z","DOIUrl":"10.1007/s10950-025-10305-z","url":null,"abstract":"<div><p>On April 3, 2024, an <i>M</i><sub>S</sub>7.3 earthquake occurred in Hualien County, Taiwan. Some strong ground motion acceleration records were captured by the relevant stations of the Central Meteorological Bureau of Taiwan and the Palert earthquake early warning system. Based on these acceleration records, this paper proposes a method for calculating and analyzing the dominant azimuth of strong ground motion. The distribution characteristics of the dominant azimuth of strong ground motion in each region are quantitatively analyzed from the perspective of epicentral distance and geological division and the characteristics of ground motion in sea area are also analyzed. The research results show that the dominant azimuth of strong ground motion is randomly distributed throughout the entire area affected by this earthquake, with no significant dominant azimuth interval. The dominant azimuth of strong ground motion in the near-field region and the eastern coastal area is clearly distributed, with most of the acceleration dominant azimuth distributed in the range 150°-210°. The dominant azimuth of strong ground motion is uniformly distributed in the western impact plain area of Taiwan without significant dominant azimuth interval, while the dominant azimuth of strong ground motion in the eastern mountainous area and the central hilly area is clearly distributed. In addition, the acceleration time history in the sea area shows the characteristic of low frequency and long period.</p></div>","PeriodicalId":16994,"journal":{"name":"Journal of Seismology","volume":"29 4","pages":"779 - 793"},"PeriodicalIF":2.0,"publicationDate":"2025-07-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145236868","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":"Evaluation of soil-structure resonance by determining local soil conditions in Gümüşhane city center","authors":"Suna Altundaş, Abdulsamed Tanış","doi":"10.1007/s10950-025-10310-2","DOIUrl":"10.1007/s10950-025-10310-2","url":null,"abstract":"<div><p>Recent earthquakes have adversely affected human life, highlighting the necessity of identifying structural and ground conditions to assess risky areas and mitigate potential damage from future seismic events. Gümüşhane, the designated study area, is susceptible to the impacts of a possible earthquake due to its proximity to the North Anatolian Fault Zone (NAFZ). Accordingly, vertical electrical sounding (VES), multi-channel analysis of surface waves (MASW), and microtremor methods were employed in the central district of Gümüşhane to characterize local ground conditions and to evaluate the resonance hazard potential arising from the interaction between structures and the underlying soil. Dominant period (T<sub>0</sub>) and frequency (f<sub>0</sub>), spectral amplification (A<sub>0</sub>), soil index of vulnerability (Kg), average S wave velocity for 30 m depth (Vs<sub>30</sub>), ground type, and soil-structure resonance hazard potential were determined. For the study area, f<sub>0</sub> values between 1.32–8.91 Hz (T<sub>0</sub> values varied between 0.11–0.76 s), A<sub>0</sub> values between 0.93–7.63 according to horizontal/vertical spectral ratio (HVSR) analysis, and Vs<sub>30</sub> values were between 390–886 m/sec to MASW. As a result of the evaluation of microtremor measurements made in all buildings using the floor spectral ratio (FSR) technique, the soil-structure potential resonance hazard was found to be low. The present study, conducted in Gümüşhane, examined the subject of soil-structure interaction. In the process of planning new residential areas, it is imperative to ensure the reliability of construction based on the local ground conditions. Furthermore, it is essential to determine the risk status of buildings in relation to potential earthquake movements.</p></div>","PeriodicalId":16994,"journal":{"name":"Journal of Seismology","volume":"29 4","pages":"875 - 902"},"PeriodicalIF":2.0,"publicationDate":"2025-07-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145236742","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":"Converting PSH estimates in terms of ground motion intensity into macroseismic intensity estimates","authors":"Dario Albarello","doi":"10.1007/s10950-025-10313-z","DOIUrl":"10.1007/s10950-025-10313-z","url":null,"abstract":"<div><p>Methodological aspects relative to the conversion of probabilistic hazard evaluations in terms of ground motion intensity into hazard in terms of macroseismic Intensity. As a first step, the probabilistic relationships between ground motion intensity and macroseismic intensity are critically re-examined to account for the different formal properties of these two dimensions of the earthquake intensity. Then these relationships are coherently considered in the conversion of hazard curves by accounting for the inherent binning of macroseismic intensity values and of their inherent probabilistic nature. It is shown that approximate procedures currently used to provide this conversion may provide biased outcomes when the dispersion of values relative to ground motion intensity corresponding to macroseismic Intensity is dismissed or not properly accounted for. The impact of this possible bias is evaluated by considering the seismic hazard at reference site condition at the city of L’Aquila in Central Italy.</p></div>","PeriodicalId":16994,"journal":{"name":"Journal of Seismology","volume":"29 4","pages":"933 - 941"},"PeriodicalIF":2.0,"publicationDate":"2025-07-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145236952","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":"Focal mechanisms and bouguer-gravity anomalies of the 2025 earthquake cluster in the Santorini-Amorgos region (Southern Aegean Sea, Greece): evidence for shallow extensional magmatism","authors":"Mustafa Toker, Evrim Yavuz, Emir Balkan","doi":"10.1007/s10950-025-10314-y","DOIUrl":"10.1007/s10950-025-10314-y","url":null,"abstract":"<div><p>Understanding clustered earthquake sequences is essential for seismic hazard assessment, as it involves constraining faulting styles and nodal planes of potential ruptures. This study investigates the nature of a dense earthquake sequence (~ 3000 events) initiated on January 27, 2025, in the Santorini-Amorgos region of the Southern Aegean Sea (SAS), a tectonically active Volcanic Island Arc (VIA). We analyzed 23 shallow crustal earthquakes (Mw ≥ 4.5, depth ≤ 10 km) that occurred between February 2–9, 2025, using full-waveform, low-frequency Centroid Moment Tensor (CMT) inversion from regional seismograms. The inversion was complemented by high-resolution Bouguer gravity anomaly data derived from the EIGEN-6C4 satellite gravity model to assess subsurface density variations. The focal mechanisms consistently indicate NE-SW striking, high-angle (≥ 45°) normal faults with NW- and SE-dipping planes and centroid depths ≤ 10 km. Integration of CMT results with gravity anomalies (90–100 mgal) suggests a migrating zone of shallower extensional magmatism (SEM) driving the sequence. These findings reveal a Precursory Seismic Cluster (PSC) and provide new constraints on the seismotectonic and magmatic processes shaping seismic hazard in the region.</p></div>","PeriodicalId":16994,"journal":{"name":"Journal of Seismology","volume":"29 4","pages":"943 - 961"},"PeriodicalIF":2.0,"publicationDate":"2025-07-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145236953","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":"Lightweight deep transfer learning for earthquake detection in resource-constrained IoT devices","authors":"Irshad Khan, Jae-Kwang Anh, Young-Woo Kwon","doi":"10.1007/s10950-025-10303-1","DOIUrl":"10.1007/s10950-025-10303-1","url":null,"abstract":"<div><p>In a natural disaster, intelligent Internet of Things (IoT) systems can be utilized to respond appropriately. Recently, the application of IoT technology in seismology, particularly in earthquake detection, has garnered much attention. This approach’s attractiveness lies in its simplicity of installation, minimal processing power requirements, cost-effectiveness, and expansive coverage, even in areas lacking Internet connectivity. However, the locality of installed sensors brings variations in seismic and noise data, making the earthquake detection task very challenging because of the false alarms. Network-based systems connecting multiple IoTs can resolve the issue by running highly computation-intensive algorithms on a powerful server or cloud and aggregating the data sent from those sensors. On the other hand, Standalone IoT devices operate independently, making decisions locally using both traditional and machine learning methods to manage false alarms. However, these techniques struggle to handle diverse noise patterns and often fail to detect low-magnitude earthquakes in noisy environments. While deep learning models can enhance earthquake detection in such conditions, their high computational cost makes them impractical for resource-constrained devices. To address these challenges, this article introduces a lightweight deep learning model incorporating a transfer learning approach for standalone devices. The proposed model outperforms traditional machine learning methods in earthquake detection using IoT sensors while significantly reducing computational demands. Designed to operate without internet connectivity, the Multi-headed Convolutional Neural Network (MCNN) model achieves 99% accuracy without incurring additional processing costs. Furthermore, it demonstrates high adaptability and the ability to update rapidly with minimal configuration changes.</p></div>","PeriodicalId":16994,"journal":{"name":"Journal of Seismology","volume":"29 4","pages":"759 - 777"},"PeriodicalIF":2.0,"publicationDate":"2025-07-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145236985","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}
Meng Gong, ChangSheng Lu, Yuyan Qi, Xiaoshan Wang, Xiao Tian, Jin Li
{"title":"CNN-ECA based classification of natural earthquakes and quarry blasting","authors":"Meng Gong, ChangSheng Lu, Yuyan Qi, Xiaoshan Wang, Xiao Tian, Jin Li","doi":"10.1007/s10950-025-10306-y","DOIUrl":"10.1007/s10950-025-10306-y","url":null,"abstract":"<div><p>The rapid and accurate identification of natural earthquakes and artificially blasted earthquakes is crucial for effective earthquake monitoring and early warning. We used waveform data from 5480 natural earthquake events and 4482 blasting events recorded by 110 seismic stations in a quarry in Utah, USA from January 2013 to August 2017, to construct a deep machine learning based CNN-ECA model and accurately and efficiently identify and verify these two types of earthquakes. Firstly, these data were preprocessed by removing mean, trend, instrument response removal, resampling (100 Hz), and bandpass filtering (1–20 Hz). Afterwards, the Fast Fourier Transform (FFT), Continuous Wavelet Transform (CWT), and Short Time Fourier Transform (STFT) methods were used to transform the time-domain data of 1453 natural earthquake events and 1103 quarry blasting events in 2013, obtaining four different types of training sample data: time-domain, frequency-domain (FFT results), and time–frequency domain (CWT and STFT results). Next, the four types of sample data were trained and tested using the Efficient Channel Attention Convolutional Network (CNN-ECA) and traditional Convolutional Neural Network (CNN). The results showed that the CNN-ECA model outperformed the CNN model in all four test samples. Especially when using time–frequency data converted through STFT and FFT as input, the recognition performance of the network model is more significant, with test set accuracies reaching 97.94% and 97.80%, respectively. Finally, the trained CNN-ECA model was used to validate and analyze the natural earthquakes and quarry blasting events recorded between 2014 and 2017. The results indicated that the combined use of FFT and STFT/CWT input data to jointly discriminate seismic events further improved the accuracy of earthquake type identification.</p></div>","PeriodicalId":16994,"journal":{"name":"Journal of Seismology","volume":"29 4","pages":"795 - 812"},"PeriodicalIF":2.0,"publicationDate":"2025-07-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145236811","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}
Richard Lewerissa, Sismanto, Jan Mrlina, Laura A. S. Lapono
{"title":"Integrating seismicity parameter analysis and gravity anomaly characterization to enhance the accuracy of earthquake hazard assessment: a case study in Cenderawasih Bay, Indonesia","authors":"Richard Lewerissa, Sismanto, Jan Mrlina, Laura A. S. Lapono","doi":"10.1007/s10950-025-10311-1","DOIUrl":"10.1007/s10950-025-10311-1","url":null,"abstract":"<div><p>In Cenderawasih Bay, Papua, we improved the assessment of earthquake hazards by integrating seismic parameter analysis with gravity anomaly data. This region, which is known for its seismic activity and complex tectonic structures, offers an ideal setting for disaster risk-mitigation studies. In eastern Indonesia, tectonic plate convergence involves several microplates, creating deformation zones with arc-continent collisions, subduction, shear faults, strike-slip faults, and extensional faults. Our study examined the relationship between earthquake activity and gravity anomalies to enhance the understanding of seismic activity distribution in the Cenderawasih Bay region. Seismicity parameters were analyzed using the IRIS Earthquake Catalog, with aftershocks and foreshocks removed through Reasenberg declustering to focus on independent events. The average values for the seismicity parameters, specifically the magnitudes of completeness (<i>M</i><sub>c</sub>), a, and b, were 4.5, 7.3, and 0.95, respectively, indicating significant tectonic activity. Gravity anomalies were analyzed using a two-dimensional radial power spectrum and three-dimensional Euler deconvolution, which resulted in an estimated depth range for the gravity anomaly source between 15 and 33 km. The results demonstrate that seismic activity is primarily concentrated in the northern region and along the main fault lines of Cenderawasih Bay. There is a correlation between the seismicity parameter <i>b</i>-value and variations in gravity anomalies, with low <i>b</i>-values associated with low to moderate gravity anomalies in regions of high crustal pressure. In contrast, elevated <i>b</i>-values are associated with significant gravity anomalies in low-pressure regions. The seismicity parameter <i>a</i>-value was low to moderate, suggesting tectonic activity. High values correlate with significant gravity anomalies, and the inverse is true. Integrating seismicity analysis with gravity anomaly characterization provides insights into the earthquake distribution in Cenderawasih Bay. This methodological approach not only improves the accuracy of earthquake hazard assessments in the region, but also aids effective disaster risk mitigation.</p></div>","PeriodicalId":16994,"journal":{"name":"Journal of Seismology","volume":"29 4","pages":"903 - 931"},"PeriodicalIF":2.0,"publicationDate":"2025-07-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145236762","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}
Adrianto Widi Kusumo, Hiroyuki Azuma, Toshiki Watanabe, Yoshiya Oda
{"title":"Seismic tomography for subsurface structures imaging beneath Hachijojima Volcanic Island, Izu-Bonin Arc, Japan","authors":"Adrianto Widi Kusumo, Hiroyuki Azuma, Toshiki Watanabe, Yoshiya Oda","doi":"10.1007/s10950-025-10309-9","DOIUrl":"10.1007/s10950-025-10309-9","url":null,"abstract":"<div><p>We present a seismic tomography study of the subsurface structure beneath Hachijojima Island, one of the volcanic fronts in the Izu-Bonin Arc, Japan. Seismic observations were conducted over two 7-month periods in 2019 and 2021, utilizing 55 densely installed stations on the island. During these periods, a total of 179 local earthquakes were recorded — 119 in 2019 and 60 in 2021 — resulting in 4671 P-wave arrival times and 3927 S-wave arrival times. The 3-D tomography, derived using the double-difference technique, revealed a shallow low-velocity region between the island’s two main volcanoes, Nishiyama and Higashiyama, suggesting the presence of volcanic sediments near the surface. Additionally, a high-velocity anomaly was identified at a depth of 4–5 km, extending vertically from deeper regions beneath Nishiyama. This feature is interpreted as a magma pathway from past volcanic activity, with high P-wave velocities and elevated Vp/Vs ratios indicating possible fluid presence. At greater depths, low P-wave velocity perturbations and elevated Vp/Vs ratios suggest a magmatic plumbing system comprising a mid-crustal magma chamber at approximately 8–12 km depth and lateral magmatic pathways at 10–20 km depth. Furthermore, a distinct zone characterized by reduced P-wave velocity and increased Vp/Vs is interpreted as a shallow magma chamber with H₂O-saturated magma accumulation. These findings provide valuable insights into the subsurface magmatic processes beneath Hachijojima Island, which are crucial for improving volcanic hazard assessment.</p></div>","PeriodicalId":16994,"journal":{"name":"Journal of Seismology","volume":"29 4","pages":"855 - 873"},"PeriodicalIF":2.0,"publicationDate":"2025-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10950-025-10309-9.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145236783","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Nicola Piana Agostinetti, Christina Dahner, Savka Dineva
{"title":"High-resolution temporal variations in rock elasticity at kiruna mine (block #30 to #34) using full 4D passive seismic tomography","authors":"Nicola Piana Agostinetti, Christina Dahner, Savka Dineva","doi":"10.1007/s10950-025-10308-w","DOIUrl":"10.1007/s10950-025-10308-w","url":null,"abstract":"<div><p>We investigate seismic velocity changes in the rock mass related to mining induced seismic events and ore exploitation by computing a one-month long 4D elastic model of Kiirunavaara mine (Sweden). We focus on a specific mine sector, where a single <span>(varvec{M_W})</span>=2.0 event occurred on May 22 (02:31 local time), damaging the infrastructure. We make use of P- and S-first-arrival times obtained from the permanent seismic system for computing the full 4D (continuous 3D volume in time) seismic velocity model of Kiruna mine using a trans-dimensional Monte Carlo sampling. The trans-dimensional approach guarantees that the resolution, both in space and in time, is strictly data-driven. Our results give the following insights into the velocity differences at the mining levels and at different time-length scales. (a) We observe a striking correlation between spatial variations of <span>(varvec{V_P})</span> and ore-body geometry, confirming the robustness of the velocity model. Clay zones appear as a low <span>(varvec{V_P/V_S})</span> ratio zones, as seen in previous tomographic studies. (b) High-frequency (hourly) fluctuations of the rock mass <span>(varvec{V_P})</span> around the ore-passes are highly correlated with seismic sequences in the same rock volumes. In particular, <span>(varvec{V_P})</span> increases rapidly when ore-passes are seismically active and <span>(varvec{V_P})</span> values keep a high value for few (1-4) hours after the end of the seismic sequence. (c) The smoothed velocity model, computed as averaged model over a 2-days moving window, suggests that low-frequency <span>(varvec{V_P})</span> fluctuations can be compared to stress cell measurements located closely.</p></div>","PeriodicalId":16994,"journal":{"name":"Journal of Seismology","volume":"29 4","pages":"835 - 853"},"PeriodicalIF":2.0,"publicationDate":"2025-07-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://link.springer.com/content/pdf/10.1007/s10950-025-10308-w.pdf","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145236782","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}