海洋模式中Boussinesq近似误差讨论

根据全球增暖的特点,设计了一个理想的数值试验方案,用Boussinesq POP海洋模式和改进的非Boussinesq POP海洋模式定量讨论了Bousiinesq近似在海洋模式计算中的误差。结果发现,在只有热力驱动的热力环流背景下,由热膨胀引起的海平面上升在水平方向上是基本均匀的,在所给的初始边界条件下,这种由Boussinesq近似引起的最大海平面误差可以达到59％,在Boussinesq POP模式中,热源中心处的海面高度要远小于由非Boussinesq模式计算的海面高度,而其周围有虚假的海面高度下降;在只有加热引起的热盐环流过程中,当模式作了Boussinesq假设以后,计算的经向和纬向垂直环流都会产生虚假的加强,虽然这种误差只是在1％左右;在Bousiinesq近似假定下,热量经向通量在赤道上垂直剖面的积分误差比质量经向通量在赤道上垂直剖面的积分误差大一个量级;非Boussinesq模式计算的气压梯度所做功的垂向分布在3000m以下是有波动的,而Boussinesq模式计算的气压梯度所做功的垂向分布在3000m以下基本上是均匀的,它的误差在10％以上。     

一个斜压海洋环流数值模式及其数值试验

为研究大尺度海洋环流和一些大尺度海洋物理过程,我们建立了一个无海底地形的三层斜压海洋环流数值模式。模式方程采用了静力近似和Boussinesq近似,在上边界采用了钢盖近似。同时,还做了两个数值试验:第一个模拟试验,海洋从静止状态开始,以4个季节的气候平均风应力和海表大气季节平均温度为强迫条件,积分区城选择45℃以北的太平洋,数值积分了5个模式年,这时海洋的上层环流和温度的年变化趋于平稳。计算结果再现了太平洋气候平均状态的海温分布和主要海流的大尺度特征以及它们的年变化特征。如西部强化流、赤道潜流等。第二个数值试验,在第一个数值模拟的基础上给定赤道西太平洋一个风应力异常区,数值积分了几个月。从计算结果中可看到类似E1 Ni■o的现象,并且初步地分析了赤道波动在E1 Ni■o形成过程中的作用。     

Stokes漂流近似公式对海洋表层流场估算的影响

本文采用波浪订正的Ekman模型,研究分析了三种Stokes漂流近似公式(单波公式、e指数公式、Phillips谱近似公式)对海洋表层流场估算的影响。海表总流场由海表面高度(SSH)数据计算的地转流和海浪模式WAVEWATCH Ⅲ输出结果计算的非地转流组成,并采用拉格朗日浮标观测数据对计算结果进行了验证。研究表明,随着Stokes漂流近似公式精度的提高,其计算的拉格朗日流速更接近于谱积分公式的计算结果,更贴近拉格朗日浮标观测数据。与谱积分公式计算的海表拉格朗日流速相比,单波公式的平均相对偏差为0.0834,e指数公式的平均相对偏差为0.0392,Phillips谱近似公式的平均相对偏差为0.0101,说明Phillips谱近似公式在不同风速下均能对谱积分公式有良好的近似效果。在低风速条件下,由Stokes漂流近似公式精度引起的海洋表层流场估算误差可以忽略不计,但随着风速增加,由近似公式精度引起的偏差逐渐变大,此时应该选择Phillips谱近似公式计算Stokes漂流,来减小误差。     

大洋环流模式中Boussinesq近似与静压近似问题的讨论

改进环流模式中的Boussinesq近似与静压近似是当前海洋环流模式的主要研究方向之一。虽然目前已提出一些在模式中包含非Boussinesq效应和静压效应的方法,但是仍有很多关键问题亟待解决。系统地回顾了改进这2个基本近似的理论分析和数值模式的研究进展,介绍了不同方法的理论基础及其数值分析结果,并分析了这些方法中所存在的缺陷和问题,指出了开展进一步研究需要首先分析的几个关键问题。     

海洋对热带气旋响应的研究 Ⅱ.不同海洋热力结构下的情形

