Relation between Upwelling Intensity and the Variability of Physical and Chemical Parameters in the Southern Benguela Upwelling System

The extent to which wind-driven seasonal upwelling cycles manifest in surface ocean temperature and nutrient variability along a monitoring line in the Southern Benguela upwelling system was investigated. Monitoring conducted monthly over a six-year period shows that surface temperature and nutrient concentrations exhibit very poor seasonality and weak correlation with the upwelling index. This is, despite clear evidence for spatial inshore-offshore gradients in temperature, nutrients, and chlorophyll, consistent with an upwelling regime.The upper ocean temperature gradient shows a much better correspondence to the upwelling index but at the same time demonstrates that surface heating, and not vertical mixing related to upwelling, controls the upper ocean temperature gradient. The results suggest that remote sensing techniques would be inadequate tools to monitor upwelling events in the Southern Benguela. Also, the incidence of phytoplankton blooms is more likely triggered by stratified conditions associated with surface heating than relaxation of upwelling winds.


Introduction
The Benguela upwelling system stretches from South Africa to Angola and is one of the ocean's four large and most productive eastern boundary upwelling areas [1][2][3].It is bound at both the equatorward and poleward ends by warm water regimes, notably the Angola-Benguela Front in the north and the Agulhas Current retroflection area in the south [4].Eastern boundary currents like the Benguela are driven by equatorward alongshore winds, which also force offshore Ekman transport in the surface boundary layer.Along the coast, this divergence results in the upwelling of cold, nutrient-rich water into the euphotic zone.A wind-driven equatorward frontal jet typically forms a boundary between the open ocean and the upwelled water, with mass balance being provided by a poleward undercurrent at depth.The enrichment of surface water inshore of the front supports high productivity and often a high frequency of harmful algal blooms.A global increase in the incidence of such blooms [5], and the role of productive upwelling areas as global carbon sinks have renewed interest in understanding upwelling processes and the possible role of global change as a driving factor behind change in upwelling dynamics.Of particular interest is the application of satellite oceanography to detect upwelling events, in order to improve spatial and temporal data coverage.
St. Helena Bay, located at the southern end of the Benguela Upwelling system, is the best studied and one of the most productive areas in this system [6].Upwelling plumes of cold nutrient-rich South Atlantic Central Water (SACW) are characterised by consistent temperature and salinity (ca 10 ∘ C and 34.80 PSU) [3,[7][8][9].Offshore transport in a northwesterly direction eventually causes mixing with or sinking below the warm, high salinity oceanic water (ca 20 ∘ C and 35.5 PSU) at a convergence front [10,11].In the summer (Dec.-Feb.), the subthermocline waters in the southern Benguela are ventilated by the ocean-shelf exchange of aerated SACW at Cape Columbine and Cape Peninsula at the southern boundary [12].In the winter (June-August), there is a disruption of the thermocline and a complete reaeration of the water column.St. Helena Bay is also a part of the Benguela Upwelling system that is characterised by expansive red tides and with significant variability spatially and temporarily.Red tides typically develop in response to upwelling relaxation and are characterised by a phytoplankton community shift from diatoms to dinoflagellates.This study applies data collected monthly over a 6-year period along the almost 200 km long St. Helena Bay Monitoring Line (SHBML) to investigate covariation between an upwelling index calculated from nearby wind records and physical and chemical properties along the transect.The extent to which the well-documented seasonal upwelling cycles is manifested at the surface along an inshore-offshore was investigated in the context of implications for the validity of remote sensing [13] as a monitoring tool in this ocean area and to improve understanding of the physical drivers of biological processes in St. Helena Bay.

Study Area.
The study area is the St. Helena Bay Monitoring Line (SHBML) in the Southern Benguela, along the west coast of South Africa (Figure 1, Table 1).Stations 1 to 12 have bottom depths of 27 to 1396 m, respectively, and are located at distances ranging from 3 to 191 km offshore (Table 1).The Berg River discharges into St.Helena Bay at Veldrift (Figure 1) and drains a catchment dominated by agriculture, wineries, canneries, and textile milling [14,15].Monthly cruises on the research vessels RV Africana, RV Algoa, and RV Ellen Khuzwayo were conducted along the SHBML from January 2004 to December 2009.The SHBML is typically sampled only once during the two-month long October/ November pelagic spawner biomass survey.Upper ocean profiles for temperature, salinity, and chlorophyll-a (Chl-a)  parameters were measured in situ using a Sea-Bird Electronics SBE 911 plus CTD (conductivity, temperature, and depth) sensor as well as fluorescence sensor.Water temperature was measured to an accuracy of ±0.001 ∘ C. The temperature gradient was used to indicate the presence or absence of a thermocline.Salinity was computed from CTD data in practical salinity units (PSUs) and calibrated using discrete samples (accuracy 0.0003 Siemens/meter).
Water samples for the study of dissolved inorganic nutrients were stored in acid-washed polyethylene bottles with pressure caps and kept frozen at −80 ∘ C until analysis ashore.No samples were kept frozen for longer than 3 months before analysis.Dissolved nitrate (NO 3 − ), phosphate (PO 4 3− ), and silicates were determined according to the methods described by [17].An Astoria Analyzer Series 300 expanded to four channels was used to perform the nutrient analysis.
Samples for Chl-a analysis were taken at four depths (near-surface, above the fluorescence maximum (F-max), Fmax, and below the F-max).Subsamples (200 mL) of seawater for Chl-a analysis were collected in prerinsed plastic measuring cylinders.Samples were filtered under vacuum onto 25 mm Whatman GF/F glass fibre filter papers [18].The filter papers were frozen in aluminium foil pouches for analysis ashore.Chl-a was measured fluorometrically on a Turner Designs 10-AU fluorometer after extraction in 90% acetone [19].The fluorometer was calibrated with Chl-a standard (Sigma Chemical Co., USA) in 90% acetone solution with a GBC Cintra 404 spectrophotometer and an extraction coefficient of 87.67 Lg −1 cm −1 .These warm water events are usually, but not always, also observed at 20 m depth (Figure 3(a)).

Relationship between the Upwelling Index (UI) and
Linear regression coefficients between temperature and the upwelling index, which has a very well defined seasonal profile (black line in Figures 2(a) and 2(b)), confirm the poor seasonality visually evident in the near-surface temperature data (Table 2).The expected trend is for surface temperature to decrease in response to upwelling (more positive UI values) and for this to be more pronounced closer to shore, where divergence is more likely to occur.However, near-surface temperature is observed to be very poorly correlated with the UI at the surface (approximately 3 to 5 m depth), except at the furthest offshore stations (11 and 12), where a weak positive correlation with the UI is observed.At 20 m depth the expected negative relationship between temperature and the UI is evident at stations 1 to 5 (inshore stations), with again a positive relationship with the UI manifesting at the offshore stations (11 and 12), with the latter slightly weaker than the correlation at the near surface.
There is a more pronounced relationship between the upwelling index and near-surface salinity at the inshore stations than observed for temperature (Table 2, Figures 2(c) and 2(d)).The data suggests a trend of decreasing salinity values at times of upwelling, at both the near-surface and 20 m depth, slightly better expressed at 20 m depth.Seasonal variations in salinity are well-defined from 2004 to 2006 (Figure 3(b)).In 2007, however, the seasonal profile is disrupted by the intrusion of a water mass with typical wintertime (June-August) temperature (13.16 ∘ C) but much lower salinity (34.41 PSU) values than is the norm in July (Figure 3), followed by a return to more typical upwelling conditions the following month, that is, colder (11.24 ∘ C) more saline (34.80 PSU) water.Similar but less pronounced intrusions of water masses with lower salinities than normally observed in winter (June-August) disrupt the seasonal salinity profiles in 2008 and 2009 (Figure 3(b)).These fresher-water intrusions are also much more evident at near-surface sampling depths than at 20 m depth and in 2008, for example, it is only evident at the near surface.

Relationship between the UI and Surface Ocean Nutrient and
Chlorophyll-a Concentrations.Near-surface ocean distribution patterns for nutrients exhibit across-shelf gradients, in Figures 4(a) and 4(b)).This inshore-offshore gradient in nutrient availability is reflected in higher inshore Chl-a levels (above 40 mg/m 3 in some instances), that is, enhanced coastal productivity fuelled by the higher availability of nutrients (Figures 4(c) and 4(d)).There is a general trend of lower nutrient and Chl-a values during the winter (June-Aug) along the SHBML and higher nutrient and Chl-a during the upwelling season (Figures 4 and 5).However, similar to the observations made for temperature, the annual cycles are disrupted by numerous events, during which nutrient and Chl-a levels are contrary to this generalized view of surface ocean nutrient dynamics.The result is variability that exhibits very poor seasonality (Figure 5).The relationship between the UI and the nutrient parameters dissolved NO 3 − , PO 4 3− , and Si is very weak at the near surface but is slightly better expressed at 20 m depth, where the relationship with the UI is positive, as expected in a upwelling area (Table 2, Figures 4(a) and 4(b)).Of the three nutrient proxies, PO 4 3− shows the strongest relationship with the UI (Table 2).Chl-a levels exhibit a weak but positive relationship with the UI at both sampling depths and all stations (Table 2, Figures 4(c) and 4(d)).
Nutrient and Chl-a variability associated with the warm water events that disrupt the seasonal temperature profile, and the fresher water events that disrupt the seasonal salinity profile, bear further scrutiny.The anomalous warm water events observed at station 1 in the summer (Dec.-Feb.) of 2005, the summer (Dec.-Feb.) of 2006, late spring (Nov.) in 2007, and in spring (Sept.-Nov.) in 2009 are all characterized by much lower than usual surface nutrient levels (0.17 to 2.41 mol/L for NO 3 − and 0.26 to 1.21 mol/L for PO 4 3− ).The most pronounced of these, during the summer (Dec-Feb) of 2005 was also associated with very high Chl-a levels (33.4 mg/m 3 ) and temperature-salinity properties ( = 17.38 ∘ C and  = 34.79PSU) that rule out the intrusion of an offshore warm-nutrient-depleted water mass (Figure 2) as an explanation for this event.The high temperature suggests thermal heating as the cause of the elevated temperature, and the low nutrient and high Chl-a levels suggest the existence of an algal bloom, caused by the existence of stratified conditions conducive to enhanced phytoplankton productivity.The combined temperature-salinity-nutrient-Chla properties that characterize the other warm water events that are evident in the temperature record (Figure 3(a)) are more ambiguous.Although conditions similar to those which appear to have caused the summer (Dec-Feb) 2005 warm event cannot be ruled out, such conditions also cannot be confirmed with the data available.The most pronounced low-salinity event observed occurred in the winter (June-August) of 2007 (Figures 2(c) and 3(b)).The temperature, nutrient, and Chl-a levels associated with this event, however, are not unusually high or low, with the exception of lower than usual dissolved Si levels (5.42 mol/L).An explanation for the origin of the lowsalinity water observed in early spring (September) in 2008 and later winter (August) 2009 is also difficult to extract from the available data.Although continental run-off from winter (June-August) rainfall over the adjacent landmass cannot be ruled out as a contributing factor, there is nothing in the data available to substantiate such a hypothesis at this stage.

