上海崇明生态农业园区花菜田的碳平衡

王适，吴永波，陈杰，丁晓叶

（1.南京林业大学森林与环境资源学院，江苏南京 210037；2.同济大学现代农业科学与工程研究院，上海 200092）

上海崇明生态农业园区花菜田的碳平衡

王适1，吴永波1，陈杰2，丁晓叶1

（1.南京林业大学森林与环境资源学院，江苏南京 210037；2.同济大学现代农业科学与工程研究院，上海 200092）

以上海崇明岛生态农业园区花菜Brassica oleracea var.botrytis田为研究对象，主要运用生物量法对花菜田植被碳储量、土壤呼吸及人工管理碳源进行测定，估算花菜田碳平衡情况，旨在通过对花菜田碳平衡的估算，探讨农田碳增汇减排措施。研究表明：①花菜在1个生育期内含碳率最小值为25.73%，最大值为41.61%，平均值为34.12%。②植被在成熟收割时植被碳储量为5.18 t·hm-2，土壤碳储量84.74 t·hm-2。③花菜田土壤呼吸速率在生长初期最高，然后逐渐减弱，土壤异养呼吸通量为2.38 t·hm-2。④生育期内人工碳源总排放量为1.81 t·hm-2。⑤花菜1个生育期净固碳量为0.99 t·hm-2，相当于固定二氧化碳4.01 t·hm-2。崇明岛花菜田碳平衡测定结果为微弱碳汇，针对地区土壤盐碱化现状，当前最主要增汇减排措施是合理配施有机肥，改善土壤，增加产量从而增加碳汇，其次还可以农林间种，调整人员管理结构，减少碳源。图3表1参18

生态学；生态农业园；碳汇；碳源；碳平衡；花菜田

Key words：ecology；eco-agriculture area；carbon sinks；carbon source；carbon balance；Brassica oleracea var.botrytisfarmland

近年来，人们发现农业生态系统中的碳库是全球碳库中最活跃的部分。在人类活动的影响下，农田土壤碳库正在迅速变化。农田碳平衡主要从2个方面研究：碳汇（系统固定碳）与碳源（系统释放碳）。现今农田碳汇多采用最大收获法（生物量法）来计算，该方法与调查森林地上植被碳蓄积量方法类似。如罗怀良[1]长期观测中丘陵地区农田植被碳储量，刘允芬等[2]对青藏高原农田植被碳蓄积进行了测算，丁晓叶[3]在对崇明岛人工林与农田碳储量研究时都用到生物量法。此外，还有赵荣钦等[4]的模型法，其主要是将作物有机碳含量与经济产量、经济系数相结合，粗略估算植被碳储量。农田碳源侧重土壤呼吸碳源研究，国内使用较多的是静态气室法和动态气室法[5]。农田系统存在特殊的人工管理投入碳排放。管理活动碳源主要包括投入农用化学品的能源碳源和农田人工管理活动（机耕、灌溉、收获）碳源。当前越来越多的学者把注意力投到农用化学品、灌溉和机耕活动碳源上，但对收获活动带来的碳源的关注较少。一直以来关于农田生态系统究竟是碳源还是碳汇的问题存在争议。一些学者如Lal等[6]认为：农田碳汇量巨大，全球耕地总的固碳潜力为0.75～1.0 Pg·a-1；而韩冰等[7]针对中国农业措施，提出中国农田固碳潜力为182.1 Tg·a-1；赵荣钦等[4]研究发现，中国沿海地区农田碳吸收总量为22 482.4万t，属于明显碳汇。而有部分学者则认为农田主要是一种碳源，如孙艳丽等[8]对华北平原冬小麦Triticum aesticum l（-1.98 t·hm-2），玉米Zeamays（-1.38 t·hm-2）轮作农田碳平衡进行测定，指出农田属于碳源，而Bouwman[9]也认为农业活动是重要的温室气体排放源。本研究运用生物量法，对上海乃至全国第1个“低碳农业园”崇明现代农业园区主要种植的花菜Brassica oleracea var.botrytis田，从植被动态含碳率和固碳量、土壤呼吸碳源释放及人工管理碳源几个方面定量分析花菜田系统碳平衡特征，通过估算判断其是碳汇还是碳源，并根据碳平衡特征探讨增汇减排技术，旨在为今后低碳农业园的建立提供理论指导与数据支持。

