晋中东山地区褐土土壤水分特征的测定与应用

利用压力膜方法测定晋中东山地区褐土的土壤水分，并分析了土壤水分及其能态特征。结果表明：本区土壤水吸力与土壤含水量之间呈明显的幂函数关系，可以用幂函数方程模拟；土壤的当量孔径分布大致随土壤水吸力的增加而减小；土壤的结构性能良好，通气孔隙、毛管孔隙、无效孔隙的当量孔径分布分别占土壤总孔隙度的１８．３％、５１．２％和３０．５％，尤其是毛管孔隙的分布占土壤总孔隙度的一半以上；土壤的比水容量随吸力的增大而迅速减小．当吸力为０．１×１０￣５～０．３×１０￣５ｐａ时，其比水容量的量级为１０￣（－１）［ｍ１／（１０￣５ｐａ·ｇ）］，而当吸力达３×１０￣５～１５×１０￣５ｐａ时，其比水容量的量级下降了１００倍，为１０￣（－３）［ｍｌ／（１０￣５Ｐａ·ｇ）］。表层土壤水分对于作物吸收利用比心土层和底土层更为有效，但深层土壤的持水性能比表层士壤要好。为此，必须提高作物对土壤深层储水的利用。     

古尔班通古特沙漠南缘丘间地土壤水分特征曲线及模拟

为研究干旱荒漠地区土壤水分特征的变化规律并探讨土壤水分特征曲线拟合模型的适宜性,以古尔班通古特沙漠南缘丘间地为研究对象,应用压力膜仪测定了丘间地剖面各土层(0—150cm)的土壤水分特征曲线,结合土壤容重、颗粒分布及饱和含水率等性质,分析了丘间地不同土层水分特征和孔隙大小的分布状况,明确了Gardner模型和van Genuchten模型拟合的适用性。结果表明:(1)丘间地剖面砂质土壤释水量伴随吸力呈规律性变化,低吸力段(0~100kPa)各土层水分特征曲线走势陡直,释水性强;中高吸力段(100~1 500kPa)曲线趋于平缓,持水性强。水分曲线特征与其物理性质一致,各土层物理性质差异较小,曲线的差异也较小。(2)丘间地剖面各土层大孔隙由上到下不断增加且20—40cm土层的小孔隙最多,有助于春季积雪融化迅速入渗到土壤里贮存,对植物吸收利用水分有重大意义。(3)两模型都能很好地拟合丘间地的土壤水分特征曲线,VG模型的拟合效果优于G模型,但G模型形式简单、参数少且易于求解,更便于实际应用。     

基于染色示踪的农膜残留农田土壤优先流特征

随着覆膜农作的发展,农膜残留日益严峻,严重影响土壤水流运动。为了探索农膜残留对水流运动的影响机制,在河套灌区九庄节水综合试验站,设计3种残膜量,分别为0,300,600 kg/hm^2,于成熟期在各试验小区选择典型区域进行20,80 mm水量的染色示踪试验,共6个处理,研究不同残膜量对土壤物理特性的影响,基于染色示踪技术分析不同入渗量和不同残膜量对土壤水流运动的影响,并评价各处理土壤剖面的优先流特征。结果表明:随着残膜量增加,土壤容重、饱和导水率和土壤初始含水率呈下降趋势,而土壤孔隙度呈增大趋势,土壤中随机分布的残膜导致残膜处理(300,600 kg/hm^2)的各物理参数有较高的变异系数。基于土壤染色剖面二值化图像分析,土壤优先流发育程度随残膜量和入渗量增加而提高,300,600 kg/hm^2残膜量处理的最大染色深度(MDSD)和土壤染色剖面的变异系数(Cv)随残膜量增加而增大,与20 mm入渗量相比,80 mm入渗量下残膜处理的MDSD和Cv分别平均增加47.23%和27.34%。不同处理的优先流特征指标差异显著,5个优先流特征指标(染色面积比、基质流深度、优先流比、长度指数和变异系数)均显示残膜处理(300,600 kg/hm^2)的优先流程度高于无残膜处理(0 kg/hm^2),可见残膜量及入渗量的增加均会增大土壤优先流的发生几率和强度。     

点源供水条件下残膜对土壤水分运移的影响

