基于生态风险评价的采煤矿区土地损毁与复垦过程分析

煤炭开采和复垦活动带来强烈的地表变化和生态环境扰动。该文以山西省平朔矿区为例,通过构建评价模型分析采煤矿区的生态风险动态变化。通过采用最小损毁累积模型测算风险源累积损毁影响值,利用遥感影像结合野外实测数据计算生态脆弱度指数,得出矿区2001年和2010年矿区生态风险值,分析不同开采年份中矿区由于采矿和复垦活动带来的生态风险变化。结果表明:随着煤炭产量的增加,2010年土地损毁的累积影响范围比2001年增长了7 095.17 hm2;但是已复垦排土场的损毁累积危害下降明显,并且随着矿区扰动区域的东移,采矿对研究区西部的影响也在减弱。经过10年的土地复垦与生态重建措施,已复垦的排土场生态风险值呈下降趋势,生态系统趋于稳定。在采矿扰动区中,2001年处于中等风险及以下的区域占比为0.02%,到2010年为16.77%。特大型采煤矿区扰动范围较大,但是采后土地复垦有助于降低局部生态风险的影响。通过研究矿区土地损毁情况、复垦过程和复垦后的状态,分析不同区域生态风险的动态变化,可以为矿区的生态环境治理和区域发展规划提供参考,为矿区管理和相关决策提供科学依据。     

流域"地-矿"土地水资源利用冲突测度确定及土地整治策略

量化区域煤炭开采与水土资源相互影响程度,有助于把握土地利用过程中的主要限制因素,提高土地整治的针对性与有效性。该文以山西省泽州县长河流域为研究区,提出煤炭开采与水土资源利用相互影响测度与土地整治策略的思路和方法。首先提出煤炭开采与水土资源相互影响分析框架,并将土地、水和煤炭作为主要的相互影响对象,将土地利用系统作为影响载体,借助数学模型和GIS空间分析技术,依次评价土地损毁累积作用指数、水资源匮乏指数和煤炭资源损失指数;然后确定综合影响值及其等级空间分布,并据此将流域划分为重点生态修复区、重点生态保护区、损毁土地复垦区、耕地资源保护区、水资源匮乏修复区、水资源与耕地资源保护区和耕地质量提升区等七大土地利用功能区,并根据主要限制因素提出相应的土地整治策略。研究结果表明综合影响高值区及中值区面积占研究区总面积的分别是20.21%和33.39%,主要分布在流域西部,是煤炭开采活动集聚区,煤炭压覆及土地损毁是冲突主要表现特征,划分为重点生态修复区、重点生态保护区、损毁土地复垦区和耕地资源保护区;综合影响低值区及(基本)无影响区分别占24.24%和22.16%,主要分布在流域东部,水资源匮乏是冲突主要表现特征,划分为水资源匮乏修复区、水资源与耕地资源保护区、耕地质量提升区,以上分区与土地利用总体规划中土地用途分区基本吻合。综合影响测度模型有利于定量分析"煤炭-水土"相互利用影响等级与主导限制因素,为制定有效的土地整治策略提供依据。     

基于多源时序NDVI的稀土矿区土地毁损与恢复过程分析

为分析稀土开采的时空分布及矿区土地毁损与恢复过程,该文以岭北稀土矿区为例,以1990—2016年的HJ-1B CCD、Landsat 5和Landsat 8遥感数据为数据源,结合回归分析法、遥感时序NDVI分析方法,对岭北稀土矿区的稀土开采状况及土地毁损与恢复情况进行分析。为减少不同数据由于传感器自身原因而带来的NDVI误差,采用回归分析法构建HJ-1B CCD、Landsat5/8数据的NDVI转换方程并利用均方根误差对转换方程的精度进行检验。结果表明:HJ-1B CCD、Landsat5/8数据的NDVI影像对应转换方程的R2值均超过0.9以及模拟影像与真实影像之间的均方根误差值均小于0.05,说明HJ-1B CCD、Landsat5/8的NDVI数据之间存在较为显著的线性正相关,所求的转换方程具有较高精度,可以消除数据来源差异对NDVI的影响;稀土开采时空分布分析表明岭北稀土矿区在2001—2004年、2006年的开采面积较大,均超过1 km~2,其中2006的开采面积最大,达到2.546 1 km~2;稀土开采在空间分布上较为分散,不仅增加了监管的困难,造成了资源的浪费,还一定程度上增加了治理困难;土地毁损与恢复分析表明受到人为扰动的区域占矿区总面积的54.59%,其中受森林砍伐扰动影响的面积最大,为86.5341 km~2;稀土开采扰动的土地面积为11.354 4 km~2,其植被恢复的平均加权时长为11年。其中复垦恢复的面积为5.004 9 km~2,仍有6.349 5 km~2的区域植被未恢复,急需治理。该研究能较好的反映稀土干扰下的土地毁损与恢复过程,为岭北稀土矿区的生态环境治理及矿区的可持续发展提供借鉴。     

