The presented data consist of a set of sediment samples obtained in this study and data from previous studies, as shown in Table 1 . The data originate from short sediment cores and surface sediment samples collected along ~140 km of Israel’s inner Mediterranean shelf from off Ashqelon in the south to off Beit Yannay in the north, representing the distal part of the Nile littoral cell at different periods ( Figure 1 ; Table 1 ). The stations were located within the ~30-50 m water depth silty belt, an area most sensitive to past Nile River discharge/flood dynamics (Hecht, 1964; Nir, 1984).

Short sediment cores were sampled in 2007 by the R/V Shikmona off Ashqelon (33 m water depth) and Beit Yannay (40 m water depth), representing the southern and northern zones of the study area. In addition, surface sediment samples were collected in 2021 by the R/V Bat-Galim at three stations (Ashkelon, Ashdod, and Beit Yannay) at ~40 m water depth. Surface sediments and the short cores were collected using a box corer (Ocean Instruments BX 700 AL, Ocean Instruments, Fall City, United States) and Perspex tubes.

The sediment cores were sliced onboard into 1 cm slices, frozen, and lyophilized in the laboratory. Sub-samples were analyzed for granulometry, major elements, CaCO₃, P speciation, total organic carbon (TOC), and δ¹³Corg. Chronology was determined in two parallel cores using lead-210 (details are below).

Duplicates of surface sediment (top 0–1 cm) were sampled for micropaleontological analysis from each box core at each station. Because the abundance of living foraminiferal assemblages is low in the study area (Avnaim-Katav et al., 2020), each duplicate included three repeated samples of the top 0–1 cm interval using a 55 mm diameter Perspex mini corer. Micropaleontological analyses followed Schonfeld et al. (Schonfeld et al., 2012). Sediment samples were stained with rose bengal solution (2 g rose bengal/l 95%- ethanol) for two weeks at the time of sampling. Rose bengal confirms the presence of cytoplasm and is widely used to distinguish between dead and presumed live foraminifera (Avnaim-Katav et al., 2020; Avnaim-Katav et al., 2021). The samples for analyses of dead foraminiferal assemblages, representing pre-AHD deposition, were collected from the bottom layer of two short cores off Beit Yannay and Ashkelon.

Grain size analysis followed Crouvi et al. (2008) and was performed on sediments <2 mm using a Malvern Mastersizer MS-2000 (Malvern Panalytical, Malvern, United Kingdom). Detailed information on the analysis of the major elements, Total Organic Carbon (TOC), δ¹³C of organic matter (δ¹³Corg), and CaCO₃ is presented in Bookman et al. (2021). In general, major elements were analyzed with an inductively coupled plasma-atomic emission spectrophotometer (ICP-OES OPTIMA 3300, Agilent Technologies, Santa Clara, United States) after the fusion of sediment sub-samples with LiBO₂. International standards (e.g., USGS SRS T-207, T-209) were used along with the known samples. The accuracy of the analyses of the major elements was better than 2%. Total Organic Carbon (TOC) content and δ¹³Corg were determined using the Thermo Flash 2000 Elemental Analyzer (manufactured by Thermo Fisher Scientific, Waltham, United States) that interfaced with a Thermo Finnigan Delta V Isotope Ratio Mass Spectrometer (IRMS) (Thermo Fisher Scientific) at the Geological Survey of Israel. The mass spectrometer measures the ¹³C/¹²C ratio in the sample and standard and calculates the δ¹³Corg in ‰. All analyses have been reported relative to VPDB for δ¹³Corg and calibrated against the international standards Urea #2, Acetanilide #1, B2151, and B2153. The CaCO₃ content was determined by gasometry.

P speciation was performed following Eckert et al. (2003) by sequential extraction of the solid phase to metal (Fe)-bound P. P bound to Fe hydroxides was released by leaching with NaOH (~0.1M) and measurement of molybdate-reactive P.

