由于经济的快速发展，我国在2002年就取代美国成为世界第一大天然橡胶消费国^[1]。2010年中国天然橡胶产量不足70万t，自给率不足20%，远低于国际公认的30%安全保障线^[2]。橡胶园的地力直接影响天然橡胶的生产量，海南农垦是中国最大的天然橡胶生产基地，自第二次全国土壤普查以来其橡胶园养分状况已发生了很大变化，干部职工并不了解当前橡胶园地力的空间变异及其影响，必需摸清橡胶园地力水平，提升其生产能力。耕地地力评价技术正快速向规范化方向发展，其评价方法日益完善和丰富^{[3, 4, 5, 6, 7, 8, 9]}。因此，运用耕地地力评价方法查清主要橡胶园的养分状况和地力水平并提出相应的措施，对橡胶园土肥资源的优化配置和海南天然橡胶产业的可持续发展具有重要的意义。

海南农垦地处18°15′35″—20°02′14″ N，108°42′09″ — 110°47′42″E，属热带季风气候，年积温高、日照充足、雨量丰富，四季不明显，但干湿季节明显，东部湿润多台风、西部干旱、南暖北凉。海南农垦现有25个国营橡胶种植基地分公司，遍布海南省18个市县(图 1)，面积和产量占全国的40%以上，橡胶职工达7.2万余人，天然橡胶种植面积约246068.7 hm²，其中已开割橡胶总株数约4747.5万株，是全国最大的天然橡胶生产基地。

图 1 海南农垦天然橡胶基地分公司分布图 Fig.1 the distribution map of China Hainan Rubber Industry Group Co.,Ltd. branch

图选项

根据20世纪80年代第二次土壤普查资料，海南垦区橡胶园可划分为花岗岩类、片麻变质岩类、玄武岩类、砂页沉积岩类4个主要土壤母质类型区。其中花岗岩类母质发育的土壤主要分布在垦区的中部并向四周辐射，是垦区天然橡胶的主要种植区，占海南垦区橡胶园面积的61.21%。花岗岩地区岩石主要由正长石、云母和石英组成，抗风化能力强，所形成土壤砂粒适中。

综合土壤类型、土地利用、管理界线、橡胶园面积和生产队数量等因素，平均约20 hm²布设一个样点，在研究区共布点采集7204个土壤样品，采集0—20 cm层土壤，记录样点坐标及所在橡胶园基本情况。将所采土壤样品在整洁无污染的室温下风干，研磨过2 mm孔径筛用于pH、有效磷、速效钾、有效态微量元素等含量及物理性状的测定，继续研磨过0.25 mm孔径筛用于有机质等含量的测定。

土壤质地采用指测法测定；pH采用电位法测定；有机质采用油浴加热重铬酸钾氧化容重法测定；有效磷采用氟化铵-盐酸浸提，钼锑抗比色法测定；速效钾采用乙酸铵浸提，火焰光度计法测定；交换性镁采用原子吸收分光光度计法测定^[10]。

一般认为月降雨量100 mm以上适宜橡胶生长，而海南岛降雨量存在明显的时空分布不均，区域上自东向西逐渐减少，一年之内有明显干湿季节^[11]；另外，在同一土壤母质类型区中橡胶园的肥力特性、土壤理化性状以及橡胶树栽培管理相类似，而不同土壤母质类型区橡胶园之间，上述条件和特点必有差异^[12]；故本研究选取降雨量和土壤母质类型作为海南农垦橡胶园类型区划分的两个因素。

以海南岛橡胶园四大土壤母质类型区划为基础，进一步考虑降雨量对橡胶树生长的影响，利用海南岛多年月平均降雨量图提取出月平均降雨量100 mm等值线，然后将该等值线与海南岛橡胶园土壤母质类型区划图叠加，对海南岛橡胶园进行分区，获取海南岛橡胶园的不同类型分区。本研究选取面积最大的花岗岩类母质土壤区中月均降雨量100 mm以上的橡胶园作为研究区的评价范围(以下称为“花岗岩类多雨区”)。

邀请相关行业专家，依据必要性、重要性、差异性、稳定性、易获取、精简性“六”原则，运用德尔菲法(Delphi Method,DM)选取pH、有机质、有效磷、速效钾、交换性镁、成土母质、质地、坡度、坡向、有效土层厚度和年降水量(共11个)作为花岗岩类多雨区橡胶园的地力评价指标，借助模糊数学(FMT)的隶属函数模型确定各评价指标的隶属度。1965年Zadehl提出模糊数学理论^[13]，以后不断得到完善和广泛应用，利用模糊子集和隶属函数确定评价指标的隶属度是其核心^{[14, 15]}。

