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.

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.

On the steadiness of separating meandering currents

The existence of inertial steady currents that separate from a coast and meander afterward is investigated. By integrating the zonal momentum equation over a suitable area, it is shown that retroflecting currents cannot be steady in a reduced gravity or in a barotropic model of the ocean. Even friction cannot negate this conclusion. Previous literature on this subject, notably the discrepancy between several articles by Nof and Pichevin on the unsteadiness of retroflecting currents and steady solutions presented in other papers, is critically discussed. For more general separating current systems, a local analysis of the zonal momentum balance shows that given a coastal current with a specific zonal momentum structure, an inertial, steady, separating current is unlikely, and the only analytical solution provided in the literature is shown to be inconsistent. In a basin-wide view of these separating current systems, a scaling analysis reveals that steady separation is impossible when the interior flow is nondissipative (e.g., linear Sverdrup-like). These findings point to the possibility that a large part of the variability in the world’s oceans is due to the separation process rather than to instability of a free jet.

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.

Observations of the inter-ocean exchange around South Africa

There is growing evidence that the interocean exchange south of Africa is an important link in the global overturning circulation of the ocean, the so‐called ocean conveyer belt. At this location, warm and salty Indian Ocean waters enter the South Atlantic and are pulled by currents that eventually reach the North Atlantic, where water cools and sinks. A major contributor to the exchange is the frequent shedding of ring eddies from the termination of the Agulhas Current south of the tip of Africa. This shedding is controlled by developments far upstream in the Indian Ocean, and variations in this ‘Agulhas Leakage’ can lead to changes in the rate and stability of the Atlantic overturning, with possible associated global climate variations [Weijer et al., 1999]. Regional climate variations in the tropical and subtropical Indian Ocean are known to affect the whole system of the Agulhas Current, including the interocean exchanges. This article reports on some of the seminal results of ongoing multinational, multidisciplinary projects that explore these issues.

Voltage-gated potassium currents within the dorsal vagal nucleus: inhibition by BDS toxin

Voltage-gated potassium (Kv) channels are essential components of neuronal excitability. The Kv3.4 channel protein is widely distributed throughout the central nervous system (CNS), where it can form heteromeric or homomeric Kv3 channels. Electrophysiological studies reported here highlight a functional role for this channel protein within neurons of the dorsal vagal nucleus (DVN). Current clamp experiments revealed that blood depressing substance (BDS) and intracellular dialysis of an anti-Kv3.4 antibody prolonged the action potential duration. In addition, a BDS sensitive, voltage-dependent, slowly inactivating outward current was observed in voltage clamp recordings from DVN neurons. Electrical stimulation of the solitary tract evoked EPSPs and IPSPs in DVN neurons and BDS increased the average amplitude and decreased the paired pulse ratio, consistent with a presynaptic site of action. This presynaptic modulation was action potential dependent as revealed by ongoing synaptic activity. Given the role of the Kv3 proteins in shaping neuronal excitability, these data highlight a role for homomeric Kv3.4 channels in spike timing and neurotransmitter release in low frequency firing neurons of the DVN.

Tracking the momentum flux of a CME and quantifying its influence on geomagnetically induced currents at Earth

We investigate a coronal mass ejection (CME) propagating toward Earth on 29 March 2011. This event is specifically chosen for its predominately northward directed magnetic field, so that the influence from the momentum flux onto Earth can be isolated. We focus our study on understanding how a small Earth-directed segment propagates. Mass images are created from the white-light cameras onboard STEREO which are also converted into mass height-time maps (mass J-maps). The mass tracks on these J-maps correspond to the sheath region between the CME and its associated shock front as detected by in situ measurements at L1. A time series of mass measurements from the STEREO COR-2A instrument is made along the Earth propagation direction. Qualitatively, this mass time series shows a remarkable resemblance to the L1 in situ density series. The in situ measurements are used as inputs into a three-dimensional (3-D) magnetospheric space weather simulation from the Community Coordinated Modeling Center. These simulations display a sudden compression of the magnetosphere from the large momentum flux at the leading edge of the CME, and predictions are made for the time derivative of the magnetic field (dB/dt) on the ground. The predicted dB/dt values were then compared with the observations from specific equatorially located ground stations and showed notable similarity. This study of the momentum of a CME from the Sun down to its influence on magnetic ground stations on Earth is presented as a preliminary proof of concept, such that future attempts may try to use remote sensing to create density and velocity time series as inputs to magnetospheric simulations.

