一种定量微藻利用源于硅酸岩风化的无机碳能力方法

摘要

本发明公开一种定量微藻利用源于硅酸岩风化的无机碳能力方法。通过测定硅酸盐岩风化过程中产生的碳酸氢根离子的碳同位素值、待测微藻二氧化碳同化过程中的稳定碳同位素分馏值、微藻的碳酸氢根离子途径份额以及微藻仅利用硅酸盐岩风化产生的碳酸氢根离子所形成的藻体碳同位素值，计算出微藻利用硅酸盐岩风化产生的无机碳源的利用份额，依据微藻利用硅酸盐岩风化产生的无机碳源的利用份额以及微藻增值倍数，最终定量出微藻利用源于硅酸岩风化的无机碳能力，为寻找遗失的碳汇提供新的途径。

说明书

技术领域

本发明涉及一种定量微藻利用源于硅酸岩风化的无机碳能力方法，属于全球变化以及生态环境系统监测领域。

背景技术

全球大气二氧化碳(CO₂)浓度逐渐升高，全球气候变暖，全球碳循环已成为国际科学界关注的热点。全球碳储存库主要有化石燃料、陆地、大气、海洋，碳酸盐岩等，随着人类对化石燃料资源的开采燃烧及其土地开荒，扰动了前农业时期的全球碳循环。地球系统的碳循环影响着全球气候变化，它对世界经济、社会和生态环境等产生了重大影响，大气CO₂浓度升高对世界各国经济的可持续发展和国家安全，等带来严峻挑战。钙镁硅酸盐岩的风化作用形成较为稳定的碳酸氢根离子(HCO₃^-)，对气候调节起决定性的作用。

硅酸盐在风化作用中，微藻作为古老的生命形式，贯穿整个地质，微藻能利用风化产生的碳酸氢根离子(HCO₃^-)，形成“生物泵”，加速岩石的风化。定量出微藻利用源于硅酸岩风化的无机碳能力，可指示微藻对硅酸岩的生物溶蚀能力，为寻找遗失的碳汇提供新的途径。

发明内容

本发明要解决的技术问题是：提供一种定量微藻利用源于硅酸岩风化的无机碳能力方法，实现微藻利用源于硅酸岩风化无机碳的光合碳汇的定量，可以补充部分被忽视的“遗失碳汇”，填补碳汇计算的空白。

本发明采取以下技术方案：它包括以下步骤：

第一，利用添加碳酸氢钠的浓度为16mM及以上、同时附加10mM乙酰唑胺也即AZ双向同位素标记培养待测微藻，获取待测微藻二氧化碳同化过程中的稳定碳同位素分馏值Δ；分步骤如下：

分步一、首先选择两种δ¹³C值差值大于10‰的碳酸氢钠作为同位素标记1和同位素标记2分别添加到检测培养液中，同时在同样的培养条件下培养起始生物量相同的待测微藻；检测培养液为添加碳酸氢钠的浓度为16mM及以上、同时附加10mM乙酰唑胺也即AZ的藻类培养液；测定待测微藻的起始生物量M_O和微藻的起始稳定碳同位素组成δ¹³C的值δ_O。同位素标记1的碳酸氢钠的δ¹³C值为δ_C1，同位素标记2的碳酸氢钠的δ¹³C值为δ_C2；

分步二、其次将待测微藻在添加同位素标记的检测培养液中培养4天后，测定待测微藻的最终生物量M_N和在两种同位素标记的检测培养液培养的被考察微藻的稳定碳同位素组成δ¹³C的值δ_S1和δ_S2；

分步三、测定此培养期间每天培养环境中大气二氧化碳的日平均稳定碳同位素组成；由此计算出在培养期间，培养环境中大气二氧化碳的平均稳定碳同位素组成δ¹³C的值δ_air；

分步四、获取待测微藻增值倍数P＝M_N/M_O；

分步五、将待测微藻的起始稳定碳同位素组成δ_O、待测微藻增值倍数P以及在两种同位素标记的检测培养液培养的被考察微藻的稳定碳同位素组成δ_S1或δ_S2代入方程：和计算出在两种同位素标记的检测培养液培养的被考察微藻的新生藻体稳定碳同位素组成δ¹³C的值δ_T1和δ_T2；

