首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 140 毫秒
1.
辽宁省人工林樟子松径向生长对水热梯度变化的响应   总被引:3,自引:0,他引:3  
基于树轮年代学的理论和方法,建立辽宁地区人工林樟子松年轮宽度年表网络,研究人工林樟子松树木年轮生长变化,探讨樟子松生长与气候变化的关系以及其对水热梯度变化的响应规律。结果表明,研究区各采样点人工林樟子松的径向生长变化呈现较为一致的显著下降趋势(P0.05),并且下降趋势由西北至东南降低。空间上水平上,研究区樟子松的径向生长与研究区降水和相对湿度的空间梯度变化吻合(P0.05),均呈西北至东南增加的特点;各采样点樟子松年轮宽度变化的平均敏感度与降水量、相对湿度显著负相关关系(P0.05),但与温度的正相关关系不显著。各样点樟子松径向生长与月降雨量和月平均相对湿度多呈正相关关系,与月平均温度多呈负相关关系,进一步表明区域水分因子对人工林樟子松生长的限制作用明显。  相似文献   

2.
河北坝上樟子松人工林径向生长及其对气候因素的响应   总被引:2,自引:0,他引:2  
利用树木年轮学方法测定了河北坝上地区樟子松逐年树轮宽度和胸高断面积增量,建立了樟子松树轮宽度差值年表并分析其对气候因素的响应以及生态弹性,为河北坝上地区利用樟子松开展退化杨树防护林更新建设的适宜性提供理论参考。研究结果表明,河北坝上地区樟子松树龄在28 a左右,其生长阶段可划分为快速生长期(0—10 a),生长下降期(11—20 a)和生长平稳期(21—28 a);樟子松树轮宽度在1992—2000年期间为快速增长期((4.49±1.15)mm/a),在2001—2010年期间呈现下降趋势((3.62±1.39)mm/a),而在2011—2019年期间樟子松生长呈平稳特征((2.21±0.68)mm/a),约为快速增长期增长速度的一半;樟子松BAI在1992—2000年间一直呈上升趋势((5.09±2.26)cm~2/a),2001—2019年期间呈平稳特征((10.46±0.67)cm~2/a),表现出稳定且持续的生长能力。樟子松差值年表与气候因素的相关性显示:樟子松径向生长主要与上一年生长季末期(9、10月)、休眠期(12月)及当年生长季(7月)平均温度和最高温度呈显著负相关,其中...  相似文献   

3.
大兴安岭是我国气候变化最为显著的地区之一,兴安落叶松和樟子松是该地区最为重要的树种,研究它们径向生长对气候变化的响应差异,可以为预测气候变化下我国北方森林动态提供科学依据。在大兴安岭地区选择6个样点共采集兴安落叶松树轮和樟子松树轮样芯451个,建立了12个标准年表。比较了1900年以来树木径向生长趋势,利用Pearson相关分析法分析各样点兴安落叶松和樟子松生长对气候因子的响应,运用线性混合模型探讨温度和降水对兴安落叶松和樟子松年径向生长的影响,通过滑动相关对比两个树种生长-气候关系的时间稳定性。结果表明: 兴安落叶松径向生长与3月平均温度呈负相关,与上一年冬季和当年7月降水呈正相关。樟子松径向生长与当年8月温度呈正相关,与当年生长季(5—9月)降水呈正相关。冬季降雪对兴安落叶松径向生长起到重要的促进作用,夏季过多降水对樟子松径向生长起到显著的限制作用。兴安落叶松和樟子松生长对气候变化的响应存在明显差异,因此,气候变化可能会影响北方森林生态系统的树木生长、物种组成以及空间分布等。  相似文献   

