Table of Contents
International Journal of Oceanography
Volume 2014 (2014), Article ID 325321, 16 pages
http://dx.doi.org/10.1155/2014/325321
Research Article

Environmental Forcing of Red Tides in the Southern Benguela

1University of Zululand, KwaDlangezwa 3886, South Africa
2Physics Department, University of Puerto Rico, Mayagüez, PR 00681, USA

Received 18 February 2014; Revised 10 April 2014; Accepted 24 April 2014; Published 16 June 2014

Academic Editor: Robert Frouin

Copyright © 2014 Mark R. Jury. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

The Southern Benguela cape upwelling plumes have inshore wind shadows prone to red tides in late summer. Their intensity and coverage are estimated by satellite fluorescence measurements in the period 1997–2012 and qualified by in situ reports. High satellite fluorescence cases are identified at daily to seasonal time scales, and characteristics of the upper ocean and lower atmosphere are studied using third generation reanalyses. A dominant feature is easterly winds over the Cape Peninsula (34°S, 18°E) induced by a ridging anticyclone-coastal low weather pattern. Over Cape Columbine (33°S), there is a wind shadow with cyclonic wind and current shear. Composite atmospheric profiles reveal a 4°C temperature inversion near 500 m that traps a coastal wind jet >6 m/s below 200 m. The composite shelf oceanography shows a relic upwelling plume below 10 m overtopped by warmer water near the coast, providing the thermal stratification needed for biotic aggregation. Data from the IPSL5 coupled climate model over the period 1980–2080 indicates that environmental conditions favoring red tides may become more frequent.

1. Introduction

The Southern Benguela (31.5–34°S, 17.5–18.5°E) is a zone of summer upwelling and equatorward flow [14]. The upwelling intensifies next to two capes at 34°S and 33°S [47]. In the north of each cape, there is wind shadow and corresponding band of warm surface water [8, 9]. Cyclonic wind shear produces a clockwise circulation in the leeward bay [912] where high chlorophyll content promotes dinoflagellate blooms or “red tides” [13]. In situ surveys have found that red tides form in ~5 × 20 km patches and accumulate over a few days [1416], mainly in February to April season.

The upwelling is pulsed by passing weather systems and amplified by coastal low pressure cells [1719]. As they transit the Southern Benguela, equatorward flow becomes shallow and sheared by the topography. Wind shadow zones grow offshore and stratified conditions lasting many days contribute to red tides [2024], which appear rapidly at high concentrations suggesting that local growth is boosted by confluent circulations [16, 25, 26]. The decay of red tides leads to low oxygen water and marine life mortalities with economic consequences [27].

Most red tides exhibit spectral radiance in the fluorescence band ±683 nm. These can be measured near the surface by satellite remote sensing [28, 29]. The sea wide-field sensor (SeaWiFS) collected data from 1997 to 2010 in spectral bands from 412 to 865 nm. Since 2002, the moderate resolution imaging spectroradiometer (MODIS) provides declouded 8-day 9 km resolution imagery in 405–877 nm spectral bands that enable a fluorescence line height (FLH) calculation [30]. MODIS FLH is found to be consistent with in situ observations of red tide [3133] as a radiometric index of relative harmful phytoplankton biomass. The European MERIS satellite also has spectral bands across the fluorescence peak capable of detecting red tides and cyanobacteria blooms [34], even in coastal turbid waters [35]. FLH is sensitive to phytoplankton nutrient status, growth rates, species, and size through variations in quantum yield. In addition to fluorescence proxies, satellite measured chlorophyll (CHL) indicates high biomass dinoflagellate blooms and sea surface temperature (SST) helps in establishing the background conditions that contribute to the formation of red tides during wind-driven coastal upwelling relaxation. Satellite altimetry can resolve the background circulation in conjunction with ocean model assimilation of in situ observations.

Multiple environmental factors can be evaluated by regression of variables onto an FLH index. Additional FLH proxies may be generated from the National Aeronautics and Space Administration (NASA) ocean biogeochemical model (NOBM). The NOBM physical [36] and chemical submodels [37] assimilate satellite data in parallel with environmental observations. The model calculates budgets and interactions between phytoplankton, nutrients, minerals and zooplankton, and consequent dissolved organic material and detritus. The NOBM employs a sophisticated Kalman ensemble filter. Its use in the coastal zone is limited by a 100 km resolution, so application here is made in the broader upwelling zone north of Cape Columbine.

Ocean-atmosphere reanalysis products are employed here to study the environmental forcing of red tides in the Southern Benguela in February–April season in 1997–2012. This is a challenge given the coastal pulsed nature of upwelling [3841]. The paper addresses the following research questions. What are the temporal and spatial characteristics of satellite fluorescence that indicate red tides? What are the environmental influences on red tides at event to seasonal time scales? Are conditions favoring red tides more likely in future?

Methods are employed to bring out the environmental forcing of red tides at time scales from daily local weather to interannual ocean influence. The spatial scale includes the two shelf zones leeward of the Cape Peninsula and Cape Columbine upwelling plumes, where red tides are most frequently observed. The study does not attempt to “predict” red tides but to understand the multivariate environmental forcing. Model fitting is done only to extend the MODIS record back to SeaWiFS and to establish a daily index corresponding with red tides.

2. Data and Methods

