首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
选取2018年中国大陆构造环境监测网络GNSS基准站及中国周边IGS站数据,基于CSR4、FES2004、GOT00、NAO99b等海潮负荷改正模型,采用单一变量的解算方式,分析不同海潮模型对我国GNSS基准站基线、三维坐标、时间序列等造成的影响。结果表明:利用不同海潮模型解算基线精度,水平方向相差不大于1 mm,垂向相差不大于3 mm;各测站NEU坐标精度均在毫米级以下,以耿马站时间序列为例,与其他模型相比,基于FES2004模型解算的各向偏差趋势波动较弱,其年序列图的精度及连续性、稳定性更佳。简言之,利用FES2004海潮模型解算GNSS站点位置,效果较好。  相似文献   

2.
Two accurately calibrated superconducting gravimeters (SGs) provide high quality tidal gravity records in three central European stations: C025 in Vienna and at Conrad observatory (A) and OSG050 in Pecný (CZ). To correct the tidal gravity factors from ocean loading effects we compared the load vectors from different ocean tides models (OTMs) computed with different software: OLFG/OLMP by the Free Ocean Tides Loading Provider (FLP), ICET and NLOADF. Even with the recent OTMs the mass conservation is critical but the methods used to correct the mass imbalance agree within 0.1 nm/s2. Although the different software agrees, FLP probably provides more accurate computations as this software has been optimised. For our final computation we used the mean load vector computed by FLP for 8 OTMs (CSR4, NAO99, GOT00, TPX07, FES04, DTU10, EOT11a and HAMTIDE). The corrected tidal factors of the 3 stations agree better than 0.04% in amplitude and 0.02° in phase. Considering the weighted mean of the three stations we get for O1 δc = 1.1535 ± 0.0001, for K1 δc = 1.1352 ± 0.0003 and for M2 δc = 1.1621 ± 0.0003. These values confirm previous ones obtained with 16 European stations. The theoretical body tides model DDW99/NH provides the best agreement for M2 (1.1620) and MATH01/NH for O1 (1.1540) and K1 (1.1350). The largest discrepancy is for O1 (0.05%). The corrected phase αc does not differ significantly from zero except for K1 and S2. The calibrations of the two SG's are consistent within 0.025% and agree with Strasbourg results within 0.05%.  相似文献   

3.
Satellite Laser Ranging (SLR) to LAGEOS has a remarkable contribution to high-precise geodesy and geodynamics through deriving and validating various global geophysical models. This paper validates ocean tide models based on the analysis of satellite altimetry data, coastal tide gauges, and hydrodynamic data, i.e., CSR3.0, TOPEX4.0, CSR4.0A, FES2004, GOT00.2, and the CSRC Schwiderski model. LAGEOS orbits and SLR observation residuals from solutions based on different ocean tide models are compared and examined. It is found that LAGEOS orbits are sensitive to tidal waves larger than 5 mm. The analysis of the aliasing periods of LAGEOS orbits and tidal waves reveals that, in particular, the tidal constituent S2 is not well established in the recent ocean tide models. Some of the models introduce spurious peaks to empirical orbit parameters, which can be associated with S2, Sa, and K2 tidal constituents, and, as a consequence, can be propagated to fundamental parameters derived from LAGEOS observations.  相似文献   

4.
This paper describes the GPS applicability for detecting the vertical displacements of ground stations caused by ocean tide loading effects. An experiment was carried out using 12 permanent GPS stations located in the coastal area of Korea using data in the period 1 July until 26 August 2003. The relative height differences were calculated from hourly DGPS data processing based on the carrier-phase observation. The power spectra of the M2 and N2 constituents of ocean tide loading were derived using the CLEAN algorithm. The differential vertical displacements generated by the ocean tide loading effect are typically 3–25 mm in coastal area of the Korea. We compared the results from GPS with those of the ocean tide models, NAO.99Jb regional model and GOT00.2, FES99 global models. The M2 (N2) amplitude differences of vertical displacements between GPS and GOT00.2 is 1.22 ± 3.61 mm (1.01 ± 1.48 mm), and that of the M2 (N2) amplitude difference between GPS and FES99 is 0.04 ± 4.64 mm (0.64 ± 1.75 mm), whereas the M2 (N2) amplitude difference between GPS and NAO.99Jb are 0.05 ± 1.03 mm (0.86 ± 1.18 mm). The highest vertical displacements at the PALM station are found for 24.5 ± 0.7 mm from GPS observation, and 22.9 mm from the regional model NAO.99Jb and 13.17 and 10.00 mm from the global models GOT00.2 and FES99, respectively. These values show that the vertical displacements derived from GPS are in good agreement with those of the regional model NAO.99Jb around Korea, more than with the global models. This result indicated that GPS is an effective tool to measure the vertical displacement caused by the ocean tide loading effect in the coastal area, and we need to use the NAO.99Jb ocean tide model rather than the global ocean tide models in and around the Korean peninsula for position determination with permanent GPS installations. This work demonstrates that vertical displacement caused by the M2 and N2 constituents of ocean tide loading can be measured by carrier-phase DGPS.  相似文献   

