农膜残留对大豆光生理特征及生物量累积的影响*

黄 珊,樊廷录,刘萌娟,陈荣桓,梁楚涛,程万莉,陈延华,薛 萐,5,杨晓梅,5**

(1.西北农林科技大学资源环境学院 杨凌 712100;2.甘肃省农业科学院旱地农业研究所 兰州 730070;3.西北农林科技大学农学院 杨凌 712100;4.北京市农林科学院植物营养与资源研究所 北京 100007;5.西北农林科技大学水土保持研究所/黄土高原土壤侵蚀与旱地农业国家重点实验室 杨凌 712100)

作为干旱半干旱地区重要的农业种植技术,农膜施用在保持土壤水分、增加土壤温度、促进种子萌发和作物生长及提高产量方面发挥着巨大的作用[1]。然而,作物收割所造成的农膜碎片化,给农膜回收及捡拾造成了极大的困难,从而累积在土壤中的农膜逐年增加[2-5]。据估计,欧洲每年输入农田土壤的塑料颗粒为6.3～43.0万t,北美农田塑料输入量约4.4～30.0万t,已严重威胁土壤质量[6]。我国农田土壤的塑料残留量为0.1～324.5 kg·hm−2,平均83.6 kg·hm−2[7];同时随着长期的施用及大面积推广,农膜残留给农田土壤及周围生态环境造成了严重的危害[8]。

残留于土壤中的农膜,一方面在耕作机械力的作用下,其破碎程度增加,并可与土壤颗粒聚合,从而可能在土体内传输;另一方面,受光照、水分以及土壤微生物的作用,进一步老化、降解,从而对植物生长及土壤环境产生影响[9-10]。研究表明:在农膜添加条件下,小麦(Triticum aestivum)的叶和根生物量积累均受到抑制[11];5～20 a农膜残留量的累积增加导致玉米(Zea mays)出苗率下降[12]。同时,由于农膜残留的积累,干扰作物的生长发育,尤其是根系发育[13]。有研究表明,随着土壤中农膜累积,马铃薯(Solanum tuberosum)根系受到严重影响,在农膜残留量为180 kg·hm−2时,根系体积减少[14];而土壤中残留的农膜不仅影响番茄(Solanum lycopersicum)幼苗根系的生长及分布,且会降低根系活力[15];当残膜量达到900 kg·hm−2时,棉花(Gossypiumspp.)根系向下伸长缓慢,长度最小[16]。农膜残留对根生长的危害也影响植物对养分的吸收与积累,进而影响植物的生长发育[17]。土壤中因农膜的存在导致马铃薯减产显著,当农膜残留量为720 kg·hm−2时,其减产幅度最大,与对照相比减产9.7%～14.5%[14]。随着土壤中农膜残留量的增加,番茄幼苗的株高、茎粗、地上部及根系鲜重逐渐降低[15],从而对农作物的生长及产量产生了明显的不良影响。此外,尽管根系是植物从土壤中获取养分的主要途径,其吸收的矿物质占总干物重的8%～9%[18],而光合作用是植物获取能量和积累养分的主要来源,光合产物累积量占总干物重的91%～92%[19-20]。因此,土壤农膜残留对农作物生长及产量的影响不能仅仅考虑萌发、根系生长及产量状况,还可能对植物光合作用特征产生影响。

光是光合作用的主导因子,植物光合作用对光的响应体现了植物净光合速率和光合有效辐射之间的关系,对研究植物光化学过程中的光化学效率具有重要的意义[21]。而CO2作为植物光合作用的原料,其浓度及其利用效率影响着植物的光合作用过程,进而影响植物的光合速率。研究发现,农作物生长情况以及干物质量积累可以通过光合作用的能力来体现,即光响应曲线[22-23]和CO2响应曲线[24]拟合参数来表征该植物的光合作用能力。叶子飘等[25]通过建立直角双曲线修正模型,将所得的光响应与CO2响应曲线拟合值与实际测量值进行对比,其所得拟合参数与实际情况更贴近[26-27],从而被用来研究植物光合作用的生理特征。尽管农膜残留对产量的研究较多[28],但农膜残留累积是否会影响作物的光合作用,从而影响其生物量累积等方面的研究还尚未清楚。基于此,本研究通过室外盆栽控制试验,研究不同类型塑料及残留累积量下,大豆(Glycine max)花期及初荚期叶片光合作用对光及CO2响应曲线以及地上、地下的生物量特征,明确塑料类型及残留量对大豆生物量累积及叶片光生理响应特征的影响,以期为农膜施用及生物可降解膜的生态安全评价提供依据。

