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

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

背景技术：

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

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

技术实现要素：

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

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

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

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

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

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

分步四、获取待测微藻增值倍数p＝mn/mo；

分步五、将待测微藻的起始稳定碳同位素组成δo、待测微藻增值倍数p以及在两种同位素标记的检测培养液培养的被考察微藻的稳定碳同位素组成δs1或δs2代入方程：和计算出在两种同位素标记的检测培养液培养的被考察微藻的新生藻体稳定碳同位素组成δ¹³c的值δt1和δt2；

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

分步七、将fb-a、δair以及δt1和δc1或δt2、δc2代入到方程：

δ＝δair+fb-aδc1-δt1-fb-aδair或δ＝δair+fb-aδc2-δt2-fb-aδair，计算出待测微藻二氧化碳同化过程中的稳定碳同位素分馏值δ；

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

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

第四，将测试组待测微藻在培养液中同样培养7天或7天以上，同样分别测定待测微藻的最初和最终生物量meo和men；微藻的起始稳定碳同位素组成δo-e；培养后的微藻的稳定碳同位素组成δ¹³c的值δs-e；

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

第六，分别计算对照组和测试组的微藻增值倍数pc和pt，pc＝mcn/mco，pt＝men/meo；

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

第八，同时，将测试组的微藻的起始稳定碳同位素组成δo-e、培养后的微藻的稳定碳同位素组成δ¹³c的值δs-e以及微藻增值倍数pt代入方程：计算出测试组中微藻的新生藻体稳定碳同位素组成δ¹³c值δt-e；

第九，将δair-e、δt-c、δ代入方程即可计算出测试组和对照组微藻的碳酸氢根离子途径份额fb；

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

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

第十二，将δt-e、δt-c、δtb代入方程即可算出微藻利用硅酸盐岩风化产生的无机碳源的利用份额fb；

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

本发明的基本原理为：

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

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

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

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

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

δti＝(1-fbi)δta+fbiδtb＝(1-fbi)[(1-fbi)(δair-δ)+fbi(δair-δ+9‰)]+fbi[(1-fbi)(δci-δ)+fbi(δci-δ+9‰)](i＝1,2)(1)

δta是指微藻只利用大气来源的无机碳时的藻体碳同位素值，包括：co2和hco3^-这两种利用途径；δtb是指微藻只利用添加的碳酸氢钠来源的无机碳时的藻体碳同位素值，同样包括：co2和hco3^-这两种利用途径。δti为添加已知δ¹³c的某一碳酸氢钠培养的微藻藻体的δ¹³c值，fbi为微藻利用添加的无机碳占总碳源的份额，(1-fbi)为微藻利用大气二氧化碳占总碳源的份额。fbi为微藻利用碳酸氢根离子途径，(1-fbi)为微藻二氧化碳利用途径份额。δair、δ以及δci分别为环境中大气二氧化碳的平均稳定碳同位素组成δ¹³c的值、二氧化碳同化过程中的稳定碳同位素分馏值以及同位素标记的碳酸氢钠的δ¹³c值。

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

δt1＝(1-fb1)[(1-fb1)(δair-δ)+fb1(δair-δ+9‰)]+fb1[(1-fb1)(δc1-δ)+fb1(δc1-δ+9‰)](2)

δt2＝(1-fb2)[(1-fb2)(δair-δ)+fb2(δair-δ+9‰)]+fb2[(1-fb2)(δc2-δ)+fb2(δc2-δ+9‰)](3)

在方程(2)和(3)中，δt1为添加第一种已知δ¹³c的碳酸氢钠培养的微藻藻体的δ¹³c值；δt2为添加第二种已知δ¹³c的碳酸氢钠培养的微藻藻体的δ¹³c值；fb1为微藻利用添加的第一种碳酸氢钠占总碳源的份额；fb2为微藻利用添加的第二种碳酸氢钠占总碳源的份额；fb1为培养在添加第一种已知δ¹³c的碳酸氢钠的培养液中的微藻利用碳酸氢根离子途径所占的份额；fb2为培养在添加第二种已知δ¹³c的碳酸氢钠的培养液中的微藻利用碳酸氢根离子途径所占的份额。δc1为同位素标记1的碳酸氢钠的δ¹³c值，δc2为同位素标记2的碳酸氢钠的δ¹³c值；δ值为微藻co2同化过程中的稳定碳同位素分馏值；δair值为培养微藻时房间空气中的碳同位素值。

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

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

