首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
Time-variable gravity data of the GRACE (Gravity Recovery And Climate Experiment) satellite mission provide global information on temporal variations of continental water storage. In this study, we incorporate GRACE data for the first time directly into the tuning process of a global hydrological model to improve simulations of the continental water cycle. For the WaterGAP Global Hydrology Model (WGHM), we adopt a multi-objective calibration framework to constrain model predictions by both measured river discharge and water storage variations from GRACE and illustrate it on the example of three large river basins: Amazon, Mississippi and Congo. The approach leads to improved simulation results with regard to both objectives. In case of monthly total water storage variations we obtained a RMSE reduction of about 25 mm for the Amazon, 6 mm for the Mississippi and 1 mm for the Congo river basin. The results highlight the valuable nature of GRACE data when merged into large-scale hydrological modeling. Furthermore, they reveal the utility of the multi-objective calibration framework for the integration of remote sensing data into hydrological models.  相似文献   

2.
Reducing aliasing effects of insufficiently modelled high-frequent, non-tidal mass variations of the atmosphere, the oceans and the hydrosphere in gravity field models derived from the Gravity Recovery and Climate Experiment (GRACE) satellite mission is the topic of this study. The signal content of the daily GRACE gravity field model series (ITG-Kalman) is compared to high-frequency bottom pressure variability and terrestrially stored water variations obtained from recent numerical simulations from an ocean circulation model (OMCT) and two hydrological models (WaterGAP Global Hydrology Model, Land Surface Discharge Model). Our results show that daily estimates of ocean bottom pressure from the most recent OMCT simulations and the daily ITG-Kalman solutions are able to explain up to 40 % of extra-tropical sea-level variability in the Southern Ocean. In contrast to this, the daily ITG-Kalman series and simulated continental total water storage variability largely disagree at periods below 30 days. Therefore, as long as no adequate hydrological model will become available, the daily ITG-Kalman series can be regarded as a good initial proxy for high-frequency mass variations at a global scale. As a second result of this study, based on monthly solutions as well as daily observation residuals, it is shown that applying this GRACE-derived de-aliasing model supports the determination of the time-variable gravity field from GRACE data and the subsequent geophysical interpretation. This leads us to the recommendation that future satellite concepts for determining mass variations in the Earth system should be capable of observing higher frequeny signals with sufficient spatial resolution.  相似文献   

3.
Better quantification of continental water storage variations is expected to improve our understanding of water flows, including evapotranspiration, runoff and river discharge as well as human water abstractions. For the first time, total water storage (TWS) on the land area of the globe as computed by the global water model WaterGAP (Water Global Assessment and Prognosis) was compared to both gravity recovery and climate experiment (GRACE) and global positioning system (GPS) observations. The GRACE satellites sense the effect of TWS on the dynamic gravity field of the Earth. GPS reference points are displaced due to crustal deformation caused by time-varying TWS. Unfortunately, the worldwide coverage of the GPS tracking network is irregular, while GRACE provides global coverage albeit with low spatial resolution. Detrended TWS time series were analyzed by determining scaling factors for mean annual amplitude (f GRACE) and time series of monthly TWS (f GPS). Both GRACE and GPS indicate that WaterGAP underestimates seasonal variations of TWS on most of the land area of the globe. In addition, seasonal maximum TWS occurs 1 month earlier according to WaterGAP than according to GRACE on most land areas. While WaterGAP TWS is sensitive to the applied climate input data, none of the two data sets result in a clearly better fit to the observations. Due to the low number of GPS sites, GPS observations are less useful for validating global hydrological models than GRACE observations, but they serve to support the validity of GRACE TWS as observational target for hydrological modeling. For unknown reasons, WaterGAP appears to fit better to GPS than to GRACE. Both GPS and GRACE data, however, are rather uncertain due to a number of reasons, in particular in dry regions. It is not possible to benefit from either GPS or GRACE observations to monitor and quantify human water abstractions if only detrended (seasonal) TWS variations are considered. Regarding GRACE, this is mainly caused by the attenuation of the TWS differences between water abstraction variants due to the filtering required for GRACE TWS. Regarding GPS, station density is too low. Only if water abstractions lead to long-term changes in TWS by depletion or restoration of water storage in groundwater or large surface water bodies, GRACE may be used to support the quantification of human water abstractions.  相似文献   

