生态学报  2014, Vol. 34 Issue (21): 6287-6296

文章信息

王恩姮, 赵雨森, 夏祥友, 陈祥伟
WANG Enheng, ZHAO Yusen, XIA Xiangyou, CHEN Xiangwei
冻融交替后不同尺度黑土结构变化特征
Effects of freeze-thaw cycles on black soil structure at different size scales
生态学报, 2014, 34(21): 6287-6296
Acta Ecologica Sinica, 2014, 34(21): 6287-6296
http://dx.doi.org/10.5846/stxb201307031829

文章历史

收稿日期:2013-7-3
修订日期:2014-5-30
冻融交替后不同尺度黑土结构变化特征
王恩姮, 赵雨森, 夏祥友, 陈祥伟     
东北林业大学 林学院, 哈尔滨 150040
摘要:冻融交替是改变黑土结构、加剧土壤侵蚀的重要因子.以典型黑土区耕作土壤为研究对象,采用野外季节性冻融循环与室内模拟冻融循环相结合、X射线计算机断层摄影(CT)与扫描电子显微镜(SEM)相结合的方法,通过水分物理性质、团聚体破坏率、孔隙数目、孔隙面积、孔隙成圆率、孔隙Feret直径的测定与分析,研究了冻融交替后0-40 cm、40-80 cm和120-160 cm 3个土层以及田间季节性冻融环刀、室内模拟冻融CT扫描和室内模拟冻融SEM3种方式下黑土结构特征的变化规律.结果表明:冻融交替能够对不同土层和不同尺度的耕地黑土结构产生不同程度的影响.季节性冻融后,表层土壤容重升高,非毛管孔隙度和持水能力显著降低(P <0.05),40-80 cm土层团聚体破坏率增加40.97%(P <0.05),土壤抗蚀性有所削弱,120-160 cm土壤没有受到季节性冻融的显著影响.CT扫描尺度上,3个土层均以1-2 mm径级的孔隙数目为最多,形状也相对规则、接近圆形;冻融循环没有对表层土壤大孔隙结构产生影响,却能够显著降低40-80 cm土层范围内大孔隙面积以及Feret直径(P <0.05).SEM扫描显示冻融后土壤表面粗糙度增加,颗粒松散、脱离,孔壁断裂,证明了冻融交替对土壤微结构的破坏作用;同时结合电子能谱的元素分析可知冻融交替能够改变土壤颗粒表面化学特征.
关键词冻融交替    土壤微结构    CT    SEM    
Effects of freeze-thaw cycles on black soil structure at different size scales
WANG Enheng, ZHAO Yusen, XIA Xiangyou, CHEN Xiangwei     
Northeast Forestry University, College of Forestry, Harbin 150040, China
Abstract:Alternating freezing and thawing is a critical factor associated with soil structure change and accelerates soil erosion in the black soil region of Northeast China. Based on the soils sampled from fields of the black soil region, the effects of freeze-thaw cycles on soil structure at different soil depths (0-40 cm, 40-80 cm, 120-160 cm) and size scales (field core sampling scale of seasonal freeze-thaw cycles, computerized tomography [CT] scale of artificial freeze-thaw cycles, and scanning electron microscope [SEM] scale of artificial freeze-thaw cycles) were studied. We measured and analyzed the hydro-physical properties, percentage of aggregate disruption (PAD), pore number, pore area, pore roundness, and Feret diameter of pores in the field and lab by using CT and SEM methods. Results showed that: freeze and thaw alternation affected black soil structure from different soil depths and different size scales both in seasonal field condition and artificial freeze and thaw cycles condition. At the scale of seasonal freeze-thaw cycles, bulk density of top soil (0-40 cm) increased, total porosity decreased significantly, with non capillary porosity decreasing dramatically. Water holding capacity decreased consequently(P <0.05). Within the soil depth of 40-80 cm, bulk density and porosity did not change significantly, but PAD increased by 40.97%(P <0.05)resulting in weak erosion resistance. Hydro-physical properties and aggregate breakdown characterization of soils in the depth range of 120-160 cm were not changed significantly by seasonal freeze and thaw cycles. At the size scale of CT, the alternation did not change the characterization of macro pores in 0-40 cm but significantly reduced the average pore area, roundness and diameter after freeze and thaw cycles. While the average area and diameter decreased significantly in the 40-80 cm depth, this might result in poor infiltration and water movement to the subsoil. Pore area of 1-2 mm and > 5 mm decreased 9.58% and 42.19% (P <0.05), respectively. Structure of soils from 120-160 cm underwent by artificial freeze-thaw cycles showed similar results as field scale did, were not affected significantly by alternate artificial freeze and thaw cycles. The number of pores with diameter between 1 and 2 mm was predominant the pore size distribution through three soil depths. At the scale of SEM, the roughness of the soil surface and alveolate cells were significantly increased; we also observed micro-aggregate displacement and rearrangement, and disruption of the pore walls. Integrating electron spectroscopy results of 40-80 cm depths, carbon and calcium did not change after 6 freezing and thawing cycles, but oxygen group elements decreased. Magnesium, silicon, aluminum, potassium and iron elements increased to different degrees, and the total increment equaled to the decrease quantity of oxygen group element. The packing arrangement of oxygen group elements and metallic oxide are the key factors to determine the type and surface chemistry of soil clay minerals, so the amount of change in oxygen group and other elements gives indirect evidence that freeze and thaw cycles can affect soil microstructure. But to determine whether or not freeze and thaw cycles would change the crystal lattice structure of clay minerals, further research will be needed and must be approached from the view of soil chemistry and minerology.
Key words: Freeze-thaw cycles    Soil microstructure    CT    SEM    

