首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
The span of coordinate time series affects the determination of an optimal noise model. We analyzed position data recorded for 10 continuous Global Positioning System(GPS) sites from 1998.0 to mid-2009 on the Australian Plate to estimate the best noise model and thereafter obtain the true uncertainties of the velocity, employing the maximum likelihood estimation(MLE) method. MLE was employed to analyze the data in four ways. In the first two analyses, the noise was assumed to be a combination of flicker noise and white noise for the raw time series and spatially filtered time series. In the final two analyses, the spectral indices and amplitudes were simultaneously estimated for a power law noise plus white noise model for the raw time series and spatially filtered time series. We conclude that the noise model of GPS time series in Australia can be best described as the combination of flicker noise and white noise. Velocity uncertainties fall below ~0.2 mm/yr when the time span exceeds ~9.5 years. A comparison of noise amplitudes and maximum likelihood estimation values between the raw and spatially filtered time series suggests that traditional spatial filtering to remove common-mode errors might not be applicable to the raw time series of this region.  相似文献   

2.
香港GPS基准站坐标序列特征分析   总被引:27,自引:3,他引:24       下载免费PDF全文
利用香港GPS连续运行参考站网络2001年1月至2007年8月的观测资料,全面深入地分析了12个基准站坐标序列特征.本文采用主成分空间滤波算法去除公共误差,来提高坐标序列的信噪比,并采用最大似然估计准则定量估计滤波后坐标序列的噪声特性,计算了地球表面质量负荷(包括大气、非潮汐海洋、积雪和土壤水)对香港GPS基准站坐标序列的影响.研究结果表明:香港GPS基准站坐标序列具有高度的空间相关性,其公共误差具有较强的季节性变化特征;地表质量负荷变化引起的香港地壳形变可以解释公共误差序列中约为3mm的垂向周年变化,经过质量负荷改正后的公共误差序列与高阶电离层误差高度相关;滤波后坐标序列的噪声特性可以用可变白噪声加闪烁噪声模型来描述,顾及闪烁噪声所计算的速度误差要比只考虑可变白噪声计算的速度误差大2~6倍;基准站间存在达1.5 mm/yr的相对水平运动,揭示香港地区存在活动断层;部分基准站坐标具有明显的振幅为1~2 mm本地季节性变化,所有测站的残差序列也表现出强烈的季节性变化.  相似文献   

3.
Most GPS time-series exhibit a seasonal signal that can have an amplitude of a few millimetres. This seasonal signal can be removed by fitting an extra sinusoidal signal with a period of one year to the GPS data during the estimation of the linear trend.However, Blewitt and Lavallée (2002) showed that including an annual signal in the estimation process still can give a larger linear trend error than the trend error estimated from data from which the annual signal has been removed by other means. They assumed that the GPS data only contained white noise and we extend their result to the case of power-law plus white noise which is known to exist in most GPS observations. For the GPS stations CASC, LAGO, PDEL and TETN the difference in trend error between having or not having an annual signal in the data is around ten times larger when a power-law plus white noise model is used instead of a pure white noise model. Next, our methodology can be used to estimate for any station how much the accuracy of the linear trend will improve when one tries to subtract the annual signal from the GPS time-series by using a physical model.Finally, we demonstrate that for short time-series the trend error is more influenced by the fact that the noise properties also need to be estimated from the data. This causes on average an underestimation of the trend error.  相似文献   