通过室内模拟试验,研究残留地膜对土壤水分运移的影响,阐明不同残膜污染水平下土壤水分入渗过程。试验共设置6个残膜梯度,采用马氏瓶恒压供水175 min,根据时间间隔记录马氏瓶读数和湿润面积,分析不同残膜污染水平下土壤水分入渗过程的差异。研究结果表明,残膜会阻碍土壤湿润锋的运移,残膜量在0~360 kg/hm2区间时,当湿润锋经过0~10 cm土层,湿润锋横向距离随着残膜量的增加呈现先减小后增加的变化,当湿润锋经过10~20 cm土层,横向距离随着残膜量的增加而逐渐降低。残膜对湿润锋的垂向运移有显著的阻碍作用,但残膜量的多少对湿润锋垂向距离变化没有显著影响。残膜的存在会提高土壤湿润比和稳定入渗率,残膜区土壤湿润体变小,水分滞留在湿润体内,影响水分在土壤中正常运移与分布。当残膜量达到720 kg/hm2时,残膜区土壤大孔隙比例增加,导致土壤优势流明显,与其他残膜处理相比,湿润锋的运移加快,湿润比和稳定入渗率降低。该研究从土壤水分运移角度阐明了残膜污染危害的过程和机理,为残膜污染的防治提供了理论的依据。     

残膜对土壤水分入渗和蒸发的影响及不确定性分析

为探究残膜对土壤水分入渗和蒸发过程的影响规律,通过室内土柱试验,设置6个残膜量水平(0、80、160、320、640和1 280 kg/hm2),研究了残膜对湿润锋运移、土壤水分分布、累积入渗量和累积蒸发量及其不确定性的影响。结果表明:随残膜量增加,湿润锋垂直运移速率和累积入渗量逐渐减小;残膜量80 kg/hm2时,湿润锋运移速率大幅下降;累积蒸发量随残膜量增加而递减而蒸发系数呈递增趋势,土壤保水能力减弱;随残膜量增加,0~10和20~45 cm含水率呈降低趋势,而土壤水分的变异系数呈增加趋势,残膜加剧了土壤水分垂直分布的变异性,残膜量320 kg/hm2的处理会出现表土层"板结"现象;基于Gibbs抽样算法分析表明,Kostiakov入渗模型和Rose蒸发模型各参数的95%后验置信区间上下限的差值和标准差均随残膜量增加而增大,累积入渗量和累积蒸发量的95%后验置信区间面积呈增大趋势,土壤累积入渗量和累积蒸发量的不确定性随残膜增多而增强。该研究可探明残膜污染区的土壤水分运动规律,并为提高Kostiakov模型、Rose模型的模拟效率和模拟精度提供参考。     

土壤水吸力对控释尿素养分释放特征的影响

为探究不同土壤水吸力对控释尿素养分释放的影响,根据供试壤质砂土的水分特征曲线设置了6种不同土壤水吸力的处理,采用25℃恒温土壤模拟实验,测定土壤中控释尿素的养分释放特征。结果表明,土壤水吸力为0 k Pa时恒温土壤培养法与静水浸提法测得的控释尿素养分释放特征无显著差异。土壤水吸力为75 k Pa、30 k Pa和0 k Pa的3个处理,测定的土壤孔隙中空气相对湿度均为95%以上,土壤水汽饱和,释放期基本相同,表明在不考虑水分流动及养分扩散状况影响时,土壤水分条件已不再是影响控释尿素养分释放的主要因素。土壤水吸力525 k Pa和260 k Pa(土壤孔隙中空气相对湿度分别为84%和91%)时释放期分别为416.4 d和120.0 d,相对于土壤水饱和时释放期(63.6 d)的相对相差分别为146.8%和59.1%,远超过《控释肥料》行业标准规定的允许差(20%),表明控释尿素的释放率和释放期受到土壤水吸力过高的抑制。土壤水吸力大小直接影响土壤孔隙空气湿度的饱和与否,土壤水吸力对控释尿素养分释放的影响通过土壤水汽作用于肥料颗粒实现。     

用MatLab确定土壤水分特征曲线参数

土壤水分特征曲线是定量研究土壤水运动的重要参数。目前,土壤物理学家普遍采用vanGenuchten模型(简称VG模型)来描述土壤水分特征曲线(h-θ曲线),本文讨论了如何确定模型参数的方法。根据中国科学院栾城试验站大田土壤剖面所采土样实测土壤负压h和土壤含水率θ数据,VG模型中的参数利用MatLab非线性拟合函数来确定,通过用四参数模型计算,计算的θ值与实测的θ值拟合较好,误差相对较小。通过长江河口地区土样检验表明,用MatLab方法确定VG模型参数,计算值与实测值拟合也较好,且计算值与实测值的残差平方和范数<0.001。     