4.
The hydrological contribution to polar motion is a major challenge in explaining the observed geodetic residual of non-atmospheric and non-oceanic excitations since hydrological models have limited input of comprehensive global direct observations. Although global terrestrial water storage (TWS) estimated from the Gravity Recovery and Climate Experiment (GRACE) provides a new opportunity to study the hydrological excitation of polar motion, the GRACE gridded data are subject to the post-processing de-striping algorithm, spatial gridded mapping and filter smoothing effects as well as aliasing errors. In this paper, the hydrological contributions to polar motion are investigated and evaluated at seasonal and intra-seasonal time scales using the recovered degree-2 harmonic coefficients from all GRACE spherical harmonic coefficients and hydrological models data with the same filter smoothing and recovering methods, including the Global Land Data Assimilation Systems (GLDAS) model, Climate Prediction Center (CPC) model, the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis products and European Center for Medium-Range Weather Forecasts (ECMWF) operational model (opECMWF). It is shown that GRACE is better in explaining the geodetic residual of non-atmospheric and non-oceanic polar motion excitations at the annual period, while the models give worse estimates with a larger phase shift or amplitude bias. At the semi-annual period, the GRACE estimates are also generally closer to the geodetic residual, but with some biases in phase or amplitude due mainly to some aliasing errors at near semi-annual period from geophysical models. For periods less than 1-year, the hydrological models and GRACE are generally worse in explaining the intraseasonal polar motion excitations.  相似文献   

5.
Since its launch in April 2002, the Gravity Recovery and Climate Experiment (GRACE) mission is recording the Earth’s time-variable gravity field with temporal and spatial resolutions of typically 7–30?days and a few hundreds of kilometers, allowing the monitoring of continental water storage variations from both continental and river-basin scales. We investigate here large scale hydrological variations in Africa using different GRACE spherical harmonic solutions, using different processing strategies (constrained and unconstrained solutions). We compare our GRACE estimates to different global hydrology models, with different land-surface schemes and also precipitation forcing. We validate GRACE observations through two different techniques: first by studying desert areas, providing an estimate of the precision. Then we compare GRACE recovered mass variations of main lakes to volume changes derived from radar altimetry measurements. We also study the differences between different publicly available precipitation datasets from both space measurements and ground rain gauges, and their impact on soil-moisture estimates.  相似文献   

6.
The increased availability of global datasets and technologies such as global hydrologic models and the Gravity Recovery and Climate Experiment (GRACE) satellites have resulted in a growing number of global‐scale assessments of water availability using simple indices of water stress. Developed initially for surface water, such indices are increasingly used to evaluate global groundwater resources. We compare indices of groundwater development stress for three major agricultural areas of the United States to information available from regional water budgets developed from detailed groundwater modeling. These comparisons illustrate the potential value of regional‐scale analyses to supplement global hydrological models and GRACE analyses of groundwater depletion. Regional‐scale analyses allow assessments of water stress that better account for scale effects, the dynamics of groundwater flow systems, the complexities of irrigated agricultural systems, and the laws, regulations, engineering, and socioeconomic factors that govern groundwater use. Strategic use of regional‐scale models with global‐scale analyses would greatly enhance knowledge of the global groundwater depletion problem.  相似文献   

7.
本文利用CSR发布的GRACE RL06时变重力场模型,结合两种水文模式、卫星测高、降雨和蒸散等多源数据,从多个角度综合系统地分析维多利亚湖流域2003-01-2017-06的陆地水储量变化.比较了正向建模方法和单一尺度因子对泄漏误差的改正效果,经对比采用正向建模方法在此流域效果更好.基于多源数据得出以下三点与此前研究...  相似文献   

8.
Seasonal and interannual changes in the Earth's gravity field are mainly due to mass exchange among the atmosphere,ocean,and continental water sources.The terrestrial water storage changes,detected as gravity changes by the Gravity Recovery and Climate Experiment(GRACE) satellites,are mainly caused by precipitation,evapotranspiration,river transportation and downward infiltration processes.In this study,a land data assimilation system LDAS-G was developed to assimilate the GRACE terrestrial water storage(TWS) data into the Community Land Model(CLM3.5) using the POD-based ensemble four-dimensional variational assimilation method PODEn4 DVar,disaggregating the GRACE large-scale terrestrial water storage changes vertically and in time,and placing constraints on the simulation of vertical hydrological variables to improve land surface hydrological simulations.The ideal experiments conducted at a single point and assimilation experiments carried out over China by the LDAS-G data assimilation system showed that the system developed in this study improved the simulation of land surface hydrological variables,indicating the potential of GRACE data assimilation in large-scale land surface hydrological research and applications.  相似文献   