1 材料与方法

1.1 试验材料及设计

本试验依托土壤侵蚀与旱地农业国家重点实验室,在干旱模拟试验基地(中国科学院水土保持研究所)进行,其位于温带半湿润半干旱气候区,具有显著的大陆性季风气候特征:年均气温12.9 ℃,无霜期211 d,年均日照时数2163.8 h,年总太阳辐射量481.26 kJ·cm−2,年均降水量635.1 mm,多集中在6—9月。试验用土采自于杨凌周边农田无农膜施用史的耕层土壤,去除杂物及砾石后,过筛(10 mm)并风干待用。农膜选取当地农户常用的聚乙烯薄膜(PE处理)及生物降解膜(98%聚己二酸丁二酯+聚乳酸,BP处理),根据田间农膜残片大小,剪裁为2 cm×2 cm、1 cm×1 cm、0.5 cm×0.5 cm的碎片并按(1∶1∶1)均匀混合,按照不同质量浓度梯度(0、0.1%、0.5%、1.0%)加入试验用土中,并反复搅拌使其充分混合均匀。大豆种子由西北农林科技大学农学院大豆育种实验室提供,品种为‘种群290’,选择饱满且大小均匀的种子(单粒重:0.23±0.02 g),进行盆栽种植(花盆高26.5 cm、底部直径20 cm、顶部直径24 cm)。种植时,每盆加入6.5 kg混有农膜残片的试验土,灌水使土壤水分含量约20%,静置隔夜后播种。每盆播种5粒大豆种子,然后在表土上添加60%含水量的含包含蛭石、珍珠岩、椰丝的基质500 g(约1 cm厚)。每个处理重复3次,且日常观测水分含量,维持正常的田间持水量60%～80%,待出苗完成1个月后每盆定株2株。在花期及豆荚初期,对植株进行光合作用测定;收集整个植株,60 ℃烘干至恒重,测定花期及收获期地上及地下生物量。

1.2 光生理指标及生物量测定

分别选择大豆花期(6月1—14日)与初荚期(7月10—27日)的晴天,于9:00—12:00和14:00—17:00,利用LI-6400XT型便携式光合测定仪测定大豆叶片的瞬时光合速率。

在大气CO2浓度条件下(约400 μmol·mol−1),仪器配备红、蓝人工光源,设定光合有效辐射强度(PAR,μmol·m−2·s−1)由高到低依次为2000、1500、1000、800、500、250、120、60、30、15、0。被测叶片在每个PAR下适应120～180 s待仪器稳定后开始测定。测定时叶片位置保持在植株的原方位。通过直角双曲线修正模型进行光响应曲线拟合,并用光滑曲线连接,获得大豆的光补偿点(LCP)、光饱和点(LSP)、初始量子效率(AQY)、暗呼吸速率(Rd)及最大光合速率(Pmax)。CO2响应作用进行时,由于生成的反应底物具有一定的滞后性,所以在观测时仪器控制的CO2浓度先从大气浓度降至50 μmol·mol−1,再从大气浓度升至2000 μmol·mol−1,设定CO2浓度(μmol·mol−1)梯度为400、300、200、100、50、400、400、800、1200、1600、2000。样品室固定光强为1200 μmol·m−2·s−1,光诱导持续约15 min,待仪器稳定后开始测定。测定时叶片位置保持在植株的原方位。每个CO2浓度的观测时长为120～180 s。通过直角双曲线修正模型对测得的数据进行CO2响应曲线拟合,并用光滑曲线连接,获得CO2补偿点(CCP)、CO2饱和点(CSP)、初始羧化效率(η)、光呼吸速率(Rp)及最大光合能力(Amax)。

