陳 浮楊寶丹馬 靜張紹良張媛媛

（1 中國礦業(yè)大學(xué)環(huán)境與測繪學(xué)院，江蘇徐州 221116）

（2 中國礦業(yè)大學(xué)低碳能源研究院，江蘇徐州 221008）

（3 中石化石油工程設(shè)計有限公司，山東東營 257026）

高濃度CO2地下泄漏對土壤微生物群落結(jié)構(gòu)的影響*

陳 浮1，2楊寶丹1馬 靜2?張紹良2張媛媛3

（1 中國礦業(yè)大學(xué)環(huán)境與測繪學(xué)院，江蘇徐州 221116）

（2 中國礦業(yè)大學(xué)低碳能源研究院，江蘇徐州 221008）

（3 中石化石油工程設(shè)計有限公司，山東東營 257026）

通過模擬高濃度CO2在農(nóng)田土壤中的地下泄漏，研究了不同濃度CO2泄漏情景下土壤微生物多樣性的變化。實驗設(shè)置了400 g m-2d-1、800 g m-2d-1、1 200 g m-2d-1和2 000 g m-2d-1持續(xù)CO2通氣60 d共計4個處理，并與對照組、恢復(fù)組（2 000 g m-2d-1組停止通氣60 d后）分期采集土壤樣品，分析土壤理化性質(zhì)、土壤閉蓄的氣體濃度、微生物多樣性指數(shù)及主要類群變化規(guī)律。結(jié)果表明，4種處理均提高了土壤中CO2濃度，分別為1.60%、4.80%、10.80%和19.60%。土壤微生物多樣性Chao指數(shù)和Shannon指數(shù)隨CO2通入量增加而減少，降幅分別達(dá)17.00%～27.80%和6.10%～9.50%。相反，非度量多維尺度（NMDS）分析顯示土壤微生物β多樣性在中、低濃度升高，在高、極端濃度表現(xiàn)為降低。擬桿菌屬（Bacteroidales）相對豐度隨CO2泄漏量增加從3.09%上升至21.20%，可作為高濃度CO2泄漏生態(tài)安全性評估的敏感性指標(biāo)?；诟咄啃蛄邢嗨贫萇TU分類的RDA分析表明土壤環(huán)境因子的變化能夠較好地解釋微生物多樣性演替。研究結(jié)果為評估和監(jiān)測地下CO2泄漏對近地表生態(tài)系統(tǒng)環(huán)境風(fēng)險提供科學(xué)依據(jù)。

碳捕集與封存；CO2泄漏；微生物群落；高通量測序；多樣性分析

溫室氣體減排已經(jīng)成為各國政府高度關(guān)注的議題和目標(biāo)，CO2捕集與封存（CCS，CO2Capture and Geological Storage）是一項重要的減排技術(shù)，因其經(jīng)濟(jì)與環(huán)境雙贏的特點(diǎn)而倍受關(guān)注［1］。盡管CCS工程一直將防止CO2泄漏列為首要條件，但現(xiàn)有的技術(shù)體系還很難保證封存的CO2完全不發(fā)生泄漏，例如，儲運(yùn)、地質(zhì)活動或人為操作失誤等均可能導(dǎo)致地下封存的CO2發(fā)生泄漏風(fēng)險［2］。1986年喀麥隆尼奧斯湖沉積相CO2暴發(fā)導(dǎo)致周邊數(shù)十平方千米范圍內(nèi)2 000多人、數(shù)萬牲畜窒息死亡［3］。目前中國已開展了16項CCS試驗，其中鄂爾多斯咸水層封存和勝利油田驅(qū)油開采均為1 Mt級CO2封存項目，一旦發(fā)現(xiàn)泄漏極可能對近地表生態(tài)系統(tǒng)產(chǎn)生嚴(yán)重的負(fù)面影響［4］?，F(xiàn)行的CCS工程建設(shè)和安全評估均側(cè)重于封存選址和CO2地下運(yùn)移［5-8］，更多關(guān)注泄漏的環(huán)境后果和評估。然而，由于CO2地下封存技術(shù)涉及的因素較多，如地質(zhì)環(huán)境、封存材質(zhì)和CO2含量等，很難采用統(tǒng)一的標(biāo)準(zhǔn)進(jìn)行CCS環(huán)境風(fēng)險評估，如近幾年歐美國家已經(jīng)多次出現(xiàn)CCS項目環(huán)境侵權(quán)訴訟［9］。因此，在我國開展高濃度CO2地下泄漏的生態(tài)環(huán)境影響研究，具有重要意義。

