二维温盐环流模式研究14C在海洋中的分布

海洋碳循环模式中常使用放射性同位素14C来检验该模式的物理模型是否较好地反映了海洋中的大尺度变化的过程。本工作使用的物理模型是包括太平洋、大西洋和南大洋在内的二维温盐环流模式，在给定的年平均海表强迫下积分4000年后，环流达到准稳态，其定常流场用来驱动14C模式。对14C模式积分5000年以上，可得到14C稳态的分布。使用不同的二氧化碳海气交换系数和垂直湍流扩散系数值进行了7个实验，发现使用Broecker等的交换系数方程和扩散系数值依深度从0.6增加到1.45 cm2/s，模拟的结果较好。使用观测到的大气中14C含量，进一步研究了工业革命后及核试验后的14C在海洋中的穿透，模拟结果与现有的一些资料进行了比较，总体上给出了合理的结果。     

全球海洋环流模式中氚分布的模拟

氚（3H）作为一种重要的被动示踪物,经常被用于研究海洋中的物理过程及评估海洋环流模式的模拟性能。使用一个全球海洋环流模式（LICOM）来研究氚在海洋中的分布、存储和输送。模拟的全球氚通量表明,1975年之前氚主要由海气交换输入海洋,特别是在1963年,氚的气体交换输入约为降水输入的2.5倍,1975年之后两种方式的氚输入通量都大幅减少。比对GEOSECS（Geochemical Ocean Sections Study,1972～1978年）和WOCE （World Ocean Circulation Experiment,1989～1995年）大洋观测计划期间的观测资料发现,我们的模式很好地模拟出了氚的海表分布、水柱总量、经向分布以及次表层的高值信号,主要缺点在于模拟的氚向深层的穿透不足,特别是在全球的两个副热带地区,表现尤为明显,氚输入函数的不确定性和模式物理场描述的不足可能是造成误差的主要原因。模式给出的海洋中氚储存总量的结果与基于观测得到的结果比较吻合,如北太平洋海区：1973～1974年模拟结果约为20.4 kg,相同期间观测估计值为21.1±4.7 kg,1989～1995年模拟结果为20.7 kg,相同期间观测估计值为23.4±2.0 kg。氚在等密度面上高低纬的侧向通风明显,模式成功模拟出氚从中高纬的海表进入,沿等密度面向低纬的次表层输送,又经大洋环流和扩散分别向南半球和高纬输送的过程。     

改进的UD模式及其在海洋碳循环中的应用

针对政府间气候变化专业委员会历次报告中"表层—深层"海水间碳通量的差异,利用改进的"上升—扩散"(upwelling-diffusion,UD)模式模拟海洋碳循环状态。除了改进模式结构外,还建立了直接计算海洋溶解无机碳(dissolved inorganic carbon,DIC)δ13C分布的控制方程。区别于目前广泛采用的深度参数的方法,直接利用沉积物捕获器的观测计算分解流量,减少了参数的引用。为更加符合海洋生物地球化学循环,采用地球化学海洋断面研究计划GEOSECS(Geochemical Ocean Sections Study)的海洋磷酸盐分布来标定两个深海动力参数:垂直扩散系数3 000 m2/a和上涌速率3.5 m/a,从而得到温跃层平均深度为860 m,这些参数都在直接观测和其他模式标定的结果之间。通过复原稳定状态下DIC分布并与GEOSECS的观测进行对比,得出自产业革命到20世纪70年代中期的海洋吸收大气CO2约为78 Gt C(Gigatonnes Carbon,1 Gt=1×1015g),与其他的模式估计及观测计算结果相吻合。模式模拟的"混合层—深层"海水间的碳通量为46 Gt C/a,远小于IPCC第四次报告(101 Gt C/a)而与第三次报告(42 Gt C/a)的估算相一致,但仍小于δ13C质量守恒法确定的范围(60~80 Gt C/a)。     

印度洋二维碳循环模式中表层CO2分压分布对物理和生化过程的敏感性试验

