城市生态--生产--生活空间功能定量识别与分析

doi:10.11821/dlxb201601004

摘要/Abstract

摘要：

土地利用多功能性识别是城市用地组织,协调与配置的基础信息源,是判定城市用地内在功能形态,功能组合模式和功能之间动态权衡的关键,具有重要的理论和实践意义,但长期以来并未构建一套可行的识别方法体系.本文从土地功能,生态系统服务和景观功能综合的视角构建城市生态--生产--生活空间功能分类体系,并以生态系统服务价值评估为基础系统整合空间功能价值量核算函数群,通过纵横对比的方法确定空间功能主导类型.研究区实证分析表明,城市生态--生产--生活功能分类体系较好反映了不同地类的功能类型;空间主导功能的识别也与不同地类的功能匹配;同时也发现三生空间的整体毗邻性较低,空间功能的互补和融合性较差的问题;三生空间功能存在一定的空间集聚性.

关键词: 三生空间, 城市, 土地功能, 生态系统功能, 景观功能, 功能识别, 价值评估

Abstract:

The identifying of land multifunctionality is a basic tool for organization, coordination and configuration of urban land, and is a key criterion for urban land functions forms, composite pattern and dynamic tradeoffs. This topic is of important theoretical and practical significance. An available identification system of urban land multifunctionality, however, had not been built for a long time. This paper develops a function classification system for urban ecological-production-living spaces from an integrated perspective of land function, ecosystem services and landscape function. We integrated a value function group of space function based on ecosystem services valuation. A comparison method of vertical and cross direction is proposed to identify dominant function type of urban land. The empirical results show that in the study area, function classification system of urban ecological-production-living spaces can reflect differentiated function types of different land use classes. The identified result of dominant function for urban space is matched with the functions of different land types. Meanwhile, we also found that there are some problems on urban land use, such as the low adjacency degree with different spaces, and poor complementarity with different space functions. The result indicates that the spatial distribution of urban ecological-production-living spaces is congregated in the study area.

Key words: ecological-production-living spaces, urban area, land function, ecosystem services, landscape function, function identifying, value assessment

李广东, 方创琳. 城市生态--生产--生活空间功能定量识别与分析[J]. 地理学报, 2016, 71(1): 49-65.

LI Guangdong,FANG Chuanglin. Quantitative function identification and analysis of urban ecological-production-living spaces[J]. Acta Geographica Sinica, 2016, 71(1): 49-65.

图/表 11

表1

城市生态--生产--生活空间功能分类体系

功能形式	一级功能	二级功能	解释	三级潜在表征指标
生态功能	调节功能	E1气体调节	大气化学成份调节	单位面积大气污染物净化能力
		E2气候调节	碳汇,温度,降水及其它由生物媒介的全球及地区性的气候调节	单位面积碳汇量/碳储量
		E3水调节	地表水体径流调节与涵养	渗透量/渗透比率(径流量/表层面积);土壤水蓄储能力(mm/m);泛滥平原水蓄储能力(mm/m)
		E4水和废物净化	水源净化沉积物,农药,病原细菌和病原体,过滤和分解水中和土壤中的有机物质	输沙量(g/l);固体溶解总量(mg/l),分解总量和种类(n/hm²);分解率(kg/hm²·a)
		E5缓和突发事件	洪水,火灾,飓风,地震和疾病控制等	洪峰强度,地震强度,飓风强度及损失量,控制的疾病等
		E6授粉	生态系统变化影响到授粉的分布,丰度和效率	植物种类,分布以及传粉者的可获得性
		E7土壤保持	土壤保持,土壤侵蚀和滑坡控制能力	植被覆盖率(%);水侵和风侵土壤流失量(kg/hm²·a);通用土壤流失方程评估滑坡频率系数(n/hm²·a)
		E8养分循环	养分循环能力(氮,磷等)	水质指标:氮(mg/l),磷(mg/l);营养流失量(kg/hm²·a);导电系数(μs/cm);固体溶解量(mg/l);营养流通率,氮,磷(y^-1)等
		E9初级生产	植物初级生产能力	NPP(净初级生产力)
