Human-driven breakdown of predator–prey interactions in the northern Adriatic Sea
Abstract
Long-term baseline data that allow tracking how predator–prey interactions have responded to intensifying human impacts are often lacking. Here, we assess temporal changes in benthic community composition and interactions between drilling predatory gastropods and their molluscan prey using the Holocene fossil record of the shallow northern Adriatic Sea, which is characterized by a long history of human transformation. Molluscan assemblages differ between the Isonzo and Po prodelta, but both show consistent temporal trends in the abundance of dominant species. Samples of mollusc prey collected at high stratigraphic resolution indicate that drilling frequencies have drastically declined in the Po prodelta since the mid-twentieth century, while a weaker trend in the more condensed sediments of the Isonzo prodelta is not statistically significant. The decrease in drilling predation intensity and the community turnover are linked to the loss of predatory gastropods and the increased relative abundance of less-preferred prey during the most recent decades. Our results align with data showing the substantial depletion of marine resources at higher trophic levels in the region and indicate that the strong simplification of the food web initiated in the late nineteenth century accelerated further since the mid-twentieth century.
1. Introduction
Global warming, ocean acidification, eutrophication and direct human interventions in marine ecosystems such as fishing, bottom trawling and species introduction markedly change ecosystem functioning and influence biotic interactions [1–3]. The removal of top vertebrate predators due to overfishing results in the loss of top-down control [4,5] and marine invertebrates at lower trophic levels often profit from this predator release [6–8]. However, for benthic invertebrates the effect of anthropogenic impacts on predator–prey interactions at lower trophic levels is poorly known over decadal to centennial timescales. This gap is driven by the lack of long-term time series pre-dating the onset of anthropogenic disturbances [9,10]. However, ichnological signatures of interspecific interactions preserved in the youngest fossil record (i.e. traces of predation) can be used to reconstruct recent shifts in predation pressure and to place those shifts in a broader historical context [9–13]. In particular, drill holes produced by gastropods in shells of other molluscs provide direct evidence of predation [14] and are one of the most widely used sources of data on predator–prey interactions in the fossil record [15,16], accounting for as much as 75% of predation occurrences documented in the palaeontological literature [17].
Here, we use the frequency of predatory traces preserved in the sedimentary record of the northern Adriatic Sea (NAS; figure 1a ) to reconstruct changes in the frequency of predatory attacks of shell-drilling gastropods on their molluscan prey since the late nineteenth century and to assess the magnitude of those changes relative to the late Holocene baseline pre-dating the onset of intense anthropogenic activities. Over the last 10 000 years, the intensity of drilling predation increased across the basin in step with the increase in organic enrichment and water depth caused by the Holocene transgression [18]. However, the NAS also ranks among the most strongly altered marine regions worldwide, with anthropogenic impact on coastal ecosystems reaching back beyond Roman times and rapidly intensifying over the last century [19–21]. Since the mid-twentieth century, the benthic community composition has changed over the entire northern Adriatic basin [22–25] in response to persistent anthropogenic pressures, including bottom trawling, eutrophication and hypoxia [26–30]. Molluscan communities were subjected to a basin-wide loss of epifaunal species, a reduction of grazers, carnivores and herbivores, and a concurrent increase in the dominance of infaunal and opportunistic species feeding on plankton or detritus [31]. A case in point is the opportunistic, infaunal, filter-feeding bivalve Varicorbula gibba, which increased in abundance and body size during the last decades [32]. This regional-scale population shift, unprecedented relative to the Holocene history of this species, was hypothesized to be driven by ecological release from predation, competition and trophic amensalism [25].
Here, we quantify changes in the composition of the molluscan fauna and drilling frequencies (DFs) in sediment cores to document temporal trends in predator–prey interactions between shell-drilling gastropods and their molluscan prey across the Holocene–Anthropocene transition in the NAS. Specifically, we evaluate the impact of human-induced loss of predatory gastropods and other community-level compositional changes on predator–prey interactions.
2. Material and methods
Sediment cores were collected at three stations in the NAS in 2013, with two stations located in the Po prodelta at 21 m water depth (Po 3, Po 4), and one station in the Isonzo prodelta in the Gulf of Trieste at 12 m water depth (Panzano) (figure 1a , electronic supplementary material, table S1). These two deltaic systems are characterized by the highest sediment accumulation rates in the northern Adriatic Sea [33,34] and thus provide high, decadal-scale resolution of stratigraphic archives [25,35]. At all three stations, the seafloor is muddy. At the Po sites, mud mounds generated by the thalassinid shrimp Jaxea nocturna and by the polychaete Sabella pavonina are abundant [35]. At Panzano, epifaunal clumps formed by sponges, ascidians, ophiurids, polychaetes and bivalves dominate [36]. Changes in community composition and DF were assessed in two 16 cm diameter replicate cores at each station [37], which were taken a few metres from each other (electronic supplementary material, table S1). These replicate cores are very similar in faunal composition [35,38] and were pooled in this study to increase sample size per increment when analysing DFs and community composition. The cores were sliced into 2 cm increments for the uppermost 20 cm, and 5 cm increments below. The smaller increments at the top were implemented to potentially resolve more finely the surface strata, which record the time interval during which anthropogenic impacts intensified. However, for the analyses presented here, the adjacent 2 cm increments were pooled to increase sample size for individual time intervals and obtain a better comparability with the 5 cm increments used for the rest of the core [35,38].
The core age models are based on radiocarbon-calibrated amino acid racemization of shells of V. gibba from cores M13 and M21 in the Po delta [35] and core M28 in the Isonzo prodelta [38]. Time averaging of molluscan assemblages in individual increments varies between 10 and 50 years at the Po prodelta and between approximately 20 years and a few centuries at the Isonzo prodelta (figure 6 in Ref. [35]). Four distinct temporal units were distinguished using the age models based on V. gibba, including the middle to late Holocene with 40 samples, the seventeenth to early nineteenth century with 21 samples, the late nineteenth to early twentieth century with 32 samples and the late twentieth to twenty-first century, with 42 samples.
Sedimentation rates at the Panzano station are approximately 2–4 mm/yr and the sedimentary record of the cores covers approximately the last 500 years of ecosystem history [39], providing baseline data on variability in DF and community composition before the onset of strong human impacts in the twentieth century. Conversely, the sediment cores from the Po prodelta provide a higher stratigraphic resolution but extend back only to the late nineteenth century due to high sedimentation rates of approximately 10–20 mm/yr [35]. Therefore, we included older samples from three sediment cores (205-S9, EM-S7 and EM-S13) drilled on the coastal Po plain near Ferrara to evaluate the significance of the recent ecological changes in the context of longer temporal trends. These older samples represent middle to late Holocene prodelta sediments (1.1 to 7.0 ky BP [40]), deposited under similar environmental conditions but subsequently buried under the prograding delta front and coastal plain deposits, and thus currently located at 19 to 26 m below the surface of the present-day Po plain.
All samples were sieved using a 1 mm mesh size and only complete molluscan shells and valves (>90% preserved) were counted and identified to the species level. The abundance of bivalves was calculated as the number of articulated specimens plus the combined number of left and right valves divided by two. To estimate DF (i.e. the per cent of prey specimens with complete drill holes representing a proxy for the frequency of successful predatory attacks), the total number of drilled shells was divided by the total number of bivalve, gastropod and scaphopod individuals per core increment. Because only one valve in bivalves is usually drilled, DF was corrected using the equation of [41]: DF = DV/[(RV + LV)/2 + A], where DV is the total number of drilled valves, RV is the total number of right valves, LV the total number of left valves and A the number of articulated individuals. Only samples with at least 20 individuals were used in the analyses. In addition to the DF of the total assemblage we also calculated DFs for bivalves, gastropods and the four most abundant species, including the bivalves V. gibba and Kurtiella bidentata and the gastropods Turritellinella tricarinata and Tritia varicosa (electronic supplementary material, figure S1, table S2). The four species, however, are not abundant enough in temporal bins to make a rigorous comparison of DFs within and between cores at the species level (electronic supplementary material, table S3). Whereas high-resolution analyses of DFs are not possible for individual taxa, we supplemented assemblage-level analyses with coarse-resolution assessments of the four common species. In that case, core samples were pooled into pre-anthropogenic (middle-late Holocene, seventeenth to early nineteenth century, late nineteenth to early twentieth century) and post-anthropogenic (late twentieth to early twenty-first century) and Fisher’s test was applied to assess for statistical differences between the two sample groups. Predatory gastropods, which make up only 1.38% of the total molluscan assemblage, are strongly dominated by two species of naticids, Euspira nitida (70%) and E. macilenta (20%), while the few muricids (9%) are mostly juveniles of Hexaplex trunculus and Bolinus brandaris. We were not able to distinguish consistently between drill holes made by these two families of shell-drilling gastropods.
We evaluated stratigraphic changes in the composition of molluscan assemblages using non-metric multidimensional scaling ordination (NMDS; k = 2 dimensions) based on square-root-transformed proportional species abundances and Bray–Curtis dissimilarities. The differences in molluscan community composition between the four temporal units were tested with non-parametric permutational multivariate analysis of variance (PERMANOVA [42]). Differences in DFs between the units were assessed with the Kruskal–Wallis test and the pairwise Wilcoxon rank-sum test, using a posteriori Bonferroni correction for multiple comparisons. Distance-based redundancy analysis was used to assess whether differences in species-level community composition are related to differences in DF [43]. The principal coordinate analysis using Bray–Curtis dissimilarity was performed on the square-root-transformed proportional species abundances, and the resulting ordination scores were subjected to redundancy analysis, using the function ‘capscale’ in the vegan R package [44].
To assess stratigraphic trends in the relative abundance of drilling predators, a rank correlation between the relative abundance of predators and depth in the core was compared with a resampling distribution of Spearman rho values under a null model in which the relative abundance of predators is constant along the core. The resampling distribution is based on a Monte Carlo simulation of vertical trends generated by random drawing of drilling predators with probability given by their overall relative frequency in cores. At each core level, the total number of simulated specimens was based on the actual sample size for this level. For each simulation run, a Spearman rho value was calculated for random data. Positive rho values indicate an upward decline and negative values indicate an upward increase in predator abundance towards the top of the core. The two-tailed significance value p is based on the number of replicate runs in which the simulated absolute values of rho exceeded the observed absolute rho value. The results are based on 10 000 iterations.
To assess trends in DF, a serial Monte Carlo model was used based on the DF data in a given core. In each iteration of the model, a random walk was simulated with an initial DF value based on the value observed at the base of the core (when the starting value was selected randomly from among all values observed in the core, the results (data not shown) are effectively the same). Then, values were shifted by adding one first difference, using a first difference value observed between two randomly selected adjacent levels. The process was repeated for all successive levels resulting in a random time series based on the actual first differences observed in the core. If the shift resulted in an impossible negative value of DF, the shift was reversed. Alternative runs with negative values converted to 0 (not shown) yielded comparable outcomes. The process was repeated 10 000 times and the resulting set of 10 000 time series was used to estimate the range of possible trends expected under a random walk. The results were also used to estimate significance of changes in DF across the core. In this approach, the core was divided into two groups of increments (upper and lower parts of the cores) starting from the third increment in the core and then successively shifting the boundary between the lower and upper group to find the dividing point that results in the maximum difference in DF. The difference was measured using t-statistic (difference between the mean DFs between the two parts of the core divided by the standard error). The resulting trend was then compared against 95% confidence bands based on the analysis of the 10 000 Monte Carlo series simulating non-directional random trends. In addition, a Spearman rank correlation was computed for DF and core depth (positive correlation indicates an upward decrease in drilling predation and negative correlation indicates an upward increase in DF). The distribution of expected values of Spearman rho was computed based on 10 000 Monte Carlo runs. Monte Carlo simulations and Spearman rank correlations were performed for the total assemblage, gastropods and bivalves at all three stations (except for gastropods at station Po 3 because of a low number of samples with a minimum of 20 individuals; electronic supplementary material, table S3). Statistical analyses were performed with the software package R (v. 4.3.1). All data used in this study and the R code are available at the Dryad Digital Repository [45].
3. Results
The overall DFs in the molluscan assemblages of the Po prodelta exceeded 10% in all but one sample from the middle-late Holocene and fluctuated around 10% in the late nineteenth to early twentieth century (figure 2). The DFs in the Isonzo prodelta were consistently higher, varying between 20% and 30% from the seventeenth to the early twentieth century, rarely exceeding 30%. Despite these differences in DFs, both regions exhibited peaks in DFs in the mid-twentieth century (above 20% in the Po prodelta and above 30% in the Isonzo prodelta), followed by a decrease during the late twentieth and early twenty-first century, with DFs reaching all-time lows in the youngest samples (figures 1b and 2). These declines were also observed when bivalves and gastropods were analysed separately (electronic supplementary material, figures S2 and S3). DFs in the total assemblage during the late nineteenth to early twentieth century and the late twentieth century were significantly lower compared with the late Holocene baseline in the Po delta when core increments are grouped into three multi-decadal time bins according to their stratigraphic position (figure 1, table 1a,b ). Differences between time bins are statistically insignificant in the Isonzo delta (figure 1, table 1a ).
a | |||
---|---|---|---|
chi-squared | d.f. | p‐value | |
Po prodelta | 15.094 | 2 | 0.0005276 |
Isonzo prodelta | 4.5819 | 2 | 0.1012 |
b | |||
Po prodelta | late Holocene | late nineteenth to early twentieth century | |
late nineteenth to early twentieth century | 0.0098 | ||
late twentieth to early twenty-first century | 0.0012 | 0.2567 | |
c | |||
chi-squared | d.f. | p‐value | |
Po prodelta | 17.713 | 2 | 0.0001425 |
Isonzo prodelta | 2.6181 | 2 | 0.2701 |
d | |||
Po prodelta | late Holocene | late nineteenth to early twentieth century | |
late nineteenth to early twentieth century | 0.01187 | ||
late twentieth to early twenty-first century | 0.00012 | 1 | |
e | |||
PERMANOVA | F-statistics | Pr(>F) | |
regions (Po versus Panzano) | 32.066 | 0.001 | |
Po prodelta temporal units | 11.207 | 0.001 | |
Isonzo prodelta temporal units | 7.2275 | 0.001 | |
f | |||
RDA | F-statistics | Pr(>F) | |
Po | 2.1656 |
|
|
Panzano | 1.5101 | 0.079 |
DFs of the total assemblage exhibit a strong positive autocorrelation among adjacent 5 cm thick levels (lag 1; electronic supplementary material, figures S4–S6). The simulated time series of DFs represent non-directional random walks, which tend to drift towards larger values. The Monte Carlo simulations suggest that the most likely transition from higher to lower DFs occurs around the early 1950s (core depth 15 cm) in Panzano, around 1980s (core depth 50 cm) at Po 3 and at the end of the twentieth century (core depth 22.5 cm) at Po 4 (electronic supplementary material, figures S4–S6). The decrease in DFs of the total assemblage towards the top of the cores is statistically significant for Po 3 and Po 4, but not for Panzano (electronic supplementary material, figure S7). Monte Carlo simulations for DFs of bivalves and gastropods alone show similar trends at the two Po stations, but at Panzano DFs of gastropods also drastically and significantly decrease. At the same time, there is no significant trend for bivalves (electronic supplementary material, figures S7–S12). Also, for the four most abundant species (V. gibba, K. bidentata, T. tricarinata, T. varicosa), the results of Fisher’s test consistently indicate that DF decreased through time and, in three out of four cases, the observed difference was statistically significant (p < 0.0001 in all three cases; electronic supplementary material, table S4).
Drilling predators (naticid and muricid gastropods) constitute only 0.3% and 0.2% of the total molluscan assemblage at the Panzano and Po stations, respectively (figures 3 and 4, electronic supplementary material, figure S13). A noticeable decline in their relative abundance can be detected in the Po cores (figure 3a ) and becomes more noticeable when the trend is smoothed using a moving average with a step of 5 (figure 3b ). Conversely, no decline in drilling predators is discernible at Panzano (figure 4a,b ). The upward decline of drilling predators is statistically significant at the Po stations when measured as rank correlation between core depth and relative abundance of predators (rho = 0.47), as indicated by a Monte Carlo simulation of a null model that assumes constant relative abundance of predators (figure 3c ), with the observed trend reflecting sampling effects only. The decline in DF in Panzano cores is not distinguishable from a Monte Carlo null model (figure 4c ). DFs and abundance of predators decline in parallel at both prodelta successions (figure 1b,c , electronic supplementary material, figure S14).
Molluscan assemblages differ between the two regions (electronic supplementary material, figure S15), but show significant temporal changes in composition accompanied by decreasing DFs in both prodelta successions (table 1e ). However, the correlation between DF and multivariate assemblage composition is significant only in the Po prodelta (figure 5, table 1f ). The changes in relative abundances of the four most abundant species (representing 37–85% of individuals per sample) are related to species-level trends in DFs, highlighting the link between compositional shifts and overall DF. In the Po prodelta, the strong increase in the relative abundance of the opportunistic bivalve V. gibba since the mid-twentieth century is paralleled by the strong decrease of the suspension-feeding gastropod T. tricarinata towards the present. The distinct increase in abundance of V. gibba coincides with a decrease in its DF, while T. tricarinata is never abundant enough in the Po prodelta to confidently evaluate temporal trends in its species-level DF (figure 6a ). Given that V. gibba dominated the molluscan assemblages in that area, the overall DFs followed the trends observed in this species. In the Isonzo prodelta, the relative abundance of the commensal bivalve K. bidentata increased steadily until the late nineteenth–early twentieth century and declined rapidly during the twentieth century. The abundance of the scavenger T. varicosa shows an opposite pattern with a strong increase during that time, while T. tricarinata decreases in relative abundance since the seventeenth century, with a drastic drop during the late twentieth century (figure 6b ). The trends in the relative abundance of K. bidentata, T. tricarinata and T. varicosa in the Isonzo prodelta affect the DFs of the total assemblage: K. bidentata and T. tricarinata are frequently drilled and disappear from the community at the expense of T. varicosa, which is rarely drilled (figure 6b ). V. gibba in the Isonzo prodelta, however, shows a strong increase in relative abundance since the early twentieth century, which comes along with a slight decrease in DFs (figure 6b ). Therefore, the community changes in the Isonzo prodelta have only a minor impact on the overall decreasing DFs of the total assemblage (figure 2, table 1f ).
4. Discussion
Predation is one of the key factors shaping marine ecosystems today [45–48]. However, the rarity or lack of baseline data prohibits assessing how predator–prey interactions may have changed under increasing anthropogenic impact on marine ecosystems (e.g. [13]). A longer term perspective, incorporating data from the Holocene fossil record is necessary to answer this question (e.g. [49]). At a regional scale and using sediment cores from areas with low sedimentation rates, DFs in the NAS were significantly lower during the early phase of the Holocene sea-level rise (approximately 10–20% per increment) than during the sea-level highstand over the last few 1000 years (approximately 20–30% per increment). DF largely increases with increasing clay content of sediments and other environmental variables correlated with water depth (organic enrichment and hydrodynamic energy) and with increasing predation pressure (as reflected in the relative abundance of naticid and muricid gastropods) [18]. However, it remains unclear how the approximately twentieth-century compositional replacement of epifaunal grazers, carnivores and herbivores by infaunal opportunistic detritus and suspension feeders [31] influenced predator–prey interactions. We demonstrate here that this mid-twentieth-century compositional shift is associated with a decline in DFs. This decline coincides with a decline in abundance of predatory gastropods in the Po prodelta but not in the Isonzo prodelta. The decline is directly preceded by distinct peaks in predation frequency (figure 2) coinciding with the onset of eutrophication [26,50]. Depending on the number of trophic levels, primary productivity driven by eutrophication can increase the abundance of predators and predator–prey interactions [51,52]. Over millennial time scales, an increase in organic enrichment of sediments (covarying with a decline in grain size) during the Holocene also had a positive effect on DFs in the NAS [18]. However, the positive effects of enhanced food supply on macrobenthic communities were apparently overwhelmed by the more frequent occurrence of seasonal hypoxia since the 1960s [26], which resulted in a loss of formerly abundant, hypoxia-sensitive species such as T. tricarinata [35]. A similarly complex relationship between the intensity of drilling predation and primary productivity, modulated by the oxygen availability, was documented in bivalve assemblages sampled along the northern Gulf of Mexico: DFs increased with net primary productivity but the development of hypoxic conditions in the most productive areas adjacent to the Mississippi River apparently limited the activity of drilling predators [53]. Although a decline in salinity might also decrease the abundance of drilling gastropods, salinity did not markedly increase or decrease in the northern Adriatic Sea between 1911 and 1982 and thus does not account for the observed decline [26,54].
Changes in assemblage-level DFs are not only dependent on the identity and abundance of the predators but also on the varying relative abundances of prey species and differences in their susceptibility to drilling predation [55]. Our data suggest that both mechanisms were likely responsible for the late twentieth-century decrease in DFs observed in the NAS. The collapse of predator–prey interactions in the Po prodelta was driven by the complete loss of predatory gastropods in the late twentieth and early twenty-first century (figure 1c ), while a prominent prey species, V. gibba, remained abundant. Conversely, the decrease in predator abundance in the late twentieth and early twenty-first century is small and the decrease of DFs in the Isonzo prodelta is only significant for gastropods (electronic supplementary material, figure S11). The decrease in DFs in the Isonzo prodelta coincides with the replacement of the preferentially drilled species (K. bidentata and T. tricarinata) by species that are only moderately affected by drilling predation such as T. varicosa. As the relative abundance of predators did not decrease, it is possible that they may have switched to a non-drilling predatory strategy [56]. Trends in the relative abundance of predators can also mask decline in their absolute abundance, if the population density of other mollusc species declined in the twentieth century as well. However, Panzano assemblages are time-averaged to centuries rather than to decades (i.e. their faunal composition is more condensed) and thus have smaller potential to record any decline in drilling predators occurring at decadal scales during the twentieth century. Our results could be biased by differences in transportability and thus accumulation of drilled and undrilled shells [57], preferential destruction of undrilled shells by crushing predators like crabs or fish [58], and preferential loss of drilled shells by post-mortem breakage due to reduced strength [59]. However, post-mortem transport of shells was not observed in the studied low-energy muddy environments of the NAS and taphonomic assessment of the fossil assemblages gives no evidence of preferential breakage of drilled shells. Moreover, given the stable water depths at Po and Panzano stations during the late Holocene, these potential biases should be similar for all samples and therefore cannot affect the trend of declining DFs. Biased preservation of drilled shells due to the activity of burrowers is also unlikely. In the Po cores, DF declines in tandem with a decrease in bioturbation intensity related to the onset of frequent hypoxia [35] and it remains low when burrowing increases again in the top of the cores following the reduction in the intensity of hypoxic events in the twenty-first century [29]. Also, the formerly and newly common prey taxa at Panzano have similar preservation potential and belong to families with an excellent Cenozoic fossil record [60].
Ecological release refers to the response of populations to the relaxation of selection acting on some ecologically important traits such as body size and abundance, providing ecological opportunity for an increase in population biomass of prey [61]. The increase in eutrophication and frequency of hypoxia since the mid-twentieth century could produce an ecological opportunity for the hypoxia-tolerant species V. gibba [62–65] by removing its competitors and predators, leading to a marked increase in its shell size and relative abundance [25]. We support this hypothesis, showing that twentieth-century ecosystem changes are associated with a decline in the abundance of the hypoxia-sensitive T. tricarinata (figure 6) [66], a decline in the abundance of drilling gastropods (figure 1c , electronic supplementary material, figure S15) and a decline in DFs in the total assemblage (figure 1b ).
Although surface death assemblages are valuable sources of historical ecological data [67], in areas with high sedimentation rates only historical layers below the mixed layer can serve as an unequivocal pre-impact ecological baseline [68]. Whereas many studies reconstructed past community states and Anthropocene impacts from sediment cores (e.g. [69–71]), only a few were evaluated for temporal changes in biotic interactions through the Holocene, and those previous analyses focused on parasite–host interactions rather than predation [49,72–74]. Although studies on marine invertebrate communities enabled the discovery of fundamental processes and mechanisms driving community structure and function (e.g. [75–77]), anthropogenic impacts on predator–prey interactions were documented for top vertebrate predators only (e.g. [4]). Our results from the NAS suggest that, although the region has been affected by humans at least since the Roman times [21], it was not until the twentieth century that human activities began to modify invertebrate predator–prey interactions strongly. This time lag is likely related to a shift in the nature and intensity of anthropogenic impacts with significant changes to water quality of the Adriatic Sea starting only in the nineteenth century [21]. The late twentieth-century increase in eutrophication-driven hypoxia resulted in a strong decline in the frequency of predator–prey interactions at lower trophic levels due to a decrease in predator abundance, a turnover towards less-preferred prey organisms and a size increase in a dominant prey species. These changes were accompanied by the strong depletion of marine resources in the study area, most notably large predators and consumers since the nineteenth century [21,30] and the basin-wide infaunalization of benthic soft-bottom communities driven by anthropogenic habitat degradation in the twentieth century [31]. All these trends indicate a strong simplification of the food web in the NAS.
Ethics
This work did not require ethical approval from a human subject or animal welfare committee.
Data accessibility
All data used in this study and the R code are available at the Dryad Digital Repository [45].
Declaration of AI use
We have not used AI-assisted technologies in creating this article.
Authors’ contributions
M.Z.: conceptualization, data curation, formal analysis, funding acquisition, investigation, methodology, project administration, writing—original draft; R.N.: conceptualization, formal analysis, investigation, methodology, writing—review and editing; M.D.: data curation, investigation, writing—review and editing; I.G.: data curation, investigation, writing—review and editing; A.H.: data curation, investigation, writing—review and editing; M.K.: conceptualization, formal analysis, funding aquisition, investigation, methodology, writing—review and editing; D.S.: conceptualization, formal analysis, funding aquisition, investigation, writing—review and editing; S.W.: data curation, investigation, writing—review and editing; A.T.: conceptualization, formal analysis, investigation, methodology, writing—review and editing.
All authors gave final approval for publication and agreed to be held accountable for the work performed therein.
Conflict of interest declaration
We declare we have no competing interests.
Funding
This study was funded by the Austrian Science Fund (FWF-project P24901), by the Italian Ministry of University and Research under the PRIN 2022 project number 2022WEZR44, by the Slovak Agency for Research and Development (APVV22-0523) and Slovak Scientific Grant Agency (VEGA2/0106/23), and by the US National Science Foundation grants EAR-0920075 and EAR-1559196.
Acknowledgements
The authors thank the captain of the coring vessel, J. Sedmak. The comments and suggestions of three anonymous reviewers improved the manuscript.
References
- 1.
Veale L, Hill A, Brand A . 2000 An in situ study of predator aggregations on scallop (Pecten maximus (L.)) dredge discards using a static time-lapse camera system. J. Exp. Mar. Biol. Ecol. 255 , 111–129. (doi:10.1016/s0022-0981(00)00295-1) Crossref, PubMed, Web of Science, Google Scholar - 2.
Trussell GC, Ewanchuk PJ, Matassa CM . 2006 Habitat effects on the relative importance of trait- and density-mediated indirect interactions. Ecol. Lett. 9 , 1245–1252. (doi:10.1111/j.1461-0248.2006.00981.x) Crossref, PubMed, Web of Science, Google Scholar - 3.
Oug E, Sundet JH, Cochrane SKJ . 2018 Structural and functional changes of soft-bottom ecosystems in northern fjords invaded by the red king crab (Paralithodes camtschaticus). J. Mar. Syst. 180 , 255–264. (doi:10.1016/j.jmarsys.2017.07.005) Crossref, Web of Science, Google Scholar - 4.
Heithaus MR, Frid A, Wirsing AJ, Worm B . 2008 Predicting ecological consequences of marine top predator declines. Trends Ecol. Evol. 23 , 202–210. (doi:10.1016/j.tree.2008.01.003) Crossref, PubMed, Web of Science, Google Scholar - 5.
Estes JA . 2011 Trophic downgrading of planet Earth. Science 333 , 301–306. (doi:10.1126/science.1205106) Crossref, PubMed, Web of Science, Google Scholar - 6.
Duggins DO, Simenstad CA, Estes JA . 1989 Magnification of secondary production by kelp detritus in coastal marine ecosystems. Science 245 , 170–173. (doi:10.1126/science.245.4914.170) Crossref, PubMed, Web of Science, Google Scholar - 7.
Worm B, Myers RA . 2003 Meta-analysis of cod–shrimp interactions reveals top-down control in oceanic food webs. Ecology 84 , 162–173. (doi:10.1890/0012-9658(2003)084[0162:MAOCSI]2.0.CO;2) Crossref, Web of Science, Google Scholar - 8.
Steneck RS, Vavrinec J, Leland AV . 2004 Accelerating trophic-level dysfunction in kelp forest ecosystems of the western North Atlantic. Ecosystems 7 , 323–332. (doi:10.1007/s10021-004-0240-6) Crossref, Web of Science, Google Scholar - 9.
Cintra-buenrostro CE, Flessa KW, Guillermo AS . 2005 Who cares about a vanishing clam? Trophic importance of Mulinia coloradoensis inferred from predatory damage. Palaios 20 , 296–302. (doi:10.2110/palo.2004.p04-21) Crossref, Google Scholar - 10.
Barclay KM, Leighton LR . 2022 Predation scars reveal declines in crab populations since the pleistocene. Front. Ecol. Evol. 10 , 810069. (doi:10.3389/fevo.2022.810069) Crossref, Web of Science, Google Scholar - 11.
Chiba T, Sato SI . 2014 Invasion of a naticid predator and associated changes in death assemblages of bivalve prey in northern Japan: implications for palaeoecological studies. Lethaia 47 , 4–14. (doi:10.1111/let.12033) Crossref, Web of Science, Google Scholar - 12.
Smith JA, Dietl GP . 2016 The value of geohistorical data in identifying a recent human‐induced range expansion of a predatory gastropod in the Colorado River delta, Mexico. J. Biogeogr. 43 , 791–800. (doi:10.1111/jbi.12644) Crossref, Web of Science, Google Scholar - 13.
Smith JA, Handley JC, Dietl GP . 2018 Effects of dams on downstream molluscan predator–prey interactions in the Colorado River estuary. Proc. R. Soc. B. 285 ,20180724 . (doi:10.1098/rspb.2018.0724) Link, Web of Science, Google Scholar - 14.
Kowalewski M . 2002 The fossil record of predation: an overview of analytical methods. Paleontol. Soc. Pap. 8 , 3–42. (doi:10.1017/S1089332600001030) Crossref, Google Scholar - 15.
Kowalewski M, Dulai A, Fürsich FT . 1998 A fossil record full of holes: the phanerozoic history of drilling predation. Geology 26 , 1091–1094. (doi:10.1130/0091-7613(1998)026<1091:AFRFOH>2.3.CO;2) Crossref, Web of Science, Google Scholar - 16.
Harper EM, Peck LS . 2016 Latitudinal and depth gradients in marine predation pressure. Glob. Ecol. Biogeogr. 25 , 670–678. (doi:10.1111/geb.12444) Crossref, Google Scholar - 17.
Klompmaker AA, Kelley PH, Chattopadhyay D, Clements JC, Huntley JW, Kowalewski M . 2019 Predation in the marine fossil record: studies, data, recognition, environmental factors, and behavior. Earth Sci. Rev. 194 , 472–520. (doi:10.1016/j.earscirev.2019.02.020) Crossref, Web of Science, Google Scholar - 18.
Zuschin M, Nawrot R, Dengg M, Gallmetzer I, Haselmair A, Wurzer S, Tomašových A . 2022 Scale dependence of drilling predation in the holocene of the northern Adriatic Sea across benthic habitats and nutrient regimes. Paleobiology 48 , 462–479. (doi:10.1017/pab.2022.6) Crossref, Web of Science, Google Scholar - 19.
Jackson JBC . 2001 Historical overfishing and the recent collapse of coastal ecosystems. Science 293 , 629–637. (doi:10.1126/science.1059199) Crossref, PubMed, Web of Science, Google Scholar - 20.
Lotze HK . 2006 Depletion, degradation, and recovery potential of estuaries and coastal seas. Science 312 , 1806–1809. (doi:10.1126/science.1128035) Crossref, PubMed, Web of Science, Google Scholar - 21.
Lotze HK, Coll M, Dunne JA . 2011 Historical changes in marine resources, food-web structure and ecosystem functioning in the Adriatic Sea, Mediterranean. Ecosystems 14 , 198–222. (doi:10.1007/s10021-010-9404-8) Crossref, Web of Science, Google Scholar - 22.
Simonini R, Ansaloni I, Bonvicini Pagliai AM, Prevedelli D . 2004 Organic enrichment and structure of the macrozoobenthic community in the northern Adriatic Sea in an area facing Adige and Po mouths. ICES J. Mar. Sci. 61 , 871–881. (doi:10.1016/j.icesjms.2004.06.018) Crossref, Web of Science, Google Scholar - 23.
Occhipinti-Ambrogi A, Savini D, Forni G . 2005 Macrobenthos community structural changes off cesenatico coast (Emilia Romagna, northern Adriatic), a six-year monitoring programme. Sci. Total Environ. 353 , 317–328. (doi:10.1016/j.scitotenv.2005.09.021) Crossref, PubMed, Web of Science, Google Scholar - 24.
Mavrič B, Orlando-Bonaca M, Bettoso N, Lipej L . 2010 Softbottom macrozoobenthos of the southern part of the Gulf of Trieste: faunistic, biocoenotic and ecological survey. Acta. Adriatica. 51 , 203–216. Google Scholar - 25.
Tomašových A . 2020 Ecological regime shift preserved in the anthropocene stratigraphic record. Proc. R. Soc. B 287 ,20200695 . (doi:10.1098/rspb.2020.0695) Link, Web of Science, Google Scholar - 26.
Justić D . 1991 Hypoxic conditions in the northern Adriatic Sea: historical development and ecological significance. Geol. Soc. Lond. Spec. Publ. 58 , 95–105. Crossref, Google Scholar - 27.
Degobbis D, Precali R, Ivancic I, Smodlaka N, Fuks D, Kveder S . 2000 Long-term changes in the northern Adriatic ecosystem related to anthropogenic eutrophication. Int. J. Environ. Pollut. 13 , 495. (doi:10.1504/IJEP.2000.002332) Crossref, Web of Science, Google Scholar - 28.
Danovaro R, Fonda Umani S, Pusceddu A . 2009 Climate change and the potential spreading of marine mucilage and microbial pathogens in the Mediterranean Sea. PLoS One 4 ,e7006 . (doi:10.1371/journal.pone.0007006) Crossref, PubMed, Web of Science, Google Scholar - 29.
Djakovac T, Supić N, Bernardi Aubry F, Degobbis D, Giani M . 2015 Mechanisms of hypoxia frequency changes in the northern Adriatic Sea during the period 1972–2012. J. Mar. Syst. 141 , 179–189. (doi:10.1016/j.jmarsys.2014.08.001) Crossref, Web of Science, Google Scholar - 30.
Fortibuoni T, Giovanardi O, Pranovi F, Raicevich S, Solidoro C, Libralato S . 2017 Analysis of long-term changes in a Mediterranean marine ecosystem based on fishery landings. Front. Mar. Sci. 4 , 33. (doi:10.3389/fmars.2017.00033) Web of Science, Google Scholar - 31.
Haselmair A, Gallmetzer I, Tomašových A, Wieser AM, Übelhör A, Zuschin M . 2021 Basin-wide infaunalisation of benthic soft-bottom communities driven by anthropogenic habitat degradation in the northern Adriatic Sea. Mar. Ecol. Prog. Ser. 671 , 45–65. (doi:10.3354/meps13759) Crossref, Web of Science, Google Scholar - 32.
Fuksi T, Tomašových A, Gallmetzer I, Haselmair A, Zuschin M . 2018 20th century increase in body size of a hypoxia-tolerant bivalve documented by sediment cores from the northern Adriatic Sea (Gulf of Trieste). Mar. Pollut. Bull. 135 , 361–375. (doi:10.1016/j.marpolbul.2018.07.004) Crossref, PubMed, Web of Science, Google Scholar - 33.
Frignani M, Langone L, Ravaioli M, Sorgente D, Alvisi F, Albertazzi S . 2005 Fine-sediment mass balance in the western Adriatic continental shelf over a century time scale. Mar. Geol. 222 , 113–133. (doi:10.1016/j.margeo.2005.06.016) Crossref, Web of Science, Google Scholar - 34.
Covelli S, Fontolan G, Faganeli J, Ogrinc N . 2006 Anthropogenic markers in the holocene stratigraphic sequence of the Gulf of Trieste (northern Adriatic Sea). Mar. Geol. 230 , 29–51. (doi:10.1016/j.margeo.2006.03.013) Crossref, Web of Science, Google Scholar - 35.
Tomašových A, Gallmetzer I, Haselmair A, Kaufman DS, Kralj M, Cassin D, Zonta R, Zuschin M . 2018 Tracing the effects of eutrophication on molluscan communities in sediment cores: outbreaks of an opportunistic species coincide with reduced bioturbation and high frequency of hypoxia in the adriatic sea. Paleobiology 44 , 575–602. (doi:10.1017/pab.2018.22) Crossref, Web of Science, Google Scholar - 36.
Zuschin M, Stachowitsch M . 2009 EPIFAUNA-dominated benthic shelf assemblages: lessons from the modern Adriatic Sea. Palaios 24 , 211–221. (doi:10.2110/palo.2008.p08-062r) Crossref, Web of Science, Google Scholar - 37.
Gallmetzer I, Haselmair A, Stachowitsch M, Zuschin M . 2016 An innovative piston corer for large-volume sediment samples. Limnol. Oceanogr. Methods 14 , 698–717. (doi:10.1002/lom3.10124) Crossref, PubMed, Web of Science, Google Scholar - 38.
Gallmetzer I, Haselmair A, Tomašových A, Stachowitsch M, Zuschin M . 2017 Responses of molluscan communities to centuries of human impact in the northern Adriatic Sea. PLoS One 12 ,e0180820 . (doi:10.1371/journal.pone.0180820) Crossref, PubMed, Web of Science, Google Scholar - 39.
Tomašových A, Gallmetzer I, Haselmair A, Kaufman DS, Vidović J, Zuschin M . 2017 Stratigraphic unmixing reveals repeated hypoxia events over the past 500 yr in the northern Adriatic Sea. Geology 45 , 363–366. (doi:10.1130/G38676.1) Crossref, Web of Science, Google Scholar - 40.
Amorosi A, Bruno L, Campo B, Morelli A, Rossi V, Scarponi D, Hong W, Bohacs KM, Drexler TM . 2017 Global sea-level control on local parasequence architecture from the holocene record of the Po Plain, Italy. Mar. Pet. Geol. 87 , 99–111. (doi:10.1016/j.marpetgeo.2017.01.020) Crossref, Web of Science, Google Scholar - 41.
Bardhan S, Chattopadhyay D, Mondal S, Das SS, Mallick S, Roy A, Chanda P . 2012 Record of intense predatory drilling from Upper Jurassic bivalves of Kutch, India: implications for the history of biotic interaction. Palaeogeogr. Palaeoclimatol. Palaeoecol. 317–318 , 153–161. (doi:10.1016/j.palaeo.2012.01.007) Crossref, Google Scholar - 42.
Anderson M . 2001 A new method for non-parametric multivariate analysis of variance. Austral Ecol. 26 , 32–46. Crossref, Web of Science, Google Scholar - 43.
Legendre P, Anderson MJ . 1999 DISTANCE-based redundancy analysis: testing multispecies responses in multifactorial ecological experiments. Ecol. Monogr. 69 , 1–24. (doi:10.1890/0012-9615(1999)069[0001:DBRATM]2.0.CO;2) Crossref, Web of Science, Google Scholar - 44.
Oksanen J . 2020Vegan: community ecology package . R package version 2.5-7. See https://CRAN.R-project.org/package=vegan. Google Scholar - 45.
Zuschin M . 2024 Data from: Human-driven breakdown of predator-prey interactions in the northern Adriatic Sea. Dryad Digital Repository. (doi:10.5061/dryad.xwdbrv1nr) Google Scholar - 46.
Connolly SR, Roughgarden J . 1999 Theory of marine communities: competition, predation, and recruitment-dependent interaction strength. Ecol. Monogr. 69 , 277–296. (doi:10.1890/0012-9615(1999)069[0277:TOMCCP]2.0.CO;2) Crossref, Web of Science, Google Scholar - 47.
Freestone AL, Osman RW, Ruiz GM, Torchin ME . 2011 Stronger predation in the tropics shapes species richness patterns in marine communities. Ecology 92 , 983–993. (doi:10.1890/09-2379.1) Crossref, PubMed, Web of Science, Google Scholar - 48.
Harper EM . 2003 Assessing the importance of drilling predation over the Palaeozoic and Mesozoic. Palaeogeogr. Palaeoclimatol. Palaeoecol. 201 , 185–198. (doi:10.1016/S0031-0182(03)00624-2) Crossref, Web of Science, Google Scholar - 49.
Huntley JW, Fürsich FT, Alberti M, Hethke M, Liu C . 2014 A complete holocene record of trematode–bivalve infection and implications for the response of parasitism to climate change. Proc. Natl Acad. Sci. 111 , 18150–18155. (doi:10.1073/pnas.1416747111) Crossref, PubMed, Web of Science, Google Scholar - 50.
Talaber I, Francé J, Flander-Putrle V, Mozetič P . 2018 Primary production and community structure of coastal phytoplankton in the Adriatic Sea: insights on taxon-specific productivity. Mar. Ecol. Prog. Ser. 604 , 65–81. (doi:10.3354/meps12721) Crossref, Web of Science, Google Scholar - 51.
Leibold MA . 1989 Resource edibility and the effects of predators and productivity on the outcome of trophic interactions. Am. Nat. 134 , 922–949. (doi:10.1086/285022) Crossref, Web of Science, Google Scholar - 52.
Wootton JT, Power ME . 1993 Productivity, consumers, and the structure of a river food chain. Proc. Natl Acad. Sci. 90 , 1384–1387. (doi:10.1073/pnas.90.4.1384) Crossref, PubMed, Web of Science, Google Scholar - 53.
Calderaro LA, Harnik PG, Rillo MC . 2023 Environmental correlates of molluscan predator–prey body size in the northern Gulf of Mexico. Paleobiology 50 , 70–84. (doi:10.1017/pab.2023.22) Crossref, Google Scholar - 54.
Benović A, Justić D, Bender A . 1987 Enigmatic changes in the hydromedusan fauna of the northern Adriatic Sea. Nature 326 , 597–600. (doi:10.1038/326597a0) Crossref, Web of Science, Google Scholar - 55.
Leighton LR . 2002 Inferring predation intensity in the marine fossil record. Paleobiology 28 , 328–342. (doi:10.1666/0094-8373(2002)028<0328:IPIITM>2.0.CO;2) Crossref, Web of Science, Google Scholar - 56.
Kabat AR . 1990 Predatory ecology of naticid gastropods with a review of shell boring predation. Malacologia 32 , 155–193. Web of Science, Google Scholar - 57.
Lever J, Kessler A, van Overbeeke AP, Thijssen R . 1961 Quantitative beach research: II. the ‘hole effect’. A second mode of sorting of lamellibranch valves on sandy beaches. Netherlands J. Sea Res. 6 , 339–358. (doi:10.1016/0077-7579(61)90008-4) Crossref, Google Scholar - 58.
Vermeij GJ, Dudley EC, Zipser E . 1989 Successful and unsuccessful drilling predation recent pelecypods. Veliger 32 , 266–273. Google Scholar - 59.
Roy K, Miller DJ, Labarbera M . 1994 Taphonomic bias in analyses of drilling predation: effects of gastropod drill holes on bivalve shell strength. Palaios 9 , 413–421. (doi:10.2307/3515059) Crossref, Web of Science, Google Scholar - 60.
Tomašových A, Dominici S, Zuschin M, Merle D . 2014 Onshore-offshore gradient in metacommunity turnover emerges only over macroevolutionary time-scales. Proc. R. Soc. B 281 ,20141533 . (doi:10.1098/rspb.2014.1533) Link, Web of Science, Google Scholar - 61.
Yoder JB . 2010 Ecological opportunity and the origin of adaptive radiations. J. Evol. Biol. 23 , 1581–1596. (doi:10.1111/j.1420-9101.2010.02029.x) Crossref, PubMed, Web of Science, Google Scholar - 62.
Christensen AM . 1970 Feeding biology of the sea-star Astropecten irregularis pennant. Ophelia 8 , 1–134. Crossref, Google Scholar - 63.
Holmes SP, Miller N . 2006 Aspects of the ecology and population genetics of the bivalve Corbula gibba. Mar. Ecol. Prog. Ser. 315 , 129–140. (doi:10.3354/meps315129) Crossref, Web of Science, Google Scholar - 64.
Hrs-Brenko M . 2006 The basket shell, Corbula gibba Olivi, 1792 (bivalve mollusks) as a species resistant to environmental disturbances: a review. Acta. Adriatica. 47 , 49–64. Google Scholar - 65.
Riedel B, Zuschin M, Stachowitsch M . 2012 Tolerance of benthic macrofauna to hypoxia and anoxia in shallow coastal seas: a realistic scenario. Mar. Ecol. Prog. Ser. 458 , 39–52. (doi:10.3354/meps09724) Crossref, Web of Science, Google Scholar - 66.
Crema R, Castelli A, Prevedelli D . 1991 Long term eutrophication effects on macrofaunal communities in northern Adriatic Sea. Mar. Pollut. Bull. 22 , 503–508. (doi:10.1016/0025-326X(91)90405-H) Crossref, Web of Science, Google Scholar - 67.
Kidwell SM . 2013 Time‐averaging and fidelity of modern death assemblages: building a taphonomic foundation for conservation palaeobiology. Palaeontology 56 , 487–522. (doi:10.1111/pala.12042) Crossref, Web of Science, Google Scholar - 68.
Tomašových A, Dominici S, Nawrot R, Zuschin M . 2023 Temporal scales, sampling designs and age distributions in marine conservation palaeobiology. Geol. Soc. Lond. Spec. Publ. 529 , SP529–2022. Crossref, Google Scholar - 69.
Yasuhara M, Yamazaki H . 2005 The impact of 150 years of anthropogenic pollution on the shallow marine ostracode fauna, Osaka Bay, Japan. Mar. Micropaleontol. 55 , 63–74. (doi:10.1016/j.marmicro.2005.02.005) Crossref, Web of Science, Google Scholar - 70.
Kowalewski M, Wittmer JM, Dexter TA, Amorosi A, Scarponi D . 2015 Differential responses of marine communities to natural and anthropogenic changes. Proc. R. Soc. B 282 ,20142990 . (doi:10.1098/rspb.2014.2990) Link, Web of Science, Google Scholar - 71.
Gallmetzer I, Haselmair A, Tomašových A, Mautner AK, Schnedl SM, Cassin D, Zonta R, Zuschin M . 2019 Tracing origin and collapse of holocene benthic baseline communities in the northern Adriatic Sea. Palaios 34 , 121–145. (doi:10.2110/palo.2018.068) Crossref, Web of Science, Google Scholar - 72.
Huntley JW, Scarponi D . 2012 Evolutionary and ecological implications of trematode parasitism of modern and fossil northern Adriatic bivalves. Paleobiology 38 , 40–51. (doi:10.1666/10051.1) Crossref, Web of Science, Google Scholar - 73.
Scarponi D, Azzarone M, Kowalewski M, Huntley JW . 2017 Surges in trematode prevalence linked to centennial-scale flooding events in the Adriatic. Sci. Rep. 7 ,5732 . (doi:10.1038/s41598-017-05979-6) Crossref, PubMed, Web of Science, Google Scholar - 74.
Fitzgerald E, Ryan D, Scarponi D, Huntley JW . 2024 A sea of change: tracing parasitic dynamics through the past millennia in the northern Adriatic, Italy. Geology 52 , 610–614. (doi:10.1130/G52187.1) Crossref, Web of Science, Google Scholar - 75.
Paine RT . 1966 Food web complexity and species diversity. Am. Nat. 100 , 65–75. (doi:10.1086/282400) Crossref, Web of Science, Google Scholar - 76.
Connell JH . 1961 The influence of interspecific competition and other factors on the distribution of the barnacle Chthamalus stellatus. Ecology 42 , 710–723. (doi:10.2307/1933500) Crossref, Web of Science, Google Scholar - 77.
Dayton PK . 1975 Experimental evaluation of ecological dominance in a rocky intertidal algal community. Ecol. Monogr. 45 , 137–159. (doi:10.2307/1942404) Crossref, Web of Science, Google Scholar