首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
贡嘎山雅家埂峨眉冷杉林线种群的时空动态   总被引:1,自引:0,他引:1  
冉飞  梁一鸣  杨燕  杨阳  王根绪 《生态学报》2014,34(23):6872-6878
通过对贡嘎山雅家埂峨眉冷杉种群林线附近6个3000 m2样地(阴阳坡各3个)中峨眉冷杉(Abies fabri Craib)种群的定位调查,分析了过去100a间该区峨眉冷杉种群的时间-空间动态。结果表明:1)雅家埂林线附近峨眉冷杉种群密度在过去100 a(主要是近50 a)有显著的升高,但树线的海拔位置并无明显的爬升;2)阴阳坡林线格局存在显著的坡向分异:阴坡林线和树线的海拔高度显著高于阳坡(分别比阳坡高152.5 m和135.8 m),阳坡林线附近峨眉冷杉早期的生长速率在大于阴坡,但后期的生长速率却低于阴坡;3)热量(温度)控制假说不能完全解释雅家埂目前的树线格局,除气候因素之外,其它因素也限制了雅家梗地区树线位置的变化。  相似文献   

2.
树线是监测森林生态系统对气候变化响应的理想区域。树线如何响应和适应气候变化是高寒区森林生态系统与气候关系研究的重点内容之一。本文综述了北半球树线位置和更新对气候变化响应的研究进展。自20世纪以来,北半球不同空间尺度上树线生态过程对气候变化的响应并不一致。变暖导致全球大部分地区的树线树木更新呈现增加趋势。同时,变暖引发的干旱也触发了部分地区树线树木的生长衰退甚至死亡。欧洲和北美洲树线位置上移的比例大致相同,均为56%,其余44%保持稳定。亚洲地区的树线位置上移的比例高达64%,而36%的树线位置保持相对稳定状态。本文讨论了造成树线波动空间差异的主要驱动因素。研究揭示,不同地区的树线波动受到不同气候因素和非气候因素的调控。这一认识对了解树线演化、预测树线在未来气候变化条件下的时空变化趋势具有一定意义。  相似文献   

3.
树线交错带是具有强烈生物交互作用的高寒生态过渡带,生物互作与树线生态过程密切相关。本研究系统综述了气候变化下植物间、动植物间和微生物与植物间互作因子对树线生态过程的影响。植物间互利或竞争作用的强度调控变暖背景下树线生态过程的变化,目前尚缺少树轮生态学证据,且亟待检验高阶互作的适用性;动物采食活动、微生物与植物间互作可通过影响土壤状况、改变树木生长和更新等生态过程动态,增强或削弱树线与气候间耦合关系。迄今为止,地下与地上过程联系如何影响气候变暖下的树线动态尚不明晰,而营养级间互作可能调制树线生态过程对气候响应。青藏高原等高寒区具有开展此类研究的优势与潜力。  相似文献   

4.
为了了解青藏高原东北部不同降水梯度下, 高山林线处的树木径向生长与气候关系是否存在差异, 在青海东北部从西北到东南沿降水梯度设置3个高山林线采样点: 乌兰县哈里哈图国家森林公园(HL, 年降水量217 mm)、都兰县曲什岗(QS, 281 mm)和同德县河北林场(HB, 470 mm), 运用树轮年轮学方法分析林线优势种祁连圆柏(Sabina przewalskii)的径向生长-气候关系随降水梯度的变化规律。结果表明: 不同降水梯度下, 降水对祁连圆柏径向生长的限制作用差异不明显, 但温度对祁连圆柏径向生长的影响存在显著差异。在低降水区域(HL), 冬、夏季最低气温主要限制祁连圆柏径向生长, 并且在不同气候特征年中无明显变化; 在中降水区域(QS), 祁连圆柏的径向生长明显受冬季最低气温影响, 与低降水区域相比, 春、夏季最低气温对祁连圆柏径向生长的限制作用减弱, 并且主要限制因子在不同气候特征年存在显著变化; 在高降水区域(HB), 冬、夏季最低气温对祁连圆柏径向生长的限制作用不显著, 而春、秋季最低气温对祁连圆柏径向生长的抑制作用显著增加, 并且主要集中在高温年和干旱年。该研究结果并未支持干旱(湿润)区高山林线树木径向生长主要由水分(温度)限制的假说, 但是林线处降水量会影响树木生长与温度的关系。随着青藏高原东北部暖湿化加剧, 不同地区林线处树木生长的气候限制因子可能存在复杂化趋势。  相似文献   

