Satellite-derived Sea Surface Height Research Articles

The Mean Upper-Layer Flow in the Central Gulf of Mexico by a New Method

AbstractPrevious studies have found a puzzling disagreement between two large datasets and the results of numerical models in the central Gulf of Mexico. The observations suggest an upper-layer mean flow to the west of order 10 cm s−1, while the numerical models find no such mean flow. A new technique is used here, using 23 yr of satellite-derived sea surface height data, to estimate the mean flow. This third, independent set of data yields the same westward flow found in previous studies. These findings require that there be sinking in the western Gulf. The details of the return flow remain an intriguing problem.

Response of micro- and mesozooplankton in the southwestern Bay of Bengal to a cyclonic eddy during the winter monsoon, 2005.

A cyclonic eddy (CE) in the southwestern Bay of Bengal (SWBoB; 10-15° N; 81-87° E) during winter monsoon 2005 and associated changes in the open ocean hydrography and productivity patterns were studied using satellite observations and in situ measurements. Analysis of the satellite-derived sea surface height anomaly (SSHA) indicated the existence of a large eddy (10-15° N; 81-87° E) from November to January, with its core centered at 13° N and 82° E. The large positive wind stress curl (~1.8 × 10(-7) N m(-2)) and resultant Ekman pumping (~3 × 10(-2) cm s(-1)) were identified as the prominent forcing mechanisms. In addition, the cyclonic storms and depressions experienced in the region during the study period seem to have served to maintain the strength of the CE through associated Ekman pumping. The cold (~26.6 °C), nutrient-enriched (NO3 > 2 μM, PO4 > 0.73 μM and SiO4 > 3 μM) upwelled waters in the upper layers of the CE enhanced the biological production (chl. a > 0.56 mg m(-3)). Dissolved oxygen in the surface waters was > 200 μM. The phytoplankton and zooplankton biomass recorded during the season was comparable or perhaps higher than the peak values reported from the northeastern Arabian Sea (chlorophyll a concentrations of 0.2 to 0.4 mg m(-3) and zooplankton biovolume 0.6 ml m(-3)) during winter. Occurrence of a higher mesozooplankton biovolume (0.8 ml m(-3)) and relatively low abundance of microzooplankton indicates the prevalence of a short food chain. In conclusion, high biological production, both at primary as well as secondary level, suggests the prevalence of an efficient food web as a result of physical forcing and subsequent ecological interactions evident up to tertiary level in an oligotrophic basin like BoB.

Variability of satellite-derived sea surface height anomaly, and its relationship with Bigeye tuna (Thunnus obesus) catch in the Eastern Indian Ocean

We analyzed the variability of sea surface height anomaly (SSHA), and its relationship with Bigeye tuna catch in the eastern Indian Ocean (EIO) off of Java Island (Indonesia). Both time series of SSHA and Bigeye tuna HR show dominant signals corresponding to the annual and inter-annual variability. During the southeast monsoon the wind blows along southern coast of Java and produces coastal upwelling. This causes sea level to drop due to an offshore Ekman transport, and thermocline becomes shallower. During El Niño and Indian Ocean Dipole (IOD) positive phase, upwelling is more intense and a large cold eddy forms in the EIO off Java. Generally, Bigeye tuna HR tends to increase during upwelling seasons and becomes even higher during El Niño and the positive phase of the IOD. The increased Bigeye tuna HR during the southeast monsoon, El Niño and the IOD positive phase can be attributed to the shallower thermocline depth and the enhancement of biological productivity due to development of eddies and strong upwelling in the EIO. The spatial distribution of SSHA indicates that Bigeye tuna catches are abundant in the frontal regions between cold and warm eddies.

Observed intraseasonal thermocline variability in the Bay of Bengal

