Using an eddy-tracking algorithm to understand the impact of assimilating altimetry data on the eddy characteristics of the Agulhas system

A complex and highly dynamical ocean region, the Agulhas Current System plays an important role in the transfer of energy, nutrients and organic material from the Indian to the Atlantic Ocean. Its dynamics are not only important locally, but affect the global ocean-atmosphere system. In working towards improved ocean reanalysis and forecasting capabilities, it is important that numerical models simulate mesoscale variability accurately—especially given the scarcity of coherent observational platforms in the region. Data assimilation makes use of scarce observations, a dynamical model and their respective error statistics to estimate a new, improved model state that minimises the distance to the observations whilst preserving dynamical consistency. Qualitatively, it is unclear whether this minimisation directly translates to an improved representation of mesoscale dynamics. In this study, the impact of assimilating along-track sea-level anomaly (SLA) data into a regional Hybrid Coordinate Ocean Model (HYCOM) is investigated with regard to the simulation of mesoscale eddy characteristics. We use an eddy-tracking algorithm and compare the derived eddy characteristics of an assimilated (ASSIM) and an unassimilated (FREE) simulation experiment in HYCOM with gridded satellite altimetry-derived SLA data. Using an eddy tracking algorithm, we are able to quantitatively evaluate whether assimilation updates the model state estimate such that simulated mesoscale eddy characteristics are improved. Additionally, the analysis revealed limitations in the dynamical model and the data assimilation scheme, as well as artefacts introduced from the eddy tracking scheme. With some exceptions, ASSIM yields improvements over FREE in eddy density distribution and dynamics. Notably, it was found that FREE significantly underestimates the number of eddies south of Madagascar compared to gridded altimetry, with only slight improvements introduced through assimilation, highlighting the models’ limitation in sustaining mesoscale activity in this region. Interestingly, it was found that the threshold for the maximum eddy propagation velocity in the eddy detection scheme is often exceeded when data assimilation relocates an eddy, causing the algorithm to interpret the discontinuity as eddy genesis, which directly influences the eddy count, lifetime and propagation velocity, and indirectly influences other metrics such as non-linearity. Finally, the analysis allowed us to separate eddy kinetic energy into contributions from detected mesoscale eddies and meandering currents, revealing that the assimilation of SLA has a greater impact on mesoscale eddies than on meandering currents.     

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.     

Three-dimensional structure of mesoscale eddies in the western tropical Pacific as revealed by a high-resolution ocean simulation

The three-dimensional structure of mesoscale eddies in the western tropical Pacific(6°S–20°N, 120°E–150°E)is investigated using a high-resolution ocean model simulation. Eddy detection and eddy tracking algorithms are applied to simulated horizontal velocity vectors, and the anticyclonic and cyclonic eddies identified are composited to obtain their three-dimensional structures. The mean lifetime of all long-lived eddies is about 52 days, and their mean diameter is 147 km. Two typical characteristics of mesoscale eddies are revealed and possible dynamic explanations are analyzed. One typical characteristic is that surface eddies are generally separated from subthermocline eddies along the bifurcation latitude(~13°N) of the North Equatorial Current in the western tropical Pacific, which may be associated with different eddy energy sources and vertical eddy energy fluxes in subtropical and tropical gyres. Surface eddies have maximum swirl velocities of 8–9 cm s~(-1) and can extend to about 1500 m depth. Subthermocline eddies occur below 200 m, with their cores at about 400–600 m depth, and their maximum swirl velocities can reach 10 cm s~(-1). The other typical characteristic is that the meridional velocity component of the eddy is much larger than the zonal component. This characteristic might be due to more zonal eddy pairs(two eddies at the same latitude),which is also supported by the zonal wavelength(about 200 km) in the high-frequency meridional velocity component of the horizontal velocity.     

A mesoscale eddy detection method of specific intensity and scale from SSH image in the South China Sea and the Northwest Pacific

