首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 21 毫秒
1.
北疆地区积雪时空变化的影响因素分析   总被引:2,自引:0,他引:2  
利用新疆北部地区2000-2007年观测的积雪资料分析北疆地区积雪开始时间、积雪结束时间、积雪日数、年最大雪深、积雪期平均雪深和年平均雪深随海拔、经纬度、坡度、坡向和植被的变化. 结果表明:随着海拔增加,积雪各变量变化明显,温度在海拔变化中起着关键作用;植被对积雪各变量有影响,但影响程度不明显. 在北疆区域范围内,纬度变化及温度差异不大,对积雪各变量影响很小;经度对积雪各变量的影响是由空间差异造成的;坡度对积雪变量的影响主要通过空间分布及坡度产生的阴影造成,进而影响太阳直射;坡向对积雪各变量的影响主要由水汽运动方向和太阳光照造成. 因此,在北疆区域内,对积雪各变量的影响程度为海拔>坡向>坡度>植被>纬度>经度.  相似文献   

2.
A comparison of the stable-isotope signatures of spring waters, snow, snowmelt, summer (July thru September) rain, and cool season (October thru June) rain indicates that the high-intensity, short-duration summer convective storms, which contribute approximately a third of the annual precipitation to the Spring Mountains, provide only a small fraction (perhaps 10%) of the recharge to this major upland in southern Nevada, USA. Late spring snowmelt is the principal means of recharging the fractured Paleozoic-age carbonate rocks comprising the central and highest portion of the Spring Mountains. Daily discharge measurements at Peak Spring Canyon Creek during the period 1978–94 show that snowpacks were greatly enhanced during El Niño events.  相似文献   

3.
闫玉娜  车涛  李弘毅  秦越 《冰川冻土》2016,38(1):211-221
随着寒区水文模拟研究的深入,春季融雪径流模拟误差大这一问题越来越明显.针对此问题,以八宝河流域为研究区,利用积雪衰减曲线将MODIS积雪面积比例产品转化为雪水当量,并用其更新分布式水文模型GBHM(Geomorphology-Based Hydrological Model)中模拟的雪水当量,达到提高春季融雪径流模拟精度的目的.利用GBHM模型对八宝河流域2005-2007年进行了模型预热,参数率定,同时选择2008年进行模型检验.结果表明:GBHM模型在八宝河流域有较好的径流模拟精度,年均纳什效率系数(NSE)达到0.64.但模型对春季融雪过程的模拟效果较差,通过引入积雪遥感数据,这一问题得到很大改善.加入积雪遥感数据后,2008年3-6月径流模拟精度NSE和相对偏差Bias分别由-1.0、-0.45改进为0.58、0.06,单点雪水当量模拟精度获得提高,流域水量平衡也更加合理.  相似文献   

4.
气候变化对中国西北地区山区融雪径流的影响   总被引:22,自引:16,他引:22  
选择祁连山黑河流域作为中国西北地区山区积雪流域的典型代表,分析了1956-1995年40a以来气候,积雪变化的状况和特点以及春季融雪径 波动趋势,利用融雪径流模型(Snowmelt Runoff Model-SRM)和卫星遥感数据模拟气温上升框架上的融雪径流变化情势,结果表明,中国西北地区山区的气候变化主要表现在年平均气温的缓慢上升而降水基本平稳,年内气温的上升幅度以1-2月份比较强烈,而3-6月融雪期的气温并没有大的变化,导致融雪期在时间尺度上的扩大,融雪径流呈慢增加趋势且受径流周期变化控制,融雪径流峰值的时间上前移。  相似文献   

5.
The Gradenbach mass movement (GMM) is an example of DGSD (deep-seated gravitational slope deformation) in crystalline rocks of the Eastern Alps (12.85°E, 47.00°N). The main body of the GMM covers an area of 1.7 km2 and its volume is about 120?×?106 m3. A reconstruction of the deformation history yields a mean displacement of?~?22 m from 1962 to 2011. In 1965/66, 1975, 2001, and 2009 high sliding velocities, exceeding several meters per year, interrupt the quasi-stationary periods of slow movement (≤0.3 m/year). Since 1999 the displacement of the main body of the GMM has been observed by GPS. Time series of extensometer readings, precipitation, snow cover water equivalent, water discharge, and hydrostatic water level observed in boreholes were re-processed and are presented in this paper. Continuous recording of seismic activity by a seismic monitoring network at the GMM began in the summer of 2006. Deformation has been monitored since 2007 by an embedded strain rosette based on fiber optics technology and a local conventional geodetic deformation network. The velocity of the GMM could be modeled to a large extent by a quantitative relation to hydro-meteorological data. During the phase of high sliding velocity in spring 2009, the seismic activity in the area increased significantly. Several types of seismic events were identified with some of them preceding the acceleration of the main body by about 6 weeks. The potential inherent in the Gradenbach Observatory data to supply early warning and hazard estimation is discussed.  相似文献   

