Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Patterns of Distribution and Spatial Indicators of Ecosystem Change Based on Key Species in the Southern Benguela

  • Katherine E. Watermeyer ,

    kwatermeyer@gmail.com

    Affiliation Marine Research Institute (Ma-Re) & Department of Biological Sciences, University of Cape Town, Cape Town, South Africa

  • Laurence Hutchings ,

    Contributed equally to this work with: Laurence Hutchings, Astrid Jarre, Lynne J. Shannon

    Affiliation Marine Research Institute (Ma-Re) & Department of Biological Sciences, University of Cape Town, Cape Town, South Africa

  • Astrid Jarre ,

    Contributed equally to this work with: Laurence Hutchings, Astrid Jarre, Lynne J. Shannon

    Affiliation Marine Research Institute (Ma-Re) & Department of Biological Sciences, University of Cape Town, Cape Town, South Africa

  • Lynne J. Shannon

    Contributed equally to this work with: Laurence Hutchings, Astrid Jarre, Lynne J. Shannon

    Affiliation Marine Research Institute (Ma-Re) & Department of Biological Sciences, University of Cape Town, Cape Town, South Africa

Abstract

Several commercially and ecologically important species in the southern Benguela have undergone southward and eastward shifts in their distributions over previous decades, most notably the small pelagic fish sardine Sardinops sagax and anchovy Engraulis encrasicolus. Understanding these changes and their implications is essential in implementing an ecosystem approach to fisheries in the southern Benguela and attempting to appreciate the potential impacts of future environmental change. To investigate possible impacts of these shifts at an ecosystem level, distribution maps for before (1985–1991), during (1997–2000) and after (2003–2008) the shift in small pelagic fish were constructed for 14 key species from catch and survey data, and used to calculate spatial indicators including proportion east and west of Cape Agulhas, relative overlap in biomass and area, index of diversity, connectivity. Potential interactions on the south and west coasts were also compared. For several species (redeye; chub mackerel; kingklip; chokka squid; yellowtail), previously unidentified increases in the proportion of biomass east of Cape Agulhas were shown to have occurred over the same period as that of small pelagic fish, although none to the same degree. On average, overlap with small pelagic fish increased over time and overall system connectivity was lowest in the intermediate period, possibly indicating a system under transition. Connectivity declined over time on the west coast while increasing on the east coast. Distributions of other species have changed over time, with the region east of Cape Agulhas becoming increasingly important in terms of potential trophic interaction. Variations in distribution of biomass and structural complexity affect the trophic structure and hence functioning of the system, and implications should be considered when attempting to identify the possible ecosystem impacts of current and future system-level change.

Introduction

An understanding of the functioning of an ecosystem is not achievable without knowledge of the interactions involved. One way in which this knowledge can be expanded is by observing any changes in the relative distributions of trophically linked species over time. Monitoring of past and current patterns becomes of even greater importance as climate change impacts become more apparent. Changes in migration and distribution ranges, as well as density distribution, are all known to be strongly influenced by environmental fluctuations [1] and have already been shown in multiple regions and species globally [26].

The southern Benguela is an eastern boundary current system extending from the Lüderitz upwelling cell at 26°S in the north around the coast of South Africa to East London at 28°E (Fig 1). The system is comprised of two physically and functionally different regions (i) the west coast, a classic wind-driven upwelling system, and (ii) the south coast which includes the Agulhas Bank and as a result, has characteristics of both upwelling and shelf systems [7,8]. Many species, particularly the small pelagic fish sardine (Sardinops sagax) and anchovy (Engraulis encrasicolus) have a migratory life history: adults spawn on the Agulhas bank and eggs and larvae are transported to the more productive west coast where recruitment and feeding takes place. Spawners then return to the south coast to complete the cycle. The system supports a relatively high fish biomass, which in turn supports a number of commercially important fisheries, the two largest in terms of catch and in economic value being the demersal trawl fishery, targeting deep and shallow water hake (Merluccius paradoxus and M. capensis), and the small pelagic purse-seine fishery.

thumbnail
Fig 1. The southern Benguela with west and south coasts divided at Cape Agulhas.

The shelf-edge is indicated by the 500m isobath [9].

https://doi.org/10.1371/journal.pone.0158734.g001

The concept of regime shifts, or a sudden and widespread change from one stable system state to another, has become more widely used in marine science as longer datasets become available for analysis [10,11]. When a number of physical and biological timeseries from the southern Benguela were analysed using methods to detect regime shifts, many displayed shifts between the mid-1990s and the early 2000s. Based on these findings, an ecosystem-level regime shift is thought to have occurred in the southern Benguela in the early 2000s [1214].

One of the most apparent changes is in the distribution of small pelagic fish (hereafter referring to sardine and anchovy unless otherwise specified), whose relative abundance declined on the west coast and increased on the south coast (east of Cape Agulhas, see Fig 1) in the late 1990s [15,16]. Small pelagic fish are trophically very important within the system and can exert both bottom-up and top-down effects on other species in what is known as wasp-waisted trophic control [17,18]. It therefore seems likely that this eastward increase in biomass, particularly when combined with the shifts in other variables that have since been identified, would have caused regional changes in the spatial interactions with trophically linked species and thus influenced the trophic functioning of the system. Observed trends in seabird abundance, distribution and breeding success since the mid-1990s have already been linked to the concurrent changes in small pelagic fish distribution [19]. Substantial declines in African penguins and gannet populations over the early 2000s, for example, have been related to the availability of small pelagics as prey, as has the increased abundance in Cape cormorants and swift terns on the south coast [1922].

Without knowledge of the interactions within a system, an effective application of an ecosystem approach to fisheries (EAF) is likely unattainable. Observing patterns in the relative distributions of trophically linked species, and monitoring change over time, can be a useful means of understanding system structure and any change it may undergo [2,4,23,24]. It can also allow spatial ecosystem indicators to be calculated. Indicators are an important component of applying an EAF and allow a large amount of data to be distilled into more easily communicable values, useful for monitoring temporal and spatial change, as well as evaluating management objectives [2528].

Distribution maps for key species in the southern Benguela were previously constructed for the 1980s and 1990s [29]. These maps aggregated over that entire period (1980s and 1990s) were then used to calculate a measure of spatial interaction between species and indicators of ecosystem state [24,30]. However, because ecosystem-level changes, including substantial distributional changes in key forage fish species (sardine and anchovy), have been identified as occurring over the late 1990s to early 2000s, previous findings are unlikely to be representative of the system beyond the mid-1990s.

