鈣同位素在高溫地質過程中的應用及前景

羅澤彬, 凌明星, 安亞軍, 劉 芳, 張兆峰

鈣同位素在高溫地質過程中的應用及前景

羅澤彬1,3, 凌明星2*, 安亞軍1, 劉 芳1, 張兆峰1

(1. 中國科學院 廣州地球化學研究所 同位素地球化學國家重點實驗室, 廣東 廣州 510640; 2. 東華理工大學 核資源與環境國家重點實驗室, 江西 南昌 330013; 3. 中國科學院大學, 北京 100049)

制約高溫過程中Ca同位素研究的因素主要有分析技術、儲庫端元和分餾機理。隨著質譜的發展和雙稀釋劑的應用, 分析技術已經能很好地滿足高溫過程同位素示蹤的要求, 但是目前對于儲庫端元的研究還比較薄弱, 而分餾機理也存在較大爭議。因此, 完善不同儲庫端元的Ca同位素組成, 厘清Ca同位素分餾機理勢在必行。全硅酸巖地球作為最重要的Ca儲庫, 前人研究表明, 其Ca同位素組成在0.94‰～1.05‰之間變化。Ca—O化學鍵能的強弱是造成礦物間Ca同位素組成差異的主要因素, 此外, 溫度和礦物組成差異也會引發Ca同位素分餾。熱擴散和化學擴散都會引起Ca同位素分餾, 但熱擴散只發生在特殊的環境中?；瘜W擴散受控于CaO化學勢梯度, 而化學勢梯度受控于體系的成分、溫度和壓力。部分熔融和熔體提取過程中輕同位素趨向于富集在熔體相中。由于不同含鈣礦物的Ca同位素組成存在差異, 因此, 在分離結晶過程中, 礦物晶出順序不同會引起Ca同位素分餾。變質作用過程中流體交代作用或碳酸鹽巖沉淀都會造成Ca同位素分餾。K-Ca衰變體系使得Ca同位素既可以對古老的富鉀巖石或礦物進行定年, 也可以利用放射性成因Ca對源區和殼幔物質循環進行示蹤。

Ca同位素; 高溫地質作用; 巖漿演化; 變質作用; 地質應用

0 引 言

Ca在元素周期表中位于第4周期第2主族, 原子序數為20, 屬于堿土金屬。廣泛分布于自然界各個角落, 既是生命體必備的營養元素, 也是組成巖石的8個主元素, 在上地殼中含量僅次于O、Si、Al和Fe, 排在第5位。此外, Ca屬于高度難熔元素, 具有極高的半凝聚溫度(1517 K), 是研究宇宙初始狀態和組成的重要元素窗口[1–2]。Ca有6個穩定同位素,40Ca (96.941%)、42Ca (0.647%)、43Ca (0.135%)、44Ca (2.086%)、46Ca (0.004%)和48Ca (0.187%)。其中40Ca占了主要的組成成分,48Ca本身屬于放射性同位素, 但由于其半衰期太長, 常被視為穩定同位素。在早期研究中, 由于技術和實驗條件的限制, 并沒有觀察到可分辨的同位素分餾作用[3]。直到1978年Russell.[4]采用雙稀釋劑法才較準確地測定了Ca同位素組成。進入21世紀以來, 隨著質譜技術的提高, Ca同位素研究取得了飛速的發展, Ca同位素相關報道也逐年增加。但主要集中在生物和低溫地質過程[5–12], 相對而言, 高溫地質過程中的Ca同位素研究報道卻相對較少。已有研究表明, Ca同位素在自然界樣品中44/40Ca((44Ca/40Ca)樣品/(44Ca/40Ca)標樣?1)變化可達約6‰[13], 其中硅酸巖中變化約2‰ (圖1)[7],近期研究發現更低的Ca同位素組成地幔橄欖巖[26], 說明高溫下Ca同位素的分餾可能更大。因此, Ca同位素在高溫地質過程存在可分辨的同位素分餾, 但就目前而言, 高溫過程中的分餾機理仍然不清楚。為了厘清Ca同位素在高溫過程中的分餾機理, 有效地利用Ca同位素示蹤高溫地質過程, 更多更詳細的工作還有待開展。本次研究擬就近年來已發表的相關Ca同位素文章, 總結Ca同位素在高溫地質過程中的同位素分餾作用并簡要展望未來發展前景。

圖1 各類火成巖的Ca同位素組成。數據來自文獻[1,2,5,9,10,14–25]

BSE為全硅酸巖地球。

BSE represents bulk silicate earth.

1 Ca的全硅酸巖地球儲庫和高溫過程中Ca同位素分餾

根據分餾的性質和原因, 同位素分餾可分為熱力學平衡分餾和動力學不平衡分餾。低溫過程中, 特別是生物參與過程, 由于元素的分布規律和遷移能力, 其同位素分餾主要是以動力學分餾為主。例如, 海水與礦物之間的Ca同位素分餾主要受控于結晶速率、溫度[27]或者含鈣礦物溶液的飽和狀態[28]。在植物和周圍無機環境中, 由于生物動力學分餾使得對Ca同位素吸收運移能力不同, 常使得植物更富集輕同位素[7]。與低溫過程不同的是, 高溫過程中的同位素分餾更小, 且既存在熱力學分餾也有動力學分餾參與。加上高溫過程中常以難熔且Ca含量低的硅酸巖等樣品為主, 導致分餾更難以被檢測。

制約高溫地質過程同位素研究的影響因素主要有3個方面: 分析技術、儲庫端元和分餾機理。就目前而言, 質譜的發展和雙稀釋劑法的應用已經可以把Ca同位素測試精度提高到約0.1‰[29–35], 能滿足利用Ca同位素示蹤大多數高溫地質過程的要求。

1.1 全硅酸巖地球的Ca同位素組成