已有的研究表明，由于深部高壓和浮力作用，CO2從斷層、裂隙、廢棄通道等逐步向上遷移，改變淺層地下水pH［10］，影響礦化物的平衡控制、遷移轉(zhuǎn)化和水環(huán)境［11］，并進(jìn)一步對上層土壤-植物生態(tài)系統(tǒng)產(chǎn)生直接或間接的負(fù)面影響［12］。英國ASGARD項目組、美國ZERT研究組、歐盟Geo-Net和中國農(nóng)業(yè)科學(xué)院等機(jī)構(gòu)已實施了多項人工控制地下CO2釋放的模擬試驗，探測農(nóng)田土壤-植被生態(tài)系統(tǒng)對地下CO2泄漏的響應(yīng)狀況，但取得的結(jié)論并不一致［13-20］。主要原因可歸于三點(diǎn)：一是試驗平臺不一，有盆栽，有小塊實驗田，土壤中CO2擴(kuò)散方式和濃度難以精確測定；二是持續(xù)暴露的時間不一，短的約為30 d，最長的則達(dá)90 d，很難明確土壤微生物群落轉(zhuǎn)變的臨界點(diǎn)；三是種植體系不一，不同根系分泌物對土壤微生物群落影響差異巨大，不同研究之間缺乏可比性。此外，以往的大氣CO2濃度升高對生態(tài)系統(tǒng)影響的研究很難作為CCS環(huán)境風(fēng)險評估的參考。主要原因包括兩點(diǎn)，首先是作用途徑不同，大氣CO2濃度倍增實驗中其濃度僅僅8.0?，從地表向土壤中滲透，土壤中CO2濃度一般不高于1.00%，因此，外源CO2對土壤本身CO2氣體濃度的影響較小，通常變化幅度小于10%，但地下CO2泄漏則由深層土壤向地表逃逸，其濃度可高達(dá)20.00%，導(dǎo)致土壤中CO2濃度增加數(shù)倍甚至數(shù)十倍；二是持續(xù)時間不同，地表CO2受空氣流動影響，極容易發(fā)生擴(kuò)散，因此多采用封閉栽培箱，對土壤氣的影響相對小，而地下泄漏CO2在土壤保存時間很長，其對土壤性質(zhì)包括微生物多樣性的影響可能更加劇烈。

據(jù)此，本研究在中國礦業(yè)大學(xué)南湖校區(qū)內(nèi)建立了大田實驗平臺，模擬高濃度CO2地下泄漏情景，研究土壤微生物群落結(jié)構(gòu)對CO2地下泄漏的響應(yīng)，為評估和監(jiān)測地下CO2泄漏對近地表生態(tài)系統(tǒng)環(huán)境風(fēng)險提供科學(xué)依據(jù)。

1 材料與方法

1.1 試驗場地

本試驗平臺在中國礦業(yè)大學(xué)南湖校區(qū)內(nèi)，地理位置為117°08′08″E、34°12′31″N，屬暖溫帶半濕潤季風(fēng)氣候，年均溫度14.2 ℃，年均降水量833 mm。土壤類型為潮褐土，陽離子交換量約為14.91 ± 3.27 cmol kg-1，質(zhì)地為輕壤，機(jī)械組成（＜0.001 mm、0.001～0.005、0.005～0.01、＞0.01 mm）為15.12∶22.34∶9.15∶53.39。2012 年9月利用鏟車分0～20 cm、20～90 cm二層取土并分別堆放，由下至上平鋪10 cm粗砂、布設(shè)內(nèi)帶PVC軟管的多孔鋼管、10 cm細(xì)沙和1 cm多孔納米板，最后再分別鋪上50 cm下層土壤和20 cm表土，剩余土鋪設(shè)田間路［21］。盡量保持原土壤構(gòu)型，簡單平整后放置1.5 a自然壓實。實驗田總面積約為1 000 m2（40 m×25 m），共建成48個2.5 m×2.5 m標(biāo)準(zhǔn)樣方（圖1）。

1.2 試驗設(shè)計

試驗平臺包含自動控制系統(tǒng)室、地下管路、田間路、水利設(shè)施和田地。本研究針對15個樣方開展，首先清除雜草和根系，種植玉米。2014年3月7日開始持續(xù)通氣60 d，氣源為工業(yè)級CO2，體積濃度為99.8%。設(shè)計低濃度（L-400，400 g m-2d-1）、中濃度（M-800，800 g m-2d-1）、高濃度（H-1200，1 200 g m-2d-1）、極端高濃度（E-2000，2 000 g m-2d-1）和對照組（C-0），各設(shè)3個重復(fù)。通氣后采3次土樣，分別在14 d、30 d 和60 d。對照組通氣前（C1）和試驗結(jié)束時（C2）各采1次?；謴?fù)組（R-0）為極端高濃度E組停止CO2通氣60 d后采樣1次。通入量由氣體調(diào)節(jié)閥和流量控制器控制，按下列公式計算注入速率：

V=FS / ρ （1）

式中，V為流量控制器上CO2注入速率（ml s-1）；F為設(shè)計注入濃度（400 g m-2d-1）；S為樣方面積（m-2）；ρ為CO2密度（g L-1）。

1.3 樣品采集與分析

采用五點(diǎn)法采集10～20 cm混合土樣。一部分室內(nèi)自然風(fēng)干，過2 mm篩剔除石礫和植物殘體等雜物，4 ℃冷藏備用。另一部分無菌袋密封，存放-20 ℃冰箱，用于微生物測試。

土壤有機(jī)質(zhì)采用重鉻酸鉀-外加熱法測定，陽離子交換量采用EDTA-銨鹽快速法測定，機(jī)械分析采用比重計法測定［22］，NO3--N采用雙波長分光光度法測定，土壤pH采用電位法測定（土：水＝1∶2.5）［23］。

土壤中CO2氣體濃度測量：每個樣方內(nèi)均垂直埋深20 cm氣體探針（長40 cm，內(nèi)徑2 cm），底端密封，最下端10 cm管壁上均勻鉆直徑0.3 cm的小孔，外貼防水生物膜，保證氣體自由擴(kuò)散又不積水。氣體探針頂端設(shè)密封閥門，連接到改制的帶自動傳輸數(shù)據(jù)功能的便攜式氣體檢測儀GT901 （TES，臺灣）上測量土壤中CO2體積濃度，量程0～30.00%，數(shù)據(jù)傳輸間隔為1 h。

圖1 試驗平臺位置與建設(shè)方案Fig. 1 Location and construction scheme of the experimental platform