5.
为了了解青藏高原东北部不同降水梯度下, 高山林线处的树木径向生长与气候关系是否存在差异, 在青海东北部从西北到东南沿降水梯度设置3个高山林线采样点: 乌兰县哈里哈图国家森林公园(HL, 年降水量217 mm)、都兰县曲什岗(QS, 281 mm)和同德县河北林场(HB, 470 mm), 运用树轮年轮学方法分析林线优势种祁连圆柏(Sabina przewalskii)的径向生长-气候关系随降水梯度的变化规律。结果表明: 不同降水梯度下, 降水对祁连圆柏径向生长的限制作用差异不明显, 但温度对祁连圆柏径向生长的影响存在显著差异。在低降水区域(HL), 冬、夏季最低气温主要限制祁连圆柏径向生长, 并且在不同气候特征年中无明显变化; 在中降水区域(QS), 祁连圆柏的径向生长明显受冬季最低气温影响, 与低降水区域相比, 春、夏季最低气温对祁连圆柏径向生长的限制作用减弱, 并且主要限制因子在不同气候特征年存在显著变化; 在高降水区域(HB), 冬、夏季最低气温对祁连圆柏径向生长的限制作用不显著, 而春、秋季最低气温对祁连圆柏径向生长的抑制作用显著增加, 并且主要集中在高温年和干旱年。该研究结果并未支持干旱(湿润)区高山林线树木径向生长主要由水分(温度)限制的假说, 但是林线处降水量会影响树木生长与温度的关系。随着青藏高原东北部暖湿化加剧, 不同地区林线处树木生长的气候限制因子可能存在复杂化趋势。  相似文献   

6.
为明晰坡向和海拔对山体土壤水热和植物分布格局的影响,以青藏高原东北缘山体高寒草甸为研究对象,利用回归分析、典型相关分析(CCA)排序和方差分解等方法,对阶地与同一山体不同坡向和海拔的189个样方内土壤温湿度和植物分布进行分析和定量分解.结果表明: 阴坡物种丰富度最高,阶地最低.随海拔升高,阴坡和山脊物种丰富度先增加后降低,而阳坡物种丰富度呈线性增加.阳坡土壤温度最高,不同坡向0~20 cm土层土壤温度随海拔升高基本不变;阴坡土壤湿度最高,不同坡向0~30 cm土层土壤湿度随海拔升高而增加.方差分解表明,在0~30 cm土层中,坡向和海拔共解释土壤温度变化的100%,土壤湿度变化的51.8%.坡向单独解释土壤温度变化的72.2%,贡献率最高;海拔单独解释土壤湿度变化的51.8%,贡献率最高.大多数植物倾向分布于中等海拔的阴坡与山脊之间.阴坡以莎草科为主,阳坡以禾本科为主,山脊为过渡地带.莎草科、禾本科和豆科主要分布于低海拔区.坡向和海拔共解释了山地植物多度变化的28.6%,坡向单独解释19.9%,贡献率最高.在山地高寒草甸生态系统,综合考虑小尺度地形造成的土壤及植物分布格局差异的基础上,在进行生产与生态恢复的分区管理时,应优先考虑坡向造成的土壤和植物差异.  相似文献   

7.
地形对亚热带山地景观尺度植被格局影响的梯度分析   总被引:35,自引:4,他引:35       下载免费PDF全文
用 7个指标反映三峡大老岭地区森林群落生境的地形特征 ,采用 DCCA排序方法定量分析不同地形因子对亚热带山地森林植被格局的影响及它们各自的生态意义 ;并分析它们的作用强度沿海拔的变化。结果表明 :1)DCCA排序的前 4个轴分别与海拔、坡位和坡度、坡向和坡面显著相关 ;山地森林植被存在多重尺度和方向的分异格局 ;2 )山地森林群落样方生境的热量、水分 (养分 )和光照条件可主要由海拔、坡位和坡度、坡向等地形特征分别反映 ;3)在垂直梯度上 :海拔的影响在海拔 6 0 0~ 80 0 m和 15 0 0~ 170 0 m最大 ;坡向影响的峰值在山顶带 ;坡度、坡位与坡形的影响在山地中部影响最大 ,其中坡度在海拔 6 0 0 m以下的影响也很显著 ;坡面的影响在山地中部最显著。  相似文献   

