Three dimensional diagnostic study of the circulation in the South China Sea during winter 1998

On the basis of hydrographic data obtained from 28 November to 27 December, 1998, the three-dimensional structure of circulation in the South China Sea (SCS) is computed using a three-dimensional diagnostic model. The combination of sea surface height anomaly from altimeter data and numerical results provides a consistent circulation pattern for the SCS, and main circulation features can be summarized as follows: in the northern SCS there are a cold and cyclonic circulation C1 with two cores C1-1 and C1-2 northwest of Luzon and an anticyclonic eddy (W1) near Dongsha Islands. In the central SCS there is a stronger cyclonic circulation C2 with two cores C2-1 and C2-2 east of Vietnam and a weaker anticyclonic eddy W2 northwest of Palawan Island. A stronger coastal southward jet presents west of the eddy C2 and turns to the southeast in the region southwest of eddy C2-2, and it then turns to flow eastward in the region south of eddy C2-2. In the southern SCS there are a weak cyclonic eddy C3 northwest of Borneo and an anti-cyclonic circulation W3 in the subsurface layer. The net westward volume transport through section CD at 119.125°E from 18.975° to 21.725°N is about 10.3 × 10⁶ m³s⁻¹ in the layer above 400 m level. The most important dynamic mechanism generating the circulation in the SCS is a joint effect of the baroclinicity and relief (JEBAR), and the second dynamical mechanism is an interaction between the wind stress and relief (IBWSR). The strong upwelling occurs off northwest Luzon.     

Upper Ocean Four-Dimensional Variational Data Assimilation in the Arabian Sea and Bay of Bengal

A regional ocean circulation model with four-dimensional variational data assimilation scheme is configured to study the ocean state of the Indian Ocean region (65°E–95°E; 5°N–20°N) covering the Arabian Sea (AS) and Bay of Bengal (BoB). The state estimation setup uses 10 km horizontal resolution and 5 m vertical resolution in the upper ocean. The in-situ temperature and salinity, satellite-derived observations of sea surface height, and blended (in-situ and satellite-derived) observations of sea surface temperature alongwith their associated uncertainties are used for data assimilation with the regionally configured ocean model. The ocean state estimation is carried out for 61 days (1 June to 31 July 2013). The assimilated fields are closer to observations compared to other global state estimates. The mixed layer depth (MLD) of the region shows deepening during the period of assimilation with AS showing higher MLD compared to the BoB. An empirical forecast equation is derived for the prediction of MLD using the air–sea forcing variables as predictors. The surface and sub-surface (50 m) heat and salt budget tendencies of the region are also investigated. It is found that at the sub-surface, only the advection and diffusion temperature and salt tendencies are important.     

Sequential forecasting of the surface and subsurface conditions in the Japan Sea

This study estimates a realistic change of the Japan Sea by assimilating satellite measurements into an eddy-resolving circulation model. Suboptimal but feasible assimilation schemes of approximate filtering and nudging play essential roles in the system. The sequential update of error covariance significantly outperforms the asymptotic covariance in the sequential assimilation due to the irregular sampling patterns from multiple altimeter satellites. The best estimates show an average rms difference of only 1.2°C from the radiometer data, and also explain about half of the sea level variance measured by the altimeter observation. The subsurface conditions associated with the mesoscale variabilities are also improved, especially in the Tsushima Warm Current region. It is demonstrated that the forecast limit strongly depends on variable, depth, and location.     

基于卫星高度计和浮标漂流轨迹的海洋涡旋特征信息对比分析

本文基于Chelton提供的涡旋数据集和浮标漂流轨迹提取的涡旋结果,对1993—2015年的全球涡旋进行特征信息对比分析。结果表明,在全球范围内高度计涡旋数据集中的欧拉涡旋和浮标漂流轨迹提取的拉格朗日涡旋的配对成功率在空间分布上并不均衡,在中纬度(20°—60°S,20°—60°N)配对成功率最高可达25%,而在20°S—20°N区域内配对成功率不到10%。由于低纬度地转效应并不显著,卫星高度计无法有效观测到涡旋,但通过浮标漂流轨迹识别出的拉格朗日涡旋却大量存在,这说明在低纬度区域内,采用漂流浮标手段对涡旋进行观测,能够有效地弥补卫星高度计识别涡旋的区域限制。进一步分析表明,总体而言,提取的欧拉涡旋半径要大于拉格朗日涡旋闭合回路半径。两种识别方法获得的涡旋(闭合回路)在20°—50°S, 20°—50°N的副热带和中纬度海区半径大致相当; 20°S—20°N度以内(特别是近赤道区域)、高纬度区域以及西边界流区域,欧拉涡旋半径是同期拉格朗日涡旋闭合回路半径的3倍或更多。此外,对配对涡旋的Rossby数分析结果显示,拉格朗日涡旋较小的闭合回路对应较大的平均相对涡度,这表明浮标在被中尺度涡俘获后,更容易在相对涡度较大的地方(如中尺度涡中心、中尺度涡边缘等)形成闭合回路。     