Distribution maps were plotted for key species in the southern Benguela and compared over relevant time periods selected based on previously identified shifts, with a focus on changes in small pelagic fish distribution. Any changes in the distribution of other species or in spatial indicators over time can be used to identify possible changes in the level of interspecific interaction and increase our understanding of the changes at an ecosystem level in the context of future effects of climate change.

Methods

Time periods

The average spatial distribution during different periods was used as a means of examining change over time, based on the concept of regime shifts. System level shifts have been previously identified across biological and physical elements of the southern Benguela, the first in the early to mid-1990s and the second in the early 2000s [12,31,32], coinciding with a peak in small pelagic fish biomass. Based on these shifts, three time periods were chosen as representing different ecosystem states to be compared. Period 1 was chosen as 1985–1991, representing the system before the increased abundance and eastward shift in sardine and anchovy; Period 2 as 1997–2000, an intermediate/transition phase; and Period 3, during which both sardine and anchovy were predominantly found east of Cape Agulhas, as 2003–2008. To create a clear snapshot of the specified periods and the ecosystem states they represent, the years between these selected periods (1992–1996 and 2001–2002) were omitted.

Distributions

Data from either scientific surveys or commercial catch records were used to plot distribution maps. When distributions have been plotted previously [29], all available data sources were combined into a single distribution map. This requires a number of assumptions to be made regarding the ability of various data sources to represent a particular species’ distribution and the possible biases involved. For example, both commercial catch data and survey data are by definition more representative of the species targeted by that fishery or for which the survey is designed, and are likely to underrepresent other species. To combine the two data sources one would either need to assume equally accurate representation of the species in question by both, or decide on a weighting to the importance or representativeness of each. The aim in the current study was to detect change over time rather than to produce the most accurate distribution map possible. Therefore, in this context, rather than combine all information, a single data source was selected for each species as the most representative and best suited to reflecting temporal change for that species (see details below). Resultant distributions, while not assumed to be complete, are representative and appropriate given the nature of the study.

A total of 14 key species were examined: anchovy; sardine; round herring (Etrumeus whiteheadi); Cape hake (Merluccius capensis) and (M. paradoxus); horse mackerel (Trachurus trachurus capensis); chub mackerel (Scomber japonicus); kingklip (Genypterus capensis); chokka squid (Loligo vulgaris reynaudi); snoek (Thyrsites atun); silver kob (Argyrosomus inodorus); yellowfin tuna (Thunnus albacares); yellowtail (Seriola lalandi) and geelbek (Atractoscion aequidens). To choose the data source most representative of each, plots for each species for Period 3 (2003–2008) were constructed from all available data sources. Based on these maps and discussions with relevant experts for each species (scientists at the relevant South African academic and government institutions), a single data source was selected as most representative and most likely to illustrate change over time for each species (Table 1). To allow comparison between data sources that have different spatial scales, all data were joined to a 10’x10’ reference grid. All distributions were plotted using ArcGIS v 10.1.

thumbnail
Table 1. Data sources and units used to plot distributions for each species.

https://doi.org/10.1371/journal.pone.0158734.t001

Data sources and plotting

Survey data.

Annual pelagic hydro-acoustic surveys of principal small pelagic fish species (anchovy, sardine and round herring), comprised of a summer (November) spawner biomass and a winter (May) recruitment survey, have produced reliable data for the southern Benguela since 1984 [33]. Survey design and methods have been thoroughly described in [33] and [34]. Although initially designed to monitor anchovy, the abundance and biomass estimates from these surveys are now an essential input into the management process and serve as a basis for total allowable catch (TAC) recommendations for both sardine and anchovy fisheries [33,35]. The area covered by the pelagic surveys initially extended from either the Orange River (May survey) or Hondeklip Bay (November survey) on the west coast to Port Alfred on the south coast (Fig 1), but in recent years has extended further east. The time series has been revised to take into account changes in survey equipment and an increased understanding of possible sources of error in earlier estimates [35,36]. Pelagic survey data from both the May and November annual surveys combined were assumed as the most representative of the distribution of the small pelagic sardine, anchovy and round herring and were used to construct maps and in analysis for these species.

Positions of May and November pelagic survey intervals and their corresponding densities of sardine, anchovy and round herring, including zero values, were plotted and joined to the 10’x10’ reference grid. Average annual density per grid cell was calculated for each period. Although May and November surveys estimate recruit numbers and spawner biomass respectively, and thus reflect different stages of the life cycle and the associated differing distributions, the data from both surveys were combined. Because the aim of this investigation was to increase understanding of ecosystem function and how changes in distribution may have affected it, the point of interest lies in the availability of food to predators. Therefore both recruits and adult fish were considered. Data from both surveys were combined assuming an equal weighting. To avoid skewing the results due to the extension of sampling effort along the east coast during the more recent periods, an easterly limit on acceptable data points was set corresponding to the most easterly extent of surveys during Period 1.

Demersal biomass surveys, as described by [37], are also conducted around the coast of South Africa. Surveys are designed and used primarily for the estimation of hake biomass to inform the management of the hake-directed trawl fishery. Cruises have been conducted on the west coast (Orange River to Cape Agulhas) biannually from 1986–1990, after which the winter cruise was discontinued [32]. On the south coast (Cape Agulhas to Port Alfred) surveys were conducted annually in spring from 1986–1990, and subsequently in both spring and autumn [38]. Demersal biomass survey data from all available cruises combined were assumed to be the best descriptor of hake, horse mackerel, and chub mackerel distribution.

Data were excluded where trawl duration was less than 15 minutes or the trawl location north of the Orange River. Survey catches, including zero values, and minutes per trawl were converted to kg/hr. Catches and trawl duration were plotted and joined to the reference grid. The average annual catch in kg/hr per reference grid cell was then calculated for each period.

Commercial data.

A demersal trawl fishery has been operating off South Africa since the early 1900s [39], initially targeting Agulhas and West Coast sole (Austroglossus pectoralis) and (A. microlepis), but only targeting hake after the First World War. Catch data are provided in the form of catch weight (kilograms) and duration of trawl per 20’x20’ grid cell. Trawl data were taken as representative of kingklip, chokka squid and snoek distribution. To allow commercial inshore and offshore trawl data already assigned to a 20’x20’ grid to be displayed on the 10’x10’ reference grid, each reference grid cell was assigned one quarter of the catch of the 20’x20’ cell within which it was situated. The average annual catch in kg/hr per reference grid cell was then calculated.