The primary index to estimate red tide intensity and coverage in the Southern Benguela is the MODIS fluorescence line height (FLH) 9 km 8-day data from July 2002 to December 2012 obtained from NASA. A SeaWiFS FLH proxy is calculated by multivariate regression of radiance data in the period of overlap with MODIS, extending the record back to September 1997. Time series are averaged over the Southern Benguela shelf zone: 34–31.5°S, 17.5–18.5°E (Figure 1(e)), to indicate red tides within the two cape upwelling shadows. In situ red tide reports from the Fisheries Service (MCM) were used to qualify the satellite estimates, and cases were screened accordingly. In addition to FLH data, satellite chlorophyll (CHL) is analyzed given its relationship to winds and currents in coastal upwelling zones [42, 43]. A proxy for FLH at daily time scale is an arithmetic average of NOBM chlorophytes and cyanobacteria. For comparison with oceanographic sets, the SeaWiFS-MODIS FLH is binned into months. Thus, we have FLH or FLH proxies at three time scales (daily, 8-day, monthly) to investigate the shelf-scale environmental forcing of red tides in the period 1997–2012. FLH case studies (qualified by in situ reports) are listed in Table 1.

tab1
Table 1: FLH cases for inclusion in composites screened by in situ reports.
fig1
Figure 1: Mean annual cycle averaged over the Southern Benguela 1997–2012 for the following: (a) 8-day SeaWiFS-MODIS chlorophyll and fluorescence; (b) monthly MODIS SST and SODA 1–100 m temperature; (c) ECMWF meridional wind stress and SODA 1–100 m salinity; (d) coupled forecast system (CFS) air-sea latent heat flux and lowest 2.5% absolute vorticity ×10−5 s−1; (e) day-time land surface temperatures for high (daily) FLH cases with coastal topography and place names. Box is the averaging area for all analyses: 34–31.5°S, 17.5–18.5°E; dot is regional airport.

Global atmospheric and oceanography datasets with 30–50 km resolution are used to understand the shelf dynamics. Monthly wind stress, surface and upper wind, wind vorticity, sea level air pressure (SLP), surface temperature (Ts), latent heat flux (evaporation), specific humidity, and precipitable water are obtained from the coupled forecast system (hereafter, CFS [44]). Daily atmospheric data are obtained from the National Center for Environmental Prediction (NCEP) operational assimilation system and time- and space-averaged to match FLH for temporal analysis. Oceanographic fields are provided by Simple Ocean Data Assimilation version 2.1.6 (SODA, [45]) and include subsurface (1–100 m) temperature, salinity, currents and vertical motion, sea surface height, and European Community Medium-Range Weather Forecasts (ECMWF) wind stress. Monthly ocean data assimilation is performed within the Geophysical Fluid Dynamics Lab Modular Ocean Model version 4 [46]. Coastal upwelling and mesoscale winds are studied using daily 25 km NOAA SST data [47], 25 km J-OFURO reanalysis [48], satellite scatterometer measurements from QuikSCAT and ASCAT, and 4 km Meteosat infrared imagery. Atmospheric mixed layer profiles are studied by composite averaging of daily aircraft (AMDAR), radiosonde, and NCEP wind and temperature profiles near Cape Town Airport for high (daily) FLH cases. In general, data from “first generation” reanalyses (at 200 km resolution) are unsuitable for shelf-scale analysis and use here.

The environmental forcing of red tides is evaluated by multivariate regression of candidate predictors onto FLH in the form of (a) mean annual cycle, (b) monthly, (c) 8-day, and (d) daily data. This is done by backward stepwise removal of less influential variables, until an optimal multivariate algorithm is found. These time scales represent a progression from the “climate envelope” to the “local weather” which induces red tides in the coastal zone. Oceanographic and atmospheric structures are analyzed as depth or height sections over the Southern Benguela shelf (34–31.5°S, 17°–19°E), based on composites of the seven highest FLH months in February to April season minus the seven lowest FLH months (hereafter “high minus low”). Red tide events in the 8-day FLH record are studied by composite-averaging maps of SLP, SST, winds, latent heat flux, and MODIS FLH imagery. FLH index-to-field correlations are calculated for monthly surface temperature, wind stress, latent heat flux, meridional upper wind, and sea level air pressure both within and around the Southern Benguela in February to April season. Pairwise correlations are analyzed between FLH and monthly environmental variables averaged over 34–31.5°S, 17.5–18.5°E, except the s-n SLP index: 10–20°E, 35–45°S minus 15–25°S. Statistical significance is evaluated according to the degrees of freedom as limited by record length, temporal resolution/averaging, and persistence.

Given the peculiar sea level air pressure (SLP) pattern found in the FLH index-to-field correlations, data from the French Institut Pierre Simon Laplace (IPSL) version 5 coupled general circulation model used for climate change projections [49] is evaluated over the period 1980–2080. Principal component clusters are calculated for February to April SLP and air temperature (). The leading mode of the IPSL5 simulated SLP resembles the FLH correlation pattern, so its trend and variance are analyzed. The IPSL5 is selected from a long list of models in the Intergovernmental Panel on Climate Change (IPCC) Intercomparison Project (CMIP5) because it exhibits low differences between simulated and observed conditions in the Benguela upwelling zone [50]. The IPSL5 model is described by Dufresne et al. [51] and has a dynamic coupling and feedback scheme that accurately simulates past climate. Its horizontal resolution is ~100 km and representative of large scale environmental forcing, but not shelf—or event—scale processes.

3. Results

3.1. Annual Cycle

