生态保护红线区生态系统服务多维驱动力的定量分离

左丽媛, 姜远, 高江波, 杜馥君, 张怡博

Quantitative separation of multi-dimensional driving forces of ecosystem services in the ecological conservation red line area

ZUO Liyuan, JIANG Yuan, GAO Jiangbo, DU Fujun, ZHANG Yibo

表1 生态系统服务模型

Tab. 1 Models used to quantify ecosystem services

生态系统服务	模型	计算公式
土壤保持	RUSLE模型^[28]	$A = R \times K \times L S \times (1 - C \times P)$ A为土壤保持量(t hm^-2 a^-1);R为降水侵蚀力因子(MJ mm hm^-2 h^-1 a^-1),考虑到研究区降水的季节性差异,本文采用月降水数据计算R因子,进而得到年尺度数据^[29];K为土壤侵蚀力因子(t hm² h hm^-2 MJ^-1 mm^-1),本文基于土壤属性数据,运用侵蚀力影响计算器(EPIC)计算获得^[30];LS为坡长坡度因子,为充分表达研究区山地海拔的复杂性特征^[31],本文选取空间分辨率为9 m的DEM计算LS因子;C是植被覆盖与管理因子,本文以NDVI为基础数据进行计算^[32];P是水土保持措施因子,本文参考华北平原的研究,根据土地利用类型进行赋值^[33];LS、C、P为无量纲。
产水量	InVEST模型产水量模块^[34]	$Y (x) = (1 - \frac{A E T (x)}{P (x)}) \times P (x)$ Y(x)为研究区栅格单元x的年产水量(mm);AET(x)为栅格单元x的年实际蒸散发量(mm);P(x)为栅格单元x的年降水量。
植被固碳	CASA模型^[35]	$N P P (x, t) = A P A R (x, t) \times ε (x, t)$ NPP(x, t)为空间位置x在时间t内积累的植被净初级生产力(g C m^-2);APAR(x, t)和ε(x, t)分别为空间位置x在时间t内植被吸收的光合有效辐射(MJ m^-2)和实际光能利用率(g C MJ^-1)。