4.
雷帅  张劲松  孟平  何子兵  程顺  孙守家 《生态学报》2020,40(13):4479-4492
通过树木年代学方法,测定了毛乌素和塞罕坝相同密度樟子松(Pinus sylvestris var.mongolica)全木(Ringwood,RW)、早材(Earlywood,EW)和晚材(Latewood,LW)宽度,计算胸高断面积增量(Basal area increment,BAI),并建立了樟子松年轮宽度年表,分析其对气候响应的差异。结果显示,毛乌素(Mu Us)樟子松轮宽随树龄呈"升-降"的曲线变化,塞罕坝(Saihanba)呈线性下降,两地樟子松BAI变化相似,呈"升-降"的曲线变化,EW占RW的65%—70%,表明EW对RW贡献较大。生长期间,毛乌素樟子松早晚材比例保持平稳,塞罕坝EW/RW值下降,LW/RW值上升,两地干旱事件均使LW/RW值下降,EW/RW值上升。差值年表(Residual chronology,RES)相关性分析显示,毛乌素樟子松径向生长主要与4、7月平均降雨,7月标准化降水蒸散发指数(Standardized precipitation evapotranspiration index,SPEI),3、8月平均温度及上年12月和当年3月最低温度呈正...  相似文献   

5.
以秦岭北麓南五台油松为样本,建立油松树轮宽度标准化年表(STD),研究油松径向生长与气候因子之间的相关性。结果表明: 秦岭北麓油松径向生长与前一年9月及当年5月水分因子呈显著正相关,与前一年11月温度因子呈显著正相关,与前一年10月、当年5月温度因子呈显著负相关。油松径向生长对不同气候因子响应均存在明显的滞后效应。油松径向生长与PDSI干旱指数具有较好的相关性,特别是与前一年9—12月、当年5月PDSI呈显著正相关。回归模型能较好地模拟树轮宽度指数与PDSI之间的关系,油松极宽窄轮的形成主要是干旱作用的结果。综合各种气候指标,PDSI能更好地反映研究区油松径向生长的特征。  相似文献   

6.
凉水自然保护区不同皮型红松径向生长对气候的响应   总被引:2,自引:2,他引:0  
运用相关函数及单年分析等树木年轮气候学方法,研究了黑龙江凉水国家自然保护区不同皮型红松径向生长与气候因子的关系、主要影响因子及这种响应关系是否长期稳定.结果表明: 细皮红松更适合做树木年轮气候学分析.两种皮型红松的径向生长对环境变化都比较敏感,对气候因子的响应无显著差异.1902—2009年生长季,尤其当年6月的气候因子是影响研究区两种皮型红松生长的主要因子.其中,温度表现为显著负相关,降水表现为显著正相关.不同时间段内红松径向生长对气候因子的响应存在差异,随着1970年后气温的快速升高及干旱的加剧,两种皮型红松径向生长对气候因子的响应更敏感,尤其表现为对生长季温度和更多季节水热复合因子变化(帕尔默干旱指数)的响应更显著.
  相似文献   

7.
根据兴安落叶松在大兴安岭山脉的分布特征,沿纬度设置了9个采样点,分析了兴安落叶松径向生长对气候因子响应的空间差异和时间动态。结果表明:总体上所有采样点的兴安落叶松径向生长与夏季(6—8月)标准化降水蒸散指数(SPEI)及降水量、2月SPEI及降水量呈正相关,与3月温度呈负相关。在空间上,南部年均温较高区域的兴安落叶松径向生长与2月SPEI呈显著正相关,在北部年均温较低的区域与3月温度呈显著负相关。在时间上,兴安落叶松生长-气候关系不稳定。随着温度升高,在年均温较高区域,夏季SPEI和降水量对树木生长的正效应以及夏季温度的负响应显著增强;在年均温较低区域,树木生长与3月温度负响应增强更明显。这表明气候变化改变了树木生长-气候关系,而且存在明显空间差异。在未来气候暖干化的背景下,大兴安岭兴安落叶松径向生长会受到抑制,在年均温较高且降水量较少的南部地区会因夏季水分亏缺和冬季干旱而衰退,温度较低的北部地区则可能因冬季干旱和冬季变暖而使生长受到抑制。  相似文献   