The time series of temperature data obtained from moored buoys deployed at 8°N, 12°N, and 15°N along 90°E in the Bay of Bengal (BoB) shows a persistent intraseasonal variability on 30-120 day time scale in three distinct periods 30-70 day, near 90 day, and near 120 day in the thermocline region. The standard deviation of moored buoy temperature data shows that half of the variability in the thermocline region is contributed from the 30-120 day variability. The relative contribution of local Ekman pumping velocity and remote wind forcing from equatorial Indian Ocean (EIO) to the intraseasonal thermocline variability in the BoB is examined using satellite-derived sea surface height anomaly (SSHA), wind and depths of 23° isotherm (D23, proxy for thermocline depth) derived from moored buoys temperature data. The analysis shows that large amplitude intraseasonal oscillations of thermocline - particularly the near 90 day and 120 day variability - could not be explained by local Ekman pumping velocity alone. The SSHA, D23, and wind fields reveal that the first and second baroclinic mode Kelvin and Rossby waves, which are generated remotely by winds from the EIO and eastern BoB, can significantly influence the thermocline variability in the BoB. The near 90 day and 120 day thermocline variability is driven primarily by the variability of equatorial zonal wind stress. While the 30-70 day thermocline variability is affected most by interior Ekman pumping over the Bay, it also appears to be influenced by zonal wind stress in the EIO and alongshore wind stress in the eastern BoB.

Episodic phytoplankton bloom events in the Bay of Bengal triggered by multiple forcings

Two episodic phytoplankton bloom events were identified in the southeastern and southwestern Bay of Bengal in January 1998 and December 2005, respectively, from SeaWiFS-derived chlorophyll-a concentration data collected between September 1997 and December 2010. Possible causes were examined using time series datasets of satellite-derived sea surface height anomalies, sea surface temperature, wind stress and Ekman pumping velocity data. It was found that strong cyclonic eddies with long residence time transported cold, nutrient-rich water from the subsurface layer to the euphotic zone, and were responsible for sustaining the continuous growth of phytoplankton that led to the two aforementioned phytoplankton blooms. Moreover, regarding the event in December 2005, three tropical cyclones passed through the region of the phytoplankton bloom and could have been additional factors that enhanced the bloom. Regarding the event in January 1998, transportation from the adjacent Andaman Sea by the geostrophic current and an anomalous zonal easterly wind contributed to the phytoplankton bloom in the southeastern Bay of Bengal. In addition, both blooms occurred during Indian Ocean Dipole (IOD) years. The absence of the second downwelling Kelvin waves in positive IOD year of 1997–1998 favouring the long residence of the strong cyclonic eddy in the eastern bay, and the strong wind-induced mixing in negative IOD year of 2005 might contribute to enhancing productivity.

A Neural Network Approach to Estimate Tropical Cyclone Heat Potential in the Indian Ocean

The tropical cyclone heat potential (TCHP) or the available upper ocean thermal energy is one of the critical factors in controlling the intensity of cyclones. Given the devastating impacts Indian Ocean cyclones could bring (e.g., the “killer cyclone” Nargis in 2008, which caused more than 130000 deaths), there is a pressing need to obtain reliable and more accurate TCHP estimates over the Indian Ocean to improve the cyclone track and intensity predictions. Using more than 25000 in situ subsurface temperature profiles during 1997-2007, this research explores the possibility of developing an artificial neural network (ANN) model to derive TCHP in the Indian Ocean using satellite-derived sea surface height anomalies, sea surface temperature, and climatological depth of 26°C isotherm. The estimations have been validated using more than 8000 independent in situ profiles during 2008-2009. The root-mean-square error and the scatter index of this validation data sets are 14.6 kJ/cm <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sup> and 0.2, respectively. Comparison of the estimations from a two-layer reduced gravity model and from a multiple regression method confirms the superiority of the ANN approach over other methods.

Analysis of horizontal and vertical processes contributing to natural iron supply in the mixed layer in southern Drake Passage

Horizontal advection, vertical mixing, and mixed-layer entrainment all affect iron concentrations and biological productivity in the Ona Basin, near the Shackleton Transverse Ridge (STR) in southern Drake Passage. Trace metal sampling in the region indicates that dissolved iron concentrations are significantly higher on the continental shelf near the Antarctic Peninsula and the South Shetland Islands than they are in the deep waters away from the shelf. Comparisons between satellite-derived sea surface height (SSH) and Chlorophyll-a (Chl-a) levels in the Ona Basin show >95% correlation between Chl-a concentrations and horizontal advection of these iron-rich shelf waters during the months of November and December for the years 1997–2010. However, no significant correlations are found for January–April, while high Chl-a concentrations in the Ona Basin persist through March. Enhanced vertical (diapycnal) mixing and mixed-layer entrainment are considered as alternative mechanisms for delivering iron into the Ona Basin mixed layer and sustaining the high Chl-a concentrations. Estimates of iron flux based on in situ measurements of dissolved iron concentrations suggest that diapycnal mixing alone can supply iron to the base of the mixed layer at a rate of 64±2nmolm−2 day−1 during the summer. In addition, the summer mixed layer in the Ona Basin deepens from January to April, allowing for iron-rich water to be steadily entrained from below. Estimates based on monthly mixed-layer climatologies produce average daily entrainment rates ranging from 5 to 25nmolm−2 day−1. While neither diapycnal mixing nor entrainment alone is always sufficient to meet the estimated iron demand for the Ona Basin bloom, numerical simulation suggests that the combined effect of the two processes can consistently supply sufficient iron to sustain the bloom.