The mean annual cycle is studied based on red tide index area-averaged 8-day FLH and CHL values (Figure 1(a)). Both variables show a gradual rising trend from December to April, following an initial surge in November. CHL drops precipitously in April while FLH remains high then. Both annual cycles bottom out in July as coastal upwelling declines. It can be noted that river inflows are near zero from January to April, and thus “mistaken” turbidity events are unlikely. The annual cycle of SST is at a maximum in January while subsurface temperatures are at a minimum in March (Figure 1(b)). The larger difference represents stratification necessary for biotic aggregation. Subsurface salinity is inversely related to SST with a minimum in January and a maximum in July. Meridional wind stress has a large annual cycle with equatorward values peaking from December to February, then falling rapidly in March to negative values in June (Figure 1(c)). Latent heat flux is at a minimum over the shelf from February to April, and its lowest 2.5% of cases falls almost to zero in February. Similarly wind vorticity is most cyclonic from February to April, when standing clockwise rotors prevail over the coast. Wind vorticity rises rapidly (anticyclonic) after May (Figure 1(d)) to inhibit upwelling and red tides. Figure 1(e) provides a map view of the Southern Benguela and its topography. Day-time land surface temperatures in high (daily) FLH cases exhibit values >40°C in wind shadow zones. Pairwise correlations of environmental variables with the annual cycle of FLH (Table 2) are high for latent heat flux (−0.95), subsurface temperature (−0.92), and mixed layer depth (−0.88). Multivariate regression onto the annual cycle of monthly FLH fits 97% of variance with negative coefficients for subsurface temperature and salinity and meridional wind stress (Table 3) with ~3 degrees of freedom.

tab2
Table 2: Pairwise cross-correlations of monthly data.
tab3
Table 3: Multiple regression statistics: mean annual cycle of monthly FLH.
3.2. Characteristics of High Fluorescence Spells

The environmental forcing of red tides in early 2002 is studied using daily FLH proxy and daily NCEP data. Multivariate regression (Table 4, Figure 2(a)) yields negative coefficients for wind (−1.06), wind (−0.44), SST (−0.49), south-north SLP (−0.67), and mixed layer depth (−0.79) and a positive coefficient for atmospheric precipitable water (+0.28). The multivariate environmental algorithm explains 52% of FLH variance (Table 4) with >50 degrees of freedom. Analyzing the normalized time series, frequent upwelling relaxations are noted (Figure 2(b)): SSTs peaked three times, meridional winds reversed, and the mixed layer was shallow. Zonal winds were persistently westward in February coincident with a positive s-n SLP index (Figure 2(c)) indicating a “ridged anticyclone-west coast trough” weather pattern. Precipitable water was high at the end of February 2002, due to an influx of tropical air. Meridional (equatorward) winds remained shallow during this period (Figure 2(d)) and thus sheared by the coastal topography. Wind stress anomalies along the coast were offshore and poleward, and satellite measured fluorescence was high in the shelf zone leeward of the two upwelling plumes (Figure 2(e)). Pitcher et al. [9] summarize in situ reports of repeated red tides, low oxygen water, and ecological mortalities in St. Helena Bay in early 2002.

tab4
Table 4: Multiple regression statistics: daily FLH index, Nov. 01–Jul. 02.
fig2
Figure 2: (a) Scatterplot of environmental model versus daily NOBM FLH proxy (cf. Table 4). (b) January to May 2002 time series of daily normalized departures of CFS meridional wind, NOAA SST, and NOBM mixed layer depth. (c) CFS zonal wind, precipitable water, and south minus north sea level air pressure index. (d) CFS meridional wind cross-section averaged 34–31.5°S, Feb–Apr 2002, indicating depth of equatorward flow in km. (e) SeaWiFS FLH proxy (10−2), 6–13 March 2002, and ECMWF wind stress anomalies (March 2002, max = 0.08 Nm−2). Arrows in (b) highlight upwelling relaxation.

Spatial features of high (daily) FLH events are studied. Wind and surface temperature maps on 18 March 2007 and MODIS FLH in the period from 14 to 21 March are given in Figures 3(a) and 3(b). Southeasterly winds prevailed off the Cape Peninsula and there was a wind shadow in St. Helena Bay. Cyclonic wind vorticity west of the two capes was ~−2 10−5 s−1. Meteosat infrared temperatures exhibited a wide zone of coastal upwelling <12°C and land temperatures above 22°C. The MODIS 8-day FLH map of 14–21 March 2007 (Figure 3(b)) suggested red tides along the coast of St. Helena Bay and seaward of the upwelling plumes at 31.3°S and 33.1°S. Scatterometer wind maps for two different high FLH cases are given in Figures 3(c) and 3(d). They exhibit a wind shadow with wind speeds <7 m/s inshore and boundary layer height ~200 m, contrasting with >15 m/s winds in the south and boundary layer height ~1 km. On 24 March 2012, there was a long shear line NW of the Cape Peninsula, initiated by friction of the Hottentots Holland Mountains and reinforced by negative sensible heat fluxes (−20 W m−2) over the upwelling zone. In St. Helena Bay, the winds curled cyclonically, and poleward wind flow was observed along the Namaqua Coast. The surface temperature and wind patterns in two further high FLH cases are analyzed in Figures 4(a)4(d). All share similar features: the upwelling zone was wider and further south than usual. SSTs were <13°C off the cape at 34°S and ~16°C in St. Helena Bay. Land surface temperatures (at 20:00 LST) were well above 22°C and the satellite imagery was cloud-free due to offshore wind flow. Scatterometer winds were >12 m/s in the south, and there was a broad wind shadow <4 m/s to the north where stress vectors turned clockwise.

fig3
Figure 3: (a) Meteosat 4 km infrared image and CFS wind vorticity (contours, 10−5 s−1) for high FLH case of 18 March 2007; (b) MODIS FLH (10−2) for 14–21 March 2007 and wind streamlines. Scatterometer wind maps (color flags, m/s) for (c) 20 March 2007 and (d) 24 March 2012 with CFS atmospheric mixed layer height (m) labeled at key points.
fig4
Figure 4: Meteosat 4 km infrared images for two high FLH cases: (a) 7 March 2002 and (b) 24 March 2012, all at 20:00 LST. (c, d) Scatterometer 25 km gridded wind speed (shaded, m/s) and CFS wind stress vectors for the same cases.

