The influence of bottom topography on the decay of modeled Agulhas rings

The influence of a large meridional submarine ridge on the decay of Agulhas rings is investigated with a 1 and 2-layer setup of the isopycnic primitive-equation ocean model MICOM. In the single-layer case we show that the SSH decay of the ring is primarily governed by bottom friction and secondly by the radiation of Rossby waves. When a topographic ridge is present, the effect of the ridge on SSH decay and loss of tracer from the ring is negligible. However, the barotropic ring cannot pass the ridge due to energy and vorticity constraints. In the case of a two-layer ring the initial SSH decay is governed by a mixed barotropic–baroclinic instability of the ring. Again, radiation of barotropic Rossby waves is present. When the ring passes the topographic ridge, it shows a small but significant stagnation of SSH decay, agreeing with satellite altimetry observations. This is found to be due to a reduction of the growth rate of the m = 2 instability, to conversions of kinetic energy to the upper layer, and to a decrease in Rossby-wave radiation. The energy transfer is related to the fact that coherent structures in the lower layer cannot pass the steep ridge due to energy constraints. Furthermore, the loss of tracer from the ring through filamentation is less than for a ring moving over a flat bottom, related to a decrease in propagation speed of the ring. We conclude that ridges like the Walvis Ridge tend to stabilize a multi-layer ring and reduce its decay.

Flux comparison of Eulerian and Lagrangian estimates of Agulhas leakage: A case study using a numerical model

Estimating the magnitude of Agulhas leakage, the volume flux of water from the Indian to the Atlantic Ocean, is difficult because of the presence of other circulation systems in the Agulhas region. Indian Ocean water in the Atlantic Ocean is vigorously mixed and diluted in the Cape Basin. Eulerian integration methods, where the velocity field perpendicular to a section is integrated to yield a flux, have to be calibrated so that only the flux by Agulhas leakage is sampled. Two Eulerian methods for estimating the magnitude of Agulhas leakage are tested within a high-resolution two-way nested model with the goal to devise a mooring-based measurement strategy. At the GoodHope line, a section halfway through the Cape Basin, the integrated velocity perpendicular to that line is compared to the magnitude of Agulhas leakage as determined from the transport carried by numerical Lagrangian floats. In the first method, integration is limited to the flux of water warmer and more saline than specific threshold values. These threshold values are determined by maximizing the correlation with the float-determined time series. By using the threshold values, approximately half of the leakage can directly be measured. The total amount of Agulhas leakage can be estimated using a linear regression, within a 90% confidence band of 12 Sv. In the second method, a subregion of the GoodHope line is sought so that integration over that subregion yields an Eulerian flux as close to the float-determined leakage as possible. It appears that when integration is limited within the model to the upper 300 m of the water column within 900 km of the African coast the time series have the smallest root-mean-square difference. This method yields a root-mean-square error of only 5.2 Sv but the 90% confidence band of the estimate is 20 Sv. It is concluded that the optimum thermohaline threshold method leads to more accurate estimates even though the directly measured transport is a factor of two lower than the actual magnitude of Agulhas leakage in this model.

On the fast decay of Agulhas rings

The Indian Ocean water that ends up in the Atlantic Ocean detaches from the Agulhas Current retroflection predominantly in the form of Agulhas rings and cyclones. Using numerical Lagrangian float trajectories in a high-resolution numerical ocean model, the fate of coherent structures near the Agulhas Current retroflection is investigated. It is shown that within the Agulhas Current, upstream of the retroflection, the spatial distributions of floats ending in the Atlantic Ocean and floats ending in the Indian Ocean are to a large extent similar. This indicates that Agulhas leakage occurs mostly through the detachment of Agulhas rings. After the floats detach from the Agulhas Current, the ambient water quickly looses its relative vorticity. The Agulhas rings thus seem to decay and loose much of their water in the Cape Basin. A cluster analysis reveals that most water in the Agulhas Current is within clusters of 180 km in diameter. Halfway in the Cape Basin there is an increase in the number of larger clusters with low relative vorticity, which carry the bulk of the Agulhas leakage transport through the Cape Basin. This upward cascade with respect to the length scales of the leakage, in combination with a power law decay of the magnitude of relative vorticity, might be an indication that the decay of Agulhas rings is somewhat comparable to the decay of two-dimensional turbulence.

Relating Agulhas leakage to the Agulhas Current retroflection location

The relation between the Agulhas Current retroflection location and the magnitude of Agulhas leakage, the transport of water from the Indian to the Atlantic Ocean, is investigated in a high-resolution numerical ocean model. Sudden eastward retreats of the Agulhas Current retroflection loop are linearly related to the shedding of Agulhas rings, where larger retreats generate larger rings. Using numerical Lagrangian floats a 37 year time series of the magnitude of Agulhas leakage in the model is constructed. The time series exhibits large amounts of variability, both on weekly and annual time scales. A linear relation is found between the magnitude of Agulhas leakage and the location of the Agulhas Current retroflection, both binned to three month averages. In the relation, a more westward location of the Agulhas Current retroflection corresponds to an increased transport from the Indian Ocean to the Atlantic Ocean. When this relation is used in a linear regression and applied to almost 20 years of altimetry data, it yields a best estimate of the mean magnitude of Agulhas leakage of 13.2 Sv. The early retroflection of 2000, when Agulhas leakage was probably halved, can be identified using the regression.

