薄荷的生境适宜性区划及品质区划研究.docx

1、薄荷的生境适宜性区划及品质区划研究薄荷的生境适宜性区划及品质区划研究 摘要 利用最大熵模型（MaxEnt）与地理信息系统（ArcGIS），结合生态因子数据、薄荷地理分布数据、薄荷挥发油含量数据，建立薄荷挥发油与生态环境因子之间的关系模型，进行我国薄荷药材品质区划研究。研究结果表明，影响薄荷分布的主要生态因子为年均降水量（45.87%）、最湿季平均温（11.92%）、最暖季平均温（7.84%）、5月平均温（6.80%）、温度季节性变化的标准差（4.42%）、最冷季平均温（3.47%）、海拔（2.92%）。最适宜薄荷生长的生态特征为：年均降水量在530，1 465 mm、最湿季平均温在24.5，2

2、9 ，最暖季平均温均在25.5，29 、5月平均降水量（prec 5）在67，133 mm、温度季节性变化的标准差在8 333，9 643、最冷季平均温为1.7，8.3 、海拔为0，165 mm。薄荷品质最佳区域主要集中分布在江苏东部及西南部、安徽中部、山东东部、浙江北部、黑龙江中东部等地区。该研究结果与江苏、安徽等地为薄荷的主产区，其中江苏太仓为薄荷的道地产区，安徽太和县是全国最大的薄荷生产基地的实际情况相吻合。该区划研究为更好地指导薄荷人工种植和生产提供科学依据。 关键词 薄荷； 最大熵模型； 地理信息系统； 生境适宜度； 品质区划 Abstract In this study， ecol

3、ogical factors， occurrence records， the essential oil components content were used to predict the potential geographical distribution and quality division of Mentha haplocalyx in China based on the MaxEnt modeling and geographic information system（GIS）. The AUC of ROC curve was above 0.950，indicatin

4、g that the predictive results with the maximum model were highly precise. The results showed that the main environmental factors determining the potential distribution were annual average precipitation （the contribution rate， 45.87%）， mean temperature of wettest quarter （11.92%）， mean temperature of

5、 warmest quarter （7.84%）， average monthly precipitation of May （6.80%）， standard deviation of seasonal temperature variation （4.42%）， mean temperature of the coldest quarter （3.47%） and altitude （2.92%）. The environmental variables in the highly potential areas were determined as annual average prec

6、ipitation around 530，1 465 mm， mean temperature of wettest quarter around 24.5，29 ， mean temperature of the warmest quarter around 25.5，29 ， average monthly precipitation of May around 67，133 mm， standard deviation of temperature seasonal change around 8 333，9 643， mean temperature of the coldest qu

7、arter around 1.7，8.3 and the altitude around 0，165 mm. The best quality distribution of M. haplocalyx was mainly located in Jiangsu， Anhui， Shandong， Zhejiang and Heilongjiang. The zoning results basically coincide with the actual situation. The quality division of M. haplocalyx can be used for prov

8、iding a scientific basis for selection of artificial planting base and guidance of its production.Key words Mentha haplocalyx； MaxEnt； ArcGIS； habitat suitability； quality division doi：10.4268/cjcmm20161709 薄荷为唇形科植物薄荷Mentha haplocalyx Briq.的干燥地上部分，是我国传统的常用中药，具有疏散风热，清利头目，利咽，透疹，疏肝行气等功效1-2。挥发油为其主要药效成分，

9、是评价薄荷质量的重要标准3。薄荷在中国大部分地区均有分布，其中以江苏、安徽、河南、江西等地为主产区，生于溪边、沟边等地，喜温暖湿润气候和阳光充足、雨量充沛的环境。薄荷用途很广，不仅作为药用，还可用于食品、化妆品、香料、烟草工业等行业，国内外需求量很大，为保证市场供求，栽培薄荷早已成为主要资源，但由于缺乏中药材引种和产地适宜性分析等相关知识，导致薄荷品质下降等诸多问题。因此，进行我国薄荷生境适宜性等级划分及品质区划研究，对于了解薄荷资源的整体概况和科学选取适宜的区域进行薄荷的人工种植具有十分重要的意义。最大熵模型（MaxEnt）是近年来用于预测生物物种潜在分布区最常用的统计软件4。1957年，J