生产功能	生存与健康物质供给功能	P1食物供给	持续生产或收获粮食,水果,野生浆果,蘑菇,坚果,牲畜,半家养动物,猎物,鱼和其它水生资源等	持续(有机的)获取粮食产品(t,hm²);持续(有机的)获取牲畜(t,hm²);持续获取渔类产品(t,hm²);野生动物可作为食物的种类;持续获取野生动植物产品(t,hm²)
		P2淡水供给	淡水提供	整体的淡水资源量(百万m³)
		P3药物供给	持续生产和获得自然药物产品,生物化学或药物产品成分	由自然药物衍生的药品种类;运用自然化合物的药物种类
		P4基因资源	当地和地方性的品种保护,品种基因库的维持	辛普森(Simpson)指数;香农-威纳(Shannon-wiener)指数;粮食种类;牲畜品种;渔类品种
	原材料生产功能	P5木材供给	木材产出	自然供给和持续产出的工业圆材量(百万m³);由自然森林获取的纸浆和纸产品量;
		P6纤维供给	纤维制品材料提供	持续的(有机的)棉花产品供给量(t,hm²);
		P7装饰资源	持续生产和获得装饰和观赏用的野生植物,木制品和贝壳等	用于工艺品制作的自然材料种类; 持续获取观赏用园艺植物的总量
	能源矿产生产功能	P8能源生产	石油,天然气,生物质能等	生物质能(百万t油当量);石油和天然气产量(百万t油当量)
	能源矿产生产功能	P9矿产生产	矿物产出	矿产,黄金等
	间接生产功能	P10商品与服务产品生产	第二,三产业生产	批发零售用地,住宿餐饮用地,商务金融用地,其它商服用地,工业用地以及商用仓储用地等
生活功能	空间承载与避难功能	L1居住承载	住房与附属设施用地	住宅用地价格(基准价格)
		L2交通承载	水陆空交通线路及设施用地	替代价格(基准地价)
		L3公共服务承载	公共服务设施用地,避难空间等	替代价格(基准地价)
	物质生活保障功能	L4基本生活保障	维持基本生活保障	征地补偿标准,最低生活保障标准
	物质生活保障功能	L5就业保障	提供就业机会	直接提供和保障的就业人数
	精神生活保障功能	L6科学和教育	为科学研究和教育提供对象	科研数量;教育课程数量,人数
		L7休闲	提供休闲旅游空间	最大可持续承载的旅游人数,营业收入
		L8文化和艺术	提供艺术欣赏对象和设计灵感	以空间作为灵感或基础产出的书籍,印刷品数量
		L9美学	提供美学基础,景观美学质量等	美学价值,例如靠近自然保护区的户数
		L10精神和历史	景观,宗教,文化,图腾,历史文物等	精神或宗教吸引的人数

表1

图1

表3

生态空间功能定量测度模型