Relationship between the UI and Vertical
Temperature and Nutrient Gradients.Evaluation of variability in the upper ocean temperature gradient (shown for station 1 in Figure 6), in comparison with seasonal variability in the upwelling index, shows that the temperature gradient exhibits better defined seasonal variability than surface ocean temperature.Along the SHBML the water column is more stratified during spring (Sept.-Nov.)and summer (Dec.-Feb.)than during the winter (June-August).Generally speaking, surface temperatures and surface-deep water temperature gradients are determined by a balance between surface heating and windinduced vertical mixing.Along the SHBML the net surface heat flux is positive and is increasing during the spring-summer upwelling season.The water column stabilization and stratification that result from this counter the tendency for vertical mixing induced by upwelling-favorable winds during the spring-summer season, during most years (Figure 6).During the winter (June-August) the water column is better mixed, despite the prevailing downwelling-favorable winds during this season.
The balance between surface heating and vertical mixing is clearly in favor of surface heating in St. Helena Bay.This has implications not only for the application of remote sensing methods to track upwelling events, but also for biological processes and spatial and temporal variability therein.Planktonic organisms tend to be distributed throughout the water column when it is well mixed, which is winter-time along the SHBML.When the water column is stratified, that is, spring-summer along the SHBML, it favors the development of algal blooms.Also, this concentration of food particles in a stratified water column is believed to be advantageous to the growth of larval fish [20].