Mesoscale eddies exist almost everywhere in the ocean and play important roles in the ocean circulation of the world. These eddies may cause sound spread singular regions and bring great influences to the upwater ship and underwater aircraft. Due to the lack of hydrographic survey datasets, study of mesoscale eddies has been greatly restricted. Fortunately, satellite altimeter provided an effective way to study mesoscale eddies. An automatic detection algorithm is introduced to detect mesoscale eddies of specific intensity and spatial/temporal scale based on satellite sea surface height (SSH) data and the algorithm is applied in a strong eddy activity region: the South China Sea and the Northwest Pacific. The algorithm includes four steps. The first step is preprocessing of the SSH image, which includes elimination of error SSH data and interpolation. The second step is to detect suspected mesoscale eddies from preprocessed SSH images by dynamic threshold adjustment and morphological method, and the suspected mesoscale eddy detection includes two procedures: suspected mesoscale eddy core region detection and suspected mesoscale eddy brim extraction. The third step is to pick out mesoscale eddies satisfied with specified criteria from suspected mesoscale eddies. The criteria include three items, that is, intensity criterion, spatial scale, criterion and temporal scale criterion. The last step is algorithm performance analysis and verification. The algorithm has the capability of adaptive parameter adjustment, and can extract mesoscale eddies of interested intensity and spatial/temporal scale. The paper can provide a basis for analyzing space-time characteristics of mesoscale eddy in the South China Sea and the Northwest Pacific.     

The mesoscale eddies and Kuroshio transport in the western North Pacific east of Taiwan from 8-year (2003–2010) model reanalysis

The relationship between the Kuroshio volume transport east of Taiwan (~24°N) and the impinging mesoscale eddies is investigated using 8-year reanalysis of a primitive equation ocean model that assimilates satellite altimetry and SST data. The mean and fluctuations of the model Kuroshio transport agree well with the available observations. Analysis of model dynamic heights and velocity fields reveals three dominant eddy modes. The first mode describes a large eddy of ~500 km in diameter, centered at ~22° N. The second mode describes a pair of the north–south counter-rotating eddies of?~?400 km in diameter each, centered at 23° and 20° N, respectively. The third mode describes a pair of the east–west counter-rotating eddies of?~?300 km in diameter each, centered at 21° N. The associated velocity fields indicate eddies extending to 600–700 m in depth with vertical shears concentrated in the upper 400 m. All three modes and the model Kuroshio transport have similar dominant timescales of 70–150 days and generally are coherent. The decreased Kuroshio volume transports typically are associated with the impinging cyclonic eddies and the increased transports with the anticyclonic eddies. Selected drifter trajectories are presented to illustrate the three eddy modes and their correspondence with the varying Kuroshio transports.     

A numerical study of the circulation and monthly-to-seasonal variability in the Caribbean Sea: the role of Caribbean eddies

This study examines the circulation and associated monthly-to-seasonal variability in the Caribbean Sea using a regional ocean circulation model. The model domain covers the region between 99.0 and 54.0°W and between 8.0 and 30.3°N, with a horizontal resolution of 1/6°. The ocean circulation model is driven by 6-hourly atmospheric reanalysis data from the National Center for Environmental Prediction and boundary forcing extracted from 5-day global ocean reanalysis data produced by Smith et al. (Mercator Newsletter 36:39–49, 2010), and integrated for 7 years. A comparison of model results with observations demonstrates that the regional ocean circulation model has skill in simulating circulation and associated variability in the study region. Analysis of the model results, as well as a companion model run that uses steady annual mean forcing, illustrates the role of Caribbean eddies for driving monthly-to-seasonal circulation variability in the model. It is found that vertically integrated transport between Nicaragua and Jamaica is influenced by the interaction between the density perturbations associated with Caribbean eddies and the Nicaraguan Ridge. The impact of Caribbean eddies squeezing through the Yucatan Channel is also discussed.     

On the variability of the flow along the Meso-American Barrier Reef system: a numerical model study of the influence of the Caribbean current and eddies

A high resolution (3–8 km grid), 3D numerical ocean model of the West Caribbean Sea (WCS) is used to investigate the variability and the forcing of flows near the Meso-American Barrier Reef System (MBRS) which runs along the coasts of Mexico, Belize, Guatemala and Honduras. Mesoscale variations in velocity and temperature along the reef were found in seasonal model simulations and in observations; these variations are associated with meandering of the Caribbean current (CC) and the propagation of Caribbean eddies. Diagnostic calculations and a simple assimilation technique are combined to infer the dynamically adjusted flow associated with particular eddies. The results demonstrate that when a cyclonic eddy (negative sea surface height anomaly (SSHA)) is found near the MBRS the CC shifts offshore, the cyclonic circulation in the Gulf of Honduras (GOH) intensifies, and a strong southward flow results along the reef. However, when an anticyclonic eddy (positive SSHA) is found near the reef, the CC moves onshore and the flow is predominantly westward across the reef. The model results help to explain how drifters are able to propagate in a direction opposite to the mean circulation when eddies cause a reversal of the coastal circulation. The effect of including the Meso-American Lagoon west of the Belize Reef in the model topography was also investigated, to show the importance of having accurate coastal topography in determining the variations of transports across the MBRS. The variations found in transports across the MBRS (on seasonal and mesoscale time scales) may have important consequences for biological activities along the reef such as spawning aggregations; better understanding the nature of these variations will help ongoing efforts in coral reef conservation and maintaining the health of the ecosystem in the region.     