分步六、将δ_T1、δ_T2、δ_C1和δ_C2代入到方程：计算出待测微藻在碳酸氢根离子利用途径完全被抑制的情况下利用添加的无机碳源的份额f_B-A；

分步七、将f_B-A、δ_air以及δ_T1和δ_C1或δ_T2、δ_C2代入到方程：

Δ＝δ_air+f_B-Aδ_C1-δ_T1-f_B-Aδ_air或Δ＝δ_air+f_B-Aδ_C2-δ_T2-f_B-Aδ_air，计算出待测微藻二氧化碳同化过程中的稳定碳同位素分馏值Δ；

第二，同时设置对照组和测试组，培养条件相同；对照组为培养液中不添加硅酸盐岩，培养待测微藻；测试组为培养液中添加硅酸盐岩，培养待测微藻；

第三，将对照组中的待测微藻在培养液中培养7天或7天以上，分别测定待测微藻的最初和最终生物量MC_O和MC_N；微藻的起始稳定碳同位素组成δ_O-C；培养后的微藻的稳定碳同位素组成δ¹³C的值δ_S-C；

第四，将测试组待测微藻在培养液中同样培养7天或7天以上，同样分别测定待测微藻的最初和最终生物量ME_O和ME_N；微藻的起始稳定碳同位素组成δ_O-E；培养后的微藻的稳定碳同位素组成δ¹³C的值δ_S-E；

第五，测定对照组和测试组每天培养环境中大气二氧化碳的日平均稳定碳同位素组成，由此计算出对照组和测试组的培养环境中大气二氧化碳的平均稳定碳同位素组成δ¹³C的值δ_air-E。

第六，分别计算对照组和测试组的微藻增值倍数P_C和P_T，P_C＝MC_N/MC_O，P_T＝ME_N/ME_O；

第七，将对照组的微藻的起始稳定碳同位素组成δ_O-C、培养后的微藻的稳定碳同位素组成δ¹³C的值δ_S-C以及微藻增值倍数P_C代入方程:计算出对照组中微藻的新生藻体稳定碳同位素组成δ¹³C值δ_T-C；

第八，同时，将测试组的微藻的起始稳定碳同位素组成δ_O-E、培养后的微藻的稳定碳同位素组成δ¹³C的值δ_S-E以及微藻增值倍数P_T代入方程：计算出测试组中微藻的新生藻体稳定碳同位素组成δ¹³C值δ_T-E；

第九，将δ_air-E、δ_T-C、Δ代入方程即可计算出测试组和对照组微藻的碳酸氢根离子途径份额f_b；

第十，确定测试组中硅酸盐岩风化过程中产生的碳酸氢根离子的碳同位素值δ_C-W；δ_C-W＝δ_air-E-1.1‰；

第十一，将δ_C-W、f_b、Δ代入方程δ_TB＝δ_C-W－Δ+9‰f_b计算出微藻仅利用硅酸盐岩风化产生的碳酸氢根离子所形成的藻体碳同位素值δ_TB；

第十二，将δ_T-E、δ_T-C、δ_TB代入方程即可算出微藻利用硅酸盐岩风化产生的无机碳源的利用份额f_B；

第十三，将微藻利用硅酸盐岩风化产生的无机碳源的利用份额f_B、测试组的微藻增值倍数P_T，代入方程D＝f_B(P_T-1),计算出待测微藻利用源于硅酸岩风化的无机碳能力D。

本发明的基本原理为：

稳定碳同位素的强烈分馏特征是识别不同无机碳来源的基础。自然界中碳元素有两种稳定同位素：¹²C和¹³C，它们的天然平均丰度分别为98.89％和1.11％。稳定碳同位素组成通常用δ¹³C(‰)表示，自然界中δ¹³C的变化为-90‰～+20‰(PDB标准)。稳定碳同位素的强烈分馏特征有利于识别不同无机碳来源。质量平衡原理以及同位素混合模型和化学计量学方法，是定量识别不同无机碳源的基础。