利用二维印度洋碳循环模式的模拟结果,集中对表层海洋中的CO2分压分布及其控制因子、海洋生物对海气CO2交换的影响、海洋营养物含量的改变和海洋环流的改变对大气CO2浓度的影响等进行了深入的分析和讨论,并与实际的GEOSECS观测数据的分析结果做比较;研究了与表层海洋CO2分压相关的海洋条件,较详细讨论了形成海洋表层CO2源与汇系统的决定因素及其相对重要性,得到了海洋热力因子和海洋环流对海洋表层的CO2化学过程起着决定性作用而生物过程仅处于次要地位的重要结论。此外,还利用建立的海洋碳模式进行了一些有意义的数值试验,详细讨论了海洋的物理化学因子改变对大气CO2浓度的可能影响。       

二维的大气CO2——大西洋碳循环模式

本文描述了一个二维（纬度×深度）的大西洋碳循环模式，模拟了大气和海洋间CO2的交换以及碳在海洋中的输送过程。模式在运行时使用了一个12层的三维动力学模拟的海洋环流的结果。大西洋被划分成397个网格箱，每个箱子中各种形式的碳的含量、总碱度、溶解的无机营养物和溶解氧的浓度以及几种14C(碳14)同位素的值分别得到求解。模式稳定状态的计算采用解大型稀疏线性方程组的直接解法。计算结果与“地球化学的海洋研究（GEOSECS）”的实际观测数据对比，表明模式较好地再现了实际大西洋中几种化学量的分布。     

海洋碳循环模式(Ⅱ)——对印度洋的模拟结果分析

把建好的海洋碳模式应用于印度洋区域，模拟得到了印度洋中与碳有关各化学量的表层分布、垂直分布和沿子午线面的等值线分布。与实测的GEOSECS(Geochemical Ocean－Section Study)数据作对比，模式较好地再现了印度洋上营养盐浓度、总碳浓度、总碱度和溶解氧的二维分布。通过模拟还发现，在稳定状态下，大气和海洋中总碳含量的分布依赖于发生在海洋中的各种物理化学过程及边界条件，水平扩散系数Kh和光合作用常数率Kg对各化学量的分布有较大影响(以前有学者认为不太重要，如Baes[1])；南印度洋中纬地区10°S至30°S是14C的重要向下渗透区域，人为排放的CO2可通过这片渗透区从海洋的表层输入海洋的深层。     

对印度洋偶极子中海洋环流异常的模拟研究

用高分辨率印度洋-太平洋区域海洋环流模式(IPOM)模拟研究印度洋偶极子(IOD)过程.用观测的1990年～1999年热带海表风应力强迫IPOM,模拟出20世纪90年代出现的两次(1994年和1997年)IOD过程中热带印度洋海温异常的一些基本特征.通过模拟的海洋环流过程,揭示出IOD过程中海洋环流异常的物理图像.发现在IOD事件时,东赤道印度洋上层出现强的向西(负)的距平流窄带,此距平流在赤道两侧向外辐散,且具有向西传的海洋Rossby波特征.IOD位相时在沿赤道的垂直剖面上,存在一个明显的距平环流圈:表层为强的向西距平流;下面为向东的补偿流;80°E以东存在着明显的涌升流,构成垂直环流圈的上升支;其下沉支主要在55°E以西的西印度洋.同时在热带东印度洋赤道两侧各有一个垂直的经向距平环流圈,其共同的上升支在赤道附近.在反IOD位相时,洋流距平分布与IOD位相截然相反,但洋流距平的绝对值较小.由上述距平洋流分布的特征发现,IOD过程中热带印度洋海温异常(东冷西暖)现象,可从水平和垂直海流的异常变化,特别是大范围异常涌升流和沉降流的出现得到解释.     

海洋对人为CO2吸收的三维模式研究

文中用包含海洋化学过程和一个简单生物过程的三维碳循环模式模拟了海洋对大气CO2 的吸收 ,并分析了碳吸收的纬度分布。模拟工业革命以来海洋对大气 CO2 的吸收表明 :海洋碳吸收再加上大气 CO2 的增加只占由化石燃料燃烧、森林砍伐和土地利用的变化而释放到大气中的 CO2 的 2 /3。1 980～ 1 989年期间海洋年平均吸收 2 .0 5Gt C。海洋人为 CO2 的吸收有明显的纬度特征。模式计算的海洋 CO2 的吸收在总量与纬度分布上与观测结果比较相符。     

太平洋大尺度环流数值模拟 Ⅲ：季节变化