A 4D-variational approach applied to an eddy-permitting North Atlantic configuration: Synthetic and real data assimilation of altimeter observations

The increasing number of oceanic observations calls for the use of synthetic methods to provide consistent analyses of the oceanic variability that will support a better understanding of the underlying mechanisms. In this study, a 1/3° eddy-permitting model of the North Atlantic (from 20°S to 70°N) is combined with a 4D-variational method to estimate the oceanic state from altimeter observations. This resolution allows a better extraction of the physical content of altimeter data since the model spatial scales are more consistent with the data than coarser assimilation exercises because of a lower error in model representativity. Several strategies for the assimilation window are tested through twin experiments carried out under the following conditions: different window lengths and either a quasi-static (also known as progressive) variational assimilation with progressive extension of the window, or a simpler direct method without prior assimilation. From our set of experiments, the most efficient strategy is the use of both a simple direct assimilation method and a 90-day window. The assimilation of synthetic altimeter data constrains the model-temperature, -salinity and -velocity fields mainly over the first 1300 m where the error is the largest. Improvements occur not only in quiescent regions, but also in more energetic meso-scale regimes. Despite the existence of model- and surface forcing-errors as well as large errors in the first guess, the assimilation of real altimeter data proves to be consistent with our twin experiments. Indeed, the analyses show a better detachment of the Gulf Stream, weaker regional biases and more accurate positions for meso-scale structures. Independent hydrographic data (Argo floats and CTD cruises) and transports estimates along the OVIDE 2002 cruise show an improvement of the analysed oceanic state with respect to the assimilation-free case though water mass properties are still incorrectly represented. After assimilation, the North Atlantic heat transport in the model is in good agreement with independent estimates based on hydrographic data.     

Stochastic study of the temperature response of the upper ocean to uncertainties in the atmospheric forcing in an Atlantic OGCM

The impact of errors in atmospheric forcing on the behaviour of ocean models is a fundamental issue for ocean modellers and data assimilation and one that has yet to be fully addressed. In this study, we use a stochastic modelling approach with 50 7-months (September–March) primitive equation eddy permitting (1/4°) integrations. We investigate the response of the oceanic circulation to atmospheric uncertainties, focusing principally on their impact on the upper oceanic temperature field. The ensemble is generated by perturbing the wind, atmospheric temperature and incoming solar radiation of the ERA40 reanalysis. Each perturbation consists of a random combination of the 20 dominant EOFs of the difference between the ERA40 and NCEP/CORE reanalysis datasets. The ensemble standard deviation of various interfacial and oceanic quantities is then examined in the upper 200 m of three distinct regions of the North Atlantic: in the Gulf Stream, in the Northern Tropical band and in the North East Atlantic. These show that even a very small perturbation of the atmospheric variables can lead to significant changes in the ocean properties and that regions of oceanic mesoscale activity are the most sensitive. The ocean response is driven by vertical diffusivity and eddy activity. The role of subsurface currents is also crucial in carrying the eddy signal away from the regions of mesoscale activity. Finally, the decorrelation time scale of the mesoscale activity is critical in determining the amplitude of the oceanic response.     

Mean sea level determination from satellite altimetry