6.
北京-张家口地区冬春季积雪特征分析   总被引:2,自引:4,他引:2  
2022年冬奥会将在北京-张家口(以下简称北-张地区)举办,揭示该地区的积雪变化特征及其在全球变暖背景下的发展趋势,对冬奥会的筹备以及当地的积雪资源的开发利用等方面都有重要意义。利用2002-2014年MODIS遥感积雪产品提取了研究区域积雪数据,结合1966-2013年台站积雪、气温和降水资料和DEM数据,分析了积雪的时空分布特征,并对冬奥会场地进行积雪资源评价。结果表明:2002-2013冬春年北张地区的整体积雪频率较小,多处于0~0.2之间,但场馆区2月的积雪频率多在0.5以上,最大值接近0.9左右,积雪的分布呈带状和点状。积雪覆盖率最大值出现在1月初,达到0.23。积雪的形成缓慢,但是消亡迅速。1966-2012冬春年冬季积雪日数的波动幅度大于春季,延庆和崇礼县的2月份积雪日数分别为4.6d和13.9d,且均呈下降状态。积雪初终日均有提前,但整体的积雪期在减少。北京和张家口整体的最大积雪深度变化平稳,在1966-1980年和2000-2012年处于高值区,波动较大,其他年份最大雪深处于低值变化平稳,延庆和崇礼县的2月份最大积雪深度分别为3.6cm和5.1cm。通过分析积雪指标与气象因子(气温、降水)的相关关系发现,在年内(年际)变化上,积雪指标与气温(降水)的关系更为密切。冬奥会场地的2月份气温在-14~2℃之间,月平均降水量仅0.2mm·d-1,积雪日数不足,预计难以形成足够深度的雪,且未来气温上升,达到0.8℃·(10a)-1,降水、积雪深度和积雪日数均呈下降趋势,可能60%~95%的赛事用雪将来自人造雪,以应对可能的积雪不足。  相似文献   

7.
During the formation and development of glacial meltwater runoff, hydrochemical erosion is abundant, especially the hydrolysis of K/Na feldspar and carbonates, which can consume H+ in the water, promote the formation of bicarbonate by dissolving atmospheric CO2, and affect the regional carbon cycle. From July 21, 2015, to July 18, 2017, the CO2 concentration and flux were observed by the eddy covariance (EC) method in the relatively flat and open moraine cover area of Koxkar Glacier in western Mt. Tianshan, China. We found that: (1) atmospheric CO2 fluxes ranged from ??408.95 to 81.58 mmol m?2 day?1 (average ? 58.68 mmol m?2 day?1), suggesting that the study area is a significant carbon sink, (2) the CO2 flux footprint contribution areas were primarily within 150 m of the EC station, averaging total contribution rates of 93.30%, 91.39%, and 90.17% of the CO2 flux in the snow accumulation, snow melting, and glacial melting periods, respectively. Therefore, the contribution areas with significant influences on CO2 flux observed at EC stations were concentrated, demonstrating that grassland CO2 flux around the glaciers had little effect at the EC stations, (3) in the predominant wind direction, under stable daytime atmospheric stratification, the measurement of CO2 flux, as interpreted by the Agroscope Reckenholz Tanikon footprint tool, was 79.09% ± 1.84% in the contribution area. This was slightly more than seen at night, but significantly lower than the average under unstable atmospheric stratification across the three periods of interest (89%). The average distance of the farthest point of the flux footprint under steady state atmospheric conditions was 202.61?±?69.33 m, markedly greater than that under non-steady state conditions (68.55?±?10.34 m). This also indicates that the CO2 flux observed using EC was affected primarily by hydrochemical erosion reactions in the glacier area, (4) a good negative correlation was found between net glacier exchange (NGE) of CO2 and air temperature on precipitation-free days. Strong ice and snow ablation could promote hydrochemical reactions of soluble substances in the debris area and accelerated sinking of atmospheric CO2. Precipitation events might reduce snow and ice melting, driven by reduced regional temperatures. However, a connection between NGE and precipitation, when less than 8.8 mm per day, was not obvious. When precipitation was greater than 8.8 mm per day, NGE decreased with increasing precipitation, (5) graphically, the slope of NGE, related to daily runoff, followed a trend: snow melting period?>?snow accumulation period?>?early glacial ablation period?>?late glacier ablation period?>?dramatic glacier ablation period. The slope was relatively large during snow melting, likely because of CO2 sinking caused by water–rock interactions. The chemical reaction during elution in the snow layer might also promote atmospheric CO2 drawdown. At the same time, the damping effect of snow cover and the almost-closed glacier hydrographic channel inhibited the formation of regional runoff, possibly providing sufficient time for the chemical reaction, thus promoting further CO2 drawdown.  相似文献   

