1. 引言

土地是人类赖以生存的基本自然资源，同时也是人类战略性的经济资源。土地资源保护问题直接关系到国家生态安全，也是土地可持续利用研究的一个重要问题 [1] 。随着全球城市化及经济快速发展和人口不断增长，对土地的需求也不断上升 [2] ，导致土地的深度和广度利用加快 [3] ，从而加重了人地矛盾，生态环境问题日益突出。因此，分析和评价土地利用现状，特别是对生态脆弱地区的可持续发展，以及土地保护，协调与经济发展的关系具有重要意义 [4] 。

蒙古国是世界最大的草原国家之一，拥有独特而珍贵的草原生态系统，具有多重生态服务功能。国内外诸多学者围绕土地利用现状做了大量研究，现阶段蒙古国在土地利用研究理论、土地利用研究方法及技术等方面已经形成了较丰富的研究成果。如Amarsaikhan D (2005) [5] 、Chinbat B (2006) [6] 、巴图赫斯格(2008)等人，主要利用遥感和地理信息系统分析土地利用变化及其驱动因子。Altansukh Ochir等人与国内外学者合作分析了蒙古国1990~2010年的土地覆盖变化 [7] 。同时，学者们也重点研究了蒙古国单一土地利用类型的动态变化，如Davydov Oyuntugs (2013) [8] 、Gantumur [9] 等人，分别采用RS、元胞自动机(CA)模型和人工神经网络(ANN)结合方法研究了不同时期的蒙古国首都乌兰巴托城市扩张的时空演变；Tuvshintogtokh [10] 、Myangan [11] 、Ganzorig [12] 、Erdenee [13] 等人对草地、水域、林地和耕地等单一土地利用类型的变化。此外，众多学者研究发现气候 [14] 及土地产权的变化，干旱、火灾等自然灾害 [15] ，加上森林砍伐、过度放牧等人为因素均是导致蒙古国土地利用、资源与环境利用发生重大变化的主导原因之一 [7] 。然而，以蒙古国单一省份为研究区的土地利用现状评价研究相对较少。因此本文立足于蒙古国国情及土地资源保护目的，利用计量地理学的多样化指数、集中化指数、土地区位指数及土地利用程度等方法，分析蒙古国中央省土地利用现状水平，明确区域土地利用结构，揭示土地资源区域组合和空间配置的合理性，探究土地开发利用的方向和重点，为区域土地合理利用提供一定的参考依据。

2. 研究区概况与数据来源

2.1. 研究区概况

蒙古国中央省(Tuv Province)位于蒙古国中部(图1)，该省将首都乌兰巴托市从四面包围，北接色楞格省，东面与肯特省，西面、南面与东戈壁省、中戈壁省、前杭爱省、布尔干省相接，地理位置46˚15'N~49˚10'N，104˚10'E~109˚2'E。中央省属于温带大陆性气候，季节变化明显，夏季炎热、冬季寒冷、雨热同期、无霜期短；海拔在834~2701 m之间，地势东北部高西北和东南部低，隶属蒙古国东部肯特山脉的高原，自然资源丰富。中央省总面积7.40万km²，其中草地面积占比为83.35%，林地和耕地面积占比分别为14.13%和1.89%。中央省下辖1个市和26个县(苏木)，省会为宗莫德市(表1)，距离乌兰巴托市43 km。

2.2. 数据来源与处理