Composite eddy structures on both sides of the Luzon Strait and influence factors

Combining Argo observations with satellite remote sensing data during the period of 2002–2014, the mean three-dimensional structures of mesoscale eddies on both sides of the Luzon Strait (LS) were obtained via a composite method and analyzed to statistically examine the influences of background marine environment and the Kuroshio current on the eddy structures. The significant signals of temperature and salinity anomalies within the composite eddies extend much deeper in the region east of the LS (zone E) than those in the region west of the strait (zone W) because of stronger eddy intensity and larger vertical gradients of background temperature and salinity in the deep layer in zone E. In the vertical structure of temperature anomaly within the eddies, two cores occur at around 200 and 400 dbar depths, respectively, in zone E and only one core is centered at about 100 dbar in zone W. There is a clear three-core sandwich pattern in the vertical structure of salinity anomaly within the eddies in zone E. The Kuroshio water trapped in the eddy is responsible for abnormally positive salinity anomaly in the surface layer of the anticyclonic eddy center in zone W. On both sides of the LS, an asymmetric dipole structure in the surface layer gradually turns into a monopole one at depths, which resulted from the competition between horizontal advection effect and eddy pumping effect. The Kuroshio current influences the distribution patterns of isotherms and isohalines and enhances background temperature and salinity horizontal gradients on both sides of the LS, determining the orientations of dipole temperature and salinity structures within eddies.     

Toward a global ocean data assimilation system based on ensemble optimum interpolation: altimetry data assimilation experiment

A global ocean data assimilation system based on the ensemble optimum interpolation (EnOI) has been under development as the Chinese contribution to the Global Ocean Data Assimilation Experiment. The system uses a global ocean general circulation model, which is eddy permitting, developed by the Institute of Atmospheric Physics of the Chinese Academy of Sciences. In this paper, the implementation of the system is described in detail. We describe the sampling strategy to generate the stationary ensembles for EnOI. In addition, technical methods are introduced to deal with the requirement of massive memory space to hold the stationary ensembles of the global ocean. The system can assimilate observations such as satellite altimetry, sea surface temperature (SST), in situ temperature and salinity from Argo, XBT, Tropical Atmosphere Ocean (TAO), and other sources in a straightforward way. As a first step, an assimilation experiment from 1997 to 2001 is carried out by assimilating the sea level anomaly (SLA) data from TOPEX/Poseidon. We evaluate the performance of the system by comparing the results with various types of observations. We find that SLA assimilation shows very positive impact on the modeled fields. The SST and sea surface height fields are clearly improved in terms of both the standard deviation and the root mean square difference. In addition, the assimilation produces some improvements in regions where mesoscale processes cannot be resolved with the horizontal resolution of this model. Comparisons with TAO profiles in the Pacific show that the temperature and salinity fields have been improved to varying degrees in the upper ocean. The biases with respect to the independent TAO profiles are reduced with a maximum magnitude of about 0.25°C and 0.1 psu for the time-averaged temperature and salinity. The improvements on temperature and salinity also lead to positive impact on the subsurface currents. The equatorial under current is enhanced in the Pacific although it is still underestimated after the assimilation.     

Effects of cyclonic mesoscale eddies on the marine ecosystem in the Kuroshio Extension region using an eddy-resolving coupled physical-biological model