8.
Heavy rainfall in June 2013 triggered flash flooding and landslides throughout the Indian Himalayan state of Uttarakhand, killing more than 6000 people. The vast majority of fatalities and destruction resulted directly from a lake outburst and debris flow disaster originating from above the village of Kedarnath on June 16 and 17. Here, we provide a systematic analysis of the contributing factors leading to the Kedarnath disaster, both in terms of hydrometeorological triggering and topographic predisposition. Topographic characteristics of the lake watershed above Kedarnath are compared with other glacial lakes across the north-western Himalayan states of Uttarakhand and Himachal Pradesh, and implications for glacier lake outburst hazard assessment in a changing climate are discussed. Our analysis suggests that the early onset of heavy monsoon rainfall (390 mm, June 10–17) immediately following a 4-week period of unusually rapid snow cover depletion and elevated streamflow was the crucial hydrometeorological factor, resulting in slope saturation and significant run-off into the small seasonal glacial lake. Between mid-May and mid-June 2013, snow-covered area above Kedarnath decreased by around 50 %. The unusual situation of the lake being dammed in a steep, unstable paraglacial environment but fed entirely from snowmelt and rainfall within a fluvial dominated watershed is important in the context of this disaster. A simple scheme enabling large-scale recognition of such an unfavourable topographic setting is introduced. In view of projected 21st century changes in monsoon timing and heavy precipitation in South Asia, more emphasis should be given to potential hydrometeorological triggering of lake outburst and debris flow disasters in the Himalaya.  相似文献   

9.
1981-2010年青藏高原积雪日数时空变化特征分析   总被引:2,自引:0,他引:2  
全球气候变暖大背景下, 作为冰冻圈最为活跃和敏感因子, 青藏高原积雪变化备受国内外关注. 本文利用青藏高原(以下简称高原)1981-2010年地面观测积雪日数资料, 较系统地分析了近30年来高原积雪日数的时空变化特点. 主要结论如下: (1) 近30年内高原平均年积雪日数出现了非常显著的减少趋势, 减少幅度达4.81 d·(10a)-1, 其中冬季减幅最为明显, 为2.36 d·(10a)-1, 其次是春季(2.05 d·(10a)-1), 而夏季最少(0.21 d·(10a)-1); (2) 30年间, 积雪日数较少的年份多数出现在本世纪初10年内, 且2010年属于异常偏少年, 高原积雪日数在1997年左右发生了由多到少的气候突变; (3) 在空间上, 北部柴达木盆地及其附件区域部分气象台站观测的年积雪日数出现了不显著的增加趋势之外, 高原91.5%的气象站年积雪日数呈减少趋势, 且高寒内陆中东部和西南喜马拉雅山脉南麓等高原历年积雪日数高值区域减少最为明显; (4) 由于受到气象台站所在地理位置、地形地貌、地表类型、海拔高度、局地气候以及大气环流等综合影响, 高原平均年积雪日数的空间差异很大, 最多达146 d, 最少的则不足1 d, 平均仅为38 d, 其中高寒内陆中东部是积雪日数最长的区域, 而东南部海拔和纬度较低的干热河谷地区积雪日数最少.  相似文献   

