A floristic comparison of three Northern Coastal Forests differing in disturbance history

Forests are closed-canopy plant communities comprising mainly woody plants more than 5 m tall (Midgley et al. 1997). Forest is the smallest and most widely fragmented biome in South Africa (Low & Rebelo 1996; Rutherford & Westfall 1994). In South Africa, forests generally occur as patches within other biomes. It is estimated that indigenous coastal forest patches cover less than 500 000 ha (Low & Rebelo 1996). Forest along the eastern coast of South Africa represents the southernmost distribution of tropical and subtropical forest fauna and flora in Africa, and contains a high level of floristic endemism, with a number of species that are at the end of their distribution range. Indigenous coastal subtropical forest types are embedded within the broader Indian Ocean Coastal Belt (IOCB) biome (Mucina & Geldenhuys 2006).


Introduction
Forests are closed-canopy plant communities comprising mainly woody plants more than 5 m tall (Midgley et al. 1997).Forest is the smallest and most widely fragmented biome in South Africa (Low & Rebelo 1996;Rutherford & Westfall 1994).In South Africa, forests generally occur as patches within other biomes.It is estimated that indigenous coastal forest patches cover less than 500 000 ha (Low & Rebelo 1996).Forest along the eastern coast of South Africa represents the southernmost distribution of tropical and subtropical forest fauna and flora in Africa, and contains a high level of floristic endemism, with a number of species that are at the end of their distribution range.Indigenous coastal subtropical forest types are embedded within the broader Indian Ocean Coastal Belt (IOCB) biome (Mucina & Geldenhuys 2006).
The narrow subtropical coastal strip of KwaZulu-Natal (KZN) contains several diverse forest types, which are patchy in distribution, but harbour relatively high biodiversity despite their limited spatial extent (Moll & White 1978;Mucina & Rutherford 2006).Previous studies have shown forest in KZN to have decreased in extent since the early 1800s owing to anthropogenic pressures such as agriculture, forestry, urbanisation, tourism and mining (Lawes et al. 2007;Midgley et al. 1997;Mucina & Geldenhuys 2006).However, more recent assessments suggest that a few conserved patches of natural evergreen forest may have increased over the last 30 years (Environmental Geographic Information Systems [EGIS] 2016).This increase can be ascribed to variable combinations of the following: urban development, controlled access and protection against fire.Coastal forests in KZN typically occur as a series of small patches (< 10 ha) (Moll & White 1978;Mucina & Rutherford 2006), accounting for only approximately 9.6% of the total land area (Rutherford et al. 2006;Von Maltitz et al. 2003).Nevertheless, coastal forests in KZN still contribute to a substantial proportion of the plant biodiversity housed in the Maputaland-Pondoland-Albany (MPA) biodiversity hotspot (Steenkamp et al. 2004).These forests are diverse and of varying types based on vegetation composition, structure and altitude (Mucina & Geldenhuys 2006), and may have been variably linked and separated owing to Quaternary glacial cycles (Griffiths & Lawes 2006).More recent anthropogenic pressures (prior to European settlement) have also influenced the forests either directly (e.g.harvesting of wood for iron production) or indirectly (e.g.altering fire regimes) (Armstrong et al. 2000;Morgenthal & Cilliers 1999).The most recent changes (post 1800s), which include high levels of urbanisation and agriculture, have had the most significant effects on the extent and integrity of these coastal forests.Historical transformation and conservation activities (e.g.rehabilitation and restoration) have collectively led to a mixture of primary and secondary forest types.
Forests are prone to both natural (e.g.fire, wind, temperature fluctuation, precipitation variability, alluvial processes, etc.) and anthropogenic (e.g.harvesting and burning) disturbances (Eeley, Lawes & Piper 1999;Fahrig 2003) that can act as exogenous agents of vegetation change (White 1979).Consequently, forest flora are rich in disturbanceadapted species and disturbance-recovery processes that follow both natural and anthropogenic disturbances, contribute to floristic and structural composition, and underlie species diversity patterns in the long term (White 1979).These disturbance-recovery processes are influenced by the type, degree and duration of disturbances, as well as the natural regenerative potential of the species within the system, which includes soil seed bank composition (Bossuyt, Heyn & Hermy 2002;Chazdon 2003).Disturbances (anthropogenic and natural) that remove large areas of tree canopies, for example, can result in excessive increases in sunlight reaching the forest floor and an associated decrease in soil moisture content.While the increased light availability can provide favourable germination conditions for lightdemanding species, the decline in soil moisture content can compromise seedling survival.Recovery following disturbances that heavily impact soils and aboveground vegetation can be considerably slow, often with long-lasting effects on species composition (Wassenaar et al. 2005).
Forest floristic composition, demography and succession stage are therefore reflections of past disturbance events and the system's ability to recover from them (Adema et al. 2005;Van Aarde & Claassens 1998).Floristic surveys are essential for biodiversity assessment, management and conservation (Hahs et al. 2009).An appraisal of current species composition in remnant patches of coastal forest in relation to their disturbance history can therefore inform their future conservation or rehabilitation.The availability of basic coastal forest floristic data (e.g.species checklists) can also help advance broader, pressing research questions related to demography, phenology, succession, alien invasions and climate change, amongst others.This motivated the present study that compared the floristics of three Northern Coastal Forest (NCF) patches in relation to their contrasting disturbance histories.
Northern Coastal Forest (designated as FOz-7 in the forest classification developed by Mucina and Geldenhuys, 2006) is composed of two forest types, viz., KwaZulu-Natal Coastal Forest and KwaZulu-Natal Dune Forest (sensu Von Maltitz et al. 2003).In the present study, we examined KwaZulu-Natal Coastal Forest specifically, which is known to be species rich and occurs as patches within the MPA (Mucina & Geldenhuys 2006).Studies that examine rehabilitation, restoration and regeneration of NCF or any other coastal subtropical forest types in South Africa are scarce (Olivier & Van Aarde 2014;Van Aarde et al. 2008).The spatial, temporal, biotic and abiotic requirements for these processes to be successful are generally unknown (Arnold & Lutzoni 2007).The three NCF patches selected for investigation here are geographically separated (by > 100 km) and included Hawaan, which is largely primary forest with a small area of regrowth (c.80 years old) that was not sampled in this study and has been subjected to low levels of disturbance and protection for c.50 years; Umdoni, which is composed of core parts of primary forest (Mucina & Geldenhuys 2006) and patches of new growth (10-30 years old) and has been exposed to high levels of disturbance for a prolonged period of time; and a secondary forest (Twinstreams) that was rehabilitated (partly with Sand Forest canopy tree species) (Moll & Cooper 1966;Untiedt 1992).The disturbance and successional histories of these forests deserve further consideration and are expanded on below.At the time of this study, there were no recent published reports on floristics of the forests examined here.Previous studies conducted at Hawaan Forest (HAF) and Twinstreams Forest (TMF) (Moll 1972;Moll & White 1978) represent a comprehensive description of the flora at that time, but are now outdated.No published data (except for an informal tree and species list) exist for Umdoni Forest (UMF).
The primary objective of the present study was to assess the impacts of disturbance on species composition and diversity in NCF.Classical vegetation survey methods were used to assess floristic composition and diversity within each forest, which was then used to assess floristic similarity amongst the three patches in relation to their contrasting disturbance histories.Important issues relate to whether a secondary forest rehabilitated with atypical species will result in these species becoming excluded from the system or successfully becoming integrated into a novel but artificial (structurally or floristically) system as reported elsewhere (Lugo & Helmer 2004;Uriarte et al. 2009).In addition, an assessment of the similarity, in terms of species composition and diversity, between UMF (high levels of disturbance) and HAF (low levels of disturbance) allowed for an investigation of the resilience of the NCF floristic signature.