隶属函数分为5种：戒上型、戒下型、峰型、直线型、概念型(散点型)，本研究中评价指标的隶属函数属于峰型、戒上型和概念型。戒上型、峰型隶属函数，是利用特尔菲法对评价指标一组实测值通过专家打分评估出相应的一组打分值(隶属度)，并根据这两组数据拟合隶属函数建立(表 1)。概念性隶属函数则是直接通过专家经验打分建立(表 2)。

通过层次分析法(AHP)确定各评价指标的权重。首先，根据研究区地力评价指标相互关系构建其层次结构图(图 2)；然后，邀请相关行业专家讨论、比较同层次每个指标与上一层次的相对重要性，并量化评估打分，形成判断矩阵；再计算矩阵最大特征根及特征向量进行一致性检验，直至通过一致性检验。

图 2 评价要素层次分析结构图 Fig.2 Hierachy plot of soil productivity and assessing indicators

图选项

将模糊数学确定的评价指标隶属度和层次分析法确定的权重，通过加权求和计算出其对应的评价单元橡胶园的地力综合指数(IFI)，计算方法如下：

式中,IFI为橡胶园地力综合指数，F_i为第i个评价指标的隶属度；C_i为第i个评价指标对应组合权重值。根据综合指数计算结果，采用累积频率曲线分级法来划分研究区橡胶园的地力等级，见表 3结果。

数据的数理统计与图表制作采用SPSS 16.0、Excel 2010和ArcGIS 10.0来完成。

由图 3可知，红线为月均降雨量100 mm等值线，该线纵贯海南岛南北，将玄武岩、花岗岩、片麻变质岩和砂页沉积岩四大土壤母质类型区分割为东西两部分。月均降雨量100 mm等值线以东地区降雨充沛，橡胶树生长很少遇到干旱的影响，而该线以西地区则会 受到干旱的影响。本研究将月均降雨量100 mm等。值线以东地区称为“多雨区”，即红线以东区域；月均降雨量100 mm等值线以西地区称为“少雨区”，即红线以西区域。本研究以土壤母质类型和降雨量为主要因素将海南农垦橡胶园划分为玄武岩类多雨区、玄武岩类少雨区、花岗岩类多雨区、花岗岩类少雨区、片麻变质岩类多雨区、片麻变质岩类少雨区、砂页沉积岩类多雨区、砂页沉积岩类少雨区等8个类型区

图 3 不同类型橡胶园分区图 Fig.3 The district map of different type rubber plantations

图选项

从图 3可以看出，花岗岩类多雨区面积占了海南岛一大半多，主要分布在琼中、儋州、白沙、屯昌、五指山、保亭、乐东、万宁、陵水、澄迈等市县，是海南农垦天然橡胶的主要种植区(图 1)，故本研究确定花岗岩类多雨区为评价范围，对海南农垦橡胶园生产具有更广泛的指导意义。

表 4表明，花岗岩类多雨区橡胶园面积140174.4 hm²，占海南农垦橡胶园面积的56.97%；其次是片麻变质岩类多雨区和砂页沉积岩类多雨区橡胶园，分别为31063.3 hm²和21211.5 hm²，分别占到海南农垦橡胶园面积的12.62%和8.62%；其余5个类型区橡胶园面积共53619.5 hm²，仅占海南农垦橡胶园面积的21.79%。

海南农垦花岗岩类多雨区橡胶园(以下简称为评价区)面积为140174.4 hm²，本研究以地力综合指数为依据，采用累积频率曲线法将橡胶园地力划分为5个等级。其中，一等地橡胶园面积8328.0 hm²，分布面积最小，占参评橡胶园面积的5.94%；二等地橡胶园面积42462.5 hm²，分布面积第二大，占参评橡胶园面积的30.29%；三等地橡胶园面积53228.6 hm²，分布面积最大，占参评橡胶园面积的37.97%；四等地橡胶园有一定面积分布，面积24661.5 hm²排第三大，占参评橡胶园面积的17.59%；五等地橡胶园面积11493.8 hm²，分布面积较小，占参评橡胶园面积的8.20%，见表 5。