4.
Time–frequency characterization is useful in understanding the nonlinear and non-stationary signals of the hydro-climatic time series. The traditional Fourier transform, and wavelet transform approaches have certain limitations in analyzing non-linear and non-stationary hydro-climatic series. This paper presents an effective approach based on the Hilbert–Huang transform to investigate time–frequency characteristics, and the changing patterns of sub-divisional rainfall series in India, and explored the possible association of monsoon seasonal rainfall with different global climate oscillations. The proposed approach integrates the complete ensemble empirical mode decomposition with adaptive noise algorithm and normalized Hilbert transform method for analyzing the spectral characteristics of two principal seasonal rainfall series over four meteorological subdivisions namely Assam-Meghalaya, Kerala, Orissa and Telangana subdivisions in India. The Hilbert spectral analysis revealed the dynamic nature of dominant time scales for two principal seasonal rainfall time series. From the trend analysis of instantaneous amplitudes of multiscale components called intrinsic mode functions (IMFs), it is found that both intra and inter decadal modes are responsible for the changes in seasonal rainfall series of different subdivisions and significant changes are noticed in the amplitudes of inter decadal modes of two seasonal rainfalls in the four subdivisions since 1970s. Further, the study investigated the links between monsoon rainfall with the global climate oscillations such as Quasi Bienniel Oscillation (QBO), El Nino Southern Oscillation (ENSO), Sunspot Number (SN), Atlantic Multidecadal Oscillation (AMO) etc. The study noticed that the multiscale components of rainfall series IMF1, IMF2, IMF3, IMF4 and IMF5 have similar periodic structure of QBO, ENSO, SN, tidal forcing and AMO respectively. As per the seasonal rainfall patterns is concerned, the results of the study indicated that for Assam-Meghalaya subdivision, there is a likelihood of extreme rare events at ~0.2 cycles per year, and both monsoon and pre-monsoon rainfall series have decreasing trends; for Kerala subdivision, extreme events can be expected during monsoon season with shorter periodicity (~2.5 years), and monsoon rainfall has statistically significant decreasing trend and post-monsoon rainfall has a statistically significant increasing trend; and for Orissa subdivision, there are chances of extremes rainfall events in monsoon season and a relatively stable rainfall pattern during post-monsoon period, but both monsoon and post-monsoon rainfall series showed an overall decreasing trend; for Telangana subdivision, there is a likelihood of extreme events during monsoon season with a periodicity of ~4 years, but both monsoon and post-monsoon rainfall series showed increasing trends. The results of correlation analysis of IMF components of monsoon rainfall and five climate indices indicated that the association is expressed well only for low frequency modes with similar evolution of trend components.  相似文献   

5.
v--vThe prototype International Data Centre (IDC) in Arlington, Virginia has been acquiring data from seismic stations at locations designated in the Comprehensive Test-Ban Treaty for the International Monitoring System (IMS) since the start of 1995. A key characteristic of these stations is their background noise levels and their seasonal and diurnal variability. Since June 1997 an automated sample selection effort has collected over 700,000 individual noise sample spectra from 39 primary and 57 auxiliary stations. Monthly median and 5 and 95 percentile estimates have been calculated for each channel of every station. Compatibility of median spectra obtained for the same station and channel in the same month for two different years confirms the consistency of the noise-sampling algorithm used. A preliminary analysis of the results shows strong (more than a factor of two) seasonal variation at a quarter of all stations. Strong diurnal variations at half of the sites indicate that many of the selected sites are poorly located with respect to cultural noise sources. The results of this study are already being used to evaluate station quality, improve those processes that require background noise values, such as automatic association and requesting auxiliary station data, and to improve the estimation of station and network detection and location thresholds.  相似文献   

6.
Water resources and soil erosion are the most important environmental concerns in the Yangtze River basin, where soil erosion and sediment yield are closely related to rainfall erosivity. The present study explores the spatial and temporal changing patterns of the rainfall erosivity in the Yangtze River basin of China during 1960–2005 at annual, seasonal and monthly scales. The Mann–Kendall test is employed to detect the trends during 1960–2005, and the T test is applied to investigate possible changes between 1991–2005 and 1960–1990. Meanwhile the Rescaled Range Analysis is used for exploring future trend of rainfall erosivity. Moreover the continuous wavelet transform technique is using studying the periodicity of the rainfall erosivity. The results show that: (1) The Yangtze River basin is an area characterized by uneven spatial distribution of rainfall erosivity in China, with the annual average rainfall erosivity range from 131.21 to 16842 MJ mm ha?1 h?1. (2) Although the directions of trends in annual rainfall erosivity at most stations are upward, only 22 stations have significant trends at the 90 % confidence level, and these stations are mainly located in the Jinshajiang River basin and Boyang Lake basin. Winter and summer are the seasons showing strong upward trends. For the monthly series, significant increasing trends are mainly found during January, June and July. (3) Generally speaking, the results detected by the T test are quite consistent with those detected by the Mann–Kendall test. (4) The rainfall erosivity of Yangtze River basin during winter and summer will maintain a detected significant increasing trend in the near future, which may bring greater risks to soil erosion. (5) The annual and seasonal erosivity of Yangtze River basin all have one significant periodicity of 2–4 years.  相似文献   