煤矿排土场复垦干扰指数构建及其时空分异

排土场复垦后受干扰情况影响重建生态系统的稳定性,直接关系煤矿区生态修复成败,针对长期恢复过程的干扰状况评价尤其重要。该文在MODIS全球干扰指数的基础上,将地表温度与增强型植被指数相结合,构建基于Landsat影像的排土场复垦干扰指数(dump reclamation disturbance index, DRDI),并计算山西平朔安太堡露天矿排土场28 a间的干扰指数及其时空分异特征。结果表明:复垦初期干扰强烈区呈现较强空间聚集特性,植被结构和功能恢复至自然水平的平均时长约为10a,单个排土场受干扰面积可下降至1.32%～17.76%之间。平朔安太堡露天煤矿4个排土场在稳定时限和干扰分布特征上呈现典型差异,受排矸方式、复垦方向、年限及地形地貌等因素驱动。DRDI可评估复垦排土场重建生态系统稳定时限,有效识别重建生态系统因干扰造成植被退化的空间范围和严重程度,应用于排土场复垦效果评价将为优化复垦模式、提高复垦效率提供参考。干扰指数对自然及人为因素响应的定量分析是未来研究重点。     

中国东北多煤层老矿区采煤沉陷地损毁特征与复垦规划

多煤层老矿区开采时间长,土地损毁程度严重,沉陷地再利用情况复杂。该文通过实地调查和搜集各煤矿地质地形图、采掘图、储量图和采矿历史等资料,借助矿区沉陷预计软件(MSPS),运用概率积分法对鹤岗矿区沉陷地近期(2014-2020年)和远景(2021年-闭矿)2个阶段进行沉陷预测,并采用ArcGIS软件对数据进行统计分析。根据有关规程中各地类沉陷地损毁程度判定标准对鹤岗采煤沉陷地损毁程度进行划分。基于鹤岗采煤沉陷地损毁特征,提出鹤岗矿区沉陷地时空复垦的规划模式。结果表明:随着煤炭不断开采,多煤层老矿区沉陷范围不断扩大,土地损毁程度加重,地表最大下沉深度达39.4 m;沉陷地在进行再利用过程中存在阶段性变化特征,可根据特定阶段沉陷地损毁特征制定土地用途。沉陷地区域稳定性特征和损毁程度是进行沉陷地时空复垦规划的重要依据。从时间安排看,应以闭矿沉陷地损毁程度为着眼点,全面考虑,确定各区域稳沉时序,鹤岗采煤沉陷划分为优先永久治理区、优先临时治理区、动态治理区和未来治理区;从空间布局看,前3个治理区功能定位依次为接续替代产业用地、短期工业用地和农业用地。     

黄土区露天煤矿排土场复垦后土壤与植被的演变规律

恢复受损的土壤和植被是矿区生态恢复的关键,植被恢复过程的实质是植被-土壤复合生态系统相互作用的过程。该文通过典型小区调查的方法,选择山西平朔州安太堡露天煤矿复垦排土场为研究区,分析了不同复垦年限（3、5、10、12和17 a）土壤环境因子和乔木林地植被生物量的动态演变规律,建立了黄土区露天煤矿排土场复垦土壤环境因子和乔木林地植被生物量Logistic演替模型,并构建了土壤-植被交互影响的偏微分方程组。相关系数及显著性检验表明所建立的土壤各环境因子演变模型和乔木林复垦地的植被生物量演变模型有效,能够很好地反映排土场的土壤因子和植被生物量的动态演变过程;随着复垦年限的增加,研究区土壤环境因子质量不断提升并逐渐接近原地貌,土壤因子和植被生物量都呈S型变化,符合Logistic生长演替模型;土壤环境因子与植被生物量二者交互作用明显,符合Kolmogorov捕食模型。该研究可为黄土区露天矿排土场土地复垦与生态恢复提供理论依据。     

煤炭开采对"矿-农"复合区农田质量影响评价

