任意正交曲线坐标系下的海洋模式动力框架的发展与评估

本文发展了一个可以适用于任意水平正交曲线坐标系的海洋模式动力框架,并将其应用于中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室发展的气候系统海洋模式LICOM2.0（LASG/IAP Climate system Ocean Model,version2.0）。在经纬网格坐标系下,新的动力框架与LICOM2.0原有的动力框架模拟结果完全一致。基于新的动力框架,海洋模式可采用能够准确描述北冰洋地形的三极网格,克服了LICOM2.0经纬网格版本必须将北极点处理为孤岛的缺陷,从而显著改进了模式对于北冰洋环流和北大西洋经圈翻转流函数（AMOC）的模拟能力。此外,引进三极网格还可以避免模式网格距随纬度增加而急剧减小带来的计算不稳定,在LICOM2.0的三极网格版本中,模式不需要采用任何空间滤波方案仍然能够保证计算的稳定性,从而与LICOM2.0的经纬网格版本相比,极大地提高了模式的并行效率,这一点在当水平分辨率提高到0.1度时表现得尤为明显,海洋模式的并行加速比可以从经纬网格版本的5.8左右提高到三极网格版本的15.0左右。     

LICOM海洋模式对南大洋环流的模拟评估

使用中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室（LASG/IAP）开发的第三代气候海洋模式（LASG/IAP Climate system Ocean Model version 3, LICOM3.0）低分辨率版本在海洋模式比较计划（Ocean Model Intercomparison Project, OMIP）试验中的模拟数据,描述了南极绕极流（Antarctic Circumpolar Current, ACC）和南大洋经向翻转环流（Meridional Overturning Circulation, MOC）在1958～2009年的平均状态及其变化,并与已有的模式模拟结果和观测结果对比以评估LICOM模式的模拟效果。通过对比已有模式模拟数据发现,LICOM3.0模式模拟的ACC和南大洋MOC在两组OMIP试验中平均状态相仿、结果在合理范围内,但OMIP1试验中海表强迫的变化趋势较OMIP2试验中的变化更大,得到的环流输送在OMIP1试验中增长趋势也更大。     

Sensitivity of the Upper Ocean Temperature and Circulation in the Equatorial Pacific to Solar Radiation Penetration Due to Phytoplankton

Solar radiation penetration in the upper ocean is strongly modulated by phytoplankton, which impacts the upper ocean temperature structure, especially in the regions abundant with phytoplankton. In the paper, a new solar radiation penetration scheme, based on the concentration of chlorophyll-a, was introduced into the LASG/IAP （State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics/Institute of Atmospheric Physics） Climate system Ocean Model （LICOM）. By comparing the simulations using this new scheme with those using the old scheme that included the constant e-folding attenuation depths in LICOM, it was found that the sea surface temperature （SST） and circulation in the central and eastern equatorial Pacific were both sensitive to the amount of phytoplankton present. Distinct from other regions, the increase of chlorophyll-a concentration would lead to SST decrease in the central and eastern equatorial Pacific. The higher chlorophyll-a concentration at the equator in comparison to the off-equator regions can enlarge the subsurface temperature gradient, which in turn strengthens the upper current near the equator and induces an enhancing upwelling. The enhancing upwelling can then lead to a decrease in the SST in the central and eastern equatorial Pacific. The results of these two sensitive experiments testify to the fact that the meridional gradient in the chlorophyll-a concentration can result in an enhancement in the upper current and a decrease in the SST, along with the observation that a high chlorophyll-a concentration at the equator is one of the predominant reasons leading to a decrease in the SST. This study points out that these results can be qualitatively different simply because of the choice of the solar radiation penetration schemes for comparison. This can help explain previously reported contradictory conclusions.     

涡分辨率模式模拟的南海中尺度涡旋

Mesoscale eddies(MEs) in the South China Sea(SCS) simulated by a quasi-global eddy-resolving ocean general circulation model are evaluated against satellite data during 1993–2007. The modeled ocean data show more activity than shown by the satellite data and reproduces more eddies in the SCS. A total of 345(428) cyclonic eddies(CEs) and 330(371) anti-cyclonic eddies(AEs) generated for satellite(model) data are identified during the study period, showing increase of ~24% and ~12% for the model data, respectively. Compared with eddies in satellite, the simulated eddies tend to have smaller radii, larger amplitudes, a slightly longer lifetime, faster movement and rotation speed, a slightly larger nonlinear properties(U/c) in the model. However, the spatial distribution of generated eddies appears to be inhomogeneous, with more CEs in the northern part of SCS and fewer AEs in the southern part. This is attributed to the exaggerated Kuroshio intrusion in the model because the small islands in the Luzon Strait are still not well resolved although the horizontal resolution reaches(1/10)°. The seasonal variability in the number and the amplitude of eddies generated is also investigated.     