1 研究区概况

崇明岛位于31°27＇00″～31°51＇15″N，121°09＇30″～121°54＇00″E，处于北亚热带南缘，气候温和湿润，夏季湿热，冬季干冷，四季分明，具有明显的亚热带海洋性季风性气候。岛内常年平均气温为15.3℃，1月份最低，7-8月份气温最高，全岛东部年平均温度略高于西部。崇明岛年平均降水量为1 055.6 mm，降水主要集中于汛期。岛上土壤类型主要为冲积母质上发育而成的滨海盐碱土，此外还有淋溶土和半水成土。全年平均日照2 104 h，日照百分率为47%，2月日照时数最少，为131.1 h。11月至2月多刮北风和西北风[10]。实验时间为2011年9月初到2012月年1月中旬，地点位于陈家镇现代农业园内，地处崇明东部，以园区中花菜田为研究对象，花菜生育期从种植到收获历时140 d左右。实验样地表层土壤氮质量分数为1.83 g·kg-1，磷1.35 g·kg-1，有机碳为12.71mg·kg-1。

2 材料与方法

根据现有研究结论[11]，农田碳平衡研究公式为净生产力=净第一生产力-土壤异养呼吸碳源。广义的农田碳平衡还包括人工管理投入碳源。由于花菜生长期间凋落物生物量几乎为0，所以在本研究中忽略不计。本研究的土壤呼吸碳源采用土壤异氧呼吸碳释放量计算，忽略作物根呼吸的影响，因为在计算根的碳储量时，直接是扣除呼吸作用消耗的碳以后的碳储量，因此，在计算作物生物量时将作物地上部分与地下部分统一计算。

2.1 样地设置

在花菜种植期间，选择相邻的3块20 m×20m的田块作为实验样地，每块样地随机设置3块1 m ×1m的小样方，隔15 d在样方内用收获法收割花菜植株，封装后带回实验室测定。

2.2 植被含碳率与碳储量测定

将每次采集的花菜样品封装带回实验室保存，最后将所有采集的样品在粉碎前均放入85℃的恒温箱中烘至恒量，再经粉碎、碾磨、过筛后，采用重铬酸外加热法[12]测定植被含碳率，再结合植被生物量计算碳储量。C=B×O。其中：C为碳储量（t·hm-2），B为生物量（t·hm-2），O为含碳率。

2.3 土壤有机碳储量测定

选择样点3个·样地-1，随机设置3个重复·样点-1。由于农田基本在土壤表层耕作，不如林地对低层土壤影响作用大，所以研究时土壤剖面挖60 cm深，按0～5，5～10，10～20，20～40，40～60 cm分层取样，相同层次混合，同时用环刀法[13]测定土壤容重。采集的土样在室内自然风干后过2mm筛，用重铬酸钾氧化-外加热法[13]测定土壤有机质质量分数，进而计算土壤有机碳储量[14]：C0=∑0.1HiBiOi。其中：C0为土壤有机碳储量（t·hm-2）；0.1为单位换算系数；Hi为第i层土壤厚度（cm）；Bi为第i层土壤的容重（g·cm-3）；Oi是第i层土壤有机碳质量分数（g·kg-1）。

2.4 土壤呼吸测定

在1m×1m样方内随机设置3个重复，采用Li-6400便携式光合作用仪测量，安装09土壤呼吸叶室测定土壤呼吸。在测定前1天将标准规格的聚氯乙烯土壤隔离圈打入土中1～2 cm深，保证每个月的测定都在同一地点进行，免除测定点移位带来的干扰。每次测定的前1天，在不破坏土壤的情况下彻底剪除隔离圈内的地表植被，从而减少土壤扰动及根系损伤对测量结果的影响。测定时间为6：00-18：00，2 h测定1次，每个土壤隔离圈重复3次。

土壤异氧呼吸采用开沟法测定，在小样方四周挖掘70～100 cm深、50 cm宽的壕沟，沟中放入玻璃纤维薄片以阻止根向样方内生长，排除根系呼吸的干扰，然后把沟填平再测定，所测得的土壤呼吸即为土壤异养呼吸。

2.5 农田管理碳排放测定