10、aynes提出了最大熵理论，该理论认为，在已知条件下，熵最大的事物最接近它的真实状态5。地理信息系统（GIS）作为空间信息分析的核心技术，以其强大的空间信息分析和管理能力，在中药产地适宜性评价中得到广泛应用6-9。研究表明，MaxEnt模型结合ArcGIS地理信息系统在预测物种潜在分布区的性能要优于同类的其他模型，特别是在物种分布数据不全的情况下，应用MaxEnt模型预测能得到满意的结果10-12。 本研究基于生态因子数据、薄荷地理分布数据、薄荷挥发油含量数据，利用MaxEnt模型计算薄荷的生境适宜度，得到薄荷在全国范围内的生境适宜度分布图，并建立薄荷挥发油与生态环境因子之间的关系模型，最终结

11、合ArcGIS软件进行我国薄荷的品质区划研究，以期为薄荷的资源利用及其他薄荷相关研究提供理论依据。 1 材料和方法 1.1 薄荷数据 1.1.1 样品采集和分布数据 通过查阅文献收集薄荷主产区信息，通过实地调查获取薄荷药材样品，于20122014年在薄荷分布较多的江苏、安徽、河南等13个省份29个县市，选取56个具有代表性的样地采集薄荷样品，采用GPS获得采样点的经度、纬度和海拔等位置信息，见图1。 1.1.2 挥发油含量数据 按照2010年版中国药典（一部）挥发油测定法提取薄荷挥发油并计算挥发油含量。 1.2 生态因子数据 气候因子数据是根据我国气象观测站19502000年的气候观测数据插值

12、而成，空间分辨率为1 km1 km；土壤类数据根据第二次全国土地调查提供的1100万中华人民共和国土壤图制成，采用的土壤分类系统主要为FAO-90；地形数据包括：高程和坡向；植被类型数据以中科院植物研究所的中华人民共和国植被图（1100万）中的植被亚类数据制成。以上所有数据来源于中国中医科学院中药资源中心中药资源空间信息网格数据库。 1.3 薄荷潜在分布区预测分析 1.3.1 生境适宜度计算 为了计算的高效和简便，根据以下原则对生态因子进行筛选：MaxEnt模型筛选：使用MaxEnt模型计算各生态因子对薄荷的贡献率，删除贡献率为0的生态因子。生物学特性筛选：基于温度对薄荷挥发油成分的积累起关键

13、作用，降水对薄荷生存有较大影响，选择具有直接影响的生态因子。最终使用的生态因子见表1。 将薄荷分布数据和生态因子数据导入MaxEnt软件（3.3.3k），随机选取25%的采样点数据为测试集，75%的采样点数据为训练集11，输出格式为 ASCII 栅格图层，计算全国范围内薄荷的生境适宜度。 1.3.2 模型预测精度评价 分析遗漏/调试曲线评价，漏报率是在“存在数据”的训练集和检验集的基础上计算得到的。基于累积阈值的定义，漏报率应该和预测的遗漏值一致。受试者工作特征曲线评价。近年来，受试者工作特征曲线（ROC）分析法在物种潜在分布预测模型评价中应用越来越广泛13。采用ROC曲线与横坐标围成的面积即

14、AUC值来评价模型预测结果的准确度，AUC 0.500.60为失败，0.600.70为较差，0.700.80为一般，0.800.90为较好，0.901.0为非常好，越接近1说明环境变量与预测的物种地理分布之间的相关性越大，即模型预测的结果越准确14。 1.3.3 薄荷潜在分布区预测分析 生长区划分类：提取采样点处的生境适宜度值，以提取结果中的最小值（Min）作为适宜生长和不适宜生长的分界线。对于适宜生长的区间，则以正态分布中的+为界分为次适宜区和最适宜区。即0，Min为不适宜区，Min，+为次适宜区，+，1为最适宜区。利用ArcGIS 10.2软件的空间分析功能，计算得出分布适宜性等级分布图。

15、 1.4 薄荷品质区划分析 1.4.1 薄荷挥发油与生态环境因子间的相关性分析 基于采样点的位置信息，应用ArcGIS的信息提取功能，提取薄荷采样点所在位置的生态环境信息。应用“R”计算薄荷挥发油与生态环境因子间的关系模型，探寻薄荷药材品质和生态环境因子之间的关系。 1.4.2 薄荷品质区划研究 基于薄荷挥发油与生态环境因子间的关系模型，利用ArcGIS 10.2软件的空间分析功能，估算全国范围内薄荷挥发油含量的空间分布。基于薄荷的潜在分布区分析结果和挥发油含量空间分布进行叠加，分析符合药典标准的分布区域。 2 结果 2.1 各地薄荷中挥发油含量 采用2010年版中国药典（一部）挥发油测定法，