Study area
Northern Coastal Forest is found embedded in the IOCB biome that extends from the Eastern Cape northwards into KZN (Mucina & Rutherford 2006).In some areas, these forests are restricted to the dunes (Acocks 1988).Beyond South Africa, these forests occur throughout the Mozambican coastline as far as southern Tanzania.Many tropical elements reach their southern limit in the NCF of South Africa (Mucina et al. 2006).The present study was conducted in three NCF patches along the Indian Ocean coastline (Figure 1).Hawaan Forest (29°05'S, 31°05'E) is located in the Umhlanga Rocks area (northern Durban), while the secondary forest (TSF) (30°23'S, 030°40'E) is situated near the town of Mtunzini, 160 km north of Durban.Umdoni Forest, the southernmost of these forests (28°59'S, 31°44'E), is located between Park Rynie and Ifafa Beach.The largest forest is UMF (200 ha), followed by HAF (55 ha) and TSF (40 ha).The details of all sites are given in Table 1.To justify the comparisons made in this study, it is instructive to consider the successional status and disturbance histories of the three forests in greater detail.
A forest with no recorded logging and which has developed under natural processes following natural disturbance is considered 'primary' (Cockle, Martin & Drever 2010).Hawaan is a primary forest that is fully fenced and access controlled for c.50 years.A narrow fringe (< 500 m wide) of secondary forest occurs along the western perimeter which is a result of this area being farmed > 80 years ago and abandoned owing to poor soil quality.Over the course of c. 80 years, this fringe has not recovered to climax forest and is clearly distinguishable (from aerial photos and in situ) from the primary forest given the low, sparse canopy and abundance of pioneer species (e.g.Albizia adianthifolia).For the purposes of this study, we only sampled the primary forest at HAF.It should also be noted that frequent clearing of alien species has taken place at this forest for many years.In light of the sampling approach adopted here (i.e.no sampling in regrowth portions), HAF was considered a primary forest with low levels of disturbance in this study.
Umdoni is composed of core parts of primary forest (Mucina & Geldenhuys 2006) that have experienced high levels of disturbance owing to the neighbouring golf course and associated activities.There are also parts of the forest that have expanded over time.These parts have been subject to uncontrolled burning and clearing, and are clearly distinguishable from the primary forest.Unlike HAF, parts of UMF are subjected to illegal access, informal pathways, clearing, illegal harvesting, erosion and only occasional alien plant removal.The anthropogenic disturbances mentioned above have not been restricted to specific parts of the forest; however, it was possible to identify the physiognomically mature part of the forest and this is the area that was sampled in the present study.With the above in place, UMF was regarded as a primary forest with high levels of disturbance.
The secondary forest (TSF) was rehabilitated from a sugar cane farm c. 70 years ago and is presently exposed to low levels of disturbance.Classical definitions of a secondary forest suggest that it is a forest that has been logged and has subsequently recovered naturally or artificially (Huang et al. 2017).The forest is presumed to be expanding out of the demarcated forest boundaries, but as in UMF, sampling at TSF was confined to the physiognomically mature part of the forest.