地膜残留量对新疆棉田蒸散及棵间蒸发的影响

为探讨残膜对干旱区农田蒸散耗水特征的影响,在新疆阿克苏典型覆膜滴灌棉田开展2 a小区试验研究,设计0、225、450 kg/hm2共3种不同的地膜残留量,采用水量平衡法,微型棵间蒸发仪法,于主要生育时期测定并计算土壤含水量、蒸散量、棵间蒸发量、作物蒸腾量、棵间蒸发占蒸散的比例。结果表明:随着地膜残留量增加棵间蒸发量、棵间蒸发占蒸散的比例呈增大趋势,而蒸散量及作物蒸腾量则逐渐减小。与无残膜处理相比,225和450 kg/hm2处理全生育期田间无效耗水的棵间蒸发量分别增加了9.27和22.20 mm,棵间蒸发占蒸散的比例增幅分别为2.6%和6.1%,作物蒸腾量降低34.89和55.94 mm。在棉花生育期内,棵间蒸发占蒸散的比例(E/ET)与叶面积指数(leaf area index,LAI)呈幂函数关系,各处理间棵间蒸发占蒸散的比例对叶面积指数的响应差异不同,450 kg/hm2处理蒸发占蒸散的比例随着LAI的增加,曲线下降趋势较明显;无残膜处理棵间蒸发占蒸散的比例与LAI的决定系数最高,均在0.9以上。土壤水分利用率也随残膜量的增加依次降低,当残膜量由0增加到450 kg/hm2时,土壤水分利用率从28.25%降至24.91%,可见,残膜增大了农田的无效耗水,不利于土壤水分的有效利用。研究可为制定缓解或克服残膜危害的应变调控技术提供依据。     

再生水水质对斥水和亲水土壤水分特征曲线的影响

为探求再生水灌溉对斥水和亲水土壤水力特性的影响,该文选用有代表性的斥水黏壤土和亲水黏壤土、斥水砂土和亲水砂土,测定其在自来水、再生水和其他4种生活污水条件下的土壤水分特征曲线,采用主成分分析法得到不同水质综合指标,分析水质综合指标对不同土壤水分特征曲线的影响,采用van Genuchten-Mualem模型对黏壤土土-水曲线的参数进行拟合,并分析水质综合指标对黏壤土累积当量孔径分布、比水容量和水分常数的影响。结果表明:在相同基质吸力情况下,黏壤土含水率随水质综合指标增加(水质变差)而减小,砂土的含水率随水质变化不大;在低吸力段,斥水和亲水黏壤土的比水容量随水质综合指标的增加而增加;土壤进气值与水质综合指标呈显著负线性相关关系(R~2分别为0.94和0.78);相同水质条件下,斥水土壤的进气值比亲水土壤小;随着水质综合指标的增加,斥水和亲水黏壤土的极微孔隙降低,而中等孔隙和大孔隙增加,小于某当量孔径的累积百分比增加;随着水质综合指标的增加,斥水和亲水黏壤土的田间持水率、凋萎系数、有效水和易利用水比例均减小,但再生水对田间持水率和易利用水比例降低作用不显著。研究结果可为大面积再生水灌溉及其管理提供一定的理论依据。     

三峡山地不同坡位土壤的水分特征曲线及水力学性质

