《近10年全国省域农机总动力的区域差异时空特征》由会员分享，可在线阅读，更多相关《近10年全国省域农机总动力的区域差异时空特征（8页珍藏版）》请在金锄头文库上搜索。
1、 近10年全国省域农机总动力的区域差异时空特征 0 引言农业机械总动力指主要用于农、林、牧、渔业的各种动力机械的动力总和，是反映农业机械化发展状况的一个重要指标，也是政府部门制定农机化发展规划重要参考数据1.在发展现代农业、建设社会主义新农村的背景下，全国农业机械化作业水平呈现稳步提升的态势，农业机械总动力也不断增长。因此，从总体上把握和认识农业机械总动力的分布和变化规律，对科学制定农业机械动力发展规划十分必要。当前，在农业机械动力方面的相关研究主要包括农机总动力需求预测、农机总动力影响因素及农机动力对粮食产量的影响等2 -7.这些研究从国家或省级的尺度出发，揭示了农机总动力的时间变化趋势及其
2、影响因素的作用规律； 但忽略了农机总动力在区域内部表现出的空间自相关及异质性。探索性空间数据分析（ Ex-ploratory Spatial Data Analysis） 是空间数据分析技术的一种，常用来描述研究对象的空间分布规律并用可视化方式进行表达8.其内容包括探测聚合数据的空间特征和空间模式、识别空间数据的异常值、拟合对应于数据地理方面的假说，以及揭示现象之间的空间相互作用机制等9.ESDA 强调发现数据的空间分布模式，注重数据的空间依赖性和空间异质性，可以用来探索农机动力的空间关联格局问题。目前，将其用于农机动力这一主要农业生产投入要素的空间格局变化的研究还鲜见报道。而我国农业机械化空
3、间区域差异显着10,东部地区耕地占全国耕地总面积的32. 1% ,而农机动力拥有量占全国农机总动力的50. 2% ; 西部地区耕地占 23. 7% ,农机动力却只有15. 3% .因此，利用探索性空间数据分析技术可以更加深入地分析我国省域农机动力投入的空间格局特征及其演变规律。基于此，本文以中国 31 个省级行政区为基本研究单元，利用基尼系数与 ESDA 空间分析技术，对近 10 年全国省域农机总动力的区域差异时空特征以及空间演变格局进行分析，以期从新的视角为我国农业机械宏观布局与规划决策提供科学依据。1 材料和方法1. 1 数据来源研究所需的全国农机总动力及农作物播种面积数据来源于中国统计年
4、鉴 2004 - 2013。研究区不包含港澳台地区。中国行政边界数据来源于中国国家基础地理信息中心发布的 1 400 万比例尺基础地理矢量数据集。1. 2 研究方法1. 2. 1 基尼系数基尼系数是由意大利经济学家 Corrado Gini 提出的一种用于刻画收入分配不均的指标11.由于该指数能非常方便地反映出总体收入差距状况而成为国际上非常流行的指标。其值介于 0 1 之间： 当值为 0时，表明分配完全平均； 当值为 1 时，表明分配极其不平衡12.本文利用基尼系数来测度全国省域单位农作物面积农机动力投入的差异情况，并在开源的科学统计分析及绘图 R 语言环境中，利用 Achim Zeilei
5、s 等开发的 ineq 软件包计算得到全国省域农机动力基尼系数。1. 2. 2 ESDA 方法ESDA 是基于数据驱动的，注重发现空间数据的分布模式，揭示数据的空间自相关性与空间异质性的可视现象13.空间自相关是指位于不同空间位置的地理事物属性值之间具有统计相关性，即空间距离临近的地理单元属性值趋于相似14.空间自相关系数是度量地理单元空间相关性的指标，一般分为全局自相关系数和局域自相关系数两大类。全局型自相关系数反映了观测变量在整个研究区域内空间相关性的总体趋势，从而判断此现象是否存在聚集特征。在众多全局统计量中，Morans I 是最常用的评测指标，其取值范围为 - 1 1 之间，正值表示
6、该地理事物属性值呈现聚集分布（ 高 - 高聚集，低 -低聚集） ,负值表示该地理事物属性值为离散分布； 指数绝对值越大表示聚集或离散的程度越深，值为 0 则表示不存在空间相关，即呈现空间随机分布。全局空间自相关统计量反应了整体上地理单元之间空间相关的平均程度，但却忽略了可能存在的局部不平稳； 而局域指标可反映一个区域单元上的某种属性值与邻近区域单元上同一属性值的相关程度15.本文采用 Local Morans I 统计量并结合 LISA 聚集图，来研究全国省域农机动力的局部空间分布规律。2 结果与分析2. 1 农机动力区域差异时空特征2. 1. 1 农机动力区域差异空间特征统计全国各省级行政区
7、 2003 - 2012 年单位面积农业机械总动力数据，并利用 ARCGIS 软件绘制新疆近 10 年省级尺度的平均单位农业总动力空间分布情况。由图 1 可知： 全国各省域单位农业机械动力差异明显，上海市的单位农业动力最低为 2. 5kW/hm2; 而位于东部地区的天津市单位农业动力高达 13kW/hm2,区域差异明显。从空间分布来看： 地处华北平原的天津、北京、山东等地的农机动力投入较高，位于9. 5kW / hm2以上； 西南地区的四川、重庆、贵州、云南等农机动力投入最低，小于 2. 7kW/hm2; 其他地区农机动力投入则处于中间水平。由图 1 还可看到： 农业机械化水平较高的黑龙江及新
8、疆等地单位农机动力投入排名并不靠前； 而农业经济发展水平相对落后的西藏、青海等地单位农机动力投入却较高。这主要是由于农业机械化水平的高低与农机作业水平、农业机械化效益等因素核心相关，而单位农机动力仅是参考指标，并不是决定性因素。如黑龙江、新疆地区拥有大量的农垦及团场区域，耕地更加规整，更加便于大型农业机械作业，使其使用较少的单位农机动力投入就能达到较高的农业机械化水平。而西藏、青海等地却呈现相反的特点。从整体上看，全国各省域农机动力投入空间分布整体呈现东部大于西部、北方大于南方、由东南沿海向西部递减的趋势。2. 1. 2 农机动力区域差异时间演变为研究全国省域单位播种面积的农机动力投入区域差异
9、的时间演变特征，计算了其 2003 - 2012 年全国、东部、中部及西部的逐年基尼系数值，结果如图 2所示。由图 2 可知： 全国省域单位农机动力投入基尼系数基本位于 0. 3 以下，表明各省域差异总体处于相对均衡水平； 从时间变化趋势来看，2003 - 2012 年全国省域单位农机动力投入基尼系数除 2007 年有小幅升高外，总体呈现明显下降趋势，由2003 年的0. 31 下降到 2012 年的 0. 25,表明近 10 年来全国省域单位农机动力投入的区域差异不断缩小。但全国不同地区的农机动力投入差异及演变特征并不一致。从基尼系数大小上看： 东部和西部地区基尼系数较高、中部地区基尼系数值
10、较小，表明东部和西部各省域农机动力投入差异较大； 而中部地区各省域农机动力投入差异较小、投入较为均衡。从基尼系数时间变化来看：东部和中部地区省域单位农机动力投入的基尼系数总体下降趋势比较明显，区域差异缩小趋势也最为显着，这与东、中部地区快速发展的经济水平以及工业对农业发展的反哺、促进作用密不可分； 西部地区则表现出与东、中部完全不同的变化特征，2003 - 2012年西部各省域单位农机投入差异则呈现相反的上升趋势，表明其西部地区各省域的单位农机动力投入差异呈扩大趋势。因此，近 10 年全国省域单位农机动力投入区域差异总体呈明显缩小趋势，但不同地区之间又存在差异性。2. 2 全局空间自相关分析利
11、用 ARCGIS 空间统计工具计算 2003 - 2012 年逐年省域农机动力投入的全局 Morans I 指数，结果如表 1 所示。由表 1 可知，全国近 10 年省域单位农机动力 Mo-rans I 指数均为正值，且通过 99% 的置信水平检验（ Z得分 2. 58） ,说明各省域单位农机动力投入高的地区和低的地区出现相对集中，即全国省域单位农机动力存在显着的空间正相关性。此外，近 10 年来，Morans I指数整体上呈现缩小的趋势，表明随着时间的推移全国省域单位农机动力空间聚集在不断减弱。因此，全国省域单位农机动力空间分布并非完全的随机分布，而是表现出明显的空间聚集，并且随着时间的推移
12、其空间聚集在不断减弱，各省农机动力投入趋向均衡发展。2. 3 局部空间自相关分析全局 Morans I 指数作为测度总体自相关特性的统计量，往往掩盖了区域内部具体地区的空间集聚特征，而局域 Morans I 指数可进一步揭示区域内部的局域自相关性16.因此，本文利用 Moran 散点图和显着性水平 LISA 图来对 2003 年和 2012 年两个横截面数据进行分析，以探索全国省域农机动力投入的局部空间集聚特性。2. 3. 1 Moran 散点图分析Moran 散点图不仅可以探索地理单元间的全局空间关联特征，还可以用来识别局部空间关联模式。Moran 散点图将全国省域农机动力投入数据点划分为4
13、 个象限，分别表征 4 种不同的局部空间关联形式17: 第 1 象限是高高聚集，该象限内的省域及其周围省域的农机动力均较高，形成高高聚集； 第 2 象限是高低聚集，该象限内的省域单位农机动力投入较高，但其周围的省域单位农机动力投入较低； 第 3 象限是低低聚集，该象限内的省域及其周围省域的单位农机动力投入均较低； 第 4 象限是低高聚集，该象限内的省域单位农机动力投入较低，但其周围的省域单位农机动力投入却较高。其中，落入第 1、第 3 象限的省域单元属性值与相邻单元属性值趋同，即存在较高的空间正相关，具有同质性； 而落入第 2、第 4 象限的空间单元属性值与相邻单元差异显着，存在空间负相关，具
14、有较高的空间异质性。利用 Geoda 软件绘制全国农机动力 2003 和 2012年的 Moran 散点图，并对各象限省域单元数量其进行统计，如表 2 所示。由表 2 可知，2003 和 2012 年全国省域单元分布在第 1、第 3 象限的比例分别为 32%、19%,第 2、第 4 象限比例分别为 68%、81%.第 2、第 4 象限省域个数明显多于第 1、第 3 象限，说明省域单位农机动力投入更多呈现空间异质性，即高低聚集和低高聚集的趋势明显。从时间变化上看，相对研究初期的 2003 年来说，2012 年位于第 1、第 3 象限的省域数量均有所下降，表明空间自相关聚集趋势趋于减弱，而呈现较多
15、的空间异质性。2. 3. 2 显着性水平 LISA 聚类图分析由于利用 Moran 散点图表示的区域分异并没有经过统计假设检验，因此只是用来作为参考。为了更直观地反映全国省域单位农机动力投入的局域空间格局演变过程，本文利用 Geoda 软件并结合 ARCGIS 绘制出在 0. 05 显着性水平下的 LISA 聚类图，如图 3 所示。由图3 可知： 在两个典型年份，全国省域农机动力空间聚集类型中的高高、低低聚类均有省域通过 0. 05的显着性检验； 而高低和低高聚集两种类型未有省域通过显着性检验。由图 3（ a） 可知： 2003 年全国省域农机动力的高高聚集区主要分布在华北平原的北京、天津及河北地区； 而低低聚集区主要分布在西南部的重庆、四川及湖南地区。随着时间的推移，2012 年的高高聚集省域单元有较大减少，只有北京市，表明全国省域农机动力高高聚集有明显的下降趋势，省域间实现了较为均衡的发展。同样，2012 年的低低聚集区空间分布位置变化并不大，还是位于西南部，省域域单元数量也有明显的下降，仅有湖南省。总体来看，由 LISA 聚类图分析得到的各聚类类型省域单元数量变化的结果，与表 2 基本一致。因此，近 10 年全国省域农机动力投入存在显着的空间自相关特征（ 高高聚类、低低聚类） .其中，高高聚集区主要分布于华北平原