功能类型	公式	解释	来源
气体调节功能	$V i, O 2 = V i, O 2_NPV × r × (1 + r) t / [(1 + r) t - 1 V O 2_NPV = (V 1 O 2 + V 2 O 2) 2 × 1 (1 + r) t V 1 O 2 = Q O 2 × S 1 O 2 V 2 O 2 = Q O 2 × S 2 O 2 Q O 2 = 1.2 × ∑ NP P i$	$V i, O 2$ 为第i种土地利用的释放氧气单位年价值量; $V i, O 2_NPV$ 为释放氧气的净现值;V₁O₂为采用造林法测算的释放氧气价值;V₂O₂为采用工业制氧法测算的释放氧气价值;r为社会折旧率,t可以视为土地利用用途不变的时间,此处r为0.07,t为10;QO₂为氧气释放物质量,一般通过NPP(净初级生产力)来计算;1.2是生产1 g植物干物质能释放氧气的数量;S₁O₂为氧气单位释放量的造林成本(元/t),由1989年造林成本为260.9元/t C (参见中国林业统计年鉴,1990年),可以推算出2009年的价格为646.29元/t C;S₂O₂为氧气单位释放量的工业制氧成本(元/a),1991年工业制氧成本为0.4元/kg (参见中国统计年鉴,1992年),可以推算出2009年价格为0.99元/kg;NPP_i是第i种植被单位面积净初级生产力.	[32-33]
气候调节功能	$V i, C = V i, C_NPV × r × (1 + i) t / [(1 + r) t - 1 V i, C_NPV = (V i, C_market + V i, C_social) × 1 (1 + r) t V i, C_market = (Q A, i + Q B, i + Q S, i) × NP P i NP P mean × C c V i, C_social = (Q A, i + Q B, i + Q S, i) × Cost (P s)$	V_i_,_C为单位碳固定和储存的年价值;V_i_,_C_{_NPV}为单位碳固定和储存的净现值;r为社会折旧率;t可以视为土地利用用途不变的时间;V_i_,_C_{_markrt}为市场价值;V_i_,_C_{_social}为社会价值;Q_A_,_i是第i类土地利用方式单位面积地表碳固定和储存量;Q_B_,_i为单位面积地下碳固定和储存量;Q_S_,_i为单位面积土壤有机碳固定和储存量;C_c为瑞典碳税率,40.5美元/t,换算为人民币为251.72元/t (2009年);Cost(P_s)为碳固定和储存的社会成本.Q_A_,_i,Q_B_,_i和Q_S_,_i的单位面积数量值参考Leh等的研究;Cost(P_s)参考InVEST的设定方法,100美元/t(2009年)左右,换算为人民币为621元/t.	[34-35]
水调节功能	$V i, a = V w × r × (1 + r) t / [(1 + r) t - 1 V w = Q t × P × 1 (1 + r) t Q t = Q c + Q l + Q s Q c = ε × P y Q l = G i × L i Q s = h i × f i$	V_i_,_a为水源涵养单位年价值(元);V_w为水源涵养总价值(元);r为社会贴现率,t为项目运行年限(a),折现率一般为0.07,水库使用年限按20a计算;P为单位水库库容造价,根据1993-1999年《中国水利年鉴》平均水库库容造价为2.17元/t,推算得到单位库容造价为7.14元/t (2009年);Q_t为水源涵养总量;Q_c为树冠截水量; $ε$ 为树冠截水率;P_y为年降水量;Q_l为枯枝落叶层持水量;G_i为枯枝落叶层干重;L_i为枯枝落叶层饱和持水率;Q_s为土壤层持水量;h_i为第i种土壤厚度;f_i为第i种土壤的粗孔隙率.	[36-37]
水和废物净化功能	$V x = Cost (p) × Retention × ∑ t = 0 T - 1 1 (1 + r) t Retention = AL V x × filtration AL V x = λ x λ w ? × po l x λ a = log (∑ Y)$	V_x为水质和废物净化年价值量;Cost(p)为污染物处理成本,以最终污染物过滤量为单位,目前全国单位面积最终污染物处理成本为300元/t;Retention为污染物过滤量;r为贴现率,一般为0.07;T为污水处理厂使用年限,此处设为20a;ALV_x为调整后的污染物载荷量;filtration为污染物过滤能力;λ_x为径流系数 $λ w ?$ 为平均径流系数,通过两者计算水文敏感分值;pol_x为污染物产出系数;λ_α为地表径流系数计算的一般公式,其中Y为产水量.相关参数设定参考Leh等的研究附录.	[19, 34]
缓和突发事件功能	$V i, e = V mean, i × f (v_inde x i, R i)$	V_i,e为第i种土地利用缓和突发事件年价值;V_mean为TEEB全球范围内计算得到的第i类土地利用方式缓和突发事件年均值;v_index为第i类土地利用方式的脆弱性指数,用于说明第i类土地利用方式的整体状况,一般取值为0~1;R为第i类土地利用方式的恢复力,用于表征缓和突发事件的效率水平,一般取值为0~1.
授粉功能	$V i_hb = ∑ i = 1 n (V i × D i × P i) D i = (Y i 0 - Y it) / Y i 0$	V_{i_hb}为每年农作物授粉产生的经济价值;V_i为农作物年产量;D_i为农作物昆虫授粉依存度;P_i为农作物产品价格;D_i为授粉依存度;Y_i₀为开放授粉区农作物产量;Y_ic为无授粉区农作物产量.	[38]
土壤保持功能	$S t = ? M S r - M S f MS (p t) = β l × μ j × σ k × MS (p s) M S f = ∑ l k MS (p l) / k M S r = ∑ j h MS (p j) / h$ $V i, S = (V S 1 + V S 2 + V S 3) × r × (1 + r) t / [(1 + r) t - 1 V S 1 = R i × (1 / n i) × P i V S 2 = Cos t i / (104 × ρ i × h i) V S 3 = 0.24 × C / ρ i$	S_t为土壤流失减少量(m³/a);MS_r和MS_f分别为无植被覆盖和有植被覆盖的土壤侵蚀模数;MS(p_t)是第i类植被--土壤--坡度复合的土壤侵蚀模数的计算公式;β_l为植被参数(l =1, 2, 3, et al.6);μ_j为土壤参数(j =1, 2, 3, et al.5);σ_k为坡度参数(k =1, 2, 3);MS(p_S)为标准植被--土壤--坡度复合状态下的侵蚀模数;k和h分别为有植被和无植被总数;V_i_,_s为第i种土地利用的单位年价值(r和t的设置同前式);V_S₁为保持土壤养分的单位价值;R_i为单位质量土壤中养分(碱解氮,速效磷,速效钾)的平均含量(t);n_i为土壤中碱解氮,速效磷,速效钾分别在磷酸二胺和氯化钾中的含量,分别为14%,15%和60%;P_i为化肥(磷酸二胺,氯化钾)价格(元/t),2009年价格分别为3500元/t,2900元/t;V_S₂为减少废弃土地的单位价值;ρ_i为土壤容重(t/m³);h_i为土壤厚度(m);Cost_i为第i种土地利用方式单位面积的机会成本,或者第i种土地利用方式的年均收益(元/(hm²·a)).V_S₃为减少泥沙淤积的单位价值;C为单位库容造价(元/m³),2009年单位库容造价为7.02元/m³.	[32, 39]
养分循环功能	$V i, n = NP P i × (N i C 1 R 1 + P i C 2 R 2 + K C i C 2 R 3 + M i C 3)$	V_i_,_n为i种土地利用类型的单位面积养分循环年价值;NPP_i为年均净初级生产力;N_i,P_i,K_i,M_i分别为不同土地利用方式下土壤中氮磷钾和有机质的含量;C₁,C₂,C₃分别为磷酸二胺,氯化钾化肥,有机质价格(元/t),2009年价格分别为3500元/t,2900元/t和500元/t;R₁,R₂,R₃分别为磷酸二胺含氮量,磷酸二胺含磷量,氯化钾含钾量(%),分别为14%,15%和60%.	[40]
初级生产功能	$NPP = RDI r R n (r 2 + R n 2 + r R n) (R n + r) (R n 2 + r 2) × exp [- (9.87 + 6.25 RDI) 0.5] × 100 R n = RDI × r × L × 2.38 × 10 - 4 RDI = 0.629 + 0.237 PER - 0.00313 PE R 2 PER = 58.93 BT TAP$ $V i_NPP = NP P i × P S × 1.474 × 10 - 6$	NPP为净初级生产力(gDWm^-2yr^-1);r为年均降水量(mm),RDI为辐射干燥指数;R_n为地表截获的净辐射量(Jcm^-1yr^-1);R_n可以通过PER计算得出,其中PER为蒸散率,通过BT年均温(0℃~30℃)和TAP年均总降水量计算,通过月数据算得;L为潜伏热(2503 Jg^-1);V_{i_NPP}为第i种土地利用类型的单位NPP年价值;NPP_i为净初级生产力;P_s为标准煤价格(元/t),2009年价格为元877元/t (1 g碳相当于2.2 g的有机物质,1 g的有机物质在燃烧释放的能量上相当于0.679 g的标准煤,而1 g碳相当于1.474 g (2.2×0.679)的标准煤).	[41-43]

表3

表4

生产空间功能定量测度模型

功能类型	公式	解释	来源
食物,纤维供给以及能源和矿产生产功能	$V i, a = ∑ Y i, a × P i, a - Cos t i, a$ $V i, f = ∑ Y i, f × P i, f - Cos t i, f$ $V i, e = ∑ Y i, e × P i, e - Cos t i, e$ $V i, e = ∑ Y i, e × P i, e - Cos t i, e$	V_i_,_a为第i种土地利用食物供给的年价值;Y_i_,_a为第i类可收获食物的单位面积年均产量(kg/hm²);P_i_,_a为第i类食物的当年市场平均价格(元/kg);Cost_i_,_a为单位面积成本;V_i_,_f为第i种土地利用纤维供给的年价值;Y_i_,_f为第i类可收获纤维的单位面积年均产量(kg/hm²);P_i_,_f为第i类纤维的当年市场平均价格(元/kg);Cost_i_,_f为单位面积成本;V_i_,_e为第i种土地利用方式能源供给的年价值;Y_i_,_e为第i类土地利用的能源产出功能量;P_i_,e为第i类能源的当年市场平均价格;Cost_i_,_e为单位面积能源的获取成本.V_i_,_e为第i种土地利用方式下矿产供给的年价值;Y_i_,_e为第i类土地利用方式下的矿产产出量;P_i_,_e为第i类矿产的当年市场平均价格;Cost_i_,_e为单位面积矿产的获取成本.
淡水供给功能	$V i, w = Y xj × WP; Y xj = (1 - AE T xj P x) × P x AE T xj P x = 1 + ω R xj 1 + ω R xj + 1 R xj; R xj = k xj × E T 0 P x ω x = Z AW C x P x$	Y_xj为产水量;WP为单位面积水价,2009年均价为2.5元/t;Y_xj为x土地覆被类型j年产水量;AET_xj为x土地覆被类型j年蒸散发量;P_x为年均降水量;ω为修正植被年可利用水量与预期降水量的比值;R_xj为x土地覆被类型j的Budyko干燥指数(实际蒸散与降水的比值),无量纲;k_xj为植被蒸散系数;ET₀为参考作物蒸散;Z为Zhang系数,此处为3.33;AWC_x为土壤有效含水量;相关参数设定参考Leh等的设定方法.	[34]
药物供给功能	$V i, m = V mean, i × f (adjustment inde x i)$	V_i_,_m为第i类土地资源产出药物资源的年价值;V_mean_,_i为TEEB统计的全球范围内的平均值;f(adjustment index_i )是根据当地调查了解到的药物资源种类,产量等而设定的调整系数.
基因资源供给功能	$V = (H × C 1 + P × 1.26 × 1025) × C 2 / A H = ∑ i = 1 s N i / Nlo g 2 N i / N P = (m / M) × (t / 12) × (s / S)$	H为Shannao-Wearer指数;C₁为bits转化率,为5.78×10¹³,C₂为能值--经济值转化率,为8.67×10¹²;P为珍稀物种支持率;A为调查总面积;N_i为第i种物种的个体数;N是S个物种的总个体数;m为生态系统中珍稀物种个体数;M为地球上珍稀物种总数;t为珍稀物种一年内在该生态系统中生活的时间(月);s为该系统面积;S为t时间内物种个数的实际活动面积.	[44]
木材供给功能	$V H i = Perc_har v i × Perc_avai l i × Pric e i × Accum_mas s i - Harv_cos t i$	VH_i为单位面积可收获木材的净值;Perc_harv_i为可收获的木材比例;Perc_harv_i为综合出材率;Price_i为木材市场价格,2009年原木价格为2800元/m³;Accum_mass_i为单位面积木材蓄积量;Harv_cost_i为单位面积木材年收获成本.	[34]
装饰资源供给功能	$V i, o = V mean, i × f (adjustment inde x i)$	V_i,o为第i类土地资源产出装饰资源的年价值;V_mean_,_i为TEEB统计的全球范围内的平均值;f(adjustment index_i )是根据当地调查了解到的装饰资源的种类,产量等而设定的调整系数.
商品与服务产品生产功能	$V i, l and = V i, b_price × (1 + ∑ K i, factors) × Y i, year × D i, data × R i, plot ratio$ $V i, l and = (E i, a + E i, d + T + R i, 1 + R i, 2 + R i, 3) × S × K$	V_i_,_land为待估宗地价格;V_i_,_{b_price}为待估宗地所处区的基准地价;K_i_,_factors为各影响因素修正系数之和,包括区域因素(城市繁华程度及区域在城市中的位置交通条件,公共及基础设施水平,区域环境条件,土地使用限制和自然条件等)和个别因素(与宗地直接有关的自然条件,市政设施条件,宗地面积,形状,长度,宽度,地质条件,使用限制和宗地位置等);Y_i_,_year为年期修正系数,一般为 $Y i, year = [1 - 1 / (1 + r) m] / 1 - 1 / 1 + r) n]$ ,m为待估宗地的使用年期;n为基准地价设定的使用年期;D_i_,_data为期日修正系数, $D i, data = Q i / Q 0$ ,Q_i为估价期日地价指数,Q₀为基准地价估价期日地价指数;R_{i,plot ratio}为容积率修正系数, $R i, plot ratio = R i / R 0$ ,R_i为待估宗地容积率,R₀为基准地价设定的容积率.E_i_,_a为土地取得费;E_i_,_d为土地开发费;T为税费;R_i_,1为利息;R_i_,2为利润;R_i_,3为土地增值;S为个别修正系数;K为年期修正系数.

表4

表5

生活空间功能定量测度模型

功能类型	公式	解释	来源
居住承载功能	$V i, l and = V i, b_price × (1 + ∑ K i, factors) × Y i, year × D i, data × R i, plot ratio$ $V i, l and = (Cos t i + R i) × S × K$ (乡村宅基地价格核算)	V_i_,_land为待估宗地价格;V_i_,_{b_price}为待估宗地所处区的基准地价;K_i_,_factors为各影响因素修正系数之和,包括区域因素和个别因素;Y_i_,_year为年期修正系数;m为待估宗地的使用年期;n为基准地价设定的使用年期;D_i_,_data为期日修正系数;R_i_,_{plot ratio}为容积率修正系数;Cost_i为土地成本价格,成本无法核算的可以运用当地征地补偿标准(每年)替代;R_i为土地增值;S为个别修正系数;K为年期修正系数.
交通与公共服务承载功能	$V i, l and = P ? i × S × D$	V_i_,_land为待估用地价格; $P ?$ 为待估用地所处区周边用地的平均价格;S为个别修正系数;D为区域修正系数.
基本生活保障功能	$V i = P i × L / K$	V_i为基本生活保障功能价值;P_i为第i种土地利用方式的人口承载力;L为研究区相应年份的城镇最低生活保障水平(元/人·a);K为相应年份城市居民人均可支配收入与农民人均纯收入的比例,即城农收入比.	[45]
就业保障功能	$V w = (P 0 - P 1) × B / K$	V_w为失业保险功能价值;P₀为人口承载力,中国土地资源在中等投入水平下人口承载力为3.68~5.51人/hm²,此处取5人/hm²;P₁为土地资源的最小承载力,此处取3.68人/hm².B为研究区相应年份的城镇失业保险金标准(元/人·a).K为相应年份城市居民人均可支配收入与农民人均纯收入的比例,即城农收入比.	[45]
休闲功能	$P travel = T cap × I i × T i - 1 × l / A i l = L × (1 + a e - bt) - 1 t = E n - 1 - 3$	P_travel为研究区休闲价值;T_cap为研究区总体旅游环境容量;I_i为研究区第i年的旅游收入,T_i为第i年的游客总量;l为社会发展阶段系数;L为极富阶段的支付意愿,取值为1;t为时间变量,表示社会发展阶段;a,b为常数,取值为1;e为自然对数;En为恩格尔系数,A_i为旅游区的总面积(hm²).	[33]
科学和教育,文化和艺术,美学功能,精神和历史功能	$V i, SE = V mean, i × f (adjustment ind e x i) V i, EA = V mean, i × f (adjustment ind e x i) V i, AE S = V mean, i × f (adjustment ind e x i) V i, SH = V mean, i × f (adjustment ind e x i)$	V_i_,_c为第i类土地资源具有的科学和教育,文化和艺术,美学功能,精神和历史的年价值;V_mean_,_i为TEEB统计的全球范围内第i类服务的平均值;f(adjustment index_i)是根据当地调查了解到的级别,功能量,文化艺术价值等信息而设定的调整系数.

表5

图2

表6

图3

图4

图5

图6

参考文献 45

[1]	Verburg P H, Van De Steeg J, Veldkamp A, et al. From land cover change to land function dynamics: A major challenge to improve land characterization. Journal of Environmental Management, 2009, 90(3): 1327-1335.
[2]	European Commission.Safeguarding the multifunctional role of agriculture: Which instruments? Brussels: DG Agriculture-European Commission, 1999.
[3]	Abler D.Multifunctionality, Agricultural policy and environmental policy. Agricultural and Resource Economics Review, 2004, 33(1): 8-17.
[4]	Andersen P S, Vejre H, Dalgaard T, et al.An indicator-based method for quantifying farm multifunctionality. Ecological Indicators, 2013, 25: 166-179.
[5]	Xie G, Zhen L, Jochen K H, et al.Assessing the multifunctionalities of land use in China. Journal of Resources and Ecology, 2010, 1(4): 311-318.
[6]	Gómez-Limón J A, Vera-Toscano E, Rico-González M. Measuring individual preferences for rural multifunctionality: The importance of demographic and residential heterogeneity. Journal of Agricultural Economics, 2012, 63(1): 1-24.
[7]	Mittenzwei K, Fjellstad W, Dramstad W, et al.Opportunities and limitations in assessing the multifunctionality of agriculture within the CAPRI model. Ecological Indicators, 2007, 7(4): 827-838.
[8]	Peng Jian, Liu Zhicong, Liu Yanxu.Research progress on assessing multi-functionalty of agriculture. Chinese Journal of Agricultural Resources and Regional Planning, 2014, 35(6): 1-8.
	[彭建, 刘志聪, 刘焱序. 农业多功能性评价研究进展. 中国农业资源与区划, 2014, 35(6): 1-8.]
[9]	Wiggering H, Müller K, Werner A, et al.The concept of multifunctionality in sustainable land development//Katharina Helming, Wiggering H. Sustainable Development of Multifunctional Landscapes. New York, Berlin, Heidelberg: Springer-Verlag, 2003.
[10]	Paracchini M L, Pacini C, Jones M L M, et al. An aggregation framework to link indicators associated with multifunctional land use to the stakeholder evaluation of policy options. Ecological Indicators, 2011, 11(1): 71-80.
[11]	Wiggering H, Dalchow C, Glemnitz M, et al.Indicators for multifunctional land use: Linking socio-economic requirements with landscape potentials. Ecological Indicators, 2006, 6(1): 238-249.
[12]	Zhang Xiaoping, Zhu Daolin, Xu Zuxue.Assessment on multi-functionality of land use in Tibet. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(6): 185-194.
	[张晓平, 朱道林, 许祖学. 西藏土地利用多功能性评价. 农业工程学报, 2014, 30(6): 185-194.]
[13]	Lv Yihe, Ma Zhimin, Fu Bojie, et al.Diversity of ecosystem services and landscape multi-functionality: From scientific concepts to integrative assessment. Acta Ecologica Sinica, 2013, 33(4): 1153-1159.
	[吕一河, 马志敏, 傅伯杰, 等. 生态系统服务多样性与景观多功能性: 从科学理念到综合评估. 生态学报, 2013, 33(4): 1153-1159.]
[14]	Zhen Lin, Wei Yunjie, Xie Gaodi, et al.Regional analysis of dynamic land use functions in China. Acta Ecologica Sinica, 2010, 30(24): 6749-6761.
	[甄霖, 魏云洁, 谢高地, 等. 中国土地利用多功能性动态的区域分析. 生态学报, 2010, 30(24): 6749-6761.]
[15]	Zhen Lin, Cao Shuyan, Wei Yunjie, et al.Land use functions. Resources Science, 2009, 31(4): 544-551.
	[甄霖, 曹淑艳, 魏云洁, 等. 土地空间多功能利用. 资源科学, 2009, 31(4): 544-551.]
[16]	Millennium Ecosystem Assessment.Ecosystems and Human Well-being Synthesis. Washington, DC: Island Press, 2005.
[17]	Costanza R, D'arge R. The value of the world's ecosystem services and natural capital. Nature, 1997, 387: 253-260.
[18]	Egoh B, Reyers B, Rouget M, et al.Mapping ecosystem services for planning and management. Agriculture Ecosystems & Environment, 2008, 127(1/2): 135-140.
[19]	Leh M D K, Matlock M D, Cummings E C, et al. Quantifying and mapping multiple ecosystem services change in West Africa. Agriculture, Ecosystems & Environment, 2013, 165: 6-18.
[20]	Gret-Regamey A, Bebi P, Bishop I D, et al.Linking GIS-based models to value ecosystem services in an Alpine region. Journal of Environmental Management, 2008, 89(3): 197-208.
[21]	Chen N, Li H, Wang L.A GIS-based approach for mapping direct use value of ecosystem services at a county scale: Management implications. Ecological Economics, 2009, 68(11): 2768-2776.
[22]	de Groot R. Function-analysis and valuation as a tool to assess land use conflicts in planning for sustainable, multi-functional landscapes. Landscape and Urban Planning, 2006, 75(3/4): 175-186.
[23]	Gulickx M M C, Verburg P H, Stoorvogel J J, et al. Mapping landscape services: A case study in a multifunctional rural landscape in the Netherlands. Ecological Indicators, 2013, 24: 273-283.
[24]	Lovell S T, Desantis S R, Nathan C A, et al.Integrating agroecology and landscape multifunctionality in Vermont: An evolving framework to evaluate the design of agroecosystems. Agricultural Systems, 2010, 103(5): 327-341.
[25]	O'farrell P J, Anderson P M L. Sustainable multifunctional landscapes: A review to implementation. Current Opinion in Environmental Sustainability, 2010, 2(1-2): 59-65.
[26]	Peng Jian, Lü Huiling, Liu Yanxu, et al.International research progress and perspectives on multifunctional landscape. Advances in Earth Science, 2015, 30(4):465-476.
	[彭建, 吕慧玲, 刘焱序, 等. 国内外多功能景观研究进展与展望. 地球科学进展, 2015, 30(4): 465-476.]
[27]	Willemen L, Hein L, Van Mensvoort M E F, et al. Space for people, plants, and livestock? Quantifying interactions among multiple landscape functions in a Dutch rural region. Ecological Indicators, 2010, 10(1): 62-73.
[28]	de Groot R, Wilson M A, Boumans R M J. A typology for the classification, description and valuation of ecosystem functions, goods and services. Ecological Economics, 2002, 41(3): 393-408.
[29]	de Groot R. Functions of Nature: Evaluation of Nature in Environmental Planning, Management and Decision Making. Groningen: Wolters-Noordhoff, 1992.
[30]	Willemen L.Spatial characterization of landscape functions. Landscape and Urban Planning, 2008, 88(1): 34-43.
[31]	de Groot R, Brander L, Van Der Ploeg S, et al. Global estimates of the value of ecosystems and their services in monetary units. Ecosystem Services, 2012, 1(1): 50-61.
[32]	Guo Z W, Xiao X M, Gan Y L, et al.Ecosystem functions, services and their values: A case study in Xingshan county of China. Ecological Economics, 2001, 38(1): 141-154.
[33]	Li Huancheng.GIS-based apporoach for regional ecosystem services valuation: Methodology and applications [D]. Hangzhou: Zhejiang University, 2010.
	[李焕承. 基于GIS的区域生态系统服务价值评估方法研究与应用. 杭州: 浙江大学, 2010.]
[34]	Tallis H, Polasky S.Mapping and valuing ecosystem services as an approach for conservation and natural-resource management. Annals of the New York Academy of Sciences, 2009, 1162(1): 265-83.
[35]	Shi Longyu, Cui Shenghui, Yin Kai, et al.The impact of land use/cover change on ecosystem service in Xiamen. Acta Geographica Sinica, 2010, 65(6): 708-714.
	[石龙宇, 崔胜辉, 尹锴, 等. 厦门市土地利用/覆被变化对生态系统服务的影响. 地理学报, 2010, 65(6): 708-714.]
[36]	Su C H, Fu B J, He C S, et al.Variation of ecosystem services and human activities. Acta Oecologica, 2012, 44: 46-57.
[37]	Lee R. Forest Hydrology.New York: Columbia University Press, 1980.
[38]	Liu Pengfei, Wu Jie.Economic values of bee pollination to China's agriculture. Scientia Agricultura Sinica, 2011, 44(24): 5117-5123.
	[刘朋飞, 吴杰. 中国农业蜜蜂授粉的经济价值评估. 中国农业科学, 2011, 44(24): 5117-5123.]
[39]	Li J, Ren Z Y, Zhou Z X.Ecosystem services and their values. Ecological Research, 2006, 21(4): 597-604.
[40]	Li Guang, Huang Gaobao,Wang Qi, et al.Spatial-temporal distribution of the ecological service value and NPP of cuItivated land on Longdong Loess Plateau . Acta Pratacultuae Sinica, 2011, 20(6): 18-25.
	[李广, 黄高宝, 王琦, 等. 陇东耕地净第一性生产力及生态服务价值的时空分异研究. 草业学报, 2011, 20(6): 18-25.]
[41]	Zhou G S, Wang Y H, Jiang Y L, et al.Estimating biomass and net primary production from forest inventory data: A case study of China's Larix forests. Forest Ecology and Management, 2002, 169(1/2): 149-157.
[42]	Zhang G G, Kang Y M, Han G D, et al.Effect of climate change over the past half century on the distribution, extent and NPP of ecosystems of Inner Mongolia. Global Change Biology, 2011, 17(1): 377-389.
[43]	Wang Yanyan, Yang Mingchuan, Pan Yaozhong.Valuation of the organic matter produced by terrestrial vegetation ecosystems in China in 1995 based on remote sensing. Ecology and Environmental Sciences, 2005, 14(4): 455-459.
	[王艳艳, 杨明川, 潘耀忠. 中国陆地植被生态系统生产有机物质价值遥感估算. 生态环境学报, 2005, 14(4): 455-459.]
[44]	Zhu Hongguang, Qin Pei, Wan Shuwen, et al.Emergy analysis of two models of wetland utilization on seashore in Jiangsu province. Chinese Journal of Ecology, 2001, 20(1): 38-44.
	[朱洪光, 钦佩, 万树文, 等. 江苏海涂两种水生利用模式的能值分析. 生态学杂志, 2001, 20(1): 38-44.]
[45]	Li Cuizhen, Kong Xiangbin, Sun Xianhai.Cultivated land resources value system and its evaluation in Beijing. Acta Geographica Sinica, 2008, 63(3): 321-329.
	[李翠珍, 孔祥斌, 孙宪海. 北京市耕地资源价值体系及价值估算方法. 地理学报, 2008, 63(3): 321-329.]

一级类		二级类		生态功能单位价值(元/hm²)	生产功能单位价值(元/hm²)	生活功能单位价值(元/hm²)	主导功能	具体功能识别			具体组合模式
01	耕地	011	水田	19289	64637	17244	2	1	2	3	生产主--生态副--生活次空间
01	耕地	013	旱地	15607	26487	15564	2	1	2	3	生产主--生态副--生活次空间
02	园地	021	果园	45270	87343	20517	2	1	2	3	生产主--生态副--生活次空间
		022	茶园	25894	36666	15967	2	1	2	3	生产主--生态副--生活次空间
		023	其他园地	29457	35157	15836	2	1	2	3	生产主--生态副--生活次空间
03	林地	031	有林地	60250	28172	21576	1	1	2	3	生态主--生产副--生活次空间
03	林地	033	其他林地	21522	18187	16664	1	1	2	3	生态主--生产副--生活次空间
04	草地	043	其他草地	10363	9332	9767	1	1	2	3	生态主--生产副--生活次空间
10	交通运输用地	102	公路用地	0	88991	153977	3		2	3	生活主--生产副空间
11	水域及水利设施用地	111	河流水面	12589	10149	8772	1	1	2	3	生态主--生产副--生活次空间
11	水域及水利设施用地	114	坑塘水面	7983	0	0	1	1			生态主空间
12	其他土地	122	设施农用地	12612	136600	0	2	1	2		生活主--生态副空间
12	其他土地	127	裸地	3217	0	0	1	1			生态主空间
20	城镇村及工矿用地	202	建制镇	0	22235	690366	3		2	3	生活主--生产副空间
		203	村庄	0	285514	217508	3		2	3	生活主--生产副空间
		204	采矿用地	0	8710	0	2		2		生产主空间
		205	风景名胜及特殊用地	18024	21600	348217	3	1	2	3	生活主--生态副--生产次空间