冻融交替现象普遍存在于北半球中、高纬度地区及高海拔地区,是影响该区土壤生态的重要因素之一[1, 2]。能够较为强烈地影响土壤物理结构[3, 4],进而改变土壤生物区系特征[5, 6],越来越多的研究也表明冻融交替区是CO2和N2O等温室气体的重要释放源[7, 8]。同时,冻融交替对改善土壤环境起着积极的作用,能够调节土壤的紧实度[9],释放有效养分,提高土壤的可耕性等[10]。但冻融交替也是改变土壤结构进而诱发土壤侵蚀的主要过程之一,尤其是在土壤下层存在冻层和滞水层的情况下[11, 12, 13, 14]。冻融循环对土壤结构与功能的影响极其复杂,迄今尚无共识的结论[15]

我国黑土区主要分布在东北松嫩平原东部及北部的山前台地,每年经历6—8个月的季节性冻融,由于传统耕作方式和理念的限制,大部分耕地在漫长的冬季几乎为裸地,与免耕系统以及秸秆还田的耕地相比,冻融强度和频度有所增强,改变土壤结构,影响土壤侵蚀的发生发展规律[16]。20世纪90年代的研究表明,黑土区侵蚀沟沟头的溯源侵蚀速率为每年1 m左右,近年在典型黑土区个别侵蚀沟的前进速度可达12 m/a[17, 18]。张永光等指出黑土区初春融雪后地表裸露,在表层土壤解冻以及存在季节性冻土层的情况下,伴随降雨的共同作用是加剧黑土区侵蚀的主要原因[19]。胡刚等人也提出了黑土区“冬季坍塌沉积,春季降雨冲刷”的侵蚀模式[18]。除了表层土壤结构受冻融交替的影响之外,土体内部水分运移以及热量的动态变化也将影响下层土壤结构的变化,进而造成潜在的侵蚀风险[20]。然而关于东北黑土区冻融交替对土壤结构尤其是对微结构影响方面的研究至今却鲜有报道。

本研究采用野外季节性冻融交替以及室内模拟冻融循环相结合的方法,借助X射线计算机断层摄影(computerized tomography,CT)和扫描电子显微镜(scanning electron microscope,SEM)手段[21, 22],以期实现冻融循环交替对不同深度和不同尺度黑土结构的研究,揭示季节性冻融对黑土结构的影响程度和规律,为黑土结构质量的演变以及土壤侵蚀机理方面的研究提供证据。

1 研究区概况与研究方法 1.1 研究区概况

试验地点位于黑龙江省西北部克山县境内的克山农场 (48°12′ — 48°23′N,125°8′ — 125°37′E),属克拜漫川漫岗地带,海拔 240 — 340 m,平均坡度 3°。区域内土壤类型以粘化湿润均腐土占绝对优势,只在部分低洼处镶嵌分布少量的草甸土,属典型黑土区[23]。该区温带大陆性季风气候特征明显,年均气温0.9 ℃,≥10 ℃ 有效积温 2296.2 ℃,年均降水量501.7 mm,年均蒸发量 1329. 4 mm,无霜期115 d,土壤冻结期从10月中旬至翌年6月末,最大冻结深度可达 2.5 m。

1.2 样品采集与制备