本文在年平均环流数值模拟基础上,用季节变化的大气强迫场又积分了12年,并对由模拟得到的太平洋海面起伏、洋流、温度、海表热通量和上层热储存率等的季节变化特征作了较为详细的阐述和与已有观测结果进行了比较,事实证明,模式基本上模拟出了观测到的太平洋大尺度环流基本特征及其季节变化,证实基于理论考虑设计的IAP OGCM模式具有较好的性能和对实际环流的模拟能力,有些结果并优于国外一些具有同等分辨率的海洋模式。     

一个大洋环流模式和相应的海气耦合模式的评估Ⅰ.热带太平洋年平均状态

评估了中国科学院大气物理研究所大气科学和地球流体动力学数值模拟国家重点实验室海洋环流模式L30T63和海气耦合模式FGCM-0模拟的热带太平洋年平均状态,资料取自L30T63由观测的大气强迫驱动的Control试验、由NCAR CCM3大气强迫驱动的Spinup试验、以及相应的海气耦合模式FGCM-0.主要的结论是:(1)在"准确"的海表强迫下,Control模拟的海面温度和温跃层与观测结果相当接近,模式的固有误差是赤道冷舌过分西伸和东南太平洋温跃层偏浅.(2)Spinup能模拟出合理的热带太平洋上层海洋环流,但存在两个问题,即:暖池区海面温度显著偏高、沿赤道的梯度过大;赤道温跃层偏浅、东西向坡度偏小,它们分别与CCM3提供的海表短波辐射通量和风应力的系统误差有关.这两个问题很可能是海气耦合模式FGCM-0运行初期误差迅速发展的重要原因.(3)FGCM-0模拟的赤道暖池区上层100 m的平均温度比观测低3℃.分析表明FGCM-0夸大了暖池区海洋动力过程的降温作用,使得模拟的"暖池"在一定程度上具有冷舌的属性.FGCM-0模拟的热带南太平洋温跃层比观测结果偏浅数十米到100 m,以致赤道两侧的上层海洋温度分布趋于对称,成为"double ITCZ"现象在上层海洋中的表现.风应力旋度的系统误差和垂直混合随深度衰减过快是温跃层偏浅的两个可能原因;FGCM-0中与北太平洋中高纬地区深厚冷偏差相关的经圈环流也有利于热带温跃层误差的维持.     

Simulation of ~(14)C in IAP/LASG L30T63 Ocean Model

1. Introduction14C is a radioactive isotope of carbon, whose halftime is about 5730 yr. Under natural circumstance,14C comes into being in the stratosphere because of thenitrogen explosion of cosmic radial. 14C is mixed inthe atmosphere, absorbed by oceans, and then trans-ported into deep ocean. 14C radioactively reduces asits age increases. The reduction process of 14C canuncover the evolvement process of sea water's age andrenewal time.14 C content is usually described with one in athousan…     

A degradation approach to accelerate simulations to steady-state in a 3-D tracer transport model of the global ocean

 We have developed a new method to accelerate tracer simulations to steady-state in a 3-D global ocean model, run off-line. Using this technique, our simulations for natural ¹⁴C ran 17 times faster when compared to those made with the standard non-accelerated approach. For maximum acceleration we wish to initialize the model with tracer fields that are as close as possible to the final equilibrium solution. Our initial tracer fields were derived by judiciously constructing a much faster, lower-resolution (degraded), off-line model from advective and turbulent fields predicted from the parent on-line model, an ocean general circulation model (OGCM). No on-line version of the degraded model exists; it is based entirely on results from the parent OGCM. Degradation was made horizontally over sets of four adjacent grid-cell squares for each vertical layer of the parent model. However, final resolution did not suffer because as a second step, after allowing the degraded model to reach equilibrium, we used its tracer output to re-initialize the parent model (at the original resolution). After re-initialization, the parent model must then be integrated only to a few hundred years before reaching equilibrium. To validate our degradation-integration technique (DEGINT), we compared ¹⁴C results from runs with and without this approach. Differences are less than 10‰ throughout 98.5% of the ocean volume. Predicted natural ¹⁴C appears reasonable over most of the ocean. In the Atlantic, modeled Δ¹⁴C indicates that as observed, the North Atlantic Deep Water (NADW) fills the deep North Atlantic, and Antartic Intermediate Water (AAIW) infiltrates northward; conversely, simulated Antarctic Bottom Water (AABW) does not penetrate northward beyond the equator as it should. In the Pacific, in surface eastern equatorial waters, the model produces a north–south assymetry similar to that observed; other global ocean models do not, because their resolution is inadequate to resolve equatorial dynamics properly, particularly the intense equatorial undercurrent. The model’s oldest water in the deep Pacific (at −239‰) is close to that observed (−248‰), but is too deep. Surface waters in the Southern Ocean are too rich in natural ¹⁴C due to inadequacies in the OGCM’s thermohaline forcing. Received: 18 March 1997 / Accepted: 27 July 1997     