9.
Time variable gravity field models derived from the satellite mission GRACE have been demonstrated to be consistent with water mass variations in the global hydrological cycle. Independent observations are provided by terrestrial measurements. In order to achieve a maximum of reliability and information gain, ground-based gravity observations may be deployed for comparison with the gravity field variations derived from the GRACE satellite mission. In this context, the data of the network of superconducting gravimeters (SG) of the ‘Global Geodynamics Project’ (GGP) are of particular interest. This study is focused on the dense SG network in Central Europe with its long-term gravity observations. It is shown that after the separation and reduction of local hydrological effects in the SG observations especially for subsurface stations, the time-variable gravity signals from GRACE agree well with the terrestrial observations from the SG station cluster.Station stability of the SG sites with respect to vertical deformations was checked by GNSS based observations. Most of the variability can be explained by loading effects due to changes in continental water storage, and, in general, the stability of all stations has been confirmed.From comparisons based on correlation and coherence analyses in combination with the root mean square (RMS) variability of the time series emerges, that the maximum correspondence between the SG and GRACE time series is achieved when filtering the GRACE data with Gaussian filters of about 1000 km filter length, which is in accordance with previous publications.Empirical Orthogonal Functions (EOF) analysis was applied to the gravity time series in order to identify common characteristic spatial and temporal patterns. The high correspondence of the first modes for GRACE and SG data implies that the first EOF mode represents a large-scale (Central European) time-variable gravity signal seen by both the GRACE satellites and the SG cluster.  相似文献   

10.
The impact of continental hydrological loading from land water, snow and ice on polar motion excitation, calculated as hydrological angular momentum (HAM), is difficult to estimate, and not as much is known about it as about atmospheric angular momentum (AAM) and oceanic angular momentum (OAM). In this paper, regional hydrological excitations to polar motion are investigated using monthly terrestrial water storage data derived from the Gravity Recovery and Climate Experiment (GRACE) mission and from the five models of land hydrology. The results show that the areas where the variance shows large variability are similar for the different models of land hydrology and for the GRACE data. Areas which have a small amplitude on the maps make an important contribution to the global hydrological excitation function of polar motion. The comparison of geodetic residuals and global hydrological excitation functions of polar motion shows that none of the hydrological excitation has enough energy to significantly improve the agreement between the observed geodetic excitation and geophysical ones.  相似文献   

11.
A Comparison of Global and Regional GRACE Models for Land Hydrology   总被引:1,自引:0,他引:1  
When using GRACE as a tool for hydrology, many different gravity field model products are now available to the end user. The traditional spherical harmonics solutions produced from GRACE are typically obtained through an optimization of the gravity field data at the global scale, and are generated by a number of processing centers around the world. Alternatives to this global approach include so-called regional techniques, for which many variants exist, but whose common trait is that they only use the gravity data collected over the area of interest to generate the solution. To determine whether these regional solutions hold any advantage over the global techniques in terms of overall accuracy, a range of comparisons were made using some of the more widely used regional and global methods currently available. The regional techniques tested made use of either spherical radial basis functions or single layer densities (i.e., mascons), with the global solutions having been obtained from the various major processing centers. The solutions were evaluated using a range of computed statistics over a selection of major river basins, which were globally distributed and ranged in size from 1 to 6 million km2. For one of the basins tested, the Zambezi, additional validation tests were conducted through comparisons against a custom designed regional hydrology model of the region. We could not prove that current regional models perform better than global ones. Monthly mean water storage variations agree at the level of 0.02 m equivalent water height. The differences in terms of monthly mean water storage variations between regional and global solutions are comparable with the differences among only global or regional solutions. Typically they reach values of 0.02 m equivalent water heights, which seems to be the level of accuracy of current GRACE solutions for river basins above 1 million km2. The amplitudes of the seasonal mass variations agree at the sub-centimetre level. Evident from all of the comparisons shown is the importance that the choice of regularization, or spatial filtering, can have on the solution quality. This was found to be true for global as well as regional techniques.  相似文献   