生物量测定采用烘干称重法,精确度为0.01 g。在大豆花期以及成熟期采样,带回室内洗净晾干,分别测定地上部分和地下部分的鲜重;然后105 ℃下烘30 min,将温度下降到65 ℃,继续烘干至恒重测定干重。

1.3 数据分析

根据叶子飘等[29]建立的直角双曲线修正模型,大豆叶片光合作用对光响应模型为:

式中:Pn(I)为净光合速率;α为光合作用光响应曲线的初始斜率,即初始量子效率(AQY);β和γ为系数;I为光合有效辐射;Rd为暗呼吸。

光饱和点(LSP)表达式为:

最大净光合速率(Pmax)表达式为:

通过拟合光合作用对光响应曲线,得到大豆叶片光补偿点(LCP)。

根据叶子飘等[30]建立的直角双曲线修正模型,大豆叶片光合作用对CO2响应模型为:

式中:An(Ci)为净光合速率,Ci为胞间CO2浓度,Rp为光呼吸速率,a为CO2响应曲线的初始羧化效率(η),b和c为系数。

CO2饱和点(CSP)表达式为:

最大光合能力(Amax)表达式为:

通过拟合光合作用对CO2响应曲线,得到大豆叶片CO2补偿点(CCP)。

处理间方差分析以及相关性分析采用SPSS 23统计分析软件进行(P＜0.05),并利用Origin 2019作图。

2 结果与分析

2.1 农膜残留对大豆叶片光合作用光响应曲线及其拟合参数的影响

聚乙烯薄膜(PE)和生物降解膜(BP)不同残留条件下大豆叶片光合作用光响应曲线拟合结果如图1所示。不同农膜残留处理下大豆叶片在PAR低于500 μmol·m−2·s−1时,净光合速率(Pn)随PAR增加急剧增加;PAR在500～1000 μmol·m−2·s−1时,Pn的增加速度减缓;PAR＞1000 μmol·m−2·s−1后,Pn基本保持平稳,达到光饱和。不同处理的拟合效果均较好(R2=0.99)。

基于光合作用光响应曲线拟合特征,大豆叶片光响应曲线拟合参数见表1。PE处理中,花期与初荚期的AQY、Pmax(除花期0.1%处理外)、LSP随农膜残留量的增加未发生显著变化;而花期LCP值随农膜残留量表现为:0.1%＞CK＞0.5%＞1.0%,且0.1%残留量显著高于0.5%、1.0%处理;Rd在残留量为0.1%处理最大,但与其他处理间无显著差异。初荚期,LCP在0.5%和1.0%农膜残留量处理下显著高于CK和0.1%农膜残留量;Rd在农膜残留量为1.0%时最大,且显著高于CK和0.1%农膜残留量处理。BP处理中,花期残留量0.1%时,AQY与Pmax最大,AQY显著高于CK和0.5%残留量处理,Pmax则显著高于其他3个处理;LSP在各残留量处理间无显著差异;而LCP在0.1%残留量时最低,显著低于0.5%和1.0%残留量,且在残留量为0.5%时达峰值。Rd在0.1%残留量处理下最高,显著高于CK和0.5%残留量处理。初荚期,BP残留对AQY无明显影响;Pmax在残留量0.1%时显著低于其他处理;LSP在残留量为0.1%时最高;LCP在残留量为0.1%和0.5%时较低,而LCP1.0%残留量时与CK无显著差异;Rd各残留量间无显著差异。

对比PE、BP两种农膜,在相同残留量条件下,AQY差异不显著;花期Pmax在0.1%残留量时BP显著高于PE,而初荚期PE显著高于BP,而在其他浓度处理中无显著差异;LSP在花期无显著差异,初荚期各残留量均为BP显著高于PE;LCP在花期,仅在0.1%残留量时PE明显大于BP,而在其他浓度时BP均大于PE;Rd花期在残留量为1.0%时BP显著高于PE。