8.
以6种不同方式对林相图中同一样带取样,采用CA、DCA和CCA 3种排序方法,研究了取样方式对排序轴的解释效果、物种与环境因子、环境因子之间、环境因子与坐标轴之间关系的影响.结果表明,样方大小和形状的变化在不同程度上改变了排序结果.大样方和长方形样方都增强了排序轴的解释效果,并对双序图中稀有种、独特种的位置有较大的影响;环境因子中土壤因子对样方的大小和形状都很敏感,坡度、经纬度只对样方大小敏感,坡位、海拔、温度和降水则对样方形状敏感;随着样方面积的增加,海拔、温度和降水的作用降低,而坡向的作用增加.  相似文献   

9.
土壤胞外酶及其化学计量比是反映土壤养分可用性和微生物底物限制的敏感指标。然而, 对全球变化敏感的高山树线过渡带土壤酶活性和化学计量比的变化特征及其关键驱动因素仍不清楚。该研究在青藏高原东南部的川西贡嘎山高山树线过渡带(森林、树线、灌丛)进行土壤采样, 测定了树线过渡带土壤中5种水解酶(β-葡萄糖苷酶(BG)、纤维素二糖水解酶(CBH)、木聚糖水解酶(XYL)、N-乙酰氨基葡萄糖苷酶(NAG)、亮氨酸氨基肽酶(LAP))和2种氧化酶(多酚氧化酶(POX)、过氧化氢酶(CAT))的活性, 并计算土壤胞外酶活性化学计量比(碳、氮(N)酶活性比和碳质量指数)。结果表明: 灌丛土壤LAP、POX、CAT活性显著低于树线和森林土壤, XYL活性在树线最低, 其他胞外酶活性在树线过渡带不同位置差异不显著。灌丛土壤lnBG/lnLAP显著高于森林和树线处土壤, lnBG/ln(NAG + LAP)在树线过渡带没有显著变化, 碳质量指数在树线处最高。非度量多维尺度分析表明, 土壤有机碳、全氮、硝态氮含量和植物叶片木质素:N是影响树线过渡带土壤酶活性差异的主要因素, 植物叶片碳氮比、木质素:N和土壤可溶性氮含量是影响树线过渡带土壤胞外酶活性化学计量比差异的主要因素。综上所述, 贡嘎山地区的部分土壤酶活性及其化学计量比沿树线过渡带会发生明显的变化, 这种变化可能是由不同植物类型下微生物群落差异导致。这表明, 未来气候变化引起的树线迁移可能会改变胞外酶活性进而影响土壤养分循环。  相似文献   

10.
通过对库姆塔格沙漠南缘植被进行实地调查,对该地区植物物种丰富度及其群落组成进行研究,并选取11个影响物种丰富度的气候因子和地形因子,利用PCA分析、方差分解等方法探讨气候、地形因子对物种丰富度的影响。结果显示:该区共有植物15科32属38种,植物群落种类匮乏,物种组成单一,植物生活型主要以灌木、多年生草本为主,占所有物种的80%以上;植物物种丰富度与植物群落组成明显受到水热条件的制约。回归分析结果表明,该区物种丰富度与能量因子呈显著负相关(P<0.05),而与水分因子呈显著正相关(P<0.05)。方差分解结果表明,水分、能量共同制约了该地区的物种丰富度,二者的共同解释率为44.3%。此外,地形因子对研究区物种丰富度也存在一定影响,能进一步提高环境因子对物种丰富度的解释率。总之,库姆塔格南缘物种组成单一、物种丰富度格局受到水热条件的共同制约,同时地形的变化也有重要影响。  相似文献   