管理碳排放主要包括能耗碳排放、灌溉碳排放和人工碳排放。能耗碳排放包括杀虫剂和肥料生产的碳排放。灌溉碳排放估算：Ci=Vi×w×h×ni。其中：Vi为煤电的碳强度系数（0.92 kg·kW-1·h-1）；w为灌溉所用的电机功率（kW）；h为每次灌溉工作时数；ni为灌溉次数。人工投入碳排放估算：Cr=Vr×N。其中：Vr为成人（60 kg体质量）每天呼出的二氧化碳容积；N为1个作物生长季投入的人工总数（人·d-1）。

根据现有中国生产肥料的化石能源碳排放研究：生产1 t杀虫剂排放碳4 931.93 kg·Mg-1；生产1 t氮肥排放1.74 t碳；1 t钾肥排放120.28 kg碳，1 t磷肥排放165.09 kg碳。

3 结果与分析

3.1 花菜含碳率变化

实验测得含碳率数据如图1所示。花菜生长期含碳率波动范围为30%～40%，生长初期最小，为25.73%，生长旺期最大，为41.61%，平均含碳率为35.12%。花菜只有在生长旺期才达到中国主要农作物平均含碳率水平[15]（水稻Oryza sativa为41.45%，玉米为47.09%，大豆Glycine soja为45.00%，小麦为48.53%等）。

3.2 花菜碳储量

将每次采集的花菜样品经烘干称量测得生物量，结合含碳率测得花菜碳储量如图2所示。花菜在生长初期生长缓慢，生物量基本无变化，此时的碳储量（0.279 4 t·hm-2）与移栽前（0.307 3 t·hm-2）相比略有降低，而在第4次测定时，也就是生长1个半月时花菜碳储量有明显增加，通过测定第4次到第8次之间为生长旺期，生物量增速最快，含碳率也相对较高，从而测得生长期最大碳储量为0.684 9 t·hm-2，一般收割均在生长末期之前，生长末期碳储量逐渐减少，但变化不大。综合计算花菜1个生育期固碳量为5.18 t·hm-2。

3.3 土壤有机碳储量及分布

如表1所示：花菜田0～5 cm土壤有机碳最高为15.77 g·kg-1，是底层土壤8.59 g·kg-1的1.84倍。土壤总碳储量为84.74 t·hm-2，其中表层土（0～20 cm）总碳储量为31.25 t·hm-2，占总碳储量的36.88%。这可能与表层土壤经常翻耕，大部分化肥等原料停留在表层有关。

表1 花菜田各土壤层有机碳储量Table 1 Organic carbon storage of the Brassica oleracea var.botrytis soil layer

图1 花菜生育期内含碳率动态变化Figure 1 Dynamics of carbon content of Brassica oleracea var.botrytis in one growth period

图2 花菜生长期间碳储量动态变化Figure 2 Dynamics of carbon storage of Brassica oleracea var.botrytis in one growth period

3.4 土壤呼吸月变化

从图3可以看出：花菜生长初期土壤呼吸速率普遍较高，第4次测得最高总呼吸速率为5.46μmol· m-2·s-1，此时花菜开始进入生长旺期。之后随着花菜的生长，总呼吸速率反而降低，逐渐减弱直至收割时最低为0.47μmol·m-2·s-1，生育期内总呼吸平均呼吸速率为3.16μmol·m-2·s-1。

花菜田异养呼吸速率动态变化与总呼吸基本一致，也是生长旺期开始时最高，成熟后逐渐降低，收获时最低，异养呼吸速率变化范围为0.43～2.84μmol·m-2·s-1，平均值为1.72μmol·m-2·s-1。

土壤异养呼吸主要受土壤温度、湿度等环境因子及微生物活动的影响，将它与自养呼吸分开讨论，对于研究未来全球碳动态变化具有十分重要的意义。

3.5 花菜田管理碳源

据统计，花菜种植期间施农药0.043 t·hm-2，化肥9.000 t·hm-2，其中磷肥2.000 t·hm-2，钾肥3.000 t·hm-2，氮肥4.000 t·hm-2，计算得能源总碳排放为0.790 t·hm-2；栽植期间共灌溉3次，6 h·次-1，收获时雇佣17人以人工方式收割，无机耕活动，计算得人工管理活动碳排放为1.02 t·hm-2。统计得花菜1个生育期内管理碳排放为1.81 t·hm-2。

3.6 花菜田碳平衡计算

