Journal of Geodesy最新文献

{"title":"The topside global broadcast ionospheric delay correction model for future LEO navigation augmentation","authors":"Yan Yang, Fei Guo, Chengpan Tang, Mengjie Wu, Kai Li, Xiaohong Zhang, Enyuan Tu","doi":"10.1007/s00190-024-01874-x","DOIUrl":"https://doi.org/10.1007/s00190-024-01874-x","url":null,"abstract":"<p>In this paper, we propose a solution of designing a topside broadcast ionospheric model to enable the future low earth orbit (LEO) navigation augmentation (LEO-NA) services. Considering the lack of global station observations to develop the LEO-NA ionosphere model, we utilize abundant global navigation satellite system (GNSS) data from LEO satellites to determine the topside global broadcast ionospheric delay. This delay can be combined with existing GNSS broadcast ionospheric delay correction models to determine LEO-NA ionospheric delay. First, the performance of the different-order spherical harmonic (SH) model is evaluated in generating a global topside ionospheric map. The results indicate that by increasing the order from 1 to 2, the internal and external accuracy of the model improves significantly. However, increasing the order from 2 to 8 leads to a decrease in accuracy of 0.10 and 0.11 TECU (total electron content unit) for the internal and external root mean square error. Taking into account compatibility with the Beidou global ionospheric delay correction model, limited data capacity in the navigation message, ionospheric model accuracy, and computational efficiency, we select the second-order SH model as the topside ionosphere broadcast model and outline the strategy for calculating broadcast coefficients. Finally, the accuracy of the topside global broadcast ionospheric delay correction model is evaluated during periods of high and low solar activity. The mean values of root mean square in 2009 and 2014 are 1.49 and 1.88 TECU, respectively. The model in 2009 and 2014 can correct for 67.30% and 72.49% of the ionospheric delay, respectively.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"25 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-07-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141557076","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":"A spatial-varying non-isotropic Gaussian-based convolution filter for smoothing GRACE-like temporal gravity fields","authors":"Fan Yang, Shuhao Liu, Ehsan Forootan","doi":"10.1007/s00190-024-01875-w","DOIUrl":"https://doi.org/10.1007/s00190-024-01875-w","url":null,"abstract":"<h3 data-test=\"abstract-sub-heading\">Abstract</h3><p>The strong noise of satellite-based Time-Variable Gravity (TVG) field is often suppressed by applying the averaging filters. However, how to appropriately compromise the data blurring and de-noising remains as a challenge. In our hypothesis, the optimum spatial averaging filter expects to contain averaging kernels that capture the same amount of orbital samples everywhere, to avoid introducing excessive data blurring. To achieve the goal, we take advantages of the spherical convolution and introduce extra spatial constraints into a Gaussian kernel: (1) its half-width radius adapts to the global inhomogeneity of satellite orbit, and (2) the kernel is reshaped as an ellipsoid to adapt to the regional anisotropy. In this way, we designed optimal filters that contain a spatially-Varying non-isotropic Gaussian-based Convolution (VGC) kernel. The VGC-based filter is compared against three most popular filters through real TVG fields and another closed-loop simulation. In both scenarios, VGC-based filters retain more realistic secular trend and seasonal characteristics, in particular at high latitudes. The spatial correlation between the VGC estimates and the simulated ground truth is found to be 0.95 and 0.86 over Greenland and Antarctica, which is found to be 10% better than other tested filters. Temporal correlations with the ground truth are also found to be considerably better than the other filters over 90% of the globally distributed river basin. Besides, the VGC-based filters provide tolerable efficiency (3.5 s per month) and sufficient accuracy (integral error less than 3%). The method can be extended to the next generation gravity mission as well.</p><h3 data-test=\"abstract-sub-heading\">Plain Language Summary</h3><p>Time-Variable Gravity (TVG) fields of the Gravity Recovery and Climate Experiment (GRACE) and its Follow-On mission (GRACE-FO) need proper filtering to suppress the noise before being applied for intended geophysical studies. Existing filters are generally designed in the spectral domain. Though they are numerically efficient, they can hardly treat the noise in fairness, globally. As a result, the TVG fields may get over-smoothed after applying those filters, particularly in regions with high-latitudes. However, it would be mathematically simple to design a filter by applying a spherical convolution, whose kernels can be easily constrained and tuned in the spatial domain. This study introduces filters with spatially-Varying non-isotropic Gaussian-based Convolution kernel (VGC) that is enforced to comply with the spatial distribution of the TVG noise. The proposed filter is found to preserve a finer spatial resolution of TVG fields, and at the same time, to be able to de-noise them at a comparable level as the existing techniques. Geophysical applications that use GRACE-like TVG fields might have benefits from this practical filtering technique.\u0000</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"18 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-07-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141553438","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":"A new spherical harmonic approach to residual terrain modeling: a case study in the central European Alps","authors":"Joachim Schwabe, Torsten Mayer-Gürr, Christian Hirt, Tobias Bauer","doi":"10.1007/s00190-024-01843-4","DOIUrl":"https://doi.org/10.1007/s00190-024-01843-4","url":null,"abstract":"<p>For decades, the residual terrain model (RTM) concept (Forsberg and Tscherning in J Geophys Res Solid Earth 86(B9):7843–7854, https://doi.org/10.1029/JB086iB09p07843, 1981) has been widely used in regional quasigeoid modeling. In the commonly used remove-compute-restore (RCR) framework, RTM provides a topographic reduction commensurate with the spectral resolution of global geopotential models. This is usually achieved by utilizing a long-wavelength (smooth) topography model known as reference topography. For computation points in valleys this neccessitates a harmonic correction (HC) which has been treated in several publications, but mainly with focus on gravity. The HC for the height anomaly only recently attracted more attention, and so far its relevance has yet to be shown also empirically in a regional case study. In this paper, the residual spherical-harmonic topographic potential (RSHTP) approach is introduced as a new technique and compared with the classic RTM. Both techniques are applied to a test region in the central European Alps including validation of the quasigeoid solutions against ground-truthing data. Hence, the practical feasibility and benefits for quasigeoid computations with the RCR technique are demonstrated. Most notably, the RSHTP avoids explicit HC in the first place, and spectral consistency of the residual topographic potential with global geopotential models is inherently achieved. Although one could conclude that thereby the problem of the HC is finally solved, there remain practical reasons for the classic RTM reduction with HC. In this regard, both intra-method comparison and ground-truthing with GNSS/leveling data confirms that the classic RTM (Forsberg and Tscherning 1981; Forsberg in A study of terrain reductions, density anomalies and geophysical inversion methods in gravity field modeling. Report 355, Department of Geodetic Sciences and Surveying, Ohio State University, Columbus, Ohio, USA, https://earthsciences.osu.edu/sites/earthsciences.osu.edu/files/report-355.pdf, 1984) provides reasonable results also for a high-resolution (degree 2160) RTM, yet neglecting the HC for the height anomaly leads to a systematic bias in deep valleys of up to 10–20 cm.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"22 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-07-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141546187","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":"Benefit of classical leveling for geoid-based vertical reference frames","authors":"Christian Gerlach, Reiner Rummel","doi":"10.1007/s00190-024-01849-y","DOIUrl":"https://doi.org/10.1007/s00190-024-01849-y","url":null,"abstract":"<p>Classically, vertical reference frames were realized as national or continent-wide networks of geopotential differences derived from geodetic leveling, i.e., from the combination of spirit leveling and gravimetry. Those networks are affected by systematic errors in leveling, leading to tilts in the order of decimeter to meter in larger networks. Today, there opens the possibility to establish a worldwide unified vertical reference frame based on a conventional (quasi)geoid model. Such a frame would be accessible through GNSS measurements, i.e., physical heights would be derived by the method of GNSS-leveling. The question arises, whether existing geodetic leveling data are abolished completely for the realization of vertical reference frames, are used for validation purposes only, or whether existing or future geodetic leveling data can still be of use for the realization of vertical reference frames. The question is mainly driven by the high quality of leveled potential differences over short distances. In the following we investigate two approaches for the combination of geopotential numbers from GNSS-leveling and potential differences from geodetic leveling. In the first approach, both data sets are combined in a common network adjustment leading to potential values at the benchmarks of the leveling network. In the second approach, potential differences from geodetic leveling are used as observable for regional gravity field modeling. This leads to a grid of geoid heights based on classical observables like gravity anomalies and now also on leveled potential differences. Based on synthetic data and a realistic stochastic model, we show that incorporating leveled potential differences improves the quality of a continent-wide network of GNSS-heights (approach 1) by about 40% and that formal and empirical errors of a regional geoid model (approach 2) are reduced by about 20% at leveling benchmarks. While these numbers strongly depend on the chosen stochastic model, the results show the benefit of using leveled potential differences for the realization of a modern geoid-based reference frame. Independent of the specific numbers of the improvement, an additional benefit is the consistency (within the error bounds of each observation type) of leveling data with vertical coordinates from GNSS and a conventional geoid model. Even though we focus on geodetic leveling, the methods proposed are independent of the specific technique used to observe potential (or equivalently height) differences and can thus be applied also to other techniques like chronometric or hydrodynamic leveling.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"18 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141462595","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":"Primal and dual mixed-integer least-squares: distributional statistics and global algorithm","authors":"P. J. G. Teunissen, L. Massarweh","doi":"10.1007/s00190-024-01862-1","DOIUrl":"https://doi.org/10.1007/s00190-024-01862-1","url":null,"abstract":"<p>In this contribution we introduce the dual mixed-integer least-squares problem and study it in relation to its primal counterpart. The dual differs from the primal formulation in the order in which the integer ambiguity vector <span>(a in {mathbb {Z}}^{n})</span> and baseline vector <span>(b in {mathbb {R}}^{p})</span> are estimated. As not the ambiguities, but rather the entries of <i>b</i> are usually the parameters of interest, the attractiveness of the dual formulation stems from its direct computation of <i>b</i>. It is shown that this potential advantage relies on the ease with which an implicit integer least-squares problem of the dual can be solved. For the convoluted cases, we introduce two methods of simplifying approximations. To be able to describe their quality, we provide a complete distributional analysis of their estimators, thus allowing users to judge whether or not the approximations are acceptable for their application. It is shown that this approach implicitly introduces a new class of admissible integer estimators of which we also determine the pull-in regions. As the dual function is shown to lack convexity, special care is required to be able to compute its global minimizer <span>({check{b}})</span>. Our proposed method, which has finite termination with a guaranteed <span>(epsilon )</span>-tolerance, is constructed from combining the branch-and-bound principle, with a special convex-relaxation of the dual, to which the projected-gradient-descent method is applied to obtain the required bounds. Each of the method’s three constituents are described, whereby special emphasis is given to the construction of the required continuously differentiable, convex lower bounding function of the dual.\u0000</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"43 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141448292","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":"Variance component adaptive estimation algorithm for coseismic slip distribution inversion using interferometric synthetic aperture radar data","authors":"Yingwen Zhao, Caijun Xu, Yangmao Wen","doi":"10.1007/s00190-024-01866-x","DOIUrl":"https://doi.org/10.1007/s00190-024-01866-x","url":null,"abstract":"<p>When conducting coseismic slip distribution inversion with interferometric synthetic aperture radar (InSAR) data, there is no universal method to objectively determine the appropriate size of InSAR data. Currently, little is also known about the computing efficiency of variance component estimation implemented in the inversion. Therefore, we develop a variance component adaptive estimation algorithm to determine the optimal sampling number of InSAR data for the slip distribution inversion. We derived more concise variation formulae than conventional simplified formulae for the variance component estimation. Based on multiple sampling data sets with different sampling numbers, the proposed algorithm determines the optimal sampling number by the changing behaviors of variance component estimates themselves. In three simulation cases, four evaluation indicators at low levels corresponding to the obtained optimal sampling number validate the feasibility and effectiveness of the proposed algorithm. Compared with the conventional slip distribution inversion strategy with the standard downsampling algorithm, the simulation cases and practical applications of five earthquakes suggest that the developed algorithm is more flexible and robust to yield appropriate size of InSAR data, thus provide a reasonable estimate of slip distribution. Computation time analyses indicate that the computational advantage of variation formulae is dependent of the ratio of the number of data to the number of fault patches and can be effectively suitable for cases with the ratio smaller than five, facilitating the rapid estimation of coseismic slip distribution inversion.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"46 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141436055","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":"Review of early ground deformation observations by electronic distance measurements (EDM) on active Sicilian volcanoes: valuable data and information for long-term analyses","authors":"Alessandro Bonforte, Salvatore Gambino, Rosanna Velardita, Laura Privitera","doi":"10.1007/s00190-024-01867-w","DOIUrl":"https://doi.org/10.1007/s00190-024-01867-w","url":null,"abstract":"<p>Electronic distance measurements (EDM) represent one of the first methods to detect ground deformation on volcanoes. Used since 1964, they enable acquiring precise distance measurements, whose time repetition may highlight changes related to volcanic activity. This technique was widely used on volcanoes from the 1970s to the early 2000s and has been used many times to model position, geometry, and volumes of magmatic and hydrothermal sources. This paper reports the EDM experiences, results and data acquired on Sicilian volcanoes (Etna, Vulcano, Stromboli and Pantelleria) from the early 1970s, which have played a major role in the birth of the volcano-geodesy for volcanic process knowledge, making the Sicilian volcanoes among those with the longest geodetic record in the world.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"26 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141430468","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":"The B-spline mapping function (BMF): representing anisotropic troposphere delays by a single self-consistent functional model","authors":"Shengping He, Thomas Hobiger, Doris Becker","doi":"10.1007/s00190-024-01864-z","DOIUrl":"https://doi.org/10.1007/s00190-024-01864-z","url":null,"abstract":"<p>Troposphere’s asymmetry can introduce errors ranging from centimeters to decimeters at low elevation angles, which cannot be ignored in high-precision positioning technology and meteorological research. The traditional two-axis gradient model, which strongly relies on an open-sky environment of the receiver, exhibits misfits at low elevation angles due to their simplistic nature. In response, we propose a directional mapping function based on cyclic B-splines named B-spline mapping function (BMF). This model replaces the conventional approach, which is based on estimating Zenith Wet Delay and gradient parameters, by estimating only four parameters which enable a continuous characterization of the troposphere delay across any directions. A simulation test, based on a numerical weather model, was conducted to validate the superiority of cyclic B-spline functions in representing tropospheric asymmetry. Based on an extensive analysis, the performance of BMF was assessed within precise point positioning using data from 45 International GNSS Service stations across Europe and Africa. It is revealed that BMF improves the coordinate repeatability by approximately <span>(10%)</span> horizontally and about <span>(5%)</span> vertically. Such improvements are particularly pronounced under heavy rainfall conditions, where the improvement of 3-dimensional root mean square error reaches up to <span>(13%)</span>.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"44 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141430396","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":"Continental and oceanic AAM contributions to Chandler Wobble with the amplitude attenuation from 2012 to 2022","authors":"Xue-Qing Xu, Ming Fang, Yong-Hong Zhou, Xin-Hao Liao","doi":"10.1007/s00190-024-01872-z","DOIUrl":"https://doi.org/10.1007/s00190-024-01872-z","url":null,"abstract":"<p>We reconstructed the Chandler Wobble (CW) from 1962 to 2022 by combining the Eigen-oscillator excited by geophysical fluids of atmospheric and oceanic angular momentums (AAM and OAM). The mass and motion terms for the AAM are further divided with respect to the land and ocean domains. Particular attention is placed on the time span from 2012 to 2022 in relation to the observable reduction in the amplitude of the CW. Our research indicates that the main contributor to the CW induced by AAM is the mass term (i.e., the pressure variations over land). Moreover, the phase of the AAM-induced CW remains relatively stable during the interval of 1962–2022. In contrast, the phase of the OAM-induced CW exhibits a periodic variation with a cycle of approximately 20 years. This cyclic variation would modulate the overall amplitude of the CW. The noticeable amplitude deduction over the past ten years can be attributed to the evolution of the CW driven by AAM and OAM, toward a state of cancellation. These findings further reveal that the variation in the phase difference between the CW forced by AAM and OAM, may be indicative of changes in the interaction between the solid Earth, atmosphere, and ocean.</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"15 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141425412","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":"An improved parameter filtering approach for processing GRACE gravity field models using first-order Gauss–Markov process","authors":"Lin Zhang, Yunzhong Shen, Qiujie Chen, Kunpu Ji","doi":"10.1007/s00190-024-01871-0","DOIUrl":"https://doi.org/10.1007/s00190-024-01871-0","url":null,"abstract":"<p>Removing stripe noise from the GRACE (Gravity Recovery and Climate Experiment) monthly gravity field model is crucial for accurately interpreting temporal gravity variations. The conventional parameter filtering (CPF) approach expresses the signal components with a harmonic model while neglecting non-periodic and interannual signals. To address this issue, we improve the CPF approach by incorporating those ignored signals using a first-order Gauss–Markov process. The improved parameter filtering (IPF) approach is used to filter the monthly spherical harmonic coefficients (SHCs) of the Tongji-Grace2018 model from April 2002 to December 2016. Compared to the CPF approach, the IPF approach exhibits stronger signals in low-degree SHCs (i.e., degrees below 20) and lower noise in high-order SHCs (i.e., orders above 40), alongside higher signal-to-noise ratios and better agreement with CSR mascon product and NOAH model in global and basin analysis. Across the 22 largest basins worldwide, the average Nash–Sutcliffe coefficients of latitude-weighted terrestrial water storage anomalies filtered by the IPF approach relative to those derived from CSR mascon product and NOAH model are 0.90 and 0.21, significantly higher than 0.17 and − 0.71, filtered by the CPF approach. Simulation experiments further demonstrate that the IPF approach yields the filtered results closest to the actual signals, reducing root-mean-square errors by 30.1%, 25.9%, 45.3%, 30.9%, 46.6%, 32.7%, 39.6%, and 38.2% over land, and 2.8%, 54.4%, 70.1%, 15.3%, 69.2%, 46.5%, 40.4%, and 23.6% over the ocean, compared to CPF, DDK3, least square, RMS, Gaussian 300, Fan 300, Gaussian 300 with P4M6, and Fan 300 with P4M6 filtering approaches, respectively</p>","PeriodicalId":54822,"journal":{"name":"Journal of Geodesy","volume":"2014 1","pages":""},"PeriodicalIF":4.4,"publicationDate":"2024-06-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141334130","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}