Research Article |
|
Corresponding author: Chiara Facca ( facca@unive.it ) Academic editor: Yuriy Kvach
© 2026 Marco Boschiero, Chiara Facca, Francesco Cavraro, Simone Redolfi Bristol, Anna Gavioli, Federico Riccato, Matteo Zucchetta, Piero Franzoi.
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:
Boschiero M, Facca C, Cavraro F, Redolfi Bristol S, Gavioli A, Riccato F, Zucchetta M, Franzoi P (2026) Ecology and biology of Callinectes sapidus in the Northern Adriatic Sea: could the small spatial scale explain its outbreak? Aquatic Invasions 21(1): 49-72. https://doi.org/10.3391/ai.2026.21.1.180751
|
The Atlantic blue crab (Callinectes sapidus) is among the 100 worst invasive species in the Mediterranean Sea, causing significant ecological and economic impacts. The aim of this study is to investigate key aspects of the species’ biology and ecology during its demographic outbreak in a Northern Adriatic area significantly affected by the species’ invasion. Year-round sampling was carried out across a short spatial gradient encompassing lagoon, estuarine, and marine habitats. This comprehensive approach aimed to elucidate the invasive success of the blue crab. Our findings shows that the species resulted widely distributed across all habitat types, with significant differences among stations, seasons, and sexes, with females being prevalent in higher salinity marine and outer lagoon waters during spawning season. These findings, along with the spatiotemporal analyses of the condition factor and the presence in the lagoon of various cohorts of juveniles over the year, highlight the completion of the complex life cycle of the blue crab on an extremely small spatial scale. Moreover, with an average of over 2 million eggs laid per female and a prolonged spawning season, the species reveals a robust reproductive potential, likely favoured by the short distance between mating and spawning habitats. In conclusion, the results of this study underscore the critical role of the short spatial environmental mosaic in facilitating the invasive success of C. sapidus, providing relevant data for managing this unprecedented demographic explosion.
Atlantic blue crab, fecundity, invasion, reproductive cycle, salinity gradient
The Atlantic blue crab Callinectes sapidus Rathbun, 1896 is a decapod crustacean of the Portunidae family (
In Europe, the blue crab was first recorded in 1900 along the French coasts (
Although the complete eradication of the species is unlikely due to its high adaptability and mobility (
Several factors can be hypothesized to have contributed to the outbreak of C. sapidus that occurred in the Northern Adriatic in 2023. Firstly, increased temperatures associated with climate change appear to favour the progressive expansion of the species to higher latitudes. Indeed, under various RCP scenarios, it has been predicted that in the near future C. sapidus will experience an increase in habitat suitability in the Northern Adriatic (
This study aims to investigate key aspects of the species’ biology and ecology during its demographic outbreak in a Northern Adriatic area. We present an annual insight about C. sapidus across various habitat types (riverine estuary, lagoon and marine coast). The investigation was carried on simultaneously with the species demographic explosion in the area and we focused on: I) Spatial and temporal distribution to determine habitat preferences and year-round migration patterns. II) Main morphological relationships to study growth, weight gain and condition factor across stations and seasons. III) Size structure within lagoon sites using two fishing gears to identify crabs at various maturation stages. IV) Species fecundity by analysing the distribution and fertility of ovigerous females across stations and seasons. This comprehensive approach aimed to elucidate the invasive success of this allochthonous species in an area characterized by different habitats along a short spatial gradient.
The present study was conducted in four locations near Chioggia (northeastern Italy), along a salinity gradient which extends from the outlet of the Brenta River to the Adriatic Sea, including the southern part of Venice lagoon (Figure
In the northwestern Adriatic, salinity gradients occur on an extremely small spatial scale, around two orders of magnitude smaller than in Atlantic environments such as the Chesapeake Bay, where the blue crab is most prevalent and females undergo post-mating migrations for 150–200 km (
The two lagoon stations exhibit salinity levels that generally range from polyhaline to euhaline, and partial macrophyte cover. Despite their geographical proximity, CHN is more confined by saltmarshes and undergoes a minor influence from both freshwater and saltwater inputs into the lagoon. Conversely, CHS is more exposed to freshwater inputs from a large canal (Canale Novissimo) and to seawater inputs from tidal flow through Chioggia inlet. PELL, a marine site, is characterized by higher and more stable salinities and presents a bare sandy substrate. BR, located at the Brenta river mouth, typically presents oligohaline waters and a bare muddy substrate. However, at this station, salinity undergoes significant fluctuations, primarily driven by meteorological conditions. The distribution of the sampling stations was arranged following the salinity gradient on a spatial transect of only a few kilometres.
Sampling activities using wire crab pots were carried out from June 2023 to June 2024. A total of 14 sampling events were performed at each of the four sampling stations during this period (9 in 2023 and 5 in 2024). During each sampling event, six wire pots of two different types (3 pots per type) were deployed at each site and retrieved 24 hours later. The first type measured 52.5×42.5×22.5 cm (L×W×H; 4×4 cm square mesh) and featured four entry funnels, no internal chambers, and three circular escape vents (6 cm diameter) for non-target species. The second type measured 40×61×40 cm (L×W×H; 5×5 cm square mesh) and was equipped with four entry funnels, an internal chamber, and no escape vents. All wire pots were baited with low-value fish species, primarily grey mullets and sardines. One sampling event in 2023 was not conducted at the PELL due to adverse sea conditions. Sampling using the fyke nets (locally named “tresse con cogolli”, a traditional fishing gear with a mesh size of 0.6 cm (
All the crabs sampled with pots and a sub-sample of crabs (representative of the sex-ratio and size distribution) collected during each fyke nets check were culled using a two-step ethical protocol as described in
The methodology for fertility estimates used was adapted from
The distribution of crabs captured with wire pots over stations and seasons by sex was analysed by performing a GLM (negative binomial family; significant p < 0.05; n = 660). Post hoc exploring relevant interactions between the three factors was conducted through Estimated Marginal Means (EMMs) with False Discovery Rate (FDR) adjustment method (significance at p < 0.05). To perform graphical representation and statistical analysis, data were organized into seasonal periods, broadly defined as: Summer (Jun-Aug 2023), Autumn (Sep-Dec 2023), Winter (Feb-Mar 2024), and Spring (Apr-Jun 2024).
The relationships between the main morphological parameters and the total weight of all sampled crabs (wire pots and fyke nets) were compared using regression models, with log-transformed variables. The following comparisons were performed (W-CL; W-CW; W-TCW; TCW-CW; TCW-CL) separately for males and females. Specimens lacking one or both chelae were excluded from the regressions involving weight and eggs weight was removed in ovigerous females. The slope of the linear relationship (parameter ‘b’) was used to determine isometry (b = 1 or b = 3) and allometry (b ≠ 1 or b ≠ 3) by means of linear F-test (significant allometry at p < 0.05) (
Population size distribution was seasonally calculated for females and males using TCW with 1 cm frequency and distinguishing fyke nets data from wire pots data, as the latter only capture adult crabs. Moreover, only data from the lagoon stations were included, as sampling data with both gear types were available. No statistical analysis was performed as the fyke net data presents a polymodal distribution that could not be fitted by a suitable model.
Spatial and temporal distribution of ovigerous females was tested using a binomial presence/absence GLM (significance at p < 0.05; n = 330). The odds ratio was calculated, and Tukey HSD post hoc for pairwise comparisons was performed (significance at p < 0.05). The tests included solely data of crabs collected by wire pots and data organized into seasonal periods as previously reported. As ovigerous females were only sampled from April to September, for further investigations three periods were identified within the egg-laying season: 1 – spring (April and May; n = 33); 2 – early summer (June and July; n = 27); 3 – late summer (August and September; n = 25). A total of 86 ovigerous females were analysed, but one was excluded from statistical analyses as an outlier using the Interquartile Range (IQR) Method. The estimated egg numbers were analysed as a function of sampling station, period, and maturation stage by means of One-Way ANOVA and Tukey post hoc test in case of significant differences (p < 0.05; n = 85). The linear equations, the adjusted coefficient of determination (Adjusted R²) and the p-value (significance at p < 0.05) were obtained. Slopes comparisons were tested by means of ANCOVA analysis (significance at p < 0.05).
Statistical analyses were performed and graphics plotted using R (v 4.5.0; R Core Team 2025) within RStudio (v 2024.12.1+563;
The seasonal and yearly average values of water environmental parameters, detected at each sampling station, are shown in Table
Seasonal and yearly average (± SD) values of water environmental parameters detected per sampling station (BR, CHN, CHS, PELL): Temperature (°C); ph; Salinity (PSU); Dissolved Oxygen (%); Turbidity (FNU).
| Station | Season | Temp. (°C) | pH | Sal. (PSU) | DO (%) | Turb. (FNU) |
|---|---|---|---|---|---|---|
| BR | Summer | 23.15 ± 2.86 | 7.83 ± 0.10 | 3.11 ± 3.61 | 85.01 ± 0.46 | 20.08 ± 6.86 |
| Autumn | 13.56 ± 8.17 | 8.04 ± 0.20 | 15.42 ± 15.46 | 104.19 ± 13.18 | 68.87 ± 97.92 | |
| Winter | 11.74 ± 1.76 | 8.02 ± 0.11 | 2.15 ± 0.91 | 106.38 ± 6.32 | 20.11 ± 0.30 | |
| Spring | 17.76 ± 2.76 | 7.83 ± 0.45 | 0.70 ± 0.80 | 116.39 ± 1.29 | 52.80 ± 59.95 | |
| Full year | 16.37 ± 6.00 | 7.92 ± 0.25 | 5.64 ± 9.54 | 102.32 ± 14.42 | 42.32 ± 56.01 | |
| CHN | Summer | 23.54 ± 1.60 | 8.38 ± 0.35 | 32.85 ± 2.20 | 74.35 ± 15.66 | 3.39 ± 2.36 |
| Autumn | 17.57 ± 8.13 | 8.25 ± 0.20 | 31.60 ± 2.47 | 105.37 ± 18.50 | 1.78 ± 1.09 | |
| Winter | 12.22 ± 1.35 | 8.47 ± 0.07 | 29.98 ± 3.87 | 109.74 ± 7.11 | 3.40 ± 2.40 | |
| Spring | 21.98 ± 4.34 | 8.50 ± 0.30 | 26.21 ± 0.59 | 108.12 ± 51.56 | 25.90 ± 41.98 | |
| Full year | 19.88 ± 5.99 | 8.38 ± 0.26 | 30.66 ± 3.29 | 95.50 ± 28.94 | 7.76 ± 19.26 | |
| CHS | Summer | 23.96 ± 1.47 | 8.39 ± 0.38 | 30.77 ± 5.46 | 81.92 ± 9.25 | 11.52 ± 6.96 |
| Autumn | 17.95 ± 7.95 | 8.30 ± 0.16 | 32.94 ± 2.08 | 104.88 ± 12.98 | 2.52 ± 0.86 | |
| Winter | 12.14 ± 1.34 | 8.43 ± 0.08 | 31.25 ± 3.22 | 110.17 ± 5.64 | 3.91 ± 1.17 | |
| Spring | 21.72 ± 4.45 | 8.43 ± 0.08 | 24.87 ± 3.98 | 102.20 ± 11.97 | 6.92 ± 3.81 | |
| Full year | 20.07 ± 6.00 | 8.38 ± 0.24 | 30.19 ± 4.75 | 96.86 ± 15.14 | 6.87 ± 5.74 | |
| PELL | Summer | 27.11 ± 2.21 | 8.14 ± 0.22 | 33.45 ± 2.16 | 122.01 ± 7.02 | 0.37 ± 0.36 |
| Autumn | 19.56 ± 6.07 | 8.35 ± 0.05 | 35.10 ± 1.56 | 112.65 ± 5.56 | 1.65 ± 1.75 | |
| Winter | 11.92 ± 1.70 | 8.37 ± 0.13 | 34.21 ± 0.54 | 133.17 ± 7.88 | 1.52 ± 0.49 | |
| Spring | 21.37 ± 2.48 | 8.52 ± 0.02 | 29.31 ± 0.83 | 146.12 ± 6.33 | 2.39 ± 2.08 | |
| Full year | 20.54 ± 6.17 | 8.34 ± 0.15 | 33.25 ± 2.68 | 123.99 ± 13.90 | 1.50 ± 1.52 |
Spatial and seasonal distribution of the blue crab, assessed through catch per unit effort (CPUE), demonstrated extreme variability and complexity as shown in Figure
The outputs of the relationships between the main morphological parameters and the total weight of crabs (W-CL; W-CW; W-TCW; TCW-CW; TCW-CL) are summarized and reported in Suppl. material
The relative condition factor (Kn) was calculated based on the regression lines for females, y = 2.249x – 1.0093 (adjusted R² = 0.7944, p < 0.001), and for males, y = 2.797x – 2.130 (adjusted R² = 0.8334, p < 0.001) (Figure
The results of the relative condition factor (Kn) for female and male crabs are shown in Figure
Population size structure analysis (Figure
The distribution of ovigerous females per station and season is shown in Figure
The ovigerous females weight ranged from 112.99 g (12.51 TCW cm) to 308.50 g (17.80 TCW cm). The egg mass weight wet ranged from 10.16 to 55.50 g (average 30.27 ± 9.12 g), representing an average of 15.16 ± 2.69% of the total female body wet weight. The number of estimated eggs ranged from 518,717 to 4,106,743 (average 2,072,450 ± 745,464), but a high of 5,982,200 was recorded and excluded from the analyses as outlier. The average egg diameter of each egg mass ranged from 217.3 to 310.3 µm (average 251.4 ± 17.6 µm). No significant relationship was found between estimated egg numbers and the three periods within the reproduction season (spring, early summer and late summer) (One-Way ANOVA, p = 0.800) (Figure
Average (± SD) egg mass weight (g) estimated egg number and egg diameter (µm) per maturation stage (early, intermediate and pre-hatching) and in total. The number of analysed crabs is also reported (n).
| Eggs Maturation Stage | n | Egg Mass Weight (g) | Estimated Egg Number | Egg Diameter (µm) |
|---|---|---|---|---|
| Early | 33 | 31.08 ± 7.61 | 2,330,196 ± 735,218 | 240.0 ± 11.7 |
| Intermediate | 36 | 30.83 ± 8.96 | 2,034,355 ± 707,586 | 254.2 ± 14.9 |
| Pre-hatching | 16 | 27.31 ± 11.99 | 1,626,563 ± 654,900 | 268.6 ± 17.8 |
| Total | 85 | 30.27 ± 9.12 | 2,072,450 ± 745,464 | 251.4 ± 17.6 |
Relations between the egg mass weight (g), female crab weight (g) and estimated egg number were also explored. Results showed a positive correlation between female weight (W) and the egg mass weight (EggW), following the equation: EggW = 0.1633 × W – 2.1937, with adjusted R2 = 0.6768 and p < 0.001 (Figure
Data of abundance and distribution of C. sapidus in four distinct aquatic habitats (two brackish, one freshwater and one coastal marine) in the area of Chioggia highlighted the magnitude of the outbreak of the species that occurred during spring and summer 2023, and its rapid and successful establishment. The blue crab was indeed found to be abundantly present in all the aquatic environments where samplings were carried out. The high significance of the three-way interaction (station:season:sex) highlights the complex distribution pattern of the species in an area where different aquatic habitats are present nearby and deeply connected between each other, showing that males and females have different space-time dynamic in the habitat use. The overall greater abundance of male specimens in the riverine station confirms their preference for lower salinity waters (
Although the blue crab was found to be abundantly present all over the study area, analysis of the relative condition factor (Kn) reveals significant differences among specimens captured in different stations and seasons. The station has a strong overall effect, as both males and females captured at the riverine site exhibit a poorer condition compared to all other stations. Despite the blue crab is slightly more abundant at the riverine site, the lower Kn might be caused by the high energy consumption for osmoregulation in a oligohaline environment (
Fecundity of the blue crab might also benefit from the longer activity period, being reproductive for almost six months. Coupling our samples with the reports of local fishermen (F. Riccato, personal communication) ovigerous females were collected from late March to mid-September. The length of egg brooding period is consistent with other recent studies conducted in the Mediterranean Sea (
In conclusion, the findings of our study highlight the massive spread of the blue crab C. sapidus in the Northern Adriatic Sea, in an area comprising diverse aquatic environments at short distances: lagoon, river estuary and marine coast, all located between 2 and 9 km from each other. Spatial and seasonal distribution of males and females appeared to be strongly connected with the reproductive cycle of the blue crab. Therefore, the proximity of all these environments represents an extremely crucial factor in favouring the rapid successful establishment of the species, which can take advantage of all the different habitats to foster a faster and more efficient life cycle compared to its native range. Based on our results, the effect of the small spatial scale “ from river to sea” gradient may be equivalent or even stronger than other factors favouring the invasive success of C. sapidus. This study represents a relevant starting point for further investigations in the area; more information are, in fact, necessary to provide a better understanding of the species biological and ecological characteristics in Northern Adriatic, such as its trophic role and importance of other drivers, such as temperature rise, larval connectivity and food availability. While the complete eradication of the species is unlikely due to its high adaptability and mobility, our findings provide one of the first insights that could be helpful for the development of targeted management strategies, such as intensive and localized fishing, and for understanding the invasive success of this allochthonous species on the Northern Adriatic ecosystems.
The authors are grateful to Giuseppe Cherubini (Veneto Region), Mauro Armelao and Walter Salvagno (Fondazione della Pesca di Chioggia) and the fishermen who supported the blue crab harvesting: Andrea Salvagno, Mattia Salvagno, Ondino Crosara, Roberto Scarpa. The Authors are grateful to the Editor and the reviewers for their constructive feedback.
The present data were collected in the framework of the projects listed below with the contribution of the Private and Public Bodies in bold. The funders had no role in data collection and analysis, and preparation of the manuscript. All the listed funders participated in the study design and approved the decision to submit the manuscript for publication.
Marco Boschiero: sample design and methodology, investigation and data collection, data analysis and interpretation, writing – original draft; Chiara Facca: sample design and methodology, writing – review & editing; Francesco Cavraro: sample design and methodology, writing – review & editing; Simone Redolfi Bristol: investigation and data collection, writing – review & editing; Anna Gavaioli: writing – review & editing; Federico Riccato: investigation and data collection, writing – review & editing; Matteo Zucchetta: data analysis and interpretation, writing – review & editing; Piero Franzoi: Research conceptualization, sample design and methodology, writing – review & editing, funding provision.
All research pertaining to this article did not require ethics approval. All research pertaining to this article was allowed by the research permit Decreto n. 50 - 24/01/2023 issued by Veneto Region (Italy).
Species georeferenced records are available at the European Alien Species Information Network: https://easin.jrc.ec.europa.eu/easin/RJD/Download/37b1faca-07d0-46d9-a05c-cd4aaf654bbf.
Suppl. table SS1
Data type: xlsx
Explanation note: Regression equations with corresponding R2 values and p-values for all pairs of morphometric parameters tested for females and males; (b) coefficients are reported in red for negative allometry, in green for positive allometry and in blue for isometry. The results of the ANCOVA test, which compared the slopes of the regression lines between sexes, are also provided.
Suppl. table SS2
Data type: xlsx
Explanation note: Min, max and mean ± SD values of blue crab’s total carapace width (TCW) caught with wire pots and fyke nets across the four sampling seasons
Residuals and normality diagnostics of the most complex GLM
Data type: png
Explanation note: figure S1: Residuals and normality diagnostics of the most complex GLM A) negative binomial family GLM for crabs’ distribution; B) gamma family GLM for the relative condition factor; C) binomial presence/absence GLM for ovigerous females’ distribution.