South Africa’s line fishery is one of its oldest, with the boat-based fishery well-established by the mid-1800s. The fishery currently operates around the whole coastline on the continental shelf [40]. Catch and effort data are in the form of catch weight (kg) and catch days reported on a 5x5 minute grid. Line fishery data were used to produce maps for silver kob, yellowtail and geelbek. The 5x5 minute grid was joined to the reference grid and the catch and number of catch days per grid cell summed and used to calculate the average annual catch per fishing day per cell during each period. Data for kob in the linefishery do not distinguish between species, however their distribution is limited to the area between Cape Point and the Kei River mouth (28.37°E) on the south-eastern coast of South Africa (C. Attwood, Marine Research Institute, University of Cape Town). Any records east of this point were excluded from analyses, and all those remaining were assumed to consist of primarily Silver Kob.

An experimental longline fishery for tuna in the southern Benguela began in 1997, and has been operating as a commercial fishery since 2005 [41]. Boats report catch weight (kg) and the number of hooks per set, as well as the start and end positions of the line. Data from this fishery was used as the basis for the yellowfin tuna distribution. Start positions were plotted and the average annual catch in kg/1000 hooks for each reference grid cell during each period calculated.

For all data sources, any data points that fell obviously outside of possible sampling areas (i.e. on land/ far outside of known sampling area) were excluded from analyses. To allow for uncertainty in areas of low density/ catch rate, cells representing the lowest 5% of all distributions in each period were removed, thus the final map used in analysis only represents the core 95% of the distribution.

Analysis

Distribution maps were used to calculate the following indicators. Indicators i–iv relate to ecosystem state, whereas v is an indicator of pressure.

  1. The proportion of biomass found east and west of Cape Agulhas for each species in each period. To test for differences in the observed proportions in each period a beta regression model with a logit link function was fitted to the observed proportions, with period as an explanatory variable. This was done using the betareg package in R [42].
  2. The relative overlapping areas and biomass (ROA and ROB) between species where a and b are the trophically related species, j is the period and ∩ symbolises the intersect between the two [24,30]. Similarly, the proportion of total biomass overlapping was calculated (ROB). When calculating overlaps between species with distributions based on demersal survey data (hake, horse mackerel and chub mackerel) and small pelagic fish (based on pelagic survey data), only the November spawner biomass pelagic survey data were used for the overlaps due to the timing of the surveys: since 1990, no winter demersal surveys have been conducted. A more accurate overlap could therefore be obtained by excluding the winter pelagic survey from the analysis. For all other species, based on year round commercial data, overlaps were calculated using combined May and November pelagic survey data for small pelagic fish species.
    The averages of ROA and ROB between all species where a trophic relationship (predation or competition) exists (Table 2) were taken as a measure of overall ecosystem connectivity for that period [30]. Although ‘connectivity’ is used in trophic models to refer to the degree of trophic linkage within the system, the term can also be applied to a range of factors relating to states of physical or trophic and literal or potential connectedness, or degree of interaction within systems or between components of those systems, e.g. [4345]. Here, ‘connectivity’ refers to the average degree of physical overlap between species and thus the potential for interaction [30].
    Relationships were based on those identified by [29], initially derived from trophic relationships in the southern Benguela [46]. Those for geelbek and yellowtail were added based on dietary literature [4749] and discussion with DAFF scientists. All relationships involving yellowfin tuna were excluded from this calculation, as there are no data for this species during Period 1 (1985–1991). Differences in the degree of overlap over time were again tested using a beta regression model with a logit link function fitted to the observed proportions, with period as an explanatory variable [42].
  3. An index of spatial biodiversity (ISBj) was calculated for each period (j), based on mapped species only, and excluding yellowfin tuna as no data are available for this species for Period 1. Maps showing the number of species out of the possible 13 found in each cell for each period were also generated. ISBj was calculated as the average proportion of total possible species S found in any grid cell during period j according to: where n is the total number of cells with observations and s is the number of species in cell g [30]. The ISBj indicator has some drawbacks, namely that different patterns can give the same overall average result, and it may not be suited to intersystem comparisons without consideration of the implications of the species included when calculating the index [30]. It does, however, allow for an exploration of possible ecosystem level change over time, and thus is considered suitable for this application (see discussion). Differences in ISB between periods were again tested using a beta regression model with a logit link function fitted to the observed proportions and period as an explanatory variable [42].
  4. The proportion of effort east and west of Cape Agulhas was calculated for each data source, in units of survey intervals (pelagic survey), trawl minutes (demersal survey), fishing hours (demersal commercial), hooks set (longline) and catch days (linefishery). These proportions were averaged overall as well as for all commercial data and all survey data, for each period.
thumbnail
Table 2. Trophic relationships between species examined, represented qualitatively as either strong (P or C) or moderate (p or c) predation and competition respectively.

Weak relationships or occasional predation are not included. Updated from [24].

https://doi.org/10.1371/journal.pone.0158734.t002

Results

The data included in analyses, after points outside of the known sampling area had been cleaned from the datasets, varied across the three time periods from several hundred survey intervals for pelagic surveys, to over three million hook sets for longline data (Table 3).

thumbnail
Table 3. Extent of the data included in analyses, averaged per year for Period 1 (1985–1991), Period 2 (1997–2000) and Period 3 (2003–2008).

Units in which the effort data are originally recorded are shown.

https://doi.org/10.1371/journal.pone.0158734.t003

Analyses, including results from beta regression GLMs, confirm the increase in proportion of sardine biomass east of Cape Agulhas when compared with Period 1 (P2: p < 0.01, P3: p < 0.001, pseudo R2 = 0.86) as previously illustrated by [50] and [51]. A similar but lesser pattern is also present in a number of other species, including kingklip (P2: p < 0.001, P3: p < 0.001, pseudo R2 = 0.79), anchovy (P2: p < 0.05, P3: p < 0.001, pseudo R2 = 0.6), chokka (P2: p < 0.001, P3: p < 0.001, pseudo R2 = 0. 57), M. paradoxus (P2: p < 0.01, P3: p < 0.001, pseudo R2 = 0.41), redeye (P3: p < 0.05, pseudo R2 = 0.32)) and chub mackerel (P2: p < 0.01, P3: p <0.001, pseudo R2 = 0.26); (Figs 25). Geelbek showed the opposite trend; the proportion found east of Cape Agulhas was significantly lower in Period 3 compared to Period 1 (p < 0.05. pseudo R2 = 0.42) (Figs 4 and 5). Not all changes over time were linear however. M. paradoxus, chub mackerel and snoek for example increased on the south coast from Period 1 to Period 2, and then declined again in Period 3, although M. paradoxus and chub mackerel both remained at significantly higher levels on the south coast in Period 3 than they had been during Period 1.

thumbnail
Fig 2.

