Introduction

In temperate lowland rivers, larval and juvenile fish in their first year of life (‘0+’) may face episodic flood events that can be deleterious to their cohorts, especially in early summer when individuals are still in early developmental stages (Nunn et al., 2002, 2007a, b). Early developmental stages have not yet developed all fins and muscle structures to withstand flow, and in channelised rivers, seeking slackwaters may not be possible (Bolland et al. 2015). The probability of over-winter survival and recruitment of 0+ fishes can also be positively correlated to their body lengths at the end of their first growth season (Kirjasniemi & Valtonen, 1997; Mills & Mann, 1985; Nunn et al., 2003).

Spawning strategies of temperate riverine cypriniform fishes vary by species (Vila-Gispert, 2002), but where each strategy is assumed to maximise the survival rates of the 0+ fish to facilitate larger numbers to survive their first year of life and subsequently recruit (Beardsley & Britton, 2012a). In species such as dace Leuciscus leuciscus (Linnaeus 1758), spawning can be a single event, usually in early spring and when no other cypriniforms are reproducing (Gutmann Roberts, 2018). This maximises the access of adults to spawning substrates as well as ensuring the 0+ fish then have an extended growth season that facilitates individual fish to attain relatively large body lengths (e.g. > 50 mm) that should increase their probability of surviving the over-wintering period (Mann, 1974; Nunn et al., 2002; Beardsley & Britton, 2012b). Accessing favourable spawning habitats is enhanced by their relatively large spawning migrations where individuals often move upstream and into tributaries (Clough & Beaumont 1998; Kottelat & Freyhof 2007). This strategy, however, also means that a stochastic event in early summer (e.g. a flood) could result in high rates of mortality and/ or downstream displacement (Nunn et al., 2003, 2007a). An alternative strategy is the use of fractional or batch spawning events (hereafter referred to as ‘multiple spawning’) (Gutmann Roberts, 2018). Utilised by species such as chub Squalius cephalus (Linnaeus 1758), multiple spawning strategies appear to involve a trade-off between prolonged spawning efforts in the adult fish (that can potentially extent over two months) versus a reduced probability of the entire cohort being exposed to the same level of mortality risk from a flood event (Mann, 1976; Nunn et al., 2002, 2007a). In the 0+ fish, there is then the potential for a trade-off between their higher survival rates in the growth season versus their likelihood of achieving low body sizes at the end of that season (Bolland et al., 2007). This can potentially limit the over-winter survival rates of the smaller individuals that were produced relatively later in the summer (Nunn et al., 2007a, b).

In the rivers of eastern England, the fish communities are generally relatively diverse due to their previous connectivity with the Rhine-Danube systems following the last glacial period (Wheeler & Jordan, 1990). Conversely, rivers in western England are of lower diversity and so, in attempts to diversify angling opportunities, some non-indigenous species have been introduced, such as European barbel Barbus barbus (Linnaeus 1758) (Wheeler & Jordan, 1990; Antognazza et al., 2016). In river basins including the Wye and Severn, these introductions have been successful, with non-indigenous B. barbus establishing abundant populations (Amat-Trigo et al., 2017), despite the propensity of these rivers to flood regularly (Marriot, 1992) that could inhibit the ability of introduced fish to establish populations (Gutmann Roberts, 2018). Indeed, B. barbus is invasive in many other European rivers following introductions to enhance angling that have expanded their natural range (Kotlik & Berrebi, 2001; Meraner et al., 2013; Zaccara et al., 2019). Nevertheless, there is a considerable knowledge gap on how their spawning strategies in their non-indigenous range influence their 0+ fish cohorts and how these potentially influence their invasion success (Gutmann Roberts, 2018). This knowledge gap is important, as hypotheses on the establishment probability of introduced species include the ‘pre-adaptation hypothesis’, which suggests when introduced fishes share similar ecological traits and behaviours with native fishes, they ought to benefit through, for example, a similar ability to acquire resources (Duncan & Williams, 2002; Ricciardi & Mottiar, 2006; Buoro et al., 2016; Gutmann Roberts, 2018). Moreover, if this strategy as per that used by the introduced species in their native range, then there is minimal requirement for their reproductive traits to adapt to the new conditions (Schlaepfer et al., 2010; van Kleunen et al., 2011; Gutmann Roberts, 2018).