11.
通过对库姆塔格沙漠南缘植被进行实地调查,对该地区植物物种丰富度及其群落组成进行研究,并选取11个影响物种丰富度的气候因子和地形因子,利用PCA分析、方差分解等方法探讨气候、地形因子对物种丰富度的影响。结果显示:该区共有植物15科32属38种,植物群落种类匮乏,物种组成单一,植物生活型主要以灌木、多年生草本为主,占所有物种的80%以上;植物物种丰富度与植物群落组成明显受到水热条件的制约。回归分析结果表明,该区物种丰富度与能量因子呈显著负相关(P 0.05),而与水分因子呈显著正相关(P 0.05)。方差分解结果表明,水分、能量共同制约了该地区的物种丰富度,二者的共同解释率为44.3%。此外,地形因子对研究区物种丰富度也存在一定影响,能进一步提高环境因子对物种丰富度的解释率。总之,库姆塔格南缘物种组成单一、物种丰富度格局受到水热条件的共同制约,同时地形的变化也有重要影响。  相似文献   

12.
与传统方法相比, 利用树木年轮学方法研究树线过渡区树木生长温度敏感性高低的问题更注重比较树木个体间的生长情况, 从各个树轮序列间的生长一致性程度和树轮序列对气候因素(气温、降水)的响应一致性程度可探讨树线过渡区树木生长的温度敏感性。为了认识高山树线过渡区内树木生长的温度敏感性问题, 选择西藏昌都地区八宿县的一条川西云杉树线过渡区为研究对象, 比较了过渡区内树木个体间的生长一致性, 分析了树木生长与气候因素的相关性及其在个体间的异同。结果显示: 树线过渡区内树轮生长在个体间的一致性较低, 树轮生长与气温的关系在树木个体间的一致性也较低, 而树轮生长与当年4-9月降水的关系相对较强。西藏八宿树线过渡区属于干旱区, 相对于气温而言, 降水对树木生长的影响更大。此外, 小生境的异质性及干扰事件的发生也有可能降低树木对温度的敏感性。在全球变暖及极端气候事件增加的背景下, 树木生长的温度敏感性被高估可能会导致对树线过渡区位置及树线过渡区内群落生产力等的预测产生偏差, 这一问题应该在区域生态模拟研究和相关林业经营与管理上得到重视。  相似文献   

13.
 与传统方法相比, 利用树木年轮学方法研究树线过渡区树木生长温度敏感性高低的问题更注重比较树木个体间的生长情况, 从各个树轮序列间的生长一致性程度和树轮序列对气候因素(气温、降水)的响应一致性程度可探讨树线过渡区树木生长的温度敏感性。为了认识高山树线过渡区内树木生长的温度敏感性问题, 选择西藏昌都地区八宿县的一条川西云杉树线过渡区为研究对象, 比较了过渡区内树木个体间的生长一致性, 分析了树木生长与气候因素的相关性及其在个体间的异同。结果显示: 树线过渡区内树轮生长在个体间的一致性较低, 树轮生长与气温的关系在树木个体间的一致性也较低, 而树轮生长与当年4-9月降水的关系相对较强。西藏八宿树线过渡区属于干旱区, 相对于气温而言, 降水对树木生长的影响更大。此外, 小生境的异质性及干扰事件的发生也有可能降低树木对温度的敏感性。在全球变暖及极端气候事件增加的背景下, 树木生长的温度敏感性被高估可能会导致对树线过渡区位置及树线过渡区内群落生产力等的预测产生偏差, 这一问题应该在区域生态模拟研究和相关林业经营与管理上得到重视。  相似文献   

14.
张谧  熊高明  陈志刚  樊大勇  谢宗强 《生态学报》2004,24(12):2686-2692
地貌单元和地形的局部变化严重制约着光照、温度、水分甚至土壤资源的再分配 ,从而对植被的分布产生影响。以永久性监测样地为基础 ,利用数字高程模型模拟 0 .96 hm2样地的地形表面 ,在地理信息系统的支持下分析了米心水青冈 -曼青冈群落的地形异质性状况 ,并结合群落调查的结果探讨了地形因素对群落组成的影响。结果表明研究范围内样方的平均坡度在 18.2~ 74 .6°之间。坡向可以分为南坡、西南坡及东南坡 3部分 ,坡度及坡向差异显著。对坡度、坡向等地形因素及群落的物种组成进行半方差分析 ,发现三者的空间异质性尺度分别为 2 85 .8m,2 80 .5 m及 6 .1m。地形因素的变化尺度远大于植被组成的变化 ,说明小尺度下地形因素不是影响群落内种群分布格局的唯一因素。坡度和坡向因子的块金值分别为 2 1.8及 4 98.0 ,说明两者在小尺度下的某些过程不能忽视。植物与地形因素的相关分析表明不同的种对地形的适应性不同 ,群落中的 2 4种主要植物可以分为 4个生态种组。因此 ,群落的物种组成变化受到多种因素的影响 ,这种影响可以在地形变化中得到反映。  相似文献   