a—f: Distribution maps for sardine, anchovy, redeye, M. capensis, M. paradoxus and horse mackerel during Periods 1–3. The inset figure represents the proportion of biomass found east and west of Cape Agulhas in each period. Asterisks indicate where a significant change from P1 was detected on further statistical analysis, results below. * = p < 0.05; ** = p < 0.01; *** = p< 0.001.

https://doi.org/10.1371/journal.pone.0158734.g002

thumbnail
Fig 3.

a—f: Distribution maps for chub mackerel, kingklip, chokka squid, silver kob, snoek and yellowtail during Periods 1–3. The inset figure represents the proportion of biomass found east and west of Cape Agulhas in each period. Asterisks indicate where a significant change from P1 was detected on further statistical analysis, results below. * = p < 0.05; ** = p < 0.01; *** = p< 0.001.

https://doi.org/10.1371/journal.pone.0158734.g003

thumbnail
Fig 4.

a and b: Distribution maps for geelbek and yellowfin tuna during Periods 1–3. The inset figure represents the proportion of biomass found east and west of Cape Agulhas in each period. Note there are no data for yellowfin tuna in Period 1 1985–1991. Asterisks indicate where a significant change from P1 was detected on further statistical analysis, results below. * = p < 0.05; ** = p < 0.01; *** = p< 0.001.

https://doi.org/10.1371/journal.pone.0158734.g004

thumbnail
Fig 5. Estimated proportion of biomass east of Cape Agulhas in each period examined, for species where a significant change over time was found.

Proportions were estimated using a beta regression model fitted to the observed proportions, with period as an explanatory variable. The median, upper quartile, lower quartile and interquartile range are shown.

https://doi.org/10.1371/journal.pone.0158734.g005

Horse mackerel, and the linefish kob, snoek and yellowtail all showed similar proportions on the west and south coasts over all periods examined. Although effort in the tuna longline fishery expanded between periods 2 and 3, only catches made in the commonly sampled area were considered, and a far higher proportion of biomass was on the south coast than the west (77%) in the later period.

Similar patterns were illustrated by ROA and ROB (Figs 6 and 7) and there was increasing overlap between small pelagic fish as a group (including sardine, anchovy and redeye), and horse mackerel, chub mackerel, chokka, kob, yellowtail, geelbek, and to a lesser extent M. capensis, over all time periods (Fig 6). M. paradoxus, although displaying a slight increase in overlap with redeye, overlapped less with sardine and anchovy over time. On closer examination of results, this is a reflection of a decrease in overlap with small pelagic fish on the west coast (see the Appendix), where the majority of M. paradoxus biomass is found (Fig 2E). Overlap of both biomass and area with small pelagic fish increased over time on the south coast for M. paradoxus however, but this is masked by the west coast trend.

thumbnail
Fig 6. Average overlap in area and biomass of all other species with small pelagic fish species (sardine, anchovy and redeye) for the three periods examined.

https://doi.org/10.1371/journal.pone.0158734.g006

thumbnail
Fig 7.

a—p: Relative overlap in area and biomass between selected species and all other species mapped during each period. 1 unit on the y-axis is equivalent to 100% overlap. For all overlaps by species and coast, see supporting information (S1S5 Figs).

https://doi.org/10.1371/journal.pone.0158734.g007

In general overlaps were lower during Period 2, increasing again in Period 3 (Fig 7), except in the case of the two hake species, which displayed similar levels in Periods 2 and 3 with regard to area overlap, while biomass overlaps declined over all three periods. All overlaps between each species, separated for south and west coasts, are illustrated in the Appendix, displayed between each species and east and west of Cape Agulhas.

The pattern over time in system connectivity was similar based on area or biomass overlap (Fig 8), with a slight declining trend west of Cape Agulhas and increasing trend east of Cape Agulhas over time. When the total system is examined however no trend is detectable (Fig 8A and 8D), with modelled ROA and ROB also showing no significant difference from Period 1. Modelled predicted overlap in both area and biomass, east and west of Cape Agulhas was, however, significantly different in Periods 2 and 3 when compared with Period 1 (p < 0.0001), but in all cases models had low explanatory power (pseudo R2 = 0.023–0.029).

thumbnail
Fig 8.

Mean and standard deviation of system connectivity, or overlap in area (ROA) and biomass (ROB) between trophically related species (see Table 2 for relationships), shown for the total system (a and d) and regions east (c and f) and west (b and e) of Cape Agulhas.

https://doi.org/10.1371/journal.pone.0158734.g008

Increasing diversity on to the east when compared with the west coast is illustrated both by the maps of the combined species distributions (excluding yellowfin tuna due to lack of data for Period 1) (Fig 9), and by the Index of Spatial Biodiversity, ISB. The index was similar around the coast during Period 1, but over time decreased in the west and increased in the east (Fig 10). When data were used to fit beta regression GLMs, overall ISB (ISBsystem) was significantly lower for both Period 2 (p < 0.01) and 3 (p < 0.001) when compared with Period 1, although the explanatory power of the model was low (pseudo R2 = 0.01). The same results were true for the west coast as well (both p < 0.001, pseudo R2 = 0.04). ISB east of Cape Agulhas was not significantly different in Periods 2 and 3, but the model had an even lower explanatory power (pseudo R2 = 0.007).

thumbnail
Fig 9. The number of the 13 mapped species found per grid cell during each period.

Yellowfin tuna were excluded from this analysis due to a lack of data in period 1.

https://doi.org/10.1371/journal.pone.0158734.g009

thumbnail
Fig 10. Index of spatial biodiversity (ISB) east and west of Cape Agulhas, and for the system as a whole, during periods 1–3.

https://doi.org/10.1371/journal.pone.0158734.g010

As in [30], ISB (ISBall) was also calculated based on all common species identified by demersal surveys during each period (161 spp.) and compared to the above findings based on 13 species. As expected due to the large number of species included, ISBall was much lower than when calculated based on the 13 mapped species included, being lowest in Period 1 (8.93%), increasing to approximately 11% in Periods 2 and 3. In all periods ISBall was approximately 20% higher in the area east of Cape Agulhas than on the west coast.

Effort from all data sources increased on the east coast over time and decreased on the west (Fig 11). Due to the initial survey-driven bias towards the west coast, this has meant effort in the most recent period (3) is the most evenly distributed by coast of the periods examined. This difference is driven by increased survey effort rather than increasing commercial effort on the east coast (for example proportion of demersal commercial effort on the east coast in Period 1: Period 3 is 56:59%, while that for demersal survey is 33:62%).

thumbnail
Fig 11. Proportion of effort east and west of Cape Agulhas averaged for all data sources used to construct distributions.

Effort is shown both combined from all sources, and for commercial and survey data separately, during periods 1–3 (P1 –P3).