微藻能且只能利用CO₂与HCO₃^-作为自身的碳源，一方面CO₂作为线性非极性，呈电中性，能以分子形式经自由扩散进入细胞双层脂膜，进入细胞内为微藻光合作用所利用，此过程称为二氧化碳途径；另一方面碳酸氢根离子(HCO₃^-)通过碳酸酐酶(CA)催化，被微藻直接利用(直接利用：HCO₃^-通过载体蛋白或阴离子交换蛋白转运进入细胞，一部分在细胞内经胞内碳酸酐酶转化为CO₂，另一部分通过叶绿体膜蛋白主动转运到叶绿体内，经叶绿体碳酸酐酶转化为CO₂参与光合作用)或间接利用(间接利用：HCO₃^-在胞外碳酸酐酶的催化作用水解，形成CO₂直接自由扩散进行光合作用)，此过程为微藻的碳酸氢根离子利用途径。

(一)微藻二氧化碳同化过程中的稳定碳同位素分馏值(Δ)的计算原理

乙酰唑胺(acetazolamide，AZ)是含1,3,4-噻二唑环的杂环磺酰胺类碳酸酐酶胞外酶抑制剂，高浓度碳酸氢钠也对碳酸酐酶胞外酶有抑制作用，在高浓度碳酸氢钠和10mMAZ的作用下，依赖于胞外碳酸酐酶的碳酸氢根离子的间接转运途径和碳酸氢根离子的直接转运途径将同时被完全抑制。通过添加的两种标记的重碳酸盐的方法可以定量出微藻同化二氧化碳过程中的稳定碳同位素分馏值，进而建立计算无机碳利用途径份额的计算方法。

微藻可以利用大气的无机碳源和添加的无机碳源，且对于每一种来源的无机碳的利用，都存在CO₂和HCO₃^-两种无机碳利用途径，为此，我们可以建立如下两端元的同位素混合模型：

δ_Ti＝(1-f_Bi)δ_TA+f_Biδ_TB＝(1-f_Bi)[(1-f_bi)(δ_air-Δ)+f_bi(δ_air-Δ+9‰)]+f_Bi[(1-f_bi)(δ_Ci-Δ)+f_bi(δ_Ci-Δ+9‰)](i＝1,2)(1)

δ_TA是指微藻只利用大气来源的无机碳时的藻体碳同位素值，包括：CO₂和HCO₃^-这两种利用途径；δ_TB是指微藻只利用添加的碳酸氢钠来源的无机碳时的藻体碳同位素值，同样包括：CO₂和HCO₃^-这两种利用途径。δ_Ti为添加已知δ¹³C的某一碳酸氢钠培养的微藻藻体的δ¹³C值，f_Bi为微藻利用添加的无机碳占总碳源的份额，(1-f_Bi)为微藻利用大气二氧化碳占总碳源的份额。f_bi为微藻利用碳酸氢根离子途径，(1-f_bi)为微藻二氧化碳利用途径份额。δ_air、Δ以及δ_Ci分别为环境中大气二氧化碳的平均稳定碳同位素组成δ¹³C的值、二氧化碳同化过程中的稳定碳同位素分馏值以及同位素标记的碳酸氢钠的δ¹³C值。

在同一条件下培养的同一种微藻，对于添加的两种标记的重碳酸盐来说，方程(1)可以分别表示如下：

δ_T1＝(1-f_B1)[(1-f_b1)(δ_air-Δ)+f_b1(δ_air-Δ+9‰)]+f_B1[(1-f_b1)(δ_C1-Δ)+f_b1(δ_C1-Δ+9‰)](2)

δ_T2＝(1-f_B2)[(1-f_b2)(δ_air-Δ)+f_b2(δ_air-Δ+9‰)]+f_B2[(1-f_b2)(δ_C2-Δ)+f_b2(δ_C2-Δ+9‰)](3)

在方程(2)和(3)中，δ_T1为添加第一种已知δ¹³C的碳酸氢钠培养的微藻藻体的δ¹³C值；δ_T2为添加第二种已知δ¹³C的碳酸氢钠培养的微藻藻体的δ¹³C值；f_B1为微藻利用添加的第一种碳酸氢钠占总碳源的份额；f_B2为微藻利用添加的第二种碳酸氢钠占总碳源的份额；f_b1为培养在添加第一种已知δ¹³C的碳酸氢钠的培养液中的微藻利用碳酸氢根离子途径所占的份额；f_b2为培养在添加第二种已知δ¹³C的碳酸氢钠的培养液中的微藻利用碳酸氢根离子途径所占的份额。δ_C1为同位素标记1的碳酸氢钠的δ¹³C值，δ_C2为同位素标记2的碳酸氢钠的δ¹³C值；Δ值为微藻CO₂同化过程中的稳定碳同位素分馏值；δ_air值为培养微藻时房间空气中的碳同位素值。