Effects of mesoscale eddies on the marine ecosystem in the Kuroshio Extension (KE) region are investigated using an eddy-resolving coupled physical-biological model. The model captures the seasonal and intra-seasonal variability of chlorophyll distribution associated with the mesoscale eddies, front variability, Kuroshio meanders, and upwelling. The model also reproduces the observed interannual variability of sea surface height anomaly (SSHA) in the KE region along a zonal band of 32–34°N from 2002 to 2006. The distribution of high surface chlorophyll corresponds to low SSHA. Cyclonic eddies are found to detach from the KE jet near 150°E and 158°E and propagate westward. The westward propagating cyclonic eddies lift the nutrient-rich thermocline into the euphotic zone and maintain high levels of chlorophyll in summer. In the subsurface layer, the pattern in chlorophyll is influenced by both lateral and vertical advection. In winter, convection inside the eddy entrains high levels of nutrients into the mixed layer, increasing production, and resulting in high chlorophyll concentration throughout the surface mixed layer. There is significant interannual variability in both the cyclonic eddy activity and the surface phytoplankton bloom south of the KE jet, although whether or not there is a causal link is unclear.     

High-resolution ensemble forecasting for the Gulf of Mexico eddies and fronts

High-resolution models and realistic boundary conditions are necessary to reproduce the mesoscale dynamics of the Gulf of Mexico (GOM). In order to achieve this, we use a nested configuration of the Hybrid Coordinate Ocean Model (HYCOM), where the Atlantic TOPAZ system provides lateral boundary conditions to a high-resolution (5 km) model of the GOM . However, such models cannot provide accurate forecasts of mesoscale variability, such as eddy shedding event, without data assimilation. Eddy shedding events involve the rapid growth of nonlinear instabilities that are difficult to forecast. The known sources of error are the initial state, the atmospheric condition, and the lateral boundary condition. We present here the benefit of using a small ensemble forecast (10 members) for providing confidence indices for the prediction, while using a data assimilation scheme based on optimal interpolation. Our set of initial states is provided by using different values of a data assimilation parameter, while the atmospheric and lateral boundary conditions are perturbed randomly. Changes in the data assimilation parameter appear to control the main position of the large features of the GOM in the initial state, whereas changes in the boundary conditions (lateral and atmospheric) appears to control the propagation of cyclonic eddies at their boundary. The ensemble forecast is tested for the shedding of Eddy Yankee (2006). The Loop Current and eddy fronts observed from ocean color and altimetry are almost always within the estimated positions from the ensemble forecast. The ensemble spread is correlated both in space and time to the forecast error, which implies that confidence indices can be provided in addition to the forecast. Finally, the ensemble forecast permits the optimization of a data assimilation parameter for best performance at a given forecast horizon.     

Existence of eddies at crossroad of the Indonesian seas

An eddy-resolving Indo-Pacific ocean circulation model was applied to highlight the behavior of eddies throughout the Indonesian seas. The complexity of the topography and coastline at the entrance of the Makassar Strait induce an eddy-type throughflow, instead of a straightforward flow. A sill and a narrow passage in the Makassar strait creates a barrier and impedes the continuation of eddies from the Pacific ocean, but the existence of a steep deep basin (>500 m depth) between the Java and Flores seas indicates a possible area for eddies. Based on our numerical results, we described the presence of a unique eddy structure north of Lombok Island, which we designated the “Lombok Eddy” and verified it by performing a drifter release field experiment and reviewing monthly mean climatology data from the World Ocean Atlas 2001 and XBT PX2 track data. NCEP/NCAR reanalysis, satellite observation data, and mixed layer depth analysis were also used to confirm these processes. By analyzing numerical simulation results and available temperature datasets, two additional eddies were found. All eddies form primarily due to eastward local winds correlated with seasonal monsoon winds during the austral summer. These eddies vary synchronously at an interannual time scale. Together, they are referred to as the Lombok Eddy Train (LET), which affects the surface layer down to a depth of 60 m, and the intensity of the eddy system is strongly affected by mixed layer depth variability from December to February.     

An assessment of the impact of surface currents on wave modeling in the Southern Ocean

This paper presents an assessment of the impact of the ocean circulation on modeled wave fields in the Southern Ocean, where a systematic positive bias of the modeled wave height against altimetry data has been reported. The inclusion of ocean currents in the wave model considerably reduces the positive bias of the simulated wave height for high southern latitudes. The decrease of wave energy in the presence of currents is almost exclusively related to the reduction of the relative wind, caused by an overall co-flowing current field associated with the Antarctic Circumpolar Current. Improvements of the model results are also found for the peak period and the mean period against a long-term moored buoy. At the mooring location, the effect of currents is greater for larger and longer waves, suggesting remotely generated swells are more influenced by the currents than local waves. However, an additional qualitative analysis using high-resolution currents in a finer grid nested to the global coarser grid shows that typical resolution of global hydrodynamic reanalysis is not sufficient to resolve mesoscale eddies, and as a consequence, the simulation of mesoscale wave patterns can be compromised. The results are also discussed in terms of the accuracy of forcing fields.     

