Abstract

Upper ocean conditions in the Caribbean Sea are studied for long-term variability and trends using filtered surface observations and ocean model reanalysis fields. A principal component analysis is made, and trends in the leading mode are extracted. Sea surface temperature shows an accelerating upward trend while air pressure exhibits quasidecadal fluctuations. Sea surface height and subsurface temperature rise linearly while subsurface salinity exhibits fresher upper and saltier lower layers. The amplitude of warming is highest in the southern Caribbean east of 75°W near 150 m and lowest near the surface, indicating little role for a top-down process such as air-sea exchange. The freshening surface layer does not appear connected to river discharge or regional rainfall, so changes in ocean advection and sources are the likely drivers. Westward currents exhibit a reduction of throughflow and an influx from the Windward Passage. The Caribbean Current has slowed ~0.06 m/s in the reanalysis era. Crop yields show little sensitivity to ocean conditions but tend to follow rainfall. Marine catch per capita in the Caribbean follows subsurface currents and vertical motion but is less affected by temperature and salinity.

1. Introduction

The accumulation of greenhouse gases and consequent absorption of radiation has resulted in an atmospheric warming that is faster in the Caribbean than elsewhere in the tropics. This is related to a local acceleration of the Hadley cell and gas plumes that drift west from Africa [1], producing a drying trend that is projected to continue in the 21st century [2]. How the sub-surface ocean contributes to global warming has recently been explored using observed temperature profiles and ocean model projections [3]. Across the central Antilles, the rate of warming at the top of the atmospheric boundary layer is triple the ocean surface, so sensible heat fluxes are declining. This has implications for the thermodynamic energy available to tropical weather systems within the region [4] and to the advection of heat, moisture and momentum out of the region.

The Caribbean Sea is bounded on the south and west by South and Central America, and fringed on the east and north by the chain of Antilles Islands (8°N–25°N, 85°W–55°W) and Atlantic Ocean. Persistent subtropical trade winds, year-round sunshine, and consistent water exchanges result in little seasonal variation. The surface warm layer is >100 m deep and the upper 1200 m is stratified [5, 6]. Most Atlantic water infiltrates the Caribbean Sea through the Grenada, St. Vincent, and St. Lucia Passages in the southeast [7, 8] sourced from the meandering North Brazil Current that bears fresh water from the Orinoco River. From there, the Caribbean Current flows westward at ~0.5 m/s in the latitudes 13–16°N [911]. It turns northwest between Nicaragua and Jamaica, with a branch forming the counterclockwise Panama Gyre [8, 12, 13]. Westward currents exit via the Yucatan Channel eventually to be drawn into the Gulf Stream, which receives contributions from the Antilles Current and its Atlantic sources [14, 15]. Schmitz and Richardson [16] and Johns et al. [7] indicate that the Caribbean Sea is well ventilated with ~28 Sv of throughflow from both the North and the South Atlantic [17, 18]. South Atlantic waters are fresher and more oxygenated than North Atlantic waters of the same density, they enter the Caribbean near Trinidad. While the mean structure of water masses and currents in the Caribbean is well known and processes underlying year-to-year fluctuations are being uncovered [19], long-term trends and influences are relatively unexplored.

Here, long-term variability and trends in the upper ocean across the Caribbean are described in terms of spatial pattern and temporal amplitude. Some of the processes underlying the trends and their biophysical consequences are studied. The key question is how is the global warming signal reflected in the Caribbean Sea? The hypothesis is that the response is top-down (greatest trend near the surface) and spatially homogeneous.

2. Data and Methods

