首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 187 毫秒
1.
王文种  黄对  刘九夫  刘宏伟  王欢 《湖泊科学》2020,32(5):1552-1563
湖泊变化是气候变化的指示器.为探索利用单一短时间尺度的卫星水位数据源估算长时间序列的湖泊水量变化的可行性,本文利用短时间尺度(2016—2018年) Sentinel-3A合成孔径雷达高度计(SRAL)作为唯一卫星水位数据源,以藏北高原内陆湖泊当惹雍错为例,结合基于Landsat光学遥感数据提取的1988—2018年的湖泊面积,综合分析2016—2018年间的非结冰期遥感湖泊面积与遥感湖泊水位变化,基于该时段范围的水位变化-面积变化关系和水量估算公式,估算1988—2018年湖泊水位水量变化与2001—2018年的年内变化,并结合GLDAS产品数据与雪线变化情况初步探讨湖泊变化的可能原因.结果表明:当惹雍错近30年湖泊面积扩张明显,湖泊水位、水量增加显著,相比1988年,2018年的湖泊面积、水位、水量分别增加21.1 km2、5.29 m、44.75亿m3.其中1988—1998年湖泊面积-水位-水量有所减少,2000—2018年间湖泊变化总体呈增加趋势.2001—2018年内湖泊面积、水位、水量变化呈现干湿季特征.1996—2014/2015年湖泊水量变化为38.3亿/35.5亿m3,水量变化趋势、变化量与以往对应时间段的研究结果具有较强的一致性.湖泊面积扩张主要发生在水下地形平缓的东南部和中西部区域.结合气候因素与雪线变化的分析表明,湖泊水量变化受降雨、气温影响复杂,长时间年际尺度上的湖泊水量增长与气温的一致性较降水量强,湖泊湿季受降水量与气温的影响都较大,其中2008—2018年的湿季降水量、气温与水量变化散点拟合的确定性系数R2分别为0.613、0.845.该研究表明Sentinel-3A合成孔径雷达数据在湖泊水量变化估算上的潜力,为利用单一且只具有短时段数据的卫星雷达数据估算长时间序列湖泊水量变化提供依据.  相似文献   

2.
周期性的淹没或出露是洪泛型湖泊湿地的重要物理特征,湖泊湿地的淹没动态对其生态过程有显著影响。三峡工程运行以来,鄱阳湖湿地淹没动态发生了显著变化并引发了剧烈生态效应,而目前研究尚未对后三峡工程时代的鄱阳湖淹没动态演变进行系统量化,也制约了对其驱动下湖泊湿地生态系统演变原因与机制的了解。鉴于此,本研究结合水文站实测数据与遥感观测资料,以淹没开始时间、淹没结束时间以及淹没历时3个变量共同表征湖泊湿地淹没动态,从站点及全湖尺度分别对三峡工程运行后鄱阳湖湿地淹没动态的变化趋势、量级及显著性进行定量分析。结果表明:(1)2000—2020年间,鄱阳湖湿地淹没开始时间在64%的湖区被推迟,推迟速率约为1.10天/年;仅在入江水道及碟形湖小幅提前,且提前趋势并不显著;(2)鄱阳湖湿地淹没结束时间在72%的湖区显著提前,提前速率约为1.46天/年;仅在有闸控工程的碟形湖因延迟泄水而有所延迟;(3)受淹没开始时间推迟而结束时间提前的影响,鄱阳湖湿地淹没历时在70%的湖区显著缩短,缩短速率约为2.19天/年,而在有闸控工程的碟形湖则有所延长。本文从站点与全湖尺度分别给出了基于实测而非模型模拟的鄱阳湖淹没动态...  相似文献   

3.
鄱阳湖水面蒸发量的计算与变化趋势分析(1955-2004年)   总被引:2,自引:0,他引:2  
闵骞  刘影 《湖泊科学》2006,18(5):452-457
利用器测折算法与气候模式法,分别计算鄱阳湖周围康山、棠荫、都昌、星子、湖口5站的单站水面蒸发量,以5站两种方法计算值的平均值代表鄱阳湖大湖面的水面蒸发量,求得鄱阳湖1955-2004年各月的水面蒸发量和蒸发水量,结果为:多年平均年蒸发量1081.2 mm.年蒸发水量27.06×10~8 m~3.对年、月水面蒸发量在近50年来的变化趋势进行了分析,表明除5月份外,其他各月蒸发量和年蒸发量均呈逐渐减少趋势,年蒸发量平均每年减小2.79 mm,年蒸发水量平均减少0.05×10~8 m~3,对湖区水资源持续利用和湖泊环境将产生明显影响.对水面蒸发量递减原因进行了初步探讨.  相似文献   

