The disposal site in question has long received sewage sludge from the Dan Region Wastewater project, known as SHAFDAN, the municipal treatment system serving the densely populated Tel Aviv metropolitan area. Sludge disposal at sea is a practice with a long and contested history worldwide, and the Israeli site has been monitored for decades precisely because the eastern Mediterranean shelf is an unusually demanding environment for such activity. The Levantine basin is warm, oligotrophic, and increasingly stressed by seawater warming, so any additional organic load delivered to the seabed lands on an ecosystem already near its tolerance limits. Earlier work by the same research group and by colleagues at the Israel Oceanographic and Limnological Research had documented the chemical and biological footprint of the sludge, but the new study spans the full transition from active disposal to its complete cessation in March 2017, giving scientists an unusually clean before-and-after comparison.

The technical heart of the study is a twelve-year benthic time series built from box-corer samples collected aboard the research vessels Shikmona and Mediterranean Explorer. Each spring and fall, the team quantified macrofaunal abundance, the composition of taxa present, taxa richness, and the total organic carbon content of the sediment, commonly abbreviated TOC. TOC serves as a chemical proxy for organic enrichment: the more organic material, whether from sludge, phytoplankton debris, or other sources, has settled into the sediment, the higher the reading. Macrofauna, the worms, mollusks, crustaceans and other small animals living within the sediment, act as the biological counterpart to that chemical signal. Because different species have very different tolerances for organic loading and low oxygen, the identity and abundance of the animals present integrate the environmental history of the site over months to years. Reading the two signals together, chemistry and biology, is the standard logic of benthic impact monitoring, and it is what allowed the researchers to disentangle the effects of sludge from those of natural seasonality and regional warming.

During the active disposal period, before March 2017, the picture was stark. The community was overwhelmingly dominated by one or more opportunistic species of capitellid polychaetes, a family of segmented worms famous in marine ecology as classic indicators of disturbed, organically enriched sediments. Capitellids thrive where other animals cannot: they tolerate the hypoxic, sulfide-rich conditions that develop when excess organic matter fuels bacterial respiration in the seabed. Their monopolization of the disposal area is textbook ecological succession in reverse, a community pushed back to the earliest, most degraded stage described in the Pearson-Rosenberg model of organic enrichment. In practical terms, the sea floor beneath the disposal zone had been converted from a diverse shelf community into a monoculture of pollution-tolerant worms, a biological signature of chronic organic loading that persisted as long as the sludge kept arriving.

Perhaps the most striking finding from the disposal years is that the impact was not static but rhythmic. The researchers identified an annual cycle driven by the interplay between waste delivery and winter hydrodynamics. Sludge accumulated on the seabed through the summer months, and then winter storms dispersed it, scouring and redistributing the organic layer across the shelf. This accumulation-and-dispersion cycle produced a matching annual cycle in macrofaunal abundance, with the fall samples, taken after the summer accumulation phase, showing the strongest suppression. By fall, the accumulated sludge overrode every other environmental factor the team could evaluate, prominently attenuating the abundance of the macrofaunal community. TOC levels, meanwhile, remained stable and relatively high year-round during the disposal period, indicating that the sediment’s organic burden never fully cleared between dumping seasons. The system, in other words, was locked into a perpetual seasonal pulse of enrichment and partial recovery, never allowed to escape the disposal footprint.

Then came a twist that makes the study more than a routine impact assessment. Between 2012 and 2017, well before the disposal actually stopped, the community began to change in complex, gradual ways. The dominance of the opportunistic capitellids attenuated, a wider variety of other taxa entered the samples, and taxa richness shifted accordingly. Crucially, these changes were not timed with the cessation of sludge disposal, which did not occur until March 2017. The researchers initially hypothesized that the driver was the warming of the eastern Mediterranean seawater, a trend well documented for the Levantine basin and increasingly implicated in ecological upheavals across the region, from multi-species collapses at the warm edge of the sea to shifts in sponge and fish populations. Whatever the precise mechanism, the episode is a cautionary tale for monitoring programs everywhere: community change in a polluted area does not necessarily mean the pollution has changed, and attributing temporal shifts to a single management action requires careful timing analysis.

The cessation itself, when it finally came, produced a response that was both delayed and decisive. In 2018, one year after disposal ended, the disposal area settled into a new steady state that was unmistakably different from the disposal-era baseline. Macrofaunal abundance settled at low levels, but taxa richness was relatively high, the taxonomic composition stabilized in a new configuration, and TOC levels dropped to background levels comparable to unaffected shelf sediments. The chemical recovery was remarkably fast by marine standards; organic enrichment that had persisted year-round during the disposal period was effectively flushed from the system within a year once the source was cut. The biological community, however, did not simply revert to some pre-disposal original state. Instead, it reorganized into a novel assemblage, a reminder that recovered ecosystems are not always restored ecosystems, and that legacies of decades of disturbance can shape community structure long after the pressure is removed.