Long-term trends at the ocean surface are characterized by ship and satellite data reanalyzed by the National Ocean and Atmosphere Administration (NOAA) for sea surface temperatures (SST; see [20]) and by the National Center for Environmental Prediction (NCEP) for wind [21, 22]. Sea level pressure (SLP) derives from ship data reanalysis by the Hadley Center [23], and tropical cyclone prevalence is from Emanuel [24]. Data were gathered from the Climate Explorer http://climexp.knmi.nl/ and Climate Library http://iridl.ldeo.columbia.edu/ websites over the period since 1854. Surface data are averaged over the Caribbean (8°N–25°N, 85°W–55°W; ~6 106 km2), where observations are relatively dense, except north of Panama (Figure 1(a)). This domain cuts off the Gulf of Mexico and Straits of Florida. Over time, ship reports of sub-surface temperature vary from near zero before 1880 to >20/1° cell/year after 1925 (Figure 1(b)). Recent years have seen a sharp upturn with profiling floats.

Upper ocean conditions are described, since 1958, by Simple Ocean Data Assimilation (SODA) fields version 2.4 [25] comprising sea surface height (SSH), temperature ( ), salinity ( ), currents ( ), and vertical motion ( ) at 50 km horizontal and ~50 m vertical resolution. Surface wind stress ( , ) is derived from European Community Medium-range Weather Forecasts (ECMWF; see [26]). The ocean fields are based on a numerical model assimilation of in situ and remotely sensed data, and include coastal station and ship hydrographic data; satellite surface temperature, altimetric height, and winds; ocean drifters and profiling floats, ingested by the Global Ocean Data Assimilation System and its predecessors. Although the Caribbean has reasonable coverage for surface elements and ocean profiles in the reanalysis era (Figure 1(b)), it is useful to minimize errors contributed by uneven observations through spatial and temporal aggregation as outlined below.

To define the pattern and trend in the upper ocean, eigenvectors within a covariance matrix are calculated via the IRI Climate Library for each variable: SSH, , , , , , , , from annual averages per grid cell and depth in the period 1958–2007. The leading principal component (mode-1) representing the dominant variability is analyzed as maps (8°N–25°N, 85°W–55°W) using depth averages: 1–100 m for and , and 1–200 m for , , . This method is repeated for E-W depth sections (85°W–55°W, 1–350 m, averaged over 8°N–25°N) and N-S depth sections (8°N–25°N, 1–350 m, averaged over 85°W–55°W). In this manner, numerous profile observations are collapsed into each analysis.

Principal component (PC) analysis is employed as a way to cluster and deduce the long-term variability of concern here. Table 1 intercompares the variance explained by the first and second modes. The leading mode containing the trend, accounts for about half the variance in all fields except the ocean currents—where observations are potentially less reliable. Trends were found to be insignificant in secondary modes. To isolate the trend signal and establish the rate of change, the slope of the mode-1 PC time score is mapped per grid cell. To determine its relative importance, the fit of the trend regression is calculated after standardizing and smoothing with a 5-yr running mean. Statistical significance is evaluated via the Pearson product-moment test with 10 degrees of freedom for sub-surface data ( at 95% confidence) and 25 degrees of freedom for the longer surface records . Comparison of single-station long-term records is limited to SSH from Puerto Rico (1966+) and the Orinoco River discharge at 8°N, 63°W (1923+) and its catchment run-off from hydrological reanalysis [27].

Searching for long-term trends has certain pitfalls. Trends in upper ocean heat content are believed to have accelerated since the 1970s [3], and the sub-surface record may reflect this signal. But the observational density changes with satellite data assimilated after 1979. Furthermore, the time series contain multidecadal oscillations that spread toward the Caribbean from higher latitudes (cf. SLP mode-1 pattern Figure 1(c); [2830]). Using the wavelet co-variance technique of Torrence and Compo [31], it is seen that Caribbean SLP and SST share spectral energy around 16 years in the reanalysis era (Figure 1(d)). Thus, it is useful to distinguish trend from variability in the analysis. Ultimately, the results here should be considered “best estimates” with uncertainty contributed by instrument sampling errors and the interpolation of uneven observations to model grid-cells.