Lagrangian validation of numerical drifter trajectories using drifting buoys: Application to the Agulhas system

The skill of numerical Lagrangian drifter trajectories in three numerical models is assessed by comparing these numerically obtained paths to the trajectories of drifting buoys in the real ocean. The skill assessment is performed using the two-sample Kolmogorov–Smirnov statistical test. To demonstrate the assessment procedure, it is applied to three different models of the Agulhas region. The test can either be performed using crossing positions of one-dimensional sections in order to test model performance in specific locations, or using the total two-dimensional data set of trajectories. The test yields four quantities: a binary decision of model skill, a confidence level which can be used as a measure of goodness-of-fit of the model, a test statistic which can be used to determine the sensitivity of the confidence level, and cumulative distribution functions that aid in the qualitative analysis. The ordering of models by their confidence levels is the same as the ordering based on the qualitative analysis, which suggests that the method is suited for model validation. Only one of the three models, a 1/10° two-way nested regional ocean model, might have skill in the Agulhas region. The other two models, a 1/2° global model and a 1/8° assimilative model, might have skill only on some sections in the region

A weaker Agulhas Current leads to more Agulhas leakage

Time series of transports in the Agulhas region have been constructed by simulating Lagrangian drifter trajectories in a 1/10 degree two-way nested ocean model. Using these 34 year long time series it is shown that smaller (larger) Agulhas Current transport leads to larger (smaller) Indian-Atlantic inter-ocean exchange. When transport is low, the Agulhas Current detaches farther downstream from the African continental slope. Moreover, the lower inertia suppresses generation of anti-cyclonic vorticity. These two effects cause the Agulhas retroflection to move westward and enhance Agulhas leakage. In the model a 1 Sv decrease in Agulhas Current transport at 32°S results in a 0.7 ± 0.2 Sv increase in Agulhas leakage.

Fast northward energy transfer in the Atlantic due to Agulhas rings

The adiabatic transit time of wave energy radiated by an Agulhas ring released in the South Atlantic Ocean to the North Atlantic Ocean is investigated in a two-layer ocean model. Of particular interest is the arrival time of baroclinic energy in the northern part of the Atlantic, because it is related to variations in the meridional overturning circulation. The influence of the Mid-Atlantic Ridge is also studied, because it allows for the conversion from barotropic to baroclinic wave energy and the generation of topographic waves. Barotropic energy from the ring is present in the northern part of the model basin within 10 days. From that time, the barotropic energy keeps rising to attain a maximum 500 days after initiation. This is independent of the presence or absence of a ridge in the model basin. Without a ridge in the model, the travel time of the baroclinic signal is 1300 days. This time is similar to the transit time of the ring from the eastern to the western coast of the model basin. In the presence of the ridge, the baroclinic signal arrives in the northern part of the model basin after approximately 10 days, which is the same time scale as that of the barotropic signal. It is apparent that the ridge can facilitate the energy conversion from barotropic to baroclinic waves and the slow baroclinic adjustment can be bypassed. The meridional overturning circulation, parameterized in two ways as either a purely barotropic or a purely baroclinic phenomenon, also responds after 1300 days. The ring temporarily increases the overturning strength. Th presence of the ridge does not alter the time scales.

Tracer Leakage from Modeled Agulhas Rings

In a numerical, isopycnal, ocean model the mixing is investigated with the environment of two idealized Agulhas rings, one that splits, and one that remains coherent. The evolution of a passive tracer , initially contained within the rings, shows that tracer leakage is associated with the for mation of filaments in the early stage of ring evolution. These filaments reach down to the ther mocline. In the deepest layers leakage occurs on a larger scale. Self-advection of the rings is ver y irregular , and it is not possible to compute a Lagrangian boundar y i n order to estimate the transport of leakage from the rings. T o describe the processes that gover n tracer leakage, in a coordinate frame moving with the ring a kinematic separatrix is defined in the streamfunction field for the nondivergent flow . Initially , filaments arise because of the elongation of the ring, which is mainly gover ned by an m 5 2 instability that is collaborating with differential rotation. Because of beta, the symmetr y i s destroyed related to the separatrix associated with a stagnation point in the flow . The filament upstream of the stagnation point grows much faster and is associated with the bulk of tracer leakage. Mixing is enhanced by time dependence of the separatrix. As a result, there are no large differences between the leakage from a coherent ring, where the m 5 2 instability equilibrates, and from a splitting ring, where the m 5 2 instability keeps growing, which confir ms that the amount of leakage is mainly gover ned by the ring’ s initial defor mation combined with unsteady self-advection of the ring and not by the splitting of the ring. The decay of tracer content in the ther mocline shows that in the first months up to 40% of the ring water can be mixed with the environment. In deeper layers the decay of tracer content may reach up to 90%.