The chronology of the two sediment cores was based on a combination of two radiotracers, ²¹⁰Pb and ¹³⁷Cs, due to their suitable half-lives of 22.2 and 30.0 years. A total of 38 slices of the sediment cores from ASH (19 upper 1 cm slices) and BY (19 upper 1 cm slices) were subjected to gamma spectroscopy. The dried sediments were sealed in cylindrical plastic dishes with a diameter of 3.5 cm using Rn-tight foil. Before measurement, they were left sealed for a minimum of 3 weeks so that the radioactive equilibrium between ²²⁶Ra and ²²²Rn (and their daughters) could be established.

For gamma spectroscopy, a coaxial HPGe detector with 50% relative efficiency (Canberra Industries/Mirion Technologies, Atlanta, United States) was used, housed in a 10 cm Pb shield with Cu and plastic inner linings, and operated under Genie 2000 software (Canberra Industries/Mirion Technologies). Measurement times varied from 77149 to 426623 seconds (long counting times of up to 5 days had to be used for samples with very low activity). Photopeak efficiencies were calculated using LabSOCS^© (Laboratory SOurceless Calibration System), the calibration tool of the Genie 2000 software (Mirion Technologies). The calculation is based on the definitions of sample geometry, composition and density, and sample-detector geometry, and the detector has been factory-characterized for use with LabSOCS. Comparison with measurements of radioactive standard materials revealed deviations of less than 5% between measured and calculated efficiencies in most cases. The activity concentrations of all estimated radioisotopes were recalculated to the date of sampling. For the determination of ²¹⁰Pb_xs, the excess ²¹⁰Pb activity and the supported ²¹⁰Pb activity (determined via the 351.9 keV line of ²¹⁴Pb) were subtracted from the total ²¹⁰Pb signal, measured via the 46.5 keV line. Additionally, the artificial isotope ¹³⁷Cs was analyzed via the 661.6 keV line of the daughter isotope ^137mBa. For several samples, no clear ¹³⁷Cs signal could be detected at the 95% confidence level. The decision threshold varied for individual samples, mainly depending on their mass and counting time, and mostly reached values between 0.4 and 1.0 Bq kg^-1. The chronological models were applied to data expressed in units of activity concentration Bq kg^-1. Comparable results (within standard error) were obtained with data recalculated to Bq cm^-3 using dry densities that do not vary significantly in the profiles. A grain size normalization procedure was applied to the ²¹⁰Pb_xs data (Kirchner and Ehlers, 1998). When ²¹⁰Pb_xs and ¹³⁷Cs enter the marine environment, they are removed from the solution by adsorption on inorganic particles or organic matter in suspension, which is later deposited on the bottom. Most of the activity is linked to the finest fraction with the highest effective surface area. The sediment cores were dated by applying a simple exponential model (constant ²¹⁰Pb_xs flux, constant sedimentation rate) to normalized data from the profiles ( Supplementary Figure 1 ). Variations in sediment particle size were accounted for using the fine fraction for grain size correction (²¹⁰Pb_xs (Bq kg^-1)/fraction <63 μm (%) * 100). When alternatively using Al content (Alvarez-Iglesias et al., 2007) for fine material normalization (²¹⁰Pb_xs (Bq kg^-1)/Al (%) * 100, data not shown), the data look very similar to those after grain size correction. We confirmed this by the appearance of the ¹³⁷Cs values and the geochemical profiles, which support this sedimentation model and the pre-1950s origin of the deeper layers. The accompanying radionuclide data and age model can be found in the PANGAEA data repository (Ransby et al., 2023).