不论添加哪种标记的碳酸氢钠，只要在同一培养条件下，同一种微藻利用添加的重碳酸盐所占总碳源的份额是相同的，即：f_B1＝f_B2＝f_B。不论添加哪种标记的碳酸氢钠，只要在同一培养条件下，同一种微藻利用碳酸氢根离子途径所占的份额是相同的，由此可以得到：f_b1＝f_b2＝f_b。此基础上，将方程(2)与方程(3)做差、化简可得：

在高浓度碳酸氢钠和10mMAZ的作用下，依赖于胞外碳酸酐酶的碳酸氢根离子的间接转运途径和碳酸氢根离子的直接转运途径将同时被完全抑制，所以：

f_b1＝f_b2＝f_b＝0(5)

这样，此时的(2)和(3)式可分别简化成：

Δ＝δ_air+f_Bδ_C1-δ_T1-f_Bδ_air(6)

Δ＝δ_air+f_Bδ_C2-δ_T2-f_Bδ_air(7)

δ_T1为添加第一种已知δ¹³C的碳酸氢钠培养的微藻藻体的δ¹³C值；δ_T2为添加第二种已知δ¹³C的碳酸氢钠培养的微藻藻体的δ¹³C值；

收获后的微藻藻体的δ¹³C值，不完全是新生成的藻体的δ¹³C，而是含有培养时的原始藻体加新生成的藻体的δ¹³C混合值，因此，为了获取新生成的藻体的δ¹³C值δ_T1或δ_T2；我们建立了如下方程：

δ_O+(P-1)×δ_Ti＝P×δ_Si(8)

P为微藻增值倍数，其值为实验处理后的微藻生物量(M_N)/初始接种时的微藻生物量(M_O)；

δ_O为初始接种的微藻藻体碳同位素组成(δ¹³C)；

δ_Si为实验处理后的微藻藻体碳同位素组成(δ¹³C)(i＝1，2)；

δ_Ti为实验处理新生成的藻体碳同位素组成(δ¹³C)(i＝1，2)。

(8)式变形后为：

因此知道了δ_Ti、δ_Ci以及培养期间培养环境中大气二氧化碳平均稳定碳同位素组成δ¹³C的值δ_air，由(6)和(7)式计算微藻二氧化碳同化过程中的稳定碳同位素分馏值Δ。

(二)微藻利用硅酸盐岩风化产生的无机碳源的利用份额的计算原理

在风化水体中微藻可以利用直接来源于大气的无机碳源和硅酸盐岩风化产生的无机碳源，建立微藻碳同位素的两端元混合模型为：

δ_T＝(1-f_B)δ_TA+f_Bδ_TB(10)

化简可得：

f_B＝(δ_T-δ_TA)/(δ_TB-δ_TA)(11)

式中：δ_TA为微藻只利用直接来源于空气二氧化碳的无机碳时的碳同位素值；

δ_TB为微藻只利用硅酸盐岩风化产生的无机碳源时的碳同位素值；

δ_T值为硅酸盐岩处理下的微藻藻体的碳同位素值；

f_B为微藻利用硅酸盐岩风化产生的无机碳源的利用份额；

1-f_B为微藻利用直接来源于空气二氧化碳的无机碳的的利用份额。

且对于每一种来源的无机碳的利用，都存在CO₂和HCO₃^-两种无机碳利用途径，在无机碳利用过程中，HCO₃^-的直接转运过程存在约10‰的碳同位素分馏，由CAex催化的HCO₃^-的间接转运过程只存在约1.1‰的碳同位素分馏，两者之间存在约9‰的稳定碳同位素分馏差异。故建立的两端元混合模型(10)式可化为：

δ_T＝(1-f_B)[(1-f_b)(δ_air-Δ)+f_b(δ_air-Δ+9‰)]+f_B[(1-f_b)(δ_C-Δ)+f_b(δ_C-Δ+9‰)](11)

化简得：

f_b＝[f_B(δ_air-δ_C)+(Δ+δ_T-δ_air)]/9‰(12)

