留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

南海神狐海域水合物发育区浅表层沉积物甲烷周转定量模拟

胡廷苍 张艳平 胡钰 罗敏 陈多福

胡廷苍, 张艳平, 胡钰, 罗敏, 陈多福. 南海神狐海域水合物发育区浅表层沉积物甲烷周转定量模拟[J]. 海洋地质与第四纪地质, 2020, 40(3): 99-108. doi: 10.16562/j.cnki.0256-1492.2019042401
引用本文: 胡廷苍, 张艳平, 胡钰, 罗敏, 陈多福. 南海神狐海域水合物发育区浅表层沉积物甲烷周转定量模拟[J]. 海洋地质与第四纪地质, 2020, 40(3): 99-108. doi: 10.16562/j.cnki.0256-1492.2019042401
HU Tingcang, ZHANG Yanping, HU Yu, LUO Min, CHEN Duofu. Quantitative assessment of methane turnover in shallow surface sediments of hydrate-bearing areas in Shenhu area of South China Sea[J]. Marine Geology & Quaternary Geology, 2020, 40(3): 99-108. doi: 10.16562/j.cnki.0256-1492.2019042401
Citation: HU Tingcang, ZHANG Yanping, HU Yu, LUO Min, CHEN Duofu. Quantitative assessment of methane turnover in shallow surface sediments of hydrate-bearing areas in Shenhu area of South China Sea[J]. Marine Geology & Quaternary Geology, 2020, 40(3): 99-108. doi: 10.16562/j.cnki.0256-1492.2019042401

南海神狐海域水合物发育区浅表层沉积物甲烷周转定量模拟


doi: 10.16562/j.cnki.0256-1492.2019042401
详细信息
    作者简介:

    胡廷苍(1993―),男,硕士研究生,主要从事海洋地质研究,E-mail:m170300611@st.shou.edu.cn

    通讯作者: 陈多福(1962―),男,教授,博士生导师,主要从事海洋地质研究,E-mail:dfchen@shou.edu.cn
  • 基金项目:  国家重点研究发展计划冷泉重点项目“中国海域冷泉系统演变过程及其机制”(2018YFC0310003);国家自然科学基金项目“南海北部冷泉和天然气水合物发育区海底浅表层沉积物碳循环数值模拟”(41730528),“冲绳海槽海底冷泉—热液系统相互作用及资源效应”(91858208)
  • 中图分类号: P736.4, P738

Quantitative assessment of methane turnover in shallow surface sediments of hydrate-bearing areas in Shenhu area of South China Sea