In total, 55 surface and sub-surface sediment samples were washed through a 63 μm sieve and dried at 50°C, and the >125 μm fraction was used for foraminiferal analyses ( Table 1 ). The living assemblage, which included individuals with bright pink and homogeneous cytoplasm, was studied in the entire duplicate sample according to the protocol of Schonfeld et al. (2012). The living assemblage in all of these samples (except for three samples with <50 specimens) had more than 100 individuals, a sufficient number for multivariate statistical methods, and relatively robust environmental interpretation (Hayward et al., 2019). The time-averaged dead assemblage in recently deposited sediments (2021) and pre-AHD dam samples was examined in samples that were split into aliquots containing generally at least 200 dead specimens (one-quarter of the samples contained 118 to 166 shells, and the remaining contained 217 to 590 shells). Specimens were identified to species level and counted to determine the total number of benthic foraminiferal individuals per gram of dry sediment (BF/g). Taxonomic identification was based on Loeblich and Tappan (1987; 1994); Cimerman and Langer (1991); Hottinger et al. (1993); Jones (1994), and the World Foraminifera Database, World Register of Marine Species (2023). The state of preservation of the specimens was evaluated and was generally very good in all samples examined.

To measure the differences in community structure between study periods and study sites in the live and dead assemblages, we performed non-metric multidimensional scaling (NMDS) ordinations after computation of a Bray-Curtis dissimilarity matrix based on foraminifera abundances using the “vegan” package (version 2.5-6) (Oksanen et al., 2019) in R v4.2.2 (R Core Team, 2021). Foraminiferal abundance data were log(x+1) transformed prior to analyses to reduce the influence of the most abundant species. As an additional test, we performed Ward.D hierarchical cluster analysis on the Bray-Curtis dissimilarity matrix using the hclust function in the “vegan” package. The Ward.D algorithm is based on minimizing variances in hierarchically identified assemblages and fits aggregated data, for which the Bray-Curtis measure is generally recommended (Singh et al., 2011). Transformed abundances and clustering were visualized using the “pheatmap” package in R (Kolde and Maintainer, 2018). The significance of the differences between live and dead foraminiferal assemblages derived from the Bray-Curtis matrix was assessed with PERMANOVA (permutational multivariate ANOVA) tests, with period and site as fixed factors, followed by post-hoc pairwise comparisons. We further calculated the relative contribution of each species to the similarities within each assemblage using a SIMPER (similarity percentage) analysis, which examines the percentage contribution of each species to the similarity within and dissimilarity between assemblages (Clarke, 1993).

Following the classification of benthic foraminifera into five categories based on their ecological quality status (EcoQS) (Alve et al., 2016), we have analyzed the temporal changes in the total relative abundance of category 1 (oligotrophic bioindicators) and category 3 (eutrophic bioindicators) using the Kruskal-Wallis rank sum test in R.

Here, we present the geochemical changes in two short sediment cores, one located at the northernmost edge of the Nile littoral cell (BY) and the other ~100 km to the south (ASH), assuming that they best represent the outcomes of the Nile damming. Indeed, the top ~10 cm, accumulated during the last ~40 years, is significantly coarser than the underlying pre-AHD sediments ( Figure 2A ).

Applying the exponential model to the grain size normalized data in cores ASH and BY, we obtained sedimentation rates of 2.3 ± 0.5 mm yr^-1 and 1.8 ± 0.3 mm yr^-1, respectively. The occurrence of ¹³⁷Cs at the depths of 10.5 cm and ~7 cm in cores ASH and BY, respectively, was thus dated to 1962 ( ± 10 yr) and 1967 ( ± 7 yr), respectively, in good agreement with the maximum of the ¹³⁷Cs bomb test fallout in 1963 ( Supplementary Figure 1 ). ¹³⁷Cs is associated with organic matter and the fine clay mineral fraction in sediments, which makes its detection analytically easier in fine-grained, organically rich sediment types. The fact that no ¹³⁷Cs could be detected in the deeper, finer layers, while it was detected with the same experimental setup in the coarser sediments in the shallower layers, strongly supports the pre-1950s origin of the deeper layers. Extrapolation of the ²¹⁰Pb-derived sedimentation rate to depths beyond 10 cm is somewhat uncertain, as the data scatter is very large. On the other hand, for both cores, the ²¹⁰Pb inventory in the upper 10 cm corresponds to about 75% of the total inventory (1305 of 1756 Bq m^-2 for ASH, 3350 of 4033 Bq m^-2 for BC), which also supports the assumption that these sections are not much older than two half-lives of the radioisotope, i.e., approximately 45 years. Although it cannot be excluded from the data that ²¹⁰Pb is present at depths below 20 cm that are not accessible in this data set, the ratio of inventories above and below 10 cm is a strong argument for the age model and its extrapolation to greater depths.

