Application of altimetry data assimilation on mesoscale eddies simulation

Mesoscale eddy plays an important role in the ocean circulation. In order to improve the simulation accuracy of the mesoscale eddies, a three-dimensional variation (3DVAR) data assimilation system called Ocean Variational Analysis System (OVALS) is coupled with a POM model to simulate the mesoscale eddies in the Northwest Pacific Ocean. In this system, the sea surface height anomaly (SSHA) data by satellite altimeters are assimilated and translated into pseudo temperature and salinity (T-S) profile data. Then, these profile data are taken as observation data to be assimilated again and produce the three-dimensional analysis T-S field. According to the characteristics of mesoscale eddy, the most appropriate assimilation parameters are set up and testified in this system. A ten years mesoscale eddies simulation and comparison experiment is made, which includes two schemes: assimilation and non-assimilation. The results of comparison between two schemes and the observation show that the simulation accuracy of the assimilation scheme is much better than that of non-assimilation, which verified that the altimetry data assimilation method can improve the simulation accuracy of the mesoscale dramatically and indicates that it is possible to use this system on the forecast of mesoscale eddies in the future.

Tectonics of the southern tip of the Parece Vela Basin, Philippine Sea Plate

We report new geophysical and petrological data collected at the southern tip of the Parece Vela Basin in the Philippine Sea. The Parece Vela Basin, which was formed as a backarc basin behind proto Mariana arc-trench system from late Oligocene to middle Miocene, provides us a good opportunity to study the nature of successive backarc basin formations in the Philippine Sea and the relationship between are and backarc magmatisms. Regional bathymetric map derived from satellite altimetry shows that the southern tip of the basin, now located just west of the Yap arc-trench system, has unique morphological and tectonic features which include: 1) the absence of spreading center or its trace, 2) shallow average depth, and 3) enigmatic curved structures. Our newly collected high-resolution bathymetric data reveal that the spreading fabric similar to the central Parece Vela Basin exists to the north of 9 degrees 20'N. Thus it appears that the present-day Yap arc and backarc region represent the western half of the seafloor that was produced by the early E-W and the following NE-SW spreading in the northern and central Parece Vela Basin, and that the eastern counterpart now lies west of the West Mariana Ridge. Unlike the northern Parece Vela Basin, there appears to be no evidence for a systematic propagation of spreading center in the southern part. Instead two rift segments, one which extends from the central Parece Vela Basin and the other which lies within the western remnant arc (Kyushu-Palau Ridge), overlap at the southern tip of the basin, producing a complex seafloor that includes curvilinear deeps and deformed topographic highs. (c) 2007 Elsevier B.V. All rights reserved.

Tectonic and geodynamic pattern of marginal seas on the West Pacific inferred from multi-satellite altewetry geoid

With the rapid development of satellite observations, we can use the altimetry geoid to study submarine tectonics and geodynamics. On the basis of the 4' x 4' geoid undulation calculated from altimeter data of Geosat, ERS-1/2 and Topex/Poseidon on the West Pacific, located between 0degreesN similar to 45degreesN, 100degreesC similar to 150degreesE, Bouguer, Glenni and isostatic geoid undulation are obtained from correction of gravitational potential of the global topography and isostacy. Moho discontinuity depth is inversed by the Glenni geoid undulation, and the stress field from small-scale mantle convection is reasonably calculated from the isostatic geoid undulation. The results show that within the Philippine Sea and the South China Sea, short-wavelength lineations of the geoid undulation are parallel or cross to magnetic lineations and rifting ridges. The Moho depth of marginal sea basins becomes shallow southward, and its values are similar to that of the Philippine Sea. These facts show that strength of tectonic activities are almost the same on the both sides of the Ryukyu-Taiwan-Philippine are. Various kinds of tectonic features with different driving mechanisms of small-middle and large-scale of mantle convection, however, display a special pattern of tectonics and geodynamics of the continental marginal seas distinguished from oceans and continents.

The variability of eddy kinetic energy in the South China Sea deduced from satellite altimeter data

