Earth and Planetary Physics最新文献

{"title":"Teleseismic waves reveal anisotropic poroelastic response of wastewater disposal reservoir","authors":"Andrew J Barbour, Nicholas M Beeler","doi":"10.26464/epp2021034","DOIUrl":"10.26464/epp2021034","url":null,"abstract":"Connecting earthquake nucleation in basement rock to fluid injection in basal, sedimentary reservoirs, depends heavily on choices related to the poroelastic properties of the fluid‐rock system, thermo‐chemical effects notwithstanding. Direct constraints on these parameters outside of laboratory settings are rare, and it is commonly assumed that the rock layers are isotropic. With the Arbuckle wastewater disposal reservoir in Osage County, Oklahoma, high‐frequency formation pressure changes and collocated broadband ground velocities measured during the passing of large teleseismic waves show a poroelastic response of the reservoir that is both azimuthally variable and anisotropic; this includes evidence of static shifts in pressure that presumably relate to changes in local permeability. The azimuthal dependence in both the static response and shear coupling appears related to tectonic stress and strain indicators such as the orientations of the maximum horizontal stress and faults and fractures. Using dynamic strains from a nearby borehole strainmeter, we show that the ratio of shear to volumetric strain coupling is ~0.41 which implies a mean Skempton's coefficient ofA = 0.24 over the plausible range of the undrained Poisson's ratio. Since these observations are made at relatively low confining pressure and differential stress, we suggest that the hydraulically conductive fracture network is a primary control on the coupling between pore pressure diffusion and elastic stresses in response to natural or anthropogenic sources.","PeriodicalId":45246,"journal":{"name":"Earth and Planetary Physics","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2021-12-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.26464/epp2021034","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41995927","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"3D geomechanical modeling of the response of the Wilzetta Fault to saltwater disposal","authors":"Behzad Hemami,&nbsp;Shahla Feizi Masouleh,&nbsp;Ahmad Ghassemi","doi":"10.26464/epp2021054","DOIUrl":"10.26464/epp2021054","url":null,"abstract":"<p>From 2009 to 2017, parts of Central America experienced marked increase in the number of small to moderate-sized earthquakes. For example, three significant earthquakes (~<i>M</i>_w 5) occurred near Prague, Oklahoma, in the U. S. in 2011. On 6 Nov 2011, an <i>M</i>_w 5.7 earthquake occurred in Prague, central Oklahoma with a sequence of aftershocks. The seismic activity has been attributed to slip on the Wilzetta fault system. This study provides a 3D fully coupled poroelastic analysis (using FLAC3D) of the Wilzetta fault system and its response to saltwater injection in the underpressured subsurface layers, especially the Arbuckle group and the basement, to evaluate the conditions that might have led to the increased seismicity. Given the data-limited nature of the problem, we have considered multiple plausible scenarios, and use the available data to evaluate the hydromechanical response of the faults of interest in the study area. Numerical simulations show that the injection of large volumes of fluid into the Arbuckle group tends to bring the part of the Wilzetta faults in Arbuckle group and basement into near-critical conditions.</p>","PeriodicalId":45246,"journal":{"name":"Earth and Planetary Physics","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2021-12-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.26464/epp2021054","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44034792","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Development of a new high resolution waveform migration location method and its applications to induced seismicity","authors":"SiYu Miao,&nbsp;HaiJiang Zhang,&nbsp;YuYang Tan,&nbsp;Ye Lin","doi":"10.26464/epp2021056","DOIUrl":"10.26464/epp2021056","url":null,"abstract":"<p>Locating seismic events is a central task for earthquake monitoring. Compared to arrival-based location methods, waveform-based location methods do not require picking phase arrivals and are more suitable for locating seismic events with noisy waveforms. Among waveform-based location methods, one approach is to stack different attributes of P and S waveforms around arrival times corresponding to potential event locations and origin times, and the maximum stacking values are assumed to indicate the correct event location and origin time. In this study, to obtain a high-resolution location image, we improve the waveform-based location method by applying a hybrid multiplicative imaging condition to characteristic functions of seismic waveforms. In our new stacking method, stations are divided into groups; characteristic functions of seismic waveforms recorded at stations in the same group are summed, and then multiplied among groups. We find that this approach can largely eliminate the cumulative effects of noise in the summation process and thus improve the resolution of location images. We test the new method and compare it to three other stacking methods, using both synthetic and real datasets that are related to induced seismicity occurring in petroleum/gas production. The test results confirm that the new stacking method can provide higher-resolution location images than those derived from currently used methods.</p>","PeriodicalId":45246,"journal":{"name":"Earth and Planetary Physics","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2021-12-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.26464/epp2021056","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49147610","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Study on electron stochastic motions in the magnetosonic wave field: Test particle simulations","authors":"Kai Fan,&nbsp;XinLiang Gao,&nbsp;QuanMing Lu,&nbsp;Shui Wang","doi":"10.26464/epp2021052","DOIUrl":"10.26464/epp2021052","url":null,"abstract":"<p>Using the test particle simulation method, we investigate the stochastic motion of electrons with energy of 300 keV in a monochromatic magnetosonic (MS) wave field. This study is motivated by the violation of the quasi-linear theory assumption, when strong MS waves (amplitude up to ~1 nT) are present in the Earth’s magnetosphere. First, electron motion can become stochastic when the wave amplitude exceeds a certain threshold. If an electron initially resonates with the MS wave via bounce resonance, as the bounce resonance order increases, the amplitude threshold of electron stochastic motion increases until it reaches the peak at about the 11th order in our study, then the amplitude threshold slowly declines. Further, we find that the coexistence of bounce and Landau resonances between electrons and MS waves will significantly reduce the amplitude threshold. In some cases, the electron motion can become stochastic in the field of an MS wave with amplitudes below 1 nT. Regardless, if neither the bounce nor Landau resonance condition is satisfied initially, then the amplitude threshold of stochastic motion shows an increasing trend for lower frequencies and a decreasing trend for higher frequencies, even though the amplitude threshold is always very large (&gt; 5 nT). Our study suggests that electron stochastic motion should also be considered when modeling electron dynamics regulated by intense MS waves in the Earth’s magnetosphere.</p>","PeriodicalId":45246,"journal":{"name":"Earth and Planetary Physics","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2021-12-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.26464/epp2021052","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41761879","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Non-storm erosion of MeV electron outer radiation belt down to L* < 4.0 associated with successive enhancements of solar wind density","authors":"Ying Xiong,&nbsp;Lun Xie,&nbsp;SuiYan Fu,&nbsp;BinBin Ni,&nbsp;ZuYin Pu","doi":"10.26464/epp2021051","DOIUrl":"10.26464/epp2021051","url":null,"abstract":"<p>We report an unusual non-storm erosion event of outer zone MeV electron distribution during three successive solar wind number density enhancements (SWDEs) on November 27−30, 2015. Loss of MeV electrons and energy-dependent narrowing of electron pitch angle distributions (PAD) first developed at <i>L</i>* = 5.5 and then moved down to <i>L</i>* &lt; 4. According to the evolution of the electron phase space density (PSD) profile, losses of electrons with small pitch angles at <i>L</i>* &gt; 4 during SWDE1 are mainly due to outward radial diffusion. However during SWDE2&amp;3, scattering loss due to EMIC waves is dominant at 4 &lt; <i>L</i>* &lt; 5. As for electrons with large pitch angles, outward radial diffusion is the primary loss mechanism throughout all SWDEs which is consistent with the incursion of the Last Closed Drift Shell (LCDS). The inner edge of EMIC wave activity moved from <i>L</i>* ~5 to <i>L</i>* ~4 and from <i>L</i> ~6.4 to <i>L</i> ~4.2 from SWDE1 to SWDE2&amp;3, respectively, observed by Van Allen Probes and by ground stations. This is consistent with the inward penetration of anisotropic energetic protons from<i>L</i>* = 4.5 to <i>L</i>* = 3.5, suggesting that the inward extension of EMIC waves may be driven by the inward injection of anisotropic energetic protons from the dense plasma sheet.</p>","PeriodicalId":45246,"journal":{"name":"Earth and Planetary Physics","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2021-12-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.26464/epp2021051","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"47822491","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Spatial-temporal evolution of injection-induced earthquakes in the Weiyuan Area determined by machine-learning phase picker and waveform cross-correlation","authors":"Wing Ching Jeremy Wong,&nbsp;JinPing Zi,&nbsp;HongFeng Yang,&nbsp;JinRong Su","doi":"10.26464/epp2021055","DOIUrl":"10.26464/epp2021055","url":null,"abstract":"<p>Anthropogenic induced seismicity has been widely reported and investigated in many regions, including the shale gas fields in the Sichuan basin, where the frequency of earthquakes has increased substantially since the commencement of fracking in late 2014. However, the details of how earthquakes are induced remain poorly understood, partly due to lack of high-resolution spatial-temporal data documenting the evolution of such seismic events. Most previous studies have been based on a diffusive earthquake catalog constructed by routine methods. Here, however, we have constructed a high resolution catalog using a machine learning detector and waveform cross-correlation. Despite limited data, this new approach has detected one-third more earthquakes and improves the magnitude completeness of the catalog, illuminating the comprehensive spatial-temporal migration of the emerging seismicity in the target area. One of the clusters clearly delineates a potential unmapped fault trace that may have led to the <i>M</i>_w 5.2 in September 2019, by far the largest earthquake recorded in the region. The migration of the seismicity also demonstrates a pore-pressure diffusion front, suggesting additional constraints on the inducing mechanism of the region. The patterns of the highly clustered seismicity reconcile the causal link between the emerging seismicity and the activity of hydraulic fracturing in the region, facilitating continued investigation of the mechanisms of seismic induction and their associated risks.</p>","PeriodicalId":45246,"journal":{"name":"Earth and Planetary Physics","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2021-12-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.26464/epp2021055","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"41947558","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Determination of the local magnitudes of small earthquakes using a dense seismic array in the Changning–Zhaotong Shale Gas Field, Southern Sichuan Basin","authors":"Wen Yang,&nbsp;GuoYi Chen,&nbsp;LingYuan Meng,&nbsp;Yang Zang,&nbsp;HaiJiang Zhang,&nbsp;JunLun Li","doi":"10.26464/epp2021026","DOIUrl":"10.26464/epp2021026","url":null,"abstract":"<p>With the development of unconventional shale gas in the southern Sichuan Basin, seismicity in the region has increased significantly in recent years. Though the existing sparse regional seismic stations can capture most earthquakes with , a great number of smaller earthquakes are often omitted due to limited detection capacity. With the advent of portable seismic nodes, many dense arrays for monitoring seismicity in the unconventional oil and gas fields have been deployed, and the magnitudes of those earthquakes are key to understand the local fault reactivation and seismic potentials. However, the current national standard for determining the local magnitudes was not specifically designed for monitoring stations in close proximity, utilizing a calibration function with a minimal resolution of 5 km in the epicentral distance. That is, the current national standard tends to overestimate the local magnitudes for stations within short epicentral distances, and can result in discrepancies for dense arrays. In this study, we propose a new local magnitude formula which corrects the overestimated magnitudes for shorter distances, yielding accurate event magnitudes for small earthquakes in the Changning−Zhaotong shale gas field in the southern Sichuan Basin, monitored by dense seismic arrays in close proximity. The formula is used to determine the local magnitudes of 7,500 events monitored by a two-phased dense array with several hundred 5 Hz 3C nodes deployed from the end of February 2019 to early May 2019 in the Changning−Zhaotong shale gas field. The magnitude of completeness () using the dense array is −0.1, compared to 1.1 by the sparser Chinese Seismic Network (CSN). In addition, using a machine learning detection and picking procedure, we successfully identify and process some 14,000 earthquakes from the continuous waveforms, a ten-fold increase over the catalog recorded by CSN for the same period, and the is further reduced to −0.3 from −0.1 compared to the catalog obtained via manual processing using the same dense array. The proposed local magnitude formula can be adopted for calculating accurate local magnitudes of future earthquakes using dense arrays in the shale gas fields of the Sichuan Basin. This will help to better characterize the local seismic risks and potentials.</p>","PeriodicalId":45246,"journal":{"name":"Earth and Planetary Physics","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2021-12-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.26464/epp2021026","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"49599359","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Multi-instrument study of longitudinal wave structures for plasma bubble seeding in the equatorial ionosphere","authors":"H. Takahashi,&nbsp;P. Essien,&nbsp;C. A. O. B. Figueiredo,&nbsp;C. M. Wrasse,&nbsp;D. Barros,&nbsp;M. A. Abdu,&nbsp;Y. Otsuka,&nbsp;K. Shiokawa,&nbsp;GuoZhu Li","doi":"10.26464/epp2021047","DOIUrl":"10.26464/epp2021047","url":null,"abstract":"<p>Large Scale Wave Structures (LSWS) in the equatorial ionospheric F-region were observed by measuring spatial and temporal variations within detrended total electron content (dTEC) data obtained by ground-based GNSS receivers over the South American continent. By using dTEC-maps, we have been able to produce, for the first-time, two-dimensional representations of LSWS. During the period from September to December, the LSWS frequently occurred starting a few hours prior to Equatorial Plasma Bubble (EPB) development. From 17 events of LSWS observed in 2014 and 2015, wave characteristics were obtained: the observed wavelengths, periods, and the phase speeds are respectively, ~900 km, ~41 min and ~399 m/s; the waves propagated from the northeast to southeast. In some cases the front of the oscillation was meridionally aligned, extending to more than 1600 km, the first time such large extension of the wavefront has been reported. From F-layer bottom height oscillation data, measured by ionosonde, LSWS exhibit two different vertical phase propagation modes, in-phase and downward phase. The former mode indicates the presence of a polarization electric field in the F-layer bottom side; the latter suggests propagation of atmospheric gravity waves. The presence of LSWS near the solar terminator, followed by the development of EPBs, suggests that the upwelling of the F-layer bottom height produces a condition favorable to the development of Rayleigh–Taylor instability.</p>","PeriodicalId":45246,"journal":{"name":"Earth and Planetary Physics","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2021-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.26464/epp2021047","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"43801889","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Case study of an Equatorial Plasma Bubble Event investigated by multiple ground-based instruments at low latitudes over China","authors":"LongChang Sun,&nbsp;JiYao Xu,&nbsp;YaJun Zhu,&nbsp;Wei Yuan,&nbsp;XiuKuan Zhao","doi":"10.26464/epp2021048","DOIUrl":"10.26464/epp2021048","url":null,"abstract":"<p>Observational evidence is insufficient to understand how equatorial plasma bubbles (EPBs) form over low latitudes. The mechanism of plasma-density enhancement (formation of “plasma blobs”) at low latitudes is in dispute. In this paper, we use data from multiple ground-based instruments (one all-sky airglow imager, five digisondes, and one Fabry–Perot interferometer) to investigate the evolution of an EPB event that occurred at low latitudes over China on the night of 06 December 2015 (06-Dec-2015). We provide observational evidence that an enhanced equatorward wind most likely induced by a substorm could have initiated the Rayleigh–Taylor instability (RTI) that destabilized several EPB depletions in an upwelling region of a large-scale wave-like structure (LSWS) in the bottomside ionosphere. Those EPB depletions were forced to surge poleward, from nearly 10° to 19° magnetic latitude, two hours before midnight. Smaller-scale bifurcations evolved rapidly from tips of airglow depletions by a secondary<i><b>E</b></i> × <i><b>B</b></i> instability when the aforementioned substorm-induced southwestward wind blew through. During the growth phase of the EPB depletions, a westward polarization electric field inside the LSWS is likely to have compressed plasma downward, inducing the two airglow-type blobs observed in the bottomside ionosphere, by a mechanism of LSWS-blob connection that we propose. We also provide observational evidence of brightness airglow depletions. We find that an enhanced poleward wind associated with a passing-by brightness wave (BW) is likely to have transported plasma to fill the airglow depletions, which finally evolved into brightness airglow structures. This study investigates the physical processes accompanied by the EPB event and those two-airglow blobs observed at low-latitudes over China.</p>","PeriodicalId":45246,"journal":{"name":"Earth and Planetary Physics","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2021-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.26464/epp2021048","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"42830378","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Plasma depletions lasting into daytime during the recovery phase of a geomagnetic storm in May 2017: Analysis and simulation of GPS total electron content observations","authors":"Yuichi Otsuka,&nbsp;Atsuki Shinbori,&nbsp;Takuya Sori,&nbsp;Takuya Tsugawa,&nbsp;Michi Nishioka,&nbsp;Joseph D. Huba","doi":"10.26464/epp2021046","DOIUrl":"10.26464/epp2021046","url":null,"abstract":"<p>This paper reports that plasma density depletions appearing at middle latitudes near sunrise survived until afternoon on 29 May 2017 during the recovery phase of a geomagnetic storm. By analyzing GPS data collected in Japan, we investigate temporal variations in the horizontal two-dimensional distribution of total electron content (TEC) during the geomagnetic storm. The SYM-<i>H</i> index reached −142 nT around 08 UT on 28 May 2017. TEC depletions extending up to approximately 38°N along the meridional direction appeared over Japan around 05 LT (LT = UT + 9 hours) on 29 May 2017, when TEC rapidly increased at sunrise due to the solar extreme ultraviolet (EUV) radiation. The TEC depletions appeared sequentially over Japan for approximately 8 hours in sunlit conditions. At 06 LT on 29 May, when the plasma depletions first appeared over Japan, the background TEC was enhanced to approximately 17 TECU, and then decreased to approximately 80% of the TEC typical of magnetically quiet conditions. We conclude that this temporal variation of background plasma density in the ionosphere was responsible for the persistence of these plasma depletions for so long in daytime. By using the Naval Research Laboratory: Sami2 is Another Model of the Ionosphere (SAMI2), we have evaluated how plasma production and ambipolar diffusion along the magnetic field may affect the rate of plasma depletion disappearance. Simulation shows that the plasma density increases at the time of plasma depletion appearance; subsequent decreases in the plasma density appear to be responsible for the long-lasting persistence of plasma depletions during daytime. The plasma density depletion in the top side ionosphere is not filled by the plasma generated by the solar EUV productions because plasma production occurs mainly at the bottom side of the ionosphere.</p>","PeriodicalId":45246,"journal":{"name":"Earth and Planetary Physics","volume":null,"pages":null},"PeriodicalIF":2.9,"publicationDate":"2021-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://agupubs.onlinelibrary.wiley.com/doi/epdf/10.26464/epp2021046","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"44140277","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}