Analysis of the normalized sand and fines (<63µm) inventories in the two cores (kg m^-2; considering the porosity and particle density of 2.65 gr cm^-3) shows a shift at depths that correspond to the mid-1960s, post-AHD, according to the ²¹⁰Pb exponential age model ( Figure 3 ). The sedimentation rates and the slope of the cumulative sands and fines ( Figure 3 ) were used to calculate the mean annual accumulation along the cores. Accordingly, for the BY core, we calculated mean sand and fines accumulation rates of 0.08 and 1.5 kg m^-2 y^-1 before the mid-1960s and 0.66 and 1.3 kg m^-2 y^-1 thereafter, respectively. For the ASH core, mean sand and fines accumulation rates were 0.2 and 1.9 kg m^-2 y^-1 before the mid-1960s and 1.6 and 1.1 kg m^-2 y^-1 thereafter, respectively. This would indicate a respective ~15% and ~40% reduction in fines accumulation rates in BY and ASH cores, respectively, and an ~8-fold increase in sand accumulation rates in both stations. A certain increase in the sand fraction (coarsening) in the top layers would occur solely from the decrease in the flux of Nile-derived fine sediments while maintaining the sand flux as in the previous AHD. Nonetheless, these unexpected calculated accumulation rates suggest that sand transport increased from the Nile Delta northward to the inner shelf band of Israel after the AHD. This estimate holds even if, for some unknown reason, our age model grossly (e.g., by a factor of four) overestimated post-AHD sedimentation rates. A possible cause for such a process is probably associated with the accelerated erosion rates of the sandy Nile Delta (tens of meters per year) after the AHD (3-5 times the rates before the AHD) and the northeastward transport and export of these sands along the coast (Frihy, 1988; Smith and Abdel-Kader, 1988; Frihy and Komar, 1991; Inman et al., 1993; Frihy and Lotfy, 1997). It has been proposed that the blockage of the blockage of the Nile sand caused by the AHD forced currents to remove sand from the delta coast, compensating for the reduction and preventing sand shortages further north in the Nile littoral cell (Zviely et al., 2007). Additionally, man-made structures such as breakwaters in ports in Egypt, Gaza, and Israel may have diverted transported sand further offshore. Other factors that may somewhat affect sediment transport and grain size distribution on the shelf may include extreme flooding of the large ephemeral El-Arish River (Klein, 2000), fluctuations in erosion rates of the coastal escarpment (Mushkin et al., 2016), and various physical processes (Van Wellen et al., 2000; Shaeri et al., 2020).

The sand content exceeds ~50% in the post-AHD sediments, compared to <20% in the pre-AHD sediments ( Figure 2A ). The grain size mode increased gradually from ~5 µm in the pre-AHD sediments to ≥85µm in the post-AHD sediments. The silt fraction, which was ~75% in the pre-AHD sediments, decreased to ~45% in the post-AHD sediments, and the clay content decreased by ~50%, from ~20% to <10%. The latter reduction in fine sediments post-AHD is reflected in the changes in the chemical composition of the sediments, which show significantly higher Al and Fe concentrations and lower Si/Al ratios in the pre-AHD sediments ( Figures 2D, F ). The significantly lower K/Al, Ti/Al, and Si/Al ratios in the pre-AHD shelf sediments ( Figures 2D, E ) coincide with these ratios in the Nile sediments (Krom et al., 1999), emphasizing their strong reduction in the post-AHD sediments. Sediment coarsening is also accompanied by a distinct but moderate increase in CaCO₃, from less than 10% to 15-20%, mainly attributed to the increase in BF abundance ( Figure 2C ). The grain size coarsening is more abrupt and rapid in the southern shelf (ASH core), closer to the Nile, than in the northern distal part (BY core; Figure 2 ). In the latter, the coarsening starts earlier, is more gradual, and increases in two steps, probably in response to an earlier stage of Nile damming (since the Aswan Low Dam in 1898 -1902) that was barely recorded in the southern shelf.