山地不同坡位的土壤水力学性质研究可为地表水文过程预测和坡地生态修复提供参考。通过采集三峡库区大老岭林区山地不同坡位的原状土样,测定饱和导水率、当量孔径、水分库容等水力学参数,利用V-G模型拟合土壤水分特征曲线,评价不同坡位土壤的渗透性能和持水能力。结果表明:坡顶土壤的入渗性能最好,平均饱和导水率为108.54 cm/d,其他坡位平均饱和导水率为34.81～47.56 cm/d。坡顶土壤的持水能力最差,吸力值从0增至100 kPa,含水量下降幅度最大(63.67%),其他坡位下降幅度在46.36%～52.07%。坡顶土壤大孔隙体积比最大可达20.95%,其他坡位在12.27%～16.58%。田间持水量库容占总库容的百分比以坡顶土壤最小(45.24%),其余在60.08%～65.33%。用V-G模型拟合土壤水分特征曲线,各样点决定系数R~2均大于0.95,拟合效果较好。饱和导水率与大孔隙体积比、砂粒含量呈正相关,与粉粒含量和黏粒含量呈负相关。V-G模型参数α与进气值相关,α值越小,持水能力越强。试验中α值与砂粒含量、大孔隙体积比呈正相关,与粉粒含量和黏粒含量呈负相关。砂粒和大孔隙越多,土壤持水性越差。参数n的取值影响拟合水分特征曲线的弯曲程度,n值越小,曲线越平缓。试验中n值与砂粒含量呈正相关,与粉粒含量呈负相关。砂粒含量越多,随吸力增大土壤排水越快。可见,山地坡面顶部土壤的持水能力最差,渗透能力最强,其他坡位差异不显著,V-G模型拟合的结果与实测结果一致。     

非饱和土双应力变量广义土水特征曲线模型验证

利用基于轴平移技术的Geo-Expert高级型应力相关土-水特征曲线压力板仪研究不同覆土压力(0、40、100、200 k Pa)对南阳膨胀土土水特征曲线(soil-water characteristic curve,SWCC)的影响;并对提出的考虑土体变形及多孔隙分布形态的双应力变量广义SWCC表征方程进行如下试验验证:1)不同覆土压力下微多孔隙分布形态的南阳膨胀土侧限固结试验及脱湿试验SWCC验证;2)零净法向应力状态双孔隙尺度硅藻土双峰SWCC试验验证;3)不同净围压状态下单孔隙尺度韩国残积土SWCC试验验证;4)多应力路径下法国Bapaume黄土,不同初始干密度下日本Edosaki砂土在脱湿-吸湿过程SWCC试验验证;并比较分析新方程与van Genuchten方程及Fredlund等方程的差异性。结果表明:1)覆土压力会显著改变膨胀土结构及孔隙通道,进而改变SWCC边界效应区、过渡区形态;也改变了双重孔隙尺度土壤的进气值;第1个波峰SWCC进气值均在1 k Pa左右;相比于零覆土压力试样,40、100、200 k Pa覆土压力试样第2个波峰SWCC进气值分别高4.74、17.58、67.23 k Pa;2)未考虑净法向应力影响的单应力状态多峰SWCC、考虑侧限双应力状态多峰SWCC、各向同性净法向应力单峰SWCC、不同脱湿-吸湿路径SWCC及不同初始干密度的SWCC试验拟合曲线均表明,双应力广义SWCC具有包容复杂影响因素的能力;3)新方程能够利用至少3个土-水数据即可拟合出具有较高的精度的整条SWCC。研究为定量描述土壤的持水、渗透及强度特性提供参考。     

残膜对土壤水分水平运动的阻滞作用

为探明不同残膜量对土壤水分水平运动过程的影响规律,设置6个残膜量水平,分别为0,80,160,320,640,1 280kg/hm~2,采用室内水平土柱吸渗法,分析残膜对水平湿润锋运移、累积入渗量、Boltzmann系数和土壤水分扩散率的影响。结果表明:水平湿润锋运移速率、土壤含水率和累积入渗量均随残膜量增加而减小,残膜量≥160kg/hm~2时,水平湿润锋运移速率和累积入渗量大幅降低;残膜对土壤水分水平运动的阻滞系数与残膜量呈正相关关系,残膜阻滞系数的最大增长率出现在80~160kg/hm~2残膜量区间,高达251.15%;随残膜量的增加,阻滞作用较强区域距离入渗点的距离增大,Boltzmann系数和水分扩散率则逐渐减小,土壤基质势对水分水平运动的驱动作用和土壤持水能力逐渐减弱;混有残膜的土壤水分水平运动过程,仍符合土壤水分扩散率单一参数模型和Kostiakov入渗模型。残膜是一种难降解持续性污染物,会降低水平湿润锋运移速率、累积入渗量和土壤水分扩散率,对土壤水分水平运动的阻滞作用随残膜增多而增强。     

紫色土埂坎根-土复合体土水特性与抗剪强度