式中：δ_C值为硅酸盐岩风化过程中产生的碳酸氢根离子的碳同位素值；

Δ值为微藻CO₂同化过程中的稳定碳同位素分馏值；

δ_air值为培养微藻时环境空气中的碳同位素值；

f_b值为微藻的碳酸氢根离子途径利用份额；

1-f_b值为微藻的二氧化碳途径利用份额。

根据(12)式，微藻在有硅酸盐岩和无硅酸盐岩处理的情况下，其碳酸氢根离子途径份额并不会改变，则f_b＝f_b1＝f_b2：

f_b1＝[f_B1(δ_air-δ_C)+(Δ+δ_T1-δ_air)]/9‰(13)

f_b2＝[f_B2(δ_air-δ_C)+(Δ+δ_T2-δ_air)]/9‰(14)

其中：f_B1值为微藻在没有硅酸盐岩处理下利用硅酸盐岩风化产生的无机碳源的利用份额；

f_B2值为微藻在有硅酸盐岩处理下利用硅酸盐岩风化产生的无机碳源的利用份额；

f_b1值为微藻在没有硅酸盐岩处理下的碳酸氢根离子途径利用份额；

f_b2值为微藻在有硅酸盐岩处理下的碳酸氢根离子途径利用份额；

δ_T1值为在没有硅酸盐岩处理下微藻藻体的碳同位素值；

δ_T2值为在有硅酸盐岩处理下微藻藻体的碳同位素值。

当不添加硅酸盐岩时f_B1值为0，计算出f_b如(15)式

f_b＝(Δ+δ_T1－δ_air)/9‰(15)

根据(11)式，可将δ_TB表示为(16)式：

δ_TB＝(1-f_b)(δ_C－Δ)+f_b(δ_C－Δ+9‰)(16)

(16)式化简得：

δ_TB＝δ_C－Δ+9‰f_b(17)

δ_C值为硅酸盐岩风化过程中产生的碳酸氢根离子的碳同位素值，其值可根据空气中的二氧化碳分馏值计算得到，如(18)式：

δ_C＝δ_air－1.1‰(18)

在25℃下，缓慢的平衡状态下约为8.5‰，但是由于藻类快速吸收利用，使得二氧化碳水解平衡一直向着碳酸氢根离子方向发展，因此，在有藻类快速利用碳酸氢根离子的情况下硅酸盐岩风化过程产生的碳同位素分馏值取为1.1‰。

(11)式中微藻利用硅酸盐岩风化产生的无机碳源的利用份额f_B＝(δ_T－δ_TA)/(δ_TB－δ_TA)中δ_T值为添加硅酸盐岩培养的藻体碳同位素值，δ_TA即为同条件下不加硅酸盐岩培养的藻体碳同位素值，δ_TB可根据(17)式计算出，即可算出微藻利用硅酸盐岩风化产生的无机碳源的利用份额f_B。最后，以微藻利用硅酸盐岩风化产生的无机碳源的利用份额f_B与微藻新增值倍数的乘积来表征微藻利用源于硅酸岩风化的无机碳能力。

本发明的有益效果：

1)本方法不仅能量化在微藻对硅酸盐的生物溶蚀作用过程中的碳酸氢根离子途径份额，而且还能量化出微藻利用硅酸盐岩风化产生的碳酸氢根离子的份额。

2)本方法可以定量出不同微藻利用源于不同的硅酸岩风化的无机碳能力，结果具有可比性。

3)本方法提供了一种微藻易被忽视的部分“碳汇”的获取方法，为寻找遗失的碳汇提供技术支撑。

4)本方法还可以通过定量出不同微藻利用源于不同的硅酸岩风化的无机碳能力来比较不同微藻对硅酸盐的生物溶蚀作用。

具体实施方式

本发明的实例：它包括以下步骤：

第一步骤，利用添加碳酸氢钠的浓度为16mM及以上、同时附加10mM乙酰唑胺也即AZ双向同位素标记培养待测微藻，获取待测微藻二氧化碳同化过程中的稳定碳同位素分馏值Δ；分步骤如下：

