李剛強,李海俠,張小凌,任二慧,李峰銳,劉大祥,陳宇東,張子豪
(昆明理工大學國土資源工程學院,昆明 650093)
湖泊沉積物中的磁性礦物是研究古環(huán)境、古氣候的重要對象[1-3]。但是受到沉積動力、物源輸入、沉積成巖等因素影響,在同一湖盆內(nèi)沉積物在空間和時間上存在相對差異[4]。與湖心穩(wěn)定的水淹環(huán)境相比,湖濱區(qū)域沉積動力多變、沉積環(huán)境復(fù)雜[5],在漫長的地質(zhì)沉積過程中能較好地記錄湖面波動的信息。因此對湖濱不同巖性沉積物進行詳細的磁學測量,是研究湖泊演化的重要途徑。盡管同一湖泊內(nèi)湖中心沉積物與湖濱沉積物在巖性、磁性礦物的含量存在差異,但其磁性參數(shù)所反映的氣候環(huán)境信息在時間上具有一致性[5-6],說明湖濱沉積物的磁學特征對氣候和環(huán)境具有更廣泛的指示意義。并且與其他類型的湖泊相比,高原類湖泊入湖河流大多短小,沉積物來源相對單一,沉積環(huán)境相對穩(wěn)定,其記錄的磁學信息更加豐富。
①~⑧為巖層序號;LFκ為低頻體積磁化率圖1 采樣位置及鉆孔柱狀圖Fig.1 Sampling position and borehole histogram
滇池是位于云貴高原腹地的典型高原大型湖泊,在近3.4 Ma連續(xù)的湖相沉積物中記錄下各類氣候與構(gòu)造事件[7]。許哲平等[7]根據(jù)植物孢粉優(yōu)勢種群的相對含量變化探討了該區(qū)古氣候的變化,雖然得出了第四紀以來昆明盆地氣候以溫暖為主的重要結(jié)論,但并未對沉積物樣品做系統(tǒng)的磁性參數(shù)測量,其沉積過程中環(huán)境變化的細節(jié)還有待進一步解讀。本研究以位于云貴高原腹地的滇池東部DC1鉆孔湖濱沉積物為研究對象,應(yīng)用環(huán)境磁學的方法,進行了系統(tǒng)的磁學參數(shù)研究,測試了質(zhì)量磁化率χ、非磁滯剩磁(anhysteretic reamanent magnetization,ARM)、飽和等溫剩磁(saturation isothermal remanence magnetization,SIRM)、低頻磁化率(χlf)、高頻磁化率(χhf)和軟剩磁(soft remanence,IRM-300 mT),以及樣品pH、總有機碳(total organic carbon,TOC)。其目的是探討磁性特征、氣候環(huán)境、湖泊演化以及總有機碳之間的相互關(guān)系,為研究區(qū)域氣候和湖泊的總體演化提供依據(jù)。
昆明盆地位于云貴高原腹地,主要受西南季風系統(tǒng)的影響[8],盆地長約70 km,寬15~23 km,面積1 070 km2。盆地內(nèi)年均溫度15.6 ℃,年均降雨 1 000 mm。盆地內(nèi)四季溫差較小,表現(xiàn)為明顯的干濕兩季。滇池位于盆地西南,湖面海拔1 886 m,面積330 km2,平均水深5 m,最大水深8 m,是云南省面積最大的湖泊,也是中國第六大淡水湖[9]。姜朝松等[10]研究表明,滇池3 000年前為圖1(a)所示的淺綠色區(qū)域,后期由于人類活動湖泊面積不斷減小,最終演化成現(xiàn)代藍色區(qū)域[圖1(a)]。DC1巖芯取自昆明市滇池東部(地理坐標為24°55′32″N,102°47′35″E),總體上可分為[圖1(b)]:①褐黃色填雜黏土層,含植物根系,厚0.3 m;②淺褐色沙質(zhì)黏土層。厚度 3.3 m;③黃褐色砂質(zhì)黏土層含礫石,厚2.72 m;④褐灰色砂質(zhì)黏土層,厚4.1 m;⑤褐灰色粉質(zhì)黏土層,夾薄層砂礫,礫石磨圓度高,厚5.4 m;⑥褐黃色粉質(zhì)黏土層,厚0.88 m;⑦蘭灰色砂礫層,含礫量約50%~80%,礫石磨圓度高,分選性好,夾薄層灰褐色粉質(zhì)黏土,厚3.7 m;⑧灰黑色粉沙層,厚1.1 m。
在去除表面30 cm的人為擾動層后,每間隔 30 cm 取樣,總共獲得69個樣品。再將樣品帶回實驗室經(jīng)自然風干、碾碎,通過0.4 mm孔徑篩后按照編號收集。再將土壤樣品依次按序號放入8 cm3塑料小盒內(nèi),壓實并密封。對樣品的磁學測量包括:低頻磁化率(κlf,976 Hz)、高頻磁化率(κhf,15 600 Hz)、ARM、IRM-300 mT和SIRM(1 T)。在此基礎(chǔ)之上,計算了質(zhì)量磁化率χ、頻率磁化率(χfd,%)、S-ratio(IRM-300 mT/SIRM)、非滯后剩磁磁化率(χARM)和硬剩磁(hard isothermal remanence magnetization,HIRM)等。所使用的儀器為MFK1-FA多頻磁化率儀、Molspin交變退磁儀、脈沖磁化儀和Minispin旋轉(zhuǎn)磁力儀。并根據(jù)Hrouda[11]給出的計算公式[式(1)],將頻率磁化率測量結(jié)果進行轉(zhuǎn)換。此外,樣品pH測量采用HANNA HI8424型酸度計測定,TOC采用油浴重鉻酸鉀容量法(簡稱油浴法)測定。
χfd=χmfdln10/(lnfmHf-lnfmLf)
(1)
式(1)中:fmLf和fmHf分別為使用儀器MFK1-FA測量低頻、高頻磁化率時所使用的頻率,本研究中fmLf=976 Hz、fmHf=15 600 Hz;χmfd為使用儀器MFK1-FA在低頻(976 Hz)、高頻(15 600 Hz)下測量的低頻磁化率和高頻磁化率的差值。
環(huán)境磁學參數(shù)綜合分析能反映環(huán)境中磁性礦物組成、含量和顆粒大小[12]。LFк、ARM、SIRM、χfd%通常對磁性礦物的濃度敏感[13]。LFκ對粗顆粒的多疇(MD,粒徑>2 μm)成分的變化靈敏,它的高值表示多疇粒子的含量較高。ARM值的大小對樣品中穩(wěn)定單疇(SSD,粒徑0.02~0.04 μm)粒子的濃度具有良好的指示作用[13]。天然針鐵礦為高矯頑力磁性礦物,當外加磁場大于6T時才會飽和,當外加磁場為1 T時,赤鐵礦對等溫剩磁(isothermal remanence magnetization,IRM)的貢獻比針鐵礦大100倍,因此,當外加磁場為1 T時SIRM可表示赤鐵礦的信息[14]。同時,SIRM也能反映磁性礦物的顆粒大小和類型的變化,在給定的濃度下,粗粒的MD顆粒比SSD顆粒有較低的SIRM值[15]。而χfd%主要用于指示樣品中超順磁顆粒(SP,粒徑<0.03 μm)的存在和相對含量[16]。當土壤樣品中小于2%時基本不含有SP顆粒,χfd%在2%~10%時表明土壤樣品中含有少量的SP顆粒且不占主導(dǎo)地位[16]。亞鐵磁性礦物是土壤質(zhì)量磁化率的主要貢獻者[17],因此質(zhì)量磁化率的大小可用作亞鐵磁性礦物的粗略量度[15]。
HIRM與S-300 mT(IRM-300 mT與SIRM的比值)能反映反鐵磁性礦物的含量。HIRM可估算反鐵磁性礦物(主要是赤鐵礦或/和針鐵礦)的絕對濃度[18]。S-300 mT值則可衡量樣品中亞鐵磁性礦物和不完整反鐵磁性礦物的相對含量。其數(shù)值接近1,表示磁性礦物以低矯頑力的磁鐵礦和磁赤鐵礦為主,相反,當S-300 mT數(shù)值越接近零或為負值時,說明硬磁性礦物(赤鐵礦+針鐵礦)的含量越高[19]。在墨西哥火山附近的湖相沉積物中該值小于0.6,可表示高矯頑力的赤鐵礦存在[3]。故當SIRM與HIRM取得峰值且S-300 mT<0.6時,表明赤鐵礦的含量較多。
另外,比值參數(shù)SIRM/χ可以指示磁性礦物顆粒的大小,其值越高則礦物顆粒越小。S-300 mT的值則代表了樣品中亞鐵磁性礦物和不完整反鐵磁性礦物的相對比例。因此,基于SIRM/χ與S-300 mT制作的SIRM/χ-S-300 mT散點圖既能反映磁性礦物顆粒的大小,又能反映亞鐵磁性礦物和不完整反鐵磁性礦物的相對含量[19]。該研究樣品的SIRM/χ-S-300 mT散點圖如圖2所示,大部分數(shù)據(jù)點都集中在右下角,表明DC1巖芯樣品的磁性礦物總體以粗粒的亞鐵磁性礦物為主,而細粒礦物與反鐵磁性礦物的含量較少。
圖2 巖芯樣品SIRM/χ-S-300 mT散點圖Fig.2 SIRM/χ-S-300 mT scatter plot of core samples
綜合磁性參數(shù)隨深度的變化特征(圖1、圖3),可分為6個特征段,大致與巖層對應(yīng),如圖3所示。
圖3 巖芯樣品磁化參數(shù)隨深度變化剖面圖Fig.3 Profile of magnetization parameters of core samples varying with depth
(1)第Ⅰ段(21.0~19.9 m):該區(qū)段上LFκ平均值達到802.0×10-5SI;ARM、χfd%、χ與SIRM平均值分別為116.3×10-5A·m2/kg、6.58%、475.7×10-8m3/kg和1 047.2×10-5A·m2/kg;S-300 mT平均值為0.62,表明磁性礦物以亞鐵磁性的粗粒多疇為主。
(2)第Ⅱ段(19.9~15.6 m):LFκ增大趨勢明顯,且平均值較高為701.0×10-5SI;ARM與χfd%的變化呈增大趨勢,平均值103.8×10-5A·m2/kg和3.75%;χ在該區(qū)段上呈現(xiàn)出峰值,且平均值為440.9×10-8m3/kg;與第Ⅰ段相比S-300 mT平均值達到0.82,更接近于1,表明磁性礦物以多疇顆粒的亞鐵磁性礦物為主。
(3)第Ⅲ段(15.6~14.7 m):LFκ與ARM減小趨勢明顯,平均值分別為85.1×10-5SI和11×10-5A·m2/kg;χfd%在該區(qū)段上取得一個小的峰值,但均值較小,為3.5%;SIRM與上區(qū)段相比也有所減小,平均值為1 132.0×10-5A·m2/kg;S-300 mT平均值為0.47,小于0.6,表明在該區(qū)段上亞鐵磁性礦物的總體數(shù)量有所減少,主要以反鐵磁性的多疇顆粒為主。
(4)第Ⅳ段(14.7~11.1 m):LFκ與ARM平均值分別為248.5×10-5SI和33.9×10-5A·m2/kg;χfd%平均值4.0%,與上區(qū)段相比數(shù)值有所增加。因此,該區(qū)段磁性礦物以較粗的多疇顆粒居多,另外含有一定量的超順次顆粒,但不占主導(dǎo)地位,SIRM平均值為1 491.6×10-5A·m2/kg,質(zhì)量磁化率曲線呈增大的趨勢;S-300 mT平均值為0.75,表明反鐵磁性礦物的含量相對減少,亞鐵磁性礦物的含量有所增加。
(5)第Ⅴ段(11.1~3.6 m):在區(qū)段LFκ與ARM明顯的增大,平均值分別為474.9×10-5SI和79.9×10-5A·m2/kg;SIRM有所變化但總體趨勢增大,平均值為1 286.7×10-5A·m2/kg;S-300 mT在該區(qū)段上平均為0.77,表明磁性礦物以亞鐵磁性礦物為主。
(6)第Ⅵ段(3.6~0.3 m):LFκ、ARM整體呈現(xiàn)峰值特征,平均值分別為237.7×10-5SI和53.5×10-5A·m2/kg;χfd%的變化趨勢比較平穩(wěn),平均值為3.45%;SIRM在該區(qū)域穩(wěn)定增加,平均值為1 519.9×10-5A·m2/kg,S-300 mT在該區(qū)段上平均為0.45,表明該區(qū)段反鐵鐵磁性礦物(以赤鐵礦為主)含量較多,磁性礦物顆粒較粗,以多疇顆粒為主。
樣品pH是反映沉積環(huán)境的指標。張澤鑫等[20]在利用熱合成法對針鐵礦的合成研究中發(fā)現(xiàn),當pH為8時,合成物主要是針鐵礦和赤鐵礦的混合體,這表明反鐵磁性礦物能在弱堿性環(huán)境中生長與富集。DC1沉積序列中pH變化幅度較小[圖4(a)],沉積環(huán)境從弱堿-弱酸-弱堿均勻變化。
圖4 pH、TOC與χ的剖面特征圖及滇池演化示意圖 Fig.4 Profile characteristics of pH,TOC and χ and evolution diagram of Dianchi Lake
沉物積的磁性參數(shù)曲線可以反映沉積物中磁性礦物在沉積過程中經(jīng)受的物理、化學及生物作用[21],在推斷沉積環(huán)境、分析古氣候的變化規(guī)律及其細節(jié)方面發(fā)揮著重要作用[1]。但值得注意的是運用沉積物的磁性參數(shù)來探討與氣候環(huán)境之間的關(guān)系時,應(yīng)考慮沉積物的來源、沉積類型以及總的沉積背景,因為這些因素將影響沉積與氣候環(huán)境之間的響應(yīng)機制[22]。在北方黃土沉積環(huán)境中,磁化率能反應(yīng)其氣候變化,低值磁化率表明黃土沉積中的鐵磁性礦物含量較低,反映沉積氣候干冷,植被覆蓋較差,區(qū)域成土作用弱;磁化率值相對較高,反映成土作用較好,成土環(huán)境較為暖濕[23]。張新榮等[24]在松嫩平原南緣現(xiàn)代沉積物磁化率特征及古氣候環(huán)境意義的研究中認為湖泊水深較大時沉積的泥質(zhì)物磁化率較高,對應(yīng)于濕潤氣候;反之,水位低時沉積的砂質(zhì)物磁化率較低,對應(yīng)于干旱氣候。除此之外,磁性礦物的類型與沉積環(huán)境之間也有著密切的聯(lián)系。如沉積物中的常見的反鐵磁性礦物(赤鐵礦和針鐵礦)也可以作為反映氣候的指標[25]。一般認為赤鐵礦形成于脫水反應(yīng),干旱環(huán)境有利于形成赤鐵礦;針鐵礦通常是從水溶液中直接沉淀形成,因此,潮濕環(huán)境利于其(針鐵礦)發(fā)育[25]。
在DC1沉積序列中,除第Ⅲ區(qū)段與第Ⅵ區(qū)段為反鐵磁性的赤鐵礦以外,其他區(qū)段均以亞鐵磁性的(磁鐵礦和磁赤鐵礦)礦物為主,這或許是整個沉積過程中(21 m的沉積過程)滇池湖泊環(huán)境發(fā)生過多次改變所導(dǎo)致。如圖4(d)所示,一開始采樣區(qū)域被滇池湖面淹沒,這一時期對應(yīng)的沉積序列為第Ⅰ、Ⅱ區(qū)段。之后,滇池湖面發(fā)生變化,采樣區(qū)域地勢較高的地方露出湖面,形成淺水沼澤,局部地區(qū)形成干旱環(huán)境,因此,在第Ⅲ區(qū)段磁性礦物表現(xiàn)為反鐵磁性的赤鐵礦行為。第Ⅳ、Ⅴ區(qū)段對應(yīng)的沉積環(huán)境為湖進時期,在這段時間滇池湖面升高,采樣區(qū)域再次被湖水淹沒。在這之后,由于人類活動的影響(如海口河的開挖、圍湖造田等)[10],滇池湖退,采樣區(qū)域再次露出水面,形成脫水環(huán)境。這也是表層(0.3~3.6 m)樣品磁性參數(shù)表現(xiàn)為赤鐵礦特征的原因。
沉積物中磁性礦物的顆粒大小受沉積環(huán)境(包括物理環(huán)境與化學環(huán)境)的影響較大[22],因此沉積物的粒度構(gòu)成對于了解研究區(qū)的沉積過程及沉積環(huán)境特征具有指示意義[26]。早期的一些研究中發(fā)現(xiàn),湖泊沉積物中磁性礦物的顆粒大小與氣候環(huán)境之間也有著良好的響應(yīng);當氣候暖濕時,沉積物中的超順磁組分含量較高;氣候干冷時,其所含的超順磁組分相對較低[27]。然而這一結(jié)論在湖相沉積中或許并不具有普遍性。因為近年來的研究發(fā)現(xiàn),在湖泊環(huán)境中磁性礦物會受到侵蝕作用、泥沙運輸、沉積作用和后期的成巖作用的影響[28-30];同時,水體和沉積物中的微生物所驅(qū)動一系列化學反應(yīng)(包括硝酸鹽還原、錳還原、鐵還原、硫酸鹽還原等)將導(dǎo)致磁性的溶解[22,29],并且,與粗粒的磁性礦物相比,顆粒較細的首先溶解[22]。因此,當沉積物中粗顆粒的磁性礦物占據(jù)主導(dǎo)地位時,不能簡單地理解為:沉積過程對應(yīng)的氣候為干冷時期。
Xu等[22]在鶴慶盆地湖相序列磁性礦物溶解與印度季風變化的關(guān)系研究中發(fā)現(xiàn),TOC含量的變化對磁性礦物的溶解具有指示作用,TOC的升高對應(yīng)著磁化率的減小,也對應(yīng)著磁性礦物的溶解過程。該研究中TOC與磁化率隨深度的變化關(guān)系如圖4(b)和圖4(c)所示,紅色區(qū)域TOC的高值恰好對應(yīng)磁化率的低值,這些特征標志著磁性礦物的溶解。除此之外,利用χ與χARM之間的關(guān)系制作的King圖能反應(yīng)磁性礦物的顆粒大小[31],如圖5所示。較為明顯的特征是樣品中80%以上的磁性顆粒分布在0.1~5 μm。χ高值樣品磁性顆粒粒徑更大,大部分分布在1~5 μm。因此滇池沉積序列樣品磁性礦物主要以粗粒為主,細粒組分較少,這一特征也是磁性礦物溶解的重要證據(jù)。
圖5 巖芯樣品King圖Fig.5 King diagram of core samples
在相關(guān)研究工作中,將沉積物中有機質(zhì)所含的碳定義為有機碳[32]。有機碳的含量可反應(yīng)沉積物中有機質(zhì)的輸入量和沉積環(huán)境對有機質(zhì)的保存能力[33]。馬龍等[33]對張北高原安固里淖湖沉積物有機碳的研究表明,TOC數(shù)值的大小可反應(yīng)氣候環(huán)境的變化,高值(均值>1.2%)與溫暖氣候相對應(yīng),而低值(均值<1.06%)與寒冷時期相對應(yīng)。TOC在各區(qū)段的平均值如圖6所示,與安固里淖湖相比,滇池樣品沉積物有機碳在各區(qū)段的平均值均大于1.2%,這表明滇池第四系表層21 m沉積序列對應(yīng)的氣候環(huán)境主要以溫暖氣候為主,并且各區(qū)段TOC均值差異明顯,這可能是由于在整個沉積過程中環(huán)境溫度出現(xiàn)了小幅度的波動。另外,TOC的高值集中在pH偏弱酸性的環(huán)境,低值分布在pH偏弱堿性的環(huán)境,這一現(xiàn)象表明與弱堿性環(huán)境相比弱酸性環(huán)境更有利于有機碳的保存。
圖6 巖芯樣品各段pH與TOC均值圖Fig.6 Mean chart of pH and TOC in each section of rock sample
(1)滇池東部第四系淺層湖濱相沉積物中磁性礦物主要以亞鐵磁性礦物為主,反鐵磁性礦物含量較少,主要集中在第Ⅲ、Ⅵ區(qū)段。磁性礦物的顆粒較粗,含有一定量的超順磁顆粒但不占主導(dǎo)地位。
(2)根據(jù)各區(qū)段磁性礦物類型,可反映出在整個21 m沉積過程中滇池經(jīng)歷了湖退—湖進—湖退的變化過程。從磁性參數(shù)和TOC所反映的環(huán)境特征來看,沉積過程對應(yīng)的氣候以溫暖為主,但略有波動。這與昆明盆地植物孢粉所反映的第四紀氣候變化特征一致。
(3)TOC的含量能反應(yīng)磁性礦物的溶解過程,在該研究中,細顆粒的磁性礦物較少可能與磁性礦物的溶解有關(guān)。TOC的高值集中在pH偏弱酸性的環(huán)境,低值分布在pH偏弱堿性的環(huán)境,這一現(xiàn)象表明與弱堿性環(huán)境相比弱酸性環(huán)境更有利于有機碳的保存。