12.
The interaction between surface water and groundwater is an important aspect of hydrological processes. Despite its importance, groundwater is not well represented in many land surface models. In this study, a groundwater module with consideration of surface water and groundwater dynamic interactions is incorporated into the distributed biosphere hydrological (DBH) model in the upstream of the Yellow River basin, China. Two numerical experiments are conducted using the DBH model: one with groundwater module active, namely, DBH_GW and the other without, namely, DBH_NGW. Simulations by two experiments are compared with observed river discharge and terrestrial water storage (TWS) variation from the Gravity Recovery and Climate Experiment (GRACE). The results show that river discharge during the low flow season that is underestimated in the DBH_NGW has been improved by incorporating the groundwater scheme. As for the TWS, simulation in DBH_GW shows better agreement with GRACE data in terms of interannual and intraseasonal variations and annual changing trend. Furthermore, compared with DBH_GW, TWS simulated in DBH_NGW shows smaller decreases during autumn and smaller increases in spring. These results suggest that consideration of groundwater dynamics enables a more reasonable representation of TWS change by increasing TWS amplitudes and signals and as a consequence, improves river discharge simulation in the low flow seasons when groundwater is a major component in runoff. Additionally, incorporation of groundwater module also leads to wetter soil moisture and higher evapotranspiration, especially in the wet seasons.  相似文献   

13.
West African countries have been exposed to changes in rainfall patterns over the last decades, including a significant negative trend. This causes adverse effects on water resources of the region, for instance, reduced freshwater availability. Assessing and predicting large-scale total water storage (TWS) variations are necessary for West Africa, due to its environmental, social, and economical impacts. Hydrological models, however, may perform poorly over West Africa due to data scarcity. This study describes a new statistical, data-driven approach for predicting West African TWS changes from (past) gravity data obtained from the gravity recovery and climate experiment (GRACE), and (concurrent) rainfall data from the tropical rainfall measuring mission (TRMM) and sea surface temperature (SST) data over the Atlantic, Pacific, and Indian Oceans. The proposed method, therefore, capitalizes on the availability of remotely sensed observations for predicting monthly TWS, a quantity which is hard to observe in the field but important for measuring regional energy balance, as well as for agricultural, and water resource management. Major teleconnections within these data sets were identified using independent component analysis and linked via low-degree autoregressive models to build a predictive framework. After a learning phase of 72 months, our approach predicted TWS from rainfall and SST data alone that fitted to the observed GRACE-TWS better than that from a global hydrological model. Our results indicated a fit of 79 % and 67 % for the first-year prediction of the two dominant annual and inter-annual modes of TWS variations. This fit reduces to 62 % and 57 % for the second year of projection. The proposed approach, therefore, represents strong potential to predict the TWS over West Africa up to 2 years. It also has the potential to bridge the present GRACE data gaps of 1 month about each 162 days as well as a—hopefully—limited gap between GRACE and the GRACE follow-on mission over West Africa. The method presented could also be used to generate a near-real-time GRACE forecast over the regions that exhibit strong teleconnections.  相似文献   

14.
Since its launch in March 2002, the Gravity Recovery and Climate Experiment (GRACE) has provided a global mapping of the time-variations of the Earth’s gravity field. Tiny variations of gravity from monthly to decadal time scales are mainly due to redistributions of water mass inside the surface fluid envelops of our planet (i.e., atmosphere, ocean and water storage on continents). In this article, we present a review of the major contributions of GRACE satellite gravimetry in global and regional hydrology. To date, many studies have focused on the ability of GRACE to detect, for the very first time, the time-variations of continental water storage (including surface waters, soil moisture, groundwater, as well as snow pack at high latitudes) at the unprecedented resolution of ~400–500 km. As no global complete network of surface hydrological observations exists, the advances of satellite gravimetry to monitor terrestrial water storage are significant and unique for determining changes in total water storage and water balance closure at regional and continental scales.  相似文献   