Surface Geostrophic Circulation of the Mediterranean Sea Derived from Drifter and Satellite Altimeter Data

AbstractDrifter observations and satellite-derived sea surface height data are used to quantitatively study the surface geostrophic circulation of the entire Mediterranean Sea for the period spanning 1992–2010. After removal of the wind-driven components from the drifter velocities and low-pass filtering in bins of 1° × 1° × 1 week, maps of surface geostrophic circulation (mean flow and kinetic energy levels) are produced using the drifter and/or satellite data. The mean currents and kinetic energy levels derived from the drifter data appear stronger/higher with respect to those obtained from satellite altimeter data. The maps of mean circulation estimated from the drifter data and from a combination of drifter and altimeter data are, however, qualitatively similar. In the western basin they show the main pathways of the surface waters flowing eastward from the Strait of Gibraltar to the Sicily Channel and the current transporting waters back westward along the Italian, French, and Spanish coasts. Intermittent and long-lived subbasin-scale eddies and gyres abound in the Tyrrhenian and Algerian Seas. In the eastern basin, the surface waters are transported eastward by several currents but recirculate in numerous eddies and gyres before reaching the northward coastal current off Israel, Lebanon, and Syria and veering westward off Turkey. In the Ionian Sea, the mean geostrophic velocity maps were also produced separately for the two extended seasons and for multiyear periods. Significant variations are confirmed, with seasonal reversals of the currents in the south and changes of the circulation from anticyclonic (prior to 1 July 2007) to cyclonic and back to anticyclonic after 31 December 2005.

Observation of first and second baroclinic mode Yanai waves in the ocean

AbstractThe structure and propagation of oceanic Yanai waves (also known as mixed Rossby–gravity waves) are investigated by analysis of satellite‐derived sea‐surface‐height (SSH) data. Significant spectral peaks along the dispersion curves of first and second baroclinic mode Yanai waves are identified by the wavenumber–frequency spectral analysis of SSH for the period 10–17 days. The spatial structure and propagation of these modes are described by an analysis of SSH time series filtered in the frequency–wavenumber domain that includes strong signals along the dispersion curves but excludes those of 17‐day tropical instability waves. The difference in meridional structure between the first and second baroclinic modes, which is consistent with theory, is evident. The slow eastward propagation of first and second baroclinic mode Yanai wave packets, which is consistent with the group velocity derived from theory, is also isolated from the SSH data. Year‐to‐year variability of Yanai wave activity is compared with that of the El Niño–Southern Oscillation (ENSO), showing a moderate correlation with ENSO Modoki. © Published in 2011 by John Wiley &amp; Sons Ltd.

The role of meso-scale eddies in mixed layer deepening and mode water formation in the western North Pacific

Distributions of mixed layer depths around the centers of anti-cyclonic and cyclonic eddies in the North Pacific Ocean were composited by using satellite-derived sea surface height anomaly data and Argo profiling float data. The composite distributions showed that in late winter, deeper mixed layers were more (less) frequently observed inside the cores of the anti-cyclonic (cyclonic) eddies than outside. This relationship was the clearest in the region of 140°E–160°W and 35°N–40°N, where the temperature and salinity of the deep mixed layers were similar to those of the lighter variety of central mode water (L-CMW). A simple one-dimensional bulk mixed layer model showed that both strong sea-surface heat and momentum fluxes and weak preexisting stratification contributed to formation of the deep mixed layer. These conditions were associated with the anti-cyclonic eddies, suggesting that these eddies are important in the formation of mode waters, particularly L-CMW.

Interannual and seasonal variations of Kuroshio transport east of Taiwan inferred from 29 years of tide-gauge data