为进一步证明和解释典型黑土区土壤结构对季节性冻融循环的相应规律,以及冻融交替对诱发沟蚀、加剧侵蚀沟发展的潜在影响,于2009年10月初,以黑龙江省克山农场13连队区内侵蚀沟为研究对象,沿着土壤剖面层次清晰的侵蚀沟壁,修整取样剖面。在0—40 cm、40—80 cm和120—160 cm土层范围内用环刀(100 cm3)采集原状土样,在取样层内再进行机械法分层,每10 cm一层,各取3个环刀,每40 cm的取样层共计12个重复;并在采样点附近取非原状土,分别以供季节性冻融前水分物理性质以及团聚体破坏率的测定。同时采用自制内径和高分别为10 cm的有机玻璃管进行原装土柱的取样,每层重复取样3次,以供模拟冻融试验所用。将取于0—40 cm、40—80 cm和120—160 cm范围内的有机玻璃管原状土柱浸泡(水面与有机玻璃管上沿平齐)在水中12 h进行饱和处理,然后将土柱取出平置于3 cm高的沙子层上自然渗水2 h,水分处理后将原状土柱样品置于自动循环低温冻融箱内(型号DRX-Ⅰ)进行冻融循环处理,循环次数为6次,每次冻融循环历时15 h,其中冻结时间为12 h,解冻时间为3 h,冻结最低温度为-20 ℃,解冻最高温度为5 ℃。从上至下,3个土层样品在进行冻融循环时的初始平均体积含水量依次为31.71%、34.36%和23.77%。于2010年6月末季节性冻融完全结束后,在侵蚀沟壁第一次采样点附近进行环刀原状土和非原状土取样,以供季节性冻融后水分物理性质以及团聚体破坏率的测定,团聚体破坏率(PAD0.25)= 100%×大于0.25 mm水稳团聚体(g)/大于0.25 mm风干团聚体(g)[24]。所采土样质地为粘壤土,pH为6.2,3个土层有机质和粘粒含量至上而下依次为6.5%、4.6%、2.4%和33.1%、39.1%、36.5%。

1.3 CT扫描与处理 1.3.1 CT扫描

试验采用的是黑龙江省中医药大学CT扫描中心,美国GE公司生产的Lightspeed16排螺旋CT扫描仪,空间分辨率为0.275 mm。扫描参数为:电压140 kV,电流60 mA,扫描间隔1.5 s,扫描层厚2 mm,窗宽(显示图像时所选用的CT值范围)和窗位(窗宽上、下限CT值的平均数)均为1300。对每个土柱样品进行横断面扫描,扫描位置为中心线1次,中心线两侧每隔10 mm扫描1次,各3次,共7次,以下数据的分析均采用7次扫描的平均结果。

1.3.2 图像处理

将原始的CT图像以*.JPG文件形式保存(图 1A),在Photoshop软件中将每一个横断面CT扫描图像切割成214×214像素的中心方块(图 1B),对应的实际边长为7.55 cm,仍以*.JPG文件形式保存。应用ImageTool3.0软件再次将图像转化为只有孔隙(黑色)和固体颗粒(白色)的灰度图像(图 1C)。在转化过程中,首先对含有已知大小的大孔隙的原状土柱进行CT扫描,然后手动设定阈值(Processing → Threshold →Manual),反复调试,直至图像中的孔隙直径与已知孔隙直径最相近为宜。在调试过程中发现,自动设定阈值(Processing → Threshold → Automatic)后的孔隙大小更接近已知孔隙,所以阈值采用自动设置法。阈值设定后,通过对已知孔隙直径的测量,设置度量单位(Setting → Calibrate Spatial Measurement),作为图像定量处理的依据。

图 1 灰度图像二值化处理 Fig. 1 Segmentation of gray-scale CT images
1.3.3 数据处理

通过ImageTool3.0软件的查找功能(Analysis→ Find Objects→Automatic)和分析功能(Analyze Analysis →Object Analysis → Analyze),获得每一个可以识别的孔隙的定量参数。以Feret直径(经过不规则物体的中心,任意方向每隔10°的36个直径平均值)为孔隙径级分类标准,选择孔隙数目、面积和成圆率(4πA/P2,式中A是孔隙面积;P是孔隙周长,成圆率为1时,孔隙为圆形)3个指标进一步对原状土柱内不同径级(<0.5 mm、0.5—1 mm、1—2 mm、2—5 mm和≥5 mm)孔隙的大小和形状进行分析。采用SPSS11.5进行统计分析。

1.4 SEM扫描与处理

模拟冻融循环前后,将土柱从有机玻璃管中小心推出,沿着自然纹理用手将土柱轻轻掰成约为1.5 cm×1.5 cm×1.5 cm的土块,每个土块至少保证一个面的纹理自然,没有受到人为掰修过程的挤压。将修整好的土块放入瑞士BAL-TEC公司生产的SCD 005型离子溅射仪进行喷金处理,采用美国FEI公司生产的Quanta200型,附带 X 射线能谱仪的扫描电子显微镜观察不同土层冻融前后微形态特征。

2 结果与分析 2.1 季节性冻融后黑土剖面结构变化特征(田间环刀尺度)

季节性冻融对不同土层结构特征影响程度不同,0—40 cm范围内的土层经过一季的冻融循环总孔隙度显著降低,通过非毛管孔隙/毛管孔隙的比值以及毛管持水量的变化进一步可知,季节性冻融循环能够显著减少表层土壤的非毛管孔隙和毛管孔隙的比例,且非孔隙减小的幅度更大,虽然团聚体破坏率有所增加,但是差异未达显著水平。与表层土壤相比,40—80 cm范围内的土壤对季节性冻融更为敏感,容重在季节性冻融后显著降低,团聚体破坏率相对冻融前提高了40.97%(P<0.05),由此可见季节性冻融削弱了该层土壤的抗蚀性。虽然总孔隙度、非毛管/毛管孔隙没有明显变化,但是饱和持水量和毛管持水量显著提高,这与该层粘粒含量高有关,增加了比表面积,持水能力增强。当土层深超过120 cm时,各水分物理性质指标以及大于0.25mm团聚体破坏率没有受到季节性冻融的显著影响。