The aim of this study was to evaluate the range of spawning strategies used by both native and non-indigenous cypriniform fishes in a spate river in Western England, the River Teme. This river forms part of the lower River Severn basin and is prone to spates throughout the year (Gutmann Roberts, 2018). The approach used also allowed testing of the pre-adaptation hypothesis of invasion biology, where the non-indigenous B. barbus were used as the model invasive species. The objectives were to: (1) analyse temporal patterns in the lengths of 0+ cypriniform fishes and length at end of summer; (2) estimate separate sub-cohorts within a spawning season to determine spawning strategies for each (single vs. multiple spawning); and (3) determine if these spawning strategies, for each species, are consistent between different sites and across different years.

Materials and methods

Sampling sites

The sampling sites (n = 3) covered the majority of the B. barbus non-indigenous range in the study river (Fig. 1). With negligible off-channel habitat throughout the river, each sampling site consisted of areas of reduced flow rates within the river channel. Site 1 was the furthest upstream, located at Tenbury Wells (52°19′ N, − 2°24′ W) (Fig. 1), where the wetted width of the channel was 26 metres (07/10/15). The sampled area was located immediately downstream of a road bridge at the downstream end of a large gravel island, near to the right-hand bank. This site is in an urbanised area with little agriculture directly nearby, with a public footpath running along the right-hand bank. Riparian vegetation included overhanging trees (Salix spp.). Within the river, there was minimal instream vegetation, with the river generally running over gravel at depths of < 1 m. The sampling area comprised an area of minimal/ negligible flow, close to the right-hand bank (Gutmann Roberts, 2018). Site 2 was located at Knightwick (52°12′ N, − 2°23′ W) (Fig. 1), where the wetted width was 27 m (07/10/15). Sampling was from the right-hand bank, with samples taken either at the downstream end of an exposed gravel beach, or upstream of the gravel beach, in shallow water of a maximum depth of 1 m. Again, instream vegetation was minimal, with the sampling area comprising relatively slack water over a gravel substrate that was contiguous with a gravel riffle (Gutmann Roberts, 2018). Site 3 was the most downstream site at Powick (52°10′ N, − 2°14′ W) (Fig. 1), where wetted width was 15 m (07/10/15), with the sampling area located at the downstream end of a gravel riffle used by spawning B. barbus (Pinder et al., 2015). Sampling was conducted from the left-hand bank in an area of shallow, wadable water. The right-hand bank was steep, incised and suffering erosion with sheep grazing pasture extending up to the river (Gutmann Roberts, 2018). The annual temperature and flow regime of the River Teme for October 2015 to October 2016 are provided in Gutmann Roberts et al. (2019a); a notable feature of the river is large flood events, especially in winter (Brown et al., 1987), but can also occur throughout the year, including early summer (Gutmann Roberts, 2018).

Fig. 1
figure 1

Map of the three survey Sites (filled circle); 1Tenbury, 2 Knightwick and 3 Powick on the River Teme which flows into the River Severn. Urban areas are shaded grey, with inset showing position of sites and river within the UK

Sampling methodology

The sampling method for the 0+ fish was micro-mesh seine netting (net dimensions: 25 × 3 m; 2.5 mm mesh size), which captures more representative samples of larger, juvenile fishes than methods such as point abundance sampling by electric fishing (but with the latter able to provide higher sampling resolution for larval fishes) (Cowx et al., 2001; Nunn et al., 2002, 2003, 2007a, b; Copp, 2010). The net was deployed in appropriate larval and juvenile fish habitats at each site, as identified by areas of habitat off the main flow of the river, where there was sufficient depth (up to 1 m) and cover (including over-hanging trees and large stones) to provide refuge for the 0+ fishes (Gutmann Roberts, 2018). Seine nets were deployed by one person staying stationary at the riverbank holding an end of the net whilst the other person walked with the opposite end of the net in a downstream direction, ensuring the bottom of the net remained on the river bed, with the net was looped round by the person then moving out into the river and upstream before returning to the stationary sampler. The net was then brought into the riverbank whilst keeping the bottom of the net on the river bed (Gutmann Roberts, 2018). In 2015, sampling commenced in early July and concluded in October (Table 1). In 2016, sampling commenced earlier to determine the arrival of larvae in nursery habitats, hence sampling commenced in late May and concluded in October (Table 1). In 2017, sampling commenced in May and continued through to September. The rationale for concluding sampling in September/ October was a series of low catches of cypriniforms in the final samples as the 0+ fishes utilised alternative, non-accessible habitats, in combination with rising water levels that severely limited safe access to the sampling sites (Gutmann Roberts, 2018). Indeed, by October the few B. barbus and S. cephalus that were caught were all juveniles in the final stages of development (scales all over body and all fins fully formed) and thus were likely to be starting to utilise deeper water that was not accessible for sampling by micro-mesh seine netting (Copp, 1992; Bischoff & Freyhof, 1999).