Low-frequency variability of a semi-closed sea induced by the circulation in an adjacent ocean in a wind-driven,quasi-geostrophic,eddy-resolving simulation

The wind-driven circulation in the northwestern Pacific and the South China Sea (SCS) is simplified as a two-layer, quasi-geostrophic model in two rectangular basins connected by an idealized strait. This model is used to investigate the impact of the western boundary current (WBC) on the adjacent marginal sea. The variability of the circulation in the two basins is investigated with a high resolution and at low viscosity, which allows the numerical solution to resolve mesoscale eddy forcing. The model ocean is driven by the time-independent asymmetric wind stress acting on the idealized Pacific (large basin) only. Under the reference parameters used in this study, the WBC can intrude into the idealized SCS (small basin) in the form of a loop current, shedding eddies regularly. The rate of eddy shedding is nearly constant throughout the entire integration time of the model; however, the intensity of the eddy-shedding exhibits multiple timescale variability ranging from quasi-biennial to decadal timescale. A set of sensitivity experiments demonstrate that our results are robust against changes to model parameters and geometry. Multivariate spectral analysis is used to extract the spatiotemporal feature of the variability. Joint analysis for the two basins shows that the circulation in the idealized SCS is significantly impacted by the variability at decadal (15-year), interannual (5–7-year and quasi-biennial), and quasi-annual timescales. The spatial structures of the modes of variablility suggest that the variance in position of the WBC, combined with mesoscale activity, act to influence the low-frequency modes of the idealized SCS. The structural differences between the modes strongly impacting the idealized SCS and those having weak influence are also presented.     

Dynamic model of mesoscale eddies. Eddy parameterization for coarse resolution ocean circulation models

In the framework of the eddy dynamic model developed in two previous papers (Dubovikov, M.S., Dynamical model of mesoscale eddies, Geophys. Astophys. Fluid Dyn., 2003, 97, 311–358; Canuto, V.M. and Dubovikov, M.S., Modeling mesoscale eddies, Ocean Modelling, 2004, 8, 1–30 referred as I–II), we compute the contribution of unresolved mesoscale eddies to the large-scale dynamic equations of the ocean. In isopycnal coordinates, in addition to the bolus velocity discussed in I–II, the mesoscale contribution to the large scale momentum equation is derived. Its form is quite different from the traditional down-gradient parameterization. The model solutions in isopycnal coordinates are transformed to level coordinates to parameterize the eddy contributions to the corresponding large scale density and momentum equations. In the former, the contributions due to the eddy induced velocity and to the residual density flux across mean isopycnals (so called Σ-term) are derived, both contributions being shown to be of the same order. As for the large scale momentum equation, as well as in isopycnal coordinates, the eddy contribution has a form which is quite different from the down-gradient expression.     

Impacts of the Kuroshio intrusion on the two eddies in the northern South China Sea in late spring 2016

During mid-May to early June 2016, a cold eddy and a warm eddy were captured on the continental slope of the northern South China Sea by the in situ measurements. A salty lens-shaped water mass in the subsurface layer existed in these two detected eddies, which indicated they had a Kuroshio water origin. The trajectories of the observed eddies from satellite altimeter data show that the cold eddy was generated in the central part of the Luzon Strait, while the warm eddy was formed southwest of Taiwan. The genesis of the cold eddy is related to a weak Kuroshio loop current, while that of the warm eddy is associated with a strong Kuroshio loop current. The warm eddy east of the Luzon Strait may trigger the Kuroshio from a leaping path to a looping path. During the evolution of these detected eddies, they had interactions with the Kuroshio and Luzon Gyre. Energy analysis from ocean reanalysis data showed that the baroclinic conversion between the cold eddy and the Kuroshio was stronger than that between the cold eddy and Luzon Gyre. During the eddy shedding stage, the warm eddy mainly acquired energy from the Kuroshio loop current through the baroclinic conversion.     

Eddy energy sources and mesoscale eddies in the Sea of Okhotsk