研究根-土复合体土-水特征曲线与抗剪强度的关系,可为紫色土埂坎根-土复合体强化机理的揭示与埂坎稳定性的维持提供科学依据。选取三峡库区典型紫色土坡耕地埂坎草本植物根—土复合体为研究对象,结合Hyprop2土壤水分特征曲线测量仪、滤纸法与直剪试验,拟合土-水特征曲线,揭示基质吸力对根-土复合体抗剪强度的影响。结果表明:(1)根-土复合体土-水特征曲线明显分为边界效应区、过渡区与非饱和残余区,3种常用模型(B-C、VG、F-X)中F-X模型拟合该曲线效果最好,根-土复合体饱和含水率、进气吸力、残余区含水率以及相同体积含水率下的基质吸力均高于素土。(2)随着体积含水率降低,根-土复合体黏聚力先增大后减小,试验范围内黏聚力最大值51.25 kPa出现在体积含水率约23%时,内摩擦角则线性增大。相同体积含水率下,根-土复合体黏聚力较素土最大增加50%,内摩擦角提升不大。(3)基质吸力对根-土复合体抗剪强度的增强作用具有阶段性特征,各阶段临界吸力值与土-水特征曲线一致,过渡区(基质吸力为3~500 kPa)土体抗剪强度提高明显,进入非饱和残余区后(基质吸力>500 kPa)由于黏聚力下降,土体抗剪强度增速减慢,根-土复合体抗剪强度随基质吸力增大而提升的幅度大于素土。通过建立埂坎根-土复合体土-水特征曲线和抗剪强度的关系,可估测实际工况下的埂坎土体抗剪强度,进而为坡耕地埂坎的建设、维护管理以及坡耕地侵蚀阻控提供理论依据。     

崩岗不同土层土壤水力学特性差异性分析

为研究崩岗不同土层土壤水力学特性的差异性,采用离心法测定不同土层土壤水分特征曲线,筛选出适合的土壤水分特征曲线拟合模型,结合统计模型,推求土壤的当量孔径分布、比水容量、非饱和导水率和扩散率,分析崩岗不同土层土壤水力学参数的变化规律。结果表明,崩岗土层从红土层到砂土层的变化过程中,土壤质地由黏土向砂土变化;Fredlund&Xing模型对崩岗土壤土水特征曲线拟合效果最好;参数θs、α、n随着质地变黏重逐渐减小;随着土层深度的增加,土壤的持水性能降低;土壤比水容量、非饱和导水率和扩散率受土壤质地和基质吸力的共同影响。在低吸力阶段,3个指标随基质吸力变化比较平缓,砂土层土壤比水容量和非饱和导水率最大,扩散率最小;而在高吸力阶段,砂土层土壤的这些指标降低较快,且低于其他土层,各层土壤间导水率和扩散率差异随着基质吸力的增加而增大。     

非饱和土双应力变量广义土水特征曲线理论模型构建

土水特征曲线(soil-water characteristic curve,SWCC)方程是非饱和土力学中最重要的土性表征手段之一。该文评价当前经典的SWCC方程,指出其未具有包容复杂因素的能力,具有灵活性的优点但却同时具有对试验数据量依赖性高的缺点,不能处理多孔隙尺度集群土体固-液-气共同运动及作用的水力-力学耦合效应问题。建立双应力变量广义SWCC概念图示并定义相对体积含水比,基于Fredlund双应力变量理论及van Genuchten土-水表征方程,构建考虑土体变形及多孔隙分布形态的双应力状态变量的广义SWCC方程。相较于2个参数的Brooks等的方程、3个参数的van Genuchten方程以及4个参数的Fredlund等的方程,广义SWCC方程仅3个参数,其中2个参数在双对数坐标系的"相对体积含水比-吸力"平面中进行最小二乘法线性拟合得到,仅1个参数需非线性最小二乘法拟合得到。该模型可利用不同应力状态下的至少3个土水试验数据点,绘制出1条具有适宜精度的单峰SWCC;方程考虑了多峰孔隙概率密度函数分布及土体变形因素,实现了从应力历史推广到应力状态的广义情况,为定量描述不同孔隙结构土体双应力状态下的持水特性、渗透特性和强度特性提供了一条途径。     

滤纸法测定干湿循环下膨胀土基质吸力变化规律