https://doi.org/10.1371/journal.pone.0158734.g011

Discussion

The pattern of eastward movement highlighted in this study reflects what has already been illustrated for both sardine [50,51] and anchovy [15,31], but is not consistently present in all higher trophic level species as one might expect. Although both hake species are strongly trophically linked to various other species that display increased abundance on the south coast over time (Table 2), e.g. small pelagic fish, horse mackerel and chokka, neither has followed the same pattern. The proportion of M. capensis on the south coast did not increase significantly in Period 2, and proceeded to decline in Period 3 while the proportion of M. paradoxus on the other hand was significantly higher on the south coast during Periods 2 and 3. The relationship to changes in sardine was not linear however, with a decline in prevalence on the south coast between Periods 2 and 3. If dietary prevalence of small pelagics or chokka is taken into account, kob, geelbek and yellowtail would also have been expected to exhibit a parallel increase on the south coast. On the contrary, geelbek decreased in abundance in the east during Period 3. Yellowtail did display an increasing but non-significant trend. While kingklip did increase proportionally east of Cape Agulhas over time, it is possible that this is related to successful implementation of management actions in the region, such as the ‘kingklip box’ seasonal closure near Port Elizabeth [52], rather than a shift in distribution.

Although various non-climate factors, such as the recovery of previously exploited stocks in the east, should be considered as possible explanations of the increase in abundance in this area, in this case they do not provide the most likely cause. Overall, fishing effort on the south coast has not declined over the periods examined (Fig 11). The main species group displaying this increase in easterly abundance–small pelagic fish–has been both relatively extensively fished and well monitored by biannual surveys, as described in the methods section. Neither survey data (used to construct these distributions) nor catch data for those species reflect any heavy exploitation and subsequent recovery in the east. They do however support the hypothesis that sustained effort on the west coast concurrent with declines in biomass in that region may have exacerbated the pattern [50]. In the line fishery, heavy depletion of silver kob and geelbek stocks along with other targeted species by the late 1990s [40] did lead to a declaration of a state of emergency in the fishery, and the development of the linefish management protocol [53] (an approach using biological reference levels for management planning). Although the resultant reduction in effort by approximately 70% in 2000 was substantial, and there is an indication that silver kob catch per unit effort may be increasing [54], this is not reflected in results presented here and the status of most line fish species remain collapsed or unknown [9]. Pressure from the line and inshore trawl fisheries remains relatively high given the low abundance of these species [41,55], and is still likely to effect the ability of line fish species to respond to increases in prey abundance. Changes in demersal assemblages over the period examined are also not thought to be linked only to changes in fishing effort [32]. As a result, and given the widespread nature of the shifts observed here as well as described elsewhere [9,14], relaxation of effort leading to a regional recovery is not a probable explanation and system-level climate forcing seems likely to have played a major role.

Top predators were not included in this study, however, as previously discussed, effects of this change in small pelagic fish distribution on seabirds have been well documented [19,21,56]. Although the prevalence of sardine in the diet of Cape fur seals, opportunistic top predators, has been shown to reflect local availability, no concurrent change in seal populations with changes in small pelagic fish abundance or distribution has been recorded in South Africa [57]. This stable distribution is largely attributed to the lack of additional suitable breeding habitat limiting any potential expansion of the population [58, 59]. The decline and subsequent expansion of populations at the northern extent of their range in Namibia and Angola during the 1990s and early 2000s have been linked to changes in abundance and distribution of prey species and supports this hypothesis, illustrating a direct effect on seals due to changes in prey availability if space limitation is not a factor [58]. An additional factor dampening possible seal distributional response is the opportunistic nature and omnivory of seals [57,59].

A number of species show peaks in abundance east of Cape Agulhas during Period 2, along with the documented shift in anchovy spawner biomass and eggs east of Cape Agulhas [15]. This abrupt change in anchovy distribution has been linked to environmental changes, possibly induced by increased coastal upwelling east of Cape Agulhas and the resultant improvement of feeding conditions for anchovy spawners relative to those on the west coast [31]. This may illustrate an environmentally-mediated, bottom-up mechanism of trophic control whereby changes in primary production affect higher trophic level abundance and distribution [60]. The same mechanism may have influenced the distribution of other species shown here to have increased on the south coast during Period 2, (e.g. chub mackerel). Because the upwelling in question is restricted to coastal regions, however, it is unlikely to have directly effected the Period 2 increase in snoek east of Cape Agulhas, where records are largely further offshore.

Although plotting only 95% of the distribution has been suggested as too low a threshold [24,29], in this context where a map based on a single data source was used to explore change over time, it was deemed a reasonable if conservative measure, allowing for meaningful interpretation of overlaps. As expected, where both species in the pair being examined for overlap exhibited an increase in proportion on the south coast, degree of overlap increased. In many cases, the trophic functioning of the system is likely to have been affected by changes in overlap since trophically related species (see Table 2) are not sharing the same proportions of their distributions as they were previously. Chokka squid, for example, is more prevalent on the south coast. The degree of overlap with a species increasingly found on the south coast, e.g. sardine, will thus have increased. As both predation and competition relationships exist between the two species (Table 2), trophic interactions within the region must be altered. However, that species pairs exhibiting increased overlap are not always coupled predator and prey species seems to imply an outside driving force, or drivers, other than trophic interactions considered here, such as environmentally favourable conditions to the east as discussed by [31].

System connectivity (Fig 8) can give an indication of the ecosystems resilience and ability to withstand change. While the connectivity of the system as a whole does not appear to have changed (Fig 8A and 8D), the decline in connectivity in the west (Fig 8B and 8E) and the increase in connectivity in the east (Fig 8C and 8F) may depict the rise in trophic importance of the south coast since Period 1 as more potential trophic interactions are located there than on the west coast. Connectivity in this study has been based on the interaction among only the 14 species examined and known to be trophically linked (Table 2) and is thus not truly a reflection of the whole ecosystem. It does however, include the most prevalent, commercially and trophically important species, allowing for comparison between time periods.

Despite its drawbacks, the ISBj indicator is an appropriate measure to have used in this analysis where the overarching, within-system change is of interest. ISBj does seem to reflect the increasing food web complexity on the south coast that has been suggested by increasing connectivity on that coast. The inclusion of a large number of less prevalent species in this calculation of ISBall (not pictured in Fig 10) gives an overview of system state from a different perspective, but also explains the lack of agreement between the outcomes of the two ISB calculations (a dip in ISBj in Period 2, and an increase in ISBall from Period 1 to Periods 2 and 3).