(1)首先选择两种δ¹³C值差值大于10‰的碳酸氢钠作为同位素标记1和同位素标记2分别添加到检测培养液中，同时在同样的培养条件下培养起始生物量相同的待测微藻；检测培养液为添加碳酸氢钠的浓度为16mM及以上、同时附加10mM乙酰唑胺也即AZ的藻类培养液。测定待测微藻的起始生物量M_O和微藻的起始稳定碳同位素组成δ¹³C的值δ_O。同位素标记1的碳酸氢钠的δ¹³C值为δ_C1，同位素标记2的碳酸氢钠的δ¹³C值为δ_C2；将待测微藻在添加同位素标记的检测培养液中培养收获后，测定待测微藻的最终生物量M_N和在两种同位素标记的检测培养液培养的被考察微藻的稳定碳同位素组成δ¹³C的值δ_S1和δ_S2。测定培养期间每天培养环境中大气二氧化碳的日平均稳定碳同位素组成，由此计算出在培养期间，培养环境中大气二氧化碳的平均稳定碳同位素组成δ¹³C的值δ_air。

(2)待测微藻增值倍数P＝M_N/M_O。将微藻的起始稳定碳同位素组成δ_O、待测微藻增值倍数P以及在两种同位素标记的检测培养液培养的被考察微藻的稳定碳同位素组成δ_S1或δ_S2代入方程：和计算出在两种同位素标记的检测培养液培养的被考察微藻的新生藻体稳定碳同位素组成δ¹³C的值δ_T1和δ_T2。将δ_T1、δ_T2、δ_C1和δ_C2代入到方程：计算出待测微藻利用添加的无机碳源的份额f_B-A。

(3)将f_B-A、δ_air以及δ_T1和δ_C1或δ_T2、δ_C2代入到方程：

Δ＝δ_air+f_B-Aδ_C1-δ_T1-f_B-Aδ_air或Δ＝δ_air+f_B-Aδ_C2-δ_T2-f_B-Aδ_air，计算出待测微藻二氧化碳同化过程中的稳定碳同位素分馏值Δ；

第二步骤，同时设置对照组和测试组，培养条件相同；对照组为培养液中不添加硅酸盐岩，培养待测微藻；测试组为培养液中添加硅酸盐岩，培养待测微藻；

第三步骤，将对照组中的待测微藻在培养液中培养7天后，分别测定待测微藻的最初和最终生物量MC_O和MC_N；微藻的起始稳定碳同位素组成δ_O-C；培养后的微藻的稳定碳同位素组成δ¹³C的值δ_S-C；

第四步骤，将测试组待测微藻在培养液中同样培养7天后，同样分别测定待测微藻的最初和最终生物量ME_O和ME_N；微藻的起始稳定碳同位素组成δ_O-E；培养后的微藻的稳定碳同位素组成δ¹³C的值δ_S-E；

第五步骤，测定对照组和测试组每天培养环境中大气二氧化碳的日平均稳定碳同位素组成，由此计算出对照组和测试组的培养环境中大气二氧化碳的平均稳定碳同位素组成δ¹³C的值δ_air-E。

第六步骤，分别计算对照组和测试组的微藻增值倍数P_C和P_T，P_C＝MC_N/MC_O，P_T＝ME_N/ME_O；

第七步骤，将对照组的微藻的起始稳定碳同位素组成δ_O-C、培养后的微藻的稳定碳同位素组成δ¹³C的值δ_S-C以及微藻增值倍数P_C代入方程:计算出对照组中微藻的新生藻体稳定碳同位素组成δ¹³C值δ_T-C；

第八步骤，同时，将测试组的微藻的起始稳定碳同位素组成δ_O-E、培养后的微藻的稳定碳同位素组成δ¹³C的值δ_S-E以及微藻增值倍数P_T代入方程：计算出测试组中微藻的新生藻体稳定碳同位素组成δ¹³C值δ_T-E；

第九步骤，将δ_air-E、δ_T-C、Δ代入方程即可计算出测试组和对照组微藻的碳酸氢根离子途径份额f_b；

第十步骤，依据δ_air-E确定测试组中硅酸盐岩风化过程中产生的碳酸氢根离子的碳同位素值δ_C-W，δ_C-W＝δ_air-E-1.1‰；

第十一步骤，将δ_C-W、f_b、Δ代入方程δ_TB＝δ_C-W－Δ+9‰f_b计算出微藻仅利用硅酸盐岩风化产生的碳酸氢根离子所形成的藻体碳同位素值δ_TB；