在核幔分異過程中Ca幾乎不進入地核, 所以全硅酸巖地球(bulk silicate earth, BSE)是地球上最初的鈣源[14–15]。為了能進一步探討高溫過程中的Ca同位素地球化學行為, 準確限定BSE的Ca同位素組成是基本前提。早期BSE的Ca同位素限定主要是調查大洋玄武巖的Ca同位素組成[16], 但這是建立在部分熔融和巖漿演化過程中Ca同位素不發生分餾的前提下。已有研究表明, Ca同位素在部分熔融過程中會發生分餾[15,17], 而且, 從已有洋島玄武巖(oceanic island basalt, OIB)數據可以看出, OIB自身存在較大的變化范圍, 暗示地幔中可能存在Ca同位素不均一性, 或者因受到后期改造而不能代表源區Ca同位素組成。因此, BSE的Ca同位素限定不能通過大洋玄武巖直接獲得。

另一種限定BSE的方法是直接測定地幔橄欖巖的Ca同位素組成。上地幔橄欖巖中主要有4種礦物, 包括橄欖石、富鋁礦物相(斜長石、尖晶石和石榴石)、斜方輝石和單斜輝石, 而橄欖石和富鋁礦物相Ca含量極低, 因此, 斜方輝石和單斜輝石Ca同位素組成能大致代表上地幔的Ca同位素組成。在這樣的條件下, Huang.[14]對來自美國亞利桑拉州和新墨西哥州的地幔橄欖巖包體中的單斜輝石和斜方輝石進行了研究, 發現斜方輝石比單斜輝石更富重Ca同位素組成, 并認為Ca—O化學鍵長短的不同是控制二者之間同位素分餾的主要因素, Ca—O鍵越短越富集重同位素組成。Huang.[14]根據上地幔橄欖巖中兩種輝石的CaO含量、44/40Ca值以及兩種輝石在上地幔中所占的比例定量模擬計算了上地幔的Ca同位素組成44/40Ca=(1.05±0.04)‰。