16、测定薄荷挥发油含量，结果见表2。 2.2 薄荷分布适宜性区分析 2.2.1 各生态因子对薄荷分布的影响 根据最大熵模型计算得到各生态因子对薄荷分布的贡献率，据此判断影响薄荷分布的主要生态因子，结果表明，年均降水量（45.87%）、最湿季平均温（11.92%）、最暖季平均温（7.84%）、prec5（6.80%）、温度季节性变化的标准差（4.42%）、最冷季平均温（3.47%）、altitude（2.92%）是影响薄荷分布的最主要的生态因子，以上7个生态因子的累计贡献率达83.24%。 年均降水量对薄荷的分布影响见图2，年均降水量对薄荷的分布影响如下：在465，530 mm，1 465，1 93

17、0 mm属于适宜区域，在530，1 465 mm属于最适宜。 最湿季平均温对薄荷的分布影响见图3，最湿季平均温对薄荷的分布影响如下：在18，24.5 ，29，34.5 属于适宜区域，在24.5，29 属于最适宜。最暖季平均温对薄荷的分布影响见图4，最暖季平均温对薄荷的分布影响如下：在19，25.5，29，33属于适宜区域，在25.5，29属于最适宜。 另外，经同上方法得出薄荷最适宜在5月降水量为67，133mm、温度季节性变化的标准差为8 333，9 643、最冷季平均温为1.7，8.3 、altitude为5，165mm的地区生长。综合以上分析表明，薄荷最适宜在年均降水量在530，1 465

18、mm、最湿季平均温在24.5，29，最暖季平均温均在25.5，29、5月平均降水量（prec5）在67，133mm、温度季节性变化的标准差在8 333，9 643、最冷季平均温为1.7，8.3、海拔（altitude）为0，165mm的地区生长。 2.2.2 薄荷潜在分布区预测 依据正态分布的参数和，结合采样点处的提取值，对生长区划和品质区划进行分类。由数据计算出生境适宜度的最小值0.152 8，最大值0.915 2，均值0.620 0，标准差0.215 985 5。 按照生境适宜度从低到高将薄荷分布区依次分为3个等级，分别为不适宜区（0，0.152 8）、次适宜区（0.152 8，0.835

19、 9）、最适宜区（0.835 9，1），结果见图5。江苏、安徽、山东、河南、湖北、辽宁、吉林、湖南、贵州、海南和甘肃等地是薄荷的主要分布区（次适宜生长区和最适宜生长区）。其中，薄荷生境最适宜生长区主要集中于江苏、安徽的部分地区。 2.2.3 模型预测精度评价 分析遗漏/调试曲线评价见图6，基于训练集的漏报率和预测的遗漏值一致，表明MaxEnt模型计算的生境适宜度的结果可信度较高。受试者工作特征（ROC）评价见图7，经受试者工作曲线检验，MaxEnt模型的ROC曲线下面积测试集的AUC为0.968，表明MaxEnt模型对薄荷潜在分布区的预测结果非常好，可信度很高。 2.3薄荷品质区划 2.3.1

20、挥发油与气候因子间的相关性分析选择 薄荷挥发油作为薄荷品种区划的指标，应用“R”进行挥发油与气候因子之间的关系分析，构建挥发油与辐射、温度、降雨量、日照时数和相对湿度年均值的回归模型7。 y1=-3.858-0.00005523x1-0.0006432x2+0.001347x3+0.06882x4（y1薄荷挥发油，x1辐射，x2降水量，x3为日照，x4相对湿度）。在显著性水平0.05条件下，由回归方程的F检验结果P=5.508e-050.4%的标准，进行重新分区，结果见图9。可知江苏东部及西南部、安徽中部、山东东部、浙江北部、黑龙江中东部等地区的薄荷的品质最佳。 3讨论 3.1薄荷生境适宜性的

21、主要影响因子 本论文采用MaxEnt模型对薄荷生境适宜性的主要影响因子分析研究表明，年均降水量、最湿季平均温、最暖季平均温是影响薄荷生长分布的最主要生态因子，其次是5月平均降水量、最冷季平均温、海拔等。薄荷生长过程中最主要的特点是喜湿润、怕干旱，所以年均降水量作为影响薄荷生长最主要生态因子这与薄荷的实际生长情况喜湿润、怕干旱相符合，且5月平均降水量对薄荷适宜等级的贡献率达6.80%，因薄荷采收期一般为7月和10月，采收期前5月份正是薄荷生长旺盛的时候，充足的雨水有利于薄荷的生长。最湿季平均温在24.529，最暖季平均温在25.529属于最适宜薄荷的生长，这与中华本草中提到薄荷喜温暖、湿润气候，