Biophysical relationships are analyzed by comparison of upper ocean PC scores with annual marine catch (including all species, ~75% fish) for the Caribbean Sea, excluding Central and North American waters, obtained from the website http://www.fao.org/fishery/ in the reanalysis era. There are known reporting problems associated with the commercial exploitation of marine resources [32]; few countries have quality controlled historical catch data. For terrestrial resources, annual crop yields per (planted) area as an aggregate for all countries in the Caribbean (excluding Central and North America) were obtained from http://faostat.fao.org/ in the reanalysis era. Agricultural outputs display a step in 1985, so departures were computed from separate means. The sensitivity of crop yield (for avocado, coconut, citrus fruit, mango, spices, and sugar cane) was tested by comparison with upper ocean PC scores and with Caribbean area-averaged rainfall from the Global Precipitation Climatology Project [33]. Variability in resources may arise from environmental effects or from human influences such as changes in effort or reporting. With this in mind, the marine catch was divided by the Caribbean population. The marine catch per capita and crop yield time series are expected to integrate and lag climate, so cross-correlations were computed on annual resource indices at a 1-yr lag.

3. Results

3.1. Surface Trends

The historical Caribbean-area surface time series are given in Figure 2. SST exhibit quasi-decadal oscillations and a cool period in the early 20th century. There is a positive second-order trend with 40% fit suggesting acceleration of warming in recent years. Tropical cyclones in the Caribbean oscillate with SST but exhibit little trend. There is a widening gap (Figure 2(a)) induced by declining sensible heat fluxes ( > SST). SLP is inversely related to SST ( ), with shared spectral energy in the bands ~30 and ~60 years before 1960 (cf. Figure 1(d)). SLP exhibits a weak downtrend and a 16-yr cycle since 1960 (Figure 2(b)). For surface wind, the reanalysis starts after 1900. Zonal wind has little trend over the 20th century while meridional wind exhibits a positive trend of 17% fit (Figures 2(c) and 2(d)).

The mode-1 loading pattern for reanalysis sea surface height (SSH; Figure 3(a)) is positive everywhere with higher values off Venezuela and around Cuba, and lower values along the central Antilles and Panama gyre. The faster rise in the vicinity of Key West/Havana relates to subsiding land. The SSH time score has a linear upward trend with 72% fit (Figure 3(c)). Station data for Puerto Rico indicate a slower rise of SSH than the regional mean. The wind stress “change vector” (from , trends) is eastward in the central Caribbean. South of Cuba and east of Trinidad NE trade winds have strengthened. The trend for is second order with 68% fit (Figure 3(d)).

3.2. Subsurface Maps

Maps of trends for , , currents, and are presented in Figure 4. There has been widespread warming of the Caribbean thermocline layer as expected. The upward slope of is greatest off Venezuela (Figure 4(a)) associated with diminished trade wind upwelling in recent decades. Salinity in the thermocline has declined particularly in the Caribbean Current 13°–17°N, east of 75°W (Figure 4(b)). The fresh tongue suggests an influx of South Atlantic waters. Currents tend to loop around Jamaica from the Windward Passage and connect with eastward trends in the Caribbean Current (Figure 4(c)). The vertical motion trend is upward around Hispanola/Puerto Rico and sinking in the south off Colombia and Venezuela (Figure 4(d)).

Roemmich [12] found that surface and thermocline waters leaving the Caribbean through the Yucatan are composed of 17 Sv from the SE Caribbean and 6 Sv from the Windward Passage. What is seen in the trends here is a reduction of inflow from the SE Caribbean and an increase from the Windward Passage. There is a trend of reduced outflow toward the Yucatan. and linear trends achieve 65% and 78% fit, respectively (Figure 4(e)). Because salinity in the SE Caribbean is affected by the Orinoco plume [34], its discharge is considered (Figure 4(f)). The trend of river input over the historical record is weak, and regional rainfall is declining, so changes in ocean current advection underlie the freshening, not local air-sea fluxes.

3.3. Depth Sections