值得注意的是, 這種方法存在一些不足之處。首先, 橄欖石雖然Ca含量很低, 但是如果其Ca同位素組成與全巖或其他礦物組分存在明顯的差異, 那么在限定BSE的時候橄欖石應該是一個不可忽略的成分。目前只存在極少橄欖石Ca同位素數據的報道[18], 很難作出判斷; 其次, Huang.[14]也提出不能否認斜方輝石和單斜輝石中存在動力學分餾, 后續研究也發現礦物的組分會影響Ca同位素分餾[18,36];再有源區的碳酸鹽巖交代會降低44/40Ca值[19], 并不能對交代做出定量評價; 另外, 其研究區只限于北美地區, 考慮到地幔Ca同位素的不均一性特征[13], 因此, 1.05‰可能更能代表上地?；蛘呔植康蒯５腃a同位素組成。不同的是, Kang.[15報道了28個來自蒙古和西伯利亞地幔橄欖巖包體的Ca同位素數據, 并挑選了其中14個飽滿的、未受交代作用影響的二輝橄欖巖對BSE進行了估計,44/40Ca位于0.90‰～0.99‰之間, 平均值為(0.94±0.05)‰, 并以此作為BSE的估計值。此外, 在經過熔體提取之后的難熔虧損地幔橄欖巖中觀察到44/40Ca與CaO呈負相關關系(圖2)[15,20,21], 指示熔融提取過程可能會導致殘余難熔組分變重, 因此, 1.05‰更能代表經過熔體提取之后的難熔虧損地幔, 而0.94‰更能代表BSE值。造成二者之間微小差異的原因可能是上地幔中含有經過熔體提取之后剩余的難熔虧損地幔。因為在橄欖巖熔體提取過程中單斜輝石優先被消耗, 導致熔體部分富集輕同位素組成而殘余難熔地幔富集重同位素組成[15。有一點需要指出的是, 就現在Ca同位素分析精度而言, 二者的差異也有可能是實驗室之間的系統誤差造成。但是同樣的問題是地幔的Ca同位素組成存在不均一性, 對于BSE的Ca同位素組成, 還需要更多的數據和研究去證明。

圖2 δ44/40Ca vs. CaO圖解, 數據來自文獻[1,2,5,9,10,14–25]

受碳酸鹽巖交代橄欖巖會富集輕同位素, 經過熔融提取過程的橄欖巖會富集重同位素; 長英質巖石中44/40Ca值變化大, 可能受控于不同含Ca礦物的分離結晶。

carbonate metasomatic peridotite will enrich light isotopes, and peridotite after melt extraction will enrich heavy isotopes;44/40Ca values are variable in feldspar rocks, and may be controlled by the crystallization of different Ca bearing minerals.

1.2 高溫過程中Ca同位素平衡分餾

Ca同位素在礦物之間存在明顯的同位素分餾, 例如, 在共生的斜方輝石(opx)和單斜輝石(cpx)中Δ44/40Caopx-cpx從?0.49‰變化至1.11‰[14–15]。高溫過程目前主要有以下幾種觀點: (1) Ca—O化學鍵的強弱。Ca在斜方輝石中配位數為6而單斜輝石配位數為8, 斜方輝石比單斜輝石具有更短的Ca—O化學鍵, Ca—O化學鍵越短,44/40Ca值越大[14]; (2) 高溫礦物間的熱力學平衡分餾。溫度越高, 同位素分餾越小, 這也是非傳統金屬穩定同位素在高溫下分餾小的原因。在不同溫度下, 不同礦物間發生同位素交換反應, 達到同位素平衡。當溫度較高時,44/40Ca同溫度的平方的倒數(1/2)成正比, Feng.[36]通過第一性原理計算了單斜輝石和斜方輝石在300～1000 ℃范圍內可分辨的Ca同位素分餾; (3) 礦物自身組分控制。已有研究發現, Δ44/40Caopx-cpx隨Ca/Mg比值降低而增大(圖3)[18,36], 并解釋為斜方輝石的組分是控制礦物間Ca同位素變化的主要因素。當溫度從1196 K到1267 K時, cpx(單斜輝石)減少, Δ44/40Caopx-cpx降低, 但這種礦物間同位素分餾的大小(1.12‰)與礦物同位素分餾系數、1/2之間的相關性不一致。相反, 它可以反映出溫度控制對斜方輝石(opx)晶體化學(即Ca含量)的影響, 進而影響Δ44/40Caopx-cpx。同時全巖的Lu、Yb、Fe等組分和全巖的44/40Ca值也表現出很好地相關性, 說明全巖組分也會影響礦物間Ca同位素分餾[18]。Wang.[38–39]和Huang.[40]等用第一性原理模擬計算了Ca同位素分餾和CaO濃度之間的變化特征, 結果表明, 在一定濃度范圍內, Ca同位素分餾表現出明顯的濃度效應, 即礦物中CaO濃度會影響Ca同位素分餾。

圖3 δ44/40Caopx-cpxvs. Ca/Mg在斜方輝石中的圖解, 數據來自文獻[14,15,19,37]

1.3 高溫過程中Ca同位素動力學分餾

前面更多討論的是高溫過程中的熱力學平衡分餾, 那么在高溫過程中動力學分餾Ca同位素又是怎樣體現的呢？實驗研究表明, 輕的同位素比重的同位素擴散速率快[13]。在多個同位素體系(Li、Ca、Fe和Mg等)中均觀察到了擴散作用引起的動力學同位素分餾[41–46]。元素或同位素的擴散作用可分為熱擴散和化學擴散, 高溫過程中, 由溫度梯度導致的濃度差異所引起的擴散作用稱為熱擴散作用; 由化學勢梯度引起的擴散作用稱為化學擴散。雖然熱擴散在實驗室條件下也能看到可分辨的同位素分餾, 但是熱擴散作用在天然物質中引起的同位素分餾大小只與溫度梯度呈線性關系[47–48], 然而天然環境中熱擴散速度比化學擴散快幾個數量級以上, 溫度梯度消散的太快了, 僅限于部分特殊的環境才能保存溫度梯度(如偉晶巖脈和圍巖, 基性巖墻和圍巖)[43,47]。而化學擴散的驅動力為化學濃度差, 只要有化學勢差異存在, 化學擴散便會持續進行, 最終到達化學平衡。例如, 夏威夷玄武巖中橄欖石和熔體中Mg-Fe相互擴散作用[49], 地幔橄欖巖和熔體之間化學擴散作用導致巖石圈地幔Fe同位素不均一[14]。在高溫Ca同位素體系中, 除了CaO自身濃度梯度, 熔體基質的性質(SiO2、Na2O和Al2O3等)也對Ca的擴散速度和同位素分餾大小有影響[45–46]。Watkins.[45]選擇不同Al2O3含量的基性熔體與流紋質熔體進行熔融擴散實驗發現, 在拉斑玄武巖-流紋巖體系中44/40Ca負漂約5‰～6‰, 玄武巖-流紋巖體系44/40Ca負漂約3‰, 而堿玄巖-流紋巖體系則出現44/40Ca正漂約2‰, 隨著Al2O3含量增加44Ca趨向于富集在酸性熔體中。不同CaO-Na2O-SiO2體系的擴散實驗表明, 在CaO和Na2O梯度較小的體系中擴散較快, Ca同位素分餾相對較大(2.5‰), 而在CaO-SiO2對擴散和Na2O均勻分布的體系中擴散較慢, Ca同位素分餾較小(1.3‰)[46]。此外, 熔體與橄欖巖反應也會導致明顯的Ca同位素動力學分餾。由于熔體和橄欖巖之間存在Ca濃度梯度, Ca從熔體中擴散到斜方輝石形成單斜輝石, 導致全巖Ca同位素組成變輕[14]。

1.4 40K衰變對Ca同位素的影響

還有一點需要說明的是,40K可以通過 ‘β-’衰變成40Ca, 這使得在討論Ca同位素分餾時需要對40K衰變的影響進行評估。評估辦法通常是選擇K含量或者K/Ca比值高的樣品進行不加稀釋劑測試, 所得40/44Ca的值大于40/44Ca標樣(47.162)表明, 樣品可能受到40K衰變的影響, 值偏離越遠衰變影響越大。雖然40K豐度很低, 但對于一些含鉀高的礦物或巖石, 例如鉀長石和A型花崗巖, 其放射性成因40Ca的累積會使得44/40Ca值偏低[22], 并偏離質量相關分餾線。對于K/Ca比值為50, 年齡為300 Ma的巖石足以引起44/40Ca值1‰的變化[7]。像Sr和Nd同位素一樣, 放射性成因的40Ca也有示蹤殼幔物質相互作用的潛力[16,50]。對于地殼長英質巖石而言, 通常具有高的Na2O和K2O含量, 而40K的衰變會累積放射性成因的40Ca, 當這些長英質物質經過俯沖帶或拆沉作用再循環進入地幔時, 會攜帶放射性成因40Ca的信息, 因此, 放射性成因Ca同位素可以指示殼幔物質循環。已有Ca同位素研究表明, 在大陸內部, 一些巖漿系統可回收70%以上的古老地殼, 有些則不到30%, 可能取決于下地殼的熱狀態或化學成分, 以及地幔提供的巖漿通量[50]。此外,40K-40Ca衰變體系還被用于古老的富鉀巖石的定年研究[16]。

2 巖漿過程中的Ca同位素分餾

巖漿從產生到成巖經歷了一系列復雜的地質過程, 包括部分熔融、分離結晶、流體交代、同化混染和巖漿混合等。哪些過程會造成Ca同位素分餾？接下來一一進行討論。

2.1 部分熔融

部分熔融是巖漿形成的重要過程, 壓力降低、溫度升高和揮發組分的加入都可能降低巖石的固相線, 從而引發部分熔融[51]。例如在俯沖帶附近, 富水的洋殼和蛇紋巖會發生脫水作用, 導致地幔楔發生部分熔融; 地幔柱和洋中脊由于軟流圈物質上涌使其溫度升高也會誘發部分熔融; 此外在大陸伸展構造背景下, 軟流圈上涌或玄武巖巖漿底侵、壓力減小和降低固相線也會引發部分熔融, 這也是A型花崗巖的形成機制之一[52]。為了研究Ca同位素在部分熔融過程中的地球化學行為, Zhu.[17]對太平洋的洋中脊玄武巖(MORB)進行了Ca同位素研究。MORB是上地幔部分熔融的產物, 受外界干擾因素較少, 是研究部分熔融過程最好的天然物質。數據表明, MORB的Ca同位素組成從(0.75±0.05)‰變化到了(0.86±0.03)‰, 低于前人估計的BSE值。為什么MORB的Ca同位素會比BSE值低？已有研究表明, 橄欖石、斜方輝石和單斜輝石的Ca同位素組成存在差異[15,18,19]。在部分熔融過程中, 富集40Ca的橄欖石和單斜輝石優先被消耗, 導致熔體相富集輕同位素, 而殘余固相變重。MgO和44/40Ca的負相關也說明熔融比例不同會造成Ca同位素分餾[20–21]。因此, 在部分熔融過程中Ca同位素會發生同位素分餾且輕同位素優先進入熔體相。MORB由地幔發生高程度部分熔融(約10%)而形成, 對基性玄武巖進行研究時應考慮對部分熔融影響進行有效扣除。但是對于大多數中酸性巖而言, 從原巖部分熔融到后期成巖的整個過程中, 部分熔融程度都相對較低(<1%), 能造成的Ca同位素分餾相對較小。

2.2 分離結晶

在巖漿形成之后, 還會經歷長時間的分離結晶過程, 那么分離結晶是否會造成Ca同位素的分餾？如果造成分餾, 那Ca同位素在分離結晶過程中的同位素分餾又如何？主要影響因素有哪些？首先, 在巖石中Ca通常以Ca2+的形式賦存于各種礦物中, 在巖漿演化過程中通常會發生斜方輝石、單斜輝石、斜長石、角閃石、磷灰石和榍石等富Ca礦物的分離結晶。Ryu.[22]對1.7 Ga的花崗閃長巖進行了全巖和單礦物的Ca同位素研究, 發現不同富Ca礦物之間存在明顯的Ca同位素組成差異(圖4), 其中鉀長石和綠泥石由于受到40K衰變的影響而富集放射性成因40Ca, 偏離了質量相關分餾線[22]。富Ca礦物中斜長石具有最輕的同位素組成, 方解石、榍石、磷灰石、角閃石和綠簾石依次變重。矯正鉀衰變影響之后的鉀長石和斜長石有著差不多的Ca同位素組成, 而黑云母具有最重的Ca同位素組成。對于造成礦物之間差異的原因目前還不清楚, 一種較為合理的解釋是, 元素在礦物中的配位數是影響化學鍵能和鍵長的主要因素, 配位數越低、鍵長越短和鍵能越強, 趨向于富集重同位素組成。但也有可能是高溫平衡分餾作用或者受到了后期差異組分的改造[13]。此外, He.[53]發現高 (Dy/Yb)N比值的巖石樣品具有相對較低的44/42Ca值, 暗示富Ca石榴子石可能富集輕Ca同位素組成。總的來說, 不同礦物之間Ca同位素組成的差異暗示分離結晶過程中可能存在Ca同位素分餾。在MORB中觀察到的44/40Ca與Eu/Eu*的負相關關系也表明斜長石可能會導致殘余相變重, 這與斜長石具有較輕的同位素組成相一致[17]。

圖4 不同含Ca單礦物δ44/40Ca值vs. δ44/42Ca值圖解, 數據來自文獻[22]

2.3 交代或混合作用

除了部分熔融和分離結晶作用, 熔流體交代、同化混染和巖漿混合作用也廣泛存在于中酸性巖的形成過程中。當源區在形成或后期受到其他差異組分的改造, 將會引起同位素分餾, 可以用Ca同位素結合其他元素或同位素進行定量模擬發生交代或混合的比例, 這也是同位素示蹤的基本原理。已有研究表明, 硅酸巖和碳酸巖熔體交代均會引起可分辨的Ca同位素分餾。因為石榴子石相對富集輕同位素組成, 在加厚的榴輝巖相下地殼部分熔融所形成的埃達克質巖(源區有石榴子石殘留)中會導致Ca同位素分餾[54]。此外, 硅酸巖下地殼與地幔橄欖巖發生熔融反應也會造成Ca同位素分餾[55]。碳酸鹽巖相對硅酸巖具有更輕的Ca、Mg同位素組成, 碳酸鹽巖交代將會引起較大的同位素分餾, 這也使得Ca、Mg同位素可以作為判斷是否有沉積碳酸巖介入俯沖帶的證據之一, 而Ca、Mg和C元素的耦合, 也使得Ca、Mg同位素顯示出了示蹤深部碳循環的潛力[7,16,56]。對中國東部大范圍的Mg同位素研究已經表明, 在晚白堊世時期, 中國東部都受到古老太平洋俯沖的影響, 其中沉積碳酸巖隨著俯沖作用被帶到深部, 再發生部分熔融經過地幔楔所形成的玄武巖具有更低的26Mg值[57–59]。大部分碳酸鹽中的Ca含量比Mg含量更高, 因此可以為碳酸鹽再循環提供更加直接的證據。Huang.[19]對夏威夷洋島玄武巖研究提出, 大約4%的古老沉積碳酸巖加入夏威夷地幔柱使得其44/40Ca值略低于BES估計值。同樣的, 俯沖帶中存在沉積碳酸巖加入的Ca同位素證據也在騰沖和西藏研究中被發現[23–24]。

3 變質作用過程中的Ca同位素分餾

Ca同位素除了巖漿演化過程和沉積物的加入, 變質作用過程中也會發生同位素分餾作用。俯沖作用過程中洋殼的脫水變質有大量的流體釋放, 而這些變質流體可能會引起Ca同位素分餾。但目前對于變質巖的Ca同位素研究甚少, 僅少數標樣和幾個天然樣品的報道。僅有的研究表明, 在未蝕變的藍片巖相巖石、周圍的榴輝巖和石英碳酸鹽巖脈中, Ca同位素更趨向于富集在石英碳酸鹽巖脈中(圖5), 指示可能有一種富重Ca同位素流體的交代作用[60]。富重Ca同位素源區可能來自于部分蝕變的巖石圈地幔, 含水量為1%～2%, 蛇紋石蝕變10%～15%。但是即使蛇紋巖具有高的44/40Ca值, 其CaO含量很低, 能否脫水產生可分辨的同位素分餾還不是很明確[60]。另一種可能是和蝕變洋殼具有相當44/40Ca值的流體在運移過程中因為碳酸鹽巖沉淀導致Ca同位素變重?？紤]到二者44/40Ca差值約0.6‰, 估算大約有40%～70%溶解態的Ca需要被沉淀才能解釋[60]。Lu.[61]進一步對俯沖帶變質過程研究表明, 在俯沖帶變質脫水過程中, 除了局部存在碳酸質流體交代作用會引起Ca同位素分餾外, Ca同位素不會發生明顯的同位素分餾, 而俯沖循環的地殼物質可能是造成上地幔Ca同位素不均一的原因之一。

圖5 變質作用過程中Ca同位素分餾圖解, 數據來自文獻[60]

4 地外樣品中的Ca同位素分餾作用

Ca同位素在行星地球化學領域的最早應用可以追溯到 Russell.[4]在 1978年發現的富鈣、鋁難熔包體(calcium-aluminum-rich inclusions, CAI))中富集48Ca同時虧損46Ca。隨后Simon.[1,62]開始分析頑火輝石球粒隕石與碳質隕石的Ca同位素組成差異以及碳質球粒隕石CV3 (carbonaceous vigarano meteorite 3)中的包裹體Ca同位素組成, 發現頑火輝石球粒隕石比碳質球粒隕石相對富集重Ca同位素組成, 對比地球、月球、火星以及其他隕石后指出, 頑火輝石球粒隕石更能代表類地行星形成的主要材料。Huang.[63]測試了Allende球粒隕石中6個CAIs、碳質球粒隕石、普通球粒隕石以及頑火輝石球粒隕石的Ca同位素組成結果指出, 在隕石中存在質量相關分餾與非質量相關分餾[64]; Valdes.[2]系統地討論了頑火輝石球粒隕石、普通球粒隕石、地球與月球樣品的Ca同位素組成差異, 進一步討論地球的初始物質、頑火輝石和普通球粒隕石均是形成地球的初始材料之一, 而碳質球粒隕石CO (CarbonaceousOrgueil)不是; Schiller.[65]通過Ca同位素的非質量相關分餾系統討論了類地行星早期的形成與演化。