农田质量的内涵应包括本底质量、经济质量和生态环境质量3部分。煤炭开采导致农田严重受损,科学合理评价煤炭开采对矿-农复合区农田质量的影响,对制定区域土地复垦规划和综合整治有着重要的意义。以山西省泽州县长河流域为研究区,该研究首先建立了煤炭开采沉陷区农田质量评价的指标体系,选取农田自然质量、农业生产条件、区位条件、农田环境条件来表征农田质量,对并依据这些指标进行空间可视化,利用多因素加权求和的方法对农田质量进行综合评价。研究结果表明:与煤炭开采非扰动区相比,位于煤矿区井田边界内的农田自然质量、农业生产条件、区位条件、农田生态环境条件评价分值均较低;煤矿井田边界范围内的农田质量(平均值为0.37)低于非矿区农田(平均值为0.55)。应结合煤炭开采扰动影响区农田质量区域分布特点,制定科学合理的矿-农复合区土地复垦规划和土地综合整治方案。     

基于遥感影像的露天煤矿区土地特征信息及分类研究

该研究利用Landsat TM数据,以安太堡大型露天煤矿为例,在对地物光谱特征深入分析的基础上,设计了大型露天矿区土地剧烈扰动下不同地物特征提取模型,提取了安太堡露天矿区植被高覆盖区、植被低覆盖区、剥离堆垫区、采煤运煤区及边坡区等不同地物信息。采用归一化植被指数(NDVI)提取植被低覆盖区与植被高覆盖区信息,采用(TM4-TM5)>0提取植被高覆盖区信息并与NDVI进行了比较,采用TM4<40提取采煤运煤区信息,采用TM4/TM7在0.99～1.01范围来提取边坡区信息, 并统计计算了各类地物所占面积和分布情况。对研究区TM影像进行主成分分析,剥离堆垫区、采煤运煤区和边坡区等反映矿区扰动特征的信息主要由第1主成分反映,植被低覆盖区和高覆盖区等反映矿区植被覆盖特征的信息主要由第2主成分反映,两个主成分的贡献率达到97.16％,并利用扰动特征和植被特征对研究区地物进行了分类。该技术与方法为露天矿地物变化动态监测以及土地复垦与生态重建均提供了准确数据支持。     

黄土丘陵沟壑区煤矿沉陷耕地复垦

探索开发适宜黄土地貌和气候特点以及黄土丘陵采煤沉陷损伤规律的土地复垦模式,对于实现该区数量质量并重管理的耕地保护目标具有重要的现实意义。该文基于开采设计、沉降观测参数和MSAS软件,预测山西兴县斜沟煤矿一采区沉陷损毁土地面积为16.54 km2,损毁土地类型以耕地和草地为主,其中轻度破坏耕地面积3.70 km2,中度破坏耕地面积2.41 km2,重度破坏耕地面积0.17 km2。结合黄土丘陵沟壑低潜水位区土地损毁特征和区域生态保育和国土整治目标,拟对轻度破坏的1.65 km2耕地采用裂缝充填和土地平整复垦措施,按照集中连片的原则对轻度、中度和重度破坏的4.63 km2耕地采用坡改梯的土地复垦工程措施进行复垦,实现增加耕地面积0.40 km2,每年减少水土流失4.94万t,提高了耕地质量等级。研究认为,结合企业复垦资金和社会资金投入,按照集中连片的原则对采煤沉陷破坏的耕地采取坡改梯复垦是该区域损毁耕地复垦的有效模式。     

矿业城市土地损毁生态风险评价与空间防范策略

为量化矿业城市土地损毁生态风险、提高土地复垦和综合整治有效性,以吉林省辽源市为研究区,提出基于GIS的矿业城市土地损毁生态风险定量评价与空间防范的思路和方法。首先将土地损毁作为直接生态风险源,将土地生态系统作为生态风险受体,借助遥感信息模型与GIS空间分析模型,依次评价土地生态系统的生态敏感性、土地损毁累积作用以及二者间的暴露关系;然后确定综合生态风险值及其等级空间分布,并据此将矿业城市土地利用划分为环境敏感特别保护区、一般生态保护区、生态保育区、生态协调区、生态复垦区、限制发展区、工矿发展区与城镇发展区等八大功能区,各区内生态风险防范措施各有侧重;最后利用2010年城市总体规划图对评价与分区结果进行有效性分析,分析结果表明,辽源市90%以上的规划建设用地属于发展主导区域、70%以上的非规划用地属于生态保护主导区,研究结果与2010年城市规划用地布局范围基本吻合。基于 GIS 定量方法进行矿业城市土地损毁生态风险分析可空间化土地利用的风险等级与主导风险因素,为制定有针对性的风险防范措施提供了空间信息;也为矿业城市土地复垦规划、土地整理规划与土地利用总体规划的有机结合提供了参考。     

