粮食单产空间格局改变.docx

粮食单产空间格局改变0引言随着中国人口的增加、社会经济的发展、耕地面积的减少粮食安全已成为一个紧迫的问题［1］。我国的粮食安全问题实质上是粮食产量问题［2－3］。而粮食产量是众多因素共同作用的结果［4］。近些年来一些学者通过对我国粮食产量进行的研究表明几十年来驱动我国粮食产量变化的直接推动力是单产其次才是播种面积［5］。我国自1952年来粮食总产的增加已经可以完全归结为单产水平的提高［6］。粮食单产水平及其波动状况直接影响粮食产量及其稳定性单产变化及其波动的区域差异已经引起我国许多学者的关注。党安荣等［5］采用GIS分析手段通过对1984年以来我国粮食单产空间变化的定量研究表明在1984－1994年这10年中我国北方地区粮食增产强度较大、南方地区下滑现象明显。程叶青［6］采用数理统计和GIS空间分析方法通过对1980－2005年东北地区183个县(市)的粮食单产变异的分析揭示了东北地区粮食单产水平空间格局变化特征为东北地区粮食可持续生产提供了科学依据。钟甫宁等［7］通过变异率和变异系数两个参数分析了全国各省市单产波动特征结果表明粮食主产区的波动系数和波动幅度比其他地区高。殷培红等［8］利用去趋势波动方法对全国2075个县(市)的20年粮食单产数据进行分析后指出吉林、辽宁东部、河南、安徽地区粮食单产波动系数较高。以往研究尽管基于国家或省级尺度较好地揭示了粮食单产宏观格局的变化规律却囿于统计分析和GIS的空间表达难以反映区域内部粮食生产的异质性。而基于GIS平台的探索性空间数据分析ESDA(ExploratorySpatialDataAnalysis)方法是通过对事物或现象空间分布格局的描述与可视化发现空间集聚和空间异常解释研究对象之间的空间相互作用机制的一种新的有效的测度方法。目前此方法已广泛应用于区域人口［9］、农业［10］、土地［11－12］等各种社会经济要素［13］的空间格局分布研究上但用于粮食单产这一主要农业经济要素的地域格局变化的研究还鲜见报道。河北省是全国重要的商品粮食生产基地之一［14］在全国粮食生产中具有重要的地位。耕地面积减少和人口增长的双重压力使河北省粮食生产只能走主攻单产之路因此如何充分挖掘耕地增产潜力稳定提高粮食单产就显得相当重要。但河北省地貌类型多样农业类型区跨度较大粮食单产水平空间分布不尽相同。基于此本文以河北省136个县(市)为研究单元基于泰尔指数模型与ESDA空间分析技术对1987－2009年河北省136个县粮食单产的区域差异特征以及动态演变格局进行了动态研究以期从一个崭新的途径上反映区域粮食单产的差异特征为该区域粮食生产宏观布局与规划决策提供科学依据。1研究区域及研究方法1．1研究区概况河北省地跨36度03'N～42度40'N113度27'E～119度50'E之间总面积为18．88万km2。该区域平均气温自南向北自东向西逐渐降低;地形地貌类型多样平原分布在东南方山地呈半环状耸于西北部高原镶嵌在西北边缘由海陆节节上升。河北省拥有良好的农业基础和农业区位优势是《全国新增1000亿斤粮食生产能力规划(2009－2020年)》中确定的重点区域之一。根据全国农业生态类型区的划分标准河北省范围内包含2个国家一级类型区和4个国家二级类型区如表1和图1所示。1．2研究方法1．2．1泰尔系数研究区域差异的定量化指标很多［15－17］如极差、标准差、变异系数和泰尔(Theil)系数等其中泰尔指数和变异系数应用较为普遍。泰尔指数(Theil)的意义在于可将区域总体差异分解成不同空间尺度的区内差异(TWR)和区间差异(TBR)便于比较它们对区域整体差异的影响和贡献。按照国家二级区的划分对泰尔指数进行一阶段分解可以将河北省总体差异分解4个生态类型区区间的差异和4个生态类型区内各县市之间的差异其计算公式为［18］式(1)中N为县市数yi为i地区的粮食产量占河北省的份额Pi为i县市的播种面积占河北省的份额。式(2)中T为泰尔指数TwR为区内差异TBR为区间差异yi为i区播种面积占河北省的总播种面积的比重pi为i区粮食产量占河北省的比重yij为i区第j个县市占i区的播种面积pij为i区第j个县市的粮食产量占i区的粮食产量密度。则其中Cv为变异系数y为河北省的平均粮食单产yi为县市i的粮食单产n为县市个数。1．2．2ESDA方法ESDA方法是一系列空间分析方法和技术的集合［19］为事物或现象的空间相关性测度提供了有效工具。根据研究目的的不同可分为全局空间自相关和局部空间自相关。1)全局自相关。全局空间自相关反映了观测变量在整个研究区域内空间相关性的总体趋势Moran‘sI指数是空间自相关分析中普遍采用的一种方法。由于粮食单产是一个比率在分析比率变量的