The 8-day FLH time series (Figure 5(a)) that is the basis for most analyses is filled with high and low spells embedded in a background seasonal rhythm. Multiple regression of surface weather variables yields an algorithm fitting 28% of variance with negative influences from PAR, zonal wind, and SST and with positive influences from meridional wind and humidity (Table 5). A higher fit could be possible with subsurface ocean variables, but these are not available at multiday time scales. Fitting continuous data provides insights into the seasonal forcing of red tide. However, peak events tend to occur in February–April, wherein seven cases (Table 1, cf. [52]) are used to calculate daily and 8-day composite maps of environmental conditions. Within each high FLH spell, there is a ridging anticyclone-coastal low event that stands out. The composite SLP map sequence (Figure 5(b)) reveals an anticyclone moving eastward along 40°S, reaching 20°E on day-0 and 33°E on day+1, while a coastal low event shifts from 28°S to 32°S. An important feature for red tides is the repetition of this weather sequence and development of a west coast trough to the north of an anticyclone belt. The composite SST anomaly map for high FLH cases (Figure 5(c)) has a cold anomaly between the two capes that represents a broadening of the upwelling zone. There is a warm anomaly offshore, so the thermal gradient is strengthened. The MODIS satellite FLH composite in February-March 2005 (Figure 5(d)) exhibits high values along the Namaqua and St. Helena Bay coasts and in a meridional axis west of Cape Columbine. Low FLH values are noted in the Cape Peninsula upwelling plume. Wind stress anomalies during February-March 2005 were poleward in the vicinity of Cape Columbine.

tab5
Table 5: Multiple regression statistics: 8-day FLH index, Sep. 97–Dec. 12.
fig5
Figure 5: (a) Time series of 8-day Sea-WiFS/MODIS fluorescence (area as in Figure 2(e)) with circles to identify the 7 high FLH cases confirmed by in situ reports. (b) Composite high day-0 (left) and day+1 sea level air pressure maps. (c) Composite high day-0 SST anomaly. (d) Example of MODIS FLH (10−2) and ECMWF wind stress anomalies (max = 0.08 Nm−2) in February-March 2005.
3.3. Structure of High Fluorescence Spells

The coastal environment accompanying red tides in the Southern Benguela is studied. Past research has identified that ocean stratification is important, but what about the atmosphere? In Figure 6, composite atmospheric profiles from NCEP model assimilation, AMDAR aircraft, and radiosonde near Cape Town are analyzed for high daily FLH cases. Winds directions tend to be southerly to varying degrees and shift to easterly with height, as expected from the Ekman spiral effect of surface friction. Model and radiosonde wind direction profiles are smoother than aircraft. Wind speeds exhibit a near-surface coastal jet of varying thickness, with strongest winds from 100 to 200 m in aircraft and radiosonde profiles, and somewhat deeper in the model representation. Dew point temperatures are 10–15°C and decrease with height consistent with dry easterly downslope flow, locally known as “berg winds.” Temperature inversions are found in all profiles, indicative of stable atmospheric conditions that inhibit turbulent momentum transfer at the top of the mixed layer. The increase is ~4°C from 10 m to 500 m in aircraft and radiosonde profiles and weaker in the model analysis. With surface air temperatures near 18°C and coastal SST below 15°C, an inversion exists over the upwelling zone which helps to maintain the wind shadow.

fig6
Figure 6: Composite (a) wind direction, (b) wind speed (m/s), (c) dew point, and (d) temperature (C) profiles with height (m) for high daily FLH cases from NCEP model assimilation, AMDAR aircraft, and radiosonde near Cape Town.

Composite ocean structure for seven high minus low FLH months (cf. Table 1) is illustrated in Figures 7(a)7(c) along coast axis. A key feature is the negative ECMWF meridional wind stress difference indicating a wind shadow over Cape Columbine (33.2–32.7°S). Zonal wind stress is neutral in the south but westward (-) over the Namaqua coast (32.3–31.5°S), indicating berg winds. Sea surface height shows a positive difference in the wind shadow and a negative slope to the north. Composite SODA currents are onshore- and poleward-directed to the south, suggesting a “pulling” action on shelf waters. Salinity differences are positive in the south with a slope toward neutral conditions elsewhere. SODA subsurface temperature differences are neutral in the south and exhibit a gradient toward cold conditions north of Cape Columbine: a broadening of the upwelling zone there. Vertical motion differences are the inverse: negative (downwelling) in the south and neutral elsewhere.

fig7
Figure 7: Composite 7 high minus 7 low FLH months averaged 17.5–18.5°E from SODA for the following: (a) wind stress and sea height, (b) currents and salinity, and (c) temperature and vertical motion differences, to study alongshore gradients.

Composite ocean structure for high minus low FLH months is analyzed as cross-shelf sections in Figures 8(a)8(d). These help to identify the “climate envelope,” within which the red tides events are formed. The temperature-salinity sections identify subsurface cold (−0.7°C) fresh differences on the shelf edge (18.0°E) denoting a relic upwelling plume, while warm differences are located offshore. Most significantly, the SODA reanalysis detects the warm salty surface layer at the coast (18.6°E). The current differences reveal several features consistent with in situ surveys during red tide [53]. A zonal downwelling circulation is found, with poleward flow differences (−0.03 m/s) in the upper layer from 17.8 to 18.4°E. The cross-shelf structure of wind stress differences is most negative at 18.2°E. The V-shape indicates shear lines that separate active upwelling and quiescent zones. The sea surface height is elevated on the coast due to onshore transport and supports surface poleward currents. It is suggested that the ability of 50 km resolution ocean reanalysis products to represent many known features underlying red tides in the Southern Benguela means that operational forecasts of conditions suitable for red tide are within reach.

fig8
Figure 8: Composite 7 high minus 7 low FLH months averaged 34–31.5°S from SODA for the following: (a) temperature, (b) salinity, (c) currents: meridional differences shaded, zonal overturning vectors max = 0.02 m/s, and (d) wind stress and sea height differences, to study cross-shelf gradients.
3.4. Environmental Influences on Seasonal Fluorescence