More Information
  • 摘要: 在天然气水合物发育区海底沉积物中甲烷厌氧氧化作用(AOM)是碳循环的重要组成部分。通过定量计算表层沉积物中甲烷迁移转化通量,可以更准确评估甲烷来源碳对沉积物碳库和海洋深部碳库影响。本文利用反应―运移模型对采集于南海神狐水合物发育区两个站位(SH-W19-PC、SH-W23-PC)采集的孔隙水SO42-、溶解无机碳(DIC)、Ca2+剖面进行拟合,同时对DIC碳同位素进行分析,确定近海底沉积物中的碳循环。研究显示两个站位孔隙水中SO42-和Ca2+浓度在剖面上随深度呈线性减少,DIC浓度随深度逐渐增加,其δ13CDIC值随深度逐渐降低至约-25‰,表明两个站位存在一定程度的AOM。模拟计算两个站位沉积物孔隙水溶解甲烷向上的通量分别为25.9和18.4 mmol·m-2 a-1,AOM作用产生的DIC分别占其总DIC量的70.7%和60%。由沉积物向海水中释放的DIC通量占DIC汇的约60%。因此,在天然气水合物发育区向海底渗漏甲烷大部分以DIC的形式进入上覆海水,这些具有极负碳同位素值的甲烷来源的DIC可能对局部深海碳库产生一定的影响。
  • 图  1  采样站位置图 (红色实点为SH-W19-PC,黑色实点为SH-W23-PC,白色实点为对比站位)

    Figure  1.  Locations of the sampling sites (red solid dot is SH-W19-PC, black solid dot is SH-W23-PC, white solid dots are contrast stations)

    图  2  SH-W19-PC和SH-W23-PC站位SO42−、Ca2+、DIC、CH4浓度在深度剖面上实测值和模型拟合结果 (红色点代表实测浓度值,黑色线为模型拟合曲线)

    Figure  2.  Measured values and simulate depth profiles of core SH-W19-PC and SH-W23-PC.Down-depth concentration of SO42−, DIC, Ca2+and CH4 are shown (red dots represent measured concentration values, black lines are model fitting curves)

    图  3  SH-W19-PC和SH-W23-PC柱状沉积物孔隙水DIC的δ13C值随深度变化剖面 (红色实心点代表数据点)

    Figure  3.  Pore water depth profiles of δ13CDIC values in sediments cores (SH-W19-PC and SH-W23-PC)

    图  4  SH-W19-PC和SH-W23-PC站位沉积表层碳转化过程及模型定量计算结果示意图

    Figure  4.  Schematic diagram of quantitative calculation results of surface carbon conversion process and model of SH-W19-PC and SH-W23-PC sites

    表  1  SH-W19-PC和SH-W23-PC站位模型参数及边界条件

    Table  1.   Model parameters and boundary conditions at sites SH-W19-PC and SH-W23-PC

    参数SH-W19-PCSH-W23-PC单位
    温度(T55
    盐度(S33.533.5%
    压力(P105105.1bar
    干燥固体密度(ρS2.52.5g/cm3
    沉积速率(ωa0.0330.033cm/a
    沉积物-水界面孔隙度(φ0b0.70.7
    POC初始年龄(a0c4040ka
    SO42−在海水中的扩散系数c191191cm2/a
    CH4在海水中的扩散系数c294294cm2/a
    DIC在海水中的扩散系数d203203cm2/a
    Ca2+在海水中的扩散系数c142142cm2/a
    POC抑制系数(KCe3535mmol/L
    POC降解米氏常数(KSO42-e0.10.1mmol/L
    AOM的动力系数(KAOMe11cm3/(a·mmol)
    SO42−的上边界条件28.228.2mmol/L
    DIC的上边界条件2.22.3mmol/L
    Ca2+的上边界条件1111mmol/L
      a为南海北部6个ODP184钻孔的平均值[52]b据Wang 等[52]c据Hu Y等[33]d基于碳酸氢根离子的分子扩散系数[27]e据Wallman 等[50]
    下载: 导出CSV

    表  2  反应深度积分速率(单位:mmol∙m−2 a−1,C)

    Table  2.   Depth-integrated rates(unit, mmol∙m−2 a−1,C)

    站位POC降解速率OSR速率ME速率AOM速率
    SH-W19-PC14.49.22.628.5
    SH-W23-PC15.411.61.920.2
    D8-1333.3330.20.1
    D17-1519.0114.030.1
    W01-1624.214.74.720.9
      表格中D8-13和D17-15为南海东沙海域两个站位数据[32-33],W01-16为神狐站位数据[33]
    下载: 导出CSV
  • [1] Kvenvolden K A. Gas hydrate-geological perspective and global change [J]. Reviews of Geophysics, 1993, 31(2): 173-187. doi: 10.1029/93RG00268
    [2] Buffett B, Archer D. Global inventory of methane clathrate: sensitivity to changes in the deep ocean [J]. Earth and Planetary Science Letters, 2004, 227(3-4): 185-199. doi: 10.1016/j.jpgl.2004.09.005
    [3] Wallmann K, Pinero E, Burwicz E, et al. The global inventory of methane hydrate in marine sediments: a theoretical approach [J]. Energies, 2012, 5(7): 2449-2498. doi: 10.3390/en5072449
    [4] Milkov A V. Global estimates of hydrate-bound gas in marine sediments: how much is really out there? [J]. Earth-Science Reviews, 2004, 66(3-4): 183-197. doi: 10.1016/j.earscirev.2003.11.002
    [5] Knittel K, Boetius A. Anaerobic oxidation of methane: progress with an unknown process [J]. Annual Review of Microbiology, 2009, 63(1): 311-334. doi: 10.1146/annurev.micro.61.080706.093130
    [6] Valentine D L, Kastner M, Wardlaw G D, et al. Biogeochemical investigations of marine methane seeps, Hydrate Ridge, Oregon [J]. Journal of Geophysical Research: Biogeosciences, 2015, 110(G2): G02005.
    [7] Stakes D S, Orange D, Paduan J B, et al. Cold-seeps and authigenic carbonate formation in Monterey Bay, California [J]. Marine Geology, 1999, 159(1-4): 93-109. doi: 10.1016/S0025-3227(98)00200-X
    [8] Komada T, Burdige D J, Li H L, et al. Organic matter cycling across the sulfate-methane transition zone of the Santa Barbara Basin, California borderland [J]. Geochimica et Cosmochimica Acta, 2016, 176: 259-278. doi: 10.1016/j.gca.2015.12.022
    [9] Pohlman J W, Bauer J E, Waite W F, et al. Methane hydrate-bearing seeps as a source of aged dissolved organic carbon to the oceans [J]. Nature Geoscience, 2010, 4(1): 37-41.
    [10] Haeckel M, Boudreau B P, Wallmann K. Bubble-induced porewater mixing: A 3-D model for deep porewater irrigation [J]. Geochimica et Cosmochimica Acta, 2007, 71(21): 5135-5154. doi: 10.1016/j.gca.2007.08.011
    [11] Dickens G R. Global change: hydrocarbon-driven warming [J]. Nature, 2004, 429(6991): 513-515. doi: 10.1038/429513a
    [12] Whiticar M J, Faber E, Schoell M. Biogenic methane formation in marine and freshwater environments: CO2 reduction vs. acetate fermentation—Isotope evidence [J]. Geochimica et Cosmochimica Acta, 1986, 50(5): 693-709. doi: 10.1016/0016-7037(86)90346-7
    [13] Chen D F, Li X X, Xia B. Distribution features of gas hydrate stable zones and resource prediction of the Qiongdongnan basin in the South China Sea [J]. Chinese Journal of Geophysics, 2004, 47(3): 548-555. doi: 10.1002/cjg2.519
    [14] Zhang H Q, Yang S X, Wu N Y, et al. Successful and surprising results for China's first gas hydrate drilling expedition [J]. Fire in the Ice: Methane Hydrate Newsletter, 2007, 7: 6-9.
    [15] Zhang G X, Yang S X, Zhang M, Liang J, et al. GMGS2 expedition investigates rich and complex gas hydrate environment in the South China Sea [J]. Fire in the Ice: Methane Hydrate Newsletter, 2014, 14(1): 1-5.
    [16] Yang S X, Zhang M, Liang J Q, et al. Preliminary results of China's third gas hydrate drilling expedition: a critical step from discovery to development in the South China Sea [J]. Fire in the Ice: Methane Hydrate Newsletter, 2015, 15(2): 1-5.
    [17] Wu S, Zhang G, Huang Y, et al. Gas hydrate occurrence on the continental slope of the northern South China Sea [J]. Marine and Petroleum Geology, 2005, 22(3): 403-412. doi: 10.1016/j.marpetgeo.2004.11.006
    [18] Deng H, Yan P, Liu H L, et al. Seismic data processing and the characterization of a gas hydrate bearing zone offshore of southwestern Taiwan [J]. Terrestrial, Atmospheric and Oceanic Sciences, 2006, 17(4): 781-797. doi: 10.3319/TAO.2006.17.4.781(GH)
    [19] Yan P, Deng H, Liu H L. The geological structure and prospect of Gas Hydrate over the Dongsha slope, South China Sea [J]. Terrestrial, Atmospheric and Oceanic Sciences, 2006, 17(4): 645-658. doi: 10.3319/TAO.2006.17.4.645(GH)
    [20] Suess E. RV SONNE Cruise Report SO177, Sino–German Cooperative Project, South China Sea Continental Margin: Geological Methane Budget and Environmental Effects of Methane Emissions and Gashydrates. IFM-GEOMAR Reports. 2005.
    [21] 李学杰. 南海北部陆坡水合物初步研究: ODP钻孔证据[J]. 南海地质研究, 2004(1):16-28. [LI Xuejie. Preliminary study of gas hydrate in the northern slope of the south china sea: evidences from ODP [J]. South China Sea Geology Research, 2004(1): 16-28.
    [22] Hui G, Li S Z, Guo L L, et al. Source and accumulation of gas hydrate in the northern margin of the South China Sea [J]. Marine and Petroleum Geology, 2006, 69: 127-145.
    [23] Feng J X, Yang S X, Liang J Q, et al. Methane seepage inferred from the porewater geochemistry of shallow sediments in the Beikang Basin of the southern South China Sea [J]. Journal of Asian Earth Sciences, 2018, 168: 77-86. doi: 10.1016/j.jseaes.2018.02.005
    [24] Chen F, Hu Y, Feng D, et al. Evidence of intense methane seepages from molybdenum enrichments in gas hydrate-bearing sediments of the northern South China Sea [J]. Chemical Geology, 2016, 443: 173-181. doi: 10.1016/j.chemgeo.2016.09.029
    [25] Hu Y, Feng D, Liang Q Y, et al. Impact of anaerobic oxidation of methane on the geochemical cycle of redox-sensitive elements at cold-seep sites of the northern South China Sea [J]. Deep Sea Research Part II: Topical Studies in Oceanography, 2015, 122: 84-94. doi: 10.1016/j.dsr2.2015.06.012
    [26] Han X Q, Suess E, Huang Y Y, et al. Jiulong methane reef: Microbial mediation of seep carbonates in the South China Sea [J]. Marine Geology, 2008, 249(3-4): 243-256. doi: 10.1016/j.margeo.2007.11.012
    [27] Luo M, Chen L Y, Wang S H, et al. Pockmark activity inferred from pore water geochemistry in shallow sediments of the pockmark field in southwestern Xisha Uplift, northwestern South China Sea [J]. Marine and Petroleum Geology, 2013, 48: 247-259. doi: 10.1016/j.marpetgeo.2013.08.018
    [28] Tong H P, Feng D, Cheng H, et al. Authigenic carbonates from seeps on the northern continental slope of the South China Sea: New insights into fluid sources and geochronology [J]. Marine and Petroleum Geology, 2013, 43: 260-271. doi: 10.1016/j.marpetgeo.2013.01.011
    [29] Feng D, Cheng M, Kiel S, et al. Using Bathymodiolus tissue stable carbon, nitrogen and sulfur isotopes to infer biogeochemical process at a cold seep in the South China Sea [J]. Deep Sea Research Part I: Oceanographic Research Papers, 2015, 104: 52-59. doi: 10.1016/j.dsr.2015.06.011
    [30] Lin Z Y, Sun X M, Peckmann J, et al. How sulfate-driven anaerobic oxidation of methane affects the sulfur isotopic composition of pyrite: A SIMS study from the South China Sea [J]. Chemical Geology, 2016, 440: 26-41. doi: 10.1016/j.chemgeo.2016.07.007
    [31] Chuang P C, Dale A W, Wallmann K, et al. Relating sulfate and methane dynamics to geology: accretionary prism offshore SW Taiwan [J]. Geochemistry, Geophysics, Geosystems, 2013, 14(7): 2523-2545. doi: 10.1002/ggge.20168
    [32] 梁华崔, 梁前勇, 胡钰, 等. 南海东沙海域浅表层柱状沉积物孔隙水地球化学特征及对冷泉流体活动的指示[J]. 地球化学, 2017, 46(4):333-344. [LIANG Huacui, LIANG Qianyong, HU Yu, et al. Pore water geochemistry of shallow surface sediments in the Dongsha area of the South China Sea and its implications for the activities of cold seep fluids [J]. Geochimica, 2017, 46(4): 333-344. doi: 10.3969/j.issn.0379-1726.2017.04.004
    [33] Hu Y, Luo M, Chen L Y, et al. Methane source linked to gas hydrate system at hydrate drilling areas of the South China Sea: Porewater geochemistry and numerical model constraints [J]. Journal of Asian Earth Sciences, 2018, 168: 87-95. doi: 10.1016/j.jseaes.2018.04.028
    [34] 吴雪停, 刘丽华, Haeckel M, et al. 南海北部深海浅层沉积物中甲烷生物地球化学过程数值模拟研究[J]. 海洋地质与第四纪地质, 2016, 36(3):81-90. [WU Xueting, LIU Lihua, Haeckel M, et al. Simulation of the biogeochemical processes in methane-bearing surface sediments of Haiyang 4 area, northern slope of south china sea [J]. Marine Geology & Quaternary Geology, 2016, 36(3): 81-90.
    [35] Luo M, Dale A W, Wallmann K, et al. Estimating the time of pockmark formation in the SW Xisha Uplift (South China Sea) using reaction-transport modeling [J]. Marine Geology, 2015, 364: 21-31. doi: 10.1016/j.margeo.2015.03.006
    [36] Hu Y, Luo M, Liang Q Y, et al. Pore fluid compositions and inferred fluid flow patterns at the Haima cold seeps of the South China sea [J]. Marine and Petroleum Geology, 2016, 103: 29-40.
    [37] Ye H, Yang T, Zhu G R, et al. Pore water geochemistry in shallow sediments from the northeastern continental slope of the South China sea [J]. Marine and Petroleum Geology, 2016, 75: 68-82. doi: 10.1016/j.marpetgeo.2016.03.010
    [38] 冯俊熙, 杨胜雄, 梁金强, 等. 南海北部神狐东南海域沉积物孔隙水地球化学特征及其对天然气水合物的指示[J]. 海洋地质前沿, 2017, 33(7):32-44. [FENG Junxi, YANG Shengxiong, LIANG Jinqiang, et al. Pore water geochemistry in shallow sediments from southeastern shenhu area of northern south china sea and their implications for gas hydrate occurrence [J]. Marine Geology Frontiers, 2017, 33(7): 32-44.
    [39] 郭依群, 杨胜雄, 梁金强, 等. 南海北部神狐海域高饱和度天然气水合物分布特征[J]. 地学前缘, 2017, 24(4):24-31. [GUO Yiqun, YANG Shengxiong, LIANG Jinqiang, et al. Characteristics of high gas hydrate distribution in the Shenhu area on the northern slope of the South China Sea [J]. Earth Science Frontiers, 2017, 24(4): 24-31.
    [40] Guo Y Q, Yang S X, Liang J Q, et al. Characteristics of high gas hydrate distribution in the Shenhu area on the northern slope of the South China Sea [J]. Earth Science Frontiers, 2017, 24(4): 24-31.
    [41] Taylor B, Hayes D E. Origin and history of the South China Sea Basin[M]//Hayes D E. The Tectonic and Geologic Evolution of Southeast Asian Seas and Islands. Washington, D.C.: American Geophysical Union, 1983, 27: 23-56.
    [42] Sun Y B, Wu F, Clemens S C, et al. Processes controlling the geochemical composition of the South China Sea sediments during the last climatic cycle [J]. Chemical Geology, 2008, 257(3-4): 240-246. doi: 10.1016/j.chemgeo.2008.10.002
    [43] 王存武, 陈红汉, 陈长民, 等. 珠江口盆地白云深水扇特征及油气成藏主控因素[J]. 地球科学—中国地质大学学报, 2007, 32(2):247-252, 266. [WANG Cunwu, CHEN Honghan, CHEN Changmin, et al. Characteristics of the Baiyun deep-water fan and main accumulation controlling factors in Pearl River Mouth Basin, South China Sea [J]. Earth Science—Journal of China University of Geosciences, 2007, 32(2): 247-252, 266.
    [44] Chen H, Savage P S, Teng F Z, et al. Zinc isotope fractionation during magmatic differentiation and the isotopic composition of the bulk Earth [J]. Earth and Planetary Science Letters, 2013, 369-370: 34-42. doi: 10.1016/j.jpgl.2013.02.037
    [45] 王家豪, 庞雄, 王存武, 等. 珠江口盆地白云凹陷中央底辟带的发现及识别[J]. 地球科学—中国地质大学学报, 2006, 31(2):209-213. [WANG Jiahao, PANG Xiong, WANG Cunwu, et al. Discovery and recognition of the central diapiric zone in Baiyun depression, Pearl River Mouth basin [J]. Earth Science-Journal of China University of Geosciences, 2006, 31(2): 209-213.
    [46] 吴庐山, 杨胜雄, 梁金强, 等. 南海北部神狐海域沉积物中孔隙水硫酸盐梯度变化特征及其对天然气水合物的指示意义[J]. 中国科学: 地球科学, 2013, 56(3):530-540. [WU Lushan, YANG Shengxiong, LIANG Jinqiang, et al. Variations of pore water sulphate gradients in sediments in indicator for underlying gas hydrate in shenhu area, the South China sea [J]. Chinese Science: Earth Science, 2013, 56(3): 530-540.
    [47] 张伟, 梁金强, 何家雄, 等. 南海北部神狐海域GMGS1和GMGS3钻探区天然气水合物运聚成藏的差异性[J]. 天然气工业, 2018, 38(3):138-149. [ZHANG Wei, LIANG Jiaqiang, HE Jiaxiong, et al. Differences in natural gas hydrate migration and accumulation between GMGS1 and GMGS3 drilling areas in the Shenhu area, northern South China Sea [J]. Natural Gas Industry, 2018, 38(3): 138-149.
    [48] 吴能友, 杨胜雄, 王宏斌, 等. 南海北部陆坡神狐海域天然气水合物成藏的流体运移体系[J]. 地球物理学报, 2009, 52(6):1641-1650. [WU Nengyou, YANG Shengxiong, WANG Hongbin, et al. Gas-bearing fluid influx sub-system for gas hydrate geological system in Shenhu Area, Northern South China Sea [J]. Chinese Journal of Geophysics, 2009, 52(6): 1641-1650.
    [49] Sayles F L. Early Diagenesis: a theoretical approach. Robert A. Berner, Princeton University Press; Princeton, N.J., 1980; XII, 241 p., Princeton Series in Geochemistry, 25 cloth; 9 paper [J]. Geochimicaet Cosmochimica Acta, 1982, 46(2): 303-303.
    [50] Wallmann K, Aloisi G, Haeckel M, et al. Kinetics of organic matter degradation, microbial methane generation, and gas hydrate formation in anoxic marine sediments [J]. Geochimica et Cosmochimica Acta, 2006, 70(15): 3905-3927. doi: 10.1016/j.gca.2006.06.003
    [51] Van Cappellen P, Wang Y F. Cycling of iron and manganese in surface sediments; a general theory for the coupled transport and reaction of carbon, oxygen, nitrogen, sulfur, iron, and manganese [J]. American Journal of Science, 1996, 296(3): 197-243. doi: 10.2475/ajs.296.3.197
    [52] Wang P, Prell W L, Blum P, et al. Proceedings of the Ocean Drilling Program, Initial Report 184, College Station, TX: Ocean Drilling Program, 2000.
    [53] Jørgensen B B, Kasten S. Sulfur cycling and methane oxidation[M]//Schulz H D, Zabel M. Marine Geochemistry. Berlin, Heidelberg: Springer, 2006: 271-309.
    [54] Reeburgh W S. Oceanic methane biogeochemistry [J]. Chemical Reviews, 2007, 107(2): 486-513. doi: 10.1021/cr050362v
    [55] Regnier P, Dale A W, Arndt S, et al. Quantitative analysis of anaerobic oxidation of methane (AOM) in marine sediments: A modeling perspective [J]. Earth-Science Reviews, 2011, 106(1-2): 105-130. doi: 10.1016/j.earscirev.2011.01.002
    [56] Middelburg J J. A simple rate model for organic matter decomposition in marine sediments [J]. Geochimica et Cosmochimica Acta, 1989, 53(7): 1577-1581. doi: 10.1016/0016-7037(89)90239-1
    [57] Reeburgh W S. Methane consumption in Cariaco Trench waters and sediments [J]. Earth and Planetary Science Letters, 1976, 28(3): 337-344. doi: 10.1016/0012-821X(76)90195-3
    [58] Boetius A, Ravenschlag K, Schubert C J, et al. A marine microbial consortium apparently mediating anaerobic oxidation of methane [J]. Nature, 2000, 407(6804): 623-626. doi: 10.1038/35036572
    [59] Zhang Y P, Luo M, Hu Y, et al. An areal assessment of subseafloor carbon cycling in cold seeps and hydrate-bearing areas in the northern South China Sea [J]. Geofluids, 2019, 2019: 2573937.
    [60] Borowski W R, Paull C K, Ussler III W. Global and local variations of interstitial sulfate gradients in deep-water, continental margin sediments: Sensitivity to underlying methane and gas hydrates [J]. Marine Geology, 1999, 159(1-4): 131-154. doi: 10.1016/S0025-3227(99)00004-3
    [61] Egger M, Riedinger N, Mogollón J M, et al. Global diffusive fluxes of methane in marine sediments [J]. Nature Geoscience, 2018, 11(6): 421-425. doi: 10.1038/s41561-018-0122-8
    [62] Wu L S, Yang S X, Liang J Q, et al. Variations of pore water sulfate gradients in sediments as indicator for underlying gas hydrate in Shenhu Area, the South China Sea [J]. Science China Earth Sciences, 2012, 56(4): 530-540.
    [63] Hong W L, Sauer S, Panieri G, et al. Removal of methane through hydrological, microbial, and geochemical processes in the shallow sediments of pockmarks along eastern Vestnesa Ridge (Svalbard) [J]. Limnology and Oceanography, 2016, 61(S1): S324-S343. doi: 10.1002/lno.10299
    [64] Hong W L, Torres M E, Kim J H, et al. Carbon cycling within the sulfate-methane-transition-zone in marine sediments from the Ulleung Basin [J]. Biogeochemistry, 2013, 115(1): 129-148.
    [65] Berner R A. The carbon cycle and carbon dioxide over Phanerozoic time: the role of land plants [J]. Philosophical Transactions of the Royal Society of London. Series B: Biological Sciences, 1998, 353(1365): 75-82. doi: 10.1098/rstb.1998.0192
    [66] Jahnke R A. The global ocean flux of particulate organic carbon: Areal distribution and magnitude [J]. Global Biogeochemical Cycles, 1996, 10(1): 71-88. doi: 10.1029/95GB03525
    [67] Krumins V, Gehlen M, Arndt S, et al. Dissolved inorganic carbon and alkalinity fluxes from coastal marine sediments: model estimates for different shelf environments and sensitivity to global change [J]. Biogeosciences, 2013, 10: 371-398.
    [68] Chuang P C, Yang T F, Wallmann K, et al. Carbon isotope exchange during anaerobic oxidation of methane (AOM) in sediments of the northeastern South China Sea [J]. Geochimica et Cosmochimica Acta, 2019, 246: 138-155. doi: 10.1016/j.gca.2018.11.003
    [69] Schubert C J, Nielsen B. Effects of decarbonation treatments on δ13C values in marine sediments [J]. Marine Chemistry, 2000, 72(1): 55-59. doi: 10.1016/S0304-4203(00)00066-9
    [70] Guan H X, Feng D, Wu N Y, et al. Fatty-acids and their δ13C characteristics of seep carbonates from the northern continental slope of Gulf of Mexico [J]. Chinese Science Bulletin, 2009, 55(8): 730-735.
    [71] Campbell K A, Farmer J D, Des Marais D. Ancient hydrocarbon seeps from the Mesozoic convergent margin of California: carbonate geochemistry, fluids and palaeoenvironments [J]. Geofluids, 2010, 2(2): 63-94.
    [72] 佟宏鹏, 冯东, 陈多福. 南海北部冷泉碳酸盐岩的矿物、岩石及地球化学研究进展[J]. 热带海洋学报, 2012, 31(5):45-56. [TONG Hongpeng, FENG Dong, CHEN Duofu. Progresses on petrology, mineralogy and geochemistry of cold seep carbonates in the northern South China Sea [J]. Journal of Tropical Oceanography, 2012, 31(5): 45-56.
    [73] Burdige D J, Komada T, Magen C, et al. Modeling studies of dissolved organic matter cycling in Santa Barbara Basin (CA, USA) sediments [J]. Geochimica et Cosmochimica Acta, 2016, 195: 100-119. doi: 10.1016/j.gca.2016.09.007
    [74] Burdige D J. Sediment pore waters[M]//Hansell D A, Carlson C A. Biogeochemistry of Marine Dissolved Organic Matter. San Diego, CA: Academic Press, 2002: 611-663.
    [75] Hung C W, Huang K H, Shih Y Y, et al. Benthic fluxes of dissolved organic carbon from gas hydrate sediments in the northern South China Sea [J]. Scientific Reports, 2016, 6: 29597. doi: 10.1038/srep29597
  • [1] 计超, 徐利强, 张一辉, 郭敏, 孔德明.  南海琼东上升流区过去1 900年海洋生产力记录 . 海洋地质与第四纪地质, 2020, 40(4): 1-10. doi: 10.16562/j.cnki.0256-1492.2019092502
    [2] 张国良, 战明君.  板块俯冲和岩浆过程中碳循环及深部碳储库 . 海洋地质与第四纪地质, 2019, 39(5): 36-45. doi: 10.16562/j.cnki.0256-1492.2019092201
    [3] 蔡观强, 李顺, 赵利, 高红芳, 钟和贤.  南海海盆中部表层沉积物地球化学特征 . 海洋地质与第四纪地质, 2018, 38(5): 90-101. doi: 10.16562/j.cnki.0256-1492.2018.05.009
    [4] 李金澜, 田军.  末次冰期南海南部暴露的巽他陆架是大气碳汇? . 海洋地质与第四纪地质, 2018, 38(4): 155-163. doi: 10.16562/j.cnki.0256-1492.2018.04.013
    [5] 张强, 吕福亮, 贺晓苏, 王彬, 毛超林.  南海成藏组合发育特征及勘探潜力评价 . 海洋地质与第四纪地质, 2017, 37(6): 158-167. doi: 10.16562/j.cnki.0256-1492.2017.06.017
    [6] 张素芳, 张向涛, 张青林, 佘清华.  南海北部白垩系发育特征及构造意义 . 海洋地质与第四纪地质, 2015, 35(6): 81-86. doi: 10.16562/j.cnki.0256-1492.2015.06.008
    [7] 傅飘儿, 庄畅, 刘坚, 陈道华, 张欣, 陈思海.  南海西沙海槽XH-CL16柱状沉积物稀土元素特征及其物源 . 海洋地质与第四纪地质, 2015, 35(4): 63-71. doi: 10.16562/j.cnki.0256-1492.2015.04.007
    [8] 刘毅, 孙立广, 罗宇涵, 孙松, 王玉宏.  南海湖泊沉积物中的陆源粉尘记录 . 海洋地质与第四纪地质, 2013, 33(3): 1-8. doi: 10.3724/SP.J.1140.2013.03001
    [9] 李三忠, 索艳慧, 刘鑫, 戴黎明, 余珊, 赵淑娟, 马云, 王霄飞, 程世秀, 薛友辰, 熊莉娟, 安慧婷.  南海的基本构造特征与成因模型:问题与进展及论争 . 海洋地质与第四纪地质, 2012, 32(6): 35-53. doi: 10.3724/SP.J.1140.2012.06035
    [10] 张向涛, 陈亮, 佘清华, 张素芳, 乔培军, 邵磊.  南海北部古韩江物源的演化特征 . 海洋地质与第四纪地质, 2012, 32(4): 41-48. doi: 10.3724/SP.J.1140.2012.04041
    [11] 陆红锋, 刘坚, 陈芳, 程思海, 廖志良.  南海东北部硫酸盐还原-甲烷厌氧氧化界面——海底强烈甲烷渗溢的记录 . 海洋地质与第四纪地质, 2012, 32(1): 93-98. doi: 10.3724/SP.J.1140.2012.01093
    [12] 陈芳, 周洋, 苏新, 刘广虎, 陆红锋, 王金莲.  南海神狐海域含水合物层粒度变化及与水合物饱和度的关系 . 海洋地质与第四纪地质, 2011, 31(5): 95-100. doi: 10.3724/SP.J.1140.2011.05095
    [13] 冯晓萍, 蔡进功.  沉积物的颗粒大小与所含有机质关系的研究进展 . 海洋地质与第四纪地质, 2010, 30(6): 141-148. doi: 10.3724/SP.J.1140.2010.06141
    [14] 蔡观强, 邱燕, 彭学超, 钟和贤.  南海西南海域表层沉积物微量和稀土元素地球化学特征及其意义 . 海洋地质与第四纪地质, 2010, 30(5): 53-62. doi: 10.3724/SP.J.1140.2010.05053
    [15] 邬黛黛, 吴能友, 付少英, 梁金强, 管红香.  南海北部东沙海域水合物区浅表层沉积物的地球化学特征 . 海洋地质与第四纪地质, 2010, 30(5): 41-51. doi: 10.3724/SP.J.1140.2010.05041
    [16] 陆红锋, 陈芳, 刘坚, 周洋, 廖志良.  南海东北部甲烷成因碳酸盐岩的矿物学及同位素组成 . 海洋地质与第四纪地质, 2010, 30(2): 51-59. doi: 10.3724/SP.J.1140.2010.02051
    [17] 陈芳, 周洋, 苏新, 陆红锋, 刘广虎, 庄畅.  南海神狐海域含水合物层底栖有孔虫群落结构与同位素组成 . 海洋地质与第四纪地质, 2010, 30(2): 1-8. doi: 10.3724/SP.J.1140.2010.02001
    [18] 周洋, 陈芳, 苏新, 刘广虎, 陈超云.  南海东沙海域HD319岩心富甲烷环境底栖有孔虫群落结构 . 海洋地质与第四纪地质, 2009, 29(3): 1-8. doi: 10.3724/SP.J.1140.2009.03001
    [19] 姚伯初, 杨木壮.  南海晚新生代构造运动与天然气水合物资源 . 海洋地质与第四纪地质, 2008, 28(4): 93-100. doi: 10.3724/SP.J.1140.2008.03093
    [20] 叶芳, 刘志飞, 拓守廷, 翦知湣.  南海北部中更新世0.78~1.0 Ma期间的陆源碎屑粒度记录 . 海洋地质与第四纪地质, 2007, 27(2): 77-83.
  • 加载中
图(4) / 表(2)
计量
  • 文章访问数:  1960
  • HTML全文浏览量:  27
  • PDF下载量:  7
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-04-24
  • 修回日期:  2019-07-08
  • 网络出版日期:  2020-04-29
  • 刊出日期:  2020-06-01

南海神狐海域水合物发育区浅表层沉积物甲烷周转定量模拟

doi: 10.16562/j.cnki.0256-1492.2019042401
    作者简介:

    胡廷苍(1993―),男,硕士研究生,主要从事海洋地质研究,E-mail:m170300611@st.shou.edu.cn

    通讯作者: 陈多福(1962―),男,教授,博士生导师,主要从事海洋地质研究,E-mail:dfchen@shou.edu.cn
基金项目:  国家重点研究发展计划冷泉重点项目“中国海域冷泉系统演变过程及其机制”(2018YFC0310003);国家自然科学基金项目“南海北部冷泉和天然气水合物发育区海底浅表层沉积物碳循环数值模拟”(41730528),“冲绳海槽海底冷泉—热液系统相互作用及资源效应”(91858208)
  • 中图分类号: P736.4, P738

摘要: 在天然气水合物发育区海底沉积物中甲烷厌氧氧化作用(AOM)是碳循环的重要组成部分。通过定量计算表层沉积物中甲烷迁移转化通量,可以更准确评估甲烷来源碳对沉积物碳库和海洋深部碳库影响。本文利用反应―运移模型对采集于南海神狐水合物发育区两个站位(SH-W19-PC、SH-W23-PC)采集的孔隙水SO42-、溶解无机碳(DIC)、Ca2+剖面进行拟合,同时对DIC碳同位素进行分析,确定近海底沉积物中的碳循环。研究显示两个站位孔隙水中SO42-和Ca2+浓度在剖面上随深度呈线性减少,DIC浓度随深度逐渐增加,其δ13CDIC值随深度逐渐降低至约-25‰,表明两个站位存在一定程度的AOM。模拟计算两个站位沉积物孔隙水溶解甲烷向上的通量分别为25.9和18.4 mmol·m-2 a-1,AOM作用产生的DIC分别占其总DIC量的70.7%和60%。由沉积物向海水中释放的DIC通量占DIC汇的约60%。因此,在天然气水合物发育区向海底渗漏甲烷大部分以DIC的形式进入上覆海水,这些具有极负碳同位素值的甲烷来源的DIC可能对局部深海碳库产生一定的影响。

English Abstract

  • 自从20世纪70年代首次发现天然气水合物后,天然气水合物作为一种潜在的洁净能源,越来越受到人们的关注。海洋沉积物中的天然气水合物(>600 Gt)是地球上最大的甲烷储库[1-2],通常分布于大陆边缘海底[3]。然而作为一个动态的甲烷储库,在温压条件改变、水合物相平衡状态遭到破坏时,水合物会发生分解,可导致甲烷进入海水在海底形成冷泉。由于甲烷是一种强烈的温室气体,在地质历史时期中由水合物分解大规模释放的甲烷可引起海水酸化和全球气候变化[1,3-4]

    大陆边缘水合物发育区通常存在着富甲烷流体的渗漏现象。一旦甲烷运移到沉积表层遇到由海水扩散下来的硫酸根时,绝大部分甲烷被甲烷厌氧氧化作用(AOM)消耗[5]。这些渗漏的甲烷作为除沉积有机质之外的有机碳源可为微生物群落提供能量,因此,冷泉渗漏区的生物量比非渗漏区的要高,且与渗漏甲烷的强度有一定的联系[6]。同时AOM也会导致孔隙水碱度升高,从而促进沉积物中自生碳酸盐岩的沉淀[7]。冷泉区孔隙水中的DIC和DOC通常具有较高的浓度,并且它们的13C和14C 出现明显亏损,表明这些碳来源于甲烷[8-9]。在甲烷渗漏强烈的地区,大量的甲烷通常是以气泡形式向上快速运移,不能充分被AOM作用所消耗,而进入上覆水体[10-11]。此外,DIC除了在AOM过程中产生外,有机质硫酸根还原作用(OSR)也是DIC的重要来源[12]。这些甲烷和有机质来源的DIC在海水中的释放可能会对局部深海碳库产生影响。

    20世纪90年代以来,我国对南海地区进行了多次高分辨率地震调查,通过似海底反射层(BSR)推断南海北部东沙、西沙、台西南和神狐地区广泛发育有高饱和度的天然气水合物[13-24]。对水合物区沉积物和孔隙水的地球化学特征显示存在与甲烷活动有关的亏损13C的碳酸盐岩和硫酸根浓度剖面异常信号[25-30]。最近运用反应运移模型对南海东沙、神狐、西沙等水合物和冷泉发育区海底表层的生物化学反应速率和甲烷通量进行了评估,认识到甲烷活动在上述区域较为活跃,并在部分站位AOM作用主导了沉积表层碳循环过程[23,27,31-38],但对甲烷和有机质参与的微生物主导的各种碳转化通量和最终归宿的定量研究并不多。

    神狐海域位于南海北部陆坡,该海域天然气水合物发育区经过水合物钻探航次调查,显示有大面积富甲烷流体活动[39-41]。本文通过对SH-W19-PC和SH-W23-PC两个站位的地球化学分析显示有较强的AOM作用,因此,选取这两个站位重力活塞柱状沉积物孔隙水的地球化学数据进行数值模拟研究,应用反应运移模型定量评估神狐海域两个站位的与甲烷周转有关的生物地球化学反应速率,计算DIC排放到海水中的通量和深部沉积层来源的甲烷通量。通过综合不同碳库之间的转化关系和转化量,评估神狐海域天然气水合物发育区浅表层沉积物碳循环对海洋碳库和海洋环境的潜在影响。

    • 南海是西太平洋最大的边缘海,在欧亚板块、印度—澳大利亚板块和太平洋板块共同作用下形成[41]。大中型新生代沉积盆地在南海北部陆坡发育,晚上新世以来北部陆坡陆源有机质大量输入,沉积物快速堆积,为南海北部陆坡气源的形成提供了有利条件[42]。ODP 184航次3个站位均发现了明显富集烃类物质的沉积层,证实南海北部发育一系列含石油和天然气水合物的沉积盆地[43]

      神狐地区位于南海北部陆坡中段,构造上处于珠江口盆地珠Ⅱ坳陷区内[37]。受珠江口新生代沉积盆地形成演化的影响,神狐海域发育了广泛的海底峡谷、中央底辟带和各种类型的底辟构造[44-45],为甲烷流体的运移提供了良好通道。在新近纪神狐海域沉积速率变大,并且有机质含量较高[46]。上述有利的构造背景和沉积环境指示神狐区域具有形成天然气水合物的良好成藏条件[46]。2007年和2015年广州海洋地质调查局在该海区进行了GMGS1和GMGS3两个航次的水合物钻探,均采集到了水合物实物样品[47],同时神狐海域多个站位的沉积物及孔隙水地球化学分析也指示该区域存在富甲烷的流体活动[48]

    • SH-W19-PC和SH-W23-PC站位沉积物样品为2016年9月由广州海洋地质调查局“海洋四号”用重力活塞柱状采样器采集。两个站位均位于南海北部陆坡神狐水合物发育区内,水深分别为约1 300 和1 500 m(图1),柱状样全长分别为820 和847 cm。柱样采集后立刻进行孔隙水抽取,使用孔径为0.2 μm的Rhizon采样器采集,采样间隔为10~40 cm,将4 h内获得15~20 mL的孔隙水装在两个10 mL的玻璃瓶内,一瓶内加入10 μL饱和HgCl2用于DIC浓度及其碳同位素分析,另一瓶用高纯浓硝酸酸化处理,用于阴离子和溶解金属元素分析。所有的样品都在4 ℃下保存直至分析。

      图  1  采样站位置图 (红色实点为SH-W19-PC,黑色实点为SH-W23-PC,白色实点为对比站位)

      Figure 1.  Locations of the sampling sites (red solid dot is SH-W19-PC, black solid dot is SH-W23-PC, white solid dots are contrast stations)

    • 硫酸根(SO42−)和钙离子(Ca2+)浓度在中国科学院南海海洋研究所用Dionex ICS-5000+离子色谱仪检测,分析精度优于2%。仪器分析之前SO42− 浓度用超纯水稀释500倍,Ca2+ 浓度用超纯水稀释100倍。阴离子浓度分析时用28 mM的KOH作为淋洗液,IonPac AS11HC作为分析柱。对于阳离子分析,20 mM甲磺酸作为淋洗液,IonPac CS12A作为分析柱。DIC的浓度和同位素的分析在路易斯安那州立大学完成,通过ThermoFinnigan Gas Bench 和ThermoFinnigan Delta V Advantage两个仪器耦合进行检测,利用浓度范围为1 mM至30 mM的标准溶液获得标准曲线。DIC浓度的分析精度优于2‰,同位素精度优于0.1‰,同位素测试结果用V-PDB标准表示。

    • 利用反应运移模型模拟1个固相组分(POC)和4个溶解组分(SO42−、DIC、Ca2+、CH4),用下面两个偏微分方程描述溶解组分(式1)和POC(式2)[49-50]

      $ \phi \frac{{\partial {{\rm{C}}_a}}}{{\partial t}} = \frac{{\partial \left( {\phi \cdot{D_s}\cdot\frac{{\partial {C_a}}}{{\partial x}}} \right)}}{{\partial x}} - \frac{{\partial \left( {\phi \cdot{v_p}\cdot{C_a}} \right)}}{{\partial x}} + \phi \cdot\mathop \sum R $

      (1)

      $ \left( {1 - \phi } \right)\frac{{\partial POC}}{{\partial t}} = - \frac{{\partial \left( {\left( {1 - \phi } \right)\cdot{v_s}\cdot POC} \right)}}{{\partial x}} + \left( {1 - \phi } \right)\cdot\mathop \sum \nolimits R $

      (2)

      式中x(cm)为POC或溶解组分所在的沉积深度;t(a)为时间;φ为孔隙度;DS(cm2 a−1)是经曲度矫正过的分子扩散系数;Ca(μmol cm-3)是溶解组分的浓度;POC是沉积物中有机质的百分含量;vp(cm a−1)是孔隙水的埋藏速率;vs(cm a−1)是固体的沉积速率;是在该模型中某一组分变化所考虑的生物化学反应速率总和。在模拟的长度范围内主要考虑的生物化学过程包括:OSR、AOM、产甲烷作用和自生碳酸盐岩沉淀。由于在重力活塞柱状样品采集过程中常导致最顶部的10~20 cm沉积物样品的缺失,故通常不考虑有机质有氧矿化、反硝化作用和铁锰还原作用[51]

      由于孔隙度和沉积速率数据没有实际测得,故参考临近站位的实测值[52]。在海底没有外部流体对流的情况下固相和孔隙水的埋藏速率相同。迂曲度矫正后的分子扩散系数DS由下面方程得到:

      $ {D_s} = \frac{{{D_m}}}{{1 - {\rm{ln}}{{\left( \phi \right)}^2}}} $

      (3)

      式中Dm是海水中溶解组分在一定温度盐度和压力下的扩散系数。

      有机质矿化为海洋沉积物中的微生物活动提供了能量,在沉积物缺氧层内硫酸根是有机质氧化的主要电子受体[53],OSR过程可被简化为:

      $ 2C{H_2}O{\left( {PO_4^{3 - }} \right)_{{r_p}}} + SO_4^{2 - } \to 2HCO_3^ - + {H_2}S + 2{r_p}PO_4^{3 - } $

      (4)

      OSR过程中每消耗1 mol硫酸根会氧化2 mol有机碳,同时释放DIC、硫化物和磷酸盐到孔隙水中,rp在是有机磷和有机碳的比例。硫酸根消耗完后有机质产甲烷开始,方程式为:

      $ 2C{H_2}O{\left( {PO_4^{3 - }} \right)_{{r_p}}} \to C{O_2} + C{H_4} + 2{r_p}PO_4^{3 - } $

      (5)

      稳态情况下该过程产生等量的甲烷和CO2,自然环境中甲烷主要由有机质发酵和CO2还原作用产生[54]。在模型中OSR和产甲烷过程的反应速率速率取决于POC的总降解速率RPOC(g C g−1a−1[55]

      $ {R_{POC}} = \frac{{{K_c}}}{{\left[ {DIC} \right] + \left[ {C{H_4}} \right] + {K_c}}}\cdot\left( {0.16\cdot{{\left( {{a_0} + \frac{x}{w}} \right)}^{ - 0.95}}} \right)\cdot POC $

      (6)

      式中KC是孔隙水中DIC和CH4积累对POC降解的抑制系数;[DIC]和[CH4]分别代表DIC和CH4的浓度;a0是有机质初始年龄,反映其活性程度[56]。OSR速率(RSR,μmol cm−3a−1,以SO42−计),产甲烷速率(RMG,μmol cm−3a−1,以CH4计)和POC产DIC速率(RDP,μmol cm−3a−1,以C计)表达式如下:

      $ {R_{OSR}} = 0.5\cdot\frac{{{\rho _s}\cdot\left( {1 - \phi } \right)\cdot{{10}^6}}}{{M{W_c}\cdot\phi }}\cdot\frac{{\left[ {SO_4^{2 - }} \right]}}{{\left[ {SO_4^{2 - }} \right] + {K_{SO_4^{2 - }}}}}\cdot{R_{POC}} $

      (7)

      $ {R_{MG}} = 0.5\cdot\frac{{{\rho _s}\cdot\left( {1 - \phi } \right)\cdot{{10}^6}}}{{M{W_c}\cdot\phi }}\cdot\frac{{{K_{SO_4^{2 - }}}}}{{\left[ {SO_4^{2 - }} \right] + {K_{SO_4^{2 - }}}}}\cdot{R_{POC}} $

      (8)

      $ {R_{DP}} = \frac{{{\rho _s}\cdot\left( {1 - \phi } \right)\cdot{{10}^6}}}{{M{W_c}\cdot\phi }}\cdot{R_{POC}} - {R_{MG}} $

      (9)

      式中ρs是沉积物干重的密度,MWc是C的相对原子质量,KSO42−是在低SO42−浓度时的SO42−还原作用的抑制系数。

      AOM反应总的方程式为:

      $ C{H_4} + SO_4^{2 - } \to HCO_3^ - + H{S^ - } + {H_2}O $

      (10)

      AOM反应速率(RAOM)可由双分子动力学方程表示[55]

      $ {R_{AOM}} = {k_{AOM}}\cdot\left[ {SO_4^{2 - }} \right]\cdot\left[ {C{H_4}} \right] $

      (11)

      式中kAOM是AOM反应速率常数,Ca2+碳酸盐岩中的沉淀(或形成)速率(RPPT),可简化为由模型模拟出的浓度和实测浓度之差计算:

      $ {R_{PPT}} = {K_{PPT}}\cdot\left( {\left[ {C{a^{2 + }}} \right] - {{\left[ {C{a^{2 + }}} \right]}_T}} \right) $

      (12)

      式中KPPT是一级速率常数,[Ca2+]和[Ca2+]T分别是模拟和实测站位(SH-W19-PC和SH-W23-PC)各层位孔隙水Ca2+浓度。

      为了求解这个模型,设定模型上边界如表1所示,SO42−、Ca2+、DIC和CH4的浓度为固定值(狄利克雷型边界条件),上边界值为海水中各组分的浓度,由于没有实测甲烷数据,故甲烷浓度设为0;POC的上边界设为通量值(罗宾型边界条件);模型的下边界除了甲烷浓度外,各组分的浓度梯度都设为0(纽曼型边界条件);下边界甲烷浓度是通过较好的拟合SO42−浓度剖面获得。为了保证有机质有较高的降解程度,模型模拟长度为20 m,在不均一的网格内(由表层到深层分辨率逐渐降低)利用有限差分和直线法将式(1)、式(2)两个偏微分方程离散化,最后用MATHMATICA V.8.0中NDSolve命令结合反应速率方程求解并计算出各个反应过程的反应速率和各组分通量及转化量。

      表 1  SH-W19-PC和SH-W23-PC站位模型参数及边界条件

      Table 1.  Model parameters and boundary conditions at sites SH-W19-PC and SH-W23-PC

      参数SH-W19-PCSH-W23-PC单位
      温度(T55
      盐度(S33.533.5%
      压力(P105105.1bar
      干燥固体密度(ρS2.52.5g/cm3
      沉积速率(ωa0.0330.033cm/a
      沉积物-水界面孔隙度(φ0b0.70.7
      POC初始年龄(a0c4040ka
      SO42−在海水中的扩散系数c191191cm2/a
      CH4在海水中的扩散系数c294294cm2/a
      DIC在海水中的扩散系数d203203cm2/a
      Ca2+在海水中的扩散系数c142142cm2/a
      POC抑制系数(KCe3535mmol/L
      POC降解米氏常数(KSO42-e0.10.1mmol/L
      AOM的动力系数(KAOMe11cm3/(a·mmol)
      SO42−的上边界条件28.228.2mmol/L
      DIC的上边界条件2.22.3mmol/L
      Ca2+的上边界条件1111mmol/L
        a为南海北部6个ODP184钻孔的平均值[52]b据Wang 等[52]c据Hu Y等[33]d基于碳酸氢根离子的分子扩散系数[27]e据Wallman 等[50]
    • 图2图3为南海北部神狐区域SH-W19-PC和SH-W23-PC两个站位的SO42−、Ca2+、DIC和δ13CDIC随深度变化的剖面图。SH-W19-PC相较SH-W23-PC站位的SO42−浓度降低更快。SH-W19-PC站位由表层到760 cm处,SO42−浓度从28.2 mM减小至0。Ca2+ 浓度从表层的12.1 mM减小到底部3.7 mM。DIC浓度随深度明显增加,变化范围为2.2~35 mM。同时δ13CDIC在深度剖面上逐渐减小,在DIC浓度最大的深度720 cm处达到最小值(−25.1‰)。SH-W23-PC站位SO42−和Ca2+浓度在深度上几乎都显示线性减小趋势,SO42−浓度变化由表层28.2 mM减小至800 cm处的7.7 mM。Ca浓度从表层10.7 mM减小至底部4.7 mM。DIC浓度随深度逐渐增大而δ13CDIC随深度减小,在800 cm处δ13CDIC达到最小值−23.9‰。

      图  2  SH-W19-PC和SH-W23-PC站位SO42−、Ca2+、DIC、CH4浓度在深度剖面上实测值和模型拟合结果 (红色点代表实测浓度值,黑色线为模型拟合曲线)

      Figure 2.  Measured values and simulate depth profiles of core SH-W19-PC and SH-W23-PC.Down-depth concentration of SO42−, DIC, Ca2+and CH4 are shown (red dots represent measured concentration values, black lines are model fitting curves)

      图  3  SH-W19-PC和SH-W23-PC柱状沉积物孔隙水DIC的δ13C值随深度变化剖面 (红色实心点代表数据点)

      Figure 3.  Pore water depth profiles of δ13CDIC values in sediments cores (SH-W19-PC and SH-W23-PC)

    • 两个站位模拟孔隙水Ca2+、DIC、SO42−剖面都较好的拟合了实际数据。两个站位POC的初始年龄a0和有机质含量参考附近站位值[33],SH-W19-PC和SH-W23-PC两个站位模拟有机质降解速率分别为14.4和15.4 mmol∙m−2a−1表2)。SH-W19-PC的OSR反应速率(9.2 mmol∙m−2a−1)小于SH-W23-PC(11.6 mmol∙m−2a−1),两个站位的产甲烷速率相当。SH-W19-PC站位下部的甲烷通量(25.9 mmol∙m−2a−1)明显高于SH-W23-PC(18.4 mmol∙m−2a−1),同时SH-W19-PC的AOM反应速率(28.5 mmol∙m−2a−1)也比SH-W23-PC(20.2 mmol∙m−2a−1)要高,这主要是由于甲烷下边界不同而导致的通量不同。利用Ca2+浓度变化计算出两个站位的自生碳酸盐岩沉淀速率分别为9.3和6.6 mmol∙m−2a−1

      表 2  反应深度积分速率(单位:mmol∙m−2 a−1,C)

      Table 2.  Depth-integrated rates(unit, mmol∙m−2 a−1,C)

      站位POC降解速率OSR速率ME速率AOM速率
      SH-W19-PC14.49.22.628.5
      SH-W23-PC15.411.61.920.2
      D8-1333.3330.20.1
      D17-1519.0114.030.1
      W01-1624.214.74.720.9
        表格中D8-13和D17-15为南海东沙海域两个站位数据[32-33],W01-16为神狐站位数据[33]
    • 在海洋缺氧沉积物中,甲烷的厌氧氧化过程主要发生在硫酸盐甲烷转换带(SMTZ)内,SMTZ内的甲烷可能为多种来源(生物成因、热成因、水合物分解等),硫酸盐主要是海水中的硫酸根在浓度梯度作用下扩散至该带内[57-58]。SMTZ的深度受多种因素的控制,其中甲烷通量的大小对该带的影响显著[59-60]。随着甲烷通量的增加,SMTZ深度逐渐变浅,硫酸盐浓度剖面呈线性变化并且斜率随之增加[61]。对于全球陆坡区,SMTZ深度的平均深度为1 280 cm[61]。在水深为1 500 m附近的神狐海域两个站位,它们的SMTZ深度明显小于平均值,可能表明具有较大的甲烷通量。两个站位孔隙水DIC浓度随深度急剧增加,同时δ13CDIC逐渐亏损至−25‰附近,说明孔隙水DIC库中有AOM来源13C亏损的DIC加入[62]。比较两个站位发现SH-W19-PC的硫酸根浓度梯度大于SH-W23-PC,SMTZ深度小于SH-W23-PC,表明SH-W19-PC站位下部甲烷通量可能更大[60]。数值模拟结果证实了上述推断,SW-19站位下部来源的甲烷通量(25.9 mmol∙m−2a−1)大于SW-23(18.4 mmol∙m−2a−1)。同时两个站位的有机质产甲烷速率分别为2.6和1.9 mmol∙m−2a−1。而且模型计算的AOM反应速率SW-19站位也高于SW-23站位(分别为28.5和20.2 mmol∙m−2a−1)(表2),总体上与先前该海区的研究结果相近(20.9 mmol∙m−2a−1[63],但均低于东沙水合物发育区的值[32-33]

      自生碳酸盐岩的沉淀主要受孔隙水总碱度和Ca2+浓度的控制。由于孔隙水碱度随AOM作用强度的增加而增加,所以自生碳酸盐岩沉淀速率SH-W19-PC站位(9.3 mmol∙m−2a−1)高于SH-W23-PC站位(6.6 mmol∙m−2a−1),但略小于神狐海区的另一冷泉活动站位(10.43 mmol∙m−2a−1[62]

      通常来讲,正常沉积物孔隙水DIC的源主要有OSR、AOM和ME过程,DIC的汇主要包括自生碳酸盐岩沉淀、向海底扩散到进入海水中和沉积埋藏部分[63]。以上6个过程直接影响SH-W19-PC和SH-W23-PC两个站位的DIC库。对于DIC的来源,两个站位AOM作用产生的DIC占总DIC产量的比例分别为70.7%和60%,每个过程对DIC产量的贡献按AOM,OSR和ME的顺序依次减少。这与南海北部东沙和西沙海域相一致[59],但是该顺序又和郁陵盆地有所不同(AOM>ME>OSR),这可能是因为有机质活性和模拟深度不同所导致[64]。对于DIC的汇,SH-W19-PC和SH-W23-PC两个站位的DIC向海底排放的通量分别为24.7和20.7 mmol∙m−2a−1,均占DIC汇的约60%。另一部分DIC在沉积物中或在海底表层与Ca2+结合形成自生碳酸盐岩沉淀。研究区两个站位形成碳酸盐的DIC分别占总DIC汇的23%和20%,较小部分的剩余的DIC随沉积埋藏保存在孔隙水中。

    • 作为海洋碳循环的重要组成部分,沉积表层有机碳埋藏、有机质降解和甲烷氧化作用共同调节了深海碳库、沉积层碳库和孔隙水碳库之间的平衡,并持续影响海洋环境[65]。由河流带来和海表真光层形成的颗粒有机碳(POC)在沉降过程中被微生物降解,到达水深大于1 000 m海底的POC只占表层初级生产力的1%~3%[66],从神狐SH-W19-PC和SH-W23-PC两个站位的数值模拟结果看(图4),海底沉积表层POC通量为100 mmol·m−2a−1,在模型下边界POC的埋藏通量分别为85.6和84.6 mmol·m−2a−1,计算获得这两个站位有机质降解速率为14.4和15.4 mmol·m−2a−1表2)。在模型中这部分有机碳被OSR和产甲烷作用利用,这些原位产生的甲烷与下部来源的甲烷在SMTZ内消耗,产生大量DIC,进而扩散到上覆海水中。由于SH-W19-PC和SH-W23-PC站位甲烷活动强度较弱,甲烷都是以扩散形式向上运移,不存在甲烷气泡的运移方式,所以全部甲烷在近海底沉积物孔隙水中被消耗,进入海水中的甲烷通量接近0。通过对比两个站位碳转化关系可知,DIC是沉积表层碳循环的重要载体。随着甲烷通量的增加,由沉积物向海水中释放的DIC也随之增加,尽管计算结果比海岸带和大陆架的DIC通量小很多[67],但这些DIC可能对冷泉区底部水的化学性质产生一定影响。

      图  4  SH-W19-PC和SH-W23-PC站位沉积表层碳转化过程及模型定量计算结果示意图

      Figure 4.  Schematic diagram of quantitative calculation results of surface carbon conversion process and model of SH-W19-PC and SH-W23-PC sites

      根据DIC碳同位素的不同,通常可以对不同来源DIC进行划分。海水中的δ13CDIC为0,海洋真光层POC来源δ13CDIC约为−20‰,陆源有机质来源δ13CDIC不超过−25‰,而甲烷来源δ13CDIC均小于−40‰[68]。研究区两个站位δ13CDIC在柱样下部都达到−25‰附近,南海沉积有机质的δ13CDIC值约为−20‰[62,69],在模型中有机质降解过程包括OSR和ME两个过程,OSR过程不考虑碳同位素的分馏效应。ME的产物甲烷和DIC发生碳同位素分馏,但是这部分甲烷通过AOM作用产生的DIC又与之前的DIC混合,所以认为有机质降解来源的DIC碳同位素值与有机质一致[33]。那么神狐海区沉积物有机质来源的DIC和海水中的DIC混合不可能达到−25‰。通过数值模拟计算得出两个站位存在不同程度的AOM作用,可以认为两个站位孔隙水DIC是3个端元混合的结果。模拟结果显示AOM作用产生的DIC占总DIC产量的比例最大。从沉积物向海水中释放的DIC通量来看,即使所有的有机质来源的DIC全部被排放到海水中,在这两个站位也仅占沉积物—水界面DIC通量的47%和65%,所以甲烷来源的DIC占沉积物—水界面DIC通量的较大比例。本研究区两个站位甲烷来源的DIC及其亏损13C的特征所造成的潜在影响可能主要包括:①在沉积物中被化能自养生物利用,导致生物体碳同位素异常[70];②AOM作用产生的DIC导致孔隙水碱度升高,促进亏损13C的自生碳酸盐岩的形成[71-72];③向海水中排放亏损13C的DIC,可能对冷泉区域局部海底碳库的组分造成一定的影响。

      此外,在地质历史时期由海底沉积物向海洋排放的甲烷对于调节海洋碳库和全球气候变化起到了重要作用[61]。由本文的模拟研究结果可知,向海底排放的甲烷有一部分通过AOM作用以DIC形式进入了海洋,因此,推测当地质历史时期海底发生大规模甲烷渗漏活动时期,甲烷的最终归宿可能是以DIC形式大量进入海水中,从而深刻影响海洋无机碳库。

      DOC作为碳循环过程中的重要组成部分,目前对于冷泉区沉积物向海水中释放的DOC已有少量研究[73-75],研究表明随着海底甲烷通量的增加,以甲烷作为能源的微生物在代谢过程中将产生更多的DOC,在一些地区甲烷来源的DOC可达到沉积物向海水释放DOC通量的28%[9],这些DOC的活性成分可以被微生物利用进而转化为DIC或惰性DOC。由于缺乏实测DOC数据,无法对SH-W19-PC和SH-W23-PC两个站位孔隙水DOC的活性及通量进行判断,但是甲烷和有机质来源的DOC与DIC的转化是不能忽视的。因此,通过冷泉区甲烷来源的DOC的研究将有助于增强对深海DOC参与的生物地球化学循环的认识[9]

    • 本文利用反应运移模型对南海北部神狐海域SH-W19-PC和SH-W23-PC两个站位沉积物浅表层物理化学过程进行模拟,定量计算出碳循环各个过程的反应速率和转化量。计算结果显示两个站位平均AOM反应速率(24.3 mmol·m−2a−1)明显高于有机质降解速率(14.9 mmol·m−2a−1),OSR和ME过程产生的DIC量约为AOM作用产生DIC的一半。两个站位DIC汇的主要形式是由沉积物释放到海水中,沉积物—水界面DIC通量分别为(24.7和20.7 mmol·m−2a−1),少部分以自生碳酸盐岩的形式保存于沉积物中。因此,在研究区两个站位甲烷活动主导了沉积物浅表层的碳循环过程,大量甲烷来源的DIC释放到海水中。这些由沉积物向海水中排放的甲烷来源亏损13C的DIC,在地质历史时期海底甲烷活动强烈时期可能会深刻影响深海碳库组成和生态环境。

      致谢:非常感谢海洋四号机组人员在海上采样的帮助。

参考文献 (75)

目录

    /

    返回文章
    返回