利用中二层非线性原始方程海洋模式，研究海洋在自身不同热力结构下对热带气旋的响应。计算结果表明，初始混合层（ML）深度和层结强度对海表温（SST）和ML深度变化起着十分重要的作用。初始ML深度对海流量值影响较大，层结强度则较小。东海陆架区特殊的海洋热力结构，极易造成SST下降。海洋对7002号台风响应的模拟结果与观测资料较一致。     

海洋对热带气旋响应的研究

利用中二层非线性原始方程海洋模式，研究海洋在自身不同热力结构下对热带气旋的响应。计算结果表明，初始混合层（ＭＬ）深度和层强强度对海表温（ＳＳＴ）和ＭＬ深度变化起着十分重要的作用。初始ＭＬ深度对海流量值影响较大，层结强度则较小。东海陆架区特殊的海洋热力结构，极易造成ＳＳＴ下降。海洋对７００２号台风响应的模拟结果与观测资料较一致。     

混合坐标模式HYCOM模拟COADS强迫下的南海平均环流

采用混合坐标模式(HYCOM)模拟南海环流,同时利用海表温度卫星资料和吕宋海峡Sb-ADCP观测海流数据来评估模式结果.地形敏感性实验表明,吕宋海峡地形数据对模拟黑潮入侵方式影响较大,在地形误差较大的情况下,模拟的黑潮可能以反气旋流套方式入侵.和Pathfinder海表温度卫星资料比较,模式输出的月平均温度在海盆区域误差较小.ERA-15资料强迫所模拟的吕宋海峡上层环流和Sb-ADCP观测一致,而COADS结果低估了吕宋海峡的体积输送.     

渤海海冰漂移数值研究

本文利用海冰－海洋动力耦合模式对渤海典型天气形势下海冰动力过程作了数值模拟，海冰模式建立在动量和质量守恒基础上，忽略了海冰变化的热力过程。海冰厚度被划分为三类：堆积冰、平整冰和开阔冰，冰的形变由一个厚度的重新分布约束条件确定。海洋模式是一个二维风暴潮模式．同时考虑了潮汐的作用。风场资料来自于沿岸气象观测站每日四次的风观测，计算网格是十分之一经纬度。主要的分析和模拟是在大气和海洋共同作用下海冰的漂流特点。观测结果比较表明，数值模拟结果基本上反映了该海域流冰的漂移性质，同时也可为短期冰情预报提供有益的参考。     

一维海洋混合层模式应用于二维数值模拟的试验研究

基于一个一维湍能海洋混合层模式,发展了一个易于与大气或海洋模式耦合的,可应用于水平二维空间的模块化海洋混合层模式。并加入垂直上翻参数化方案,引入日平均海面温度及气候态松弛项,进行二维海洋模拟理想气旋实验和飓风实例模拟实验。理想气旋强迫实验表明,垂直上翻过程的参数化方案可有效地弥补原有一维海洋混合层模式无法形成气旋中心动力上翻的不足,消除了虚假暖水核心。卡特里娜飓风个例实验的结果表明,加入垂直上翻后,飓风中心附近的海表面温度误差明显减小。海洋混合层模式对海洋表面温度日变化模拟能力在二维应用中依然表现良好。经过上述改进,发展的海洋混合层模式可以较为真实地模拟平常海表温度高频日变化的同时,还对剧烈天气过程也有一定的模拟能力,具有广泛的应用前景。     

渤海海水漂移数值研究

本文利用海冰－海洋动力耦合模式对渤海典型天气形势下海冰动力过程作了数值模拟，海冰模式建立在动量和质量守恒基础上，忽略了海冰变化的热力过程。海冰厚度补划分为三类：堆积冰、平整冰和开阔冰，冰的形变由一个厚度的重新分布约束条件确定，海洋模式是一个二维风暴潮模式，同时考虑了潮汐的作用。风场资料来自于沿岸气象观测站每日四次的风观测。计算网络是十分之一经纬度。主要的分析和模拟是在大气和海洋共同作用下海洋的漂流     

Validation of a decadal OGCM simulation for the tropical Pacific