[1] Twenty-nine years of tide-gauge data are analyzed in conjunction with wind and satellite-derived sea-surface height and ocean velocity data to study the interannual and seasonal variations of the Kuroshio transport off the northeastern coast of Taiwan. The data reveals an interannual variation of ±0.1 m (transport-variation of approximately ±3.5 Sv; 1 Sv = 106 m3 s−1), and a much weaker (5–10 times weaker) seasonal fluctuation that is minimum in May and maximum in November. The interannual fluctuations are not directly wind-driven by linear dynamics; rather, the Kuroshio strengthens in years of abundant eddies of the Subtropical Counter Current, which is related to the current's instability state driven by the slow fluctuations of the large-scale wind stress curl in the western Pacific. The seasonal transport fluctuation is also eddy-forced, but has weaker amplitude because the seasonal time scale is of the same order as the eddy-propagation time scale, and transport-producing eddy signals tend to overlap east of Taiwan.

Oceanic Heat Content Variability in the Eastern Pacific Ocean for Hurricane Intensity Forecasting

Abstract Recent evidence supports the premise that the subsurface ocean structure plays an important role in modulating air–sea fluxes during hurricane passage, which in turn, affects intensity change. Given the generally sparse in situ data, it has been difficult to provide region-to-basin-wide estimates of isotherm depths and upper-ocean heat content (OHC). In this broader context, satellite-derived sea surface height anomalies (SSHAs) from multiple platforms carrying radar altimeters are blended, objectively analyzed, and combined with a hurricane-season climatology to estimate isotherm depths and OHC within the context of a reduced gravity model at 0.25° spatial intervals in the eastern Pacific Ocean where tropical cyclone intensity change occurs. Measurements from the Eastern Pacific Investigation of Climate in 2001, long-term tropical ocean atmosphere mooring network, and volunteer observing ship deploying expendable bathythermograph (XBT) profilers are used to carefully evaluate satellite-based measurements of upper-ocean variability. Regression statistics reveal small biases with slopes of 0.8–0.9 between the subsurface measurements compared with isotherm depths (20° and 26°C), and OHC fields derived from objectively analyzed SSHA field. Root-mean-square differences in OHC range between 10 and 15 kJ cm−2 or roughly 10%–15% of the mean signals. Similar values are found for isotherm depth differences between in situ and inferred satellite-derived values. Blended daily values are used in the Statistical Hurricane Intensity Prediction Scheme (SHIPS) forecasts as are OHC estimates for the Atlantic Ocean basin. An equivalent OHC variable is introduced that incorporates the strength of the thermocline at the base of the oceanic mixed layer using a climatological stratification parameter /No, which seems better correlated to hurricane intensity change than just anomalies as observed in Hurricane Juliette in 2001.

Interannual variability of Kelvin wave propagation in the wave guides of the equatorial Indian Ocean, the coastal Bay of Bengal and the southeastern Arabian Sea during 1993–2006

The observed variability of the Kelvin waves and their propagation in the equatorial wave guide of the Indian Ocean and in the coastal wave guides of the Bay of Bengal (BoB) and the southeastern Arabian Sea (AS) on seasonal to interannual time scales during years 1993–2006 is examined utilizing all the available satellite and in-situ measurements. The Kelvin wave regime inferred from the satellite-derived sea surface height anomalies (SSHA) shows a distinct annual cycle composed of two pairs of alternate upwelling (first one occurring during January–March and the second one occurring during August–September) and downwelling (first one occurring during April–June and the second one occurring during October–December) Kelvin waves that propagate eastward along the equator and hit the Sumatra coast and bifurcate. The northern branches propagate counterclockwise over varied distances along the coastal wave guide of the BoB. The potential mechanisms that contribute to the mid-way termination of the first upwelling and the first downwelling Kelvin waves in the wave guide of the BoB are hypothesized. The second downwelling Kelvin wave alone reaches the southeastern AS, and it shows large interannual variability caused primarily by similar variability in the equatorial westerly winds during boreal fall. The westward propagating downwelling Rossby waves triggered by the second downwelling Kelvin wave off the eastern rim of the BoB also shows large interannual variability in the near surface thermal structure derived from SODA analysis. The strength of the equatorial westerlies driven by the east–west gradient of the heat sources in the troposphere appears to be a critical factor in determining the observed interannual variability of the second downwelling Kelvin wave in the wave guides of the equatorial Indian Ocean, the coastal BoB, and the southeastern AS.

