Risk Characterization of Microplastics in San Francisco Bay, California

doi:10.21203/rs.3.rs-1507100/v1

Download PDF

Research Article

Risk Characterization of Microplastics in San Francisco Bay, California

https://doi.org/10.21203/rs.3.rs-1507100/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Assessing microplastics risk to aquatic ecosystems has been limited by lack of holistic exposure data and poor understanding of biological response thresholds. Here we take advantage of two recent advances, a toxicological meta-analysis that produced biotic response thresholds and a method to quantitatively correct exposure data for sampling methodology biases, to assess microplastic exposure risk in San Francisco Bay, California, USA. Using compartment-specific particle size abundance data, we rescaled empirical surface water monitoring data obtained from manta trawls (> 333 µm) to a broader size (1 to 5,000 µm) range, corrected for biases in fiber undercounting and spectroscopic subsampling, and assessed the introduced uncertainty using probabilistic methods. We then compared these rescaled concentrations to four risk thresholds developed to inform risk management for California for each of two effect categories/mechanisms - tissue translocation-mediated effects and food dilution - each aligned to ecologically relevant dose metrics of surface area and volume, respectively. More than three-quarters of samples exceeded the most conservative food dilution threshold, which rose to 85% when considering just the Central Bay. Within the Central Bay, 38% of the samples exceeded a higher threshold associated with management planning, which was statistically significant at the 95% confidence interval. For tissue translocation-mediated effects, no samples exceeded any threshold with statistical significance. The risk associated with food dilution is higher than that found in other systems, which likely reflects this study having been conducted for an enclosed water body. A sensitivity analysis indicated that the largest contributor to assessment variability was associated with estimation of ambient concentration exposure due to correcting for fiber undercounting. Even after compensating for biases associated with fibers and other small particles, concentrations from the trawl samples were still significantly lower than the 1-L grab samples taken at the same time, suggesting our SFB risk estimates are an underestimate. We chose to rely on the trawl data because the 1-L grab sample volume was too small to provide accurate spatial representation, but future risk characterization studies would be improved by using in-line filtration pumps that sample larger volumes while capturing a fuller range of particle size than a towed net.

Microplastics

environmental risk assessment

probability density functions

marine ecosystems

estuarine ecosystems

management framework

Microplastics have been found in a wide array of aquatic environments, from pristine mountain streams to the Arctic (Gonzalez-Pleiter et al. 2020) to deep undersea habitats (Abel et al. 2021). Toxicological studies have determined microplastics can cause adverse effects, such as tissue inflammation (Pirsaheb et al., 2020), impaired growth (Zimmerman et al. 2020, Redondo-Hasselerharm et al. 2018), feeding disruption (Wang et al. 2020), developmental anomalies (Gardon et al. 2020), and changes in gene expression (Zhang et al. 2019). However, the prevalence of those biotic effects in natural aquatic ecosystems is not well understood (Gouin et al. 2019).

Quantifying the risk of microplastics in aquatic ecosystems is challenging for two reasons. First, the concentrations at which those effects manifest in biota are not well understood. That uncertainty arises because of shortcomings in existing toxicological studies (de Ruijter et al. 2020), with fewer than half of the studies conducted to date having included more than two exposure concentrations and many of those exposures at extreme concentrations well beyond what is typically encountered in the natural environment (Burns and Boxall 2018, Hampton et al. 2022). Although this testing provides useful insights into potential effects and mechanisms of toxicity, testing at multiple relevant concentrations to establish a dose-response relationship is necessary to quantitatively characterize risk. Exacerbating this problem is that microplastics have a diversity of properties, such as size, shape, and polymer type, that can affect toxicity, and few studies have quantified the relative importance of these factors (Hampton et al.,2022; Bucci and Rochman, 2022). Most (72% out of 163) toxicity studies have been conducted using single-sized beads of a single polymer type (Mehinto et al. 2022) which is a poor representation of mixtures encountered in the ambient environment (Rubin et al., 2021).

The second challenge is a lack of holistic exposure data to compare directly to toxicologically derived response thresholds. Most ambient microplastics data are collected by towing ~ 330 µm mesh nets, which underestimates the abundance of microplastics smaller than the mesh size (Brander et al., 2020). Studies that have sampled from the environment and report broader size distributions find that the smaller sized particles are exponentially more abundant (Covernton et al., 2019; Kooi et al. 2021), suggesting the need for sampling regimes and/or estimation methods that capture a more complete size range of particles. Additionally, field monitoring particle data often suffers from unquantified biases due to self-contamination (Scopetani et al. 2020), difficulties associated with sampling and analyzing fibers (Miller et al. 2021), spectroscopic interferences and library mischaracterizations (Cowger et al. 2020; Werbowski et al. 2021), spatial and temporal heterogeneity (Koelmans et al. 2019), as well as spectroscopy subsampling regimes performed to ensure feasibility when particles are counted manually (Zhu et al. 2021).

Here we take advantage of two recent advances that address these challenges. The first is a meta-analysis in which a broad array of toxicological studies were incorporated into a single risk assessment framework (Mehinto et al. 2022), which produced thresholds for a range of biotic responses and recommended management actions. This meta-analysis applied critical quality criteria to screen reliable toxicity studies and integrated the results into a combined outcome that transcended shortcomings of the underlying individual studies. The second is the use of probability density functions (PDFs) to quantitatively correct exposure data for biases due to sampling methodologies (Koelmans et al., 2020, Kooi and Koelmans 2019). Size abundance microplastic particle data can be used to derive probability density functions (PDFs) that allow the rescaling of field monitoring data restricted to a given size range (e.g., > 330 µm) to a more holistic size range (e.g., 1 to 5,000 µm), enabling direct comparison to toxicity thresholds from laboratory studies aligned to the same size range (Koelmans et al 2022).

Combining these two advances, we assess the risk to aquatic ecosystems from microplastic exposure in San Francisco Bay, California (SFB) where a comprehensive study of ambient exposure was conducted (Zhu et al 2021) (Fig. 1). After rescaling to a common size distribution, we compare the monitoring data to aligned risk and management thresholds from Mehinto et al. (2022) to estimate the likelihood and pervasiveness of a local biological response. Using PDFs and Monte-Carlo modeling, we quantify the uncertainty of the rescaling methods to determine where the greatest uncertainties in this risk characterization lie, pointing to the science advancements needed to improve risk assessments in the future.

Figure 1. Flowchart of general steps involved in microplastics risk characterization as employed in this study. Data obtained/derived from respective studies are annotated using colors; Zhu et al. (2021) is green, Mehinto et al.(2022) is red, and this study is blue.

Data quality

Crucial to assessing risks is selecting data fit for that assessment purpose. Surface water monitoring data for microplastics in SFB, California, USA reported by Zhu et al. (2021) were quantitatively assessed for quality according to the criteria defined for water sampling in Koelmans et al. (2019). Data reported for other matrices (e.g., stormwater effluent, fish tissue, sediment) were not scored due to a lack of established quality criteria for such matrices. Briefly, nine criteria which relate to the reproducibility and reliability of aqueous microplastics sampling, contamination mitigation, sample processing/handling, and chemical analysis were applied. For each criterion, a score of 0, 1, or 2 was applied and a total accumulated score was calculated by adding scores for individual criteria (maximum 18 points). Samples that received a ‘zero’ value for any individual score were not considered sufficiently reliable (Hermsen et al. 2018).

Microplastics characterization

Blank-corrected environmental microdebris occurrence data from SFB and outside of SFB, in National Marine Sanctuaries, were obtained from Zhu et al. (2021). Sampling details can be found in Zhu et al., 2021. In short, sampling sites were picked to represent each region within SFB. Regions are characterized by differences in population sizes and point sources upstream, such as wastewater and stormwater. Surface water concentrations were spectroscopy corrected. Spectroscopy was only performed on a subset of particles, i.e., roughly 10% of each morphology within each size fraction (see Zhu et al., 2021 for more detail). This methodology was chosen in an attempt to be more representative. This led to spectroscopy conducted on 23% of all particles from surface manta trawls. The systematic removal of all fiber particle counts from manta trawl data in Zhu et al (2021) was corrected for using a subset of manta samples in which all fibers were counted and reported in Hung et al (2020). While Zhu et al. (2021) presented a novel method to correct manta trawl samples for size, a different size rescaling method was used here (i.e., Koelmans et al. 2020). Blank-, fiber-, and plastic polymer-corrected particle concentrations were rescaled to a common size distribution (1 to 5,000 µm) to allow direct comparison to hazard thresholds (Mehinto et al., 2022) according to the methods in Koelmans et al (2020; described below) using marine surface water size distribution data from Kooi et al. (2022). Additionally, a statistically significant outlier was identified based on four-times the mean Cook’s distance and was removed from the Zhu et al. (2021) dataset, which was a sample collected from a tidal front and was highly contaminated with microplastics and other debris (sample identification: CB9-Manta-11 Jan 18).