TOC contents of 0.8-1wt.% with δ¹³Corg of approximately -19.5‰ in pre-AHD sediments decreased to less than 0.5wt.% TOC and approximately -22-23‰ δ¹³Corg in the overlying younger sediments ( Figure 2B ). Metal (mainly iron oxides)-bound phosphorus concentrations are lower in the post-AHD sediments, probably reflecting reduced scavenged phosphate in iron oxides as compared to the enhanced discharge of nutrients and iron in the pre-AHD period, especially during the seasonal floods ( Figure 2F ).

The bulk geochemical data from the two sediment cores are presented in Supplementary Table 1 . Processing the geochemical data for principal component analysis (PCA) reveals two principal components that account for approximately 94% of the total variance in both sediment cores. Supplementary Figure 2 shows that the variance in elemental composition corresponds to the relative contribution of two major sources: terrestrial and marine-biogenic. The marine component contributes mainly to Ca or CaCO₃ as a major constituent of micro- and macro-faunal shells and to Si as a major constituent of quartz grains, both of which affect the proportion of sand grain size. Al, Fe, Mg, K, and Ti mostly correspond to terrestrial aluminosilicate minerals. Metal (mainly iron oxides)-bound phosphorus (P-metal) concentrations are linked to terrestrial sources of both dissolved phosphate and iron. The variance of TOC and its isotopic composition seem predominantly affected by the change in the terrestrial discharge of dissolved nutrients that trigger the marine-born TOC.

The relationships between δ¹³Corg and TOC ( Figure 4 ) show a mixing curve between depleted TOC with relatively more negative δ¹³Corg in the post-AHD period and a more enriched TOC with a marine/algae δ¹³Corg signature in the pre-AHD sediment. The ~3 ‰ decrease (-19.2 to -22.6 ‰; one measurement of -24.6 ‰; Figure 4 ) in δ¹³Corg may reflect a shift in the relative contribution of autochthonous versus allochthonous organic carbon. Marine planktonic δ¹³C in the eastern Mediterranean and other marine areas at similar latitudes ranges between -18 and -22 ‰ (Goericke and Fry, 1994; Meyers, 1994; Harmelin-Vivien et al., 2008; Li et al., 2016), and terrestrial δ¹³C ranges from -25 to -28 ‰ (Hedges et al., 1997; Li et al., 2016). The ~50% decrease in TOC may indicate a drastic reduction in dissolved nutrient supply post-AHD. This observation suggests that prior to the AHD, the Nile system supplied significant amounts of dissolved nutrients to its distal cell, which enhanced algal biomass and the sedimentation of autochthonous organic carbon, but only an insignificant amount of allochthonous organic carbon. While the pre-AHD transported significant amounts of silty-clay sediments, it is speculated that the Nile-born detrital organic matter was decomposed or retained in the Nile delta, while the dissolved nutrients, representing a much larger stock, were transported northward along the eastern Levantine coast, fueling the primary producers, especially during the seasonal floods. The floods peaked between August and October/November (Sharaf El Din, 1977) and triggered a large phytoplankton bloom in the nearshore waters off the Egyptian and Israeli coasts (Oren and Komarovsky, 1961; Dowidar, 1984). At present, the southeastern Mediterranean Levantine Basin and coastal waters off Israel display oligotrophic conditions (Herut et al., 2000; Kress et al., 2019). Nonetheless, partial, but probably not dominant, changes in the TOC content can be attributed to the significant increase in grain size after the AHD (Hedges and Oades, 1997; Magill et al., 2018).

