留言板

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

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

白洋淀大清河流域油松精准适宜性空间分布

郭虹扬 史明昌 杨建英 陈春阳

郭虹扬, 史明昌, 杨建英, 陈春阳. 白洋淀大清河流域油松精准适宜性空间分布[J]. 浙江农林大学学报, 2021, 38(6): 1100-1108. doi: 10.11833/j.issn.2095-0756.20200599
引用本文: 郭虹扬, 史明昌, 杨建英, 陈春阳. 白洋淀大清河流域油松精准适宜性空间分布[J]. 浙江农林大学学报, 2021, 38(6): 1100-1108. doi: 10.11833/j.issn.2095-0756.20200599
GUO Hongyang, SHI Mingchang, YANG Jianying, CHEN Chunyang. Precise spatial distribution of suitability of Pinus tabulaeformis in Daqing River Basin, Baiyangdian[J]. Journal of Zhejiang A&F University, 2021, 38(6): 1100-1108. doi: 10.11833/j.issn.2095-0756.20200599
Citation: GUO Hongyang, SHI Mingchang, YANG Jianying, CHEN Chunyang. Precise spatial distribution of suitability of Pinus tabulaeformis in Daqing River Basin, Baiyangdian[J]. Journal of Zhejiang A&F University, 2021, 38(6): 1100-1108. doi: 10.11833/j.issn.2095-0756.20200599

白洋淀大清河流域油松精准适宜性空间分布

doi: 10.11833/j.issn.2095-0756.20200599
基金项目: 国家水体污染控制与治理科技重大专项(2018ZX07110001)
详细信息
    作者简介: 郭虹扬(ORCID: 0000-0002-8938-4789),从事“3S”技术集成开发与应用研究。E-mail: hongyangguo@yeah.net
    通信作者: 史明昌(ORCID: 0000-0001-8289-552X),教授,从事“3S”技术集成开发与应用研究。E-mail: shimc@bjfu.edu.cn
  • 中图分类号: S791.254;S718.5

