Daniel B ShawpersName, Qi LipersName, Janine K NunespersName, Luc DeikepersName (2023date). “Ocean emission of microplastic.” PNAS Nexustitle·j.
Object of study — transport of microplasticterm·dom from the ocean into the air by jet dropsterm·dom of bursting bubbles, and its global emissionterm·key
Research question — how much microplasticterm·dom can be emitted from the ocean by bubble-bursting sea sprayterm·key
Contribution — a mechanistic, experimentally calibrated global model of ocean microplasticterm·dom emissionterm·key
phenomenon · claim-type · cross-domain-transfer · mechanism · claim-relation · replication-status · limitation — grounded in the full text| phenomenon | microplastic emission from the ocean by bursting bubbles — “We demonstrate and quantify that jet drops are efficient at emitting microplastics up to 280 μ m in diameter and are thus expected to dominate the emitted mass of microplastic.” | |
| claim-type | mechanism — “We demonstrate and quantify that jet drops are efficient at emitting microplastics up to 280 μ m in diameter and are thus expected to dominate the emitted mass of microplastic.” | |
| method — evidence | — “By varying the liquid properties, bubble size, depth of bubble rise, microplastic size and concentration, equations for jet drop capture of microplastic are developed.” | |
| object of study — evidence | — “We present a mechanistic study of bursting bubbles transporting microplastics.” | |
| research question — evidence | — “Sea spray from bubble bursting can transport salt and biological material from the ocean into the atmosphere, and there is a need to quantify the amount of microplastic that can be emitted from the ocean by this mechanism.” | |
| contribution — evidence | — “Our study provides an independent estimate of the microplastic emission by the ocean and reduces the associated uncertainties by identifying and quantifying the leading emission mechanism.” | |
| mechanism | loss of surface mobility A rising bubble loses its free-slip condition as it scavenges particles and surfactants, slowing collection. — “The N ∝ ( H / R b ) 1 / 2 scaling suggests that when considering large variations of H / R b , the bubble’s surface looses its free-slip condition as it scavenges microplastic particles and surface-active molecules ( 43 ), transitioning to a nonslip condition ( 54–57 ).” | mechanisms[0] |
| mechanism | breaking waves Ocean bubbles are created by breaking waves entraining air. — “Bubbles on the ocean are created by breaking waves as air is entrained under the surface.” | mechanisms[1] |
| sense | ejectable microplastic - particles no larger than the jet drops that can carry them — “We demonstrate and quantify that jet drops are efficient at emitting microplastics up to 280 μ m in diameter and are thus expected to dominate the emitted mass of microplastic.” | senses[0] |
| relation | nonslip surface —contrasted-with→ free-slip surface — always less than — “Indeed the collection rate for a sphere with a nonslip surface is always less than that of one with a free-slip surface for a fixed particle and bubble size ( 43 , 58 ), which provides a rationale for the observed N ∝ ( H / R b ) 1 / 2 scaling.” — “the collection rate for a sphere with a nonslip surface is always less than that of one with a free-slip surface” | concept_edges[0] |
| weight | nonslip surface —contrasted-with→ free-slip surface — 0.85 — “Indeed the collection rate for a sphere with a nonslip surface is always less than that of one with a free-slip surface for a fixed particle and bubble size ( 43 , 58 ), which provides a rationale for the observed N ∝ ( H / R b ) 1 / 2 scaling.” | concept_edges[0] |
| effect direction | bubble rise height —causes→ scavenging rate — reduces — “Because the scavenging rate is found to decrease with bubble rise height, the total amount of microplastic ejected depends on the ocean conditions through both F d and the bubble scavenging.” — “the scavenging rate is found to decrease with bubble rise height” | concept_edges[1] |
| relation | bubble rise height —causes→ scavenging rate — decrease with — “Because the scavenging rate is found to decrease with bubble rise height, the total amount of microplastic ejected depends on the ocean conditions through both F d and the bubble scavenging.” — “the scavenging rate is found to decrease with bubble rise height” | concept_edges[1] |
| weight | bubble rise height —causes→ scavenging rate — 0.75 — “Because the scavenging rate is found to decrease with bubble rise height, the total amount of microplastic ejected depends on the ocean conditions through both F d and the bubble scavenging.” | concept_edges[1] |
| relation | SDS —same-as→ sodium dodecyl sulfate — (SDS) — “A small amount of sodium dodecyl sulfate (SDS) is added to the liquid at a concentration of 6.6 μ m which is several orders of magnitude below the critical micelle concentration of SDS of 8.2 mM .” — “sodium dodecyl sulfate (SDS)” | concept_edges[2] |
| weight | SDS —same-as→ sodium dodecyl sulfate — 0.95 — “A small amount of sodium dodecyl sulfate (SDS) is added to the liquid at a concentration of 6.6 μ m which is several orders of magnitude below the critical micelle concentration of SDS of 8.2 mM .” | concept_edges[2] |
| evidential status | mechanism — suggests — “The N ∝ ( H / R b ) 1 / 2 scaling suggests that when considering large variations of H / R b , the bubble’s surface looses its free-slip condition as it scavenges microplastic particles and surface-active molecules ( 43 ), transitioning to a nonslip condition ( 54–57 ).” | mechanisms[0] |
| replication | original | |
| limitation | The dominant uncertainty is scarce, coarse data on ocean microplastic concentration. — “However, the dominant uncertainty of our global emission estimate is the limited knowledge for ocean microplastic concentration; it is both scarce in temporal and spatial coverage and has a coarse size resolution ( 3 , 8 ).” | limitations[0] |
| limitation | The model gives emission from the surface, not the flux reaching the upper atmosphere. — “As such, our model describes the emission of microplastic from the ocean surface, not the flux that is actually transported to the upper atmosphere.” | limitations[1] |
| limitation | Annual-mean concentrations against hourly emission may miss nonlinear temporal effects. — “As such, nonlinearities between temporal variations in the microplastic concentration and jet drop production flux could affect the final result due to this difference in timescale sampling.” | limitations[2] |
| limitation | Biofouling was not considered. — “Biofouling is another effect which was not considered in this work, and it has been shown to change the effective density of microplastic particles ( 75 ).” | limitations[3] |
| limitation | Only spherical particles were used. — “While only spherical particles were used in this study, the emission of other shapes has been studied previously.” | limitations[4] |
finding · quantity · variable · definition · key-term — grounded in the full text| finding | Jet drops efficiently emit microplastics up to 280 μm and should dominate the emitted mass. — “We demonstrate and quantify that jet drops are efficient at emitting microplastics up to 280 μ m in diameter and are thus expected to dominate the emitted mass of microplastic.” | content_tags.findings[0] |
| finding | No microplastic particles were seen in the thin-film cap. — “Throughout the entire study, no microplastic particles are observed in the thin-film cap which is O ( 0.1 − 10 μ m ) thick and smaller than the size of particles studied.” | content_tags.findings[1] |
| finding | Jet drop sizes with microplastics follow scaling laws derived for drops without them. — “The size of jet drops containing microplastics is found to be well described by theoretical scaling laws derived for drops without microplastics.” | content_tags.findings[2] |
| finding | High coastal emissions follow from high coastal microplastic concentration. — “High emissions near the coast are the result of the high concentration of microplastic there.” | content_tags.findings[3] |
| finding | The remaining uncertainty lies in limited knowledge of ocean microplastic. — “We show that the remaining uncertainty resides in the limited knowledge of ocean microplastic.” | content_tags.findings[4] |
| quantity | 0.02 to 7.4 Mt/yr — — “Results are then integrated globally, accounting for how air bubbles are produced and ocean measurements of microplastic, to estimate the global ocean microplastic emission which ranges from 0.02 to 7.4 Mt/yr with a best guess of 0.1 Mt/yr.” | content_tags.quantities[0] |
| quantity | 19 to 23 Mt — — “It is estimated that 19 to 23 Mt of plastic is currently flowing into the ocean annually ( 12 , 13 ).” | content_tags.quantities[1] |
| quantity | 8.6 Mt/yr — — “The upper bound is similar in magnitude to the best guess value of 8.6 Mt/yr from Ref. ( 5 ).” | content_tags.quantities[2] |
| variable | scavenging rate — “Because the scavenging rate is found to decrease with bubble rise height, the total amount of microplastic ejected depends on the ocean conditions through both F d and the bubble scavenging.” | content_tags.variables[0] |
| definition | Film drops — “Film drops come from the liquid that was in the bubble’s thin-film cap at the moment of bursting.” | content_tags.definitions[0] |
| definition | Jet drops — “Jet drops are formed from the collapse of the bubble’s underwater cavity and are responsible for most super-micron drops ( 16 , 19 , 29 );” | content_tags.definitions[1] |
| key-terms | microplastic, jet drops, bubble bursting, sea spray, scavenging, emission, film drops, efficiency factor |
| concept | Microplastic — typically defined as plastic particles of size between 1 μ m and 5 mm ( 3 ) — “Microplastic, typically defined as plastic particles of size between 1 μ m and 5 mm ( 3 ), is found in the atmosphere, and recent papers suggest that ocean is a potentially significant source of atmospheric microplastic ( 4–6 ).” | concepts[0] |
| concept | Film drops — come from the liquid that was in the bubble’s thin-film cap at the moment of bursting — “Film drops come from the liquid that was in the bubble’s thin-film cap at the moment of bursting.” | concepts[1] |
| concept | Jet drops — are formed from the collapse of the bubble’s underwater cavity — “Jet drops are formed from the collapse of the bubble’s underwater cavity and are responsible for most super-micron drops ( 16 , 19 , 29 );” | concepts[2] |
| concept | efficiency or enrichment factor — It physically represents how much more concentrated the particles are in the jet drop compared to the bulk liquid — “where E is defined as an efficiency or enrichment factor ( 35 , 36 ). It physically represents how much more concentrated the particles are in the jet drop compared to the bulk liquid.” | concepts[3] |
| concept | nonslip surface — not defined in this paper's own words; the bubble surface state reached as it scavenges surface-active molecules — “Indeed the collection rate for a sphere with a nonslip surface is always less than that of one with a free-slip surface for a fixed particle and bubble size ( 43 , 58 ), which provides a rationale for the observed N ∝ ( H / R b ) 1 / 2 scaling.” | concepts[4] |
| concept | free-slip surface — not defined in this paper's own words; a bubble surface that is completely mobile — “Indeed the collection rate for a sphere with a nonslip surface is always less than that of one with a free-slip surface for a fixed particle and bubble size ( 43 , 58 ), which provides a rationale for the observed N ∝ ( H / R b ) 1 / 2 scaling.” | concepts[5] |
| concept | jet drop emission function — provides the size distribution per unit ocean area per unit time — “The jet drop emission function provides the size distribution per unit ocean area per unit time and is sensitive to wind, waves and sea surface temperature ( 20 ).” | concepts[6] |
| concept | scavenging rate — not defined in this paper's own words; the rate at which a rising bubble collects particles — “Because the scavenging rate is found to decrease with bubble rise height, the total amount of microplastic ejected depends on the ocean conditions through both F d and the bubble scavenging.” | concepts[7] |
| concept | bubble rise height — not defined in this paper's own words; the height a bubble rises to the surface — “Because the scavenging rate is found to decrease with bubble rise height, the total amount of microplastic ejected depends on the ocean conditions through both F d and the bubble scavenging.” | concepts[8] |
| concept | sodium dodecyl sulfate — not defined in this paper's own words; a surfactant added to weakly contaminate the surface — “A small amount of sodium dodecyl sulfate (SDS) is added to the liquid at a concentration of 6.6 μ m which is several orders of magnitude below the critical micelle concentration of SDS of 8.2 mM .” | concepts[9] |
| concept | SDS — sodium dodecyl sulfate — “A small amount of sodium dodecyl sulfate (SDS) is added to the liquid at a concentration of 6.6 μ m which is several orders of magnitude below the critical micelle concentration of SDS of 8.2 mM .” | concepts[10] |
Full text — tagged
The reference list is not part of the text we hold for this paper. Our copy was extracted from the publisher's page and the bibliography did not survive that extraction, so the article below ends before its references. Nothing has been removed by us.
Ocean emissionterm·key of microplasticterm·dom Abstract Microplasticsterm·dom are globally ubiquitous in marine environments, and their concentration is expected to continue rising at significant rates as a result of human activity. They present a major ecological problem with well-documented environmental harm. Sea sprayterm·key from bubble burstingterm·key can transport salt and biological material from the ocean into the atmosphere, and there is a need to quantify the amount of microplasticterm·dom that can be emitted from the ocean by this mechanism. We present a mechanistic study of bursting bubbles transporting microplasticsterm·dom. We demonstrate and quantify that jet dropsterm·dom are efficient at emitting microplasticsterm·dom up to 280 μ m in diameter and are thus expected to dominate the emitted mass of microplasticterm·dom. The results are integrated to provide a global microplasticterm·dom emissionterm·key model which depends on bubble scavengingterm·key and bursting physics; local wind and sea state; and oceanic microplasticterm·dom concentration. We test multiple possible microplasticterm·dom concentration maps to find annual emissionsterm·key ranging from 0.02 to 7.4—with a best guess of 0.1—mega metric tons per year and demonstrate that while we significantly reduce the uncertainty associated with the bursting physics, the limited knowledge and measurements on the mass concentration and size distribution of microplasticterm·dom at the ocean surface leaves large uncertainties on the amount of microplasticterm·dom ejected.
Microplasticterm·dom particles are increasingly prevalent in the ocean. Bursting bubbles are known to transmit salt crystals and organic materials from the ocean into the atmosphere: influencing the radiative balance and serving as cloud condensation nuclei. We demonstrate through laboratory experimentsterm·mth that bursting bubble jet dropsterm·dom can scavenge and eject microplasticterm·dom (10 to 280 μ m in size) to the air. Results are then integrated globally, accounting for how air bubbles are produced and ocean measurements of microplasticterm·dom, to estimate the global ocean microplasticterm·dom emissionterm·key which ranges from 0.02 to 7.4 Mt/yr with a best guess of 0.1 Mt/yr. We show that the remaining uncertainty resides in the limited knowledge of ocean microplasticterm·dom.
Introduction The increasing prevalence of plastics in the ocean is a global-scale issue with wide-ranging impacts. Considerable scientific focus has been placed on ocean plastic debris since the 1970s ( 1 ) with an estimated 10% of all plastic produced eventually being deposited into the ocean ( 2 ). Microplasticterm·dom, typically defined as plastic particles of size between 1 μ m and 5 mm ( 3 ), is found in the atmosphere, and recent papers suggest that ocean is a potentially significant source of atmospheric microplasticterm·dom ( 4–6 ). Estimates range from 0 to 22 mega metric tons, Mt, per year ( 5 , 7 ); at the upper range, oceans would be one of the largest sources of atmospheric microplasticterm·dom ( 5 ), while other studies estimate that microplasticterm·dom emissionsterm·key by the ocean are negligible ( 7 ). As such, understanding microplasticterm·dom emissionterm·key by the ocean is an urgent unmet need to close global plastic budgets ( 3 , 5 ).
Plastic debris can be found in all marine environments ( 8 ) and are being transported by ocean currents ( 9 , 10 ) and waves ( 11 ). It is estimated that 19 to 23 Mt of plastic is currently flowing into the ocean annually ( 12 , 13 ). If current trends continue, it is expected that by 2040 the annual rate of microplasticterm·dom entering aquatic environments from land will have increased by 260% from 2016 ( 14 ) and a “peak plastic waste” is not expected to be reached until 2100 ( 15 ). Their increasing pervasiveness in the world’s environments presents a serious issue motivating the accurate quantification of the microplasticsterm·dom cycling in and out of the ocean.
Material such as water, salts, and biological material have long been known to be transported from the ocean to the atmosphere via sea sprayterm·key droplets ( 16–20 ). The rate of transport is high enough to affect global climate dynamics; salt crystals and organic aerosols influence the radiative balance of the atmosphere and serve as cloud condensation nuclei ( 21–24 ).
Sea sprayterm·key aerosols are generated by two pathways: spume drops resulting from high wind shearing of wave crests ( 19 , 25 ) and surface bubbles bursting ( 16 , 26 )—itself decomposed into film and jet dropsterm·dom ( 16 , 17 ). Film dropsterm·dom come from the liquid that was in the bubble’s thin-film cap at the moment of bursting. Their size is controlled by the bubble’s radius and the cap thickness at burst, ( 17 , 27 ) and they are responsible for most submicron spray drops ( 16 , 19 , 28 ). Jet dropsterm·dom are formed from the collapse of the bubble’s underwater cavity and are responsible for most super-micron drops ( 16 , 19 , 29 ); capillary waves travel down the bubble’s empty cavity to focus and form a jet which destabilizes into droplets ( 30–34 ).
As bubbles rise to the surface, they scavenge soluble and insoluble material ( 18 , 35–40 ), leading to an enhancement of the material’s concentration in the ejected droplet. The efficiency factorterm·key E , by which the concentration in the drop changes, increases with the height a bubble rises to the surface, H ; a linear scaling is predicted by a simple interceptor model ( 41 , 42 ). However, laboratory experimentsterm·mth have reported a decrease in the rate of material collection with increasing H ( 18 , 43 , 44 ), which can be related to saturation of the bubble’s surface area with particles ( 35 ) and decreasing surface mobility of the bubble as it rises due to the scavengingterm·key of surface-active molecules ( 43 ). Correction factors have been proposed for the scavengingterm·key of nonspherical bacteria ( 35 ), and the origin of liquid in the jet dropterm·dom has been modeled numerically ( 36 ).
In this work, we characterize and quantify how microplasticterm·dom is ejected from the ocean by bubble burstingterm·key as shown in Fig. 1 . An experimental study of microplasticterm·dom ejection by individual bubbles is first presented. By varying the liquid properties, bubble size, depth of bubble rise, microplasticterm·dom size and concentration, equations for jet dropterm·dom capture of microplasticterm·dom are developed. Subsequently, we integrate our findings on the individual transport mechanism into a global estimation of microplasticterm·dom emissionsterm·key, by considering a physics-based sea sprayterm·key generation function for jet dropsterm·dom which is a function of wind and waves at the ocean surface together with estimates of the ocean microplasticterm·dom concentration ( 8 ). Finally, an estimate of the global emissionterm·key of microplasticterm·dom from the ocean for multiple possible microplasticterm·dom concentrations is discussed. We argue that remaining uncertainties reside predominantly in the limited knowledge of ocean microplasticterm·dom concentration maps.
Sketch of the relevant processes ejecting microplasticterm·dom out of the ocean adapted from Refs. ( 4 , 16 ). A) Microplasticterm·dom (red/darker color) in the ocean is transported into the atmosphere by sea sprayterm·key drops. B) Bursting bubbles create small drops or aerosols such as jet dropsterm·dom. C) Microplasticsterm·dom present in the liquid can be carried up into the jet dropsterm·dom produced. The arrows point to 100 μ m microplasticterm·dom pieces. Drops produced can be picked up by wind and carry microplasticterm·dom material up into the atmosphere. The liquid eventually evaporates leaving behind microplasticterm·dom pieces. D) The relevant physical processes for bubble-bursting ejection of microplasticterm·dom start with the scavengingterm·key of particles as the bubble rises Di). After arriving at the surface Dii), the bubble eventually settles into its equilibrium shape which—upon bursting—focuses capillary waves at its base to form jet dropsterm·dom Diii) which carry microplasticterm·dom material.
Single bubble microplasticterm·dom ejection We observe the transport of microplasticsterm·dom by individual bubble burstingterm·key through laboratory experimentsterm·mth using high-speed photography as shown in Fig. 2 : the top row is an above-surface view, while the two lower sequences (A and B) show an underwater view of two different bursting events at the same conditions (liquid, bubble size, particle size, and particle number concentration). The red arrows in each sequence indicate microplasticterm·dom particles of interest. The above-water sequence shows the formation of a drop from the jet which carries a clearly visible microplasticterm·dom particle. In sequences A and B, the highlighted particles are captured even though the particle in B is not at the bottom center of the cavity. Throughout the entire study, no microplasticterm·dom particles are observed in the thin-film cap which is O ( 0.1 − 10 μ m ) thick and smaller than the size of particles studied.
An R b = 1.9 mm bubble burstingterm·key in deionized water ( R b / l c = 0.70 and R b / l μ = 150,000 ) to produce a jet dropterm·dom, r d = 380 μ m , which can transport 100 μ m diameter polyethylene ( ρ MP = 1.00 g / cm 3 ) microplasticterm·dom pieces. There are χ = 3.38 × 10 8 m − 3 pieces of plastic per unit volume of liquid in the bulk surrounding the bubble. In this regime, microplasticterm·dom pieces are able to be captured and transported by the jet dropterm·dom, as shown in sequences A and B.
Size of jet dropsterm·dom carrying microplasticterm·dom The radius, r d , (and velocity) of the first jet dropterm·dom produced by a bursting bubble has been shown to be controlled by the ratio of the bubble radius, R b , to the visco-capillary length scale, l μ = μ l 2 / ( ρ l σ ) ( 45 , 46 ), where μ l is the liquid viscosity, ρ l is the liquid density, and σ is the surface tension. Measurements of the jet dropterm·dom radius as a function of bubble size from multiple prior (experimental and numerical) studies ( 30 , 34 , 46–50 ) as well as proposed relationships ( 46 , 47 , 51 ) (see formulas in online supplementary material ) are shown in Fig. 3 . The radii of jet dropsterm·dom carrying microplasticterm·dom in the present study are shown in Fig. 3 as solid circles. They agree well with existing data which indicates that the established relations for r d can be applied to the microplasticterm·dom transport process.
Jet dropterm·dom size, r d , as a function of bubble size, R b , both nondimensionalized by the visco-capillary length scale l μ = μ l 2 / ( ρ l σ ) . The results from our experimentsterm·mth where microplasticterm·dom is transported by jet dropsterm·dom (large solid circles) are compared to multiple experimental and numerical studies measuring the size of jet dropsterm·dom ( 30 , 34 , 46–50 ). The shape of a surface bubble changes with size, R b / l c , as shown by the inset. The shape of the largest and smallest bubbles in this study are shown in red (darker color) and orange (lighter color), respectively.
The ratio of the bubble size to the capillary length, l c = σ / ( Δ ρ g ) , (where Δ ρ is the density difference between the two fluids and g is gravity) determines the bubble shape. The range of bubble sizes in this study is shown in the inset of Fig. 3 , and size is a secondary parameter controlling jet dropterm·dom existence and the speed at which they are ejected ( 50 ) (see online supplementary material ). The inset of Fig. 3 shows the outline of the largest and smallest bubbles in this study: R b / l c = 0.70 (red) and R b / l c = 0.31 (orange).
Bursting at the surface Having determined the jet dropterm·dom radius, we measure the number of microplasticterm·dom particles captured by each bursting bubble. This is obtained by catching each individual drop on a flat plate suspended above the free surface. The drops are then dried and any particles left behind are counted in a microscope. We measure about 50–100 bubble-bursting events at each condition to account for statistical variability in the system (see Tables S3 and S4 in the online supplementary material ).
Microplasticterm·dom ejection is first studied by releasing bubbles at a fixed height, H = 1.5 cm, below the free surface. At this fixed value of small H (comparable to the bubble size), the number of particles ejected, N , is a function of the particle size, r MP , the particle concentration in the liquid, χ, and the drop size, r d , which is determined by R b / l μ (with l μ being the visco-capillary length scale accounting for the properties of the liquid, i.e. viscosity, surface tension, and density following Fig. 3 ). The Stokes number defined as St = ( ρ MP − ρ l ) r MP 2 u b / ( 9 μ l R b ) with ρ MP being the particle density and u b the bubble’s rise velocity, describes the particle’s response time to drag forces versus the characteristic timescale of the flow. For all conditions, St ≪ 1 (in fact, for some conditions ρ MP = ρ l = 1.00 g / cm 3 such that St ≈ 0 ). As such, two dimensionless groups are expected to control the number of microplasticterm·dom particles ejected: Figure 4 A shows the number of particles per drop, N , versus the first dimensionless group, r MP / r d at constant H = 1.5 cm. The liquids used include deionized (DI) water, ethanol–water mixtures, and salt water with various microplasticterm·dom sizes and bubble sizes. The data are colored by χ 4 π r d 3 / 3 , the second dimensionless group. The average number of particles transported ranges from about 0.1 to 300. While O ( 10 2 ) particles can be captured for r MP / r d ≪ 1 , the vertical extent of the data with decreasing χ (and the vertical striation by color) indicates that N is a stronger function of χ 4 π r d 3 / 3 . r MP / r d describes a cutoff as r MP / r d ≈ 1 above which no particles are transported via jet dropsterm·dom: the size of the jet dropterm·dom being the maximum particle size that can be transported [consistent with results of other studies ( 7 )]. Bubbles are created in salt water, DI water and ethanol–water mixtures to test the influence of the liquid properties on r d and the number of microplasticsterm·dom particles emitted, and we demonstrate a universal behavior in the emissionsterm·key of microplasticterm·dom by a bursting bubble jet dropterm·dom (Figs. 3 and 4 ).
Number of microplasticterm·dom particles transported by a jet dropterm·dom, N . Data for deionized water, salt water (at ocean salinity), and 20% ethanol are plotted as circles, crosses, and diamonds, respectively, and the error bars are the standard error of N . A) Data at H = 1.5 cm for N versus r MP / r d color coded by χ 4 3 π r d 3 , the bulk liquid’s number concentration times drop volume. The gray dashed line represents a theoretical maximum number of particles in the drop. B) Same data with N as a function of χ 4 3 π r d 3 with r MP / r d color coded. The solid black line is Eq. 2 with a prefactor of 11 found by best fit. At a constant rise height, the number of microplasticsterm·dom per bubble is well described by the particle concentration and jet dropterm·dom size. C) N as a function of ( H / R b ) . Grouping the data by consistent color (which denotes χ 4 π r d 3 / 3 of each point) shows an N ∝ ( H / R b ) 1 / 2 relationship. The dashed line shows the linear scaling from the interceptor model ( 41 ) and does not describe the data with large H / R b values, suggesting that as a bubble rises it becomes increasingly less able to scavenge particles, consistent with the transition from free-slip to nonslip on the bubble’s surface as it collects surfactant. D) Scaling of N accounting for the bursting physics and particle scavengingterm·key as a function of ( H / R b ) 1 / 2 χ 4 3 π r d 3 ( Eq. 3 ), with r MP / r d color coded. The solid line is a fit with dimensionless prefactor of 4.5. The inset shows E = N / ( χ 4 3 π r d 3 ) as a function of ( H / R b ) 1 / 2 as well as the data from ( 35 ) [reprinted with permission from Ref. ( 35 ). Copyright 2023 American Chemical Society]. As in C, each set of experimentsterm·mth varying only H (data with a consistent color) shows the E ∝ ( H / R b ) 1 / 2 scaling.
The maximum possible number of particles that can be transported is found by replacing the total volume of liquid in the jet dropterm·dom by densely packed spheres of size r MP , yielding N max = π 3 2 ( r MP r d ) − 3 , where π / ( 3 2 ) is Carl Gauss’s dense-packing coefficient for spheres ( 52 ), shown by the gray dashed line in Fig. 4 A, and all of the data lies well beneath it.
Assuming that the particles are well mixed, any volume of liquid, V , contains an average of χ V randomly dispersed particles ( 53 ). As such, the expected number of microplasticterm·dom pieces per jet dropterm·dom can be written as where E is defined as an efficiency or enrichment factorterm·mth ( 35 , 36 ). It physically represents how much more concentrated the particles are in the jet dropterm·dom compared to the bulk liquid. Figure 4 B shows N as a function of χ 4 3 π r d 3 for the same data as in A, with the data color coded by r MP / r d . The data are shown to be well described by Eq. 2 as the black line shows. At the constant value of H = 1.5 cm, the efficiency factorterm·key is well approximated as a constant value of E = E c ≈ 11 obtained by least square fit. The strong effect of changes in χ can be seen for the data with the smallest values of r MP / r d where a large range of χ spreads the data over almost all values of N observed.
Particle scavengingterm·key After analyzing microplasticterm·dom transport at a constant rise height, bubbles are released from increasing depth (from 1.5 to 80 cm) to study their scavengingterm·key of particles. The effect of increasing H is shown in Fig. 4 C. Four different particle-drop-bubble size combinations are shown (see colorbar). Only the needle depth was changed across each set of points with a consistent color, and the vertical spread of the data is due to variations in r d and χ. Figure 4 C shows that the number of particles N increases with H and can be described by a ( H / R b ) 1 / 2 scaling.
The dashed line represents the N ∝ H / R b scaling from the interceptor model ( 41 ) which assumes that the surface is completely mobile and captures particles at a constant rate. Nonconstant rate of particle collection as the bubble rises has been previously reported ( 18 , 43 , 44 ). The N ∝ ( H / R b ) 1 / 2 scaling suggests that when considering large variations of H / R b , the bubble’s surface looses its free-slip condition as it scavenges microplasticterm·dom particles and surface-active molecules ( 43 ), transitioning to a nonslip condition ( 54–57 ). Indeed the collection rate for a sphere with a nonslip surfaceterm·dom is always less than that of one with a free-slip surfaceterm·dom for a fixed particle and bubble size ( 43 , 58 ), which provides a rationale for the observed N ∝ ( H / R b ) 1 / 2 scaling.
Finally, Fig. 4 D shows the final scaling for N accounting for both the scavengingterm·key ( E = E 0 ( H / R b ) 1 / 2 ) and surface-bursting ( Eq. 2 ), leading to where E 0 is a dimensionless prefactor found to be E 0 ≈ 4.5 by best fit shown by the solid line. The inset shows E = N / ( χ π r d 3 4 / 3 ) as a function of ( H / R b ) 1 / 2 and includes data from Ref. ( 35 ). The E ∝ ( H / R b ) 1 / 2 scaling is clear for groups of consistent colors. While correction factors for E have been proposed accounting for shape variation of the particles ( 35 ) and transfer from the bubble to the jet dropterm·dom ( 36 ), we show that the the average number of microplasticsterm·dom per jet dropterm·dom can be described by a simple model only accounting for the bursting and scavengingterm·key processes.
Modeling global emissionterm·key of oceanic microplasticterm·dom Having elucidated and quantified the microplasticterm·dom emissionsterm·key by individual bubble burstingterm·key, we now aim to quantify microplasticsterm·dom emissionsterm·key at the ocean surface. For drop production at the ocean surface, we leverage the mechanistic approach proposed by Ref. ( 20 ) that specifically introduces a jet drop emission functionterm·mth F d ( r d ) , coherent with sea sprayterm·key aerosols field observations and previously proposed sea sprayterm·key generation functions ( 17 , 19 , 21 , 60 , 61 ). The jet drop emission functionterm·mth provides the size distribution per unit ocean area per unit time and is sensitive to wind, waves and sea surface temperature ( 20 ). With data of the microplasticterm·dom concentration at the surface of the ocean, α MP (mass concentration per unit ocean surface area) ( 8 ), the oceanic microplasticterm·dom emissionterm·key by bubble burstingterm·key can be written as which has units of mass of microplasticterm·dom per unit area of ocean per unit time. The ratio α MP / z MP represents the mass concentration per unit volume over a well mixed microplasticterm·dom layer of depth z MP . As shown in Fig. 4 D, E = E 0 ( H / R b ) 1 / 2 , and since r d is a unique function of R b for a given l μ ( 47 ) (see Fig. 3 ), the efficiency factorterm·key can be represented as a function of H and r d : E ( H , R b ) ≡ E ( H , r d ) = E 1 [ H / r d ( R b / l μ ) ] 1 / 2 , with E 1 ≈ 1.7 fitted to the data (Fig. S7 in the online supplementary material ).
Bubbles on the ocean are created by breaking waves as air is entrained under the surface. The characteristic depth of bubble entrainment has been found to scale with the significant wave height, H s ( 62–65 ) so that we consider the height over which the bubbles are rising H = H s / 2 . While there are multiple mechanisms which disperse microplasticterm·dom in the ocean ( 66 ), the turbulent velocity field created by breaking waves—which entrain air to create bubbles—extends with the significant wave height ( 67 ), so we consider that the microplasticterm·dom layer can be written as z MP = H s / 2 . This is consistent with microplasticterm·dom transport models which have a vertical length scale proportional to wave height ( 11 , 68 , 69 ). With these assumptions on H and z MP , Eq. 4 becomes where ζ ≈ 2.4 is the combined prefactor from the efficiency function, depth of entrainment and microplasticterm·dom layer. Because the scavenging rateterm·dom is found to decrease with bubble rise heightterm·dom, the total amount of microplasticterm·dom ejected depends on the ocean conditions through both F d and the bubble scavengingterm·key.
The number and mass of microplasticterm·dom particles emitted can be evaluated using the sea sprayterm·key generation function ( 20 ), ocean-surface microplasticterm·dom concentration ( 8 ), and estimates of the global sea state. Figure 5 A shows the annual mean of the jet drop emission functionterm·mth [expressed in drop volume emitted per year per ocean surface area, computed for a representative year (2014) using realistic wind forcing, see online supplementary material and Ref. ( 20 )] with strong production at high latitude corresponding to high winds. Figure 5 B shows the concentration of “ejectable” microplasticterm·dom estimated from Ref. ( 59 ). Particles ranging in size from 100 μ m < 2 r MP < 800 μ m within 5 m of the ocean’s surface in 2020 are shown (expressed in mass per ocean surface area). Figure 5 C shows the total amount of microplasticterm·dom ejected out of the ocean accounting for the wind and waves through the sea sprayterm·key generation function (Figure 5 A) and the microplasticterm·dom concentration map (Figure 5 B) shown to produce an annual ocean microplasticterm·dom emissionterm·key map. High emissionsterm·key near the coast are the result of the high concentration of microplasticterm·dom there. The annual emissionterm·key is then obtained by global integration. Considering various possible microplasticterm·dom concentration maps leads to an annual emissionterm·key range of 0.02 to 7.4 Mt/yr with a best guess estimate of 0.1 Mt/yr.
Global microplasticterm·dom ejection from the ocean. Terms in Eq. 4 are illustrated, as yearly average. A) Volume of liquid contained in jet dropsterm·dom ejected from the ocean by bursting bubbles per unit area per unit time in the ocean, constrained by surface wind and waves following Ref. ( 20 ), for the year of 2014. The strong production at high latitude corresponds to high winds and waves. B) Microplasticterm·dom mass concentration within 5 m of the ocean surface and ranging in size from 100 μ m < 2 r MP < 800 μ m in 2020 from Ref. ( 59 ) C) Microplasticterm·dom emissionsterm·key following Eq. 5 using the microplasticterm·dom distribution shown in panel B. Integrating over the globe, the emissionterm·key is ≈ 0.1 Mt/yr.
Caveats and limitations As shown in Fig. 4 , Eq. 3 describes well the amount of microplasticterm·dom ejected by jet dropsterm·dom across a wide range of microplasticterm·dom concentrations, bubble sizes, and rise distance. The uncertainties in the bursting physics is well constrained in the laboratory as shown in Fig. 4 , while upscaling uncertainties are estimated to be within a factor of 2 to 4 even considering the complex multi-scale processes ( 20 ). However, the dominant uncertainty of our global emissionterm·key estimate is the limited knowledge for ocean microplasticterm·dom concentration; it is both scarce in temporal and spatial coverage and has a coarse size resolution ( 3 , 8 ). Furthermore, they are interpolated from limited data sets across multiple years standardized in time to what would be observed in 2014 from Ref. ( 8 ), or 2020 in Ref. ( 59 ). Table 1 outlines the sensitivity of the global emissionterm·key to the choice of microplasticterm·dom concentration coverage map. An upper-limit of concentration is provided by Ref. ( 8 ) which bins together particles sizing from 330 μ m to 200 mm in diameter. The average particle radius in any given area however is 1 ± 0.5 mm , which indicates that a considerable portion of particles are of an “ejectable” size, and with it the annual emissionterm·key rate for 2014 is found to be 7.4 Mt. To estimate only the ejectable amount of microplasticterm·dom from this dataset, a size distribution P ( r MP ) must be assumed (see methods) and be limited to 1 mm.
Sensitivity of global ocean microplasticterm·dom emissionterm·key calculated by Eq. 5 to the microplasticterm·dom concentration map employed ( α MP ).
We test two different size distributions, an exponential and a piecewise distribution which both preserve the total number for every datum of the original dataset from Ref. ( 8 ) integrating over the full range of sizes (see methods). Once determined, the distributions are then integrated within the range of ejectable microplasticsterm·dom from 330 μ m < 2 r MP < 1 mm to describe the amount of ejectable microplasticsterm·dom. The global emissionterm·key assuming the exponential and piecewise size distributions are shown in Table 1 to be 0.06 Mt/yr and 0.1 Mt/yr respectively; demonstrating a high sensitivity to the choice of distribution and size cutoff of microplasticsterm·dom. We note that in all these calculations, the drop diameter being emitted is the same, going up to 1 mm, and argue that the large drops ejected are essential to accurately estimating the mass of microplasticterm·dom emissionsterm·key. Using the recent dataset from Ref. ( 59 ) and an upper size of 800 microns leads to emissionterm·key of 0.1 Mt/yr.
Some works calculating microplasticterm·dom emissionsterm·key chose to assume a constant microplasticterm·dom concentration across all oceans ( 70 ). Others ( 7 ) employed a concentration distribution from Ref. ( 5 ) which features a maximum particle size of 70 μ m and was calibrated with Ref. ( 8 ). Using the coverage proposed by Ref. ( 5 ) in Eq. 4 , we find a microplasticterm·dom emissionterm·key rate of 0.02 Mt/yr, significantly lower due to the maximum size of 70 μ m considered by Ref. ( 5 ) which is well below the maximum particle size able to be transported by jet dropsterm·dom ( ≈ 500 μ m ). The sensitivity of the global emissionterm·key to the concentration map motivates improved concentration and size-distribution measurements of microplasticterm·dom ranging in size from 50 μ m < r MP < 500 μ m .
While we demonstrate in this work that microplasticterm·dom can be ejected if it is the same size or smaller than the jet dropsterm·dom (which have a maximum size of 2 r MP ≈ 1 mm for sea water), our model does not describe their subsequent transport into the turbulent atmospheric boundary layer. Even without microplasticterm·dom, a large volume of liquid droplets redeposit to the ocean after being ejected to the region just above the surface: see Fig. 5 of Ref. ( 19 ). As such, our model describes the emissionterm·key of microplasticterm·dom from the ocean surface, not the flux that is actually transported to the upper atmosphere.
Another limitation of the global emissionterm·key function is that the microplasticterm·dom concentration data used is a constant annualized mean for each spatial position. The emissionterm·key function (and sea state) however are calculated on a much shorter timescale of O (hours). As such, nonlinearities between temporal variations in the microplasticterm·dom concentration and jet dropterm·dom production flux could affect the final result due to this difference in timescale sampling.
Conclusion The amount of microplasticterm·dom ejected from the ocean is determined by studying individual drops created by bursting bubbles to create a global emissionterm·key model which depends on local sea state. We demonstrate experimentally by direct high-speed video visualization and measurement that microplasticterm·dom particles (with diameters from 10 to 280 μ m ) are transported out of a bulk liquid such as the ocean and into the air above by jet dropsterm·dom resulting from bursting bubbles. The size of jet dropsterm·dom containing microplasticsterm·dom is found to be well described by theoretical scaling laws derived for drops without microplasticsterm·dom. The number concentration of microplasticterm·dom particles in emitted drops is found to be proportional to the concentration of particles in the bulk and jet dropterm·dom volume, with an efficiency or enrichment factorterm·mth E = E 0 ( H / R b ) 1 / 2 (we find E 0 ≈ 4.5 ), suggesting a reduction in surface mobility as the bubbles rise over large H / R b . The number of emitted microplasticterm·dom particles for a single bubble burstingterm·key is thus given by N = E 0 ( H / R b ) 1 / 2 χ 4 π 3 r d 3 .
Given that jet dropsterm·dom produced by bubble burstingterm·key at the ocean surface dominate sea sprayterm·key aerosols emissionsterm·key of size 2 μ m to 1 mm ( 17 , 20 ) and that microplasticterm·dom pieces are carried by jet dropsterm·dom equal to or larger in size, jet dropsterm·dom can be effective at emitting microplasticterm·dom up to O (1 mm) in size and are hence responsible for most of the emitted mass of microplasticterm·dom from the ocean to the atmosphere. To estimate the global microplasticterm·dom emissionsterm·key by the ocean, the individual emissionterm·key efficiency is combined with a jet drop emission functionterm·mth [from Ref. ( 20 )] and microplasticterm·dom concentration in the ocean [from Ref. ( 8 ), see the Materials and methods section for a discussion of particle size compared to that of the emissionterm·key function]. Thus, we provide an independent bottom-up estimate of microplasticterm·dom emissionsterm·key from the ocean that can be compared to previous estimates obtained from atmospheric observations and inverse modeling ( 3 , 5 ) and other independent estimates which did not consider any efficiency factorterm·key or scavengingterm·key dynamics ( 7 ). We obtain a range of annual microplasticterm·dom ocean emissionterm·key from 0.02 to 7.4 Mt/yr with the upper bound accounting for sizes larger than what can be ejected by jet dropsterm·dom. The upper bound is similar in magnitude to the best guess value of 8.6 Mt/yr from Ref. ( 5 ). Our best guess using the most recent concentration data set from Ref. ( 59 ) gives an estimate of 0.1 Mt/yr. The emissionterm·key employing the microplasticterm·dom maps of Ref. ( 5 ) are much larger than those from Refs. ( 7 , 70 ) as we account for larger droplets being emitted. While the upper bound of our emissionterm·key estimate supports the hypothesis that ocean emissionterm·key of microplasticterm·dom plays a significant role in global-scale microplasticterm·dom transport as proposed by Ref. ( 5 ), the precise emissionterm·key remains difficult to estimate better without improved concentration maps that include accurate size distributions.
Our study provides an independent estimate of the microplasticterm·dom emissionterm·key by the ocean and reduces the associated uncertainties by identifying and quantifying the leading emissionterm·key mechanism. As a consequence, the remaining large part of the uncertainty comes from the limited oceanic observations of microplasticterm·dom concentration and the lack of high resolution size distributions of microplasticterm·dom pieces present at the surface of the ocean (particularly at scales of less than 50 μ m ) confirming previous assessments ( 3 , 5 ). A significant effort in observational data of marine microplasticterm·dom data is necessary to reduce these present uncertainties.
Materials and methods Experimental methods We analyze the ejection dynamics visually, using two high-speed cameras, a phantom V2012 and a phantom 4 K, to capture videos (up to 22,000 frames per second) of a single bubble burstingterm·key. An example of bubble burstingterm·key in water is shown in Fig. 2 . At each experimental condition, approximately 10 videos are taken of the bursting event from both above and below the free surface.
The radius of the bubble, R b , is varied by employing different sized needles; inner diameters ranging from 0.4 to 2.7 mm to produce bubbles of size R b = 0.88 to 1.90 mm by slowly pushing air through the needle with a syringe pump so that only a single bubble is created at a time. Two different liquid containers were employed. For the H = 1.5 cm conditions, a shallow dish of size 2 × 10 × 10 cm was used to hold the liquid, while the H > 1.5 cm trials were conducted in a tank of size 80 × 20 × 20 cm . In both containers, the bubbles did not touch the bottom or sides during the rising or bursting process.
In each trial, spherical microplasticterm·dom particles ranging in radius, r MP , from 5 to 140 μ m were added to the liquid. All particles were polyethylene with a density of 1.00 g / cm 3 except for one set of trials with r MP = 5 μ m hollow borosilicate glass with a density of 1.1 g / cm 3 . We have considered microplasticterm·dom of near neutral density compared to the liquid, with slight buoyancy variations, the density ratio ρ MP / ρ l varying from 0.97 to 1.14 (see Table S2 of the online supplementary material ). We do not observe any significant change due to density in the transport over this range of density variation. The concentration of particles added ranges from χ = 1.02 × 10 8 to 1.166 × 10 11 or equivalently α = 0.1 gm − 3 to 1166 gm − 3 ; [ χ ] = L − 3 is a number concentration of pieces per unit volume of liquid while [ α ] = M L − 3 is the sum of microplasticterm·dom mass per unit volume of liquid. During measurements, the mixture is mixed routinely to ensure that χ is constant for the duration of each trial. The size of the tank is much larger (at least a factor of 10 6 ) than the volume of jet dropsterm·dom, so the drop ejection does not change χ over the duration of each experimentterm·mth.
While only spherical particles were used in this study, the emissionterm·key of other shapes has been studied previously. Some works find that emissionterm·key does not vary significantly between spheres and fibers ( 71 ), and others find that fibers are not ejected as well ( 7 ); particle shape could be an area of future research. This work does not consider the transport of particles at the ocean’s surface which is complex, and has been shown to be a function of both the particle’s shape and the flow field ( 11 , 72–74 ). This study used both polyethylene and borosilicate glass particles, and no significant difference was observed in the capture and emissionterm·key between the different materials, which is consistent with other work ( 7 ). Biofouling is another effect which was not considered in this work, and it has been shown to change the effective density of microplasticterm·dom particles ( 75 ). As the concentration measurements used in this model are taken in the ocean, they have already accounted for any effect biofouling has on the presence of particles at the ocean’s surface. A recent work ( 76 ) showed that there is a correlation between the presence of microplasticsterm·dom and surfactant concentrations in the ocean; this could be due to biofouling which would have the potential to alter the liquid properties (such as surface tension).
Multiple liquids were used in the present study: DI water, salt water with a salinity of 42 g/kg similar to sea water (“Sea Salt” ASTM D1141-98), and ethanol–water solution (20% by mass). A small amount of sodium dodecyl sulfateterm·dom (SDSterm·dom) is added to the liquid at a concentration of 6.6 μ m which is several orders of magnitude below the critical micelle concentration of SDSterm·dom of 8.2 mM . As a surface-active (surfactant) chemical, SDSterm·dom is able to provide surface forces and is a common additive to weakly contaminate the surface—making it more like the marine environment. In this concentration range, the difference in surface tension due to the addition of surfactant is within measurement error using the pendant drop method ( 77 ). Liquid properties (density, viscosity, and surface tension) are provided in Table S1 of the online supplementary material . Varying the liquid properties enables us to understand the role of the physico-chemical properties and derive universal scalings that can be applied to sea water or any other liquid.
Together with the direct images obtained by high-speed videos (see Fig. S1 of the online supplementary material for a diagram of the experimental setup), the jet dropsterm·dom produced are also captured to count the number of microplasticterm·dom particles they contain. A flat petri dish is positioned above the free surface: ranging from 1 to 6 cm above the surface so that the jet dropsterm·dom impact at a low enough velocity to stick. The drops are allowed to dry, leaving behind the microplasticterm·dom that was captured by the jet dropterm·dom. The deposited particles are counted using a Leica DMI4000 B microscope. For each condition, O ( 50 − − 100 ) drops are collected and the number of microplasticterm·dom particles in each drop are counted and their average is reported. Examples of the microscope view of particles on the collection plate are shown in Figs. S2 and S3 of the online supplementary material , and examples of microplasticterm·dom number distributions are shown in Figs. S4–S6 of the online supplementary material .
Global model parameters and uncertainties Number and mass coverage of microplasticsterm·dom have been proposed with interpolation models varying by a factor of 2 (shown in Figs. S8–S11 of the online supplementary material ) ( 5 , 8 ), but the uncertainty on the data is probably larger and related to the binning of most microplasticterm·dom size measurements. The concentration map of particles sized 100 μ m < 2 r MP < 800 μ m located within 5 m of the ocean’s surface in 2020 from a recent study ( 59 ) is also shown in Fig. S12 of the online supplementary material . The uncertainty of the resulting output depending on mass coverage map is shown in Figs. S13–S17 of the online supplementary material . Furthermore, the limited measurements in the ocean are taken in a variety of conditions which affects their spatial distribution in the ocean ( 9–11 ). In Ref. ( 8 ), the bin size ranged from 0.15 mm < r MP < 100 mm, and over 90% of the data was collected with a 0.33 mm net mesh. As such, it is treated as an upper bound of microplasticterm·dom concentration. Other studies [such as Refs. ( 5 , 78 )] suggest a somewhat smaller number could also be appropriate. Estimates of ejectable microplasticterm·dom are determined from Ref. ( 8 ) by assuming a particle size distribution. As such, the total number of particles is ∫ min ( r MP ) max ( r MP ) P ( r MP ) d r MP and the total mass is ∫ min ( r MP ) max ( r MP ) P ( r MP ) 4 π 3 ρ MP r MP 3 d r MP for any given latitude and longitude. The full range of sizes from Ref. ( 8 ) goes from 330 μ m < 2 r MP < 200 mm . We consider two options for the distribution in order to account for the fact that particles above 1 mm are not ejected by jet dropsterm·dom. The distributions we consider are an exponential P ( r MP ) = A e − r MP / ⟨ r MP ⟩ (where ⟨ r MP ⟩ is the average size for each datum) and a piecewise distribution, P ( r MP ) = { B for r MP < ⟨ r MP ⟩ and C for r MP > ⟨ r MP ⟩ } , where A , B , and C are determined at each point. Concentration maps of each assumption are shown in Figs. S18 and S19 of the online supplementary material , respectively. The piecewise distribution maintains both the mass and number for each datum. Once determined, the distributions are then integrated from 330 μ m < 2 r MP < 1 mm to describe the amount of ejectable microplasticterm·dom, and summarized in Table 1 . The ejection rate for each size distribution is shown in Fig. S20 and S21 of the online supplementary material , respectively. A comparison between the concentration and subsequent emissionterm·key assuming a piecewise distribution of the data in ( 8 ) and the most recently published concentration dataset in Ref. ( 59 ) is shown in Fig. S22 of the online supplementary material . While the assumptions that produce each are different, the estimated concentration of ejectable microplasticterm·dom and the subsequent emissionterm·key maps are strikingly similar and the total annual emissionterm·key of each is 0.1 Mt/yr.
Global emissionterm·key is also computed with the dataset from ( 5 ) which is treated as a lower bound due to its relatively small maximum particle size of 70 μ m . Measurements ( 79 ) of oceanic plastic with detailed size distinctions below 330 μ m found that particle number concentration was approximately constant for all size bins between 50 and > 1 , 000 μ m , and as much as 50% of the number of microplasticterm·dom particles in the ocean are smaller than 50 μ m . Another study ( 74 ) focusing on microplasticterm·dom concentration from the surface to 10 m of depth found that the total amount of microplasticterm·dom near the ocean’s surface may be 10 times larger than what would be estimated by a surface net alone as most ocean-surface microplasticterm·dom studies use. While precise (both in spatial resolution and the size-binning resolution) data for the size-range necessary are lacking, we have estimated the amount of ejectable microplasticterm·dom from Ref. ( 8 ) and find a range spanning about two orders of magnitude that depends solely on the microplasticterm·dom concentration coverage chosen.
Microplasticterm·dom concentration data are reported as a number per unit ocean surface area. The raw data was collected at the ocean’s surface and standardized to account for mixing due to ocean conditions ( 8 ). The depth over which this microplasticterm·dom is spread, z MP needs to be estimated in our emissionterm·key model. Following measurements and modeling of ocean-surface microplasticterm·dom transport, we assume that the characteristic depth scales with significant wave height to obtain the volumetric concentration α MP / z MP . However, the relationship between z MP and ocean conditions remains an open question as there are very limited observations of underwater microplasticterm·dom concentration profiles ( 66 , 68 , 69 , 79 ), and there is some degree of uncertainty with this estimate. The prefactor by which z MP and H scale with significant wave height is estimated to range from 0.25 to 0.75 . The uncertainty of the prefactor of the jet drop emission functionterm·mth is estimated to be a factor of 2 in Ref. ( 20 ). These uncertainties of the emissionterm·key process are estimated to be less significant than the uncertainties of the microplasticterm·dom concentration data.
Supplementary Material Click here for additional data file.