The primary experiment on the Geodynamics Experimental Ocean Satellite‐3 (GEOS‐3) is the radar altimeter. This experiment's major objective is to demonstrate the utility of measuring the geometry of the ocean surface, i.e., the geoid. Results obtained from this experiment so far indicate that the planned objectives of measuring the topography of the ocean surface with an absolute accuracy of ±5 m can be met and perhaps exceeded. The GEOS‐3 satellite altimeter measurements have an instrument precision in the range of ±25 cm to ±50 cm when the altimeter is operating in the “short pulse”; mode. After one year's operations of the altimeter, data from over 5,000 altimeter passes have been collected. With the mathematical models developed and the altimeter data presently available, mapping of local areas of ocean topography has been realized to the planned accuracy levels and better. This paper presents the basic data processing methods employed and some interesting results achieved with the early data. Plots of mean sea surface heights as inferred by the altimeter measurements are compared with a detailed 1^o × 1° gravimetric geoid.     

Eddy Detection Over Southern Indian Ocean Using TOPEX/POSEIDON Altimeter Data

Mesoscale eddies constitute the most energetic component of the variability of ocean currents. An attempt has been made for the detection of oceanic mesoscale eddy signatures over the Southern Indian Oceanic (SIO) regions using the dynamic topography derived from TOPEX/POSEIDON (T/P) altimeter data, by the signal processing technique, called matched filtering. After applying all the ocean and atmospheric corrections, data of a complete cycle of T/P over SIO has been used for detection of eddy signatures. The geoid undulations are removed from the data of corrected sea surface height from T/P and the resulting dynamic topographic data are passed through a matched filter designed to detect a generic eddy signature of Gaussian signal embedded in noise. The filter is optimized to detect eddies with amplitude 20 to 30 cm and diameters roughly 100?250 km. Out of all the analyzed data of T/P orbits over SIO a few examples are presented for brevity. Qualitative verification of eddies is done with some independent T/P sea level anomaly data over the region. The analysis shows that the matched filtering technique is most suitable for monitoring eddy signatures along the subsatellite track instantly over the remote and most hostile regions of the southern global oceans.     

Comparison of EnOI and EnKF regional ocean reanalysis systems

This study compares two regional eddy resolving ocean reanalysis systems, based on the ensemble Kalman filter (EnKF) and ensemble optimal interpolation (EnOI), focusing on data assimilation aspects. Both systems are configured for the Tasman Sea using the same ocean model with 0.1° resolution and commonly available observations of satellite altimetry, sea surface temperature and subsurface temperature and salinity. The primary goals are to quantify the difference in performance of the EnKF and EnOI and investigate how important this difference might be from an oceanographic perspective. We find that both systems generally constrain mesoscale circulation in the region, with some exceptions for the East Australian Current separation region, the most energetic and chaotic part of the domain. Overall, the EnKF is found to consistently outperform the EnOI, producing on average 9–21% smaller innovations. The EnKF also has better forecast skill relative to the persisted analysis than the EnOI. For SST the EnKF forecast outperforms persisted analysis by about 17%, which indicates that the surface circulation is mainly constrained. The EnKF and EnOI are shown to produce qualitatively different increments of unobserved or sparsely observed variables; however, we find only moderate improvements of the EnKF over EnOI in subsurface temperature fields when compared against withheld XBT observations. We attribute this lack of a major improvement in subsurface reconstruction to the inability of the EnKF to linearly constrain the system due to initialisation shock, model error caused by open boundaries, and possibly insufficient observations.     

卫星高度计资料在三维海温和盐度数值预报中的应用

随着卫星遥感观测技术的发展,越来越多的卫星观测资料被应用于数值模式的同化研究中.基于国家海洋环境预报中心西北太平洋三维湿盐流预报系统,利用法国CLS中心的沿轨高度计资料的海表面高度异常的融合数据,结合基于三维变分的OVALS(ocean variational analysis system)同化系统,在垂向将海面高度...     

北太平洋涡旋对基于细尺度参数化的海洋内部混合的影响

基于Vector Geometry方法对2016—2018年的高度计资料进行涡旋识别,并使用细尺度参数化方法和Argo数据计算了涡旋附近的海洋内部扩散率,分析了北太平洋的涡旋对海洋内部混合的影响。结果显示,研究区域在涡旋影响下的平均扩散率比无涡旋影响下的值大6%,并且气旋涡增强了600—1200m深度的混合,对600—900m深度的混合影响最大,可达18%;反气旋涡明显增强了300—900m深度的混合,但对900—1200m深度的混合没有明显影响。随着与涡旋中心距离的增大,涡旋外围混合扩散率缓慢减小,涡旋内部混合扩散率变化不明显,此结果与2014年3—10月在24°—36°N、132°—152°E区域的一个个例分析结果一致。此外,随着涡旋强度的增大,海洋内部混合明显增强。统计结果表明,在研究区域, 90%的扩散率值在10~(-5.5)—10~(-4)m~2/s范围内。     