2.2 模拟冻融循环条件下黑土孔隙结构变化特征(CT扫描尺度)

冻融循环后,表层0—40 cm范围土层的大孔隙结构没有发生明显变化,虽然孔隙数目、面积和Feret直径略有降低但并均未达到显著水平,成圆率略有增加,孔隙的形状更加不规则(表 2)。与表层土壤不同,40—80 cm范围土层于冻融前后大孔隙数量几乎没有增减,孔隙形状也没有显著变化,但孔隙大小的变化差异明显,面积从冻融前的12.22 mm2减小至6.51 mm2,Feret直径则相应从2.31 mm降低至1.91 mm(P<0.05)。对于120—160 cm范围土层而言,冻融循环降低了大孔隙数目(P<0.05),但对孔隙的大小和形状均无显著影响,与田间季节性冻融的结果一致(表 1)。

表1 模拟冻融循环前后水分物理性质及团聚体破坏率 Table 1 Hydro-physical property and percentage of aggregate disruption before and after freeze-thaw cycles
土层深 Soil depth/cm处理 Treatment容重 Bulk density/ (g/cm3)总孔隙度 Total porosity /%非毛管/毛管 NCP/CP饱和持水量 Saturated water holding capacity/%毛管持水量 Capillary water holding capacity/%团聚体破坏率 PAD0.25/%
不同大写字母表示同土层冻融前后差异显著,P<0.05. NCP / CP = Non capillary porosity / capillary porosity,PAD0.25=Percentage of aggregate (> 0.25mm) disruption.
0—40冻融前1.06±0.06A53.88±3.01A0.12±0.01A53.83±2.53A47.57±2.61A48.28±5.28A
冻融后1.11±0.06A42.99±4.30B0.03±0.02B39.10±4.85B37.83±4.37B53.29±9.96A
40—80冻融前1.25±0.16A48.66±8.19A0.06±0.06A39.05±13.97A36.91±10.79A35.73±9.42A
冻融后1.14±0.12B48.52±6.15A0.05±0.03A42.79±8.64B40.51±7.35B50.37±10.62B
120—160冻融前1.20±0.42A35.17±16.70A0.02±0.02A25.99±14.46A25.57±13.80A43.73±10.25A
冻融后1.31±0.05A39.89±5.59A0.02±0.06A30.63±6.46A30.17±4.82A59.10±15.06A
表2 冻融前后孔隙平均参数的变化 Table 2 Average change of pore characterization before and after freeze-thaw cycles
土层深度 Soil depth /cm处理 Treatment数目 ± sd Number面积(mm2) ± sd Area成圆率 ± sd RoundnessFeret直径 ± sd Feret diameter
不同大写字母表示同土层冻融前后差异显著,P<0.05.
0—40冻融前120±22 A3.71±1.08 A1.19±0.14 A1.68±0.17 A
冻融后110±20 A2.93±0.56 A1.22±0.11 A1.56±0.13 A
40—80冻融前52±3 A12.22±3.69 A1.16±0.11 A2.31±0.22 A
冻融后50±9 A6.51±2.12 B1.23±0.10 A1.91±0.20 B
120—160冻融前78±7 A2.60±0.68 A1.28±0.19 A1.43±0.10 A
冻融后68±3 B2.51±0.69 A1.19±0.11 A1.48±0.15 A

表 3所示,冻融循环对40—80 cm土层范围内各个径级的孔隙数目没有显著影响,但0—40 cm土层内≥5 mm的孔隙和120—160 cm土层内<1 mm的孔隙数目则有所减少(P<0.05)。1—5 mm的孔隙数目在全剖面范围内冻融循环前后均差异显著。同时发现,各个土层在冻融循环前后均是中间径级的孔隙数目较多,以1—2 mm的孔隙最多,>5 mm的孔隙最少。

表3 冻融前后各径级孔隙数目的变化 Table 3 Pore number change of different diameter sizes before and after freeze-thaw cycles
土层深 Soil depth/cm处理 Treatment<0.5 mm0.5—1 mm1—2 mm2—5 mm≥5 mm
不同大写字母表示同土层冻融前后差异显著,P<0.05
0—40冻融前15±6 A26±5A46±10A29±8 A3±1 A
冻融后13±3 A26±9 A43±10A25±4 A2±1 B
40—80冻融前6±2 A14±2 A16±4 A12±2 A4±1 A
冻融后7±3 A11±3A19±2 A11±3 A3±1 A
120—160冻融前12±4 A18±4 A34±3 A12±2 A2±1 A
冻融后7±3 B14±3 B35±5 A11±4 A1±1 A

