Abstract

This study considers environmental factors influencing aggregate fish catch in the South Coast transition of South Africa. The environmental forcing is studied via (i) seasonal analysis of SeaWifs chlorophyll and related variables, (ii) composite analysis of atmospheric and oceanographic reanalysis data, (iii) statistical analysis of annual FAO fish catch with climatic indices, and (iv) analysis of depth-latitude hydrographic sections over the shelf (33–36S, 22–26E). In years of higher fish catch there is a northward shift of the subtropical anticyclones and upwelling that is partially related to Pacific El Nino. Westerly troughs skirt the Agulhas Bank creating onshore Ekman transport. Higher sea surface height inshore, and cooler sea temperatures and lower salinity offshore induce a gradient that weakens the Agulhas Current. These environmental conditions favour the southeastward migration of juvenile fish from west to south coast. A multivariate model of aggregate fish stocks, using four environmental variables: salinity and zonal currents in the Agulhas Current, sea temperature in the Agulhas source region, and geopotential height over the Cape, accounts for 53% of variance at 0-1 year lead. Freshening of the boundary current is a factor influencing aggregate fish catch in the South Coast transition.

1. Introduction

Marine fish catch fluctuates in relation to human pressure, climatic conditions, and physiological requirements [110]. Long-term records on annual catch from the Food and Agricultural Organization [11, 12] exhibit quasidecadal cycles often with short bursts of high catch interspersed with longer quiescent periods [1317]. Although fisheries often show climate dependence, the linkages are indirect because of complex feedbacks, lags between fishing effort and natural production, and weakness of the data [18, 19].

Retrospective studies that account for slow variations of the ocean and the life history and migration patterns of the fish have been conducted in the past 50 years using physical oceanographic fields [20, 21]. The alternation of sardine and anchovy landings [1, 22] has been attributed to environmental conditions that affect spawning and feeding location [2325]. Because of the worldwide extent and quasidecadal time scale of fishery fluctuations, global ocean-atmosphere coupling has been proposed as a driver [26] through features such as the Pacific decadal oscillation (PDO), the North Atlantic oscillation (NAO) [27, 28], and the El Niño southern oscillation (ENSO) [21, 29, 30]. Resolving environmental influences on fisheries will contribute to better projections and management of fish stock fluctuations. Following from this discussion, key questions include: to what extent are regional fisheries affected by the surrounding environmental conditions? Are environmental factors within the catch region/during the catch year more important than remote/precursor teleconnections? Which atmosphere/ocean elements are most influential?

The marine fish catch along the southern tip of Africa (Figure 1) is comprised of short-lived shelf-zone pelagic and longer-lived deep-water species. Much is known about the region’s oceanography [31] characterized by seasonal upwelling and abundant populations of sardine and anchovy. The fishery responds to the annual cycle [5] anchovy spawn on the western Agulhas Bank in early summer when high wind speeds are prevalent, while sardines spawn over a broader season and area [32, 33]. Spawning over the eastern Agulhas Bank can result in losses by strong currents. Success rates are highest when shelf temperatures are 16–19°C. Eggs are transported from the Agulhas Bank to the St Helena Bay nursery area (chlorophyll max., Figure 1) by frontal currents [34]. Nursery area conditions are also important: upwelling is needed to maintain primary productivity, but not so much that larvae are exported. Juveniles then migrate back to the Agulhas Bank for recruitment, following various paths southeastward along the shelf. Roy et al. [35] suggest that the pelagic fishery shifted eastward after 1995 due to cooling over the shelf from enhanced easterly winds and seasonal upwelling. Here an alternative hypothesis is explored that the fishery is modulated by changes in the boundary current [36]; specifically, when the Agulhas Current is cool, fresh, and slow and westerly winds produce onshore Ekman transport, larvae retention and juvenile recruitment are favoured. With ocean-atmosphere reanalyses and historical catch information this hypothesis is tested for the Agulhas Bank. Section 2 covers the data and methods, while Section 3 presents the results divided into seasonal and fishery variability, large-scale forcing, local conditions, and climatic change. Section 4 provides conclusions and some recommendations.

2. Data and Methods