5.
Introduction With the development of science and technology, the accuracy of gravity measurement is im-proved. The gravity observation with FG5 gravimeter has reached accuracy of μ magnitude. And the high accurate superconducting gravimeter can detect the tiny signal of 10?2 μ magnitude in frequency domain. With the high-accuracy gravity observation on Earth′s surface, the Earth′s tidal parameters can be determined precisely. And the observations can also be used to invert the struc-ture…  相似文献   

6.
本文探索了海潮负荷特征值的反演方法.基于中国近岸海岛GPS站(平潭与闸坡)观测数据,采用FFT方法提取了海潮负荷特征值,分析了反演特征值与全球海潮模型FES2004、NAO.99b和GOT4.7计算出的相关特征值之间的差异,评估了反演海潮分潮频率项的精度,并利用平潭站得到的海潮负荷特征值对附近的三沙站进行海潮改正以评价反演效果.结果表明:(1)频谱分析可精确提取4个半日分潮和4个全日分潮负荷的频率信息;与已知频率相比,P1分潮的反演频率误差为1.4%,其他7种分潮负荷反演频率误差均小于1%.(2)两个海岛分属不同的潮波入侵通道,反演分潮振幅和初相存在差异,但反演分潮频率几乎一致,间接证实它们属于同一潮波系统,也表明其潮差有别.(3)反演振幅与三种全球模型具有较好的一致性;其中S2、O1、P1、Q1四个分潮在水平方向互差为1~2 mm,高程方向上的互差均小于3 mm;K2、K1、M2、N2在水平方向振幅互差多数小于2 mm,个别差异高达4 mm,高程方向互差多为5~6 mm,个别超过10 mm.(4)反演得到海潮负荷改正模型相对于3种全球模型在三沙站的改正效果略佳,间接表明反演结果有效、可靠.(5)动态PPP结果中虽然存在多种误差,其时间序列仍可分离并提取海潮负荷的影响.  相似文献   

7.
The seasonal cycle of the main lunar tidal constituent M 2 is studied globally by an analysis of a high-resolution ocean circulation and tide model (STORMTIDE) simulation, of 19 years of satellite altimeter data, and of multiyear tide-gauge records. The barotropic seasonal tidal variability is dominant in coastal and polar regions with relative changes of the tidal amplitude of 5–10 %. A comparison with the observations shows that the ocean circulation and tide model captures the seasonal pattern of the M 2 tide reasonably well. There are two main processes leading to the seasonal variability in the barotropic tide: First, seasonal changes in stratification on the continental shelf affect the vertical profile of eddy viscosity and, in turn, the vertical current profile. Second, the frictional effect between sea-ice and the surface ocean layer leads to seasonally varying tidal transport. We estimate from the model simulation that the M 2 tidal energy dissipation at the sea surface varies seasonally in the Arctic (ocean regions north of 60°N) between 2 and 34 GW, whereas in the Southern Ocean, it varies between 0.5 and 2 GW. The M 2 internal tide is mainly affected by stratification, and the induced modified phase speed of the internal waves leads to amplitude differences in the surface tide signal of 0.005–0.0150 m. The seasonal signals of the M 2 surface tide are large compared to the accuracy demands of satellite altimetry and gravity observations and emphasize the importance to consider seasonal tidal variability in the correction processes of satellite data.  相似文献   