冻融循环对各个径级孔隙平均面积的影响如表 4所示,仅有40—80 cm土层内1—2 mm和>5 mm的孔隙面积变化显著,分别较冻融循环以前降低了9.58%和42.19%(P<0.05),对其他土层、其他径级的孔隙面积没有明显的影响,由此也进一步说明40—80 cm土层对季节性冻融交替的敏感程度与其上下土层相比,更为明显。

表4 冻融前后各径级孔隙平均面积的变化 Table 4 Average pore area change of different diameter sizes before and after freeze-thaw cycles/mm2
土层深 Soil depth/cm处理 Treatment<0.5 mm0.5—1 mm1—2 mm2—5 mm≥5 mm
不同大写字母表示同土层冻融前后差异显著,P<0.05
0—40冻融前0.15±0.02 A0.45±0.06 A1.69±0.09 A7.32±0.95 A37.80±15.67 A
冻融后0.15±0.01 A0.48±0.02 A1.54±0.2 A7.15±0.69 A23.84±17.00 A
40—80冻融前0.15±0.02 A0.49±0.03 A1.67±0.07 A7.82±1.02 A127.19±5200 A
冻融后0.16±0.01 A0.49±0.05 A1.51±0.14 B7.81±0.67 A73.53±26.47 B
120—160冻融前0.14±0.02 A0.48±0.05 A1.61±0.08 A6.22±0.84 A28.00±21.00 A
冻融后0.14±0.01 A0.50±0.03 A1.68±0.09 A4.79±2.2 A22.64±22.15A

成圆率是表示不规则图形接近标准圆形的程度,成圆率越接近1说明不规则图形越近似圆形。由表 5可知,冻融循环处理能够小幅度改变各个径级孔隙的成圆率,但是都没有达到显著水平,仅有表层1—2 mm的孔隙成圆率显著增加。同时发现,中间径级的孔隙成圆率较为接近1,以1—2 mm的孔隙形状最接近圆形,其次是0.5—1 mm径级的孔隙,较大径级和较小径级的孔隙成圆率与1的差值逐渐增加,说明其孔隙形状越来越不规则。

表5 冻融前后径级孔隙成圆率的变化 Table 5 Roundness change of different diameter sizes before and after freeze-thaw cycles
土层深 Soil depth/cm处理 Treatment<0.5 mm0.5—1 mm1—2 mm2—5 mm≥5 mm
不同大写字母表示冻融前后差异显著,P<0.05
0—40冻融前3.32±0.91 A1.38±0.14 A0.85±0.02 A0.62±0.05 A0.36±0.06 A
冻融后3.37±0.51 A1.31±0.11 A0.92±0.05 B0.65±0.05 A0.27±0.19 A
40—80冻融前3.45±0.84 A1.27±0.15 A0.83±0.07 A0.53±0.05 A0.25±0.08 A
冻融后3.84±0.80 A1.27±0.22 A0.89±0.06 A0.57±0.06 A0.34±0.07 A
120—160冻融前2.81±0.82 A1.31±0.12 A0.95±0.01 A0.72±0.08 A0.22±0.13 A
冻融后3.16±0.37 A1.20±0.11 A0.95±0.06 A0.77±0.12 A0.20±0.19 A
2.3 模拟冻融循环条件下黑土微形态变化特征(SEM扫描尺度)

由冻融前的SEM图像可以直观看出黑土剖面不同层次土壤的微结构,并且能够粗略判断黑土剖面粘土矿物的类型。0—40 cm表层土壤中清晰可见腐殖质与絮片状粘粒物质连成松散的团聚体和残存着的半分解有机质,腐殖质化较强(图 2a);次生粘土矿物主要以云母为主,并伴有原生石英的残存体(图 2b,c)。冻融前纹理清晰,结构性、完整性好,模拟冻融循环6次以后,则出现明显的裂隙(图 2d,e),出现蜂窝状突起(图 2e,f),表面粗糙度增加,小骨骼颗粒较冻融前松散,并表现出脱离的迹象。

图 2 冻融循环6次后0—40 cm 土层SEM 图像 (a—c 为冻融前;d—f为冻融后) Fig. 2 SEM images from 0—40 cm before (a—c) and after (d—f) freeze and thaw cycles

图 3所示,40—80 cm土层细骨骼颗粒排列紧密有序,细粒物质定向性不明显(图 3b,c),并伴有非常明显的片层状特征(图 3a),由此也印证了该层粘粒含量相对较高,比表面积大的特性。冻融循环后,土体表面出现纵横交错的裂隙(图 3d,f),同时也有疑似孔道的蜂窝状突起出现(图 3f)并伴随骨骼颗粒的定向性趋势,有较细颗粒松散脱落或重新排列(图 3)。

