地球信息科学学报 ›› 2022, Vol. 24 ›› Issue (9): 1835-1852.doi: 10.12082/dqxxkx.2022.220015

• 遥感科学与应用技术 • 上一篇    

基于时序影像的农业活动因子提取与闽西耕地SOC数字制图

聂祥琴1,2(), 陈瀚阅1,2,3,*(), 牛铮3, 张黎明1,2, 刘炜1,2, 邢世和1,2, 范协裕1,2, 李家国4   

  1. 1.福建农林大学资源与环境学院,福州 350002
    2.福建省土壤生态系统健康与调控重点实验室,福州 350002
    3.中国科学院空天信息创新研究院 遥感科学国家重点实验室,北京 100101
    4.中国科学院空天信息创新研究院,北京 100101
  • 收稿日期:2022-01-10 修回日期:2022-04-08 出版日期:2022-09-25 发布日期:2022-11-25
  • 通讯作者: *陈瀚阅(1985—),女,福建莆田人,讲师,主要研究方向为环境遥感。E-mail: Chenhanyue.420@163.com
  • 作者简介:聂祥琴(1994—),女,贵州毕节人,硕士生,研究方向为土壤属性制图。E-mail: niexiangqin2022@163.com
  • 基金资助:
    遥感科学国家重点实验室开放基金项目(OFSLRSS202112);高分辨率对地观测系统重大专项(06-Y30F04-9001-20/22);福建省自然科学基金项目(2021J01117);福建省自然科学基金项目(2019J01660);国家自然科学基金项目(41971050)

Digital SOC Mapping in Croplands Using Agricultural Activity Factors Derived from Time-Series Data in Western Fujian

NIE Xiangqin1,2(), CHEN Hanyue1,2,3,*(), NIU Zheng3, ZHANG Liming1,2, LIU Wei1,2, XING Shihe1,2, FAN Xieyu1,2, LI JiaGuo4   

  1. 1. College of Resources and Environment, Fujian Agriculture and Forestry University, Fuzhou 350002, China
    2. Key Laboratory of Soil Ecosystem Health and Regulation, Fujian Agriculture and Forestry University, Fuzhou 350002, China
    3. State Key Laboratory of Remote Sensing Science, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100101, China
    4. Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100101, China
  • Received:2022-01-10 Revised:2022-04-08 Online:2022-09-25 Published:2022-11-25
  • Contact: CHEN Hanyue
  • Supported by:
    Open Fund of Key Laboratory of Remote Sensing Science(OFSLRSS202112);Major Project of High Resolution Earth Observation(06-Y30F04-9001-20/22);Natural Science Foundation of Fujian Province(2021J01117);Natural Science Foundation of Fujian Province(2019J01660);Natural Science Foundation of China(41971050)

摘要:

人类活动对表层耕地土壤有机碳(Soil Organic Carbon, SOC)影响强烈,但目前大范围复杂地貌地形区的耕地SOC数字制图对人为因素的空间刻画不足。本文以福建省西部耕地为研究对象,基于Sentinel-2/MSI时间序列数据提取轮作模式分类信息(Crop Rotation, CR),以及可反映轮作模式信息的植被特征变换变量(Harmonic Analysis of Time Series, HANTS),分别作为农业活动定性和定量因子,将常规气候和地形因子作为自然环境因子,并对不同类型环境变量进行组合(气候+地形、气候+地形+轮作模式、气候+地形+HANTS变量、气候+地形+轮作模式+HANTS变量)。基于随机森林模型(Random Forest, RF)对不同环境变量组合驱动的耕地表层SOC空间预测精度进行对比分析,探索以轮作模式为例的农业活动因子提高耕地表层SOC数字制图精度的可能性。结果表明,同时加入两种农业活动因子的RF模型表现最佳,其模型预测精度相较于纯自然环境变量驱动的模型有明显提高(R2提高了89.47%,RMSEMAE分别下降了10.66%和12.05%)。轮作模式类型(CR)和HANTS变量两种农业活动因子均被保留参与建模,尤其是轮作模式类型显著影响耕地SOC,在最佳模型的环境变量重要性中排序第四。由此可见,轮作模式相关农业活动因子可有效提高耕地SOC空间预测精度。在所有RF模型中,年降水量(Annual Rainfall, Rainfall)的重要性排名都是第一位。通过最佳模型反演得出该区耕地土壤有机碳均值为18.22±2.99 g/kg,范围为8.25~30.69 g/kg,双季稻和烟稻种植区域SOC含量高于稻菜种植区域。研究结果为复杂地貌地形区耕地土壤有机碳协同变量的更新提供了新的思路。

关键词: 土壤有机碳, HANTS, 轮作模式, 农业活动因子, 空间预测, Sentinel-2, 随机森林, 变量组合

Abstract:

Human activities significantly affect the amount and spatial variation of top Soil Organic Carbon (SOC) in croplands. However, the spatial distribution of agricultural management practices has not been carefully considered in SOC mapping in croplands, especially for croplands in large-scale complex landforms. A case study was conducted in agricultural area in western Fujian Province. Sentinel-2/MSI NDVI time series data were used to derive two types of variables that contains crop rotation information. One is the Crop Rotation (CR) pattern type, which was regarded as qualitative factors of agricultural activities. The other are variables generated using Harmonic Analysis based on sentinel-2 NDVI time series data (HANTS), which were regarded as quantitative factors of agricultural activities. Two types of agricultural activities factors, as well as natural environmental variables were adopted as predictive environmental variables. Four different combinations of above variables according to different categories were formed respectively (i.e., climate factor + terrain factor, climate factor + terrain factor + crop rotation pattern, climate factor + terrain factor + HANTS variables, and climate factor + terrain factor + crop rotation pattern + HANTS variables). Random Forest (RF) models were developed based on four different combinations of above variables for predicting SOC. These RF models were compared to explore whether incorporating agricultural activity factors could improve the SOC mapping accuracy in croplands. Results showed that the combination of natural environment variables with both crop rotation type and variables derived through HANTS yielded the highest accuracy. Compared with the combination of natural environment variables, the prediction accuracy of the optimal model was significantly improved (R2 increased by 89.47%, RMSE and MAE decreased by 10.66% and 12.05%, respectively). Two types of agricultural activity factors were both adopted in optimal model, especially CR significantly affected the SOC in croplands, ranking fourth in the importance of environmental variables of the optimal model. In all RF models, annual rainfall (Rainfall) ranked first in the importance of environmental variables. This indicated that climate factors play a dominant role in soil organic carbon digital soil mapping. The SOC content in croplands of the region predicted from the optimal model was (18.22±2.99) g/kg on average and varied in the range of 8.25~30.69 g/kg. The SOC content in double cropping rice and tobacco-rice planting area were higher than that in rice-vegetable planting area. The results provide a new vision for updating the environmental variables of SOC mapping in complex landform areas.

Key words: soil organic carbon, harmonic analysis of time series, crop rotation mode, factors of agricultural activities, spatial prediction, Sentinel-2, random forest model, combination of variables