鄂尔多斯位于内蒙古自治区西南部，地处鄂尔多斯高原腹地，介于北纬37˚35′24″~40˚51′40″，东经106˚42′40″~111˚27′20″之间，东西长约400 km，南北宽约340 km，总面积8.7万km²。研究区西、北、东三面被黄河环绕，境内湖泊、草原、沙漠、丘陵、高山等地质景观齐全。鄂尔多斯地区有两大沙漠，位于北部的库布齐沙漠和横跨内蒙古鄂尔多斯和陕西榆林的毛乌素沙地，大部分位于鄂尔多斯境内(图1)。研究区属于半干旱大陆性气候，年平均气温在5.3℃~8.7℃之间，极端最高气温可达37.4℃~38℃ [ 12 ] 。南部地区比中北部气温更高。此外，研究区东部年平均降水量为300至400毫米，而西部为200至300毫米，降水主要集中在7月和9月。鄂尔多斯市独特的地理位置与地形条件构成了我国北部生态屏障的重要组成部分。

图1. 研究区示意图

本文采用的CLCD LUCC数据集是基于Google Earth Engine平台的多传感器观测数据，该数据由武汉大学杨杰等人利用随机森林分类器处理得到，空间分辨率为30米，包含1985 + 1990~2020中国逐年土地覆盖信息，总体准确率达80% [ 13 ] 。为分析鄂尔多斯2000年以来草地退化和恢复状况，本研究基于CLCD LUCC数据集2000~2020研究时段，详细刻画研究区LUCC和草地退化回复时空规律。