An ocean general circulation model is forced with NCEP reanalysis over the 1948–1999 period. The simulated dynamic height and sea level are compared respectively to the dynamic height computed from hydrological data and to the sea level measured by tide gauges in the tropical Pacific. The model is shown to capture important features of the temporal structure of variability in the tide gauge data over the last several decades. However, the comparison reveals a largely artificial trend in the simulation, which consists of a decreasing dynamic height and sea level in the southwest and northwest of the tropical Pacific. Model sensitivity experiments show this trend is controlled by the NCEP surface wind stresses and more precisely by a weakening in the trade winds and a trend in the off-equatorial wind curl, with this trend existing mainly before the mid 1970s. For studies of decadal variability, the simple removal of a linear trend is an inadequate way to solve the problem, due to the inhomogeneities in the data used in reanalysis products and the non-linearity of models.     

西北太平洋中国沿海海面倾斜及其机理初探

海面倾斜与高程基准密切相关,它已经受到大地测量学界和海洋学界的重视.进一步讨论了大地水准的测量精度,指出了大地(几何)水准与沿海验潮资料的不符是由于海面的倾斜;分析了海面倾斜的机理,其中包括海流、海水密度分布、气压、台风引起的破碎波等的作用;提出了用海水异常密度的三维分布计算大地水准面的扰动,实质上它反映了海面的倾斜,该量占大地水准测量结果的73%,这进一步说明密度的异常乃是我国沿海海面倾斜的主要原因.     

Comparison of sea surface heights observed by TOPEX altimeter with sea level data at Chichijima

The sea surface heights (SSHs) observed by the TOPEX altimeter are compared with tide gauge data at Chichijima in Ogasawara (Bonin) Islands and hydrographic data taken around the islands, in order to quantitatively verify the altimeter observations and oceanic tide corrections by three tide models proposed by Cartwright and Ray (1991), Rayet al. (1994), and Maet al. (1994). First, performance of the new tide models is assessed by comparing tidal variations consisting of diurnal and semi-diurnal constituents with the tide gauge data at Chichijima. The tide model proposed by Rayet al. gives the smallest root-mean-squared (rms) difference of 2.61 cm. Errors in amplitude and phase in each tide model are evaluated by spectral analysis. The TOPEX SSHs corrected by the tide models are compared with sea level data at Chichijima. A long-term variation of a period of about 1 year is found in the residual between the SSHs and the Chichijima sea levels. This variation is also found in the difference between the dynamic height anomalies calculated from hydrographic data around the island and the Chichijima sea levels. By subtracting the long-term variation, the rms difference between the TOPEX SSHs and the Chichijima sea levels is reduced to about 4 cm and the slope of the regression line is improved to unity. The residual shows variations related to aliasing caused by incompleteness of the ocean tide correction with the repeat cycle of the altimeter observation.     

Accuracy Assessment of Jason-1 and TOPEX/POSEIDON Along-Track Sea Surface Slope Special Issue: Jason-1 Calibration/Validation

Sea surface slope computed from along-track Jason-1 and TOPEX/POSEIDON (T/P) altimeter data at ocean mesoscale wavelengths are compared to determine the equivalent 1 Hz instrument height noise of the Poseidon-2 and TOPEX altimeters. This geophysical evaluation shows that the Ku-band 1-Hz range noise for both instruments is better than 1.7 cm at 2 m significant wave heights (H_1/3), exceeding error budget requirements for both missions. Furthermore, we show that the quality of these instruments allows optimal filtering of the 1-Hz along-track sea surface height data for sea surface slopes that can be used to calculate cross track geostrophic velocity anomalies at the baroclinic Rossby radius of deformation to better than 5 cm/sec precision along 87.5% of the satellite ground track between 2 and 60 degrees absolute latitude over the deep abyssal ocean (depths greater than 1000 m). This level of precision will facilitate scientific studies of surface geostrophic velocity variability using data from the Jason-1 and T/P Tandem Mission.     

Accuracy Assessment of Jason-1 and TOPEX/POSEIDON Along-Track Sea Surface Slope