表1 聚乙烯薄膜(PE)和生物降解膜(BP)不同残留条件下大豆叶片光合作用光响应曲线拟合参数Table 1 Parameters of light response curves of soybean photosynthesis under different treatments of polyethylene film(PE)and biodegradable plastic film(BP)residues in soils

2.2 农膜残留对大豆叶片光合作用CO2响应曲线及其拟合参数的影响

聚乙烯薄膜(PE)和生物降解膜(BP)不同残留条件下大豆光合作用CO2响应曲线拟合结果如图2所示。不同类型农膜残留下,大豆叶片在胞间CO2浓度＜200 μmol·mol−1时,Pn随胞间CO2浓度的增加急剧增加,胞间CO2浓度在200～800 μmol·mol−1时,Pn随胞间CO2浓度的增加而增加速度减缓,胞间CO2浓度＞800 μmol·mol−1时,Pn基本保持平稳,达到饱和,且拟合效果较好(R2=0.90)。

基于光合作用CO2曲线拟合特征,大豆叶片CO2响应曲线拟合参数见表2。结果表明:PE处理中,在花期,η、Amax、CSP、CCP、Rp均在各农膜残留量之间无显著差异。初荚期η值表现为:CK＞0.1%＞0.5%＞1.0%,0.5%和1.0%残留量显著低于CK;Amax在残留量0.5%时最低,且与其他残留量间和CK间差异显著;CSP各农膜残留量间无显著差异;CCP随残留量提高而降低,表现为:CK＞0.1%＞0.5%＞1.0%,且0.5%和1.0%残留量处理显著低于CK和0.1%残留量处理,0.1%残留量处理与CK间无显著差异;Rp值在农膜残留量为0.1%时最大,并与其余残留量和CK间差异显著。

BP处理中,花期η在残留量为0.1%时最大,与其他残留量间和CK间差异显著;农膜残留量对Amax和CCP无显著影响;CSP在残留量为0.1%时最大,并与0.5%残留量间差异显著,而其他各残留量间则无显著性差异;Rp在农膜残留量为0.1%时显著最大,而后随着残留量的增加不断减小。初荚期,η、Amax的值随农膜残留量的变化为:CK＞0.1%＞0.5%＞1.0%;不同处理间CSP无显著差异;CCP在低浓度处理中逐渐降低,而在残留量为1.0%时,又恢复到与CK同一水平,且无显著差异;Rp在各残留量处理间无明显差异。

对比PE、BP两种农膜残留,η仅在花期残留量为0.1%时,BP显著高于PE。Amax在花期残留量为1.0%时PE显著高于BP;而在初荚期残留量为0.5%时,BP显著高于PE,而在残留量为1.0%时,PE显著高于BP。CCP仅在初荚期残留量为1.0%时,BP显著大于PE。Rp在花期残留量0.1%时,BP显著大于PE;而在初荚期残留量为0.1%时,PE显著大于BP。

表2 聚乙烯薄膜(PE)和生物降解膜(BP)不同残留条件下大豆叶片光合作用CO2响应曲线拟合参数Table 2 Parameters of CO2 response curves of soybean photosynthesis under different treatments of polyethylene film(PE)and biodegradable plastic film(BP)residues in soils

2.3 农膜残留对大豆生物量的影响

大豆花期及收获期生物量的测定结果(图3)表明:大豆地上、地下生物量在不同生长阶段及处理间的差异不同。总体来说,随土壤中农膜残留量增加,大豆地上及地下生物量均有所降低。特别是在花期,不同处理间地上、地下生物量差异较大;而在收获期,生物量(地上/地下)差异不明显。不同类型塑料残留处理中,农膜残留量为0.1%和0.5%的BP处理花期地下生物量显著降低,但当农膜残留量1.0%时,生物量与对照并无显著差异。而在PE处理中,其花期地上生物量在残留量为0.5%时达最大值,但不同残留量之间差异不显著。而在收获期,随着残留浓度的增加,地上和地下生物量逐渐减小,但不同处理间差异不显著。