Floating marine debris surface drift: Convergence and accumulation toward the South Pacific subtropical gyre

Whatever its origin is, a floating particle at the sea surface is advected by ocean currents. Surface currents could be derived from in situ observations or combined with satellite data. For a better resolution in time and space, we use satellite-derived sea-surface height and wind stress fields with a 1/3° grid from 1993 to 2001 to determine the surface circulation of the South Pacific Ocean. Surface currents are then used to compute the Lagrangian trajectories of floating debris. Results show an accumulation of the debris in the eastern-centre region of the South Pacific subtropical gyre ([120°W; 80°W]–[20°S; 40°S]), resulting from a three-step process: in the first two years, mostly forced by Ekman drift, the debris drift towards the tropical convergence zone (∼30°S). Then they are advected eastward mostly forced by geostrophic currents. They finally reach the eastern-centre region of the South Pacific subtropical gyre from where they could not escape.

Indian Ocean Rossby waves detected in HYCOM sea surface salinity

Rossby waves have been well identified in satellite derived sea surface height (SSH), sea surface temperature (SST) and ocean color observations. Studies of Rossby waves have yet to include sea surface salinity (SSS) as a parameter, largely because presently available in situ measurements of salinity lack sufficient spatial and temporal coverage, and as of now no methods are available for measuring salinity from a satellite. In this paper, we demonstrate that Rossby waves can be observed in SSS in the Indian Ocean by using simulations of the 1/12° global Hybrid Coordinate Ocean Model (HYCOM). HYCOM results compared favorably to SSS data provided by Argo floats and the World Ocean Atlas 2005 (WOA05) in selected grid boxes in the Indian Ocean. Hovmöller diagrams of HYCOM SSS anomalies and gradient show the distinct westward propagating signature of Rossby waves, with a steeper slope in longitude/time plots further from the equator. The propagation speeds, calculated from a 2D Radon Transform are comparable with new theoretical speeds for Rossby waves. Annual westward propagating signals in the SSS simulations at most of the latitudes in the Indian Ocean coincide with previous studies. We hope that future studies of Rossby waves in SSS using model results and eventually satellite measurements of salinity data will allow a better understanding of Rossby wave dynamics.

Eddy-Induced Heat Transport in the Subtropical North Pacific from Argo, TMI, and Altimetry Measurements

Abstract Basin-scale heat transport induced by mesoscale oceanic eddies is estimated by combining satellite-derived sea surface height and temperature [temperature data are from the TRMM Microwave Imager (TMI)] data with Argo float temperature–salinity data. In the North Pacific Ocean subtropical gyre, warm (cold) temperature anomalies of mesoscale eddies are found to be consistently located to the west of high (low) SSH anomalies. The phase misalignment between the temperature and velocity anomalies, however, is largely confined to the seasonal thermocline, causing most of the eddy-induced heat transport to be carried in the surface 200-m layer. By establishing a statistical relationship between the surface and depth-integrated values of the eddy heat transport, the basin-scale eddy heat transport is derived from the concurrent satellite SSH/SST data of the past six years. In the Kuroshio Extension region, the meandering zonal jet is found to generate oppositely signed eddy heat fluxes. As a result, the zonally integrated poleward heat transport associated with the Kuroshio Extension is at a level O(0.1 PW), smaller than the previous estimates based on turbulent closure schemes. Large poleward eddy heat transport is also found in the subtropical North Pacific along a southwest–northeast-tilting band between Taiwan and the Midway Islands. This band corresponds to the region of the subtropical front, and it is argued that the relevant temperature field for identifying this band in the turbulent closure scheme models should be that averaged over the seasonal thermocline.

Hydrographic Features off the Southeast Coast of Kyushu during the Kuroshio Small Meanders: A Case Study for Small Meanders that Occurred in 1994 and 1995 Spring