7.
Trend identification is a substantial issue in hydrologic series analysis, but it is also a difficult task in practice due to the confusing concept of trend and disadvantages of methods. In this article, an improved definition of trend was given as follows: ‘a trend is the deterministic component in the analysed data and corresponds to the biggest temporal scale on the condition of giving the concerned temporal scale’. It emphasizes the intrinsic and deterministic properties of trend, can clearly distinguish trend from periodicities and points out the prerequisite of the concerned temporal scale only by giving which the trend has its specific meaning. Correspondingly, the discrete wavelet‐based method for trend identification was improved. Differing from those methods used presently, the improved method is to identify trend by comparing the energy difference between hydrologic data and noise, and it can simultaneously separate periodicities and noise. Furthermore, the improved method can quantitatively estimate the statistical significance of the identified trend by using proper confidence interval. Analyses of both synthetic and observed series indicated the identical power of the improved method as the Mann–Kendall test in assessing the statistical significance of the trend in hydrologic data, and by using the former, the identified trend can adaptively reflect the nonlinear and nonstationary variability of hydrologic data. Besides, the results also showed the influences of three key factors (wavelet choice, decomposition level choice and noise content) on discrete wavelet‐based trend identification; hence, they should be carefully considered in practice. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

8.
利用17年的SABER(Sounding of the Atmosphere using Broadband Emission Radiometry)Level2C数据研究了中间层与低热层大气(MLT, Mesosphere and Lower Thermosphere) CO2 VMR(Volume Mixing Ratio)的年际变化特征.使用多元线性回归模型对双月平均时间序列拟合,定量地提取各变化特征.结果表明,SABER CO2 VMR长期趋势在中间层保持在5.5%/decade左右,在中间层顶和低热层降低至4.5%/decade左右;结果与模式预测在统计意义上相符.长期趋势没有显著的纬度差异,但在各纬度上都具有明显的季节依赖,MLT CO2 VMR长期趋势的季节性改变源自低层大气长期趋势季节性改变.SABER CO2 VMR对QBO (Quasi-Biannual Oscillation)和ENSO (El Nino-Southern Oscillation)在绝大多数区域没有统计显著的响应;对太阳活动11年循环以负响应为主,在部分区域出现的微弱正响应目前没有合适的物理机制解释.  相似文献   

9.
Series of gravity recordings at the stations Medicina (Italy) and Wettzell (Germany) are investigated to separate seasonal gravity variations from long-term trends in gravity. The findings are compared to height variations monitored by continuous GPS observations. To study the origin of these variations in height and gravity the environmental parameters at the stations are included in the fact finding. In Medicina, a clear seasonal signal is visible in the gravity and height data series, caused by seasonal fluctuations in the atmosphere including mass redistribution, the ocean, groundwater but also by geo-mechanical effects such as soil consolidation and thermal expansion of the structure supporting the GPS antenna. In Wettzell, no seasonal effect could be clearly identified, and the long-term trend in gravity is mainly caused by ground water variations. The successful combination of height and gravity series with the derived ratio of gravity to height changes indicates that the long-term trends in height and gravity are most likely due to mass changes rather than to tectonic movements.  相似文献   

10.
In this paper,we analyze the time series of site coordinates of 27 continuously monitoring GPS sites covered bythe Crustal Movement Observation Network of China over the whole country.The data are obtained in the periodfrom the beginning of the observation to the November of 2005.On the basis of data processing,we analyze thepower spectrum density of coordinate component noise at each site and calculate the spectral indexes manifestingthe noise property of each component.The spectral indexes indicate that for most sites,the noise of time series ofeach coordinate component can be addressed by the model of white noise flicker noise;and for a small amountof sites,it can be described by the model of white noise flicker noise random walk noise.We also quantita-tively estimate each noise component in the model by using the criterion of maximum likelihood estimation.Theresult shows that the white noise in the time series of GPS site coordinates does not constitute the main part ofnoise.Therefore,the error estimation of site movement parameters is usually too small,or too optimistic if weconsider the white noise only.Correspondingly,if this factor is not fully considered in explaining these movementparameters,it might mislead the readers.  相似文献   