图 3 冻融循环6次后40—80 cm土层 SEM 图像 (a—c 为冻融前;d—f为冻融后) Fig. 3 SEM images from 40—80 cm before (a—c) and after (d—f) freeze and thaw cycles

图 4所示,120—160 cm 土层土壤骨骼颗粒排列更加紧密(图 4a),絮凝状特征明显并伴有团聚形体的块体出现(图 4c),由于更靠近发育母质,次生矿物和原生矿物结构清晰可见(图 4b)。冻融循环6次以后,絮凝状表面被破坏(图 4d),孔道增加,孔壁粗糙(图 4e),有明显断裂和重新排列的特征(图 4f)。

图 4 冻融循环6次后120—160 cm 土层SEM 图像 (a—c 为冻融前;d—f为冻融后) Fig. 4 SEM images from 120—160 cm before (a—c) and after (d—f) freeze and thaw cycles

通过进一步分析对季节性冻融交替最为敏感的40—80 cm土层范围内各元素组成比例的变化可知(图 5),冻融循环并没有对土壤颗粒表面C元素的含量产生影响,由于典型黑土剖面通体Ca含量较少,也没有在冻融前后发生变化,但其他阳离子含量在冻融循环6次以后均有不同程度的增加,其中以Fe元素增加的幅度为最大,各种阳离子总增加量等于氧元素含量的降低量。氧元素与Si、Al、Fe、Mg等元素在土壤层状硅酸盐粘土矿物和粘粒氧化物中的组合、结构决定着粘粒矿物的类型、表面活性等一系列主要性质。在各种阳离子与O2-阴离子键结合成的硅酸盐和氧化物矿物的结构里,阳离子都是填充在O2-阴离子的间隙之中[25],因此,氧元素以及其他元素含量的变化可以作为冻融循环改变土壤微结构和土壤颗粒表面化学特征的间接证据。

图 5 模拟冻融前后40—80 cm土壤元素组成的变化 Fig. 5 Element change from 40—80 cm before and after freeze and thaw cycles
3 讨论与结论

冻融交替能够对不同深度和不同尺度的耕地黑土结构产生不同程度的影响。季节性冻融后,田间表层(0—40 cm)土壤总孔隙度显著减少,且非毛管孔隙减小的幅度较大,持水能力迅速减弱;这与温美丽以及笔者前期的研究结果一致,即多次冻融循环能够使疏松的黑土土壤结构变紧实,同时使紧实的黑土土壤结构变得疏松,使土壤容重趋于一致[26, 27]。40—80 cm土层的土壤孔隙没有显著变化,但是团聚体破坏率大幅降低,土壤抗蚀性有所削弱;这可能与该层粘粒含量有关系,Zhao等在相同区域的研究结果表明,相对表层土壤和1 m以下的土层而言,50—70 cm粘化层土壤对季节性冻融响应最为强烈,无论是坡上还是坡下,团聚体破坏率均在冻融前后发生了显著的变化[28]。而120—160 cm土层的孔隙结构以及团聚体破坏率没有受到季节性冻融的显著影响,这与农田黑土水分动态的垂直变化有关。由于冻结过程中土壤水分在垂直方向上的再分配,1 m以内的土壤含水量和持水能力会发生较大的变化,而120—200 cm或者更深的土层,常年保持湿润状态,土壤体积含水量的增减不超过6%—7%,主要受土壤母质的物理作用力控制,受生物-气候因素的影响较小[29, 30]

模拟冻融循环后,CT扫描尺度上,40—80 cm土层的平均孔隙面积显著降低,这与田间尺度的研究结果相吻合,也是导致季节性冻融后黑土垂直剖面心土层滞水的原因[30, 31]。120—160 cm土层的平均孔隙大小与形状同样没有受到冻融循环的显著影响,这可能与样品本身的孔隙度小、含水量低有关系;在同一冻融温差下,低含水量土壤经过冻融循环后孔隙结构的变化较小[32]。3个土层范围内虽然随着孔隙径级的增大,相应的孔隙面积也逐渐增加,但是却以中间径级1—2 mm孔隙数目最多,形状也最规则、接近圆形,随着孔隙径级的增大或者减小,孔隙数目和面积逐渐减少,这个径级的孔隙可能是冻融过程中水分垂直运移的主要通道。SEM扫描尺度上,冻融前3个土层的土壤微结构纹理清晰,结构性、完整性好,冻融后表面粗糙度增加,出现蜂窝状突起,并伴有小骨骼颗粒松散、脱离的现象,也有明显的孔壁断裂和重新排列的特征;虽然不能量化SEM扫描的结果,但是结合40—80 cm土层的能谱分析中氧元素以及其他阳离子含量的变化,可以间接证明冻融循环对土壤微结构和土壤颗粒表面化学特征的影响。黑土区粘土矿物以2 ∶ 1型的蒙脱石为主,胀缩性大,冻融过程中随着水分的进出,晶层胀缩、晶格开放,粘土矿物吸附的阳离子暴露,致使在能谱分析时更易捕获,表现出较高的含量。但冻融循环是否能够破坏粘土矿物晶格的基本构型,还需从土壤化学和矿物学的角度进行深入研究。