Transports and budgets of volume,heat, and salt from a global eddy-resolving ocean model

The results from an integration of a global ocean circulation model have been condensed into an analysis of the volume, heat, and salt transports among the major ocean basins. Transports are also broken down between the model's Ekman, thermocline, and deep layers. Overall, the model does well. Horizontal exchanges of mass, heat, and salt between ocean basins have reasonable values; and the volume of North Atlantic Deep Water (NADW) transport is in general agreement with what limited observations exist. On a global basis the zonally integrated meridional heat transport is poleward at all latitudes except for the latitude band 30°S to 45°S. This anomalous transport is most likely a signature of the model's inability to form Antarctic Intermediate (AAIW) and Antarctic bottom water (AABW) properly. Eddy heat transport is strong at the equator where its convergence heats the equatorial Pacific about twice as much as it heats the equatorial Atlantic. The greater heating in the Pacific suggests that mesoscale eddies may be a vital mechanism for warming and maintaining an upwelling portion of the global conveyor-belt circulation. The model's fresh water transport compares well with observations. However, in the Atlantic there is an excessive southward transport of fresh water due to the absence of the Mediterranean outflow and weak northward flow of AAIW. Eddies in the mid-latitudes act to redistribute heat and salt down the mean gradients. Residual fluxes calculated from a sum of the computed advective (including eddies), forced, and stored fluxes of heat and salt represent transport mostly due to vertical sub-grid scale mixing processes. Perhaps the model's greatest weakness is the lack of strong AAIW and AABW circulation cells. Accurate thermohaline forcing in the North Atlantic (based on numerous hydrographic observations) helps the model adequately produce NADW. In contrast, the southern ocean is an area of sparse observation. Better thermohaline observations in this area may be needed if models such as this are to produce the deep convection that will achieve more accurate simulations of the global 3-dimensional circulation.     

The impact of a downslope water-transport parametrization in a global ocean general circulation model

The impact of a downslope water-transport parametrization on the circulation and water mass characteristics of a global depth-level ocean general circulation model is investigated. The spreading of dense water from the formation regions into the deep ocean is known to be poorly represented in depth-level models with no bottom boundary layer resolved or attached. The new scheme is simple and intends to parametrize the effects of various oceanographic processes (rather than the processes themselves) that help dense water to descend topographic slopes by which the formation regions are separated from the world ocean. The new scheme significantly improves the large scale properties of the North Atlantic Deep Water. Changes in the North Atlantic circulation, however, are rather small. In the Southern Ocean, the exchange between the dense water formation regions on the continental shelves and the deep ocean is strengthened at the expense of deep water mass formation by open ocean convection. In all three ocean basins, the density of the deep and bottom water is higher with the new parametrization, which brings the simulations closer to observations in the Atlantic and Indian Oceans. In the Pacific Ocean, however, where the density has already been well reproduced without the downslope transport, it becomes slightly too high. The results are in agreement with those from other model studies.     

Stability dependent parameterisations of vertical mixing in ocean general circulation models

Summary Parameterisations of mixing induced through shear instability, internal wave breaking, and double diffusion are investigated in simulations of ocean climate using a global ocean general circulation model (OGCM). Focus is placed on the sensitivity of the large scale circulation, water mass formation and transport of heat as measures of the model's ability to represent current climate. The model resolution is typical of OGCMs being coupled to atmospheric. GCMs in climate models and the parameterisations investigated are all computationally inexpensive enough to allow for integrations on long time scales. Under the assumption of constant vertical eddy coefficients (the control case), the model climatology displays acceptable values of North Atlantic Deep Water formation, Antarctic Circumpolar Current (ACC) transport, and Indonesian through-flow but an excessively deep and diffuse pycnocline structure with weak stratification in the deep ocean. It is found that various circulation and water mass properties are sensitive to the choice of parameterisation of vertical mixing and that determining a scheme which works satisfactorily over all regions (tropical, mid-latitude, and polar) of the domain is not straightforward. Parameterisations of internal wave breaking or upper ocean shear instability lead to some improvements in the model water mass formation. ACC and poleward heat transport when compared to the control case whereas parameterisations of double diffusive processes did not. Based on these and other results, various recommendations are made for mixing parameterisations in ocean climate models.With 8 Figures     