从图 4可以看出，评价区主要分布在琼中、屯昌、澄迈、儋州、万宁、保亭、乐东、三亚、陵水、定安、琼海、五指山、昌江、文昌和东方等15个县市。一等地橡胶园分布较散，而且面积很少，在白沙的北部和西部、琼中的中部和西部、定安的南部、万宁的中北部和西南部、三亚和乐东交界线上的中间部等均有少面积分布；二等地橡胶园主要分布于白沙、琼中和三亚；三等地橡胶园主要分布于儋州、白沙、屯昌、澄迈、保亭的中部、琼中的中北部；四等地橡胶园主要分布于屯昌、乐东、保亭、琼中的东南部；五等地橡胶园分布比较集中，主要分布于琼中的西北和东南部、儋州的中西部、保亭的中部、陵水的中部、乐东的东北部。

图 4 评价区橡胶园地力水平分布图 Fig.4 The distribution map of soil productivity level of rubber plantations in assessment area

图选项

从表 6中可以看出，评价区土壤速效钾、有机 质、有效磷和pH的标准误都很小，表明该区橡胶园

土壤采样测试的pH、有机质、有效磷和速效钾含量均具有很好的代表性；而该区pH平均值为4.68，说明其橡胶园土壤整体呈酸性；其中橡胶园土壤速效钾、有机质和有效磷含量的平均值分别为42.39 mg/kg、10.51 g/kg和4.62 mg/kg，故此区有机质、有效磷和速效钾含量都不高。

表 6表明，评价区不同地力等级橡胶园之间土壤的pH和有机质含量均达到显著差异，但只有pH和速效钾是随着地力水平升高而含量增加，在地力水平高的橡胶园土壤有效磷含量较高，有机质含量高低与地力水平高低不相关；另外，土壤有效磷含量只有二等地与四等地之间没有达到显著差异，土壤速效钾含量只有三等地与四等地之间没有达到显著差异，其余各地力水平橡胶园之间的土壤有效磷或速效钾含量均达到了显著差异；说明上述4个指标对该区橡胶园的地力水平具有重要影响，尤其是pH和有机质的影响较大。

耕地地力是指由特定气候区域地形、地貌、成土母质、土壤理化性状、农田基础设施、培肥水平等要素综合构成、影响并决定的耕地生产能力^[16]。耕地地力评价是以利用方式为目的，估价耕地生产潜力和土地适宜性的过程，揭示生物生产力的高低和潜在生产力^[17]。美国国际环境与发展研究所和世界资源研究所在《世界资源》中认定种植天然橡胶的橡胶园是耕地的一种类型。1984年全国农业区划委员会和原国家土地管理局在《土地利用现状调查技术规程》中把种植农作物的土地确定为耕地，《辞海》(经济分册)把耕地定义为经过开垦用以种植农作物并经常耕耘的土地。海南省在2010年以公文《琼农字[2010]73号》确定天然橡胶为海南省的主要农作物。因此，橡胶园地力评价的对象是橡胶园生产能力，属于耕地地力评价的一种，可以看成是土地评价学的一个新分支。橡胶园是海南垦区农业生产最重要的资源，其地力的好坏直接影响到天然橡胶的可持续发展，重视橡胶园的质量建设、落实橡胶园地力评价、实现橡胶园土壤肥料科学管理是海南垦区的重要课题。

近些年，很多学者借助新兴的“3S”技术，采用系统聚类法、层次分析法和模糊评价法等建立了耕地地力评价指标体系，开展了大量的耕地地力评价方面的研究^{[3, 4, 5, 6, 7, 8, 9]}。其中，跨行政区域范围的土地评价研究也有报道，如新疆天然草地类型质量评价^[17]、湖北省清江流域旅游景区开发潜力评价^[18]及何毓蓉等对川江流域及其周边典型农业生态区的耕地地力评价^[19]。另外，欧阳进良等依托GIS还评价了种植不同作物的土地适宜性，并划分出作物种植区^[20]。但是，国内大量的地力评价研究几乎全部集中于大田作物，而对诸如橡胶树等此类多年生高大乔木种植园区的地力评价研究很少。关于橡胶园的地力评价仅见于数篇文献^{[21, 22, 23]}，又局限于农场区域小范围、调查采样量少或指标内容片面等不足。而在跨行政区域范围的分区基础上开展耕地地力评价，尤其对天然橡胶这类特殊作物的园地，进行大量采样、系统、定量、全面的跨行政区域范围的地力评价研究，目前还未见报道，本次研究尚属首次，进一步完善了耕地地力评价方法和技术在各类作物耕地上的应用，尤其是积累了耕地地力评价方法在橡胶园这类特殊而又重要作物上管理养分等方面的运用经验，这对国内外橡胶园的地力评价与养分管理工作具有重要的借鉴意义。