FGOALS高分辨率气候模式系统模式研制与应用综述

当今气候系统模式发展的重要趋势之一,是通过提高模式的空间分辨率,改进对气候系统中多尺度相互作用过程和极端事件的模拟能力。过去5年里,中国科学院大气物理研究所发展并完善了25 km分辨率大气环流分量模式FAMIL2.2、1/10°分辨率海洋环流分量模式LICOM3.0,并以此为基础建立了高分辨率气候系统模式FGOALS-f3-H。利用上述高分辨率模式,开展了大量的数值模拟试验和预报/预测研究,其中包括国际耦合模式比较计划第六阶段（CMIP6）的高分辨率模式比较子计划（HighResMIP）,建立了海洋环流预测系统（LFS）等。初步评估分析表明,相对于低分辨率模式,高分辨率模式对气候平均态和气候变率的模拟能力均有明显改进。其中高分辨率大气环流模式可以更好地模拟台风、极端降水事件,高分辨率海洋模式可以更好地模拟海洋中尺度涡旋和西边界流,而高分辨率耦合模式则可以更好重现中尺度海气相互作用过程、热带不稳定波动（TIW）等事件。     

黑潮入侵优化对南海北部中尺度涡旋模拟的影响

基于高分辨率海洋环流模式，通过比较吕宋海峡处地形优化后的黑潮入侵形态和强度不同的试验，我们研究了黑潮入侵优化后对南海中尺度涡模拟的影响。我们发现黑潮入侵的减弱导致了涡旋活动的减弱，这使得模式结果与观测结果更为相近。在这种情况下，模式模拟的吕宋海峡西部及北部陆坡区域的涡动动能明显减弱。模式涡动动能的减弱与模式反气旋式涡数量的减少和气旋式涡强度的减弱有关。涡动动能收支的分析进一步表明，黑潮入侵的优化将通过改变水平速度切变和温跃层斜率来改变涡动动能，而这两个参数分别与正压和斜压不稳定性有关。前者在模式涡动动能减弱中起着更为重要的作用，而黑潮入侵导致的涡动动能的水平输送对吕宋海峡西部区域的能量收支同样起着重要的作用。     

一个高分辨率太平洋-印度洋海盆环流模式的初步结果

利用LASG/IAP发展的一个0.25°×0.25°高分辨率太平洋-印度洋海盆环流模式，初步分析了模式在太平洋区域的模拟结果，并与海洋同化资料以及前人的研究结果作比较，检验此模式对该区域平均气候态、年际变化的模拟能力。分析表明，模式较好地再现了海表温度（SST）分布、赤道温跃层和纬向流结构、赤道流系分布形态、海表高度以及正压流函数空间分布特征；同时，对显著的El Ni?o和La Ni?a事件的模拟等方面与Simple Ocean Data Assimilation（SODA）2.0.2版本结果相近。此外，模式模拟北赤道流（NEC）分叉点位置的季节和年际变化以及吕宋海峡流量的年际变化与已有研究结果基本一致。进一步分析还发现，在年际尺度上，NEC分叉点位置和吕宋海峡流量与ENSO密切相关。     

时间滑动相关方法在SST可预报性及可信计算时间研究中的应用