Precise spatial distribution of suitability of Pinus tabulaeformis in Daqing River Basin, Baiyangdian

  • 摘要:   目的  油松Pinus tabulaeformis作为白洋淀大清河流域上游山区主要造林树种,研究其精准适宜性空间分布对流域内水源涵养林的建设与提升具有生产指导作用。  方法  依据植物生态学、植物地理学原理,采用空间限制性因子叠加及最大熵模型(MaxEnt),以白洋淀大清河流域油松为研究对象,基于野外190个样方数据,结合24个地理环境变量(高程、坡度、坡向、气象、岩石类型、土壤类型等),其中气象因子通过多元线性回归克里金插值方法(multiple linear regression Kriging, MLRK)得到,对白洋淀大清河流域油松的适宜性空间分布进行了研究。  结果  油松适宜性预测模型效果达到准确水平(受试者工作曲线ROC值为0.869),积温、年均最大日降水量、高程、年均气温标准差为流域油松适宜性主导地理环境变量,油松高适宜区主要分布在灵丘县、涞源县、涞水县、易县、房山区等。  结论  白洋淀大清河流域油松的适宜性条件为:积温2 800~3 800 ℃、海拔500~1 200 m、年均日最大降水量70~110 mm、年均气温标准差0.31~0.42 ℃,按照该适宜性条件进行油松林营造,可充分发挥其有益机能。图8表2参30
  • 图  1  白洋淀大清河流域油松样方分布示意图

    Figure  1  Distribution of P. tabulaeformis quadrats in Baiyangdian Daqing River Basin

    图  2  地形变量示意图

    Figure  2  Terrain variables

    图  3  土壤地质变量示意图

    Figure  3  Soil geological variables

    图  4  气象变量示意图

    Figure  4  Meteorological variables

    图  5  ROC曲线图

    Figure  5  ROC curve

    图  6  刀切法地理环境变量贡献图

    Figure  6  Jackknife environmental variable contribution graph

    图  7  主导地理环境变量反馈曲线

    Figure  7  Leading environmental variable feedback curve

    图  8  油松精准适宜性分布及细节放大图

    Figure  8  Accurate suitability distribution and detailed enlarged drawing of P. tabulaeformis

    表  1  地理环境变量数据

    Table  1.   Geographical environment data

    变量类型变量名称变量代码单位变量类型变量名称变量代码单位
    地形  高程          1m 气象多年平均气温     13
    坡度          2(°)多年年极端最高气温  14
    坡向          3多年年极端最低气温  15
    土壤地质土壤类型        4多年年平均气温标准差 16
    岩石类型        5多年平均气温年较差  17
    气象  多年日平均降水量    6mm多年年平均气温日较差 18
    多年年平均最大日降水量 7mm多年最冷月平均气温  19
    多年年平均降水量标准差 8mm多年最热月平均气温  20
    多年平均最湿季降水量  9mm多年平均地面温度   21
    多年平均最干季降水量 10mm多年年平均积温    22
    多年平均相对湿度   11%多年平均日照时数   23h
    多年年平均无霜期日数 12d多年平均年蒸发量   24mm
      说明:多年平均数据均为1981−2019年的数据(共39 a)
    下载: 导出CSV

    表  2  各县高适宜区面积及比例

    Table  2.   Area and proportion of high suitable area in each county

    适生区面积/km2占总面积比例/%
    涞源县462.601.00
    灵丘县402.650.87
    易县 366.420.80
    涞水县181.620.39
    清苑县74.380.16
    涿鹿县64.290.14
    下载: 导出CSV
  • [1] 夏军, 张永勇. 雄安新区建设水安全保障面临的问题与挑战[J]. 中国科学院院刊, 2017, 32(11): 1199 − 1205.

    XIA Jun, ZHANG Yongyong. Problems and challenges faced by the construction of water safety security in Xiong’an New Area [J]. J Chin Acad Sci, 2017, 32(11): 1199 − 1205.
    [2] 姜凤岐, 曾德慧, 于占源. 从恢复生态学视角透析防护林衰退及其防治对策: 以章古台地区樟子松林为例[J]. 应用生态学报, 2006, 17(12): 2229 − 2235.

    JIANG Fengqi, ZENG Dehui, YU Zhanyuan. Decline of protective forest and its prevention strategies from viewpoint of restoration ecology: taking Pinus sylvestris var. mongolica pantation in Zhanggutai as an example [J]. Chin J Appl Ecol, 2006, 17(12): 2229 − 2235.
    [3] 黄森旺, 李晓松, 吴炳方, 等. 近25年三北防护林工程区土地退化及驱动力分析[J]. 地理学报, 2012, 67(5): 589 − 598.

    HUANG Senwang, LI Xiaosong, WU Bingfang, et al. The distribution and drivers of land degradation in the Three-North Shelter Forest region of China during 1982−2006 [J]. Acta Geogr Sin, 2012, 67(5): 589 − 598.
    [4] 王耀, 张昌顺, 刘春兰, 等. 三北防护林体系建设工程区森林水源涵养格局变化研究[J]. 生态学报, 2019, 39(16): 5847 − 5856.

    WANG Yao, ZHANG Changshun, LIU Chunlan, et al. Research on the pattern and change of forest water conservation in Three-North Shelterbelt Forest program region, China [J]. Acta Ecol Sin, 2019, 39(16): 5847 − 5856.
    [5] 胡耀升, 翟洪波, 田野. 太行山绿化工程建设可持续性分析[J]. 林业经济, 2017, 39(9): 48 − 52.

    HU Yaosheng, ZHAI Hongbo, TIAN Ye. Analysis on sustainability of Taihang Mountain Greening Program construction [J]. For Econ, 2017, 39(9): 48 − 52.
    [6] 张春霞, 冯自茂, 李文鑫, 等. 陕西黄陵油松人工林立地类型划分及评价[J]. 西北农林科技大学学报(自然科学版), 2021, 49(1): 55 − 63.

    ZHANG Chunxia, FENG Zimao, LI Wenxin, et al. Site classification and evaluation of Pinus tabulaeformis plantation in Huangling, Shaanxi [J]. J Northwest A&F Univ Nat Sci Ed, 2021, 49(1): 55 − 63.
    [7] 赖文豪, 席沁, 武海龙, 等. 内蒙古兴和县低山丘陵立地类型划分与林草适宜性评价[J]. 浙江农林大学学报, 2018, 35(2): 331 − 339.

    LAI Wenhao, XI Qin, WU Hailong, et al. Site classification type and vegetation suitability evaluation for hilly land in Xinghe, Inner Mongolia [J]. J Zhejiang A&F Univ, 2018, 35(2): 331 − 339.
    [8] 武吉华. 植物地理学[M]. 4版. 北京: 高等教育出版社, 2004.
    [9] 苏忠. 城市峡谷中GPS/BeiDou伪距单点定位性能分析[J]. 测绘通报, 2019(增刊): 30 − 35.

    SU Zhong. Positioning performance analysis of GPS/BeiDou pseudorange single point positioning in the urban canyon[J]. 2019 (suppl): 30 − 35.
    [10] JIN Qiutong, ZHANG Jutao, SHI Mingchang, et al. Estimating loess plateau average annual precipitation with Multiple Linear Regression Kriging and Geographically Weighted Regression Kriging [J]. Water, 2016, 6(8): 266 − 286.
    [11] HENGL T. A practical guide to geostatistical mapping of environmental variables [J]. Geoderma, 2007, 140: 417 − 427.
    [12] WANG Ku, ZHANG Chuanrong, LI Weidong. Comparison of geographically weighted regression and regression Kriging for estimating the spatial distribution of soil organic matter [J]. GISci Remote Sensing, 2013, 49(6): 915 − 932.
    [13] 张年国, 王娜. 基于大数据的沈阳国土空间规划用地适宜性评价探讨[J]. 价值工程, 2019, 38(14): 10 − 13.

    ZHANG Nianguo, WANG Na. Discussion on suitability evaluation of shenyang land use planning based on big data [J]. Value Eng, 2019, 38(14): 10 − 13.
    [14] 陈有民. 园林树木学[M]. 2版. 北京: 中国林业出版社, 2011.
    [15] 任宪威. 树木学(北方本)[M]. 北京: 中国林业出版社, 1997.
    [16] PHILLIPS S J, ANDERSON R P, SCHAPIRE R E. Maximum entropy modeling of species geographic distributions [J]. Ecol Modelling, 2006, 190(3/4): 231 − 259.
    [17] ELITH J, PHILLIPS S J, HASTIE T, et al. A statistical explanation of MaxEnt for ecologists [J]. Diversity Distrib, 2011, 17(1): 43 − 57.
    [18] QIN Aili, JIN Kun, BATSAIKHAN M E, et al. Predicting the current and future suitable habitats of the main dietary plants of the Gobi Bear using MaxEnt modeling [J]. Global Ecol Conserv, 2020, 22: e1032. doi: 10.1016/j.gecco.2020.e01032.
    [19] 邱浩杰, 孙杰杰, 徐达, 等. 基于MaxEnt模型预测鹅掌楸在中国的潜在分布区[J]. 浙江农林大学学报, 2020, 37(1): 1 − 8.

    QIU Haojie, SUN Jiejie, XU Da, et al. MaxEnt model-based prediction of potential distribution of Liriodendron chinense in China [J]. J Zhejiang A&F Univ, 2020, 37(1): 1 − 8.
    [20] HANLEY J A, MCNEIL B J. The meaning and use of the area under a receiver operating characteristic (ROC) curve [J]. Radiology, 1982, 143(1): 29 − 36.
    [21] 邓飞, 李晓兵, 王宏, 等. 基于MaxEnt模型评价紫花苜蓿在锡林郭勒盟的分布适宜性及主导因子[J]. 草业科学, 2014, 31(10): 1840 − 1847.

    DENG Fei, LI Xiaobing, WANG Hong, et al. The suitability of geographic distribution and the dominant factors of alfalfa based on MaxEnt model in Xilin Gol [J]. Grassl Sci, 2014, 31(10): 1840 − 1847.
    [22] PEARSON R G, RAXWORTHY C J, NAKAMURA M, et al. Predicting species distributions from small numbers of occurrence records: a test case using cryptic geckos in Madagascar [J]. J Biogeogr, 2007, 34(1): 102 − 117.
    [23] 李培琳. 浙江省森林立地分类与杉木适宜性研究[D]. 杭州: 浙江农林大学, 2018.

    LI Peilin. Forest Site Classification and Suitability of Cunninghamia lanceolate (Lamb.) Hook in Zhejiang Province[D]. Hangzhou: Zhejiang A&F University, 2018.
    [24] 麻亚鸿. 基于最大熵模型(MaxEnt)和地理信息系统(ArcGIS)预测藓类植物的地理分布范围[D]. 上海: 上海师范大学, 2013.

    MA Yahong. Applying MaxEnt and ArcGIS to Predict Mosses Geographic Distribution Range: A Case Study of Huaping Nature Reserve, Guangxi[D]. Shanghai: Shanghai Normal University, 2013.
    [25] 周立江. 低效林评判与改造途径的探讨[J]. 四川林业科技, 2004(1): 16 − 21.

    ZHOU Lijiang. Discussion on judgment and rebuilding approaches of low-efficiency forest [J]. J Sichuan For Sci Technol, 2004(1): 16 − 21.
    [26] 李铁民. 太行山低质低效林判定标准[J]. 山西林业科技, 2000(3): 1 − 8.

    LI Tiemin. The criterion of low-quality-and-benefit forests in Taihang Mt. [J]. Shanxi For Sci Technol, 2000(3): 1 − 8.
    [27] 魏安琪. 陕西吴起退耕还林植被恢复及常见造林树种适生性研究[D]. 北京: 北京林业大学, 2019.

    WEI Anqi. Study on Vegetation Restoration and Adaptability of Common Afforestation Tree Species in Converted Farmland to Forestry Region: A Case of Wuqi County, Shaanxi Province[D]. Beijing: Beijing Forestry University, 2019.
    [28] 王志明. 基于GIS的晋西黄土区立地因子分析及类型划分研究[D]. 北京: 北京林业大学, 2012.

    WANG Zhiming. Based on GIS Technology of Research of Forest Site Classification in the Semiarid Loss Hilly-Gully Region[D]. Beijing: Beijing Forestry University, 2012.
    [29] 李超, 裴顺祥, 张连金, 等. 北京油松人工林竞争指数的适用性评价[J]. 浙江农林大学学报, 2019, 36(6): 1115 − 1124.

    LI Chao, PEI Shunxiang, ZHANG Lianjin, et al. Applicability evaluation of competition indexes for Pinus tabuliformis plantations in Beijing [J]. J Zhejiang A&F Univ, 2019, 36(6): 1115 − 1124.
    [30] 侯淑艳. 北京市低山区低效人工林结构特征与评价研究[D]. 北京: 北京林业大学, 2013.

    HOU Shuyan. Study on the Assessment and Structure Feature of Low Function Plantation in the Low Mountain of Beijing[D]. Beijing: Beijing Forestry University, 2013.
  • [1] 王祥福, 李愿会, 王维枫, 孙杰杰, 王倩, 董文婷, 王荣女, 杨娅琪.  土地利用和气候变化对青海省雪豹潜在适宜生境的影响 . 浙江农林大学学报, doi: 10.11833/j.issn.2095-0756.20230259
    [2] 张世林, 高润红, 高明龙, 韩淑敏, 张文英, 赵静.  气候变化背景下中国樟子松潜在分布预测 . 浙江农林大学学报, 2023, 40(3): 560-568. doi: 10.11833/j.issn.2095-0756.20220451
    [3] 汪雁楠, 赖国桢, 黄建建, 刘丽婷, 余良森, 温强, 龚春.  浙江红花油茶潜在适生区分布及其对未来气候变化的响应 . 浙江农林大学学报, 2022, 39(5): 989-997. doi: 10.11833/j.issn.2095-0756.20210641
    [4] 靳孟理, 胡俊, 齐实, 逯进生, 李月, 蒋九华.  北京市侧柏低效林物种多样性改造策略 . 浙江农林大学学报, 2020, 37(1): 27-35. doi: 10.11833/j.issn.2095-0756.2020.01.004
    [5] 杜雨菲, 吴保国, 陈玉玲.  基于机器学习算法的广西桉树适宜性研究 . 浙江农林大学学报, 2020, 37(1): 122-128. doi: 10.11833/j.issn.2095-0756.2020.01.016
    [6] 李超, 裴顺祥, 张连金, 郭嘉, 辛学兵.  北京油松人工林竞争指数的适用性评价 . 浙江农林大学学报, 2019, 36(6): 1115-1124. doi: 10.11833/j.issn.2095-0756.2019.06.008
    [7] 刘崴, 魏天兴, 朱清科.  半干旱黄土丘陵区河北杨和油松生长季树干液流特征 . 浙江农林大学学报, 2018, 35(6): 1045-1053. doi: 10.11833/j.issn.2095-0756.2018.06.007
    [8] 王金池, 冉啟香, 邓华锋, 黄国胜, 王雪军.  基于度量误差方法的油松林分生长模型 . 浙江农林大学学报, 2018, 35(1): 68-74. doi: 10.11833/j.issn.2095-0756.2018.01.009
    [9] 滕飞, 刘勇, 王琰, 胡嘉伟, 孙巧玉, 万芳芳, 张劲.  底部渗灌下容器类型及规格对油松根系结构的影响 . 浙江农林大学学报, 2017, 34(3): 449-458. doi: 10.11833/j.issn.2095-0756.2017.03.010
    [10] 王巧, 聂鑫, 孙德浩, 王华田, 孟先鹏, 曹桂萍, 李健, 黄昌豹.  基于AHP-模糊综合评价法的泰山油松古树树势评价 . 浙江农林大学学报, 2016, 33(1): 137-146. doi: 10.11833/j.issn.2095-0756.2016.01.019
    [11] 武亚敬, 毕君, 李秋艳, 高红真.  油松容器育苗香菇渣基质粒径配方筛选 . 浙江农林大学学报, 2015, 32(3): 483-487. doi: 10.11833/j.issn.2095-0756.2015.03.023
    [12] 史锋厚, 朱灿灿, 沈永宝, 施季森.  超声波与渗透调节处理对油松种子萌发的影响 . 浙江农林大学学报, 2011, 28(4): 545-549. doi: 10.11833/j.issn.2095-0756.2011.04.004
    [13] 邵东华, 任琴, 宁心哲, 白淑兰.  油松和虎榛子不同林型根系分泌物组分及化感效应 . 浙江农林大学学报, 2011, 28(2): 333-338. doi: 10.11833/j.issn.2095-0756.2011.02.026
    [14] 马丰丰, 贾黎明.  北京地区侧柏、油松带皮胸径与去皮胸径的关系 . 浙江农林大学学报, 2009, 26(1): 13-16.
    [15] 刘永红, 杨培华, 韩创举, 樊军锋, 李新会, 李安平, 杨世荣.  油松不同种源种实性状的变异分析 . 浙江农林大学学报, 2008, 25(2): 163-168.
    [16] 景丽, 朱志红, 王孝安, 郭华.  秦岭油松人工林与次生林群落特征比较 . 浙江农林大学学报, 2008, 25(6): 711-717.
    [17] 逯军峰, 王辉, 曹靖, 袁宏波.  不同林龄油松人工林枯枝落叶层持水性及养分含量 . 浙江农林大学学报, 2007, 24(3): 319-325.
    [18] 杨培华, 樊军锋, 刘永红, 韩创举, 段乐, 吕晓锋, 肖非, 郭宏杰.  油松第2 代无性系种子园营建技术 . 浙江农林大学学报, 2005, 22(2): 157-160.
    [19] 刘永红, 樊军锋, 杨培华, 韩创举.  油松单亲子代苗期生长性状遗传分析 . 浙江农林大学学报, 2005, 22(5): 513-517.
    [20] 周永学, 樊军锋, 杨培华, 高建社, 刘永红.  奥地利黑松与油松1 年生苗生长和生物量对比分析 . 浙江农林大学学报, 2003, 20(4): 438-441.
  • 加载中
  • 链接本文:

    https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.20200599

    https://zlxb.zafu.edu.cn/article/zjnldxxb/2021/6/1100