We used fifteen years (1993-2007) of altimetric data, combined from different missions (ERS-1/2, TOPEX/Poseidon, Jason-1, and Envisat), to analyze the variability of the eddy kinetic energy (EKE) in the South China Sea (SCS). We found that the EKE ranged from 64 cm(2)/s(2) to 1 390 cm(2)/s(2) with a mean value of 314 cm(2)/s(2). The highest EKE center was observed to the east of Vietnam (with a mean value of 509 cm(2)/s(2)) and the second highest EKE region was located to the southwest of Taiwan Island (with a mean value of 319 cm(2)/s(2)). We also found that the EKE structure is the consequence of the superposition of different variability components. First, interannual variability is important in the SCS. Spectral analysis of the EKE interannual signal (IA-EKE) shows that the main periodicities of the IA-EKE to the east of Vietnam, to the southwest of Taiwan Island, and in the SCS are 3.75, 1.87, and 3.75 years, respectively. It is to the south of Taiwan Island that the IA-EKE signal has the most obvious impact on EKE variability. In addition, the IA-EKE exhibit different trends in different regions. An obvious positive trend is observed along the east coast of Vietnam, while a negative trend is found to the southwest of Taiwan Island and in the east basin of Vietnam. Correlation analysis shows that the IA-EKE has an obvious negative correlation with the SSTA in Nio3 (5A degrees S-5A degrees N, 90A degrees W-150A degrees W). El Nio-Southern Oscillation (ENSO) affects the IA-EKE variability in the SCS through an atmospheric bridge-wind stress curl over the SCS. Second, the seasonal cycle is the most obvious timescale affecting EKE variability. The locations of the most remarkable EKE seasonal variabilities in the SCS are to the east of Vietnam, to the southwest of Taiwan, and to the west of Philippines. To the east of Vietnam, the seasonal cycle is the dominant mechanism controlling EKE variability, which is attributed primarily to the annual cycle there of wind stress curl. In this area, the maximum EKE is observed in autumn. To the southwest of Taiwan Island, the EKE is enlarged by the stronger SCS circulation, which is caused by the intrusion branch from the Kuroshio in winter. Finally, intra-annual and mesoscale variability, although less important than the former, cannot be neglected. The most obvious intra-annual and mesoscale variability, which may be the result of baroclinic instability of the background flow, are observed to the southwest of Taiwan Island. Sporadic events can have an important effect on EKE variability.

Bifurcation of Pacific North Equatorial Current at the surface

The grid altimetry data between 1993 and 2006 near the Philippines were analyzed by the method of Empirical Orthogonal Function (EOF) to study the variation of bifurcation of the North Equatorial Current at the surface of the Pacific. The relatively short-term signals with periods of about 6 months, 4 months, 3 months and 2 months are found besides seasonal and interannual variations mentioned in previous studies. Local wind stress curl plays an important role in controlling variation of bifurcation latitude except in the interannual timescale. The bifurcation latitude is about 13.3A degrees N in annual mean state and it lies at the northernmost position (14.0A degrees N) in January, at the southernmost position (12.5A degrees N) in July. The amplitude of variation of bifurcation latitude in a year is 1.5A degrees, which can mainly be explained as the contributions of the signals with periods of about 1 year (1.2A degrees) and 0.5 year (0.3A degrees).

Eddy generation and evolution in the North Pacific Subtropical Countercurrent (NPSC) zone

The seasonal generation and evolution of eddies in the region of the North Pacific Subtropical Countercurrent remain poorly understood due to the scarcity of available data. We used TOPEX/POSEIDON altimetry data from 1992 to 2007 to study the eddy field in this zone. We found that velocity shear between this region and the neighboring North Equatorial Current contributes greatly to the eddy generation. Furthermore, the eddy kinetic energy level (EKE) shows an annual cycle, maximum in April/May and minimum in December/January. Analyses of the temporal and spatial distributions of the eddy field revealed clearly that the velocity shear closely related to baroclinic instability processes. The eddy field seems to be more zonal than meridional, and the energy containing length scale shows a surprising lag of 2-3 months in comparison with the 1-D and 2-D EKE level. A similar phenomenon is observed in individual eddies in this zone. The results show that in this eddy field band, the velocity shear may drive the EKE level change so that the eddy field takes another 2-3 months to grow and interact to reach a relatively stable state. This explains the seasonal evolution of identifiable eddies.

