李緒龍,張 霞,林春明,黃舒雅,李 鑫
內(nèi)生金屬礦床成礦機(jī)制研究國(guó)家重點(diǎn)實(shí)驗(yàn)室,南京大學(xué) 地球科學(xué)與工程學(xué)院,南京 210023
化學(xué)風(fēng)化作用是影響全球氣候變化、地表物質(zhì)循環(huán)和海水化學(xué)組成的重要因素,有利于調(diào)節(jié)全球碳循環(huán)(Walker et al., 1981; Berner et al., 1983)。河流連接陸地和海洋,從陸地向海洋搬運(yùn)大量物質(zhì),其所攜帶的風(fēng)化物質(zhì)具有重要的大陸風(fēng)化、侵蝕和環(huán)境等地質(zhì)信息(McLennan, 1993; Milliman and Farnsworth, 2011)。目前,許多研究利用大量的礦物學(xué)和地球化學(xué)指標(biāo)分析河流沉積物的化學(xué)風(fēng)化強(qiáng)度,進(jìn)而評(píng)價(jià)其流域的化學(xué)風(fēng)化特征,重建古氣候變化(Nesbitt et al., 1996; Ding et al., 2009; Garzanti et al., 2013),但因沉積物化學(xué)風(fēng)化的復(fù)雜性,常用地球化學(xué)指標(biāo)的使用都存在一定局限性,受多種因素控制,包括源巖類(lèi)型、氣候狀況、構(gòu)造和地形、植被、土壤發(fā)育、粒徑、水動(dòng)力過(guò)程和人類(lèi)活動(dòng)等(Gibbs, 1970; Meybeck, 1987; Grantham and Velbel,1988; Berner, 1992; Stallard, 1995; Gaillardet et al.,1999a; Oliva et al., 2003; Garzanti et al., 2008; Bouchez et al., 2011)。并且,化學(xué)風(fēng)化指標(biāo)所反映的是特定流域的綜合風(fēng)化情況及其繼承性,并不是沉積物在搬運(yùn)和沉積過(guò)程中的瞬時(shí)風(fēng)化情況(Gaillardet et al., 1999b; Li and Yang, 2010)。因此,在利用各種風(fēng)化指標(biāo)定量評(píng)價(jià)化學(xué)風(fēng)化強(qiáng)度和氣候變化時(shí),必須深入理解這些化學(xué)風(fēng)化指標(biāo)的適用范圍和控制因素。本文結(jié)合前人研究成果,分析了常用化學(xué)風(fēng)化指標(biāo)的基本原理、計(jì)算方法、不同化學(xué)風(fēng)化指標(biāo)所存在的局限性以及受到的各種影響因素?;诟骰瘜W(xué)風(fēng)化指標(biāo)原理和適用性,作者進(jìn)一步提出綜合運(yùn)用化學(xué)風(fēng)化指標(biāo)進(jìn)行化學(xué)風(fēng)化強(qiáng)度評(píng)價(jià)的流程,希望所提出的評(píng)價(jià)流程對(duì)實(shí)際工作具有借鑒意義。
其中應(yīng)用較多的化學(xué)風(fēng)化指標(biāo)有化學(xué)蝕變指數(shù)CIA(Chemical Index of Alteration)、帕克風(fēng)化指數(shù)WIP(Weathering Index of Parker)、斜長(zhǎng)石蝕變指數(shù)PIA(Plagioclase Index of Alteration)、化學(xué)風(fēng)化指數(shù)CIW(Chemical Index of Weathering)、成分變異指數(shù) ICV(Index of Compositional Variability) 、化學(xué)風(fēng)化指數(shù)CIX(Chemical Index of Weathering excluded CaO)和αAlE 指數(shù)等 (表1)。
Parker(1970)最先提出通過(guò)計(jì)算堿金屬和堿土金屬元素在風(fēng)化產(chǎn)物中所占的比例來(lái)評(píng)價(jià)硅酸鹽巖的風(fēng)化強(qiáng)度,將其定義為帕克風(fēng)化指數(shù)(WIP;表 1),后經(jīng) Hamdan 和 Burnham(1996)得以推廣,廣泛用于判別化學(xué)風(fēng)化強(qiáng)度。該指標(biāo)采用元素(K、Ca、Na、Mg)與氧結(jié)合的鍵強(qiáng)作為加權(quán)因子來(lái)反映巖石的風(fēng)化情況以及對(duì)進(jìn)一步風(fēng)化的敏感程度,因鍵強(qiáng)既可衡量破壞化學(xué)鍵所需能量,也能評(píng)估元素參與風(fēng)化反應(yīng)的相對(duì)可能性(Nicholls,1963)。因SiO2在風(fēng)化過(guò)程中移動(dòng)不規(guī)律,且損失比例小,所以指數(shù)未考慮SiO2。WIP值越小,風(fēng)化作用越強(qiáng)。WIP指數(shù)不適用于強(qiáng)風(fēng)化情況,因?yàn)閴A金屬和堿土金屬元素在強(qiáng)風(fēng)化情況下所剩無(wú)幾(Eswaran et al., 1973)。
表1 常用化學(xué)風(fēng)化指標(biāo)Table 1 Common chemical weathering indices
化學(xué)蝕變指數(shù)(CIA; 表1)用于定量評(píng)價(jià)長(zhǎng)石的化學(xué)風(fēng)化強(qiáng)度,為Al2O3與不穩(wěn)定氧化物的比例,反映長(zhǎng)石轉(zhuǎn)化為黏土礦物的風(fēng)化程度(Nesbitt and Young, 1982)。CIA值越高,指示硅酸鹽礦物中K、Na和Ca元素從母巖中淋失越多,長(zhǎng)石的化學(xué)風(fēng)化越強(qiáng)(Nesbitt and Young, 1984, 1989)。Fedo等(1995)總結(jié)得出: CIA=50~60,反映了弱風(fēng)化程度;CIA=60~80,為中等風(fēng)化程度;CIA=80~100,反映了強(qiáng)烈風(fēng)化程度。不同礦物和巖石具有不同的CIA值,平均上地殼的CIA值為50,長(zhǎng)石容易風(fēng)化,所以長(zhǎng)石的CIA值(50)一般較低,諸如更新世冰川黏土、黃土、平均頁(yè)巖的CIA值間于60~80,伊利石、蒙脫石、亞馬遜泥巖和殘留黏土等物質(zhì)的CIA值較高,而徹底風(fēng)化能達(dá)到綠泥石、高嶺石的CIA值(100)(表2; Nesbitt and Young, 1982)。
表2 部分礦物和巖石的CIA值 (Nesbitt and Young, 1982)Table 2 CIA values of some minerals and rocks(Nesbitt and Young, 1982)
在進(jìn)行CIA值計(jì)算的同時(shí)還應(yīng)進(jìn)行成分變異指數(shù)(ICV; Index of Compositional Variability; 表1)的檢查,其表示為風(fēng)化過(guò)程容易遷移的元素之和與不易遷移的鋁元素的比值。隨著風(fēng)化的不斷進(jìn)行,長(zhǎng)石遭受風(fēng)化逐漸轉(zhuǎn)化為黏土,ICV值也隨之逐漸降低,非黏土的硅酸鹽礦物具有比黏土更高的ICV值(表3),因此,ICV值可以表征礦物成分的成熟度。一般認(rèn)為,若ICV>1,表明細(xì)碎屑巖含很少黏土物質(zhì),成分成熟度較低;若ICV值<1,則表明細(xì)碎屑巖含較多黏土成分,指示強(qiáng)烈的化學(xué)風(fēng)化作用(Cox et al., 1995)。
表3 不同硅酸鹽礦物的ICV值 (Cox et al., 1995)Table 3 ICV values of different silicate minerals (Cox et al., 1995)
即使CIA值計(jì)算中除去碳酸鹽和磷酸鹽中的CaO值,但碳酸鹽的殘余會(huì)導(dǎo)致CIA值評(píng)價(jià)化學(xué)風(fēng)化強(qiáng)度存在誤差,所以在碳酸鹽顆粒存在的地方,為避免風(fēng)化強(qiáng)度計(jì)算誤差,通常使用化學(xué)風(fēng)化指數(shù)CIX (表1) 評(píng)估化學(xué)風(fēng)化情況,該指標(biāo)不用考慮CaO(Garzanti et al., 2014)。
除CIX的其他化學(xué)風(fēng)化指標(biāo)中均計(jì)算了硅酸鹽中的CaO*,對(duì)于CaO*的計(jì)算和校正,一般用公式 CaO*=n(CaO) -(10/3)×n(P2O5)-r×n(CO2),n 為摩爾含量,當(dāng)碳酸鹽礦物為方解石時(shí)r=1,為白云巖時(shí)r=0.5。因?qū)嶋H操作中不易區(qū)分獲得的CO2是來(lái)自于方解石還是白云石,通常采用McLennan (1993)提出的方法進(jìn)行校正,若n(CaO) < n(Na2O),則采用n(CaO)作為樣品CaO*,相反則采用n(Na2O)作為CaO*。
Gaillardet等(1999b)利用具有相似巖漿相容性的易遷移元素(Na、K、Mg、Ca、Sr、Ba)和不易遷移元素(Al、Ti、Sm、Nd、Th)的濃度比值估計(jì)風(fēng)化程度。這一系列風(fēng)化指標(biāo)被定義為αE,E為所研究元素,αE中的不易遷移元素一般為T(mén)h,Ti和 REEs(Rare Earth Elements),但由于 REEs、Th和Ti受到水動(dòng)力作用的強(qiáng)烈影響而優(yōu)先賦存在重礦物中,所以對(duì)于富含重礦物的樣品,利用αE評(píng)估風(fēng)化強(qiáng)度可能會(huì)產(chǎn)生很大誤差。為了降低水動(dòng)力造成的偏差,Garzanti等(2013)提出αAlE=(Al/E)sample/(Al/E)UCC,將Al作為不易遷移元素,因?yàn)锳l不在重礦物中富集。通常河流沉積物中的αAlE值變化可以反映出易遷移元素相對(duì)于大陸上地殼(the Upper Continental Crust; UCC)的遷移順序(Garzanti et al., 2014)。對(duì)于大多數(shù)沉積物而言,其遷移順序?yàn)?Na>Ca>Sr>Mg>K>Ba>Rb,說(shuō)明 Na和Ca的淋濾極多,Sr的淋濾強(qiáng),Mg和K的淋濾適度,Ba和Rb淋濾少量(Garzanti et al., 2013)。硅酸鹽風(fēng)化開(kāi)始階段為Na、Ca和少量K的淋濾,元素淋濾越強(qiáng),αAlE值越大。在硅酸鹽巖風(fēng)化為黏土過(guò)程中,同一風(fēng)化階段的αAlNa、αAlCa和αAlK與Al/Si呈正比(Garzanti et al., 2013)。
化學(xué)風(fēng)化是地球科學(xué)研究中的熱點(diǎn)話題,利用化學(xué)風(fēng)化指標(biāo)指示風(fēng)化情況在前人研究中不斷深化,但缺乏對(duì)控制化學(xué)風(fēng)化強(qiáng)度的因素和機(jī)制的深入探討。針對(duì)這一科學(xué)問(wèn)題,前人對(duì)化學(xué)風(fēng)化的控制因素已開(kāi)展大量研究,包括物源巖性、粒度、水動(dòng)力作用、氣候和植被等(Taylor and Mclennan, 1985; Oliva et al., 2003; Houben et al., 2012;Eynatten et al., 2012; Schneider et al., 2016; Garzanti and Resentini, 2016; Guo et al., 2018)。
較早研究中利用美國(guó)、波多黎各和冰島101條流域的水化學(xué)數(shù)據(jù),發(fā)現(xiàn)不同巖性因其化學(xué)成分差異有其獨(dú)特的風(fēng)化速率,從砂巖、花崗巖、玄武巖、頁(yè)巖至碳酸鹽巖,化學(xué)風(fēng)化速率逐漸增加 (Bluth et al., 1994; Sun et al., 2018)。但以上研究?jī)?nèi)容并未得出物源對(duì)化學(xué)風(fēng)化程度和化學(xué)風(fēng)化指標(biāo)的控制結(jié)果,而了解物源信息是解釋化學(xué)風(fēng)化變化和合理使用化學(xué)風(fēng)化指標(biāo)的前提(Gaillardet et al., 1999a;Chetelat et al., 2013; Garzanti and Resentini, 2016; Sun et al., 2018; 徐小濤和邵龍義, 2018)。
Sc和Th作為惰性元素,在化學(xué)風(fēng)化和沉積物形成過(guò)程中受到水動(dòng)力作用影響較小,能夠很好記錄沉積物的源巖性質(zhì)。通常,Sc/Th<1為花崗巖,反之為鎂鐵質(zhì)巖石(Taylor and Mclennan, 1985)。除指示物源信息外,還能反映粒徑對(duì)源巖成分信息的影響,如果同一流域的懸浮物和岸灘沉積物的Sc/Th近似相等,表明受粒徑控制作用較?。▓D1;Guo et al., 2018)??衫肧c/Th反映化學(xué)風(fēng)化指標(biāo)與火成巖源巖之間的關(guān)系,如果CIA與Sc/Th比值之間具明顯相關(guān)性,則表明源區(qū)風(fēng)化強(qiáng)度與鎂鐵質(zhì)或長(zhǎng)英質(zhì)物源有關(guān),反之則表明源區(qū)化學(xué)風(fēng)化強(qiáng)度受物源巖性影響較弱。前人大量研究證明,巖性對(duì)中國(guó)大陸流域硅酸鹽風(fēng)化的控制作用并不明顯(Li and Yang, 2010; Shao et al., 2012; Bi et al., 2015)。Chetelat等(2013)對(duì)長(zhǎng)江流域的懸浮物質(zhì)進(jìn)行了化學(xué)風(fēng)化研究,長(zhǎng)江上游為峨眉山大火成巖省有關(guān)的鐵鎂質(zhì)物源,中下游為長(zhǎng)英質(zhì)物源,通過(guò)化學(xué)蝕變指數(shù)(CIA)與不同火成巖分異指數(shù)(Cr/Th、Eu/Eu*)的相關(guān)性分析顯示長(zhǎng)英質(zhì)比鎂鐵質(zhì)的化學(xué)風(fēng)化程度高,研究表明氣候變化對(duì)沉積物源產(chǎn)生影響,進(jìn)而導(dǎo)致化學(xué)風(fēng)化強(qiáng)度的差異。邵菁清和楊守業(yè)(2012)通過(guò)研究長(zhǎng)江干支流懸浮物的CIA值在時(shí)空上的變化特征發(fā)現(xiàn)長(zhǎng)江下游(南京和南通)干流枯季懸浮物的CIA值多高于洪季(圖2),這主要是因?yàn)榱饔蚣竟?jié)性的降雨區(qū)遷移使得進(jìn)入下游干流的懸浮物來(lái)源,即上游和中下游的相對(duì)貢獻(xiàn)發(fā)生明顯變化。洪季5~6月,強(qiáng)降雨帶位于中下游地區(qū),中下游流域物理風(fēng)化加強(qiáng),對(duì)下游干流的懸浮細(xì)顆粒物質(zhì)貢獻(xiàn)增多,導(dǎo)致 5~6月入海懸浮物的CIA值較高;而6月中旬到7月上旬為長(zhǎng)江中下游地區(qū)的梅雨季節(jié),CIA值呈下降趨勢(shì),這可能與中下游河流物質(zhì)貢獻(xiàn)有關(guān);7~9月強(qiáng)降雨區(qū)向上游地區(qū)遷移,上游強(qiáng)烈的物理剝蝕輸運(yùn)大量碎屑物質(zhì)進(jìn)入下游干流,使得這期間入海懸浮物中來(lái)自上游的組分占主導(dǎo),CIA值較低;從10月開(kāi)始降雨帶又逐漸往中下游遷移,另外三峽水庫(kù)的枯季蓄水也攔截了大量上游侵蝕來(lái)沙,使得中下游流域的物質(zhì)輸入對(duì)河口區(qū)的懸浮物貢獻(xiàn)比例又逐漸增大,使得CIA值總體又逐漸升高(圖2)。以上研究?jī)?nèi)容證明CIA與研究區(qū)沉積物質(zhì)來(lái)源有關(guān),沉積物地球化學(xué)特征也會(huì)因此而改變。
圖1 長(zhǎng)江、浙閩和臺(tái)灣河流沉積物的CIA-Sc/Th 圖解Fig. 1 Plot of CIA vs. Sc/Th ratios of fluvial sediments along the East China Sea coast
圖2 長(zhǎng)江流域懸浮物CIA值季節(jié)性變化 (a) 南京樣品; (b) 南通樣品 (邵菁清和楊守業(yè), 2012)Fig. 2 Seasonal variation of CIA values of suspended solids in the Yangtze River basin (a) samples of Nanjing; (b) samples of Nantong
Nesbitt和Young(1982)最初提出CIA時(shí)以泥質(zhì)巖為研究對(duì)象,之后也被用來(lái)進(jìn)行砂巖常量元素的CIA計(jì)算并推測(cè)物源區(qū)風(fēng)化程度。碎屑沉積物在產(chǎn)生過(guò)程中元素組成受物理和化學(xué)過(guò)程控制,物理風(fēng)化過(guò)程根據(jù)沉積物顆粒大小、密度和形狀對(duì)其進(jìn)行水動(dòng)力分選,產(chǎn)生的粗細(xì)顆粒會(huì)對(duì)沉積物的礦物組成產(chǎn)生影響(粒度效應(yīng)),造成粗細(xì)顆粒中不同元素的富集,使得CIA等化學(xué)風(fēng)化指標(biāo)不能指示其真實(shí)的化學(xué)風(fēng)化強(qiáng)度,表明化學(xué)風(fēng)化指標(biāo)會(huì)受到粒徑影響(Eynatten et al., 2012; Guo et al.,2018)。如砂巖中保留的黏土礦物相對(duì)較少,CIA計(jì)算值偏低,因此,用砂巖計(jì)算出的CIA值不能錯(cuò)誤的采用泥質(zhì)巖的判斷標(biāo)準(zhǔn)來(lái)分析物源區(qū)風(fēng)化程度(Armstrong-Altrin et al., 2004; Garzanti et al., 2011)。
Shao等(2012)研究了長(zhǎng)江中下游河漫灘沉積物和懸浮物質(zhì),發(fā)現(xiàn)流域內(nèi)沉積物CIA與粒度存在較好相關(guān)性(R2=0.79),但不同樣品性質(zhì)和流域不同地區(qū)的粒度與CIA相關(guān)性不同,較粗粒沉積物(Ф<7.0)的平均粒徑與CIA相關(guān)性好(R2=0.71),而較細(xì)的懸浮物質(zhì)(Ф>7.0)平均粒徑與CIA相關(guān)性差(R2=0.23),各相關(guān)系數(shù)均在0.01水平(雙側(cè))上進(jìn)行相關(guān)性檢驗(yàn) (圖3)。長(zhǎng)江上游懸浮物粒度粗,CIA值總體較低,中下游懸浮物粒度相對(duì)較細(xì),CIA值高??梢?jiàn)粒度越細(xì),黏土組分相對(duì)富集,CIA值就越高。因此,CIA在實(shí)際應(yīng)用中需要考慮粒度影響(Shao et al., 2012)。長(zhǎng)江下游干流和入海懸浮物的CIA受懸浮物粒度的影響非常弱,且下游河道寬坡度小,水動(dòng)力弱,導(dǎo)致懸浮物的平均粒徑較細(xì)且變化不大,可反映流域綜合風(fēng)化歷史(Potter et al., 2001)。
圖3 長(zhǎng)江沉積物平均粒徑 (Mz) 與CIA 的相關(guān)性 (a: 長(zhǎng)江下游河漫灘與懸浮物CIA與粒度的相關(guān)性;虛線為總相關(guān)性;b: 長(zhǎng)江流域懸浮物CIA與粒度的相關(guān)性. 其中N代表樣品數(shù)量, Shao et al., 2012)Fig. 3 Correlations of the CIA values with mean grain size (Mz) from the Changjiang River sediments (a, Correlation between CIA and mean grain size of the lower Yangtze River floodplain and SPM samples; The dotted line is the total correlation; b, Correlation between CIA and mean grain size of SPM of the Yangtze River. The N value means the number of samples, Shao et al., 2012)
毋庸置疑,選擇不同粒徑沉積物指示同一物源區(qū)的風(fēng)化情況未必能取得風(fēng)化強(qiáng)度真實(shí)值。由于泥質(zhì)巖具有較好的均質(zhì)性和沉積后的低滲透性,受粒度和再風(fēng)化作用的影響較小,更好地保留了物源區(qū)信息,泥質(zhì)巖比其他碎屑巖更適合進(jìn)行物源區(qū)風(fēng)化程度及古氣候的研究(Garzanti et al., 2013)。Bao等(2018)考慮到化學(xué)風(fēng)化指標(biāo)的適用性,認(rèn)為選取沉積物中的黏土組分(<2 μm)進(jìn)行化學(xué)風(fēng)化指標(biāo)的計(jì)算,進(jìn)而可以更加準(zhǔn)確的反映源區(qū)的化學(xué)風(fēng)化強(qiáng)度。Xiong等(2010)發(fā)現(xiàn)同一位置的中粗粒沉積物(5~20 μm、20~63 μm)與細(xì)粒沉積物 (<5 μm)有不同的CIA值,表明風(fēng)化程度不一致,認(rèn)為選?。? μm細(xì)粒組分即可很好反映源區(qū)的風(fēng)化強(qiáng)度。因此,在利用CIA評(píng)價(jià)化學(xué)風(fēng)化時(shí),最好選用能反映物源成分的泥巖或細(xì)碎屑巖、具代表性的懸浮物,以此可降低粒徑對(duì)化學(xué)風(fēng)化指標(biāo)的控制作用(王自強(qiáng)等, 2006; Shao et al., 2012)。但如何精確界定最適合用于化學(xué)風(fēng)化分析的沉積物粒度范圍仍須進(jìn)一步研究。
沉積物再旋回效應(yīng)對(duì)化學(xué)風(fēng)化的影響不容忽視(Gaillardet et al., 1999a)。前人在研究青藏高原東部和俄羅斯遠(yuǎn)東河床沉積物的物源及風(fēng)化作用時(shí),發(fā)現(xiàn)沉積作用的不斷旋回以及對(duì)前一個(gè)沉積旋回物質(zhì)的繼承性會(huì)對(duì)化學(xué)風(fēng)化指標(biāo)的應(yīng)用產(chǎn)生較大影響(Borges et al., 2008)。再旋回過(guò)程中母巖物質(zhì)經(jīng)歷再次風(fēng)化,CIA值偏大,因此不能準(zhǔn)確評(píng)價(jià)源區(qū)風(fēng)化程度和古氣候特征,而ICV可用來(lái)判斷源區(qū)物質(zhì)是否遭受再旋回作用。一般ICV<1說(shuō)明樣品成熟度較高,富黏土礦物成分(高嶺石、蒙脫石、絹云母等),可能經(jīng)歷了再旋回作用或首次沉積時(shí)經(jīng)歷強(qiáng)風(fēng)化作用;ICV>1表明樣品成熟度較低,黏土成分不多(Barshad, 1966)。在利用CIA評(píng)價(jià)風(fēng)化時(shí)應(yīng)選取ICV>1的樣品排除再旋回影響。利用ICV確定是否遭受再旋回的基礎(chǔ)上,利用Th/Sc-Zr/Sc圖區(qū)分沉積物分選和再旋回作用的影響(圖4;Kasanzu et al., 2008),Th/Sc和Zr/Sc在沉積初期明顯正相關(guān),經(jīng)歷再旋回作用后,Zr/Sc大幅增加,Th/Sc無(wú)明顯變化。樣品分布呈趨勢(shì)線1表明沉積物并未受再旋回影響,成分接近于源巖;若樣品分布呈趨勢(shì)線2,則表明沉積物受到再旋回影響,使得Zr元素富集(圖4)。
圖4 Th/Sc vs. Zr/Sc圖解(藍(lán)色和紅色點(diǎn)分別代表坦桑尼亞?wèn)|北部新元古代Ikorongo群的頁(yè)巖和粉砂巖;Kasanzu et al., 2008)Fig. 4 Diagram of Th/Sc vs. Zr/Sc. Blue and read points represent the shales and siltstones of the Neoproterozoic Ikorongo group in northeastern Tanzania (Kasanzu et al., 2008)
石英作為物理性質(zhì)和化學(xué)性質(zhì)均十分穩(wěn)定的礦物,在再旋回過(guò)程中不斷匯聚,特別是在化學(xué)侵蝕過(guò)程中,通過(guò)多次沉積旋回而累積下來(lái),但石英的富集并不能簡(jiǎn)單認(rèn)為就是化學(xué)風(fēng)化強(qiáng)度的增加(Garzanti et al., 2019)。沉積物富含石英顆粒很難被界定是由于廣泛的風(fēng)化作用和長(zhǎng)石等不穩(wěn)定礦物的破碎,還是石英砂巖中的石英顆粒在再旋回過(guò)程中的不斷匯聚(Basu, 2017)。由于石英富集而導(dǎo)致的石英稀釋作用解釋了除Si以外的大多數(shù)元素的消耗,在探討再旋回作用對(duì)化學(xué)風(fēng)化指標(biāo)的控制作用中,也同樣要考慮到石英稀釋作用的影響,即化學(xué)風(fēng)化指標(biāo)隨SiO2相對(duì)含量增加的變化特征(Garzanti et al., 2013)。Guo等 (2018) 認(rèn)為CIA和αAlNa均不受石英稀釋作用的控制,而WIP隨石英稀釋作用的增強(qiáng)呈現(xiàn)減小趨勢(shì)(圖5)。
圖5 CIA、αAlNa 和 WIP受石英稀釋作用的影響情況 (改自Guo et al., 2018)Fig. 5 Diagrams illustrating the variations of CIA, αAlNa and WIP with increasing quartz dilution (modified from Guo et al., 2018)
CIA值常與A-CN-K [Al2O3-(CaO*+Na2O) -K2O]三角圖相結(jié)合來(lái)反映碎屑沉積物的化學(xué)風(fēng)化程度和母巖成分(Nesbitt and Young, 1989; Fedo et al., 1995;Nesbitt et al., 1996; Eynatten et al., 2003)。根據(jù)元素的遷移順序,風(fēng)化作用可分為以下三個(gè)階段:早期的去Na、Ca階段,主要為斜長(zhǎng)石的風(fēng)化,在A-CN-K圖上,風(fēng)化趨勢(shì)線平行于A-CN連線并指向A;中期的風(fēng)化為去K階段,風(fēng)化趨勢(shì)線平行于A-K連線,此時(shí)以鉀長(zhǎng)石和伊利石的風(fēng)化為主,斜長(zhǎng)石幾乎完全消失;晚期階段風(fēng)化殘余物為高嶺石、石英、三水鋁石以及鐵氧化物,落于A-CN-K圖的A點(diǎn)(圖6;陳駿等, 2001)。
圖6 Al2O3-CaO*+Na2O-K2O (A-CN-K) 三角圖 (改自Fedo et al., 1995; 徐小濤等, 2018)Fig. 6 Plotted as molar proportions on Al2O3-(CaO*+Na2O) -K2O (A-CN-K) diagram (modified from Fedo et al., 1995)
前人通過(guò)采集粉砂巖和砂質(zhì)粉砂巖進(jìn)行分析,發(fā)現(xiàn)細(xì)碎屑巖在成巖過(guò)程中的鉀交代會(huì)改變?cè)瓗r成分,需要進(jìn)行鉀交代作用校正以獲得原巖真實(shí)風(fēng)化程度。在使用CIA之前通過(guò)計(jì)算ICV進(jìn)行成熟度判別,選取ICV>1的樣品進(jìn)行CIA計(jì)算,再利用A-CN-K圖提出的公式對(duì)鉀交代進(jìn)行校正,經(jīng)校正后的CIA計(jì)算值可判斷源巖的風(fēng)化程度(徐小濤和邵龍義, 2018;林春明等, 2021)。除了可以利用A-CN-K圖校正鉀交代作用,還可利用 Panahi(2000)提出的公式進(jìn)行校正,計(jì)算得出未發(fā)生鉀交代的CIA值(CIAcorr),元素氧化物取摩爾含量,相關(guān)表達(dá)式如下:
Fedo等(1995)根據(jù)動(dòng)力學(xué)浸出速率繪制出古元古代的平均頁(yè)巖和晚太古代上地殼的風(fēng)化趨勢(shì)線 (圖6a中箭頭b所示; Nesbitt and Young, 1984),但實(shí)測(cè)的古元古代平均頁(yè)巖數(shù)據(jù)竟位于該趨勢(shì)線之下,表明已經(jīng)歷了鉀交代作用。鉀交代一般有兩種途徑:(1)鋁土礦向伊利石轉(zhuǎn)化(圖 6b);(2) 斜長(zhǎng)石向鉀長(zhǎng)石轉(zhuǎn)化,斜長(zhǎng)石被自生鉀長(zhǎng)石取代,即K對(duì)Na和Ca的取代(圖6c)。兩個(gè)過(guò)程都會(huì)導(dǎo)致K的富集,從而使得化學(xué)風(fēng)化趨勢(shì)偏移。
通過(guò)A-CN-K圖也可確定鉀的富集量和風(fēng)化前后對(duì)應(yīng)的CIA值(Fedo et al., 1995),圖6a中的鉀交代為第一種,即向鋁土礦中加入K2O,K的頂點(diǎn)與平均頁(yè)巖的連線與預(yù)測(cè)的風(fēng)化趨勢(shì)線相交于一點(diǎn),該點(diǎn)為沉積物發(fā)生鉀交代之前的成分點(diǎn),d線所指即為其對(duì)應(yīng)的CIA值(圖6a),而平均頁(yè)巖的CIA值為a線所指 (圖6a),表明頁(yè)巖在交代過(guò)程中平均增加了約10%的K2O。
在黏土形成過(guò)程中K2O被淋濾形成鉀離子,鉀離子進(jìn)入溶液后既可以形成鉀礦物也可以通過(guò)離子交換吸附在黏土表面,CIA值不能準(zhǔn)確反映化學(xué)風(fēng)化作用強(qiáng)度,而使用化學(xué)風(fēng)化指數(shù)CIW(Chemical Index of Weathering)則可避免該點(diǎn)。CIW在CIA的基礎(chǔ)上去除了鉀,若CIW變化趨勢(shì)與CIA一致,既可排除鉀交代作用對(duì)CIA指數(shù)演化趨勢(shì)的影響,又可證明研究數(shù)據(jù)的可靠性(表1; Harnois,1988; 齊靚等, 2015)。但CIW的計(jì)算可能忽略了鉀長(zhǎng)石中的鋁(Fedo et al., 1995),會(huì)導(dǎo)致計(jì)算結(jié)果偏高,而斜長(zhǎng)石蝕變指數(shù)PIA(Plagioclase Index of Alteration; 表1)將鉀長(zhǎng)石中的鋁去除,可以作為定量分析斜長(zhǎng)石化學(xué)風(fēng)化強(qiáng)度的指標(biāo)。
Young和Nesbitt(1999)在研究古元古代化學(xué)風(fēng)化時(shí)得出,泥質(zhì)巖在沉積區(qū)經(jīng)歷進(jìn)一步風(fēng)化會(huì)導(dǎo)致CIA值偏大,因此CIA并不能真實(shí)反映源區(qū)的風(fēng)化程度。沉積物暴露于地表經(jīng)歷風(fēng)化作用形成古土壤和松散的殘積物。古土壤層表明巖石在沉積區(qū)經(jīng)歷成土作用并遭受淋濾作用,這種淋濾作用會(huì)導(dǎo)致 Ca、Na、K 等元素進(jìn)一步流失,但并不是物源區(qū)風(fēng)化作用的結(jié)果,因此使用經(jīng)歷過(guò)淋濾作用的古土壤泥質(zhì)巖化學(xué)成分計(jì)算的CIA 值不能反映物源區(qū)風(fēng)化作用(徐小濤和邵龍義, 2018)。成土作用過(guò)程中生物活動(dòng)和植物根系分泌的有機(jī)酸會(huì)影響到古土壤的風(fēng)化,對(duì)物源區(qū)風(fēng)化作用結(jié)果產(chǎn)生干擾 (Egli et al., 2001; Wilson,2004; Skiba, 2007)。利用叢枝菌根真菌(AMF)死亡或菌根裂解時(shí)向土壤釋放的球囊霉素研究球囊霉素蛋白(glomalin-related soil protein, GRSP)與地球地表風(fēng)化作用的相關(guān)性,該蛋白與地表植被初級(jí)生產(chǎn)力呈正相關(guān)。GRSP含量的升高表明該區(qū)域植被生長(zhǎng)旺盛,植物根系生長(zhǎng)過(guò)程中對(duì)沉積物中 K、Na 等礦物元素的吸收和對(duì)巖石的物理破壞將導(dǎo)致 CIA 值升高。此外,物源區(qū)的巖漿巖與上覆古土壤層的 CIA 會(huì)有明顯不同,反映出同一風(fēng)化殼不同深度的古土壤所代表的風(fēng)化作用程度明顯不同,例如 Szymański 和 Szkaradek (2018)所研究的安山巖及其上覆的古土壤層。因此,在選取樣品時(shí)應(yīng)避開(kāi)沉積區(qū)風(fēng)化殼或根土巖,否則會(huì)導(dǎo)致物源區(qū)風(fēng)化程度的判斷不準(zhǔn)確(Yang et al.,2014; Awasthi, 2017; Babeesh et al., 2017)。
一般來(lái)說(shuō),在濕潤(rùn)地區(qū)堿金屬和堿土金屬等化學(xué)元素的流動(dòng)性最高(Schneider et al., 2016),在寒冷地區(qū)化學(xué)風(fēng)化強(qiáng)度會(huì)由于低溫和較低濕度而表現(xiàn)較弱(Nesbitt and Young, 1996; Potter et al., 2001;Hall et al., 2002)。在濕熱的氣候條件下,風(fēng)化作用導(dǎo)致源巖中化學(xué)不穩(wěn)定礦物中流動(dòng)元素的淋濾,從而富集黏土礦物(如高嶺石; Garzanti et al., 2019)。研究發(fā)現(xiàn)黏土礦物的分布與緯度有關(guān):高嶺石、鋁土礦和蒙脫石向南北兩極逐漸減少,主要分布在濕熱地區(qū);綠泥石和伊利石主要分布在中、高緯度的寒冷和濕潤(rùn)地區(qū)(Fagel, 2007)。Dinis等(2017)以非洲南部大西洋邊緣地區(qū)為研究對(duì)象探討氣候?qū)Τ练e物組成的影響,并驗(yàn)證不同化學(xué)風(fēng)化指標(biāo)指示氣候的可靠性。通過(guò)研究非洲南部西海岸河流砂和黏土的化學(xué)風(fēng)化情況,結(jié)果顯示化學(xué)元素的遷移明顯受到強(qiáng)烈緯度控制和氣候梯度影響。
Yang 等(2004)在討論長(zhǎng)江和黃河流域硅酸鹽礦物的化學(xué)風(fēng)化時(shí),認(rèn)為氣候是控制化學(xué)風(fēng)化的主要因素,源巖等其他因素影響不大。Li和Yang(2010)比較了全球主要河流的CIA組成,認(rèn)為相比全球流域而言,中國(guó)緯向分布的主要河流流域的硅酸鹽巖化學(xué)風(fēng)化對(duì)氣候的響應(yīng)比對(duì)源巖、流域面積和物理侵蝕更為敏感。從構(gòu)造角度講,長(zhǎng)江的青藏高原源區(qū)及上游地區(qū)受新構(gòu)造運(yùn)動(dòng)影響顯著,新生代抬升剝蝕快,新鮮巖石暴露多而有利于化學(xué)風(fēng)化進(jìn)行,但由于上游源區(qū)氣候條件惡劣且快速剝蝕使得土壤層并不發(fā)育,硅酸鹽巖風(fēng)化難以深入進(jìn)行,而中下游的揚(yáng)子克拉通相對(duì)穩(wěn)定,地形坡度小,適宜的氣候條件使得化學(xué)風(fēng)化更為完全(邵菁清和楊守業(yè), 2012)。在利用化學(xué)風(fēng)化指標(biāo)評(píng)價(jià)化學(xué)風(fēng)化強(qiáng)度時(shí),要同時(shí)考慮研究區(qū)氣候與構(gòu)造特征對(duì)物源區(qū)化學(xué)風(fēng)化的控制作用。
Dinis等(2017)利用Hijmans 等(2005)的氣候數(shù)據(jù)計(jì)算了各流域的降雨量,結(jié)果顯示年平均降雨量與大多數(shù)風(fēng)化指標(biāo)(CIA、WIP、CIX)的相關(guān)性較差,說(shuō)明這些參數(shù)不能真實(shí)反映當(dāng)?shù)貧夂驐l件。在很大程度上,由于再旋回作用和對(duì)過(guò)去地質(zhì)歷史的繼承,河流黏土質(zhì)和砂質(zhì)的CIA、CIX和WIP指標(biāo)只能部分用于推斷集水區(qū)的降雨。在排除沉積物以長(zhǎng)英質(zhì)為主的小流域后,河泥的αAlMg與降雨相關(guān)性最好,高降水區(qū)域(>600 mm)河流砂中的αAlNa與降雨相關(guān)性最好(圖7)。鐵鎂質(zhì)物源沉積物中的鎂含量通常比長(zhǎng)英質(zhì)物源沉積物高,但是這種由物源成分造成的差異在細(xì)粒沉積物中會(huì)很大程度上減弱(Eynatten et al., 2012, 2016)。黏土的αAlNa與降雨缺乏相關(guān)性,可能是主要受再旋回作用的影響,在旋回過(guò)程中沉積物中的Na含量不斷被損耗。關(guān)于黏土礦物組合與降水的相關(guān)性特征,結(jié)果顯示低降水(<600 mm)時(shí)蒙脫石與降水呈負(fù)相關(guān),說(shuō)明高緯度干旱區(qū)河流沉積物中蒙脫石與降水之間的關(guān)系更為明顯(Dinis et al., 2017)。高嶺石含量和蒙脫石/高嶺石與降雨量的相關(guān)性較差。高嶺石的形成與地形有關(guān),在風(fēng)化長(zhǎng)期進(jìn)行的平坦地區(qū)更容易形成和富集,因此高嶺石并不能很好的用于指示氣候特征。除了利用泥質(zhì)沉積物的αAlMg、砂質(zhì)沉積物的αAlNa和蒙脫石反映氣候特征以外,黏土和沙中的CIA、WIP與CIX均與降水相關(guān)性較差,利用這些參數(shù)可在排除氣候(主要是降水)控制的情況下更加準(zhǔn)確反映源區(qū)的風(fēng)化特征,且采用細(xì)粒沉積物作為化學(xué)風(fēng)化樣品可能減小降雨變化對(duì)化學(xué)風(fēng)化程度的影響(圖7)。
趙占侖等(2018)在研究位于甘肅省古浪縣土門(mén)鎮(zhèn)的一處黃土—古土壤—湖相沉積—沙丘砂沉積旋回的地球化學(xué)特征時(shí),通過(guò)計(jì)算該剖面中位于末次冰消期的各沉積相的CIA值,發(fā)現(xiàn)其化學(xué)風(fēng)化強(qiáng)度與該時(shí)期內(nèi)的氣候變化特征相矛盾,認(rèn)為沉積物質(zhì)來(lái)源的不同,導(dǎo)致不同沉積相在堆積初始階段主量元素便存在差異,不同物源的不同沉積相計(jì)算的CIA指數(shù)失去了對(duì)比的前提。同一剖面的不同沉積相在堆積過(guò)程中因作用條件存在差異,可能會(huì)受到不同干擾因素的影響(如物源差異、粒度效應(yīng)、碳酸鹽沉積、水動(dòng)力等),是否化學(xué)蝕變指數(shù)(CIA)能夠指示同一剖面出現(xiàn)的不同沉積相的相對(duì)化學(xué)風(fēng)化強(qiáng)度,進(jìn)而揭示該地區(qū)古沉積環(huán)境及氣候變化?該問(wèn)題值得研究與探討。
沉積物在化學(xué)風(fēng)化過(guò)程中經(jīng)歷了復(fù)雜地質(zhì)過(guò)程,如何精確判別物源區(qū)化學(xué)風(fēng)化強(qiáng)度依舊是值得深入研究的科學(xué)問(wèn)題。利用化學(xué)風(fēng)化指標(biāo)評(píng)價(jià)源巖化學(xué)風(fēng)化強(qiáng)度要考慮各指標(biāo)的受控因素,否則其評(píng)價(jià)結(jié)果必然會(huì)失真。本文以化學(xué)蝕變指數(shù)CIA為主要評(píng)價(jià)指標(biāo),輔以WIP、CIW、PIA、ICV、CIX和αAlE等指標(biāo)共同指示化學(xué)風(fēng)化強(qiáng)度,理清化學(xué)蝕變指數(shù)CIA的控制因素及其影響因子,進(jìn)而達(dá)到定量化解釋研究區(qū)化學(xué)風(fēng)化強(qiáng)度的目的(圖8)。
圖8 物源區(qū)化學(xué)風(fēng)化強(qiáng)度判別流程圖Fig. 8 Diagram showing how to judge chemical weathering intensity in provenance area
掌握研究區(qū)相關(guān)地質(zhì)資料是研究物源區(qū)化學(xué)風(fēng)化的首要前提,需了解源區(qū)巖石、土壤特征和沉積物類(lèi)型等地質(zhì)信息。選擇細(xì)粒沉積物或懸浮物質(zhì)且避開(kāi)根土巖和古土壤等再風(fēng)化產(chǎn)物,可有效避免粒度效應(yīng)對(duì)化學(xué)風(fēng)化指標(biāo)的控制作用(Nesbitt and Young, 1982, 1984; 邵菁清和楊守業(yè), 2012; Yang et al., 2014; Guo et al., 2018)。主微量元素測(cè)試前要對(duì)樣品進(jìn)行前處理和去除其中的非硅酸鹽成分,樣品稱(chēng)量后放入特氟龍燒杯中,加入一定濃度的稀鹽酸和H2O2溶液去除樣品中的碳酸鹽和有機(jī)質(zhì)組分,然后做超聲震蕩,使得在分析之前沉積物完全解聚。一般采用自然沉降法或離心處理獲取特定粒徑的沉積物組分,向樣品中加入1 mol/L的至純鹽酸在60℃恒溫震蕩箱中進(jìn)行淋濾,從而去除自生的和吸附的磷酸鹽/碳酸鹽成分(Bao et al., 2018;Guo et al., 2018)?;瘜W(xué)風(fēng)化指標(biāo)計(jì)算過(guò)程中常利用McLennan(1993)提出的方法對(duì)硅酸鹽中的CaO進(jìn)行校正(Garzanti et al., 2014; Dinis et al., 2017)。通過(guò)獲取處理后樣品的地球化學(xué)信息進(jìn)行化學(xué)風(fēng)化控制因素分析,利用Sc/Th-CIA 判別圖反映原巖信息,Sc/Th比值與CIA之間相關(guān)性強(qiáng)弱指示化學(xué)風(fēng)化是否受源區(qū)原巖性質(zhì)的影響。利用Th/Sc-Zr/Sc判別圖進(jìn)一步確認(rèn)沉積分異或沉積再旋回的控制作用,若樣品點(diǎn)呈圖4中的趨勢(shì)線1則表明該樣品經(jīng)歷沉積分異作用,趨勢(shì)線2表明受沉積再旋回作用控制,通過(guò)選取ICV>1的樣品排除再旋回影響。鉀交代作用導(dǎo)致CIA值不能準(zhǔn)確反映硅酸鹽化學(xué)風(fēng)化強(qiáng)度,用A-CN-K圖或Panahi(2000)提出的公式進(jìn)行鉀交代作用的校正,最終通過(guò)CIA計(jì)算得出源巖化學(xué)風(fēng)化強(qiáng)度。
為驗(yàn)證和確保CIA受控因素以及真實(shí)反映源區(qū)化學(xué)風(fēng)化,有學(xué)者提出研究源區(qū)化學(xué)風(fēng)化的新思路。為剝離諸多因子的復(fù)雜影響,有研究者采集青藏高原東南緣的雅礱江下游不同地貌部位和植物群落的表土樣品做CIA的影響因子分析。發(fā)現(xiàn)其CIA與母巖、植被、海拔、坡度、坡向、土壤含水量、土壤有機(jī)碳(SOC)、全氮(TN)含量等環(huán)境因子呈顯著相關(guān)性(Feng et al., 2020)。各環(huán)境因子之間也具有一定共線性,運(yùn)用嶺回歸法(Zhang et al., 2013)可以解決自變量共線問(wèn)題并定量估算各環(huán)境因子對(duì)CIA值的控制作用。研究結(jié)果表明母巖性質(zhì)是控制化學(xué)風(fēng)化過(guò)程的主要因素, 地貌要素(海拔和坡度)是驅(qū)動(dòng)化學(xué)風(fēng)化過(guò)程的主要外部因素且海拔高度的影響大于坡度; 生物作用對(duì) CIA值有一定的貢獻(xiàn)。王亞暐等(2018)使用SPSS.22軟件將風(fēng)化前后巖樣的主量元素組成進(jìn)行單因素方差分析(ANOVA),得出其風(fēng)化過(guò)程中的控制性元素后進(jìn)行判別方程式擬合,與前人已有的風(fēng)化指標(biāo)進(jìn)行對(duì)比,得出適用于研究區(qū)千枚巖類(lèi)風(fēng)化的判別式 Fi=42.161×ln(SiO2)+0.062×ln(TiO2)-4.100×ln(A l2O3)+1.033×ln(Fe2O3)+0.005×ln(MgO)-0.230×ln(K2O)-1.179×ln(Na2O)-87.97。在進(jìn)行研究區(qū)的化學(xué)風(fēng)化判別式擬合之前,對(duì)源區(qū)物質(zhì)和風(fēng)化后的沉積物樣品的主量元素組成進(jìn)行單因素方差分析,進(jìn)而從中篩選出Sig.<0.05的元素組分,顯著性小于0.05意味著這些組分是控制化學(xué)風(fēng)化的關(guān)鍵組分,并且可以得知各元素組分影響化學(xué)風(fēng)化的比重。挑選出主量元素中的關(guān)鍵組分用于化學(xué)風(fēng)化指數(shù)擬合,使用Fisher判別法或主成分分析,通過(guò)數(shù)據(jù)降維處理,依據(jù)上述挑選出的關(guān)鍵組分的貢獻(xiàn)率確定特征化學(xué)風(fēng)化判別函數(shù)。為避免使用CIA評(píng)價(jià)化學(xué)風(fēng)化的單一性和誤差,結(jié)合WIP、CIW、PIA、ICV、CIX和αAlE等指標(biāo)共同指示源區(qū)元素遷移和淋濾特征,硅酸鹽化學(xué)風(fēng)化強(qiáng)度。通過(guò)以上研究思路和方法,定量解釋源區(qū)化學(xué)風(fēng)化強(qiáng)度。
各化學(xué)風(fēng)化指標(biāo)適用條件有所差異,受物源巖性、沉積分異、沉積再旋回、鉀交代作用沉積區(qū)進(jìn)一步風(fēng)化作用和環(huán)境因素差異等不同因素控制,厘清化學(xué)風(fēng)化指標(biāo)的受控因素是正確評(píng)價(jià)化學(xué)風(fēng)化強(qiáng)度的必要前提。再旋回影響可通過(guò)選擇細(xì)粒沉積物或懸浮物質(zhì)有效減弱。樣品地球化學(xué)分析前處理包括去除雜質(zhì)和硅酸鹽中CaO的校正。為明確化學(xué)風(fēng)化指標(biāo)受控因素特征,分別利用Sc/Th-CIA判別源區(qū)原巖特征,Th/Sc-Zr/Sc進(jìn)一步區(qū)分沉積分異和沉積再旋回影響,選取不受再旋回影響的樣品判斷并校正鉀交代影響。研究區(qū)的差異決定了化學(xué)風(fēng)化指標(biāo)所受控制因素不同,CIA受控因子分析可定量評(píng)價(jià)化學(xué)風(fēng)化指標(biāo)受控因素的影響因子,并以此得出研究區(qū)特征風(fēng)化指標(biāo)。在利用化學(xué)風(fēng)化指標(biāo)評(píng)價(jià)源區(qū)化學(xué)風(fēng)化情況時(shí),為避免使用單一化學(xué)風(fēng)化指標(biāo)導(dǎo)致的評(píng)價(jià)結(jié)果失真,往往結(jié)合WIP、CIA、PIA、ICV、CIX、αAlE等多個(gè)指標(biāo)共同指示源區(qū)化學(xué)風(fēng)化情況。
致謝:南京大學(xué)地球科學(xué)與工程學(xué)院凌洪飛教授、地理與海洋科學(xué)學(xué)院李徐生副教授給予悉心指導(dǎo)與熱情幫助;同濟(jì)大學(xué)海洋與地球科學(xué)學(xué)院郭玉龍老師、中國(guó)礦業(yè)大學(xué)(北京)地球科學(xué)與測(cè)繪工程學(xué)院徐小濤老師在成文過(guò)程中給予有益的建議和討論;兩位匿名審稿專(zhuān)家提出了建設(shè)性修改意見(jiàn),在此一并向他們致以衷心的感謝!