8.
王婷  李聪  张弘  任思远  李鹿鑫  潘娜  袁志良  叶永忠 《生态学报》2016,36(17):5324-5332
树木年轮资料能够提供区域内过去长时间的环境和气候信息,成为获取过去气候变化信息的重要手段之一。利用采自宝天曼自然保护区的两种针叶树种油松和华山松树木年轮样本分别建立了油松和华山松树轮宽度标准年表PT和PA,并将油松和华山松样本合并建立了联合树种的区域年表(RC)。3种年表分别与不同气候要素(月平均气温、月平均最高气温、月平均最低气温、月降水量)及其不同月份组合进行相关分析,结果表明,油松年表PT和华山松年表PA都包含较高的气候信息,且都和生长季不同月份温度显著负相关和降水显著正相关。其中,油松和华山松都与当年4月和5月降水显著正相关,油松还与当年5月的平均最高温度和上年11月最低温显著负相关,与当年3月平均最低气温显著正相关;华山松与上年10月和当年4月的平均最高温度显著负相关,与上年12月和当年7月平均最低温度显著正相关;联合年表RC包含了单个年表PT和PA共同的气候信息,与当年4—5月降水和3月最低温显著正相关,与当年7月最低温和4、5月最高温显著负相关。不同年表与生长季(3—8月)内气候要素月份组合的相关分析也表明联合年表RC包含和单物种年表PT、PA相似的气候信息并加强了PT、PA受当年生长季气候变化影响的公共信号。可见,同一地区的多树种联合年表一定程度上能体现出区域性树木生长对气候变化响应的生态生理特征,为同地区建成多树种联合年表来探讨当地区域性气候变化提供了可行性的理论基础和一定的参考作用.  相似文献   

9.
研究人工林径向生长与气候变化的关系对全球气候变暖背景下人工林合理经营有着重要的意义。该文对在辽东山区广泛栽培的黄花落叶松(Larix olgensis)和日本落叶松(Larix kaempferi)人工林, 运用树木年轮气候学方法建立了辽宁草河口和湾甸子林场落叶松人工林年表, 分析了落叶松径向生长对气候变化的响应以及气候条件、树种、立地条件和林分因子(林龄、密度、蓄积量等)的相对影响程度。结果发现在影响年轮-气候关系的因素中, 气象因子的潜在蒸发散(PET)的影响力最大; 林龄、密度和蓄积量同时也具有重要的影响作用。中龄落叶松人工林径向生长主要与气温呈正相关关系, 成熟落叶松人工林径向生长主要与气温呈负相关关系; 而其他因素, 如树种、立地条件等的影响作用不大。这表明在气候变暖背景下随着林龄增加, 林分会逐渐受到气温升高导致的水分亏缺的限制, 导致明显的生长下降趋势, 因而气候变暖对成熟落叶松人工林威胁更为严重, 所以要注重对成熟林的优先保护, 同时可以预测, 随着东北地区今后气候进一步变暖, 可能将逐步影响到林龄更小的林分的生长, 因此需要进一步研究如何在落叶松人工林经营中采取科学的措施来更好地应对未来气候变化。  相似文献   

10.
为研究滇西北高原树木径向生长与气候关系随海拔的变化规律, 分别在玉龙雪山低、中、高海拔采集丽江云杉(Picea likiangensis)年轮样本, 建立了不同海拔丽江云杉树轮宽度残差年表, 将年轮指数与气候因子进行响应分析、冗余分析以及滑动响应分析。结果表明: 玉龙雪山丽江云杉径向生长受气温和降水共同影响, 但不同海拔径向生长响应模式存在差异。其中当年1-3月降水与不同海拔丽江云杉径向生长均呈显著正相关关系; 当年生长季后期降水与中、低海拔树木生长呈显著负相关关系, 与高海拔树木生长呈显著正相关关系; 中、低海拔树木生长还受当年春季干旱胁迫; 而当年7月气温升高促进高海拔丽江云杉生长。冗余分析与响应分析结果基本一致, 说明冗余分析能够有效量化树轮宽度指数与气候因子的关系。滑动响应分析显示气温和降水在小时间尺度上的变化也会影响树木生长。结合不同海拔丽江云杉生长对气候因子的响应模式及未来气候预测, 玉龙雪山高海拔丽江云杉生长将得到加强, 而中、低海拔丽江云杉生长则表现出不确定性。  相似文献   