The changes in δ¹³Corg composition in shallow continental shelf sediments are a first-order estimate due to mixing between marine/autochthonous and terrestrial/allochthonous organic carbon (Shultz and Calder, 1976). The isotopic mixing equation then represents a conservative mixing of these two end members (Shultz and Calder, 1976):

where Fm and Ft are the fractions of marine and terrestrial organic carbon (Fm + Ft = 1), and δ¹³Corg-m and δ¹³Corg-t are the isotopic compositions of the terrestrial and marine source end members.

Rearranging Equation 1,

Based on the above literature data, we calculated Fm (eq. 2) using a terrestrial δ¹³C-t of -27 ‰ and δ¹³Corg-m -19 ‰ as end members. Fm ranges from 0.4 to ~0.9 before the AHD, corresponding to higher TOC concentrations. We assumed that the Si/Al or Al ratios between the sample and the Nile sediment end member (Si/Al ratio of 2.483 wt./wt. and Al concentration of ~11 wt.%; Krom et al., 1999) could serve as a proxy for the intensity of dissolved nutrient supply, as they coincide with the transport of Nile-derived fine sediments. These ratios may represent the fraction of marine-born organic carbon attributable to the supply of Nile-derived nutrients, as they show a similar fractional range and a significant linear correlation with Fm calculated from carbon isotopic composition ( Figure 4 ).

The decline in these variables reflects a major change in the regional nutrient budget and a large decrease in primary production that coincides with the damming of the Nile (Halim et al., 1995). The trend of increasing oligotrophy on the inner shelf of Israel differs from reports of a large anthropogenic contribution of nutrients supporting increasing fisheries in the Mediterranean coastal waters off Egypt (Nixon, 2003). This increase in fertility seems to be on a local scale, mainly restricted to the delta area, unlike the pre-AHD summer floods, which were of significant and large magnitude and affected the Israeli shelf annually. The decrease in TOC and δ¹³Corg, indicators of nutrient supply and primary production, predates the sharp sediment coarsening. This may indicate that earlier phases of Nile damming had already contributed to the increasing oligotrophy of the EM, which accelerated after the operation of the AHD and continues at present.

The current study traces the temporal variations in the BF assemblage composition linked to changes in the sedimentary regime in the distal part of the Nile littoral cell following the damming of the Nile. We presented the BF assemblages in sedimentary layers representing pre (~1910) and post (1997; 2011; 2021) AHD conditions. The integrated database of surface and subsurface sediment samples ( Table 1 ) contains a total of 170 taxa, of which 120 were identified at the species level after taxonomic refinement ( Supplementary Table 2 ). The comparison of species accumulation curves of 10 surface sediment samples collected in 2021 between BY and ASH stations (three versus five sites) indicates that the three chosen sites represent 92% of the total species richness (56/61 species) ( Supplementary Table 3 ; Supplementary Figure 3 ).