2.4 大豆生物量与光合作用光响应和CO2响应参数的相关性分析

将大豆光合作用光响应与CO2响应曲线拟合参数与生物量进行相关性分析(表3),结果表明各参数与生物量相关性明显不同。生物量与光合作用光响应拟合参数结果表明:地上、地下生物量与各参数的正负相关性不同。AQY与地上生物量呈正相关关系,而与地下生物量呈负相关关系,且塑料类型之间并无差异(收获期PE地下生物量除外)。而在不同处理的花期与收获期,Pmax与生物量关系变化较大,但均相关不显著。对于生物量与LCP分析发现,PE处理均呈负相关关系;且在收获期与地上生物量相关性达到极显著;而在BP处理中,在不同时期,LCP与地上、地下生物量之间的正负关系变化较大。不同处理间生物量与光合作用CO2响应参数Amax呈正相关关系,η、CSP、CCP与生物量的相关关系差异显著,但是不同生育期及地上、地下生物量相关性不同。BP处理下,大豆地下生物量花期与Rp呈显著负相关,而收获期与η呈正相关。

表3 聚乙烯薄膜(PE)和生物降解膜(BP)残留下大豆地上部分和地下部分生物量与光合作物光响应和CO2响应曲线参数的相关性Table 3 Pearson correlation analysis between soybean biomass(aboveground and underground)and photosynthesis curve fitting parameters under different treatments of polyethylene film(PE)and biodegradable plastic film(BP)residues in soils

3 讨论与分析

3.1 土壤中农膜残留特征对光响应曲线特征的影响

环境是影响植物光合作用首要因素,而生长环境的改变,尤其是土壤环境的改变,将对植物光合作用产生较大的影响[31-32]。随着农膜的长期施用,其残留累积量逐年增加,直接影响土壤性质,从而影响农作物生长及生物量累积。光合作用光响应曲线是植物在不同光强下光合作用的集中反映,光响应曲线拟合所得光合参数可反映植物光合能力与生长状况,同时可判断植物抗逆性[29-30]。本研究结果表明:不同类型塑料不同残留量下光合作用光响应曲线拟合效果较好,决定系数均高于0.99,不同生育期表现出相对一致的特征;且各处理的拟合曲线相关参数差异显著。其中,AQY的变化范围为0.043～0.095,符合植物自然条件下的阈值[33]。AQY值越大表明植物利用弱光的能力越强[34]。花期BP低残留量(0.1%)AQY值显著高于其他处理,然而当期残留量达到1.0%时AQY值与其他处理间并无显著差异;而在初荚期,BP各处理的AQY值均大于PE处理。这说明生物降解膜的残留会影响植物对弱光的利用能力,可能具有增强大豆对弱光利用能力的作用。同时,LCP反映植物叶片光合作用过程中光合同化作用与呼吸消耗相当时的光强,其值越低,意味着植物利用弱光的能力越强[35]。本研究中PE残留花期LCP降低而初荚期LCP升高,说明其增强了花期大豆叶片对于弱光利用能力,而减弱了初荚期叶片对于弱光的利用能力,而BP残留则完全相反,这可能与降解膜的生物降解有关。而LSP和Pmax体现出植物在强光下的生长状况[36],植物拥有较高的LSP和Pmax时植物在强光下生长不易受到抑制。而本研究中不同农膜残留条件下,大豆叶片的Pmax并无明显差异;随着BP残留量的增加,初荚期叶片的LSP增大,增强了叶片对强光的适应能力。因此,在初荚期BP残留增大了LSP,同时降低了LCP,从而使得大豆叶片光合作用可利用的光强范围增大[37],光的生态适应能力更强,这与林萌萌等[38]、张占琴等[39]研究结果一致。此外,Rd反映植物在黑暗条件下的呼吸速率,是在无光照条件下植物呼吸释放的CO2的量,呼吸作用分解植物积累的淀粉以及糖类,Rd越大表明在不进行光合作用时植物消耗自身养分越多[34],而Rd值越小越有利于植物在低Pn条件下保持碳平衡[31],从而应对干旱的胁迫。本研究拟合结果表明,随着农膜残留量的增加,大豆叶片自身消耗营养成分增加,且在BP处理的消耗量略高于PE处理,这可能与BP残留降解使土壤特征发生变化有关,特别是其降解过程中刺激了土壤微生物的活性及其释放的有机物质,从而使得其消耗量大于普通PE残留。由此可见,不同类型农膜及其在土壤中的残留量可能直接影响着植物的光合作用过程,特别是生物降解膜,其对光响应特征的影响及其生理机制还应进一步探讨。