Based on eddy-permitting ocean circulation model outputs, the mesoscale variability is studied in the Sea of Okhotsk. We confirmed that the simulated circulation reproduces the main features of the general circulation in the Sea of Okhotsk. In particular, it reproduced a complex structure of the East-Sakhalin current and the pronounced seasonal variability of this current. We established that the maximum of mean kinetic energy was associated with the East-Sakhalin Current. In order to uncover causes and mechanisms of the mesoscale variability, we studied the budget of eddy kinetic energy (EKE) in the Sea of Okhotsk. Spatial distribution of the EKE showed that intensive mesoscale variability occurs along the western boundary of the Sea of Okhotsk, where the East-Sakhalin Current extends. We revealed a pronounced seasonal variability of EKE with its maximum intensity in winter and its minimum intensity in summer. Analysis of EKE sources and rates of energy conversion revealed a leading role of time-varying (turbulent) wind stress in the generation of mesoscale variability along the western boundary of the Sea of Okhotsk in winter and spring. We established that a contribution of baroclinic instability predominates over that of barotropic instability in the generation of mesoscale variability along the western boundary of the Sea of Okhotsk. To demonstrate the mechanism of baroclinic instability, the simulated circulation was considered along the western boundary of the Sea of Okhotsk from January to April 2005. In April, the mesoscale anticyclonic eddies are observed along the western boundary of the Sea of Okhotsk. The role of the sea ice cover in the intensification of the mesoscale variability in the Sea of Okhotsk was discussed.     

Mesoscale eddies and baroclinic instability over the eastern Sakhalin shelf of the Sea of Okhotsk: a model-based analysis

Based on an eddy-permitting numerical model, the mesoscale variability in the East-Sakhalin Current is investigated during the winter-spring period. Analysis of necessary conditions for the development of baroclinic instability showed that the nearshore component of the East-Sakhalin Current is potentially baroclinic unstable in the first half-year. The simulated circulation uncovered a generation of anticyclonic eddies on the eastern Sakhalin shelf. It was established that a spatial scale of these eddies and the first baroclinic Rossby radius of deformation are values of the same order; a lifetime of these eddies varies from 4 to 6 weeks, given the Rossby number varies from 0.05 to 0.2. Analysis of the rate of eddy energy conversion on the eastern Sakhalin shelf showed that the generation of the revealed mesoscale eddies results from, mainly, baroclinic instability, whereas barotropic instability can be both favoring and preventing to the generation of these eddies.     

Numerical study of surface water circulation around Sekisei Lagoon, southwest Japan

This paper discusses the variability of surface currents around Sekisei Lagoon using a nested grid ocean circulation model. We developed a triple-nested grid system that consists of a coarse-resolution (1/60° or ∼1.85 km) model off Taiwan, an intermediate-resolution (1/300° or ∼370 m) model around the Yaeyama Islands, and a fine-resolution (1/900° or ∼123 m) model of Sekisei Lagoon. The nested grid system was forced by wind and heat flux calculated from six-hourly atmospheric reanalysis data and integrated over the period from May to July 2003. The coarse-resolution model was driven by lateral boundary conditions calculated from daily ocean reanalysis data to include realistic variation of the Kuroshio and mesoscale eddies with spatial scales of ∼500–700 km in the open ocean. The tidal forcing was included in the intermediate-resolution model by interpolating sea level data obtained from a data-assimilative tidal model. The results were then used to drive the fine-resolution model to simulate the surface water circulation around Sekisei lagoon. Model results show that (1) currents inside the lagoon are mainly driven by tide and wind; (2) there exists a strong southwestward current along the bottom slope in the southeast portion of the lagoon; the current is mainly driven by remote mesoscale eddies and at times intensified by the local wind; (3) the flow relaxation scheme is effective in reducing biases along the open boundaries. The simulated currents were used to examine the retention and dispersion of passive particles in the surface layer. Results show that the surface dispersion in the strong open ocean current region is significantly higher than that inside the lagoon.     

An eddy tracking algorithm based on dynamical systems theory

This work introduces a new method for ocean eddy detection that applies concepts from stationary dynamical systems theory. The method is composed of three steps: first, the centers of eddies are obtained from fixed points and their linear stability analysis; second, the size of the eddies is estimated from the vorticity between the eddy center and its neighboring fixed points, and, third, a tracking algorithm connects the different time frames. The tracking algorithm has been designed to avoid mismatching connections between eddies at different frames. Eddies are detected for the period between 1992 and 2012 using geostrophic velocities derived from AVISO altimetry and a new database is provided for the global ocean.