图(8) / 表(2)
计量
  • 文章访问数:  568
  • HTML全文浏览量:  146
  • PDF下载量:  77
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-09-17
  • 修回日期:  2021-05-07
  • 网络出版日期:  2021-06-29
  • 刊出日期:  2021-12-08

白洋淀大清河流域油松精准适宜性空间分布

doi: 10.11833/j.issn.2095-0756.20200599
    基金项目:  国家水体污染控制与治理科技重大专项(2018ZX07110001)
    作者简介:

    郭虹扬(ORCID: 0000-0002-8938-4789),从事“3S”技术集成开发与应用研究。E-mail: hongyangguo@yeah.net

    通信作者: 史明昌(ORCID: 0000-0001-8289-552X),教授,从事“3S”技术集成开发与应用研究。E-mail: shimc@bjfu.edu.cn
  • 中图分类号: S791.254;S718.5

摘要:   目的  油松Pinus tabulaeformis作为白洋淀大清河流域上游山区主要造林树种,研究其精准适宜性空间分布对流域内水源涵养林的建设与提升具有生产指导作用。  方法  依据植物生态学、植物地理学原理,采用空间限制性因子叠加及最大熵模型(MaxEnt),以白洋淀大清河流域油松为研究对象,基于野外190个样方数据,结合24个地理环境变量(高程、坡度、坡向、气象、岩石类型、土壤类型等),其中气象因子通过多元线性回归克里金插值方法(multiple linear regression Kriging, MLRK)得到,对白洋淀大清河流域油松的适宜性空间分布进行了研究。  结果  油松适宜性预测模型效果达到准确水平(受试者工作曲线ROC值为0.869),积温、年均最大日降水量、高程、年均气温标准差为流域油松适宜性主导地理环境变量,油松高适宜区主要分布在灵丘县、涞源县、涞水县、易县、房山区等。  结论  白洋淀大清河流域油松的适宜性条件为:积温2 800~3 800 ℃、海拔500~1 200 m、年均日最大降水量70~110 mm、年均气温标准差0.31~0.42 ℃,按照该适宜性条件进行油松林营造,可充分发挥其有益机能。图8表2参30