8.
We have estimated the parameters of fluid core resonance (FCR) due to the nearly diurnal free wobble of the Earth's core based on the superconducting gravimeter (SG) data obtained at the following four observation sites; Esashi and Matsushiro in Japan, Canberra in Australia and Membach in Belgium. By fitting the tidal admittances normalized with the O1 wave at each site to a model of the damped harmonic oscillator, we obtained values of 429.66 ± 1.43 sidereal days, 9350–10,835, −4.828E−4 ± 3.4E−6, −3.0E−5 ± 4.5E−6 for the eigenperiod, the Q-value and the real and imaginary parts of the resonance strength, respectively. Our values obtained from only using the gravity data are very consistent with those inferred from the VLBI nutation data. Our study strongly indicates that the systematic difference between two estimations from the gravity and the nutation in particular for the Q-value, which has been shown in previous works, is mainly caused by the inaccurate correction for the ocean tide effects. The error in the ocean tide correction is discussed based on the comparison among four global ocean tide models; Schwiderski model (1980), NAO.99b (Matsumoto et al., 2000), CSR4.0 (Eanes and Bettadpur, 1994) and GOT99.2b (Ray, 1999).  相似文献   

9.
The daily samples of the hourly measurements of the foF2 critical frequency, obtained on January 5–21, 1995, at the midlatitude and high-latitude automated ionospheric stations (geographic latitude higher than 60°), are considered. The {fo} sets are transformed into the {δfoF2} sets of relative variations, for which asymmetry (A) and excess (E) are calculated. The selected stations are grouped into 20 pairs of automated ionospheric stations (AISs) located at distances of 200–10 000 km from one another. Sign estimates of the cross-correlation coefficients between the sets of 16 A and E values for different pairs of stations are used. Two types of structures of the statistical invariant spatial distribution are established: the structures with a scale of about 300 km, invariant with respect to latitude, and with a scale of about 6000 km (for only high latitudes).  相似文献   

10.
The deep structure of the upper mantle is determined from data on phase velocities of Love and Rayleigh waves measured by a differential method on traces between two stations in central Western Europe. One-dimensional velocity structures are first constructed from data of each pair of stations, after which two-dimensional distributions of SH and SV velocities are calculated by the method of two-dimensional tomography from S wave velocities at fixed depths. The results are presented in the form of 2-D vertical structures of the average S wave velocity (S = (SV + SH)/2) constructed along profiles crossing the region in directions of the best resolution. The main structural features are a higher velocity zone at depths of 60–80 km in the area (48°–50°N, 9°–11°E) and a lower velocity zone in the western part of the region at depths of 100–150 km, probably extending farther beyond the studied area.  相似文献   

11.
The present study reports the analysis of GPS TEC prior to 3 earthquakes (M > 6.0). The earthquakes are: (1) Loyalty Island (22°36′S, 170°54′E) on 19 January 2009 (M = 6.6), (2) Samoa Island (15°29′S, 172°5′W) on 30 August 2009 (M = 6.6), and (3) Tohoku (38°19′N, 142°22′E) on 11 March 2011 (M = 9.0). In an effort to search for a precursory signature we analysed the land and ocean parameters prior to the earthquakes, namely SLHF (Land) and SST (Ocean). The GPS TEC data indicate an anomalous behaviour from 1–13 days prior to earthquakes. The main purpose of this study was to explore and demonstrate the possibility of any changes in TEC, SST, and SLHF before, during and after the earthquakes which occurred near or beneath an ocean. This study may lead to better understanding of response of land, ocean, and ionosphere parameters prior to seismic activities.  相似文献   

12.
武汉超导重力仪观测最新结果和海潮模型研究   总被引:18,自引:5,他引:13       下载免费PDF全文
利用武汉台站GWR_C032超导重力仪观测资料,在对原始数据进行有效预处理的基础上作调和分析,获得反映地球内部介质特征的重力潮汐参数.基于卫星测高技术和有限元方法同时考虑验潮站数据作约束条件获得的多个全球海潮模型,利用负荷理论和数值褶积积分技术计算了重力负荷,对周日和半日频段内的重力潮汐参数实施负荷改正,提出了“负荷改正有效性”概念,研究了全球海潮模型适应性.数值结果说明,海潮改正的有效性高达91%(O1,NAO99)和92%(M2,ORI96).基于11个海潮模型对主波(O1,K1,M2和S2)的负荷改正说明平均有效性为(86%,70%,73%和84%),振幅因子与理论模型间的差异分别从(212%,155%,116%和080%)降到(031%,039%,034%和008%),同时还说明利用NAO99和ORI96全球海潮模型能获得比其他模型更佳的负荷改正效果.文章还利用国际地球动力学计划网络其他7个台站的超导重力仪观测研究了全球海潮模型的适定性问题,结果说明不同模型中不同潮波具有明显的区域特点,早期构制的SCW80全球海潮模型仍可作为大地测量研究中的重要参考模型.  相似文献   