卫星高度计资料在海洋潮汐研究中的应用

本文依据收集到的392个地面验潮站8个主要分潮（M2、S2、K1、O1、N2、K2、P1及Q1）的调和常数，对现有7个全球大洋潮汐模式的准确度进行了检验，结果显示各模式在深海区域均达到了比较高的准确度，相互之间差别也不大。经验模式GOT00和CSR4.0、同化模式NAO99、反演同化模式TPXO7.0、数值同化模式FES2002和FES2004的M2分潮均方根偏差在3 cm左右，其它分潮（S2、K1、O1、N2、K2、P1及Q1）大约在1～2 cm。本文还依据中国近海18个岛屿的调和常数对其中的5个大洋潮汐模式的准确度进行了检验，结果表明，M2分潮均方根偏差在6～14 cm，明显高于大洋部分的偏差，其中日本国家天文台的潮汐模式NAO99在中国近海的结果相对较准确。 我们利用1992年8月至2008年8月的TOPEX/POSEIDON和JASON-1(T/P-J)卫星高度计资料，对沿卫星轨道的302816个站点进行了14个分潮的潮汐调和分析，得到了全球大洋潮汐的8个主要分潮以及2个气象分潮Sa、Ssa的经验同潮图。主要结果有：（1）各分潮在卫星上升轨道与下降轨道的交叉点(约7000个)相关性分析表明：M2分潮的振幅和迟角的相关系数很高（分别为0.9965和0.9961）；S2，K1，O1和Sa分潮也有较好的相关性（相关系数为0.94～0.99）；（2）该结果与392地面个验潮站吻合较好，其中M2分潮的振幅、迟角和向量的均方根偏差分别为：1.73 cm，2.340和2.93 cm；S2，K1和O1分潮的振幅、迟角和向量的均方根偏差为1 cm左右，5.250～7.270和1.5～2.1 cm，该精度与最近几年国际上的主要大洋潮汐模式的准确度相近；（3）首次通过卫星资料获得了Sa、Ssa分潮的同潮图。周期为1年的Sa分潮与大洋105个地面站相比，振幅、迟角和向量的均方根偏差分别为1.50 cm、18.360和2.16 cm。在此基础上，进一步分析了构成Sa、Ssa气象分潮的两个主要因素（海水密度以及海面气压）在全球的分布。 在T/P-J等卫星资料无法覆盖到南大洋和北冰洋，本文利用Princeton Ocean Model（POM）进行了数值模拟，模拟结果与162个地面实测站（其中南大洋30个，北冰洋132个）的观测比较一致。基于卫星资料分析的结果和数值模拟结果合并得到了全球大洋的8个主要分潮同潮图。在此基础上通过全球潮汐能量耗散的计算得到潮能通量的分布，并得到全球M2、S2、K1和O1分潮的潮汐能量耗散率为2.431TW、0.401TW、0.336TW和0.176TW。 本文还利用卫星资料对南海潮汐进行了研究，在中国南海，获得了主要的半日潮、全日潮、四分日分潮和长周期分潮（M2，S2，N2，K2，K1，O1，P1，Q1，M4, MS4，Sa, Ssa）的经验同潮图。与南海沿岸94个地面验潮站的数据符合得比较好，M2，S2，K1及O1等4个主要分潮的平均振幅差为2～4 cm，均方根偏差分别是9～11 cm.其它4个主要分潮N2，K2，P1，Q1的平均振幅差为1～2 cm，均方根偏差为2～4 cm。此外，本文还利用卫星高度计资料潮汐分析结果沿卫星轨道进行高通滤波，分离得出中国近海的M2，S2，K1及O1分潮的内潮信息。