10.
H. P. Hong  W. Ye 《Natural Hazards》2014,73(2):355-371
Snow depth records from daily measurements at climatological stations were obtained from Environment Canada and were processed and analyzed. It was identified that there are 549 stations, each with at least 20 years of useable annual maximum snow depth data. Both the Gumbel distribution and generalized extreme value distribution were used to fit the annual maximum snow depth, considering several distribution fitting methods. Statistical analysis results indicated that, according to the Akaike information criterion, the Gumbel distribution is preferred for 72 % stations. The estimated return period value of annual maximum snow depth at stations was used to calculate their corresponding ground snow load. The at-site analysis results were used as the basis to spatially interpolate the ground snow loads for locations tabulated in the National Building Code of Canada (NBCC) since a code location and a climatological site are usually not co-located. For the interpolation, the ordinary co-kriging method with elevation as co-variate was used because a cross-validation analysis by using several deterministic and probabilistic spatial interpolation techniques indicated that the ordinary co-kriging method is preferred. A comparison of the newly estimated ground snow loads to those locations tabulated in the 1995 edition and 2010 edition of the NBCC was also presented.  相似文献   

11.

The gully systems develop a compact network dissecting the loess cover of the NW part of the Lublin Upland. They are overgrown by dense and multi-storey forest assemblages. It is very difficult to apply teledetection techniques, GPS positioning, and geodesic tools, such as a laser rangefinder or total station, for monitoring of gully development. The detailed research, involving the comparison of traditional and modern methods, was conducted in a gully with a secondary erosional dissection, developing in the gully bottom with varied intensity for the last 15 years. In the period of comparative research (2012/2013), the rate of development of the erosional dissection was measured by means of traditional methods as well as a Leica HDS 7000 laser scanner. 3D laser scanning permitted precise terrain modelling and the assessment of qualitative and quantitative spatial changes. This paper compares the results of three measurement campaigns. The analysis showed intensive development of the dissection, manifested in an increase in its volume (37 %), retreat of the erosion scarp (32 m), deep erosion (19 %), and lateral erosion (27 %). Registered moderate snowmelt in 2013 caused the volume dissection to increase by 57 m3, mainly as a result of the gully bottom scouring. During heavy rainfalls in June, the landforms were widened by a further 217 m3 and observed deepening and widening of erosive potholes. The prevalence of deep erosion resulted in the transformation back into a V-shaped cross section. A total of 37.3 Mg km−2 of material was discharged from the gully catchment during snowmelt runoffs in 2013. A precipitation event in June 2013 caused the discharge of 856.9 Mg km−2 of material from the catchment. The obtained results confirm that snowmelt erosion determines the spatial development of perennial gullies, but the volume of material discharged from the catchment depends on the occurrence of extreme precipitation during the warm season.

  相似文献   

12.
亚洲中部干旱区的大尺度遥感积雪信息研究,可在跨界河流水资源分配利用方面提供数据支持,对国家重大战略的生态安全保障有重要作用。采用数据融合方法,将MOD10A2和MYD10A2数据进行融合去云处理,结合气象站点积雪数据评估去云后的积雪识别精度;提取积雪覆盖率(SCP)与积雪日数(SCD)信息,分析SCP与SCD年际、年内变化差异;结合数字高程模型,分析不同高程带下SCP的时空变化规律。结果表明:(1)MOD10A2与MYD10A2融合去云处理,可有效去除云的干扰,准确提取亚洲中部干旱区积雪变化信息。(2)年内SCP最大值范围为55.7%~77.4%,最小值范围为1.6%~2.9%,融雪期SCP下降速率具有明显地域差异,总体SCP呈缓慢增加趋势。(3)总体SCD呈略微下降趋势,32.2%的区域呈下降趋势,30.9%的区域呈增加趋势,36.9%的区域保持稳定不变。(4)海拔1 000 m以下,SCP年内随季节变化呈U型,年际变化显著;1 000~4 000 m区域,SCP年内均随季节的变化呈现出V型,年际变化呈现出稳定性波动;6 000 m以上为永久性积雪,季节、时空变化差异性均不明显。  相似文献   

