留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

海南岛近20年GPP变化格局及驱动因素分析

雷济舟 崔嵬 朱济帅 张润卿 赵俊福 章杰 张翔 孙仲益

雷济舟,崔嵬,朱济帅,等. 海南岛近20年GPP变化格局及驱动因素分析[J]. 热带生物学报,2024, 15(1):42−51. doi:10.15886/j.cnki.rdswxb.20230006 doi:  10.15886/j.cnki.rdswxb.20230006
引用本文: 雷济舟,崔嵬,朱济帅,等. 海南岛近20年GPP变化格局及驱动因素分析[J]. 热带生物学报,2024, 15(1):42−51. doi:10.15886/j.cnki.rdswxb.20230006 doi:  10.15886/j.cnki.rdswxb.20230006
LEI Jizhou, CUI Wei, ZHU Jishuai, ZHANG Runqing, ZHAO Junfu, ZHANG Jie, ZHANG Xiang, SUN Zhongyi. Analysis of varying patterns and driving factors of GPP in Hainan Island during the past 20 years[J]. Journal of Tropical Biology, 2024, 15(1): 42-51. doi: 10.15886/j.cnki.rdswxb.20230006
Citation: LEI Jizhou, CUI Wei, ZHU Jishuai, ZHANG Runqing, ZHAO Junfu, ZHANG Jie, ZHANG Xiang, SUN Zhongyi. Analysis of varying patterns and driving factors of GPP in Hainan Island during the past 20 years[J]. Journal of Tropical Biology, 2024, 15(1): 42-51. doi: 10.15886/j.cnki.rdswxb.20230006

海南岛近20年GPP变化格局及驱动因素分析

doi: 10.15886/j.cnki.rdswxb.20230006
基金项目: 海南大学科研启动项目(KYQD(ZR)20069);海南省自然科学基金青年基金项目(320QN202);国家重点研发项目(2021YFD2200404)
详细信息
    第一作者:

    雷济舟(1996−),男,海南大学生态与环境学院2020级生态学硕士研究生。E-mail:hnhn604@163.com

    通信作者:

    孙仲益(1989−),男,副教授,博导。研究方向:生态信息学。E-mail: gis.rs@hainanu.edu.cn

  • 中图分类号: X 82.826