综合季相节律和特征光谱的冬小麦种植面积遥感估算

及时准确地获取区域和国家尺度的作物种植面积和空间分布具有重要意义。针对目前中低分辨率遥感数据相结合方法的局限,提出一种新的作物类型识别方法。首先基于MODIS NDVI数据的时间优势,提取研究区各类植被的NDVI时间序列曲线,从而分析冬小麦在季相节律上的识别特征,构建冬小麦识别模型。再将MODIS像元分类处理,纯耕地像元利用冬小麦的季相节律特征识别;耕地与其他植被的混合像元利用混合像元分解的思想提取耕地组分的NDVI时间序列,从而进行识别,进一步根据空间关系将识别结果重新定位到中分辨率尺度上;冬小麦与其他作物的混合像元覆盖区则利用TM遥感影像的光谱差异加以区分。在伊洛河流域主要农业区,以冬小麦为识别对象,结果表明识别精度达到96.3%。该方法为作物种植信息的提取提供了新的解决问题的途径,也对其他类型作物的识别也具有重要的参考价值。     

Impact of fertilization on cover crops and microbial community on a bauxite‐mined soil undergoing reclamation

The reclamation of bauxite‐mined areas can be favored by the application of organic and/or chemical fertilization to restore the vegetation. Otherwise, the impact of fertilizations on soil microbiota or plant–microbe interactions as land reclamation progresses is less understood. To address this issue, we evaluated the impact of organic and chemical fertilization on plants and soil microbial community within the first 36 months of land reclamation in a bauxite‐mined site. The experiment was set up according to a split‐plot design in which the main plots received fertilizations [non‐fertilized control (NF), chemical fertilization (CF; NPK and rock phosphate), organic fertilization (OF; poultry litter), and CF+OF combined], and the subplots received cover crops [no cover crops (NC), grass (B; Brachiaria), legume (S, Stylosanthes), and B+S combined]. Cover crops biomass yield was assessed annually with five field campaigns per year. We used phospholipid fatty acids (PLFAs) to infer the impacts of mining and restoration practices on actinobacteria, Gram‐negative and Gram‐positive bacteria, arbuscular mycorrhizal fungi, and fungi. Accordingly, PLFAs were determined before bauxite mining (pre‐mining), six months after topsoil reconfiguration (post‐mining), and after 14 and 36 months following the application of the fertilizations and cover crops. PLFAs results indicated that in post‐mining, the living microbiota was significantly lower than in pre‐mining. Cover crops biomass yield was highest for B and B+S fertilized with CF+OF at 14 and 36 months. Both parametric and non‐parametric statistics showed a temporal variation in the response of living microbes to the treatments applied. After 14 months, living microbes showed greatest response to OF, while at 36 months their response was strongest in the treatments with highest plant biomass production (B and B+S). These results suggest that in the early stages of land reclamation, living microbial biomass benefit the most from organic fertilizers. As this initial boost decline, living microbes are more likely to thrive in areas undergoing reclamation where they can develop synergistic interactions with plants.     

融合光谱混合分解与面向对象的土地利用/覆被分类

错综复杂的土地利用模式和破碎的地物斑块制约了土地利用/覆被分类的精度和效率。一方面,混合像元模糊了地物的光谱信息,影响了分类精度。另一方面,如何高效利用地物的光谱、形状和纹理特征是当前土地利用/覆被分类的研究热点。为了提高基于遥感技术的土地利用/覆被分类精度,该研究基于Sentinel-2A遥感影像,开展融合光谱混合分解与面向对象的土地利用/覆被分类研究。首先,基于地物的光谱、形状和纹理特征,在3个分割尺度通过NDWI（Normalized Difference Water Index)、NDVI（Normalized Difference Vegetation Index）、SBL（Soil Background Level）等8个特征参数构建了不同地物信息的提取规则。其次,利用光谱混合分解模型提取研究区基质（SL;岩石和土壤）、植被（GV;光合作用叶片）和暗色物质（DA;阴影和水）3类通用端元。最后,尝试融合3端元光谱特征优化地物信息提取规则。研究结果表明：1）基于构建的光谱、形状和纹理的地物信息提取规则,使用模糊函数、阈值法进行土地利用/覆被分类,获得了较高的分类精度,总体精度为80.83%,Kappa系数为0.76。2）融合3端元的光谱特征的提取规则将分类精度提升至90.00%,Kappa系数提升至0.88。3）具有明确物理意义的3端元的融入增强了像元内各组分信息的差异性,弥补了传统光谱指数对植被与土壤间的亮度信息解析度不足的缺陷。该方法能充分利用影像的光谱信息,是一种由易到难、对不确定因素进行逐层剥离的土地利用/覆被信息提取技术。因此,对中高分辨率的多光谱遥感影像十分友好,在土地利用/覆被的精细化分类中有较大应用潜力。     