4.
鄱阳湖水体夏季气温效应   总被引:3,自引:2,他引:1  
本文讨论了定量计算湖泊纯水体气温效应的方法,在计算鄱阳湖水体夏季纯水体气温效应值的基础上,重点分析鄱阳湖水体夏季气温效应的特征、强度和分布,其结果是鄱阳湖水体夏季对湖区气温影响呈热源效应。表现在提高了湖区夏季平均气温,湖中心棠荫站比湖外高O.1—1.0℃,纯水体气温效应为0.4℃。随着距水体距离的增大,气温效应也逐渐减弱,水体影响气温的范围约为10—40km。由于湖区夏季受到副热带高压控制,陆风比湖风旺盛,气温效应夜间比白天强,加上地形的作用,影响范围南部比北部大。1980—1987年间的7、8月考察资料也和上述结果相吻合。  相似文献   

5.
博斯腾湖水盐动态变化(1951-2011年)及对气候变化的响应   总被引:5,自引:1,他引:4  
分析了1951-2011年博斯腾湖历史水位和湖水矿化度的动态变化特征,解析了博斯腾湖水量与水质对气候变化的响应及未来变化趋势.结果表明,博斯腾湖水位在60年内经历了两个突变时期,突变时间分别为1974年和1994年,湖水矿化度也相应地呈现了三个动态变化阶段,水量与水位呈极显著负相关,但水质变化滞后于水位变化1年;流域气温呈显著增加趋势,气温升高的突变时间为1993年,与开都河出山口径流突变时间一致,但降水变化不显著;1993年前,博斯腾湖水量主要受气温和人类活动双重影响,1993年后博斯腾湖水量主要受气温的显著影响,气温主要通过改变入湖水量及湖区蒸发损耗来调控湖泊水位和水质;未来气温持续升高情景下,博斯腾湖水位将面临降低趋势,水质也将有恶化趋势.因此,为合理开发利用博斯腾湖水资源,减少水资源无效损耗,抑制水质恶化趋势,确保流域可持续发展,建议将博斯腾湖调水时间集中在5-9月,并严格控制孔雀河流域工农业用水量及工农业、生活污染源,减少污水排放量,减少周边地下水开发量.  相似文献   

6.
1973-2013年红碱淖水域水质变化及驱动力分析   总被引:3,自引:1,他引:2  
赵宁  马超  杨亚莉 《湖泊科学》2016,28(5):982-993
对1973-2013年8期Landsat MSS、TM、ETM+、OLI影像进行了辐射定标、大气校正、辐射归一化和波段运算等处理;利用归一化差值植被指数(NDVI),分期提取了红碱淖水域面积,分析了湖水水质及红碱淖周围植被变化.研究表明:40 a间红碱淖水域面积呈阶段性萎缩趋势,1990s后萎缩速率加剧,水域面积总体缩小45.7%;湖区周边NDVI波动性增加显示了水退草进的变化趋势;湖面NDVI值的骤增,暗示叶绿素a或悬浮物浓度增加,间接表明湖区水质变差.根据40 a来水域面积变化,红碱淖的演变进程可依据湖泊面积动态度划分为稳定期(1973-1994年)和萎缩期(1994-2013年)两个阶段,气候暖干化是影响稳定期(1973-1994年)湖泊变化的主要因素,在萎缩期(1994-2013年)气候暖干化叠加人类活动是湖泊水量减少和水质变差的诱因,高强度的人为干扰如人工筑坝、灌溉耗水和煤炭开采是红碱淖水域面积锐减的主要原因.  相似文献   

