FPAR(x,t)?(NDVI(x,t)?NDVIi,min)(NDVIi,max?NDVIi,min)?(FPARmax?FPARmin)?FPARmin
式中：NDVIi,max和NDVIi,min分别对应第i种植被类型的NDVI最大和最小值。FPAR与比值植被指数（SR）也存在着较好的线性关系，可由以下公式表示：
FPAR(x,t)?(SR(x,t)?SRi,min)(SRi,max?SRi,min)?(FPARmax?FPARmin)?FPARmin
式中，FPARmin和FPARmax的取值与植被类型无关，分别为0.001和0.95； SRi,max和SRi,min分别对应第i种植被类型NDVI的95%和5%下侧百分位数， SR(x,t)由以下公式表示：
SR(x,t)?1?NDVI(x,t)
1?NDVI(x,t)通过对FPAR-NDVI和FPAR-SR所估算结果的比较发现，由NDVI所估算的FPAR比实测值高，而由SR所估算的FPAR则低于实测值，但其误差小于直接由NDVI所估算的结果，因此我们可以将二者结合起来，取其加权平均或平均值作为估算FPAR的估算值：
FPAR(x,t)??FPARNDVI?(1??)FPARSR
光能利用率的估算：光能利用率是在一定时期单位面积上生产的干物质中所包含的化学潜能与同一时间投射到该面积上的光合有效辐射能之比。环境因子如气温、土壤水分状况以及大气水汽压差等会通过影响植物的光合能力从而调节植被的NPP。
?(x,t)?T?1(x,t)?T?2(x,t)?W?(x,t)??max
式中：Tε1(x,t)和Tε2(x,t)表示低温和高温对光能利用率的胁迫作用；Wε(x,t)为水分胁迫影响系数，反映水分条件的影响；εmax是
— 38 —
理想条件下的最大光能利用率（g C/MJ）。
温度胁迫因子的估算：Tε1(x, t)的估算：其反映在低温和高温时植物内在的生化作用对光合的限制而降低第一性生产力。
Tε1(x,t)?0.8?0.02?Topt(x)?0.0005?[Topt(x)]2
式中：Topt(x)为植物生长的最适温度，定义为某一区域一年内NDVI值达到最高时的当月平均气温（℃）；当某一月平均温度小于或等于-10℃时，其值取0。 Tε2(x, t)的估算：表示环境温度从最适温度Topt(x)向高温或低温变化时植物光能利用率逐渐变小的趋势，这是因为低温和高温时高的呼吸消耗必将会降低光能利用率，生长在偏离最适温度的条件下，其光能利用率也一定会降低。
T?2(x,t)?1.184/{1?exp[0.2?(Topt(x)?10?T(x,t))]}?1/{1?exp[0.3?(?Topt(x)?10?T(x,t))]}
当某一月平均温度T(x, t)比最适温度Topt(x)高10℃或低13℃时，该月的Tε2(x, t)值等于月平均温度T(x, t)为最适温度Topt(x)时Tε2(x, t)值的一半。
水分胁迫因子的估算：水分胁迫影响系数Wε(x,t)反映了植物所能利用的有效水分条件对光能利用率的影响，随着环境中有效水分的增加，Wε(x,t)逐渐增大,它的取值范围为0.5(在极端干旱条件下)到1(非常湿润条件下)。
W?(x,t)?0.5?0.5?EET(x,t)/EPT(x,t)
式中：EET为区域实际蒸散量（mm）； EPT为区域潜在蒸散量（mm） 最大光能利用率的确定：月最大光能利用率εmax的取值因不同的植被类型而有所不同，在CASA模型中全球植被的最大光能利用率为
— 39 —
0.389 g C〃MJ。
全国尺度上利用CASA模型计算NPP时主要采用MODIS250米每16天合成的NDVI数据产品和反射率数据产品，以及MODIS1千米每8天合成的陆地表面温度(LST)数据产品，并利用气象站点太阳辐射量数据作为遥感估测模型的补充，通过克里金插值得到250米分辨率日照时数数据，以实现NPP估算由点向面的推演和扩展。在数据条件允许的情况下，可使用更高时空分辨率数据进行计算。
-1
— 40 —
附录B
修正通用水土流失方程（USLE）的土壤保持功能算法
模型结构：Ac = Ap - A r ＝ R× K× L×S×（1- C） 式中： Ac为土壤保持量，Ap为潜在土壤侵蚀量，Ar为实际土壤侵蚀量，R为降水因子，K为土壤侵蚀因子，L、S为地形因子，C为植被覆盖因子。
R－降水因子
?122?R??（）/P??Pi??i?1??
式中：Pi为月均降水量，P为年均降水量, α= 0.3589，β= 1.9462。
K－土壤侵蚀因子
K = fcsand ×fcl-si ×forgc ×fhisand fcsand = 0.2 + 0.3 exp [-0.0256 ms (1-msilt/100)] fcl-si = [msilt/(mc + msilt)]0.3 forgc = 1- 0.25orgC / [orgC + exp(3.72 - 2.95 orgC)] fhisand = 1-0.7(1-ms/100)/{(1-ms/100)+exp[-5.51+22.9(1-ms/100)]} 式中：ms为土壤粗砂含量，msilt为土壤粉砂含量，mc为土壤粘粒含量，orgC为有机碳含量。
— 41 —
L、S—地形因子 L = (λ/22.13)m S = 10.8sinθ + 0.03 θ＜5° 16.8sinθ — 0.5 5°≤ θ ＜10° 式中：
— 42 —
21.91sinθ — 0.96 为坡长（m），m为坡长指数，θθ ≥ 10°为坡度（°）。
λ