The influence of the Bering Strait on the circulation in a coarse resolution global ocean model

An ocean general circulation model of global domain, full continental geometry and bottom topography, is used to study the influence of the Bering Strait on the general circulation by comparing equilibrium solutions obtained with and without a land-bridge between Siberia and Alaska. The model is integrated with restoring boundary conditions (BC) on temperature and salinity, and later, with mixed BC in which a restoring BC on temperature is maintained but a specified flux condition on salinity is imposed. In both cases, the effect of the Bering Strait is to allow a flow of about 1.25–1.5 Sv from the North Pacific to the Arctic Ocean and, ultimately, back to the North Pacific along the western boundary current regions of the Atlantic and Indian Oceans. When a restoring BC on salinity is used, the overturning associated with North Atlantic Deep Water and Antarctic Intermediate Water formation are increased if the Bering Strait is present in the model geometry. The result of switching to a specified flux BC on salinity is to cause a transition in the THC in which the overturning associated with North Atlantic Deep Water formation increases from about 12 Sv to about 22 Sv. This transition occurs in an essentially smooth fashion with no significant variability and is about 12% smaller in magnitude if the Bering Strait is present in the model geometry. Because the Bering Strait appears to exert some influence on the general circulation and the formation of deep water masses, it is recommended that this Strait be included in the geometry of similar resolution models designed to study the deep ocean and potential changes in climate. Correspondence to: CJC Reason     

A carbon cycle model with latitude dependence

A two-dimensional carbon cycle model is divided into three zones representing equatorial, middle and high latitude regions. The three zones are coupled together by a deep ocean meridional convective cell and atmospheric transport terms. The model is applied to the calculation of the dispersion of radiocarbon and tritium from nuclear weapons tests, to the calculation of the atmospheric record of bomb radiocarbon and to the calculation of the Mauna Loa record of atmospheric CO₂. Calibrating on the basis of the Northern hemisphere bomb test data yields a model which has approximately twice the CO₂ ocean uptake of the one-dimension box diffusion models calibrated on the basis of deep water equilibrium carbon 14.Work performed under the auspices of the U.S. Department of Energy by the Lawrence Livermore Laboratory under contract No. W-7405-ENG-78.     

Formation and pathways of North Atlantic Deep Water in a coupled ice–ocean model of the Arctic–North Atlantic Oceans

We investigate the formation process and pathways of deep water masses in a coupled ice–ocean model of the Arctic and North Atlantic Oceans. The intent is to determine the relative roles of these water masses from the different source regions (Arctic Ocean, Nordic Seas, and Subpolar Atlantic) in the meridional overturning circulation. The model exhibits significant decadal variability in the deep western boundary current and the overturning circulation. We use detailed diagnostics to understand the process of water mass formation in the model and the resulting effects on the North Atlantic overturning circulation. Particular emphasis is given to the multiple sources of North Atlantic Deep Water, the dominant deep water masses of the world ocean. The correct balance of Labrador Sea, Greenland Sea and Norwegian Sea sources is difficult to achieve in climate models, owing to small-scale sinking and convection processes. The global overturning circulation is described as a function of potential temperature and salinity, which more clearly signifies dynamical processes and clarifies resolution problems inherent to the high latitude oceans. We find that fluxes of deep water masses through various passages in the model are higher than observed estimates. Despite the excessive volume flux, the Nordic Seas overflow waters are diluted by strong mixing and enter the Labrador Sea at a lighter density. Through strong subpolar convection, these waters along with other North Atlantic water masses are converted into the densest waters [similar density to Antarctic Bottom Water (AABW)] in the North Atlantic. We describe the diminished role of salinity in the Labrador Sea, where a shortage of buoyant surface water (or excess of high salinity water) leads to overly strong convection. The result is that the Atlantic overturning circulation in the model is very sensitive to the surface heat flux in the Labrador Sea and hence is correlated with the North Atlantic Oscillation. As strong subpolar convection is found in other models, we discuss broader implications.