由于花菜生育期内凋落物量几乎为0，所以在计算中忽略不计。通过对各子系统碳汇源计算发现，花菜1个生育期内植被固碳量为5.18 t·hm-2，花菜田土壤释放碳量为2.38 t·hm-2，农田管理碳排放量为1.81 t·hm-2，花菜田表现为弱碳汇，测得净碳汇为0.99 t·hm-2，相当于固定二氧化碳4.01 t·hm-2。

4 讨论

图3 花菜田土壤总呼吸与异养呼吸月变化Figure 3 Month dynamics of soil total and heterotrophic respiration for Brassica oleracea var.botrytis farmland

实验测得花菜含碳率基本为30%～40%，平均值为35.12%，结合生物量计算得花菜植被生育期内固碳量为5.18 t·hm-2。与前人研究相比，低于河南1997年平均固碳量（5.65 t·hm-2）和2007年平均固碳量（9.77 t·hm-2）[4]，但高于江西2003年和2007年平均固碳量（1.12 t·hm-2和1.51 t·hm-2）[16]，同时介于李洁静等[17]测算的太湖地区水稻（4.7 t·hm-2）和油菜（8.0 t·hm-2）固碳量之间。花菜含碳率较低，固碳能力较弱，但花菜植被固碳是花菜田系统的主要碳汇方式，因此，增加植被碳储量是最好的碳增汇措施，而要增加碳储量主要是通过增加花菜产量实现。由于崇明岛地区土壤盐碱化较严重、土地贫瘠，目前，最适宜的增汇方式就是增加有机肥配施，但有机肥施用也不能无限制增加，如何找到一个最利于作物生长的比例还需要长期实验研究。此外，还可以在同一块田地上进行轮作，通过不同品种的更替减缓土壤流失，增加土壤碳储量，从而促进作物生长。

花菜田1个生育期内异养呼吸碳排放量为2.38 t·hm-2。呼吸作用是农田最主要碳源，因此，要增加花菜田碳汇能力就要减少碳源释放，也就是减小土壤呼吸作用。农田土壤异养呼吸速率一般比人工林土壤异养呼吸速率大，因为农田郁闭度相对较低，阳光能更多的到达地面，从而引起农田土温快速上升，再加上农田施用肥料，使得土壤微生物活性较强，进而释放出更多二氧化碳。针对这一特点，在构建低碳园区时可以采取农林间作的方式，降低农田郁闭度，进而降低土壤表面温度，以减弱农田土壤呼吸。

花菜在1个生育期内管理碳源总排放量为1.81 t·hm-2。由于农田对管理的特殊依赖性，人工管理碳源也是碳释放的重要途径，要减少人工管理碳源就要求园区进行不断管理改革，改变农田管理方式，减少人员或机械投入，或者引进节能机器等从而减少人工管理碳源。

花菜田生育期净碳汇为0.99 t·hm-2，表现为弱碳汇。这一数值与中国其他学者关于农田碳汇研究结果相比，低于朱咏莉等[18]测得的亚热带稻田净固定碳量，但与郝庆菊等[5]测定的水旱轮作农田的净碳汇值相当。本研究主要根据现今崇明岛农业园区耕作方式，选择园区主要种植品种花菜进行研究，虽然园区目前已将部分田地退田还林，但人工林种植面积相对于农田较小，且种植时间较短，固碳能力较弱。园区主要碳流动载体依然以农田为主，如何根据农田碳流动方式采取碳增汇减排措施是当前最主要问题，也是构建低碳农业园区最主要环节。

研究中将作物生长期的生物产量作为碳吸收，实际上已包括了作物凋落物所带来的碳，并扣除了作物呼吸作用消耗的碳。如果需要进一步细化不同时期的碳流动方式，还需进一步测定各时期作物光合作用与呼吸作用对碳变化量的影响，再结合土壤呼吸与根呼吸释放碳量进行分析。本研究中对人工管理活动碳排放的计算只适用于崇明地区，因为不同地区、不同耕作方式所需要的人力物力不同，从而计算结果也会有所不同。其次是燃料标准，由于目前对燃料没有较为严格的要求和规范，所以在计算时也存在一定的不确定性。同一种燃料碳的含量可能有很大差别，而在选择碳排放系数时，还要根据生产化肥、杀虫剂的过程不同，以及灌溉耗电对化石燃料的利用不同来进行选取。