5 總 結

(1) Ca同位素不僅在不同礦物之間存在可分辨的Ca同位素分餾, 而且在巖漿演化和變質作用過程中會發生同位素分餾。不同礦物間的分餾可能受控于Ca—O化學鍵的強度、溫度和礦物自身的組分控制, 如Ca/Mg比值。40K衰變也是造成不同K/Ca比值礦物之間Ca同位素組成差異的原因之一。在擴散作用中由于CaO濃度差異, 也會造成Ca同位素分餾, 同時Ca同位素分餾也與熔體的性質有關(SiO2、Na2O和Al2O3)。部分熔融過程中輕的Ca同位素優先進入熔體相, 導致殘余固相Ca同位素組成變重。不同含Ca礦物之間的Ca同位素組成差異說明在分離結晶過程中由于礦物分離結晶的先后順序不同會引起Ca同位素分餾。變質作用過程中變質流體的交代作用和流體運移過程中碳酸鈣沉淀都可能引起Ca同位素分餾。

(2) K-Ca衰變體系既可以對某些富K巖石或礦物進行同位素定年, 放射性成因40Ca又能示蹤殼幔物質循環比例, 在示蹤長英質巖漿中與超爆發、大陸裂谷有關的地殼貢獻方面具有廣闊的應用前景。