DNA提取、擴(kuò)增與測序：按DNA快速提取試劑盒（Felix bio-tech，美國）上操作說明從45個（每個0.5 g）土樣中提取總DNA，利用核酸定量儀（NanoDrop ND-1000，美國）檢測提取的總DNA濃度和純度。采用特異引物338F （5′-ACTCCTACGGGAGGCAGCAG-3′）和806R （5′-GGACTACVSGGGTATCTAAT-3′）對16S rDNA基因的V4區(qū)進(jìn)行PCR擴(kuò)增（Pyrobest DNA Polymerase，TaKaRa，DR500A）。25 μlPCR體系包含：5 μl 5×Q5 反應(yīng)緩沖劑，5 μl 5×Q5 GC 強(qiáng)化劑，2 μl 2.5 mmol L-1dNTPs，兩個特異引物各1 μl（10 μmol L-1），0.25 μl Q5 聚合酶，8.75 μl超純無菌水和2 μl 模板DNA（20 ng μl-1）。PCR擴(kuò)增具體程序如下：98 ℃預(yù)變性5 min，98 ℃變性30 s，50 ℃退火30 s，72 ℃延伸30 s，共持續(xù)25個循環(huán)周期。最后72 ℃持續(xù)10 min結(jié)束程序。進(jìn)行1%瓊脂糖凝膠電泳檢測，割膠回收目標(biāo)條帶（Axygen Biosciences，Union City，CA，USA），Tris-HCL洗脫。最終利用AxyPrep DNA凝膠抽取儀對PCR進(jìn)行定量（Axygen Biosciences，Union City，CA，USA）。利用PCR選擇性地富集兩端連有接頭的DNA片段，同時擴(kuò)增DNA文庫。利用PicoGreen和熒光分光光度計方法定量文庫。此外，使用Agilent 2100對PCR富集片段進(jìn)行質(zhì)量控制，驗證DNA文庫的片段大小及分布。將所有樣品DNA文庫均一化至10 nmol L-1后等體積混合。再將混合好的文庫（10 nmol L-1）逐步稀釋定量至4～5 pmol L-1后進(jìn)行上機(jī)測序。

1.4 數(shù)據(jù)處理

首先需要去除引物接頭序列，再根據(jù)PE reads之間的overlap關(guān)系，將成對的reads 拼接（merge）成一條序列，然后按照barcode標(biāo)簽序列識別并區(qū)分樣品得到各樣本數(shù)據(jù)，最后對各樣本數(shù)據(jù)的質(zhì)量進(jìn)行質(zhì)控過濾，得到各樣本有效數(shù)據(jù)。利用QIIME（Version 1.7.0，http：//qiime.org/）對原始FASTQ文件進(jìn)行多路復(fù)用剔除和質(zhì)量篩選，調(diào)用uclust對優(yōu)質(zhì)序列按序列相似度97%進(jìn)行聚類分析，選取每個類中最長的序列作代表序列。再將優(yōu)化序列map比對回OTU代表序列，獲得每個OTU代表序列的分類學(xué)信息。選取相似度在97%條件下OTU生成預(yù)期的稀釋曲線，利用Chao（The Chao1 estimator，http：//www.mothur.org/wiki/Chao）和Shannon（The Shannon index，http：//www.mothur. org/ wiki/Shannon）指數(shù)反映單樣本（一個特定區(qū)域或生態(tài)系統(tǒng)）群落的豐富度和均勻性，即α多樣性；利用非度量多維標(biāo)度指數(shù)（NMDS index）計算各樣本間的進(jìn)化及豐度信息的距離來反映是否具有顯著的群落差異，即β多樣性；利用SPSS17.0 和Canoco windows for 4.5（http：//www.canodraw. com/）軟件進(jìn)行數(shù)據(jù)統(tǒng)計和冗余分析（RDA）。多重相關(guān)比較顯著性檢驗采用最小顯著差法（LSD）。

2 結(jié) 果

2.1 高濃度CO2地下泄漏對土壤理化性質(zhì)的影響

從表1可以看出，持續(xù)通氣60 d后土壤化學(xué)特性發(fā)生顯著變化。土壤pH、電導(dǎo)率和NO3--N隨CO2泄漏量增加而下降，pH和NO3--N尤為顯著。停止通氣后pH則恢復(fù)明顯，NO3--N略有回升，但電導(dǎo)率則繼續(xù)下降，有機(jī)質(zhì)含量卻隨CO2泄漏量呈上升趨勢，停止通入后明確回落。統(tǒng)計分析發(fā)現(xiàn)CO2泄漏量與主要土壤因子變化存在顯著性相關(guān)關(guān)系（表1）。

2.2 泄漏對土壤中CO2濃度的影響

根據(jù)上述公式計算，L-400、M-800、H-1200、E-2000對應(yīng)的注入速率分別為16.64 ml s-1、29.27 ml s-1、43.91 ml s-1和73.18 ml s-1。樣品主要采集10～20 cm表層土壤，因此CO2傳感器探針設(shè)置監(jiān)測位置為埋深10～20 cm處。從圖2可以看出，E-2000組通氣6 h后土壤中CO2濃度趨向穩(wěn)定，約為20.00%，H-1200組通氣8 h后趨向穩(wěn)定，約為10.06%，24 h內(nèi)各組CO2濃度均呈現(xiàn)飽和，C-0組一直保持約0.60%，組內(nèi)差異不顯著。其中的主要原因是，土壤中CO2濃度受通氣量、氣體來源、地表氣候因子、通氣時間、土壤孔隙度、土壤含水率和pH等多因素的影響，本模擬平臺設(shè)置的CO2泄漏呈近似均勻擴(kuò)散，與先前模擬試驗呈同心圓點(diǎn)源擴(kuò)散完全不同，保證土壤中CO2濃度只受土壤自身狀況的影響［24］。當(dāng)通氣一段時間后土壤中CO2濃度達(dá)到飽和，持續(xù)通氣則向地表逸散，并不能持續(xù)增加土壤中CO2濃度。

