Research Article |
Corresponding author: Kenneth J. Killgore ( jack.killgore@usace.army.mil ) Academic editor: Michal Janáč
© 2024 Kenneth J. Killgore, Jan J. Hoover, William T. Slack, James P. Kirk, Bradley R. Lewis, Steven G. George, Leandro E. Miranda.
This is an open access article distributed under the terms of the CC0 Public Domain Dedication.
Citation:
Killgore KJ, Hoover JJ, Slack WT, Kirk JP, Lewis BR, George SG, Miranda LE (2024) Population characteristics of silver carp from the source of their North American introduction in the Lower Mississippi River. Aquatic Invasions 19(3): 329-343. https://doi.org/10.3391/ai.2024.19.3.133295
|
Silver carp, Hypophthalmichthys molitrix, escaped into the Lower Mississippi River (LMR) over 50 years ago, established reproductive populations, and spread across much of the Mississippi River Basin. Demographic rates of silver carp are needed to inform decisions on control and management of this invasive species, but have not been published for the LMR. The purpose of this paper is to report silver carp age and growth estimates from fish collected in riverine (mainstem) and backwater (lake) habitats in the LMR during the period 2011–2019, to compare our results with populations from other geographic areas in the Upper Mississippi River drainage, and to evaluate latitudinal and habitat differences in demographic parameters. Silver carp gained weight with increasing length similarly throughout the lower and upper basin. However, annual growth rates were higher in the LMR compared to northern rivers including the Illinois, Wabash, Missouri, and Middle Mississippi rivers. In the LMR, regression analyses demonstrated that females were heavier in lakes than males or females in the mainstem and that females in lakes had the lowest instantaneous mortality (-0.186). Maximum age was 8 and 10 years for females and males, respectively. The largest male weighed 13.8 kg with a total length of 1022 mm, and was 7 years old. The largest female weighed 16.0 kg with a total length of 1034 mm TL, and was 7 years old. Rapid growth rates, larger sizes, and lower mortality in the LMR, in combination with limited commercial fishing, extensive river-floodplain connectivity, and vast amounts of spawning areas, ensure that LMR silver carp will continue to act as a source of fast-growing invasive individuals for other reaches and other rivers throughout the Mississippi River Basin.
Latitudinal gradients, recruitment, mortality, age, growth, mainstem channel, oxbow lakes, backwater
Silver carp, Hypophthalmichthys molitrix, are one of the most widely cultured fish in the world and are now found in 88 countries and territories including the United States (
Silver carp can weigh over 27 kg, grow to over 1 m total length, and live up to 20 years in their native range of Southeast Asia (
Growth rates of silver carp in the Lower Mississippi River mainstem (LMR), which extends from the mouth of the Ohio River to the Gulf of Mexico, have not been previously reported despite the fact that they originated in Arkansas and the LMR represents the center of population abundance and a primary dispersal route. The LMR provides unlimited spawning sites in flowing reaches and extensive backwaters and lakes with high productivity for rearing (
The LMR mainstem sites included six locations extending from the mouth of the Ohio River (river kilometer, Rkm 1534) downstream to near the mouth of the Mississippi River at the Gulf of Mexico (Rkm 32) (Fig.
Capture locations of silver carp in the Lower Mississippi River (LMR) from 2011-2019. Boxed river kilometers (Rkm) indicate mainstem sampling reaches. Lakes indicated by boxed numbers 1-4. Samples were collected from Fort Jackson (Rkm 32, 29.2130N, -89.2722W) upstream to the near the mouth of the Ohio River (37.0053N, -89.1765W).
To characterize the large-scale variation in silver carp length at age and their growth performance throughout the vast LMR, we relied on samples collected opportunistically as part of various investigations at the U.S. Army Engineer Research and Development Center in Vicksburg, MS. Four gear types were used to collect silver carp. Gill nets produced 65% of the total catch, electrofishing 23%, hoop nets 3%, and the remaining 9% were jumpers that landed in one of our survey boats. Most gill nets were multifilament with 10.2 cm bar mesh, 91 m long, and 2.4 m deep. A 46 by 2.4 m trammel net with 6.3 cm square mesh was also used. Silver carp from Horn Lake were by-catch of commercial fishing during the 2011 paddlefish Polyodon spathula harvest which used 12.7 cm bar mesh. Boat electrofishing was generally operated as pulsed direct current, 60 Hz, and approximately 5 to 7 mean amps read from a Smith-Root GPP unit®. Hoop nets were 1.2 m diameter opening, 7 hoops, 4.6 m long, and 7.6 cm square mesh.
All silver carp collected were euthanized, and immediately after, fish were measured for total length to the nearest 1 mm and for weight to the nearest 10 g. Gonads were examined to determine sex. The leading pectoral fin rays were removed, labeled, and stored for subsequent aging. To prepare for aging, fin rays were sectioned along the basal portion of the ray to a thickness of approximately 0.70 mm using a Buehler Isomet® saw with a diamond wafering blade. Aging involved counting annuli on the sectioned rays while viewed under transmitted light with an Olympus SZX16® zoom stereomicroscope, equipped with an Olympus DP72® camera system and cellSens® imaging software. Two readers agreed in advance about what constituted an annular mark and independently aged each fish by counting annual marks. Opaque zones viewed under reflected light were counted as annuli. When a discrepancy occurred (about 45% of the time, 36% of those with a difference of 1 year), a consensus age was assigned if agreement could be reached; otherwise, the rays were discarded. The percent agreement and a coefficient of variation (CV) were calculated to document precision differences among agers (
We examined weight-length and mortality curves using linear methods, and growth curves using nonlinear methods. The relation between total body weight (W; g) and total length (L; mm) is represented by a logarithmic curve typically applied with a log10 transformation of W and L:
log 10 W = b0 + b1 ∙ log10L (1)
where b0 is an intercept parameter and b1 is a slope parameter that describes weight increase relative to length increase. To test for potential differences in weight-length relations between habitats and between genders, equation 1 was expanded into an analysis of covariance to include habitat type (H: main channel = 0, lakes = 1), and fish gender (G: male = 0, female = 1):
log 10 W = b0 + b1 ∙ log10L+ b2 ∙ log10L ∙ H + b3∙ log10L ∙ G + b4 ∙ log10L ∙ H ∙ G (2)
where b0 is an intercept parameter, b1 tests for an effect of length on weight, b2 tests for an effect of habitat on the length-weight model, b3 tests for an effect of gender on the length-weight model, and b4 tests if an effect of habitat on the length-weight model depended on gender, or if an effect of gender depended on habitat.
The relation between number of fish (N) and age (years) is represented by a logarithmic curve that is typically linearized with a loge transformation of N (i.e., catch curve):
logeN = b0 + b1 ∙ age (3)
where b0 is an intercept parameter and b1 is a slope parameter that represents the instantaneous rate of annual mortality. To test for potential differences in mortality between habitats and between gender equation 3 was expanded to include habitat type and fish gender:
logeN = b0 + b1 ∙ age + b2 ∙ age ∙ H + b3 ∙ age ∙ G + b4 ∙ age ∙ H ∙ G (4)
where b0 is an intercept parameter, b1 tests for an effect of age on number, b2 tests if b1 is influenced by habitat, b3 tests if b1 is influenced by gender, and b4 tests if any effect of gender on b1 depends on habitat, or any effect of habitat depends on gender. Equations 1–4 were fit with the GLM procedure (
Length at age data were used to estimate von Bertalanffy growth coefficients (
Lt = L∞ – (L∞ – L0) ∙ e–Kt (5)
where Lt is the length at age t in years, L∞ is the mean asymptotic length, L0 is the length at time zero (i.e., hatching), K is a growth coefficient, and e is the exponent for natural logarithms. L0 was fixed as 5.1 mm based on estimates published by
Lt = [L∞ + (b1 ∙ H) + (b2 ∙ G) + (b3 ∙ H ∙ G)] – [L∞ + (b1 ∙ H) + (b2 ∙ G) + (b3 ∙ H ∙ G) – L0] ∙ e[–K + (b4 ∙ H) + (b5 ∙ G) + (b6 ∙ H ∙ G)] ∙ t (6)
where b1 tests if L∞ depends on habitat, b2 tests if L∞ depends on gender, and b3 tests if L∞ depends on habitat and gender. In a similar manner, the tests conducted by b4, b5, and b6 are analogous to those performed by b1 - b3, with the distinction that they are applied to the variable K rather than L∞ Equations 5–6 were fit with the NLIN (nonlinear) procedure (
A total of 360 silver carp comprised of 209 females and 151 males were collected in the LMR mainstem and lakes from 2011 to 2019. The weight-length regression for all fish combined (equation 1; Table
Relations between log10 weight (g) and log10 total length (mm) for male and female silver carp collected in mainstem and lake habitats in the Lower Mississippi River, 2011-2019. Equation 1 was a regression model fit to the full data set regardless of gender or habitat; equation 2 was an analysis of covariance that included the effect of gender, habitat, and their interaction on the weight-length relation.
Parameter | Estimate | SE | t-value | P > t |
---|---|---|---|---|
Equation 1 (R2 = 0.95) | ||||
b0 | -9.88 | 0.120 | -82.1 | <0.01 |
b1 | 3.665 | 0.041 | 88.6 | <0.01 |
Equation 2 (R2 = 0.96) | ||||
b0 | -9.596 | 0.159 | -69.3 | <0.01 |
b1 | 3.572 | 0.047 | 75.8 | <0.01 |
b2: mainstem | -0.012 | 0.003 | -4.3 | <0.01 |
b2: lakes | 0 | |||
b3: male | -0.010 | 0.002 | -4.6 | <0.01 |
b3: female | 0 | |||
b4: mainstem, male | 0.013 | 0.004 | 3.7 | <0.01 |
b4: mainstem, female | 0 | |||
b4: lakes, male | 0 | |||
b4: lakes, female | 0 |
Regional comparisons of weight-length regressions indicated that the slope of LMR fish weight gains were greater than upper basin fish (Table
Length-weight regression for silver carp based on
River | L-W Regression Equation | R2 | Intercept 95% CL | Slope 95% CL | 450 mm | 800 mm | Reference | ||
---|---|---|---|---|---|---|---|---|---|
South Dakota Tributaries, n=469 | |||||||||
James | log10 weight = -5.26+3.11(log10 length) | 0.96 | -5.43 | -5.1 | 3.05 | 3.17 | 981 | 5869 |
|
Vermillion | log10 weight = -4.82+2.90(log10 length) | 0.98 | -5.16 | -4.47 | 2.82 | 3.07 | 748 | 3971 |
|
Big Sioux | log10 weight = -5.53+3.21(log10 length) | 0.98 | -5.91 | -5.14 | 3.07 | 3.36 | 970 | 6150 |
|
Missouri River | |||||||||
Gavins Point, n=7 | log10 weight= -6.92+3.70(log10 length) | 0.97 | 0.21 | 0.59 | 788 | 6628 |
|
||
Interior Highlands, n=68 | log10 weight = -5.35+3.13(log10 length) | 0.93 | 0.76 | 0.21 | 900 | 5453 |
|
||
Middle Mississippi River/Illinois River | |||||||||
Middle, n=145 | log10 weight = -5.29+3.11(log10 length) | 0.81 | 915 | 5477 |
|
||||
Illinois River, n=452 | log10 weight = -5.29+3.12(log10 length) | 0.99 | 972 | 5856 |
|
||||
Lower Mississippi River, n=360 | log10weight= -6.88+3.66(log10length) | 0.95 | -6.647 | -7.121 | 3.58 | 3.75 | 691 | 5695 | This Study |
Ages were estimated for 360 silver carp ranging from 2 to 10 years. Of this total, 59% of the pectoral fin rays were read by two readers independently. Percent agreement was 77% for ages 1–2, decreased to 51% for ages 3–5, and 53% for older (ages 6–10) fish. Overall, there was a 56% agreement on the same age, a 95% agreement for ±1 year, and a coefficient of variation of 122% (n=211).
A total of 151 males and 209 females were aged with an overall sex ratio between males and females of 1:1.3. However, sex ratio differed between habitats with females constituting 66% of the individuals in lakes compared to 42% females in mainstem. Males and females had similar age distributions (Table
Age (years) frequency of silver carp in the Lower Mississippi River, 2011-2019.
Age | Females (N = 209) | Males (N = 151) | ||||||
---|---|---|---|---|---|---|---|---|
Frequency | Percent | Cumulative Frequency | Cumulative Percent | Frequency | Percent | Cumulative Frequency | Cumulative Percent | |
2 | 10 | 4.8 | 10 | 4.8 | 17 | 11.3 | 17 | 11.3 |
3 | 27 | 12.9 | 37 | 17.7 | 19 | 12.6 | 36 | 23.8 |
4 | 37 | 17.7 | 74 | 35.4 | 35 | 23.2 | 71 | 47.0 |
5 | 45 | 21.5 | 119 | 56.9 | 31 | 20.5 | 102 | 67.6 |
6 | 30 | 14.4 | 149 | 71.3 | 24 | 15.9 | 126 | 83.4 |
7 | 49 | 23.4 | 198 | 94.7 | 20 | 13.3 | 146 | 96.7 |
8 | 11 | 5.3 | 209 | 100.0 | 3 | 2.0 | 149 | 98.7 |
9 | 0 | 0 | ||||||
10 | 0 | 2 | 1.3 | 151 | 100.0 |
The catch curve for all fish combined (equation 3; Table
Relation between loge number of fish and age (years) for male and female silver carp collected in mainstem and lake habitats in the Lower Mississippi River, 2011-2019. Equation 3 was a regression model fit to the full data set regardless of gender or habitat; equation 4 was an analysis of covariance that included the effect of gender, habitat, and their interaction on the relation between number and age.
Parameter | Estimate | SE | t-value | P > t |
---|---|---|---|---|
Equation 3 (R2 = 0.45) | ||||
b0 | 4.317 | 0.37 | -7.9 | <0.01 |
b1 | -0.372 | 0.088 | -4.2 | <0.01 |
Equation 4 (R2 = 0.80) | ||||
b0 | 4.411 | 0.159 | 11.9 | <0.01 |
b1 | -0.454 | 0.064 | -7.1 | <0.01 |
b2: mainstem | 0 | |||
b2: lakes | 0.067 | 0.048 | 1.4 | 0.18 |
b3: male | 0 | |||
b3: female | -0.055 | 0.061 | -0.9 | 0.38 |
b4: mainstem, male | 0 | |||
b4: mainstem, female | 0 | |||
b4: lakes, male | 0 | |||
b4: lakes, female | 0.256 | 0.080 | 3.2 | <0.01 |
The growth model estimated that silver carp in the LMR attained an L∞ of 911 mm TL, with 95% confidence limits stretching 891–930 (equation 5, Table
Von Bertalanffy growth coefficients of silver carp collected in the Lower Mississippi River from 2011-2019 using a nonlinear model where L∞ is the mean asymptotic length and K is a growth coefficient (Equation 5). Equation 5 was expanded to test for potential differences in growth characteristics between habitats and between fish gender (Equation 6).
Parameter | Estimate | Approximate SE | Lower 95% CL | Upper 95% CL |
---|---|---|---|---|
Equation 5 | ||||
L∞ | 911 | 9.789 | 891 | 930 |
K | 0.548 | 0.025 | 0.498 | 0.598 |
Equation 6 | ||||
L∞ | 776 | 15.573 | 745 | 806 |
K | 0.925 | 0.102 | 0.726 | 1.125 |
b1: lake | 137 | 29 | 80 | 194 |
b1: mainstem | 0 | |||
b2: female | 48 | 27 | -5 | 102 |
b2: male | 0 | |||
b3: lake, female | 29 | 41 | -52 | 111 |
b3: lake, male | 0 | |||
b3: mainstem, female | 0 | |||
b3: mainstem, male | 0 | |||
b4: lake | -0.425 | 0.118 | -0.657 | -0.194 |
b4: mainstem | 0 | |||
b5: female | -0.119 | 0.145 | -0.404 | 0.165 |
b5: male | 0 | |||
b6: lake, female | 0.037 | 0.159 | -0.276 | 0.350 |
b6: lake, male | 0 | |||
b6: mainstem, female | 0 | |||
b6: mainstem, male | 0 |
The model predicted that silver carp in the LMR reached total lengths of 600 mm or greater by age 2, and almost 900 mm total length by age 5 (Fig.
Von Bertalanffy growth curves for silver carp (combined sexes) in LMR mainstem and lakes compared to the Tennessee River (
Since silver carp were introduced into the USA over 50 years ago, they have become widespread in the LMR potentially providing source populations for range expansion and establishment.
Silver carp gained weight with increasing length at larger sizes similarly throughout the lower and upper basin, whereas annual growth rates were faster in the lower basin including the Tennessee River (
Spatial variation in growth, weight, lifespan and mortality of silver carp have been noted in several studies. In addition to this study, a latitudinal decrease in silver carp growth rates have been reported (
Lock and dams on the Upper Mississippi, Illinois, and Ohio rivers affect passage (
Our study contributes to a better understanding of key demographic rates of silver carp in the LMR needed for three distinct goals: i.) evaluating the probability of local establishment based on an age-structured population model, predicting the risk of establishing a population, and identifying the most vulnerable stage to control (
KJK, JJH, WTS, JPK, LEM contributed to research conceptualization, sampling, and data analysis, and BRL and SGG led field and laboratory data collection activities
All research was performed in accordance with the Guidelines for the Use of Fishes in Research by the American Fisheries Society, and under the Institute of Animal Care Committee at the Engineer Research and Development Center, and accredited by the Care and Use of Laboratory Animals in DOD Programs.
Funding was provided from multiple sources within U. S. Army Corps of Engineers including the Aquatic Nuisance Species Research Program, the U. S. Army Corps of Engineers New Orleans District, and the U. S. Army Corps of Engineers Mississippi Valley Division.
We thank K. Boysen, J. Collins, A. Katzenmeyer, W. Lancaster, B. Maus, and C. Murphy with the U.S. Army Corps of Engineers team for their assistance in collecting silver carp. Amanda Oliver created the map. Funding was provided from multiple sources within U. S. Army Corps of Engineers including the Aquatic Nuisance Species Research Program, the U. S. Army Corps of Engineers New Orleans District, and the Mississippi River Geomorphology and Potamology Program at the U. S. Army Corps of Engineers Mississippi Valley Division. We thank three anonymous reviewers for excellent comments and suggestions. This study was partially performed under the auspices of the Institute and Animal Care Committee protocol # EL-FR-2022-1. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.