Table 1 Micromesh seine net sampling dates at sites on the River Teme and day length based on location of S2 and mean daily river level from Bransford located between sites 2 and 3 from riverlevels.uk and mean daily water temperature from S2 Environment Agency data

Following netting, all fish were removed from the net and, where catches were sufficiently high at the species level, sub-samples were taken randomly from the main sample with a small hand net, euthanized (anaesthetic overdose, MS-222) and then preserved in 70% IMS. They were then kept in chilled conditions (approximately 5°C) until their processing in the laboratory (Gutmann Roberts, 2018).

Data collection and analysis

In the laboratory, for each sampling date and site, the fish were identified to species (Pinder, 2001) and measured using digital calipers (standard length (Ls), nearest 0.1 mm). Assessment of these data involved calculating their length distributions in 1 mm class intervals for B. barbus, S. cephalus, Phoxinus phoxinus (Linnaeus 1758) and L. leuciscus (Gutmann Roberts, 2018). These length distributions were plotted temporally by species and site to identify whether there was the appearance of ‘new’ fish into the cohort throughout the summer, i.e. whether fish of < 20 mm were regularly appearing in samples collected in July and August that would suggest multiple spawning activities (Nunn et al. 2002, 2007a).

These length distributions were then used to identify the presence of length modes in the samples per site and to assess their growth through each growth season using modal progression analysis (MPA) (Gutmann Roberts, 2018). This method was used here on the assumption that each mode identified by MPA represented a discrete spawning event and that each mode could be tracked through subsequent samples. For each species, site and sampling date, the length distributions were analysed for MPA by decomposition assessment using Bhattacharya’s method in FiSAT (Bhattacharya, 1967; Bolland et al., 2007; Hamidan & Britton, 2014). This analysis identifies the presence of modes in each length distribution by separating them into a series of normal distributions (King, 2013). For each mode, the output was the number of individuals, their mean length and standard deviation (SD) (Bolland et al., 2007), with modes separated by application of a separation index (SI), calculated as the ratio of the difference between successive means and the difference between the SD of their modes. Values of the SI > 2.0 indicate significant length differences with the other identified modes (Bhattacharya, 1967; Bolland et al., 2007). The overall output of MPA for each site per sample and species was thus the number of modes in the cohort and their mean length (± SD), plus their SI with adjacent modes. These outputs were then plotted for mean length (± 95% confidence interval) per mode and per sample for each species and site to identify and track the modes over time. This enabled visual assessment of the different modes in the cohort over time. Due to low sample sizes of B. barbus and S. cephalus being caught at Site 3 in 2016, these data were omitted from analyses. Site 1 also had a low sample size of B. barbus in 2016 and was not included for analysis of length modes (Gutmann Roberts, 2018).

Results

Overview of 0+ fish samples

The cypriniform fish samples from Sites 1 and 2 were dominated by P. phoxinus, with S. cephalus and B. barbus also prominent (Supplementary Material Tables S1 to S3). Other fish species occasionally present in samples but at numbers insufficient for further analyses included stone loach Barbatula barbatula (Linnaeus 1758), bullhead Cottus gobio (Linnaeus 1758), gudgeon Gobio gobio (Linnaeus 1758), and three-spined stickleback Gasterosteus aculeatus (Linnaeus 1758). The species composition of samples from Site 3 were similar to Sites 1 and 2, but with L. leuciscus also present (Tables S1 to S3).

Barbus barbus