7.
城市湖泊时空演变的遥感分析——以武汉市为例   总被引:8,自引:4,他引:4  
曾忠平  卢新海 《湖泊科学》2008,20(5):648-654
在GIS技术支持下,以多时相Landsat TM/ETM 遥感影像为数据源,利用NDWI指数提取武汉市1991年、1995年、2000年、2002年的湖泊水域面积信息,引入湖泊萎缩强度等测度指标,从时空两方面对主城区湖泊水域变化和26个主要湖泊的发展演化特征进行了系统地分析.结果表明,从1991年到2002年,武汉市主城区湖泊水域面积急剧减少,达38.67km2;主要湖泊水域的面积变化具有时空分异特征,湖泊萎缩变化与湖泊所处地区的城市建设水平、道路交通规划、政策导向等因素密切相关.分析认为,城市湖泊转变为公园后,萎缩现象得到明显遏制,通过立法加强湖泊保护力度并尽早制定城市湖泊的规划管理措施,是保护武汉幸存湖泊的有力手段.  相似文献   

8.
2000年以来青藏高原湖泊面积变化与气候要素的响应关系   总被引:1,自引:0,他引:1  
青藏高原星罗密布的湖泊对气候变化十分敏感,在自然界水循环和水平衡中发挥着重要作用.以MODIS MOD09A1和SRTM DEM为数据源,提取了2000-2016年青藏高原丰水期面积大于50 km2的湖泊边界,从内外流分区、湖泊主要补给来源和湖水矿化度三个方面对2000年以来湖泊面积变化进行分析,并结合青藏高原近36年气象数据,根据气象要素变化趋势分区,初步探讨青藏高原湖泊面积变化与气候要素的关系.结果表明:青藏高原面积大于50 km2的138个湖泊整体扩张趋势显著,总面积增加2340.67 km2,增长率为235.52 km2/a.其中,扩张型湖泊占67.39%,萎缩型湖泊占12.32%,稳定型湖泊占20.29%.内流湖扩张趋势显著,外流湖扩张趋势较明显;以冰雪融水为主要补给来源的湖泊整体扩张趋势明显,以地表径流和河流补给为主要补给源的湖泊也呈扩张趋势;盐湖和咸水湖以扩张为主,淡水湖的扩张、萎缩和稳定三种类型较均衡.在青藏高原气候暖湿化方向发展背景下,湖泊面积变化与气候要素具有显著的区域相关性.气温和降水变化趋势分区结果表明,气温增加、降水增加强趋势的高原Ⅰ区湖泊扩张程度(78.18%)依次大于气温降低、降水量呈增加趋势的Ⅴ区(66.67%),气温、降水量呈增加趋势的Ⅱ区(60.78%),气温呈降低、降水量呈增加强趋势的Ⅳ区(58.83%)和气温呈增加、降水量呈减少趋势的Ⅲ区(50.00%).湖泊面积变化对气候变化响应研究表明,升温引起的冰雪融水补给对Ⅰ区、Ⅱ区和Ⅲ区湖泊面积扩张的影响显著,加之降水量的增加,湖泊扩张速率明显;Ⅳ区和Ⅴ区湖泊面积扩张主要受降水量增加影响显著.整体而言,气温主要影响以冰雪融水为主要补给来源的湖泊,降水量主要影响以降水和地表径流为主要补给来源的湖泊.  相似文献   

9.
近50a中国东部和西部湖区水量空间变化趋势   总被引:3,自引:1,他引:2  
黄智华  薛滨  逄勇 《湖泊科学》2007,19(5):497-503
基于近代湖泊水量(水位、水深和面积)记录,对近50 a中国东部和西部湖区水量的空间变化趋势进行半定量研究.结果表明:整体上看,20世纪80年代以前,东部湖区水量呈现波动变化而西部湖区水量基本趋于减少;80年代以后东部和西部湖区水量均表现为不同程度的增加.其中,20世纪50年代,东部和西部湖区水量状况基本相似,均保持高湖面特征.20世纪60年代,东部湖区水量有所减少但幅度不明显,而西部湖区以中高水量为主.20世纪70年代,不同区域湖泊基本都保持中湖面特征.20世纪80年代,东部湖区水量呈增加趋势而西部湖区水量相对70年代略有减少.20世纪90年代,东部和西部湖区水量均表现为增加趋势.进一步研究发现,最近50 a来中国不同区域湖泊水量发生变化除了气候因素的驱动外,人类活动的影响也是不可忽视的.  相似文献   