表1 不同處理各時期土壤理化特征變化Table 1 Soil physicochemical properties relative to treatment and timing of sampling

2.3 土壤微生物群落多樣性和豐富度變化

全部45個樣品采用Illumina Miseq測序技術(shù)共獲取了2 121 986個序列，所有樣品的稀釋曲線均傾向于飽和平穩(wěn)狀態(tài)，表明測序數(shù)據(jù)量合理。圖3a顯示Chao和Shannon指數(shù)隨CO2泄漏量和時間的變化，2個指數(shù)均隨CO2泄漏量增加而減小，L-400、M-800組與C-0對照組變化小，說明中、低濃度CO2通入對土壤微生物影響較小。E-2000組與C-0對照組相比二個指數(shù)最大分別下降了17.00%～27.80%和6.10%～9.50%，說明高、極高濃度CO2通入對微生物影響顯著。CO2泄漏量800 g m-2d-1可能是微生物的一個重要閾值，更高的濃度可能徹底改變了微生物對土壤的環(huán)境適應(yīng)性。與此同時，停止CO2通入60 d，微生物群落α多樣性和豐富度均呈現(xiàn)大幅度反彈，說明土壤微生物對CO2脅迫具備較強(qiáng)的耐受性。

圖2 土壤中CO2濃度隨通氣時間的變化Fig. 2 Variation of CO2concentration in the soil with aeration going on

在O T U分類基礎(chǔ)上采用非度量多維尺度（NMDS）指數(shù)作β多樣性分析。NMDS指數(shù)顯示C-0、L-400、和R-0樣品和其他樣品群落區(qū)分開，第2軸又將M-800和H-1200、E-2000組樣品分開（圖3b）。NMDS指數(shù)更好地解釋了H-1200、 E-2000組高濃度下“抑制狀態(tài)”與L、M組低濃度下“促進(jìn)狀態(tài)”下土壤微生物群落結(jié)構(gòu)的變化。

2.4 土壤微生物結(jié)構(gòu)與組成變化

所有樣品中95%以上的序列可以明確分類至不同類群，包括酸桿菌門（Acidobacteria）、放線菌門（Actinobacteria）、擬桿菌門（Bacteroidetes）、綠彎菌門（Chloroflexi）、厚壁菌門（Firmicutes）、浮霉菌門（Planctomycetes）、變形菌門（Proteobacteria）、芽單胞菌門（Gemmatimonadetes）和疣微菌門（Verrucomicrobia）等9大菌門相對豐度占90%以上，其中放線菌門（Actinobacteria）最為豐富，約占35.51%（圖4）。C O2泄漏造成了厭氧環(huán)境，酸桿菌門（Acidobacteria）、綠彎菌門（Chloroflexi）和芽單胞菌門（Gemmatimonadetes）等好氧菌門相對豐度隨C O2泄漏量呈下降趨勢，擬桿菌門（Bacteroidetes）和厚壁菌門（Firmicutes）等厭氧菌門則呈相反的趨勢。其中擬桿菌屬（Bacteroidales）變化最為顯著，隨CO2通入從C-0處理開始時相對豐度3.0 9%上升至E-2000處理60 d后相對豐度21.20%，乳酸菌屬（Lactobacillus）次之，從C-0處理開始時相對豐度1.10%增加至E-2000處理60 d后相對豐度6.95%。熱微菌屬（Thermomicrobia）則隨CO2通入從C-0處理開始時相對豐度16.70%下降至E-2000處理60 d后相對豐度3.20%。

圖3 不同處理微生物群落α多樣性（a）與β多樣性（b）分析Fig. 3 Alpha diversity（a）and beta diversity（b）analysis of microbial community relative to treatment

圖4 不同處理的土壤微生物菌門相對豐度變化Fig. 4 Phylum composition of the soil microbes relative to treatment

為進(jìn)一步檢測土壤環(huán)境因子與微生物群落變化的對應(yīng)關(guān)系，引入了RDA分析。從圖5可以看出：不同處理各時期土壤微生物群落與pH、EC值和NO3--N等環(huán)境因子呈現(xiàn)顯著相關(guān)（p＜0.005），在第1軸上解釋了C-0、L-400組與M-800、H-1200、E-2000組之間的差異，pH是最長的箭頭，表示其是引起群落分異的最重要因子。有機(jī)質(zhì)在第2軸上較好地解釋了L-400、M-800組與H-1200、E-2000組之間群落差異。

圖5 土壤微生物群落結(jié)構(gòu)與土壤因子RDA分析Fig. 5 RDA analysis of soil microbial community with soil factors

3 討 論