FLH index-to-field correlations are analyzed for the February–April season in Figures 9(a)9(e). This analysis promotes an understanding of the interannual forcing of red tides. The CFS correlation maps exhibit a pattern of increased upwelling and equatorward wind flow in respect of high FLH seasons. Surface temperature correlations are negative in the upwelling zone () and positive inland. Wind stress correlations are positive and upwelling favorable, but negative correlations with latent heat flux suggest light winds between the two capes. Further inland there are positive correlations referring to increased moisture flux. Taking a wider perspective (Figure 9(d)), the FLH index-to-field correlations are positive (equatorward) for meridional upper winds and refer to an anticyclonic ridge in the midlatitude jet stream. The sea level air pressure map is most important (Figure 9(e)). It exhibits south-high/north-low values: a ridged anticyclone-west coast trough weather pattern. Thus, the climate “envelope” contains event-like features that could enable forecasts at long lead times.

fig9
Figure 9: Correlation maps using February–April FLH index with CFS reanalysis: (a) surface temperature , (b) meridional wind stress tau , (c) latent heat flux , (d) 200 mb meridional wind V200, and (e) sea level air pressure SLP. Values of significant at 90% confidence in the period 1998–2010.

Coupled model simulations that form part of the IPCC CMIP5 assessment are available for analysis here. Use is made of the IPSL5 model outputs for SLP and with an rcp 4.5 W/m2 scenario, considering the principal component time scores and loading patterns (Figures 10(a)10(c)). The leading SLP mode is the ridged anticyclone-west coast trough pattern explaining 38% of variance in February to April season over the period 1980–2080. The IPSL5 model, which represents the Benguela upwelling zone better than most, predicts a poleward shift of subtropical anticyclones as CO2 rises. With more easterly winds, the west coast trough deepens. The SLP1 loading pattern of IPSL5 (Figure 10(b)) resembles the FLH-SLP map (Figure 9(e)). The upward trend of its time score (18% of variance, Figure 10(a)) suggests that the frequency of red tide events could grow in the Southern Benguela, as the wind shadow widens. Concurrently, the easterly winds raise air temperatures over the upwelling zone (Figure 10(c)), giving rise to stronger inversions that promote cyclonic wind shear. These changes are underway, drawing the fisheries southward along the coast [54].

fig10
Figure 10: (a) February–April principal component scores for SLP1 (38%) and (11%) based on IPSL simulation and (b, c) their loading patterns. Linear trend for SLP1 refers to increasing ridging high, coastal trough pattern, contributing to warm easterly winds over the Southern Benguela.

4. Summarizing Discussion

This exploration of environmental influences on red tide initially considered a wide variety of spatial and temporal scales and satellite proxies. Following statistical tests, the 8-day fluorescence line height averaged 34–31.5°S, 17.5–18.5°E appeared to be the best indicator and yields the most in situ confirmations. Local ship and aerial surveys have found that red tides form in ~5 × 20 km patches and accumulate over a few days [1416], mainly in February to April season. Efforts to detect and predict red tides have focused on local observations and nested models over multiday events. Here, 30–50 km monthly reanalysis fields have identified the key environmental forcing. The large scale weather is comprised of ridging anticyclone and west coast trough which generates sustained easterly winds over the Cape Peninsula and cold upwelling plumes that extend >100 km seaward (Figures 3(b) and 4(a)4(c)). As the weather evolves into a quiescent phase, there is poleward Ekman transport and inshore warming. The composite SLP anomaly south of Cape Town reaches +9 mb during high FLH cases. The climate “envelope” contains repeated red tide weather events (compare Figures 5(b) and 9(e)) involving a retreat of the circumpolar westerlies and an atmospheric ridge south of Cape Town (Figure 9(d)). These features can be forecast by global models at seasonal lead times.

The windward cape-leeward bay geography creates upwelling shadows where increased residence time builds phytoplankton biomass and harmful dinoflagellate blooms [5559]. Alongshore winds and currents accelerated by the cape sweep clockwise into the leeward bay (Figure 3(f)). The upwelling shadow promotes onshore poleward currents (Figure 8(c)) and thermal stability (Figures 8(a) and 8(b)) which transports and concentrates red tides near the coast. As they dissipate and decay, anoxic water stresses the marine life. The economic consequences of red tides could be mitigated through improved forecasts of its environmental forcing as outlined here. Further progress can be made by employing higher resolution daily products such as the 8 km HYCOM-NCODA model assimilation from http://apdrc.soest.hawaii.edu/las8.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The data suppliers are acknowledged: MCM (red-tide cases), NCEP-NOMAD (CFS daily), NASA-Giovanni (SeaWiFS/MODIS color, Meteosat IR), IRI Climate Library (SODA, ECMWF, and MODIS Ts), Climate Explorer KNMI (CFS monthly, correlation, and PC analysis), APDRC-Hawaii (J-OFURO winds), NOAA-AMDAR (aircraft profiles), and NESDIS-MANATI (scatterometer winds).