Analysis of varying patterns and driving factors of GPP in Hainan Island during the past 20 years

  • 摘要: 为探究气象因素与人类活动对近20年海南岛总初级生产力(Gross Primary Production,GPP)变化的相对贡献,首先利用Theil-Sen及Mann-Kendall方法获取海南岛GPP整体时空分布特征,在此基础上,分别以土地利用及覆被变化(Land Use and Cover Change,LUCC)指示人类活动,以气温(Air Temperature,Ta)、饱和水汽压差(Vapor Pressure Deficit,VPD)以及光合有效辐射( Photosynthetically Active Radiation,PAR)作为气象要素指标,通过空间统计与机器学习手段,构建海南岛GPP变化格局归因模型,量化驱动因素的相对贡献。结果表明:研究期间内,时间上,海南岛GPP呈现0.44 Tg ·a−1的显著增加趋势(P=0.024);空间上,海南岛87.8%面积的GPP表现增加趋势,海口及三亚周边等小部分区域(约9%)则表现为下降;海南岛土地利用共计转移15 528.40 km2,主要发生于林地地区,林地净增加642.88 km2,草地地区转移面积为4 759.28 km2,耕地地区转移的面积为4 051.23 km2;相较于人类活动,研究期年际间气象条件的差异是海南岛GPP变化的主导因素,但特殊年份中LUCC对海南省市县的影响则有所凸显。
  • 图  2  海南岛GPP时空间变化趋势

    图  3  各土地利用类型在LUCC所驱动GPP变化

    图  4  相对贡献率变化图

    图  5  土地分类对相对贡献率的影响及各市(县)相对贡献率变化

    图  6  海南岛各气象要素变化图

    表  1  根据IGBP分类系统对土地利用类型的重分类

    土地利用
    类型编号
    重分类后土地
    利用类型
    IGBP分类
    1 林地 常绿阔叶林
    常绿针叶林
    落叶阔叶林
    混交林
    有林草地
    稀树草地
    2 灌丛 稀疏灌木丛及
    封闭灌木丛
    3 草地 草地
    4 耕地 耕地、
    耕地或自然植被
    5 城市 城市以及建筑用地
    6 裸地 裸地
    7 湿地 永久湿地及水域
    下载: 导出CSV

    表  2  模型评价精度表

    土地利用
    类型
    SE/
    (g·m−2·mon−1
    R2模型
    算法
    用地类型
    未变栅格
    数量/个
    林地9.600.60RF97534
    草地6.350.66RF2199
    耕地4.880.81RF7762
    城市5.410.62RF1568
    裸地3.390.29GPR33
    湿地10.380.12GPR1583
    下载: 导出CSV

    表  3  海南岛近20年间土地利用每5年的转移情况

    转移路径2001—2005
    面积/km2
    2005—2010
    面积/km2
    2010—2015
    面积/km2
    2015—2019
    面积/km2
    林地→草地517.43327.40325.38658.50
    林地→耕地695.82826.56794.801684.88
    林地→城市19.8549.4939.1740.23
    林地→裸地1.060.260.790.79
    林地→湿地10.0633.0827.5354.52
    灌丛→草地2.121.590.790.53
    草地→林地574.331237.59810.42873.14
    草地→耕地351.48215.97246.67358.63
    草地→城市2.127.414.233.44
    草地→裸地0.791.851.061.32
    草地→湿地3.9725.6711.6526.20
    耕地→林地329.25933.22742.131143.64
    耕地→草地181.30476.14328.19361.54
    耕地→城市1.593.440.530.26
    裸地→林地0.532.381.590.79
    裸地→草地0.531.592.651.85
    裸地→湿地0.792.121.321.59
    湿地→林地30.1713.2316.4140.49
    湿地→草地11.123.713.7112.17
    湿地→裸地1.320.530.531.06
    下载: 导出CSV
  • [1] CAMPBELL J E, BERRY J A, SEIBT U, et al. Large historical growth in global terrestrial gross primary production[J]. Nature, 2017, 544(7648): 84 − 87. doi:  10.1038/nature22030
    [2] SUN Z, WANG X, ZHANG X, et al. Evaluating and comparing remote sensing terrestrial GPP models for their response to climate variability and CO2 trends[J]. Science of the Total Environment, 2019, 668: 696 − 713. doi:  10.1016/j.scitotenv.2019.03.025
    [3] BO Y, LI X, LIU K, et al. Three decades of gross primary production (GPP) in China: variations, trends, attributions, and prediction inferred from multiple datasets and time series modeling[J]. Remote Sensing, 2022, 14(11): 2564. doi:  10.3390/rs14112564
    [4] 丁仲礼, 段晓男, 葛全胜, 等. 国际温室气体减排方案评估及中国长期排放权讨论[J]. 中国科学(D辑: 地球科学), 2009, 39(12): 1659 − 1671.
    [5] JANSSENS I A, FREIBAUER A, CIAIS P, et al. Europe's terrestrial biosphere absorbs 7 to 12% of European anthropogenic CO2 emissions[J]. Science, 2003, 300(5625): 1538 − 1542. doi:  10.1126/science.1083592
    [6] SUN Z, WANG X, YAMAMOTO H, et al. Spatial pattern of GPP variations in terrestrial ecosystems and its drivers: Climatic factors, CO2 concentration and land-cover change, 1982–2015[J]. Ecological informatics, 2018, 46: 156 − 165. doi:  10.1016/j.ecoinf.2018.06.006
    [7] KNAPP A K, CIAIS P, SMITH M D. Reconciling inconsistencies in precipitation-productivity relationships: implications for climate change[J]. New Phytologist, 2017, 214(1): 41 − 47. doi:  10.1111/nph.14381
    [8] HU L, FAN W, LIU S, et al. Temporal and spatial distribution and variation of GPP in MOHE, China[C]. Fort Worth, TX, USA: IEEE, 2017.
    [9] MCCARTHY J K, DWYER J M, MOKANY K. Direct climate effects are more influential than functional composition in determining future gross primary productivity[J]. Landscape Ecology, 2020, 35(4): 969 − 984. doi:  10.1007/s10980-020-00994-x
    [10] QIU R, HAN G, MA X, et al. CO2 concentration, A critical factor influencing the relationship between solar-induced chlorophyll fluorescence and gross primary productivity[J]. Remote Sensing, 2020, 12(9): 1377. doi:  10.3390/rs12091377
    [11] WU C, NIU Z, GAO S. Gross primary production estimation from MODIS data with vegetation index and photosynthetically active radiation in maize[J]. Journal of Geophysical Research: Atmospheres, 2010, 115: D12.
    [12] SUN Z, WANG X, YAMAMOTO H, et al. The effects of spatiotemporal patterns of atmospheric CO2 concentration on terrestrial gross primary productivity estimation[J]. Climatic Change, 2020, 163(2): 913 − 930. doi:  10.1007/s10584-020-02903-2
    [13] HUANG X, XIAO J, WANG X, et al. Improving the global MODIS GPP model by optimizing parameters with FLUXNET data[J]. Agricultural and Forest Meteorology, 2021, 300: 108314. doi:  10.1016/j.agrformet.2020.108314
    [14] VERMA M, FRIEDL M A, LAW B E, et al. Improving the performance of remote sensing models for capturing intra-and inter-annual variations in daily GPP: An analysis using global FLUXNET tower data[J]. Agricultural and Forest Meteorology, 2015, 214: 416 − 429.
    [15] JOINER J, YOSHIDA Y, ZHANG Y, et al. Estimation of terrestrial global gross primary production (GPP) with satellite data-driven models and eddy covariance flux data[J]. Remote Sensing, 2018, 10(9): 1346. doi:  10.3390/rs10091346
    [16] CHEN Y, FENG X, TIAN H, et al. Accelerated increase in vegetation carbon sequestration in China after 2010: A turning point resulting from climate and human interaction[J]. Global Change Biology, 2021, 27(22): 5848 − 5864. doi:  10.1111/gcb.15854
    [17] DING Z, ZHENG H, LIU Y, et al. Spatiotemporal patterns of ecosystem restoration activities and their effects on changes in terrestrial gross primary production in Southwest China[J]. Remote Sensing, 2021, 13(6): 1209. doi:  10.3390/rs13061209
    [18] REZENDE LUIZ F C, ANDERSON C A, CELSO V R, et al. Impacts of land use change and atmospheric CO2 on gross primary productivity (GPP), evaporation, and climate in southern Amazon[J]. Journal of Geophysical Research: Atmospheres, 2022, 127(8): e2021JD034608.
    [19] MA J, XIAO X, MIAO R, et al. Trends and controls of terrestrial gross primary productivity of China during 2000–2016[J]. Environmental Research Letters, 2019, 14(8): 084032. doi:  10.1088/1748-9326/ab31e4
    [20] ZHANG Y, ZHANG C, WANG Z, et al. Vegetation dynamics and its driving forces from climate change and human activities in the Three-River Source Region, China from 1982 to 2012[J]. Science of the Total Environment, 2016, 563/564: 210 − 220.
    [21] YOU N, MENG J, ZHU L, et al. Isolating the impacts of land use/cover change and climate change on the GPP in the Heihe River Basin of China[J]. Journal of Geophysical Research: Biogeosciences, 2020, 125(10): e2020JG005734.
    [22] ICHII K, HASHIMOTO H, NEMANI R, et al. Modeling the interannual variability and trends in gross and net primary productivity of tropical forests from 1982 to 1999[J]. Global and Planetary Change, 2005, 48(4): 274 − 286. doi:  10.1016/j.gloplacha.2005.02.005
    [23] ZHANG Y, XIAO X, GUANTER L, et al. Precipitation and carbon-water coupling jointly control the interannual variability of global land gross primary production[J]. Scientific reports, 2016, 6: 39748. doi:  10.1038/s41598-016-0001-8
    [24] PAU S, DETTO M, KIM Y, et al. Tropical forest temperature thresholds for gross primary productivity[J]. Ecosphere, 2018, 9(7): e02311.
    [25] ZENG Z, ESTES L, ZIEGLER A D, et al. Highland cropland expansion and forest loss in Southeast Asia in the twenty-first century[J]. Nature Geoscience, 2018, 11(8): 556 − 562. doi:  10.1038/s41561-018-0166-9
    [26] QUESADA B, ARNETH A, ROBERTSON E, et al. Potential strong contribution of future anthropogenic land-use and land-cover change to the terrestrial carbon cycle[J]. Environmental Research Letters, 2018, 13(6): 064023. doi:  10.1088/1748-9326/aac4c3
    [27] 覃艳. 自贸港建设背景下海南金融与外贸的关系研究[J]. 中国商论, 2021(19): 4 − 6.
    [28] 熊坚, 冯婉吉, 朱罗娜. 海南自贸港建设的理论探索与实践创新——基于“双循环”新发展格局视角[J]. 管理现代化, 2021, 41(5): 61 − 67. doi:  10.19634/j.cnki.11-1403/c.2021.05.015
    [29] SULLA-MENASHE D, FRIEDL M A. User guide to collection 6 MODIS land cover (MCD12Q1 and MCD12C1) product[J]. Usgs:Reston, Va, Usa, 2018, 1: 18.
    [30] ZHANG Y, WANG Q, WANG Z, et al. Impact of human activities and climate change on the grassland dynamics under different regime policies in the Mongolian Plateau[J]. Science of the Total Environment, 2020, 698: 134304. doi:  10.1016/j.scitotenv.2019.134304
    [31] RUNNING S W, ZHAO M. Daily GPP and annual NPP (MOD17A2/A3) products NASA Earth Observing System MODIS land algorithm[J]. MOD17 User′s Guide, 2015: 1-28.
    [32] LI Y, ZHANG Y, LV J. Interannual variations in GPP in forest ecosystems in Southwest China and regional differences in the climatic contributions[J]. Ecological Informatics, 2022, 69: 101591. doi:  10.1016/j.ecoinf.2022.101591
    [33] SUN W, MU X, SONG X, et al. Changes in extreme temperature and precipitation events in the Loess Plateau (China) during 1960–2013 under global warming[J]. Atmospheric Research, 2016, 168: 33 − 48. doi:  10.1016/j.atmosres.2015.09.001
    [34] SUN W, SHAO Q, LIU J, et al. Assessing the effects of land use and topography on soil erosion on the Loess Plateau in China[J]. CATENA, 2014, 121: 151 − 163. doi:  10.1016/j.catena.2014.05.009
    [35] ALDUCHOV O A, ESKRIDGE R E. Improved magnus form approximation of saturation vapor pressure[J]. Journal of Applied Meteorology, 1996, 35(4): 601 − 609. doi:  10.1175/1520-0450(1996)035<0601:IMFAOS>2.0.CO;2
    [36] LIANG S, ZHAO X, LIU S, et al. A long-term global land surface satellite (GLASS) data-set for environmental studies[J]. 国际数字地球学报(英文), 2013, 6(S1): 5 − 33. doi:  10.1080/17538947.2013.805262

    LIANG S, ZHAO X, LIU S, et al. A long-term global land surface satellite (GLASS) data-set for environmental studies[J]. International Journal of Digital Earth, 2013, 6(sup1): 5-33. doi:  10.1080/17538947.2013.805262
    [37] DING Z, ZHENG H, LI H, et al. Afforestation-driven increases in terrestrial gross primary productivity are partly offset by urban expansion in Southwest China[J]. Ecological Indicators, 2021, 127: 107641. doi:  10.1016/j.ecolind.2021.107641
    [38] ZHENG J, WU W, YU L, et al. Coconut trees detection on the tenarunga using high-resolution satellite images and deep learning[C]//2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS. IEEE, 2021: 6512-6515.
    [39] YI Z F, CANNON C H, CHEN J, et al. Developing indicators of economic value and biodiversity loss for rubber plantations in Xishuangbanna, southwest China: a case study from Menglun township[J]. Ecological Indicators, 2014, 36: 788 − 797. doi:  10.1016/j.ecolind.2013.03.016
    [40] 王雪芳. 海南省霸王岭国有林场生态保护对策研究[D]. 长沙: 中南林业科技大学, 2016
    [41] 胡会峰, 刘国华. 中国天然林保护工程的固碳能力估算[J]. 生态学报, 2006, 26(1): 291 − 296. doi:  10.3321/j.issn:1000-0933.2006.01.035
    [42] 许涵, 李意德, 骆土寿, 等. 达维台风对海南尖峰岭热带山地雨林群落的影响[J]. 植物生态学报, 2008, 32(6): 1323. doi:  10.3773/j.issn.1005-264x.2008.06.013
    [43] DELGADO R C, PEREIRA M G, TEODORO P E, et al. Seasonality of gross primary production in the Atlantic Forest of Brazil[J]. Global Ecology and Conservation, 2018, 14: e00392. doi:  10.1016/j.gecco.2018.e00392
    [44] CUI W, XIONG Q, ZHENG Y, et al. A study on the vulnerability of the gross primary production of rubber plantations to regional short-term flash drought over Hainan Island[J]. Forests, 2022, 13(6): 893. doi:  10.3390/f13060893
    [45] NIU Z, YAN H, LIU F. Decreasing cropping intensity dominated the negative trend of cropland productivity in southern China in 2000–2015[J]. Sustainability, 2020, 12(23): 10070. doi:  10.3390/su122310070
    [46] YAN H, LIU F, QIN Y, et al. Tracking the spatio-temporal change of cropping intensity in China during 2000–2015[J]. Environmental Research Letters, 2019, 14(3): 035008. doi:  10.1088/1748-9326/aaf9c7
    [47] GUO Z, ZHANG Y, DEEGEN P, et al. Economic analyses of rubber and tea plantations and rubber-tea intercropping in Hainan, China[J]. Agroforestry systems, 2006, 66(2): 117 − 127. doi:  10.1007/s10457-005-4676-2
    [48] YU B, CHAO X, ZHANG J, et al. Effectiveness of nature reserves for natural forests protection in tropical Hainan: a 20 year analysis[J]. Chinese Geographical Science, 2016, 26(2): 208 − 215. doi:  10.1007/s11769-016-0800-7
    [49] OLSSON P O, HELIASZ M, JIN H, et al. Mapping the reduction in gross primary productivity in subarctic birch forests due to insect outbreaks[J]. Biogeosciences, 2017, 14(6): 1703 − 1719. doi:  10.5194/bg-14-1703-2017
    [50] ZHAO M, HEINSCH F A, NEMANI R R, et al. Improvements of the MODIS terrestrial gross and net primary production global data set[J]. Remote sensing of Environment, 2005, 95(2): 164 − 176. doi:  10.1016/j.rse.2004.12.011
    [51] SULLA-MENASHE D, GRAY J M, ABERCROMBIE S P, et al. Hierarchical mapping of annual global land cover 2001 to present: The MODIS Collection 6 Land Cover product[J]. Remote Sensing of Environment, 2019, 222: 183 − 194. doi:  10.1016/j.rse.2018.12.013
    [52] LI X, LIANG S, YU G, et al. Estimation of gross primary production over the terrestrial ecosystems in China[J]. Ecological Modelling, 2013, 261/262: 80 − 92.
    [53] LI G, ZHANG F, JING Y, et al. Response of evapotranspiration to changes in land use and land cover and climate in China during 2001–2013[J]. Science of the Total Environment, 2017, 596/597: 256 − 265.
  • [1] 刘嘉慧, 余瑞, 王有, 李欣, 陈荣昊.  基于多光谱无人机的不同放牧策略对草地叶面积指数变化动态解析 . 热带生物学报, 2024, 15(1): 60-72. doi: 10.15886/j.cnki.rdswxb.20230047
    [2] 谢建引, 曾安, 张洪亮.  中国农作物生产及种业现状与分析 . 热带生物学报, 2023, 14(3): 259-267. doi: 10.15886/j.cnki.rdswxb.2023.03.003
    [3] 钟王月, 邹伟, 徐诗琴, 林积泉.  海南岛农用地土壤重金属的生态风险评价 . 热带生物学报, 2023, 14(6): 668-674. doi: 10.15886/j.cnki.rdswxb.20220104
    [4] 郭沈沈, 吴志祥, 刘文杰, 赵豪, 李运帅, 王誉博, 王鹏.  热带北缘地区橡胶林水分利用效率的季节差异及归因分析 . 热带生物学报, 2023, 14(4): 424-432. doi: 10.15886/j.cnki.rdswxb.2023.04.010
    [5] 史娴, 聂堂哲, 熊千, 刘兆鑫, 张嘉怡, 刘文杰, 乌兰, 崔嵬, 孙仲益.  基于InVEST与MaxEnt模型的海南岛红树林生态系统碳储量增量预估 . 热带生物学报, 2023, 14(3): 298-306. doi: 10.15886/j.cnki.rdswxb.2023.03.007
    [6] 柯宏基, 陈傅晓, 李向民, 樊佳伟, 王永波, 符书源, 谭围.  基于线粒体控制区的海南岛3种弹涂鱼的遗传多样性 . 热带生物学报, 2023, 14(2): 203-213. doi: 10.15886/j.cnki.rdswxb.2023.02.010
    [7] 安映荷, 张润卿, 刘文杰, 王鹏, 章杰, 乌兰, 孙仲益.  海南岛橡胶林区域不同SIF产品的差异性分析及其对GPP估算的影响 . 热带生物学报, 2023, 14(4): 412-423. doi: 10.15886/j.cnki.rdswxb.2023.04.009
    [8] 邹海平, 张京红, 李伟光, 陈小敏, 白蕤, 吕润.  海南岛水稻需水量与缺水量的时空变化特征 . 热带生物学报, 2023, 14(5): 569-576. doi: 10.15886/j.cnki.rdswxb.20230074
    [9] 陈红, 郭冬艳, 吴俞, 李勋.  两类季风槽特征及海南岛的强降水机制 . 热带生物学报, 2022, 13(4): 348-357. doi: 10.15886/j.cnki.rdswxb.2022.04.005
    [10] 刘应帅, 余瑞, 郑彬彬, 刘嘉慧, 宋奇, 陈荣昊, 严哲.  海南岛森林植被NEP季节性时空变化规律及气候驱动因素分析 . 热带生物学报, 2022, 13(2): 166-176. doi: 10.15886/j.cnki.rdswxb.2022.02.008
    [11] 杨静, 张亚杰, 张京红, 张明洁.  基于MODIS NDVI的海南岛植被覆盖面积的反演技术 . 热带生物学报, 2022, 13(4): 397-403. doi: 10.15886/j.cnki.rdswxb.2022.04.011
    [12] 吴佳妮, 樊晶, 李勋, 邢彩盈.  影响海南岛的热带气旋强度和路径的主要特征(1949—2020年) . 热带生物学报, 2022, 13(4): 339-347. doi: 10.15886/j.cnki.rdswxb.2022.04.004
    [13] 符钉辉, 丁伟品, 龙军桥, 黄赞慧, 张小海, 陈泽恒.  海南岛湿地鸟类的资源调查及多样性研究 . 热带生物学报, 2022, 13(2): 120-126. doi: 10.15886/j.cnki.rdswxb.2022.02.003
    [14] 李光伟, 邢峰华, 敖杰, 毛志远.  海南岛30年空中云水资源时空特征分析 . 热带生物学报, 2022, 13(4): 331-338. doi: 10.15886/j.cnki.rdswxb.2022.04.003
    [15] 符传博, 佟金鹤, 徐文帅, 刘丽君.  海南岛臭氧污染时空特征及其成因分析 . 热带生物学报, 2022, 13(4): 404-409. doi: 10.15886/j.cnki.rdswxb.2022.04.012
    [16] 吴胜安, 邢彩盈, 朱晶晶.  海南岛气候特征分析 . 热带生物学报, 2022, 13(4): 315-323. doi: 10.15886/j.cnki.rdswxb.2022.04.001
    [17] 李伟光, 张京红, 刘少军, 车秀芬, 陈小敏, 邹海平.  海南岛干旱的气象特征及监测指标 . 热带生物学报, 2022, 13(4): 324-330. doi: 10.15886/j.cnki.rdswxb.2022.04.002
    [18] 陈小敏, 李伟光, 梁彩红, 白蕤, 吴慧.  海南岛主要农业气象灾害特征及防御措施分析 . 热带生物学报, 2022, 13(4): 416-421. doi: 10.15886/j.cnki.rdswxb.2022.04.014
    [19] 侯祥文, 徐诗涛, 王德立, 梁靖雯, 司更花.  海南岛鸟巢蕨自然种群遗传多样性SRAP分析 . 热带生物学报, 2021, 12(1): 25-32. doi: 10.15886/j.cnki.rdswxb.2021.01.004
    [20] 李亚军, 王先明, 程贤松, 魏盟智, 李江月, 朱为菊, 邓晓东.  海南岛海尾湾浮游植物群落结构及水质的调查 . 热带生物学报, 2020, 11(3): 257-265. doi: 10.15886/j.cnki.rdswxb.2020.03.001
  • 加载中
图(6) / 表 (3)
计量
  • 文章访问数:  244
  • HTML全文浏览量:  270
  • PDF下载量:  12
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-01-20
  • 录用日期:  2023-05-07
  • 修回日期:  2023-05-01
  • 网络出版日期:  2023-12-09
  • 刊出日期:  2024-01-16

海南岛近20年GPP变化格局及驱动因素分析

doi: 10.15886/j.cnki.rdswxb.20230006
    基金项目:  海南大学科研启动项目(KYQD(ZR)20069);海南省自然科学基金青年基金项目(320QN202);国家重点研发项目(2021YFD2200404)
    作者简介:

    雷济舟(1996−),男,海南大学生态与环境学院2020级生态学硕士研究生。E-mail:hnhn604@163.com

    通讯作者: 孙仲益(1989−),男,副教授,博导。研究方向:生态信息学。E-mail: gis.rs@hainanu.edu.cn
  • 中图分类号: X 82.826

摘要: 为探究气象因素与人类活动对近20年海南岛总初级生产力(Gross Primary Production,GPP)变化的相对贡献,首先利用Theil-Sen及Mann-Kendall方法获取海南岛GPP整体时空分布特征,在此基础上,分别以土地利用及覆被变化(Land Use and Cover Change,LUCC)指示人类活动,以气温(Air Temperature,Ta)、饱和水汽压差(Vapor Pressure Deficit,VPD)以及光合有效辐射( Photosynthetically Active Radiation,PAR)作为气象要素指标,通过空间统计与机器学习手段,构建海南岛GPP变化格局归因模型,量化驱动因素的相对贡献。结果表明:研究期间内,时间上,海南岛GPP呈现0.44 Tg ·a−1的显著增加趋势(P=0.024);空间上,海南岛87.8%面积的GPP表现增加趋势,海口及三亚周边等小部分区域(约9%)则表现为下降;海南岛土地利用共计转移15 528.40 km2,主要发生于林地地区,林地净增加642.88 km2,草地地区转移面积为4 759.28 km2,耕地地区转移的面积为4 051.23 km2;相较于人类活动,研究期年际间气象条件的差异是海南岛GPP变化的主导因素,但特殊年份中LUCC对海南省市县的影响则有所凸显。

English Abstract

雷济舟,崔嵬,朱济帅,等. 海南岛近20年GPP变化格局及驱动因素分析[J]. 热带生物学报,2024, 15(1):42−51. doi:10.15886/j.cnki.rdswxb.20230006 doi:  10.15886/j.cnki.rdswxb.20230006
引用本文: 雷济舟,崔嵬,朱济帅,等. 海南岛近20年GPP变化格局及驱动因素分析[J]. 热带生物学报,2024, 15(1):42−51. doi:10.15886/j.cnki.rdswxb.20230006 doi:  10.15886/j.cnki.rdswxb.20230006
LEI Jizhou, CUI Wei, ZHU Jishuai, ZHANG Runqing, ZHAO Junfu, ZHANG Jie, ZHANG Xiang, SUN Zhongyi. Analysis of varying patterns and driving factors of GPP in Hainan Island during the past 20 years[J]. Journal of Tropical Biology, 2024, 15(1): 42-51. doi: 10.15886/j.cnki.rdswxb.20230006
Citation: LEI Jizhou, CUI Wei, ZHU Jishuai, ZHANG Runqing, ZHAO Junfu, ZHANG Jie, ZHANG Xiang, SUN Zhongyi. Analysis of varying patterns and driving factors of GPP in Hainan Island during the past 20 years[J]. Journal of Tropical Biology, 2024, 15(1): 42-51. doi: 10.15886/j.cnki.rdswxb.20230006
  • 陆地生态系统(以下称为陆生系统)总初级生产力(Gross Primary Production,GPP)是指大气中的CO2通过植被光合作用进入陆生系统并转化为有机碳的总量[1],在一定程度上决定着碳汇时空分布格局[2],影响着生态系统物质及能量流动[3];为实现区域碳平衡,减缓气候变化起到关键作用[4],具有重要意义[5]。但受到生态系统及其与环境要素间相互作用的影响,GPP的时空分布具有较大的变异性与异质性特征[6],因此,对GPP分布格局的驱动因素定量分析将会有助于提升区域碳循环的理解,也更有助于揭示区域生态系统对人类活动与区域气候变化的响应与反馈。气候变化以及人类活动作为影响热带地区GPP分布格局的主要驱动因素,相关研究已取得丰富成果,对于前者主要集中于水分[7]、温度[8-9]、CO2浓度[10]以及太阳辐射[11]等气象因素对GPP时空格局的影响[12],已通过机理模型[13]、遥感模型[14]及统计方法[15]进行了量化分析与机理解释。而后者,相关研究也证实了诸如城市建设[16]、退耕还林还草[17]、刀耕火种[18]等人类活动会显著改变区域GPP分布[19],同时土地利用及覆被变化(Land Use and Cover Change,LUCC)作为人类活动的直接体现[20],是最为经典的表征指标[21]。热带地区水热条件充足,各种环境因素都会对GPP产生影响,许多研究对此进行了讨论,其中在亚马逊热带地区通过机理模型揭示出辐射是生产力变化的主导因素之一[22];热带地区水热条件虽然相对充足,但时空间的不均匀分配使得年际降水对GPP变化产生较为关键影响[23],同时GPP对温度的变化也有着较高的敏感性[24]。在全球尺度的研究中已有结果表明,温度、辐射以及水分对于GPP的影响程度分别为13.07%、−7.24%和11.74%,具体表现为低纬度地区水分对于GPP变化起到主导作用,中高纬度地区温度起到主导作用[6]。随着人类在热带地区的活动愈发频繁[25],越来越多的研究者发现开荒等农业生产行为致使的土地利用变化才是导致GPP改变的主要驱动因素,这一观点得到在亚马逊雨林、东南亚热带雨林等地的相关研究支持[25-26]。综上所述,热带地区植被结构复杂,环境因素多变,虽然科学界对GPP时空分布格局的驱动因素及响应机理已基本达成共识,但主导因素及其相对贡献大小依旧为争论焦点。

    海南岛是我国第一大热带岛屿,作为相对独立的地理单元是进行热带地区GPP研究的理想切入点;并且近年来随着海南自由贸易港等政策的落实[27-28],海南岛建设程度加深,人才引进以及城市扩张等一系列活动导致土地覆盖变化剧烈,为探究人类活动与气候变化对GPP的影响提供了契机。因此,本研究选择海南岛作为研究靶区,对其近20年GPP动态变化趋势进行归因分析,具体目标为:探究2001—2019年间海南岛GPP时空变化格局;量化气象要素与人类活动对于年际间GPP变化的贡献率,揭示主导因素。旨在为解明热带地区生态系统响应区域气候变化提供科学参考也为海南生态文明试验区的政策制定提供理论依据。

    • 海南岛(18°10′—20°10′N,108°37′—111°03′E),地处热带北缘,属热带季风海洋性气候,温差较低,全年高温,年平均气温22.5~25.6 ℃。年平均降雨约为1 640 mm,降水充足但分配不均,雨旱两季分明。海南岛地势为四周低平,中间高耸,呈穹隆山地形,中部山区以五指山、鹦哥岭为隆起核心,向外围逐级下降,由山地、丘陵、台地、平原构成环形层状地貌,梯级结构明显,土地覆盖类型丰富,以森林、农田以及草原为主,其中,森林生态系统包括常绿阔叶林、常绿针叶林、落叶阔叶林以及混交林等。

    • 海南岛土地利用类型数据集来源于MODIS(Moderate-resolution Imaging Spectroradiometer)Land Cover产品(MCD12Q1),时间分辨率为1 a,空间分辨率为500 m,本研究时间范围为2001—2019年[29]。本研究采用了IGBP(International Geosphere Biosphere Programme)分类方案所确定的17个土地覆盖方案。结合相关研究[30]以及本研究实际情况,对土地覆盖类型进行重分类,重分类情况见表1

      表 1  根据IGBP分类系统对土地利用类型的重分类

      土地利用
      类型编号
      重分类后土地
      利用类型
      IGBP分类
      1 林地 常绿阔叶林
      常绿针叶林
      落叶阔叶林
      混交林
      有林草地
      稀树草地
      2 灌丛 稀疏灌木丛及
      封闭灌木丛
      3 草地 草地
      4 耕地 耕地、
      耕地或自然植被
      5 城市 城市以及建筑用地
      6 裸地 裸地
      7 湿地 永久湿地及水域
    • 2001—2019年海南岛GPP数据集为MODIS产品MOD17A2HV006。该产品为4级标准产品[31],主要通过光能利用率模型进行计算,逐8天合成空间分辨率为500 m的GPP产品,从区域尺度到全球尺度已经被广泛应用[32]。本研究将GPP数据求和处理成月值尺度后转化为年尺度进行后续分析。

    • 本研究所使用2001—2019年逐月Ta、PAR以及相对湿度(Relative Humidity,RH)气象数据,其中Ta、RH数据均源自于国家地球系统科学数据中心(http://www.geodata.cn/),其空间分辨率均为1 km,是利用全国基本气象站观测数据进行空间插值并综合地形数据所绘制的空间分布数据[33-34]。根据Ta与RH数据并结合Tetens经验公式[35]对海南岛逐月VPD(饱和水汽压差)的数据进行计算;PAR(光合有效辐射)的数据为Global Land Surface Satellite(GLASS)产品(http://www.glass.umd.edu/index.html),其空间分辨率为0.01°,时间分辨率为1 d是基于多源遥感数据和地面实测数据进一步反演得到的长时间序列、高精度的全球地表遥感产品[36]

    • 首先以海南岛年际间用地类型未改变的栅格作为数据样本,利用其多年VPD、Ta及PAR年距平值作为解释变量,GPP年距平均值作为目标变量,通过高斯过程回归(Gaussian Process Regression, GPR)与随机森林(Random Forest,RF)两种机器学习算法构建预测模型,GPR是常用的监督分类学习方法,是一种基于贝叶斯方法的非参数概率模型,回归的目的是通过学习样本,经过训练得到输入变量与输出变量之间的函数关系,常常用于小样本回归分析。RF是经典的基于分类和回归树的集成学习算法,可以解释若干自变量对因变量的作用;其次,选择较好的机器学习算法后,利用年际间VPD、PAR以及Ta的差异作为输入数据,对其年际间气象要素差异所引起的GPP差异进行估算,完成气象因素对GPP影响的相对贡献率的计算;最后,以年际间用地类型发生改变的栅格作为研究目标,计算相邻年间LUCC改变的栅格其GPP的变化值,进而完成LUCC对GPP影响的相对贡献率的计算;最终进行海南岛GPP变化主导因素的探讨(图1)。

      图  1  研究流程图

    • 本研究通过GPR及RF机器学习算法构建模型以分析气象因素对GPP的相对贡献率。如表2所示,海南岛近20年间各用地类型均有不同程度改变,裸地与湿地类型未变栅格所占比例较低,用于训练模型的样本数量较少,因此选择GPR算法;而林地、草地、耕地等其他类型用地样本数量充足,利用RF能够取得较好效果。所有用地利用的模型训练均采用五折交叉验证进行精度的评估,以决定系数(Coefficient Of Determination,R2)、均方根误差(Root Mean Square Error,RMSE,记作SE)作为模型精度评价指标。

      表 2  模型评价精度表

      土地利用
      类型
      SE/
      (g·m−2·mon−1
      R2模型
      算法
      用地类型
      未变栅格
      数量/个
      林地9.600.60RF97534
      草地6.350.66RF2199
      耕地4.880.81RF7762
      城市5.410.62RF1568
      裸地3.390.29GPR33
      湿地10.380.12GPR1583
    • 为探究LUCC以及气象因素对于海南岛GPP的影响程度,本研究中LUCC以及气象因素相对贡献度计算公式如下:

      $$ {\text{Con}}{({\text{LUCC}})_{(i,i + 1)}} = \frac{{{\text{GPP}}{{({\text{LUCC}})}_{(i,i + 1)}}}}{{{\text{GP}}{{\text{P}}_i}}} $$ (1)
      $$ {\text{Con}}{({\text{Climate}})_{(i,i + 1)}} = \frac{{\displaystyle \sum\limits_{i = 1}^n {f({\text{VP}}{{\text{D}}_{(i,i + 1)}}{\text{T}}{{\text{a}}_{(i,i + 1)}}{\text{PA}}{{\text{R}}_{(i,i + 1)}})} }}{{{\text{GP}}{{\text{P}}_i}}} $$ (2)

      式中,${\text{Con}}{({\text{LUCC}})_{(i,i + 1)}}$为第i年至第i+1年间LUCC对于GPP变化的相对贡献率,$ {\text{GP}}{{\text{P}}_i} $i年海南岛GPP总值,${\text{GPP}}{({\text{LUCC}})_{(i,i + 1)}}$ii+1年间LUCC所驱动的GPP变化值;${\text{Con}}{({\text{Climate}})_{(i,i + 1)}}$i年至第i+1年间气象因素对GPP变化的相对贡献率,$ f $为1.3.1所介绍的机器学习模型。${\text{VP}}{{\text{D}}_{(i,i + 1)}}$${\text{T}}{{\text{a}}_{(i,i + 1)}}$ 以及${\text{PA}}{{\text{R}}_{(i,i + 1)}}$ 分别为第i年至第i+1年间VPD、Ta以及PAR的变化量,其中,i为年份;n为各土地利用类型,n=1,2……7。

    • 海南岛时间上,近20年期间GPP(以C表示,下同)时间变化如图2-a所示,GPP年际变化明显,总体上研究时期GPP呈现上升趋势(0.44 Tg·a−1)。2005年海南岛GPP值最低仅为56.46 Tg C,2017年后GPP值持续增加,在2019年达到峰值74.12 Tg C。其GPP主要由3种土地利用类型构成,林地为GPP值最高的土地利用类型,约占到全年GPP总值的85%,19年期间GPP值以0.46 Tg·a−1增加;耕地作为GPP值仅次于林地的土地利用类型,其年GPP值最大占比达到10.7%,耕地19年期间变化幅度不明显,其年GPP值以0.07 Tg·a−1增加;草地GPP值在19年期间持续降低(−0.11 Tg·a−1);空间上,近20年间海南岛87.8%的面积呈现出GPP增加趋势(图2-b);其中,海南岛北部以及东北部是GPP极显著增加的集中区,中部山区部分地区GPP无明显变化,而海口市与三亚市周边小部分区域的GPP表现为减少趋势,且碎片化程度较高,占比为8.9%。

      图  2  海南岛GPP时空间变化趋势

    • 对近20年海南岛土地利用变化数据进行分析(表3),研究时期内转移面积小于1 km2的分支不进行展示。结果显示,2001—2005年,海南岛土地利用转移总面积为2 736.94 km2。其中,林地的转出面积最大,净转出面积为309.13 km2,主要转化为耕地;2005—2010年土地转移总面积为4 164.04 km2,草地转化成林地以及耕地转化为林地为最主要的转出方式;2010—2015年期间,海南岛的土地覆盖变化主要发生在林地、草地以及耕地的地区面积分别为1 188.10 km2、1 074.56 km2以及1 070.85 km2;2015—2019年最主要的转移方式是林地转化为耕地和草地转化为林地,其转化净面积分别为541.25 km2和214.65 km2

      表 3  海南岛近20年间土地利用每5年的转移情况

      转移路径2001—2005
      面积/km2
      2005—2010
      面积/km2
      2010—2015
      面积/km2
      2015—2019
      面积/km2
      林地→草地517.43327.40325.38658.50
      林地→耕地695.82826.56794.801684.88
      林地→城市19.8549.4939.1740.23
      林地→裸地1.060.260.790.79
      林地→湿地10.0633.0827.5354.52
      灌丛→草地2.121.590.790.53
      草地→林地574.331237.59810.42873.14
      草地→耕地351.48215.97246.67358.63
      草地→城市2.127.414.233.44
      草地→裸地0.791.851.061.32
      草地→湿地3.9725.6711.6526.20
      耕地→林地329.25933.22742.131143.64
      耕地→草地181.30476.14328.19361.54
      耕地→城市1.593.440.530.26
      裸地→林地0.532.381.590.79
      裸地→草地0.531.592.651.85
      裸地→湿地0.792.121.321.59
      湿地→林地30.1713.2316.4140.49
      湿地→草地11.123.713.7112.17
      湿地→裸地1.320.530.531.06
    • 研究时期,海南岛LUCC所驱动的GPP变化主要发生在林地、草地以及耕地地区(图3)LUCC所驱动的GPP年际变化明显。其中,林地地区LUCC所驱动的GPP呈增加趋势,在2018—2019年林地地区LUCC所驱动的GPP达到最大值8.41×10−2 Tg C,占比57%;草地地区19年期间呈显著微弱增加趋势,2005—2006年草地地区达到LUCC所驱动的GPP峰值7.56×10−2 Tg C,2004—2005年草地地区LUCC所驱动的GPP仅为−5.99×10−2 Tg C;耕地地区在近20年内也是呈现出微弱增加趋势,在2008—2009年LUCC所驱动的GPP值最大。

      图  3  各土地利用类型在LUCC所驱动GPP变化

    • LUCC以及气象因素对于海南岛GPP相对贡献率的年际变化如图4所示,年际间LUCC与气象要素的相对贡献率同频率较高,但程度相差较大,气象因素主导海南岛年际间GPP变化。LUCC的相对贡献率近20年期间波动较为剧烈,其中2008—2009年期间相对贡献率达0.31%,而在2004—2005年LUCC对于GPP的相对贡献率最小为−0.24%;气象因素对年际间GPP变化的相对贡献率逐年差异明显,2005—2006年其正向相对贡献率达到最大值10.79%,而在2010—2011年气象要素对于GPP的影响程度为历年来最低值−14.99%。

      图  4  相对贡献率变化图

    • LUCC一直是陆生系统GPP的主要驱动因素。林地、草地以及耕地地区是海南岛LUCC发生的主要区域。研究时期LUCC共变化15 528.40 km2。其中,植树造林6 749.33 km2,累积GPP达到0.23 Tg C,这与相关研究结论[37]基本一致即植树造林对于GPP 增长有着较为良好的促进作用。这一现象的原因可能是椰林[38]、橡胶林[39]等林种具有较高的经济价值。此外,近52%的新林地资源主要来源于草地,改善了海南岛植被覆盖,强化了林地生态系统的防线保护作用;耕地也在此过程中净增加679.41 km2,GPP在此过程中增加1.63×10−2 Tg C。耕地的快速扩张一定程度改变原有的生态环境,由单一经济作物所构成的生态系统抵抗力稳定性降低可能会在未来引发水土流失,所以如何合理规划耕地资源是要深入考虑的。

      LUCC对GPP年际间变化的相对贡献率远低于气象因素,其一在于,海南岛林地面积高(79.1%),并可占全岛GPP的84.8%,并且主要集中在中部山区,由于天然林保护等政策[40-41],中部山区人类活动强度较低,所以在近20年来,林地并未发生较大的面积变化(近20年变化不足2%);二是由于MCD12Q1产品的将海南岛土地利用与覆盖分为14类,而为方便分析本研究将其重分类为7类,不同分类系统所造成的相对贡献率影响(图5-a)最高可相差3倍(2004—2005年)。尺度效应也是重要因素,如图5-b所示,将海南岛近20年逐年逐市县进行拆解,并使用IGBP原始分类系统,LUCC在各市县中对年际间GPP变化的相对贡献率普遍高于全岛平均值,LUCC相较于气象因素而言作用范围小,因此在进行全岛LUCC相对贡献率计算时不同地区间差异相互抵消,因此也造成LUCC的相对贡献率较低。LUCC相较于气象因素,对于小尺度范围的GPP年际间变化起到主要作用。LUCC与气象因素对于年际间GPP变化的相对贡献率具有较高的同频性,其原因在于相对贡献率的算法以前一年GPP为基准,而非年际间GPP差异为基准,好处在于避免出现年际间GPP相差极小而导致的极端贡献;但负效应变为同频效应以2005—2006年为例,极端事件第18号热带风暴“达维”对全岛GPP产生显著负面影响[42],GPP为近20年最低值,因此2005—2006年会产生极高的LUCC与气象因素的正向相对贡献,各市县贡献率受整体变化所影响。

      图  5  土地分类对相对贡献率的影响及各市(县)相对贡献率变化

    • 本研究仅选取Ta、VPD以及PAR[43]作为影响GPP年际间变化的主要气象驱动因素,未考虑如大气CO2浓度、氮沉降等环境因素,这为本研究带来一定不确定性,但根据Sun等人研究表明[6],环境要素中气象因素是环境因素中主导GPP年际间变化的关键因素,因此不确定性能够最小程度降低。此外,本研究选取VPD而非降水作为水分条件,主要由于像海南岛这类低纬度湿润地区,大气的湿润程度是限制植物生长的主要因素[44]。其次,对于气象因素间的交互作用,本研究所采用机器学习算法(GPR与RF),属于数据驱动模式,能够减弱自变量间的交互关系[15]

      气象因素相对贡献较高,分别选取正负相对贡献率最高的2008—2009年(2005—2006年为最高正相对贡献,但主要受极端事件“达维”影响,因此选取2008—2009年)与2010—2011年为例;2008—2009年,PAR、Ta及VPD均呈现超距平的转好态势;而2010—2011年PAR与Ta呈现出远超平均年际间变化的降幅,同时水分制约也略有所增加(图6);可见本研究所提出的基于机器学习的数据驱动算法能够捕捉气象要素对GPP的影响。除气象要素与LUCC外,人类的管理与经营活动,如耕地农业管理及政策支持[45-46]、橡胶林下经营[47]、天然林的保护[48]、极端事件以及病虫害等异常事件也会对海南岛GPP产生影响[49]

      图  6  海南岛各气象要素变化图

    • 首先本研究所使用的GPP与LUCC数据来源于MODIS产品,MODIS作为光学遥感产品应用于热带地区时不可避免精度会降低[50],这是传统光学遥感产品共同面临问题。MODIS的LUCC产品在全球尺度的精度为73.6%[51],本研究将可能出现异物同谱的相似类进行了重分类,共7大类,这能够大大降低其分类精度所带来的不确定性(图5-a),这里不能够忽视。其次,MODIS的GPP产品基于光能利用率算法,是现阶段主流大区域尺度的遥感GPP算法,包括EC-LUE[52]、CASA以及VPM[6]等多种形式,其中可吸收光合有效辐射比率(fAPAR,fraction of Absorbed PAR)是模型核心部分,无论何种形式模型均将fAPAR考虑为关于NDVI等植被指数的函数;长时间序列植被指数主要源于MODIS产品,因此本研究中直接使用了MODIS的GPP产品,避免二次计算时增加不确定性。此外关于计算方案引入的不确定性主要为,当年际间LUCC发生改变时,对GPP的影响则全部归于LUCC的贡献;但实则气象要素依旧对该区域产生影响,可此部分的影响并未考虑到本研究中,低估了气象因素造成的总体影响。综上,鉴于本研究的时间范围与尺度,结合各遥感产品的精度,最终确定MCD12Q1为研究数据,该数据产品在国内的研究已广泛应用于国内相关研究[53],并取得较为可靠的成果。考虑到海南岛全年较高的云覆盖,产品精度有限。未来要进一步研究LUCC对GPP影响等问题时,更高精确度土地利用数据产品是不可或缺的;本研究使用机器学习模型虽然可以在一定程度上将气象因素与人类活动进行解耦并单独分析其对GPP的影响,实际环境中GPP的影响是多元驱动的,无法对影响因素的共同作用进行驱动分析。本研究中假设LUCC发生改变时GPP的变化归因为LUCC所引起,将气象因素作用归于LUCC贡献中,这也带来了一定不确定性。

参考文献 (53)

目录

    /

    返回文章
    返回