Diagnostic calculation of the circulation in the South China Sea during summer 1998

On the basis of hydrographic data obtained from 12 June to 6 July, 1998, the three-dimensional structure of circulation in the South China Sea (SCS) is computed using a three-dimensional diagnostic model. The combination of sea surface height anomaly from altimeter data and numerical results provides a consistent circulation pattern for the SCS, and the main circulation features can be summarized as follows: In the northern SCS there are a cyclonic eddy C1 near Dongsha Islands and an anti-cyclonic eddy W1 west of Luzon Island. In the central SCS a strong anti-cyclonic eddy W3 and a cyclonic eddy C3 compose a quasi-dipole southeast of Vietnam. A coastal northward jet is present at the western boundary near the Vietnam coast above 300 m level. This northward coastal jet flows northward and turns eastward at about 14°N, and then flows southeastward into the area between eddies W3 and C3. In the southern SCS the current is weaker. The most important dynamic mechanism underlying the circulation in the SCS is the joint effect of the baroclinicity and relief (JEBAR), and the second dynamical mechanism is the interaction between the wind stress and relief (IBWSR). Comparison of the characters of circulation in the SCS during summer 2000 with that during summer 1998 reveals no obvious variability of the main characteristics.     

Observed Variability of Thermohaline Fields,Currents and Eddies in the Western Bay of Bengal During BOBMEX-99

Field measurements during the Bay of Bengal Monsoon Experiment (BOBMEX-99), data from a deep sea moored buoy, and satellite altimeter were used to describe variability in the hydrographic and meso-scale features in the Bay of Bengal (BoB) during the summer monsoon of 1999. The thermohaline fields showed two regions of upsloping of isopleths centered at 82°E and 84.75°E, ～110 km and 450 km away from the coast, respectively, followed by downsloping. The upsloping/downsloping of isopleths and the alternating currents was part of cyclonic and anti-cyclonic circulation patterns in the western BoB. In this region, both wind and current were important in the dynamics of coastal upwelling. The observations showed a relationship between the propagating waves and eddy on variability of thermohaline fields. On an annual cycle, four Kelvin waves were observed in the BoB, but only the downwelling Kelvin wave formed during October entered the Arabian Sea. During the monsoon season, four eddies were formed in the western BoB, of which the anticyclonic eddy centered at 15°N, 84°E and the cyclonic eddy centered at 17.5°N, 84.5°E were prominent. The baroclinic instability caused by the opposing currents along the east coast and the wind stress curl favored the formation of eddies. Okhubo-Weiss and Isern-Fontanet parameter confirmed the presence of eddies in the BoB.     

Evolution and propagation of amesoscale eddy in the northern South China Sea during winter

In situ observations,satellite data,and the output from an eddy-resolving ocean circulation model were used to study the generations and propagations of an anticyclonic eddy in the northern South China Sea (NSCS) during the winter of 2009-2010.In the NSCS,the anticyclonic eddy firstly appeared to the west of the Luzon Strait,migrated generally along the continental slope and dissipated around the Xisha Archipelago. The evolution of the warm eddy contains three phases:development,maturation,and decay.The eddy mainly stayed near 119.7°E in December and then gradually moved to 118.7°E until January 15,when its intensity,as indicated by the thermocline temperature and salinity anomalies,increased significantly,reflecting the growth of the eddy.The eddy reached its peak on January 15 and persisted until February 23. During this period,the eddy propagated westward to 116.4°E.After,the warm eddy weakened significantly and dissipated finally near the Xisha Archipelago.     

气候模式中海洋数据同化对热带降水偏差的影响