11.
连续观测站的噪声分析   总被引:18,自引:0,他引:18       下载免费PDF全文
黄立人  符养 《地震学报》2007,29(2):197-202
分析了中国地壳运动观测网络在全国布设的27个GPS连续站开始运行以来至2005年11月的站坐标时间序列. 在对数据进行清理的基础上, 分析了各站坐标分量噪声的功率谱密度, 计算了表征各分量的噪声特性的谱指数. 谱指数显示, 大部分站的各坐标分量时间序列的噪声可以用白噪声+闪烁噪声的模型来描述, 少部分则可用白噪声+闪烁噪声+随机漫步噪声的模型来描述. 采用最大似然估计准则, 定量估计了模型中的各噪声分量. 结果表明, GPS站坐标时间序列中白噪声甚至不是噪声的主要成分. 因此, 仅顾及白噪声得到的测站运动参数的误差估计, 实际上是过高的, 或者说是过于乐观的. 相应地, 对这些运动参数作出解释时, 如果不充分考虑这一因素, 有可能会对读者产生误导.   相似文献   

12.
Temporal and spatial variations of stable oxygen (18O) and hydrogen (2H) isotope measurements in precipitation act as important proxies for changing hydro‐meteorological and regional and global climate patterns. Temporal trends in time series of the stable isotope composition in precipitation were rarely observed, and they are poorly understood. These might be a result of a lack of proper trend detection tools and effort for exploring trend processes. Here, we investigate temporal trends of δ18O in precipitation at 17 observation stations in Germany between 1978 and 2009. We test if significant trends in the isotope time series from different models can be observed. Mann–Kendall trend tests are applied on the isotope series, using general multiplicative seasonal autoregressive integrate moving average (ARIMA) models, which account for first and higher order serial correlations. Effects of temperature, precipitation, and geographic parameters on isotope trends are also investigated in the proposed models. To benchmark our proposed approach, the ARIMA results are compared with a trend‐free pre‐whitening procedure, the state of the art method for removing the first order autocorrelation in environmental trend studies. Moreover, we further explore whether higher order serial correlations in isotope series affects our trend results. Overall, three out of the 17 stations show significant changes when higher order autocorrelation are adjusted, and four show a significant trend when temperature and precipitation effects are considered. The significant trends in the isotope time series generally occur only at low elevation stations. Higher order autoregressive processes are shown to be important in the isotope time series analysis. Results suggest that the widely used trend analysis with only the first order autocorrelation adjustment may not adequately take account of the high order autocorrelated processes in the stable isotope series. The investigated time series analysis method including higher autocorrelation and external climate variable adjustments is shown to be a better alternative. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

13.
Ezer  Tal 《Ocean Dynamics》2022,72(11):741-759

The long-term variability of sea level and surface flows in the Gulf of Mexico (GOM) is studied using global monthly sea level reconstruction (RecSL) for 1900–2015. The study explored the long-term relation between the dynamics of the GOM and inflows/outflows through the Yucatan Channel (YC) and the Florida Straits (FS). The results show a century-long trend of increased mean velocity and variability in the Loop Current (LC); however, no significant upward trend was found in the YC and FS flows, only increased variability. Empirical orthogonal function (EOF) analysis of sea surface height found spatial patterns dominated by variations in the LC and temporal variations on time scales ranging from a few months to multidecadal. The time evolution of each EOF mode of sea level is correlated with the velocity of either the LC, the YC, or the FS or some combination of the different flows. The mean sea level difference between the GOM and the northwestern Caribbean Sea was found to be influenced by the North Atlantic Oscillation (NAO), with unusually high differences during the 1970s when the NAO index was low and the Atlantic Ocean circulation was weak. Extreme peaks in SL difference coincide with the extension of the LC and the seasonal eddy shedding pattern. The observed seasonal cycle in the extension area of the LC as obtained from 20 years of altimeter data is significantly correlated (R = 0.63; confidence level = 98%) with the seasonal YC flow obtained from 116 years of the RecSL data. However, the same LC extension record had lower correlation (R = 0.45; confidence level = 90%) with the observed YC transport obtained from direct moored measurements over ~ 5 years, indicating the need for much longer measurements, since the LC extension and the YC flow are strongly affected by interannual and decadal variations. The study demonstrates the usefulness of even a coarse-resolution reconstruction for studies of regional ocean variability and climate change over longer time scales than current direct observations allow.

  相似文献   

