Long-Term Change of Demersal Fish Assemblages on the Inshore Agulhas Bank Between 1904 and 2015

Without baseline data from near pristine assemblages, measures of ecosystem change may be significantly underestimated. A unique historical dataset provided an opportunity to investigate long-term change in demersal fish assemblages of South Africa’s inshore trawl grounds. Three sites surveyed over a period from 1903 to 1904 were re-surveyed in 2015 using replicated historical gear and methods. Catch composition was contrasted between historical and modern periods using unconstrained ordination, permutational multivariate analysis of variance, permutational tests of the homogeneity of multivariate group dispersions and similarity percentage analyses. After 111 years, the re-survey revealed a drastically transformed demersal assemblage, with the period effect explaining nearly half of the variance among samples. Historical catches were dominated by kob (Argyrosomus spp.), panga (Pterogymnus laniarius) and east coast sole (Austroglossus pectoralis), jointly contributing 70–84% of the catch. The same taxa made up a minor component (1.5–5.5%) of modern assemblages. Instead, the re-survey catches consisted predominantly of gurnards (Chelidonichthys spp.), Cape horse mackerel (Trachurus capensis), spiny dogfish (Squalus spp.), shallow-water hake (Merluccius capensis), and white sea catfish (Galeichthys feliceps), with their summed contribution rising to 85% from the historical 3%. These results suggest that a century of trawling may have altered benthic habitats, indirectly contributing to changes in the fish community. Historical assemblages included a substantial proportion of taxa that associate with reef habitats, whereas the re-survey assemblages were characterized by species that inhabit unconsolidated sediments or both reef and non-reef habitats. The unique historical context and data, comparable gear and methods and long temporal period revealed striking baseline changes that may be overlooked in most fisheries. Reconstructing this important historical context improves our ability to assess, interpret and manage changing marine ecosystems.