基于GIS与RS的大理河流域植被格局分形维数时空变化特征

为掌握大理河流域植被格局时间和空间分布特征,构建了基于分形布朗运动理论的流域植被格局量化模型。该文利用自1990年至2006年共计5期TM/ETM影像为基础信息源,在ARCGIS平台下计算得到大理河流域各时相的像元归一化植被指数(normalized difference vegetation index,NDVI),建立了流域像元尺度的栅格结构数字植被模型(digital vegetation model,DVM),并逐个计算得到各子流域植被格局分形维数。结果显示分形维数都大于2.5,介于2.7311~2.8499之间,各下一级子流域植被覆盖分形维数的大小都没有超过其所在的更大流域的分形维数。植被格局分形维数(fractional brownian motion,FBM)从下游到上游随着离出口点距离的增大而逐渐变小,植被格局自下游到上游逐渐趋于破碎复杂。各子流域从1990年至2006年植被格局分形维数基本随着时间的变化呈现了先减小后增大的趋势。该文建立的流域植被格局分形量化模型在综合量化流域植被格局破碎程度等方面具有一定的优越性。     

基于多源遥感数据融合和LSTM算法的作物分类研究

准确、及时地获取农作物的空间分布信息,对于指导农业生产、制定农业政策具有重要意义。为了检验长短时记忆网络(long short-term memory,LSTM)算法在基于时序遥感数据进行作物分类中的优势,该文以临汾盆地为研究区域,利用Savitzky-Golay滤波对MODIS NDVI进行平滑处理,并采用ESTARFM(enhanced spatial and temporal adaptive reflectance fusion model)算法对滤波后的MODIS NDVI和Landsat NDVI进行融合,生成空间分辨率为30 m、时间分辨率为8天的时序NDVI。基于Landsat NDVI利用LSTM算法进行作物分类,同时,基于融合NDVI分别利用LSTM算法和神经网络(neuralnetwork,NN)算法进行作物分类,并对比3种方法的分类精度。结果表明,Savitzky-Golay滤波后的时序MODISNDVI能够反映不同作物的物候特征;基于融合NDVI的分类精度明显高于基于LandsatNDVI的分类精度,表明融合后的时序NDVI由于具有更高的时间分辨率,能够更加突出不同作物的物候特征,显著提高作物分类精度;基于融合NDVI和LSTM算法的分类精度高于基于融合NDVI和NN算法的分类精度,前者的冬小麦面积估测精度高于后者的估测精度,表明LSTM算法的分类精度高于NN算法。该文可为基于遥感影像进行不同作物种植区域提取的研究提供重要的方法参考。     

Application of spectral mixture modelling to the regional assessment of land degradation: A case study from basilicata,Italy

Mapping and monitoring of land degradation processes such as soil erosion has become an important task for both agricultural and environmental planners. The potential of using SPOT-HRV data for assessing the types of soil erosion and land degradation is obvious. The spectral information and the spatial resolution of the multispectral data allow a high accuracy in local mapping and rapid regional assessment. In Basilicata, southern Italy, spectral mixture modelling (end-member techniques) has been applied to extract relevant information for assessing soil erosion. This method allows an estimate of the proportion of surface cover types (end-members) within the pixel. From the data sets average spectral responses were extracted for growing vegetation, non-green sclerophylous vegetation mixed with dry grasses and bare soil. Using these spectra as end-members in the model, it was possible to decompose the spectral information of all pixels into the three surface components, giving the percentage cover of the three types within every pixel. Classification of the pixels according to the percentage of surface cover or bare soil allows an assessment of erosion or erosion risk.     

Bird utilisation of rehabilitated bauxite minesites in Western Australia