(3) 高溫過程中Ca同位素存在明顯的同位素分餾, 對于分餾的機理和影響因素還不是很清楚, 高溫過程中的應用也相對較少。Ca同位素在地質過程中的分餾展現出了廣泛的應用潛力, 需要更多的天然樣品結合實驗和理論計算等多方向多手段去完善Ca同位素體系, 用Ca同位素這把“鑰匙”去認識地球的未知領域。

感謝審稿人對文章提出的寶貴意見和建議。此外, 文章寫作過程得到李賀博士以及吳韋同學的支持和幫助, 在此一并感謝。

[1] Simon J I, Depaolo D J. Stable calcium isotopic composition of meteorites and rocky planets[J]. Earth Planet Sci Lett, 2010, 289(3): 457–466.

[2] Valdes M C, Moreira M, Foriel J, Moynier F. The nature of Earth’s building blocks as revealed by calcium isotopes[J]. Earth Planet Sci Lett, 2014, 394: 135–145.

[3] Heumann K G, Luecke W. Calcium isotope ratios in carbonate rocks[J]. Earth Planet Sci Lett, 1973, 20(3): 341–346.

[4] Russell W A, Papanastassiou D A. Calcium isotope fractionation in ion-exchange chromatography[J]. Anal Chem, 1978, 50(8): 1151–1154.

[5] Bagard M L, Schmitt A D, Chabaux F, Pokrovsky O S, Viers J, Stille P, Labolle F, Prokushkin A S. Biogeochemistry of stable Ca and radiogenic Sr isotopes in a larch-covered permafrost- dominated watershed of Central Siberia[J]. Geochim Cosmochim Acta, 2013, 114: 169–187.

[6] Fantle M S, DePaolo D J. Ca isotopes in carbonate sediment and pore fluid from ODP Site 807A: The Ca2+(aq)–calcite equilibrium fractionation factor and calcite recrystallization rates in Pleistocene sediments[J]. Geochim Cosmochim Acta, 2007, 71(10): 2524–2546.

[7] Fantle M S, Tipper E T. Calcium isotopes in the global biogeochemical Ca cycle: Implications for development of a Ca isotope proxy[J]. Earth-Science Reviews, 2014, 129: 148– 177.

[8] Farka? J, Buhl D, Blenkinsop J, Veizer J. Evolution of the oceanic calcium cycle during the late Mesozoic: Evidence from44/40Ca of marine skeletal carbonates[J]. Earth Planet Sci Lett, 2007, 253(1/2): 96–111.

[9] Holmden C, Bélanger N. Ca isotope cycling in a forested ecosystem[J]. Geochim Cosmochim Acta, 2010, 74(3): 995– 1015.

[10] Jacobson A D, Grace Andrews M, Lehn G O, Holmden C. Silicate versus carbonate weathering in Iceland: New insights from Ca isotopes[J]. Earth Planet Sci Lett, 2015, 416: 132– 142.