10.
于2014年4、7和10月以及2015年1月(分别代表春、夏、秋和冬季)对鄱阳湖13个常规监测点表层水体中氧化亚氮(N_2O)浓度进行测定,并选择合适的模型估算其释放量.结果表明,鄱阳湖全年N_2O平均浓度为32.57±17.35 nmol/L,总体处于过饱和状态,平均饱和度为256.83%±129.05%.鄱阳湖N_2O年平均交换通量为0.83±0.69μmol/(m2·h).鄱阳湖水体N_2O季节性释放规律为春季最高,平均交换通量为1.71μmol/(m2·h),其次是夏季和冬季,秋季最低.从空间上来看,春季北部湖区交换通量显著高于南部湖区.相关性分析表明,铵态氮浓度是影响夏季和冬季鄱阳湖水体N_2O产生的主要因素.结合水域面积初步估算出全年释放N_2O约1.29×107mol,其中春季和夏季是鄱阳湖水体N_2O释放的高峰期,总释放量约占全年的80.40%.全年通过N_2O输出氮素约为361.93 t,对鄱阳湖流域内N_2O分布及质量平衡具有一定影响.  相似文献   

11.
Urbanization influences hydrologic cycle significantly on local, regional even global scale. With urbanization the water resources demand for dense population sharpened, thus it is a great challenge to ensure water supply for some metropolises such as Beijing. Urban area is traditionally considered as the area with lower evapotranspiration (ET) on account of the impervious surface and the lower wind speed. For most remote sensing models, the ET, defined as latent heat in energy budget, is estimated as the difference between net radiation and sensible heat. The sensible heat is generally higher in urban area due to the high surface temperature caused by heat island, therefore the latent heat (i.e. the ET) in urban area is lower than that in other region. We estimated water consumption from 2003 to 2012 in Beijing based on water balance method and found that the annual mean ET in urban area was about 654 mm. However, using Surface Energy Balance System (SEBS) model, the annual mean ET in urban area was only 348 mm. We attributed this inconsistence to the impact of anthropogenic heat and quantified this impact on the basis of the night-light maps. Therefore, a new model SEBS-Urban, coupling SEBS model and anthropogenic heat was developed to estimate the ET in urban area. The ET in urban area of Beijing estimated by SEBS-Urban showed a good agreement with the ET from water balance method. The findings from this study highlighted that anthropogenic heat should be included in the surface energy budget for a highly urbanized area.  相似文献   

12.
鄱阳湖地区大气边界层特征的数值模拟   总被引:5,自引:1,他引:4       下载免费PDF全文
采用WRFV2.2中尺度数值模式对鄱阳湖地区200 km×200 km范围内,2009年11月5日00∶00至2009年11月6日12∶00不同高度的气象要素进行了数值模拟,得到了水平分辨率为1 km的鄱阳湖地区大气边界层风、温、湿度场和廓线分布的大气边界层物理特征.模拟结果发现:白天鄱阳湖面上空存在着冷岛效应并伴随湖风,而夜间湖面上空存在着热岛效应并伴随陆风,湖面与陆地之间最大温差可达6 ℃;同时地形以及下垫面类型对鄱阳湖区风场的分布具有很大影响,夜间存在一条东北西南走向的低空辐合带,白天逐渐消失;此外受风场和地形作用湖面上空的湿度分布也不均匀,白天湿度层厚度低而夜晚湿度层厚度高,湖中心右侧湿度值大于左侧湿度值.模拟结果能较好地反映鄱阳湖的大气边界层物理特征,有助于了解鄱阳湖地区区域气候的特点,以及由于地形、地理环境、地表特征所形成的不同高度上的风、温、湿的分布规律和大气边界层物理特征,为鄱阳湖地区局地天气预报、风能资源开发、环境保护等提供了科学依据.  相似文献   