Due to time-constraints of spectroscopically confirming the polymer identity of all particles in samples, Zhu et al. (2021) subsampled particles from samples based on the number of particles of a particular morphology and size class within a given sample. For each site in the SFB and compartment (i.e., stormwater, wastewater, fish tissue, sediment, surface water) the proportion of particles that were spectroscopically determined to be a specific polymer (e.g., polyester, polyethylene, etc.) were divided by the total number of particles spectroscopically characterized for that compartment-site combination (Fig. 1). Zhu et. al. (2021) reported that interferences such as dyes and carbon black prevented the spectroscopic confirmation of all particles and reported some polymers using suspected terms such as “anthropogenic (synthetic)” or “anthropogenic (unknown base)”, etc. To be conservative, particles that could not be polymerically confirmed were excluded from the proportion of microdebris particles considered to be plastic, with a hierarchical schema developed here and employed to classify particle types (Figure S1). To determine if a single polymer correction factor should be applied to all matrices or should be performed separately for each matrix, a one-way analysis of variance (ANOVA) was conducted (Table S3). Then, for surface water data obtained using manta trawl, a one-way ANOVA was conducted to determine if there are site-specific differences in plastic proportions of microdebris particles, with site-specific correction factors being applied only if the ANOVA determined statistically significant differences. Proportions of microdebris confirmed to be plastic values were multiplied by the total number of microdebris particles reported by Zhu et al. (2021) for each compartment to obtain microplastic occurrence data in a probabilistic manner as part of the Monte-Carlo modeling method (described below).

Due to the mesh size of a manta trawl net (333 µm), fibers are known to pass through the sampling apparatus, leading to a systematic shape-based undercounting bias. Zhu et al (2021) did not include fibers in total particle counts in their blank-corrected concentrations for samples collected using manta trawls to remove that uncertainty with the impact of systematically undercounting microplastic particles. Hung et al (2020) reported particle count data for all shapes - including fibers - in 9 manta trawl samples taken from various sites across the SFB and Marine Sanctuaries. This manta trawl fiber subsampling dataset was used to derive a correction factor to estimate the amount of microplastics that would be present in the other manta trawl samples from SFB if plastic fibers had been counted in those samples, with uncertainties propagated probabilistically using the Monte-Carlo method described below. For each manta trawl sample with fiber counts, the proportion of particles that were fibers was calculated and a fiber correction factor was derived as the inverse of one minus the median fraction of particles that were fibers. Due to the small sample size, site-specific differences in fiber proportions were not considered.

Rescaling of environmental concentrations

Environmental microdebris occurrence data from SFB reported in Zhu et al. (2021) included various size ranges of particles based on each sampling technique (e.g., > 333µm for manta trawl; > 50 µm for grab samples: see Table S8) which were rescaled to a common size distribution of 1–5,000 µm to compare to ecotoxicity thresholds aligned to the same size distribution in Mehinto et al (2020) based on the methods described in Koelmans et al. (2020; 2022). Environmental concentrations were multiplied by a correction factor derived for each sampling technique based on their particle size limits (Eq. 1) (Koelmans et al. 2020).

Equation 1\({C}_{env}={CF}_{meas}* {C}_{meas}\)

In Eq. 1, C_env is the environmentally realistic occurrence concentration in particles · L^− 1 (adjusted for non-alignment of mesh sizes), CF_meas is a dimensionless correction factor for the environmentally monitored concentration (meas); and C_meas is the measured environmental concentration, expressed in particles · L^− 1 (Koelmans et al. 2020). Environmental concentrations are rescaled to an upper (UL,D; µm) and lower default size range (LL,D; µm) (here 5000 and 1 µm respectively), using the power law slope of microplastic particle abundance in the environment based on size (a, unitless), with the upper limit (UL, meas; µm) and lower limit (LL, meas; µm) defined by the size limits of quantification of the monitoring method employed (Eq. 2) (Koelmans et al. 2020).

Equation 2\({CF}_{meas}=\frac{{L}_{UL,D}^{1-a}-{L}_{LL,D}^{1-a}}{{L}_{UL, meas}^{1-a}-{L}_{LLmeas}^{1-a}}\)

Kooi et al. (2021) derived power slope exponents (a) based on size for freshwater and marine environments across several locations in Europe using individual particle datasets obtained using state-of-the-science Fourier-transform infrared imaging coupled with a focal-plane array detector and automated image analysis. When measured particles length data was pooled across all samples for each compartment across distinct locations (e.g., Rhine and Dommel rivers), power law exponent values contained low variability within each compartment, but were significantly different between compartments (e.g., a = 2.64 ± 0.01 for marine and a = 2.07 ± 0.03 for freshwater surface waters), implying microplastic size relationships are highly conserved within compartments (Kooi et al. 2022). While within-compartment particle size distributions are not expected to deviate significantly across regions, site-specific data would be preferable to rescale environmental concentrations so long as the data is high resolution and is reliable (Kooi et al. 2022). Lacking site-specific high-resolution particle size distribution data, compartment-specific a values for length were used from Kooi et al. (2022) to rescale environmental concentrations in SFB (Table S7). To account for within-compartment variability in rescaling environmental concentrations, uncertainties in the derived a value were propagated probabilistically using Monte-Carlo methods described below.