3.2 土壤中农膜残留特征对CO2响应曲线特征的影响

光合作用CO2响应曲线反映植物利用不同浓度CO2作为原料进行光合作用的能力,光合作用CO2响应曲线拟合光合参数可用来分析植物的光合能力[40]。其中,η反映植物对CO2的利用能力[22],η值越高,植物对于CO2的利用能力越强。本研究中,花期各处理间η的变化并不显著;而初荚期,PE和BP残留条件下η值相比花期均略有降低。一般来说,进入初荚期,叶片活力达到顶峰[41],光合能力最强,然而由于农膜残留,直接导致大豆叶片η不升反降,从而说明农膜残留抑制了光合速率。而Amax反映植物叶片光合电子传递能力、磷酸化的活性以及对高浓度CO2的利用能力[24]。而本研究中各处理与对照相比,大豆叶片的Amax均有所降低,但其CSP无显著性变化,从而可排除气孔因素限制,由此表明:农膜残留导致叶肉细胞光合活性变低,即光合电子传递能力及磷酸化活性减弱,从而降低了叶片Amax和叶片对高浓度CO2的利用能力[42]。此外,CCP是植物叶片的光合同化作用与呼吸消耗相同时CO2的浓度,反映了植物对低浓度CO2利用的能力[43-44],CCP越高,对低浓度CO2利用的能力越小,而大豆属于C3植物,其CCP一般大于30 μmol·mol−1[24]。本试验结果表明,大豆的CCP范围为38～120 μmol·mol−1,与上述研究结果一致。初荚期叶片CCP随PE农膜残留量的增加而降低,而BP处理中CCP在残留量为0.5%时显著降低,说明农膜残留对叶片利用低浓度CO2的能力有提升作用,且PE＞BP,这可能因为大豆叶片羧化酶的活性增强,促使大豆叶片光合作用过程做出调整来适应土壤中农膜残留的累积。此外,随农膜残留量增加和大豆生长,Rp有所降低,但差异不显著,这表明在有光的条件下,土壤农膜残留对于大豆生长过程中的干物质消耗起到抑制作用,从而有利于大豆干物质的累积,这与低Rp的作物品种具有产量高这一特点相一致[24]。综上所述,土壤中塑料残留特征对光合作用CO2响应曲线各拟合参数的影响较大,特别是生物可降解塑料及其残留量的特征,直接影响着大豆干物质累积及生物量,故后续在可降解塑料的应用及其对植物应对气候变化方面,还应进一步深入追踪。

3.3 土壤中农膜残留特征对大豆生物量累积的影响