13.
利用已被广泛使用的MODIS积雪数据, 获得了塔里木河源区之一的托什干河流域积雪变化信息. 结果表明: 流域积雪覆盖率时空差异显著, 在积雪丰富的年份, 1月积雪覆盖率可达90%以上, 但在积雪少的年份, 则只有50%; 2000年以来流域积雪呈现微弱增加的趋势, 积雪变化趋势呈现明显的时空差异. 相对于其他季节, 流域冬季积雪增加更为明显; 与其他高度带相比, 作为主要积雪覆盖区的海拔3 000~4 000 m高度带积雪的增加趋势也更为明显. 以流域所在的气象格网数据和积雪覆盖率变化曲线作为输入数据, 应用融雪径流模型(SRM)模拟了流域春季融雪径流过程, 率定了模型主要参数, 获得了较好的结果. 以CMIP5的3种RCP情景为驱动数据, 应用模型预估了流域2021-2050年的融雪径流状况, 结果显示 4月之前径流变化不明显, 之后径流峰值增大显著, 不同气候情景对径流的影响不明显.  相似文献   

14.
For the sustainable utilization of rivers in the mid and downstream regions, it is essential that land surface hydrological processes are quantified in high cold mountains regions, as it is in these regions where most of the larger rivers in China acquire their headstreams. Glaciers are one of the most important water resources of north-west China. However, they are seldom explicitly considered within hydrological models, and climate-change effects on glaciers, permafrost and snow cover will have increasingly important consequences for runoff. In this study, an energy-balance ice-melt model was integrated within the Variable Infiltration Capacity (VIC) macroscale hydrological model. The extended VIC model was applied to simulate the hydrological processes in the Aksu River basin, a large mountainous and glaciered catchment in north-west China. The runoff components and their response to climate change were analyzed based on the simulated and observed data. The model showed an acceptable performance, and achieved an efficiency coefficient R 2 ≈ 0.8 for the complete simulation period. The results indicate that a large proportion of the catchment runoff is derived from ice meltwater and snowmelt water. In addition, over the last 38 years, rising temperature caused an extension in the snow/ice melting period and a reduction in the seasonality signal of runoff. Due to teh increased precipitation runoff, the Aksu catchment annual runoff had a positive trend, increasing by about 40.00 × 106 m3 per year, or 25.7 %.  相似文献   

15.
Snowmelt runoff is an important source of water resources in the arid mountain area. Modelling snowmelt runoff for cold regions remains a problematic aspect because of the lack of data by gauges in large basins. In order to overcome the shortage of measured data in the snowmelt runoff modelling, the temperature interpolation method would greatly help in improving the simulation accuracy and describing the snow-hydrological behaviours of the study catchments. In this study, the temperature is the principal variable used to estimate the importance of the melting of snow cover using the snowmelt runoff model. Five different temperature interpolation attempts were performed over the Kaidu River Basin for the snowmelt season of the year 2000. Three temperature inputs were taken directly from the individual weather stations in or near the study area, and the other two temperature inputs were interpolated from the three weather stations. The results indicated that the temperature estimated from different methods could result in quite a difference in runoffs in comparison with the observed ones. The simulation results using average temperature from the three stations showed good results; the simulation run with the weighted average temperature generated a lower R 2 than the average temperature of three stations and using temperature directly adopted from three individual stations gave various results. The weather stations used to perform the snowmelt runoff simulation should be located in the place which is most representative of the mountain weather conditions, and the land cover and topography that those stations represented also play an important role in the snowmelt runoff simulation.  相似文献   

16.
Shallow landslides are fairly frequent natural processes which emerge as a result of both rainfall and rapid snowmelt in the Flysch Belt of the Outer Western Carpathians. We estimated the total water content thresholds for the previously defined seven phases of increased landsliding which took place between 1939 and 2010 around the Napajedla meteorological station. The time series were reconstructed on the basis of data from surrounding stations. Rainfalls with the highest intensities (>1 mm/min) were removed from the set. Rainfall of such an intensity primarily causes overland flow and soil erosion and does not contribute to landslide threshold. The snow water equivalent was computed on the basis of the snow height, and possible errors were evaluated as interval estimations. An interval of 10 days before a landslide phase was selected for the total water content threshold. The resulting lower boundary (67.0 mm/10 days) and upper boundary (163.3 mm/10 days) thresholds of water infiltrated into soil during an event shall be part of the prepared online warning system in this area.  相似文献   