13.
In this study, we collected 1 156 broadband vertical components records at 22 digital seismic stations in Xinjiang region, Ürümqi station, and 7 stations in the adjacent regions during the period of 1999–2003. The records were firstly processed by the stacked spectral ratio method to obtain Q 0 (Q at 1 Hz) and the frequency correlation factor η corresponding to each path. Based on the results, the distribution images of Q 0 and η in 1°×1° grids for Xinjiang region were gained by the back-projection technique. The results indicate that Q 0 is high (300–450) in the Tarim platform and marginal Siberian platform, while Q 0 is low (150–250) in the southern regions as west Kunlun fold system and Songpan-Ganzi fold system. In the northern regions as Junggar fold system and Tianshan fold system, Q 0 is also low (250–300) and η varies between 0.5 and 0.9.  相似文献   

14.
A three-dimensional regional ocean model is used to examine the impact of positive Indian ocean dipole (pIOD) events on the coastal upwelling features at the southwest coast of India (SWCI). Two model experiments are carried out with different surface boundary conditions that prevailed in the normal and pIOD years from 1982 to 2010. Model experiments demonstrate the weakening of coastal upwelling at the SWCI in the pIOD years. The reduced southward meridional wind stress off the SWCI leads to comparatively lower offshore Ekman transport during August–October in the pIOD years to that in normal years. The suppressed coastal upwelling results in warmer sea surface temperature and deeper thermocline in the pIOD years during June–September. The offshore spatial extent of upwelled colder (<?22 °C) water was up to 75.5° E in August–September in normal years that was limited up to 76.2° E in pIOD years. The heat budget analysis reveals the decreased contribution of vertical entrainment process to the mixed layer cooling in pIOD years which is almost half of that of normal years in October. The net heat flux term shows warming tendency during May–November with a higher magnitude (+?0.4 °C day?1) in normal years than pIOD years (+?0.28 °C day?1). The biological productivity is found to reduce during the pIOD years as the concentration of phytoplankton and zooplankton decreases over the region of coastal upwelling at SWCI. Nitrate concentration in the pIOD years dropped by half during August–September and dropped by an order of magnitude in October as compared to its ambient concentration of 13 μmol L?1 in normal years.  相似文献   

15.
Concentrations of lead were measured in a surface transect and at two vertical profile stations (15°N and 20°S) in the Central Pacific. These measurements complement similar measurements made earlier in the North Pacific at 33°N and in the Northwest Atlantic at 34°N [1,2], as well as recent measurements of eolian lead input fluxes near each of these locations [3]. The new transect of surface water concentrations of lead corroborates previous measurements, which decrease from 13 ng/kg at 30°N to 4 ng/kg at 17°S in the Central Pacific [4]. This transect gradient is shown to overlie a similar geographic gradient of subsurface maximum concentrations of lead in the three Pacific vertical profile stations, decreasing from 14 ng/kg at 33°N to 11 ng/kg at 14°N to 2.5 ng/kg at 20°S. Lead concentrations at each of those locations exhibit maxima at 400 m, decreasing concentrations to 2500 m and approximately concentrations of 0.8–1.1 ng/kg below that depth. The subsurface maximum at the northwest Atlantic profile station (36 ng/kg at 34°N) is also congruent with surface water lead concentrations which decrease from 806 ng/kg to 32 ng/kg in an offshore transect from Rhode Island to 34°N, 66°W [5], and the shape of the Atlantic profile is congruent with those in the Pacific. There is a positive correlation between the magnitudes of eolian lead input fluxes and the magnitudes of the upper water maxima in lead concentration profiles at corresponding locations as follows: South Pacific easterlies 3 ng/cm2 yr vs. 2.5 ng/kg; North Pacific easterlies 6 ng/cm2 yr vs. 11 ng/kg; North Pacific westerlies 50 ng/cm2 yr vs. 14 ng/kg; and North Atlantic westerlies 170 ng/cm2 yr vs. 36 ng/kg.This relationship enables one to view the anthropogenic perturbations of the marine lead cycle on a global scale, since the industrial origin of eolian and seawater lead has been established by correlations between geographic patterns of industrial lead emissions to the atmosphere and isotopic ratios of industrial leads [3] and by geographic patterns of Pb/silicate-dust ratios and lead isotopic ratios in ocean surface waters [3–5]. These new data coupled with earlier biogeochemical data indicate that surface water concentrations of lead in the North Pacific and North Atlantic are now conservatively estimated to be 8 to 20-fold greater and those in the South Pacific are 2-fold greater than natural concentrations because of industrial emissions of lead to the atmosphere.  相似文献   