The study’s spatial dimension adds further texture to the story. By sampling multiple stations across the disposal area and its surroundings, the team could distinguish the core impact zone, where sludge effects were strongest, from peripheral areas where the annual accumulation-dispersion cycle and the post-cessation recovery played out differently. This spatial gradient is what makes the dataset valuable beyond the Israeli coast. Similar disposal sites and wastewater outfalls around the world, from the New York Bight to Australian waters to the Gulf of Lions, have produced comparable patterns of opportunistic dominance and slow recovery, and the Israeli series confirms that the same ecological principles operate in the ultra-oligotrophic, rapidly warming eastern Mediterranean. It also demonstrates the value of sustained, twice-yearly sampling with consistent methods; a shorter or less frequent program could easily have missed either the pre-cessation community shift or the one-year chemical rebound, and drawn the wrong conclusion about cause and effect.

For marine managers, the takeaways are concrete. First, cessation works: cutting off the sludge source allowed the sediment chemistry to reset within roughly a year and released the macrofaunal community from its capitellid monoculture. Second, recovery is not instantaneous or symmetrical; the biological reorganization took years and produced a new steady state rather than a return to the past. Third, disentangling multiple stressors requires long time series, because climate-driven change can mimic, mask, or interact with local pollution signals in ways that short-term monitoring cannot resolve. As coastal nations worldwide phase out sea disposal of sewage sludge under international conventions, the Israeli shelf offers one of the clearest documented examples of what the other side of that decision looks like: a sea floor that, given the chance, reassembles itself into something more diverse, more stable, and far more representative of what a healthy Mediterranean shelf should be.

Subject of Research: Long-term effects of sewage-sludge disposal and its cessation on benthic macrofaunal communities of the Israeli Mediterranean shelf

Article Title: Effects of sewage-sludge disposal and its cessation on benthic macrofaunal community on the Israeli Mediterranean shelf

Article References: Lubinevsky, H., & Tom, M. (2026). Effects of sewage-sludge disposal and its cessation on benthic macrofaunal community on the Israeli Mediterranean shelf. Environmental Monitoring and Assessment, 198(10), Article 1122. https://doi.org/10.1007/s10661-026-15967-x

Image Credits: AI Generated

DOI: 10.1007/s10661-026-15967-x

Keywords: sewage sludge, benthic macrofauna, Mediterranean Sea, Israeli shelf, polychaetes, capitellids, total organic carbon, sediment monitoring, marine pollution, ecological recovery, Levantine basin, environmental monitoring

The research team, led by Amanda Martins of the Federal University of Bahia together with Marcos Krull and Francisco Barros, combined two rarely united sources of evidence: a long-term benthic monitoring dataset stretching back to 2004, and controlled laboratory experiments that measured exactly how much organic matter individual worms eat and how much sediment they displace. The monitoring data came from ten sampling stations along the estuary’s salinity gradient, surveyed on six separate occasions between 2004 and 2022. At each station, researchers collected cores of sediment, sieved the samples through a 0.5 millimeter mesh, and identified the invertebrates to the family level. Polychaetes made up more than half of all benthic macroinvertebrates in the system, making them the natural focal group for the study.

To translate counts of worms into estimates of ecosystem-wide function, the team turned to Bayesian generalized linear mixed models, a statistical framework that allows researchers to predict abundance across unsampled areas while quantifying uncertainty. The models incorporated environmental predictors including salinity and sediment grain size fractions—ranging from coarse sand upstream to fine sand and mud downstream—drawn from high-resolution environmental layers covering the modeled estuarine portion of the system. The researchers tested both Poisson and negative binomial distributions, with and without zero-inflation parameters, and compared generalized additive and generalized linear mixed models. Model selection relied on leave-one-out cross-validation using Pareto-smoothed importance sampling, and the winning models were fitted with Markov Chain Monte Carlo methods using a Hamiltonian sampler, with four independent chains and weakly informative priors to regularize the estimates.

The result of this statistical machinery was a population estimate of roughly 2,062.5 million individual polychaetes across the estuarine system, with a 95 percent credible interval spanning about 1.87 to 2.63 billion individuals. Seven families dominated the analysis: Orbiniidae, Spionidae, Eunicidae, Pilargidae, Lumbrineridae, Onuphidae, and Goniadidae. Orbiniidae alone accounted for an estimated 1,045.9 million individuals, followed by Onuphidae with 671.6 million. Densities were highest in the marine-influenced regions of the estuary, particularly for these two dominant families, setting the stage for a strongly spatially patterned delivery of ecological functions.