14.
This paper addresses the role of meteorological forcing on mean sea level (MSL) variability at the tide gauge of Cuxhaven over a period from 1871 to 2008. It is found that seasonal sea level differs significantly from annual means in both variability and trends. The causes for the observed differences are investigated by comparing to changes in wind stress, sea level pressure and precipitation. Stepwise regression is used to estimate the contribution of the different forcing factors to sea level variability. The model validation and sensitivity analyses showed that a robust and timely independent estimation of regression coefficients becomes possible if at least 60 to 80 years of data are available. Depending on the season, the models are able to explain between 54 % (spring, April to June) and 90 % (winter, January to March) of the observed variability. Most parts of the observed variability are attributed to changes in zonal wind stress, whereby the contribution of sea level pressure, precipitation and meridional wind stress is rather small but still significant. On decadal timescales, the explanatory power of local meteorological forcing is considerable weaker, suggesting that the remaining variability is attributed to remote forcing over the North Atlantic. Although meteorological forcing contributes to linear trends in some sub-periods of seasonal time series, the annual long-term trend is less affected. However, the uncertainties of trend estimation can be considerably reduced, when removing the meteorological influences. A standard error smaller than 0.5 mm/year requires 55 years of data when using observed MSL at Cuxhaven tide gauge. In contrast, a similar standard error in the meteorologically corrected residuals is reached after 32 years.  相似文献   

15.
Trends in the timing and magnitude of floods in Canada   总被引:2,自引:0,他引:2  
This study investigates trends in the timing and magnitude of seasonal maximum flood events across Canada. A new methodology for analyzing trends in the timing of flood events is developed that takes into account the directional character and multi-modality of flood occurrences. The methodology transforms the directional series of flood occurrences into new series by defining a new location of the origin. A test of flood seasonality (multi-modality) is then applied to identify dominant flood seasons. Floods from the dominant seasons are analyzed separately by a seasonal trend analysis. The Mann–Kendall test in conjunction with the method of pre-whitening is used in the trend analysis. Over 160 streamflow records from one common observation period are analyzed in watersheds with relatively pristine and stable land-use conditions. The results show weak signals of climate variability and/or change present in the timing of floods in Canada during the last three decades. Most of the significant trends in the timing of spring snowmelt floods are negative trends (earlier flood occurrence) found in the southern part of Canada. There are no significant trends identified in the timing of fall rainfall floods. However, the significance of the fall, rainfall-dominated flood season has been increasing in several analyzed watersheds. This may indicate increasing intensity of rainfall events during the recent years. Trends in the magnitude of floods are more pronounced than the trends in the timing of floods. Almost one fifth of all the analyzed stations show significant trends in the magnitude of snowmelt floods. Most of the significant trends are negative trends, suggesting decreasing magnitudes of snowmelt floods in Canada over the last three decades. Significant negative trends are found particularly in southern Ontario, northern Saskatchewan, Alberta and British Columbia. There are no significant trends in the magnitude of rainfall floods found in the analyzed streamflow records. The results support the outcomes of previous streamflow trend studies conducted in Canada.  相似文献   

16.
张风霜  畅柳 《地震》2017,37(1):92-102
利用2010—2015年陆态网络GPS连续观测站的三维站坐标时序资料, 在采用最大似然估计MLE方法计算得到云南地区GPS基准站三维站坐标分量最优噪声模型的基础上, 进一步定量对比分析了有色噪声和白噪声模型下陆态网络GPS连续站速度场稳定性与时间序列长度的关系。 结果表明: 在闪烁噪声模型下, 中国大陆大部分地区GPS连续站1 a时段的速度不确定度大。 而随着时间的逐年累积, 闪烁噪声模型下的速度不确定度与同时段白噪声模型下的差异逐步减小, 闪烁噪声模型下3~4 a时段的速度不确定度与同时段白噪声模型的差异小于0.2 mm。 研究说明在有色噪声模型下要获得高精度的速度场需要3~4 a的观测数据, 白噪声模型下GPS速度场精度可靠可信的前提是采用不少于3~4 a的数据。  相似文献   