15.
局部Slepian函数是将局部区域内的地球物理信号转化为空间谱的一种方法,其可以保证在球面上局部范围内获得最优谱平滑解,非常适用于局部范围地球物理信号的研究.本文利用中国陆态网西南地区72个测站的连续GPS观测资料分析川云渝地区陆地水负荷形变特征,并基于Slepian函数方法解算60阶的空间谱基函数,结合弹性质量负荷理论研究了川云渝地区2011年至2015年陆地水储量变化的时空分布模式.针对Slepian函数的边界效应问题,本文使用GLDAS格网数据计算得到站点处垂直负荷位移时间序列,然后利用该位移数据来进行水储量变化恢复实验,结果表明当边界扩充为3°时能较好地恢复GLDAS模型输出的陆地水储量变化.通过对比区域内GPS、GRACE、GLDAS得到的等效水高以及降雨数据,发现季节性降水是陆地水变化的一个重要驱动因子,GPS反演结果与GRACE和GLDAS数据具有较强的空间一致性.云南地区周年变化要强于川渝地区,其中云南西部的山区陆地水变化最大,约为30 cm,最小为川北以及重庆地区仅为7 cm.相较于GPS反演结果,GRACE与GLDAS明显低估了陆地水储量的季节性变化,分别达到24%和47%.比较分析地区内平均等效水高时间序列的相位发现,GPS得到的陆地水变化与降雨数据一致性较好,而GRACE与GLDAS存在一到两个月左右的时延.同时GPS能较好的探测出2015年1月左右南方地区大范围的强降水,而GRACE与GLDAS并没有体现出该现象,说明GPS能更为灵敏地探测到局部地区陆地水的变化.在站点等效水高时间序列上,GPS与GRACE的相关性总体上要优于GPS与GLDAS,陆地水周年变化较大的云南和四川西部地区站点三种数据间相关性较好,而其他季节性信号不明显的地区则相关性较差.本文的研究表明运用GPS-Slepian方法能够独立地监测高时空分辨率的陆地水储量变化,是作为当前补充GRACE观测资料空缺期的有益尝试.  相似文献   

16.
本文利用GRACE (Gravity Recovery and Climate Experiment) 卫星重力资料研究了亚马逊流域2002-2010年的陆地水变化,并与水文模式和降雨资料进行了比较分析.在年际尺度上,GRACE结果表明:2002-2003年和2005年,亚马逊流域发生明显的干旱现象;2007年至2009年,陆地水呈逐年增加的趋势,并在2009年6月变化值达到最大,为772±181 km3;自2009年6月至2010年12月,陆地水总量又急剧减少了1139±262 km3,这相当于全球海平面上升3.2±0.7 mm所需的水量.水文模式得到的亚马逊流域陆地水在2010年也表现出明显的减少.降雨资料与GRACE观测资料有很好的一致性.在2005年和2010年的干旱期,亚马逊流域的降雨显著减少,说明降雨是亚马逊流域陆地水变化的重要因素.此外,本文采用的尺度因子的方法有效地降低了GRACE后处理误差的影响.  相似文献   

17.
华北平原作为我国重要的工农业基地和政治经济中心,面临着严重的水资源危机.因此,开展对华北平原地下水储量变化的监测工作具有重要现实意义与科学价值.本文基于GRACE重力卫星的空间约束方法,研究了华北平原地下水储量变化的时空分布规律,并与地面水井实测与地下水模型结果进行了综合比较和分析.结果表明:2002-2014年,华北平原地下水存在明显的长期亏损,GRACE估计的亏损速率为-7.4±0.9 km~3·a~(-1),而地面水井资料估计的浅层地下水亏损速率为-1.2 km~3·a~1,对比两者之间的差异可以发现,华北平原的地下水亏损以深层地下水为主.2002-2008年,GRACE估计的华北平原地下水亏损速率为-5.3±2.2 km~3·a~(-1),这与华北平原两个地下水模型得到的平均亏损速率-5.4 km~3·a~(-1)十分吻合.通过华北平原区域地下水模型的独立验证,说明GRACE可以有效评估华北平原的地下水储量变化趋势.除了长期亏损的趋势项之外,华北平原地下水还存在明显的年际变化特征,并与该地区年降雨量变化特征一致.在降雨偏少的2002年、2005-2009年和2014年,华北平原地下水储量显著减少.在空间分布上,GRACE结果表明,华北平原的地下水储量减少主要发生在山前平原和中部平原区,这也与水井实测资料和区域地下水模型结果较为吻合.与GRACE和区域地下水模型相比,目前的全球水文模型仍无法准确估计华北平原地下水变化的空间分布和亏损速率.上述研究表明,GRACE提供了评估华北平原地下水储量变化的重要监测手段.  相似文献   