Conclusions
It is clear from the results that all physical and chemical parameters exhibit pronounced across-shelf variability along the SHBML in the surface ocean, but very poor seasonality and no significant relationship with the calculated upwelling index.This does not mean that the observed weak relationships do not provide valuable insight into upwelling processes in St. Helena Bay.It does mean, however, that the weak relationship between wind-induced upwelling and inner shelf temperature cast doubt on the usefulness of satellite observations in detecting upwelling events in this part of the southern Benguela Upwelling system, as reflected either in decreased sea-surface temperatures or in surface Chl-a distributions.It also poses a problem for using moored offshore buoys as detectors of upwelling events.The tabulated regression coefficients for the UI and temperature (Table 2) suggest that station 2 (located approximately 7 km offshore) would be the best location to moor such a buoy and that the focus should be on temperature variability at 20 m depth rather than at the surface.
The observation that surface heating and water column stratification resulting from it outbalance wind-induced vertical mixing during the upwelling season in St. Helena Bay has important implications.First of all, it explains the poor correlation between sea surface temperature variability and the upwelling index, as well as the poor correlation between the UI and nutrient parameters.Secondly, it implies that remote sensing techniques will be poor indicators of upwellingrelated variability in surface temperature and chlorophyll distributions in St. Helena Bay.It is also a critical observation in regard to understanding the factors that contribute to harmful algal blooms and the increasing incidence of such events.