Alcoa of Australia Limited rehabilitate bauxite mined areas in the southwest of Western Australia. One of the aims of the revegetation programme is to promote the return of fauna species which inhabited areas prior to mining.This paper discusses the extent to which bird species utilise revegetated bauxite minesites. It was shown that revegetated areas as young as 4–5 years can support similar bird species numbers, densities and diversities as unmined forest. However, the avifaunal populations occupying any rehabilitated area depend largely on the rehabilitation technique used. Factors which may be used to increase further avifaunal utilisation of minesites are discussed.The effects of jarrah dieback and fire on forest birds are compared with results for bauxite mined areas. From this, it is possible to assess the combined impact of dieback, fire and mining on the longer term conservation of jarrah forest bird species.     

基于归一化植被指数的流域植被覆盖分形维数研究

量化和表征植被覆盖空间分布的复杂性是研究植被与地表物质迁移关系的重要问题。为研究植被覆盖分形维数的表征特点及其在不同尺度上的变化特征,该文基于地理信息系统平台,利用分形布朗运动理论,结合像元归一化植被指数(normalized difference vegetation index,NDVI)的空间分布,提出并计算了植被覆盖分形布朗运动(fractional brownian motion,FBM)分形维数。结果显示,流域植被覆盖的空间分布具有统计自相似性,可用分形维数表征,其值在2.5~3之间,越接近2.5,表示植被覆盖空间分布越复杂。植被覆盖FBM分形维数与流域内均值化NDVI值和NDVI值的变异系数无直接关系,与单位面积不同NDVI值的像元数呈极显著负相关(R=0.66,P0.01)。植被覆盖FBM分形维数具有尺度效应,随流域面积增大而增大,到一定尺度后趋于平稳。流域植被覆盖FBM分形维数能够克服NDVI像元点奇异值等对植被空间分布量化表征计算的干扰,相对传统表征植被覆盖的指数,其在水文、土壤侵蚀等模型中具有更广泛的应用意义。     

潞安矿区净初级生产力和土地覆被变化及成因

净初级生产力和土地覆被变化是指示区域生态环境变化的重要指标.基于EOS/MODIS卫星遥感数据,运用地理信息系统软件ArcGIS 9.3的统计与分析模块对2001-2006年潞安矿区净初级生产力、土地覆被变化及其成因进行了研究和分析.结果表明:(1) 2001-2003年的矿区绝大部分区域净初级生产力大于0.2 kg/(m2·a),而2004-2006年的净初级生产力低于0.2 kg/(m2·a)即研究时段内矿区净初级生产力呈明显减少趋势;(2)灌丛向稀树草地,稀树草地向农用地,地表水域向裸地转换是土地覆被变化的主要特征;(3)煤炭开采,农用地开垦等人类活动和降水量减少,气温升高为特征的气候变化是研究区植被覆盖度和净初级生产力降低的两个主要原因.总之,煤炭开采和土地垦殖为主的人类活动所引起的土地覆被变化直接导致了潞安矿区净初级生产力的降低.     

基于地块级时序遥感的喀斯特石漠化地区撂荒地时空演变

撂荒地在全球呈扩散趋势,传统的调查方法已无法满足及时、精确掌握撂荒情况的需求。针对喀斯特石漠化耕地的撂荒因光学影像的不足和严重的混合像元不易精确识别和提取,以贵州省息烽县为例,运用高精度影像精确识别提取耕地地块,再叠加通过MVC方法计算Landsat数据得到的2003-2018年时序NDVI数据,得到撂荒和复垦时空分布。结果表明:(1)地块和时序NDVI组合能够精准地识别提取撂荒地,该方法对于光学数据缺乏的和耕地破碎的多云雨山区有很好的适用效果;(2)2003-2018年息烽县撂荒地呈随机散点状分布,主要分布在有岩石裸露潜在、轻度、中度石漠化地区,2008年撂荒面积最大2 545.6 hm^2,撂荒率最高7.7%;有7块耕地持续撂荒时间最长12年,有8块地累计撂荒最多15年。(3)2003-2018年息烽县复垦率与撂荒率波动趋势基本一致,复垦率的响应延后1年,较高的石漠化等级负面影响复垦行为。研究结果为喀斯特多云雨山区撂荒地的精准识别提取提供高效、可行的思路和方法,为撂荒地的驱动因子分析、趋势监测和风险评价、效应评估及政策设计提供精确的基础数据。