11.
曹恭祥  王云霓  季蒙  郭中  李艳慧  耿庆春  包瑞  刘佳 《生态学报》2021,41(13):5352-5359
以呼伦贝尔沙地樟子松人工林为研究对象,利用树木年轮学和解析木法,推算过去41年樟子松林的生物量和年净初级生产力,并分析年净初级生产力与气温、降水、湿度等气象因子的关系。研究表明:随着林龄的增加,樟子松林生物量变化符合逻辑斯蒂生长方程;樟子松林生物量从1977年的4.83 t/hm2到2017年的167.95 t/hm2。樟子松林净初级生产力年际差异较大,呈先增加后减小、最后趋于稳定的趋势,多年平均净初级生产力为4.08 t hm-2 a-1,最高达到6.13 t hm-2 a-1,最低为2.63 t hm-2 a-1。气候因子与樟子松林净初级生产力关系较为密切,并具有一定的"滞后效应";樟子松林净初级生产力与上一年8月、12月和当年3月的降水量显著相关;净初级生产力与上一年和当年8-9月的月均湿度极显著正相关;净初级生产力与各月温度均呈不同程度的负相关,尤其是上一年和当年的6-9月温度显著影响着净初级生产力。研究表明,樟子松林净初级生产力受温度和降水的综合影响,净初级生产力与气候因子的关系属于温度敏感型,未来该地区持续升温的条件下樟子松林净初级生产力可能会降低。  相似文献   

12.
外生菌根菌对樟子松苗木生长的影响   总被引:1,自引:0,他引:1  
张文泉  闫伟 《西北植物学报》2013,33(5):998-1003
以内蒙古大兴安林樟子松分布区较典型的6个乡土菌种与本实验室从法国引进的宿主范围较广的2个菌种及内蒙古地区分布广泛的土生空团菌为研究对象,采用固体和液体两种菌剂类型,对出苗30d左右的樟子松实生苗进行外生菌根菌接种试验,接种130d后观察菌根侵染情况及生长特性,同时对菌根化的樟子松苗木进行自然干旱胁迫,研究外生菌根菌接种对樟子松苗木生长及抗逆性的影响。结果显示:(1)5个乡土菌种和1个法国引进菌种均能够与樟子松实生苗共生菌根,菌根化苗木较对照的苗高增长10.5%~111.4%、地茎增长4.55%~92.8%、苗木干重增加6.8%~184.5%,且固体菌剂效果优于液体菌剂。(2)在干旱胁迫下,菌根化苗木叶绿素含量较对照增加5.82%~36.92%、净光合速率提高120.2%~217.03%,叶片推迟10.3~25h出现萎蔫,临界致死时间推迟13.8~38.2h。(3)接种菌种厚环乳牛肝菌和点柄乳牛肝菌的苗木菌根侵染率、苗木生长量、干旱胁迫下各项生理指标均显著优于其他菌种和对照。研究表明,樟子松实生苗菌根化对苗木生长有明显的促进作用,能够提高苗木对干旱胁迫的抵御能力;厚环乳牛肝菌和点柄乳牛肝菌是有潜在研究和应用价值的优良菌株,这两个菌株的固体菌剂可用于樟子松菌苗生产。  相似文献   

13.
吕刚  王婷  李叶鑫  魏忠平  王凯 《生态学报》2017,37(24):8294-8303
以科尔沁沙地南缘的彰武县章古台万亩林为研究对象,野外取样调查和实验室检测分析相结合,以樟子松固沙林采伐迹地为对照,研究不同植被恢复类型下樟子松固沙林更新迹地生境的改良效果,探讨植被重建后植物多样性及土壤理化性质的响应。结果表明:1)樟子松(Pinus sylvestris)固沙林更新迹地植被重建后,草本植物物种多样性增加,林下植物偶见种数目增多,13种草本植物在8个样地中未重复出现;2)与采伐迹地相比,不同植被重建类型土壤理化性质均有所改善,土壤容重、田间持水量、有机质、全氮、碱解氮、速效钾在彰武小钻杨(Populus xiaozhuanica)、五角枫(Acer mono)林地改良效果较好,全钾在红刺榆林(Hemiptelea davidii)地含量明显提高、山杏(Armeniaca sibirica)林地的有效磷相比采伐迹地明显提高,且均表现为上层改良效果优于下层;3)不同植被恢复类型土壤理化性质间存在显著差异,随着物种多样性的增加,土壤理化性质逐渐改善,土壤容重、田间持水量、有机质、全磷与植物多样性具有显著的相关性,土壤理化性质的与植物多样性相互作用,共同促进生态系统正向演替。研究结果为控制科尔沁沙地土壤沙化,加速该区生态系统的恢复与重建提供理论依据。  相似文献   