第十二步骤，将δ_T-E、δ_T-C、δ_TB代入方程即可算出微藻利用硅酸盐岩风化产生的无机碳源的利用份额f_B；

第十三步骤，将微藻利用硅酸盐岩风化产生的无机碳源的利用份额f_B、测试组的微藻增值倍数P_T，代入方程D＝f_B(P_T-1),计算出待测微藻利用源于待测硅酸岩风化的无机碳能力D。

实施例1：在不同营养下莱茵衣藻利用源于硅酸岩风化的无机碳能力

实施例1的实施效果如下：

实验1：莱茵衣藻在二氧化碳同化过程中的稳定碳同位素分馏值(Δ)的计算；

培养材料为：莱茵衣藻。培养基为在基本培养液(SE培养基)添加16mM碳酸氢钠和10mMAZ，基本培养条件为：光周期L/D：12h/12h；温度25℃；光照强度为100μmol·m^-2·s^-1，pH值8.3。添加的碳酸氢钠的δ¹³C分别为-17.4‰(PDB)(δ_C1)和-28.4‰(PDB)(δ_C2)，表1表示莱茵衣藻有关实验数据。

表1莱茵衣藻新生藻体稳定碳同位素组成δ¹³C的值

依据表1可以获得莱茵衣藻的二氧化碳同化过程中的稳定碳同位素分馏值Δ如表2。

表2莱茵衣藻的二氧化碳同化过程中的稳定碳同位素分馏值Δ(‰)

实验2：莱茵衣藻利用硅酸岩风化产生的无机碳源的利用份额以及莱茵衣藻对硅酸岩的溶蚀作用能力的计算；

培养材料为：莱茵衣藻。培养基为SE培养基与剔镁SE培养基(用等物质量的硫酸钠代替硫酸镁)，其pH值6.3。基本培养条件为：光周期L/D：12h/12h；温度25℃；光照强度为100μmol·m^-2·s^-1。SE培养基条件件下培养莱茵衣藻，一组为添加橄榄岩的，一组不添加橄榄岩；剔镁SE培养基条件下培养莱茵衣藻，一组为添加橄榄岩的，一组不添加橄榄岩；每一组处理均设置三个实验平行，培养7天后测定各个处理的藻体碳同位素值，测定空气的稳定碳同位素值δ_air，结果见表3。

表3莱茵衣藻的碳同位素值、空气稳定碳同位素值以及微藻增值倍数

注：C.R+L表示剔镁SE培养基、添加莱茵衣藻的处理组；

C.R+R+L表示剔镁SE培养基、添加莱茵衣藻与橄榄岩样品的处理组；

C.R+F表示SE培养基、添加莱茵衣藻的处理组；

C.R+R+F表示SE培养基、添加莱茵衣藻与橄榄岩样品的处理组。

根据表3的相关实验数据，根据本发明的计算方法，可计算出测试组和对照组微藻的碳酸氢根离子途径份额f_b；硅酸盐岩风化过程中产生的碳酸氢根离子的碳同位素值δ_C-W、微藻仅利用硅酸盐岩风化产生的碳酸氢根离子所形成的藻体碳同位素值δ_TB、微藻利用硅酸盐岩风化产生的无机碳源的利用份额f_B以及微藻利用源于硅酸岩风化无机碳能力D(如表4)

表4莱茵衣藻利用的橄榄岩风化过程中的碳酸氢根的份额和微藻利用源于硅酸岩风化无机碳能力

*f_b不能超过1，超过1为测定精度造成的，所以取1。

从表4中可以看出，莱茵衣藻在缺镁培养基与不缺镁培养基的两种培养条件下微藻利用源于硅酸岩风化的无机碳能力明显不同。在缺镁条件下，微藻能利用的橄榄岩风化过程产生的碳酸氢根离子份额较大，促进了溶蚀反应方程式向右进行，产生的碳酸氢根离子多，因此，微藻利用源于硅酸岩风化无机碳能力强，这与实际相符。上述测得的模式藻类利用的橄榄岩风化过程中的碳酸氢根的份额一方面可用于估算野外环境下微藻利用硅酸岩风化产生的无机碳形成的碳汇，另一方面还可以用于不同培养环境下无机碳利用途径的份额、无机碳源利用份额的确定。