The different data sources used have their own advantages and disadvantages. Survey data are by definition more suited to produce an accurate index of biomass, although it should be kept in mind that demersal surveys have been designed around hake and the pelagic surveys around anchovy and sardine, thus are not as representative for other species. For example, both redeye and chub mackerel are likely to have ranges extending further offshore than is captured by surveys., and horse mackerel have been recorded by acoustic survey over the shelf-break area when demersal survey methods failed to detect them there [61]. Commercial data on the other hand, while providing far greater sampling effort than surveys can, are unable to provide unbiased data as the effort is not random. Commercially viable concentrations of target species are actively sought out, for example data for the pelagic fishery will reflect only those fish that are accessible in terms of port and processing facilities. Commercial data can potentially provide more accurate information and better coverage for species that are not the object of survey data collection, but the limitations must be kept in mind when assessing results. Unfortunately, any single data source is very unlikely to represent the full extent of a species distribution. The method can, however, still provide a useful means of comparing change over time, as has been done here, rather than reporting absolute distributions.

Results could be improved if the disparate distributions and trophic roles of juveniles and adults of trophically important species could be taken into account. These data are only available at an ecologically satisfactory level of detail for sardine and anchovy (although data for redeye are also available, they do not cover the entire distribution). The decision to combine the May and November recruit and spawner biomass survey data for the small pelagic fish species was based on the goal of investigating ecosystem structure and functioning, and as such the need to understand the availability and distribution of small pelagic as potential food represented by all life-stages of a predatory species within the ecosystem. However, the differences between the recruit and spawner distributions identified by each survey should be kept in mind. For example, the eastward distributional shift in anchovy [15] is evident only in the spawner biomass, while recruits remain almost entirely on the traditional west coast nursery grounds (see supporting information S6 Fig), with greater implications for the seasonality of prey availability to predators. On the other hand, sardine and redeye display similar trends in both recruit and spawner biomass distributions. Both increasingly are found east of Cape Agulhas over time, with redeye recruits actually showing a more pronounced distributional trend towards the east than displayed by the redeye spawner biomass data (37% vs 18% increase in proportion found east of Cape Agulhas from Period 1 to Period 3). Thus, while sardine and redeye seem to have undergone a distributional shift, anchovy have experienced a shift in spawning area, and this must be taken into consideration when interpreting results.

In many cases distributions and patterns of interaction have changed over time. These changes and changes reported in top predator species such as seabirds support the idea that an understanding of trophic relationships is important if potential system-wide changes are to be understood or anticipated. Although previous studies based on ecosystem models have described changes to overall system functioning [62,63], patterns of change identified were not as strong when considering the average state of the ecosystem as a whole. Regional differences between the areas east and west of Cape Agulhas were evident, however, with notable shifts in species distributions and the potential interactions between species (predators, prey and competition for common prey). During Period 1, a number of species displayed a far greater prevalence on the west coast compared to their Period 3 distributions. Connectivity was highest on the west coast during Period 1, and to the east in Period 3, with an overall dip during the intermediate Period 2, a pattern largely echoed by the index of spatial biodiversity. Although unfortunately beyond the scope of this study, future expansion of this work to investigate links with physical variables would be an important step in understanding system function. Spatio-temporal variations in biomass and structural complexity clearly affect the structure and functioning of the system, and further understanding of the implications of changes such as those described here is important when attempting to appreciate the possible implications of current and future system-level change for the ecosystem and management of the fisheries it supports.

Supporting Information

S1 Fig. Overlap in area (ROA) and biomass (ROB) between sardine, anchovy and redeye and all other species, east and west of Cape Agulhas.

https://doi.org/10.1371/journal.pone.0158734.s001

(TIF)

S2 Fig. Overlap in area (ROA) and biomass (ROB) between M. capensis, M. paradoxus and horse mackerel, and all other species, east and west of Cape Agulhas.

https://doi.org/10.1371/journal.pone.0158734.s002

(TIF)

S3 Fig. Overlap in area (ROA) and biomass (ROB) between chub mackerel, kingklip and chokka squid, and all other species, east and west of Cape Agulhas.

https://doi.org/10.1371/journal.pone.0158734.s003

(TIF)

S4 Fig. Overlap in area (ROA) and biomass (ROB) between Silver kob, snoek and yellowfin tuna, and all other species, east and west of Cape Agulhas.

https://doi.org/10.1371/journal.pone.0158734.s004

(TIF)

S5 Fig. Overlap in area (ROA) and biomass (ROB) between yellowtail and geelbek, and all other species, east and west of Cape Agulhas.

https://doi.org/10.1371/journal.pone.0158734.s005

(TIF)

S6 Fig. Average distribution of sardine, anchovy and redeye recruits (May) and spawner biomass (Nov) during each period, based on pelagic survey data.

https://doi.org/10.1371/journal.pone.0158734.s006

(TIF)

Acknowledgments

We thank staff at the South African governmental Departments of Agriculture, Forestry and Fisheries (DAFF), and Environmental Affairs (DEA) for generously assisting with both data and expert opinions when needed: Dr Carl van der Lingen, Janet Coetzee, Dr Tracey Fairweather, Dagmar Merkle, Jan van der Westhuizen, Dr Sven Kerwath, Dr Rob Crawford, Sobahle Somhlaba, Charlene Da Silva, as well as Prof. Colin Attwood (Ma-RE, UCT) and Sarah Wilkinson (CapFish). We also thank the anonymous reviewers for their helpful comments.

Author Contributions

Conceived and designed the experiments: KEW LH AJ LS. Performed the experiments: KEW. Analyzed the data: KEW. Wrote the paper: KEW. Formulation of approach and interpretation of data: KEW LH AJ LJS. Critical assessment and final approval of version to be published: KEW LH AJ LJS.