22、植株生长适宜温度为2030较一致15。最冷季平均温直接影响薄荷能否安全过冬，极寒冷地区不适宜其生存，因此该因子也成为影响薄荷生长分布的重要因子之一。薄荷在海拔为0165m属于最适宜生长区，这与薄荷在海拔2100m以下地区都可以生长，而以低海拔栽培，其精油和薄荷脑含量较高的生物学特性相一致15。 3.2薄荷的生境适宜性区划 对薄荷生境适宜性区划的研究表明，薄荷生境适宜区主要集中于江苏、安徽、河南、湖北、山东等地区，这些地区主要集中在低山丘陵、太湖平原、沿海、大别山脉、鄱阳湖等气候温湿多雨地区，这与薄荷的生长习性相一致。比较薄荷采样点地理分布示意图（图1）和MaxEnt预测出的薄荷生境适宜分布图（

23、图5）可以看出，薄荷样本记录的实际分布均位于潜在生境适宜分布区，说明预测的潜在分布区与实际分布高度一致。 3.3薄荷的品质区划分析 本文首次将薄荷中有效成分含量和生态环境因子相结合，建立了薄荷挥发油、主要成分与生态环境因子之间的关系模型对薄荷进行品质区划研究，区别于以往仅根据环境因子对物种生境适宜性的区划。研究结果表明，薄荷品质最佳区域主要集中分布在江苏东部及西南部、安徽中部、山东东部、浙江北部、黑龙江中东部等地区。这与江苏、安徽等地为薄荷的主产区，其中江苏太仓为薄荷的道地产区，安徽太和县是全国最大的薄荷生产基地的实际情况相吻合。本课题以品质结合生态因子为基础的薄荷药材品质区划研究，比仅以生态

24、因子为指标的区划结果能更合理地指导薄荷药材人工栽培区域的选取，并进一步提高薄荷的产量和品质。薄荷品质区划研究为薄荷在全国范围内大面积种植及推广提供了理论基础，为进一步进行薄荷药材资源发展规划提供了科学依据。 注本文地图插图审图号：GS（2016）2073号。 参考文献 1中国药典.一部S.2010：354. 2中国科学院中国植物志编辑委员会.中国植物志.第66卷M.北京：科学出版社，1977：264. 3郭晓恒.不同地区野生薄荷挥发油的比较J.天然产物研究与开发，2011，23（6）：1139. 4ElithJ，GrahamCH，AndersonRP，etal.Novelmethodsimpr

25、ovepredictionofspeciesdistributionsfromoccurrencedataJ.Ecography，2006，29：129. 5PhillipsSJ，AndersonRP，SchapireRE.MaximumentropymodelingofspeciesgeographicdistributionsJ.EcolModel，2006，190：231. 6郭兰萍，黄璐琦，蒋有绪，等.影响苍术挥发油组分的气候主导因子及气候适宜性区划研究J.中国中药杂志，2007，32（10）：888. 7张小波，郭兰萍，黄璐琦.我国黄花蒿中青蒿素含量的气候适宜性等级划分J.药学学报，

26、2011，46（4）：472. 8张小波，周涛，郭兰萍，等.基于地形因子的贵州省头花蓼生态适宜性等级划分J.中国中药杂志，2011，36（3）：311. 9张小波，郭兰萍，韦霄，等.广西青蒿种植气候适宜性等级区划研究J.中国中药杂志，2008，33（15）：1794. 10RobertoM，ZamoraR，RamnJ，etal.PredictivemodelingofmicrohabitatsforendemicbirdsinsouthChileantemperateforestsusingMaximumentropyJ.EcolInform，2011，6：364. 11YangXQ，KushwahaSPS，SaranS，etal.Maxentmodelingforpredictingthepotentialdistributionofmedicinalplant，JusticiaadhatodaL.inLesserHimalayanfoothillsJ.EcolEng，2013，51：83. 12车乐，曹博，白成科，等.基于MaxEnt和ArcGIS对太白米的潜在分布预测及适宜性评价J.生态学杂志，2014，33（6）：1623.