Research Article |
Corresponding author: Matthew S. Bird ( mbird@uj.ac.za ) Academic editor: Christopher Majka
© 2023 Matthew S. Bird, David T. Bilton, Musa C. Mlambo, Renzo Perissinotto.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Bird MS, Bilton DT, Mlambo MC, Perissinotto R (2023) Water beetles (Coleoptera) associated with Afrotemperate Forest patches in the Garden Route National Park, South Africa. ZooKeys 1182: 237-258. https://doi.org/10.3897/zookeys.1182.102866
|
Southern Afrotemperate Forest is concentrated in the southern Cape region of South Africa and whilst it is relatively well known botanically, the fauna, specifically the aquatic invertebrate fauna, is poorly documented. The majority of remaining intact forest habitat is contained within the Garden Route National Park (GRNP), which straddles the provincial boundary between the Western and Eastern Cape. This study undertakes a survey of the water beetle fauna inhabiting the GRNP. The aquatic ecosystems within temperate forests of the region are poorly researched from an ecological and biodiversity perspective, despite being known to harbour endemic invertebrate elements. We collected water beetles and in situ physico-chemical data from a total of 31 waterbodies across the park over two seasons (summer and late winter) in 2017. The waterbodies sampled were mostly small freshwater perennial streams and isolated forest ponds. A total of 61 beetle taxa was recorded (29 Adephaga, 32 Polyphaga) from these waterbodies. The water beetle fauna of these forests appears to be diverse and contains many species endemic to the fynbos-dominated Cape Floristic Region, but very few of the species appear to be forest specialists. This is in contrast to the fynbos heathland habitat of the region, which harbours a high number of water beetle species endemic to this habitat, often with Gondwanan affinity. Our study is the first to document the water beetles of Afrotemperate Forests in the southern Cape region and provides an important baseline for future work on such habitats in the region and in other parts of southern Africa.
aquatic Coleoptera, aquatic invertebrates, biodiversity census, forest conservation, freshwater biodiversity, southern Cape, temperate forests
Closed-canopy evergreen indigenous forest is a relatively scarce biome in southern Africa, most of this vegetation type in South Africa being located in the east and north of the country (
Our understanding of the biodiversity of Afrotemperate Forests in southern Africa remains limited and patchy, both taxonomically and geographically. Floristic composition and endemism are relatively well understood for most groups; Southern Afrotemperate forests being dominated by palaeoendemic trees such as the podocarps (Afrocarpus and Podocarpus) and Cunonia, Ocotea, and Olea, but with relatively few strictly endemic plant taxa (
Very little work has been conducted on the freshwaters of Southern Afrotemperate Forests to date.
Samples were collected from a total of 31 waterbodies spanning the length of the Garden Route National Park (GRNP) along the southern Cape coast of South Africa between the towns of Storms River in the east and George in the west (Fig.
Location of the study sites along the southern Cape coastline of South Africa. Sites were grouped according to five main areas that represent different forest fragments within the Garden Route National Park (GRNP): Storms River; Nature’s Valley; Harkerville; Diepwalle and Wilderness (a). The position of the study sites in relation to the remaining core Afrotemperate Forest habitat in the region is depicted, as well the boundaries of the GRNP (b).
Sample collection was performed across two seasons, with 20 sites being sampled in early February 2017 (mid-summer, hereafter ‘summer’) and 14 sites sampled in mid-September 2017 (late winter, hereafter ‘winter’). Three of the sites sampled in summer were sampled again in winter. In terms of waterbody types sampled in this study, 15 of the sites sampled in summer were small perennial streams, whilst four such sites were sampled in winter. Five sites sampled in summer were ponds, whilst nine ponds were sampled in winter. One of the sites sampled in winter was a seepage wetland. The waterbodies were low lying, all occurring at less than 400 m altitude. Several clusters of sites were sampled, with 13 sites occurring in the vicinity of Storms River, 12 sites at Nature’s Valley, three at Harkerville, two at Diepwalle and one site was sampled at Wilderness. All sampled sites occurred within patches of Southern Afrotemperate Forest and were in a relatively pristine condition, being located inside the GRNP. The site locality information for all sampled waterbodies is provided in Table
Examples of the waterbodies and habitat types sampled in the Garden Route National Park (pictures taken during the summer survey in mid-February 2017). a Stream in the Plaatbos forest, Storms River (site 4) b main channel of the Storms River where the bridge crosses near the public picnic site (site 6) c the main channel of the Groot River at Nature’s Valley (site 13) d marshy pond at Plaatbos, Storms River (site 3) e stream crossing a hiking trail in the Plaatbos forest, Storms River (site 2) f small stream on the Kalanderkloof hiking trail, Nature’s Valley (site 9) g DTB examining water beetles at a pond in the Harkerville forest (site 15) h pond in the Diepwalle forest (site 19) i the authors hard at work sampling a pond at Nature’s Valley (site 12) adjacent to the Groot River j typical Southern Afrotemperate Forest habitat at Harkerville.
Site locality information for waterbodies sampled during this study. Two collection trips were undertaken, the first being during February 2017 (mid-summer) and the second during September 2017 (late winter). SR: Storms River; NV: Nature’s Valley; HV: Harkerville; DW: Diepwalle; WN: Wilderness.
Site | Date sampled (dd/mm/yyyy) | GPS (DD) | Altitude (m) | Region | Waterbody type | February (summer) | September (winter) |
---|---|---|---|---|---|---|---|
1 | 07/02/2017 | -34.02138, 23.88472 | 41 | SR | Stream | X | |
2 | 07/02/2017; 14/09/2017 | -33.97711, 23.89476 | 237 | SR | Stream | X | X |
3 | 07/02/2017; 14/09/2017 | -33.97541, 23.90689 | 239 | SR | Pond | X | X |
4 | 07/02/2017 | -33.98300, 23.90829 | 195 | SR | Stream | X | |
5 | 07/02/2017 | -33.98180, 23.91132 | 233 | SR | Stream | X | |
6 | 07/02/2017 | -33.98871, 23.91929 | 78 | SR | Stream | X | |
7 | 08/02/2017 | -33.97638, 23.88886 | 228 | SR | Stream | X | |
8 | 08/02/2017 | -33.97800, 23.88846 | 224 | SR | Stream | X | |
9 | 09/02/2017 | -33.97403, 23.55288 | 64 | NV | Stream | X | |
10 | 09/02/2017 | -33.97158, 23.54332 | 137 | NV | Stream | X | |
11 | 09/02/2017 | -33.96859, 23.55978 | 3 | NV | Stream | X | |
12 | 09/02/2017; 15/09/2017 | -33.96860, 23.55861 | 9 | NV | Pond | X | X |
13 | 09/02/2017 | -33.97605, 23.56169 | 2 | NV | Pond | X | |
14 | 10/02/2017 | -33.97428, 23.51926 | 39 | NV | Stream | X | |
15 | 11/02/2017 | -34.05024, 23.22491 | 240 | HV | Pond | X | |
16 | 11/02/2017 | -34.07092, 23.20679 | 189 | HV | Stream | X | |
17 | 11/02/2017 | -34.07839, 23.22742 | 177 | HV | Stream | X | |
18 | 11/02/2017 | -33.96131, 23.15123 | 392 | DW | Stream | X | |
19 | 11/02/2017 | -33.96436, 23.14399 | 381 | DW | Pond | X | |
20 | 12/02/2017 | -33.98355, 22.65148 | 5 | WN | Stream | X | |
21 | 14/09/2017 | -34.02140, 23.8886 | 25 | SR | Stream | X | |
22 | 15/09/2017 | -34.01701, 23.88892 | 101 | SR | Stream | X | |
23 | 15/09/2017 | -34.02209, 23.89196 | 68 | SR | Pond | X | |
24 | 15/09/2017 | -33.98311, 23.90889 | 195 | SR | Seep | X | |
25 | 15/09/2017 | -33.97967, 23.90582 | 217 | SR | Stream | X | |
26 | 15/09/2017 | -33.96713, 23.56006 | 3 | NV | Pond | X | |
27 | 16/09/2017 | -33.96937, 23.53168 | 218 | NV | Pond | X | |
28 | 16/09/2017 | -33.96966, 23.52587 | 223 | NV | Pond | X | |
29 | 16/09/2017 | -33.97414, 23.52207 | 36 | NV | Pond | X | |
30 | 16/09/2017 | -33.97509, 23.52778 | 87 | NV | Pond | X | |
31 | 16/09/2017 | -33.98408, 23.53546 | 4 | NV | Pond | X |
Water beetles were collected during both seasons using sweep netting. A long-handled square-framed pond net with a 30-cm mouth and 1-mm mesh was used for this purpose. With each sweep the net would be swept from the water surface to the bottom substrate and back to the surface again, in similar fashion to the protocols of
To provide baseline information on the freshwater habitats of GRNP, and an environmental context for the water beetle assemblages, basic in situ physico-chemical parameters were measured at each site. Temperature, conductivity, pH, turbidity, and dissolved oxygen were recorded using a YSI 6600-V2 multi-system probe. Physico-chemical measurements could not be taken from two of the sites during the summer survey due to logistical constraints.
All identifications were conducted by DTB, using a wide range of literature and, in some cases, comparison with reference/voucher material. All identifications were based, at least in part, on the study of male genitalia, unless otherwise stated.
Spatio-temporal patterns in the physico-chemistry of the waterbodies were assessed to determine whether the beetle assemblages mirrored physico-chemical patterns. Differences in physico-chemistry amongst sampled waterbodies were depicted using Principal Components Analysis (PCA), after first normalising the variables in the matrix. The variables constituting each matrix were temperature, conductivity, dissolved oxygen, pH, depth, and turbidity. Physico-chemical differences were compared across three factors of interest, which were overlaid on the PCA plots: season (summer, winter); region (Storms River, Nature’s Valley, Harkerville, Diepwalle, Wilderness); and waterbody type (streams, ponds, seeps). Permutational MANOVA (PERMANOVA,
Spatio-temporal patterns in beetle assemblage composition were depicted using non-metric multidimensional scaling (MDS). The MDS plots were overlaid by the same factors as per the physico-chemical data (seasons, regions, and waterbody types). Beetle presence-absence data were converted to a Bray-Curtis dissimilarity matrix in order to construct the MDS plots. PERMANOVA was used to test for differences in beetle assemblage composition (represented by a Bray-Curtis dissimilarity matrix) between the two seasons sampled (February 2017 – mid-summer vs. September 2017 – late winter) and amongst the different regions (i.e., separate forest patches) sampled along the Tsitsikamma coast (Storms River, Nature’s Valley, Harkerville, Diepwalle), as well as between waterbody types (streams vs. ponds). For regional comparisons using PERMANOVA, Wilderness was not included in the test due to only one site being sampled in that region, and similarly seeps were excluded from the comparisons of waterbody types due to only one seep being sampled. Species richness (number of species recorded per waterbody) was similarly compared amongst seasons, regions, and waterbody types. Richness patterns were visually assessed using boxplots and comparisons between seasons, regions and waterbody types were performed using t-tests (two-group comparisons) or one-way ANOVA (three-group comparisons), given that the richness data followed a Gaussian distribution and significant heteroscedasticity was not evident for any of the comparisons (
All tests were performed using an a priori significance level of α = 0.05. P values for PERMANOVA models were tested using 999 unrestricted permutations of the raw data. The PCA, MDS, PERMANOVA and DISTLM procedures were implemented with PRIMER v. 7.0.21 software (
The waterbodies encountered in the forests of the GRNP were predominantly small perennial rocky streams, although a small proportion of these streams (e.g., sites 9 and 10) are expected to dry up intermittently. There are several larger running waters in the park, such as the Groot, Storms, and Salt rivers, which were also sampled in this study. The second-most abundant waterbody type encountered was ponds (or depression wetlands according to the South African wetland classification system;
Table
Physico-chemical variables recorded at each waterbody during the February and September 2017 surveys. Median, minimum, and maximum values are reported for each survey. Readings were not recorded at sites 13 and 20.
Survey date | Site | Temperature (°C) | Conductivity (mS.cm-1) | pH | Dissolved O2 (mg.L-1) | Turbidity (NTU) | Depth (m) |
---|---|---|---|---|---|---|---|
February 2017 | 1 | 20.41 | 0.980 | 7.80 | 8.45 | 1.6 | 0.10 |
2 | 18.54 | 0.166 | 4.64 | 8.19 | 0.5 | 0.16 | |
3 | 20.41 | 0.221 | 6.41 | 2.53 | 8.0 | 0.17 | |
4 | 18.71 | 0.175 | 5.65 | 7.63 | 0.0 | 0.17 | |
5 | 18.64 | 0.163 | 5.31 | 2.35 | 1.4 | 0.21 | |
6 | 20.61 | 0.100 | 4.49 | 8.76 | 1.8 | 0.17 | |
7 | 19.95 | 1.138 | 6.74 | 7.53 | 1.8 | 0.10 | |
8 | 18.51 | 0.350 | 7.07 | 3.48 | 0.7 | 0.09 | |
9 | 17.51 | 0.811 | 6.91 | 6.22 | 0.2 | 0.16 | |
10 | 20.38 | 0.425 | 7.01 | 5.68 | 9.6 | 0.12 | |
11 | 21.44 | 0.178 | 6.89 | 3.83 | 13.2 | 0.11 | |
12 | 19.57 | 0.757 | 7.02 | 5.90 | 52.5 | 0.23 | |
14 | 20.93 | 0.250 | 5.51 | 9.09 | 210.0 | 0.07 | |
15 | 18.90 | 0.240 | 6.70 | 0.56 | 12.0 | 0.22 | |
16 | 18.50 | 0.373 | 6.87 | 8.82 | 0.0 | 0.60 | |
17 | 18.71 | 0.415 | 6.78 | 7.56 | 0.8 | 0.18 | |
18 | 17.52 | 0.193 | 6.35 | 3.31 | 11.6 | 0.30 | |
19 | 19.38 | 0.311 | 6.81 | 0.72 | 86.5 | 0.45 | |
Median | 19.14 | 0.281 | 6.76 | 6.06 | 1.8 | 0.17 | |
Minimum | 17.51 | 0.100 | 4.49 | 0.56 | 0.0 | 0.07 | |
Maximum | 21.44 | 1.138 | 7.80 | 9.09 | 210.0 | 0.60 | |
September 2017 | 2 | 13.48 | 0.955 | 8.16 | 9.10 | 2.1 | 0.30 |
3 | 12.08 | 0.140 | 8.53 | 1.82 | 7.1 | 0.07 | |
12 | 14.93 | 0.802 | 7.32 | 1.53 | 12.9 | 0.05 | |
21 | 13.92 | 0.761 | 9.70 | 9.75 | 2.2 | 0.45 | |
22 | 13.74 | 1.724 | 8.18 | 2.46 | 28.6 | 0.60 | |
23 | 14.52 | 2.245 | 7.70 | 0.83 | 431.0 | 0.50 | |
24 | 12.72 | 0.454 | 8.15 | 3.41 | 12.5 | 0.05 | |
25 | 14.24 | 0.138 | 7.15 | 7.83 | 5.1 | 0.05 | |
26 | 14.54 | 4.605 | 7.78 | 5.93 | 26.2 | 0.05 | |
27 | 14.20 | 0.159 | 7.64 | 7.71 | 20.1 | 0.15 | |
28 | 15.28 | 0.334 | 6.70 | 6.51 | 5.2 | 0.07 | |
29 | 15.00 | 0.359 | 9.65 | 3.17 | 28.4 | 0.15 | |
30 | 17.54 | 0.279 | 8.77 | 6.88 | 31.5 | 0.08 | |
31 | 15.63 | 0.370 | 7.98 | 4.88 | 107.0 | 0.07 | |
Median | 14.37 | 0.412 | 8.06 | 5.41 | 16.5 | 0.08 | |
Minimum | 12.08 | 0.138 | 6.70 | 0.83 | 2.1 | 0.05 | |
Maximum | 17.54 | 4.605 | 9.70 | 9.75 | 431.0 | 0.60 |
Sites were generally shallow, being < 1 m in depth (the deepest recording was 0.60 m for sites 16 and 22). However, this does not reflect the true depth of some of the larger rivers such as the Groot River, where water beetles were targeted in the shallow marginal vegetation at the edges (e.g., site 13) rather than the deeper middle section of the channel. Although some of the waterbodies were well oxygenated (dissolved oxygen concentrations > 7 mg.L-1), a large proportion of the sites had low dissolved oxygen concentrations, with some of the streams and ponds recording remarkably low values (< 2 mg.L-1, see Table
According to the PCA plot in Fig.
Principal components analysis depicting multivariate differences in the physico-chemistry of the various waterbodies sampled in this study. The site numbers are indicated above the symbol for each site and symbols have been differentiated according to season (summer vs. winter trips). The physico-chemical variables measured at each site have been overlaid as vectors on the plot.
Non-parametric permutational MANOVA (PERMANOVA) results for models comparing the physico-chemistry of the waterbodies between (a) seasons, (b) regions and (c) waterbody types. The multivariate models tested for differences between group centroids in multivariate space, represented by Euclidean distance. An asterisk indicates significant P values at α = 0.05.
(a) | df | SS | MS | F | P |
---|---|---|---|---|---|
Season | 1 | 45.64 | 45.64 | 9.75 | 0.001* |
Residual | 30 | 140.35 | 4.67 | – | – |
Total | 31 | 186 | – | – | – |
(b) | df | SS | MS | F | P |
Region | 3 | 12.85 | 4.28 | 0.69 | 0.762 |
Residual | 28 | 173.15 | 6.18 | – | – |
Total | 31 | 186 | – | – | |
(c) | df | SS | MS | F | P |
Waterbody type | 1 | 18.344 | 18.34 | 3.26 | 0.005* |
Residual | 29 | 163.04 | 5.62 | – | – |
Total | 30 | 181.39 | – | – | – |
The list of water beetles recorded in this study is reported in Table
Water beetles collected from the Garden Route National Park during the course of this study. The sites are listed from which each taxon was collected on each of the two sampling trips (February and September 2017). Site numbers 1 – 31 correspond to those listed in Table
Taxa | Sampling date | Region | |||||
---|---|---|---|---|---|---|---|
February | September | SR | NV | HV | DW | WN | |
Gyrinidae: | |||||||
Dineutus grossus (Modeer, 1776) | 1, 6, 8, 14 | 23 | X | X | |||
+Aulonogyrus formosus knysnanus Brinck, 1955 | 13, 14, 16, 17 | X | X | ||||
Aulonogyrus varians Brinck, 1955 | 6 | 25 | X | ||||
+Orectogyrus capicola Brinck, 1955 | 14 | X | |||||
Haliplidae: | |||||||
+Haliplus exsecratus Guignot, 1936 | 11, 20 | X | X | ||||
Noteridae: | |||||||
Synchortus simplex Sharp, 1882 | 3 | X | |||||
Dytiscidae: | |||||||
+Agabus austellus Englund, Bilton & Bergsten, 2020 | 15 | X | |||||
+Copelatus caffer Balfour-Browne, 1939 | 1, 2, 3, 4, 10, 11, 12, 13, 15, 18, 19 | 2, 3, 12, 21, 23, 24, 25, 26, 30 | X | X | X | X | |
+Copelatus capensis Sharp, 1882 | 1, 3, 9, 10, 11, 12, 13, 15, 18 | 3, 12, 21, 23, 24, 26, 28, 30 | X | X | X | X | |
Copelatus erichsoni Guérin-Méneville, 1849 | 10, 11, 12 | 3, 12, 23, 24, 30 | X | X | |||
+Copelatus notius Omer-Cooper, 1965 | 11 | X | |||||
Aethionectes apicalis (Boheman, 1848) | 12 | 12 | X | ||||
Hydaticus capicola Aubé, 1838 | 10, 11, 13, 14, 15, 17, 18, 20 | 12, 23, 27 | X | X | X | X | X |
Hydaticus dregei Aubé, 1838 | 8 | X | |||||
Hydaticus galla Guérin-Méneville, 1849 | 1, 3, 4, 5, 7, 8, 10, 11, 12, 13, 14, 15, 16, 18, 19 | 3, 12, 21, 22, 23, 24, 27, 30 | X | X | X | X | |
+Bidessus mundulus Omer-Cooper, 1965 | 28 | X | |||||
Clypeodytes meridionalis Régimbart, 1895 | 1, 3, 6, 8, 13, | 25 | X | X | |||
Hydroglyphus lineolatus (Boheman, 1848) | 27 | X | |||||
Uvarus opacus (Gschwendtner, 1935) | 3 | X | |||||
Yola frontalis Régimbart, 1906 | 4, 6, 8, 11, 4 | X | X | ||||
+Canthyporus fluviatilis Omer-Cooper, 1956 | 3, 15 | X | X | ||||
+Canthyporus hottentottus (Gemminger & Harold, 1868) | 3, 8 | 26, 27 | X | X | |||
+Hydrovatus amplicornis Régimbart, 1895 | 3 | 28 | X | X | |||
+Darwinhydrus solidus Sharp, 1882 | 15 | 27, 28, 29, 31 | X | X | |||
+Hydropeplus trimaculatus (Laporte, 1835) | 27 | X | |||||
+Hyphydrus soni Biström, 1982 | 1, 3, 6, 7, 10, 11, 12, 13, 15, 16, 20 | 12, 22, 25, 26, 27 | X | X | X | X | |
+Africophilus jansei Omer-Cooper & Omer-Cooper, 1957 | 14 | X | |||||
Laccophilus lineatus Aubé, 1838 | 3, 6, 7, 11, 12, 13, 14, 20 | 22, 25 | X | X | X | ||
Hydrochidae | |||||||
Hydrochus sp. | 27, 28, 29 | X | |||||
Spercheidae | |||||||
Spercheus cerisyi Guérin-Méneville, 1842 | 27, 28 | X | |||||
Hydrophilidae | |||||||
Amphiops globus Erichson, 1843 | 1, 11, 12, 14, 19 | 12, 22 | X | X | X | ||
Amphiops senegalensis (Laporte, 1840) | 13 | X | |||||
+Anacaena capensis Hebauer, 1999 | 25 | X | |||||
+Anacaena glabriventris Komarek, 2004 | 10, 14 | 27 | X | ||||
Agraphydrus albescens (Régimbart, 1903) | 6, 13 | X | X | ||||
+Enochrus hartmanni Hebauer, 1998 | 27, 28, 29 | X | |||||
Enochrus (Methydrus) sp. | 1, 3, 4, 8, 9, 10, 12, 13, 15, 17, 19 | 23, 24, 30 | X | X | X | X | |
Helochares longipalpis (Murray, 1859) | 3 | X | |||||
Helochares sp. | 6 | X | |||||
+Limnoxenus sjoestedti Knisch, 1924 | 27 | X | |||||
Hydrochara elliptica (Fabricius, 1801) | 27 | X | |||||
Sternolophus mundus (Boheman, 1851) | 1, 11, 12, 13 | 12 | X | X | |||
Laccobius praecipuus Kuwert, 1890 | 14 | X | |||||
Coelostoma sp. | 14 | X | |||||
Hydraenidae | |||||||
Hydraena cooperi Balfour-Browne, 1954 | 3, 13 | X | X | ||||
+Mesoceration apicalum Perkins & Balfour-Browne, 1994 | 2, 4, 16, 17 | X | X | ||||
+Mesoceration barriotum Perkins, 2008 | 17 | X | |||||
+Mesoceration dissonum Perkins & Balfour-Browne, 1994 | 2, 4, 5 | X | |||||
+Mesoceration distinctum Perkins & Balfour-Browne, 1994 | 6 | X | |||||
+Mesoceration integrum Perkins, 2008 | 17 | X | |||||
+Nucleotops interceps Perkins, 2004 | 29 | X | |||||
+Parhydraena asperita Perkins, 2009 | 1, 2, 4, 15, 17 | 31 | X | X | X | ||
+Parhydraena seriata Balfour-Browne, 1954 | 22, 26, 29 | X | X | ||||
Dryopidae | |||||||
+Strina sp. | 6, 17 | 25 | X | X | |||
Elmidae | |||||||
Stenelmis sp. | 2 | X | |||||
+Elpidelmis capensis (Grouvelle, 1890) | 2, 4, 6, 17 | 25 | X | X | |||
+Elpidelmis fossicollis Delève, 1966 | 25 | X | |||||
+Peloriolus sp. 1 | 6 | X | |||||
+Peloriolus sp. 2 | 2, 5, 6, 17 | 25 | X | X | |||
Ptilodactylidae | |||||||
Ptilodactylidae (larvae) | 13 | X |
Mean taxon richness across all sites and sampling trips was 7.1±3.7 (±SD) taxa per site. The most taxa recorded at a single site was 14, recorded at sites 3 and 13, which were both ponds. This was followed by sites 6 (stream) and 27 (pond), where 13 taxa were collected at each of these sites. Therefore, three out of the four richest sites were ponds. The boxplots in Fig.
Boxplots comparing the median and spread of water beetle taxon richness (number of taxa per site) between a seasons b regions and c waterbody types at GRNP. The middle line represents the median, whilst the boxes demarcate the interquartile range and the ‘whiskers’ extend to the maximum and minimum values. The black circles on the graphs represent individual data points (number of taxa) for each site sampled. Unpaired t-tests reported no significant difference in richness between the two seasons (t32 = 1.604, p = 0.119) and between the waterbody types (streams vs. ponds, t31 = 0.959, p = 0.345). One-way ANOVA reported no significant difference in richness between the regions (F3,29 = 0.809, p = 0.499). ‘Seeps’ was excluded as a factor from the waterbody comparisons due to only one sample being taken from this habitat and ‘Wilderness’ was similarly excluded from the regional comparison due to only one sample being collected in this region.
Water beetle assemblage composition differed between seasons, regions, and waterbody types at GRNP, as depicted visually in the MDS plots in Fig.
Multidimensional scaling (MDS) plots depicting the similarity of sites sampled at GRNP in terms of their water beetle assemblages. Symbols on the plot have been coded in terms of a season b region and c waterbody type. Convex hulls (dashed lines) have been overlaid on each plot to clarify groupings according to season, region, or waterbody type.
Non-parametric permutational MANOVA (PERMANOVA) results for models comparing beetle assemblage composition across (a) seasons, (b) regions and (c) waterbody types. The multivariate models tested for differences between group centroids in Bray-Curtis dissimilarity space. SR: Storms River; NV: Nature’s Valley; HV: Harkerville; DW: Diepwalle. For the regional comparison, Wilderness was not included due to only one site being sampled there on one occasion and for the comparison of waterbody types, ‘seeps’ was excluded as a factor because only one seep was sampled on one occasion (i.e., streams were compared with ponds). An asterisk indicates significant P values at α = 0.05.
(a) | df | SS | MS | F | P | |||
Season | 1 | 6241.6 | 6241.6 | 2.20 | 0.018* | |||
Residual | 32 | 90669 | 2833.4 | – | – | |||
Total | 33 | 96910 | – | – | – | |||
(b) | df | SS | MS | F | P | Post hoc pairwise comparisons | ||
Groups | t | P | ||||||
Region | 3 | 12497 | 4165.6 | 1.51 | 0.048* | SR, NV | 1.477 | 0.021* |
Residual | 29 | 80134 | 2763.2 | – | – | SR, HV | 1.238 | 0.140 |
Total | 32 | 92630 | – | – | – | SR, DW | 0.840 | 0.661 |
NV, HV | 1.265 | 0.103 | ||||||
NV, DW | 0.979 | 0.329 | ||||||
HV, DW | 1.282 | 0.155 | ||||||
(c) | df | SS | MS | F | P | |||
Waterbody type | 1 | 6152.4 | 6152.4 | 2.136 | 0.033* | |||
Residual | 31 | 89270 | 2879.7 | – | – | |||
Total | 32 | 95422 | – | – | – |
The measured environmental variables in this study were together able to explain approximately 78.5% of the variation in beetle assemblage composition among the waterbodies sampled in the GRNP (Table
Results of the dbRDA multivariate regression tests of environmental variables against beetle assemblage composition. Independent marginal tests are first presented (a), followed by variables selected by the step-wise procedure using the AICc selection criterion (b) and the ‘best’ (most parsimonious, considering all combinations of variables) overall model according to the AICc criterion (c). ‘% Var’: the percentage of variation in each Bray-Curtis similarity matrix that is explained by the respective predictor variable in each test; ‘Cum. % var’: the cumulative percentage variation across all tests; ‘Res. df’: residual degrees of freedom associated with each test. An asterisk indicates significant variables at α = 0.05.
(a) Marginal tests: | |||||
Variable | F | P | % Var | ||
Latitude | 1.24 | 0.240 | 3.99 | ||
Longitude | 1.20 | 0.294 | 3.86 | ||
Season | 2.01 | 0.036* | 6.31 | ||
Altitude | 1.05 | 0.405 | 3.40 | ||
Region: ‘Storms River’ | 1.83 | 0.065 | 5.77 | ||
Region: ‘Nature’s Valley’ | 2.31 | 0.019* | 7.16 | ||
Region: ‘Harkerville’ | 1.48 | 0.136 | 4.72 | ||
Region: ‘Diepwalle’ | 0.70 | 0.691 | 2.29 | ||
Waterbody type: ‘Stream’ | 2.21 | 0.021* | 6.89 | ||
Waterbody type: ‘Pond’ | 2.01 | 0.042* | 6.29 | ||
Waterbody type: ‘Seep’ | 0.45 | 0.842 | 1.49 | ||
Temperature | 1.85 | 0.053 | 5.84 | ||
Conductivity | 1.17 | 0.33 | 3.78 | ||
DO | 1.34 | 0.198 | 4.29 | ||
pH | 2.64 | 0.006* | 8.11 | ||
Depth | 0.68 | 0.717 | 2.24 | ||
Turbidity | 0.66 | 0.702 | 2.18 | ||
Total: | 78.59 | ||||
(b) Sequential tests: | |||||
Variable | AICc | F | P | % Var. | Res. df |
pH | 256.1 | 2.64 | 0.01* | 8.11 | 30 |
(c) Best solution: | |||||
Variable | AICc | F | P | % Var. | Res. df |
pH | 256.1 | 2.64 | 0.01* | 8.11 | 30 |
Our work demonstrates that the waterbodies of forests of the Garden Route National Park support a diverse water beetle fauna, including a number of South African endemics. The total of 61 taxa recorded from the region is, however, considerably lower than the 116 reported from similar surveys by the same team in the subtropical iSimangaliso Wetland Park, further north on the KwaZulu-Natal coast (
Of the species recorded here, 32 are endemic to South Africa. The vast majority of these are Cape endemics, more widespread in the fynbos biome to the west, and not tied to forest waterbodies. Such species include the dytiscids Darwinhydrus solidus Sharp, 1884 and Hydropeplus trimaculatus (Laporte, 1835), both of which are widespread and often abundant in lentic waters in fynbos in the far southwestern Cape, a number of the stream-dwelling Mesoceration (Hydraenidae) found in GRNP and the two lotic Elpidelmis species (Elmidae). Very few water beetle species found in these forests are either local endemics or forest specialists, the suite of taxa recorded during our surveys being dominated by species typical of fynbos waterbodies of the southern Cape (DTB, pers. obs.). Taxa which appear to be genuinely restricted to this region are Aulonogyrus formosus knysnanus Brinck, 1955 (Gyrinidae) and Parhydraena asperita Perkins, 2009 (Hydraenidae). Of these two, only the latter appears to be predominantly a forest species, which is particularly abundant in the margins of small standing waters filled with decaying leaf litter, although it has also been reported from stream margins in the nearby Little Karoo (
Our study demonstrates that there are clear, measurable, differences between the aquatic beetle assemblages in different forested sections of the Garden Route National Park, as revealed by nMDS and PERMANOVA analyses, but no significant differences in species richness. Clearly, despite these forested catchments being close geographically, there is significant spatial variation in aquatic habitats, reflected in the different beetle faunas. Interestingly, the relatively few environmental parameters recorded during our study are able to explain almost 79% of the variation in beetle assemblage composition across sites, suggesting that these measures capture the main environmental drivers of species composition in the region. In most studies, even with many more environmental parameters, the proportion of explained variation is typically much lower (e.g.,
In summary, our study documents the aquatic beetle faunas of southern Cape Afrotemperate Forests for the first time, providing an important baseline for future work in the area and similar habitats in other parts of southern Africa. We show that these systems support a wide range of water beetle species, including a number of South African endemics, but do not, apparently, harbour any truly local endemics, even in running waters. This observation is in marked contrast to streams draining fynbos catchments, particularly further west in the Cape, where high concentrations of locally endemic water beetles are known, many with Gondwanan affinities. Whilst
Sampling in the Garden Route National Park was authorised by a research collection permit from South African National Parks (permit number: PERRI-R/2016-023). We are very grateful to Nerina Kruger (Science Liaison & GIS at South African National Parks Scientific Services, Sedgefield) and the managers of the various sections of the GRNP for providing general assistance with the project and logistic support during the two expeditions of this study. We would like to thank Lynette Clennell for assisting with field collections. Albrecht Komarek (Vienna) kindly examined the Agraphydrus species. This work is based on the research supported by the South African Research Chairs Initiative of the Department of Science and Technology (DST) and National Research Foundation (NRF) of South Africa. Any opinion, finding and conclusion or recommendation expressed in this material is that of the author(s) and the NRF does not accept any liability in this regard.
No conflict of interest was declared.
No ethical statement was reported.
National Research Foundation (NRF); Department of Science and Technology (DST).
Conceptualization: DTTB, MCM, MSB, RP. Data curation: MCM, DTTB. Formal analysis: MSB, MCM. Funding acquisition: RP. Investigation: DTTB, MCM, MSB, RP. Methodology: DTTB, MCM, RP, MSB. Project administration: RP. Resources: DTTB, RP. Software: MSB. Validation: MCM, DTTB. Visualization: MSB. Writing – original draft: MCM, MSB, RP, DTTB. Writing – review and editing: RP, DTTB, MSB, MCM.
Matthew S. Bird https://orcid.org/0000-0001-9163-1008
David T. Bilton https://orcid.org/0000-0003-1136-0848
Musa C. Mlambo https://orcid.org/0000-0001-7624-5686
Renzo Perissinotto https://orcid.org/0000-0002-9224-3573
All of the data that support the findings of this study are available in the main text or Supplementary Information.