In 2015, the lengths of the 0+ B. barbus in samples ranged between 13 and 37 mm across the three sites. Fish of the smallest SLs were recorded in July and fish present > 30 mm were present at all sites in August samples (Fig. S1). Similarly, in 2016, the length range of 0+ B. barbus was 12 to 37 mm, but with no B. barbus in samples until 29/06/16 (Fig. S2). In 2017, sampling commenced in May, with 0+ B. barbus first detected on 05/06/17 (Fig. 2).

Fig. 2
figure 2

Standard length (LS mm) distributions of 0+ Barbus barbus at Site 1, 2 and 3, River Teme 2017. Note differences in values on the Y axis for comparative purposes

In 2015, a relatively large size range of 0+ B. barbus was present in each sample, with the largest length range in the 20/08/15 sample (18 mm; Fig. S1). Length frequency distributions revealed fish of < 15 mm were still appearing in samples collected in August, with MPA consistently identifying three or four modes in the cohorts over the sampling period and in the sites where the SI exceeded 2.0 (Figs. S1, S3). In 2016, B. barbus of 11 to 12 mm were in samples collected on 28/07/16, and with fish of 9 to 12 mm also present in samples collected on 09/08/16 when other individuals were present to lengths > 30 mm (Fig. S2). Modal progression analysis on the 2016 samples for Site 2 consistently revealed three modes in the cohorts (Fig. S4). In 2017, up to five modes were detected at all three sites, with the majority able to tracked across the sampling period (Fig. 3).

Fig. 3
figure 3

Mean length (mm, ± 95% CI) per mode of Barbus barbus from Site 1–3, River Teme from 2017, as identified by Modal Progression Analysis. The symbols represent each mode

Squalius cephalus

In 2015, 0+ S. cephalus were present in all samples from all sites and at lengths of 13 to 37 mm (Fig. S5), with fish of < 20 mm were always present in samples until late August when fish > 30 mm were also present (Fig. S5). 0+ S. cephalus were in always samples from early June in 2016 and 2017 (Fig. 4; S6).

Fig. 4
figure 4

Standard length (LS mm) distributions of 0+ Squalius cephalus at Sites 1, 2 and 3, River Teme 2017. Note differences in values on the Y axis for comparative purposes

Similar to B. barbus, there was considerable variability in the S. cephalus size ranges per sample, providing evidence that there were extended spawning periods in each year (Figs. S5, S6 and 4). MPA was also able to consistently identify and track up to five modes through the cohorts (Figs. S7, S8 and 5); for example, fish of < 15 mm were present in samples collected from S1 in 2015 between early July and late August (Fig. S7).

Fig. 5
figure 5

Mean length (± 95% CI) per mode of Squalius cephalus from Sites 1, 2 and 3, River Teme from 2017, as identified by Modal Progression Analysis. The symbols represent each mode

Phoxinus phoxinus

For 0+ P. phoxinus, their lengths in samples varied between 7 and 50 mm (Fig. 6, S9, S10), although fish present in May samples of > 27 mm were considered to be age 1+ and were not included in subsequent analyses. In all years, smaller fish (< 20 mm) were apparent in samples collected throughout the summer, and even in October, although fish < 10 mm were generally only present in May to early August (Figs. S9, S10 and 6). As up to six modes could be tracked through the growth season, multiple spawning events were apparent in all sites and years (Fig.s 7, S11, S12).

Fig. 6
figure 6

Standard length (LS mm) distributions of Phoxinus phoxinus at Sites 1, 2 and 3, River Teme 2017. Note differences in values on the Y axis for comparative purposes

Fig. 7
figure 7

Mean length (± 95% CI) per mode of Phoxinus phoxinus from Sites 1, 2 and 3, River Teme from 2017, as identified by Modal Progression Analysis. The symbols represent each mode

Leuciscus leuciscus

The only site where L. leuciscus were sampled consistently was S3, where fish of < 20 mm were never present in samples collected after 1st July of each year, suggesting there was a discrete annual spawning period relatively early in spring, although the length range of the cohorts was relatively large by October 2015 (32 to 50 mm) (Figs. 8, 9, S13, S14). In 2016, 0+ L. leuciscus ranged from 18 to 43 mm, with 0+ fish in samples as early as 26/06/16 (Fig. S14). In 2017, sample sizes were reduced from previous years, but the fish maintained their pattern of a single spawning event (Figs. 8 and 9).