CO2濃度、根系分泌物和時間通常被認(rèn)為是影響土壤微生物多樣性的3個重要變量［25］，本研究在考慮上述變量基礎(chǔ)上設(shè)計了4個通氣梯度研究高濃度CO2地下泄漏對土壤微生物的脅迫，結(jié)果發(fā)現(xiàn)多樣性、豐富度和群落結(jié)構(gòu)均發(fā)生了變化，這與前人CO2濃度倍增實驗取得的結(jié)論相對一致［26-29］，但變化特征卻具有明顯的差異。理論而言，地下CO2泄漏對土壤和農(nóng)作物的影響完全不同于以往的大氣CO2濃度倍增實驗，前者首先影響土壤環(huán)境，土壤CO2濃度飽和后再逸散入大氣，同時在近地表受氣流影響很難形成高濃度的聚集［30］。因此，實際狀況下地下CO2泄漏對農(nóng)作物幾乎不可能產(chǎn)生“施肥”效應(yīng)，為此精準(zhǔn)的實驗?zāi)M平臺至關(guān)重要。以往模擬試驗均采用點(diǎn)源，尤其是盆栽實驗幾乎不考慮泄漏與土壤中CO2濃度的關(guān)系，一般情況下點(diǎn)源擴(kuò)散在土壤中CO2濃度呈同心圓狀［16-18］，但環(huán)狀條帶極不穩(wěn)定，易造成多重濃度反復(fù)交叉干擾。這與地下CO2泄漏以及大氣CO2濃度倍增實驗緩慢向上或向下均勻的影響特征完全不同，每次采樣均需要重復(fù)測定土壤氣濃度，耗時費(fèi)力精度差。后續(xù)植被實驗更難，往往一個濃度條帶上僅能生長一株玉米，樣品難以滿足實驗要求。本模擬實驗平臺解決了泄漏與土壤中CO2濃度的關(guān)聯(lián)性，獲取均勻態(tài)擴(kuò)散，極大地方便了研究［24］。CO2注入加大了土壤中CO2濃度，可能改變了土壤物理化學(xué)性質(zhì)，導(dǎo)致微生物群落結(jié)構(gòu)的變化。一般認(rèn)為CO2等酸氣對土壤pH屬性影響最為明顯影響，本實驗結(jié)果顯示土壤pH隨CO2通入濃度梯度降低而逐漸下降，這與Beaubien［26］和Oppermann［31］等關(guān)于泄漏點(diǎn)pH呈下降趨勢的結(jié)論一致，但McFarland等［32］并未發(fā)現(xiàn)極端CO2濃度變化會引起pH的明顯變化。造成這一差異的原因可能來自兩個方面：一是與土壤本底pH有關(guān)，CO2泄漏不可能造成土壤無限的酸化；二是當(dāng)土壤中CO2達(dá)到飽和后，大量CO2外排并逸散入大氣，換言之，大氣中CO2濃度即使再高，也不可能無限擴(kuò)大土壤中CO2濃度，這與地下泄漏直接作用于土壤環(huán)境的影響完全不同。此外，本實驗證實了Luis等［33］關(guān)于土壤有機(jī)質(zhì)隨CO2通入呈上升的趨勢，表明60 d高濃度CO2泄漏能改變土壤有機(jī)碳庫，這可能與還原條件下有機(jī)質(zhì)轉(zhuǎn)化率更低密切相關(guān)。同時發(fā)現(xiàn)NO3--N也是隨著CO2泄漏量急劇下降的，可能是CO2通入促進(jìn)反硝化菌生長，利用硝態(tài)氮參加反硝化反應(yīng)有關(guān)［34］，然而，已知的絕大多數(shù)反硝化菌均為異養(yǎng)微生物，利用有機(jī)碳而非CO2生長，因此，也有可能是其他土壤性質(zhì)的變化，例如土壤酸化導(dǎo)致大量有機(jī)碳溶出，進(jìn)而刺激了反硝化微生物生長，或者是土壤理化性質(zhì)的變化，導(dǎo)致土壤硝態(tài)氮快速淋溶至深層土壤。但也有相反的研究結(jié)果，如McFarland等［32］認(rèn)為NO3--N會隨CO2增加而增大，其原因可能是CO2通入可能會抑制生物的生長，減少植被和土壤微生物對氮的需求，但具體的機(jī)制仍需進(jìn)一步研究，并在不同CO2封存條件下具有不同的響應(yīng)規(guī)律。

根系分泌物是影響微生物群落結(jié)構(gòu)及多樣性的重要因子［27，35］，但相互作用關(guān)系極為復(fù)雜。本實驗前剔除了雜草和根系，再統(tǒng)一種植玉米，目的即是消除不同植物根系分泌物對微生物影響的不確定性，同一植物根系且初始土壤因子基本無差異，CO2注入則成為影響微生物的單一因子。已有報道表明高濃度CO2泄漏會降低土壤微生物多樣性和豐富度，本研究也得到了類似結(jié)果［3，25，36］。不同處理各時期群落特征存在差異，說明CO2泄漏已經(jīng)對土壤微生物多樣性產(chǎn)生了重要影響。2014年5月勝利油田G89-1驅(qū)油開采區(qū)生態(tài)本底值調(diào)查時發(fā)現(xiàn)了數(shù)處地下CO2泄漏點(diǎn)，連續(xù)觀察了3 a，最高一處近地表CO2濃度高達(dá)0.20%，是正常大氣CO2濃度值的5倍，泄漏點(diǎn)周邊1.5 m范圍內(nèi)植物均無法成活。與天然CO2通道周邊觀察結(jié)果極為相似［26］，影響范圍小主要原因是CO2在風(fēng)速大于1 m s-1時就能迅速擴(kuò)散［30］，一旦發(fā)生大量泄漏則破壞力極大［3］。其外，本研究發(fā)現(xiàn)玉米根系下擬桿菌目（Bacteroidales）相對豐度隨CO2通氣量持續(xù)顯著增加，但這類群的新陳代謝及功能多樣性還需要更充分的信息加以分析，這是未來研究的一個重要方向。