Zhu et al. (2021) reported particle length data for SFB for all compartments (except surface water obtained using 1-L grab samples) which was measured manually and was not intended to provide high-resolution information on particle size distributions. Nonetheless, this dataset was used to derive compartment-specific size a values according to the methods described in Kooi & Koelmans (2019) as part of a sensitivity analysis only and were not used to rescale concentrations for risk characterization purposes. Briefly, all particle monitoring data from SFB reported on an online open data repository managed by the California Natural Resources Agency (https://data.cnra.ca.gov/dataset/microplastic-sf-bay) were aggregated, and only the measured lengths of individual particles spectroscopically confirmed to be plastic were used in estimating size distributions. Plastic particle data were grouped by compartment (i.e., surface water, sediment, fish tissue, wastewater effluent, stormwater runoff), and abundance particle length-based data were plotted on a log-log scale with relative abundance on the y-axis and particle size on the x-axis and a linear trendline fit to the data was used to derive the a exponent value. Since power laws usually only apply to values greater than some minimum value (in this case, particle length), both the minimum applicable value and the final a value were determined using a maximum likelihood estimation method (Clauset et al., 2009; Newman, 2005) using the poweRlaw package (Gillespie, 2014), and bootstrap 100 times. As part of the sensitivity analysis only, the a value derived using manta trawl data (> 333 µm) was used to rescale the manta trawl-derived surface water monitoring data.

Risk Characterization

Aquatic ecotoxicological thresholds were used to characterize risk by deriving the ratio of predicted no-effect concentration (PNEC) thresholds to predicted environmental concentrations (PEC), with exceedances greater than one indicative of risk. While data were available for microplastics concentrations in various matrices (e.g., surface water, stormwater, wastewater effluent, sediment, fish tissue) in the SFB, only surface water concentration data were used for risk characterization as direct comparisons of undiluted stormwater or wastewater are not representative of environmental exposures, and due to the lack of hazard thresholds for marine sediment or fish tissue.

Surface water concentration data were compared to ecotoxicological thresholds derived by Mehinto et al.(2022) using species sensitivity distributions (SSDs) based on chronic no-observed-adverse-effect concentrations for 14–16 freshwater and marine species from 6 to 7 taxonomic groups. Prior to derivation, Mehinto et al. screened 162 peer-reviewed laboratory toxicity studies for a set of pre-defined quality criteria based on the standards defined by de Ruijter et al. (2020). A total of 290 threshold data points were extracted from 21 studies that met the minimum pre-defined criteria. These thresholds were aligned to a common size distribution of 1–5,000 µm using environmental PDFs and based on mechanisms of action as described in Koelmans et al. (2020) and Kooi et al. (2021). Following Kooi et al (2021), Mehinto et al (2022) derived thresholds for two effect mechanisms/pathways – food dilution and tissue translocation, which were aligned by volume and surface area ecologically relevant metrics, respectively. The food dilution-based effect considered particles small enough to be ingested by the organism of interest to be accessible (i.e., exclude non-accessible particles), then aligned (both monodisperse and polydisperse) laboratory effect concentrations to environmentally realistic concentrations based on particle volume (Mehinto et al.2022). The tissue translocation mechanism of action considered particles wide enough to translocate across tissues (83 µm) following ingestion to be accessible and aligned laboratory effect concentrations to environmentally realistic concentrations based on translocatable surface area (Mehinto et al.2022). Although the methodology used to align thresholds was identical for the mechanism of the food dilution effect, additional studies were used in Mehinto et al. (2022) compared to Koelmans et al. (2020), which only used effect thresholds data for studies in which the authors confirmed that a food dilution mechanism was demonstrated or was plausible. For both effect mechanisms/pathways, four PNEC thresholds were derived which correspond to different levels of confidence that microplastics can cause adverse effects to aquatic organisms and call for varying levels of management action - ranging from increasing monitoring to implementing Source control measures (Table S9) (Mehinto et al.2022). At the time of writing, the risk management actions associated with thresholds from Mehinto et al. (2022) carry no regulatory or legal authority in California or any other jurisdiction and are only suggestions.

Sensitivity Analysis

Uncertainties were evaluated probabilistically using Monte Carlo methods based on PDFs derived for each correction factor, including: manta trawl fiber correction (shape under-counting bias), plastic fraction of total microdebris particles (spectroscopic subsampling variability), and rescaling concentrations to a common size range (a variability) (Fig. 1). To obtain the combined correction factor with probabilistic propagation of uncertainties, a data frame of 10,000 values was generated for each correction factor based on their modeled distribution, and each of these three data frames was multiplied by one another row-wise. The 50th percentile value from this combined correction factor distribution was used to correct the manta trawl surface water monitoring data, and the 5th and 95th percentile values were used to calculate uncertainty. This methodology accurately accounts for the underlying distributions of the correction factor data and is preferable to error propagation techniques that rely on assumptions of normality.

For each correction factor, a theoretical distribution was fit to the data based on the shape of the underlying distribution (see Table S10) - which was evaluated visually using Cullen and Frey graphs (Figure S8) and using a maximum likelihood estimation approach with the fitdistrplus R package (Delignette-Muller & Dutang 2015). Due to the highly skewed distribution in the manta trawl fiber correction data (Figure S6), these data were log-transformed prior to derivation of a PDF, then back-transformed into linear space following the Monte Carlo simulation. Following log-transformation, fiber correction data (unitless; >1) were well-fit by a normal distribution (Figure S8a). Plastic proportion data (unitless; 0–1) were well-fit by a two-shape beta distribution (Figure S8b) that was truncated such that values greater than one was not produced during the Monte Carlo simulated due to their theoretical implausibility (i.e., more than 100% of particles cannot be plastic). Because the size-based a values used to rescale concentrations (i.e., from Kooi et al. 2022) were derived using maximum likelihood estimation based on a log-log linear distribution (Figure S7; Table SS7), a normal distribution was assumed, and the PDF was approximated using a normal distribution based on the mean and standard deviation of the linear regression slope (Table S7). Correction factors were derived from the Monte-Carlo simulated distribution of alpha values using Eq. 2.

To quantify and compare the relative sensitivity of correction and rescaling factors applied here on the resulting exposure assessment, variability for each parameter was held constant while uncertainty in the other variables was calculated. Finally, to assess the relative uncertainty between rescaled and corrected environmental occurrences with modeled risk thresholds, the 95th percentile of the Monte Carlo simulated occurrence data water was compared to the 95% confidence intervals for microplastics hazard thresholds reported in Mehinto et al.(2022) based on the SSD model. An additional sensitivity analysis was performed using site-specific a values derived for SFB using manta trawl particle length data as described above.

Statistics

All statistical analyses were conducted in R (version 4.1.1; R Core Team, 2021) and figures were produced using the package ggplot2 (Wickham, 2009). Base maps sourced from Google were used for mapping using the ggmap package (Kahle and Wickham, 2013). One-way ANOVAs were used to determine if plastic correction factors should be separated for each matrix and site (for manta trawl only). To determine if bias corrections (i.e., rescaling, fiber correction, plastic correction) resulted in comparable concentrations between water matrices, one-way ANOVAs were run for both raw and rescaled/corrected concentrations. For all hypothesis tests, statistical significance was determined at an alpha level of 0.05, and multiple comparisons were performed using Tukey’s Honest Significant Difference post-hoc test, when applicable. All Monte-Carlo simulations were performed with 10,000 iterations with a seed set for reproducibility.

Data Quality

Microplastics monitoring data reported in Zhu et al. (2021) received total accumulated scores of 13, 10, and 14 for manta trawl, grab samples, and wastewater treatment plant effluent samples respectively according to criteria defined in Koelmans et al. (2019) (Table S1). While manta trawl and wastewater treatment plant effluent data received a score of at least one for each quality criteria, grab samples received “zero” scores for several criteria (sample size and sample treatment) (Table S1). Accordingly, manta trawl and wastewater treatment plant effluent data from Zhu et al. (2021) are considered sufficiently reliable for the purposes of risk characterization, while the grab sample data are not, however only manta trawl data (surface water) were used for risk characterization due to the non-applicability of wastewater data for estimating exposure. Because blank-corrections were applied based on color-morphology combinations instead of polymer identification, all matrices received a score of “1” for negative controls instead of “2”. The blank correction procedure applied in Zhu et al. (2021) may lead to an underestimation of concentrations if microplastic particles of the same color have a different polymer identity - an uncertainty which is not accounted for in this probabilistic assessment. The grab, wastewater treatment plant effluent, and manta trawl data scores were higher than the average score for surface waters (7.9; range 4 to 15; n = 55) reported in Koelmans et al. (2019). While additional quality criteria are available for biota and sediment (Hermsen et al. 2018; Redondo-Hasselerharm et al. in press; Bäuerlein et al. Submitted), risk thresholds are unavailable for these compartments and these data were not quality scored here. Additionally, stormwater was not scored due to a lack of established quality criteria.

Microplastics characterization

The percentage of analyzed particles spectroscopically determined to be plastic was significantly different between matrices according to a one-way ANOVA (p < 1x10^− 16; Table S3), so matrix-specific plastic correction factors were derived accordingly using PDFs (Figures S3 and S9; Tables S2 - S4; Table S10). Tukey’s post-hoc test for significance revealed significant differences in plastic proportions of total particles between manta trawl and sediment, fish tissue, wastewater treatment plant (WWTP) effluent, and surface water collected with 1-L grab (Table S4). Surface water samples collected with manta trawl contained the highest percentages of confirmed microplastics (72% ± 24%), followed by surface water collected by 1-L grab (42% ±24), sediment (37% ± 14%), wastewater treatment plant effluent (31% ± 18%), and fish tissue (24% ±12%) (Table S2; Figure S2). Additional significant differences were found between sediment and fish tissue; sediment and 1-L grab surface water; and fish tissue and 1-L grab surface water (Table S4). As surface water data obtained using manta trawl were the only monitoring data used for risk characterization here, site-specific differences were tested using a one-way ANOVA, which demonstrated no significant difference in proportions of plastic particles relative to all spectroscopically characterized particles by location (p = 0.12; Table S5; Figures S3 – S4). Accordingly, a single correction factor for plastic percentages was applied to all manta trawl data regardless of location, which was the median value of 0.63 (0.31 to 0.95: 95% CI) that was derived from a two-shape beta distribution PDF (Table S10; Figure S9).

To correct for the systematic removal of fiber particle data from the blank-corrected dataset reported in Zhu et al. (2021), data were used from Hung et al. (2020) for 9 manta trawl samples from SFB in which all fibers were counted. On average, fibers constituted 78% (± 28% sd) of particles in the manta trawl samples in which they were counted (Table S6). Other aqueous matrices and sediment contained lower percentages of fibers, while fish tissue contained a higher percentage (Table S6, Figure S5). Based on the PDF of the fiber proportions, a fiber correction factor of 8.87 (95% CI: 1.29 to 50.89; Table S10) was calculated and applied to manta trawl monitoring data as part of the Monte-Carlo analysis. Due to the relatively small sample size (n = 9) and skewed nature of the fiber proportion distribution (Figure S6), the fiber correction factor contains relatively high uncertainty compared to the plastic spectroscopy correction factor and size rescaling correction factor (Figure S13).

Rescaled Environmental Occurrence Data

Size abundance distributions of microplastics in SFB were fit by linear regression on a log(10)-log(10) scale using a maximum likelihood estimation approach (Kooi et al. 2021), with a exponent values ranging from 2.15 to 3.02 (Figure S7, Table S7). Length-based power law exponent values (a) derived for microplastics in SFB were comparable to values derived from various locations in Europe reported by Kooi et al. (2021) (Table S7). Notably, the a values for marine surface waters were 2.15 ± 0.48 and 2.07 ± 0.03 (mean ± sd) for SFB and in Europe, respectively, thus representing less than a 5% difference and are not statistically significant from one another (Table S7). Additionally, a law exponents followed the same rank order by matrix between Kooi et al. (2021) and those derived here for comparable matrices (i.e., marine surface water < wastewater effluent < marine sediment). The greatest difference between a law exponents was for marine sediment (2.90 ± 0.41 and 2.57 ± 0.20 for SFB and Europe, respectively) (Table S7). Direct comparisons for power law exponents derived for SFB to other studies/locations were not possible for stormwater runoff (which has not been reported elsewhere) or fish tissue (“biota” reported in Kooi et al. 2021 corresponds to benthic invertebrates).

Size-based correction factors for matrices ranged from 58 to 9,774 depending on matrix (a value) and mesh size (Table S8). Fish tissue data had the smallest mesh size (25 µm) and had the smallest correction factor accordingly (58; 95% CI: 53 to 63) (Table S8). Manta trawl data had the largest mesh size (333 µm) and had the second largest correction factor (529; 95% CI: 401 to 704 based an a value of 2.07 form Kooi et al. 2021) (Table S8). Stormwater data had a smaller mesh size than manta trawl (106 µm), however the correction factor was over 10x higher (9,774; 95% CI: 5,614 to 17,019) due to the high a value (2.97 ± 0.83) (Table S8).

In theory, rescaling data and correcting for systematic biases (i.e., fiber correction, spectroscopic subsampling) should reduce differences in monitoring concentrations taken at similar times and locations within a given matrix due to size-differences in mesh sizes of sampling apparatus (Koelmans et al. 2020). Before rescaling and correcting, surface water concentrations collected using manta trawl as well as effluent from wastewater treatment plants were not significantly different from one another according to one-way ANOVA with Tukey’s post-hoc (p > 0.05; Tables S11 – S12) but were both significantly lower than surface water collected through other means (stormwater, 1-L grab surface water; p < 0.001) (Table S12). Additionally, 1-L grab surface water concentrations were significantly higher than stormwater concentrations collected with a depth-integrated peristaltic pump (p = 0.03; Table S12). Following rescaling and correcting, manta trawl-collected surface water concentrations were still significantly lower than other surface water concentrations collected via other methods (i.e., stormwater, 1-L grab surface water) (p < 0.001; Tables S13 – S14), however wastewater concentrations were no longer significantly different from both 1-L grab and manta trawl-collected surface water concentrations (p > 0.05) (Fig. 2; Table S14). Despite the combined correction factor to account for systematic under-counting based on size and fibers as well as fractions of particle spectroscopically confirmed to be plastics, rescaled and corrected manta trawl surface water data were still significantly lower (p = 3.4 x 10^− 14; Table S14) than rescaled surface water 1-L grab samples, of which the majority were taken at similar locations and times. These results suggest that additional systematic biases are present in either the manta trawl (likely undercounting) or the 1-L grab samples (potentially overcounting). Undercounting in manta trawl samples may be due in part or in whole to imprecise blank corrections based on shape-color combinations as opposed to polymer-based corrections.

Risk characterization

Depending on the postulated effect mechanism/pathway (i.e., food-dilution or tissue translocation), risk exceedances of microplastics in SFB vary significantly. For all comparisons of PNECs (i.e., hazard thresholds from Mehinto et al. 2022) with PECs (i.e., corrected, and rescaled surface water concentrations in SFB) stated throughout this manuscript, only those in which the 95% CI does not include ‘0%’ represent statistically significant exceedances. Accordingly, only food-dilution thresholds one, two, and three have statistically significant exceedances in the SFB, while all other thresholds (i.e., food-dilution threshold four, and all tissue translocation thresholds) do not.

Comparison of corrected and rescaled manta-trawl collected surface water samples with food-dilution thresholds derived by Mehinto et al. (2022) resulted in 82% (95% CI: 27–100%) of samples exceeding the most conservative risk threshold (i.e. “Investigative monitoring” threshold one), 27% (95% CI: 3–73%) of samples exceeding threshold two (“Discharge monitoring”), 21% (95% CI: 3–58%) of samples exceeding threshold three (“Management planning”), and 3% (95% CI: 0–18%) of samples exceeding threshold four (“Source control measures”) (Table S15).

Comparison of surface water samples with tissue translocation-based thresholds derived by Mehinto et al. (2022) resulted in 3% (95% CI: 0–9%) of samples exceeding the most conservative risk threshold (i.e., “Investigative monitoring” threshold one), 0% (95% CI: 0 to 3%) of samples exceeding threshold two (“Discharge monitoring”), 0% (95% CI: 0 to 3%) of samples exceeding threshold three (“Management planning”), and 0% (95% CI: 0 to 0%) of samples exceeding threshold four (“Source control measures”) (Table S15).

Risk exceedances were higher during the rainy season, with 94% (95% CI: 41–100%) of surface water samples collected following a storm event exceeding food dilution threshold one compared with 71% (95% CI: 12–100%) of samples collected during the dry season (Figure S10). Rainy season samples exceeded food dilution threshold three within confidence limits (29%; 95% CI: 6–71%), however dry season samples did not (12%; 95% CI: 0–47%) (Figure S10).

Risk exceedances varied by location within the SFB (Fig. 4, Table S16). The Central Bay had the highest proportion of samples exceeding risk thresholds, with 85% (95% CI: 38–100%) exceeding Mehinto et al. (2022)’s most conservative food dilution threshold one (“Investigate monitoring”), 38% (95% CI: 8–85%) exceeding food-dilution threshold two (“Discharge monitoring”), 38% (95% CI: 8–77%) exceeding threshold three (“Management planning”), and 8% exceeding food dilution threshold four (“Source control measures”), however exceedances of threshold four were not statistically significant (95% CI: 0 to 31%) (Table S16). Additionally, the Central Bay was the only location with any samples exceeding a tissue translocation-based threshold at the 50th percentile, with 8% exceeding threshold one – however these exceedances were not statistically significant (95% CI: 0 to 23%) (Table S16).

Comparison of SFB samples to samples taken from outside of the bay demonstrated substantially higher risk within the bay. Samples taken from the National Marine Sanctuaries - which is an open-ocean location with minimal inputs from wastewater discharge or stormwater runoff and was selected as a reference location as part of the study design (Zhu et al. 2021) did not have any samples exceeding the most conservative threshold (i.e., food-dilution threshold one) with statistical significance (i.e., 35%; 95% CI: 0–91%) (Table S16; Fig. 4). Of the samples in the National Marine Sanctuaries exceeding food dilution threshold one at 50th percentile, the three highest were at the mouth of the bay just West of the Golden Gate Bridge, suggesting rapid dilution of microplastic particle concentrations outside of the SFB (Fig. 4).

Sensitivity Analysis

Comparison of the influence of factors in estimating environmental occurrence from manta trawl data reveals that the fiber correction factor contains the highest relative uncertainty compared to the spectroscopic sub-sampling correction factor for plastics and the size-based alignment correction factor (Figure S13). Holding variability for all correction/rescaling factors constant except for the fiber correction, the 95% confidence interval for percentage of samples in SFB exceeding Mehinto et al. (2022)’s food-dilution threshold one is (29–100%), compared with (76–88%) for size rescaling, (65–94%) for the plastic-proportion due to spectroscopic subsampling correction, and (26–100%) for combined rescaling and corrections (Figure S13). If the fiber correction factor is omitted from the analysis entirely, uncertainty decreases substantially in the risk characterization, and the overall number of statistically significant risk exceedances decreases as well (Figure S11). If the fiber correction factor is not applied, 27% of SFB samples would exceed food-dilution threshold one (95% CI: 18–39%) compared to 82% of samples when the fiber correction is applied (95% CI: 27–100%) (Figure S11 and Table S15).

While the size distribution value (a) has a substantial impact on the outcome of the risk characterization due to the high correction factor values derived for manta trawls (529; 95% CI: 401 to 704, Figure S9), the site-specific values for marine surface waters for SFB were of minimal difference from those derived elsewhere and applied for this risk characterization (i.e. Kooi et al. 2021) and therefore had limited uncertainty in this assessment (Figures S13 - S14). However, larger mesh sizes correspond to exponentially larger correction factors (Eq. 2) and are therefore highly influential in the case of manta trawl data (333 µm mesh). For example, the correction factor for 1-L grab samples would be 66 (95% CI: 55 to 80), using the same a value and uncertainty applied for manta trawl here, indicating the higher uncertainty and influence of rescaling manta trawl data compared to grab samples.

Comparison of the total uncertainties associated with estimating environmental surface water concentrations with the uncertainties in risk thresholds from Mehinto et al (2022) reveals comparable levels of uncertainty, with food dilution thresholds spanning ~ 2 to 5 orders of magnitude between 95% confidence intervals depending on the tier (Table S10) while estimated environmental concentrations for manta trawl surface samples span ~ 2.5 orders of magnitude between 95th percentiles based on the combined correction and rescaling uncertainties (Table S11).

Here, we combine occurrence data from SFB, California with a risk assessment framework to estimate the risk to local aquatic ecosystems. The risk framework includes hazard thresholds for two ecologically relevant categories of effect mechanisms - food dilution and mechanisms triggered upon tissue translocation. Based on the best available toxicological evidence and monitoring data, our results suggest that microplastic exposure in SFB in 2017 was high enough to cause biological perturbation through the food dilution effect mechanism. Eighty-two percent of the SFB had concentrations that exceeded Mehinto et al. 's (2022) tier one food-dilution threshold (“Investigative monitoring”) with statistical significance, with the highest percentages of statistical exceedances occurring within the Central Bay. Furthermore, the Central Bay was the only region within the SFB with any samples exceeding the third food dilution threshold (“Management planning”) with statistical significance, however no samples exceeded the highest food dilution threshold (need of immediate source control measures) with statistical significance. Because samples were not taken with the goal of being spatially or temporally representative of the SFB (Zhu et al. 2021), additional monitoring is suggested to improve confidence in risk characterizations.

Our analysis suggests that the risk associated with tissue translocation-mediated effects is substantially less than that for the food-dilution endpoint in SFB. While there were a few samples with concentrations greater than the first threshold at the 50th percentile, the exceedances were not statistically significant, and there were no samples above any of the other three thresholds at the 50th percentile. This lesser effect likely reflects tissue translocation-mediated effects being initiated by the subset of particles that are small enough (< 83 µm) to permeate the intestinal wall (Jeong et al. 2016; Qiao et al. 2019; Mehinto et al., 2022), whereas food-dilution is caused by a wider spectrum of particle sizes - based on ingestibility - that artificially fill the gut and lead to reduced food assimilation by blocked food passage (Lee et al. 2013; Murphy et al.2018; Blarer et al.2016; Cole et al.2015). Still, our analysis suggests that early stages of tissue translocation-mediated toxicities (e.g., oxidative stress, inflammation; Limonta et al., 2019; Xia et al. 2020) are possible in SFB, a finding that would not have been apparent without the rescaling procedures used to correct the underestimation of small particles captured in trawl nets (Koelmans et al. 2020). Rescaling particle counts based on size to correct for sampling bias in combination with toxicological thresholds aligned to ecologically relevant metrics provides the opportunity to compare exposure and hazard more appropriately (Koelmans et al. 2022).

Based on the species sensitivity distributions used to derive the food-dilution thresholds in Mehinto et al. (2022), the most sensitive species are the black-lip pearl oyster (Pinctada margaritifera), the marine medaka fish (Oryzias melastigma), and a water flea (Ceriodaphnia dubia), thus representing a diversity of taxonomic groups (mollusk, fish, and crustacea). Like many productive marine/freshwater estuarine systems, these three taxa are present in the SFB, with some similar species for which these laboratory model organisms may be suitable indicators. For instance, the Olympia oyster (Ostrea conchaphila) is native to the SFB and has experienced declining abundance which has been primarily attributed to loss of habitat and other factors (Parker and Boyer 2019). Experiments using model species within the Pinctada genus (e.g., Pinctada mazatlanica) have been used to inform risk management of the Olympia oyster in SFB (Sumudu 2008).

The ecological risk that we found for SFB was larger than that for several previous risk characterizations conducted for other geographies (e.g., Everaert et al. 2020; Koelmans et al. 2020; Adam et al. 2021). Key reasons for these differences are the use of different hazard threshold values and alignment procedures (or lack thereof), and that this sampling effort focused inside an urban enclosed water body with limited circulation. The SFB area has a large population of over 7 million people (MTC-ABAG Library 2022), and 39 WWTPs feeding into the bay (Hummel et al. 2018). This difference is confirmed by the application of our risk characterization methodology in the National Marine Sanctuary areas in the open ocean. Here, we did not find any samples exceeding any thresholds with statistical significance (Table S16), with the three highest samples taken at the mouth of the SFB (Fig. 4, interactive map in SI). Comparison of SFB concentrations to global marine surface water concentrations rescaled by Everaert et al. (2020) to the same size range used here (1 to 5,000 µm) using similar methods, reveals higher concentrations within SFB than in ~ 75% of global marine locations (Everaert et al 2020). While most marine monitoring data has been conducted in the open ocean, measurements in enclosed areas near urban centers indicate higher contamination. Everaert et al. (2020) reported 50 microplastics·L^− 1 in the Yellow Sea near China - an enclosed water body adjacent to a population of ~ 600 million people (Wang et al. 2016) - which is higher than ~ 98% of rescaled surface water samples in SFB. The high concentrations in the SFB and other enclosed water bodies demonstrates the importance of targeted monitoring to protect coastal resources globally.

While our probabilistic 95th percentiles based on Monte Carlo modeling and PDFs were sufficiently small to confidently state SFB contains microplastics at concentrations of biological concern, there is room for improvement in future risk characterizations to reduce uncertainties. Despite the research advancements that allowed us to compare exposure data to hazard data with higher certainty than what has been previously possible without the rescaling and realignment procedures developed in Koelmans et al. (2020), our risk estimates for SFB still contain substantial uncertainty and understanding the factors that contribute to that variability will help focus advancements needed to improve future estimates. Here we employed four key factors to calculate risk: ambient concentration measurements, critical thresholds at which biological effects manifest, size rescaling and other correction procedures used to correct for data collection biases, and the alignments of the biological thresholds (Koelmans et al. 2022). Understanding the relative contribution of these factors to variability will help improve precision of future studies.

Variability associated with the fiber correction factor represented the largest uncertainty in this risk characterization (figure S13). The ambient concentrations used to characterize risk were collected using surface manta trawls (Zhu et al., 2021), which are known to significantly undercount particles smaller than the mesh size (~ 333 µm), which includes the width-measurement of fibers allowing them to pass through the mesh like spaghetti (Dris et al., 2018; Zhu et al., 2020). For this reason, even after correcting for the uncounted fibers and rescaling concentrations to a common size range using compartment-specific size distributions, estimates generated by the 1-L grab samples taken at the same time as the trawl remained several orders of magnitude higher (Fig. 2), suggesting our SFB risk estimates based on corrected and rescaled manta trawl data were either an underestimate, or that the 1-L grab samples overcounted particles. Another possible cause for underestimation of number concentrations is the fact that the blank corrections were made based on color and morphology combinations, rather than polymer identity due to the use of spectroscopic subsampling as opposed to identification of all particles. This leads to underestimation of concentrations if microplastic particles of the same or similar color have a different polymer identity. We chose to rely on the manta trawl data because the 1-L grab sample volume was too small to provide accurate representation (Koelmans et al. 2019), with the higher percent relative standard deviations for duplicate manta samples (up to 13.3) being much lower than 4 duplicate grab samples (up to 46.9; Hung et al., 2020). This variability was so great that the manta and grab samples were not even significantly correlated (R² = 0.04), even though they were paired closely in space and time (Figure S12). Future risk characterization studies would be improved by applying polymer-specific data and blank corrections and using pumps with in-line filtration that include small mesh sizes (< 300 µm) to reduce small-scale spatial variability while capturing a fuller range of particle sizes than a towed net can (Tamminga et al., 2019; Lorenz et al. 2019; Mintenig et al. 2020).

The rescaling procedures described in Koelmans et al (2020) help with addressing some of the size undercounting biases associated with use of 333 µm mesh manta trawl nets, but these corrections yield additional sources of uncertainty, and complicate comparisons to other risk characterizations that do not perform rescaling procedures. The large mesh size of a manta trawl (333 µm) means that the size-based correction factor is relatively sensitive to the size exponent value used (see Table S8). For instance, the manta-trawl collected rescaled concentrations in SFB are nearly an order of magnitude higher than similarly rescaled data from manta-trawls reported in Everaert et al (2020) due to the larger size power exponent value used here. Everaert et al (2020) used a power exponent of 1.6 (± 0.5) corresponding to a correction factor of 40x for a 333 µm mesh, while a power exponent of 2.07 (± 0.03) was used here and corresponds to a correction factor of 530x for a 333 µm mesh. The power exponent used in Everaert et al (2020) was derived by Kooi and Koelmans (2019) using the best available data at the time - which was arbitrarily size binned data extracted from tables and graphs from other studies - and is therefore less accurate than the values derived in Kooi et al (2021) which used high-resolution datasets at the individual particle level. Using data from five studies that used state-of-the-art Fourier-transform infrared imaging and automated analysis, Kooi et al (2021) derived a length-based power exponent value of 2.07 (± 0.03) for marine surface waters, which is slightly smaller and substantially more certain than the site-specific value derived here using low-resolution manta trawl data measured using manual techniques (i.e., 2.15 ± 0.48). To reduce uncertainty due to size rescaling in SFB, small mesh-size samples could be obtained using in-line filtration and be analyzed using state-of-the-science analytical techniques to derive local data (e.g., Mintenig et al. 2020, Primpke et al. 2020).

A minor point of uncertainty was the spectroscopy correction factor applied to concentration data to ensure that only plastic particles were used to characterize risk (Figure S13). We did not consider particles that were clearly anthropogenic but were not unequivocally plastic, due to spectroscopic interferences from dyes and/or plastic additives (Zhu et al. 2021; see Figure S13). As such, the concentrations used here are also an underestimation for this reason. Additionally, unquantified inaccuracies in blank corrections based on color-shape combinations are expected to result in further underestimations of exposure and therefore risk. To reduce this uncertainty, future studies should use microplastic-specific spectral libraries to reduce the proportion of spectra that are less polymer-specific (Munno et al., 2019; Cowger et al. 2021) as well as automation to allow for chemical confirmation of all particles and polymer-specific blank corrections (Primpke et al. 2020).

Uncertainties in this risk characterization due to the selection of concentration thresholds at which effects manifest are both clearly illustrated and robust, as Mehinto et al. (2022) identified four thresholds that bracket the severity of response whereas other risk assessments have relied on single thresholds and tested the sensitivity based on relevant factors. While the variability associated with the distribution modeling component of the SSDs in Mehinto et al (2022) spanned ~ 2 to 5 orders of magnitude, the median threshold values were insensitive to individual studies - which had only approximately a two-fold influence. Furthermore, the eight threshold values used in this study from Mehinto et al. (2022) bracket the published range for microplastics (see Koelmans et al. 2022), with 95% confidence intervals based on SSDs being smaller in the food dilution thresholds in Mehinto et al. (2022) than those of thresholds derived in previous thresholds. Koelmans et. al (2020) developed and applied the same rescaling and alignment methodology to obtain an SSD using studies which demonstrated ingestion and suggested food dilution as the effect mechanism and derived a hazard concentration for 5% of species (HC5) of 76 particles/L (95% CI: 11 to 521 particles/L). Additional microplastics HC5 values which have not been rescaled or aligned span several orders of magnitude, however six out of nine published values are within confidence intervals of 76 particles/L (Everaert et al. 2018; Skåre et al. 2019; Zhang et al. 2020; Everaert et al. 2020; Jung et al. 2021), with three exceeding the range due to their inclusion of nanoplastics data (Adam et al. 2019; Besseling et al. 2019; Adam et al. 2021) as demonstrated in a review by Koelmans et al. (2022). There is room for additional studies to improve the threshold values (Hampton et al. this volume), with specific attention to experimental design to assess risk more accurately (de Rujiter et al., 2020) and additional studies using fibers (which were highly abundant in the SFB), however the uncertainty associated with the threshold values is still smaller than that of the exposure data.

While the Zhu et al. (2021) monitoring study in SFB included more than just surface water samples, here our risk assessment was limited to surface water exposure. At present, quantitative risk thresholds have not yet been developed for marine sediment exposure. However, microplastics are known to accumulate in sediments (Koelmans et al., 2017), and the SFB had relatively high concentrations of microplastics in the sediments. Moreover, we only considered the particle-induced effects of microplastics and did not account for additional potential risks resulting from pathogens (Zettler et al. 2013; Bowley et al., 2021), or potential chemical-mediated effects (see Koelmans et al. 2021 for summary of sorbed contaminants), including the leaching of chemical additives (Zimmerman et al. 2020) that preferentially desorb following ingestion (Coffin et al. 2019). In particular, tire wear particles were a large portion of the microplastics found in SFB sediments and chemical derivatives from those tire products have been found to cause acute mortality in salmon (Tian et al. 2021). Furthermore, microplastics are anticipated to interact with and exacerbate effects from additional stressors such as thermal stress due to climate change (Weber et al. 2020; Kolomijeca et al., 2020). Future work should seek to look holistically at risk in the SFB, including comparing relative risks from microplastics particles and other stressors such as dissolved and sorbed chemicals, for which risk assessment frameworks already exist.

Overall, our results indicate that SFB has regions where present exposure concentrations are above thresholds for risk with statistical confidence based on the best available ecotoxicological hazard thresholds. If we continue business as usual, inputs of microplastics to coastal environments are anticipated to triple over the next twenty years (Lebreton and Andrady 2020). As such, the region might consider management actions now to prevent greater risk in the future. Beyond San Francisco Bay, our study can inform risk assessments in local regions across the globe. Combined these methods can be used globally to inform management locally.

ANOVA - analysis of variance
PDF - probability density function
SFB - San Francisco Bay, California

PNEC - predicted no-effect concentration

PEC - predicted environmental concentration

SSD - species sensitivity distribution

Ethics Declaration statement

Not applicable.

Consent to Participate statement

Not applicable.

Consent for publication

All authors have reviewed the final manuscript and consent to its publication.

Availability of data and material

Code used for modeling is available in https://github.com/ScottCoffin/SFBayMPRiskCharacterization/tree/Manuscript

Size-rescaled, blank-corrected, fiber fraction-corrected, and plastic spectroscopic sub-sampling corrected microplastics occurrence data are available in the supplemental information as a comma-separated value file (.csv).

An interactive map showing concentrations of microplastics and risk exceedances is available in the supplemental information as a HyperText Markup Language (.html) file.

Competing interests

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Funding

Scott Coffin was supported through funding from the California State Water Resources Control Board. Chelsea Rochman and Stephen Weisberg did not receive funding for this work. Albert Koelmans and Merel Kooi were supported through funding from Wageningen University & Research.

Author Contributions

S.C. wrote the main manuscript text, conceptualized, developed, and performed the methodology, conducted the formal analysis, curated the data, produced visualizations, and administered the project. S.B.W. wrote and edited the main manuscript text and conceptualized the project. C.R. wrote and edited the main manuscript text, curated the data, and conceptualized the project. M.K. edited the main manuscript and conducted formal analysis for power laws. A.A.K. wrote and edited the main manuscript and conceptualized the project. All authors reviewed the manuscript.

Acknowledgments

The authors are grateful to Dr. Diana Lin and other staff at the San Francisco Estuary Institute for assisting in accessing data, advising on methods, and reviewing the manuscript.

Author’s information

The views and opinions expressed in this article are those of the authors and do not necessarily reflect the official policy or position of any government agency or institution.

Abel SM, Primpke S, Int-Veen I, Brandt A, Gerdts G. Systematic identification of microplastics in abyssal and hadal sediments of the Kuril Kamchatka trench. Environmental Pollution. 2021;269:116095.
Adam V, von Wyl A, Nowack B. Probabilistic environmental risk assessment of microplastics in marine habitats. Aquatic Toxicology. 2021 Jan;230:105689.
Adam V, Yang T, Nowack B. Toward an ecotoxicological risk assessment of microplastics: Comparison of available hazard and exposure data in freshwaters. Environ Toxicol Chem. 2019 Feb;38(2):436–47.
Bäuerlein PS, Erich MW, van Loon W, Bakker I, Mintenig SM, Koelmans AA. A monitoring and data analysis method for microplastics in marine sediments for OSPAR and MSFD. Mar Poll Bull. Submitted.;
Besseling E, Redondo-Hasselerharm P, Foekema EM, Koelmans AA. Quantifying ecological risks of aquatic micro- and nanoplastic. Critical Reviews in Environmental Science and Technology. 2019 Jan 2;49(1):32–80.
Blarer P, Burkhardt-Holm P. Microplastics affect assimilation efficiency in the freshwater amphipod Gammarus fossarum. Environmental science and pollution research. 2016;23(23):23522–32.
Bowley J, Baker-Austin C, Porter A, Hartnell R, Lewis C. Oceanic hitchhikers–assessing pathogen risks from marine microplastic. Trends in microbiology. 2021;29(2):107–16.
Brander SM, Renick VC, Foley MM, Steele C, Woo M, Lusher A, et al. Sampling and QA/QC: A guide for scientists investigating the occurrence of microplastics across matrices.:52.
Bucci K, Rochman CM. Microplastics: a multidimensional contaminant requires a multidimensional framework for assessing risk. Micropl&Nanopl. 2022 Dec;2(1):7.
Burns EE, Boxall ABA. Microplastics in the aquatic environment: Evidence for or against adverse impacts and major knowledge gaps: Microplastics in the environment. Environ Toxicol Chem. 2018 Nov;37(11):2776–96.
Clauset A, Shalizi CR, Newman ME. Power-law distributions in empirical data. SIAM review. 2009;51(4):661–703.
Coffin S, Huang G-Y, Lee I, Schlenk D. Fish and Seabird Gut Conditions Enhance Desorption of Estrogenic Chemicals from Commonly-Ingested Plastic Items. Environ Sci Technol. 2019 Apr 16;53(8):4588–99.
Cole M, Galloway TS. Ingestion of nanoplastics and microplastics by Pacific oyster larvae. Environmental science & technology. 2015;49(24):14625–32.
Cook RD. Detection of influential observation in linear regression. Technometrics. 1977;19(1):15–8.
Covernton GA, Pearce CM, Gurney-Smith HJ, Chastain SG, Ross PS, Dower JF, et al. Size and shape matter: A preliminary analysis of microplastic sampling technique in seawater studies with implications for ecological risk assessment. Science of The Total Environment. 2019 Jun;667:124–32.
Cowger W, Gray A, Christiansen SH, DeFrond H, Deshpande AD, Hemabessiere L, et al. Critical Review of Processing and Classification Techniques for Images and Spectra in Microplastic Research. Appl Spectrosc. 2020 Sep;74(9):989–1010.
Cowger W, Steinmetz Z, Gray A, Munno K, Lynch J, Hapich H, et al. Microplastic Spectral Classification Needs an Open Source Community: Open Specy to the Rescue! Anal Chem [Internet]. 2021 May 19; Available from: https://doi.org/10.1021/acs.analchem.1c00123
de Ruijter VN, Redondo Hasselerharm PE, Gouin T, Koelmans AA. Quality criteria for microplastic effect studies in the context of risk assessment: A critical review. Environ Sci Technol. 2020 Aug 28;acs.est.0c03057.
Delignette-Muller ML, Dutang C. fitdistrplus: An R package for fitting distributions. Journal of statistical software. 2015;64:1–34.
Everaert G, De Rijcke M, Lonneville B, Janssen CR, Backhaus T, Mees J, et al. Risks of floating microplastic in the global ocean. Environmental Pollution. 2020 Dec;267:115499.
Everaert G, Van Cauwenberghe L, De Rijcke M, Koelmans AA, Mees J, Vandegehuchte M, et al. Risk assessment of microplastics in the ocean: Modelling approach and first conclusions. Environmental Pollution. 2018 Nov;242:1930–8.
Gardon T, Huvet A, Paul-Pont I, Cassone A-L, Koua MS, Soyez C, et al. Toxic effects of leachates from plastic pearl-farming gear on embryo-larval development in the pearl oyster Pinctada margaritifera. Water research. 2020;179:115890.
Gillespie CS. Fitting heavy tailed distributions: the poweRlaw package. arXiv preprint arXiv:14073492. 2014;
González-Pleiter M, Velázquez D, Edo C, Carretero O, Gago J, Barón-Sola Á, et al. Fibers spreading worldwide: Microplastics and other anthropogenic litter in an Arctic freshwater lake. Science of the Total Environment. 2020;722:137904.
Gouin T, Becker RA, Collot A, Davis JW, Howard B, Inawaka K, et al. Toward the development and application of an environmental risk assessment framework for microplastic. Environmental toxicology and chemistry. 2019;38(10):2087–100.
Hampton LMT, Brander S, Coffin S, Cole M, Hermabessiere L, Koelmans A, et al. Characterizing microplastic hazards: Which concentration metrics and particle characteristics are most informative for understanding toxicity in aquatic organisms? Microplastics and Nanoplastics. 2022.
Hermsen E, Mintenig SM, Besseling E, Koelmans AA. Quality criteria for the analysis of microplastic in biota samples: a critical review. Environmental science & technology. 2018;52(18):10230–40.
Hummel MA, Berry MS, Stacey MT. Sea Level Rise Impacts on Wastewater Treatment Systems Along the U.S. Coasts. Earth’s Future. 2018 Apr;6(4):622–33.
Hung C, Klasios N, Zhu X, Sedlak M, Sutton R, Rochman CM. Methods Matter: Methods for Sampling Microplastic and Other Anthropogenic Particles and Their Implications for Monitoring and Ecological Risk Assessment. Integr Environ Assess Manag. 2021 Jan;17(1):282–91.
Jeong C-B, Won E-J, Kang H-M, Lee M-C, Hwang D-S, Hwang U-K, et al. Microplastic size-dependent toxicity, oxidative stress induction, and p-JNK and p-p38 activation in the monogonont rotifer (Brachionus koreanus). Environmental science & technology. 2016;50(16):8849–57.
Jung J-W, Park J-W, Eo S, Choi J, Song YK, Cho Y, et al. Ecological risk assessment of microplastics in coastal, shelf, and deep sea waters with a consideration of environmentally relevant size and shape. Environmental Pollution. 2021 Feb 1;270:116217.
Kahle DJ, Wickham H. ggmap: spatial visualization with ggplot2. R J. 2013;5(1):144.
Kirchner J. Data analysis toolkit# 5: uncertainty analysis and error propagation. University of California Berkeley Seismological Laboratory Available online at: http://seismo berkeley edu/~ kirchner/eps_120/Toolkits/Toolkit_05 pdf. 2001;
Klasios N, De Frond H, Miller E, Sedlak M, Rochman CM. Microplastics and other anthropogenic particles are prevalent in mussels from San Francisco Bay, and show no correlation with PAHs. Environmental Pollution. 2021 Feb;271:116260.
Koelmans AA, Kooi M, Law KL, Van Sebille E. All is not lost: deriving a top-down mass budget of plastic at sea. Environmental Research Letters. 2017;12(11):114028.
Koelmans AA, Mohamed Nor NH, Hermsen E, Kooi M, Mintenig SM, De France J. Microplastics in freshwaters and drinking water: Critical review and assessment of data quality. Water Research. 2019 May;155:410–22.
Koelmans AA, Redondo Hasselerharm PE, Mohamed Nor NH, Kooi M. Solving the non-alignment of methods and approaches used in microplastic research in order to consistently characterize risk. Environ Sci Technol. 2020 Sep 4;acs.est.0c02982.
Koelmans AA, Redondo-Hasselerharm PE, Nor NHM, de Ruijter VN, Mintenig SM, Kooi M. Risk assessment of microplastic particles. Nat Rev Mater [Internet]. 2022 Jan 21 [cited 2022 Jan 24]; Available from: https://www.nature.com/articles/s41578-021-00411-y
Koelmans A, Diepens N, Hazimah N. Weight of Evidence for the Microplastic Vector Effect in the Context of Chemical Risk Assessment. In: Microplastic in the Environment: Pattern and Process [Internet]. Springer; 2021. Available from: https://doi.org/10.1007/978-3-030-78627-4_6
Kolomijeca A, Parrott J, Khan H, Shires K, Clarence S, Sullivan C, et al. Increased temperature and turbulence alter the effects of leachates from tire particles on fathead minnow (Pimephales promelas). Environmental science & technology. 2020;54(3):1750–9.
Kooi M, Koelmans AA. Simplifying Microplastic via Continuous Probability Distributions for Size, Shape, and Density. Environ Sci Technol Lett. 2019 Sep 10;6(9):551–7.
Kooi M, Primpke S, Mintenig SM, Lorenz C, Gerdts G. Characterizing the multidimensionality of microplastics across environmental compartments. Water Research. 2021;24.
Lee K-W, Shim WJ, Kwon OY, Kang J-H. Size-dependent effects of micro polystyrene particles in the marine copepod Tigriopus japonicus. Environmental science & technology. 2013;47(19):11278–83.
Limonta G, Mancia A, Benkhalqui A, Bertolucci C, Abelli L, Fossi MC, et al. Microplastics induce transcriptional changes, immune response and behavioral alterations in adult zebrafish. Scientific reports. 2019;9(1):1–11.
Mintenig S, Kooi M, Erich M, Primpke S, Redondo-Hasselerharm P, Dekker S, et al. A systems approach to understand microplastic occurrence and variability in Dutch riverine surface waters. Water research. 2020;176:115723.
MTC-ABAG Library. San Francisco Bay Area 2010 Decennial Census Data [Internet]. MTC-ABAG Library; 2022 Feb. Available from: http://www.bayareacensus.ca.gov/bayarea.htm
Murphy F, Quinn B. The effects of microplastic on freshwater Hydra attenuata feeding, morphology & reproduction. Environmental pollution. 2018;234:487–94.
Newcombe RG. Propagating Imprecision: Combining Confidence Intervals from Independent Sources. Communications in Statistics - Theory and Methods. 2011 Sep;40(17):3154–80.
Newman ME. Power laws, Pareto distributions and Zipf’s law. Contemporary physics. 2005;46(5):323–51.
Parker VT, Boyer KE. Sea-level rise and climate change impacts on an urbanized Pacific Coast estuary. Wetlands. 2019;39(6):1219–32.
Paula Elisa R-H, A. R, Koelmans, AA. Risk assessment of microplastics for freshwater benthic ecosystems guided by strict quality criteria and data alignment methods. In Press;
Pirsaheb M, Hossini H, Makhdoumi P. Review of microplastic occurrence and toxicological effects in marine environment: Experimental evidence of inflammation. Process Safety and Environmental Protection. 2020;142:1–14.
Primpke S, Christiansen SH, Cowger W, De Frond H, Deshpande A, Fischer M, et al. Critical Assessment of Analytical Methods for the Harmonized and Cost-Efficient Analysis of Microplastics. Appl Spectrosc. 2020 Sep;74(9):1012–47.
Qiao R, Sheng C, Lu Y, Zhang Y, Ren H, Lemos B. Microplastics induce intestinal inflammation, oxidative stress, and disorders of metabolome and microbiome in zebrafish. Science of the Total Environment. 2019;662:246–53.
Redondo-Hasselerharm PE, Falahudin D, Peeters ET, Koelmans AA. Microplastic effect thresholds for freshwater benthic macroinvertebrates. Environmental science & technology. 2018;52(4):2278–86.
Rubin AE, Sarkar AK, Zucker I. Questioning the suitability of available microplastics models for risk assessment – A critical review. Science of The Total Environment. 2021 Sep;788:147670.
Scopetani C, Esterhuizen-Londt M, Chelazzi D, Cincinelli A, Setälä H, Pflugmacher S. Self-contamination from clothing in microplastics research. Ecotoxicology and Environmental Safety. 2020 02;189:110036.
Skåre JU, Alexander J, Haave M, Jakubowicz I, Knutsen HK, Lusher A, et al. Microplastics; occurrence, levels and implications for environment and human health related to food. Scientific opinion of the Scientific Steering Committee of the Norwegian Scientific Committee for Food and Environment. VKM Report. 2019;
Tamminga M, Stoewer S-C, Fischer EK. On the representativeness of pump water samples versus manta sampling in microplastic analysis. Environmental pollution. 2019;254:112970.
Wang Q, Wangjin X, Zhang Y, Wang N, Wang Y, Meng G, et al. The toxicity of virgin and UV-aged PVC microplastics on the growth of freshwater algae Chlamydomonas reinhardtii. Science of The Total Environment. 2020;749:141603.
Wang Q, Song J, Zhou J, Zhao W, Liu H, Tang X. Temporal Evolution of the Yellow Sea Ecosystem Services (1980–2010). Heliyon. 2016 Mar 1;2(3):e00084.
Weber A, Jeckel N, Wagner M. Combined effects of polystyrene microplastics and thermal stress on the freshwater mussel Dreissena polymorpha. Science of The Total Environment. 2020 May;718:137253.
Welaratna S. The native oyster recruitment study in central and south San Francisco Bay 2006–2007. San Jose State University; 2008.
Werbowski LM, Gilbreath AN, Munno K, Zhu X, Grbic J, Wu T, et al. Urban Stormwater Runoff: A Major Pathway for Anthropogenic Particles, Black Rubbery Fragments, and Other Types of Microplastics to Urban Receiving Waters. ACS EST Water. 2021 Jun 11;1(6):1420–8.
Xia X, Sun M, Zhou M, Chang Z, Li L. Polyvinyl chloride microplastics induce growth inhibition and oxidative stress in Cyprinus carpio var. larvae. Science of The Total Environment. 2020;716:136479.
Zettler ER, Mincer TJ, Amaral-Zettler LA. Life in the “plastisphere”: microbial communities on plastic marine debris. Environmental science & technology. 2013;47(13):7137–46.
Zhang C, Jeong C-B, Lee J-S, Wang D, Wang M. Transgenerational Proteome Plasticity in Resilience of a Marine Copepod in Response to Environmentally Relevant Concentrations of Microplastics. Environ Sci Technol. 2019 Jul 16;53(14):8426–36.
Zhang X, Leng Y, Liu X, Huang K, Wang J. Microplastics’ pollution and risk assessment in an urban river: A case study in the Yongjiang River, Nanning City, South China. Exposure and Health. 2020;12(2):141–51.
Zhu X, Munno K, Grbic J, Werbowski LM, Bikker J, Ho A, et al. Holistic Assessment of Microplastics and Other Anthropogenic Microdebris in an Urban Bay Sheds Light on Their Sources and Fate. ACS EST Water. 2021 Jun 11;1(6):1401–10.
Zimmermann L, Göttlich S, Oehlmann J, Wagner M, Völker C. What are the drivers of microplastic toxicity? Comparing the toxicity of plastic chemicals and particles to Daphnia magna. Environmental Pollution. 2020 Dec;267:115392.

No competing interests reported.

SupplementalInformation.pdf

Download PDF

Editorial decision: Major revision
27 Apr, 2022
Reviews received at journal
24 Apr, 2022
Reviewers agreed at journal
11 Apr, 2022
Reviewers invited by journal
07 Apr, 2022
Editor assigned by journal
07 Apr, 2022
Submission checks completed at journal
06 Apr, 2022
First submitted to journal
30 Mar, 2022

You are reading this latest preprint version

Risk Characterization of Microplastics in San Francisco Bay, California

Status:

Version 1

Abstract

Figures

Introduction

Materials And Methods

Results

Discussion

Conclusion

Abbreviations

Declarations

Ethics Declaration statement

Consent to Participate statement

Consent for publication

Availability of data and material

Competing interests

Funding

Author Contributions

Acknowledgments

Author’s information

References

Additional Declarations

Supplementary Files

Status:

Version 1