本发明涉及海洋微生物生态技术领域，具体的说是一种基于细胞碳生物量的海洋超微型浮游植物固碳量计算的模型方法。
背景技术：
海洋超微型浮游植物(0.2-2μm)是海洋中数量最多的初级生产者，在全球海洋中广泛分布，在某些寡营养大洋海区可以贡献初级生产力的80％以上。海洋超微型浮游植物包括原核(原绿球藻Pro、聚球藻Syn)和真核(真核超微型浮游植物Euk)类群。原绿球藻集中分布在南北纬40度之间寡营养的大洋海区。聚球藻和真核超微型浮游植物广布在全球海洋中，并在近岸、陆架海、高纬度开阔大洋海区对初级生产力有重要贡献。
目前对于超微型浮游植物丰度的研究报道较多，利用流式细胞仪可以精确测定样品中原核(原绿球藻Pro、聚球藻Syn)和真核(真核超微型浮游植物Euk)的丰度，但关于超微型浮游植物初级生产力的数据较少。原因之一是超微型浮游植物初级生产力的测定方法复杂、费时费力，且涉及放射性同位素的使用，对研究人员身体伤害较大。如果能够建立一种数学模型，利用数据较多的超微型浮游植物丰度的数据来估算超微型浮游生物初级生产力，将有望大大提高我们对超微型浮游植物固碳能力的认识。
1998年，Field等人在《Science》发表论文《Primaryproductionofthebiosphere:integratingterrestrialandoceaniccomponents》，针对整个浮游植物群落做了初级生产了模型的估算，而对于不同粒级浮游植物初级生产力估算的工作较少。Uitz等人分别在2006年、2008年和2010年，利用高效液相色谱得到的不同类群特征色素的组成进行分析，得到了对于不同粒级浮游植物初级生产力的估算模型。但是，因为高效液相色谱无法根据色素特征将真核超微型浮游植物与其他粒级真核浮游植物区分，该模型对于超微型浮游植物初级生产力的计算存在不足，低估了超微型浮游植物对初级生产力的贡献。
目前还没有一种方便有效的对于海洋超微型浮游植物固碳量计算的模型方法。
另一方面，由于聚球藻Syn一般都属于远海区域，很难采集到样本，因此，发明估算模型只针对于原绿球藻Pro和真核超微型浮游植物Euk。
技术实现要素：
根据上述不足之处，本发明提供一种海洋超微型浮游植物固碳量的估算方法。
为实现上述目的，本发明的技术方案在于：一种海洋超微型浮游植物固碳量的估算方法，其特征在于：包括如下步骤：
(1)利用流式细胞仪分别测定采集样品中聚球藻Syn的丰度ASyn和超微型真核藻类Euk的丰度AEuk；
(2)分别计算聚球藻Syn和超微型真核藻类Euk的碳生物量，
CSyn＝ASyn×CFSyn×106/1012，
其中，CSyn为聚球藻碳生物量，单位：毫克·C/立方米；ASyn为聚球藻丰度，单位：细胞/毫升；CFSyn为聚球藻平均单细胞碳量，为255飞克·C/细胞；1立方米＝106毫升；1毫克＝1012飞克；
CEuk＝AEuk×CFEuk×106/1012，
其中，CEuk为超微型真核藻类碳生物量，单位：毫克·C/立方米；AEuk为超微型真核藻类丰度，单位：细胞/毫升；CFEuk为超微型真核藻类平均单细胞碳量，为2590飞克·C/细胞；1立方米＝106毫升；1毫克＝1012飞克；
(3)分别计算聚球藻Syn和超微型真核藻类Euk的生长速率，
LnμSyn＝-0.55/(kT)+21.5，其中μSyn为聚球藻生长速率，单位：/天；Ln为自然对数；k为玻尔兹曼常数，8.62×10-5eV·K-1，1eV＝96.49kJ·mol-1，T为样品采样点的绝对温度，摄氏温度+273.15；
LnμEuk＝-0.57/(kT)+22.4，其中μEuk为超微型真核藻类生长速率，单位：/天，Ln为自然对数；k为玻尔兹曼常数，8.62×10-5eV·K-1，1eV＝96.49kJ·mol-1，T为样品采样点的绝对温度，摄氏温度+273.15；
(4)分别得到聚球藻Syn采样点的真光层深度ZeuSyn和超微型真核藻类Euk采样点的真光层深度ZeuEuk；
(5)分别计算聚球藻Syn和超微型真核藻类Euk固碳量基于水深的变化对表层光强的反应，
h(I0)Syn＝0.66125I0Syn/(I0Syn+4.1)，其中I0Syn是聚球藻Syn根据表层表面光合作用有效辐射，I0Syn，单位molesphotonsm-2h-1；
h(I0)Euk＝0.66125I0Euk/(I0Euk+4.1)，其中I0Syn是聚球藻Syn根据表层表面光合作用有效辐射，I0Syn，单位molesphotonsm-2h-1；
(6)分别计算聚球藻Syn和超微型真核藻类Euk的初级生产力NPP，
NPPSyn＝CSyn×μSyn×ZeuSyn×h(I0)Syn；
NPPEuk＝CEuk×μEuk×ZeuEuk×h(I0)Euk；
(7)将聚球藻Syn的初级生产力NPPSyn和超微型真核藻类Euk的初级生产力NPPEuk相加，得到采样点超微型浮游植物的初级生产力NPPPico，NPPPico＝NPPSyn+NPPEuk。
优选的是：所述步骤(3)中样品采集点的原位温度是实测数据或者根据采样点的时间和经纬度从NASA卫星对应的遥感数据月平均值获得。
优选的是：所述步骤(3)中聚球藻Syn和超微型真核藻类Euk的生长速率的计算方法是当叶绿素浓度超过0.5毫克/立方米时，基于稀释试验得到的原位聚球藻和超微型真核藻类生长速率与绝对温度之间具有良好的线性关系，利用这两个线性关系和原位温度计算聚球藻Syn和超微型真核藻类Euk的生长速率，所述的原位聚球藻的生长速率与绝对温度之间的线性关系为y＝-0.55x+21.5，r2＝0.47，n＝123，所述的超微型真核藻类的生长速率与绝对温度之间的线性关系为y＝-0.57x+22.4，r2＝0.41，n＝113。
优选的是：所述步骤(4)中聚球藻Syn采样点的真光层深度ZeuSyn和超微型真核藻类Euk采样点的真光层深度ZeuEuk是实测数据或利用采样点对应的漫衰减系数k490计算真光层深度。
优选的是：所述利用采样点对应的漫衰减系数k490计算真光层深度的公式为：Zeu＝Ln(0.01)/k490。
优选的是：当应用于不同海域的真光层计算时，计算真光层深度的公式中加入统计系数m，Zeu＝Ln(0.01)/(m×k490)。
本发明的有益效果在于：通过本方法可以计算近海、陆架海、高纬度开阔大洋(纬度高于40度)超微型浮游植物的固碳量(初级生产力)，与基于14C原位测定超微型浮游植物初级生产力相比，本方法省时省力，并且可以在较大的研究区域估算超微型浮游植物的初级生产力。而且经验证，利用该方法计算的初级生产力与基于14C测定的超微型浮游植物初级生产力之间具有较好的线性关系，估算可靠性强，为后期实际研究提供理论基础。
附图说明
图1是本发明的操作流程图。
图2是本发明聚球藻的生长速率与绝对温度之间的线性关系图。
图3是本发明超微型真核藻类的生长速率与绝对温度之间的线性关系图。
图4是利用本发明模型估算的超微型浮游植物初级生产力与基于14C测定的超微型浮游植物初级生产力之间的关系图，采样地点为：南大洋；北大西洋；南海。
具体实施例
下面结合具体实施例对本发明做进一步说明。
模型的建立
一种海洋超微型浮游植物固碳量的估算方法，其特征在于：包括如下步骤：
(1)利用流式细胞仪分别测定采集样品中聚球藻Syn的丰度ASyn和超微型真核藻类Euk的丰度AEuk；
(2)分别计算聚球藻Syn和超微型真核藻类Euk的碳生物量，
CSyn＝ASyn×CFSyn×106/1012，
其中，CSyn为聚球藻碳生物量，单位：毫克·C/立方米；ASyn为聚球藻丰度，单位：细胞/毫升；CFSyn为聚球藻平均单细胞碳量，为255飞克·C/细胞；1立方米＝106毫升；1毫克＝1012飞克；
CEuk＝AEuk×CFEuk×106/1012，
其中，CEuk为超微型真核藻类碳生物量，单位：毫克·C/立方米；AEuk为超微型真核藻类丰度，单位：细胞/毫升；CFEuk为超微型真核藻类平均单细胞碳量，为2590飞克·C/细胞；1立方米＝106毫升；1毫克＝1012飞克；
(3)分别计算聚球藻Syn和超微型真核藻类Euk的生长速率，
LnμSyn＝-0.55/(kT)+21.5，其中μSyn为聚球藻生长速率，单位：/天；Ln为自然对数；k为玻尔兹曼常数，8.62×10-5eV·K-1，1eV＝96.49kJ·mol-1，T为样品采样点的绝对温度，摄氏温度+273.15；
LnμEuk＝-0.57/(kT)+22.4，其中μEuk为超微型真核藻类生长速率，单位：/天，Ln为自然对数；k为玻尔兹曼常数，8.62×10-5eV·K-1，1eV＝96.49kJ·mol-1，T为样品采样点的绝对温度，摄氏温度+273.15；
聚球藻Syn和超微型真核藻类Euk的生长速率的计算方法是当叶绿素浓度超过0.5毫克/立方米时，基于稀释试验得到的原位聚球藻和超微型真核藻类生长速率与绝对温度之间具有良好的线性关系，利用这两个线性关系和原位温度计算聚球藻Syn和超微型真核藻类Euk的生长速率，所述的原位聚球藻的生长速率与绝对温度之间的线性关系为y＝-0.55x+21.5，r2＝0.47，n＝123，所述的超微型真核藻类的生长速率与绝对温度之间的线性关系为y＝-0.57x+22.4，r2＝0.41，n＝113。
(4)分别得到聚球藻Syn采样点的真光层深度ZeuSyn和超微型真核藻类Euk采样点的真光层深度ZeuEuk；聚球藻Syn采样点的真光层深度ZeuSyn和超微型真核藻类Euk采样点的真光层深度ZeuEuk是实测数据或利用采样点对应的漫衰减系数k490计算真光层深度。
其中，样品采集点的原位温度是实测数据或者根据采样点的时间和经纬度从NASA卫星对应的遥感数据月平均值获得。
利用采样点对应的漫衰减系数k490计算真光层深度的公式为：Zeu＝Ln(0.01)/k490。当应用于不同海域的真光层计算时，计算真光层深度的公式中加入统计系数m，Zeu＝Ln(0.01)/(m×k490)。
(5)分别计算聚球藻Syn和超微型真核藻类Euk固碳量基于水深的变化对表层光强的反应，
h(I0)Syn＝0.66125I0Syn/(I0Syn+4.1)，其中I0Syn是聚球藻Syn根据表层表面光合作用有效辐射，I0Syn，单位molesphotonsm-2h-1；
h(I0)Euk＝0.66125I0Euk/(I0Euk+4.1)，其中I0Syn是聚球藻Syn根据表层表面光合作用有效辐射，I0Syn，单位molesphotonsm-2h-1；
(6)分别计算聚球藻Syn和超微型真核藻类Euk的初级生产力NPP，
NPPSyn＝CSyn×μSyn×ZeuSyn×h(I0)Syn；
NPPEuk＝CEuk×μEuk×ZeuEuk×h(I0)Euk；
(7)将聚球藻Syn的初级生产力NPPSyn和超微型真核藻类Euk的初级生产力NPPEuk相加，得到采样点超微型浮游植物的初级生产力NPPPico，NPPPico＝NPPSyn+NPPEuk。
实施例1
如图2所示取值，即ASyn＝1922，AEuk＝278，T＝4.9839473℃，I0＝40.77893，Chl＝4.1153时(T、I0、Chl根据样品采集时间和地点：2005年3月，东经121.4833°，北纬40.45°，从NASAAqua/MODIS卫星的2005年3月平均值得到)，模型计算得到聚球藻的初级生产力为0.81毫克·C/立方米，超微型真核藻类的初级生产力为1.27毫克·C/立方米，超微型浮游植物初级生产力为2.08毫克·C/立方米。
实施例2：
取值为ASyn＝5420，AEuk＝1472，T＝16.34129℃，I0＝49.78449，Chl＝3.8289时(T、I0、Chl根据样品采集时间和地点：2005年6月，东经121.03°，北纬38.70°，从NASAAqua/MODIS卫星的2005年6月平均值得到)，模型计算得到聚球藻的初级生产力为5.4毫克·C/立方米，超微型真核藻类的初级生产力为16.6毫克·C/立方米，超微型浮游植物初级生产力为22.0毫克·C/立方米。
实施例3：
取值为ASyn＝24593，AEuk＝2165，T＝22.85763℃，I0＝33.25993，Chl＝6.72709时(T、I0、Chl根据样品采集时间和地点：2005年9月，东经121.07°，北纬40.73°，从NASAAqua/MODIS卫星的2005年9月平均值得到)，模型计算得到聚球藻的初级生产力为24.2毫克·C/立方米，超微型真核藻类的初级生产力为24.3毫克·C/立方米，超微型浮游植物初级生产力为48.5毫克·C/立方米。
实施例4：
取值为ASyn＝10566，AEuk＝3774，T＝4.824732℃，I0＝16.64843，Chl＝4.72373时(T、I0、Chl根据样品采集时间和地点：2005年12月，东经119.64°，北纬39.90°，从NASAAqua/MODIS卫星的2005年12月平均值得到)，模型计算得到聚球藻的初级生产力为4.3毫克·C/立方米，超微型真核藻类的初级生产力为16.7毫克·C/立方米，超微型浮游植物初级生产力为21.0毫克·C/立方米。