盡管很難解釋高濃度CO2泄漏條件下一些微生物的增加或減少，但高通量測序技術(shù)提供了不同類群微生物的動態(tài)變化趨勢［18，37］，可在一定程度較為明確地解釋與CO2通入量相關(guān)的土壤環(huán)境下微生物群落的差異性。此外，不同處理組類群中序列分布是不同的，它表明時間和CO2通入量均產(chǎn)生了顯著的影響，R-0組群落α多樣性和豐富度更為接近C-0對照組，但反映群落結(jié)構(gòu)和組成的β多樣性則更接近H-1200組和E-2000組，它顯示土壤微生物盡管可以適應(yīng)新的環(huán)境體系，但其功能的恢復(fù)需要更長的時間，不同于α多樣性和豐富度的恢復(fù)。這些研究結(jié)果表明，高濃度CO2地下泄漏極可能對土壤微生物群落及其功能產(chǎn)生重要影響，未來仍需對土壤微生物群落的底物利用方式，生長代謝速率及其恢復(fù)力對CO2泄漏的適應(yīng)規(guī)律進(jìn)一步研究。

4 結(jié) 論

不同CO2泄漏情景下農(nóng)田土壤微生物群落結(jié)構(gòu)、多樣性和豐富度存在顯著性差異。Chao指數(shù)和Shannon指數(shù)隨CO2泄漏量增加而減少，降幅分別達(dá)17.00%～27.80%和6.10%～9.50%。同時，微生物類群隨CO2泄漏梯度增加發(fā)生了明顯的異化。β多樣性分析表明不同CO2泄漏情景影響了農(nóng)田土壤微生物群落組成，土壤pH和CO2通入量可能是影響微生物多樣性的最重要環(huán)境因子。擬桿菌屬（Bacteroidales）相對豐度隨CO2泄漏量增加從3.09%上升至21.20%，可作為監(jiān)測和評估CO2地下泄漏生態(tài)環(huán)境風(fēng)險的敏感性指標(biāo)。盡管60 d恢復(fù)實驗表明土壤微生物群落多樣性和豐富度恢復(fù)較好，但其功能恢復(fù)力尚未可知，未來仍需重點(diǎn)研究高濃度CO2地下泄漏對土壤生態(tài)系統(tǒng)功能的潛在風(fēng)險。

參 考 文 獻(xiàn)

［1］Dai Z X，Stauffer P H，Carey J W，et al. Presite characterization risk analysis for commercialscale carbon sequestration. Environmental Science & Technology，2014，48（7）：3908—3915

［2］任韶然，李德祥，張亮，等. 地質(zhì)封存過程中泄漏途徑及風(fēng)險分析.石油學(xué)報，2014，35（3）：591—601

Ren S R，Li D X，Zhang L，et al. Leakage pathways and risk analysis of carbon dioxide in geological storage （In Chinese）. Acta Petrolei Sinica，2014，35 （3）：591—601

［3］Patil R H，Colls J J，Steven M D. Effects of CO2gas as leaks from geological storage sites on agro-ecosystems. Energy，2010，35（12）：4587—4591

［4］IPCC. Climate change 2014：Mitigation of climate change//Contribution of working group III to the fifth assessment report of the intergovernmental panel on climate change. Cambridge，United Kingdom and New York，USA：Cambridge University Press，2014：18—28

［5］Bachu S，Bonijoly D，Bradshaw J，et al. CO2storage capacity estimation：Methodology and gaps. International Journal of Greenhouse Gas Control，2007，1（4）：430—443

［6］Gibson-Pool C M，Svendsen L，Underschultz J，et al. Site characterization of a basin-scale CO2geological storage system：Gippsland basin，Southeast Australia. Environmental Geology，2008，54（1）：1583—1606

［7］Okwen R T，Stewart M T，Cunninghan J A. Analytical model for screening potential CO2repositories. Computational Geosciences，2011，15（4）：755—770

［8］Wielopolski L，Mitra S. Near-surface soil carbon detection for monitoring CO2seepage from a geological reservoir. Environmental Earth Sciences，2010，60 （2）：307—312

［9］Wilson E J，Mark A，F(xiàn)igueiredo D. Geologic carbon dioxide sequestration：An analysis of subsurface property law. Journal of Environmental Law Report，2006，36：l14—121

［10］Lemieux J M. Review：The potential impact of underground geological storage of carbon dioxide in deep saline aquifers on shallow groundwater resources. Hydrogeology Journal，2011，19（4）：757—778

［11］Kharaka Y K，Thordsen J J，Kakouros E，et al. Changes in the chemistry of shallow groundwater related to the 2008 injection of CO2at the ZERT field site，Bozeman，Montana. Environmental Earth Sciences，2010，60（2）：273—284

［12］李琦，劉桂臻，張建，等. 二氧化碳地質(zhì)封存環(huán)境監(jiān)測現(xiàn)狀及建議. 地球科學(xué)進(jìn)展，2013，28（6）：718—727 Li Q，Liu G Z，Zhang J，et al. Status and suggestion of environmental monitoring for CO2geological storage （In Chinese）. Advances in Earth Science，2013，28 （6）：718—727

［13］West J M，Jones D G，Annunziatellis A，et al. Comparison of the impacts of elevated CO2soil gas concentrations on selected European terrestrial environments. International Journal of Greenhouse Gas Control，2015，42：357—371

［14］Fessen J E，Clegg S M，Rahn T A，et al. Novel MVA tools to track CO2seepage tested at the ZERT controlled release site in Bozeman，MT. Environmental Earth Sciences，2010，60（2）：325—334

［15］Wells A，Strazisar B，Diehl J R，et al. Atmospheric tracer monitoring and surface plume development at the ZERT pilot test in Bozeman，Montana，USA. Environmental Earth Sciences，2010，60（2）：299—305