This study considers oceanic and atmospheric variability at annual time scales around Southern Africa (0–40S, 4–52E), the upstream Agulhas Current, and the Agulhas Bank (33–36S, 20–28E, Figure 1). The methods seek to understand what environmental conditions influence fish catch, in addition to the effects of abundance, management, and effort. Annual fish catch statistics in the period 1976–2006 are drawn from the FAO fisheries data base. The major species oscillate together as indicated by cross-correlations between an aggregate fish index and the individual species: hake comprising half the catch has , pelagic comprising 17%   , sardines and anchovy comprising 12%   , and demersal comprising 12%   , where is significant at 90% confidence for . The sardine and anchovy catch has risen since 1998, while the demersal fishery has declined. Hake remained steady, and the pelagic catch (including sardine and anchovy before 1998) tends to oscillate. Since the various fisheries fluctuate together and details of individual species are uncertain in the early years, an aggregate fish catch index is developed for the Agulhas Bank. The annual fish catch is divided by the normalized population as a way of accounting for increasing fishing effort. Analyses were also made per species, but these support the aggregate results and were not further considered.

Atmospheric conditions are represented primarily by the National Center for Environmental Prediction (NCEP) reanalysis product [37]. These 2° gridded fields are based on surface ship and coastal data, optimally interpolated using a numerical weather prediction model and enhanced by satellite remote sensing since 1979. Because the south coast of South Africa has a busy shipping lane, the comprehensive ocean atmosphere data set (COADS) is used for local wind time series. With the ocean affected by fresh water inputs from the land and atmosphere, 1° gridded rainfall analyses from the Global Precipitation Climatology Project (GPCP) are considered that involve gauges over land and satellite estimates since 1979. River run-off data are extracted from the South African hydrological services for the Agulhas-facing coast, including the Tugela River at 29S.

The ocean environment is described by the SODAv2.4 reanalysis product of the University of Maryland at 0.5° horizontal and ~30 m vertical resolution [38]. These make use of hydrographic data (from SADCO and GTS), optimally interpolated using a numerical ocean model that employs European Community Medium-range Forecast (ECMWF) wind stress and surface fluxes, satellite thermometry since 1979, and satellite altimetry since 1985. A principal component analysis (PCA) was made to determine the dominant spatial clusters and detrended time scores over the period 1976–2006 for 0–200 m depth-averaged temperature, salinity, currents, vertical motion, and sea surface height in the southern African domain. The top three PCA modes were retained, and their annual time scores were cross-correlated with the FAO fish catch and other climatic indices. Because of the potential ambiguities generated by PCA, environmental time series were extracted for smaller domains, averaged over the Agulhas Current (26–33S, 29–35E) and the Agulhas Bank for comparison with the catch statistics. With the advent of SeaWifs chlorophyll estimates in 1997 that reflect primary productivity [39], an analysis of the seasonal cycle over the Agulhas Bank was done, and comparisons were made with the SODA2.4 ocean reanalysis variables, COADS winds, ECMWF wind stress, and GPCP rainfall. The analysis reveals seasonal amplitudes and processes underlying the Agulhas Bank shelf hydrography and biochemical productivity.

To understand the spatial pattern favoring fish catch, the atmosphere and ocean reanalysis data were composite-averaged to produce high minus low ( ) catch fields. These were done as maps and 0–600 m depth sections south of the Agulhas Bank in the longitudes 22–26E. To construct the composite, the catch record was divided into years of high and low catch in each 6-7-year cycle: 1979, 1985, 1992, 1997, 2005; 1977, 1981, 1988, 1995, 1999. Two domains were mapped: large-scale: 0–50S, 10W–52E and local: 31–38S, 16–28E. Contrasts between the pattern of years with high and low catch were analyzed by field subtraction, and differences were evaluated for: SST, 200 hPa upper wind, GPCP rainfall, ECMWF wind stress, 0–200 m depth-averaged temperature, salinity, currents and 0–500 m vertical motion. Supplementary analyses were made for seasonal periods to understand how the environmental signals vary from summer to winter. Composite H-L differences exceeding one standard deviation are interpreted as significant.

The statistical analysis explores cycles in the fish catch using wavelet spectral methods [40]. The annual atmospheric and oceanic data averaged over the Agulhas Current and Agulhas Bank, the PCA mode time scores, and global climate indices: PDO, NAO, ENSO (Nino3 SST), and river run-off were compared with the annual fish catch time series to determine relationships. Pairwise cross-correlations were evaluated, and a two-year period: catch year + previous (0-1 yr lead) was found to be optimal. Using a stepwise regression procedure to insert successively influential (normalized) environmental variables, a multivariate model was formulated to account for fluctuations in aggregate fish catch. After adjustment for multiple predictors, 26 degrees of freedom requires for significance at 90% confidence.

3. Results

3.1. Ocean Colour and Fish Catch Variability