将时间滑动相关方法 STC(sliding temporal correlation)用于研究混沌系统和海洋环流模式的可信计算时间RCT(reliable computation time),Lorenz混沌系统的数值试验表明用STC求得的可预报时间和可信计算时间,与使用传统误差限方法所得结果一致,证明了其有效性。对海洋环流模式LICOM和NEMO的研究发现:1.当海洋模式以非耦合的方式运行时,试验的结果表明其海表温度SST的可信计算时间较长,平均达到6个月以上,这主要是由于海洋模式的运行过程中,采用恢复性边界条件使模拟结果不会太过偏离观测值。对于强迫场从1月开始的试验,LICOM模式的SST可信计算时间在赤道东太平洋和西北太平洋地区存在RCT低值区,其数值不超过2个月。而NEMO模式在赤道太平洋地区全是RCT高值区,NEMO模式的RCT低值区域出现在赤道外的太平洋和大西洋中纬度地区,强迫场从7月开始的试验,RCT纬向平均分布与1月有相反的形式。2.海洋模式以耦合方式运行时,由于去掉了恢复边界条件作用,海洋模式预报的SST可信计算时间明显减小,年平均RCT为1个月左右。按季节平均得到的RCT变化不大,在30~40天之间,RCT的大值区春季位于南半球,而秋季位于北半球,可达2个月以上。耦合模式中所模拟的500 h Pa高度场的RCT与单独运行的大气模式所得结果相差不大,仍在2周以内。3.无论是按季节平均还是按海区平均所得到的RCT分布,都在30~60天左右,只有极少数区域在特定季节可以达到80天以上,这说明在海气耦合模式中,由于计算不确定造成的可预报上限一般不超过2~3个月,这比使用资料分析得到可预报期限短很多,因此根据木桶原理,RCT可能是制约海气耦合模式SST预报能力的一个重要因素。     

南海中尺度涡的时空变化规律Ⅰ.统计特征分析

利用TOPEX/Poseidon混合ERS1/2高度计资料对1993～2002年间南海中尺度涡进行辨认和动态追踪,并对其时空变化规律进行统计分析,结果表明:南海平均每年出现18个涡旋,出现数量的年际变化与El Niño/La Niña有一定关系.多数涡旋的生命周期在180 d以内,近半数为30～60 d;半径大致在100 km到250 km间,其经向变化与斜压罗斯贝变形半径的经向变化趋势一致;80%的涡旋向西移动,纬向移速大致为在-8 cm/s到3 cm/s间,随纬度变化呈“∑”型分布.涡旋发生的区域主要位于越南南部以东至台湾西南一线海域,呈东北-西南向分布,其中,吕宋海峡以西海域和越南南部以东海域涡旋的出现概率相对较大,约为23%和25%,是涡旋的高发区.涡旋的存在对水深大于200 m海域的海面高度变化的平均贡献约为36%,在涡旋高发区的贡献高达80%.从涡旋出现地点、传播路径和速度、出现概率及对海面高度变化的贡献综合来看,涡旋大体可以沿17°N分成南北两个相对独立系统,一般不会有涡旋跨系统传播.     

CMIP5多模式模拟两类El Niño海表盐度分布及与降水的关系

利用CMIP5提供的25个工业革命前控制试验（piControl）模拟数据评估了热带太平洋两类El Ni?o（即东部EP和中部CP型El Ni?o）的海表盐度（SSS）空间结构差异及其与海表温度（SST）和降水的关系。结果表明:（1）大部分模式能够模拟出EP和CP型空间结构,两类El Ni?o中的SST、降水和SSS的空间技巧评分依次减小,其中,EP型SST和降水水平分布的模拟能力强于CP型,SSS则为CP型强于EP型,CP型模拟的SST、SSS和降水异常中心位置较EP型偏西且强度偏弱;（2）CP型SST、降水和SSS三者空间分布的线性一致性比EP型好,即在CP型中,SST影响降水,进而影响SSS,同时SSS对SST调制的反馈机制较显著,而对于EP型,由于海洋水平平流和非局地效应等因素,使得SST与SSS空间对应较差;（3）依据多模式模拟的SSS空间技巧评分高低将CMIP5模式分为两类,技巧评分低（高）的模式模拟的SST、SSS和降水异常值的中心位置偏西（偏东）,引起中心位置偏移的原因与模式模拟赤道太平洋冷舌的位置有关,即赤道太平洋冷舌西伸显著,导致发生El Ni?o时SST异常变暖西伸显著,进而使得降水异常和SSS异常位置偏西。同时,技巧评分低的模式还易出现向东南延伸的负SSS异常,原因是双赤道辐合带的东南分支过于明显,即降水偏多,导致SSS偏淡。SSS变化会影响ENSO的发生发展。因此,探讨两类El Ni?o盐度分布的差异及相关物理场的关系,为提高模式的气候模拟和预测提供有益的借鉴。