Sediment grain size determines microplastic exposure landscapes for sandy beach macroinfauna.

Despite the global occurrence of microplastic contamination on sandy beaches, evidence of microplastic distribution within beaches remains contradictory. When conflicting evidence is used to inform sampling surveys, it increases uncertainty in resulting data. Moreover, it hampers spatially explicit risk characterization of microplastic pollution to intertidal fauna. We aimed to guide sampling designs for microplastic monitoring on beaches, and to quantify macroinfauna exposure to microplastics. Microplastic abundance, quantified between 5 mm-66 µm, lacked a significant zonation across the top sediment layer of sub-terrestrial, upper and lower midlittoral, and swash zones at two sites with varying anthropogenic influence on a microtidal dissipative beach in Uruguay. Microplastic abundance decreased exponentially with increasing grain size, as revealed by Bayesian Poisson regression, although the decrease was less steep compared to prior knowledge regarding sediment - plastic interactions obtained for large (millimeter-sized) industrial pellets. Significant differences in microplastic contamination between the two sites with varying anthropogenic influence likely related to their proximity to a freshwater canal. Corresponding field measurements of body burdens of fibers and irregular particles were significantly lower for the polychaete Euzonus (Thoracophelia) furcifera, despite its preference for finer sediments with higher microplastic loads, compared to the isopods Excirolana braziliensis and Excirolana armata. Results provide critical insights toward representative sampling of microplastics within beach sites. Specifically, we caution against sampling limited to the drift line, and instead recommend: 1) reporting beach morphodynamic characteristics; 2) using clearly defined, ecologically-informed zonation schemes; and 3) accounting for sediment grain size as a covariate to normalize among reported contamination levels. The results contribute valuable baseline data toward realistic exposure landscapes relative to the sediment grain size preferences of macroinfauna, needed to inform laboratory experiments.

Consequences of a freshwater canal discharge on the benthic community and its habitat on an exposed sandy beach.

This paper analyses the spatial and temporal effects of a freshwater discharge (Canal Andreoni) on the macroinfauna community and its habitat in a sandy beach of Uruguay. Bimonthly, we examined 17 environmental variables plus macroinfauna abundance, biomass, richness, evenness and diversity of three sites: Andreoni, at the canal mouth, Coronilla, at 1 km, and Barra, at 13 km from the mouth. Both univariate and multivariate techniques showed an increasing degree of perturbation towards the canal. This was reflected by abiotic and biotic differences between sites and by a consistent two-dimensional ordination of the samples. A clear seasonal pattern was found, specially accentuated at Coronilla, where the effects of the canal were stronger in winter and weaker in summer. Multivariate linking between macroinfauna and its habitat highlighted the role of salinity as explanatory variable of the observed trends.

Modified whole effluent toxicity test to assess and decouple wastewater effects from environmental gradients.

Environmental gradients and wastewater discharges produce aggregated effects on marine populations, obscuring the detection of human impact. Classical assessment methods do not include environmental effects in toxicity tests designs, which could lead to incorrect conclusions. We proposed a modified Whole Effluent Toxicity test (mWET) that includes environmental gradients in addition to effluent dilutions, together with the application of Generalized Linear Mixed Models (GLMM) to assess and decouple those effects. We tested this approach, analyzing the lethal effects of wastewater on a marine sandy beach bivalve affected by an artificial canal freshwater discharge used for rice crops irrigation. To this end, we compared bivalve mortality between canal water dilutions (CWd) and salinity controls (SC: without canal water). CWd were prepared by diluting the water effluent (sampled during the pesticide application period) with artificial marine water. The salinity gradient was included in the design by achieving the same final salinities in both CWd and SC, allowing us to account for the effects of salinity by including this variable as a random factor in the GLMM. Our approach detected significantly higher mortalities in CWd, indicating potential toxic effects of the effluent discharge. mWET represents an improvement over the internationally standardized WET tests, since it considers environmental variability and uses appropriate statistical analyses.

Disentangling diversity patterns in sandy beaches along environmental gradients.

Species richness in sandy beaches is strongly affected by concurrent variations in morphodynamics and salinity. However, as in other ecosystems, different groups of species may exhibit contrasting patterns in response to these environmental variables, which would be obscured if only aggregate richness is considered. Deconstructing biodiversity, i.e. considering richness patterns separately for different groups of species according to their taxonomic affiliation, dispersal mode or mobility, could provide a more complete understanding about factors that drive species richness patterns. This study analyzed macroscale variations in species richness at 16 Uruguayan sandy beaches with different morphodynamics, distributed along the estuarine gradient generated by the Rio de la Plata over a 2 year period. Species richness estimates were deconstructed to discriminate among taxonomic groups, supralittoral and intertidal forms, and groups with different feeding habits and development modes. Species richness was lowest at intermediate salinities, increasing towards oceanic and inner estuarine conditions, mainly following the patterns shown for intertidal forms. Moreover, there was a differential tolerance to salinity changes according to the habitat occupied and development mode, which determines the degree of sensitivity of faunal groups to osmotic stress. Generalized (additive and linear) mixed models showed a clear increase of species richness towards dissipative beaches. All taxonomic categories exhibited the same trend, even though responses to grain size and beach slope were less marked for crustaceans and insects than for molluscs or polychaetes. However, supralittoral crustaceans exhibited the opposite trend. Feeding groups decreased from dissipative to reflective systems, deposit feeders being virtually absent in the latter. This deconstructive approach highlights the relevance of life history strategies in structuring communities, highlighting the relative importance that salinity and morphodynamic gradients have on macroscale diversity patterns in sandy beaches.