[11] Jacobson A D, Holmden C.44Ca evolution in a carbonate aquifer and its bearing on the equilibrium isotope fractionation factor for calcite[J]. Earth Planet Sci Lett, 2008, 270(3/4): 349–353.

[12] Tipper E T, Galy A, Bickle M J. Riverine evidence for a fractionated reservoir of Ca and Mg on the continents: Implicationsfor the oceanic Ca cycle[J]. Earth Planet Sci Lett, 2006, 247(3/4): 267–279.

[13] Gussone N, Schmitt A D, Heuser A, Wombacher F, Dietzel M, Tipper E, Schiller M. Calcium Stable Isotope Geochemistry[M]. Springer: Heidelberg, 2016: 1–272.

[14] Huang S C, Farka? J, Jacobsen S B. Calcium isotopic fractionation between clinopyroxene and orthopyroxene from mantle peridotites[J]. Earth Planet Sci Lett, 2010, 292(3/4): 337–344.

[15] Kang J T, Ionov D A, Liu F, Zhang C L, Golovin A V, Qin L P, Zhang Z F, Huang F. Calcium isotopic fractionation in mantle peridotites by melting and metasomatism and Ca isotope composition of the Bulk Silicate Earth[J]. Earth Planet Sci Lett, 2017, 474: 128–137.

[16] Depaolo D J. Calcium isotopic variations produced by biological, kinetic, radiogenic and nucleosynthetic processes[J]. Rev Mineral Geochem, 2004, 55(1): 255–288.

[17] Zhu H L, Liu F, Li X, Wang G Q, Zhang Z F, Sun W D. Calcium isotopic compositions of normal mid-ocean ridge basalts from the southern Juan de Fuca Ridge[J]. J Geophys Res Solid Earth, 2018, 123(2): 1303–1313.

[18] Kang J T, Zhu H L, Liu Y F, Liu F, Wu F, Hao Y T, Zhi X C, Zhang Z F, Huang F. Calcium isotopic composition of mantle xenoliths and minerals from Eastern China[J]. Geochim Cosmochim Acta, 2016, 174: 335–344.

[19] Huang S C, Farka? J, Jacobsen S B. Stable calcium isotopic compositions of Hawaiian shield lavas: Evidence for recycling of ancient marine carbonates into the mantle[J]. Geochim Cosmochim Acta, 2011, 75(17): 4987–4997.

[20] Amini M, Eisenhauer A, B?hm F, Holmden C, Kreissig K, Hauff F, Jochum K P. Calcium isotopes (44/40Ca) in MPI-DING reference glasses, USGS rock powders and various rocks: Evidence for Ca isotope fractionation in terrestrial silicates[J]. Geost Geoanal Res, 2009, 33(2): 231–247.

[21] Liu F, Zhu H L, Li X, Wang G Q, Zhang Z F. Calcium isotopic fractionation and compositions of geochemical reference materials[J]. Geost Geoanal Res, 2017, 41(4): 675–688.

[22] Ryu J S, Jacobson A D, Holmden C, Lundstrom C, Zhang Z F. The major ion,44/40Ca,44/42Ca, and26/24Mg geochemistry of granite weathering at pH=1 and= 25 ℃: Power-law processes and the relative reactivity of minerals[J]. Geochim Cosmochim Acta, 2011, 75(20): 6004–6026.

[23] Liu F, Li X, Wang G Q, Liu Y F, Zhu H L, Kang J T, Huang F, Sun W D, Xia X P, Zhang Z. Marine carbonate component in the mantle beneath the southeastern Tibetan Plateau: Evidence from magnesium and calcium isotopes[J]. J Geophys Res Solid Earth, 2017, 122(12): 9729–9744.

[24] 劉峪菲, 許繼峰, 張兆峰, 王桂琴, 陳建林, 黃豐, 祝紅麗, 劉芳. 青藏高原拉薩地塊中西部超鉀質巖Ca-Mg同位素特征及其地質意義[J]. 地質學報, 2018, 92(3): 545–559.

Liu Yu-fei, Xu Ji-feng, Zhang Zhao-feng, Wang Gui-qin, Chen Jian-lin, Huang Feng, Zhu Hong-li, Liu Fang. Ca-Mg isotopic compositions of ultra potassic volcanic rocks in the Lhasa Terrane, southern Tibet and their geological implications[J]. Acta Geol Sinica, 2018, 92(3): 545–559 (in Chinese with English abstract).

[25] Hindshaw R S, Reynolds B C, Wiederhold J G, Kretzschmar R, Bourdon B. Calcium isotopes in a proglacial weathering environment: Damma glacier, Switzerland[J]. Geochim Cosmochim Acta, 2011, 75(1): 106–118.

[26] Zhao X M, Zhang Z F, Huang S C, Liu Y F, Li X, Zhang H F. Coupled extremely light Ca and Fe isotopes in peridotites[J]. Geochim Cosmochim Acta, 2017, 208: 368-380.

[27] Tang J W, Dietzel M, B?hm F, K?hler S J, Eisenhauer A. Sr2+/ Ca2+and44Ca/40Ca fractionation during inorganic calcite formation: II. Ca isotopes[J]. Geochim Cosmochim Acta, 2008, 72(15): 3733–3745.

[28] Lemarchand D, Wasserburg G J, Papanastassiou D A. Rate-controlled calcium isotope fractionation in synthetic calcite[J]. Geochim Cosmochim Acta, 2004, 68(22): 4665– 4678.

[29] Feng L P, Zhou L, Yang L, Depaolo D J, Tong S Y, Liu Y S, Owens T L, Gao S. Calcium isotopic compositions of sixteen USGS reference materials[J]. Geost Geoanal Res, 2017, 41(1): 93–106.

[30] Zhu H L, Liu F, Li X, An Y J, Wang G Q, Zhang Z F. A “peak cut” procedure of column separation for calcium isotope measurement using the double spike technique and thermal ionization mass spectrometry (TIMS) [J]. J Anal Atom Spectr, 2018, 33(4): 547–554.

[31] Zhu H L, Zhang Z F, Wang G Q, Liu Y F, Liu F, Li X, Sun W D. Calcium isotopic fractionation during ion-exchange column chemistry and thermal ionisation mass spectrometry (TIMS) determination[J]. Geost Geoanal Res, 2016, 40(2): 185–194.