14.
沙质草地营造樟子松林后土壤容重的变化及其影响因子   总被引:2,自引:0,他引:2  
采用野外调查和室内试验相结合,以辽宁省章古台地区不同生长阶段(包括幼龄林、中龄林、成熟林和过熟林)的20块樟子松人工固沙林样地(以临近的7块天然草地为对照)为研究对象,研究了沙质草地营造樟子松人工林后不同生长阶段0—100 cm土层土壤容重的变化及其影响因子。结果表明:天然沙质草地营造樟子松人工固沙林后,不同生长阶段樟子松林在0—10 cm土层土壤容重变化的变异系数为78%,其他土层变异系数范围为1.08%—4.35%。随着樟子松人工林林龄的增加,土壤容重变化量在0—20 cm和60—100 cm层逐渐降低,在20—60 cm层先降低,到37年左右后逐渐升高,过熟林较成熟林显著增大。林龄对不同土层容重变化的决定系数由大到小依次为40—60、60—80、20—40、10—20、0—10、80—100 cm层。土壤容重变化在60—80 cm层与土壤粗颗粒(粒径0.05 mm)含量、在0—10、20—40 cm和60—80 cm层与土壤全氮含量、在0—10、20—60 cm和80—100 cm层与土壤全磷含量、在20—40 cm和80—100 cm层与土壤全钾含量显著负相关,且土壤全氮和全磷含量对土壤容重的影响效果随土层深度的增加逐渐降低,土壤容重变化在10—20 cm层与土壤含水率、在20—40 cm层与土壤有机碳含量呈显著的正相关。总体上,沙质草地营造樟子松人工林可以改善土壤结构,提高土壤质量,建议采取封育禁牧等营林措施增加樟子松林下枯落物积累,提高土壤养分含量,同时对37年樟子松人工林逐渐进行更新。  相似文献   

15.
樟子松人工林营建对土壤颗粒组成变化的影响   总被引:1,自引:0,他引:1  
植被恢复是退化生态系统的主要恢复措施,也是人类改善区域生态环境较为重要和直接的活动。目前,针对不同植被恢复方式对干旱半干旱地区土壤理化性质及生物特征开展了大量研究。然而,关于科尔沁沙地樟子松人工林营建对土壤颗粒组成变化的影响却鲜有报道。因此,以辽宁省章古台地区不同生长阶段(包括幼龄林、中龄林、成熟林和过熟林)的20块樟子松人工林样地为研究对象(以临近的7块天然草地为对照),研究了沙地樟子松人工林营建对0—100 cm土层土壤颗粒组成变化的影响。结果表明:沙质草地营建樟子松人工林后,不同土层土壤细颗粒(<0.05 mm)含量均呈增加趋势,并且在0—10 cm层增加趋势明显,随土层深度增加土壤细颗粒增加量逐渐降低(除幼龄林外),但樟子松林地土壤颗粒组成仍以砂粒为主,土壤粘粒和粉粒含量极低(仅占5%左右)。随着樟子松人工林林龄的增加,土壤细颗粒变化量在0—10 cm层逐渐升高,而在10—100 cm层并无显著变化趋势。土壤细颗粒含量的变化在10—100 cm层与土壤含水量呈显著正相关,在0—10、20—40 cm和80—100 cm层与土壤全钾极显著负相关,在20—60 cm层与土壤有...  相似文献   

16.
以3龄樟子松幼树为材料,2013年在科尔沁沙地研究了不同沙埋深度下其株高、叶片膜透性、渗透调节物质含量及保护酶活性变化,以揭示沙埋条件下樟子松幼树生长及其对逆境的生理响应特征。结果显示:(1)在沙埋深度低于株高以上2cm时被埋樟子松幼树能够正常生长,其株高和芽长均明显高于非沙埋对照,并以沙埋深度为株高的50%时增长幅度最大;当沙埋深度超过株高2cm以上时,虽然植株高度和芽长也较埋前有一定增长,但均低于对照,且所有处理植株均未破土,后来全部死亡。(2)所有沙埋处理的叶片可溶性糖含量均显著低于对照,而POD活性显著高于对照,可溶性蛋白质和脯氨酸含量也高于对照。(3)随沙埋深度增加,叶片相对含水量总体呈增加趋势,但大多数处理与对照差异不显著;丙二醛含量基本呈显著下降趋势,可溶性蛋白和脯氨酸含量先增加后下降,而大多数处理的膜透性与对照差异不显著;随着沙埋深度增加,叶片可溶性糖含量显著下降,SOD和POD活性均先增加后下降。(4)相关分析显示,樟子松幼树叶片膜透性变化与MDA含量变化相关性几乎为零,可溶性蛋白与脯氨酸含量呈显著正相关关系,可溶性糖含量与脯氨酸含量呈显著负相关关系。研究表明,沙埋深度低于樟子松株高以上2cm能够促进其幼树生长;沙埋并没有导致樟子松幼树体内的膜脂过氧化,也没有引起细胞膜的损伤,在受到沙埋胁迫时,樟子松幼树体内SOD、POD以及可溶性蛋白和脯氨酸分别在防止其膜脂过氧化和维持细胞膨压中起到重要作用,而可溶性糖含量在沙埋过程中没有起到渗透调节作用。  相似文献   

