1. 前言
土壤有机碳的含量及其组分的分布在很大程度上影响着土壤碳库的稳定[1],土壤碳库的相对稳定性不仅影响土壤的固碳能力,并且会对土壤有机碳的排放过程产生影响,同时也是缓解温室效应的途径之一[2]。土壤有机碳组分的研究是近年来的研究热点[3]。土壤有机碳受耕作措施方面的影响最为直接有效,适宜的物料添加和耕作方式可以有效提高土壤的有机碳含量,有机碳的分组是以某一特性作为分组标准[4],划分出来时某一类的化合物的混合物,不同的分类方法间的组分在组成上会有重叠,并不是组成结构单一的化合物[5]。目前有机碳的分组方法主要有化学、物理和生物三种方法[6]。按物理性质分为颗粒碳有机碳,矿物态结合碳,主要存在于土壤团聚体的内部或外部,周转速度快,对环境条件变化颇为敏感[7]。按密度分为轻组有机碳和重组有机碳,土壤肥力的大小对土壤轻重组有机碳的影响显著,按氧化的程度可分为易氧化有机碳,能敏锐感知土壤内部有机碳的变化[8]。按有机碳的稳定性可分为活性有机碳、慢性有机碳和惰性有机碳,其特点是稳定性差,迁移能力强,可以对农田耕作措施的改变做出快速响应,因此常被用来评估农田耕作措施对土壤质量的影响。惰性有机碳相对稳定,活性低,其含量高说明有助于有机碳的积聚[9] [10] [11]。不同耕作方式的改变导致了各种有机碳组分的变化,从而影响土壤碳库的稳定性[12] [13],但是目前对于有机碳组分的转化机制及稳定性机制还缺乏相关的定论,不利于合理耕作方式的探究,当前土壤有机碳越来越受到关注,土壤有机碳作为土壤质量高低的标准,已经成为重要的参考依据之一。
当前我国大力推进各种保护性耕作措施,来改善土壤有机碳质量,旋耕是引黄罐区以往以来的传统耕作方式,由于长期以来对土壤进行扰动性耕作,不仅破坏了土壤结构,而且还加速了有机碳的损失,土壤有机碳是土壤肥力的核心,其周转过程涉及一系列物理、化学过程,是影响土壤生产力的决定性因素,使得土壤肥力下降,作物产量停滞不前。目前有关保护性耕作对有机碳的研究大多数停留在有机碳含量变化,因此本文采用免耕、旋耕和深松与秸秆还田相结合的耕作手段,通过对有机碳物理组分和化学组分指标的变化情况分析,为引黄灌区培肥地力提供理论依据和技术支撑。
2. 材料方法
2.1. 试验区概况
试验设在甘肃省农业科学院白银沿黄灌区农业试验站(37˚05'N, 104˚40'E),是由黄河水提升而形成的灌区,属于干旱区,土壤类型为灰钙土,质地为中壤。该试验开始于2017年,采用玉米连作的种植方式。试验区气候状况见表1,土壤基本理化性质见表2。
Table 1. Indicators of climatic basic conditions in the test area
表1. 试验区气候基础状况指标
气候指标 |
年平均降水量 |
蒸发量 |
年平均气温 |
>0℃积温 |
>10℃积温 |
年日照时数 |
辐射量 |
数值 |
259 mm |
2369 mm |
6.6˚C |
3208˚C |
2622˚C |
2919 h |
616.2 kJ/cm2 |
Table 2. Basic physical and chemical properties of soil in experimental area
表2. 试验区土壤基本理化性质
养分 0~20 cm |
有机质(g/kg) |
全氮 (g/kg) |
全磷 (g/kg) |
全钾 (g/kg) |
碱解氮 (mg/kg) |
速效钾 (mg/kg) |
速效磷 (mg/kg) |
pH |
含量 |
12.58 |
1.22 |
1.09 |
1.35 |
45.4 |
193 |
11.5 |
8.25 |
2.2. 试验设计
试验于2017年开始,设4个处理:(1) 旋耕(RT,机械旋耕);(2) 免耕(NT,土壤不耕作);(3) 深松(ST,全方位深松机深松35 cm);(4) 深松 + 秸秆还田(STS,全方位深松机深松35 cm,秋季将收获后的秸秆全部翻压还田),以旋耕处理为对照(CK),小区面积900 m2;每个处理重复3次,随机区组排列。旋耕采用东方红1GQN-125型旋耕机,耕作深度平均15 cm,深松采用沃野ISQ-340型全方位深松机,在秋季玉米收获后深松1次。小区面积900 m2,种植作物采用玉米连作(先玉335),密度为6.75万株/hm2,灌溉技术采用当地成熟的灌溉施肥技术,施肥量N375 kg/hm2、P2O5150 kg/hm2;40%的氮肥和全部磷肥作为基肥,于播种前结合整地施入耕层,剩余60%氮肥于玉米拔节期结合灌水追施;玉米全生育期灌水4次,灌溉定额5400 m3/hm2,灌水分配比例为出苗~拔节16%、拔节~抽雄28%、抽雄~乳熟31%、乳熟~成熟25%。不同耕作处理实施后的耙耱、镇压、开沟覆膜等措施及田间管理同当地大田。2023年4月中旬播种,10月上旬收获玉米。
2.3. 取样与测定方法
2.3.1. 采样方法
玉米收获后,在每个试验小区采用五点法取样法,分别取0~20 cm、20~40 cm,两个层次的原状土约2 kg,带回试验站,掰成1 cm的小块,风干之后去除粗根和小石块,装袋带回实验室,根据不同类型有机碳实验指标的土样要求,处理土样。
2.3.2. 测定方法
有机碳[14]:称取过0.25 mm孔筛的风干土样0.5 g (精确到0.0001 g),放入一干燥的硬质试管中,用移液管准确加入标准溶液5 mL,加H2SO4 5 mL充分摇匀,管口盖上小漏斗,以冷凝蒸出只水汽。将8~10个试管放入自动控温的铝块管座中,放入温度为185~190℃的石蜡油中,必须使温度在170℃~180℃之间,待试管液体沸腾发生气泡开始计时,煮沸5分钟,取出试管,待冷却之后,移入250 mL三角瓶中,总体积为60~70 mL,然后加入3滴指示剂,用0.2 mol/L的硫酸亚铁溶液滴定,颜色变为暗绿,即为滴定终点。
易氧化有机碳[15]:称取过0.5 mm土壤筛的含碳15 mg的风于土;将上述土放入100 mL离心管中,加入25 mL的333 mol/LKMnO4溶液,拧紧盖子,然后在200 r/pm下振荡1 h空白实验同时进行(不含土样,其它操作样):震荡后在4000 r/pm离心5 min,每一离心管中用微量移液器吸取0.4 mL滤液转入100 mL容量瓶,用蒸馏水定容至100 mL (稀释250倍):稀释液在565 nm波长处进行比色,同时标准系列浓度的KMnO4:溶液也比色,使样品浓度在标准系列范围内。
颗粒碳有机碳[16]:土壤颗粒有机碳的测定,部分采用六偏磷酸钠提取法进行提取测定。将采集土样剔除石砾与植株残茬在自然条件下风干,过2 mm孔径筛。称取待测土样10 g置于50 ml塑料瓶中,并加5 g/L的六偏膦酸钠溶30 mL,摇匀后置于摇床上振荡18 h,将分散液和土壤全部置于0.053 mm和0.25 mm筛上,清水冲至水流清晰不含细土颗粒为止。筛上大于53 um的土壤转移到事先称好质量的铝盒中,放置于60℃恒温烘箱中烘干并称重,计算烘干后土样重占整个土样的比例,即为土壤颗粒组分比例。颗粒态有机碳(0.053~0.25 mm),然后在60℃下烘干称重。测定各粒级有机碳含量,乘以各粒级所占土壤的百分比计算细颗粒有机碳的含量(g/kg全土)。矿物结合态有机碳(MAOC, <0.053 mm):在测定颗粒态有机碳的过程中,将筛下部分(<0.053 μm)的有机质在60℃条件下烘干,称重,测定有机碳含量,乘以各粒级所占土壤的百分比计算出矿物结合态有机碳含量。
轻重组有机碳[17]:采回的土样剔除肉眼可见的植物根系,在室温下全部风干后,研磨并过2 mm筛。称取25 g土壤,放入150 mL三角瓶中,加入比重1.70 g/mL的NaI溶液50 mL,在200 r/min条件下振荡1 h,将振荡后的样品置于100 mL离心管中,然后在2500 r/min离心20 min,将浮在NaI表面的轻组有机质倾倒于装有0.45 mm纤维滤膜的容器过滤器中进行抽滤,依次用总量100 mL 0.01 mol/L CaCl2和100 mL蒸馏水冲洗纤维滤膜上轻组有机质至少3次。再向离心管中加入50 mL比重为1.70 g/mL的NaI溶液进行第2次浸提,步骤同上。将两次浸提的轻组有机质转移至烧杯,在60℃下烘干、称重。烘干的轻组有机质用玛瑙研钵研细,过0.25 mm筛。离心管中的剩余土为重组,土样依次用0.01 mol/L CaCl2洗净土样中的NaI,再用蒸馏水冲洗土样若干次,洗净土样中的CaCl2。在60℃下烘干、称重,将烘干的重组土样研细,用于测定重组有机碳。
惰性碳:上述剩余残渣用去离子水洗涤土壤成近中性,烘干、研磨、过0.15 mm筛子,用重铬酸钾外加热法测定土壤中的碳。
实验仪器:分光光度计、离心机、摇床尼龙筛子油浴锅烘箱比重计离心机,真空泵,容器过滤器,油浴锅,滴定仪,80/100消煮管、TOC仪。
2.4. 计算方法
(1)
式中:c——0.8000 mol∙L−1 (1/6K2Cr2O7)标准溶液的浓度;
5——重铬酸钾标准溶液加入的体积(mL);
V0——空白滴定用去FeSO4体积(mL);
V——样品滴定用去FeSO4体积(mL);
3.0——碳原子的摩尔质量(g·mol·L−1);
10−3——将mL换算为L;1.1为氧化校正系数;
m——烘干土样质量(g);
k——将风干土换算成烘干土的系数(即为1-W%,W%:风干土样含水量)。
(2)
式中:空白、样品——机器测定的读数mM;
25——待测液总体积ml;
9——氧化过程中1 mol高锰酸钾消耗9 mg碳;
M——称取样品的质量g。
(3)
式中:WPOC——颗粒有机物所占土壤质量的比例(%);
M0——待测土样的重量(g);
M1——过筛后的烘干土壤重量(g)。
(4)
式中,Ct-POC——供试土壤中颗粒态有机碳的含量(g·kg−1);
CPOC——土壤颗粒组分中有机碳的含量(g/kg)。
(5)
式中,WMAOC——矿物结合有机碳所占土壤质量的比例(%)。
(6)
式中,Ct-MAOC——供试土壤中矿物结合有机碳的含量(g/kg);
CMAOC——土壤矿物结合有机碳组分中有机碳的含量(g/kg)。
(7)
(8)
(9)
(10)
2.5. 数据处理
采用Excel 2019、IBM SPSS Statistics 23.0和Origin 2021软件进行试验数据处理和统计分析。采用单因素分析不同耕作处理之间的差异显著性,显著性水平p < 0.05,图表中数据为平均值。
3. 结果分析
3.1. 不同耕作方式对灌溉农田土壤总有机碳的影响
由图1可以看出,在0~40 cm土层,总有机碳含量由大到小依次为STS > NT > ST > RT,随着土层的加深,总有机碳含量表现为减小的趋势,处理RT、NT、ST、STS在0~20 cm土层的总有机碳含量较20~40 cm土层的有机碳,分别减少了20%、14%、5%、11%,各个处理的总有机碳含量的降低幅度由大到小依次为RT > NT > STS > ST。从不同处理来看,在0~20 cm土层,总有机碳含量由大到小依次为STS > NT > ST > RT,其中处理STS、ST较对照差异显著,分别提高了23%、6%,处理NT较对照RT差异不显著。在20~40 cm土层,总有机碳含量由大到小依次表现为STS > ST > NT > RT,其中处理STS较对照RT差异显著,提高了33%,处理ST、NT较对照RT差异不显著。
注:不同小写字母表示相同土层数据间的差异显著性(P < 0.05),下同。
Figure 1. Total organic carbon content of irrigated farmland soils under different tillage methods (g/kg)
图1. 不同耕作方式下灌溉农田土壤总有机碳含量(g/kg)
3.2. 不同耕作方式对土壤有机碳物理分组的影响
3.2.1. 不同耕作方式对轻、重组有机碳含量的的影响
由图2可以看出,在0~40 cm土层,STS处理的轻组有机碳含量较其它三个处理呈现出较大差异,随着土层的加深呈现出轻组有机碳含量减小的趋势,处理RT、NT、ST、STS的轻组有机碳含量在0~20 cm土层较20~40 cm土层分别降低了46%、32%、90%、92%,轻组有机碳含量降低的幅度由大到小依次为STS > ST > RT > NT,免耕处理的降低幅度最小,可能是因为免耕处理未受到剧烈的扰动,从而导致降低幅度小。从不同土层来看,在0~20 cm土层,轻组有机碳含量表现为STS > ST > NT > RT,其中处理STS、ST较RT差异显著,分别提高了117%、54%,秸秆还田之后对轻组有机碳含量的提高幅度较大。处理NT较RT差异不显著。在20~40 cm土层,各处理轻组有机碳含量表现为STS > NT > ST > RT,只有处理STS较RT差异显著,处理ST、NT较RT差异不显著。由图3可以看出,在0~40 cm土层,随土层的加深,有机碳含量表现为降低,处理RT、NT、ST、STS在20~40 cm土层的重组有机碳含量较0~20 cm土层有机碳含量分别降低了26%、32%、37%、35%,降幅由大到小依次为ST > STS > NT > RT。从不同处理来看,在0~20 cm土层,处理STS、NT、ST的重组有机碳含量较RT差异显著,分别增加了37%、24%、13%。在20~40 cm土层,处理STS、ST的重组有机碳含量较RT差异显著,分别提高了29%、15%。处理NT的重组有机碳含量较RT差异不显著。结合图3、图4综合来看,轻组有机碳和重组有机碳含量表现为较大的差异,重组有机碳的含量远远高于轻组有机碳的含量。
Figure 2. Organic carbon content of light groups in irrigated farmland under different tillage methods
图2. 不同耕作方式下灌溉农田轻组有机碳含量(g/kg)
Figure 3. The content of heavy fraction organic carbon in irrigated farmland under different tillage methods (g/kg)
图3. 不同耕作方式下灌溉农田重组有机碳含量(g/kg)
3.2.2. 不同耕作方式对颗粒态、矿质结合态有机碳的影响
从图4可以看出,随着土层的加深,颗粒态有机碳随着土层的加深,呈减小的趋势,土层加深导致植物残枝落叶减少,并且随着深度的加深土壤温度降低也会影响土壤颗粒态有机碳。处理RT、NT、ST、STS随土层增加分别减小了19%、88%、9%、61%,变化幅度由大到小的顺序为NT > STS > RT > ST。从不同处理来看,在0~20 cm土层,各处理的颗粒态有机碳含量由大到小依次为STS > NT > ST > RT,STS、ST、NT较RT差异显著,分别提高了132%、52%、82%,ST处理提高的最为明显。在20~40 cm土层,各处理的颗粒态有机碳含量由大到小依次为STS > ST > NT > RT,ST、STS较RT差异显著,分别提高了70%、67%,NT较RT差异不显著。不同层次之间,均是对照RT含量最小。耕作模式在0~20 cm土层对颗粒态有机碳的影响较20~40 cm土层强烈。
Figure 4. Particulate organic carbon content (g/kg) in irrigated farmland under different tillage methods
图4. 不同耕作方式下灌溉农田颗粒态有机碳含量(g/kg)
Figure 5. Mineral bonded organic carbon content (g/kg) in irrigated farmland under different tillage methods
图5. 不同耕作方式下灌溉农田矿质结合态有机碳含量(g/kg)
由图5可以看出,随着土层的加深,矿质结合态有机碳含量呈减少趋势,土层加深土壤中的矿物质含量减少,从而使得土壤矿质结合态有机碳含量降低。处理RT、NT、ST、STS随着土层增加,分别较少了42%、26%、71%、26%,减小幅度由大到小依次为ST > RT > NT = STS。从不同层次来看,在0~20 cm土层,STS、ST较RT差异显著,分别提高了38%/29%,处理NT较RT差异不显著。在20~40 cm土层,处理STS较RT差异限制,提高了57,处理ST、NT较RT差异不显著,随着土层加深,深松处理对矿物质结合态有机碳的影响有所减缓,秸秆还田对有机质与矿物质的结合有促进作用。
3.3. 不同耕作方式对农田土壤有机碳化学组分影响
3.3.1. 不同耕作方式对惰性有机碳含量的影响
由图6可以看出,从不同土层来看,随着土层的增加,各个处理的惰性有机碳含量均表现为下降,在0~40 cm土层,各处理的惰性有机碳含量由大到小依次为STS > ST > NT > RT,处理RT、NT、ST、STS在0~20 cm土层较20~40 cm土层分别增加了24%、30%、17%、16%,降低幅度依次为NT > RT > ST > STS,免耕处理的惰性有机碳随着土层的加深,减小的幅度最为明显,而深松加秸秆还田的方式减小幅度最小。从不同处理来看,在0~20 cm土层,各处理的惰性有机碳含量由大到小依次为STS > NT > ST > RT,处理STS的惰性有机碳含量较对照RT差异显著,提高了20%,处理ST、NT的惰性有机碳含量较RT差异不显著,在20~40 cm土层,惰性有机碳含量由大到小依次表现为STS > ST > NT > RT,处理STS的惰性有机碳含量较RT差异显著,提高了30%,处理NT、ST的惰性有机碳含量较RT差异不显著。
Figure 6. Inert organic carbon content (g/kg) in irrigated farmland under different tillage methods
图6. 不同耕作方式下灌溉农田惰性有机碳含量(g/kg)
3.3.2. 不同耕作方式对易氧化有机碳含量的影响
从图7可以看出,随着土层的加深,易氧化有机碳含量呈现减小的趋势,NT两个层次间差异较大,20~40 cm土层相较于0~20 cm土层减少了45%。RT、ST、STS分别减小了21%、7%、8%,减小幅度大小为NT > RT > ST > STS。在0~20 cm土层,含量大小顺序为STS > ST > NT > RT。处理NT、ST、STS较RT差异均显著,易氧化有机碳含量分别提高了24%、35%、49%,在20~40 cm土层,含量大小顺序为STS > ST > NT > RT。处理STS、ST较RT差异显著,易氧化有机碳含量分别提高了97%、68%、处理NT较RT差异不显著。
Figure 7. Content of readily oxidized organic carbon in irrigated soil under different tillage methods (g/kg)
图7. 不同耕作方式下灌溉农田土壤易氧化有机碳含量(g/kg)
3.4. 不同耕作方式下各组分之间的相关性分析
注:SOC表示总有机碳,ROOC表示易氧化有机碳,LFOC表示轻组有机碳,HFOC表示重组有机碳,POC表示颗粒态有机碳,MOC表示矿质结合态有机碳,RPC表示惰性有机碳,*表示在0.05水平上显著相关,**表示在0.01水平上显著相关,右向红色椭圆表示正相关,椭圆大小表示相关系数大小,下同。
Figure 8. Pearson correlation between 0~20 cm soil total organic carbon and different organic carbon components
图8. 0~20 cm土壤总有机碳和不同有机碳组分之间的Pearson相关性
将土壤总有机碳含量与土壤各有机碳组分含量进行Pearson相关性分析发现(图8、图9),在0~40 cm土层,SOC与ROOC、MOOC、差异不显著。随着土层加深ROOC与LFOC不相关,而在0~20 cm呈显著正相关(p < 0.01),LFOC与POC在0~20 cm呈显著正相关(p < 0.05),20~40 cm不相关。在20~40 cm MOC与RPC不相关,但在20~40 cm呈显著正相关(p < 0.05),以上结果反映出随着土层的加深相同组分之间存在着差异性。
Figure 9. Pearson correlation between 20~40 cm soil total organic carbon and different organic carbon components
图9. 20~40 cm土壤总有机碳和不同有机碳组分之间的Pearson相关性
4. 讨论
土壤有机碳是土壤肥力的重要指标[18],其含量和组分随着耕作等管理措施的变化而变化,而有机碳根据其性质又可以分为多种类型,因而对土壤有机碳各组分变化的探究,有助于制定土壤管理策略[19],土壤中的活性有机碳,惰性有机碳,重组有机碳主要来源是植物通过一系列的生物化学过程中分解形成的有机质,可以在土壤中长期保留,表征土壤有机碳的稳定性[20]。本研究通过对各个组分的分析研究,得出来不同组分的有机碳在不同耕作方式下的变化情况,其中易氧化碳作为活性碳,是有机碳的重要组成部分,可以在土壤有机碳发生变化前反映土壤的微小变化,并且直接参与土壤生物化学转化过程,其有效性高,并且容易被微生物利用和分解,其数量受耕作方式等其它因素的影响较为显著[21]。从本实验结果来看,随着土层的加深,易氧化有机碳的含量表现为随土层的加深含量减少,其中秸秆还田能明显增加0~40 cm土层的易氧化有机碳含量,除此之外,免耕也能提高易氧化有机碳含量,这与周欢等[22]人的研究结果一致,其原因可能是免耕能增加其通透性,而秸秆的输入使得微生物的物料充足。土壤颗粒态有机碳是由未完全分解的植物残体以及根系残体组成,对于玉米作物而言,表层土壤受其影响较为敏感,并且周转速度较快,能敏感的反映作物残体归还量和有机碳的周转情况,本研究中,颗粒碳有机碳随着土壤深度的增加,颗粒态有机碳含量减少,在0~20 cm土层,处理深松加秸秆还田较其它三个处理显著的增加了颗粒碳有机碳的含量,这与于建光等[23]人的研究结果一致,然而免耕处理的颗粒态有机碳含量高于深松处理,可能是深松处理连年耕作在一定程度上损失的根系有机质,而免耕处理根系有较多的根系有机质。土壤中的矿物态有机质被认为是有机碳固持的重要机制之一,这部分有机碳其含量相对较高,属于惰性碳库,李欢等[24]研究结果表明,添加秸秆物质,可显著增加稻田土壤的矿质态有机碳的含量,在本研究中,在0~20 cm土层,秸秆还田处理矿质结合态有机碳含量最高,在20~40 cm土层,深松加秸秆还田的耕作方式使得土壤中的矿质结合态有机碳含量显著高于其他处理,随着土层的加深,秸秆还田增加了深层土壤的矿质结合态有机碳含量。轻组有机碳是属于土壤碳库的活性部分,轻组有机碳对耕作方式,秸秆还田等人类活动较为敏感,是衡量土壤碳库质量的重要指标,在本研究中,深松处理和深松加秸秆还田处理免耕处理均在一定程度上增加了轻组有机碳的含量,在0~20 cm土层,处理深松和深松加秸秆还田处理较对照差异显著,在20~40 cm土层,只有深松加秸秆还田较对照差异显著,秸秆还田提高了整体的轻组有机碳含量,而深松可以提高表层的含量,这与胡月华等人[25]的研究结果一致。而土壤重组有机碳是土壤碳库中较为稳定的有机碳,可以指示评价土壤有机碳的固存性能,在本研究中,就0~40 cm土壤整体而言,重组有机碳含量大小相差较小,重组有机碳含量表现为表层高,相较于轻组有机碳含量远大于轻组有机碳含量,深松加秸秆还田显著提高两个层次的有机碳含量,深松提高了0~20 cm土层的重组有机碳含量,这与高雪莹等人[26]的研究结果相似。土壤有机碳可以被划分为活性碳库和惰性碳库,在对有机碳的分组研究基础之上,在本研究中,在0~40 cm土层,深松加秸秆还田处理重组有机碳含量较对照差异显著,其它处理之间差异均不显著。
综上所述,通过进一步比较不同耕作方式下有机碳各组分的含量发现,土壤易氧化有机碳、轻重组有机碳、颗粒、矿物结合态有机碳以及活性碳库、惰性碳库的含量以及变化和分布规律基本趋于一致,秸秆还田的保护性耕作方式显著提高其含量。其次通过对两个层次的各组分碳含量分析可以得出,各处理均有表层富集的现象。所以对于引黄灌区的灌溉农田,采用秸秆全量还田且深松的耕作方式有利于碳的固持。
5. 结论
(1) 除慢性有机碳含量外,保护性耕作较传统耕作均提高了有机碳及有机碳各组分的含量,其中秸秆还田提高最显著。秸秆还田能显著提升土壤固碳能力。
(2) 总有机碳含量和各有机碳组分含量之间存在相关性,可见土壤有机碳及各组分之间相互促进和转化。随着土层的加深,各组分含量呈现出降低的趋势,旋耕和免耕的降低幅度小。
NOTES
*通讯作者。