18.
Seasonal water storage change of the Yangtze River basin detected by GRACE   总被引:13,自引:0,他引:13  
1 Introduction Large-scale mass redistribution, or temporal varia- tion of mass within the Earth system, the driving force of interactions between solid Earth and geophysical fluids envelope (i.e., atmosphere, ocean, and hydro- sphere), is an important geophysical process critical to human life. Most of the interactions between solid Earth and the atmosphere/oceans happen at seasonal and inter-annual time scales. One important contribu- tor of mass redistribution at seasonal and inter-annual …  相似文献   

19.
《Journal of Geodynamics》2010,49(3-5):157-165
Since 2002 the Earth’s gravity field is globally observed by the Gravity Recovery and Climate Experiment (GRACE) satellite mission. The GRACE monthly gravity field solutions, available from several analysis centres, reflect mass variations in the atmosphere, hydrosphere and geosphere. Due to correlated noise contained in these solutions, it is, however, first necessary to apply an appropriate filtering technique. The resulting, smoothed time series are applied not only to determine variations with different periodic signatures (e.g., seasonal, short and medium-term), but to derive long-periodic mass variations and secular trends as well. As the GRACE monthly solutions always show the integral effect of all mass variations, for separation of single processes, like the GIA (Glacial isostatic adjustment)-related mass increase in Fennoscandia, appropriate reduction models (e.g. from hydrology) are necessary.In this study we show for the example of the Fennoscandian uplift area that GRACE solutions from different analysis centres yield considerably different secular trends. Furthermore, it turns out that the inevitable filtering of the monthly gravity field models affects not only the amplitudes of the signals, but also their spatial resolution and distribution such as the spatial form of the detected signals. It also becomes evident that the determination of trends has to be performed together with the determination of periodic components. All periodic terms which are really contained in the data, and only such, have to be included. The restricted time span of the available GRACE measurements, however, limits the separation of long-periodic and secular signals. It is shown that varying the analysis time span affects the results considerably. Finally, a reduction of hydrological signals from the detected integral secular trends using global hydrological models (WGHM, LaDWorld, GLDAS) is attempted. The differences among the trends resulting from different models illustrate that the state-of-the-art hydrology models are not suitable for this purpose as yet. Consequently, taking the GRACE monthly gravity field solutions from one centre, choosing a single filter and applying an insufficiently reliable reduction model leads sometimes to a misinterpretation of considered geophysical processes. Therefore, one has to be cautious with the final interpretation of the results.  相似文献   

20.
Global Terrestrial Water Storage Changes and Connections to ENSO Events   总被引:1,自引:0,他引:1  
Improved data quality of extended record of the Gravity Recovery and Climate Experiment (GRACE) satellite gravity solutions enables better understanding of terrestrial water storage (TWS) variations. Connections of TWS and climate change are critical to investigate regional and global water cycles. In this study, we provide a comprehensive analysis of global connections between interannual TWS changes and El Niño Southern Oscillation (ENSO) events, using multiple sources of data, including GRACE measurements, land surface model (LSM) predictions and precipitation observations. We use cross-correlation and coherence spectrum analysis to examine global connections between interannual TWS changes and the Niño 3.4 index, and select four river basins (Amazon, Orinoco, Colorado, and Lena) for more detailed analysis. The results indicate that interannual TWS changes are strongly correlated with ENSO over much of the globe, with maximum cross-correlation coefficients up to ~0.70, well above the 95% significance level (~0.29) derived by the Monte Carlo experiments. The strongest correlations are found in tropical and subtropical regions, especially in the Amazon, Orinoco, and La Plata basins. While both GRACE and LSM TWS estimates show reasonably good correlations with ENSO and generally consistent spatial correlation patterns, notably higher correlations are found between GRACE TWS and ENSO. The existence of significant correlations in middle–high latitudes shows the large-scale impact of ENSO on the global water cycle.  相似文献   

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

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