Patterns of seagrass macrobenthic biodiversity in the warm-temperate Knysna estuarine bay, Western Cape: a review

Knysna estuarine bay in South Africa's Garden Route National Park is that country's most significant estuarine system for biodiversity and conservation value. One outstanding feature is support of 40% of South Africa's—and maybe 20% of the world's—remaining vulnerable and decreasing dwarf-eelgrass, Zostera capensis, whose associated benthic macrofauna has been studied since 2009. For these invertebrates, Knysna comprises several significantly different compartments: sandy mouth; well-flushed marine embayment; poorly flushed central sea-water 'lagoon'; and two disjunct but faunistically similar peripheral regions–marine backwater channels, and low-salinity upper estuary. Although macrofauna ranges from dilute brackish to fully marine, its abundance, local patchiness, and over considerable stretches, species density remains remarkably constant; further, one-third of species occur throughout. Intertidally, all but peripheral compartments are low density and infaunally dominated, while some peripheral areas, and much of the subtidal, are higher density and epifaunally dominated. Overall, seagrass macrobenthos appears maintained below carrying capacity (e.g., by abundant juvenile fish) and of random species composition within a site. Two further characteristics are notable: Unusually, seagrass supports fewer animals than adjacent unvegetated areas, probably because of lack of bioturbatory disturbance in them, and the vegetation cover may ameliorate ambient habitat conditions. Unfortunately, continual heavy and effectively unpreventable exploitation for bait occurs, and chlorophyte blooms have developed because of high nutrient input. Knysna presents a microcosm of problems facing biodiverse and high-value habitats set within areas of high unemployment where subsistence fishing provides the main source of protein and seagrass provides the only source of bait.