鄂尔多斯地区CLCD LUCC土地利用数据共有8个一级类，包括农田、森林、草地、裸地、湿地、水体、冰雪、不透水面。根据研究区景观特征，本研究将CLCD LUCC数据按照一级地类重新分类如下：农田、森林、草地、裸地、水域和不透水面(见表1)。社会经济指标数据来自2000~2020年鄂尔多斯地区统计年鉴(http://www.zgtjcbs.com)，部分数据来自如内蒙古统计年鉴(http://tj.nmg.gov.cn/)。

表1. CLCD LUCC土地利用分类

在2000~2020年期间，鄂尔多斯市共有17473.02 km²的国土面积发生变化，约占研究区国土面积的20% (包含6种土地利用/覆被类型) (图2(d))。研究区LUCC变化较为明显，但地类种类和空间分布并无明显演替。裸地和草地仍是鄂尔多斯主要土地利用类型，分别约占研究区总面积的13.46%和75.65% (2020年)，农田约占少于9.53%，是研究区的第三大地类，不透水面、森林和水域约占研究区总面的1.35% (表2)。各个地类具有明显的空间分布特征：其中，裸地主要分布在杭锦旗；草地则主要分布在研究区东部和东南部旗县；农田空间分布较为集中，主要分布在“黄河几字弯”(图2(a)~(c))。

表2. 2000、2010、2020年鄂尔多斯市土地利用/覆被面积统计表

图2. 鄂尔多斯市土地利用现状及变化图

如图3所示，2000年以来的三个阶段(2000~2010，2010~2020，2000~2020)，鄂尔多斯土地利用变化主要包括四个过程：分别为草地面积增加、裸地面积减少、农田面积增加和不透水面面积扩张。同时，在过去二十年期间，研究区森林面积持续增加，水域面积轻微增加。

在2000~2020年期间，研究区草地净面积增加最多，为5855 km²，占研究区总面积的6.75% (图3)，增加的草地面积主要来源于裸地和农田，面积分别为8419.5 km²和2882.3 km²(表4)，增加的草地在整个研究区均有大面积分布。在2000~2010和2010~2020年两个研究时段期间，转入的草地同样来源于裸地和农田，变化趋势保持一致(见表3)。

裸地是2000年以来变化最大的土地类型，2000~2020年期间裸地净面积减少为6947.2 km²，占总面积的8.01% (图3)，减少的裸地主要是转变为草地和农田，面积分别为8419.5 km²和157.7 km²，其次是林地、水域和不透水面(表3)，减少的裸地面积在研究区各个旗县均有分布。在2000~2010和2010~2020年期间，裸地变化趋势与2000~2020年整体变化趋势保持一致(见表3)。

2000~2020年期间，不透水面和农田净面积分别增加534 km²和453 km²(图3)，占研究区总面积的0.62%和0.52%。其中不透水面面积不断扩张，2010~2020年农田面积增幅比例大于2000~2010年的农田增幅(图3，表3)。

森林面积在三个研究时段(2000~2010, 2010~2020, 2000~2020)增长趋势保持一致，增幅差别较小(图3)。2000~2020年，森林净面积增加21 km²，增加的森林面积主要来源于草地和农田，面积分别为11.4 km²和8.3 km²，其空间分布主要在杭锦旗、达拉特旗等北部地区。水域净面积为84 km²，占总面积的0.10% (图3)。

图3. 鄂尔多斯市土地利用/覆被净变化(2000~2020, 2000~2010, 2010~2020)

表3. 鄂尔多斯市土地利用/覆被转移矩阵(2000~2010, 2010~2020)

2000~2020年期间，研究区草地生态状况趋于良好，草地恢复面积(11346.1 km²)大于草地退化面积(5490.8 km²)。2000~2010和2010~2020两个研究时段内，草地恢复趋势与2000~2020年整体保持一致，草地恢复态势明显。恢复面积分别为7899.7 km²和6162.9 km²，退化面积为3614.7 km²和4592.9 km²。2000~2010年研究时段内，在空间分布特征是草地恢复主要集中在鄂托克旗，鄂托克前旗，杭锦旗，草地退化在鄂托克旗和鄂托克前旗；2010~2020年研究时段内，草地恢复与退化空间分布与2000~2010年保持一致(见图4)。

土地利用是自然生态和社会经济等多种因素影响下的产物，融合了人类活动和生态系统相互作用的结果。本文主要从社会经济的方面来分析鄂尔多斯市土地利用的变化驱动力，以找出对其土地利用产生影响的主要驱动因素。

由鄂尔多斯市2000~2020年的统计年鉴中提供的数据和鄂尔多斯市在2000~2020年的土地利用变化特点，运用主成分分析的相关性筛选出以下15个因素X₁生产总值(亿元)、X₂第一产业总值(亿元)、X₃第三产业总值(亿元)、X₄人均GDP(万元)、X₅城市化率(%)、X₆农牧民纯收入(万元)、X₇城镇居民人均收入(元)、X₈年末人口(万人)、X₉年末牲畜总量(万只)、X₁₀城镇总人口(万人)、X₁₁建筑业(万元)、X₁₂牛肉产量(吨)、X₁₃羊肉产量(吨)、X₁₄粮食产量(吨)、X₁₅玉米(吨)。对15个影响因素的标准化后使用SPSS分析软件进行主成分分析处理数据，获得影响区域土地利用变化的主要驱动力。

在因子分析前，首先要进行 KMO和巴特利特检验。检验结果(见表4)显示，KMO值为0.830，显著性远小于0.05。这表明变量之间存在强相关性，适合进行因子分析，为接下来的主成分分析提供了基础。

图4. 2000~2010、2010~2020、2000~2020年草地退化与恢复

表4. KMO和巴特利特检验结果

由原始数据得到相关系数矩阵如表5所示，可知X₁、X₂、X₃、X₄、X₅、X₆、X₇、X₈、X₁₄、X₁₅和X₁、X₂、X₃、X₄、X₅、X₆、X₇、X₈、X₁₄、X₁₅之间的相关系数大于 91%，表明这些因素之间存在极强的正相关关系，同时X₁₀和X₅，X₁₁和X₁、X₂、X₃、X₄、X₅、X₆、X₇、X₈、X₁₄，X₁₂和X₁₃，X₁₃和X₅、X₈、X₁₂、X₁₄、X₁₅之间的相关系数大于 90%，表明这些因素之间有较强的相关性。它们的相关性可以满足主成分分析的要求。

表5. 相关性矩阵

主成分特征值及贡献率的结果如表6所示，取主成分特征值大于1，第一主成分的特征值为13.212、贡献率88.078%、累计贡献率88.078%，说明鄂尔多斯市土地利用结构的变化与社会经济发展密切相关；第二主成分的特征值为1.138、贡献率7.589%、累计贡献率95.668%，说明牲畜数量变化、牛肉产量等民生发展是影响鄂尔多斯市土地利用结构变化的一个重要因素。因此，影响鄂尔多斯市的主要驱动力因子为：经济社会发展、民生发展，同时，草地相关政策对研究区域的土地利用有重要的影响作用。

表6. 主成分特征值与主成分载荷表

根据因子得分系数和原始变量的标准化值，可得第一主成分与第二主成分的表达式：

Y 1 = 0. 267X 1 + 0. 268X 2 + 0. 267X 3 + 0. 27X 4 + 0. 271X 5 + 0. 265X 6 + 0. 269X 7 + 0. 27X 8 + 0. 191X 9 + 0. 247X 1 0 + 0. 267X 11 + 0. 226X 12 + 0. 248X 13 + 0. 267X 14 + 0. 267X 15

Y 2 = − 0. 211X 1 − 0. 187X 2 − 0. 2 0 7X 3 − 0. 18X 4 + 0.0 26X 5 − 0. 185X 6 − 0. 183X 7 + 0.00 7X 8 + 0. 636X 9 − 0. 112X 1 0 − 0. 158X 11 + 0. 43X 12 + 0. 347X 13 + 0. 126X 14 + 0. 12X 15

以对每个主成分的特征值所占主成分的总特征值之和的比例作为权重，计算出主成分综合得分，即计算主成分的得分公式为：

Y = λ 1 Y 1 + λ 2 Y 2 = 0. 921Y 1 + 0.0 79Y 2

经过计算得到社会经济驱动力的综合得分表7。

表7. 2000~2020年驱动力综合得分表

由表7，在2000年到2020年间，社会经济的驱动力综合得分为正数，逐年升高，这是由于鄂尔多斯市的经济结构趋于合理，对土地利用的影响逐渐变大。这反映了社会经济因素对鄂尔多斯市土地利用的驱动力不断变大。

过去几十年，鄂尔多斯经历了经济迅速发展和自然生态保护的双重成就，土地利用政策也在多年社会发展过程中不断调整和变化，从大面积开垦“以粮为纲”到现在坚持“生态优先，绿色发展”的发展之路，研究区社会经济发生深刻变化，土地利用强度和方式都发生了显著变化。

2000年以来，覆盖到研究区的各项重大生态工程措施见表8。其中季节性轮牧、禁牧、休牧、舍饲、科学饲养等措施促使年末牲畜总量和乡村人口总量下降，减少草地承载压力，从而促进草地自然恢复。生态移民等措施对于严重退化的地区的居民，通过提高生态补偿的方式鼓励牧民迁移，既减轻严重退化区域草地压力，又提升农牧民收入，鼓励牧民向城镇迁移。该措施一方面遏制了草地退化态势，提高草场质量；另一方面改变了人口结构，鼓励牧民迁入城市，也拉大了研究区农村人口和城镇人口数量差距，促进研究了城镇化进程，间接导致草地退化。植树造林、沙地治理的各项措施不但直接增加了植被面积，也改善改变了农村发展理念、发展模式及生产生活方式，增加了农村就业岗位，提高了农牧民收入，加快了民生发展。综上所述，各项生态功能过程措施不仅直接作用于草地退化区域，更通过作用于人类活动，如改变人口结构、减少牲畜压力、促进民生发展等遏制草地退化态势，提高草地质量。

2000年以来，研究区草地退化态势仍在持续。西部大开发战略的实施，极大带动了区域经济发展，推动了研究区城镇化/工业化进程。伴随着城镇化/工业化发展，大量劳动力涌向城镇，导致了人口结构和数量的改变。截至2014年鄂尔多斯市城镇化速率已经达到72 %，高于全国平均水平23%，鄂尔多斯市城镇化水平达到高速发展水平阶段 [ 18 ] 。尤其2000年以后，鄂尔多斯市煤的开采量呈现幂次方增长 [ 19 ] ，煤矿的开采使得大量人口迁移入城市，与之对应的GDP总量、第二产业总值、城市人口等趋势。工业化和产业结构的改变促进了鄂尔多斯市城市的扩展，城镇化、工业化发展又促进了道路的建设，以及与矿业发展相关的次生产业的发展。鄂尔多斯市典型能源型城市，约70%的鄂尔多斯国土面积下有煤炭资源。大规模的矿业发展导致土地资源直接被破坏，煤炭产业发展过程中固体肥料占用土地资源，地面塌陷破坏土地资源。同时，煤炭发展的一系列衍生产业，例如为了煤炭运输需要加大发展道路建设，这直接破坏道路草地，促使草地退化或者破碎化。煤炭发电产业、煤炭化工产业等污染空气，同时也消耗大量水资源，进而导致内蒙古地区地下水位下降，水体面积减少等一系列土地利用/覆被变化 [ 20 ] [ 21 ] 。因此，在带动地方经发展的同时，更多的生态环境保护措施，例如提高煤炭工业用水效率，发展绿色能源等减少资源消耗的政策应该吸引更加广泛的关注。

表8. 基于关键土地利用管理策略