16.
The increasing resolution of ground based gravity measurements (e.g. by superconducting gravimeters) as well as satellite based gravity field studies allows to study very small signals, globally as well as local. On the other hand, this requires the correction of such signals to uncover others. To study the Earth’s deep interior and the on-going dynamic processes requires the correction of disturbing signals, and one of these signals is related to ocean tidal loading. Although new ocean tide models are being derived from current satellite missions, there are still uncertainties.In this paper we present an intercomparison ocean tide models to test their fit to world-wide observations. Therefore, three TOPEX/POSEIDON (T/P) satellite derived models (CSR3.0, FES95.2 and TPXO.2) beside the classical SCHW80 model were selected for an accuracy assessment study. The selected models have been subjected to an intercomparison test, tide gauge validation test and comparison to 59 tidal gravity stations.The intercomparison test shows a good agreement between the T/P-based models for the open ocean and remarkable disagreement between the selected models in the coastal regions indicating that such models are still problematic in these regions. The tide gauge validation shows that the T/P derived models fit tide gauges better than SCHW80, with a better fit for the semidiurnal constituents than for the diurnal constituents. Comparing the gravimetric ocean-tide loading computed from the selected models with the residuals from a set of 59 tidal gravity stations shows that there is an improvement of the T/P derived models with respect to the Schwiderski model, especially in M2. However, this improvement is not as significant as the result of the comparison with the pelagic data. The procedure developed for the comparison of T/P derived models with SCHW80 is presented. The results provide not only information and improvement with regard to SCHW80, but also information about the properties of the new models. It is intended to continue this work applying the very recent models to see how they perform compared to this study.With this study we provide boundary conditions for the improvement of new ocean-tide models in order to benefit from the gravity measurements now possible regarding the evaluation of Earth structures and dynamic processes.  相似文献   

17.
In this paper a composite analysis was used to assess the influence of the North Atlantic Oscillation (NAO) on the winter daily rainfall and seasonal runoff at 28 stations of the Abruzzo region (Central Italy) during the period 1951–2012. Compositing was based on NAO? and NAO+ phases, identified by mean winter values of the normalized NAO index (NAOI) ≤?0.75 and ≥+0.75, respectively. In accordance with previous studies, it was found that NAO? phases determine, in general, a greater number of wet days (N w ) and (consequently) higher seasonal rainfall amounts in comparison to NAO+ phases. However, the NAO influence is characterized by a certain spatial variability, that can mostly be explained by orographic effects due to the Apennine Mountains. This is particularly evident for the mean rainfall depth per event (P e ) that, during NAO? phases, increases for the stations to the west of the Apennines, while it decreases for most of the stations to the east. The structure of winter daily rainfall of NAO+ and NAO? type, was described by a simple but effective first-order Markov process, determining the transition probabilities P01 (dry to wet) and P10 (wet to dry) and modelling the rainfall depth on wet days by a Weibull distribution. The most significant influence of NAO concerns P01 and the shape parameter of the Weibull distribution that are both higher during the NAO? phase. This means that NAO? phases are characterized by less persistent dry periods and less variable daily rainfall depths, in comparison to NAO+ phases. The effect of these differences on the winter seasonal runoff was explored by applying a Curve Number rainfall-runoff model. Significant increments of the mean seasonal runoff during NAO? phases were observed only for few stations (mainly on the west), characterized by corresponding increments of N w , P tot and P e .). NAO+ phases, instead, are characterized by relevant increments of the seasonal runoff variability, particularly on the eastern areas. In this context, the important regulating function of the watershed conditions was also discussed.  相似文献   