Sea surface slope computed from along-track Jason-1 and TOPEX/POSEIDON (T/P) altimeter data at ocean mesoscale wavelengths are compared to determine the equivalent 1 Hz instrument height noise of the Poseidon-2 and TOPEX altimeters. This geophysical evaluation shows that the Ku-band 1-Hz range noise for both instruments is better than 1.7 cm at 2 m significant wave heights (H_1/3), exceeding error budget requirements for both missions. Furthermore, we show that the quality of these instruments allows optimal filtering of the 1-Hz along-track sea surface height data for sea surface slopes that can be used to calculate cross track geostrophic velocity anomalies at the baroclinic Rossby radius of deformation to better than 5 cm/sec precision along 87.5% of the satellite ground track between 2 and 60 degrees absolute latitude over the deep abyssal ocean (depths greater than 1000 m). This level of precision will facilitate scientific studies of surface geostrophic velocity variability using data from the Jason-1 and T/P Tandem Mission.     

海浪搅拌混合对北太平洋海表面温度模拟的影响

利用NCEP再分析风场驱动WAVEWATCH III海浪模式对北太平洋海域的海浪过程进行模拟,利用浮标观测资料对模拟出的海浪要素有效波高进行验证,发现他们之间具有很好的一致性。基于模式输出的有效波高等波浪要素,利用特征波参数化理论,在海洋环流模式中引入海浪搅拌混合作用,分析其对北太平洋海表面温度模拟的影响,初步数值模拟结果表明,sbPOM模式在考虑海浪搅拌混合作用以后,模拟精度进一步提升,这对提供一个准确的大气模式下边界条件具有重要作用。     

Applications of satellite altimetry to oceanography and geophysics

Satellite-borne altimeters have had a profound impact on geodesy, geophysics, and physical oceanography. To first order approximation, profiles of sea surface height are equivalent to the geoid and are highly correlated with seafloor topography for wavelengths less than 1000 km. Using all available Geos-3 and Seasat altimeter data, mean sea surfaces and geoid gradient maps have been computed for the Bering Sea and the South Pacific. When enhanced using hill-shading techniques, these images reveal in graphic detail the surface expression of seamounts, ridges, trenches, and fracture zones. Such maps are invaluable in oceanic regions where bathymetric data are sparse. Superimposed on the static geoid topography is dynamic topography due to ocean circulation. Temporal variability of dynamic height due to oceanic eddies can be determined from time series of repeated altimeter profiles. Maps of sea height variability and eddy kinetic energy derived from Geos-3 and Seasat altimetry in some cases represent improvements over those derived from standard oceanographic observations. Measurement of absolute dynamic height imposes stringent requirements on geoid and orbit accuracies, although existing models and data have been used to derive surprisingly realistic global circulation solutions. Further improvement will only be made when advances are made in geoid modeling and precision orbit determination. In contrast, it appears that use of altimeter data to correct satellite orbits will enable observation of basin-scale sea level variations of the type associated with climatic phenomena.     

Improving the Coastal Mean Dynamic Topography by Geodetic Combination of Tide Gauge and Satellite Altimetry

Abstract

The ocean mean dynamic topography (MDT) is the surface representation of the ocean circulation. The MDT may be determined by the ocean approach, which involves temporal averaging of numerical ocean circulation model information, or by the geodetic approach, wherein the MDT is derived using the ellipsoidal height of the mean sea surface (MSS), or mean sea level (MSL) minus the geoid as the geoid. The ellipsoidal height of the MSS might be estimated either by satellite or coastal tide gauges by connecting the tide gauge datum to the Earth-centred reference frame. In this article we present a novel approach to improve the coastal MDT, where the solution is based on both satellite altimetry and tide gauge data using new set of 302 tide gauges with ellipsoidal heights through the SONEL network. The approach was evaluated for the Northeast Atlantic coast where a dense network of GNSS-surveyed tide gauges is available. The typical misfit between tide gauge and satellite or oceanographic MDT was found to be around 9?cm. This misfit was found to be mainly due to small scale geoid errors. Similarly, we found, that a single tide gauge places only weak constraints on the coastal dynamic topography.