本文采用海洋卫星观测海表温度(SST)和海面高度异常(SLA)数据,对国家海洋局第一海洋研究所地球系统模式FIO-ESM(First Institute of Oceanography Earth System Model version 1.0)中海洋模式分量进行了集合调整卡尔曼滤波(EAKF)同化,对比分析了大气环流、湿度和云量对海洋数据同化的响应,探讨了海洋同化对热带降水模拟偏差的影响。结果表明:海洋数据同化能有效改善海表温度和上层海洋热含量的模拟,30°S~30°N纬度带内年平均SST的绝均差降低60%。同化后大气模式模拟的赤道两侧信风得到明显改善,上升气流在赤道以北热带地区增强而在赤道以南热带地区减弱,热带降水模拟的动力结构更为合理,水汽和云量分布也更切合实际。热带年平均降水的空间分布和强度在同化后均得到改善,赤道以南的纬向年平均降水峰值显著降低,降水偏差明显减小,同化后30°S~30°N纬度带内年平均降水绝均差降低35%。     

Impact of Ocean Mean Dynamic Topography on Satellite Data Assimilation

The response of an eddy-permitting ocean model to changes imposed by the use of different mean dynamic topographies (MDT) is analyzed in a multivariate assimilation context, allowing the evaluation of this impact, not only on the surface circulation, but also on the interior ocean representation. The assimilation scheme is a reduced-order sequential Kalman filter (SEEK). In a first set of experiments, high resolution sea surface temperature, along-track sea surface height and sea surface salinity from climatology are assimilated into a 1/3° resolution North and Tropical Atlantic version of the HYCOM model. In a second experiment, in situ profile data are assimilated in addition to the surface measurements.

The first set of experiments illustrates that important differences in the representation of the horizontal model circulation pattern are related to differences in the MDT used. The objective of assimilation is to improve the representation of the 3D ocean state. However, the imperfect representation of the mean dynamic topography appears to be an important limiting factor with regard to the degree of realism obtained in the simulated flow.

Vertical temperature and salinity profiles are key observations to drive a general circulation ocean model toward a more realistic state. The second set of experiments shows that assimilating them in addition to sea surface measurements is a far from trivial exercise. A specific difficulty is due to inconsistencies between the dynamic topography diagnosed from in situ observations and that diagnosed from sea surface height. These two fields obtained from different data sources do not contain exactly the same information. In order to overcome this difficulty, a strategy is proposed and validated.     

Observation and parameterization of small icebergs: Drifting breakwaters in the southern ocean

The variability of small-size iceberg distributions is revealed from a novel analysis of satellite altimeter data. A strong annual cycle is modulated by pulse-like events confined to single ocean basins, with dense iceberg populations in the South Atlantic in 2004–2005, and in the South Pacific in 2008. Anomalies in sea surface temperatures of the order of 1 °C may be related to the iceberg distribution. Icebergs also appear very strongly associated with anomalies in the heights of ocean waves. A preliminary parameterization of wave blocking by icebergs significantly reduces wave model errors in the region south of 45° South, and has a perceptible influence on all the west coasts of the Southern hemisphere.     

Numerical Study of the Upper-Layer Circulation in the South China Sea

Upper-layer circulation in the South China Sea has been investigated using a three-dimensional primitive equation eddy-resolving model. The model domain covers the region from 99° to 122°E and from 3° to 23°N. The model is forced by the monthly averaged European Centre for Medium-Range Weather Forecasts (ECMWF) model winds and the climatological monthly sea surface temperature data from National Oceanographic Data Center (NODC). Inflow and outflow through the Taiwan Strait and the Sunda shelf are prescribed monthly from the Wyrtki estimates. Inflow of the Kuroshio branch current in the Luzon Strait is assumed to have a constant volume transport of 12 Sv (1 Sv = 10⁶ m³/s), and the outflow from the open boundary to the east of Taiwan is adjusted to ensure the net volume transport through all open boundaries is zero at any instant. The model reveals that a cyclonic circulation exists all year round in the northern South China Sea. During the winter time this cyclonic eddy is located off the northwest of Luzon, coinciding with the region of positive wind stress curl in this season. This cyclonic eddy moves northward in spring due to the weakening of the northeast winds. The cyclonic circulation becomes weak and stays in the continental slope region in the northern South China Sea in the summer period. The southwest wind can raise the water level along the west coast of Luzon, but there is no anticyclonic circulation in the northern South China Sea. After the onset of the northeast monsoon winds in fall, the cyclonic eddy moves back to the region off the west coast of Luzon. In the southern South China Sea and off the Vietnam coast, the model predicts a similar flow structure as in the previous related studies. This revised version was published online in July 2006 with corrections to the Cover Date.