随着农膜的长期施用,土壤中农膜累积量不断增加,其对土壤质量和农作物产量的影响日益突出。赵素荣等[45]研究发现:农膜残留不仅影响土壤的物理性质,破坏土壤结构,而且妨碍了作物扎根和对水肥的吸收,影响作物的生长发育和产量。随着土壤中塑料残片的累积,作物幼小的根系生长被残膜阻挠,导致根系发育受损[46]。本研究发现,与对照相比,土壤中塑料残留量导致花期大豆地下生物量降低,其中随PE残留量增加,大豆地下生物量降低,但差异并不显著,这与前人的相关研究结果一致[14-16]。相应地,大豆花期地下生物量随BP残留量的增加显著减少,且在残留量0.1%和0.5%时更为明显,与Meng等[47]的研究结果一致。可能由于降解膜的生物降解所产生的有机物质,影响了土壤生物区系和土壤生物、物理特性,从而间接地影响了作物的生物量、养分循环和植物生产力[48-49]。Qi等[2]的研究表明:生物可降解膜对小麦根系生物量的影响也体现在其生长的初期(前两个月),但在收获期小麦的根系生物量与普通PE膜处理并无差异,这与本研究所观测的结果一致,即在收获期,随BP残留量的增加,BP处理中大豆地下生物量并无明显差异。这表明随大豆生长,其抵抗环境变化的能力加强,而农膜残留对大豆根系生长造成的影响逐渐减弱。此外,尽管与对照相比,各塑料处理下的大豆生物量均有所下降,但对地上生物量的影响并不显著,这可能与大豆的光合作用有关。光合作用光响应曲线和CO2响应曲线参数与生物量间的相关性分析表明:光响应曲线参数中仅LCP与初荚期PE处理的大豆地上生物量有明显的负相关性(P＜0.01),而LCP的增加表明大豆叶片对弱光利用能力减弱,从而导致地上生物量的降低。CO2响应曲线参数Rp与BP处理中大豆地下生物量呈显著负相关,且BP残留量为0.1%、0.5%时Rp增高,说明这两种处理条件下,干物质加快分解,从而导致地下生物量降低。初荚期,PE处理的η和CCP与大豆地上生物量呈显著正相关,BP处理η与大豆地下生物量呈显著正相关,随着η的降低,大豆叶片对于CO2利用能力以及光合速率降低,使得大豆叶片生物量也随之降低。由此表明,CO2响应曲线的相关参数与大豆生物量的相关性更强,且在农膜残留影响下,大豆光合作用对CO2浓度的变化更加敏感,从而在气候变化与极端天气状况下,土壤中农膜残留可能间接影响作物的干物质累积。尽管光响应与CO2响应曲线拟合参数与生物量相关性变化较大,然而各处理间地上生物量的差异并不显著,与Qi等[2]的研究结果略有不同,这可能与塑料类型及颗粒大小、作物种类及种植环境的差异有关。因此,后续研究还应进一步探讨不同塑料类型、残留量大小、作物种类以及植物-土壤之间的互作响应关系,为农膜的可持续应用、可降解塑料的评估以及生态环境安全提供依据。

4 结论

本试验通过设置聚乙烯薄膜(PE处理)及生物降解膜(BP)不同残留量,研究了农膜残留对大豆叶片光合作用光响应与CO2响应曲线及参数的影响。结果表明:随着聚乙烯薄膜(PE)残留量的增加,花期叶片对弱光的利用能力增强,但初荚期叶片对弱光的利用能力减弱;而随着生物降解膜(BP)残留量的增加,花期叶片对弱光的利用差异较小,但初荚期叶片对弱光的利用能力增强,同时也提升了叶片适应强光的能力,从而扩展了叶片光能的利用范围,增强了叶片光合作用的能力。随着大豆生长(初荚期),尽管农膜残留降低了叶片对CO2的利用能力,而叶片对低浓度CO2的利用能力有所提升,同时也降低呼吸作用对于大豆干物质的消耗。基于光响应与CO2响应曲线各拟合参数的变化,虽然各处理中生物量与对照相比均有所下降,但差异并不显著,特别是地上生物量差异较小。然而,随着土壤中塑料残留特征的变化(种类及残留量),大豆在生长初期地下生物量差异明显。随BP残留量的增加,根系生物量显著降低,这可能与BP残留降解过程影响了根系土壤的特征,从而导致地下生物量的变化。此外,通过光合参数与生物量间相关分析可以得出,物质累积所涉及的影响因素众多,不同类型农膜不同残留量的参数间表现出不同的正/负相关关系,而这些仅仅表征了其对大豆光生理特征及生物量的宏观影响,机理及影响机制还有待深入。因此,基于本研究及前人的研究结论,农田土壤中农膜残留对作物生理特征的影响还需进一步的研究,特别是随着可降解农膜的推广,其对作物及土壤生态环境的影响评价应得到更多的关注。