［16］Oldenburg C M，Lewicki J L，Pan L，et al. Origin of the patchy emission pattern at the ZERT CO2release test. Environmental Earth Sciences，2010，60（2）：241—250

［17］伍洋，馬欣，李玉娥，等. 地質(zhì)封存CO2泄漏對農(nóng)田生態(tài)系統(tǒng)的影響評估及耐受閾值.農(nóng)業(yè)工程學(xué)報，2012，28（2）：196—205

Wu Y，Ma X，Li Y E，et al. Impact assessment and tolerable threshold value of CO2leakage from geological storage on agro-ecosystem（In Chinese）. Transactions of the CSAE，2012，28（2）：196—205

［18］田地，馬欣，李玉娥，等.利用高通量測序?qū)Ψ獯鍯O2泄漏情景下土壤微生物的研究. 環(huán)境科學(xué)，2013，34 （10）：4096—4104

Tian D，Ma X，Li Y E，et al. Research on soil bacteria under the impact of sealed CO2leakage by high-throughput sequencing technology（In Chinese）. Environmental Science，2013，34（10）：4096—4104

［19］蔣金豹，Michael D S，何汝艷，等. 利用大豆光譜特征判定地下封存CO2泄漏. 農(nóng)業(yè)工程學(xué)報，2013，29 （12）：163—169

Jiang J B，Michael D S，He R Y，et al. Judgment of CO2leaking in underground storage using spectral characteristics of soybean（In Chinese）. Transactions of the CSAE，2013，29（12）：163—169

［20］聶莉娟，馬俊杰，趙雪峰，等. 模擬CCS技術(shù)CO2泄漏對C3、C4作物土壤化學(xué)性質(zhì)的影響. 水土保持學(xué)報，2015，29（5）：200—205Nie L J，Ma J J，Zhao X F，et al. Effects of simulation CO2leakage of CCS on soil chemical properties under C3 and C4 crops（In Chinese）. Journal of Soil and Water Conservation，2015，29（5）：200—205

［21］Chen F，Zhang S L，Ma J，et al. A set of efficient device simulating the potential impact of CO2leakage from CO2enhanced oil recovery on farmland ecosystem// The 14th annual carbon capture，utilization & storage conference. Pittsburgh：Access Intelligence Press，2015：81

［22］張紅，呂家瓏，曹瑩菲，等. 不同植物秸稈腐解特性與土壤微生物功能多樣性研究. 土壤學(xué)報，2014，51 （4）：743—752

Zhang H，Lü J L，Cao Y F，et al. Decomposition characteristics of different plant straws and soil microbial functional diversity（In Chinese）. Acta Pedologica Sinica，2014，51（4）：743—752

［23］楊帆，黃來明，李德成，等. 高寒山區(qū)地形序列土壤有機(jī)碳和無機(jī)碳垂直分布特征及其影響因素. 土壤學(xué)報，2015，52（6）：1226—1236

Yang F，Huang L M，Li D C，et al. Vertical distributions of soil organic and inorganic carbon and their controls along toposequences in an alpine region （In Chinese）. Acta Pedologica Sinica，2015，52 （6）：1226—1236

［24］陳浮，譚敏，馬靜，等. 一種模擬CO2在土壤中均勻擴(kuò)散的試驗裝置與方法.中國專利：201510375957.9. 2015-08-06

Chen F，Tan M，Ma J，et al. A Set of devices and methods simulating CO2homogeneous diffusion in soil （In Chinese）. China Patent：201510375957.9. 2015-08-06

［25］Fernández M I，Touceda M，Pedescoll A，et al. Shortterm effects of simulated below-ground carbon dioxide leakage on a soil microbial community. International Journal of Greenhouse Gas Control，2015，36：51—59

［26］Beaubien S E，Ciotoli G，Coombs P，et al. The impact of a naturally occurring CO2gas vent on the shallow ecosystem and soil chemistry of a Mediterranean pasture （Latera，Italy）. International Journal of Greenhouse Gas Control，2008，3（2）：373—387

［27］West J M，Pearce J M，Coombs P，et al. The impact of controlled injection of CO2on the soil ecosystem and chemistry of an English lowland pasture. Energy Procedia，2009，1（1）：1863—1870

［28］Antruck H，Imek M. Effect of soil CO2concentration on microbial biomass. Biology and Fertility of Soils，1997，25（3）：269—273

［29］Morales S E，Holben W E. Simulated geologic carbon storage leak reduces bacterial richness and alters bacterial community composition in surface soil. Soil Biology & Biochemistry，2014，76：286—296

［30］Xu Y Q，Zhang S L，Hou H P，et al. Influence of CO2leakage from oil-producing wells on crop growth based on improved CASA model. International Journal of Remote Sensing，2016，37（2）：290—308

［31］Oppermann B I，Michaelis W，Blumenberg M，et al. Soil microbial community changes as a result of longterm exposure to a natural CO2vent. Geochimica et Cosmochimica Acta，2010，74（9）：2697—2716

［32］McFarland J W，Waldrop M P，Haw M. Extreme CO2disturbance and the resilience of soil microbial communities. Soil Biology & Biochemistry，2013，65：274—286

［33］Luis E S M，Paula A，Estanislao L C，et al. Highthroughput sequencing of 16S RNA genes of soil bacterial communities from a naturally occurring CO2gas vent. International Journal of Greenhouse Gas Control，2014，29：176—184

［34］周漢昌，劉文釗，劉毅，等. 土壤團(tuán)聚體N2O釋放與反硝化微生物豐度和組成的關(guān)系.土壤學(xué)報，2015，52 （5）：1144—1152

