本发明专利技术公开了一种基于水量的湖泊沉积物埋藏有机碳的估算方法。步骤包括采集湖泊沉积物样本；测试分析沉积物年代、沉积速率、沉积物干密度以及有机碳含量；计算湖泊的有机碳埋藏速率；利用卫星高度计数据提取湖泊正高水位；剔除异常值，获取单天水位，并计算年均水位；获取年均水位的高值、低值和均值；获取水位极值年份的秋季无云Landsat数据，根据MNDWI直方图设定阈值提取湖泊的水体面积；假定湖泊为圆锥体，根据水位极高年份和极低年份的水体面积确定湖泊形状，估算多年平均水位对应的湖泊水量；根据水量换算等体积圆锥体面积，结合沉积物年代和有机碳埋藏速率，获取基于水量的湖泊沉积物埋藏有机碳总量。沉积物埋藏有机碳总量。沉积物埋藏有机碳总量。
【技术实现步骤摘要】
一种基于水量的湖泊沉积物埋藏有机碳估算方法
[0001]本专利技术涉及一种基于水量的湖泊沉积物埋藏有机碳估算方法，属于湖泊水文领域。
技术背景
[0002]湖泊是重要的碳汇。早期由于湖泊水体占地球陆地表面积较小(湖泊约占陆地面积的2.8％)，而忽视了湖泊在碳循环中的作用。尽管从湖泊水库等内陆水体排放到大气中的二氧化碳和甲烷会成为温室气体，是碳源，但是保存在沉积物中的碳会进入长期地质循环，也是永久的、重要的碳汇。局地和区域尺度的研究表明，即使湖泊和水库沉积物中的碳埋藏通量通常小于碳排放通量，但仍是一种重要的长期碳汇，在区域乃至全球碳循环中发挥着重要作用。
[0003]此外，与海洋相比，湖泊水库等内陆水体埋藏碳的效率更高，这是因为内陆水体的沉积速率高、供氧量低，陆源有机碳比例高，所以有更高比例的沉积埋藏有机碳未被矿化，而是保留在沉积物中。最近也有研究表明自人类世以来，由于水土流失、拦河筑坝和富营养化，内陆水体的有机碳埋藏速率正在增加。
[0004]湖泊面积及流域面积是影响湖泊碳埋藏的重要因素，湖泊面积与湖泊流域面积比越大，单位面积接收的陆源有机质越多，湖泊碳通量越高。湖泊水深也是湖泊埋藏碳的重要影响因素。但限于数据不足，目前研究大多还是直接将有机碳埋藏速率直接与湖泊面积相乘获得有机碳埋藏总量。而湖泊面积通常是某一年的瞬时结果，这增加了估算结果的不确定性。因此，有必要结合湖泊实际水量估算埋藏有机碳总量。
技术实现思路
[0005]本专利技术要解决的问题是：克服现有技术的上述缺陷，提供一种基于水量的湖泊沉积物埋藏有机碳估算方法，该方法采集湖泊沉积物样本，获取湖泊有机碳埋藏速率；利用遥感手段获取湖泊水位和面积，并估算湖泊水量。从而获取基于湖泊水量的湖泊埋藏有机碳总量。
[0006]为了解决上述技术问题，本专利技术提出的技术方案是：基于水量的湖泊沉积物埋藏有机碳估算方法，包括以下步骤：
[0007]第一步、采集湖泊沉积物样本：每个湖泊采集3
‑
10个湖底柱状沉积物样品，并记录采样点经纬度坐标、海拔、湖泊深度和柱状沉积物长度；
[0008]第二步、测试分析得到沉积物年代t，沉积物干密度DBD，以及有机碳TOC 含量；
[0009]第三步、根据柱状沉积物的长度L和沉积物年代t计算沉积速率SR，计算公式如下：
[0010]SR＝L/t
[0011]结合有机碳含量TOC及沉积物干密度DBD计算每个湖泊岩芯近百年来有机碳埋藏速率OCBR，计算公式如下：
[0012]OCBR＝TOC
×
SR
×
DBD
[0013]计算采样湖泊多个岩芯样品有机碳埋藏速率的均值，代表这个湖泊的年平均有机碳埋藏速率；
[0014]第四步、根据雷达高度计和激光高度计数据，计算湖泊正高水位，获得湖泊的年均水位最高值h
max
和年均水位最低值h
min
的及其对应年份，并据此计算多年平均水位h
ave
；
[0015]第五步、获取水位极高年份和水位极低年份的秋季无云Landsat影像，根据改进的归一化水体指数MNDW直方图设定阈值提取湖泊的水体面积，分别得到湖泊最大面积S
max
和湖泊最小面积S
min
；
[0016]第六步、假定湖泊为圆锥体，根据最高水位h
max
、最低水位h
min、
以及最高水位、最低水位对应的湖泊最大面积S
max
和湖泊最小面积S
min
，估算多年平均水位对应的湖泊水量，并根据水量换算等体积圆锥体表面积S，以此近似湖泊表面积，等体积圆锥体表面积S的计算公式如下：
[0017][0018]第七步、结合沉积物年代t、有机碳埋藏速率OCBR和换算后的表面积S，获取基于水量的湖泊埋藏有机碳总量TOC
lake
，公式如下：
[0019]TOC
lake
＝OCBR
×
S
×
t。
[0020]综上，本专利技术基于水量估算湖泊沉积物埋藏有机碳的总量，流程清晰、结果可靠。目前湖泊沉积物埋藏有机碳的估算都采用埋藏速率乘以沉积年代和某一年的水体面积获得，但湖泊的深度、形状也会影响湖泊有机碳的埋藏，基于水量估算埋藏有机碳更严谨合理。且本方法中的水量估算完全依靠遥感数据，减少野外工作量，方便快捷。
附图说明
[0021]下面结合附图对本专利技术作进一步的说明。
[0022]图1为湖泊埋藏有机碳估算方法的流程图。
[0023]图2为利用卫星高度计提取的湖泊年均水位。
[0024]图3为水位极高年份的Landsat MNDWI结果以及提取的水体范围。
[0025]图4为水位极低年份的Landsat MNDWI结果以及提取的水体范围。
具体实施方式
[0026]下面以太湖为例，根据附图详细阐述本专利技术，使本专利技术的技术路线和操作步骤更加清晰。
[0027]如图1为估算湖泊埋藏有机碳的方法流程图，具体步骤如下：
[0028]第一步、采集湖泊沉积物样本：每个湖泊采集3
‑
10个湖底柱状沉积物样品，并记录采样点经纬度坐标、海拔、湖泊深度和柱状沉积物长度等详细参数；同时采集湖泊所在流域内的典型植被、土壤及水的样品。根据需要，在室内实验室将所有沉积物样品按照0.5或1.0cm等间距分样，并放置于密封塑料袋内冷藏保存。
[0029]第二步、测试分析沉积物年代t和理化指标。将沉积物样品低温烘干后称重，采用γ分析方法测试天然放射性核素
210
Pb及人工放射性核素
137
Cs。用于年代分析的过剩
210
Pb
(210Pb
ex
)为
210
Pb
tot
与
226
Ra比活度的差值。以662kev处的特征峰值作为
137
Cs的比活度。利用
137
Cs比活度随钻孔深度变化的峰值确定其特征年代，并与210Pb比活度随钻孔深度变化进行对比验证，并最终确定沉积物样品的年代。理化指标测试分析包括沉积物含水量、干容重与烧失量测试分析，总碳TC和总氮TN以及有机碳TOC、无机碳TIC等元素测试分析等。得到沉积物干密度DBD，以及有机碳TOC含量。
[0030]第三步、湖泊沉积物有机碳埋藏速率估算：根据柱状沉积物的长度L和沉积物年代t计算沉积速率SR，计算公式如下：
[0031]SR＝L/t
[0032]结合TOC含量及干密度DBD计算每个湖泊岩芯近百年来有机碳埋藏速率 OCBR，计算公式如下：
[0033]OCBR＝TOC
×
SR
×
DBD
[0034]然后计算采样湖泊多个岩芯样品有机碳埋藏速率的均值，代本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种基于水量的湖泊沉积物埋藏有机碳估算方法，包括以下步骤：第一步、采集湖泊沉积物样本：每个湖泊采集3
–
10个湖底柱状沉积物样品，并记录采样点经纬度坐标、海拔、湖泊深度和柱状沉积物长度；第二步、测试分析得到沉积物年代t，沉积物干密度DBD，以及有机碳TOC含量；第三步、根据柱状沉积物的长度L和沉积物年代t计算沉积速率SR，计算公式如下：SR＝L/t结合有机碳含量TOC及沉积物干密度DBD计算每个湖泊岩芯近百年来有机碳埋藏速率OCBR，计算公式如下：OCBR＝TOC
×
SR
×
DBD计算采样湖泊多个岩芯样品有机碳埋藏速率的均值，代表这个湖泊的年平均有机碳埋藏速率；第四步、根据雷达高度计和激光高度计数据，计算湖泊正高水位，获得湖泊的年均水位最高值h
max
和年均水位最低值h
min
的及其对应年份，并据此计算多年平均水位h
ave
；第五步、获取水位极高年份和水位极低年份的秋季无云Landsat影像，根据改进的归一化水体指数MNDW直方图设定阈值提取湖泊的水体面积，分别得到湖泊最大面积S
max
和湖泊最小面积S
min
；第六步、假定湖泊为圆锥体，根据最高水位h
max
、最低水位h
min
、以及最高水位、最低水位对应的湖泊最大面积S
max
和湖泊最小面积S
min
，估算多年平均水位对应的湖泊水量，并根据水量换算等体积圆锥体表面积S，以此近似湖泊表面积，等体积圆锥体表面积S的计算公式如下：第七步、结合沉积物年代t、有机碳埋藏速率OCBR和换算后的表面积S，获取基于水量的湖泊埋藏有机碳总量TOC
lake
，公式如下：TOC
lake
＝OCBR
×
S
×
t。2.根据权利要求1所述估算湖泊埋藏有机碳的方法，其特征在于：第二步中，将沉积物样品低温烘干后称重，采用γ分析方法测试天然放射性核素
210
Pb及人工放射性核素
137
Cs；用于年代...
【专利技术属性】
技术研发人员：柯长青，蔡宇，
申请(专利权)人：南京大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。