在海南岛大区域尺度下，土壤母质类型与气候往往是影响橡胶树生长的主要因素，其中土壤母质类型对橡胶树生长的影响最明显^[11]，不同地力水平橡胶园土壤母质类型的养分状况差异很大(表 6)。气候条件中，降雨和温度又是影响橡胶树生长的关键因子。海南岛不同区域之间温度的差异不大，且年均温都在23—25℃之间，满足橡胶树生长对温度(20—30℃)的要求。与温度的均匀性不同，海南岛不同区域尤其是东部和西部之间降雨量的差异很明显，降雨的差异显著地影响到橡胶树的生长和产胶量。海南岛月平均降雨量都在75 mm以上，都是可以满足种植橡胶的需求。但是，一般月平均降雨量在100 mm以下的区域一年中都会发生不同程度的旱害，对橡胶树的生长造成不利影响^[11]。因此，依据对橡胶树影响大的土壤母质和降雨量划分评价区域，其评价结果与橡胶树生长环境更吻合，更适宜为橡胶园养分的合理管理与利用以及配方施肥提供科学依据。

相对于王晶晶和韦娉婷针对单个农场橡胶园地力的评价^{[22, 23]}，或袁秀杰等^[24]对平原和丘陵区耕地地力评价及其两套指标衔接的研究来讲，本研究覆盖了整个海南农垦花岗岩类多雨区的橡胶园，地力评价结果能代表整个海南岛此类型区橡胶园的地力状况，跨农场区域且同地力评价体系的研究使得各农场间橡胶园地力等级具有统一可比性，其评价结果无需通过复杂的线性模型转换后才具可比性，对橡胶园生产实践具有更广泛的指导作用。吴敏等^[21]虽然评估了海南岛整个橡胶园的土壤肥力状况，但由于其评价指标仅限于土壤化学属性(如pH、有机质、全氮、有效磷等)，而未涉及诸如立地条件、气候等因素，其评价指标并不能很准确地反映整个海南岛橡胶园的综合地力；同时，其采样点个数也较少(共202个土样)；相比而言，本研究采集了7204个土壤样品，评价指标涵盖了理化性状、立地条件、养分状况、剖面性状等4个方面11项指标，其评估结果更系统、更全面，能更准确地反映海南岛同类型区橡胶园的地力状况。

本研究以占海南农垦橡胶园面积56.97%的花岗岩类多雨区为代表，其研究结果可以对海南农垦主要天然橡胶种植园区的管理提供理论支持。研究结果发现，评价区地力以中等水平为主(二等地、三等地、四等地)，占到其总面积的85.85%，中等地力橡胶园遍布海南省琼中、屯昌、澄迈、儋州、万宁、保亭、乐东、三亚、陵水、定安、琼海、五指山、昌江、文昌和东方等15个县市。

研究结果表明评价区土壤属于酸性，pH位于适合橡胶树生长的4.5—6.5(pH)之间；但表 6中数据显示pH下降时地力等级同步下降，且各等级橡胶园间土壤pH均达到显著差异，说明该区橡胶园土壤进一步酸化将会降低其地力水平，从而可能影响到橡胶树的正常生长，应该重视生产中对该类橡胶园土壤的酸碱度改良，适当提高其pH。

研究结果还发现，一等地和五等地橡胶园土壤的有机质、有效磷和速效钾含量，与中等水平(二等地、三等地、四等地)橡胶园土壤中对应养分含量达到了显著差异，说明土壤有机质、有效磷和速效钾的含量是影响该区橡胶园地力水平的主要因素，应该要加强其上述主要养分的管理。根据适宜橡胶树正常生长所需土壤养分含量^[25](有机质为20—25 g/kg，有效磷为5—8 mg/kg，速效钾为40—60 mg/kg)来判断，表 6中数据说明评价区土壤的有机质含量远低于橡胶树正常生长的要求水平，其有效磷含量也不能满足橡胶树正常生长；而该区一等地和二等地橡胶园土壤速效钾含量水平达到了橡胶树正常生长所要求的范围，其三等地和四等地的才刚刚达到，五等地橡胶园土壤速效钾含量水平甚至低于橡胶树正常生长所要求的范围。因此，提高该区橡胶园土壤有机质、有效磷和速效钾等养分的含量，将会提升其地力水平；建议在大力增施有机肥的基础上增施磷肥，尤其是增施钙镁磷肥等碱性肥调节土壤酸碱度；同时，加强氮钾养分的管理，提高其利用率，在橡胶树抽叶等关键期补施氮钾肥，平衡氮钾的流失。此外，要在橡胶园间作豆科绿肥作物，诸如爪哇葛藤、兰花毛蔓豆和无刺含羞草等，以增加有机质，培肥土壤。