The functional measurements themselves came from a previous experimental study in which worms collected from the Jaguaripe estuary were kept under controlled laboratory conditions. Organic matter decomposition was quantified as the consumption rate of standardized pieces of shrimp tissue over 96 hours, corrected for natural decomposition using controls. Bioturbation—the reworking of sediment—was measured as the volume of sediment displaced by each species, reconstructed from computed tomography scans of the burrow systems the worms produced. This CT-based approach allowed the researchers to visualize and measure the three-dimensional architecture of worm burrows without disturbing them, providing an unusually precise metric of a notoriously difficult-to-quantify process.

When the per-individual rates were multiplied by the model-predicted abundances across the estuary, the numbers were striking. Polychaetes consumed an estimated 232.56 tons of organic material over 96 hours, with Orbiniidae responsible for 123.16 tons—about 53 percent of the total—and Onuphidae contributing another 89.81 tons, or nearly 39 percent. Bioturbation totaled an estimated 6,098.3 cubic meters of sediment over the same period, again dominated by Orbiniidae at 4,495.1 cubic meters and Onuphidae at 1,070.7 cubic meters, together accounting for more than 91 percent of the total. Smaller families such as Eunicidae, Spionidae, and Lumbrineridae contributed comparatively little to bioturbation, while Eunicidae and Pilargidae were minor players in organic matter consumption as well.

The spatial patterns were equally revealing. Both functions peaked in regions of elevated salinity, in the euhaline and polyhaline zones near the estuary’s mouth. Orbiniidae and Onuphidae delivered their largest contributions at salinities between 25 and 40, whereas Spionidae showed a strikingly different profile, contributing up to 80 percent of decomposition and bioturbation in the lower-salinity upstream reaches, between roughly 5 and 20 on the salinity scale. Sediment grain size also mattered: the dominant families responded to different granulometric fractions, with Orbiniidae and Onuphidae key in coarse sand and mud regions respectively, and Spionidae playing its characteristic upstream role in very fine sand. No single sediment variable was consistently selected across all taxa, underscoring that each family responds to its own combination of environmental conditions.

Correlation analyses of the functional contributions revealed another important pattern: functional redundancy. In the marine lower estuary, the contributions of Orbiniidae, Onuphidae, and Lumbrineridae rose and fell together, with correlation coefficients for bioturbation reaching 0.93 between Orbiniidae and Onuphidae and 0.92 between Lumbrineridae and Onuphidae. This overlap means that if one family were lost, others performing similar roles could partially compensate, buffering the ecosystem against disturbance. In the oligohaline upstream zones, by contrast, fewer families contributed to the functions, and redundancy was correspondingly lower—making those reaches potentially more vulnerable to species loss. Spionidae, with low correlations to all other taxa, stands out as functionally irreplaceable in the upstream areas where it dominates.

The study’s authors are candid about its limitations. Only seven families and two functions were analyzed, other taxa in the estuary may contribute substantially, and variables such as food resources beyond the measured organic matter could not be included. The decomposition measurements represent early-stage, macrofauna-mediated breakdown driven by direct consumption at the sediment surface, not the full complexity of benthic metabolism. Body size, which likely influences per-individual function, could not be incorporated because weight data were unavailable for several taxa. Still, the researchers argue, the demonstration that laboratory experiments can be combined with long-term monitoring and Bayesian modeling to produce ecosystem-level estimates is itself a methodological advance, offering a template for other systems where functional data remain scarce—particularly in the tropics, where most species-specific functional measurements to date have been conducted in temperate terrestrial environments.

The broader implications extend to conservation and climate adaptation. Estuaries are hotspots of organic matter input and nutrient cycling, disproportionately productive relative to their area, yet they face mounting pressures from eutrophication, metal pollution, shrimp farming, and sea-level rise. Because polychaete bioturbation aerates deeper sediment layers, reduces compaction, and accelerates organic matter processing, the loss of these worms could tip sediments toward anoxia, with cascading effects on entire benthic communities. The functional maps generated by the study highlight regions of high activity that could guide conservation priorities, and the identification of low-redundancy upstream zones flags areas where monitoring should be intensified. As climate change reshapes salinity gradients and species distributions, the authors suggest that modeling the distribution and functioning of organisms under future scenarios will be essential for predicting trends—and for protecting the billions of small engineers whose invisible labor keeps coastal ecosystems running.

Subject of Research: Ecosystem-level estimation of polychaete-mediated organic matter decomposition and bioturbation in a tropical estuary

Article Title: Scaling up polychaete contributions to estuarine ecosystem functions

Article References: Martins, A., Krull, M., & Barros, F. (2026). Scaling up polychaete contributions to estuarine ecosystem functions. Discover Ecology, 2(1), Article 6. https://doi.org/10.1007/s44396-026-00023-2

Image Credits: AI Generated

DOI: 10.1007/s44396-026-00023-2

Keywords: polychaetes, estuarine ecology, bioturbation, organic matter decomposition, ecosystem functioning, Bayesian modeling, functional redundancy, salinity gradient, macrofauna, tropical estuary, sediment reworking, conservation