全球海洋碳循环模式模拟研究.docxVIP

全球海洋碳循环模式模拟研究 1 海洋碳循环模式 工业以来，婚煤气中的coa浓度继续增加，从1800年的28010-6增加到2005年的37910-6，增长近10010-6，并有连续增长趋势。大气CO2浓度不断增长引起的温室效应是全球变暖的重要原因之一。由全球变暖引发的一系列环境、生态问题,直接威胁到人类社会与经济的可持续发展。为了应对CO2,CH4等温室气体诱发的全球环境问题,有关碳循环的生物地球化学研究已引起了各国科学家和政府的关注,并在全球范围内开展了以国际地圈生物圈计划(IGBP)、国际人文因素计划(IHDP)、世界气候研究计划(WCRP)和生物多样性计划(DIVERSITES)为核心的国际合作研究。 占地球表面71%的海洋是一个巨大的碳库,其碳贮量约为38 000Gt,约是大气碳库的60倍,陆地生物土壤碳库的20倍,海洋可吸收大量的人为碳排放,从而减缓气候的变化。因而,海洋碳循环在全球碳循环中具有极其重要的作用,海洋碳循环研究对全球气候变化研究具有重要意义。然而,由于海洋观测的时空局限性,很难覆盖整个研究区域,长期观测资料更显匮乏,这限制了人类对海洋碳循环的直接认识,但也促进了海洋碳循环模式的不断发展。海洋碳循环模式一方面补充了观测的不足,另一方面也为我们预测未来海洋碳循环的变化趋势及其在全球气候变化中的作用提供了途径。在过去的几十年中,海洋碳循环模式的发展经历了由简单到复杂的发展过程,由最初的箱式模型,发展到现在基于海洋动力学的三维海洋碳循环模式,海洋碳循环研究取得了很大的进展。 海洋碳循环模式强烈的依赖海洋环流模式,环流模式的不断发展,反过来也会改进碳循环模式。本文基于海洋环流模式POP[10—12]和生物地球化学模型OCMIP-2[13—14]建立了一个全球海洋碳循环模式,对全球海洋碳循环进行了模拟,分析了全球海洋对人为碳吸收,以及人为碳在海洋中的分布。 2 类型介绍 2.1 海洋生物地球化学 本文所采用的海洋碳循环模式基于海洋碳循环模式比对计划第二阶段的生物地球化学模型OC-MIP-2(Ocean Carbon Model Inter-comparison Project phase 2)[13—14]。CO2在海气界面通过气体交换过程进入或离开海洋,其海气通量的计算为 其中,k为液相传输速度,是风速和海表面温度(SST)的函数;β是CO2的溶解度,是SST和海表面盐度(SSS)的函数;pCO2atm和pCO2ocn分别是大气底层和海洋表层的CO2分压,fice是海冰覆盖率。pCO2ocn是由表层海水的CO2体系确定,表层海水的CO2体系受表层溶解无机碳(DIC)、总碱度(ALK)、SST、SSS等影响。如果大气CO2分压大于海洋CO2分压,海洋将吸收CO2;反之,海洋会释放CO2。 海洋中碳的生物地球化学过程包括:浮游植物通过光合作用吸收溶解在海水中的CO2,将其转化为生物体内的有机物;通过食物链、生命过程等将碳以有机物的形态向深海输送,并在深海中部分被矿化再转化成无机碳,以生物泵形式实现碳从表层向深层的输运。OCMIP-2模型将海洋分成上层生产区和下层消耗区,其分界线补偿深度zc取为75m。上层生物的生产将CO2转化为有机碳,有机碳在下沉过程中将碳带入深层海洋,并在深层海洋发生矿化,又转化为溶解无机碳。根据Doney等,生产力计算方案由营养盐驱动方案替代原来的营养盐恢复方案,考虑了温度、营养盐、太阳辐射强度、生物量以及混合层深度的共同作用;并在生物地球化学循环中引入了铁循环,将铁的限制作用考虑到生产力计算方案中。海洋上层的生产力计算公式为: 其中,FT是温度限制因子;FN是营养盐限制因子,包括了磷酸盐和铁的共同限制作用;FI是光限制因子;B是生物量,是磷酸盐和铁的浓度共同决定的;zc为补偿深度;zml是混合层深度,该项的作用是在混合层深度大于补偿深度zc的海区,将补偿浓度等价到混合层深度;τ为生物生产的时间尺度,模式中取15d。 考虑碳的生物地球化学过程后,海洋中碳的控制方程为 其中,L(C)为平流项;HDiff(C)和VDiff(C)分别为水平和垂向扩散项;JC为生物地球化学过程产生的源汇项。 2.2 海洋模式的积分 本文所用海洋环流模式POP(Parallel Ocean Program)[10—12]是美国Los Alamos国家实验室发展的用于高性能计算的海洋环流模式,被广泛应用于海洋研究,不仅用于全球和区域的海洋环流的模拟[17—18],还作为气候耦合模式CCSM3.0中的海洋模式用于气候研究。 由于海洋碳循环达到工业革命前平衡态的积分时间较长,我们的海洋模式采用较低分辨率网格,其水平分辨率在经向上为0.8°~1.8°,纬向上约为3.6°,垂向为25层。在气候态大气场强迫下,大气CO2浓度为工业