成渝经济区植被生态质量演变特征及驱动力分析.pdf

第 卷第期 年月物探化探计算技术 收稿日期：基金项目：中国博士后科学基金第 批特别项目（）第一作者：王丽芝（），女，硕士，主要研究方向为“”技术与数字国土，：。通信作者：杨鑫（），男，博士，副教授，主要研究方向为“”技术在地学中的应用，：。文章编号：（）成渝经济区植被生态质量演变特征及驱动力分析王丽芝，杨鑫，郝利娜（成都理工大学地球科学学院，成都 ；成都理工大学地球勘探与信息技术教育部重点实验室，成都 ）摘要：利用 数据反演成渝经济区 年 年植被覆盖度（犞 犉 犆）、植被净初级生产力（犖犘 犘），构建植被生态质量指数（犙）。通过趋势线性分析、相关性分析、残差分析，研究 年 年成渝经济区植被生态质量时空变化特征，探索成渝经济区植被生态质量与气候和人类活动等影响因素的关系。结果表明：年 年成渝经济区植被覆盖度达，犖犘 犘显著升高地区达，植被生态质量指数（犙）显著增加地区达，三者之间空间分布相对一致，植被生态质量整体提升；通过趋势分析，降水、温度增加的空间分布和植被生态质量指数（犙）增加的空间分布范围相对一致，表明气候对改善植被生态质量有重要的驱动作用；相关性分析显示，植被生态质量指数（犙）与降水呈正相关，与温度呈负相关，且降水与植被生态质量指数（犙）的相关性强于温度；驱动分区显示，成渝经济区植被生态质量变化约 地区受气候因素驱动，约 地区受非气候因素（自然灾害和人类活动等）的驱动；通过残差分析，年 年成渝经济区人类活动对植被的正面作用大于负面作用，说明成渝经济区植被生态质量的时空变化是气候因素和人类活动共同影响的结果。关键词：成渝经济区；植被生态质量；时空变化；驱动力中图分类号：文献标志码：犇 犗 犐：引言植被作为生态系统不可或缺的组成成分，不仅有联结土壤、水文及环境的作用，并且植被在区域初级净生产力、气候系统以及陆地碳平衡中发挥着重要作 用，对 生 态 自 然 环 境 有 较 为 直 观 的 指 示 作用。同时，植被作为生态系统中的重要组成部分，其具有丰富的生态过程、复杂的结构和生物多样性，植被的变化影响到许多重要的生态系统服务功能，在涵养水源、防风固沙、固碳释氧、维护生物多样性、调节气候和保持生态平衡等方面具有不可替代的作用。因此，植被的生态环境质量（植被生态质量）是评价生态系统平衡状况和环境质量的一项重要指标。植被生态质量的好坏决定着一个地区的环境质量及生态系统平衡状态。随着遥感技术的快速发展和广泛使用，利用遥感对植被特征进行研究的技术与方法也越来越成熟。目前，常用遥感手段，采用单一指标法和多指标法监测植被生态质量。徐涵秋等采用单一指标法，利用不透水面盖度评价城市生态环境等；伊博力等采用单一指标法，利用植被指数间接地对草地生态环境系统进行评价监测。单一指标只能片面的反应植被生态质量某一维度的变化特征，而植被生态系统是一个错综复杂的生态系统，需要综合各种指标才能较为完整的解释其变化特征。吴 宜 进 等，采 用 多 指 标 法，引 入 改 进 的 ，以 数据反演得到植被总初级生产力、叶面积指数、区域热度、改进的湿度与植被覆盖度项指标为基础，提出基于主成分分析法的西藏植被生态质 量 评 价 方 法；李 超 等采 用 多 指 标 法，以 数据反演得到中国总初级生产力、叶面积指数、植被覆盖度、温度植被干旱指数和陆面温度项指标，综合评价了中国植被生态质量。但是利用遥感数据构建多指标植被生态质量评价模型，会受到地区地势、时间及尺度上的限制，难以实现精确的评价，这就需要探索一种更适用的植被生态质量指数构建模型。近年来，诸多学者对西南地区的植被动态变化研究诸多，单独针对成渝经济区植被动态变化的研究相对较少，多采用单一指标法对成渝经济区植被动态变化进行了分析。但成渝经济区自然环境复杂，单一指标无法精确的反应成渝经济区植被生态质量的变化特征。笔者利用 数据反演得到的成渝经济区 年 年植被覆盖度（犞 犉 犆）、植被净初级生产力（犖犘 犘）指标为基础，构建植被生态质量指数（犙），通过趋势线性分析、相关性分析、残差分析，分析了此 间成渝经济区植被生态质量时空变化特征，探究了气候变化和人为活动等驱动因素对成渝经济区植被生态质量时空变化的制约关系，为成渝经济区植被生态保护与防护提供理论基础和科学依据。研究区及数据来源 研究区概况成渝经济区由区域内的两大核心城市而得名，范围包括四川省 个市，以及重庆市 个区县（图）。东部有三峡库区这一国家重点工程，西部有汶川地震灾后重建区，是长江中上游生态屏障的重要组成部分，。成渝经济区位于西南腹地、四川盆地核心地区，地貌以丘陵为主，中心为平原，四周山体环绕，由于其特殊的盆地地形，其内部热量不易扩散，从而使得成渝经济区的温度高于同纬度的其他地区，并且东部和南部气温整体高于西部和北部，降水量则表现为东多西少，盆地多高原少。成渝经济区属亚热带湿润季风气候，区内植被以亚热带常绿阔叶林和高原山地高寒草甸为主。图 成渝经济区地理位置及概况图 数据来源及预处理研究区 和行政边界数据源自中科院资源环境科学数据中心（：），气象数据来源于中国气象数据网（：），选用成渝经济区 个气象站 年月 年 月的中国地面气候资料日值数据集。对 数据进行投影转换、镶嵌等处理，利用波段计算等方法反演得到植被覆盖度（犞 犉 犆）和植被净初级生产力（犖犘 犘）。统计每个站点年降水量和年均温度，把气象站点坐标生成站点值数据，利用薄盘光滑样条函数分别生成研究区降雨和气温的空间分布图。将所有影像数据统一重采样为 分辨率数据，以便后续进行分析。研究方法 植被覆盖度（犞 犉 犆）植 被 覆 盖 度（，犞 犉 犆）指植被（包括叶、茎、枝）在地面的垂直投影面积占研究区总面积的百分比，它量化了植被的繁茂程度。目前已有很多利用遥感监测手段估算植被覆盖度的方法，其中最为常见的是利用植被指数犖犇犞 犐进行估算，其计算公式为式（）。犞 犆 犉犖 犇 犞 犐犖 犇犞 犐 （犖 犇犞 犐 犖 犇犞 犐 ）（）式中：犖犇犞 犐 为完全是裸土或无植被覆盖区域的犖犇犞 犐值；犖犇犞 犐 为完全被植被所覆盖的像元的期王丽芝，等：成渝经济区植被生态质量演变特征及驱动力分析犖犇 犞 犐值，即纯植被像元的犖犇 犞 犐值。植被生态质量指数（犙）模型利用植被净初级生产力犖犘 犘和植被覆盖度构成植被生态质量指数，综合监测评估植被生态质量的优劣，其值越大，表明植被生态质量越好。植被生态质量指数（犙）模型为式（）。犙犻（犞 犉 犆犻犖 犘 犘犻犖 犘 犘犿 ）（）其中：犙犻为第犻年植被生态质量指数；犞 犉 犆犻为第犻年植被覆盖度；犖犘 犘犻为第犻年植被犖犘 犘；犖犘 犘犿为监测年份中犖犘 犘最大值，即当地最好气象条件下植被犖犘 犘。趋势线分析法线性回归方程是最广泛应用于分析不同因子和研究对象之间的变化趋势方式。计算公式为式（）。犛 犾 狅 狆 犲狀狀犻犻 狓狀犻犻狀犻狓狀狀犻犻（狀犻犻）（）式中：犻为年份；狓为第犻年的研究对象值；狀为时间长度；犛 犾 狅 狆 犲为狀年时间内，影响因子的变化对应的研究对象总体变化趋势。相关性分析法相关性分析是指对两个或多个具备相关性的变量元素进行分析，从而衡量两个变量因素的相关密切程度。线性相关性分析的前提条件是相关性的元素之间需要存在一定的联系或者概率。线性相关性分析犚的计算公式为式（）。犚狓 狔狀犻（狓犻狓）（狔犻狔）狀犻（狓犻狓）槡狀犻（狔犻狔）槡（）其中：犚狓 狔为狓、狔之间的线性相关系数；狓犻、狔犻为两个影响因素在第犻年的值；狓、狔为两个影响因素狀年的平均值。偏相关分析是在线性相关的基础上，排除其他影响因素的干扰，计算某两个影响因素的相关性。偏相关系数计算公式为式（）。犚狓狔，狕犚狓狔犚狓狕犚狔狕犚槡狓狕犚狔槡狕（）式中：犚狓 狔，狕为将自变量狕固定后，因变量狓与自变量狔的偏相关系数。用狋检验法进行显著性检验，公式为式（）。狋犚狓狔，狕犚狓狔，槡狕狀犿槡（）式中：狀为样本数；犿为自变量个数。为了综合考虑多个相关影响因素的共同作用，将两个或多个影响因素采用复相关分析法并入计算，计算公式为式（）。犚狓，狔狕（犚狓狔）（犚狓狕，狔槡）（）式中：犚狓狔表示狓与狔的线性相关系数；犚狓狕，狔为固定自变量狔；因变量狓与自变量狕之间的偏相关系数。用犉检验法进行显著性检验，公式为式（）。犉犚狓，狔狕犚狓，狔狕狀犽犽（）式中：狀为样本数；犽为自变量个数。残差分析法残差趋势法（）是一种逐像元建立植被响应气候的线性回归模型，分别以预测值和残差值（实际值与预测值之差）随年份的变化趋势来表达自然和人为因素对植被生态质量的独立作用。）利用研究区植被生态质量指数犙和气象数据建立回归方程计算模拟犙值。犙犪狓犜狓犘（）式中：犪表示回归系数；狓、狓为常数；犜为温度；犘为降水。）求解回归系数和常数为式（）式（）。狀犻（犜犜）犻狓狀犻（犜犜）犻（犘犘）犻狓狀犻（犜犜）犻（犙犙）犻（）狀犻（犘犘）犻狓狀犻（犜犜）犻（犘犘）犻狓狀犻（犘犘）犻（犙犙）犻（）犪犙狓犜狓犘（）将计算出的回归系数和常数带入回归方程得到模拟植被生态质量指数犙值，用实际犙值减去模拟犙值得到犙残差值。实验与分析 成渝经济区 时空变化利用 趋势分析方法对犞 犉 犆变化的趋势进行分析计算，对成渝经济区植被覆盖度变化趋势进行分析，结合研究区实际情况对犞 犉 犆进行分级，如表所示。物探化探计算技术 卷表植被覆盖度（犞 犉 犆）退化分级表 （犞 犉 犆）退化等级 值 表植被净初级生产力（犖犘 犘）退化分级表 犖犘 犘 退化等级 值 图 年 年犞 犉 犆年均变化曲线及时序变化趋势图 犞 犉 犆 年 年成渝经济区植被覆盖度整体呈波动式上升（图），年年均犞 犉 犆为 ，年 年 均犞 犉 犆为 ，平 均 增 加 速 率 为 。年 年成渝经济区植被覆盖度达 以上，植被覆盖度整体较高。成渝经济区 年 年植被覆盖度从西向东呈增加趋势，其中犞 犉 犆变化趋势呈增加的区域占全域的 ，犞 犉 犆呈 减 少 趋 势 的 区 域 占 全 域 的 ；犞 犉 犆退化的区域占全域的 ，主要集中在成渝经济区的城镇集中区；犞 犉 犆基本不变的区域占全域的 ，主要分布在以遂宁市、资阳市为中心的中心圈地区；犞 犉 犆改善的区域占全域面积的 ，主要分布在遂宁市、资阳市以及成渝经济区边缘地带等。整体来看，成渝经济区植被覆盖度（犞 犉 犆）从西向东呈改善退化改善趋势，植被覆盖整体呈现改善状态。成渝经济区 时空变化植被净初级生产力（犖犘 犘）是指绿色植物在单位时间、单位面积内所积累的有机物数量，是由光合作用所产生的有机质总量（，）中扣除自养呼吸（，）消耗量后的剩余部分，用以表征绿色植物的固碳 能 力 。利 用 趋 势 分 析 方 法 对犖犘 犘变化的趋势进行分析计算，结合研究区实际情况对犖犘 犘进行分级，如表所示。从图（）来看，年 年成渝经济区犖犘 犘总体呈波动式增长，年均 值从 年的 增长到 年的 ，平均增长速率为 。从图（）看，犖犘 犘从西向东呈现退化改善基本不变的趋势，其中 的地区犖犘 犘减少，的地区犖犘 犘增加，犖犘 犘增加速率介于 的区域占研究区面积的 ，增加速率大于 的区域占研究区面积的 。整体来看，成渝经济区中心区域（如资阳、内江、潼南、大足、铜梁等）犖犘 犘数值呈逐年上升趋势，植被净初级生产力呈改善状态；成渝经济区边缘地区（如雅安市等）犖犘 犘数值呈逐年下降趋势，植被净初级生产力呈退化趋势。表明成渝经济区在城市化集中地区大力发展的同时，也在增强植被生态保护。植被生态质量指数（犙）时空变化利用 趋势分析法，对成渝经济区植被生态质量指数（犙）进行年际趋势变化分析，结合研究区实际情况对植被生态质量指数（犙）进行分级（表），摸清了成渝经济区植被生态质量的退化改善情况。期王丽芝，等：成渝经济区植被生态质量演变特征及驱动力分析图 年 年犖犘 犘年均变化曲线及时序变化趋势图 犖犘 犘 图 年 年植被生态质量指数（犙）年均变化曲线及时序变化趋势图 （犙）植被生态质量是衡量一个地区自然生态状态和环境质量的关键指标。从图（）来看，成渝经济区植被生态质量指数（犙）整体呈明显的增加趋势，但存在一定的波动，表明成渝经济区植被生态质量整体处于变好的趋势。年 年成渝经济区植被生态质量指数（犙）介于 之间，平均值为 ，平均年际增长速率为 。年植被生态质量指数（犙）最小，与气候灾害和人为因素（城市基础建设等）有关，年植被生态质量指数（犙）最大，与人类活动（退耕还林还草、植树造林等）有关。从图（）看，成渝经济区植被生态质量指数（犙）从西向东呈退化改善基本不变的空物探化探计算技术 卷间格局。植被生态质量退化的地区主要分布在成渝经济区西部（绵阳、成都、雅安、乐山）；植被生态质量明显改善的地区主要分布在成渝经济区中心地区（眉山市、资阳市、内江市和重庆市以西）；植被生态质量基本不变的地区主要分布在成渝经济区东部地区。整体上，成渝经济区植被生态质量呈退化状态的地区占全区的 ，呈改善和不变状态的地区占全区的 。年际变化趋势表明，成渝经济区 年 年植被生态质量总体呈变好趋势。表植被生态质量指数（犙）退化分级表 （犙）退化等级 值 图 年 年降水与温度时序数据变化趋势分析图 图 年 年植被生态质量指数（犙）与降水、温度的偏相关系数空间分布 （犙）期王丽芝，等：成渝经济区植被生态质量演变特征及驱动力分析 植被生态质量与气候因素的相关性分析气候是植被生长的自然条件，植被是气候变化的反馈者。从图可知，成渝经济区年降水量变化趋势从北向南呈减少增加减少趋势，年均温度变化趋势呈明显增加减少趋势。成渝经济区内 的地区年降水量呈明显增长趋势，尤其是成渝经济区城镇集中地区年降水量增加尤为明显；成渝经济区内 的地区年降水量呈明显下降趋势，雅安市、绵阳市和达州市年降水量下降的最为明显。成渝经济区温度上升地区占 ，温度下降的地区占 ，其中成都市及其周边地区温度逐年上升，重庆市温度逐年下降。年 年成渝经济区年降水与植被生态质量指数（犙）的偏相关系数介于 之间，平均值为 （图（），表明成渝经济区植被生态质量与降水呈显著正相关关系。其中呈负相关的地区主要集中在雅安市、绵阳市等，约占全区面积的 。呈正相关的地区主要集中在成都市、资阳市以及重庆市等城镇集聚地，约占全区面积的 。狋双侧检验表明，通过 显著性检验水平的地区主要分布在成渝经济区中心地区（如成都市、资阳市及其重庆市城市群等），约占正相关区域面积的 。温度与植被生态质量指数（犙）的偏相关系数介于 之间，平均值为 （图（），表明成渝经济区植被生态质量与温度呈不显著负相关关系。其 中 呈 负 相 关 的 区 域 约 占 全 区 面 积 的 ，主要集中在成都市、德阳市、资阳市以及重庆市等城市群；其余地区均呈正相关，约占全区面积的 。狋双侧检验表明，通过 显著性检验水平的地区主要集中在成渝经济区东北部和西北部地区，约占正相关区域面积的 。综上可知，成渝经济区植被生态质量指数（犙）与降水的相关性显著强于温度，说明降水对植被生态质量的影响程度大于温度。植被生态质量指数（犙）与气候因素的驱动分析植被生态质量指数（犙）与降水、温度的复相关系数介于 之间，平均值为 ，大部分地区的植被生态质量指数（犙）与降水、温度的复相关性较为明显（图）。资阳市、遂宁市和重庆市西北部等地复相关性较强，绵阳市、泸州市及雅安市等地复相关性最弱。总体上，成渝经济区从北向南复相关性呈弱强弱的趋势。检验结果表明，植被生态质量通过 显著性水平检验的地区约占全区图 年 年植被生态质量指数（犙）与降水、温度的复相关 （犙）面积的 。采用驱动分区方法 （表）对成渝经济区植被生 态 质 量 的 变 化 进行 驱 动 力分 析。年 年成渝经济区植被生态质量驱动分区（图）结果显示，气候对成渝经济区植被生态质量的影响较为明显。年 年成渝经济区植被生态质量变化受降水、温度共同影响的面积约占全区面积的 ；温 度 主 驱 动 区 域 约 占 全 区 面 积 的 ；降水主驱动区域约占全区面积的 ；其余非气候因素（自然灾害和人类活动等）影响区域约占全区面积的 。综上所述，成渝经济区植被生态质量变化约 的地区受气候因素驱动，约 的地区受非气候因素（自然灾害和人类活动等）驱动。整体上，年 年成渝经济区植被生态质量驱动机制表现为非气候因素驱动。人类活动对植被生态质量的影响残差分析是分离自然和人类活动对植被覆盖变化的影响。从图成渝经济区植被生态质量指数犙的残差值线性变化趋势可知，在研究时段内，残差值呈波动上升趋势，并且残差值由负逐渐转为正，说明人类活动对成渝经济区的植被生态质量影响由负面逐渐转为正面。年残差值最低，说明这一年人类活动对植被生态质量的负面作用最为明显；从 年开始残差值为正并持续上升，到 年达到最大，说明在这段时间内人类活动对植被生态质物探化探计算技术 卷图 年 年植被生态质量指数（犙）与降水、温度的驱动分区 （犙）表植被生态质量变化影响因子分区规则 项目驱动因子驱动准则狋检验（降水）狋检验（温度）犉检验 犜犘狋狋 狋狋 犉犉 犜狋狋 犉犉 犘狋狋 犉犉 犜犘狋狋 狋狋 犉犉 犖 犆犉犉 犜犘温度、降水强驱动；犜温度主驱动；犘降水主驱动；犜犘温度、降水弱驱动；犖犆非气候因子驱动量的正面作用在不断加强，对成渝经济区的生态保护有了全新的认识；年 年残差值虽有下降，但残差值仍为正，对植被生态质量仍起正面作用。整体看，人类活动在时间序列上对研究区内植被生态质量的正面作用大于负面作用，植被生态质量呈改善状态。期王丽芝，等：成渝经济区植被生态质量演变特征及驱动力分析图成渝经济区植被生态质量指数犙的残差值线性变化趋势 犙 结论笔者利用 年 年 数据反演得到的犞 犉 犆、犖犘 犘构建植被生态质量指数（犙），结合成渝经济区 个气象站的气温、降水数据，分析了成渝经济区植被生态质量时空变化特征，探讨了植被生态质量与气温降水和人类活动等影响因素的制约关系，主要结论如下：）成渝经济区 年 年植被覆盖度达 以上，植被覆盖度较高。犞 犉 犆增加的区域占全区的 ，犞 犉 犆减少的区域占全区的 。的地区犖犘 犘减少，的地区犖犘 犘增加。植被 生 态 质 量 呈 退 化 状 态 的 区 域 占 全 区 的 ，呈 改 善 和 不 变 状 态 的 区 域 占 全 区 的 。犞 犉 犆、犖犘 犘和植被生态质量指数（犙）的年度趋势变化空间分布相对一致。整体来看，成渝经济区 年 年植被生态质量良好并在逐年改善。）趋势分析显示，降水和温度升高的空间分布区域和植被生态质量指数（犙）升高的空间分布区域相对一致。定量分析了植被生态质量指数（犙）与气温和降水的相关性，结果表明植被生态质量与降水呈正相关、与温度呈负相关，说明成渝经济区植被生态质量受降水和气温的共同制约，但降水的相关性强于温度。）驱动分区表明，植被生态质量受气候因素驱动的区域约占 ，其中降水主驱动占比最大，约占 ；约 的地区受非气候因素（自然灾害和人类活动等）的驱动，说明 年 年成渝经济区植被生态质量驱动机制表现为非气候因素驱动。）由于驱动分区的局限性，不能准确地说明人类活动对植被生态质量的正负面影响，所以通过残差分析，研究发现 年 年，人类活动对成渝经济区的植被生态质量起负面作用；年 年，人类活动对研究区植被生态质量起正面作用。整体来看，人类活动对植被生态质量的正面作用强于负面作用，说明在成渝经济区大力发展、扩张城市的同时，对植被生态的保护意识逐渐加强。成渝经济区植被生态质量逐年改善，其时空变化是气候因素和人类活动共同影响的结果。综上，成渝经济区的植被生态质量正在逐年改善，人们对植被生态的保护意识在逐渐加强。但成渝经济区植被生态质量仍然存在空间分布不均的现象，有些地方植被生态质量相对较低，各市县政府还需有针对性地出台一些政策加强植被生态质量相对偏差区域的植被生态保护和治理力度。气候条件是植被生态质量好坏的前置条件，成渝经济区的气候条件相对较好，比较有利于植被的生长，但旱灾洪灾仍有发生，对植被生态质量有着不利的影响，应适时做好防洪防旱工作，保护植被生长环境。参考文献：赵杰，杜自强，武志涛，等中国温带昼夜增温的季节性变化及其对植被动态的影响地理学报，（）：，（）：（），（）：梁浩近 年来克拉地峡东西海岸带红树林变化遥感监测与对比分析呼和浩特：内蒙古大学，：，（）徐涵秋基于城市地表参数变化的城市热岛效应分析生态学报，（）：物探化探计算技术 卷 ，（）：（）伊博力基于神经网络的牧户尺度草地产草量监测研究呼和浩特：内蒙古师范大学，：，（）吴宜进，赵行双，奚悦，等基于 的 年西藏生态质量综合评价及其时空变化地理学报，（）：，（）：（）李超，李雪梅 年中国植被生态质量时空变化特征长江流域资源与环境，（）：，（）：（）陈付国，孙荣，梁少民成渝经济区陆地植被碳储量动态及碳 汇 战 略 研 究 环 境 科 学 与 管 理，（）：，（）：（）李茜荣，杨东，冯磊，等成渝经济圈 年植被 的动态变化生态学杂志，（）：，（）：（）金凡琦，程志刚，靳立亚，等成渝城市群热环境效应与植被覆盖度关系研究成都信息工程大学学报，（）：，（）：（）汲玉河，周广胜，王树东，等 年秦岭地区植被生态质量演变特征及驱动力分析植物生态学报，（）：，（）：（）黄林峰，田鹏举，帅士章，等 年赤水河流域植被生态质量变化分析中低纬山地气象，（）：，（）：（），：，（）：穆少杰，李建龙，周伟，等 一 年内蒙古植被净初级生产力的时空格局及其与气候的关系生态学报，（）：，（）：（），（）：，：杨婷婷，韦丹，马祥庆，等基于随机森林模型长汀县植被覆盖度的影响因子福建农林大学学报（自然科学版），（）：，（），（）：（）方精云，唐艳鸿，林俊达全球生态学：气候变化与生态响应北京：高等教育出版社，：期王丽芝，等：成渝经济区植被生态质量演变特征及驱动力分析 ：，（）叶红，张廷斌，易桂花，等 年黄河源区 时空特征及其与气候因子关系地理学报，（）：，（）：（）艾则孜提约麦尔麦麦提，玉素甫江如素力，姜红，等 年博斯腾湖流域 时空变化特征及影响因子分析草业科学，（）：，（）：（）犈 狏 狅 犾 狌 狋 犻 狅 狀犮 犺 犪 狉 犪 犮 狋 犲 狉 犻 狊 狋 犻 犮 狊犪 狀 犱犱 狉 犻 狏 犻 狀 犵犳 狅 狉 犮 犲 狊狅 犳 狏 犲 犵 犲 狋 犪 狋 犻 狅 狀犲 犮 狅 犾 狅 犵 犻 犮 犪 犾狇 狌 犪 犾 犻 狋 狔 犻 狀狋 犺 犲犆 犺 犲 狀 犵 犱 狌 犆 犺 狅 狀 犵 狇 犻 狀 犵犲 犮 狅 狀 狅 犿 犻 犮狕 狅 狀 犲 ，（，；，）犃 犫 狊 狋 狉 犪 犮 狋：（）（），（），：（），（），（），（），；（），（）（），（，）；，犓 犲 狔 狑 狅 狉 犱 狊：；物探化探计算技术 卷