为了获得干湿循环作用下膨胀土基质吸力的变化规律,首先采用人工模拟降雨和蒸发的方法开展了膨胀土室内干湿循环试验,然后利用滤纸法进行了不同含水率下试样的基质吸力测定试验,获得了干湿循环条件下膨胀土的土水特征曲线,求出了相应的进气值与残余值,结合Fredlund土水特征曲线模型对经历不同干湿循环次数下的土壤土水特征曲线进行了拟合,最终建立了考虑干湿循环效应的膨胀土土水特征曲线模型。结果表明:1)随着干湿循环次数的增加,土壤的进气值呈下降趋势,从循环1次时的134.5 k Pa降至循环4次时的58.5 k Pa,降幅达56.5%。从循环1次至2次的进气值下降较大,往后降幅明显减小,趋于基本稳定,这表明对土壤进气值的影响以初次干湿循环为主。2)残余值亦呈下降趋势,从循环1次时的1 040.5降至循环4次时的528.5 k Pa,降幅达49.2%。每经历一次干湿循环,残余值降幅均较大,尚未趋于稳定,这表明干湿循环效应对土壤残余值的影响比对土壤进气值的影响要大。3)新建土水特征曲线模型中的拟合参数与干湿循环次数成较好线性关系,表明随着干湿循环次数的增加,土壤进气值逐渐减小,水分变化速率有所降低,而残余含水率逐渐增加。该成果可为深入研究土壤基质吸力及其工程应用提供参考。     

非常规水水质对红壤水分运移参数的影响及其过程模拟

为了探究非常规水水质对亚热带红壤水分运移参数的影响,采用蒸馏水(CK)、再生水(RW)、稀释2,4,6倍再生水(RW-D2、RW-D4、RW-D6)和生活污水(_(WW))等6种不同水质的非常规水模拟红壤水分入渗过程,利用水平土柱吸渗法测定红壤湿润锋运移速率和水分扩散度,应用压力膜仪法测定了土壤水分特征曲线并且借助van Genuchten(VG)模型提取了相应的土壤残余含水量θ_r、土壤饱和含水量θs、α和n等特征参数值。结果表明:不同水质入渗下红壤湿润锋运移速率、水分扩散度均存在着一定差异,湿润锋运移速率由快到慢分别为V_(CK)≈V_(RW-D2)V_(RW-D6)≈V_(RW-D4)V_(WW)V_(RW);在同一含水率下D_(RW-D6)(θ)D_(RW-D4)(θ)D_(WW)(θ)D_(RW-D2)(θ)D_(CK)(θ)D_(RW)(θ)。不同水质处理对VG模型参数变化率的影响程度不同,对α的影响最大,对θ_r的影响最小。因此,非常规水灌溉会降低红壤入渗速率、持水性能以及抗旱能力。     

新疆棉田地膜残留对棉花产量及土壤理化性质的影响

为了探讨地膜残留对棉花产量的影响,该文在全面调查新疆北疆短绒棉主栽区地膜残留状况的基础上,根据地膜残留趋势,人工设置地膜残留密度梯度（0、250、500、1 000、1 500、2 000 kg/hm2）,选择新陆早33号（发达根系）和新陆早13号（不发达根系）2个棉花主栽品种,研究了地膜残留对棉花产量性状和土壤理化性质的影响。结果表明,随着地膜残留密度的升高,2个品种产量和成苗率均呈下降趋势,新陆早13号减产趋势更加明显,早熟度明显提前。如果按照现有的地膜残留趋势,则覆膜68 a左右,即现在往后38 a,残膜密度将达到1 000 kg/hm2,新陆早33号和新陆早13号产量分别下降18.1%、13.5%,成苗率分别下降8.4%、16.5%,早熟度分别上升17.0%、26.7%。当残膜密度为2 000 kg/hm2（覆膜141 a左右）时,新陆早33号和新陆早13号产量分别下降38.3%、45.2%,成苗率分别下降13.0%、21.1%,早熟度分别上升32.2%、27.7%;生物量、根冠比和根系表面积等生长指标都与没有残膜下生长的棉株存在显著差异。地膜残留可导致土壤理化性质恶化,水分分布不均,土壤营养下降。其中在残膜密度为2 000 kg/hm2时,碱解氮、速效磷分别下降55.0%、60.3%。研究表明新疆地膜残留对棉花产量和土地产生明显影响,如不及时采取措施,将严重危害棉花可持续发展与土地的可持续利用,必须引起足够重视。同时发现,应用、培育发达根系的棉花品种可以一定程度上缓解地膜残留对产量的不利影响。