Research Article |
|
Corresponding author: Michaela Palmieri ( michaelapalmieri10@gmail.com ) Academic editor: Charles Martin
© 2026 Michaela Palmieri, Leandro E. Miranda, Melanie R. Boudreau, Corey G. Dunn, Leslie M. Burger, Dennis Riecke.
This is an open access article distributed under the terms of the CC0 Public Domain Dedication.
Citation:
Palmieri M, Miranda LE, Boudreau MR, Dunn CG, Burger LM, Riecke D (2026) Forecasting spread of invasive fish over a largescale network of lakes using local expert knowledge. Aquatic Invasions 21(2): 127-146. https://doi.org/10.3391/ai.2026.21.2.190069
|
Understanding spatial distribution patterns is essential to management of invasive species. Aquatic invasive species can be notably challenging to detect due to the substantial effort required to locate them underwater. This limitation has resulted in a lack of timely distribution maps, particularly over vast regions, and hindered efforts to understand, forecast, and manage the proliferation of invasive bigheaded carps (Hypophthalmichthys spp.). Much of the Mississippi River basin, particularly the Lower Mississippi Alluvial Valley, provides access to a massive network of interconnected floodplain lakes. In the absence of lake-specific monitoring data on carp occurrence status, we used local expert knowledge, provided by fish managers interviewed virtually, in conjunction with Maximum Entropy (MaxEnt) modeling, to predict bigheaded carps distribution in relation to lake physical characteristics. We predicted widespread carp invasion in more than 60% of over one thousand floodplain lakes, with lake size, inundation, and proximity to rivers closely related to carp presence. The resultant distribution map may be imprecise given the swift proliferation of bigheaded carps and sparse monitoring data, but it offers a baseline upon which presence data and range can be compared. This assessment method is also a resource for identifying priority management and conservation areas and can serve as a first step in conservation planning.
Bighead carp, floodplain lakes, MaxEnt, silver carp, species distribution model
Invasive organisms can be a major threat to ecosystems by disrupting biodiversity and ecosystem functions (
Understanding the geographical patterns of species distributions is key for the successful management of invasive species (
Bighead carp, Hypophthalmichthys nobilis (Richardson, 1845), and silver carp, H. molitrix (Valenciennes, 1884), collectively known as bigheaded carps, are invasive fish species in North America. They originated from East Asia and were moved to the Mississippi River basin in the 1970s for biocontrol in aquaculture operations. Likely due to flooding events, by the 1990s they had escaped captivity and began dispersing and reproducing across the Mississippi River basin (
Much of the Mississippi River basin, particularly the Lower Mississippi Alluvial Valley (LMAV), provides access to massive networks of interconnected floodplain lakes. Floodplain lakes are identified as the primary natural habitat for the development of juvenile bigheaded carps and for rapid growth to adulthood (
A key challenge in managing rapidly expanding invasive bigheaded carps in this extensive aquatic network is having information on their spatial distribution. Mapping the known and predicted locations of bigheaded carps is central to guiding effective management and assessing the effect of management actions (
Conventional lake monitoring used for fish populations can be effective at relatively small scales but becomes expensive at scales as large as the hundreds of lakes scattered throughout the LMAV. One low-cost option for gathering information on carp distribution could be by collecting local expert knowledge, which is accumulated through personal observations by people working in a region. This type of information has been used for documenting species occurrences in previous studies (
The LMAV is a large area in the United States that stretches from Cape Girardeau, Missouri until its terminus south of New Orleans, Louisiana. The region is roughly 1,000 km long and up to 150 km wide, spanning an area approximately 100,000 km2 (
The LMAV was essentially covered with bottomland hardwood forests before the European arrival (
To gather bigheaded carp presence data, we relied on the local knowledge available from agency fish biologists. The LMAV includes seven states, and each of these states is organized into multiple contiguous and non-overlapping geographic management areas. Depending on the state, these areas have various names (e.g., management districts, service areas, management regions), but we hereafter refer to these areas as fish management districts. Thus, the LMAV was geographically stratified into the existing fish management districts, and each of the 1,350 permanent lakes in the LMAV identified by
Fish biologists from each fish management district were identified and contacted to schedule interviews regarding bigheaded carp occurrence. Interviews were conducted virtually from January to July of 2023. Interviews began by reciting a script that provided a brief overview of the project and interview procedure (Appendix
We used Maximum Entropy (MaxEnt) presence-only modeling to predict the distribution of bigheaded carps because absence data were deemed unreliable. In place of absences, MaxEnt uses background points sampled from the study area for which presence is unknown, i.e., pseudoabsences (
Hypothesized lake characteristics that could be important in determining bigheaded carp (Hypophthalmichthys spp.) distribution in lakes of the Lower Mississippi Alluvial Valley (LMAV). CV = coefficient of variation; * = retained by Lasso procedure. Sources:
| Variable | Definition | Mean | CV | Rationale |
|---|---|---|---|---|
| Area* | Surface area of lake polygon feature (ha). | 76 | 362 | Surface area are correlated to fish assemblages in the LMAV ( |
| Mean width | Mean of several width measurements taken approximately every 100 m along the length of the lake (m). | 172 | 117 | Lake size and shape can influence fish species distributions ( |
| Length* | The length of a line drawn through the center of the lake along its longest axis (km). | 3.0 | 138 | Elongated lakes tend to have strong interactions with the terrestrial ecosystem, allowing more nutrient input and possibly access corridors ( |
| Length-width ratio | Length divided by mean width. | 27 | 158 | Lake shape can influence fish species distributions ( |
| Perimeter | Length of shoreline surrounding a lake (km). | 7.7 | 165 | Lake size and shape can influence fish species distributions ( |
| Inundation index* | Relative frequency of lake inundation from 1983 to 2011, with 100 indicating the entire lake is always fully covered by water and values closer to 0 indicating that the area covered by water is reduced to nearly nil. | 42 | 45 | Flooding is strongly linked to presence of young-of-year bigheaded carp. Inundation can also increase connectivity and provide corridors for carp movement ( |
| Shoreline development index | Ratio of lake shore length to the circumference of a circle with lake area. | 3.0 | 54 | Lake size and shape can influence fish species distributions ( |
| Agriculture cover* | The percent of area within a 1-km buffer around the lake polygon that is comprised of agricultural cover (%). | 43 | 70 | Eutrophication due to excessive agricultural runoff can lead to anoxia and fish death ( |
| Forest and wetland cover | The percent of area within a 1-km buffer around the lake polygon that is comprised of forest or wetland cover (%). | 41 | 70 | Wetlands may provide a corridor between rivers and lakes, and reduce excess nutrient concentrations of through flowing water, influencing water quality and possibly carp distribution ( |
| Batture* | Whether or not the lake is protected from major floods by a levee or floodwall, designated by a 0 if protected and a 1 if not protected. | 0.53 | 94 | Lakes regularly exposed to floods, like those between levees and the Mississippi River (batture) are more likely to be connected to waterways ( |
| Disconnectedness* | Maximum number of consecutive months that a lake remained disconnected from a river during 2000–2022. | 156 | 59 | Carp cannot complete their lifecycle in lakes and must travel from lake to river to spawn. Greater disconnection may make it more difficult and less likely for invasion to occur in lakes that are disconnected ( |
| Longitude | The east-west position of a lake. | -91.04 | 0.01 | Geographic location can influence water quality ( |
| Latitude | The north-south position of a lake. | 33.79 | 4 | Geographic location has been shown to influence water quality ( |
| Nearest neighbor* | Straight line distance to the nearest lake (km). | 3.0 | 97 | Lakes located close to other lakes that have been invaded may be more likely to be invaded ( |
| Distance to river* | Shortest network distance between a lake and a river with mean annual flow of at least 170 m2 s-1 (km). | 26 | 129 | Bigheaded carps have been documented in rivers throughout the Lower Mississippi Alluvial Valley ( |
| Probability depth* | Probability of the lake maximum depth being > 2 m. | 0.63 | 29 | Shallow depth and backwaters have been associated with higher bigheaded carp density ( |
Highly correlated environmental variables can affect model performance (
The MaxEnt model was applied to 1,086 lakes in the LMAV that had complete environmental variable data available. We used presence records obtained from the interviews, and we randomly created 10,000 pseudoabsences from lakes that had no known presence records (
To assess model performance, we used the ROCR R package (
To identify lakes potentially suitable for bigheaded carps (objective 1), we transformed the continuous MaxEnt output into a binary map of predicted suitable/unsuitable lakes using a maximum model sensitivity plus specificity threshold (
A total of 26 interviews were completed with agency fish biologists from 13 fish management districts throughout the LMAV. We administered the questionnaire for 131 of the 1,086 lakes and recorded 64 carp presences (48.9% of sampled lakes) from lakes within the study area (Fig.
The LASSO model implemented with logistic regression identified nine environmental variables suitable for inclusion as covariates in the MaxEnt model. These variables included lake area, lake length, probability depth, inundation index, agriculture cover, batture, disconnectedness, distance to river, and nearest neighbor (Table
Percent contribution and permutation importance of covariates when the final MaxEnt model was run to predict bigheaded carp (Hypophthalmichthys spp.) presence throughout the Lower Mississippi Alluvial Valley floodplain (LMAV) lakes. Variables are as defined in Table
| Variables | Percent contribution | Permutation importance |
|---|---|---|
| Area | 36.0 | 54.80 |
| Inundation index | 19.3 | 16.60 |
| Distance to river | 18.4 | 3.50 |
| Nearest neighbor | 9.7 | 6.70 |
| Agriculture cover | 7.0 | 5.30 |
| Length | 5.6 | 8.60 |
| Probability depth | 2.6 | 1.00 |
| Disconnectedness | 1.3 | 3.60 |
| Batture | 0.1 | 0.00 |
Adding the maximum model sensitivity to the model specificity threshold resulted in a threshold of 0.24, which was used to transform continuous habitat suitability values to binary predicted presences (>0.24) and absences (<0.24). This transformation forecasted that roughly 62% of the lakes in the LMAV supported bigheaded carps (Fig.
We observed a geographic pattern in predicted carp habitat suitability across the LMAV. The Moran’s I index value was 0.12, revealing a positive and significant (P < 0.001) spatial autocorrelation, indicating that lakes that may support carp tend to cluster. In addition, nearest neighbor was moderately important in predicting suitability in the MaxEnt model. Clustering predominantly occurred along major rivers (Fig.
Bigheaded carps are widespread and continue to expand their distribution in the Mississippi River basin, creating significant risks to the sustainability of native freshwater ecosystems (
We acknowledge that while the map may be imprecise given the swift proliferation of bigheaded carps and sparse monitoring data, the map is still an informative management tool. Its value lies in revealing broad patterns, core areas, prioritizing surveillance and monitoring, guiding interventions, and supporting communication and planning. The map shows that this expansion can occur throughout this extensive floodplain given the lake characteristics available to support occurrence of carps in this region. Considering that these lakes offer bigheaded carp preferred habitats and that carps have access to a diversity of river systems to facilitate reproduction and dispersal, we anticipate that dispersal throughout the Mississippi basin is likely, barring obstacles (e.g., dams) or infrequent connections (
In the absence of lake-specific monitoring data on carp occurrence status, we used local expert knowledge in conjunction with MaxEnt to predict bigheaded carp potential occurrences based on lake physical characteristics. We predicted widespread carp invasion throughout LMAV lakes with lake size, inundation potential, and proximity to rivers closely related to carp distribution in the LMAV. Most lakes likely to have carps were < 25 ha in size and were <5 km from a river. We suspect that these small lakes may have a higher likelihood of carp occurrence because their river proximity exacerbates the chance of recurring flooding that allows for carp dispersal and occupancy. A potential bias is that carp may be more visible in smaller lakes because these waters are easier to scan and have a higher proportion of littoral area, which is their preferred environment (
The relationships between inundation index, distance to rivers, and predicted habitat suitability are intuitive. Inundation index is a measure of long-term fluctuation in lake volume. Predicted habitat suitability was generally positive, where lakes with higher inundation index were predicted to be more likely to be suitable for bigheaded carp. A high inundation index means it is more likely that a lake remains close to full volume over time. Bigheaded carps are tolerant of poor water-quality (
We determined a significant geographic pattern existed in predicted carp occurrence across the LMAV, with nearly two thirds of lakes predicted to be suitable for carp often clustering in the northwest LMAV along the White River. Similarly, nearest neighbor was moderately important in predicting lake suitability. Because bigheaded carp were first released near this region, they have had a longer period to invade. Moreover, fish scientists in this region may have developed more programs and competence for surveillance of bigheaded carp populations. Spatial clustering of predicted carp occurrences within the river system of initial introduction could foreshadow that bigheaded carps will eventually invade higher densities of lakes throughout the LMAV.
Using local expert knowledge to gather invasive fish presence data was cost-effective and relatively low effort, but there were challenges. First, local expert knowledge is based on personal observations, which can be influenced by individuals’ capacity to recall information and by the time elapsed since the waterbody was last visited. Moreover, respondents likely had varying tenures in their position, which also could have affected their knowledge of presence. We did not attempt to account for any tenure effects, as the relationship between tenure and knowledge, although presumably positive, was not documented. Second, because lakes were randomly selected, fish biologists were unable to provide information for roughly half of the lakes selected because they were on private land or otherwise not managed by their agency. While we could have obtained a larger sample size by collecting information only for lakes that biologists had knowledge about, it was important to select lakes at random to ensure that the sample of study lakes was representative of the breadth of lakes in the LMAV, reducing bias associated with agency management foci, and increasing the reliability of our results. Future studies may incorporate private landowners to obtain presence data for inaccessible lakes. Another drawback of using local expert knowledge is that it does not allow for the collection of absence data. Presence-absence data have been shown to be more reliable than presence-only data because presence-only runs the risk of not detecting fish that are actually present, such as where invasive species occur at low densities (
The potential for widespread carp distribution identified in the LMAV highlights the need for rapid and aggressive management action. The map we produced offers a baseline for which data collection on presences and range expansion can be compared, and a resource for identifying priority management and conservation areas. Our results could be applied in many ways, but we see them used primarily as a first step in conservation planning. For instance, agencies may sample lakes predicted to not support carp to confirm their absence. They may then develop preventative measures against invasion such as barriers or other exclusionary constructions. Alternatively, agencies may target lakes predicted to support carp to confirm presence and then decide what management action would be most beneficial to minimize their effects.
Regardless of carp predicted status, continued monitoring through on-site surveys or local expert knowledge surveys such as ours are necessary to validate and update our estimates before and after conservation strategies are implemented. Modeling bigheaded carp distribution in the hundreds of LMAV oxbow lakes was an important step that can facilitate early detection and can guide on-site monitoring. This study offered a framework for an economical instrument that facilitates an initial assessment of large-scale invasive species monitoring.
We thank Neal Jackson and three anonymous reviewers for constructive evaluations. This publication is a contribution of the Forest and Wildlife Research Center at Mississippi State University. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government. The survey described in this article was organized and implemented by the authors and was not conducted on behalf of the U.S. Geological Survey.
MP contributed to the research conceptualization, sample design and methodology, carried out data collection, performed data analysis and interpretation, ethics approval, and writing both the original draft and reviews and editing. LEM contributed to the research conceptualization, sample design and methodology, data collection, interpretation, ethics approval, funding provision and writing both the original draft and reviews and editing. MRB contributed to data analysis and interpretation, ethics approval, and review and editing. CGD contributed to sample design and methodology, data analysis and interpretation, ethics approval, and review and editing. LMB contributed to sample design and methodology, ethics approval, and review and editing. DR contributed to sample design and methodology, ethics approval, and review and editing.
Funding was provided by the U.S. Fish and Wildlife Service, Grant F21AP02996-00, to the Mississippi Department of Wildlife, Fisheries and Parks. The Mississippi Cooperative Fish and Wildlife Research Unit is jointly sponsored by the Mississippi Department of Wildlife, Fisheries and Parks; Mississippi State University; the U.S. Geological Survey; the U.S. Fish and Wildlife Service; and the Wildlife Management Institute.
All of the data that support the findings of this study are available in the main text or Supplementary material.
Interview introduction script.
I am a graduate student at Mississippi State University. My thesis project is on bigheaded carp distribution patterns in the Lower Mississippi Alluvial Valley. Through a short survey I am gathering data on carp presence and lake characteristics that influence carp distribution. I will recommend management strategies based on the results.
The questions I ask you will be mostly closed ended. For each lake, I will read out the questions and record your responses. I will share a map of the lake to be sure we are talking about the same place. After each round of questions, I will ask if there are any adjacent lakes within a certain area, shown on the map, that you have information about. I will then administer the questionnaire for those lakes as well. You may ask for clarification at any time. I will repeat the questionnaire for as many lakes as we can within our allotted interview time. Do you have any questions before we begin?
Interview questionnaire.
Lake ID:
State and management district:
Effort/Detectability
Carp Presence
Current Management Efforts (to be asked at the end of the interview)
Geo-referenced bigheaded carp presence in lakes throughout the Lower Mississippi Alluvial Valley recorded during interviews with fish biologists
Data type: docx
Pearson’s product-moment correlation coefficients between environmental covariates considered for evaluating bigheaded carp presence in lakes throughout the Lower Mississippi Alluvial Valley
Data type: docx
Relationships between MaxEnt habitat suitability values (y axis) and covariates (x axis; Table 2.1) included in the species distribution model of bigheaded carp presence in lakes throughout the Lower Mississippi Alluvial Valley
Data type: docx