15.
近40年青藏高原植被动态变化对水热条件的响应   总被引:4,自引:0,他引:4  
植被对水热条件的响应因生态系统的空间异质性而具有显著的差异。基于GIMMS NDVI3g和MODIS NDVI逐旬数据集,通过数据融合构建的青藏高原1982—2020年的植被时间序列,利用Mann-Kendall趋势法分析近40年植被动态变化及其对温度、降水和辐射等水热条件的响应,并划定了植被动态的主要水热驱动因子分区。结果表明:(1)近40年青藏高原植被生长季平均NDVI呈现显著上升趋势,增速为0.006/10a,植被NDVI显著增加和减少的区域分别占青藏高原总面积的73.97%和18.38%;(2)青藏高原植被对水热条件的响应在静态上表现为高原腹地较高原边缘更加明显;在动态上表现为不同植被类型区对水热因子的响应关系、方向、程度均有所不同;整体上除森林和灌丛外,所有高寒植被类型与降水的响应程度要优于温度和辐射;(3)青藏高原植被生长受水热因子驱动的区域占高原总面积的55.95%,其中42.72%以上的区域气温、降水和太阳辐射的驱动作用是互补的,13.23%的区域由多个水热因子联合驱动;44.05%的区域为非气候驱动区。  相似文献   

16.
地形因子对森林景观格局多尺度效应分析   总被引:4,自引:3,他引:4  
以TM影像和野外调查为数据源,3S技术为研究手段,分析了泰山的森林景观结构动态变化以及景观格局沿高程的分异。选择相对高差较大的区域设计了8个地形方位,采用3种由连续样方组成的辐射状样带,在遥感影像上进行信息采集,研究地形对森林景观格局的多尺度效应,分析了不同地形因子对森林景观分异的影响。结果表明,在景观尺度上,地形方位、海拔、山地类型是影响森林景观镶嵌格局的控制因素,坡向、坡度是重要因素,坡形、坡位是不明显因素。以TM影像为信息源,从森林景观分异和梯度分析上,首先要考虑海拔坡向指数、坡形坡位指数和海拔,并立足于地形方位。地形主要因子间存在着稳定的显著正相关,为森林分异多尺度格局提供了较强的综合解释能力。  相似文献   

17.
印度特莱东部林地-草地系统中非生物和空间变量对木本和草本物种丰度的影响 目前尚不清楚哪些环境因素决定了热带稀树草原特别是在潮湿地带的林地和草地镶嵌处的林地和草地的物种多度。基于此,本研究探究了非生物和空间变量对印度东北部的台拉河生态系统木本和草本物种分布的影响,评估了气候和非气候因素在整个景观中保持可变的树草比和空间连通性和分散性的相对重要性。在519 km2的受保护的特莱栖息地中随机建立了134个30 m × 30 m的抽样样方,并调查了每个样方的木本和草本植物的物种多度和气候,以及非气候环境因素。基于不同的地点空间连通性模型,通过构建变量检验气候和非气候环境因素对物种多度的影响。使用冗余分析和方差分解定量解析环境变量和空间结构对林地和草地物种多度的相对重要性。研究结果表明,降雨、火灾、水分胁迫、地形和土壤养分在内的环境变量对物种多度和林草比有显著的影响。空间结构显著,最佳空间模型为反距离加权模型(inverse distance-weighted model), 而且显示最大的空间扩散距离可以达到23.5 km,表明扩散限制较弱。约21%的物种多度变化能够被环境和空间因素解释。这些结果揭示了植物群落动态的决定因素,即环境因子的时空变化可能驱动物种分布和多度的随机性,并对植被镶嵌产生主导影响。  相似文献   