Climate and altitude
The coastal belt of KZN has a seasonally moist subtropical climate.Rainfall is usually in excess of 1 000 mm per annum, becoming drier inland with less than 600 mm per annum (Kirkwood & Midgley 1999;Mucina & Rutherford 2006).The highest monthly precipitation falls between September and April, resulting in hot, humid summers and cool, dry winters.Altitude is the main factor determining both the rainfall and temperature, and ranges from 10 m to 160 m above sea level (a.s.l.) with the topography being relatively flat (Table 1).All NCFs selected in this study occur at low altitudes (10 m -160 m a.s.l.).

Vegetation
Floristically, coastal KZN is diverse and complex with many different vegetation types (Mucina & Rutherford 2006;Scott-Shaw & Escott 2011).More recently, all NCF patches were placed within the IOCB biome, but these forests are also a distinct vegetation type of the forest biome (Von Maltitz et al. 2003;Mucina & Rutherford 2006).Northern Coastal Forest is described as species-rich, tall to medium-height subtropical coastal forest on coastal (rolling) plains that stabilise the sand dunes on which they occur (Moll & White 1978;Mucina & Rutherford 2006;Mucina et al. 2006).These forests are dense, structured into several tree layers and exhibit a well-developed shrub layer (Mucina & Geldenhuys 2006).Interestingly, in False Bay there is Eastern Sand Forest ecotonal to what may be NCF (Kirkwood & Midgley 1999).