17.
1979-2014年东北地区雪深时空变化与大气环流的关系   总被引:2,自引:2,他引:0  
基于被动微波遥感反演的雪深数据集(1979-2014年),利用Mann-Kendall检验、R/S分析、相关分析和小波分析等方法研究了东北地区雪深时空变化特征及其与大气环流的关系。结果表明:1979-2014年,东北地区年均雪深总体呈减小趋势,减小速率为-0.084 cm·(10a)-1。其中,春季雪深减小速率最大,为-0.19 cm·(10a)-1P<0.05),其次是冬季[-0.17 cm·(10a)-1],而秋季雪深减小速率最小,仅为-0.05 cm·(10a)-1。空间上,平原区(东北平原和三江平原)与少部分高原区(呼伦贝尔高原西南部)年均雪深呈增大趋势,山地(大、小兴安岭和长白山)与高原大部(内蒙古高原)雪深呈减小趋势,而且雪深增大区域的面积和变化速率均小于雪深减小的地区。东北地区年均雪深变化的Hurst指数为0.85,表明雪深未来减小的持续性很强;同时雪深变化具有22 a的主周期。春秋季雪深变化与东亚槽强度及北半球极涡面积呈显著负相关性,而冬季雪深与北半球副高强度关系密切。  相似文献   

18.
我国东北低山区不同坡位积雪特性研究   总被引:4,自引:3,他引:1  
曹志  范昊明 《冰川冻土》2017,39(5):989-996
为研究坡位对积雪性质的影响,利用Snow Fork雪特性分析仪等采集积雪物理性质(积雪深度、密度、液态水含量)数据,分析了不同坡位积雪特性的变化趋势、差异及成因。结果表明:坡位的差异可显著影响积雪,特别是积雪表层的特性,积雪表层温度与积雪反辐射强度呈显著正相关关系;阴、阳坡各坡位积雪液态水含量变化趋势一致,但阳坡液态水含量最大值出现在上坡位,阴坡则是中坡位最大;阴、阳坡各坡位在液态水含量增加的影响下雪密度也在逐渐增大,并且随着积雪深度的增加雪密度逐渐减小;试验区融雪期阳坡积雪液态水含量最先对环境变化做出响应,进而影响雪密度,深度随之响应;阴坡则是液态水含量首先响应,积雪深度次之,密度最后响应。研究结果将为融雪径流形成、融雪侵蚀防治以及季节性积雪区生态系统评估提供理论依据。  相似文献   

19.
利用新疆89个地面站逐日积雪深度观测资料,研究探讨了1961—2017年新疆区域积雪期、积雪初日、积雪终日的时空变化规律,分析了北疆和天山山区积雪期的年代际和周期变化特征及其与气温、降水的关系.结果表明:新疆各地积雪期、积雪初日和终日存在明显的差异,积雪期以天山为界北多南少;从空间分布看,天山山区和新疆北部阿勒泰、塔城...  相似文献   

20.
青藏高原高寒草地植被指数变化与地表温度的相互关系   总被引:3,自引:1,他引:2  
为了解脆弱的高原生态环境对升温过程的响应, 利用1982-2006年国家标准地面气象站地表温度和GIMMS-NDVI数据集, 探讨了青藏高原高寒草地植被指数和地表温度的变化特征及其相互关系. 结果表明:1982-2006年, 高寒草地NDVI、地表温度整体均呈现增加趋势, 年均NDVI、生长季NDVI、年最大NDVI(NDVImax)与年均地表温度、生长季地表温度的上升趋势分别为0.007 (10a)-1、0.011 (10a)-1、0.007 (10a)-1与0.60 ℃·(10a)-1、0.43 ℃·(10a)-1; NDVImax与地表温度显著相关的地区达70.49%. 但是高原地形、气候、水文环境的空间差异性导致高寒草地NDVI与地表温度的相关关系十分复杂. NDVImax与年均地表温度的相关性最为显著; 在返青期和枯萎期, NDVI与地表温度均为显著正相关. 不同的植被覆盖条件下, NDVI对地表温度的响应不同:植被覆盖差以及退化严重的地区, NDVImax与地表温度呈负相关性; 反之, NDVImax与地表温度主要表现为正相关.  相似文献   

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

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