fb1＝fb2＝fb＝0(5)

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

δ＝δair+fbδc1-δt1-fbδair(6)

δ＝δair+fbδc2-δt2-fbδair(7)

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

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

δo+(p-1)×δti＝p×δsi(8)

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

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

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

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

(8)式变形后为：

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

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

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

δt＝(1-fb)δta+fbδtb(10)

化简可得：

fb＝(δt-δta)/(δtb-δta)(11)

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

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

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

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

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

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

δt＝(1-fb)[(1-fb)(δair-δ)+fb(δair-δ+9‰)]+fb[(1-fb)(δc-δ)+fb(δc-δ+9‰)](11)

化简得：

fb＝[fb(δair-δc)+(δ+δt-δair)]/9‰(12)

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

δ值为微藻co2同化过程中的稳定碳同位素分馏值；

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

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

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

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

fb1＝[fb1(δair-δc)+(δ+δt1-δair)]/9‰(13)

fb2＝[fb2(δair-δc)+(δ+δt2-δair)]/9‰(14)

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

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

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

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

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

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

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

fb＝(δ+δt1－δair)/9‰(15)

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

δtb＝(1-fb)(δc－δ)+fb(δc－δ+9‰)(16)

(16)式化简得：

δtb＝δc－δ+9‰fb(17)

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

δc＝δair－1.1‰(18)

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

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

本发明的有益效果：

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

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

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

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

具体实施方式

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

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

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

(2)待测微藻增值倍数p＝mn/mo。将微藻的起始稳定碳同位素组成δo、待测微藻增值倍数p以及在两种同位素标记的检测培养液培养的被考察微藻的稳定碳同位素组成δs1或δs2代入方程：和计算出在两种同位素标记的检测培养液培养的被考察微藻的新生藻体稳定碳同位素组成δ¹³c的值δt1和δt2。将δt1、δt2、δc1和δc2代入到方程：计算出待测微藻利用添加的无机碳源的份额fb-a。

(3)将fb-a、δair以及δt1和δc1或δt2、δc2代入到方程：

δ＝δair+fb-aδc1-δt1-fb-aδair或δ＝δair+fb-aδc2-δt2-fb-aδair，计算出待测微藻二氧化碳同化过程中的稳定碳同位素分馏值δ；

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

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

第四步骤，将测试组待测微藻在培养液中同样培养7天后，同样分别测定待测微藻的最初和最终生物量meo和men；微藻的起始稳定碳同位素组成δo-e；培养后的微藻的稳定碳同位素组成δ¹³c的值δs-e；

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

第六步骤，分别计算对照组和测试组的微藻增值倍数pc和pt，pc＝mcn/mco，pt＝men/meo；

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

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

第九步骤，将δair-e、δt-c、δ代入方程即可计算出测试组和对照组微藻的碳酸氢根离子途径份额fb；

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

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

第十二步骤，将δt-e、δt-c、δtb代入方程即可算出微藻利用硅酸盐岩风化产生的无机碳源的利用份额fb；

第十三步骤，将微藻利用硅酸盐岩风化产生的无机碳源的利用份额fb、测试组的微藻增值倍数pt，代入方程d＝fb(pt-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的相关实验数据，根据本发明的计算方法，可计算出测试组和对照组微藻的碳酸氢根离子途径份额fb；硅酸盐岩风化过程中产生的碳酸氢根离子的碳同位素值δc-w、微藻仅利用硅酸盐岩风化产生的碳酸氢根离子所形成的藻体碳同位素值δtb、微藻利用硅酸盐岩风化产生的无机碳源的利用份额fb以及微藻利用源于硅酸岩风化无机碳能力d(如表4)

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

*fb不能超过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的相关实验数据，根据本发明的计算方法，可计算出测试组和对照组微藻的碳酸氢根离子途径份额fb；硅酸盐岩风化过程中产生的碳酸氢根离子的碳同位素值δc-w、微藻仅利用硅酸盐岩风化产生的碳酸氢根离子所形成的藻体碳同位素值δtb、微藻利用硅酸盐岩风化产生的无机碳源的利用份额fb以及微藻利用源于硅酸岩风化无机碳能力d(如表8)

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

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

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