研究区行政界线数据从GADM (Database of Global Administrative Areas)全球县级尺度矢量(shp)底图数据中提取；土地利用分类数据来源于欧洲航天局2021年全球土地利用数据 (https://worldcover2021.esa.int/downloader)，空间分辨率为10米。

本文参考中国科学院资源环境科学与数据中心土地利用遥感监测数据集分类，按研究需要将土地利用类型划分为耕地、林地、草地、水域、建设用地、未利用地等6类。利用EXCEL计算计量地理模型，并借用ArcGIS10.8软件制作土地利用空间分布图。

3. 研究方法

本文采用土地利用多样化指数、土地利用集中化指数、土地利用区位意义和土地利用程度等方法评价中央省土地利用现状。

3.1. 土地利用多样化指数

计算土地利用多样化指数的目的是为分析区域内各类土地利用类型的齐全程度或多样化状况 [16] ，可计算出区域各类土地利用类型的均衡程度。本研究采用吉布斯–马丁多样化指数(Gibbs-Martin)来计算中央省土地利用多样化程度 [17] ，其中G-M的值为{0, 1}，计算结果接近于1，表明区域内土地利用类型越多样，反之计算结果越接近0，土地利用类型越单一，其计算公式如(1)所示：

$\text{G-M}=1-\frac{{\displaystyle {\sum }_{i=1}^n{x}_i^2}}{{\left({\displaystyle {\sum }_{i=1}^n{x}_i}\right)}^2}$(1)

式中，G-M为土地利用多样化指数，n为土地利用类型，x_i代表i类土地面积。按公式(1)计算获得中央省各苏木土地利用多样化指数，并利用ArcGIS10.8自然间断点分级法，将土地利用多样化指数按照G-M > 0.30为高多样化区、0.04 < G-M ≤ 0.30为中多样化区和G-M ≤ 0.04为低多样化区3个等级进行分级归类。

3.2. 土地利用集中化指数

集中化指数的计算有利于揭示区域土地利用类型的集中化程度 [18] ，其计算公式如下：

$I_i=\frac{A_i-R}{M-R}$(2)

式中，I_i为集中化指数；A_i为第i个区域各地类累计面积百分比之和；R为高一级区域各地类累计面积百分比之和，本文中中央省为高一级区域；M为土地集中分布时的累计面积百分比之和。本文土地利用类型分为6类，因此M值为600，经计算R值为579.98，R值是衡量集中化程度的基准 [19] 。集中化指数越大，表明该区域土地利用的集中化程度越高。

利用公式(2)计算出中央省下辖各苏木土地利用集中化指数，并利用ArcGIS10.8重分类法，对中央省土地利用集中化指数进行分级归类，I > 0.40为高集中化区、0 < I ≤ 0.40为中集中化、I ≤ 0为低集中化区3个等级。

3.3. 土地数量结构区位意义

土地数量结构区位意义可用土地区位指数来计算。土地区位指数是指某一区域的某种土地利用类型在高一级区域内的相对集中程度 [20] [21] ，如果Q_i> 1，则表示该区域第i类地类区位意义显著，否则相反。土地区位指数公式为：

$Q_i=\frac{f_i}{{{\displaystyle \sum }}^{}{f}_i}/\frac{F_i}{{{\displaystyle \sum }}^{}{F}_i}$(3)

式中，Q_i为区位指数；f_i代表第i类土地利用类型面积；F_i代表高一级区域第i类土地利用类型面积。

3.4. 土地利用程度综合分析

土地利用程度综合指数是针对一个具体年份内所有土地利用类型构成的综合反映，它是通过对土地利用程度的分级，量化人类活动对土地的影响程度，定量描述该区域土地利用的综合水平和变化趋势，是衡量区域土地利用情况的一个重要指标。土地利用程度指数的计算需要将土地利用类型按照人类活动影响下的自然平衡状态分级 [22] ，计算公式为

$L=100\times {\displaystyle {\sum }_{i=1}^n\left(B_i\times C_i\right)}$(4)

式中，L为土地利用程度综合指数，L的值为 $\left\{100,400\right\}$；n为分级数量；B_i为区域第i级土地利用程度分级指数；C_i为第i级地类占土地总面积的百分比。本文参考并借鉴现有研究成果 [21] [22] [23] [24] ，将土地利用程度指数划分为4等级(表2)。

4. 结果分析

4.1. 土地利用多样化分析

中央省土地利用多样化指数平均值为0.28，与理论最大值0.83相差很大，说明中央省土地利用类型总体多样化程度较低，土地利用类型较为单一。中央省下辖各苏木多样化指数较低，介于0~0.5之间，巴图孙部尔、巴彦庆大毛尼、宝日淖尔、鄂尔顿、吉日嘎朗图、孟根毛日图、孙布尔、车格勒等苏木等8个苏木的多样化指数高于全省总体水平，其中，巴图孙部尔和鄂尔顿苏木多样化指数最高，为0.50，表明该两个苏木土地利用类型相对较丰富，其余18个苏木土地利用多样化指数均低于平均指数，阿拉坦宝力格、巴彦和巴彦吉日嘎朗3个苏木为土地利用多样化指数最低，为0.01，说明这些苏木土地利用类型较单一，用地主要集中在某一用地类型(表3)。

按土地利用多样化等级划分标准，中央省高多样化区的苏木有7个，占研究区苏木总数的26.92%；中多样化区的苏木有9个，占34.62%，低多样化区的苏木有10个，占38.46% (图2)。研究区土地利用多样化情况在空间上呈现出东北部高于西北部、北部高于南部的分布规律。高多样化区主要分布于东北部、西北部少部分区域；中多样化区主要分布于西北部及东端的巴彦德力格尔和阿日胡斯图、西端的鄂尔顿桑土苏木；南部区域均为低多样化区。分析土地利用多样化状况的这种空间变化特点，与研究区地形地貌特点有一定联系，全省地势从东北向西南逐渐降低，东北部海拔最高，以肯特山脉为主，以森林和草原资源较为丰富，土地利用类型主要以林地和草地为主；西北部海拔相对较低，地势平坦，距离乌兰巴托市较近，市郊少面积的耕地分布于草地之中；南部区域以平原草地为主，草地占总面积的90%以上，其他用地类型较少。基于以上自然本底条件及市场距离的影响，土地利用多样性呈现出东北部高于西北部，北部高于南部的特点。

4.2. 土地利用集中化分析

中央省各苏木土地利用集中化指数介于−1.50~0.97之间，下辖苏木间集中化指数差异较小(表4)。其中孙部尔、吉日嘎朗图、孟根毛日图、宝日淖尔、车格勒、巴彦庆大毛尼和孙布尔等8个苏木的土地集中化指数为负值，鄂尔顿苏木的最低，为−1.5，这8个苏木的土地利用集中化程度低于全省平均水平；其余18个苏木的土地集中化指数为正值，其中巴彦吉日嘎朗苏木的土地利用集中化指数最高，为0.97。

中央省各苏木土地利用集中化程度差异性明显，集中化程度总体上呈现为较高的状态。按高中低集中化区三个等级划分，中央省高集中化区的苏木16个，占研究区苏木总数的61.54%；中集中化区的苏木为2个，占7.69%，低集中化区的苏木有8个，占30.77% (图3)。研究区土地利用集中化程度在空间上呈现出南部高于北部的特点。高集中化区主要分布于南部区域；中集中化区分布于西北部乌戈塔拉和巴彦朝歌图2个苏木；北部区域为低集中化区。究其原因，草地占南部各苏木绝大部分土地面积，土地利用类型单一；因自然、社会原因使得北部土地类型变得更为多样化，促使土地集中程度降低。

4.3. 土地区位指数分析

中央省下辖各苏木土地类型的区位指数均不同(表5)。耕地区位指数大于1的苏木共有12个，其中阿日噶狼图、阿日胡斯图、巴彦庆大毛尼、巴彦朝歌图、宝日淖尔位于首都乌兰巴托市周围，满足首都城市发展中对于农副产品的需求，鄂尔顿桑土、吉日嘎朗图、孙布尔、车格勒、乌戈塔拉、扎姆尔和巴彦杭盖主要分布在全省地势较低的西和西北部。林地具有区位意义的苏木共5个，孟根毛日图、鄂尔顿和巴图孙部尔等3个苏木区位意义最显著，Q_i值分别为4.30、3.58和2.81，宝日淖尔和吉日嘎朗图苏木区位意义较显著，Q_i值为1.49和1.35，这5个苏木是中央省重要的森林资源区。草地区位意义显著的苏木有18个，占研究区苏木总数的69.23%，其中巴彦吉日嘎朗苏木的草地区位意义最显著，Q_i值为1.20，其他17个苏木区位意义较显著，Q_i值为1.04~1.18，剩余8个苏木介于0.46~0.99，草地区位意义不明显。草地占中央省总面积83.35%，分布范围广，畜牧业是国民经济主要来源，因此草地区位意义显著苏木数量为最多。伦格、扎姆尔和宝日淖尔苏木的水域区位意义最显著，Q_i值 > 2，阿日胡斯图、巴彦、巴彦查干、布仁、孟根毛日图、温杜尔席勒图等6个苏木的水域区位指数Q_i分布于1.14~1.56之间，较为显著。建设用地具有区位意义的苏木共9个，巴彦、鄂尔顿、伦格和温杜尔苏木区位意义指数介于1.35~1.74，具有较显著的区位意义，其他5个苏木Q_i> 2.90，区位意义很显著。斯日古楞苏木建设用地区位条件最优势，这与中央省省会宗莫德市位于该苏木有关系。未利用地具有区位意义的苏木共6个，占研究区苏木总数的23.08%，其中巴彦乌力吉、巴彦查干、德力格尔汗、扎姆尔4个苏木区位意义最为显著，Q_i均 > 2，布仁和伦格2个苏木区位意义较为显著，Q_i分别为1.60和1.59，具有一定的后备资源。

4.4. 土地利用程度综合分析

中央省土地利用程度不高，综合指数为201.53，高于全省平均水平的苏木有14个，占研究区总苏木总数的53.85%，低于全省平均水平的苏木有12个。下辖各苏木的差异性不明显，车格勒苏木土地利用程度最高，综合指数为216.95，德力格尔汗苏木土地利用程度最低，综合指数为198.74。本文根据土地利用程度综合指数计算结果(表6)，以L > 300为高土地利用程度区，200 < L ≤ 300为中土地利用程度区，100 < L ≤ 200为低土地利用程度区，将全省土地利用综合程度划分为3等级区(图4)。中央省农牧区，牧户居所为便于迁移的蒙古包，住宅等建设用地少，土地利用受人类活动干扰程度低，因此该省没有高土地利用程度区的分布。由于耕地主要分布在西北地区，人类活动比其余地区剧烈，所以中土地利用程度区主要分布于西北部地区；低土地利用程度区主要分布在人类生产活动较少的南部和东北部林草水用地区域。

5. 结论

1) 蒙古国中央省土地利用多样化指数较低，土地利用类型单一，高多样化区主要分布于东北部、西北部少部分林地草地分布区域，西北部为中多样化区，南部地区为低多样化区，草地分布广泛。

2) 中央省土地利用类型的集中化程度较高，各苏木集中化差异性较为明显，空间上呈现出南部高于北部的特点。

3) 中央省各苏木各类土地类型的区位意义指数均不同，其中耕地区位意义较突出的苏木主要分布在首都乌兰巴托周边和全省地势较低的西和西北部。全省草地面积最大，畜牧业是国民经济支柱产业之一，因此草地区位意义显著苏木数量为最多。

4) 中央省土地利用程度不高，全省土地利用程度综合指数为201.53，下辖各苏木的差异性不明显，空间分布上西北部土地利用程度高于南部和东北部。

基金项目

国家自然科学基金项目(42161044)；内蒙古自治区自然科学基金项目(2020MS04019)和(2019MS04014)；内蒙古师范大学高层次人才科研启动项目(2018YJRC039)。

NOTES

^*第一作者。

^#通讯作者。

参考文献