18.
《Journal of Geodynamics》2007,43(4-5):159-174
The increasing resolution of ground based gravity measurements (e.g. by superconducting gravimeters) as well as satellite based gravity field studies allows to study very small signals, globally as well as local. On the other hand, this requires the correction of such signals to uncover others. To study the Earth’s deep interior and the on-going dynamic processes requires the correction of disturbing signals, and one of these signals is related to ocean tidal loading. Although new ocean tide models are being derived from current satellite missions, there are still uncertainties.In this paper we present an intercomparison ocean tide models to test their fit to world-wide observations. Therefore, three TOPEX/POSEIDON (T/P) satellite derived models (CSR3.0, FES95.2 and TPXO.2) beside the classical SCHW80 model were selected for an accuracy assessment study. The selected models have been subjected to an intercomparison test, tide gauge validation test and comparison to 59 tidal gravity stations.The intercomparison test shows a good agreement between the T/P-based models for the open ocean and remarkable disagreement between the selected models in the coastal regions indicating that such models are still problematic in these regions. The tide gauge validation shows that the T/P derived models fit tide gauges better than SCHW80, with a better fit for the semidiurnal constituents than for the diurnal constituents. Comparing the gravimetric ocean-tide loading computed from the selected models with the residuals from a set of 59 tidal gravity stations shows that there is an improvement of the T/P derived models with respect to the Schwiderski model, especially in M2. However, this improvement is not as significant as the result of the comparison with the pelagic data. The procedure developed for the comparison of T/P derived models with SCHW80 is presented. The results provide not only information and improvement with regard to SCHW80, but also information about the properties of the new models. It is intended to continue this work applying the very recent models to see how they perform compared to this study.With this study we provide boundary conditions for the improvement of new ocean-tide models in order to benefit from the gravity measurements now possible regarding the evaluation of Earth structures and dynamic processes.  相似文献   

19.
Two cases when Pc5 geomagnetic pulsations were registered at the IMAGE Scandinavian network of stations and with STARE radars in the afternoon sector (1700–1800 MLT) during the recovery phase of the moderate magnetic storm are analyzed in detail. Using the ground-based observations, it has been indicated that classical quasimonochromatic resonance Pc5 pulsations were observed in the first case (on October 12, 1999; Kp = 5); in this case the maximal amplitude of the spectral maximum at a frequency of 2.5 mHz was registered at Φ ~ 65°. Two maximums were observed in the spectrum in the second case (on October 13, 1999; Kp = 4): ~2.5 mHz (the same maximum) and 2.9 mHz; in this case the maximal oscillation amplitude (2.5 mHz) shifted to Φ > 67°. These results were compared with the echo signal intensity simultaneously registered with the STARE Finland radar on a beam oriented along the 105° geomagnetic meridian. The spatial-temporal maps of the Pc5 pulsation amplitude latitudinal distribution (“keograms”), constructed based on the radar measurements in the wide range of geomagnetic latitudes (63°–70°) where the resolution was substantially higher than that of the ground-based observations, made it possible to detect two regions spaced in latitude (Φ ~ 65° and Φ ~ 67°–68°) with the simultaneous excitation of oscillations (double resonance?), between which the plasmapause projection was supposedly located.  相似文献   

20.
In this study, pre-seismic and post-seismic total electron content (TEC) anomalies of 63 Mw?≥?5.0 earthquakes in Turkey (36°–42°N, 26°–45°E) were statistically investigated. The largest earthquake that occurred in Turkey during 2003–2016 is the Mw 7.1 Van earthquake on October 23, 2011. The TEC data of epicenters is obtained from CODE-GIM using a simple 4-point bivariate interpolation. The anomalies of TEC variations were determined by using a quartile-based running median process. In order to validate GIM results, we used the GPS-TEC data of available four IGS stations within the size of the Van earthquake preparation area. The anomalies that are detected by GIM and GPS-TEC show a similar pattern. Accordingly, the results obtained with CODE-GIM are reliable. The statistical results show that there are not prominent earthquake precursors for Mw?≤?6.0 earthquakes in Turkey.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号