Data collection and analyses
Studies investigating effects of disturbance within forests generally adopt a core to edge sampling approach (e.g.Matlack 1994) in order to compare vegetation across different parts of the same forest (intra-site comparisons).However, where the  effects of disturbances on vegetation are to be compared across different forests (inter-site comparisons), quadrats are most often distributed across the physiognomically mature forest portions to accommodate the effects of disturbance on the system as a whole (Denslow 1995;Frelich & Lorimer 1991;Gentry 1982).Fifteen quadrats were randomly laid out within the physiognomically mature part of the forest, more than 100 m apart from each other and any edge, and sampled between June 2014 and June 2015, in each of the four seasons.
Quadrat sizes for forest vegetation surveys vary greatly across studies, ranging from 10 m × 10 m (Matthews et al. 2001) to 50 m × 50 m (Zhang & Cao 1995) and even greater (Gentry 1982), depending on the site characteristics and research question(s).For the purpose of this study, it was imperative to reach a sampling effort of ≥ 80%, irrespective of quadrat size.Given the fact that three forests had to be sampled over multiple seasons, a quadrat size of 15 m × 15 m was deemed manageable and repeated sufficiently at all sites to achieve the acceptable sampling effort.
Cover (%) and abundance were recorded for all non-graminoid taxa within each quadrat, while only cover was recorded for graminoids, given their clonal habit.To produce a comprehensive species checklist, species occurring outside of the quadrats were also recorded (but not in terms of cover and abundance) using transects once every three months (i.e.within each season).A series of ten 2 m × 100 m transects were laid out in directions that avoided trails and clearings.
The transects were roughly parallel to each other and separated by 20 m or more.
Flowering or fruiting material was collected for voucher specimens (deposited in the Ward Herbarium [UDW], University of KwaZulu-Natal).Specimens were identified in the herbarium using field guides (e.g.(WIP 2006).In some cases (especially naturalised exotics), the alien status was determined using SANBI (2018b).Aliens not assessed in the above literature and websites were treated as non-categorised alien taxa.
Sampling effort was based on rarefaction curves.As the number of species is highly dependent on sample size, rarefaction curves (Gotelli & Colwell 2001) were plotted in EstimateS (Version 9.1.0)(Colwell 2013) using quadrat data of all three forests to determine whether sampling effort was sufficient (Chao et al. 2006(Chao et al. , 2009)).Graminoids (viz.grasses) were excluded from the above-mentioned analyses because abundance data for grasses were not available.Two nonparametric estimators (Chao1 and Chao2) were used to estimate total species richness.The percentage sampling effort was calculated using the number of species found in situ and the projected number of species based on the above two estimators.

Floristic similarity
In the present study, cluster and ordination analyses were applied to quadrat (n = 45) data to assess similarity amongst the three forest sites (Albernaz et al. 2012).The presence/ absence data were used to generate a similarity matrix based on the Bray-Curtis index.Clustering was performed with the Unweighted Pair-Group Method using the Arithmetic Average (UPGMA) in PAST v3.The package PAST v3 was also used to conduct a non-metric multidimensional scaling (NMDS) ordination in two dimensions, with the Bray-Curtis index applied to quadrat abundance data (n = 45).The ordination was evaluated using the 'stress' value (a measure of the mismatch between instance measures and the distance in ordination space).Stress values smaller than 20 generally lead to usable interpretations (Jackson et al. 2012;Kruskal 1964).
The ordination was conducted with the abundance data with default settings (two-dimensional solution, maximum number of random starts = 100).Grasses were excluded from the NMDS, but not from the cluster analyses.Previous studies (e.g.Jackson et al. 2012) have adopted a similar approach.However, the groups identified by the UPGMA analysis were compared with the ordination plot to gain a better understanding of similarity amongst the study sites.

Sampling effort and floristics
A total of 152, 177 and 179 species were found in quadrats at HAF (primary forest), UMF (disturbed primary forest) and TSF (secondary forest), respectively.Sampling effort based on the Chao1 and Chao2 estimators ranged from 76.7% to 83.2%.Values > 80% are generally regarded as an indication of adequate sampling (e.g.Moro, De Sousa & Matias 2012), and because the relatively high values obtained here (76.7% -83.2%) were based exclusively on quadrat data, this suggested that sampling was sufficient.When transect, quadrat and existing species list data were combined, the lowest number of species were found at the undisturbed primary forest (HAF; n = 312), followed by TSF (n = 390) and UMF (n = 490) (Figure 2).Collectively, 688 species belonging to 400 genera and 113 families were identified across all three forest sites, when all the data were considered (Online Appendix 1).A total of 151 species (c.21.95%) were common (i.e.shared) to the three forests (Figure 2).Eighty-three species were shared between the two primary forests (HAF and UMF), 24 between HAF and TSF and 95 between UMF and TSF (Figure 2).The primary forest (HAF) exhibited the lowest number of unique (i.e.unshared) taxa (n = 54), while 161 and 120 unique taxa were found at UMF and TSF, respectively.These dominant families contributed to 50.73% (349 species) of the total species found across the three forests.The 14 most dominant genera in terms of species richness were:

Taxa of conservation concern
Twenty-eight species (4.07% of the total number of species found across the three forests) are of conservation concern: seven Endangered, nine Vulnerable, seven Near Threatened, two Rare and three Declining (Table 2).The lowest number of species of conservation concern occurred at HAF (n = 8), while TSF and UMF had 17 and 18 of these species, respectively.The species of conservation concern were classified according to the SANBI National Red List (SANBI 2018b), which provides fairly detailed information on species habitats and vegetation types.All 28 species of conservation concern occur across more than one vegetation, but only 12 of these are considered typical forest elements (  (declining).Hawaan Forest and TSF shared no taxa of conservation concern, while HAF and UMF shared three species of conservation concern.Twinstreams Forest and UMF shared six species of conservation concern (Table 2).

Alien taxa
Alien taxa accounted for 8.72% (60 species) of the total number of species identified across the three forests (Online Appendix 2).These included Category 1-3 aliens (Category 1 [n = 31], Category 2 [n = 3], Category 3 [n = 2]), followed by not declared aliens (n = 8), naturalised exotics (n = 9) and non-categorised aliens (n = 7).In summary, there were 36 Category 1-3 aliens which are considered as invasive alien plant species (IAPs).The highest number of aliens were found in the disturbed primary forest (n = 43), followed by the secondary forest (n = 30).The primary relatively undisturbed forest (HAF) housed far fewer aliens (n = 17).A few aliens (e.g.Chromolaena odorata, Lantana camara, Litsea sebifera, Rivina humilis, Solanum seaforthianum, Achyranthes aspera and Carex sylvatica) were present in all three forests.Nine alien species were shared between the disturbed primary forest and the secondary forest, three between the two primary forests and just four between HAF and TSF (Online Appendix 2).

Diversity indices
The Shannon's evenness index was comparable across the three forests: 0.87, 0.80 and 0.86 for HAF, UMF and TSF, respectively (Table 3).These results suggest that all three forests have equivalent evenness.The Shannon diversity (H') index did, however, differ across forests, with values being highest for the primary forest (4.37), followed by the secondary forest (3.79) and disturbed primary forest (3.69).The Simpson's index of diversity (1-D) was almost the same for all three forests.Whittaker's β diversity (βw) index was highest for the disturbed primary forests (2.65), followed by the secondary forest (2.38), and relatively lower in the primary forest (1.84) (Table 3).

Floristic similarity
The phenogram of the UPGMA analysis shown in Figure 3 was based on the presence/absence data of 282 species (including grasses) found in the quadrats.There are three main clusters (A-C) at c. 20% similarity.Clusters A, B and C are composed of quadrats from HAF, TSF and UMF, respectively.Each cluster is forest specific, that is, contains  quadrats from a single forest site only.Clusters B and C grouped together, indicating that the disturbed primary forest (UMF) and secondary forest (TSF) are more similar to each other than to Cluster A, the primary forest (HAF).Furthermore, it is worth noting that the two primary forests are geographically more proximal to each other than either are to the secondary forest (TSF).
The NMDS analyses showed a separation of three forest specific groups with low similarity (Figure 4).The primary   forest (HAF) was clearly separated from (i.e.dissimilar to) the other two forests along axis 1 in the ordination space (stress = 22.97).These results suggest that UMF and TSF are more similar to each other than either are to HAF, a pattern also found in the UPGMA analysis above.These results are assumed to reflect the differences in species composition and abundances of the three forests.

Sampling approach
Some studies assessing the effects of disturbance on forests have sampled along a disturbance gradient (e.g.Matlack 1994), from core to edge, within individual forests.The major advantage of this approach is that it allows one to assess the effects of disturbance within parts of the same system.However, differences in species composition arising from core to edge sampling are largely a reflection of differences in vegetation age.Disturbances such as illegal harvesting and tree felling, and even natural disturbances, are also not confined to the edge.In fact, our observations suggest that in many forest patches this occurs within the core of the forest.
In addition, dispersers, plant propagules and alien plants move freely between disturbed and undisturbed portions of the same forest, making it a very fluid system.For these reasons, core to edge sampling may not allow one to assess the effects of disturbance on the system as a whole.Randomly positioning quadrats within the physiognomically mature part of the forest is a common approach when examining disturbance effects on forests (e.g.Denslow 1995; Frelich & Lorimer 1991) and may be an alternative way to capture heterogeneous distribution of disturbances within forests.Furthermore, this approach facilitates floristic comparisons across different forests.

Floristic composition and structure
Studies on other tropical and subtropical forests have shown species richness and diversity to be influenced by a range of biotic (e.g.competition, regeneration, growth, migration and succession [Newbold et al. 2014;Whittaker, Willis & Field 2001]) and abiotic (e.g.temperature, humidity, precipitation and soil moisture [Arnold & Lutzoni 2007;Pimm et al. 2014]) factors.In lowland coastal forests of KZN, succession can occur rapidly (30-40 years), and large-scale natural disturbances (e.g.periodic fires and herbivory) are essential for forest processes (Van Wyk et al. 1996).However, these coastal forests (including NCF) are also subjected to anthropogenic related pressures and disturbances (e.g.urban development, fragmentation, pollution and mining [Midgley et al. 1997]).Of the three forests compared here, the primary forest (HAF) has received the best protection over the last 50 years and has been maintained in a relatively undisturbed state (except for some parts of its margin) through fencing, controlled access and alien plant control.The disturbed primary forest (UMF) and the secondary forest (TSF) have more recently received some protection, but both are still open to human disturbance owing to poor controlled access and protective measures.Furthermore, as mentioned earlier, TSF is a secondary forest that was rehabilitated using a number of Sand Forest taxa.
When the quadrat data were considered, the number of species found in the three forests (152-179) was in the range reported for subtropical forests in eastern KZN (e.g.104 woody species [Gaugris & Van Rooyen 2008]); 195 tree species [Griffiths & Lawes 2006]).However, the number of species was lowest in the primary forest (HAF) and comparably higher in UMF and TSF, despite the latter being a secondary forest.When the quadrat data were supplemented with transect and existing species list data, the total number of species was again lowest in the undisturbed primary forest (HAF; n = 312) and higher in the secondary forest (TSF; n = 390) and disturbed primary forest (UMF; n = 490) (Figure 2).It was expected that UMF would have the highest species richness based on size (200 ha) (MacArthur & Wilson 1967).However, TSF, the smallest forest (40 ha), had far more species than HAF (55 ha) (Table 1).As in other studies (e.g.Pimm et al. 2014), these data suggest that habitat loss and disturbance influence forest species richness (i.e.alpha diversity).Studies on forests in others parts of the world have also shown that certain disturbances can lead to an increase in forest alpha diversity (Hobbs et al. 2006;Macdonald et al. 1989).
Despite the differences in alpha diversity noted above, the three forest patches exhibited a number of similarities in terms of dominant families.For instance, Fabaceae was the most species-rich family in TSF and ranked second in HAF and UMF (where Rubiaceae and Asteraceae were the most speciose families, respectively).The dominance of Fabaceae (8.58% of the total flora found) in all three forests is not surprising as it is amongst the most important subtropical forest families (Terborgh & Andresen 1998).Northern Coastal Forest is known to have nutrient-deficient, sandy soil substrates (Ceruti 1999) and legumes are particularly successful on nutrient-poor substrates (Dickinson 2003) owing to their ability to fix atmospheric nitrogen, and produce large seeds, which provides them with an advantage (Ter Steege et al. 2006).
Furthermore, the 13 most dominant families (≥ 13 species) in this study comprised 349 species (50.73% of the total flora).The Shannon's evenness index was comparable across the three forests irrespective of disturbance history.A threshold value of 2 for the Shannon diversity (H') index is considered as a minimum value, above which an ecosystem may be regarded as having medium to high diversity (Whittaker 1972).Values for this index for all forests studied were higher than 2 and in the range of values reported for other KZN coastal forests (2.3-3.6 [Olivier & Van Aarde 2014]).
Interestingly, values for this index were highest for the primary forest (4.37), followed by the secondary forest (3.79) and disturbed primary forest (3.69).Values for the Simpson's index of diversity index (1-D) obtained here were slightly higher than that reported for KZN coastal forests (0.58-0.88 [Olivier & Van Aarde 2014]).Unlike species richness, diversity indices take into account species abundance as well.Despite this, the patterns observed for diversity appear to contradict those observed for species richness generated for each forest via the compilation of quadrat, transect and species lists, which showed species richness to be the greatest in the secondary forest and lowest in the primary forest (Figure 2).A comparison of the quadrat data exclusively also showed species richness to be the lowest at HAF (n = 152) and higher at both UMF (n = 177) and TSF (n = 179) (Table 3).
This contradiction may be a manifestation of a disturbanceinduced increase in species richness at UMF and TSF.Disturbances associated with anthropogenic activities were notably higher at UMF and TSF (Kambaj, Sershen & Ramdhani 2016), which could have accommodated opportunistic species and, hence, increased richness as reported for forests elsewhere (Gibson et al. 2011;Macdonald et al. 1989).The explanation offered above is supported by the fact that alien plant species in UMF and TSF were higher than that in HAF.These aliens are known to be disturbance opportunists (Richardson et al. 2005), which contributed to the higher species richness within the disturbed forests.It is also worth noting that the secondary and disturbed primary forest exhibited higher numbers of unique species (121 and 161, respectively) than the primary forest (54) (Figure 2).The low number of unique flora, together with the relatively lower species richness at this forest, suggests that high levels of protection in HAF may have facilitated higher levels of stability (despite its size) compared with the disturbed patches.A number of unique species in the disturbed forests were also abundant (e.g.Tricalysia lanceolata and Oricia bachmannii in UMF; Microsorum scolopendria and Millettia grandis in TSF [data not shown]) and included aliens (Ageratum houstonianum, Rubus apetalus and Solanum mauritianum in UMF; Mangifera indica in TSF) as well as taxa of conservation concern (Ocotea bullata in TSF; Cryptocarya latifolia in UMF).
In addition to the above, βw values for the three forests (1.84 for HAF; 2.65 for UMF; and 2.38 for TSF; Table 3) suggest that these sites are dissimilar.Beta diversity is the spatial turnover or change in the identities of species.It measures the difference in species composition either between two or more local assemblages or between local and regional assemblages (Koleff, Gaston & Lennon 2003).A βw of 0 indicates that sample units have all the same species, while a high βw value indicates the opposite (Whittaker 1972).These results are supported by the fact that of the 688 species collectively found across all three forests, only 21.95% were shared by all three forests (Figure 2), with HAF sharing only 24 species with TSF.As alluded to earlier, high βw at UMF and TSF may be related to the disturbance and rehabilitation history of these two sites.
A shortcoming of the present study is that data on stem density and diameter data were not collected.These data would have allowed us to relate differences in species composition and diversity across the three forests to possible differences in vegetation age.Studies have shown species richness to be significantly correlated with stem density (Denslow 1995).However, Denslow (1995) also showed that the wellestablished relationship between abundance and diversity offers a suitable null model for the effects of disturbances on forest tree species diversity.Species abundance data were collected in the present study and differences in species composition and abundance across the three forests are possibly the product of variable combinations of disturbance, rehabilitation, protection and natural succession.Twinstreams Forest, for example, housed several species that are atypical of NCF and Sand Forest elements (e.g.Elaeodendron transvaalense and Newtonia hildebrandtii) (Boon 2010).In addition, Ocotea bullata and Sterculia murex were recorded in TSF.Ocotea bullata is a mist-belt forest element, while S. murex is naturally distributed in northern Swaziland and Mpumalanga province (Boon 2010).Their presence and more importantly persistence at TSF is based on the fact that this patch of NCF was reforested with Sand Forest and other elements.
The cluster analysis and NMDS plot reflect the site differences described above.Quadrats of TSF and UMF are placed closer, while HAF quadrats were more dissimilar to those of the disturbed sites (Figure 4).It is generally expected that fragmented landscape patches closer in geographic proximity have similar species compositions, that is, they are more similar (MacArthur & Wilson 1967).Similarities amongst forest patches have in the past been attributed to the regional species pool being subjected to the same selection pressures during glacial maxima (Griffiths & Lawes 2006).Species characteristic of NCF (Gaugris & Van Rooyen 2008;Mucina & Geldenhuys 2006) were present in both primary forests (e.g.Allophylus natalensis, Brachylaena discolour, Crotalaria natalensis, Dracaena aletriformis, Drypetes arguta, Ficus natalensis, Harpephyllum caffrum, Protorhus longifolia, Teclea natalensis and Zanthoxylum capense) and also in the secondary forest (TSF).However, HAF and UMF, which are closer in geographic proximity to each other, than either are to TSF, did not cluster together.An ecosystem's stability is based on two separate but interrelated properties: resistance and resilience (Lehman & Tilman 2000).Resilience refers to the speed of recovery from disturbance, while resistance is associated with an ecosystem's ability to retain structure and composition following a perturbation and/or to regenerate.The comparisons of the floristics between a primary forest (HAF) and primary disturbed (UMF) forest conducted here suggest that certain elements of the NCF signature may be more resistant to disturbance than others.
Disturbance and the presence of Sand Forest and other elements in TSF could account for its intermediate positioning in the ordination and suggest that TSF (and UMF to some extent) may represent novel ecosystems in which disturbance (and rehabilitation) has facilitated the introduction of species outside the NCF floristic signature.However, these disturbances and the introduction of new species have not led to a change in vegetation structure and exclusion of typical NCF taxa.Human disturbances can drastically affect ecosystem trajectories, and, in changing species composition and relative abundance, can create 'novel ecosystems', commonly referred to as 'emerging ecosystems' (Hobbs, Higgs & Harris 2009;Pimm et al. 2014).Novel or emerging ecosystems have not been reported for NCF previously, but have been reported for other forest types (Milberg & Lamont 1995).Communities that develop on these modified patches are generally persistent and characterised by new species combinations (e.g.Milberg & Lamont 1995), which may be the case at TSF (and UMF).