参考文献
[1] Robinson D A, Dewey K F, Heim R R. Global snow-cover monitoring: An update. Bulletin of the American Meteorological Society, 1993, 74: 1689-1696.
[2] Wang L F, Cai Y J, Xie H T. Relationships of soil physical and microbial properties with nitrous oxide emission under effects of freezing-thawing cycles. Chinese Journal of Applied Ecology, 2007, 18(10): 2361-2366.
[3] Oztas T, Fayetorbay F. Effect of freezing and thawing processes on soil aggregate stability. Catena, 2003, 52(1): 1-8.
[4] Sahin U, Angin I, Kiziloglu F M. Effect of freezing and thawing processes on some physical properties of saline-sodic soils mixed with sewage sludge or fly ash. Soil & Tillage Research, 2008, 99(2): 254-260.
[5] Walker V K, Palmer G R, Voordouw G. Freeze-thaw tolerance and clues to the winter survival of a soil community. Applied & Environmental Microbiology, 2006, 72(3): 1784-1792.
[6] Yang S Z, Jin H J. Physiological and ecological effects of freezing and thawing processes on microorganism in seasonally-froze ground and in permafrost. Acta Ecologica Snica, 2008, 28(10): 5065-5074.
[7] Ludwig B, Wolf I, Teepe R. Contribution of nitrification and denitrification to the emission of N2O in a freeze-thaw event in an agricultural soil. Journal of Plant Nutrition & Soil Science, 2004, 167(6): 678-684.
[8] Teepe R, Ludwig B. Variability of CO2 and N2O emissions during freeze-thaw cycles: results of model experiments on undisturbed forest-soil cores. Journal of Plant Nutrition & Soil Science, 2004, 167(2): 153-159.
[9] Unger P W. Overwinter changes in physical-properties of no-tillage soil. Soil Science Society of America Journal, 1991, 55(3): 778-782.
[10] Herrmann A, Witter E. Sources of C and N contributing to the flush in mineralization upon freeze-thaw cycles in soils. Soil Biology & Biochemistry, 2002, 34(10): 1495-1505.
[11] Qi J L, Vermeer P A, Cheng G D. A review of the influence of freeze-thaw cycles on soil geotechnical properties. Permafrost & Periglacial Processes, 2006, 17(3): 245-252.
[12] Sharma S, Szele Z, Schilling R, Munch C J, Schloter M. Influence of freeze-thaw stress on the structure and function of microbial communities and denitrifying populations in soil. Applied and Environmental Microbiology, 2006, 72(3): 2148-2154.
[13] Luo L F, Robock A, Vinnikov K Y, Schlosser C A, Slater A G, Boone A, Etchevers P, Habets F, Noilhan J, Braden H. Effects of frozen soil on soil temperature, spring infiltration, and runoff: Results from the PILPS 2 (d) experiment at Valdai, Russia. Journal of Hydrometeorology, 2003, 4(2): 334-351.
[14] Froese J C, Cruse R M, Ghaffarzadeh M. Erosion mechanics of soils with an impermeable subsurface layer. Soil Science Society of America Journal, 1999, 63(6): 1836-1841.
[15] Henry H A L. Soil freeze-thaw cycle experiments: Trends, methodological weaknesses and suggested improvements. Soil Biology & Biochemistry, 2007, 39(5): 977-986.
[16] Cruse R M, Mier R, Mize C W. Surface residue effects on erosion of thawing soils. Soil Science Society of America Journal, 2001, 65(1): 178-184.
[17] Liu X J, Jing G C, Qi H Y. The main characteristics of Kebai black soil freezing-thawing erosion. Science & Technology Information of Soil Water Conservation, 1999, (1): 28-30.
[18] Hu G, Wu Y Q, Liu B Y, Zheng Q H, Zhang Y G, Liu H H. The growth characteristics of gully erosion over rolling hilly Black Soil Area of Northeast China. Acta Geographica Sinica, 2007, 62(11): 1165-1173.
[19] Zhang Y G, Wu Y Q, Liu B Y. Ephemeral gully erosion during the spring freeze-thawing period in black soil rolling-hill regions of Northeast. Journal of Mountain Science, 2006, 24(3): 306-311.
[20] Gatto L W. Soil freeze/thaw-induced change to a simulated rill: potential impacts on soil erosion. Geomorphology, 2000, 32(1/2): 147-160.
[21] He J, Liu J L, Lü F. Characterizing fractal characteristics of soil pores based on high-resolution digital CT images. Soils, 2008, 40(4): 662-666.
[22] Jun Y H, Wu Z J, Chen L J, Wang R, Gong L, Bao H J, Liu Y. Preliminary characterization of clay minerals from four typical soils of northeast China. Spectroscopy and Spectral Analysis, 2010, 30(7): 1918-1921.
[23] Gong Z T. Classification System of Chinese Soil Theories, Methods and Practice. Beijing: Science Press, 1999: 474-479.
[24] Chen L X. Experiments and Practices of Soil Course. Harbin: Northeast Forestry University Press, 2005: 43-55.
[25] Qin Y D. Soil Physics. Beijing: Higher Education Press, 2003: 10-20.
[26] Wen M L, Liu B Y, Wei X, Liu H H. Impact of freezing and thawing on bulk density of black soil. Chinese Journal of Soil Science, 2009, 40(3): 492-495.
[27] Wang E H, Zhao Y S, Chen X W. Responses of soil structure to seasonal freeze-thawing cycles in a typical black soil region. Chinese Journal of Applied Ecology, 2010, 21(7): 1744-1750.
[28] Zhao Y S, Wang E H, Cruse R M, Chen X W. Characterization of seasonal freeze-thaw and potential impacts on soil erosion in Northeast China. Canadian Journal of Soil Science, 2012, 92(3): 567-571.
[29] Qiao Q, Shen S M, Zhou S Q. Study on water status of black soil in North part of Northeast China. ⅠBasic characteristics of black soil water status and its relationship with soil genesis. Acta Pedologica Sinica, 1963, 11(2): 143-159.
[30] Qiao Q, Shen S M, Zeng Z S. Study on water status of black soil in North part of North-east China. Ⅱ Agricultural water status and water cycle of black soils. Acta Pedologica Sinica, 1979, 16(4): 329-338.
[31] Meng K, Zhang X Y, Sui Y Y, Zhao J. Black soil water characteristic in Hailun, Heilongjiang. Chinese Journal of Soil Science, 2003, 34(1): 11-14.
[32] Liu J, Fan H M, Zhou L L, Wu M, Chai Y, Liu Y H. Study on effects of freeze-thaw cycle on bulk density and porosity of black soil. Journal of Soil and Water Conservation, 2009, 23(6): 186-189.
[2] 王连峰, 蔡延江, 解宏图. 冻融作用下土壤物理和微生物性状变化与氧化亚氮排放的关系. 应用生态学报, 2007, 18(10): 2361-2366.
[6] 杨思忠, 金会军. 冻融作用对冻土区微生物生理和生态的影响. 生态学报, 2008, 28(10): 5065-5074.
[17] 刘绪军, 景国臣, 齐恒玉. 克拜黑土区沟壑冻融侵蚀主要形态特征初探. 水土保持科技情报, 1999, (1): 28-30.
[18] 胡刚, 伍永秋, 刘宝元, 郑秋红, 张永光, 刘洪鹄. 东北漫岗黑土区切沟侵蚀发育特征. 地理学报, 2007, 62(11): 1165-1173.
[19] 张永光, 伍永秋, 刘宝元. 东北漫岗黑土区春季冻融期浅沟侵蚀. 山地学报, 2006, 24(3): 306-311.
[21] 何娟, 刘建立, 吕菲. 基于CT数字图像的土壤孔隙分形特征研究. 土壤, 2008, 40(4): 662-666.
[22] 隽英华, 武志杰, 陈利军, 汪仁, 宫亮, 包红静, 刘艳. 东北4种典型土壤粘粒矿物的初步表征. 光谱学与光谱分析, 2010, 30(7): 1918-1921.
[23] 龚子同. 中国土壤系统分类: 理论、方法和实践. 北京: 科学出版社, 1999: 474-479.
[24] 陈立新. 土壤实验实习教程. 哈尔滨: 东北林业大学出版社, 2005: 43-55.
[25] 秦耀东. 土壤物理学. 北京: 高等教育出版社, 2003: 10-20.
[26] 温美丽, 刘宝元, 魏欣, 刘洪鹄. 冻融作用对东北黑土容重的影响. 土壤通报, 2009, 40(3): 492-495.
[27] 王恩姮, 赵雨森, 陈祥伟. 典型黑土耕作区土壤结构对季节性冻融的响应. 应用生态学报, 2010, 21(7): 1744-1750.
[29] 乔樵, 沈善敏, 周绍权. 东北北部黑土水分状况之研究 Ⅰ.黑土水分状况的基本特征及其与成土过程的关系. 土壤学报, 1963, 11(2): 143-159.
[30] 乔樵, 沈善敏, 曾昭顺. 东北北部黑土水分状况之研究 Ⅱ.黑土农业水分状况及水分循环. 土壤学报, 1979, 16(4): 329-338.
[31] 孟凯, 张兴义, 随跃宇, 赵军. 黑龙江海伦农田黑土水分特征. 土壤通报, 2003, 34(1): 11-14.
[32] 刘佳, 范昊明, 周丽丽, 武敏, 柴宇, 刘艳华. 冻融循环对黑土容重和孔隙度影响的试验研究. 水土保持学报, 2009, 23(6): 186-189.