Fig. 8
figure 8

Standard length (LS) distributions of Leuciscus leuciscus at Site 3, River Teme 2017. Note variable y-axis

Fig. 9
figure 9

Mean length (± 95% CI) per mode of 0+ Leuciscus leuciscus from Site 3, River Teme from 2015 (white), 2016 (grey) and 2017 (black) with 0+ as squares, as identified by Modal Progression Analysis

Discussion

There was strong evidence that the majority of the cypriniform fishes in the River Teme had extended spawning periods that suggested a multiple spawning strategy. The native S. cephalus and P. phoxinus both had new cohorts of 0+ fish appear within the annual samples at each site, with multiple spawning also observed in the other populations in England (Nunn et al., 2002, 2007a). Conversely, native L. leuciscus appeared to have a single spawning event each year, which resulted in relatively larger 0+ fish appearing in samples, at least when compared to the other species. The presence of L. leuciscus solely at Site 3, the most downstream site, is likely due to adult dace residing in the River Severn and migrating into the River Teme tributary for reproduction (Kottelat & Freyhof 2007). The non-indigenous B. barbus also had a multiple spawning strategy, but where 0+ S. cephalus were present in samples from May, 0+ B. barbus were only present from June onwards.

The introduction of non-indigenous B. barbus into the River Severn and their subsequent dispersal into the River Teme thus provided the opportunity, via their spawning strategies, to test the pre-adaptation hypothesis of invasion biology (Ricciardi & Mottiar, 2006; Schlaepfer et al., 2010; van Kleunen et al., 2011). It was apparent from across the three sampling years, there was a consistent pattern of B. barbus having 0+ fish of < 20 mm regularly appearing in samples collected between June and late August, suggesting their spawning period extended over several weeks (e.g. between May and July). This multiple spawning was also detected in the native populations of S. cephalus and P. phoxinus. Whilst this multiple spawning strategy could have been mediated in B. barbus by plasticity in their reproductive traits, they utilise similar spawning strategies in their native range (Ovidio et al., 2009). For example, individual adult B. barbus in the River Ourthe, Belgium, were detected as spawning as least twice per year (Baras, 1995) and in captivity, B. barbus can spawn up to 15 times in 1 year under constant photoperiod and high thermal regimes (Poncin, 1992). Moreover, this apparent pre-adaption is not just limited to their spawning strategies, but also includes their somatic growth rates, as there were no significant differences in adult B. barbus growth rates between populations in their British indigenous and non-indigenous ranges (Britton et al., 2012). In entirety, these results suggest that the invasion success of B. barbus in the study river, and potentially elsewhere in their non-indigenous range, is at least partially related to their ability to express their life history traits and behaviours in a similar manner to their indigenous range.

The differences in the spawning strategies between L. leuciscus (single spawning, spring) and the other fishes (multiple spawning, early to mid-summer) at least partially relate to differences in the water temperatures and photoperiod required for initiating their spawning (Gutmann Roberts, 2018). For example, spawning in L. leuciscus usually starts when water temperatures exceed 10°C (Kennedy, 1969); in British rivers, this generally occurs between March and April (Mills, 1981) when no other cypriniform fishes are spawning, so they have minimal competition for spawning substrates and then nursery habitats. Consequently, this single spawning event could produce sufficient numbers of 0+ fish to ensure some survive any subsequent deleterious event, especially if they have achieved a relatively large body lengths by that time. Conversely, the spawning of B. barbus, S. cephalus and P. phoxinus tends to be initiated when water temperatures are higher, usually at least 11 to 12°C (Baras, 1995, Koç et al., 2007, Mills, 1987). As these three species all reproduce at similar times (i.e. May to July), and with anadromous fishes, such as shad Alosa spp. and sea lamprey Petromyzon marinus (Linnaeus 1758) also spawning in the river then and on the same spawning locations (Pinder et al., 2015), then there is potential for relatively high competition for spawning substrates and nursery habitats (Gutmann Roberts, 2018; Gutmann Roberts et al., 2019b). Consequently, in addition to increasing the resilience of 0+ fish to deleterious stochastic events that would subsequently impact recruitment success (Nunn et al., 2002), multiple spawning events could also act to reduce intra- and inter-specific competition for spawning and nursery habitats after June (Gutmann Roberts, 2018).