实施例2：在不同营养下蛋白核小球藻利用源于硅酸岩风化的无机碳能力

实施例2的实施效果如下：

实验1：蛋白核小球藻二氧化碳同化过程中的稳定碳同位素分馏值(Δ)的计算；

培养材料为：蛋白核小球藻。培养基为在基本培养液(SE培养基)添加16mM碳酸氢钠和10mMAZ，基本培养条件为：光周期L/D：12h/12h；温度25℃；光照强度为100μmol·m^-2·s^-1，pH值8.3。添加的碳酸氢钠的δ¹³C分别为-17.4‰(PDB)(δ_C1)和-28.4‰(PDB)(δ_C2)，表5表示蛋白核小球藻的有关实验数据。

表5蛋白核小球藻新生藻体稳定碳同位素组成δ¹³C的值

依据表5可以获得蛋白核小球藻的二氧化碳同化过程中的稳定碳同位素分馏值Δ如表6。

表6蛋白核小球藻的二氧化碳同化过程中的稳定碳同位素分馏值Δ(‰)

实验2：蛋白核小球藻利用硅酸岩风化产生的无机碳源的利用份额以及蛋白核小球藻利用源于硅酸岩风化无机碳能力的计算；

培养材料为：蛋白核小球藻。培养基为SE培养基与剔镁SE培养基(用等物质量的硫酸钠代替硫酸镁)，其pH值6.3。基本培养条件为：光周期L/D：12h/12h；温度25℃；光照强度为100μmol·m^-2·s^-1。SE培养基条件件下培养莱茵衣藻，一组为添加橄榄岩的，一组不添加橄榄岩；剔镁SE培养基条件件下培养蛋白核小球藻，一组为添加橄榄岩的，一组不添加橄榄岩；每一组处理均设置三个实验平行，培养七天后测定各个处理的藻体碳同位素值，测定空气的稳定碳同位素值δ_air，其结果见表7。

表7蛋白核小球藻的碳同位素值、空气稳定碳同位素值以及微藻增值倍数

注：C.P+L表示剔镁SE培养基、添加蛋白核小球藻的处理组；

C.P+R+L表示剔镁SE培养基、添加小球藻与橄榄岩样品的处理组；

C.P+F表示SE培养基、蛋白核小球藻的处理组；

C.P+R+F表示SE培养基、蛋白核小球藻、橄榄岩样品的处理组；

根据表7的相关实验数据，根据本发明的计算方法，可计算出测试组和对照组微藻的碳酸氢根离子途径份额f_b；硅酸盐岩风化过程中产生的碳酸氢根离子的碳同位素值δ_C-W、微藻仅利用硅酸盐岩风化产生的碳酸氢根离子所形成的藻体碳同位素值δ_TB、微藻利用硅酸盐岩风化产生的无机碳源的利用份额f_B以及微藻利用源于硅酸岩风化无机碳能力D(如表8)

表8蛋白核小球藻利用的橄榄岩风化过程中的碳酸氢根的份额和蛋白核小球藻利用源于硅酸岩风化无机碳能力

从表8中可以看出，在缺镁培养基与不缺镁培养基的两种培养条件下蛋白核小球藻利用源于硅酸岩风化的无机碳能力差异也较大。同样地，在缺镁条件下，蛋白核小球藻能利用的橄榄岩风化过程产生的碳酸氢根离子份额也大于不缺镁培养基中蛋白核小球藻能利用的橄榄岩风化过程产生的碳酸氢根离子份额，最终导致蛋白核小球藻在缺镁培养基中利用源于硅酸岩风化的无机碳能力大于不缺镁培养基中的蛋白核小球藻利用源于硅酸岩风化的无机碳能力。

综合实施例1和实施例2的结果可以看出，蛋白核小球藻利用源于橄榄岩风化的无机碳能力远远大于莱茵衣藻利用源于橄榄岩风化的无机碳能力，这与蛋白核小球藻的细胞显著小于莱茵衣藻的细胞有关。具有同样生物量的蛋白核小球藻因细胞小，比表面积大，接触硅酸岩的面积大，因此具有较强的溶蚀能力，产生的碳酸氢根离子量大，最终导致微藻利用源于硅酸岩风化的无机碳能力强。