13.
鄱阳湖流域径流模型   总被引:8,自引:4,他引:4  
流域径流是鄱阳湖主要来水,建立鄱阳湖流域径流模型对揭示湖泊水量平衡及其受流域自然和人类活动的影响具有重要意义.针对鄱阳湖-流域系统的特点:流域面积大(16.22×104km2)、多条入湖河流、湖滨区坡面入湖径流等,研究了相应的模拟方法,建立了考虑流域土壤属性和土地利用空间变化的鄱阳湖流域分布式径流模型.采用6个水文站1991-2001年的实测河道径流对模型进行了率定和验证.结果显示,模型整体模拟精度较高.其中,赣江、信江和饶河均取得了较好的模拟结果,月效率系数为0.82-0.95;抚河和修水模拟精度略低,为0.65-0.78.模型揭示了研究时段内年平均入湖径流总量为1623×108m3,其中,赣江最多,占47%,其次为信江和抚河,分别占13%和12%,湖滨区坡面入湖径流约占4%,其余24%来自饶河、修水以及其它入湖支流.模型将用于评估流域下垫面或气候变化引起的入湖水量变化,为湖泊水量平衡计算提供依据.  相似文献   

14.
Monthly evapotranspiration (ET) rates (2000 to 2009) across Nebraska at about 1‐km resolution were obtained by linear transformations of the MODIS (MODerate resolution Imaging Spectroradiometer) daytime surface temperature values with the help of the Priestley–Taylor equation and the complementary relationship of evaporation. For positive values of the mean annual precipitation and ET differences, the mean annual net recharge was found by an additional multiplication of the power‐function‐transformed groundwater vulnerability DRASTIC‐code values. Statewide mean annual net recharge became about 29 mm (i.e., 5% of mean annual precipitation) with the largest recharge rates (in excess of 100 mm/year) found in the eastern Sand Hills and eastern Nebraska. Areas with the largest negative net recharge rates caused by declining groundwater levels due to large‐scale irrigation are found in the south‐western region of the state. Error bounds of the estimated values are within 10% to 15% of the corresponding precipitation rates and the estimated net recharge rates are sensitive to errors in the precipitation and ET values. This study largely confirms earlier base‐flow analysis‐based statewide groundwater recharge estimates when considerations are made for differences in the recharge definitions. The current approach not only provides better spatial resolution than available earlier studies for the region but also quantifies negative net recharge rates that become especially important in numerical modeling of shallow groundwater systems.  相似文献   

15.
As a critical water discharge term in basin‐scale water balance, accurate estimation of evapotranspiration (ET) is therefore important for sustainable water resources management. The understanding of the relationship between ET and groundwater storage change can improve our knowledge on the hydrological cycle in such regions with intensive agricultural land usage. Since the 1960s, the North China Plain (NCP) has experienced groundwater depletion because of overexploitation of groundwater for agriculture and urban development. Using meteorological data from 23 stations, the complementary relationship areal evapotranspiration model is evaluated against estimates of ET derived from regional water balance in the NCP during the period 1993–2008. The discrepancies between calculated ET and that derived by basin water balance indicate seasonal and interannual variations in model parameters. The monthly actual ET variations during the period from 1960 to 2008 are investigated by the calibrated model and then are used to derive groundwater storage change. The estimated actual ET is positively correlated with precipitation, and the general higher ET than precipitation indicates the contributions of groundwater irrigation to the total water supply. The long term decreasing trend in the actual ET can be explained by declining in precipitation, sunshine duration and wind speed. Over the past ~50 years, the calculated average annual water storage change, represented by the difference between actual ET and precipitation, was approximately 36 mm, or 4.8 km3; and the cumulative groundwater storage depletion was approximately 1700 mm, or 220 km3 in the NCP. The significantly groundwater storage depletion conversely affects the seasonal and interannual variations of ET. Irrigation especially during spring cause a marked increase in seasonal ET, whereas the rapid increasing of agricultural coverage over the NCP reduces the annual ET and is the primary control factor of the strong linear relationship between actual and potential ET. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

16.
本文分析了太湖湖区降水的时空分布特征;估算和讨论了太湖的降水效应;揭示了太湖的降水效应不太显著;太湖湖区降水量的变化是该湖水位变化的决定性因素之一。  相似文献   