Introduction
The permanently open Knysna estuarine bay (34803 0 S, 23803 0 E) is a drowned river valley in South Africa's Western Cape separated from the adjacent Indian Ocean by a narrow gorge (300 m wide, 700 m long, and 4 m deep at low tide) carved through the coastal quartzite ridge by the Knysna River during times of lower sea level. The bay forms part of the open-access Garden Route National Park, and on a basket of criteria, including its size, diversity of habitat, zonal rarity, and biodiversity, is ranked South Africa's most significant estuarine system in terms of conservation importance (Turpie and Clark 2007;van Niekerk et al. 2019). Known locally as the Knysna Lagoon, the system has an area of some 10 km 2 at low tide and 16 km 2 at high tide and receives the inflow of the Knysna River at its head and a large number of smaller streams along its northern and eastern shores. Nevertheless, it is marine-dominated, consequent on low average rates of freshwater inflow and a very large tidal prism, during spring tides equaling 19 9 10 6 m 3 [the largest of any South African estuary (Grindley 1985)], causing semi-diurnal flushing of its main channel. The estuarine bay can be divided hydrologically into three linear compartments, which vary in areal extent and precise geographical position with the tidal cycle and magnitude of river flow: An outer marine bay tidally flushed with cool water from the Indian Ocean and with salinities usually [ 34; a middle, more isolated lagoonal water body also of high salinity (30-34) but with long residence times (c. 4 weeks); and an inner, stratified and well-flushed estuarine region with low and variable salinity (0-30) as a result of freshwater input from the Knysna River (Largier et al. 2000).
This system is also one of the most thoroughly researched of any South African estuary with more than 100 published articles (Russell et al. 2012;Whitfield and Baliwe 2013), work there beginning in 1947 (Day et al. 1951) (see the summaries of Day 1967;Grindley 1985;Russell et al. 2012;and the articles in Hodgson and Allanson 2000). Among other important features, Knysna supports 40% of South Africa's dwarf-eelgrass, Zostera (Zosterella) capensis [or Nanozostera capensis in the recent revision of the Zosteraceae of Coyer et al. (2013)] which may equate to 20% of its world area (Adams 2016;Wasserman et al. 2020). It also forms the only known African locality of the unusual marine valvatoidean gastropod Cornirostra (GBIF 2020), as well as being the main habitat of several other rare seagrass-associated species, including the endangered Knysna seahorse (Hippocampus capensis), and it is one of the only two localities that support the critically endangered seagrass false-limpet (Siphonaria compressa) and a seagrass population of the dwarf cushion-star Parvulastra exigua. Except in the immediate vicinity of the mouth, Z. capensis, together with some mixed Halophila ovalis, occurs virtually throughout the intertidal zone of the system as one continuous bed (Maree 2000), and it also occurs subtidally though more patchily (Wasserman et al. 2020;Barnes and Claassens 2020). Such meadows support well-developed invertebrate macrofaunas that serve the vital functions of consuming epiphytic algal growths and providing the trophic link between microphytobenthic production and that of the larger, more mobile nekton (Murphy et al. 2021). This article synthesizes the main findings of the disparate series of researches conducted on these invertebrate faunal assemblages at Knysna since 2009, re-analyzing the original data where appropriate, with particular emphasis on patterns of macrofaunal assemblage composition, abundance, species richness, and patchiness along the bay's main axial transitional gradient, as well as along the gradient of shelter, located perpendicular to that axis across its eastern section.

General methodology
Patterns described in this review are based mainly on a series of 23 sites positioned to represent the whole area over which Zostera capensis is present intertidally (Fig. 1), with additional comparison between the macrofaunal eelgrass assemblages at some of those sites and the equivalent assemblages in immediately adjacent areas of bare sediment and/or in subtidal seagrass. Sampling was conducted each year between 2009 and 2020 during the austral summer, the research being approved by SANParks and conducted in accordance with their scientific research regulations and requirements. A standard procedure was used, involving series of core samples, each of 0.0027 m 2 diameter prior to 2013 and of 0.0054 m 2 diameter thereafter and of 100 mm depth, taken from continuous stretches of seagrass while still covered by [ 10 cm of water. Cores were gently sieved ('puddled') through 710-lm mesh on site. This sampling procedure collects the smaller and more numerous members of the benthic and epibenthic macrofauna that constitute the large majority of invertebrate biodiversity (Bouchet et al. 2002;Albano et al. 2011), though not the meiofauna nor much scarcer megafauna nor sessile animals attached to the seagrass leaves.
In the laboratory, retained animals from each core were identified to species level wherever possible, with all organismal nomenclature here being as listed in the World Register of Marine Species (www. marinespecies.org), accessed November 2020, except in respect of the currently genus-less microgastropods 'Assiminea' capensis and 'A'. globulus (see Barnes 2017). It should be noted, however, that the specific identity of several animals, especially among the Polychaeta, is questionable because of lack of recent revision; those of South African taxa of Polycladida, Oligochaeta, and Nemertini, and many members of other groups less than 3 mm-4 mm in largest dimension are virtually unknown. Such animals were treated as morphospecies, an operationally appropriate procedure to detect spatial patterns of numbers of species and their differential abundance (Dethier and Schoch 2006;Gerwing et al. 2020  Steenbok Channels separating the two large bay islands from the mainland. The approximate geographical extents of the faunally distinct regions suggested in Fig. 2 are also indicated statistically significant departures from random being determined by Monte Carlo simulation using 9999 iterations. Correlations were assessed using Spearman's rank coefficient S q or the Pearson productmoment coefficient P q as appropriate; number of species per site unit of 30 cores and per region unit of 180 cores was determined by Mao tau rarefaction; curves were fitted using KaleidaGraph 4.5.4; and, where not known, information on life style of individual species was derived from that of close relatives in compendia such as Macdonald et al. (2010).
Multivariate comparison of macrofaunal assemblage composition used hierarchical clustering analysis of S 17 Bray-Curtis similarity, ANOSIM, ANCOVA, SIMPER, and ordination by non-metric multidimensional scaling (nMDS), with 9999 permutations. For such comparison, all data sets were standardized for overall species density (by dividing all ranks by the total number of species in the set) and for sample size (by dividing each species total by the overall number of individuals in the set) to reflect solely differential taxonomic composition and to permit comparison of curve slopes (Passy 2016), a measure of equitability in individual species contribution to the total (Whittaker 1972). Overlaps in quantitative assemblage composition between adjacent regions were measured by the Bray-Curtis similarity index. All multivariate analyses were based on sample sizes of [ 250 animals, well above the minimum number recommended by Forcino et al. (2015).

Patterns in assemblage composition
In total, some 67,000 individual macrofauna, representing 160 species, were examined in 2,100 core samples from the Zostera beds during the study. These ranged from typical freshwater/dilute-brackish species such as Afrochiltonia capensis, Corallana africana and Melanoides tuberculata through to fully marine backwaters and axial channel conditions in the marine section would appear to lie along the Steenbok Channel in that site D falls within one grouping and E in the other. Envelopes enclose sites with the stated levels of Bray-Curtis similarity, for which the approximate geographical locations are shown in Fig. 1 forms such as Gibbula cicer, Limaria tuberculata, Nebalia capensis and Parechinus angulosus. Ordination by nMDS of Bray-Curtis similarity data from the 23 intertidal sites suggested that four significantly different faunal clusters occurred in the system (ANOSIM R = 0.88; P \ 0.0001) (Figs. 1and 2), an essentially similar pattern to that derived earlier using non-standardized (but fourth-root transformed) abundance data (Barnes 2013a). These represented: (i) the sandy mouth region immediately adjacent to the true mouth, (ii) the outer marine embayment, and (iii) the lagoon plus lower-estuary divisions of the main axial channel, and (iv) the fringing backwater-creek system of the smaller, saltmarsh-enclosed creeks and channels that separate the bay's two large islands (each c. 82-84 ha) from the mainland, together with sites in the upper estuary. Separation of the backwaters/upperestuary sites from those along the main axial channel was the most marked, with a Bray-Curtis similarity between the two blocks of sites of only 20%, and the mouth region was an outlier within the axial channel. The three points of change along the longitudinal axis of the bay, however, were not marked by sharp faunal contrasts (Fig. 3). Indeed, SIMPER indicates that most ([ 50%) of the differences are brought about by the relative abundances of just eight common and widespread species, the gastropod molluscs Hydrobia knysnaensis and 'Assiminea' capensis (dominant in region iv), Turritella and Alaba (dominant in i), and Nassarius (dominant in iii), and the polychaetes Prionospio (dominant in iii), and Caulleriella and Simplisetia (dominant in ii). Despite statistically significant regionalization, 32% of the species (representing [ 75% of the total individuals) occurred in all four regions in more than token quantities (Table 1 lists the more numerous of these shared taxa, and Table 2 displays those characteristic of each region). Patterns of relative species abundance within the four regions did not differ (ANCOVA equality of means F = 0.17, P = 0.92; equality of slopes F = 0.59, P = 0.62) (Fig. 4), further indicating similarity between the different local assemblages. Number of species per sample did not vary across test areas of up to 1.5 ha at a given site, whether in the bay or in the lagoon (Table 3). The observation that the seagrass macrofauna of the brackish upper estuary did not differ from that in the fully saline, saltmarsh-enclosed backwater channels of the marine embayment is noteworthy and reinforces the earlier comments of Day (1959) and Barnes (1989) that so-called estuarine faunas may be as characteristic of sheltered areas of fully marine soft sediment as they are of regions subject to low salinity.
Major differences, however, did occur in the relative importance of infauna versus epifauna. Except at the lagoonal site 9, where the small biofilm-feeding cushion star Parvulastra exigua occurs in large numbers, the intertidal zone of the whole axial channel apart from the upper estuary is dominated by infaunal species (Fig. 5a), principally by polychaetes. From sites 1 to 15, the infauna comprised 68% (SE 3.7) of animals with no significant trend in their relative importance along the gradient (S q = 0.31; P = 0.26) (Fig. 5B). In contrast, the shores of the upper estuary and the marine backwater channels were dominated by epifaunal truncatelloid microgastropods, especially by 'Assiminea' capensis and Hydrobia knysnaensis, epifauna here comprising 64.4% of individuals. Only a few subtidal Z. capensis sites have so far been examined, but such areas are also overwhelmingly dominated by an epifaunal microgastropod, here by the cerithioid Alaba pinnae, although the importance of this species and hence of the subtidal epifauna in general decreases upstream so that epifauna and infauna contribute equally in the upper estuary ( Fig. 5a). Thus in the bay region there is a transition at some LWS between a burrowing polychaete infauna and a seagrass-leaf-associated gastropod epifauna, and although upstream sub-and intertidal faunas are relatively similar, downstream in the bay they are markedly different (Fig. 6a). Few data are available to help explain the great downstream subtidal abundance of the epifaunal Alaba (a mean density of 28,000 m -2 ), although various studies have suggested that few fish consume significant numbers of shelled gastropods, even relatively small ones (McCormick 1998;Reynolds et al. 2018), not least because of their low nutritive value per unit intake (Vinson and Baker 2008). It is known that in South Africa, mugilids will  (Pollard et al. 2017). Several equivalent subtidal areas of seagrass in other continents are also dominated by species of Alaba, although Knysna is the only known such locality outside the tropics (Barnes and Claassens 2020). These other areas are of relatively high salinity which may help to account for the lesser importance of this gastropod in and near the upper estuary. Why the same suite of truncatelloid microgastropods dominates the otherwise contrasting habitats of the intertidal backwaters and upper estuary is not known for certain, but their common shelter (see paragraph above) is likely to be an important component. With one exception, no evidence of any strong species interactions within any given site was forthcoming. The exception was the positive correlation between numbers of the ectoparasitic pyramidellid snail Sayella sp. and those of its probable host Hydrobia at the backwater site 'A' in Fig. 1 (R s-= 0.78; P \ 0.00001). Such a parasite/host association is known from the western Atlantic (e.g., Hershler and Davis 1980), but although the pyramidellid concerned is a widely distributed animal, it is otherwise not recorded from Africa (GBIF 2020). That exception apart, however, in a large sample (325 cores) from the Kingfisher Creek seagrass (site 2 in Fig. 1), for example, Barnes (2013b) recorded 75 macrofaunal species at overall and mean densities of 2581 and 34 m -2 , respectively. Considering the 34 relatively common species there that each attained a mean density of at least 10 m -2 (and together comprised 96% of the total individuals), all pairwise correlations of species abundance were very weak to non-existent (sensu Moore et al. 2018), positives averaging only P q = 0.069 (± 0.060 SD) and negatives P q = 0.047 (± 0.035 SD); and allowing for the familywise errors inherent in such a large correlation matrix (via Bonferroni correction), no negative correlations and only three positive ones were significant at a critical a of \ 0.05 (between the polychaetes Simplisetia and Caulleriella, Glycera and Cirriformia, and between the polychaete Prionospio and the gastropod Nassarius). Equivalently, although qualitative co-occurrence patterns across the whole of the marine-influenced embayment at Knysna show deterministic structuring (Barnes and Elwood 2011), as indeed might be expected granted the location of the sampled sites in three distinct faunal regions (sandy mouth, marine bay, and backwater system), syntopic species within a single one of those regions did not differ from random co-occurrences (Barnes and Barnes 2014b).
In the absence of strong bioturbators such as Kraussillichirus kraussi (Callianassa kraussi in the older literature) from the majority of the system, faunal relationships between seagrass and bare sediment at Knysna are not the classic one of seagrass supporting the greater number of species and of individuals per unit area (Hemminga and Duarte 2000; Pillay et al. 2011;Hyman et al. 2019, etc.). To date studies have only concerned the outer marine embayment, but there seagrass macrofauna at a given site is more similar to those occurring in adjacent areas of bare sediment than either habitat is to other areas of the same type in the general region [Bray-Curtis faunal similarity between the two contiguous habitat types being a mean 0.58, whereas within-habitat-type similarity averaged 0.26 for the seagrass and 0.25 for the bare sediment (ANOVA F 1,14 = 5.05; P \ 0.05)] (see   Fig. 6b). In general, seagrass beds supported lower, not higher, levels in half the metric comparisons in which there was a significant difference (Barnes and Barnes 2014a). Overall, faunal abundance was lower in seagrass in the ratio of 0.64: 1, while species density was indeed higher, but only by 1.13 to 1, with in large measure the higher numbers in the unvegetated sediments resulting from a quadrupled abundance of infaunal polychaetes, maybe because of the greater volume of available sedimentary habitat in the absence of eelgrass rootmass, although numbers of epifaunal crustaceans were 15 times less there (from a much smaller base). The same overall effect was not the case, however, in bare areas created by the death of seagrass following blanketing by the chlorophyte blooms described by Allanson et al. (2016) and Human et al. (2016). In these circumstances, the former seagrass sites clustered together, as did the same areas when de-vegetated, although macrofaunal abundance was again significantly lower in the former seagrass than it was in the replacement bare sediment (in a ratio of 0.62: 1) and again largely because of an increased number of polychaetes and decreased number of crustaceans in the unvegetated sediment (Barnes 2019a).
Knysna's marine embayment forms a natural harbor, has been in the past a busy port (Grindley 1985), and today supports several marinas, and hence it is one of the centers of ship-borne alien immigrant species in South Africa (Griffiths et al. 2009). Alien species of Boccardia, Polydora, Dipolydora, Pseudopolydora, Diopatra, Capitella, Desdemona, Ericthonius, Jassa, Monocorophium, Paracerceis, Elysia, Favorinus and Indothais all form part of its seagrass fauna, as do amphipods such as Cymadusa filosa, Melita zeylanica and Americorophium triaeonyx that are regarded by Robinson et al. (2005) and Mead et al.
(2011) as being cryptogenic-to which could presumably be added Victoriopisa chilkensis. Relatively recently, these aliens have been joined by more northerly species spreading southward probably as a result of global warming. Smaragdia souverbiana, for example, is now a member of the subtidal seagrass fauna (Barnes and Claassens 2020). In the Knysna intertidal, Melanoides tuberculata has arrived and joined Cerithidea decollata (Hodgson and Dickens 2012) and Austruca occidentalis (formerly Uca annulipes) (Peer et al. 2015), the latter two in the adjacent saltmarsh or at the seagrass/saltmarsh interface.

Patterns in assemblage metrics along the axial gradient
As would be expected, the number of species at the 17 sites that were spaced along the system's longitudinal axis decreased with distance upstream (S q = -0.82; P \ 0.0001; Fig. 7a), but the form of the decrease in species density suggests the occurrence of a step change within the general area of the lower estuary, with the downstream sites showing a considerable degree of uniformity of species density ( Fig. 7b; Table 4). The points in Fig. 7 are based on the whole available 12-year dataset, and hence, the location of faunal and regional boundaries will have been blurred showed an almost identical (and sharper) feature (Barnes and Ellwood 2012) in the same general location. Comparison of data across different spatial scales shows that decline upstream in number of species when assessed per site (Fig. 7) is greater than when assessed per region (Table 4): Clearly, the bay and lagoon ? lower-estuary regions are large and, as also would be expected, the total fauna contained in each is considerably in excess of that at any individual site.
Assemblage abundance per unit area (S q = -0.35; P = 0.16; Fig. 8a) and patchiness in assemblage abundance (S q = -0.20; P = 0.45; Fig. 8b), however, showed no significant change with distance upstream; indeed, degree of patchiness along the axial gradient was significantly unchanging (Barnes 2019b). Neither Fig. 6 nMDS plots based on Bray-Curtis similarities between percentage species composition data: a at four sites in the Knysna estuarine bay, one in each faunal compartment, showing the similarity between macrofaunal assemblages of seagrass at LWS in the intertidal zone and in the adjacent subtidal area (from data in Barnes and Claassens 2020), and b at three sites in the marine embayment forming a transect from the main channel (site 1) into the backwater Steenbok Channel and its tributaries, showing similarity between macrofaunal assemblages in seagrass and in adjacent areas of bare sediment (from data in Barnes and Barnes 2014a). Envelopes enclose sites at the stated levels of Bray-Curtis similarity, and site codes are those given in Fig. 1, plus in (A) an additional site, 'X', located between sites 1 and 2 Fig. 7 Change in number of seagrass macrofaunal species along the longitudinal axis of the estuarine bay a. Note in b the apparent break between sites 10 and 11, corresponding to the one in the same area shown by Barnes and Ellwood (2012) was there any significant relationship between number of species per site and overall assemblage abundance there (S q = 0.43; P = 0.08). However, significant relationships have been found between how patchy an individual species is and its occupancy and, to a lesser extent, its abundance: The more abundant and widespread the species, the less its patchiness, both in subtidal and in intertidal seagrass (Barnes 2019c, 2020), and both in interspecific comparisons (Barnes 2020) and intraspecifically (Barnes, in prep.) (Fig. 9). This suggests that the well-known macroecological abundance-occupancy pattern (e.g., He and Gaston 2003) can be extended into a patchinessabundance-occupancy one, at least in this habitat type. As can be seen in Fig. 9, the slopes of the power laws relating logit occupancy to log patchiness in individual species are much more variable than those interspecifically in the different faunal regions; thus, while the interspecific occupancy-patchiness slopes representing different regions do not differ (ANCOVA F = 1.3, P = 0.3), the equivalent intraspecific slopes are heterogeneous (ANCOVA F = 4.9, P \ 0.0001) with a further six of the dominants (including the epifaunal Alaba and Cymadusa, and infaunal Caulleriella and Salmacoma) not showing significant occupancypatchiness relationships at all. This also indicates that disparate species together form assemblages with similar properties in the various regions. There were no discernable trends in either metric upstream, although the upper estuary did display the largest value of both b and R 2 .

General discussion and conclusions
The most striking feature of the Knysna intertidal seagrass-associated macrobenthos is its relative spatial uniformity. Macrofaunal abundance does not vary markedly along the longitudinal axis, neither does patchiness of macrofaunal density. Number of species per unit area at a given site is a constant, while species density along whole sections of the gradient can be relatively uniform, and its fauna appears to form a single assemblage with only local variation in relative frequency of its dominant components. Even the subtidal fauna does not differ qualitatively from the intertidal one, although there are marked quantitative differences and overall it is much more abundant especially in the most marine-influenced regions where Alaba dominates (Barnes and Claassens 2020). Admittedly, being a marine-dominated system with under normal circumstances relatively little freshwater input (Day et al. 1951), for an estuarine system salinity is relatively constant; however, occasional episodes of severe freshwater flooding do occur once every 10-12 years or so, rendering most or all of the system temporarily fresh (see, e.g., Korringa 1956; Blake and Chimboza 2010). But, although sea water may penetrate far upstream and dominate most areas, there is much change along the Knysna axis in other features of direct relevance to macrobenthos, as in transitional paralic systems in general (Tagliapietra et al. 2012;Pérez-Ruzafa et al. 2019). Sediment changes from clean sand at the mouth, to soft organic mud in the lagoonal and lower estuarine regions, and Table 4 Biodiversity metrics of the various intertidal faunal regions of Knysna estuarine bay: Mao tau species density (N 0 ), with Chao 2 estimations; N 1 and N 2 species diversity; equitability of species abundance (J); taxonomic diversity (D) and distinctness (D * ); and N 2 functional diversity (F d ) [with functional categories as per Macdonald et al. (2010) and Barnes and Hendy (2015a)]. Each axial regional metric is based on the common sample size of 180 9 0.0054 m 2 cores, yielding some 90% of the likely total species; backwaters metric based on 148 9 0.0054 m 2 samples. Peak values in bold  et al. (2020) as that displaying peak intensity of animal food consumption and hence potential top-down control of prey species. Like other South African estuarine areas supporting dwarf-eelgrass (Whitfield et al. 1989;Nel et al. 2018), it is a nursery area for many nektonic species (Whitfield and Kok 1992), Fig. 8 Relative constancy of intertidal seagrass macrofaunal assemblage metrics along the longitudinal axis of the estuarine bay: a assemblage abundance, with an inset showing abundances at the six backwater-channel sites; and b(i) assemblage patchiness (data from Barnes 2019b). b(ii) illustrates spatial variation in macrofaunal density 0.01 m -2 across an area of site 2 (from data in Barnes 2016) schools of juvenile fish being a common sight in its seagrass beds. Indeed, seagrass beds have been regarded as one the most important types of coastal nursery (Whitfield 2017; Lefcheck et al. 2019), both because of the food they provide (Whitfield 2017) and as refuge from larger fish predators (Whitfield 2020a). If, as is generally held to be the case in seagrass (Moksnes et al. 2008;Lewis and Anderson 2012;Duffey et al. 2015, etc.), the effect of this consumption, together with that exerted by adult fish (Pollard et al. 2017) and predatory members of the invertebrate macrobenthos, is top-down control of the seagrass microphytobenthic-biofilm consumers that dominate both epifauna and infauna, then many of features of Knysna's seagrass macrobenthos become more understandable.
Three lines of evidence suggest that across the whole system seagrass macrofaunal abundance is below carrying capacity and not structured by densitydependent factors. First, the prevailing intertidal density along the longitudinal axis of some 4,000 m -2 is very low compared to the [ 40,000 m 2 animals occurring in similar intertidal dwarf-eelgrass beds in cool-temperate Europe (Blanchet et al. 2004;Barnes and Ellwood 2011;etc.) where predator rates are almost certainly lower on a fauna of similarly sized animals that are often members of the same families as represented at Knysna (Barnes and Hendy 2015b). Secondly, constancy of number of species per unit area, as demonstrated at the Kingfisher Creek site at Knysna (Barnes 2013b), is exactly what would be expected were the various species to be distributed independently of each other (granted their overall frequencies of occurrence) (Barnes and Barnes 2014b). Such independence of distribution is likely only if the whole assemblage is being maintained below the level at which species would otherwise interact. Thirdly, the large quantitative dataset of Barnes (2013b) from the same Kingfisher Creek site also showed that there were very few significant correlations (0.5%) between the abundances of pairs of species and all those were very weak. Moreover, as seen elsewhere, for example within Pacific Canadian Zostera marina meadows (Stark et al. 2020), weak and very weak positive relationships greatly out-numbered negative ones, further suggesting the lack of competitive interspecific interactions. In such overall circumstances of low and unpredictable density of potential invertebrate prey species, predators could thus be expected to have to forage optimally (Beseres and Feller 2007) and to graze down local prey stocks to threshold levels before moving and repeating the process elsewhere, and having reduced their food stocks to low levels over wide areas could themselves then experience food limitation (Saulnier et al. 2000). The actions of such predators roaming widely over the surface might also help explain the uniform levels of macrofaunal patchiness characterizing large areas (Barnes and Hamylton 2019).
A second process that could help account for the observed features of the Knysna macrobenthic assemblages, and particularly their relative uniformity along the longitudinal axis, is if the presence of the seagrass Within the relatively uniform and stable conditions provided by the Z. capensis bed, species can penetrate upstream further than they might otherwise be able to do (Barnes and Ellwood 2012). Seagrass beds and other structurally complex systems (Hyman et al. 2019) generally appear to support macrofaunal assemblages that display spatial and temporal stability of abundance and composition, and high levels of resilience (Whanpetch et al. 2010;Blake et al. 2014;Gartner et al. 2015), and this may have important knock-on effects on the whole local coastal food web (Jankowska et al. 2019).
Granted the current high loss rates of seagrass throughout the world (Waycott et al. 2009;Short et al. 2011) and of Zostera capensis in southern Africa (Adams 2016), an important question is the implication of the operation of such potential structuring factors at Knysna for the future of its important seagrass system. Threats to the health of Knysna estuarine bay have recently been reviewed by Claassens et al. (2020), and the two that might particularly affect the seagrass and its inhabitants are uncontrolled bait harvesting by destructive means and chlorophyte blooms. Bait harvesting by 'pushing' for mudprawns (Upogebia capensis) and trenching for worms (Marphysa, Polybrachiorhynchus and Arenicola spp.) is rife in Knysna (Simon et al. 2019), even (and arguably especially) in the formally protected bait-reserve area (Fig. 10) that also supports 79% of the estuarine bay's seagrass-associated species. Unfortunately, there are immense logistic and social problems associated with preventing illegal and restricting legal bait harvesting in southern Africa, especially in areas of high unemployment where subsistence fishing provides the main or only source of protein (Bandeira and Gell 2003;Napier et al. 2009) and the local intertidal provides the only source of bait (Barnes and Claassens 2020). Subsistence harvesting at Knysna is worth some ZAR 1 9 10 6 (Turpie 2007). Its precise effect on the seagrass macrofauna is unknown; however, basically because the extent of harvesting means that a like-with-like, harvested versus unharvested, comparison is not possible: No area of intertidal seagrass remains un-pushed or un-pumped. At other localities, however, it and the associated trampling are known to have severe consequences (e.g., Pillay et al. 2010;Garmendia et al. 2017;Short et al. 2011;Adams 2016). Nevertheless, the inherent resilience and spatial uniformity of the Knysna seagrass beds referred to above, together with the fact that they form one large interconnected system with dispersal possible between all sections, does offer hope.
In respect of the second major threat, at the moment the problems of eutrophication, to which Z. capensis is known to be sensitive (Mvungi and Pillay 2019), and consequent algal blooms are only local, affecting mainly the backwater channels into which the municipal sewage treatment plant discharges and in which there is a legacy of organic matter retention (Human et al. 2020). Knysna's large tidal prism proves invaluable insofar as minimizing blooms in the main axial channel is concerned. But the local effect of the chlorophyte blanket is dramatic and destructive. Animal numbers may increase on its seasonal dieback, except in the very local areas of anoxia, although this increase is almost entirely confined to densities of the dominant polychaete groups (except, for some reason, cirratulids) (Barnes 2019a). Crustaceans do not bounce back so readily, and as they provide the food for most fish in the nearby Swartvlei Estuary (Whitfield 1988, 2020b) and presumably therefore do so at Knysna, the fish populations may suffer the consequences: There is some evidence that this is indeed the case (Pollard et al. 2018). The chlorophyte problem is soluble (no pun intended) but at considerable cost (see Human et al. 2020). The financial reward, however, might also be large, not least because such blooms impact on tourism (Boesch et al. 1996) and tourists contribute some ZAR 1 9 10 9 per annum to the Knysna economy (Turpie 2007).
Knysna is perhaps not a typical South African estuary in being permanently open (see van Niekerk et al. 2020), without a significant presence of Kraussillichirus, and with a subtidal dominated by the littleknown Alaba pinnae, but nevertheless its fauna generally appears to be the classic Cape estuarine one (Day 1981;de Villiers et al. 1999); indeed, it supports [ 40% of South African estuarine biodiversity. It also seems a microcosm of the country's estuarine problems, both environmental and socioeconomic, including loss of area to reclamation ([ 30% of the saltmarsh has already gone, and 25% of that remaining is under threat). Thus far, at least, it has managed to fare better than many other South African estuarine systems (van Niekerk et al. 2013), and it may yet continue to do so.
Acknowledgements I am most grateful to SANParks area office, Knysna, and SANParks Scientific Services, Rondevlei, for permission to conduct research on the Knysna system over the last decade or so; to Rhodes University Research Committee for their support; to my collaborators at various times, sequentially Farnon Ellwood, Morvan Barnes, Ian Hendy and Louw Claassens, for their contribution to the study; to Brian Allanson for his copious encouragement and support; and to an anonymous reviewer for saving me from error.
Author contributions RSKB conceived, designed, and executed this study and wrote the manuscript. No other person is entitled to authorship.
Funding No funding was received for conducting this study.
Data availability and material Datasets on which these analyses are based have been lodged in electronic format in the Rondevlei Office of SANParks Scientific Services (http:// dataknp.sanparks.org/sanparks/metacat/Nerinak.23.11/ sanparks) and are available on request.

Compliance with ethical standards
Conflict of interest The author declares no conflict of interest.

Consent for publication
The author consents to publication.
Ethical approval All applicable international, national, and/ or institutional guidelines for sampling, care, and experimental use of organisms for the study were followed, and all necessary permissions and approvals were obtained in respect of the original collections of the data.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.