18.
为探究广西金钟山细叶云南松径向生长及其对气候的响应,该文建立了细叶云南松树木年轮标准年表,采用响应分析探讨了径向生长与气候因子的关系,并使用逐步回归和方差分解量化了不同气候因子对径向生长的影响。结果表明:(1)响应分析结果显示,上一年1月、8月、9月、11月的日照时数以及上一年10月的平均最低气温和平均气温与标准年表年轮宽度均呈显著正相关,而上一年6月降雨量>10 mm的天数与年轮宽度呈显著负相关; 当年2月的降雨量、3月的空气相对湿度和平均最低气温、6月的日照时数以及9月的平均气温和平均最高气温均与年轮宽度呈显著正相关。(2)逐步回归最终模型的方差分解结果显示,上一年10月的平均最低气温对细叶云南松径向生长的影响最大(方差解释量达23.35%),其次是当年9月的平均最高气温(方差解释量为10.39%); 上一年1月和11月的日照时数分别解释了径向生长变异的3.94%和6.58%。综上表明,秋季的温度和冬季的光照条件是限制细叶云南松径向生长最主要的气候因子,早春干旱和雨季大降雨量会降低细叶云南松的径向生长。  相似文献   

19.
郑娇  李东  袁旭东  赵小祥  刘峰  田秋香 《生态学报》2023,43(21):8704-8715
树线过渡带作为高山地区重要的生态过渡带之一,是响应温度变化的敏感区域。树线过渡带内土壤碳储量丰富,其碳周转在全球碳循环方面扮演着重要角色。探究树线过渡带土壤有机碳矿化及其温度敏感性,对于预测气候变化背景下高山地区土壤碳循环过程具有重要的指导意义。为此以青藏高原东南部贡嘎山树线过渡带(森林、树线、灌丛)的土壤为对象,在室内开展90 d不同温度(15℃和20℃)的培养实验,测定土壤有机碳矿化速率,计算单位土壤有机碳累积矿化量、温度敏感性,并分析影响它们的相关因素。结果表明:土壤有机碳矿化速率受温度和样地类型的显著影响。升温显著增加土壤有机碳矿化速率,而不同样地类型间矿化速率差异显著,矿化速率大小表现为森林>树线>灌丛。本研究用单位土壤有机碳累积矿化量表征土壤有机碳的稳定性,经90 d的培养,15℃下树线过渡带从森林、树线到灌丛单位土壤有机碳累积矿化量分别为12.33 mg/g、12.99 mg/g和10.53 mg/g, 20℃下则分别为19.16 mg/g、21.14 mg/g和16.26 mg/g,灌丛土壤单位土壤有机碳累积矿化量显著低于森林和树线土壤,这表明灌丛土壤具备更...  相似文献   

20.
山地森林样带植被-环境关系的多尺度研究   总被引:28,自引:6,他引:28  
沈泽昊 《生态学报》2002,22(4):461-470
为了评价和比较地形和其它“直接”环境因子对植被变化影响的大小 ,通过调查一条亚热带山地森林样带及其生境因子 ,分析不同层次植被的水平格局 ;并利用多元回归和 DCCA方法比较了诸生境因子对植被格局变化的影响及其与分析尺度相关的变化。结果表明 :( 1 )影响乔木、灌木和草本层植被格局变化的首要因子都是地形特征 ;最重要的土壤因子是土壤发生层厚度。在对林下物种构成的影响方面 ,坡向和坡度的作用大于林冠空隙度。 ( 2 )除坡向外 ,其它地形因子之间存在不同程度的相关性 ;林冠孔隙度和土壤理化属性分别与不同的地形特征显著相关 ,显示了地形对林隙干扰和土壤侵蚀过程的控制。乔、灌、草的覆盖度是生境条件的综合反映 ;它们与非生物因子的相关性显示了成熟林下灌木和草本层物种生态习性的整体差异。 ( 3)各分析尺度上诸生境因子间的相关关系基本稳定。但随尺度增大 ,生境的小尺度异质性被平滑掉。诸因子之间的相关显著性程度逐渐降低。 ( 4 )从 5× 5 m2到 5× 4 0 m2 ,立地因子对乔、灌、草 3层植被变化的解释量随尺度的倍增而接近成倍增长 ;但其主要原因是样方数量增加 ;分析尺度变化的影响小得多。 ( 5 )地形与土壤、光照等直接生境因子对植被格局解释量之和相当。但在不同尺度上 ,地形特征都是  相似文献   

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

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