[32] 劉芳. 鈣同位素分析方法及其地質應用——以騰沖火山巖和粵北碳酸鹽風化殼為例[D]. 廣州: 中國科學院廣州地球化學研究所, 2018.

Liu Fang. Calcium isotope analytical method and geological applications: Case studies on Tengchong volcanic rocks and the weathering crust in northern Guangdong. Guangzhou: Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, 2018 (in Chinese with English abstract)

[33] 劉芳, 祝紅麗, 譚德燦, 劉峪菲, 康晉霆, 朱建明, 王桂琴, 張兆峰. 熱電離質譜測定鈣同位素過程中雙稀釋劑的選擇[J]. 質譜學報, 2016, 37(4): 310–318.

Liu F, Zhu H L, Tan D C, Liu Y F, Kang J T, Zhu J M, Wang G Q, Zhang Z F. Optimization of calcium double-spike on measuring calcium isotopic compositions by thermal ionization mass spectrometer (TIMS) [J]. J Chinese Mass Spectr Soc, 2016, 37(4): 310–318 (in Chinese with English abstract).

[34] 劉峪菲, 祝紅麗, 劉芳, 王桂琴, 許繼峰, 張兆峰. 鈣同位素化學分離方法研究[J]. 地球化學, 2015, 44(5): 469–476.

Liu Y F, Zhu H L, Liu F, Wang G Q, Xu J F, Zhang Z F. Methodological study of chemical separation of calcium for TIMS measurements [J]. Geochimica, 2015, 44(5): 469–476 (in Chinese with English abstract).

[35] 祝紅麗. 鈣同位素分析方法及應用[D]. 廣州: 中國科學院廣州地球化學研究所, 2016.

Zhu Hong-li. Calcium isotope analytical method and its geological applications[D]. Guangzhou: Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, 2016 (in Chinese with English abstract).

[36] Feng C Q, Qin T, Huang S C, Wu Z Q, Huang F. First-principles investigations of equilibrium calcium isotope fractionation between clinopyroxene and Ca-doped orthopyroxene[J]. Geochim Cosmochim Acta, 2014, 143: 132–142.

[37] Chen C, Liu Y, Feng L, Foley S F, Zhou L, Ducea, M N, Hu Z. Calcium isotope evidence for subduction-enriched lithospheric mantle under the northern North China Craton[J]. Geochim Cosmochim Acta, 2018, 238: 55–67.

[38] Wang W Z, Zhou C, Qin T, Kang J T, Huang S C, Wu Z Q, Huang F. Effect of Ca content on equilibrium Ca isotope fractionation between orthopyroxene and clinopyroxene[J]. Geochim Cosmochim Acta, 2017, 219: 44–56.

[39] Wang W Z, Qin T, Zhou C, Huang S C, Wu Z Q, Huang F. Concentration effect on equilibrium fractionation of Mg-Ca isotopes in carbonate minerals: Insights from first-principles calculations[J]. Geochim Cosmochim Acta, 2017, 208: 185– 197.

[40] Huang F, Zhou C, Wang W Z, Kang J T, Wu Z Q. First-principlescalculations of equilibrium Ca isotope fractionation: Implicationsfor oldhamite formation and evolution of lunar magma ocean[J]. Earth Planet Sci Lett, 2019, 510: 153–160.

[41] Richter F M, Dauphas N, Teng F Z. Non-traditional fractionation of non-traditional isotopes: Evaporation, chemical diffusion and Soret diffusion[J]. Chem Geol, 2009, 258(1/2): 92–103.

[42] Richter F M, Davis A M, Depaolo D J, Watson E B. Isotope fractionation by chemical diffusion between molten basalt and rhyolite[J]. Geochim Cosmochim Acta, 2003, 67(20): 3905– 3923.

[43] Richter F M, Watson E B, Mendybaev R, Dauphas N, Georg B, Watkins J, Valley J. Isotopic fractionation of the major elementsof molten basalt by chemical and thermal diffusion[J]. Geochim Cosmochim Acta, 2009, 73(14): 4250–4263.

[44] Richter F M, Bruce Watson E, Mendybaev R A, Teng F Z, Janney P E. Magnesium isotope fractionation in silicate melts by chemical and thermal diffusion[J]. Geochim Cosmochim Acta, 2008, 72(1): 206–220.

[45] Watkins J M, Depaolo D J, Huber C, Ryerson F J. Liquid composition-dependence of calcium isotope fractionation during diffusion in molten silicates[J]. Geochim Cosmochim Acta, 2009, 73(24): 7341–7359.

[46] Watkins J M, Depaolo D J, Ryerson F J, Peterson B T. Influence of liquid structure on diffusive isotope separation in molten silicates and aqueous solutions[J]. Geochim Cosmochim Acta, 2011, 75(11): 3103–3118.

[47] Huang F, Chakraborty P, Lundstrom C C, Holmden C, Glessner J J G, Kieffer S W, Lesher C E. Isotope fractionation in silicate melts by thermal diffusion[J]. Nature, 2012, 464(7287): 396–400.

[48] Li X F, Liu Y. A theoretical model of isotopic fractionation by thermal diffusion and its implementation on silicate melts[J]. Geochim Cosmochim Acta, 2015, 154: 18–27.

[49] Teng F Z, Dauphas N, Helz R T, Gao S, Huang S C. Diffusion- driven magnesium and iron isotope fractionation in Hawaiian olivine[J]. Earth Planet Sci Lett, 2011, 308(3/4): 317–324.

[50] Mills R D, Simon J I, Depaolo D J. Calcium and neodymium radiogenic isotopes of igneous rocks: Tracing crustal contributions in felsic magmas related to super-eruptions and continental rifting[J]. Earth Planet Sci Lett, 2018, 495: 242–250.

[51] 孫鼎, 彭亞嗚. 火成巖石學[M]. 北京: 地質出版社, 1985: 1–524.

Sun Ding, Peng Ya-wu. Igneous Lithography[M]. Beijing: Geological Publishing House, 1985: 1–524 (in Chinese).

[52] Bonin B. A-type granites and related rocks: Evolution of a concept, problems and prospects[J]. Lithos, 2007, 97(1/2): 1–29.