5 结论

研究表明：崇明农业园区花菜田生态系统表现出微弱碳汇能力，花菜田系统固碳主要是通过植被固碳，主要碳释放途径为土壤呼吸与人工管理碳释放。现阶段关于崇明岛农业园区农田碳增汇减排技术最可行的措施就是改善施肥制度，比如等氮投入、不同氮源施肥、有机-无机肥配施处理等。建议通过不同配比实验，寻找最适合当地土壤条件的有机肥配施标准，提高作物生物量，并最大限度地减少土壤碳排放，在耕作强度很大的农田，最有效的办法是提高土壤中稳定态氮含量及比例。

[1]罗怀良.中丘陵地区近55年来农田植被碳储量动态研究：以四川省盐亭县为例[J].自然资源学报，2009，24（2）：251-258.

LUO Huailiang.Dynamic of vegetation carbon storage of farm land ecosystem in hilly area of central Sichuan Basin during the last 55 years：a case study of Yanting County，Sichuan Province[J].JNat Resour，2009，24（2）：251 -258.

[2]刘允芬，欧阳华，曹广明，等.青藏高原东部生态系统土壤碳排放[J].自然资源学报，2001，16（2）：152-l60.

LIU Yunfen，OU Yanghua，CAO Guangmin，et al.Soil carbon emission from ecosystems of eastern Qinghai-Tibet Plateau[J].JNat Resour，2001，16（2）：152-160.

[3]丁晓叶.崇明现代农业园区人工林和农田碳平衡研究[D].南京：南京林业大学，2011.

DING Xiaoye.Research on Carbon Balance of Plantation and Farm land in Chongming Modern Agricultural Zone[D]. Nanjing：Nanjing Forestry University，2011.

[4]赵荣钦，刘英，丁明磊，等.河南省农田生态系统碳源/汇研究[J].河南农业科学，2010（7）：40-44.

ZHAO Rongqin，LIU Ying，DING Minglei，et al.research on carbon source and sink of farm land ecosystem in Henan Province[J].JHenan Agric Sci，2010（7）：40-44.

[5]郝庆菊，王跃思，宋长春，等.三江平原农田CO2收支研究[J].农业环境科学学报，2007，26（4）：1556-1560.

HAO Qingju，WANG Yuesi，SONG Changchun，et al.CO2budget in agro-ecosystems in the Sanjiang Plain[J].JA-gro-Environ Sci，2007，26（4）：1556-1560.

[6]LAL R，BRUCE JP.The potential ofworld cropland soil to sequester C andmitigate the greenhouse effect[J].Environ Sci&Polic，1999，2：177-185.

[7]韩冰，王效科，逯非，等.中国农田土壤生态系统固碳现状与潜力[J].生态学报，2008，28（2）：612-619.

HAN Bing，WANG Xiaoke，LU Fei，et al.Soil carbon sequestration and its potential by cropland ecosystems in China[J].Acta Ecol Sin，2008，28（2）：612-619.

[8]孙艳丽，陆佩玲，李俊，等.华北平原冬小麦/夏玉米轮作田土壤N2O通量特征及影响因素[J].中国农业气象，2008，29（1）：1-5.

SUN Yanli，LU Peiling，LI Jun，et al.Characteristics of soil N2O flux in a winter wheat-summermaize rotation system in north China Plain and analysis of influencing factors.[J].Chin JAgrometeorol，2008，29（1）：1-5.

[9]BOUWMAN A F.Soils and the Greenhouse Effect[M].Chichester：John Wiley&Sons，1990.

[10]崇明县人民政府.崇明岛概况（地理环境）[EB/OL].上海崇明岛政府网.[2013-03-04].http：//www.cmx.gov.cn.

[11]BALDOCCHID D.Asscessing carbon balance：problems and prospects of the eddy covariance technique[J].Global Change Biol，2003，9：478-492.

[12]黄斌，王敬国，龚元石，等.冬小麦夏玉米农田土壤呼吸与碳平衡研究[J].农业环境科学学报，2006，25（1）：156-160.

HUANG Bin，WANG Jingguo，GONG Yuanshi，et al.Soil respiration and carbon balance in winter wheat and summermaize fields[J].JAgro-Environ Sci，2006，25（1）：156-160.