References

  1. G. Nelson and L. Hutchings, “The Benguela upwelling area,” Progress in Oceanography, vol. 12, no. 3, pp. 333–356, 1983. View at Google Scholar · View at Scopus
  2. J. Rogers and J. M. Bremner, “The Benguela ecosystem part I. Evolution of the Benguela, physical features and processes,” Oceanography and Marine Biology: An Annual Review, vol. 23, pp. 105–182, 1995. View at Google Scholar · View at Scopus
  3. S. L. Garzoli and A. L. Gordon, “Origins and variability of the Benguela current,” Journal of Geophysical Research, vol. 101, no. 1, pp. 897–906, 1996. View at Publisher · View at Google Scholar · View at Scopus
  4. L. V. Shannon and G. Nelson, “The Benguela: large scale features and processes and system variability,” in The South Atlantic: Present and Past Circulation, G. Wefer, W. H. Berger, G. Siedler, and D. J. Webb, Eds., pp. 163–210, Springer, Berlin, Germany, 1996. View at Google Scholar
  5. M. R. Jury, “The spatial distribution of wind driven upwelling off the West coast of Africa,” Continental Shelf Research, vol. 8, no. 11, pp. 1257–1271, 1988. View at Google Scholar · View at Scopus
  6. A. Bakun and C. S. Nelson, “The seasonal cycle of wind-stress curl in subtropical eastern boundary current regions,” Journal of Physical Oceanography, vol. 21, no. 12, pp. 1815–1834, 1991. View at Google Scholar · View at Scopus
  7. W. Fennel, “Theory of the Benguela upwelling system,” Journal of Physical Oceanography, vol. 29, no. 2, pp. 177–190, 1999. View at Google Scholar · View at Scopus
  8. M. R. Jury, “Mesoscale variations in summer wind over the Cape Columbine-St Helena Bay region,” South African Journal of Marine Science, vol. 3, no. 1, pp. 77–88, 1985. View at Google Scholar
  9. G. C. Pitcher, P. Monteiro, and A. Kemp, “The potential use of a hydrodynamic model in the prediction of harmful algal blooms in the Southern Benguela,” in Harmful Algae 2002, K. A. Steidinger, J. H. Landsberg, C. R. Tomas, and G. A. Vargo, Eds., pp. 11–13, Florida Fish and Wildlife Conservation Commission Florida Institute of Oceanography, Intergovernmental Oceanographic Commission of UNESCO, 2004. View at Google Scholar
  10. C. J. Holden, “Currents in St Helena Bay inferred from radio-tracked drifters,” in South African Ocean Colour and Upwelling Experiment, L. V. Shannon, Ed., pp. 97–109, Sea Fisheries Research, 1985. View at Google Scholar
  11. R. Lamberth and G. Nelson, “Field and analytical drogue studies applicable to the St Helena Bay area off South Africa’s West coast,” in The Benguela and Comparable Ecosystems, A. L. Payne, J. A. Gulland, and K. H. Brink, Eds., vol. 5, pp. 163–169, African Journal of Marine Science, 1987. View at Google Scholar
  12. P. Penven, C. Roy, A. C. de Verdière, and J. Largier, “Simulation of a coastal jet retention process using a barotropic model,” Oceanologica Acta, vol. 23, no. 5, pp. 615–634, 2000. View at Publisher · View at Google Scholar · View at Scopus
  13. G. C. Pitcher and D. Calder, “Harmful algal blooms of the Southern Benguela current: a review and appraisal of monitoring from 1989–1997,” South African Journal of Marine Science, no. 22, pp. 255–271, 2000. View at Google Scholar · View at Scopus
  14. P. L. Donaghay and T. R. Osborn, “Toward a theory of biological-physical control of harmful algal bloom dynamics and impacts,” Limnology and Oceanography, vol. 42, no. 5, pp. 1283–1296, 1998. View at Google Scholar · View at Scopus
  15. G. C. Pitcher, A. J. Richardson, and J. L. Korrubel, “The use of sea temperature in characterizing the mesoscale heterogeneity of phytoplankton in an embayment of the Southern Benguela upwelling system,” Journal of Plankton Research, vol. 18, no. 5, pp. 643–657, 1996. View at Google Scholar · View at Scopus
  16. G. C. Pitcher, A. J. Boyd, D. A. Horstman, and B. A. Mitchell-Innes, “Subsurface dinoflagellate populations, frontal blooms and the formation of red tide in the Southern Benguela upwelling system,” Marine Ecology Progress Series, vol. 172, pp. 253–264, 1998. View at Google Scholar · View at Scopus
  17. C. J. Reason and M. R. Jury, “On the generation and propagation of the Southern African coastal low,” Quarterly Journal of the Royal Meteorological Society, vol. 116, no. 495, pp. 1133–1151, 1990. View at Google Scholar · View at Scopus
  18. M. R. Jury and G. B. Brundrit, “Temporal organization of upwelling in the Southern Benguela ecosystem by resonant coastal trapped waves in the ocean and atmosphere,” South African Journal of Marine Science, vol. 12, pp. 219–224, 1992. View at Google Scholar · View at Scopus
  19. T. A. Probyn, G. C. Pitcher, P. Monteiro, A. J. Boyd, and G. Nelson, “Physical processes contributing to harmful algal blooms in Saldanha Bay, South Africa,” South African Journal of Marine Science, no. 22, pp. 285–297, 2000. View at Google Scholar · View at Scopus
  20. U. Send, R. C. Beardsley, and C. D. Winant, “Relaxation from upwelling in the Coastal Ocean Dynamics Experiment,” Journal of Geophysical Research, vol. 92, no. C2, pp. 1683–1698, 1987. View at Publisher · View at Google Scholar
  21. G. C. Pitcher, J. Agenbag, D. Calder, D. Horstman, M. R. Jury, and J. Taunton-Clark, “Red tides in relation to the meteorology of the Southern Benguela upwelling systemp,” in Harmful Marine Algal Blooms, P. Lassus, G. Arzul, E. Erard, P. Gentien, and C. Marcaillou, Eds., pp. 657–662, Intercept, 1995. View at Google Scholar
  22. J. J. Cullen and J. G. MacIntyre, “Behaviour, physiology and the niche of depth-regulating phytoplankton,” in Physiological Ecology of Harmful Algal Blooms, D. M. Anderson, A. D. Cembella, and G. M. Hallegraeff, Eds., vol. G 41 of NATO ASI, pp. 559–579, Springer, Berlin, Germany, 1998. View at Google Scholar
  23. M. Estrada and E. Berdalet, “Effects of turbulence on phytoplankton,” in Physiological Ecology of Harmful Algal Blooms, D. M. Anderson, A. D. Cembella, and G. M. Hallegraeff, Eds., vol. G 41 of NATO ASI, pp. 601–618, Springer, 1998. View at Google Scholar
  24. T. J. Smayda and C. S. Reynolds, “Community assembly in marine phytoplankton: application of recent models to harmful dinoflagellate blooms,” Journal of Plankton Research, vol. 23, no. 5, pp. 447–461, 2001. View at Google Scholar · View at Scopus
  25. G. C. Pitcher and A. J. Boyd, “Across-shelf and alongshore dinoflagellate distributions and the mechanisms of red tide formation within the Southern Benguela upwelling system,” in Harmful and Toxic Algal Blooms, T. Yasumoto, Y. Oshima, and Y. Fukuyo, Eds., pp. 243–246, Intergovernmental Oceanographic Commission of UNESCO, 1996. View at Google Scholar
  26. A. C. Cockcroft, D. S. Schoeman, G. C. Pitcher, G. W. Bailey, and D. C. VanZyl, “A mass stranding, or “walkout” of West coast rock lobster Jasus lalandii in Elands Bay, South Africa: causes, results and implications,” in Biodiversity Crises and Crustacea, J. C. von Kaupel Klein and F. R. Schram, Eds., Crustacean Issues, pp. 362–688, 2000. View at Google Scholar
  27. G. C. Pitcher and T. A. Probyn, “Anoxia in Southern Benguela during the autumn of 2009 and its linkage to a bloom of the dinoflagellate Ceratium balechii,” Harmful Algae, vol. 11, pp. 23–32, 2011. View at Publisher · View at Google Scholar · View at Scopus
  28. S. J. Weeks, G. C. Pitcher, and S. Bernard, “Satellite monitoring of the evolution of a coccolithophorid bloom in the Southern Benguela upwelling system,” Oceanography, vol. 17, no. 1, pp. 83–89, 2004. View at Publisher · View at Google Scholar · View at Scopus
  29. L. Shen, X. Huiping, and X. Guo, “Satellite remote sensing of harmful algal blooms and a potential synthesized framework,” Sensors, vol. 12, no. 6, pp. 7778–7803, 2012. View at Publisher · View at Google Scholar · View at Scopus
  30. R. M. Letelier and M. R. Abbott, “An analysis of chlorophyll fluorescence algorithms for the moderate resolution imaging spectrometer (MODIS),” Remote Sensing of Environment, vol. 58, no. 2, pp. 215–223, 1996. View at Publisher · View at Google Scholar · View at Scopus
  31. C. M. Hu, F. E. Muller-Karger, C. J. Taylor et al., “Red tide detection and tracing using MODIS fluorescence data: a regional example in SW Florida coastal waters,” Remote Sensing of Environment, vol. 97, no. 3, pp. 311–321, 2005. View at Publisher · View at Google Scholar · View at Scopus
  32. M. C. Tomlinson, T. T. Wynne, and R. P. Stumpf, “An evaluation of remote sensing techniques for enhanced detection of the toxic dinoflagellate, Karenia brevis,” Remote Sensing of Environment, vol. 113, no. 3, pp. 598–609, 2009. View at Publisher · View at Google Scholar · View at Scopus
  33. D. Z. Zhao, X. G. Xing, Y. G. Liu, J. H. Yang, and L. Wang, “The relation of chlorophyll-a concentration with the reflectance peak near 700 nm in algae-dominated waters and sensitivity of fluorescence algorithms for detecting algal bloom,” International Journal of Remote Sensing, vol. 31, no. 1, pp. 39–48, 2010. View at Publisher · View at Google Scholar · View at Scopus
  34. T. Kuster, L. Metsamaa, E. Vahtmae, and N. Strömbeck, “On suitability of MODIS 250 m resolution band data for quantitative mapping cyanobacterial blooms,” Journal of Experimental Marine Biology and Ecology, vol. 55, pp. 318–328, 2006. View at Google Scholar
  35. Y. H. Ahn and P. Shanmugam, “Detecting the red tide algal blooms from satellite ocean color observations in optically complex Northeast-Asia coastal waters,” Remote Sensing of Environment, vol. 103, no. 4, pp. 419–437, 2006. View at Publisher · View at Google Scholar · View at Scopus
  36. P. S. Schopf and A. Loughe, “A reduced gravity isopycnal ocean model: Hindcasts of El Nino,” Monthly Weather Review, vol. 123, no. 9, pp. 2839–2863, 1995. View at Google Scholar
  37. W. W. Gregg, “Assimilation of SeaWiFS ocean chlorophyll data into a three-dimensional global ocean model,” Journal of Marine Systems, vol. 69, no. 3-4, pp. 205–225, 2008. View at Publisher · View at Google Scholar
  38. P. Cury and C. Roy, “Optimal environmental window and pelagic fish recruitment success in upwelling areas,” Canadian Journal of Fisheries and Aquatic Sciences, vol. 46, no. 4, pp. 670–680, 1989. View at Google Scholar · View at Scopus
  39. L. Hutchings, M. R. Roberts, and H. M. Verheye, “Marine environmental monitoring programmes in South Africa: a review,” South African Journal of Science, vol. 105, no. 3-4, pp. 94–102, 2009. View at Google Scholar · View at Scopus
  40. S. E. Nicholson, “A low-level jet along the Benguela coast, an integral part of the Benguela current ecosystem,” Climatic Change, vol. 99, no. 3, pp. 613–624, 2010. View at Publisher · View at Google Scholar · View at Scopus
  41. M. R. Jury, “Representation of coastal upwelling in the Southern Benguela by satellite era reanalysis,” International Journal of Marine Science, vol. 3, no. 34, pp. 267–277, 2013. View at Publisher · View at Google Scholar
  42. H. Demarcq, “Trends in primary production, sea surface temperature and wind in upwelling systems (1998–2007),” Progress in Oceanography, vol. 83, no. 1–4, pp. 376–385, 2009. View at Publisher · View at Google Scholar · View at Scopus
  43. P. Fréon, M. Barange, and J. Arístegui, “Eastern boundary upwelling ecosystems: integrative and comparative approaches,” Progress in Oceanography, vol. 83, no. 1–4, pp. 1–14, 2009. View at Publisher · View at Google Scholar · View at Scopus
  44. S. Saha, S. Moorthi, and P. Hua-Lu, “The NCEP climate forecast system reanalysis,” Bulletin of the American Meteorological Society, vol. 91, pp. 1015–1057, 2010. View at Publisher · View at Google Scholar
  45. J. A. Carton and B. S. Giese, “A reanalysis of ocean climate using simple ocean data assimilation (SODA),” Monthly Weather Review, vol. 136, no. 8, pp. 2999–3017, 2008. View at Publisher · View at Google Scholar · View at Scopus
  46. S. M. Griffies, M. J. Harrison, R. C. Pacanowski, and A. Rosati, “Technical guide to MOM4,” Technical Report 5, GFDL Ocean Group, Princeton, NJ, USA, 2004. View at Google Scholar
  47. R. W. Reynolds, T. M. Smith, C. Liu, D. B. Chelton, K. S. Casey, and M. G. Schlax, “Daily high-resolution-blended analyses for sea surface temperature,” Journal of Climate, vol. 20, no. 22, pp. 5473–5496, 2007. View at Publisher · View at Google Scholar · View at Scopus
  48. M. Kubota, N. Iwasaka, S. Kizu, M. Konda, and K. Kutsuwada, “Japanese ocean flux data sets with use of remote sensing observations (J-OFURO),” Journal of Oceanography, vol. 58, no. 1, pp. 213–225, 2002. View at Publisher · View at Google Scholar · View at Scopus
  49. K. E. Taylor, R. J. Stouffer, and G. A. Meehl, “An overview of CMIP5 and the experiment design,” Bulletin of the American Meteorological Society, vol. 93, no. 4, pp. 485–498, 2012. View at Publisher · View at Google Scholar
  50. M. R. Jury, “Climate trends in Southern Africa,” South African Journal of Science, vol. 109, no. 1-2, article 980, 11 pages, 2013. View at Publisher · View at Google Scholar · View at Scopus
  51. J. L. Dufresne, M. A. Foujols, S. Denvil et al., “Climate change projections using the IPSL-CM5 earth system model: from CMIP3 to CMIP5,” Climate Dynamics, vol. 40, no. 9-10, pp. 2123–2165, 2013. View at Publisher · View at Google Scholar · View at Scopus
  52. C. Roy, S. Weeks, M. Rouault, G. Nelson, R. Barlow, and C. van der Lingen, “Extreme oceanographic events recorded in the Southern Benguela during the 1999-2000 summer season,” South African Journal of Science, vol. 97, no. 11-12, pp. 465–471, 2001. View at Google Scholar · View at Scopus
  53. G. C. Pitcher and G. Nelson, “Characteristics of the surface boundary layer important to the development of red tide on the Southern Namaqua shelf of the Benguela upwelling system,” Limnology and Oceanography, vol. 51, no. 6, pp. 2660–2674, 2006. View at Publisher · View at Google Scholar
  54. C. Roy, C. van der Lingen, J. C. Coetzee, and J. Lutjeharms, “Abrupt environmental shift associated with changes in the distribution of Cape anchovy Engraulis encrasicolus spawners in the Southern Benguela,” African Journal of Marine Science, vol. 29, no. 3, pp. 309–319, 2007. View at Publisher · View at Google Scholar · View at Scopus
  55. W. M. Graham and J. L. Largier, “Upwelling shadows as nearshore retention sites: the example of Northern Monterey Bay,” Continental Shelf Research, vol. 17, no. 5, pp. 509–532, 1997. View at Publisher · View at Google Scholar · View at Scopus
  56. L. B. Joyce, G. C. Pitcher, A. D. Randt, and P. M. S. Monteiro, “Dinoflagellate cysts from surface sediments of Saldanha Bay, South Africa: an indication of the potential risk of harmful algal blooms,” Harmful Algae, vol. 4, no. 2, pp. 309–318, 2005. View at Publisher · View at Google Scholar · View at Scopus
  57. M. Roughan, A. J. Mace, J. L. Largier, S. G. Morgan, J. L. Fisher, and M. L. Carter, “Subsurface recirculation and larval retention in the lee of a small headland: a variation on the upwelling shadow theme,” Journal of Geophysical Research, vol. 110, no. 10, Article ID C10027, pp. 1–18, 2005. View at Publisher · View at Google Scholar · View at Scopus
  58. G. C. Pitcher and S. J. Weeks, “The variability and potential for prediction of harmful algal blooms in the Southern Benguela ecosystem,” in The Benguela: Predicting a Large Marine Ecosystem, V. Shannon, G. Hempel, P. Malanotte-Rizzoli, C. Moloney, and J. Woods, Eds., pp. 125–146, Elsevier, 2006. View at Google Scholar
  59. G. C. Pitcher, F. G. Figueiras, B. M. Hickey, and M. T. Moita, “The physical oceanography of upwelling systems and the development of harmful algal blooms,” Progress in Oceanography, vol. 85, no. 1-2, pp. 5–32, 2010. View at Publisher · View at Google Scholar · View at Scopus