Zhou H C，Liu W Z，Liu Y，et al. Relationships of N2O emission with abundance and composition of denitrifying microorganisms in soil aggregates（In Chinese）. Acta Pedologica Sinica，2015，52（5）：1144—1152

［35］Shen C C，Xiong J B，Zhang H Y，et al. Soil pH drives the spatial distribution of bacterial communities along elevation on Changbai Mountain. Soil Biology & Biochemistry，2013，57（2）：204—211

［36］Liu Y，Zhou H M，Pan G X，et al. Short-term response of nitrifier communities and potential nitrification activity to elevated CO2and temperature interaction in a Chinese paddy field. Applied Soil Ecology，2015，96：88—98

［37］Singh K M，Shah T，Deshpange S，et al. High throughput 16S rRNA gene-based pyrosequencing analysis of the fecal microbiota of high FCR and low FCR broiler growers. Molecular Biology Reports，2012，39（12）：10595—10602

Effects of Underground Leakage of High Concentration CO2on Soil Microbial Community Structure

CHEN Fu1，2YANG Baodan1MA Jing2?ZHANG Shaoliang2ZHANG Yuanyuan3
（1 School of Environment Science and Spatial Informatics，China University of Mining and Technology，Xuzhou，Jiangsu 221116，China）
（2 Low Carbon Energy Institute，China University of Mining and Technology，Xuzhou，Jiangsu 221008，China）
（3 Sinopec Petroleum Engineering Design Company Limited，Dongying，Shandong 257026，China）

【Objective】 To cope with the trend of global warming，CO2capture and storage （CCS）is one of the major technologies for reduction of CO2emission. And then the captured CO2，nil in commercial value，is injected underground to raise the output of petroleum or coalbed methane，so as to maximize its economic profit. However，the prevailing technologies are far from being capable of guaranteeing zero leakage of the stored CO2during the processes of CO2storage and transportation，geological activities and human misoperation. Once the risk of potential CO2leakage becomes real，the leakage will pose an enormous threat to the near surface ecosystem. It is，therefore，essential to explore in depth effects of underground CO2leakage on farmland soil environment，especially tolerance and sensitivity of soil microbial communities to different concentrations of CO2. So the study was conducted. 【Method】 An experimental platform to simulate underground CO2leakage was constructed in an idle farmland，in the South Lake Campus of the China University of Mining and Technology. The platform was used to simulate leakage of CO2varying in intensity，i.e. 400 g m-2d-1，800 g m-2d-1，1 200 g m-2d-1and 2 000 g m-2d-1for 60 days，thus forming four treatments，i.e. L-400，M-800，H-1200 and E-2000. Besides，the experiment also had a control group and a recovery group. Soil samples were collected from the four treatments and their 3 replicates on the 14th，30th and 60th day after the start of simulated leakage，from the control group on the day before the start of leakage （C1）and at the end of the experiment（C2），and from the recovery group，which was actually Treatment E2000，60 days after the stop of the leakage. The soil samples were analyzed for soil physical and chemical properties，concentration of occluded soil，structure and α and β diversities of soil microbial community with the conventional physicochemical analysis method and the Illumna second generation gene sequencing method based on the Miseq platform.【Result】 Results show that CO2leakage decreased soil pH，electrical conductivity and nitrate nitrogen content and the effect was enhanced with rising CO2concentration whereas it had an opposite effect on soil organic matter content. In all the four treatments. Soil CO2concentration increased till it reached saturation in 24 h，when soil CO2gas concentration leveled off at 1.60%，4.80%，10.80% and 19.60%，respectively. Along with increasing CO2flux，soil microbial community decreased in diversity，Chao index and Shannon index，by 17.00%～27.80% and 6.10%～9.50%，respectively. In contrast，soil microbial community increased in β diversity（NMDS index）in Treatments L-400 and M-800（low CO2concentration）but decreased in Treatments H-1200 and E-2000. Some bacteria，like Bacteroidales，varied extremely，either rising up or falling down in relative abundance with increasing CO2leakage. 【Conclusion】 The structure，diversity and abundance of soil microbial community varied significantly from treatment to treatment. Soil pH and CO2flux were the two most important environmental factors affecting soil microbial diversity. Bacteroidales was very sensitive to CO2stress so that it can beused as a key indicator in monitoring and evaluating ecological risk of underground CO2leakage. The 60 d short term recovery experiment indicates that the soil microbial community recovered well in diversity and richness，but it is still unclear whether it did in function. Therefore，in future studies，focuses should be laid on impacts of underground high concentration CO2leakage on functions of the soil ecosystem.

Carbon capture and storage；CO2leakage；Microbial community；High-throughput sequencing；Diversity analysis

X171

A

10.11766/trxb201607110265

（責(zé)任編輯：盧 萍）

* 國家科技支撐計劃項目（2012BAC24B05）和江蘇省煤基CO2捕集與地質(zhì)儲存重點(diǎn)實驗室開放基金項目（2015A01、2015B02）資助 Supported by the National Key Technology R&D Program of China（No. 2012BAC24B05）and the Project of Key Laboratory of Coal-based CO2Capture and Geological Storage，Jiangsu Province（Nos. 2015A01，2015B02）

? 通訊作者 Corresponding author，E-mail：jingma2013@cumt.edu.cn

陳 ?。?974—），男，江蘇射陽人，博士，教授，主要從事生態(tài)安全監(jiān)測、評價與修復(fù)研究。E-mail：chenfu@cumt.edu.cn

1；

2016-09-26；優(yōu)先數(shù)字出版日期（www.cnki.net）：2016-10-20