Cusp ion steps, field-aligned currents and poleward moving auroral forms

We predict the field-aligned currents around cusp ion steps produced by pulsed reconnection between the geomagnetic field and an interplanetary magnetic field (IMF) with a B-Y component that is large in magnitude. For B-Y > 0, patches of newly opened flux move westward and eastward in the Northern and Southern Hemispheres, respectively, under the influence of the magnetic curvature force. These flow directions are reversed for B-Y < 0. The speed of this longitudinal motion initially grows with elapsed time since reconnection, but then decays as the newly opened field lines straighten. We predict sheets of field-aligned current on the boundaries between the patches produced by successive reconnection pulses, associated with the difference in the speeds of their longitudinal motion. For low elapsed times since reconnection, near the equatorward edge of the cusp region where the field lines are accelerating, the field-aligned current sheets will be downward or upward in both hemispheres for positive or negative IMF B-Y, respectively. At larger elapsed times since reconnection, as events slow and evolve from the cusp into the mantle region, these field-aligned current directions will be reversed. Observations by the Polar spacecraft on August 26,1998, show the predicted upward current sheets at steps seen in the mantle region for IMF B-Y > 0. Mapped into the ionosphere, the steps coincide with poleward moving events seen by the CUTLASS HF radar. The mapped location of the largest step also coincides with a poleward moving arc seen by the UVI imager on Polar. We show that the arc is consistent with a region of upward field-aligned current that has become unstable, such that a potential drop of about 1 kV formed below the spacecraft. The importance of these observations is that they confirm that the poleward moving events, as seen by the HF radar and the UV imager, are due to pulsed magnetopause reconnection. Milan et al. [2000] noted that the great longitudinal extent of these events means that the required reconnection pulses would have contributed almost all the voltage placed across the magnetosphere at this time. The observations also show that auroral arcs can form on open field lines in response to the pulsed application of voltage at the magnetopause.

Coherent EISCAT Svalbard Radar spectra from the dayside cusp/cleft and their implications for transient field-aligned currents

Naturally enhanced incoherent scatter spectra from the vicinity of the dayside cusp/cleft, interpreted as being due to plasma turbulence driven by short bursts of intense field-aligned current, are compared with high-resolution narrow-angle auroral images and meridian scanning photometer data. Enhanced spectra have been observed on many occasions in association with nightside aurora, but there has been only one report of such spectra seen in the cusp/cleft region. Narrow-angle images show considerable change in the aurora on timescales shorter than the 10-s radar integration period, which could explain spectra observed with both ion lines simultaneously enhanced. Enhanced radar spectra are generally seen inside or beside regions of 630-nm auroral emission, indicative of sharp F region conductivity gradients, but there appears also to be a correlation with dynamic, small-scale auroral forms of order 100 m and less in width.

DE-2 observations of filamentary currents at ionospheric altitudes

Conjunctive measurements made by the Dynamics Explorer 1 and 2 spacecraft on October 22, 1981, under conditions of southward IMF, suggest the existence of a cusp ion injection from a region at the magnetopause with a scale size of ∼ 1/2 to 1 R E . Current signatures observed by the LAPI and MAGB instruments on board DE-2 indicate the existence of a rotation in the magnetic field that is consistent with a filamentary current system. The observed current structure can be interpreted as the ionospheric signature of a flux transfer event (FTE). In addition to this large-scale current structure there exist three small-scale filamentary current pairs. These current pairs close locally and thus, if our interpretation of this event as an FTE is correct, represent the first reported observations of FTE interior structure at low-altitudes.

Deforestation, leakage and avoided deforestation policies: a spatial analysis

This paper analyses the impact of several avoided deforestation policies within a patchy forested landscape. Central is the idea that deforestation choices in one area influence deforestation decisions in nearby patches. We explore the interplay between forest landscapes comprising heterogeneous patches, localised spatial displacement, and avoided deforestation policies. Avoided deforestation policies at a landscape level are respectively: two Payments for Environmental Services (PES) policies, one focused on deforestation hotspots, the second being equally available to all agents; a conservation area; and, an agglomeration bonus. We demonstrate how the "best" policy, in terms of reduced leakage, depends on landscape heterogeneity. Agglomeration bonuses are shown to be more effective where there is less landscape heterogeneity, whilst conservation areas are most effective where there is more spatial heterogeneity.

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%.