References

  1. 1. Lehodey P, Alheit J, Barange M, Baumgartner T, Beaugrand G, Drinkwater KF, et al. Climate Variability, Fish, and Fisheries. J Clim. 2006;19: 5009–5030.
  2. 2. Murawski SA. Climate Change and Marine Fish Distributions: Forecasting from Historical Analogy. Trans Am Fsiheries Soc. 1993;122: 647–658.
  3. 3. Perry AL, Low PJ, Ellis JR, Reynolds JD. Climate Change and Distribution Shifts in Marine Fishes. Science (80-). 2005;308: 1912–1915. pmid:15890845
  4. 4. Hiddink JG, ter Hofstede R. Climate induced increases in species richness of marine fishes. Glob Chang Biol. 2008;14: 453–460.
  5. 5. Rijnsdorp AD, Peck MA, Engelhard GH, Mo C, Pinnegar JK. Resolving the effect of climate change on fish populations. 2009;
  6. 6. Simpson SD, Blanchard JL, Genner M. Impacts of climate change on fish. 2013; 113–124.
  7. 7. Hutchings L, van der Lingen CD, Shannon LJ, Crawford RJM, Verheye HMS, Bartholomae CH, et al. The Benguela Current: An ecosystem of four components. Prog Oceanogr. 2009;83: 15–32.
  8. 8. Shannon L V. The Benguela ecosystem: Part I. Evolution of the Benguela, physical features and processes. Oceanogr Mar Biol An Annu Rev. 1985;23: 105–182.
  9. 9. Blamey LK, Shannon LJ, Bolton JJ, Crawford RJM, Dufois F, Evers-King H, et al. Ecosystem change in the southern Benguela and the underlying processes. J Mar Syst. 2014;144: 9–29.
  10. 10. de Young B, Harris R, Alheit J, Beaugrand G, Mantua N, Shannon LJ. Detecting regime shifts in the ocean: Data considerations. Prog Oceanogr. 2004;60: 143–164.
  11. 11. Jarre A, Moloney CL, Shannon LJ, Fréon P, van der Lingen CD, Verheye HMS, et al. Developing a Basis for Detecting and Predicting Long-Term Ecosystem Changes. In: Shannon L V, Hempel G, Malanotte Rizzoli P, Moloney CL, Woods J, editors. Benguela: Predicting a Large Marine Ecosystem. Amsterdam: Elsevier; 2006. pp. 239–273.
  12. 12. Howard JAE, Jarre A, Clark AE, Moloney CL. Application of the sequential t-test algorithm for analysing regime shifts to the southern Benguela ecosystem. African J Mar Sci. 2007;29: 437–451.
  13. 13. Atkinson LJ, Jarre A, Shannon LJ, Field JG. Evidence for Shifts in Demersal Fish Populations on the West Coast of South Africa: 1986 to 2009. In: Kruse GH, Browman HI, Cochrane KL, Evans D, Jamieson GS, Livingston PA, et al., editors. Alaska Sea Grant, University of Alaska Fairbanks; 2012. p. 20.
  14. 14. Blamey LK, Howard JAE, Agenbag JJ, Jarre A. Regime-shifts in the southern Benguela shelf and inshore region. Prog Oceanogr. Elsevier Ltd; 2012;106: 80–95.
  15. 15. van der Lingen CD, Coetzee JC, Hutchings L. Temporal shifts in the spatial distribution of anchovy spawners and their eggs in the Southern Benguela: Implications for recruitment. GLOBEC Rep. 2002; 46–48.
  16. 16. Fairweather TP, van der Lingen CD, Booth AJ, Drapeau L, Westhuizen JJ Van Der. Indicators of sustainable fishing for South African sardine Sardinops sagax and anchovy Engraulis encrasicolus. African J Mar Sci. 2006;28: 661–680.
  17. 17. Cury PM, Bakun A, Crawford RJM, Jarre A, Quin RA, Shannon LJ, et al. Small pelagics in upwelling systems: patterns of interaction and structural changes in “‘ wasp-waist ‘” ecosystems. ICES J Mar Sci. 2000;57: 603–618.
  18. 18. Shannon LJ, Coll M, Neira S, Cury PM, Roux J-P. The role of small pelagic fish in the ecosystem. 2007;
  19. 19. Crawford RJM, Sabarros PS, Fairweather TP, Underhill LG, Wolfaardt AC. Implications for seabirds off South Africa of a long-term change in the distribution of sardine. African J Mar Sci. 2008;30: 177–184.
  20. 20. Crawford RJM, Dyer BM, Kemper J, Simmons RE, Upfold L. Trends in numbers of Cape Cormorants (Phalacrocorax capensis) over a 50-year period, 1956–57 to 2006–07. Emu. 2007;107: 253–261.
  21. 21. Crawford RJM. A recent increase of swift terns Thalasseus bergii off South Africa–The possible influence of an altered abundance and distribution of prey. Prog Oceanogr. 2009;83: 398–403.
  22. 22. Crawford RJM, Altwegg R, Barham PJ, Durant JM, Dyer BM, Geldenhuys D, et al. Collapse of South Africa ‘ s penguins in the early 21st century. African J Mar Sci. 2011;33: 139–156.
  23. 23. Beare DJ, Burns F, Greig A, Jones EG, Peach K, Kienzle M, et al. Long-term increases in prevalence of North Sea fishes having southern biogeographic affinities. Mar Ecol Prog Ser. 2004;284: 269–278.
  24. 24. Drapeau L, Pecquerie L, Fréon P, Shannon LJ. Quantification and representation of potential spatial interactions in the Southrn Benguela ecosystem. African J Mar Sci. 2004;26: 141–159.
  25. 25. Garcia SM, Staples DJ. Sustainability reference systems and indicators for responsible marine capture fisheries: a review of concepts and elements for a set of guidelines. Mar Freshw Res. 2000;51: 385–426.
  26. 26. Cury PM, Christensen V. Quantitative Ecosystem Indicators for Fisheries Management. ICES J Mar Sci. 2005;62: 307–310.
  27. 27. Babcock EA, Pikitch EK, McAllister K, Apostolaki P, Santora C. A perspective on the use of spatialized indicators for ecosystem-based fishery management through spatial zoning. ICES J Mar Sci. 2005;62: 469–476.
  28. 28. Shin Y-J, Bundy A, Shannon LJ, Blanchard JL, Chuenpagdee R, Coll M, et al. Global in scope and regionally rich: an IndiSeas workshop helps shape the future of marine ecosystem indicators. Rev Fish Biol Fish. 2012;22: 835–845.
  29. 29. Pecquerie L, Drapeau L, Fréon P, Coetzee JC, Leslie RW, Griffiths MH. Distribution patterns of key fish species of the southern benguela ecosystem: an approach combining fishery-dependent and fishery-independent data. African J Mar Sci. 2004;26: 115–139.
  30. 30. Fréon P, Drapeau L, David JHM, Moreno A, Leslie RW, Oosthuizen WH, et al. Spatialized ecosystem indicators in the southern Benguela. ICES J Mar Sci. 2005;62: 459–468.
  31. 31. Roy C, van der Lingen CD, Coetzee JC, Lutjeharms JRE. Abrupt environmental shift associated with changes in the distribution of Cape anchovy Engraulis encrasicolus spawners in the southern Benguela. African J Mar Sci. 2007;29: 309–319.
  32. 32. Atkinson LJ, Leslie RW, Field JG, Jarre A. Changes in demersal fish assemblages on the west coast of South Africa, 1986–2009. African J Mar Sci. 2011;33: 157–170.
  33. 33. Hampton I. The role of acoustic surveys in the assessment of pelagic fish resources on the South African continental shelf. South African J Mar Sci. 1992;12: 1031–1050.
  34. 34. Barange M, Hampton I, Roel BA. Trends in the abundance and distribution of anchovy and sardine on the South African continental shelf in the 1990s, deduced from acoustic surveys. South African J Mar Sci. 1999;21: 367–391.
  35. 35. de Moor CL, Butterworth DS, Coetzee JC. Revised estimates of abundance of South African sardine and anchovy from acoustic surveys adjusting for echosounder saturation in earlier surveys and attenuation effects for sardine. African J Mar Sci. 2008;30: 219–232.
  36. 36. Coetzee JC, Merkle D, de Moor CL, Twatwa NM, Barange M, Butterworth DS. Refined estimates of South African pelagic fish biomass from hydro-acoustic surveys: quantifying the effects of target strength, signal attenuation and receiver saturation. African J Mar Sci. 2008;30: 205–217.
  37. 37. Badenhorst A, Smale MJ. The distribution and abundance of seven commercial trawlfish from the Cape south coast of South Africa, 1986–1990. South African J Mar Sci. 1991;11: 377–393.
  38. 38. Yemane D, Field JG, Leslie RW. Indicators of change in the size structure of fish communities: A case study from the south coast of South Africa. Fish Res. 2008;93: 163–172.
  39. 39. Lees R. The Story of the Fishing Industry in Southern Africa–and the Men who Made It. Cape Town: Purnell; 1969.
  40. 40. Griffiths MH. Long-term trends in catch and effort of commercial linefish off South Africa’s Cape Province: snapshots of the 20th century. South African J Mar Sci. 2000;22: 81–110.
  41. 41. DAFF. Status of the South African Marine Fishery Resources 2012. Cape Town; 2012.
  42. 42. Cribari-neto F. Beta Regression in R. J Stat Softw. 2010;34: 1–24.
  43. 43. With KA, Gardner RH, Turner MG. Landscape Connectivity and Population Distributions in Heterogeneous Environments. Oikos. 1997;78: 151–169.
  44. 44. Link JS. Does food web theory work for marine ecosystems? Mar Ecol Prog Ser. 2002;230: 1–9.
  45. 45. Cadenasso ML, Pickett ST a., Grove JM. Dimensions of ecosystem complexity: Heterogeneity, connectivity, and history. Ecol Complex. 2006;3: 1–12.
  46. 46. Shannon LJ, Moloney CL, Jarre A, Field JG. Trophic flows in the southern Benguela during the 1980s and 1990s. J Mar Syst. 2003;39: 83–116.
  47. 47. Nepgen CS de V. Diet of predatory and reef fish in False Bay and possible effects of pelagic purse-seining on their food supply. Fish Bull. 1982;16: 75–93.
  48. 48. Smale MJ. The feeding biology of four predatory reef fishes off the South-East Cape coast, South Africa. South African J Zool. 1986;21: 111–130.
  49. 49. Griffiths MH, Hecht T. On the life-history of Atructoscion aeyuidens, a migratory sciaenid off the east coast of southern Africa. J Fish Biol. 1995;47: 962–985.
  50. 50. Coetzee JC, van der Lingen CD, Hutchings L, Fairweather TP. Has the fishery contributed to a major shift in the distribution of South African sardine? ICES J Mar Sci. 2008;65: 1676–1688.
  51. 51. van der Lingen CD, Coetzee JC, Demarcq H, Drapeau L, Fairweather TP, Hutchings L. An eastward shift in the distribution of southern Benguela sardine. GLOBEC Newsl. 2005;11: 17–22.
  52. 52. Field JG, Attwood CG, Jarre A, Sink K, Atkinson LJ, Petersen S. Cooperation between scientists, NGOs and industry in support of sustainable fisheries: The South African hake Merluccius spp. trawl fishery experiencea. J Fish Biol. 2013;83: 1019–1034. pmid:24090560
  53. 53. Griffiths MH. The application of per-recruit models to Argyrosomus inodorus, an important South African sciaenid fish. Fish Res. 1997;30: 103–115.
  54. 54. Winker H, Kerwath S, Attwood CG. Comparison of two approaches to standardize catch-per-unit-effort for targeting behaviour in a multispecies hand-line fishery. Fish Res. Elsevier B.V.; 2013;139: 118–131.
  55. 55. Winker H, Booth AJ, Kerwath S, Attwood CG. Fitting Bayesian state-space biomass dynamics models to standardized CPUE for carpenter and silver kob stocks. Internal Working Document. MARAM IWS/Nov12/LF/P1. 2012.
  56. 56. Crawford RJM, Underhill LG, Coetzee JC, Fairweather TP, Shannon LJ, Wolfaardt AC. Influences of the abundance and distribution of prey on African penguins Spheniscus demersus off western South Africa. African J Mar Sci. 2008;30: 167–175.
  57. 57. Kirkman SP, Oosthuizen WH, Meÿer MA, Kotze PGH, Boucher M, Underhill LG. Chapter 7 Ecological responses of Cape Fur Seals in South Africa to temporal shifts in pelagic prey availability. In: Kirkman SP, editor. Final report of BCLME (Benguela Current Large Marine Ecosystem) Project on Top Predators as Biological Indicators of Ecosystem Change in the BCLME. Cape Town: Avian Demography Unit; 2006.
  58. 58. Kirkman SP, Yemane D, Oosthuizen WH, Meÿer MA, Kotze PGH, Skrypzeck H, et al. Spatio-temporal shifts of the dynamic Cape fur seal population in southern Africa, based on aerial censuses (1972–2009). Mar Mammal Sci. 2013;29: 497–524.
  59. 59. David JHM. Diet of the South African fur seal (1974–1985) and an assessment of competition with fisheries in southern Africa. South African J Mar Sci. 1987;5: 693–713.
  60. 60. Cury PM, Shannon LJ, Shin Y-J. The functioning of marine ecosystems. Responsible Fisheries in the Marine Ecosystem. 2003. pp. 103–123.
  61. 61. Barange M, Pillar SC, Hampton I. Distribution patterns, stock size and life-history strategies of Cape horse mackerel Trachurus trachurus capensis, based on bottom trawl and acoustic surveys. South African J Mar Sci. 1998;19: 433–447.
  62. 62. Watermeyer KE, Shannon LJ, Griffiths CL. Changes in the trophic structure of the southern Benguela before and after the onset of industrial fishing. African J Mar Sci. 2008;30: 351–382.
  63. 63. Osman W. Trophic model-generated indicators of the southern Benguela ecosystem for communicating with fisheries managers. Production. University of Cape Town. 2010.