17.
Accurate estimates of seasonal evapotranspiration (ET) at different temporal and spatial scales are essential for understanding the biological and environmental determinants of ecosystem water balance in arid regions and the patterns of water utilization by the vegetation. For this purpose, remote sensing ET estimates of a Patagonian desert in Southern Argentina were verified with field measurements of soil evaporation and plant transpiration using an open top chamber. Root distribution and seasonal variation in soil volumetric water content were also analysed. There was a high correlation between remote sensing and field measurements of ecosystem water fluxes. A substantial amount of the annual ET occurred in spring and early summer (73.4 mm) using winter rain stored in the soil profile and resulting in water content depletion of the upper soil layers. A smaller amount of annual ET was derived from few rainfall events occurring during the mid or late summer (41.4 mm). According to remote sensing, the 92.9% of the mean annual precipitation returns to the atmosphere by transpiration or evaporation from the bare soil and by canopy interception. Only 7.1% infiltrates to soil layers deeper than 200 cm contributing to the water table recharge. Fourier time series analysis, cross‐correlation methods and multiple linear regression models were used to analyse 11 years of remote sensing data to assess determinants of water fluxes. A linear model predicts well the variables that drive complex ecosystem processes such as ET. Leaf area index and air temperature were not linearly correlated to ET because of the multiple interaction among variables resulting in time lags with ET variations and thus these two variables were not included in the linear model. Soil water content, the fraction of photosynthetic active radiation and precipitation explained 86% of the ET monthly variations. The high volumetric water content and the small seasonal variations at 200‐cm depth were probably the result of little water uptake from deeper soil horizons by roots with low hydraulic conductivity. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

18.
鄱阳湖区域极端降水异常的特征及成因   总被引:3,自引:1,他引:2  
闵屾  刘健 《湖泊科学》2011,23(3):435-444
鄱阳湖区域4-7月极端降水总量存在显著的2-4年周期变化,且鄱阳湖4-7月极端降水总量的变率存在比较显著的增加趋势.1998年6月鄱阳湖地区极端降水总量异常偏多,而2001年6月异常偏少,仅为1998年的13.5%左右.通过分析大气环流场、水汽输送等方面的差异,发现1998年6月500hPa高度场上乌拉尔山高脊和鄂霍次...  相似文献   

19.
The Bosten Lake watershed investigated in this study has seen significant land cover and climate change. The spatiotemporal relationship between evapotranspiration (ET) and environmental factors remain unclear. In this study, trend analysis and correlation methods are applied to analyze the spatiotemporal characteristics of ET and the relationship between ET and its driving factors using remotely sensed ET data and measured climate data between 2001 and 2018. During the study period, high values of ET primarily occurr in the wetlands of the plain area and the mid‐elevation mountain areas. The ET values show a significantly increasing trend in the different vegetation types due to climate change and other factors. The ET change trend in the study area is in the range of ?13.4 to ≈35.9 mm per year; the desert area exhibits a significant decrease and most of the mountain areas show a significantly increasing trend. ET is significantly correlated with land surface temperature, normalized difference vegetation index (NDVI), and solar radiation. The dominant factor affecting ET is NDVI, accounting for 15.2% of the study area. The results of this study highlight the need for appropriate land‐use strategies for managing water resources in arid land ecosystems.  相似文献   

20.
Pinyon‐juniper (PJ) cover has increased up to 10‐fold in many parts of the western U.S. in the last 140+ years. The impacts of these changes on streamflows are unclear and may vary depending on the intra‐annual distribution and amount of precipitation. Given the importance of streamflow in the western U.S., it is important to understand how shifts in PJ woodland cover may produce changes in streamflow across the region's diverse hydroclimates. To this end, we simulated the land surface water balance with contrasting woodland and grassland cover with the Hydrologiska Byråns Vattenbalansavdelning (HBV) model at a 4‐km resolution across the distribution of PJ woodlands in the western U.S. We used shifts in evapotranspiration (ET) between woodland and grassland cover as a proxy for potential changes in streamflows. Comparison of HBV model results with paired catchment studies indicated the model reasonably simulated annual decreases in ET with changes from woodland to grassland cover. For the northern and western ecoregions of the PJ distribution in the western U.S. where precipitation predominantly occurs in the winter, HBV simulated a 25 mm (37%) annual decrease in ET with conversion to grassland from woodland. Conversely, in southern ecoregions of PJ distribution with prominent summer monsoons, annual differences in ET were only 6 mm (19%). Our results suggest that only 29% of the PJ distribution, compared to an estimated 45% based on precipitation amount alone, has the potential for meaningful increases in streamflow with land cover change from woodland to grassland.  相似文献   

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

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