skip to main content

Title: Sinking Versus Suspended Particle Size Distributions in the North Pacific Subtropical Gyre

The particle size distribution (PSD) is a fundamental property that influences all aspects of phytoplankton ecology. In particular, the size (e.g., diameterd[μm]) and sinking speedw(m/day) of individual particles are inextricable, but much remains unknown about howdandware related quantitatively for bulk particulate matter. There is significant interest in inferring sinking mass fluxes from PSDs, but doing so requires knowing how both mass andwscale withd. To this end, using both laser diffraction and imaging, we characterized for the first time both sinking and suspended PSDs in the oligotrophic North Pacific subtropical gyre. Comparing these PSDs via a power law parameterization indicates an approximately linearw‐to‐dscaling, suggesting particles are more fractal‐like than sphere‐like in this respect. This result is robust across multiple instruments, depths, and sediment trap deployments and is made comparatively precise by a high degree of replication.

more » « less
Author(s) / Creator(s):
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The abundance and size distribution of marine particles control a range of biogeochemical and ecological processes in the ocean, including carbon sequestration. These quantities are the result of complex physical‐biological interactions that are difficult to observe, and their spatial and temporal patterns remain uncertain. Here, we present a novel analysis of particle size distributions (PSDs) from a global compilation of in situ Underwater Vision Profiler 5 (UVP5) optical measurements. Using a machine learning algorithm, we extrapolate sparse UVP5 observations to the global ocean from well‐sampled oceanographic variables. We reconstruct global maps of PSD parameters (biovolume [BV] and slope) for particles at the base of the euphotic zone. These reconstructions reveal consistent global patterns, with high chlorophyll regions generally characterized by high particle BV and flatter PSD slope, that is, a high relative abundance of large versus small particles. The resulting negative correlations between particle BV and slope further suggests synergistic effects on size‐dependent processes such as sinking particle fluxes. Our approach and estimates provide a baseline for an improved understanding of particle cycles in the ocean, and pave the way to global, three‐dimensional reconstructions of PSD and sinking particle fluxes from the growing body of UVP5 observations.

    more » « less
  2. Abstract

    The sinking of organic particles produced in the upper sunlit layers of the ocean forms an important limb of the oceanic biological pump, which impacts the sequestration of carbon and resupply of nutrients in the mesopelagic ocean. Particles raining out from the upper ocean undergo remineralization by bacteria colonized on their surface and interior, leading to an attenuation in the sinking flux of organic matter with depth. Here, we formulate a mechanistic model for the depth-dependent, sinking, particulate mass flux constituted by a range of sinking, remineralizing particles. Like previous studies, we find that the model does not achieve the characteristic ‘Martin curve’ flux profile with a single type of particle, but instead requires a distribution of particle sizes and/or properties. We consider various functional forms of remineralization appropriate for solid/compact particles, and aggregates with an anoxic or oxic interior. We explore the sensitivity of the shape of the flux vs. depth profile to the choice of remineralization function, relative particle density, particle size distribution, and water column density stratification, and find that neither a power-law nor exponential function provides a definitively superior fit to the modeled profiles. The profiles are also sensitive to the time history of the particle source. Varying surface particle size distribution (via the slope of the particle number spectrum) over 3 days to represent a transient phytoplankton bloom results in transient subsurface maxima or pulses in the sinking mass flux. This work contributes to a growing body of mechanistic export flux models that offer scope to incorporate underlying dynamical and biological processes into global carbon cycle models.

    more » « less
  3. Abstract

    Total ice water content (IWC) derived from an isokinetic evaporator probe and ice crystal particle size distributions (PSDs) measured by a two-dimensional stereo probe and precipitation imaging probe installed on an aircraft during the 2014 European High Altitude Ice Crystals–North American High IWC field campaign (HAIC/HIWC) were used to characterize regions of high IWC consisting mainly of small ice crystals (HIWC_S) with IWC ≥ 1.0 g m−3and median mass diameter (MMD) < 0.5 mm. A novel fitting routine developed to automatically determine whether a unimodal, bimodal, or trimodal gamma distribution best fits a PSD was used to compare characteristics of HIWC_S and other PSDs (e.g., multimodality, gamma fit parameters) for HIWC_S simulations. The variation of these characteristics and bulk properties (MMD, IWC) was regressed with temperature, IWC, and vertical velocity. HIWC_S regions were most pronounced in updraft cores. The three modes of the PSD reveal different dominant processes contributing to ice growth: nucleation for maximum dimensionD< 0.15 mm, diffusion for 0.15 <D< 1.0 mm, and aggregation forD> 1.0 mm. The frequency of trimodal distributions increased with temperature. The volumes of equally plausible parameters derived in the phase space of gamma fit parameters increased with temperature for unimodal distributions and, for temperatures less than −27°C, for multimodal distributions. Bimodal distributions with 0.4 mm in the larger mode were most common in updraft cores and HIWC_S regions; bimodal distributions with 0.4 mm in the smaller mode were least common in convective cores.

    more » « less
  4. Abstract

    We compile full ocean‐depth size‐fractionated (1–51 and >51 μm) particle concentration and composition of suspended particulate matter from three recent U.S. GEOTRACES cruises, and exploit detailed information of particle characteristics measured to give insights into controls on sinking velocity and mass flux. Our model integrates the concept of fractal scaling into Stokes' Law by incorporating one of two porosity‐size power law relationships that result in fractal dimensions of 1.4 and 2.1. The medians of pump‐derived total (>1 μm) mass flux in the upper 100 m of gyre stations are 285.1, 609.2, and 99.3 mg/m2/d in the North Atlantic, Eastern Tropical South Pacific, and Western Arctic Ocean cruises, respectively. In this data set, variations in particle concentration were generally more important than sinking velocity in controlling variations in mass flux. We examine different terms in a Stokes' Law model to explore how variations in particle and water column characteristics from these three cruises affect mass flux. The decomposition of different aspects of the Stokes' relationship sheds light on the lowest total mass flux of the three cruises in the Western Arctic, which could be explained by the Arctic having the lowest particle concentrations as well as the lowest sinking velocities due to having the smallest particle sizes and the most viscous water. This work shows the importance of both particle characteristics and size distribution for mass fluxes, and similar methods can be applied to existing and future size‐fractionated filtered particulate measurements to improve our understanding of the biological pump elsewhere.

    more » « less
  5. Abstract

    This study derives simple analytical expressions for the theoretical height profiles of particle number concentrations (Nt) and mean volume diameters (Dm) during the steady-state balance of vapor growth and collision–coalescence with sedimentation. These equations are general for both rain and snow gamma size distributions with size-dependent power-law functions that dictate particle fall speeds and masses. For collision–coalescence only,Nt(Dm) decreases (increases) as an exponential function of the radar reflectivity difference between two height layers. For vapor deposition only,Dmincreases as a generalized power law of this reflectivity difference. Simultaneous vapor deposition and collision–coalescence under steady-state conditions with conservation of number, mass, and reflectivity fluxes lead to a coupled set of first-order, nonlinear ordinary differential equations forNtandDm. The solutions to these coupled equations are generalized power-law functions of heightzforDm(z) andNt(z) whereby each variable is related to one another with an exponent that is independent of collision–coalescence efficiency. Compared to observed profiles derived from descending in situ aircraft Lagrangian spiral profiles from the CRYSTAL-FACE field campaign, these analytical solutions can on average capture the height profiles ofNtandDmwithin 8% and 4% of observations, respectively. Steady-state model projections of radar retrievals aloft are shown to produce the correct rapid enhancement of surface snowfall compared to the lowest-available radar retrievals from 500 m MSL. Future studies can utilize these equations alongside radar measurements to estimateNtandDmbelow radar tilt elevations and to estimate uncertain microphysical parameters such as collision–coalescence efficiencies.

    Significance Statement

    While complex numerical models are often used to describe weather phenomenon, sometimes simple equations can instead provide equally good or comparable results. Thus, these simple equations can be used in place of more complicated models in certain situations and this replacement can allow for computationally efficient and elegant solutions. This study derives such simple equations in terms of exponential and power-law mathematical functions that describe how the average size and total number of snow or rain particles change at different atmospheric height levels due to growth from the vapor phase and aggregation (the sticking together) of these particles balanced with their fallout from clouds. We catalog these mathematical equations for different assumptions of particle characteristics and we then test these equations using spirally descending aircraft observations and ground-based measurements. Overall, we show that these mathematical equations, despite their simplicity, are capable of accurately describing the magnitude and shape of observed height and time series profiles of particle sizes and numbers. These equations can be used by researchers and forecasters along with radar measurements to improve the understanding of precipitation and the estimation of its properties.

    more » « less