0 引言
【研究意义】耕地是人类赖以生存和发展的基本资源和条件[1],其数量和空间分布直接关乎粮食安全乃至全球生态格局[2]。大洋洲虽然其耕地面积总量远小于其他大洲,但仍然是全球主要农产品的重要生产和供给基地。尤其是占大洋洲陆地面积85.7%的澳大利亚一直十分重视推动农业现代化和耕地集约化利用,采用先进的水权管理制度[3,4],成为目前世界重要的粮仓和牧场。已有部分****对全球[5]和APEC地区[6]开展了粮食安全评价,初步估算分析了大洋洲土地生产潜力和粮食生产潜力。随着农业技术的发展,抗旱节水粮食作物品种的开发,澳大利亚成为粮食生产和出口大国,对中国粮食出口不断加大。此外,大洋洲频发的极端天气气候事件对当地农业土地利用影响巨大,2001—2008年澳大利亚经历了长达8年的大旱,导致产量锐减,影响全球粮食供给。因此,科学掌握大洋洲的耕地利用格局及其时空变化规律,有利于准确分析大洋洲的粮食生产格局和粮食供给状况,更好制定我国农产品国际贸易政策,服务国家粮食安全战略。【前人研究进展】一般来说,在全球或大洲尺度开展耕地利用格局变化分析多采用遥感手段和统计资料[7,8]。如YAO利用CG-LTDR（全球及中国区域长时间序列卫星数据集）对全球1982—2011年耕地变化格局开展研究[9],但该地表覆盖产品为5 km,较粗的空间分辨率难以很好揭示大洋洲内部更为细致的耕地利用变化特征,研究结果不确定性也较大。虽然有的研究为得到更为精细的土地利用变化数据,将遥感影像与社会经济数据结合或直接采用参与式调查方法[10],但这种方法多局限于小区域研究,难以进行大尺度分析。【本研究切入点】最新研制的GlobeLand30是全球首套30 m分辨率的全球地表覆盖数据,包括2000和2010年两个基准年,为大区域耕地利用格局时空变化研究提供了新的数据资料[11,12,13]。GlobeLand30耕地制图产品精度高达83.06%[14],相比FROM-GLC、GlobCover、MODIS、Collection 5和MODIS Cropland等低分辨率的全球地表覆盖数据,在耕地数量计算和耕地空间位置描述方面更具优势[15]。可见,采用GlobeLand30数据开展区域耕地资源时空变化研究具有可行性和可靠性[16],能兼顾大尺度和微观尺度精细化分析需求。【拟解决的关键问题】本研究针对目前大洋洲耕地利用格局变化研究几乎空白的现状,采用中国自主开发的地表覆盖遥感制图数据产品（GlobeLand30）开展大洋洲2000—2010年耕地利用格局特征变化研究,重点揭示国家尺度、10 km网格和30 m像元尺度耕地数量、位置、利用强度变化规律以及耕地转入转出动态变化特征。1 材料与方法
1.1 研究区概况
大洋洲是世界上最小一个洲,陆地总面积约897×104 km2,约占世界陆地总面积的6%。众多岛屿散落在辽阔的太平洋海域上,气候差异明显,土地利用类型多样。大洋洲有一半以上陆地面积为干早地区。从降水量分布看,大洋洲东部群岛部分降水量远远大于西部大陆地区降水量,且有自东向西和由赤道向南北两侧减少的特点。受气候和地理条件制约,大洋洲耕地主要分布在澳大利亚和新西兰两国。作为约占大洋洲陆地面积85.7%的澳大利亚,土壤相对贫瘠,大部分地区属于干旱半干旱地带,全国31%地区年降雨量不足500 mm,39%地区年降雨量不足250 mm[17]。受恩索现象（厄尔尼诺与南方涛动）影响,降雨量低且难以预测,农作物产量对气候变化响应极为敏感,农业发展困难。早期移民开荒建城,毁掉了适应干旱环境的植被,破坏了原本还算稳固的脆弱水循环系统[18]。之后引进羊、牛和大量高需水农作物,加剧了大陆干旱且相对贫瘠的生态系统矛盾[19]。1.2 数据源与预处理
GlobeLand30 （http://www.globeland30.com）是全球首套30 m空间分辨率的地表覆盖数据产品,利用美国陆地资源卫星（Landsat）TM5、ETM+多光谱影像和国产环境减灾星（HJ-1）多光谱影像,采用像元分类、对象分类与知识规则等综合集成方法（POK, Pixel-Object-Knowledge）研制[20,21]。该数据产品采用UTM投影和WGS-84坐标系,覆盖全球80°S—80°N的陆地范围,包括耕地、水体、森林、苔原、草地、人造地表、灌木地、裸地、湿地、冰川或永久积雪等10类一级土地覆盖类别[22]。GlobeLand30产品定义的耕地主要指用来种植农作物的土地,是通过播种耕作生产粮食和纤维的地表覆盖,包括开荒地、休闲土地、轮歇地和草田轮作地;以种植农作物为主的间有零星果树、桑树或其他树木的土地;耕种三年以上的滩地和滩涂[23]。GlobeLand30耕地定义与已有的国家或国际土地利用分类系统存在差异[24]。本研究选取大洋洲为研究区,涉及2000和2010两个基准年的GlobeLand30数据。根据大洋洲的行政区划及GlobeLand30全球接图表,筛选出两期共206幅覆盖大洋洲范围的全要素数据。由于研究的数据量较大,因此采用Python语言调用ArcGIS进行批量预处理。首先利用接图表对原始数据进行批量裁剪,保证图幅之间没有重叠。其次,进行批量属性提取,得到两期大洋洲范围的所有耕地数据。再次,为了准确计算各个地类的面积,将全要素数据和耕地数据的原始UTM投影批量转换为圆柱等面积投影（Cylindrical Equal Area）,WGS-84坐标系统保持不变。
1.3 分析指标与方法
本研究采用耕地面积数量、利用强度格局和类型转换比率[25]等3个指标群,分析2000—2010年大洋洲耕地利用格局的动态变化特征。面积数量指标包括耕地面积变化率和人均耕地面积变化率。利用强度格局包括表征集约化的复种指数变化率[26]、耕地-复种指数协调度、耕地破碎度变化率[27]。类型转换比率描述了耕地和其他地表覆盖类型之间的转换,包括耕地面积转入率和面积转出率两个指标。具体指标含义及计算方法见表1。Table 1
表1
表1耕地分析指标及含义
Table 1Analysis indexes and its definitions
| 分析指标群 Analysis index group | 分析指标 Indexes | 指标含义 Index meaning | 计算公式 Calculation formula |
|---|---|---|---|
| 面积数量指标 Cultivated land quantity index | 面积变化率 Area change percentage | 指一个时间段中一个特定区域内耕地动态变化的剧烈程度, ΔCAR The intensity of the dynamic change of cultivated land in a particular area in a period of time, ΔCAR | ${{\bigtriangleup{CAR}=}{\frac{{CPCt_2}-{CPCt_1}}{CPCt_1}{\times}100\%}}$ |
| 人均耕地面积变化率Cultivated area per capita percentage | 指一个时间段人均耕地面积动态变化程度, ΔCPCR The change percentage of cultivated land per capita is the range of change of cultivated land per capita in the period t1 to t2, ΔCPCR | ${{\bigtriangleup{CPCR}=}{\frac{{CPCt_2}-{CPCt_1}}{CPCt_1}{\times}100\%}}$ | |
| 利用强度指标 Cultivated land use pattern index | 复种指数变化率 Multi-cropping index change percentage | t1-t2时间段内复种指数变化幅度, ΔMCIR ΔMCIR is the change percentage of multiple cropping index from t1 to t2 | ${{\bigtriangleup{MCIR}=}{\frac{{MCIt_2}-{MCIt_1}}{MCIt_1}{\times}100\%}}$ |
| 耕地破碎度变化率 Cultivated land fragmentation change percentage | 指t1-t2时间段内破碎度变化幅度, ΔFIR The change percentage of cultivated land fragmentation is the range of change of cultivated land fragmentation in the period, ΔFIR | ${{\bigtriangleup{FIR}=}{\frac{{FIt_2}-{FIt_1}}{FIt_1}{\times}100}\%}$ | |
| 耕地-复种指数变化协调度 Cultivated land& Multi- cropping index match degree | 指耕地面积变化率与复种指数变化率的比值, D The ratio of the dynamic change of cultivated land and dynamic change of multiple cropping index in statistical units, D | $D={\frac{\bigtriangleup{CAR}}{\bigtriangleup{MCIR}}}$ | |
| 转换比率指标 Cultivated land conversion index | 耕地转入率 Ratio of cultivated land in | 指t1-t2时间段内,其他地物类型转入成耕地的面积InCA占t2年耕地面积CAt2的比重 The ratio of other land types into cultivated land and cultivated land area in t2 during t1 to t2 | $CIn={\frac{InCA}{CAt_2}}$ |
| 耕地转出率 Ratio of cultivated land out | 指t1-t2时间段内,耕地转出成其他地物类型的面积OutCA占t1年耕地面积CAt2的比重 The ratio of the area of cultivated land changed into other land types and cultivated land area in t1 during t1 to t2 | $COut={\frac{OutCA}{CAt_1}}$ |
新窗口打开
进行耕地面积数量指标计算时,利用ArcGIS10.2中的面积制表（TabulateArea）工具完成,统计出不同尺度下的各单元耕地面积。复种指数计算在国家尺度上进行,播种面积数据来自FAOSTAT的收获面积。为与GlobeLand30耕地定义相吻合,收获面积的统计不仅包括了农作物收获面积,也包括了水果、蔬菜、茶园等的收获面积。上述分析指标计算通过Python语言编程调用ArcGIS10.2实现。
2 结果
2.1 2000—2010年大洋洲耕地面积变化
2000年大洋洲耕地面积为5 832.10×104 hm2,2010年大洋洲耕地面积为6 053.17×104 hm2,10年间耕地面积增涨220.96×104 hm2,增幅为3.79%。图1描述了2000—2010年大洋洲30 m像元尺度耕地变化空间分布。大洋洲耕地主要分布于澳大利亚从昆士兰州北部海岸沿着海岸线延伸维多利亚州和南澳大利亚州南部的狭长海岸地带,西澳大利亚西南部与塔斯马尼亚等雨水较为充沛的地区,以及新西兰沿海平原地区和巴布亚新几内亚中部区域。
显示原图|下载原图ZIP|生成PPT图12000—2010年大洋洲国家耕地变化空间分布
-->Fig. 1Distribution of cultivated land change in Oceanian countries during 2000-2010
-->
图1-b、1-c、1-d、1-e为大洋洲耕地分布集中区域,同时也是变化显著地区。其中图1-c、1-d、1-e均在澳大利亚境内,1-b为新西兰全境。图1-d处于澳大利亚西南角,耕地增加和减少无明显空间分异特征,增减相对均衡;图1-e位于南澳大利亚州南部,该地区山地与平原交错,耕地以减少为主,减少区域主要分布在南部濒临印度洋地区;图1-c地处大分水岭山脉以东,气候湿润,土壤条件相对较好,该区域10年来耕地面积以增加为主;图1-b为新西兰,耕地以减少为主,以地处北岛的新西兰最大城市奥克兰附近耕地减少尤为明显。
大洋洲有14个独立国家,其余十几个地区尚在美、英、法等国管辖之下。各国经济发展水平差异显著,耕地面积大于10 000 hm2的国家和地区仅有6个,从高到低分别为澳大利亚、新西兰、巴布亚新几内亚、斐济、新喀里多尼亚和瓦努阿图（表2）。2000—2010年,耕地面积增加国家为澳大利亚和巴布亚新几内亚,增幅分别为5.39%和2.34%。耕地面积减少的国家和地区按减幅绝对值由高到低排列为新喀里多尼亚、新西兰、斐济和瓦努阿图,变化幅度分别为-27.28%、-11.61%、-8.28%和-7.65%。
Table 2
表2
表2大洋洲主要国家2010年耕地面积、人均耕地面积、复种指数和破碎度及其10年变化率
Table 2Changing rate of cultivated land area,per capita cultivated land,multiple cropping index and fragmentation index over 2000-2010 in Oceania
| 国家 Country | 2010年耕地 面积 Cultivated land area in 2010 (×104 hm2) | 10年耕地面积变化率 Cultivated land change percentage (%) | 2010年人均耕地面积 Cultivated land per capita in 2010（hm2/person） | 10年人均耕地变化率 Cultivated land per capita change percentage (%) | 2010年耕地复种指数 Multiple cropping index in 2010 (%) | 10年耕地复种指数变化率 Multiple cropping index change percentage (%) | 2010年耕地破碎度指数 Cultivated land fragmentation in 2010 (%) | 10年耕地破碎度变化率 Cultivated land fragmentation change percentage (%) |
|---|---|---|---|---|---|---|---|---|
| 澳大利亚 Australia | 5521.24 | 5.39 | 2.51 | -8.38 | 88.52 | -0.66 | 0.11 | -47.35 |
| 新西兰 New Zealand | 458.00 | -11.61 | 1.05 | -21.63 | 11.27 | 15.25 | 0. 20 | -70.50 |
| 巴布亚新几内亚 Papua New Guinea | 57.07 | 2.34 | 0.08 | -19.68 | 183.01 | 9.99 | 0.41 | 19.07 |
| 斐济 Fiji | 11.30 | -8.28 | 0.13 | -13.47 | 122.30 | 8.69 | 1.34 | -12.18 |
| 新喀里多尼亚 New Caledonia | 3.00 | -27.28 | 0.12 | -37.98 | 40.10 | 52.60 | 0. 91 | -18.78 |
| 瓦努阿图 Vanuatu | 2.52 | -7.65 | 0.11 | -27.68 | 441.63 | 37.91 | 0.16 | -7.53 |
| 平均 Average | 1008.86 | -7.85 | 0.67 | -21.47 | 147.80 | 20.63 | 0.52 | -22.88 |
新窗口打开
10年间,大洋洲主要耕地国家人均耕地面积均有不同程度减少,其中澳大利亚人均耕地面积由2000年的2.74 hm2,减少为2010年的2.51 hm2,减少8.38%。新西兰人均耕地面积由2000年的 1.34 hm2减少为2010年的1.05 hm2,减幅为21.63%。巴布亚新几内亚、斐济、新喀里多尼亚和瓦努阿图的人均耕地面积分别减少19.68%、13.47%、37.98%和27.68%。
2.2 2000—2010年大洋洲耕地利用强度变化
10年间,大洋洲平均复种指数提高20.63%。按照复种指数变化率由高到低排列为新喀里多尼亚提高52.60%,瓦努阿图增加37.91%,新西兰增加15.25%,巴布亚新几内亚提高9.99%,斐济增加8.69%,澳大利亚略微降低0.66%。耕地-复种指数变化协调度指标描述10年耕地面积变化率与复种指数变化率的协调程度（表1）。协调度在-1和0范围中,值越接近0表明耕地面积变化与复种指数变化协调度越低,耕地利用强度越高,即收获面积增长对耕地面积增长的依赖度越低。将六国按照象限图分为耕地复种双增型、耕增复减型、耕减复增型3种类型,本大洲无耕地复种双减型。耕地复种双增型国家仅为巴布亚新几内亚,协调度为0.23,表明该国耕地利用强度呈现温和增长。澳大利亚为唯一耕增复减型国家,协调度为-8.17,表明该国耕地利用强度降低,单位耕地新增未带来相应的收获面积提升。其余四国均为耕减复增型国家,协调度由高到低为瓦努阿图-0.20,新喀里多尼亚-0.52,新西兰-0.76,斐济-0.95。大洋洲岛国受地形、社会经济与气候影响,农业生物技术与机械化程度对耕地利用强度影响极小,耕地面积与收获面积变化弹性较大。
耕地破碎度对耕地利用效率和格局影响较大[28]。10年间,大洋洲耕地破碎度平均减少了22.88%。新西兰耕地破碎度减少70.50%,耕地规模化利用程度提高迅速。巴布亚新几内亚是唯一破碎度增加的国家,增幅19.07%。由图2可见,耕地破碎度增加明显的区域主要集中在耕地资源较为集中的澳大利亚西南角和南澳大利亚州南部临海区域。这部分耕地大多处于地中海气候区和常湿温暖气候区,10年间耕地有缓慢减少趋势。耕地破碎度降低明显的区域主要集中在澳大利亚东部墨累达令河流域上游,10年间该区域耕地增加明显。新西兰北岛北部耕地破碎度降低,同时耕地面积也减少较快。整体而言,澳大利亚耕地破碎度变化主要是由于土地利用方式改变导致的耕地面积变化,耕地面积增加的区域破碎度降低,耕地面积减少的区域破碎度提高,这种变化对区域单位面积耕地的集约化程度影响较小。新西兰北岛在耕地面积减少的同时促进了地块间的合并,扩大了农场规模,提高了集约化程度,促使破碎度迅速降低。
显示原图|下载原图ZIP|生成PPT图22000—2010年大洋洲国家10 km×10 km网格耕地破碎度变化分布
-->Fig. 2Change ratio of fragmentation index of cropland in 10 km×10 km grid in Oceanian countries during 2000-2010
-->
2.3 2000—2010年大洋洲耕地转换特征
2000—2010年,大洋洲耕地转出面积最大的地类为草地,共计转出528.21×104 hm2,占2000年大洋洲耕地总面积的9.11%（图3）。耕地转出总面积最高的国家是澳大利亚,共计630.25×104 hm2耕地转为其他地类,占本国耕地总面积12.03%,其中转为草地471.27×104 hm2,占所有转出地类的74.77%。新西兰共有12.77×104 hm2耕地转出,占本国耕地总面积的14.04%,其中转为草地最多,约57.03×104 hm2,占所有转出地类的78.37%。巴布亚新几内亚的转出耕地面积为11.96×104 hm2,占本国耕地面积的21.44%,其中转为森林最多,约8.93×104 hm2,占所有转出地类面积的74.65%。澳大利亚、新西兰、新喀里多尼亚和斐济耕地主要转为草地,巴布亚新几内亚和瓦努阿图耕地主要转为森林。大洋洲热带雨林气候区和地中海式气候区耕地转为森林比重最高,其他气候区主要转出目标为草地。
显示原图|下载原图ZIP|生成PPT图32000—2010年国家尺度耕地转出和转入类型比例
-->Fig. 3Conversion ratio of cultivated land in 2000 and 2010 in national scale
-->
10年间,大洋洲耕地转入面积最大的地类仍为草地,共计转入555.58×104 hm2,占2010年大洋洲耕地总面积的9.23%（图3）。耕地转入总面积最高的国家是澳大利亚,共计912.47×104 hm2其他地类转为耕地,占本国耕地总面积的16.53%,其中由草地转为耕地最多,约544.94×104 hm2,占所有转入地类的59.72%。新西兰共有12.59×104 hm2其他地类转为耕地,仅占本国耕地总面积的2.75%。巴布亚新几内亚转入耕地面积13.26×104 hm2,占本国耕地总面积的23.24%,其中由森林转为耕地最多,约7.76×104 hm2,占所有转入地类的58.53%。澳大利亚和新西兰转入耕地的主要来源为草地,巴布亚新几内亚、新喀里多尼亚、斐济和瓦努阿图的耕地主要由森林转入。大洋洲所有气候区由草地转为耕地占比最大。
由前文可知10年间澳大利亚和巴布亚新几内亚耕地增加,其他四国耕地减少。通过10年间耕地增减平衡测算可知,澳大利亚耕地净增加贡献最多地类为灌木、草地和森林,分别净增加166.45×104 hm2、73.67×104 hm2和31.01×104 hm2。巴布亚新几内亚净增加耕地主要来源草地,约为3.32×104 hm2。新喀里多尼亚耕地净减少贡献最多地类为草地、森林和灌木,分别为0.64×104 hm2、0.34×104 hm2和0.15×104 hm2。新西兰和斐济耕地净减少主要由于耕地转草地,分别净减少50.24×104 hm2和1.18×104 hm2。瓦努阿图耕地净减少主要由于耕地转为森林和灌木,分别为0.16×104 hm2和0.06×104 hm2。从大洋洲整体来看,尽管耕地与草地之间转换面积最大,但基本保持增减平衡,而灌木地转为耕地对耕地净增贡献最大,共净增耕地约165.03×104 hm2。
3 讨论
土地利用/覆被变化影响气候变化,同时人类为了适应气候变化也会改变土地利用方式[29]。在大洋洲,2000—2010年耕地与草地之间的转换是最常见的地类转换,特别在常湿温暖区和草原气候区交界的地方转换尤为频繁。10年间有143.01×104 hm2退耕还林,占全洲耕地总面积的2.47%,但528.21×104 hm2耕地转换为草地,占全洲耕地总面积的9.11%。2010年草地转为耕地占耕地总面积约为9.23%。10年里耕地与草地间的增减基本平衡,少量净增耕地贡献主要来源于灌木、森林和草地。耕地与草地之间相互转换频繁可能原因是澳大利亚推行将饲养牲畜和谷物生产有机结合起来的混合农业,农场内的土地交替种植小麦、牧草或休耕,充分保持麦田的肥力。澳大利亚的混合农业主要分布在国土东南部的墨累-达令盆地和西南角,处于常湿温暖气候区和地中海式气候区。除采取混合农业模式保护土壤肥力,2009—2010年,澳洲有1 900×104 hm2的农场土地采用免耕耕作[30]的方法种植农作物。保护性农业的目标是提高长期生产率、增加收益产出和保证粮食安全,但短期而言对产量增加有一定负面影响[31]。澳大利亚按自然区划可分为东部山地、中部平原和西部高原3个区域。东部山地北起约克半岛, 南到塔斯马尼亚岛;中部平原北起卡奔塔利亚湾, 纵贯大陆中部, 向南延伸到墨累河河口;西部高原大部分为沙漠和半沙漠地区。西澳大利亚寒流和南回归线上的副热带高压稳定的控制了中部西部大部分地区,空气干燥,降水稀少。东部狭长的暖湿地带则是东澳大利亚暖流带来的福泽,东西中间通过大分水岭山脉隔绝开来。因此,澳大利亚耕地主要集中在水资源相对充沛的大分水岭山脉以东、南澳大利亚州和新南威尔士州南部、维多利亚州、塔斯马尼亚州北部和东部、西澳大利亚州西南角。10年间澳大利亚耕地新增区域主要集中在大分水岭山脉以东,墨累达令河流域上游。受2001—2008年大旱影响,澳大利亚在墨累达令河流域附近（图 1-c和图1-e）加大了垦殖力度,促使新增耕地朝可能水源地集中。2010年全国垦殖率约为7.94%,相比2000年提高了5.39%。墨累河流域水能资源主要集中在干流上游及其支流。由于河流流经的大部分地区为干旱地区,流域水资源开发的主要目的是灌溉和供水[32],在此流域种植着澳大利亚90%的灌溉作物。筑坝,大量灌溉用水和城市供水让澳洲大陆本就脆弱的水源不堪重负。灌溉导致土壤含盐量过高[33],破坏了湿地,大片土地变得不适合耕种。可见,澳大利亚耕地分布受到气候及水资源分布影响十分明显。
新西兰位于太平洋西南部,山地和丘陵占其总面积75%以上,属温带海洋性气候,四季温差不大,平均气温夏季20°,冬季12°。畜牧业产品出口收入占出口总收入60%以上。农业高度机械化,主要农作物有小麦、大麦、燕麦、水果等。粮食不能自给,需要从澳大利亚进口。全大洲耕地减少比较集中的区域主要在新西兰北岛,全国最大城市奥克兰附近。2000年全国垦殖率为17.66%,2010年下降了11.61%。近10年来,新西兰政府加大了奥克兰城市开发力度,并出台了《奥克兰统一规划草案》（Proposed Auckland Unitary Plan）,预计到2030年,奥克兰将开发相当于两个汉密尔顿规模的住宅区域。10年间新西兰耕地空间格局的变化对气候变化的响应不如澳大利亚那么敏感,对城市化政策的响应更为明显。
大洋洲不同国家和地区在耕地利用格局变化驱动机制上存在差异,这种差异不仅体现在气候变化、区域政策和社会经济等因素上,全球化背景下的国际农产品贸易也会影响耕地利用格局变化。未来将进一步基于GlobeLand30对大洋洲主要国家耕地格局变化驱动机制展开研究,并基于远程耦合综合框架[34]分析中国与澳大利亚的远距离农业社会经济活动和环境相互作用机制,帮助理解和权衡本地粮食需求与远程粮食生产供给的复杂关系。
需要说明的是,虽然GlobeLand30数据产品精度总体高,但该套数据产品在使用中也存在一定局限性。一方面,GlobeLand30耕地定义遵循FAO标准,不仅包括常规农作物,也包括牧草种植地,果树、茶园、咖啡园等灌木类经济作物,其耕地定义不同于已有国家土地利用分类系统或者中低空间分辨率的全球地表覆盖数据产品,如GlobeCover、UMD-GLC和BU-MODIS,使得GlobeLand30和已有的数据产品之间不具有可比性。另一方面,澳洲不同区域、不同国家的农业土地利用多样化,农田地块大小、种植制度和种植模式存在较大差异性,基于遥感影像的光谱和时相信息提取的耕地精度不可避免地存在区域差异,在一定程度上影响本研究的分析结果。
4 结论
从耕地面积变化来看,大洋洲2000—2010年耕地面积增加3.79%,耕地面积增幅最大的国家为澳大利亚,增幅5.39%。新增最多的区域主要集中在澳大利亚大分水岭山脉以东墨累-达令河流域上游。耕地面积减少的区域主要在新西兰北部岛屿,澳大利亚东部沿海和巴布亚新几内亚东部岛屿。从利用格局变化来看,主要国家复种指数平均增加20.63%,耕地破碎度平均减少22.88%。耕地面积-复种指数协调度弹性较大。
从转换类型来看,全大洲耕地与草地之间转换面积最大,但净增加耕地贡献最大的是灌木地,净增165.03×104 hm2。
The authors have declared that no competing interests exist.