[13]张树人，赵协哲，宋钢，等.吉林省土壤微量元素含量与分布规律研究：土壤中的铜和锌[J].吉林农业科学，1987，12（3）：50-57.

ZHANG Shuren，ZHAO Xiezhe，SONG Gang，et al.Trace elements content and distribution in soils of Jilin Province（I）copper and zinc in soil[J].JJilin Agric Sci，987，12（3）：50-57.

[14]方晰，田大伦，项文化，等.不同密度湿地松人工林中碳的积累和分配[J].浙江林学院学报，2003，20（4）：374-379.

FANG Xi，TIAN Dalun，XIANGWenhua，et al.On carbon accumulation，distribution of different densities in slash pine plantation[J].JZhejiang For Coll，2003，20（4）：374-379.

[15]赵荣钦，秦明周.中国沿海地区农田生态系统部分碳源/汇时空差异[J].生态与农村环境学极，2007，23（2）：1-6，11.

ZHAO Rongqin，QIN Mingzhou.Temporospatial variation of partial carbon source/sink of farmland ecolystem in coastal China[J].JEcol Rural Environ，2007，32（2）：1-6，11.

[16]金姝兰，杨芳英.江西农作物碳储量估算与分析[J].广东农业科学，2011（2）：216-219.

JIN Shulan，YANG Fangying.Carbon stock estimation and analysis of crops in Jiangxi[J].Guangdong Agric Sci，2011（2）：216-219.

[17]李洁静，潘根兴.太湖地区长期施肥条件下水稻-油菜轮作净碳汇效应及收益评估[J].应用生态学报，2009，20（7）：1664-1670.

LI Jiejing，PAN Genxing.An evaluation of net carbon sink effect and cost/benefits of a rice-rape rotation ecosystem under long-term fertilization from Tai Lake region of China[J].Chin JAppl Ecol，2009，20（7）：1664-1670.

[18]朱咏莉，吴金水，周卫军，等.亚热带稻田生态系统CO2排放及影响因素[J].中国环境科学2005，25（2）：151 -154.

ZHU Yongli，WU Jinshui，ZHOUWeijun，et al.CO2emission from the paddy in subtropical region and its influence factors[J].China Environ Sci，2005，25（2）：151-154.

Carbon balance of Brassica oleracea var.botrytis from Chongming Ecological Agriculture Park in Shanghai

WANG Shi1，WU Yongbo1，CHEN Jie2，DING Xiaoye1
（1.College of Forest Resources&Environment，Nanjing Forestry University，Nanjing 210037，Jiangsu，China；2. Modern Agricultural Science&Engineering Institute，Tongji University，Shanghai200092，China）

To contribute tomeasurement techniques for carbon sink additions and emission reductions，and to provide data to support the development of a low-carbon agricultural park，the carbon balance ofBrassica oleraceavar.botrytiswas estimated.To evaluate carbon balance，themain carbon storage components of B.oleracea var.botrytis，soil respiration，and carbon sources of artificialmanagementwere determined bymeasuring biomass in an eco-agricultural area on Chongming Island in Shanghai.Themain results showed that（1）the minimum carbon content ofB.oleracea var.botrytiswas 25.7%，the maximum was 41.6%，and the average was 34.1%.（2）At harvest，carbon storage of the vegetable layer was 5.18 t·hm-2，and the soil layer was 84.74 t·hm-2.（3）The respiration rate of the soilwas highest in the initial stages of growth，and then gradually weakened.For one growing period，soil heterotrophic respiration effluxes for theB.oleracea var.botrytisfarmland was 2.38 t·hm-2.（4）Total carbon emission was 1.81 t·hm-2，and（5）net carbon sequestration in one growing period was 0.99 t·hm-2fixing CO24.01 t·hm-2.Thus，to improve the weak carbon sink，reasonable fertilizer including organic fertilizer should be applied to increase production that will increase carbon sinks；as production increases，agroforestry intercropping and personnel management structure adjustments can be used to reduce the carbon source.[Ch，3 fig.1 tab.18 ref.]

S181

A

2095-0756（2014）02-0190-06

2013-03-26；

2013-06-04

“十一五”国家科技支撑计划项目（2009BAC62B04）

王适，从事环境生态学研究。E-mail：wangshi_198803@163.com。通信作者：吴永波，副教授，从事生态学研究。E-mail：yongbowu0920@163.com