Conservation of Northern Coastal Forest
Northern Coastal Forests fall within the MPA biodiversity hotspot (Steenkamp et al. 2004), contributing to its richness and endemism.There are strong indications that NCF has harboured species of tropical forest origin (e.g.Baphia racemosa, Celtis africana, Cola natalensis and Vepris lanceolata) (Kirkwood & Midgley 1999;Moll & Cooper 1966) for most of its history, forming a relatively stable vegetation type in which extinction levels have been relatively low (Olivier & Van Aarde 2014).However, NCF is presently fragmented in eastern KZN and is under anthropogenic pressure and disturbance (from agriculture, illegal clearing, unsustainable harvesting of medicinal plants and alien plant invasions [Von Maltitz et al. 2003]) in various coastal urban matrices where remnant patches still occur (Mucina et al. 2006).These forests are naturally fragmented even without human influence.However, anthropogenic disturbances related to the built environment prevent expansion and reconnection as is the case in HAF, which is surrounded almost entirely by urban development.A few patches of NCF may be expanding (e.g.Twinstreams [Authors unpubl. observ.]),but the nature and extent are uncertain and require quantitative confirmation.
Conserving biodiversity in areas of high human density or activity is a major conservation challenge (Dearborn & Kark 2010).This is particularly apparent along the KZN coastline, where large-scale habitat transformation has led to the loss of approximately 50% of native vegetation (Midgley et al. 1997).
Northern Coastal Forest presently faces a number of threats, particularly on coastal dunes (Olivier & Van Aarde 2014).
Research conducted in the late 1990s indicated that anthropogenic activities led to the expansion of IAPs, which have, in turn, had negative effects on coastal KZN vegetation (Midgley et al. 1997).However, more recent studies have shown that facilitative interactions between IAPs and native species occur in a wide range of habitats, including forests (Rodriguez 2006).According to Rodriguez (2006), habitat modification (e.g.shading) is most often the mechanism through which IAPs aid the emergence and subsequent growth of native species, but a combination of different mechanisms (trophic subsidy, pollination, competitive release, predatory release, etc.) can also apply.This implies that IAPs may facilitate expansion of the forests investigated here, but this would require a focused study on alien-native plant interactions along forest edges.This was beyond the scope of the present study, but certainly a recommendation for future studies, given reports that the facilitative effects of alien species can aid forest restoration efforts by modifying habitats and promoting the establishment of native species (D'Antonio & Meyerson 2002).It should be noted, however, that a high number of IAPs were present within mature parts of all three forests investigated, including the one with low levels of disturbance (HAF).These species included Category 1-3 IAPs (Online Appendix 2), many of which are significant invaders in South Africa (Nel et al. 2004).Studies have shown that the general assumption that undisturbed forests are highly resistant to plant invasions may not hold under all circumstances; for example, 139 exotic plant species have been shown to invade relatively undisturbed deeply shaded forest understories (Martin, Canham & Marks 2009).These authors go on to state why exotics present a unique management challenge: (1) exotics often increase in abundance during succession; (2) their invasion rate may be comparatively slow under natural disturbance regimes; (3) their rate of invasion may be accelerated by anthropogenic processes, including the spread of exotic pests and pathogens.Many IAPs were reported to be prominent in NCFs as far back as 1985 (Cooper 1985) and our results suggest that many of these aliens have persisted while new ones have also emerged.
It was encouraging to note that several species of conservation concern were still present in the patches investigated.All 28 species of conservation concern occur across more than one vegetation type, while only 12 of these are considered typical forest elements (SANBI 2018b), which suggests that these species are not confined to NCF, but also occur in adjacent coastal vegetation types.The lack of species of conservation concern should reflect disturbance (i.e.fewer species of conservation concern in increasingly disturbed sites).However, in the present study, the lowest number of species of conservation concern occurred at the forest with the lowest levels of disturbance (HAF; n = 8).The high number of such species in TSF (n = 17) may be explained by the deliberate planting of these species during rehabilitation.The relatively large size of UMF (200 ha) could explain why it housed the highest number of species of conservation concern (n = 18).Endangered species were most prevalent (in terms of the number of species) in TSF followed by UMF and HAF, while UMF harboured the most Vulnerable and Declining species.Near Threatened species were also most prevalent in TSF, followed by UMF and HAF.An equal number of rare species were found in HAF and TSF, with UMF having no rare species.
These results emphasise the conservation value of both undisturbed and disturbed NCF patches, the benefits of rehabilitating transformed patches and the need to increase levels of protection within disturbed patches.The success of conservation plans for these NCF patches and other forests that occur within urban matrices depends on the combined and coordinated efforts of researchers, land owners, corporates, government, local communities and non-governmental stakeholders (Kambaj, Sershen & Ramdhani 2016).

Conclusion
The comparison of three NCF patches differing in disturbance history carried out here suggests that these forests all exhibit high levels of alpha and beta diversities.Disturbance and rehabilitation in the case of TSF has led to differences in both alpha and beta diversities across the forests, with the primary forest composition being distinct from the two disturbed forests.Disturbance and/or rehabilitation has also led to the presence of species atypical of NCF, many of which are aliens.
Nevertheless, the patches shared some similarities in species composition and physical structure, and a number of typical NCF taxa have persisted in all three forests.This together with the high number of species of conservation concern within all three patches emphasises their conservation value.Most importantly, the results suggest that even though disturbance and rehabilitation of NCF can lead to the creation of novel ecosystems, NCF signature taxa occurred at both the disturbed and restored sites.

FIGURE 1 :
FIGURE 1: Location of the three Northern Coastal Forest study sites in KwaZulu-Natal.

FIGURE 2 :
FIGURE 2: Venn diagram (based on quadrat, transect and species list data) showing species found within each forest site.

FIGURE 3 :
FIGURE 3: Phenogram generated from the Unweighted Pair-Group Method using the Arithmetic Average analyses of the three sites (45 quadrats).The letters (A-C) and the coloured bars denote the three Northern Coastal Forest study sites (A, white = Hawaan Forest [HAF]; B, grey = Twinstreams Forest [TSF]; C, black = Umdoni Forest [UMF]).

FIGURE 4 :
FIGURE 4: Non-parametric multidimensional scaling ordination plot of the three study sites.

TABLE 1 :
Summary of the features of the three Northern Coastal Forest sites.

Table 2
).Only three species of conservation concern were common to all forests: Dioscorea sylvatica (Vulnerable), Atalaya natalensis (Near Threatened) and Adenia gummiferaThe total number of species recorded at each site is indicated in parentheses.HAF, Hawaan Forest; UMF, Umdoni Forest; TSF, Twinstreams Forest.

TABLE 2 :
Taxa of conservation concern found in the three Northern Coastal Forest sites.
Note: The cross mark 'X' denotes the presence of the species in the forest while dash '-' denotes absence.HAF, Hawaan Forest; UMF, Umdoni Forest; TSF, Twinstreams Forest; SFE , Sand Forest Element.

TABLE 3 :
Species richness and various diversity indices of the three Northern Coastal Forest sites (all diversity indices based on quadrat data).
(Van Wyk & Smith 2001)rat sampling.The dominant families mentioned above are also well represented in the South African flora in general(Van Wyk & Smith 2001).
Orchidaceae, which is one of the largest Angiosperm families with c. 26 000 species globally (APW 2018), was poorly represented in this study (10 species, Online Appendix 1).Some epiphytic orchids (e.g.Cyrtorchis arcuata, Cyrtorchis praetermissa and Mystacidium capense) were recorded, but those which reside in the upper reaches of the canopy may http://www.abcjournal.orgOpen Access have been