Three Kuroshio small meanders off the southeast coast of Kyushu that occurred during 1994 to 1995 were investigated by using satellite-derived sea surface temperature (SST) and sea surface height (SSH) maps, World Ocean Circulation Experiment (WOCE) Hydrographic Program (WHP) repeat section and Japan Meteorological Agency (JMA) hydrographic observations. Based on the satellite data, we observed that the three small meanders are formed by different processes: the triggering and growth of these meanders are caused by a cyclonic eddy propagating from the Kuroshio recirculation region or Kuroshio front meanders traveling from the East China Sea. Investigation of the two small meanders in 1994 and 1995 spring that are captured by the WHP observation showed quite consistent hydrographic features. On the nearshore side of the meandering Kuroshio, a countercurrent appears, associated with vertically uniform upward lifts of the isopycnals from sea surface to bottom at the boundary between the countercurrent and the Kuroshio. In the countercurrent region, the waters in the density ranges of the North Pacific subtropical mode water (NPSTMW) and the North Pacific Intermediate Water (NPIW) are more saline and less saline than typical waters that the Kuroshio carries in a non-small meander state, respectively. There are indications that high-salinity NPSTMW and low-salinity NPIW distributed off the Kuroshio was supplied to the countercurrent region. In the meandering Kuroshio flow, while there is no notable change in properties around the NPSTMW density range, salinity of the NPIW is significantly higher than that carried by the Kuroshio in a non-small meander state, but not higher than that in the Kuroshio at the Tokara Strait, which suggests that saline NPIW from the Tokara Strait, less mixed with low-salinity NPIW off the Kuroshio, may be carried by the meandering Kuroshio.

EDDY-INDUCED HEAT TRANSPORT IN THE SUBTROPICAL NORTH PACIFIC FROM ARGO, TMI AND ALTIMETRY MEASUREMENTS

Abstract Basin-scale heat transport induced by mesoscale oceanic eddies is estimated by combining satellite-derived sea surface height and temperature [temperature data are from the TRMM Microwave Imager (TMI)] data with Argo float temperature–salinity data. In the North Pacific Ocean subtropical gyre, warm (cold) temperature anomalies of mesoscale eddies are found to be consistently located to the west of high (low) SSH anomalies. The phase misalignment between the temperature and velocity anomalies, however, is largely confined to the seasonal thermocline, causing most of the eddy-induced heat transport to be carried in the surface 200-m layer. By establishing a statistical relationship between the surface and depth-integrated values of the eddy heat transport, the basin-scale eddy heat transport is derived from the concurrent satellite SSH/SST data of the past six years. In the Kuroshio Extension region, the meandering zonal jet is found to generate oppositely signed eddy heat fluxes. As a result, the ...

Satellite-derived sea surface height and sea surface wind data fusion for spilled oil tracking

An attempt is made to estimate the trajectory of the spilled oil from the sunken tanker Nakhodka occurred on January 2, 1997 in the Japan Sea by fusing two microwave sensor data, namely ERS-2 altimeter and ADEOS/NSCAT scatterometer data. In this study ‘fusion’ is defined as the method of more reliable prediction for the trajectory of spilled oil than before. Geostrophic current vectors are derived from ERS-2 altimeter and wind-induced drift vectors are derived from ADEOS/NSCAT scatterometer data These two different satellite-derived vectors are ‘fused’ together in the surface current model to estimate and evaluate the trajectory of spilled oil from the sunken tanker Nakhodka. The distribution of component of spill vector is mostly accounted for by the distribution of geostrophic velocity component during the study period with some discrepancies during March, 1997.

Seasonal and interannual variation in the hydrography of the Cariaco Basin: implications for basin ventilation

The hydrography of the Cariaco Basin (temperature, salinity, density, dissolved oxygen concentration) was studied using monthly observations collected between November 1995 and August 1998 at the CArbon Retention In A Colored Ocean (CARIACO) time-series station (10.5°N, 64.66°W). Satellite scatterometer wind estimates showed that changes in the wind preceded changes in hydrography by 1–2 weeks. Upward migration of isopleths within the upper 150 m was observed between November and May each year, when the Trade Wind was more intense. A seasonal deepening of the isopleths was observed when winds relaxed. A secondary upwelling event was observed every year between July and August, in response to an intensification of the southward component of the Trade Wind. Interannual variation in the upwelling cycle was driven in part by variations in wind intensity and in part by other events at time scales of 1–3 months. The latter were associated with 90–140 m deep intrusions of Caribbean Sea water that forced waters above them to the surface. Satellite-derived sea surface height anomaly maps demonstrated that these events were related to cyclonic and anticyclonic eddies moving along the continental shelf. Waters deeper than 1200 m showed small temperature and salinity increases of 0.0075°C yr −1 and 0.0016 yr −1, consistent with previous estimates.