English Abstract

郭虹扬, 史明昌, 杨建英, 陈春阳. 白洋淀大清河流域油松精准适宜性空间分布[J]. 浙江农林大学学报, 2021, 38(6): 1100-1108. doi: 10.11833/j.issn.2095-0756.20200599
引用本文: 郭虹扬, 史明昌, 杨建英, 陈春阳. 白洋淀大清河流域油松精准适宜性空间分布[J]. 浙江农林大学学报, 2021, 38(6): 1100-1108. doi: 10.11833/j.issn.2095-0756.20200599
GUO Hongyang, SHI Mingchang, YANG Jianying, CHEN Chunyang. Precise spatial distribution of suitability of Pinus tabulaeformis in Daqing River Basin, Baiyangdian[J]. Journal of Zhejiang A&F University, 2021, 38(6): 1100-1108. doi: 10.11833/j.issn.2095-0756.20200599
Citation: GUO Hongyang, SHI Mingchang, YANG Jianying, CHEN Chunyang. Precise spatial distribution of suitability of Pinus tabulaeformis in Daqing River Basin, Baiyangdian[J]. Journal of Zhejiang A&F University, 2021, 38(6): 1100-1108. doi: 10.11833/j.issn.2095-0756.20200599
  • 白洋淀及大清河水系是雄安新区水资源安全保障体系的重要部分,但气候干旱、上游土壤流失严重、森林水源涵养功能低下,导致白洋淀干淀、上游断流,难以为新区发展提供保障[1]。土地或岩石的裸露、造林树种的不适宜以及森林空间结构不良等是森林水源涵养功能低下的主要原因[2-3]。早期植树造林时,主要根据区划及经验进行判断,导致油松Pinus tabulaeformis等部分防护林违背适宜性原则,致使其逐渐衰退为低效林,未达预期效果[4-5]。张春霞等[6]以海拔、坡度、坡向为主要立地因子对陕西黄陵油松进行立地质量评价,赖文豪等[7]通过地形、土壤、气象等立地因子分析了油松的适宜性情况。可见,已有研究主要集中在小流域等小尺度范围内,难以为当地区域防护林的营造提供支撑,且各立地因子精度较低,如仅用个别站点数据直接代表整个区域,或克里金插值(Kriging)的气象因子站点数据较少,不够精准化,进而导致适宜性分析结果无法贴合生产实际。本研究基于植物生态学特性,考虑了地形、土壤地质、气象等因子,且将各因子落实到30 m分辨率,采用了空间限制性因子叠加和最大熵模型(MaxEnt),对较大尺度下油松的精准适宜性空间分布进行了研究,以期为白洋淀大清河流域造林规划提供科学依据。

    • 大清河地处海河流域中部(38°10′~40°10′N,113°39′~117°34′E),西起太行山,东临渤海湾,北临永定河,南至子牙河。河系跨山西、河北、北京、天津4省市,总面积43 060 km2,其中山区18 659 km2,丘陵平原24 401 km2,分别占河系总面积的43.33%和56.67%。大清河流域地处温带半干旱大陆性季风气候区,多年平均气温为6.58~13.26 ℃,多年平均降水量为365~626 mm,年内降水量极为不均,多年平均积温为858~4 664 ℃,多年平均无霜期为73~217 d。研究区内优势植物有华北落叶松Larix principis-rupprechtii、油松、侧柏Platycladus orientalis、杨树Populus spp.、刺槐Robinia pseudoacacia、臭椿Ailanthus altissima、绣线菊Spiraea salicifolia、荆条Vitex negundo var. heterophylla、胡枝子Lespedeza bicolor等。

    • 本研究数据分为地形、土壤地质、气象等3类,其中:地形数据来源于地理空间数据云(http://www.gscloud.cn/);土壤地质数据来源于中国寒区旱区科学数据中心(http://westdc.westgis.ac.cn)和全国地质资料馆(http://ngac.org.cn/Map/List);气象数据为站点数据,来源于中国气象数据网(http://data.cma.cn/)各数据集。站点数量为79个,其中流域内36个,流域外43个,各地理环境变量如表1

      表 1  地理环境变量数据

      Table 1.  Geographical environment data

      变量类型变量名称变量代码单位变量类型变量名称变量代码单位
      地形  高程          1m 气象多年平均气温     13
      坡度          2(°)多年年极端最高气温  14
      坡向          3多年年极端最低气温  15
      土壤地质土壤类型        4多年年平均气温标准差 16
      岩石类型        5多年平均气温年较差  17
      气象  多年日平均降水量    6mm多年年平均气温日较差 18
      多年年平均最大日降水量 7mm多年最冷月平均气温  19
      多年年平均降水量标准差 8mm多年最热月平均气温  20
      多年平均最湿季降水量  9mm多年平均地面温度   21
      多年平均最干季降水量 10mm多年年平均积温    22
      多年平均相对湿度   11%多年平均日照时数   23h
      多年年平均无霜期日数 12d多年平均年蒸发量   24mm
        说明:多年平均数据均为1981−2019年的数据(共39 a)

      在白洋淀大清河流域上游山区,以油松为优势种,随机选取不同高程、坡度、坡向、坡位、土壤类型、龄组的样方进行调查,高程按照中山、低山区间结合太行山植物垂直地带性进行考虑[8],记录树高、胸径、郁闭度、植被覆盖度、林下灌草种类及盖度和地理位置等信息。在样方调查时,山脉、高大植物会遮挡全球定位系统(GPS)卫星信息接收[9],故地理位置信息通过手持GPS仪结合卫星地图比对位置并标注后自动获取,其精确度要高于手持GPS仪直接获取的位置数据,样方分布如图1所示。根据植物学和植物地理学等原理,在筛选油松可能适宜性空间单元时,考虑人为干预较少的样方和海拔100 m以上地区,因此,本研究调查样方主要分布于流域西部山区,共计190个30 m×30 m油松样方数据(图1)。

      图  1  白洋淀大清河流域油松样方分布示意图

      Figure 1.  Distribution of P. tabulaeformis quadrats in Baiyangdian Daqing River Basin

    • 采用多元线性回归克里金法(multiple linear regression Kriging, MLRK)对气象站点数据进行插值,得到连续气象空间分布数据,其主要公式如下[10-11]

      $$ \hat {{z}}\left( {{{{s}}_{\bf{0}}}} \right) = \mathop \sum \limits_{{{k}} = {\bf{0}}}^{{\rho }} {\hat {{\beta }}_{{k}}}{{{q}}_{{k}}}\left( {{{{s}}_{\bf{0}}}} \right) + \mathop \sum \limits_{{{i}} = 1}^{{n}} {{{\lambda }}_{{i}}}{{e}}\left( {{{{s}}_{{i}}}} \right)\text{。} $$ (1)

      式(1)中:$\hat z\left( {{s_0}} \right)$表示预测位置$ {{{{s}}_{\bf{0}}}}$点的插值结果,$ \displaystyle\sum \limits_{{{k}} = {\bf{0}}}^{{\rho }} {\hat {{\beta }}_{{k}}}{{{q}}_{{k}}}\left( {{{{s}}_{\bf{0}}}} \right) $是通过回归拟合的确定性部分;$ \displaystyle\sum \limits_{{{i}} = 1}^{{n}} {{{\lambda }}_{{i}}}{{e}}\left( {{{{s}}_{{i}}}} \right)$是普通克里金(OK)对回归残差的插值结果部分;$ {{k}}$表示回归拟合时的位置序号,$ \; \rho$表示空间位置总数,$ \; {\hat {{\beta }}_{{k}}}$是回归模型的系数,当k=0时,$ \; {\hat{{\beta }}_{\bf{0}}}$是截距;$i $表示回归残差插值时的位置序号,$ n$表示空间位置总数,${{{q}}_{{k}}}\left( {{{{s}}_{\bf{0}}}} \right) $是预测位置$ {{{{s}}_{\bf{0}}}}$点辅助变量的值,$ {{{\lambda }}_{{i}}}$是由回归残差的空间相关性结构确定的普通克里金插值的权重;$ {{e}}\left( {{{{s}}_{{i}}}} \right)$是位置$ {{{{s}}_{{i}}}}$处的残差。

      在进行气象数据插值时,将高程、坡度、坡向数据标准化后作为辅助变量,利用逐步回归法筛选后进行回归拟合,计算多元线性回归在气象站点位置的残差,再利用普通克里金法对回归残差进行插值,将每个预测点的确定部分和残差插值结果部分相加得到各点的预测值[12]。本研究在插值时考虑到研究区面积较大且跨越不同气候区,故将研究区按照气候区划分为2个区,即半湿润大区与半干旱大区,分别扩宽30 km进行插值。为提升流域边缘结果的准确性,站点除选取流域内的36个站点外,另选取流域周围43个站点纳入插值过程,插值完成后再进行拼接得到连续的气象要素空间分布数据,分辨率为30 m,根据精度检验,各气象数据MLRK插值结果交叉验证值R2均大于0.9,插值结果达到准确水平。

    • 根据油松的生态习性,利用限制性因子叠加[13],得到油松生长的可能适宜单元。油松是中国特有树种,主产华北、西北;生于海拔100~2 600 m,为研究区中、低山区重要的造林树种。喜光,适干冷气候,产区极端低温−25 ℃,年降水量400~750 mm,耐干旱,也能生长在年降水量250 mm的部分山区;在酸性、中性、钙质土壤中均能生长,极耐干瘠,石质荒山亦能生长,不耐水涝和盐碱[14-15]。考虑到以上生态习性,空间叠加模型为:

      $$ {P_{{\rm{SG}}}} = D \times T \times S \times R\text{。} $$ (2)

      式(2)中:$ {P_{{\rm{SG}}}}$表示油松可能适宜单元,当$ {P_{{\rm{SG}}}}=1$时,为油松可能适宜空间单元,当$ {P_{{\rm{SG}}}}=0$,为不适宜;当高程≥100 m时,D取值为1,反之则取值为0;当极端最低气温≥−25 ℃时,T取值为1,反之则取值为0;S表示土壤类型,S为草甸土、沼泽土及城区、岛屿、盐土时,取值为0,其余土壤类型为1;R表示岩石类型,R为水域时,取值为0,其余类型取值为1。由于大清河流域降水量符合油松降水量限制条件,故在叠加模型中暂不予考虑。

    • 油松适宜性计算主要采用最大熵模型(MaxEnt)。PHILLIPS等[16]基于最大熵原理编写了MaxEnt生态位模型软件,被广泛用来预测各物种的潜在适宜分布区[17-19]。本研究使用MaxEnt 3.4.1,将分辨率为30 m的24个地理环境变量作为模型的输入参数,地形变量如图2,土壤地质变量如图3,气象变量如图4

      图  2  地形变量示意图

      Figure 2.  Terrain variables

      图  3  土壤地质变量示意图

      Figure 3.  Soil geological variables

      图  4  气象变量示意图

      Figure 4.  Meteorological variables

    • 经MaxEnt软件计算,得到受试者工作曲线(ROC)图,该图反映模型灵敏度,可评价预测模型的准确程度,ROC曲线下面积(AUC)越大,评价模型越准确[20]。由图5可见:基于MaxEnt模型和地理环境变量构建的油松适宜性预测模型的训练数据集AUC为0.869,大于0.8,达到很准确水平;检验数据集AUC为0.817,大于0.8,达到很准确水平。综合模型准确性,说明油松适宜性预测模型有效性较高,可用于白洋淀大清河流域油松精准适宜性空间分布研究。

      图  5  ROC曲线图

      Figure 5.  ROC curve

      刀切法(Jackknife)[21-22]可预测各地理环境变量对油松适宜性空间分布影响的得分情况,可确定各个环境变量的重要程度。从图6可见:多年年平均积温、多年年平均最大日降水量、高程、多年年平均气温标准差依次是对油松预测贡献率最高的地理环境变量,为流域内油松精准适宜性预测提供了丰富信息。此外,当忽略岩石类型变量时,会最大程度上降低地理环境变量总得分,因此,它具有最多其他变量中没有的信息,也需要综合考虑。

      图  6  刀切法地理环境变量贡献图

      Figure 6.  Jackknife environmental variable contribution graph

      通过MaxEnt模型得到各地理环境变量与发生概率之间的反馈曲线,该曲线可反映单个地理环境变量值与发生概率之间的响应关系。从图7A可以看出:随着多年年平均积温的升高,分布概率呈先增大后减小的趋势,以分布概率0.5为界,积温适宜区间为2 800~3 800 ℃。由图7B看出:随着高程升高,分布概率先增大后减小,以分布概率0.5为界,高程适宜区间为500~1 200 m。由图7C看出:随着多年年平均最大日降水量升高,分布概率先增大后减小,以分布概率0.5为界,年平均最大日降水量适宜区间为70~110 mm。由图7D看出:随着多年年平均气温标准差升高,分布概率先增大后减小,以分布概率0.5为界,年平均气温标准差适宜区间为0.31~0.42 ℃。

      图  7  主导地理环境变量反馈曲线

      Figure 7.  Leading environmental variable feedback curve

    • 参考李培琳[23]、麻亚鸿[24]的适宜性分级标准,确定概率在0~0.05为不适宜区,0.05~0.30为低适宜区,0.30~0.50为中适宜区,0.50~1.00为高适宜区。白洋淀大清河流域油松高适宜区面积为1 924 km2,中适宜区为3 477 km2,低适宜区为7 621 km2,不适宜区为30 038 km2,中、高适宜区面积达到流域总面积的12.54%,中、高适宜区主要分布在大清河流域上游山区的西北部,在灵丘县、涞源县、涞水县、易县、房山区、阜平县等区县。油松精准适宜性空间分布及其细节放大图如图8,分辨率为30 m。表2为油松高适宜区面积大于20 km2的县,可见,高适宜区主要集中在涞源县、灵丘县、易县、涞水县、清苑县、涿鹿县,其中涞源县面积最大,为462.60 km2,占流域总面积的1%,这是由于这些县符合油松适宜条件的区域较多,在利用油松林进行绿化时可着重考虑。

      图  8  油松精准适宜性分布及细节放大图

      Figure 8.  Accurate suitability distribution and detailed enlarged drawing of P. tabulaeformis

      表 2  各县高适宜区面积及比例

      Table 2.  Area and proportion of high suitable area in each county

      适生区面积/km2占总面积比例/%
      涞源县462.601.00
      灵丘县402.650.87
      易县 366.420.80
      涞水县181.620.39
      清苑县74.380.16
      涿鹿县64.290.14
    • 油松是华北地区主要造林树种,由于早期造林时,没有考虑油松适宜性或仅根据气候区划种植了大面积的人工油松林,使得油松生长缓慢,或衰退现象明显[25-26]。近年来,针对油松适宜性的研究仅停留在小尺度,且数据和研究成果精度不高,难以为油松防护林营造提供科学支撑。

      本研究从油松生态习性的角度,得出影响油松适宜分布的主导地理环境变量有积温、最大日降水量、高程、年平均气温标准差,贡献率均大于20%,高适宜区分布在积温2 800~3 800 ℃、高程500~1 200 m、年平均最大日降水量70~110 mm、年平均气温标准差0.31~0.42 ℃的区域内。张春霞等[6]认为:油松中低海拔(900~1 100 m)立地质量为中等以上水平,与本研究高程最大值接近,最小值有差异,是由于其调查样方集中于900 m以上。魏安琪[27]认为:坡向是影响油松适宜性最大的因子,但在本研究中坡向不是主导因子,这是由于该研究仅考虑了地形因子未考虑气象因子,且本研究分析后认为坡向通过影响区域局部水热条件而对气象因子产生影响,这在气候因子中有所体现。

      在油松林管护过程中,针对不适宜导致的油松低效林,可间伐及间种其他适宜植物,通过逐渐演替由新植物种代替油松成为优势种;其次,在油松防护林营造时,需根据上述条件判断该地块是否适宜油松生长。本研究由于研究区范围较大,仅考虑了土壤类型对油松生长的影响,没有考虑土壤理化性质[28]、土壤厚度[21]等,在防护林建设中需综合考虑。此外,在以地理环境为角度进行油松适宜性研究时,未考虑林业经营管理、经济因素、其他人为因素等。在防护林建设中,整地方式、苗木品质、森林采伐等经营管理技术均对油松有较大影响;从经济角度看,造林成本、管护投入等因素也会造成一定影响[29-30],也需综合考虑。

参考文献 (30)

目录

    /

    返回文章
    返回