Figure 1 :
Figure 1: The SHBML station positions superimposed on bathymetry.Shallow depths are shown in green and deeper depths in blue tones.

Figure 2 :
Figure 2: Inshore-offshore near-surface (a and c) and 20 m depth (b and d) variability over time for temperature (a and b) and salinity (c and d).The thick black line represents the upwelling index (UI).Temperature ( ∘ C) and salinity (PSU) values are indicated as labels on respective isobars.The temperature ranged between <11 ∘ C (dark blue) and >18 ∘ C (red).Salinity ranged between <34.7 PSU (dark blue) and >35.4 (red).

Figure 3 :
Figure 3: Temporal changes in (a) temperature and (b) salinity at station 1, near surface and at 20 m depth.

Figure 4 :
Figure 4: Inshore-offshore near-surface (a and c) and 20 m depth (b and d) variability over time for dissolved phosphate (a and b) and Chl-a (c and d).Dissolved phosphate (mol/L) and Chl-a (mg/m 3 ) values are indicated as labels on respective isobars.The dissolved phosphate ranged between <0.5 mol/L (grey) and >2 mol/L (dark orange).Chl-a ranged between <4 mg/m 3 (white) and >20 mg/m 3 (dark green).

Table 1 :
St. Helena Bay Monitoring Line sampling station details.

Table 2 :
Linear regression coefficients (slopes) for physical and chemical variables as a function of the upwelling index at the surface (3 to 5 m) and at 20 m sampling depths, for stations 1 to 12. Bold values indicate slopes that are significant at  < 0.01.