[53] He Y S, Wang Y, Zhu C W, Huang S C, Li S G. Mass-independent and mass-dependent Ca isotopic compositions of thirteen geological reference materials measured by thermal ionisation mass spectrometry[J]. Geost Geoanal Res, 2017, 41(2): 283– 302.

[54] Wang Y, He Y S, Wu H J, Zhu C W, Huang S C, Huang J. Calcium isotope fractionation during crustal melting and magma differentiation: Granitoid and mineral-pair perspectives[J]. Geochim Cosmochim Acta, 2019, 259: 37–52.

[55] Kang J T, Ionov D A, Zhu H L, Liu F, Zhang Z F, Liu Z, Huang F. Calcium isotope sources and fractionation during melt-rock interaction in the lithospheric mantle: Evidence from pyroxenites, wehrlites, and eclogites[J]. Chem Geol, 2019, 524: 272–282.

[56] 祝紅麗, 張兆峰, 劉峪菲, 劉芳, 康晉霆. 鈣同位素地球化學綜述[J]. 地學前緣, 2015, 22(5): 44–53.

Zhu Hong-li, Zhang Zhao-feng, Liu Yu-fei, Liu Fang, Kang Jin-ting. Calcium isotope geochemistry review[J]. Earth Sci Front, 2015, 22(5): 44–53 (in Chinese with English abstract).

[57] Huang J, Li S G, Xiao Y L, Ke S, Li W Y, Tian Y. Origin of low26Mg Cenozoic basalts from South China Block and their geodynamic implications[J]. Geochim Cosmochim Acta, 2015, 164: 298–317.

[58] Li S G, Yang W, Ke S, Meng X N, Tian H C, Xu L J, He Y S, Huang J, Wang X C, Xia Q, Sun W D, Yang X Y, Ren Z Y, Wei H Q, Liu Y S, Meng F C, Yan J. Deep carbon cycles constrained by a large-scale mantle Mg isotope anomaly in eastern China[J]. Nat Sci Rev, 2017, 4(1): 111–120.

[59] Yang W, Teng F Z, Zhang H F, Li S G. Magnesium isotopic systematics of continental basalts from the North China craton: Implications for tracing subducted carbonate in the mantle[J]. Chem Geol, 2012, 328: 185–194.

[60] John T, Gussone N, Podladchikov Y Y, Bebout G E, Dohmen R, Halama R, Klemd R, Magna T, Seitz H M. Volcanic arcs fed by rapid pulsed fluid flow through subducting slabs[J]. Nat Geosci, 2012, 5(7): 489–492.

[61] Lu W N, He Y S, Wang Y, Ke S. Behavior of calcium isotopes during continental subduction recorded in meta-basaltic rocks[J]. Geochim Cosmochim Acta, 2020, 278: 392–404.

[62] Simon J I, Depaolo D J, Moynier F. Calcium isotope compositionof meteorites, earth, and mars[J]. Astrophys J, 2009, 702: 707–715.

[63] Huang S C, Farka? J, Yu G, Petaev M I, Jacobsen S B. Calcium isotopic ratios and rare earth element abundances in refractory inclusions from the Allende CV3 chondrite[J]. Geochim Cosmochim Acta, 2012, 77: 252–265.

[64] Huang S C, Jacobsen S B. Calcium isotopic compositions of chondrites[J]. Geochim Cosmochim Acta, 2017, 201: 364–376.

[65] Schiller M, Bizzarro M, Fernandes V A S M. Isotopic evolution of the protoplanetary disk and the building blocks of Earth and the Moon[J]. Nature, 2018, 555(7697): 507–510.

Application and prospects of calcium isotopes in high temperature processes

LUO Ze-bin1,3, LING Ming-xing2*, AN Ya-jun1, LIU Fang1and ZHANG Zhao-feng1

1.State Key Laboratory of Isotope Geochemistry, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou 510640, China; 2. State Key Laboratory of Nuclear Resources and Environment, East China University of Technology, Nanchang 330013, China; 3. University of Chinese Academy of Sciences, Beijing 100049, China

The main factors restricting the development of Ca isotopes in high temperature processes are analytical techniques, Ca isotopic composition of the end members, and fractionation mechanisms. With the application of mass spectrometry and double spikes, the analytical method can meet the requirements of high temperature process tracing. However, research on the Ca isotopic composition of the end members is insufficient, and the fractionation mechanisms are controversial. Therefore, it is imperative to identify the Ca isotopic composition of the end members and to clarify the mechanisms of Ca isotope fractionation. As the most important end member, the Ca isotopic composition of bulk silicate Earth varies from 0.94‰ to 1.05‰. The inter-mineral fractionation of Ca isotopes is primarily controlled by the strength of the Ca—O bond. In addition, the difference in temperature and mineral composition also leads to the fractionation of Ca isotopes. Both thermal and chemical diffusion can cause Ca isotope fractionation, but thermal diffusion only occurs in special environments. Chemical diffusion is controlled by the chemical potential of CaO, and the chemical potential is controlled by the composition of the system, temperature, and pressure. Light isotopes tend to be enriched in the melt during partial melting and melt extraction. The Ca isotopic heterogeneity of felsic minerals leads to Ca isotopic fractionation in the process of fractional crystallization. Both fluid metasomatism and carbonate precipitation during metamorphism lead to Ca isotopic fractionation. The K-Ca decay system makes it possible to date the ancient K-rich rocks or minerals. Radioactive Ca can also be used to trace the source region and crust-mantle material cycle.

calcium isotopes; high temperature geological processes; magmatic evolution; metamorphism; geological application

P597

A

0379-1726(2021)05-0433-09

10.19700/j.0379-1726.2021.05.001

2019-12-31;

2020-04-19;

2020-05-09

東華理工大學核資源與環境國家重點實驗室自主基金重點項目(2020Z03); 江西省“雙千計劃”創新領軍人才(青年類)項目; 東華理工大學高層次人才啟動基金

羅澤彬(1994–), 男, 博士研究生, 地球化學專業。E-mail: zebin_l@163.com

通訊作者(Corresponding author):LING Ming-xing, E-mail: mxling@ecut.edu.cn; Tel: +86-18665000415