The seasonal cycle of SeaWifs chlorophyll is shown in Figure 2(a). Primary production peaks during the autumn months March–May with cyclonic wind vorticity (curl). Wind stress is favourable for upwelling ( ) in January–March (Figure 2(b)), yet the surface layer is stratified by high solar radiation. Marked downwelling and warmer sea temperatures are experienced in May-June (Figure 2(c)) as westerly winds sweep over the area [41, 42] inhibiting phytoplankton blooms. There is a weak secondary peak in chlorophyll concentration in September-October, and a broad minimum from November to February (Figure 2(a)). As in most of the subtropics, productivity is confined to the shelf (cf. Figure 1). In the Agulhas Bank, the area of high chlorophyll widens toward the west with the axis of the Agulhas Current forming the outer edge. Cross-correlations of environmental variables (Table 1) indicate that spells of easterly winds induce vertical uplift and cooler fresher conditions, but simultaneous associations with chlorophyll are weak. When chlorophyll is lagged by 2 months, the relationship with vertical motion strengthens. Agulhas Bank shelf waters require stratification following upwelling to generate higher primary productivity. Fluctuations of chlorophyll are small in the food-poor, spring, and summer season and greatest in autumn and tend to follow the trend of mean values (cf. Figure 2(a)).

In Figure 3(a) the time series of FAO marine fish catch in the Agulhas Bank area is presented. Although the annual time series is normalized by population, a 2nd-order rising trend remains that accounts for 42% of variance. The increase in fish catch is attributed to a southward shift of the anticyclonic wind belt and more frequent upwelling favourable winds [35], an issue that is investigated below. There is a clear cycle, confirmed by wavelet spectral analysis to be significant (Figure 3(b)): 6 years in the early part of the record, lengthening to 7 years toward the end. The amplitude of the cycle also grows. In the following sections, the regional and local environmental factors that contribute to changes in Agulhas Bank fish catch are explored.

3.2. Large-Scale Influences

The PCA of ocean reanalysis fields in the South Africa domain reveals a number of distinct patterns, and those significantly correlated with the catch index are described. Figure 4 includes depth-averaged temperature mode-1, salinity mode-1, and atmospheric geopotential height mode-3. The patterns and temporal relationships suggest that cooling and freshening of Agulhas source waters at 0-1 year lead time promote higher fish catch. The mean currents that connect the source waters in the northern Mozambique Channel to the Agulhas Bank are such that it would take a few months for the water mass to arrive. The upper atmospheric geopotential PCA mode-3 pattern and its temporal correlation indicate that lower pressure over southern Africa associates with higher fish catch. The PCA time scores exhibit spectral energy in the range of 3–11 years. No single variable matches fish catch, so multivariate analysis is explored (cf. Section 4).

Composite high minus low catch ( ) atmospheric reanalysis and SST maps are given in Figure 5. A remote signal is found in southern Angola waters. The Cape Frio upwelling plume is colder and extends further north than usual before and during years with high fish catch over the Agulhas Bank. SST differences are −1.2 C around 15S, 10E (Figure 5(a)). More vigorous coastal upwelling in the southeast Atlantic cold tongue is related to an equatorward-displaced subtropical anticyclone. In the upper atmosphere, 200 hPa wind differences are from the west in the tropics and from the east over the SW Indian Ocean (Figure 5(b)). These wind patterns are consistent with Pacific El Nino [4345] that bring dry weather to South Africa, as seen in the rainfall composite map (Figure 5(c)). The zone of reduced rainfall stretches from Angola to the Agulhas Current, causing Tugela River run-off to decline ( with respect to catch, cf. Table 2). Cool SST in the east Atlantic (cf. Figure 5(a)) often precedes Pacific El Nino. The near surface geopotential height map indicates the presence of a low/high pattern in the midlatitudes (L: 50S, 10E/H: 45S, 55E, Figure 5(d)). Analysis of fields at different lags indicates that higher pressure in the south Indian Ocean “blocks” atmospheric troughs over South Africa. Composite easterly winds from south of Madagascar meet westerly winds over the eastern Agulhas Bank. The regional environment thus exhibits key signals in composite analysis.

Time series in the upstream Agulhas Current region (26–33S, 29–35E) exhibits correlations with fish catch: salinity , followed by rainfall −0.40, and zonal currents +0.37 (Table 2). Hence lower salinity and a weakened Agulhas Current benefit the fisheries. Reduced rainfall is associated with higher fish catch and is consistent with the presence of low pressure over the Cape (cf. Figure 4(c)), increased westerly winds, and reduced river run-off along the east coast.