Generation and Evolution of Natal Pulses: Solitary Meanders in the Agulhas Current

Solitar y meanders of the Agulhas Current, so-called Natal pulses, may play an important role in the overall dynamics of this current system. Several hypotheses concer ning the triggering of these pulses are tested using sea sur face height and temperature data from satellites. The data show the for mation of pulses in the Natal Bight area at irregular inter vals ranging from 50 to 240 days. Moving downstream at speeds between 10 and 20 km day 2 1 they sometimes reach sizes of up to 300 km. They seem to play a role in the shedding of Agulhas rings that penetrate the South Atlantic. The inter mittent for mation of these solitar y meanders is argued to be most probably related to barotropic instability of the strongly baroclinic Agulhas Current in the Natal Bight. The vorticity structure of the obser ved basic flow is argued to be stable anywhere along its path. However , a proper perturbation of the jet in the Natal Bight area will allow barotropic instability , because the bottom slope there is considerably less steep than elsewhere along the South African east coast. Using satellite altimetr y these perturbations seem to be related to the inter mittent presence of offshore anticyclonic anomalies, both upstream and eastward of the Natal Bight.

Assimilation of Geosat Altimeter Data for the Agulhas Current Using the Ensemble Kalman Filter with a Quasigeostrophic Model

The ring-shedding process in the Agulhas Current is studied using the ensemble Kalman filter to assimilate geosat altimeter data into a two-layer quasigeostrophic ocean model. The properties of the ensemble Kalman filter are further explored with focus on the analysis scheme and the use of gridded data. The Geosat data consist of 10 fields of gridded sea-surface height anomalies separated 10 days apart that are added to a climatic mean field. This corresponds to a huge number of data values, and a data reduction scheme must be applied to increase the efficiency of the analysis procedure. Further, it is illustrated how one can resolve the rank problem occurring when a too large dataset or a small ensemble is used.

Translation, decay and splitting of Agulhas rings in the southeastern Atlantic Ocean

All Agulhas rings that were spawned at the Agulhas retrofiec- tion between 1993 and 1996 (a total of 21 rings) have been monitored using TOPEX/Poseidon satellite altimetry and followed as they moved through the southeastern Atlantic Ocean, decayed, interacted with bottom topography and each other, or dissipated completely. Rings preferentially crossed the Walvis Ridge at its deepest parts. After having crossed this ridge they have lower translational speeds, and their decay rate decreases markedly. Half the decay of long-lived rings takes place in the first 5 months of their lifetimes. In addition to the strong decay of rings in the Cape Basin, about one third of the observed rings do not seem to leave this region at all but totally disintegrate here. The interaction of rings with bottom topography, in particular with the Verna Seamount, is shown frequently to cause splitting of rings. This will enhance mixing of the rings' Indian Ocean water into that of the southern Atlantic. This localized mixing may well provide a considerable source of warm and salty Indian Ocean water into the Atlantic overturning circulation.

Upstream control of Agulhas Ring shedding

Rings shed in the Agulhas Retroflection region play an important role in the global thermohaline circulation. The shedding of these rings has been considered very irregular. In this paper, we present evidence for remote control of the timing and frequency of the ring shedding events. This turns out to be a far more regular process, at a frequency of 4–5 cycles per year. The movement of the Agulhas Retroflection, and thereby the shedding of rings, is timed by incoming eddies from the upstream regions. Eddies from the Mozambique Channel, and from the East Madagascar current reach the retroflection region at the frequency of 4–5 times per year. The existence of these eddies can be related to incoming Rossby waves that cross the Indian Ocean and reach the Agulhas Current system. These may in turn be part of a basin-wide oscillation. The irregularity found in ring shedding statistics can be ascribed to processes occurring between the actual shedding and the first unambiguous observation of a separated ring.

Observation and origin of an interannual salinity anomaly in the Mozambique Channel

A positive salinity anomaly of 0.2 PSU was observed between 50 and 200 m over the years 2000–2001 across the Mozambique Channel at a section at 17°S which was repeated in 2003, 2005, 2006, and 2008. Meanwhile, a moored array is continued from 2003 to 2008. This anomaly was most distinct showing an interannual but nonseasonal variation. The possible origin of the anomaly is investigated using output from three ocean general circulation models (Estimating the Circulation and Climate of the Ocean, Ocean Circulation and Climate Advanced Modeling, and Parallel Ocean Program). The most probable mechanism for the salinity anomaly is the anomalous inflow of subtropical waters caused by a weakening of the northern part of the South Equatorial Current by weaker trade winds. This mechanism was found in all three numerical models. In addition, the numerical models indicate a possible salinization of one of the source water masses to the Mozambique Channel as an additional cause of the anomaly. The anomaly propagated southward into the Agulhas Current and northward along the African coast.