Abstract
The coastal intertidal ecosystem of the Bijagós Archipelago, Guinea-Bissau, one of the largest and most important in West Africa, sustains a considerable proportion of the migratory shorebird populations of the East Atlantic Flyway and operates as a nursery area for benthic fish in the region. The macrozoobenthos in these mudflats constitute the main food source for both groups so that spatial and temporal variation in their abundance and community composition is likely to influence the abundance and distribution of fish and birds. In this study we described the spatial and temporal dynamics in the density, biomass, and community composition of macrozoobenthos across six intertidal flats in three islands of the Bijagós Archipelago. Overall, the Bijagós Archipelago was characterised by a highly species-rich macrozoobenthic community, with ca. 88 taxa identified across all sites, reaching a mean density of 1871 ± 58.3 ind.m-2 (mean ± SE) and mean biomass of 5.65 ± 0.41 g of AFDM.m-2 (ash-free dry mass per m2), values much lower than what was described for nearby intertidal areas, namely the Band d’Arguin, Mauritania. Density and biomass of major macrozoobenthos classes (Bivalvia, Polychaeta, Malacostraca and Gastropoda) differed across sites and months, displaying an overall increase in density towards the final months of the dry season (March and April). Similarly, community composition also differed significantly between sites and throughout the season. The site with most distinct community composition (Adonga) supported low diversity and high abundance of a few bivalve species, whilst other two sites that hosted the most diverse communities, were also the most similar between them (Anrumai and Abu). These spatial and temporal patterns constitute an important baseline to improve knowledge of this intertidal ecosystem and will contribute towards a better understanding of the spatial and temporal distribution patterns of their consumers.
Citation: Coelho AP, Henriques M, Rocha AD, Paulino J, Schaars LK, Ramos C, et al. (2022) Spatial and seasonal variation in macrozoobenthic density, biomass and community composition in a major tropical intertidal area, the Bijagós Archipelago, West-Africa. PLoS ONE 17(11): e0277861. https://doi.org/10.1371/journal.pone.0277861
Editor: Marco Fusi, KAUST University, SAUDI ARABIA
Received: July 5, 2022; Accepted: November 4, 2022; Published: November 28, 2022
Copyright: © 2022 Coelho et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting information files.
Funding: This research was conducted within the framework of the project “Waders of the Bijagós” funded by MAVA Foundation. Fundação para a Ciência e Tecnologia (FCT Portugal) provided additional financial support to CESAM (UIDP/50017/2020+UIDB/50017/2020+LA/P/0094/2020), and through the project “MigraWebs: Migrants as a seasonal ecological force shaping communities and ecosystem functions in temperate and tropical coastal wetlands” (PTDC/BIA/ECO/28205/2017). MH was supported by FCT through doctoral grants (SFRH/BD/131148/2017 and COVID/BD/151924/2021). TC was supported by FCT trough contract 2021.00573.CEECIND. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
1. Introduction
Intertidal systems tend to be highly productive and provide numerous ecosystem services, from nutrient recycling to shoreline stabilisation and food production [1, 2]. These ecosystems are among the most widespread coastal habitats globally [3], playing an important role in the support, protection and resilience of coastal areas and their inhabitants [4–6]. Intertidal flats are threatened by global change, with sea level rise, global warming and other climatic, chemical and physical changes among the major threats, aggravated by direct anthropogenic action, like coastal development [7, 8]. These have an impact over the species assemblages in intertidal ecosystems, affecting their distribution, food and habitat availability, community composition and ultimately the entire ecosystem functioning [7, 9, 10]. It is therefore urgent to describe the current state of these systems across biological organization levels, in order to better understand future changes, particularly in understudied areas of the globe.
Intertidal flats host macrozoobenthos communities, a diverse group of species adapted to the variable environmental conditions associated with tidal regimes [11]. These communities play a fundamental role on intertidal ecosystem functioning and structure being involved in trophic and non-trophic interactions and energy transfer (e.g. [12–16], with a particularly relevant role as food source for other species in the intertidal system, such as fish and shorebirds [17–19]. The distribution and abundance of macrozoobenthos is linked to environmental conditions, such as temperature, hydrodynamic regime, sediment type and composition, primary food sources and nutrients [20–23], and to ecological characteristics, such as the presence of ecological engineers or predators [24–27]. For example, the presence of the West African fiddler crab (Afruca tangeri), an ecosystem engineer in intertidal sediments, has been shown to strongly limit the abundance and diversity of macrozoobenthos communities in the Bijagós Archipelago, Guinea-Bissau [28]. Similarly, exclusion experiments in a tropical intertidal flat have described the cascading effects that the removal of a predator species can have on different levels of the ecosystem, from changes in the abundance of and growth rates of prey species eventually leading to population declines, to altering the chemical processes in the mudflats [14]. The distribution and community composition of macrozoobenthos can in turn influence the distribution of its predators [28–31], which in turn can have interacting effects on the rest of the food chain [14, 32].
Compared with temperate systems, tropical intertidal flats are often characterised by a relatively higher species richness but lower density and biomass of macrozoobenthos [11, 20, 33]. Community composition in tropical regions varies greatly between sites, being highly influenced by the action of ecosystem engineers and generally having only a few abundant species, with the remaining occurring at low densities [20, 24]. Even though temperature fluctuation tends to be lower than in temperate sites, heavy rainfall during the wet season may produce substantial changes in environmental parameters, leading to marked seasonal variation in macroinvertebrate distribution and abundance [34]. However, there is currently limited information for many tropical intertidal ecosystems, particularly in West African areas, and establishing baseline knowledge on macrozoobenthos communities is important to improve our global understanding of intertidal ecosystems.
The Bijagós Archipelago, located off the West African coast, is an inverse deltaic estuary system (i.e. an estuary where salinity increases upstream), with more than 80 islands, and is surrounded by vast intertidal flats that are mostly fringed by dense mangrove forests [35, 36]. The Archipelago is one of the most important non-breeding areas for migratory shorebirds in the East Atlantic Flyway, being used by up to 10% of all the shorebirds migrating along this flyway [37–39]. It also provides refuge for important fish populations, being considered one of the major nursing areas in West Africa, and home to several species of international conservation concern [36, 40–42]. Both shorebirds and benthic fish rely on the productivity of the intertidal flats, feeding on its macrozoobenthic communities to sustain their populations. Despite the national and international recognition of the archipelago, and even though recent studies have been shedding light on the structure and functioning of its intertidal ecosystem (e.g. [28, 32, 43–45]), the spatio-temporal patterns and composition of the macrozoobenthic community remain poorly understood, with significant knowledge gaps on many groups of species and on the spatial and temporal representativity of existing studies [45, 46].
Studies on the macrozoobenthos of the Bijagós Archipelago are scarce, but recent papers reported overall low densities and biomass, with considerable variation in species composition between sites [46] and habitat partitions [28]. Nevertheless, these studies included a limited number of sites (three and one, respectively), hampering wider conclusions to the whole Bijagós Archipelago. Furthermore, little is currently known regarding how these macrozoobenthic communities vary within and between seasons, particularly during periods of increased predation pressure caused by migratory shorebirds. These shorebirds congregate in the Bijagós Archipelago from October to April [17, 37, 47], and were shown to have the potential to deplete their food stocks during their temporary stay at another West-African intertidal area, the Banc d’Arguin, in Mauritania [29]. In this study we aim to describe the overall spatial and temporal variation in the density, biomass, and community composition of macrozoobenthos in the intertidal flats of the Bijagós Archipelago, from the end of the wet season to the end of the dry season (October to April). Specifically, we investigate if the density and biomass of major macrozoobenthic classes, and the overall community composition are different between sites across the months. Furthermore, we describe which species are driving the differences observed between sites and periods (by aggregating months in three climatic and ecologically distinct periods, see methods), and also assess the differences in community diversity between sites within each period. In addition, we contribute with a set of undocumented length-biomass regressions for important macrozoobenthic taxa in the study area.
2. Methods
2.1. Study area
The Bijagós Archipelago is a group of 88 islands and islets off the coast of Guinea-Bissau, in West Africa, encompassing ca. 450 km2 of intertidal flats [48, 49]. Its topographical characteristics promotes large tidal amplitudes reaching up to 4.5 m in spring tides [50]. The climate is tropical and bi-seasonal, with a wet season from May to October, and a dry season from November to April. The islands are mostly surrounded by dense mangrove forests (covering ca. 699 km2) [36, 51] and the intertidal flats are mainly composed by bare sediments, with some areas covered by large macroalgal mats (ca. 12% of the intertidal flats) and few rocky reefs (4% of the intertidal flats) [35, 36, 49, 52]. The vast adjacent mangrove forests, the continental runoff during the rainy season, and the upwelling from the northern coast [36], contribute to biological productivity of this system in ways that remain to be discovered (but see [44, 45]).
2.2. Macrozoobenthos sampling
Sampling was carried out in six sites across three islands (Fig 1): the island of Formosa, in a Marine Community Protected Area (Anrumai and Abu); the island of Bubaque, the most populated of the archipelago (Escadinhas, Bijante and Bruce); and in the islet of Adonga, inside Orango National Park (Adonga). Sites were sampled from the end of the wet season to the end of the dry season (October-April; Table 1), starting in January 2018 until April 2020 (months were not sampled consecutively within each year: S1 Table). This time of year was chosen as it overlaps with the presence of hundreds of thousands of migratory shorebirds in the archipelago, for which macroinvertebrates constitute a major part of their diet and are thus very likely to influence macrozoobenthos densities.
[Figure omitted. See PDF.]
Black circles in each intertidal site (insets) represent the locations where sediment cores were taken. Note the differences in size between each intertidal site when considering the extent and distribution of sampling. Grey areas represent land, dark grey areas represent mangroves, and digital elevation model is represented for the intertidal area (ranging from yellow–orange–red showing lower, mid, and higher elevation areas).
[Figure omitted. See PDF.]
Macrozoobenthos from four Classes (Polychaeta, Bivalvia, Gastropoda and Malacostraca) were collected using sediment cores of 0.00785 m2 of area for all sites except for Adonga, where cores had an area of 0.0113 m2, and max depth of 15 to 20 cm (often limited by rocks or shell banks). Sampling at all sites was conducted in areas where shorebirds where observed feeding. Except in one site (Adonga), polychaetes were not quantified in 2020, only bivalves, gastropods and malacostracans, as fieldwork came to a sudden halt due to the Covid-19 pandemic (Table 1, S1 Table). Core samples were divided in two fractions: the top 5 cm layer, where the majority of individuals are concentrated and therefore most of the harvestable biomass for intertidal predators, which was sieved through a 0.5 mm mesh size; and the bottom layer, which was sieved with a 1-mm mesh size. All individuals collected were stored in ethanol (70% or 96%), followed by later identification to the lowest taxonomic level possible with a magnifier lens, and different structures (see below) were measured to the closest mm. Identification was conducted with the support of a photographic reference collection built in this study (S1 Appendix), and also using identification keys for other areas published in books or journals [53–56]. Taxonomic names followed WoRMS (World Register of Marine Species; www.marinespecies.org).
2.3. Macrozoobenthos biomass estimation
A total of 26 taxa were considered for the biomass estimation, by including all taxa with densities above 0.8% of the total density of macrozoobenthos. The biomass of each individual found in almost all the core samples, was estimated from their size through taxa specific regression equations (Table 2). But for a small number of the collected cores (n = 20) the biomass of each individual was measured directly by following the methods described in Lourenço et al. (2017) [18], in which individuals were measured, dried to constant weight for 48h at 60°C and weighed (dry weight), and then incinerated in a muffle furnace for 2h at 500°C and weighted again (ash weight). The AFDM was calculated as the difference between dry weight and ash weight.
[Figure omitted. See PDF.]
Previously published equations relating size to biomass for the study area were used whenever available [18], and for the remaining taxa new equations were calculated from individuals of different sizes sampled in the field, following the previously described procedure for obtaining the ADFM of each individual. For some species, a group of specimens had to be pooled together by size-class (Table 2), in order to obtain a detectable amount of biomass after incineration, and others were grouped into a pre-defined size class (e.g., size between 10 and 15mm), for which the body size attributed to the individual was the average value of that size class (in this case 12.5mm). New regression equations were calculated by determining the nonlinear least-squares estimates of the parameters of the model AFDM = a*BSb [57–59] where BS is the Body Size in mm (antero-posterior length in bivalves and gastropods, total length in polychaetes and most malacostracans, and carapace width in crabs), using the nls function in R Software version 3.5.0 [60]. Whenever significant regression equations were not obtained for a certain taxon (i.e. size variation had no effect on biomass, in some cases potently due to limited sample size), an average AFDM value was used (regardless of the individual’s size) obtained from previously published studies or from incinerated individuals (Table 2).
2.4. Data analysis
Temporal variation was analysed using two variables: 1) month, a continuous variable that consisted of a pool of the samples of each sampled month in any year; and 2) period, a categorical variable that consisted in aggregating months into three main periods: end of wet season (October and November), early dry season (December, January and February) and late dry season (March and April; Table 1). Month was used for the univariate analysis and period for the multivariate analysis (see below). Aggregation into periods served the main purpose of increasing the sample size for the community analysis, allowing a more robust application of multivariate statistics. Those periods were selected because they represent three different environmental moments: the transitional period just after the rainy season (end of wet), with a stronger prevalence of the influence from rainfall (especially continental runoff) during the preceding months, which may promote inter-seasonal differences in the macrozoobenthic community [61, 62]; and the first (early dry) and last (late dry) months of the dry season, with the climate becoming increasingly drier and hotter. Furthermore, these three periods also represent three different moments in relation to the predation pressure exerted by migratory shorebirds on macrozoobenthos: the end of wet season, corresponding to the progressive arrival of shorebirds [37, 63] and thus relatively lower predation pressure on macrozoobenthos; the early dry season, corresponding to the period when shorebirds forage only to meet their daily energy requirements; and the late dry season, when shorebirds typically increase their daily intake in preparation for the migration towards the breeding areas [64–66]. In order to check that pooling of months of different years was acceptable, we tested for differences in density on samples collected in the same month (and site) in different years, using two-sample Wilcoxon tests. This showed that in the large majority (21 out 26 comparisons) of cases, there were no significant differences in the densities of macrozoobenthos between years (S2 Table).
In order to describe the spatial and temporal variation in overall density (ind.m-2) and biomass (mg of AFDM.m-2) of macrozoobenthos, and to assess if temporal patterns varied between sites, we conducted a Generalised Linear Model with negative binomial family (GLMnb), to account for overdispersion, and used the log link function. To model the variation in density, we used the number of individuals (all macrozoobenthos) counted in each core as the dependent variable, and to model variation in the biomass we used mg of AFDM in each core. The log of the area of each sampling core was used as an offset in all models, as this varied between sites, and the explanatory variables were site (with six levels) and month (continuous variable, coded sequentially from 1 to 7, starting in October). We also explored the interaction between the two explanatory variables. Model selection was undertaken using Akaike Information Criterion (AIC), by fitting all possible model combinations and selecting the model with the lowest AIC. When AIC differ less than two, a common “rule of thumb”, indicating no support for differences between competing models (see [67]), we selected the most parsimonious model. The same approach was used to describe the variation in spatial and temporal patterns for each of the major macrozoobenthic (sub)classes (Bivalvia, Polychaeta sedentaria, Polychaeta errantia, Gastropoda and Malacostraca). When differences between sites were significant, we performed post-hoc Tukey HSD pairwise contrasts to determine which sites were different from each other. When the interaction between site and month was significant for a given macrozoobenthos (sub)class, we computed individual GLMnb within each site to assess the significance of the variation of density and biomass along months.
We computed Shannon-Wiener diversity index and species richness for each site and overall to compare with previous studies in the Bijagós Archipelago and with macrozoobenthos communities in other tropical and temperate intertidal systems elsewhere.
In order to explore spatial variation in the community composition (i.e. the identity and relative abundances of all taxa in the community) and assess if temporal variation is similar between sites, we first calculated zero-adjusted Bray-Curtis dissimilarities [68], which summarizes differences between densities of taxa (ind.m-2) among samples (cores) with many zeros, using a log (x+1) transformation to reduce the weight of the most abundant species. We then used this dissimilarity matrix as the dependent variable on a permutational analysis of variance (PERMANOVA), with site (six levels) and period (three levels) as explanatory variables. We also tested the interaction between site and period and when the interaction was significant, we further inspected differences among sites separately for each level of period using PERMANOVAs, having site as the explanatory variable. Post-hoc PERMANOVAs with Bonferroni’s p-value adjustment for multiple comparisons were used to assess significant differences between sites. PERMANOVA’s assumption of homogeneity of multivariate dispersal was tested using betadisper and permutest functions of the VEGAN package, which simultaneously tested differences in community diversity between sites. In the event of heterogeneity of multivariate dispersal, PERMANOVA remains a robust test in the presence of balanced designs [69]. Nonetheless, in these cases we regarded the results of PERMANOVA with care, as the reasons for the differences may not be straightforward in some cases.
To visualize the differences in the centroids and in multivariate dispersal of the community composition among sites for each period, we computed Canonical Analysis of Principal Coordinates based on discriminant analysis (with three axis; hereafter referred to as CAP), using the function CAPdiscrim, from the BiodiverityR package [70] following the methods described in Anderson and Willis (2003) [71]. We used zero-adjusted Bray-Curtis dissimilarity matrix (based on densities) and the variable site for the constrained ordinations. Finally, to assess which species best characterized the differences in the community composition, we performed multiple linear regressions with a permutational test of significance of the fitted vectors, between the densities of each species (with a log (x+1) transformation) and the axis of the CAP ordination. For clarity, only the species with significant fits and R2≥ 0.15 were plotted as vector arrows over the CAP ordinations, with arrow length proportional to R2 value of the regression, thus informing on the importance of each species in the macrozoobenthic community. Plots showing all the species with significant fits, regardless of their R2, were also performed and presented in the Supporting Information. All the analysis were performed in R Software version 3.5.0 [60] and results are shown as mean ± SE, unless otherwise stated.
2.5. Inclusivity in global research
Permits to access field sites and sample were given by Instituto da Biodiversidade e das Áreas Protegidas Dr. Alfredo Simão da Silva, under the guardianship of the Ministry of Environment. Additional information regarding the ethical, cultural, and scientific considerations specific to inclusivity in global research is included in the Supporting Information (S1 Checklist).
3. Results
3.1. Spatial and temporal variation in density and biomass of macrozoobenthos
From a total of 9907 individuals sampled in 609 sediment cores, the Bijagós intertidal flats had an average macrozoobenthos density of 1871 ± 58.3 ind.m-2, corresponding to a mean biomass of 5.65 ± 0.41 gAFDM.m-2 (Table 1; S3 Table). Sedentary polychaetes were the most abundant (sub)class, representing 41% of total density, followed by bivalves (24%), errant polychaetes (17%), malacostracans (8%) and gastropods (8%). Regarding biomass, bivalves were the predominant class, representing 72% of the total biomass, followed by sedentary polychaetes (17%), malacostracans (7%), errant polychaetes (4%) and lastly gastropods, which represented less than 1% of the total (Fig 2). When aggregating all months, mean densities and biomass varied markedly across sites for each macrozoobenthic (sub)class (Fig 2).
[Figure omitted. See PDF.]
Map of the Bijagós archipelago (inset) is shown, with indication of the spatial distribution of the study sites, represented by coloured squares with matching colour scheme.
3.1.1. Variation in macrozoobenthos density.
Overall density of macrozoobenthos varied significantly among sites (p<0.001), with Adonga having the highest mean density, significantly higher than all sites except Bijante, and followed by Abu, Anrumai, Bruce and Escadinhas, the later having significantly lower densities than all other sites (Table 3; S3 and S4 Tables). Overall density increased significantly throughout months (p<0.001; Table 3), and the interaction between month and site was not significant (interaction term, p = 0.07). Removing Adonga, the site with only three months sampled, from the GLM analysis presented in Table 3, yielded the same overall results.
[Figure omitted. See PDF.]
At the (sub)class level, Adonga had highest mean density of bivalves and sedentary polychaetes (the later together with Bijante), but lowest of gastropods and malacostracans (S3 and S4 Tables). Bruce had the highest densities of malacostracans, but the lowest of bivalves (together with Escadinhas) and errant polychaetes (S3 and S4 Tables). This polychaete subclass had its highest density in Abu, where gastropods were also found to be the highest (together with Bijante), while Abu had the lowest densities of sedentary polychaetes, together with Anrumai (S3 and S4 Tables).
Bivalve density increased significantly across months regardless of site (Table 3, Fig 3). Sedentary and errant polychaetes also increased significantly in Escadinhas (Fig 3) with the former showing the same pattern on Bruce and the latter on Bijante. Malacostracans and gastropods increased in Bruce and Bijante, with the later also increasing in Adonga (Fig 3).
[Figure omitted. See PDF.]
Black circles and vertical lines represent mean with associated standard errors. Black lines represent Generalized linear models (family negative binomial, with log link) fitted to data points, with grey areas showing 95% confidence intervals. Significant relationships are represented by solid lines, and non-significant by dashed lines. Estimates for the variation of density for each unit of month and p-values are shown in each panel. Bivalves are the only (sub)class for which no significant interaction was found between site and month (Table 3), and therefore has only one estimate of month and p-value for all sites. Regressions lines for Adonga were not computed due to only having 3 sampled months.
3.1.2. Variation in macrozoobenthos biomass.
The biomass of macrozoobenthos varied significantly among sites (p<0.001), with Abu, Adonga, and Anrumai having significantly higher biomass than all other sites, followed by Bijante, with Bruce and Escadinhas having significantly less biomass than the remaining sites (Table 3; S3 and S4 Tables). Overall, biomass did not vary significantly throughout months (p = 0.997), and the interaction term was also not significant (interaction term, p = 0.16), even when removing Adonga from the analysis.
At the (sub)class level, Abu, Adonga and Anrumai had the highest biomass of bivalves, with the last two sites having the lowest biomass values of gastropod and malacostracans, respectively (S3 and S4 Tables). The highest biomass of malacostracans was found in Escadinhas and of errant polychaete in Abu (S3 and S4 Tables).
Bivalve biomass did not vary across months, whereas the remaining (sub)classes showed different patterns of monthly variation, the majority of which were positive (Table 3, Fig 4). Both sedentary and errant polychaete biomass increased significantly in Bijante, but decreased in Adonga (sedentary) and Anrumai (errantia; Table 3, Fig 4). The later site also showed a decrease in malacostracan biomass, which in turn increased in three other sites (Abu, Bruce and Adonga). Gastropod was the only class with an increase in biomass regardless of sites (Table 3, Fig 4).
[Figure omitted. See PDF.]
Black circles and vertical lines represent mean and standard errors, respectively. Black lines represent Generalized linear models (family negative binomial, with log link) fitted to data points, with grey areas showing 95% confidence intervals. Significant relationships are represented by solid lines, and non-significant by dashed lines. Estimates for the variation of density for each unit of month and p-values are shown in each panel. Bivalves and Gastropods are the only (sub)classes for which no significant interaction was found between site and month (Table 3), and therefore has only one estimate of month and p-value for all sites. Regressions lines for Adonga were not computed due to only having 3 sampled months.
3.2. Variation in macrozoobenthos community composition and diversity
Overall, at least 88 taxa were identified in all sites, at different taxonomical levels (S3 Table). To assess the diversity, species richness and community composition (analysed based on the density of each species), these 88 taxa were aggregated into 82 to exclude taxa repetitions along different taxonomical levels (e.g. when some individuals of a given family were possible to identify at the species level, but most were not, all were aggregated at the family level). Overall, we found a Shannon-Wiener of 3.3 when aggregating samples from all sites. Site specific species richness and diversity are presented in Supporting Information (Table 1 in S2 Appendix).
3.2.1. Overall spatial and temporal differences.
The macrozoobenthic community composition showed overall significant differences among sites (PERMANOVA F5,507 = 22.76, R2 = 0.17, p<0.001) and periods (but with a small amount of explained variance; PERMANOVA F2,507 = 6.4, R2 = 0.02, p<0.001), and the variation between periods was significantly different among the sites (PERMANOVA interaction term, F10,507 = 2.6, R2 = 0.04, p<0.001). CAP ordinations of sites suggest that Adonga has the most different community composition, mostly due to the high densities of the bivalves P. isocardia, A. nymphalis, Keletistes sp. and T. adansonii (Fig 5A). The errant polychaete Sigambra sp. appeared to be especially abundant in Abu and Escadinhas, while a higher density of the sedentary polychaete Orbiniidae seemed to be responsible for the community composition differences in Bruce. Multivariate dispersions (community diversity) were significantly different across the sites (Permutation test, F5,519 = 20.2, p = 0.001), with Adonga as the only site with significantly smaller dispersion than all the others (S5 Table). Temporal patterns are unclear in the overall CAP ordinations and difficult to interpret at such wide level (Fig 5B). Furthermore, there were no significant differences in multivariate dispersions between periods (Permutation test, F2,522 = 2.57, p = 0.087).
[Figure omitted. See PDF.]
Ellipses represent 60% of the data for each group. Arrows show species with the most influence in the distribution of data points, and the size of the arrows are scaled according to the R2 of each regression between the densities of each species and the axis of the CAP ordination (i.e., longer arrows represent more influent species). Species with significant fits but with an R2 smaller than 0.15 are omitted.
3.2.2. Differences between sites in each period.
When comparing the community composition in each period, there were significant differences between sites in all periods (PERMANOVAs, end of wet season: F5,146 = 10.71, R2 = 0.27, p = 0.001; early dry season: F5,204 = 6.32, R2 = 0.13, p = 0.001; late dry season: F5,157 = 11.48, R2 = 0.27, p = 0.001). The same holds for differences between sites in the multivariate dispersion (permutation tests, end of wet season: F5,146 = 8.85, p = 0.001; early dry season: F5,204 = 8.00, p = 0.001; late dry season: F5,157 = 8.41, p = 0.001).
In the end of the wet season (October and November) the most similar communities were those in Anrumai and Abu (S6 Table). Adonga was the site with the most different community composition, influenced by the importance of the bivalves P. isocardia, Keletistes sp. and A. nymphalis. In Abu the sedentary polychaetes Paraonidae and Orbiniidae were important families to differentiate this community from the others, while the Orbiniidae were also important in Bruce and Bijante (Fig 6A). The significant differences in multivariate dispersal at the end of wet period were mostly due to Bruce and Adonga having significantly smaller dispersion in their communities in comparison to the other sites (S5 Table; Fig 6A).
[Figure omitted. See PDF.]
Ellipses represent 60% of the data for each group. Arrows show species with the most influence in the distribution of data points, and the size of the arrows is scaled according to the R2 of each regression (i.e., longer arrows represent more influential species). Species with a significant fit but with an R2 smaller than 0.15 are omitted. A map of the Bijagós Archipelago is shown, with indications of the study site locations represented as squares with a matching colour scheme.
In early dry season (December, January and February), the communities in Adonga remained the most distinct (S6 Table, Fig 6B). In this period, however, it is noticeable that the sedentary polychaete Maldanidae replaced the bivalves that dominated the community composition of Adonga in the previous period, especially P. isocardia (Fig 6B). This bivalve was abundant in this period in Anrumai instead, increasing markedly its importance at this site, together with A. nymphalis, although to a lesser extent. In Bruce, Orbiniidae maintained their importance in comparison to the previous period, while increasing in Escadinhas. In the early dry season, the multivariate dispersal of the macrozoobenthic community in Adonga was also significantly smaller than in all other sites, in which dispersal increased instead (S5 Table, Fig 6B).
Finally, in late dry season, most sites remain significantly different from each other, except for Anrumai and Abu, as in the end of wet season. Adonga also remained the most different site (S6 Table, Fig 6C), and similarly to the end of wet season, is strongly influenced by the bivalves P. isocardia and A. nymphalis, which increased markedly from early dry season to late dry season, as well as the bivalve T. adansonii, that became important in Adonga only in late dry season. In Bruce, Orbiniidae maintained its importance from early to late dry season, while in Bijante Orbiniidae’s importance increased to a similar importance as at the end of the wet season. The community in Escadinhas shifted to be dominated by sedentary polychaetes (including Capitellidae) and the errant polychaete Sigambra sp. The multivariate dispersion also differs significantly, mostly because of the significantly higher dispersion in Anrumai and Abu in this final period (S5 Table).
4. Discussion
In this study we undertook a large sampling effort over the intertidal flats of the Bijagós Archipelago, covering six sites and seven months (n = 609 cores), to describe the spatial and temporal variation of macrozoobenthos density, biomass, and community diversity and composition. We found a highly diverse community, with more than 88 taxa, but with typically modest densities and biomass. With both univariate and multivariate approaches, we showed significant differences in the macrozoobenthos density and community structure among sites. Interestingly, overall density increased linearly throughout months, with all significant relationships being positive across sites and (sub)classes. However, contrarily to the expected reduction in biomass due to shorebird predation in this period, we found no clear overall trends for biomass when aggregating all sites, and at the site level most of the significant relationships across sub(classes) were positive (although two were not). These patterns dispute the expected depletion in the macrozoobenthos community by the end of the dry season due to intense predation pressure from migratory shorebirds, which adds to the local fish predation over the studied period. The findings are at odds with the observation at Banc d’Arguin, Mauritania, where depletions over the shorebird predation season did seem to occur [29].
4.1. Overall patterns of macrozoobenthos communities
On average, macrozoobenthos of the Bijagós Archipelago had a density of 1871 ± 58 ind.m-2, corresponding to an average biomass of 5.65 ± 0.41 gAFDM.m-2. These values are in accordance with previous studies done locally at smaller spatial and temporal frames [46, 72]. But while density values are within the expected range for tropical sites [20], our biomass values were overall much lower when compared to the nearby intertidal area of the Banc d’Arguin, which wielded 17.0 gAFDM.m-2 [57]. Even larger differences are observed when compared to biomass values reported for temperate intertidal areas within the East Atlantic Flyway, as the Tagus estuary, or the Wadden Sea (with 21.9 g of Dry Mass.m-2 and 26.6 gAFDM.m-2 respectively) [73, 74].
We found that sedentary polychaetes, the most abundant taxonomic group (41% of total density), accounted for only 17% of the total macrozoobenthos biomass, whereas the second most common group, bivalves (27% of total density), had a disproportionally high contribution towards biomass, representing 72% of the total. This pattern of bivalves occurring at low densities but high biomass is typical in other intertidal systems, both temperate, such as the Tagus estuary where the bivalve Scrobicularia plana is responsible for most of the biomass [75, 76], and tropical, like in South Sumatra in Indonesia where the Arcidae bivalve Anadara granosa takes that role [77]. In our study, a large proportion of bivalve biomass (46%) is composed of individuals of S. senilis (Arcidae family), which is in line with previous studies in the Bijagós Archipelago [45, 46] and is similar to the assemblages in Banc d’Arguin [29, 78]. This species is important for example in the diet of Oystercatchers [79] and some specialized rays [80], and although it is not consumed by most shorebirds [18], on some occasions it can also be part of the diet of some of these species, like for Whimbrel [81], while also being commercially harvested and one of the major protein sources for local people [82].
Tropical intertidal flats are known to often show higher macrozoobenthos diversity than their temperate counterparts [11, 20, 33], and the Bijagós Archipelago is no exception. Overall, we found that species richness (88 taxa) and diversity (Shannon-Wiener index of 3.3) were considerably higher than previously reported values for this area. Lourenço et al. (2018) [46] reported a species richness of 40 taxa with a Shannon-Wiener index of 2.89, whilst Meijer et al. (2021) [45] found 48 taxa. These differences are likely due to the more restricted sampling presented in those studies, which were also more constrained in time and/or space, even though in both studies all sites sampled in the present study were included. It should nevertheless be noted that the accumulation curves produced in our study (Fig 1 of S2 Appendix) show an overall rapid increase in the number of taxa found for each additional core and despite our larger sampling effort none has reached the asymptote. This suggests that further sampling may render even higher values of diversity. Furthermore, species richness in another intertidal system in the region, the Banc d’Arguin, in Mauritania, is reported to be higher (111 taxa) [83] despite intrinsic differences between systems, it could be expected that more species may be found in the Bijagós, particularly in a sampling scheme that covers other periods of the annual cycle.
The macrozoobenthic community composition was shown to be overall well segregated in the CAP space, with the community structure gradients particularly dominated by the differences in the occurrence and density between eight key taxa. The first axis of the CAP ordination showed mainly the distribution of sampling cores along the gradients of abundance of the bivalve species, while the second axis mainly reflected the gradient along the abundance of the polychaete taxa. This suggests a diverse community that depends on several important taxa, an important trait of stable and resilient ecosystems [84, 85]. These results partially agree with those reported by Lourenço et al. (2018) [46], in which also several (ca. 9) polychaete worms and bivalves were statistically explaining the community structure gradients. However, a notable difference is the presence of the fiddler crab Afruca tangeri as a key species in Lourenço et al. (2018) [46], which was not the case in our work. This discrepancy is likely related to methodological differences: we sampled macrozoobenthos using exclusively core samplers (which tend to miss crabs for being considerably mobile) whereas Lourenço et al. (2018) [46] conducted a separate method specific for sampling fiddler crabs. This likely resulted in an underestimation of their importance in our study. Nevertheless, fiddler crabs are ecosystem engineers within the intertidal flats of the Bijagós Archipelago, and while not being direct predators of macrozoobenthos [86], they might negatively impact macrozoobenthos abundance and structure the community composition of both benthic invertebrates and shorebird consumers [28]. Thus, their inclusion within the lot of key species is warranted, as previously suggested by Lourenço et al. (2018) [46].
4.2. Spatial variation of macrozoobenthos communities
We showed overall significant spatial variation in both density and biomass of macrozoobenthos when comparing sites, and this was also true at the (sub)class level for the majority of macrozoobenthos. Considering density, when aggregating all macrozoobenthos (sub)classes, four sites stand out, Adonga and Bijante with comparatively high mean densities (2507 ± 157 ind.m-2 and 2321 ± 149 ind.m-2, respectively), and Escadinhas and Bruce with lower densities (1195 ± 100 in.m-2 and 1633 ± 119 in.m-2, respectively), than all other sites. These patterns were the same for Adonga, Escadinhas and Bruce when considering mean biomass (8.37 ± 0.83 g AFDM m-2, 2.21 ± 0.4 g AFDM m-2 and 2.57 ± 0.47 g AFDM m-2, respectively), but the high density found in Bijante did not translate into high mean biomass (4.68 ± 0.67 g AFDM m-2). This is likely due to the comparatively lower prevalence of bivalves in this site (which represented the main source of biomass in the whole system), in comparison to other sites like Anrumai, Abu and Adonga, and which in Bijante were replaced by higher density of sedentary polychaetes. Additionally, the mean size of bivalves in Bijante was the lowest of all sites (5.0 ± 0.2 cm in Bijante vs 7.6 ± 0.3 cm in Adonga and 5.5 cm for both Anrumai and Abu, with SE of 0.4 and 0.3 respectively), which also contributes to lower overall biomass. Lourenço et al. (2018) [46] also found similar patterns in biomass with Bijante having significantly higher values than Bruce and Escadinhas. However, the mean biomass values reported in that study for Bijante and Bruce in January were considerably higher than our findings for the same month (8.5 vs 3.71 g AFDM.m-2 in Bijante; 5.5 vs 2.08 g AFDM m-2 in Bruce), while for Escadinhas, values were similar. These differences may derive from the different methods used to estimate biomass (direct measure vs. estimates from equations in our study) and/or a lower contribution from fiddler crabs in our study because of the sampling method used.
The spatial variation at the (sub)class level in terms of density and biomass was translated in marked differences in the most important (sub)classes between the sites. Furthermore, spatial differences were also apparent in the results of the multivariate analysis, with overall significant differences between sites in the community composition and in community diversity (multivariate dispersal). The aggregated results of all analysis clearly indicate that spatial variation needs to be taken into account to understand the intertidal ecosystem of the Bijagós Archipelago. Spatially heterogeneous intertidal landscapes have been reported in tropical systems elsewhere [83, 87, 88], and in other systems like in the Yellow Sea [89], and in the Wadden Sea [21]. But in temperate and northern intertidal areas, these spatial differences in the macrozoobenthic community structure are typically driven by a few abundant species (e.g. [21, 75, 89, 90]) unlike what was found in the Bijagós Archipelago.
Various macrozoobenthos taxa were found to have high importance in characterizing the community composition of each site. In Adonga, the bivalves P. isocardia, Keletistes sp., A. nymphalis and T. adansonii were the main taxa driving the marked differences to the community compositions of the other sites. Pelecyora isocardia is a small sized bivalve that constitutes the most important prey item for Red knots Calidris canutus in the Bijagós Archipelago, totalling more than 80% of the proportion of biomass consumed by this shorebird [18]. Pelecyora isocardia is also important for Bar-tailed godwits Limosa lapponica, although to a lesser extent, for which it was found as a prey item in 30% of the droppings analysed in the Bijagós Archipelago [18]. Also, T. adansonii and T. nymphalis are consumed to some extent by Red knots and Bar-tailed godwits [18]. Sites like Escadinhas, Bruce and Abu were shown to have macrozoobenthic communities with high influence from different polychaete taxa, which comprise important prey items for several shorebird species [18, 72, 91] and for both juvenile and adult benthic fish [19, 92]. Due to the differences in food preference between different species of these intertidal consumers, the spatial variation in the macrozoobenthic communities found in the Bijagós Archipelago is likely to promote differential space use between the species of these high-level consumers [93–95].
In the CAP analysis we compared the community composition and diversity among sites for each period (i.e. end of wet, early dry and late dry), with significant differences between most sites emerging in all periods. It is noteworthy that Adonga was the site showing the largest separation from the other sites, with the macrozoobenthic community being characterized by high densities of bivalve taxa, and a significantly lower community variability (multivariate dispersal). Moreover, when considering the weighted mean Shannon-Wiener diversity index and species richness in each site, Adonga was also the least diverse and poorest site (although these results should be interpreted cautiously, as sampling there had limited temporal representation; Table 1 in S2 Appendix). On the other hand, Anrumai and Abu were the most similar in the CAP analysis, both in community structure and multivariate dispersal, and also displayed virtually identical species richness and Shannon-Wiener index, having only showed significant differences between each other in the community structure during the early dry season. The similarities between these two sites are also supported by the comparisons between their density and biomass, which never differed significantly for any macrozoobenthos (sub)class. This suggests that these two sites may have similarities in key biotic and abiotic habitat traits and, therefore, in shorebird and benthic fish species assemblages preying on their macrozoobenthic communities. These sites are also the two closest in our sampling, which can be an important factor determining the similarities between their macrozoobenthic communities.
4.3. Temporal variation of macrozoobenthos communities
Interestingly, the overall density of macrozoobenthos increased from the end of the wet season to the end of the dry season. After the wet season, the number of individuals started to increase despite the presence of their seasonal predators (shorebirds), suggesting recruitment at this stage [34]. This recruitment hypothesis is further supported by the fact that some species had the highest average body sizes in the late wet season, suggesting the appearance of smaller sized individuals in the population during the course of the dry season. For example, the largest shell sizes of the bivalve Pelecyora isocardia, were found at the end of the wet season (7.0 ± 0.2 mm), and the smallest in the dry season (2.5 ± 0.1 mm in the beginning and 3.9 ± 0.1 mm in the end of the season) and the same happened with body sizes of the Orbiniidae polychaetes (20.3 ± 1.3 mm in the end of the wet season vs 17.7 ± 1.1 and 17.5 ± 1.2 in early and late dry season, respectively). These results contrast with the findings from Banc d’Arguin, where the timing of the wet season is similar to that of the Bijagós, but where, next to the Sahara, the magnitude of the wet is of course much smaller. At the Banc d’Arguin intertidal flats most bivalve species densities peaked in October and declined until March, declines which were argued to be explained by depletion by molluscivore shorebirds [29]. In the Bijagós mudflats no significant decreases in density in the course of the shorebird predation season (October-March) were detected, with no evidence for depletion. We propose that the heavy rainfall during the wet season decreases water salinity to an extent that it increases the physiological stress on macrozoobenthos and induces mortality [96]. The subsequent recruitment and high growth rates over the following dry season would, however, ‘overrule’ any depletion effects shown by the (rather more abundant) shorebirds at Banc d’Arguin.
Overall biomass did not vary significantly through time. However, these results must be interpreted with caution, as they were extremely influenced by the dominance of bivalve biomass (72% of total biomass), which remained stable throughout the sampling season. Regardless, no clear pattern in the variation of the majority of (sub)classes emerged, with some sites showing a significant increase in biomass over time, while others showed a decrease. Only gastropod biomass increased consistently across sites, even though values were overall very low. This is likely due to the sampling strategy, as the larger gastropod species have very sparse distribution (e.g. Pugilina morio), while some additionally show a preference for mudflat channels (e.g. Cymbium sp). Biomass at temperate intertidal sites is typically higher in spring, when macrozoobenthic species’ growth rates are high due to increasing temperatures, and is followed by a decrease until winter [76, 97]. But at tropical sites temperature tends to be more stable throughout the year, and rainfall patterns seem to influence temporal variation the most [96]. Nonetheless, in this study no clear temporal variation of biomass was detected, which may indicate that macrozoobenthos biomass is rather stable throughout the year in the Bijagós Archipelago, as it has been suggested for other tropical sites [98].
We found no clear patterns when assessing the temporal variation in the community composition overall and within each site. While in some sites there was some variation in the most important taxa between periods, the absence of clear patterns between periods is an indication that communities at each site remained rather stable throughout time, and that the spatial variability was higher than the temporal variability, maybe due to site-specific abiotic factors that were not explored in this study, as demonstrated for macroinvertebrate communities elsewhere [76, 99].
The observed increase in density (potentially due to recruitment) and the absence of a decrease in biomass throughout the study period, which coincides with the presence of large numbers of migratory shorebirds wintering in the Bijagós Archipelago, together with the fact that communities remain rather stable between periods, suggests that despite their large numbers, shorebirds do not cause depletion of macroinvertebrates at this site.
5. Conclusion
In this work we tackled an evident knowledge gap, by providing a detailed description of the composition, density and biomass of macrozoobenthic communities at a wide spatial-temporal scale within the Bijagós Archipelago. This constitutes an important contribution for the understanding of the intertidal ecosystem of this key biodiversity site, which harbours important benthic fish communities and is a key non-breeding area for Palearctic migratory shorebirds, for which macrozoobenthos constitutes an important prey group. This baseline knowledge will be relevant for subsequent ecological studies in this area focused on the populations of these high-level consumers and on further assessing the interactions within the whole intertidal food web.
Given the high sample size and wide spatial and temporal representativeness, the main findings of this study allow for the conclusion that the macrozoobenthic community of the Bijagós Archipelago is diverse but has relatively low biomass values, is highly variable between sites, and presents densities that increase towards the end of the dry season, despite the high predation pressure by migratory shorebirds during this period. Furthermore, our results suggested that there is a high number of important macrozoobenthic taxa, possibly indicating a resilient intertidal ecosystem. These traits may well be among the important services that the Bijagós Archipelago offers and that may help explain the regional and international importance of this site, for benthic fish in the West African region, and within the East Atlantic Flyway for migratory shorebirds.
Supporting information
S1 Table. Number of samples (sediment cores) collected in each month across sites.
https://doi.org/10.1371/journal.pone.0277861.s001
(DOCX)
S2 Table. Inter-annual comparisons of density and biomass.
https://doi.org/10.1371/journal.pone.0277861.s002
(DOCX)
S3 Table. Mean density (ind.m2) and biomass (mg of AFDW.m2) and standard error, of the macrozoobenthos taxa found in 609 sediment core samples collected in six intertidal flats of the Bijagós Archipelago.
https://doi.org/10.1371/journal.pone.0277861.s003
(DOCX)
S4 Table. Post-hoc pairwise contrasts of density and biomass between in each site.
https://doi.org/10.1371/journal.pone.0277861.s004
(DOCX)
S5 Table. Post-hoc Tukey HSD tests between sites for multivariate dispersion.
https://doi.org/10.1371/journal.pone.0277861.s005
(DOCX)
S6 Table. Post-hoc PERMANOVAs comparing sites within each period.
https://doi.org/10.1371/journal.pone.0277861.s006
(DOCX)
S1 Fig. Canonical analysis of principal components based on discriminant analysis with all species displayed.
https://doi.org/10.1371/journal.pone.0277861.s007
(DOCX)
S1 Appendix. Macrozoobenthos photographic reference collection.
https://doi.org/10.1371/journal.pone.0277861.s008
(DOCX)
S2 Appendix. Spatial variation in species richness and diversity.
https://doi.org/10.1371/journal.pone.0277861.s009
(DOCX)
S1 Checklist. Questionnaire on inclusivity in global research.
https://doi.org/10.1371/journal.pone.0277861.s010
(DOCX)
Acknowledgments
We are grateful to Instituto da Biodiversidade e das Áreas Protegidas Dr. Alfredo Simão da Silva for their permission to work in the Bijagós and all the support with the logistics and transportation within and between the islands, in particular to Justino Biai, Quintino Tchantchalam, Emanuel Dias, Samuel Pontes and Saído Gomes. We are also grateful to NGO Tiniguena for the logistical support on the island of Formosa, particularly to Emanuel Ramos and Nfali. We thank Paula Lopes, Filipe Moniz, Camilo Carneiro, Josh Nightingale, João Belo, Jorge Gutierrez, Edna Correia, Rui Rebelo, Maio Onemo, Seco Leão, Baba Germano, Tchandú and Agostinho for help with fieldwork, while Lotte Neimeijer, Bianka Rasch, Pedro Lourenço and Serge Gofas contributed for the benthos identification.
Citation: Coelho AP, Henriques M, Rocha AD, Paulino J, Schaars LK, Ramos C, et al. (2022) Spatial and seasonal variation in macrozoobenthic density, biomass and community composition in a major tropical intertidal area, the Bijagós Archipelago, West-Africa. PLoS ONE 17(11): e0277861. https://doi.org/10.1371/journal.pone.0277861
About the Authors:
Ana Pinto Coelho
Contributed equally to this work with: Ana Pinto Coelho, Mohamed Henriques
Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft
E-mail: [email protected]
Affiliation: Centre for Environmental and Marine Studies (CESAM), Department of Biology, University of Aveiro, Aveiro, Portugal
https://orcid.org/0000-0002-3337-0310
Mohamed Henriques
Contributed equally to this work with: Ana Pinto Coelho, Mohamed Henriques
Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft
Affiliations Centre for Environmental and Marine Studies (CESAM), Department of Animal Biology, Faculdade de Ciências da Universidade de Lisboa, Lisbon, Portugal, Conservation Ecology Group, Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, The Netherlands
https://orcid.org/0000-0002-3061-5423
Afonso Duarte Rocha
Roles: Investigation, Methodology, Writing – review & editing
Affiliations Centre for Environmental and Marine Studies (CESAM), Department of Biology, University of Aveiro, Aveiro, Portugal, Conservation Biology Research Group, Department of Anatomy, Cell Biology and Zoology, University of Extremadura, Badajoz, Spain
João Paulino
Roles: Data curation, Investigation, Methodology, Writing – review & editing
Affiliation: Centre for Environmental and Marine Studies (CESAM), Department of Animal Biology, Faculdade de Ciências da Universidade de Lisboa, Lisbon, Portugal
https://orcid.org/0000-0003-0779-1960
Loran Kleine Schaars
Roles: Data curation, Investigation, Methodology, Resources, Writing – review & editing
Affiliation: Department of Coastal Systems, NIOZ Royal Netherlands Institute for Sea Research, Texel, The Netherlands
Catarina Ramos
Roles: Data curation, Investigation, Writing – review & editing
Affiliation: Centre for Environmental and Marine Studies (CESAM), Department of Biology, University of Aveiro, Aveiro, Portugal
Aissa Regalla de Barros
Roles: Funding acquisition, Project administration, Resources, Writing – review & editing
Affiliation: Instituto da Biodiversidade e das Áreas Protegidas Dr. Alfredo Simão da Silva (IBAP), Bissau, Guiné-Bissau
Teresa Catry
Roles: Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing
Affiliation: Centre for Environmental and Marine Studies (CESAM), Department of Animal Biology, Faculdade de Ciências da Universidade de Lisboa, Lisbon, Portugal
José Pedro Granadeiro
Roles: Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing
Affiliation: Centre for Environmental and Marine Studies (CESAM), Department of Animal Biology, Faculdade de Ciências da Universidade de Lisboa, Lisbon, Portugal
https://orcid.org/0000-0002-7207-3474
Theunis Piersma
Roles: Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing
Affiliations Conservation Ecology Group, Groningen Institute for Evolutionary Life Sciences, University of Groningen, Groningen, The Netherlands, Department of Coastal Systems, NIOZ Royal Netherlands Institute for Sea Research, Texel, The Netherlands
José Augusto Alves
Roles: Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing
Affiliations Centre for Environmental and Marine Studies (CESAM), Department of Biology, University of Aveiro, Aveiro, Portugal, South Iceland Research Centre, University of Iceland, Laugarvatn, Iceland
1. Finlayson C, Milton G, Prentice R, Davidson N. The wetland book. II. Distribution, description and conservation. Dordrecht: Springer; 2017.
2. Beninger PG. Mudflat Ecology. Aquatic Ec. Beninger PG, editor. Cham: Springer International Publishing; 2018. https://doi.org/10.1007/978-3-319-99194-8
3. Murray NJ, Phinn SR, DeWitt M, Ferrari R, Johnston R, Lyons MB, et al. The global distribution and trajectory of tidal flats. Nature. 2019;565: 222–225. pmid:30568300
4. Arkema KK, Guannel G, Verutes G, Wood SA, Guerry A, Ruckelshaus M, et al. Coastal habitats shield people and property from sea-level rise and storms. 2013;3. https://doi.org/10.1038/NCLIMATE1944
5. Bernhardt JR, Leslie HM. Resilience to climate change in coastal marine ecosystems. https://doi.org/10.1146/annurev-marine-121211-172411 pmid:22809195
6. Spalding MD, Ruffo S, Lacambra C, Meliane I, Zeitlin L, Shepard CC, et al. The role of ecosystems in coastal protection: Adapting to climate change and coastal hazards. Ocean Coast Manag. 2014;90: 50–57.
7. Harley CDG, Hughes AR, Hultgren KM, Miner BG, Sorte CJB, Thornber CS, et al. The impacts of climate change in coastal marine systems. Ecol Lett. 2006;9: 228–241. pmid:16958887
8. Piersma T. Threats to intertidal soft-sediment ecosystems. In: Reinhard S, Folmer H, editors. Water Policy in the Netherlands Integrated Management in a Densely Populated Delta. Washington, D.C: Resources for the Future; 2009. pp. 57–69.
9. Galbraith AH, Jones R, Park R, Clough J, Harrington B, Page G, et al. Global climate change and sea level rise: potential losses of intertidal habitat for shorebirds. 25;2:173–183.
10. Eriksson BK, van der Heide T, van de Koppel J, Piersma T, van der Veer HW, Olff H. Major changes in the ecology of the wadden sea: Human impacts, ecosystem engineering and sediment dynamics. Ecosystems. 2010;13: 752–764.
11. Piersma T, de Goeij P, Tulp I. An evaluation of intertidal feeding habitats from a shorebird perspective: towards relevant comparisons between temperate and tropical mudflats. Netherlands J Sea Res. 1993;31: 503–512.
12. Hamilton DJ, Diamond AW, Wells PG. Shorebirds, snails, and the amphipod (Corophium volutator) in the upper Bay of Fundy: top–down vs. bottom–up factors, and the influence of compensatory interactions on mudflat ecology. 2006; 285–306.
13. Thrush SF, Hewitt JE, Gibbs M, Lundquist C, Norkko A. Functional role of large organisms in intertidal communities: community effects and ecosystem function. 2006; 1029–1040.
14. Van Gils JA, Van Der Geest M, Jansen EJ, Govers LL, De Fouw J, Piersma T. Trophic cascade induced by molluscivore predator alters pore-water biogeochemistry via competitive release of prey. Ecology. 2012;93: 1143–1152. pmid:22764500
15. Cheverie AV, Hamilton DJ, Coffin MRS, Barbeau MA. Effects of shorebird predation and snail abundance on an intertidal mudflat community. J Sea Res. 2014;92: 102–114.
16. Christianen MJA, Middelburg JJ, Holthuijsen SJ, Jouta J, Compton TJ, van der Heide T, et al. Benthic primary producers are key to sustain the Wadden Sea food web: stable carbon isotope analysis at landscape scale. Ecology. 2017;98: 1498–1512. pmid:28369845
17. van de Kam J, Ens B, Piersma T, Zwarts L. Shorebirds: an illustrated behavioural ecology. Utrecht, The Netherlands: KNNV Publishers; 2004.
18. Lourenço PM, Catry T, Granadeiro JP. Diet and feeding ecology of the wintering shorebird assemblage in the Bijagós archipelago, Guinea-Bissau. J Sea Res. 2017;128: 52–60.
19. Dolbeth M, Martinho F, Leitão R, Cabral H, Pardal MA. Feeding patterns of the dominant benthic and demersal fish community in a temperate estuary. J Fish Biol. 2008;72: 2500–2517.
20. Dittmann S. Benthic fauna in tropical tidal flats—A comparative perspective. Wetl Ecol Manag. 2002;10: 189–195.
21. Compton TJ, Holthuijsen S, Koolhaas A, Dekinga A, ten Horn J, Smith J, et al. Distinctly variable mudscapes: Distribution gradients of intertidal macrofauna across the Dutch Wadden Sea. J Sea Res. 2013;82: 103–116.
22. Matias MG, Arenas F, Rubal M, Pinto IS. Macroalgal composition determines the structure of benthic assemblages colonizing fragmented habitats. PLoS One. 2015;10. pmid:26554924
23. Li K, Liu X, Zhou Y, Xu Y, Lv Q, Ouyang S, et al. Temporal and spatial changes in macrozoobenthos diversity in Poyang Lake Basin, China. Ecol Evol. 2019; 6353–6365. pmid:31236226
24. Dittmann S, Vargas JA. Tropical tidal flat benthos compared between Australia and central America. 2001;151: 275–293.
25. Rosa S, Granadeiro JP, Vinagre C, França S, Cabral HN, Palmeirim JM. Impact of predation on the polychaete Hediste diversicolor in estuarine intertidal flats. 2008;78: 655–664.
26. Passarelli C, Olivier F, Paterson DM, Meziane T, Hubas C. Organisms as cooperative ecosystem engineers in intertidal flats. J Sea Res. 2014;92: 92–101.
27. Umanzor S, Ladah L, Calderon-Aguilera L, Zertuche-González J. Intertidal macroalgae influence macroinvertebrate distribution across stress scenarios. Mar Ecol Prog Ser. 2017;584: 67–77.
28. Paulino J, Granadeiro JP, Henriques M, Belo J, Catry T. Composition and abundance of shorebird and macroinvertebrate communities differ according to densities of burrowing fiddler crabs in tropical intertidal flats. Hydrobiologia. 2021;848: 3905–3919.
29. Ahmedou Salem MV, van der Geest M, Piersma T, Saoud Y, van Gils JA. Seasonal changes in mollusc abundance in a tropical intertidal ecosystem, Banc d’Arguin (Mauritania): Testing the “depletion by shorebirds” hypothesis. Estuar Coast Shelf Sci. 2014;136: 26–34.
30. Oudman T, Piersma T, Ahmedou Salem MV., Feis ME, Dekinga A, Holthuijsen S, et al. Resource landscapes explain contrasting patterns of aggregation and site fidelity by red knots at two wintering sites. Mov Ecol. 2018;6: 24. pmid:30598823
31. Alves JA, Gunnarsson TG, Potts PM, Sutherland WJ, Gill JA. Sex-biases in distribution and resource use at different spatial scales in a migratory shorebird. Ecol Evol. 2013;3: 1079–1090. pmid:23610645
32. Catry T, Lourenço PM, Lopes RJ, Carneiro C, Alves JA, Costa J, et al. Structure and functioning of intertidal food webs along an avian flyway: A comparative approach using stable isotopes. Funct Ecol. 2016; 468–478.
33. Ricciardi A, Bourget E. Global patterns of macroinvertebrate biomass in marine intertidal communities. Mar Ecol Prog Ser. 1999;185: 21–35.
34. Mandal S, Harkantra SN. Changes in the soft-bottom macrobenthic diversity and community structure from the ports of Mumbai, India. Environ Monit Assess. 2013;185: 653–672. pmid:22396068
35. Pennober G. Analyse spatiale de l’environnement côtier de l’archipel des Bijagos (Guinée-Bissau). Université de Bretagne occidentale—Brest. 1999.
36. Campredon P, Catry P. Bijagos Archipelago (Guinea-Bissau). In: Finlayson CM, Milton GR, Prentice RC, Davidson NC, editors. The Wetland Book II: Distribution, Description, and Conservation. Dordrecht: Springer; 2017. pp. 1333–1340.
37. Zwarts L. Numbers and distribution of coastal waders in Guinea-Bissau. Ardea. 1988;76: 42–55.
38. Delany S, Scott D, Dodman T, Stroud D. An Atlas of wader populations in Africa and Western Eurasia. Wageningen, The Netherlands: Wetlands International and International Wader Study Group; 2009.
39. Henriques M, Belo JR, Sá J, Monteiro H, Alves JA, Piersma T, et al. The Bijagós Archipelago: a key area for waterbirds of the East Atlantic Flyway. In: Roomen M van, Citegetse G, Crowe O, Dodman T, Hagemeijer W, Meise K, et al., editors. East Atlantic Flyway Assessment 2020: The status of coastal waterbird populations and their sites. Wilhelmshaven, Germany; Wageningen, The Netherlands; Cambridge, United Kingdom.: Wadden Sea Flyway Initiative, Wetlands International, BirdLife International; 2022. pp. 81–93.
40. Campredon P, Cuq F. Artisanal fishing and coastal conservation in West Africa. J Coast Conserv. 2001;7: 91–100.
41. Correia E, Granadeiro J, Regalla A, Dias E, Almeida A, Catry P. Predatory pelagic fishes of the Bijagós Archipelago (Guinea-Bissau) show high overlap in diets dominated by sardinella. African J Mar Sci. 2017;39: 389–396.
42. Correia E, Granadeiro JP, Regalla A, Catry P. Coastal waters of a marine protected area of the Bijagós Archipelago, West Africa, shelter juvenile fishes of economic and ecological importance. Reg Stud Mar Sci. 2021;46: 101892.
43. Lourenço P, Catry T, Lopes R, Piersma T, Granadeiro J. Invisible trophic links? Quantifying the importance of non-standard food sources for key intertidal avian predators in the Eastern Atlantic. Mar Ecol Prog Ser. 2017;563: 219–232.
44. Henriques M, Granadeiro JP, Piersma T, Leão S, Pontes S, Catry T. Assessing the contribution of mangrove carbon and of other basal sources to intertidal flats adjacent to one of the largest West African mangrove forests. Mar Environ Res. 2021;169: 105331. pmid:33878552
45. Meijer KJ, El-Hacen EHM, Govers LL, Lavaleye M, Piersma T, Olff H. Mangrove-mudflat connectivity shapes benthic communities in a tropical intertidal system. Ecol Indic. 2021;130: 108030.
46. Lourenço PM, Granadeiro JP, Catry T. Low macroinvertebrate biomass suggests limited food availability for shorebird communities in intertidal areas of the Bijagós archipelago (Guinea-Bissau). Hydrobiologia. 2018;816: 197–212.
47. Salvig JC, Asbirk S, Kjeldsen JP, Rasmussen PAF. Wintering waders in the Bijagos Archipelago, Guinee-Bissau 1992–1993. Ardea. 1994;82: 137–142.
48. Hill NK, Woodworth BK, Phinn SR, Murray NJ, Fuller RA. Global protected-area coverage and human pressure on tidal flats. Conserv Biol. 2021;35: 933–943. pmid:32969049
49. Henriques M, Catry T, Belo JR, Piersma T, Pontes S, Granadeiro JP. Combining multispectral and radar imagery with machine learning techniques to map intertidal habitats for migratory shorebirds. Remote Sens. 2022;14: 3260.
50. Granadeiro JP, Belo J, Henriques M, Catalão J, Catry T. Using Sentinel-2 images to estimate topography, tidal-stage lags and exposure periods over large intertidal areas. Remote Sens. 2021;13.
51. Temudo MP, Cabral AI. The social dynamics of mangrove forests in Guinea-Bissau, West Africa. Hum Ecol. 2017;45: 307–320.
52. Cardoso P. Os mangais da Guiné-Bissau: análise a 40 anos de evolução da sua extensão. Departamento de Biologia Animal Faculdade de Ciências da Universidade de Lisboa, Lisbon. 2017.
53. Day JH. A monograph on the polychaeta of South Africa. Trust Br MUSEUM (NATURAL Hist LONDON. 1967.
54. Cossignani R, T. A&. West African seashells (including Azores, Madeira and Canary. Ancona, Italia: L’Infonnatore Piceno; 2004.
55. Muñoz I, García-Isarch E, Sobrino I, Burgos C, Funny R, González-Porto M. Distribution, abundance and assemblages of decapod crustaceans in waters off Guinea-Bissau (north-west Africa). J Mar Biol Assoc United Kingdom. 2012;92: 475–494.
56. von Cosel R, Gofas S. Marine bivalves of tropical West Africa. Paris, Marseille, France: Muséum National d’Histoire Naturelle de Paris; IRD Marseille; 2020.
57. Wolff WJ, Duiven AG, Duiven P, Esselink P, Gueye A, Meijboom A, et al. Biomass of macrobenthic tidal flat fauna of the Banc d’Arguin, Mauritania. Hydrobiologia. 1993;258: 151–163.
58. Zwarts L, Wanink JH. How the food supply harvestable by waders in the Wadden Sea depends on the variation in energy density, body weight, biomass, burying depth and behaviour of tidal-flat invertebrates. Netherlands J Sea Res. 1993;31: 441–476.
59. Carneiro C, Gunnarsson TG, Méndez V, Soares AMVM, Alves JA. Linking range wide energetic tradeoffs to breeding performance in a long-distance migrant. Ecography (Cop). 2021;44: 512–524.
60. R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2018.
61. Imoobe TOT. Variation in benthic macroinvertebrate assemblages in Ologe Lagoon, Nigeria. African J Aquat Sci. 2008;33: 45–50.
62. Kon K, Kurokura H, Tongnunui P. Influence of a microhabitat on the structuring of the benthic macrofaunal community in a mangrove forest. Hydrobiologia. 2011;671: 205–216.
63. Salvig JC, Asbirk S, Kjeldsen JP, Palle AF, Quade A, Frikke J, et al. Coastal waders in Guinea-Bissau—aerial survey results and. 1997; 33–38.
64. Zwarts L, Blomert A-M, Hupkes R. Increase of feeding time in waders preparing for spring migration from the Banc d’Arguin, Mauritania. Ardea. 1990;78: 237–256.
65. Martins RC, Catry T, Santos CD, Palmeirim JM, Granadeiro JP. Seasonal variations in the diet and foraging behaviour of dunlins Calidris alpina in a South European estuary: Improved feeding conditions for northward migrants. PLoS One. 2013;8. pmid:24312532
66. Zwarts L, Dirksen S. Digestive bottleneck limits the increase in food intake of Whimbrels preparing for spring migration from the Banc d’Arguin Mauritania. Ardea. 1990;78: 257–278.
67. Burnham KP, Anderson DR. Multimodel inference: Understanding AIC and BIC in model selection. Sociological Methods and Research. 2004. pp. 261–304.
68. Clarke KR, Somerfield PJ, Chapman MG. On resemblance measures for ecological studies, including taxonomic dissimilarities and a zero-adjusted Bray-Curtis coefficient for denuded assemblages. J Exp Mar Bio Ecol. 2006;330: 55–80.
69. Anderson MJ, Walsh DCI. PERMANOVA, ANOSIM, and the Mantel test in the face of heterogeneous dispersions: What null hypothesis are you testing? Ecol Monogr. 2013;83: 557–574.
70. Kindt R, Coe R. Tree diversity analysis. A manual and software for common statistical methods for ecological and biodiversity studies. Nairobi: World Agroforestry Centre (ICRAF); 2005.
71. Anderson MJ, Willis TJ. Canonical analysis of principal coordinates: A useful method of constrained ordination for ecology. Ecology. 2003;84: 511–525.
72. Zwarts L. The winter exploitation of Fiddler Crabs Uca tangeri by waders in Guinea-Bissau. Ardea. 1985;73: 3–12.
73. Beukema J. Biomass and species richness of the macro-benthic animals living on the tidal flats of the Dutch Wadden Sea. Netherlands J Sea Res. 1976;10: 236–261.
74. Chainho P, Silva G, Lane MF, Costa JL, Pereira T, Azeda C, et al. Long-term trends inintertidal and subtidal benthic communities in response to water quality improvement measures. Estuaries and Coasts. 2010;33: 1314–1326.
75. Rodrigues AM, Meireles S, Pereira T, Gama A, Quintino V. Spatial patterns of benthic macroinvertebrates in intertidal areas of a Southern European estuary: The Tagus, Portugal. Hydrobiologia. 2006;555: 99–113.
76. França S, Vinagre C, Pardal MA, Cabral HN. Spatial and temporal patterns of benthic invertebrates in the Tagus estuary, Portugal: Comparison between subtidal and an intertidal mudflat. Sci Mar. 2009;73: 307–318.
77. Purwoko A, Wolff WJ. Low biomass of macrobenthic fauna at a tropical mudflat: An effect of latitude? Estuar Coast Shelf Sci. 2008;76: 869–875.
78. Wolff WJ, Michaelis H. Do shorebirds graze down zoobenthic biomass at the Banc d’Arguin tidal flats in Mauritania? Estuar Coast Shelf Sci. 2008;79: 491–495.
79. Swennen C. Oystercatchers (Haematopus ostralegus) feeding on giant bloody cockles (Anadara senilis) on the Banc d’Arguin in Mauritania. Ardea. 1990;78: 53–62.
80. Yoloye V. The ecology of the West African bloody cockle, Anadara (Senilia) senilis. BulI inst Français d’Afrique Noire, Dakar. 1976;38: 25–56.
81. Carneiro C, Gunnarsson TG, Alves JA. Bloody cockles: A novel and important food item for whimbrels in the banc d’arguin. Wader Study. 2017;124: 161–162.
82. Maretti CC. The Bijagós Islands; culture, resistance and conservation. IUCN-CEESP. 2015.
83. Wijnsma G, Wolff WJ, Meijboom A, Duiven P, De Vlas J. Species richness and distribution of benthic tidal flat fauna of the Banc d’Arguin, Mauritania. Oceanol Acta. 1999;22: 233–243.
84. Douglas EJ, Pilditch CA, Kraan C, Schipper LA, Lohrer AM, Thrush SF. Macrofaunal functional diversity provides resilience to nutrient enrichment in coastal sediments. Ecosystems. 2017;20: 1324–1336.
85. Greenfield B, Kraan C, Pilditch C, Thrush S. Mapping functional groups can provide insight into ecosystem functioning and potential resilience of intertidal sandflats. Mar Ecol Prog Ser. 2016;548: 1–10.
86. Ens BJ, Klaassen M, Zwarts L. Flocking and feeding in the fiddler crab (Uca tangeri): Prey availability as risk-taking behaviour. Netherlands J Sea Res. 1993;31: 477–494.
87. Silva RF, Filho JSR, Souza SR. Spatial and temporal changes in the structure of soft-bottom benthic communities in an Amazon estuary Caeté estuary, Brazil. J Coast Res. 2011;SI 64: 440–444.
88. Zabbey N, Hart AI. Spatial variability of macrozoobenthic diversity on tidal flats of the Niger Delta, Nigeria: the role of substratum. African J Aquat Sci. 2014;39: 67–76.
89. Choi CY, Battley PF, Potter MA, Ma Z, Liu W. Factors affecting the distribution patterns of benthic invertebrates at a major shorebird staging site in the Yellow Sea, China. Wetlands. 2014;34: 1085–1096.
90. Silva G, Costa JL, de Almeida PR, Costa MJ. Structure and dynamics of a benthic invertebrate community in an intertidal area of the Tagus estuary, western Portugal: a six year data series. Marine Biodiversity. Dordrecht: Springer Netherlands; 2006. pp. 115–128. https://doi.org/10.1007/1-4020-4697-9_10
91. Lourenço PM, Catry T, Piersma T, Granadeiro JP. Comparative feeding ecology of shorebirds wintering at Banc d’Arguin, Mauritania. Estuaries and Coasts. 2016;39: 855–865.
92. Gning N, Vidy G, Thiaw OT. Feeding ecology and ontogenic diet shifts of juvenile fish species in an inverse estuary: The Sine-Saloum, Senegal. Estuar Coast Shelf Sci. 2008;76: 395–403.
93. Piersma T. What is habitat quality? Dissecting a research portfolio on shorebirds. In: Fuller R, editor. Birds and Habitat: Relationships in Changing Landscapes. 2012. pp. 383–407. https://doi.org/10.1017/cbo9781139021654.019
94. Folmer EO, Olff H, Piersma T. How well do food distributions predict spatial distributions of shorebirds with different degrees of self-organization? J Anim Ecol. 2010;79: 747–756. pmid:20337758
95. Gill JA, Sutherland WJ, Norris K. Depletion models can predict shorebird distribution at different spatial scales. Proc R Soc B Biol Sci. 2001;268: 369–376. pmid:11270433
96. Alongi DM. The ecology of tropical soft-bottom benthic ecosystems. Oceanogr Mar Biol Annu Rev. 1990;28: 381–496.
97. Ysebaert T, Herman PMJ, Meire P, Craeymeersch J, Verbeek H, Heip CHR. Large-scale spatial patterns in estuaries: Estuarine macrobenthic communities in the Schelde estuary, NW Europe. Estuar Coast Shelf Sci. 2003;57: 335–355.
98. Baron J, Clavier J. Structure and seasonal variations in a sandy beach macrofauna on the South-West coast of New Caledonia (SW Pacific Ocean). Asian Mar Biol. 1994;11: 1–7.
99. Mucha AP, Vasconcelos MT, Bordalo A. Spatial and seasonal variations of the macrobenthic community and metal contamination in the douro estuary (Portugal). Mar Environ Res. 2005;60: 531–550. pmid:15919109
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2022 Coelho et al. This is an open access article distributed under the terms of the Creative Commons Attribution License: http://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
The coastal intertidal ecosystem of the Bijagós Archipelago, Guinea-Bissau, one of the largest and most important in West Africa, sustains a considerable proportion of the migratory shorebird populations of the East Atlantic Flyway and operates as a nursery area for benthic fish in the region. The macrozoobenthos in these mudflats constitute the main food source for both groups so that spatial and temporal variation in their abundance and community composition is likely to influence the abundance and distribution of fish and birds. In this study we described the spatial and temporal dynamics in the density, biomass, and community composition of macrozoobenthos across six intertidal flats in three islands of the Bijagós Archipelago. Overall, the Bijagós Archipelago was characterised by a highly species-rich macrozoobenthic community, with ca. 88 taxa identified across all sites, reaching a mean density of 1871 ± 58.3 ind.m-2 (mean ± SE) and mean biomass of 5.65 ± 0.41 g of AFDM.m-2 (ash-free dry mass per m2), values much lower than what was described for nearby intertidal areas, namely the Band d’Arguin, Mauritania. Density and biomass of major macrozoobenthos classes (Bivalvia, Polychaeta, Malacostraca and Gastropoda) differed across sites and months, displaying an overall increase in density towards the final months of the dry season (March and April). Similarly, community composition also differed significantly between sites and throughout the season. The site with most distinct community composition (Adonga) supported low diversity and high abundance of a few bivalve species, whilst other two sites that hosted the most diverse communities, were also the most similar between them (Anrumai and Abu). These spatial and temporal patterns constitute an important baseline to improve knowledge of this intertidal ecosystem and will contribute towards a better understanding of the spatial and temporal distribution patterns of their consumers.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details
; Paulino, João
; Loran Kleine Schaars; Ramos, Catarina; Aissa Regalla de Barros; Catry, Teresa; José Pedro Granadeiro
; Piersma, Theunis; Alves, José Augusto