Cluster analysis and nonmetric multidimensional scaling (NMDS) ordination revealed a clear separation between dead (cluster A) and live (cluster B) BF assemblages and time periods, as detailed below ( Figure 5 ; Supplementary Figure 4 ; Supplementary Table 4 , Permanova analysis). We focused mainly on the dead BF assemblage, which provides the long-term record since before the AHD and integrates time-averaged information. The dead BF (Cluster A, Supplementary Figure 4 ) are divided into two major sub-clusters, A1 and A2, which are further subdivided by periods. A1 includes the separation between pre-AHD and 1997 samples, and A2 shows the division between 1997, 2011, and 2021 ( Figure 5 ). Pairwise analysis between periods showed a significant difference (p<0.02) between pre-AHD, 2011 and 2021, and between 1997, 2011, and 2021 ( Supplementary Table 5 ). The opportunistic species Ammonia tepida, Porosononion subgranosus, and Cribroelphidium poeyanum are among the most significant species in the pre-AHD samples (in sub-cluster A1), contributing 35% of the similarity between these samples (SIMPER, Supplementary Table 6 ). Their relative abundances ranged between 21-42%, 4-22%, and 8-18% of the assemblage composition, respectively. While these species may inhibit sediments poor in TOC, their highest abundances are favored by organic carbon enrichment (e.g., Bouchet et al., 2018; Jorissen et al., 2018). These species are also part of the 1997 samples; however, their contribution is relatively much lower, while other species such as Ammonia parkinsoniana, Asterigerinata mamilla, and Textularia agglutinans show increased relative abundances, ranging between 3-7%, 4-14% and 2-3%, respectively. These latter species are known to be sensitive to organic carbon enrichment and occur mainly in natural oligotrophic ecosystems (Alve et al., 2016; Jorissen et al., 2018).

We used the classification of ecological groups following Bouchet et al. (2021) as a proxy for changes in sedimentary total organic carbon levels. Based on the ecological quality status (EcoQS) BF groups, the two most relevant categories for changes in the Nile damming sediments are category 1 (sensitive species - oligotrophic bioindicators) and category 3 (tolerant species to organic carbon enrichment – elevated trophic state bioindicators). EcoQS category 1 showed an increasing trend over time (p<0.02; Kruskal-Wallis rank sum test), while category 3 showed significantly lower levels during 1997-2021 than in the pre-AHD ( Figure 6 ; Supplementary Table 7 ). These trends are consistent with the geochemical observations showing lower TOC and fine sediment post-AHD.

EcoQS category 2 includes “indifferent species” to organic carbon enrichment, which are usually observed at relatively low abundances (Bouchet et al., 2021). It is thus expected that this group will be less responsive to the inferred damming effects and the decrease in TOC, as indeed observed along the periods ( Figure 6 ). They occur in low abundance, with no changes between the different periods.

Additionally, based on pairwise analyses, the similarity between sites within periods showed that in pre-AHD and 1997, the BY station was the most dissimilar to the ASH and ASD stations (p<0.04, Permanova, Supplementary Table 8 ), while in 2011-2021, a greater dissimilarity developed between the two southern stations (ASH and AHD). This spatial dissimilarity trend suggests a further ongoing process from the northernmost distal part of the Nile cell southward.

Although live BF populations represent a snapshot of seasonal distribution that may record taphonomic processes compared to the subsurface assemblages, they are also responsive to anthropogenic impacts (e.g., Kidwell, 2007; Goineau et al., 2015). Indications for the ongoing changes in the BF composition associated with the damming of the Nile are also recorded in the live BF communities collected in 2011 and 2021. Eggerelloides scaber, Paratrochammina madeira, and the sensitive species Bolivina striatula, Cycloforina quinquecarinata, and Textularia agglutinans are among the living BF that contribute most (60% contribution) to the similarity between the samples of cluster B, based on SIMPER ( Figure 5 ; Supplementary Table 6 ). EcoQS category 1 showed a slight decrease (of about 10%) in their abundance compared to their appearance in the dead assemblages of 2020 ( Figure 6 ), while category 3 of opportunistic species continued to decrease in their density, representing 4% of the assemblage. The ongoing changes in sediment properties demonstrate the continuous impact of the Nile damming and the subsequent enhanced benthic oligotrophication process. Changes in the composition of the live BF assemblage have been observed over the last two decades on the southeastern Mediterranean shelf, where the Nilotic biotope has retreated southward (Avnaim-Katav et al., 2020; Avnaim-Katav et al., 2021). We speculate that the ongoing environmental changes in this unique ecosystem will continue to affect the distribution of benthic communities along the shelves of the Nile domain and thus require further study.