17.
樟子松幼苗生长及光合特性对强风沙流吹袭的响应   总被引:1,自引:0,他引:1  
为了解樟子松幼苗对不同时间强风沙流吹袭的生理生态响应,2013年春季在内蒙古科尔沁沙地研究了8级大风风沙流(风速18m·s-1,风沙流强度173g·cm-1·min-1)吹袭10、20和30min下樟子松幼苗生长与光合特性的变化。结果显示:(1)随着风吹时间的增加,樟子松的株高生长量减少、茎粗生长加快,落叶数量增加,其中30min处理与CK相比的株高生长量下降52.63%,茎粗生长量增加233.30%,落叶指数增加466.70%。(2)风沙流吹袭没有改变樟子松幼苗的净光合速率、蒸腾速率、气孔导度和胞间CO2浓度的日变化规律,但日光合峰值下降,日最大蒸腾速率增加;与CK相比,30min处理的日最大光合速率下降22.69%,日最大蒸腾速率增加11.89%。(3)随风吹时间增加,其叶片温度、叶片相对含水量、日均光合速率、水分利用效率下降,30min处理较CK依次下降0.60%、4.37%、28.57%和31.58%,且日均蒸腾速率、气孔导度和胞间CO2浓度增加,30min处理较CK依次增加6.25%、6.67%和12.60%。研究表明,在持续风沙流胁迫下,樟子松幼苗光合作用受到抑制,蒸腾耗水增加,水分利用效率降低;樟子松幼苗生长速率降低主要源于光合面积减少和光合作用受到抑制,而其光合速率下降主要因幼苗叶片叶温和叶片含水量下降所致,蒸腾速率的增加主要源于气孔导度的增加;为了适应风沙流持续吹袭的胁迫,樟子松幼苗采取了降低株高生长速率,加快茎粗生长速率的适应策略。  相似文献   

18.
刘明慧  柳叶  任悦  高广磊  丁国栋  张英  赵珮杉  刘轩 《生态学报》2023,43(23):9912-9924
为揭示科尔沁沙地樟子松(Pinus sylvestris var.mongolica)人工林土壤真菌间相互作用,以不同林龄(中龄、近熟和成熟)、不同土层(0—10和10—20 cm)樟子松人工林土壤为研究对象,以未造林地为对照样地,采用分子生态网络分析法,基于随机矩阵理比较分析樟子松人工林土壤真菌网络特征及其影响因素。研究结果表明:(1)中龄林到近熟林阶段,科尔沁沙地樟子松人工林土壤真菌网络平均连通度增大,平均路径长减小,网络愈加复杂,而近熟林到成熟林阶段则相反,网络愈加简单;不同土层间,10—20 cm土层土壤真菌平均路径长度小于0—10 cm土层,平均连通度和聚类系数大于0—10 cm土层,10—20 cm较0—10 cm土壤真菌间相互联系更为密切;与樟子松人工林相比,未造林地土壤真菌网络中的平均连通度略低于近熟林,平均路径长低于樟子松人工林,未造林地土壤真菌群落更不稳定。(2)中龄、近熟和成熟樟子松人工林土壤真菌关键菌种分别为Guehomyces sp.、Oidiodendron sp.和Pseudeurotium sp.,未造林地关键菌种为Alphamyces sp.。在0—1...  相似文献   

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

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