3.3. Composite Shelf Conditions

Ocean reanalysis H–L composites maps and sections over the Agulhas Bank are given in Figures 6 and 7, respectively. The maps cover the recruitment pathway and include the west coast nursery area. There is a significant westerly wind difference over the Agulhas Bank, and offshore winds are evident over the west coast (Figure 6(a)). Ekman transport tends to be poleward on the west coast and onshore along the south coast. Thus juvenile fish seeking recruitment are “helped” by favourable winds and currents. Westerly wind differences are consistent with El Nino [46, 47]; Pacific SST with respect to fish catch (cf. Table 2). The onshore Ekman transport lifts the sea surface height nearshore (Figure 6(d)) inducing eastward current differences over the shelf edge 35–37S, 20–26E (Figure 6(b)). Although wind differences are easterly over the Agulhas Bank during summer (not shown) in agreement with Roy et al. [35], the annual composite (Figure 6(a)) indicates that westerly winds prevail.

The composite temperature and sea surface height patterns indicate that coastal upwelling is reduced (Figures 6(c) and 6(d)), while cooler waters are found offshore. This gradient works against the mean state and weakens the Agulhas Current. Salinity is neutral in the coastal zone, but there is a fresh signal offshore (Figure 6(e)). These features are analyzed in vertical section in Figures 7(a) and 7(b). Cooler temperatures (−1.2 C from 100–300 m) and lower salinity (–0.12 g kg−1 from 250–450 m) are evident beyond the shelf edge (36S). There is anomalous onshore flow and vertical uplift (>0.5 m/day) on the shelf (35S, 25E) in the 200–600 m layer (Figures 6(f) and 7(c)), despite the anomalous westerly winds and onshore Ekman transport near the surface. The uplift appears related to anomalous cyclonic wind stress curl created by westerlies on the Agulhas Bank meeting easterlies from the SW Indian Ocean (Figure 6(a)). Weakening of the Agulhas Current is most evident in the 0–200 m layer around 36S (Figure 7(d)). Prevailing before and during high catch years, these features favour the migration of pelagic species from west to south coast.

4. Discussion

The relative role of local and remote atmosphere and ocean forcing is explored by multivariate stepwise regression of key environmental variables onto the aggregate catch index. The scatterplot is given in Figure 8 and the statistical summary in Table 2. Two PCA time scores and two area-averaged time series are selected in stepwise regression to account for 53% of variance at 0-1 year lead. An important variable is salinity in the upstream Agulhas Current (coefficient −0.53), followed by temperature mode-1 (+0.47), upper atmospheric geopotential height mode-3 (+0.36), and upstream Agulhas zonal current (+0.34). This mix of “predictors” has wavelet spectral energy of 4 yr in the 1976–1994 period and 6 yr since 1995 and a weak upward trend contributed by declining salinity and currents in the Agulhas Current. The scatterplot is relatively symmetrical, so the environmental forcing that induces high catch is opposed for low catch, albeit with half the variance still unexplained.

The ocean-atmosphere environment around the Agulhas Bank has been analyzed in the context of aggregate fish catch, with a focus on contrasts between high and low years. There is an equatorward shift of the atmospheric circulation and upwelling in high catch years. Westerly troughs skirt the Agulhas Bank around a midlatitude low-pressure anomaly (cf. Figure 5(d)) that suppresses rainfall and river discharges over South Africa. The anomalous wind and current regime (cf. Figures 6(a) and 6(b)) favours the migration of juvenile fish from west to south coast, with a cycle of 6-7 years. Higher sea surface height inshore and cool fresh waters offshore induce a gradient that weakens the Agulhas Current. The composite difference maps reveal links between Agulhas Bank fish catch and ocean climate at 0-1 year lead.

A multivariate regression model used regional environmental variables to “hindcast fit” aggregate fish catch and suggested upstream conditions in the Agulhas Current as a key influence. The composite fields gave evidence of an anomalous conveyor belt that favours the migration of juvenile fish from west to south coast. Further work can be directed towards biochemical forcing using ocean reanalysis data on water chemistry and chlorophyll, intraseasonal analyses of recent years with better data, and separation of environmental responses per fish species.

Acknowledgments

The data used here were extracted from websites of the Food and Agricultural Organization, International Research Institute for Climate, the National Center for Environmental Prediction, and the Climate Explorer of the Netherlands Meteorological Institute.