The spawning behaviour of the adult fish was not assessed directly in the study and so the inferences being made from the 0+ fish samples could not be verified. Nevertheless, the length distributions of the 0+ B. barbus were generally consistent with the production of progeny over prolonged periods (Nunn et al., 2002, 2007a, b), with multiple spawning also strongly suggested in the length data of S. cephalus. This reproductive behaviour has also been detected in other S. cephalus populations, including the River Spree, Germany (Fredrich et al., 2003), and the River Rokytna, Czech Republic (Lelek, 1980), where individual adults were observed spawning twice in one spawning season. In England, multiple spawning events have been apparent in S. cephalus and P. phoxinus from length data from 0+ fish (Nunn et al., 2007a). Although larger 0+ B. barbus (> 40 mm) were not present in the River Teme samples, this was assumed to be the result of these individuals attaining body sizes that enabled their utilisation of mid-channel habitats of higher flows, as in their native range, individuals can only withstand water velocities in excess of 10 cm s−1 when they are of body lengths over 50 mm (Bischoff and Freyhof, 1999). By late summer, minnow and dace prefer higher water velocities, although minnow have strong preferences for shallow water (Copp 2010), whereas chub prefer deeper waters (Copp 1992). Although there are both similarities and differences in these larval and juvenile fishes in their microhabitat use, the dietary niche of barbel can overlap with these species, although the extent to which this occurs varies spatially (Gutmann Roberts & Britton 2018). Competition from non-indigenous barbel will depend on densities of fish in comparison with native confamilials and prey availability (Gutmann Roberts, 2018). The densities of larval and juvenile cypriniform fish in similar habitats is also dependent, at least initially, on the fecundity of adults, although the relationship between adult fecundity and spawner density versus the subsequent density of larvae remains uncertain (Gutmann Roberts, 2018; Gutmann Roberts et al., 2019b).

The consequence of the multiple spawning of B. barbus and S. cephalus in the River Teme was that in samples collected in late summer and autumn when the 0+ fish growth year was almost complete, some fish were still present in samples at lengths below 20 mm. This indicates that the fish produced relatively late in spawning cycle were unable to compensate for their short growth season through faster growth rates. For S. cephalus in other English rivers, this had a life-long effect on the growth of individuals, where those fish of smaller lengths at the end of their first growth year then remaining relatively small for their age throughout life, but with a tendency to be longer lived than larger, faster growing fish (Bolland et al., 2007). These small 0+ fish body sizes are also potentially important for recruitment, as length at the end of the first growth season is a significant correlate of recruitment strength in some riverine cypriniforms, including S. cephalus (Nunn et al., 2007a, b). In the River Teme, winter spates can result in river levels increasing by over 4 m in several hours, with concomitant increases in flows (e.g. to over 60 m3 s−1) (Gutmann Roberts, 2018). In conjunction with negligible off-channel refugia and the ability of fish to hold position in flows as a function of their body size (Müller et al., 1996), this suggests that over-winter survival of 0+ B. barbus at lengths below 20 mm might be limited. This, however, remains speculative in the absence of knowledge of how B. barbus late spawning events translates into lifetime consequences for individuals and cohorts (Gutmann Roberts, 2018). Nevertheless, Nunn et al. (2010) revealed that 0+ S. cephalus in the River Trent were able to survive the winter period at lengths as small as 13 mm.

In conclusion, multiple spawning events were detected here in this non-indigenous B. barbus population, with this spawning strategy also utilised by two native fishes, as well as B. barbus in their native range. Consequently, it was argued that this aspect of the population’s invasion biology had consistency with the pre-adaption hypothesis. The potential implications of these spawning behaviours are increasing the resilience of 0+ cohorts to deleterious and stochastic events, such as floods, as well as potentially, reducing inter-specific competition for spawning and nursery resources. However, for fish produced later in the summer there is a consequence of these individuals attaining relatively small body sizes at the end of their first growth season that could impact their ability to over-winter.