INTRODUCTION
Effective policy advice on the management of biodiversity, fisheries and the ecosystems that support them requires information about the scale and drivers of past changes in the ocean (Schwerdtner Máñez et al., 2014;Engelhard et al., 2016). Without knowledge of previous change and baseline reference points, recent variability within a fishery or ecosystem may provide little context for interpretation (Thurstan et al., 2010;Moloney et al., 2013), which hinders informed decision making.
Anthropogenic impacts such as fishing and climate change have altered marine ecosystems (Jackson et al., 2001;Lotze et al., 2006;Brander, 2010;Cheung et al., 2013;Moloney et al., 2013). Most of our knowledge of these systems is based on research and monitoring initiated long after such impacts started altering them (Lotze and Worm, 2009). For example, some of the most intense trawl fishing pressures in South Africa occurred between the late 1950s and 1977, when South Africa declared its exclusive fishing zone (Griffiths et al., 2004;Sink et al., 2012b). Yet investigations into periods prior to modern time-series of research trawl surveys (1983 to present), commercial catches (1978 to present) and fishery observer data (1995 to present), are few and typically hampered by the availability of comparable data. As a result, knowledge of near pristine ecosystems and the reference points that describe them are broadly lacking, both in South Africa and across most of the globe.
Valuable insight has been gained from research into longterm assemblage changes in demersal catches. By comparing trawl surveys from 1909 and 1910 to the 1980s, Novaglio et al. (2020) showed that substantial population declines and change in fish community structure occurred during the early decades of commercial trawling in south-east Australia. Periods of warmer water have been associated with greater abundances of certain taxa, specifically smaller, faster-growing species (Genner et al., 2004(Genner et al., , 2010ter Hofstede and Rijnsdorp, 2011). Investigating demersal assemblage changes in the North Sea over three decades, Engelhard et al. (2011) concluded that both climatic changes and fishing impacts explained observed patterns, which included higher abundances of warmwater, smaller, lower-trophic-level taxa and declines of mostly larger, cooler-water, higher-trophic-level species. Fishing impacts are commonly blamed for decreased abundances observed in demersal populations or assemblages (Silvestre et al., 1986;Rijnsdorp et al., 1996;Jackson et al., 2001;Kongprom et al., 2003;Ferretti et al., 2013;Novaglio et al., 2020), as well as declines in larger species or larger sizes within a population (Rijnsdorp et al., 1996;Greenstreet et al., 1999;Rogers and Ellis, 2000;Genner et al., 2010). Decreased diversity and greater species dominance in trawl-caught assemblages have also been attributed to the effects of fishing (Greenstreet and Hall, 1996;Rijnsdorp et al., 1996;Greenstreet et al., 1999). Heath and Speirs (2012) investigated fishing impacts on the demersal community following the resumption of trawling after it had been banned for 73 years from the Firth of Clyde in Scotland. They documented substantial change in the demersal assemblage in the decades that followed, including reduced evenness, marked reductions in fish sizes and an increase in small-sized fish, notably whiting (Merlangius merlangus).

Development and Change in South African Demersal Communities
Pre-colonial exploitation of marine resources in South Africa has a long and important history (Parkington, 2001), but it was focused on intertidal reaches (Griffiths et al., 2004). Relative to most temperate shelf regions around Europe and North America, industrialized fishing was late to reach South African ecosystems, similar to the developments in south-east Australia (Novaglio et al., 2018). Commercial fishing of finfish began near Cape Town in the late 1700s and spread with the growth of colonial settlements along the coasts. Fish were netted from the shore, but greater volumes were taken by handline fishers from small boats that caught most of their fish close to ports (Lees, 1969). By 1895 the line fishery had spread east to Port Elizabeth, and reported impressive catches across a broad spectrum of coastal fishes (Department of Agriculture, 1897).
Following substantial catches made by the government research vessel, SS Pieter Faure, industrial fishing arrived, with the first steam trawler registered at Simon's Town (Cape Town) in 1898 (Department of Agriculture, 1899). The inshore Agulhas Bank sole grounds were initially the most important trawl fishing grounds (Scott, 1949). The industry expanded with developing transport networks and markets. In 1922 there were 17 steam trawlers registered in the Cape Province (Lees, 1969). At this time, large hake resources (Merluccius capensis and M. paradoxus) were discovered in deeper waters on the west coast (Scott, 1949;Payne and Punt, 1995), after which their catches dominated regional trawl landings. Catches of hake have remained the most valuable and largest part of South Africa's trawl catches to the present.
Growth of the trawl industry accelerated after the Second World War, with hake landings growing from 50,000 tons in 1950 to 160,000 tons by 1960 (Payne and Punt, 1995). Foreign fleets arrived in the region in the 1960s, rapidly escalating fishing effort (Payne and Punt, 1995;Durholtz et al., 2015). Hake landings reached a peak in the early 1970s (∼300,000 tons), after which they fell due to over-exploitation (Payne and Punt, 1995). In 1972 the International Commission for the Southeast Atlantic Fisheries was formed to regulate regional fishing pressures and in 1977 the South African government proclaimed its exclusive fishing zone that progressively removed international effort (Durholtz et al., 2015). Hake populations started to recover in the 1980s and their catches have subsequently remained relatively stable at ∼150,000 tons per year (Griffiths et al., 2004;DAFF, 2015). An overview of the development, gear and vessel technology, management strategies and fleet structure of South Africa's trawl fishery are provided in Sink et al. (2012b) and Durholtz et al. (2015).
The line fishery also grew substantially during the 20th century and its efficiency greatly increased with technological progress. Although the inshore trawlers are confined to softsediments and line fishers exploit mostly reef habitat, conflict arose early on and there remains substantial overlap in fishing grounds and in catch composition. Line fishers first experienced catch restrictions as late as 1985, in the way of bag and size limits and a cap on commercial effort. These provisions were insufficient and likely too late. Griffiths (2000) showed substantial changes in the composition of line catches from similar reef sites sampled on the Agulhas Bank in 1931-1933 and 1987-1993. Other investigations of line catches during the 20th century similarly documented concerning depletion of species (Penney et al., 1989;Van der Elst, 1989;Attwood and Farquhar, 1999;Griffiths, 2000). Such evidence contributed toward the government's decision to legislate emergency provisions in the management of line-fish resources in 2000, leading to 70% curtailment of commercial line fishing effort (Blamey et al., 2015). Confounding the effects of fishing regulations was the creation of several coastal marine protected areas (MPAs) between 1970 and 2008, some of which improved the status of line-caught species (Bennett and Attwood, 1991;Kerwath et al., 2013).
Investigations of trawl species composition changes in South Africa have focused predominantly on recent decades. On the west coast, Atkinson et al. (2011b) assessed spatial and temporal changes in the demersal fish composition from trawl surveys conducted between 1986 and 2009, while Kirkman et al. (2015) looked for regime shifts on the west coast of South Africa, Namibia, and Angola over a similar period. The Agulhas Bank demersal community appears to have increased in diversity and evenness over the period of 1986 to 2003 (Yemane et al., 2010), even though fish size indices decreased (Yemane et al., 2008) and abundance-biomass-comparison curves pointed to greater disturbance during the same period (Yemane et al., 2005).
Retrospective analysis of historical trawl survey data (prior to the 1980s) has received little attention in South Africa. Mussgnug (2013) investigated long-term changes in the Agulhas Bank trawl assemblage, by comparing a variety of datasets from historical (1898-1904; 1922-1948) and recent (1985-2010) trawl surveys and a recent scientific observer program (2003)(2004)(2005)(2006). Despite the uncertainty of comparing among different gears and methods, his results reveal notable assemblage changes, showing that several taxonomic groups decreased markedly while others had increased. Mussgnug (2013) concluded that small, fast-growing generalists withstood expanding fishing pressure better than large, slow-growing, low-fecundity taxa.

Standardization of Fishing Gear and Repeat Surveys
The materials, construction and methods used in trawling operations influence the fishing performance, specifically the relative selectivity of fish that are caught due to their size, shape and behavior (Engås and Godø, 1986;Dealteris et al., 1989;Johnson et al., 2008). Researchers typically try to minimize fishing performance effects by standardizing survey gear and methods throughout a survey series. If that is not possible, estimation of selectivity factors, usually for each species, is a way of addressing the problem when comparing among gear types (Sparre and Venema, 1998). Studies contrasting trawl records from different periods generally have to grapple with this potential bias of varying (and unquantified) selectivity among gear technologies (Rijnsdorp et al., 1996;Greenstreet et al., 1999;Rogers and Ellis, 2000;Mussgnug, 2013;Novaglio et al., 2020), which limits certainty in those results. A few examples exist where researchers set out to repeat historical trawl surveys and thus had the opportunity to imitate trawl gear used previously.
On the south-east coast of England, McHugh et al. (2011) repeated trawl surveys in 2008-2009 at two sites that had been surveyed during 1913-1922. They used an otter trawl net that was the same size and had similar mesh to the historical gear. Their comparison documented notable changes in the assemblage, predominantly in elasmobranch species, which included dramatic declines and an apparent extirpation of angel shark (Squatina squatina).
Employing data from 1964 trawl surveys conducted before initiation of a prawn trawl fishery in the Gulf of Carpentaria in Australia, Harris and Poiner (1991) compared the fish assemblage with 1985/1986 repeat surveys that used the same fishing gear. They documented substantial changes in the fish assemblage. Of the 82 taxa analyzed, 18 had decreased in abundance and 12 had increased significantly. Trawling pressure was considered the likely cause of these changes although environmental variability due to river input and sediment changes also might have had an influence (Harris and Poiner, 1991).
In New South Wales, Australia, Andrew et al. (1997) and Graham et al. (2001) compare upper continental slope trawl survey catches from 1976/1977, a time when commercial trawling was expanding into those areas, with a repeat survey using the same vessel and similar gear 20 years later. Although they did not investigate the assemblage structure, their results showed that both target and non-target species decreased substantially between survey periods and the authors concluded that fishing pressure was likely responsible for such declines.

Unique Pre-exploitation Data From South Africa
The development of a South African trawl industry was prompted by impressive catch records from the government-led exploration of trawl grounds between 1897 and 1906 (Scott, 1949). Recent digitization of detailed survey records from those early research efforts provide an unprecedented view of the demersal ecosystem at a time when human pressures on those grounds are considered to have been negligible (Currie, 2017). As trawl catch records collected prior to significant exploitation pressures are rare globally (though see Silvestre et al., 1986;Harris and Poiner, 1991;Andrew et al., 1997;Klaer, 2001;Kongprom et al., 2003;Novaglio et al., 2020), our knowledge of demersal change stems predominantly from previously impacted systems. The early South African trawl survey data provided a valuable opportunity to investigate baselines and subsequent changes from a minimally disturbed temperate demersal ecosystem.
Despite their obvious value, quantitative contrasts between the historical records and more recent survey data are challenging due to the substantial changes in trawling technology over time (Rijnsdorp et al., 1996;Rogers and Ellis, 2000;Mussgnug, 2013). We addressed this challenge by replicating the trawl surveys conducted by the Pieter Faure survey vessel more than a century previously. The aim was to assess long-term changes in the assemblage composition of trawl-caught fish on the inshore Agulhas Bank. Survey sites were revisited in 2015 with reconstructed trawl gear and methods that imitated the fishing power of 1903/1904 surveys. Changes in the fish assemblage were investigated and key taxa driving such changes were considered. Besides detection and description of long-term change, this experiment provided a measure of current demersal fish assemblages relative to their near-pristine state. Results reveal the scale of change, reiterate the value of pre-disturbed baselines and provide informative context for management of demersal ecosystems.

Selection of Data and Re-survey Sites
Three sites were identified as suitable for repeat surveys, named here as Cape Infanta, Mossel Bay and Bird Island (west to east). They were chosen due to spatial aggregations of historical trawl samples within the modern commercial trawl footprint. Once sites were identified, trawl samples were selected for inclusion in analyses based on the following criteria: (1) they included a start and end time or trawl duration, of which the latter was limited to <200 min; (2) they commenced at least 30 min after sunrise and ended before sunset; (3) their start locations were positioned outside current MPAs or trawl exclusion zones; (4) they were not associated with comments suggesting problems or deviation from normal trawling practices; and (5) their survey dates were restricted within a 3-month period at each site. The location, selection grids and trawl samples are illustrated in Figure 1.
Qualifying trawl samples were plotted in geographic information software (QGIS Development Team, 2015). A sampling grid was created, the width and height of each cell equal to the approximate length (1.3 nautical miles, ∼2.4 km) of a standard 30-min research trawl at a towing speed of 2.6 knots. The survey grid thus defined sampling units from which re-survey cells could be selected. Re-survey targets were chosen by random selection, without replacement, from those grid cells that contained the start location of at least one historical trawl.
Prior to fieldwork, the target sample sizes for repeat surveys were considered using univariate and multivariate approaches (Currie, 2017). To assess the impact of increasing sample sizes on multivariate precision, the pseudo multivariate dissimilaritybased standard error (MultSE) was calculated (Anderson and Santana-Garcon, 2015). Specifically, the 'MSEgroup.d' function (Anderson and Santana-Garcon, 2015) was applied to the historical survey data with site as a grouping factor and using a Bray-Curtis dissimilarity matrix based on square-root transformed catch data that had been standardized by their trawl distance (see below). Profiles of decreasing MultSE with increasing sample size were visually examined and indicated minimal precision increases beyond 20 samples per site (Supplementary Figure S1). A target sample size of 25 trawls per site was set. Unless specified differently, data preparation and analyses were conducted in the R programming environment (R Core Team, 2016).
Although recent decades of government-led research trawl surveys were not contrasted directly with historical catches (due to differences in gear technology), they were employed to assess the magnitude of seasonal or interannual variability. These recent  demersal trawl survey data were obtained from the Department of Agriculture, Forestry and Fisheries and are referred to as DAFF data hereafter. They were collected during annual or bi-annual stratified random surveys of the south coast continental shelf between 1986 and 2015. The sampling protocol for the DAFF data is detailed by Yemane et al. (2010) and is similar to that described below for the repeat surveys. The trawl gear used to collect these data changed in the mid-2000s, as detailed by Atkinson et al. (2011b). Both sets of gear consisted of a 55 m (180 ft) German otter trawl and are similar in many aspects, yet the fishing performance (and selectivity among taxa) may vary between the two. Therefore, analyses in this study treated the DAFF 'old gear' (1986-2003; 2006) and 'new gear' (2004-2005; 2007-2015) separately. DAFF survey samples were selected if they fell within the spatial grids described above, if they had a duration between 20 and 35 min, a towing speed of 3-3.5 knots and did not suffer net damage or other disruptions that may have influenced the catch.

Trawl Gear
A historical Granton trawl net and flat wooden trawl doors were constructed specifically to imitate those used on the Pieter Faure between 1897 and 1906. Details of the gear are provided in Currie et al. (2019a) and only a salient summary is provided here.
The trawl net was based on the plan and descriptions provided by Kyle (1903) and was made from Manila hemp and a mixed polyester-polyethylene alternative with similar properties. The headline consisted of a 28-mm Manila rope, 27.4 m long. The ground-rope consisted of a 20-mm fiber-core cable served with Manila rope to make up a total diameter of 96 mm. Mesh sizes varied from 150 mm stretched mesh in the forward parts of the net, to 100 mm in the trailing parts of the net. The trawl doors were 395-kg, flat, wooden doors with a heavy steel frame and shoe, which were shackled directly onto the net headline and footrope.

Repeat Fieldwork
In February, March and May 2015, repeat trawl surveys were conducted at Mossel Bay, Cape Infanta and Bird Island, accomplishing 25, 25, and 23 trawl samples respectively (Figure 1). These were compared with 42, 60, and 54 historical trawls. The repeat surveys were conducted from a commercial trawler, the MFV Leeukop, owned by the inshore division of Viking Fishing Company (Pty) Ltd. One of the last remaining side-trawlers in the local industry, the vessel was built in 1963, is 20.6 m in length, has a gross tonnage of 96.26 tons and an indicated engine shaft power of 189.41 kW.
Survey trawls commenced half an hour after sunrise and were completed prior to sunset. If conditions were overcast, those time limits were shifted 30 min later and/or earlier, respectively. The tow period, measured between settling of the net on the sea-floor (as observed from vessel speed) and the start of hauling, was approximately 30 min. The targeted towing speed was 2.5 knots. The length of trawl warps was set at ∼5.5 times the depth of water after initial attempts with shorter ratios failed to provide good ground contact of the doors and net.
After each trawl, the content of the net was sorted to species level where possible. Reference samples of unidentified fauna and voucher specimens were stored for examination ashore. The catch weight of each species was measured using calibrated marine scales (Marel M1100). A representative sub-sample was selected for enumeration when large catches were encountered. In such cases the number of fish caught was calculated as the mean unit weight from the sub-sample multiplied by the weight of that species' catch.
A field identification error of Polysteganus coeruleopunctatus (blueskin) meant that this taxon was recorded with Pterogymnus laniarius (panga) at Bird Island. P. coeruleopunctatus normally occurs further east than Bird Island (Smith and Heemstra, 1986) and was not expected in the trawl catches. While their catches are not differentiated in analyses, the great majority are expected to have consisted of P. laniarius and are referred to as such hereafter.
To monitor depth and temperature, a compact CTD probe was attached to the headline of the trawl net, about 1.5 m behind the starboard trawl door. During the first survey (Mossel Bay) a RBR TDR-2050 device was employed and during the following surveys a RBR Concerto unit was used. In both cases, the device was set to record the pressure (depth) and temperature every second.

Data Preparation
This study was restricted to classes Actinopteri, Elasmobranchii, and Holocephali, which are broadly included in the term 'fish' hereafter. Historical records resolved 27 taxonomic groups that could inform analyses (Table 1), the majority of which were at species-level but included five genera, one family, three orders and one mixed group of remaining 'other' fish. Grouped taxa were assigned the lowest taxonomic classification that could be assumed with confidence. In comparison with subsequent datasets, catches of taxa were summed within the same groups and contrasts were restricted to taxa that were recorded throughout each of the periods.
Data from the RBR CTD sensors were processed using RBR Ruskin software (version 1.8.14). Profiles of measured variables were exported as text files for further analysis.
The primary unit of measurement in the DAFF trawl survey data is mass (kg). The associated length-frequency records were Targeting indicated for trawl (T) and line (L) fishing sectors when the taxon is retained or indicated as discarded (-d) if not. Identification issues are outlined in the Notes column. *Most sharks discarded other than certain sizes of Galeorhinus galeus and Mustelus mustelus.
employed to convert catches into abundances of fish. If a subsample had been taken for length measurements, a mean weight per fish was calculated from the sub-sample and scaled up to the catch weight. A few cases existed where no length-frequencies were recorded for a particular species' catch. In such cases, the mean weight per fish was calculated for the relevant site, using remaining trawl samples that did contain weight and length data.

Comparison of Historical and Re-survey Catches
Analyses used catches standardized by the distance of the trawl tow, assuming a constant net width across the study: where C st is the standardized catch (count · nautical mile −1 ); C is the count for the specified taxon; d is the trawl distance; r is the duration of the haul (hours); andv is the average trawl speed (knots). The Pieter Faure trawl speed was estimated to be 1.35 m · s −1 (2.62 knots; Currie et al., 2019a). The re-survey trawl speed was calculated as the GPS-logged distance over the trawl duration. Trawl distance (duration * speed) was considered a marginally more accurate measure of trawl effort than duration (at least for the re-surveys) and was therefore used in analyses.
Due to logistical constraints, the re-survey trawl durations (∼30 min) had to be shorter than those of the historical surveys (60-200 min). The trawl duration was recorded historically as 'time trawl down, ' which was the difference between the 'trawl hauled' time and the 'trawl shot' time. For the re-survey trawls, 'fishing duration' (period that the net is settled on the sea-floor) was estimated using profiles from the temperature-depth sensors attached to the net. The latter method is a more accurate measure of the fishing effort, as it excludes the sinking time of the gear (mean = 3.5 min ± 0.3 SE across re-survey records). Inclusion of sinking duration would have minimal effect on the longerduration historical samples but proportionally greater influence on the catch rates of shorter re-survey trawls (Battaglia et al., 2006), hence the more accurate 'fishing duration' was used for standardization of the re-survey samples.
Because the DAFF survey net is a different size, analyses including those survey data standardized catches by swept area (rather than distance): where C sw is the swept-area standardized catch (count · nautical mile −2 ); a is the area of sea-floor swept by the trawl (nautical mile 2 ); w is the width of the mouth of the trawl net (m). Further analyses were based on a dissimilarity matrix of the community catch data. To down-weight the dominant influence of abundant taxa, dissimilarity matrices were based on square-root transformed standardized catches. Ordination of trawl samples was achieved by non-metric multidimensional scaling (NMDS) of a Bray-Curtis dissimilarity matrix. This was applied using the 'metaMDS' function from package 'vegan' (Oksanen et al., 2016).
The null hypothesis that there was no difference in the multivariate community composition between periods was assessed using a permutational multivariate analysis of variance using distance matrices (Anderson, 2001;hereafter PERMANOVA). A Bray-Curtis dissimilarity matrix was regressed against covariate depth and categorical variables site, period and a site-period interaction to account for potential depth and between-site influences while testing the period effect: where D is the dissimilarity matrix, β 1 represents the intercept, β 2 the estimated coefficient for depth, β 3 the estimated coefficient for the site effect, β 4 the estimated coefficient for the period effect, β 5 the estimated coefficient for the site-period interaction, and ∈ the residual error term.
The 'adonis' function in 'vegan' (Oksanen et al., 2016) was used to run the PERMANOVA models, with 9,999 permutations restricted to within-site levels. Because the period effect was of main interest, after having accounted for depth and site effects, we used sequential sums of squares and added the period effect last in the model structure. Maximum variability is attributed to the remaining main effects in this way, thereby providing a conservative estimate of the period effect.
To assess multivariate differences among pairs of surveys, a pairwise PERMANOVA analysis was also conducted among the levels of a combined period-site factor, including depth as a covariate. The Holm adjustment for multiple comparisons (Holm, 1979) was applied to avoid an inflated likelihood of obtaining a significant result.
PERMANOVA models are robust in the face of multivariate heterogeneity, but not if the experimental design is unbalanced (Anderson and Walsh, 2013). Therefore, the results from the unbalanced design were augmented with parameter ranges obtained from 1,000 randomly re-sampled (without replacement), balanced datasets containing 23 trawl samples per site per period.
To ascertain whether there may have been a change in multivariate heterogeneity between periods, a permutation test of the homogeneity of multivariate group dispersions was employed (PERMDISP2; Anderson, 2006). The same dissimilarity matrix was used as in the multivariate model and ordination above. Pairwise tests among period-site level combinations were assessed using 9,999 permutations and employing the Holm adjustment of p-values for multiple comparisons (Holm, 1979). The multivariate group dispersions were implemented using the 'betadisper' function in package 'vegan' (Oksanen et al., 2016).
A two-way (period * site) similarity percentage analysis (SIMPER; Clarke, 1993) was used to investigate taxa that discriminate between periods across sites. Similarly, a pairwise SIMPER was used to identify the taxa discriminating periods for each site. These analyses were based on the same Bray-Curtis matrix described above and conducted in Primer-E (Clarke and Gorley, 2006).
To contrast the magnitude of variability between historical and re-survey samples with that among modern DAFF research trawl surveys, a multivariate measure of pseudo variance (MultV) was employed as per Anderson and Santana-Garcon (2015). Consistent with other analyses, MultV was computed from Bray-Curtis dissimilarity matrices derived from square-root transformed, swept-area standardized catches (count · nm −2 ). At each site, 500 pairs of historical and re-survey samples were randomly selected with replacement. DAFF samples were stratified by their temporal proximity to each other, so that the difference between pairs of seasonal (≥2 months and <12 months), interannual (≥12 months and <18 months) and multi-year (≥18 months but ≤60 months) samples could be calculated. As the sample sizes were relatively small for the DAFF pairs, all possible pairwise contrasts were included. The magnitude of MultV values were visually contrasted on scatterplots using their mean and bootstrapped non-parametric 95% confidence intervals.
As some uncertainties remained in the size of the historical net and the speed at which it was towed (Currie et al., 2019a), their potential effect on results was examined. This was accomplished by repeating the full set of analyses on swept-area standardized catches, using the extreme combinations of plausible bounds for headline length and trawl speed. This resulted in two additional sets of results, the one based on the assumption of a 26.8 m (88 ft) headline length and trawl speed of 1.10 m · s −1 (2.15 knots); the other on a 36.6 m (120 ft) headline and a speed of 1.61 m · s −1 (3.13 knots; Currie et al., 2019a).
In a similar manner, the impact on results of potential biases of depth were examined. Again the entire set of results was re-generated, excluding re-survey Bird Island samples with a depth ≥ 100 m. Differences between these additional results and the main set presented below were assessed.

RESULTS
The targeted sites were successfully re-surveyed in 2015, revealing catches similar in abundance but substantially different in taxonomic composition to those of the historical surveys. The seasonality, time-of-day and trawl depths of samples were similar between periods, but trawl durations (and therefore distances) were greater in the historical surveys (Supplementary Figure S2) and their potential effect on results is discussed below.
Ordination of the standardized catches showed a clear separation of historical from re-survey trawl samples (Figure 2), predominantly along the primary x-axis. Site groupings were also discernible in the ordination, partially separating along the secondary y-axis. The large period effect was confirmed by the PERMANOVA results, accounting for 41-47% of the variance explained by the model ( Table 2). The highly significant site effect explained 10-13% of the model variance. Together with their interaction, these two variables jointly explained 64-69% of the model variance. The depth of trawls also affected multivariate distances among trawls, but to a lesser extent (6-10% variance). The R 2 estimates from the full (unbalanced) dataset fell within ranges of those estimated by the under-sampled (balanced) models, except for the depth FIGURE 2 | Non-metric multidimensional scaling (NMDS) plot based on Bray-Curtis dissimilarities of square-root transformed standardized catch. Historical and re-survey samples are indicated by circles and triangles respectively, while sites are separated by color (Cape Infanta blue, Mossel Bay green, Bird Island red). Bounding ellipses represent 95% probability contours estimated for each period-site level. effect, which was underestimated in the full design ( Table 2). All terms combined, accounted for 72-76% of total model variance. Results of a pairwise PERMANOVA suggested that all period-site levels were significantly different from each other (p ≤ 0.01), except for the re-survey samples from Cape Infanta and Mossel Bay. Shifting focus to the dispersion among samples, historical trawls tended to cluster closer together than the repeat surveys (Figure 2). The PERMDISP results confirmed that dispersion was different between periods (F = 101.3, p ≤ 0.0001). This difference appears to be caused mainly by divergence between the modern Bird Island samples and those of the other two sites, and greater variance among Mossel Bay re-survey samples (Figure 2). Dispersion was significantly different among periodsite levels (F = 14.7, p ≤ 0.0001), but a pairwise test revealed that only Mossel Bay differed between historical and re-survey periods ( Table 3). The modern survey at Mossel Bay was significantly more heterogeneous than all other period-site levels (Figure 2 and Table 3).
A large dissimilarity (81.7%) of historical and modern assemblages was supported by SIMPER analysis. Focusing on the distinguishing taxa that contributed 90% to between-period dissimilarity, six teleost taxa and Torpediniformes (electric rays) were more abundant historically, whereas four teleost taxa and Squalus spp. (spiny dogfish) were more numerous in the re-surveys (Figure 3). Argyrosomus spp. (kob), Pterogymnus laniarius (panga), Austroglossus pectoralis (east coast sole) and Chelidonichthys spp. (gurnards) contributed over 50% of the dissimilarity between periods and were relatively consistent in their distinguishing contribution (dissimilarity/SD > 1; Supplementary Table S1). The first three of those taxa dominated the composition of historical catches at each of the three sites, jointly contributing 82, 70, and 84% of Cape Infanta, Mossel Bay and Bird Island catches respectively. Contrasting the repeat surveys, Argyrosomus spp. were absent and the summed contribution of the further two species had dropped to between 1.5 and 5.5%. Other species that were nearly absent in repeat surveys included Argyrozona argyrozona (carpenter seabream; three individuals) and Rhabdosargus globiceps (white stumpnose; one individual), whereas they made up 8 and 4% of historical catches, respectively. Standardized abundances of Chelidonichthys spp., Trachurus capensis (Cape horse mackerel), Squalus spp., Merluccius capensis (shallow-water hake) and Galeichthys feliceps (white sea catfish) increased between periods, their joint contribution rising from 3% of historical catches to 85% of re-survey catches. Inspection of DAFF trawl survey data showed that by the late 1980s (1986)(1987)(1988)(1989)(1990), the proportional contribution of the same suite of taxa matched more closely the re-survey than the historical assemblage structure: Argyrosomus spp., P. laniarius, and A. pectoralis contributed a cumulative average of 11% (by weight), whereas the five taxa that were most numerous in re-survey catches contributed 75% of the DAFF catch biomass. The longterm differences were substantially greater than shorter-term (seasonal, interannual and 5-year) variability characterized in the DAFF dataset (Supplementary Figure S3).
Due to the survey design, depths were expected to be similar between periods. However, at Bird Island, a bias toward deeper sampling depths was noted in the re-survey trawls (mean = 90 m) relative to historical surveys (mean = 78 m; Supplementary Figure S2). Repetition of analyses on a dataset that excluded resurvey Bird Island samples with a depth ≥ 100 m showed no material impact on results other than reducing statistical power.

DISCUSSION
This research revealed dramatic change in demersal fish communities on South Africa's Agulhas Bank over a 111-year period. The contrast between historical and re-survey records convincingly rejected a null hypothesis of no difference in catch composition between survey periods. The substantial temporal separation of assemblage composition was evident from an unconstrained ordination plot and supported by the period effect explaining almost half of the variance among trawl samples.
The striking transformation of the community assemblage was driven by previously dominant species that were absent or rare in the re-survey, as well as a suite of taxa that were numerically sparse in historical surveys but increased to make up the majority of 2015 assemblages (Figure 3). The most severe catch declines were on average detected in Argyrosomus spp., P. laniarius, A. pectoralis, A. argyrozona and, specifically at Mossel Bay, R. globiceps. The relatively consistent proportions across sites (roughly a quarter each) of Argyrosomus spp., P. laniarius, and A. pectoralis, suggests that these taxa may have historically dominated similar habitats in this region and are expected to have played a substantial role in the trophic ecosystem. Their large size and varied diet makes this especially relevant for Argyrosomus spp. (Smale and Bruton, 1985), which we presume exerted substantial predation pressure on the surrounding community.
The drastic decreases of the above taxa were accompanied by large increases of Chelidonichthys spp., Squalus spp., M. capensis, and T. capensis caught in repeat surveys. Whereas the joint historical contribution of these taxa was 5% at Bird Island and 2% at Cape Infanta and Mossel Bay, their numbers dominated the modern assemblages (73-87% of catches). Abundances of these taxa were less consistent across re-survey sites, with M. capensis more numerous at the two western sites and Squalus spp. exceptionally abundant at the eastern site of Bird Island. Inspection of the earliest years (1986)(1987)(1988)(1989)(1990) from current government trawl survey records indicated that the above transition of dominant taxa had taken place before the start of the modern-era trawl surveys. Without the historical context, this substantial shift would be overlooked, highlighting the value of re-establishing these baselines.

Factors Affecting Assemblage Changes Between Periods
During the 20th century, fishing impacts are likely to have been a dominant driver of ecological change in marine communities in South Africa (Moloney et al., 2013) and many other regions (Jackson et al., 2001). Several of the taxa that declined in our study are fishery targets and we assume much of the change observed is associated with fishing pressures. However, some targeted taxa increased or remained similar in abundance, signaling complexity in responses. Additional pressures no doubt played a role, which could including climate change, changes in riverine inputs, pollution and hydrocarbon exploration and extraction (Moloney et al., 2013). Besides evolving fishery regulations, the closure of bays to trawling (which spanned 1928-1986) and proclamation of new MPAs (1970-2019) may have mitigated some of the impacts. Unraveling past effects of the various pressures on demersal fish assemblages will require substantial additional data, much of it not currently available. As they may indicate potential drivers of change, we instead considered factors that might explain why certain taxa had declined and others had increased or maintained numbers, following 111 years of mostly increasing pressures.
Life-history characteristics are frequently cited as determinants of the sensitivity of populations to fishing pressures, including fecundity, age at maturity, spawning strategies, growth rate, longevity and body size (Rogers and Ellis, 2000;Stevens et al., 2000;Graham et al., 2001;Genner et al., 2010;McHugh et al., 2011). Declines of elasmobranch taxa, for example, are often linked to their low fecundity, slow growth and/or old age at maturity (Stevens et al., 2000;Graham et al., 2001;Ferretti et al., 2010;McHugh et al., 2011). Yet in the current study, some elasmobranch catches appeared to increase between periods. While abundances of Torpediniformes declined across sites, the numbers of Squalus spp., Myliobatiformes and Carcharhiniformes increased overall (Supplementary Table S1). It is worth noting that those increases were predominantly due to large catches at only one site (Bird Island; Supplementary  Table S2) and therefore may not be representative of the region. Rajidae appeared to decline at the two western sites but were also caught in far greater abundances in the Bird Island re-survey. The reasons for these high elasmobranch catches at Bird Island in this study are not clear.
Habitat preferences may provide insight into the potential drivers of observed assemblage changes. Multiple taxa that declined heavily are associated with reef habitats to varying degrees, including Argyrosomus spp. (Griffiths, 1997), P. laniarius (Booth and Buxton, 1997), A. argyrozona (Brouwer and Griffiths, 2005), and R. globiceps (Mann, 2013). On the contrary, the majority of species that increased between surveys tend to be those associated with soft substrates (e.g., Chelidonichthys spp. (Booth, 1997), Squalus spp. (Stehlik, 2007), G. feliceps (Tilney and Hecht, 1990) and Rajidae (Serra-Pereira et al., 2014;De Vos et al., 2015), or inhabit both soft and hard benthic substrates (M. capensis; Fairweather et al., 2006;Mann, 2013). Trawling has continued on these grounds for over a century and reef-like habitat or structure-forming benthic communities that may have been present historically are likely to have been degraded or removed (Engel and Kvitek, 1998;Auster and Langton, 1999;McConnaughey et al., 2000;Moran and Stephenson, 2000). Trawling can reduce habitat complexity in a number of ways (Auster and Langton, 1999), including the removal and damage of emergent and highbiomass fauna, which may become replaced by smaller-bodied and infaunal animals . Initial trawl passes may cause the greatest impact in sensitive epifaunal communities that include three-dimensional structure (Collie et al., 2000;Hiddink et al., 2006). Sainsbury et al. (1993) document a case in which trawl-related habitat changes appeared to impact the fish assemblage: Two decades of escalating trawl pressure led to marked changes in the fish catch composition in north-western Australia. Examination of the composition changes showed that previously dominant taxa (Lutjanus spp. and Lethrinus spp.), which associated with dense epibenthic cover, had been replaced by two commercially less-valuable fish genera (Nemipterus spp. and Saurida spp.), which preferred areas of open sandy substrate (Sainsbury et al., 1993). The catches of epibenthic fauna (mainly sponges, alcyonarians and gorgonians) had declined over the same period.
It is conceivable that habitat-related assemblage changes similar to those documented in north-west Australia (Sainsbury et al., 1993) played out during the early years of trawling, both in South Africa and elsewhere. Trawl intensity has been linked to changes in epifaunal and infaunal communities on South Africa's west coast (Atkinson et al., 2011a), but those changes were from an environment that had been extensively fished prior to its study. Knowledge of benthic invertebrate baselines and their subsequent changes are difficult to reconstruct (Callaway et al., 2007), as fishers and scientists tended to pay little attention to documenting catches of non-target invertebrate fauna until relatively recently. Our investigation lacked data on the epifaunal biogenic habitat that occurred historically and we cannot make predictions about how long its recovery might take if fishing impacts were to cease (though see e.g., Collie et al., 2000;Kaiser et al., 2006;Hiddink et al., 2019). Further work is needed to evaluate the impact of trawling on Agulhas Bank habitats, their recovery potential and link with observed fish assemblage changes.
Geographic and depth ranges are additional factors that might contribute toward the sensitivity of species to exploitation. Argyrosomus spp., P. laniarius, A. pectoralis, and R. globiceps are restricted to relatively shallow depths (≤120 m), generally close to the coastline. Sought-after fishing targets limited to shallower coastal areas and in reach of the majority of fishers and fishing gears will inherently suffer greater exposure to fishing pressures. Most taxa that increased (Chelidonichthys spp., Squalus spp., T. capensis, M. capensis) have a substantially wider range in terms of distance from the coastline and depth. Their greater distribution range may enable larger populations, but such populations are also more likely to have maintained refuges that are less accessible or uneconomical for fishers to exploit.

Increased Variability and Spatial Differences in Catch
The variability (dispersion) of Mossel Bay samples was inflated compared to all other period-site combinations (Figure 2). One explanation might be an impacted demersal community due to persistent anthropogenic pressures: Caswell and Cohen (1991) predicted that disturbance would increase heterogeneity within a community and Warwick and Clarke (1993) demonstrated that increased variability was associated with stress or degradation across four examples of widely varying marine communities.
Of the three survey sites, Mossel Bay appears to be exposed to the highest level of anthropogenic pressures (Sink et al., 2012a), which include fishing effort and its proximity to an industrialized town and harbor servicing offshore hydrocarbon extraction.
Summed across taxa, catches were historically lowest at Bird Island and increased to highest levels at Cape Infanta. This contrasts with the re-survey data, which showed Bird Island catches to be larger than at the two western sites, due to considerable abundances of Chelidonichthys spp., G. feliceps and a host of elasmobranchs (Squalus spp., Rajidae, Myliobatiformes and Carcharhiniformes). Inspection of the last decade (2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015) of government trawl survey data suggested that some of these catch contributions were unusual (Currie, 2017), but that catches of Chelidonichthys spp., Squalus spp., Myliobatiformes and all taxa combined have recently been higher at Bird Island than at the other two sites.
A measure of multivariate variability between historical and re-survey samples suggests that demersal communities may have transformed less at Bird Island than at Cape Infanta and Mossel Bay. Multiple reasons might contribute to a less impacted Bird Island site. Relative to the other sites, trawl intensity has been lower at Bird Island in recent decades (Currie et al., 2019b). In addition, the site is adjacent to the Bird Island MPA and in relative proximity to the eastern-most boundary of commercial trawling activities during the last four decades (Booth and Hecht, 1998). Its demersal communities are, therefore, likely under less direct fishing pressure and may benefit from spillover from untrawled areas nearby. The Cape Infanta site is also adjacent to an MPA, which has benefited surf-zone fish species (Bennett and Attwood, 1991). MPAs can offer protection against multiple pressures and safeguard resilient assemblages (Bennett and Attwood, 1991;Dayton et al., 2000;Micheli et al., 2012;Mellin et al., 2016), together with their critical habitats (Lindholm et al., 2001;Kritzer et al., 2016). The benefits of such MPAs to demersal trawl communities have not been tested in South Africa. Neither has the effect of closing bays to trawling, some of which were implemented as early as 1928 (Scott, 1949). Contrasting habitats and demersal communities inside and outside these protected areas could provide valuable insight on trawling-induced habitat impacts, their recovery and the efficacy of protection.

Pre-disturbance Baselines
Few other studies have had access to demersal survey records prior to, or early in the development of a trawl fishery. Of those that did (Silvestre et al., 1986;Harris and Poiner, 1991;Andrew et al., 1997;Graham et al., 2001;Klaer, 2001;Kongprom et al., 2003;Mussgnug, 2013;Novaglio et al., 2020), most did not assess assemblage composition changes in the same way that we did. Nonetheless, some commonalities with their results were identified. Examining commercial catch rates between 1918and 1957, Klaer (2001 documented strong declines in the south-east Australia trawl fishery, despite substantial technological improvements during that time. Marked changes in the composition of catches included the virtual disappearance of two species that had been abundant during early years of the fishery. Before the same fishery expanded into deeper waters, trawl surveys were conducted on the upper continental slope in 1976/1977. Twenty years later, when a repeat survey was conducted, the previously dominant dogsharks (Centrophorus spp.) and redfish (Centroberyx affinis) had declined to <1% and <5% of initial catch rates, respectively (Andrew et al., 1997). Comparing south-east Australian trawl surveys from the 1980s with those conducted prior to commercial trawling in 1909-1910, Novaglio et al. (2020 documented substantial (83-99%) declines for several demersal fish families and notable shifts in catch composition, which they attributed predominantly to fishing. Unlike our study, they found a marked catch rate decline of 77% for the entire assemblage. On the northern coast of Australia, Harris and Poiner (1991) showed that two decades after the development of prawn trawling in the Gulf of Carpentaria, the combined catch rate of all fish had declined to less than half of initial values. Their results revealed substantial changes in the fish assemblage within 21 years of the fishery, including collapse of a previously dominant taxon (Paramonacanthus spp.).
Mussgnug (2013) included the same historical data and investigated comparable geographic areas to those of the current study, but with a focus on chondrichthyan taxa. He showed similar notable differences between historical (1898-1933) and recent  assemblages, including decreased catch contributions of P. laniarius, A. inodorus, A. pectoralis, A. argyrozona, Umbrina canariensis and R. globiceps. He also reported increases in catch contributions for T. capensis, Chelidonichthys spp. and M. capensis. The main discordance between the two studies pertains to certain elasmobranch taxa: Whereas Mussgnug (2013) emphasized decreases of all elasmobranch taxa, multiple elasmobranch groups increased between periods in our results (Supplementary Table S1). This may partly be explained by the unusually high catches of elasmobranch taxa encountered at Bird Island. A benefit of current results over those of Mussgnug (2013), is that differences in trawl gear and trawl speeds were explicitly examined and accounted for, reducing uncertainty and strengthening confidence in results.
In the North Sea, where commercial fishing (including trawling) was extensive prior to early trawl surveys, demersal assemblage changes in the 20th century appear to have been less severe. McHugh et al. (2011) investigated long-term changes in the demersal fish community of the western English Channel. Although multivariate differences between 1913-1922 and 2008-2009 did separate statistically, their ordination did not show strikingly disparate periods (as seen in Figure 2). Differences were driven predominantly by alterations in the elasmobranch assemblage, which was interpreted as evidence of the loss of larger, slow-growing taxa due to fishing pressures (McHugh et al., 2011). Other studies showed relatively subtle abundance changes of mostly less common demersal fish species in the North Sea (Greenstreet and Hall, 1996;Greenstreet et al., 1999;Rogers and Ellis, 2000). To illustrate, Greenstreet et al. (1999) noted that the same five dominant species contributed over 50% of similarity across all four periods they examined between 1925 and 1996.
It appears that the few studies with access to records from the earliest stages of a trawl fishery demonstrate more drastic assemblage changes than those based on baseline data collected many years after the development of fishing impacts. This implies that the magnitude of change measured in long-term comparisons depends largely on the state of communities first sampled. If so, it illustrates the erosive effect of shifting baselines (Pauly, 1995) and the importance of documenting pre-disturbed reference points to counter them. The majority of research into long-term demersal assemblage changes has been based in Northern Hemisphere regions, where decades, if not centuries of trawling preceded the earliest survey records (e.g., Greenstreet et al., 1999;Genner et al., 2010;e.g., Greenstreet and Hall, 1996;Rijnsdorp et al., 1996;ter Hofstede and Rijnsdorp, 2011;McHugh et al., 2011). Such studies provide useful insight on change and previous ecological states, yet the historical systems they document had suffered substantial impacts and do not represent near-pristine ecosystems (Jackson et al., 2001;Thurstan et al., 2010;McHugh et al., 2011). This study, however, had access to trawl survey data from almost pristine demersal communities. It contributes valuable knowledge of the baseline fish assemblages on the inshore Agulhas Bank and how these changed during the intervening century. Such information is pertinent to species and ecosystem assessments, could help estimate the productivity potential of undisturbed ecosystems and provide benchmarks against which change or the efficacy of management strategies can be evaluated. Estimates of unfished species abundances provide useful diagnostics and informative priors to assess or constrain fishery stock assessment models (Hilborn and Walters, 2001) and inform the choice of baseline conditions used in ecological models (e.g., Mackinson, 2001;Watermeyer et al., 2008).

Examination of Assumptions
Factors that might bias comparisons between historical and re-survey trawl catches were accounted for as far as possible. Season, time of day, geographical locations and replication of gear were standardized. Nonetheless, contrasts of the historical and re-survey catches required certain assumptions, which evoked their examination. These included that (1) the comparisons were not biased by remaining uncertainty in the historical gear dimensions and trawling speeds; (2) differences in trawl durations did not skew catch rates; (3) a depth bias at Bird Island did not impact results; (4) shorter-term interannual variability could not explain the changes observed; and (5) the taxonomic resolution imposed by historical records did not confound results or inferences drawn. These assumptions are unpacked in detail in the Supplementary Material. In summary, none of these factors appeared likely to have caused or substantially biased the observed differences between periods: (1) The effect on standardized catches that might result from uncertainties of gear size and trawling speed had negligible effect on results and interpretations reported; (2) examination of the literature suggested that a trawl duration impact would likely be marginal; (3) omitting historical Bird Island stations that caused a slight depth bias did not cause material differences in results; (4) seasonal, interannual and multi-year variability in trawl surveys from recent decades (expected to include the effects of environmental variability at these time scales) could not account for the magnitude of differences measured between periods; (5) the reduced taxonomic resolution recorded in the historical surveys would likely act to reduce (hide) rather than magnify differences observed between periods.

CONCLUSION
This quantitative comparison of historical and repeat trawl survey data provided evidence of dramatic long-term change in the demersal fish assemblages of South Africa's inshore Agulhas Bank. Although historical data have previously indicated severe declines of line-caught species (Griffiths, 2000), our results highlight the scale of ecosystem change and the remarkable historic abundance of three species that previously dominated these demersal habitats. With this new information, scientists and managers are better equipped to understand change and consider historical reference points in developing strategies to manage and rebuild stocks.
Although we did not anticipate it, the observed fish assemblage changes suggest altered benthic habitats, most likely caused by the cumulative demersal trawling impacts at the studied sites. If substantial trawling pre-dates baseline data collection, such potential habitat impacts are likely to be missed in similar research globally. If trawling has altered the sediment structure (Auster et al., 1996), the benthic habitats, together with invertebrate and fish communities dependent on them, may fail to recover to their previous states even if trawling and other pressures were to be removed.
As noted, additional human impacts, including other fisheries, climate change, pollution and river flow reduction may have contributed to the state of demersal assemblages investigated here. Improved management of these communities may require addressing impacts from multiple sources. The recently declared Agulhas Bank, Agulhas Muds and Addo Elephant National Park MPAs are expected to contribute to the protection and recovery of demersal communities in the studied region. Monitoring and dedicated research efforts are required to track current and future states of those demersal ecosystems.
The severity of the assemblage reorganization, relative to similar investigations elsewhere, is interpreted to be due to the near-pristine state of the historical community surveyed. The majority of previous studies may miss the initial, and perhaps most dramatic, shifts in assemblage composition if they lack baseline data collected prior to substantial exploitation. Researchers and managers need to be cognizant of such missing baselines when interpreting change, estimating reference points in fished ecosystems, and in advancing plans for resource management and protection of seabed habitats.

DATA AVAILABILITY STATEMENT
The historical data are available from SAEON (https://doi.org/ 10.15493/HCDB-8A05). The repeat survey data are available upon request from JC. The DAFF survey data can be requested from the director of Resources Research in the Fisheries Branch of the Department of Fisheries, Forestry and Environment in South Africa.

ETHICS STATEMENT
The animal study was reviewed and approved by Faculty of Science Animal Ethics Committee (University of Cape Town).

AUTHOR CONTRIBUTIONS
This manuscript was based on the Ph.D. work of JC. CA, KS, and LA co-supervized the work. JC conducted analyses and wrote the draft. LA, KS, and CA contributed text, comments and editorial input. All authors made a substantial, direct and intellectual contribution to the work, and approved it for publication.

ACKNOWLEDGMENTS
Fieldwork assistance from Stewart Norman, Stephen Lamberth, and Denham Parker is gratefully acknowledged, as well as the guiding expertise from Noel Less. The Fisheries Management branch at the Department of Fisheries, Forestry and Environment are acknowledged for allowing the use of their data. This work benefited from JC's involvement in the ICES Working Group on the History of Fish and Fisheries.