The temperature trend in N-S depth section (Figure 5(a)) exhibits a maximum from 11°–16°N at 100–200 m depths. The amplitude of warming is weakest at the surface and in the south below 300 m. The salinity depth section (Figure 5(b)) exhibits negative (positive) loading above (below) ~150 m, hence the halocline has weakened. This interface slopes downward to the north with greatest freshening from 12°–17°N above 50 m. The N-S depth section for zonal currents is given in Figure 5(c). Alternating bands are evident; the maximum eastward trend is near 15°N in the Caribbean Current. The rate of change for currents expressed as vectors in depth sections is southeastward. For the N-S section of current (Figure 5(d)) trends are greatest from 18–23°N above 300 m reflecting inflow from the Antilles Current. Vertical components are largely upward (downward) north (south) of 15°N, reflecting a trend toward rising (sinking) motions in the central Antilles (off South America). Weaker upwelling in the southern Caribbean relates to a weakening of trade winds since 1990 (cf. Figure 3(b)). The current trend is generally eastward over the Nicaragua Rise and in the central Caribbean 75°–60°W above 200 m (Figure 5(e)). Trends and slow fluctuations of the Caribbean thermocline circulation appear modulated by inflow from Grenada ( ) and Windward ( ) Passages. The time series for currents has a positive second-order trend with 94% fit (Figure 5(f)). Multidecadal oscillations appear to influence the circulation (second-order trends) more than thermodynamic properties (linear) in the reanalysis era.

3.4. Resources

The marine catch per capita is relatively stationary in the 1960s (Figure 6(a)) and rises in the 1980s as tropical cyclones decline (cf. Figure 2(a)). Thereafter, catch becomes depleted, with a negative second-order trend of 59% fit. Population growth follows global trends, doubling in the reanalysis era and stressing resources. Caribbean crop yield followed rainfall until recently; a weak drying trend is noted (Figure 6(b)). Crop yield is less sensitive to upper ocean conditions; cross-correlations at 1-yr lag are generally negative with declining yields and moderate with respect to ocean currents (cf. Table 2).

Table 2 summarizes the PC score mode-1 cross-correlations with linear trend. Upward, significant values are seen for all variables except wind stress. Thus, global warming signals (and reduced dS/dz) emerge from the analysis. Cross-correlations between environmental PC scores and marine fish catch lagged by 1-year are also summarized in Table 2. Surprisingly, the values are insignificant, suggesting that patterns of advection and thermodynamic properties may be less important than factors affecting catch effort and efficiency (e.g., storm frequency, sea state). The moderate correlation with suggests that increased catch follows northward winds originating from the Venezuela upwelling zone.

4. Summary

Conditions in the Caribbean Sea have been studied for long-term trends in the period of historical surface observations using area-averages and in sub-surface fields using principal component analysis and trend extraction. Continuous data coverage started in the 1850s for SLP and SST, the 1900s for winds, and in the late 1950s for model reanalyzed sea surface height and sub-surface , , currents and vertical motion. Data aggregation methodologies were employed to reduce variability in time, space, and depth, leaving the long-term trend which was analyzed by regression to first or second order. SST reflected an accelerating upward trend while SLP exhibited quasi-decadal fluctuations. SSH and sub-surface rose linearly at near-global rates while sub-surface exhibited fresher upper and saltier lower layers. The amplitude of warming was spatially homogeneous but weak in the upper 100 m, so the top-down hypothesis related to air-sea exchange was unsupported. Averaged over the Caribbean, temperatures have risen at twice the rate (.015 C/yr) at 150 m than in other layers. Sea surface height trends were found to be highest (+.0014 m/yr) at 14° and 23°N and lower (+.0008 m/yr) at 10° and 18°N. Freshening of the surface layer was unrelated to regional rainfall or river discharge, suggesting changes in ocean advection and sources. Salinity has declined fastest (−.004 ppt/yr) near 30 m. Crop yields were insensitive to ocean conditions, and followed rainfall as expected. Marine catch per capita in the Caribbean was unrelated to and , and followed enhanced northward winds that spread productivity.

The Caribbean Sea has experienced a rate of climate change that is faster than most other coral reef zones of the world according to a recent study [35]. This is a cause for concern—given declining food resources, receding beaches, and the dependence of local economies on coastal tourism.