17.
Based on the data of the rocket sounding of the middle atmosphere (25–75 km) in the Northern and Southern hemispheres during 1964–1994, the average monthly vertical distributions of the long-term temperature trend have been obtained for the low, middle, and high latitudes. These distributions have been approximated by the series of the harmonic functions describing the latitudinal, seasonal, and vertical variations in the temperature trend. The characteristics of the obtained harmonics are presented.  相似文献   

18.
This study investigated using Monte Carlo simulation the interaction between a linear trend and a lag‐one autoregressive (AR(1)) process when both exist in a time series. Simulation experiments demonstrated that the existence of serial correlation alters the variance of the estimate of the Mann–Kendall (MK) statistic; and the presence of a trend alters the estimate of the magnitude of serial correlation. Furthermore, it was shown that removal of a positive serial correlation component from time series by pre‐whitening resulted in a reduction in the magnitude of the existing trend; and the removal of a trend component from a time series as a first step prior to pre‐whitening eliminates the influence of the trend on the serial correlation and does not seriously affect the estimate of the true AR(1). These results indicate that the commonly used pre‐whitening procedure for eliminating the effect of serial correlation on the MK test leads to potentially inaccurate assessments of the significance of a trend; and certain procedures will be more appropriate for eliminating the impact of serial correlation on the MK test. In essence, it was advocated that a trend first be removed in a series prior to ascertaining the magnitude of serial correlation. This alternative approach and the previously existing approaches were employed to assess the significance of a trend in serially correlated annual mean and annual minimum streamflow data of some pristine river basins in Ontario, Canada. Results indicate that, with the previously existing procedures, researchers and practitioners may have incorrectly identified the possibility of significant trends. Copyright © Environment Canada. Published by John Wiley & Sons, Ltd.  相似文献   

19.
连续GPS观测中的相关噪声分析   总被引:1,自引:0,他引:1       下载免费PDF全文
田云锋  沈正康  李鹏 《地震学报》2010,32(6):696-704
分析了中国地壳运动观测网络GPS基准站位置时间序列中相关噪声的特性.利用最大似然估计反演了幂指数噪声的谱指数;考察了多种噪声模型,发现除闪烁噪声之外,约50%台站的东分量中的相关噪声可用"随机漫步噪声+一阶高斯-马尔可夫噪声"模型来描述,与其它分析机构产出的全球解并不一致,显示出处理策略的不同能够造成相关噪声的差异.白噪声和闪烁噪声的大小存在明显的地理分异,东南沿海台站的噪声要大于西北内陆台站.考虑相关噪声时,各站的速率估计方差要增大一个数量级以上,水平分量多在1mm/a以内,而垂向分量则较大.  相似文献   

20.
The relationship between air (Ta) and water temperature (Tw) is very important because it shows how the temperature of a water body might respond to future changes in surface Ta. Mean monthly Tw records of three gauging stations (Bezdan, Bogojevo i Veliko Gradi?te) were analysed alongside mean monthly discharge (Q) for the same stations. Additionally, Ta series from two meteorological stations (Sombor and Veliko Gradi?te) were correlated with Tw variations over the period 1950–2012. Locally weighted scatter point smoothing (LOWESS) was used to investigate long‐term trends in the raw data, alongside the Mann–Kendall (MK) trend test. Trend significance was established using Yue–Pilon's pre‐whitening approaches to determine trends in climate data. Also, the rescaled adjusted partial sums (RAPS) method was used to detect dates of possible changes in the time series. Statistically significant warming trends were observed for annual and seasonal minimum and maximum Tw at all investigated sites. The strongest warming was observed at Bogojevo gauging station for seasonal maximum Tw, with +0.05 °C per year on average. RAPS established that the trend began in the 1980s. This behaviour is linked to climate patterns in the North and East Atlantic which determine the amount of heat advected onto mainland Europe. Statistically significant correlations were found for all Tw on an annual basis. Overall, the strongest correlations (p < 0.01) between Tw residuals and the North Atlantic Oscillation (NAO) were recorded for the winter period. These findings suggest possible predictability of Tw over seasonal time‐scales. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

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

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