Onset of coupled atmosphere-ocean oxygenation 2.3 billion years ago.

The initial rise of molecular oxygen (O2) shortly after the Archaean-Proterozoic transition 2.5 billion years ago was more complex than the single step-change once envisioned. Sulfur mass-independent fractionation records suggest that the rise of atmospheric O2 was oscillatory, with multiple returns to an anoxic state until perhaps 2.2 billion years ago1-3. Yet few constraints exist for contemporaneous marine oxygenation dynamics, precluding a holistic understanding of planetary oxygenation. Here we report thallium (Tl) isotope ratio and redox-sensitive element data for marine shales from the Transvaal Supergroup, South Africa. Synchronous with sulfur isotope evidence of atmospheric oxygenation in the same shales3, we found lower authigenic 205Tl/203Tl ratios indicative of widespread manganese oxide burial on an oxygenated seafloor and higher redox-sensitive element abundances consistent with expanded oxygenated waters. Both signatures disappear when the sulfur isotope data indicate a brief return to an anoxic atmospheric state. Our data connect recently identified atmospheric O2 dynamics on early Earth with the marine realm, marking an important turning point in Earth's redox history away from heterogeneous and highly localized 'oasis'-style oxygenation.

Neoarchaean oxygen-based nitrogen cycle en route to the Great Oxidation Event.

The nitrogen isotopic composition of sedimentary rocks (δ15N) can trace redox-dependent biological pathways and early Earth oxygenation1,2. However, there is no substantial change in the sedimentary δ15N record across the Great Oxidation Event about 2.45 billion years ago (Ga)3, a prominent redox change. This argues for a temporal decoupling between the emergence of the first oxygen-based oxidative pathways of the nitrogen cycle and the accumulation of atmospheric oxygen after 2.45 Ga (ref. 3). The transition between both states shows strongly positive δ15N values (10-50‰) in rocks deposited between 2.8 Ga and 2.6 Ga, but their origin and spatial extent remain uncertain4,5. Here we report strongly positive δ15N values (>30‰) in the 2.68-Gyr-old shallow to deep marine sedimentary deposit of the Serra Sul Formation6, Amazonian Craton, Brazil. Our findings are best explained by regionally variable extents of ammonium oxidation to N2 or N2O tied to a cryptic oxygen cycle, implying that oxygenic photosynthesis was operating at 2.7 Ga. Molecular oxygen production probably shifted the redox potential so that an intermediate N cycle based on ammonium oxidation developed before nitrate accumulation in surface waters. We propose to name this period, when strongly positive nitrogen isotopic compositions are superimposed on the usual range of Precambrian δ15N values, the Nitrogen Isotope Event. We suggest that it marks the earliest steps of the biogeochemical reorganizations that led to the Great Oxidation Event.

A climate threshold for ocean deoxygenation during the Early Cretaceous.

Oceanic anoxic events (OAEs) are historical intervals of global-scale ocean deoxygenation associated with hyperthermal climate states and biological crises1,2. Massive volcanic carbon dioxide (CO2) emissions frequently associated with these events are thought to be a common driver of ocean deoxygenation through several climate-warming-related mechanisms1,3,4. The Early Cretaceous OAE1a is one of the most intense ocean deoxygenation events, persisting for more than 1 Myr (refs. 5,6). However, existing records of marine chemistry and climate across OAE1a are insufficient to fully resolve the timing and dynamics of the underlying processes, thus obscuring cause-and-effect relationships between climate forcing and ocean oxygenation states. Here we show that rapid ocean deoxygenation during OAE1a is linked to volcanic CO2 emissions and the crossing of an associated climate threshold, after which the sluggish pace of the silicate-weathering feedback and climate recovery delayed reoxygenation for >1 Myr. At the end of OAE1a, recrossing this threshold allowed for ocean reoxygenation. Following OAE1a, however, the Earth system remained sufficiently warm such that orbitally forced climate dynamics led to continued cyclic ocean deoxygenation on approximately 100-kyr timescales for another 1 Myr. Our results thus imply a tight coupling between volcanism, weathering and ocean oxygen content that is characterized by a climate threshold.

Neogene burial of organic carbon in the global ocean.

Organic carbon buried in marine sediment serves as a net sink for atmospheric carbon dioxide and a source of oxygen1,2. The rate of organic carbon burial through geologic history is conventionally established by using the mass balance between inorganic and organic carbon, each with distinct carbon isotopic values (δ13C)3,4. This method is complicated by large uncertainties, however, and has not been tested with organic carbon accumulation data5,6. Here we report a 'bottom-up' approach for calculating the rate of organic carbon burial that is independent from mass balance calculations. We use data from 81 globally distributed sites to establish the history of organic carbon burial during the Neogene (roughly 23-3 Ma). Our results show larger spatiotemporal variability of organic carbon burial than previously estimated7-9. Globally, the burial rate is high towards the early Miocene and Pliocene and lowest during the mid-Miocene, with the latter period characterized by the lowest ratio of organic-to-carbonate burial rates. This is in contrast to earlier work that interpreted enriched carbonate 13C values of the mid-Miocene as massive organic carbon burial (that is, the Monterey Hypothesis)10,11. Suppressed organic carbon burial during the warm mid-Miocene is probably related to temperature-dependent bacterial degradation of organic matter12,13, suggesting that the organic carbon cycle acted as positive feedback of past global warming.

De novo evolution of macroscopic multicellularity.

While early multicellular lineages necessarily started out as relatively simple groups of cells, little is known about how they became Darwinian entities capable of sustained multicellular evolution1-3. Here we investigate this with a multicellularity long-term evolution experiment, selecting for larger group size in the snowflake yeast (Saccharomyces cerevisiae) model system. Given the historical importance of oxygen limitation4, our ongoing experiment consists of three metabolic treatments5-anaerobic, obligately aerobic and mixotrophic yeast. After 600 rounds of selection, snowflake yeast in the anaerobic treatment group evolved to be macroscopic, becoming around 2 × 104 times larger (approximately mm scale) and about 104-fold more biophysically tough, while retaining a clonal multicellular life cycle. This occurred through biophysical adaptation-evolution of increasingly elongate cells that initially reduced the strain of cellular packing and then facilitated branch entanglements that enabled groups of cells to stay together even after many cellular bonds fracture. By contrast, snowflake yeast competing for low oxygen5 remained microscopic, evolving to be only around sixfold larger, underscoring the critical role of oxygen levels in the evolution of multicellular size. Together, this research provides unique insights into an ongoing evolutionary transition in individuality, showing how simple groups of cells overcome fundamental biophysical limitations through gradual, yet sustained, multicellular evolution.

A well-oxygenated eastern tropical Pacific during the warm Miocene.

The oxygen content of the oceans is susceptible to climate change and has declined in recent decades1, with the largest effect in oxygen-deficient zones (ODZs)2, that is, mid-depth ocean regions with oxygen concentrations <5 µmol kg-1 (ref. 3). Earth-system-model simulations of climate warming predict that ODZs will expand until at least 2100. The response on timescales of hundreds to thousands of years, however, remains uncertain3-5. Here we investigate changes in the response of ocean oxygenation during the warmer-than-present Miocene Climatic Optimum (MCO; 17.0-14.8 million years ago (Ma)). Our planktic foraminifera I/Ca and δ15N data, palaeoceanographic proxies sensitive to ODZ extent and intensity, indicate that dissolved-oxygen concentrations in the eastern tropical Pacific (ETP) exceeded 100 µmol kg-1 during the MCO. Paired Mg/Ca-derived temperature data suggest that an ODZ developed in response to an increased west-to-east temperature gradient and shoaling of the ETP thermocline. Our records align with model simulations of data from recent decades to centuries6,7, suggesting that weaker equatorial Pacific trade winds during warm periods may lead to decreased upwelling in the ETP, causing equatorial productivity and subsurface oxygen demand to be less concentrated in the east. These findings shed light on how warm-climate states such as during the MCO may affect ocean oxygenation. If the MCO is considered as a possible analogue for future warming, our findings seem to support models suggesting that the recent deoxygenation trend and expansion of the ETP ODZ may eventually reverse3,4.

Basin-scale reconstruction of euxinia and Late Devonian mass extinctions.

The Devonian-Carboniferous transition marks a fundamental shift in the surface environment primarily related to changes in ocean-atmosphere oxidation states1,2, resulting from the continued proliferation of vascular land plants that stimulated the hydrological cycle and continental weathering3,4, glacioeustasy5,6, eutrophication and anoxic expansion in epicontinental seas3,4, and mass extinction events2,7,8. Here we present a comprehensive spatial and temporal compilation of geochemical data from 90 cores across the entire Bakken Shale (Williston Basin, North America). Our dataset allows for the detailed documentation of stepwise transgressions of toxic euxinic waters into the shallow oceans that drove a series of Late Devonian extinction events. Other Phanerozoic extinctions have also been related to the expansion of shallow-water euxinia, indicating that hydrogen sulfide toxicity was a key driver of Phanerozoic biodiversity.

Uncovering the Ediacaran phosphorus cycle.

Phosphorus is a limiting nutrient that is thought to control oceanic oxygen levels to a large extent1-3. A possible increase in marine phosphorus concentrations during the Ediacaran Period (about 635-539 million years ago) has been proposed as a driver for increasing oxygen levels4-6. However, little is known about the nature and evolution of phosphorus cycling during this time4. Here we use carbonate-associated phosphate (CAP) from six globally distributed sections to reconstruct oceanic phosphorus concentrations during a large negative carbon-isotope excursion-the Shuram excursion (SE)-which co-occurred with global oceanic oxygenation7-9. Our data suggest pulsed increases in oceanic phosphorus concentrations during the falling and rising limbs of the SE. Using a quantitative biogeochemical model, we propose that this observation could be explained by carbon dioxide and phosphorus release from marine organic-matter oxidation primarily by sulfate, with further phosphorus release from carbon-dioxide-driven weathering on land. Collectively, this may have resulted in elevated organic-pyrite burial and ocean oxygenation. Our CAP data also seem to suggest equivalent oceanic phosphorus concentrations under maximum and minimum extents of ocean anoxia across the SE. This observation may reflect decoupled phosphorus and ocean anoxia cycles, as opposed to their coupled nature in the modern ocean. Our findings point to external stimuli such as sulfate weathering rather than internal oceanic phosphorus-oxygen cycling alone as a possible control on oceanic oxygenation in the Ediacaran. In turn, this may help explain the prolonged rise of atmospheric oxygen levels.

Continental configuration controls ocean oxygenation during the Phanerozoic.

The early evolutionary and much of the extinction history of marine animals is thought to be driven by changes in dissolved oxygen concentrations ([O2]) in the ocean1-3. In turn, [O2] is widely assumed to be dominated by the geological history of atmospheric oxygen (pO2)4,5. Here, by contrast, we show by means of a series of Earth system model experiments how continental rearrangement during the Phanerozoic Eon drives profound variations in ocean oxygenation and induces a fundamental decoupling in time between upper-ocean and benthic [O2]. We further identify the presence of state transitions in the global ocean circulation, which lead to extensive deep-ocean anoxia developing in the early Phanerozoic even under modern pO2. Our finding that ocean oxygenation oscillates over stable thousand-year (kyr) periods also provides a causal mechanism that might explain elevated rates of metazoan radiation and extinction during the early Palaeozoic Era6. The absence, in our modelling, of any simple correlation between global climate and ocean ventilation, and the occurrence of profound variations in ocean oxygenation independent of atmospheric pO2, presents a challenge to the interpretation of marine redox proxies, but also points to a hitherto unrecognized role for continental configuration in the evolution of the biosphere.

Enhanced ocean oxygenation during Cenozoic warm periods.

Dissolved oxygen (O2) is essential for most ocean ecosystems, fuelling organisms' respiration and facilitating the cycling of carbon and nutrients. Oxygen measurements have been interpreted to indicate that the ocean's oxygen-deficient zones (ODZs) are expanding under global warming1,2. However, models provide an unclear picture of future ODZ change in both the near term and the long term3-6. The paleoclimate record can help explore the possible range of ODZ changes in warmer-than-modern periods. Here we use foraminifera-bound nitrogen (N) isotopes to show that water-column denitrification in the eastern tropical North Pacific was greatly reduced during the Middle Miocene Climatic Optimum (MMCO) and the Early Eocene Climatic Optimum (EECO). Because denitrification is restricted to oxygen-poor waters, our results indicate that, in these two Cenozoic periods of sustained warmth, ODZs were contracted, not expanded. ODZ contraction may have arisen from a decrease in upwelling-fuelled biological productivity in the tropical Pacific, which would have reduced oxygen demand in the subsurface. Alternatively, invigoration of deep-water ventilation by the Southern Ocean may have weakened the ocean's 'biological carbon pump', which would have increased deep-ocean oxygen. The mechanism at play would have determined whether the ODZ contractions occurred in step with the warming or took centuries or millennia to develop. Thus, although our results from the Cenozoic do not necessarily apply to the near-term future, they might imply that global warming may eventually cause ODZ contraction.

Biogeochemical extremes and compound events in the ocean.

The ocean is warming, losing oxygen and being acidified, primarily as a result of anthropogenic carbon emissions. With ocean warming, acidification and deoxygenation projected to increase for decades, extreme events, such as marine heatwaves, will intensify, occur more often, persist for longer periods of time and extend over larger regions. Nevertheless, our understanding of oceanic extreme events that are associated with warming, low oxygen concentrations or high acidity, as well as their impacts on marine ecosystems, remains limited. Compound events-that is, multiple extreme events that occur simultaneously or in close sequence-are of particular concern, as their individual effects may interact synergistically. Here we assess patterns and trends in open ocean extremes based on the existing literature as well as global and regional model simulations. Furthermore, we discuss the potential impacts of individual and compound extremes on marine organisms and ecosystems. We propose a pathway to improve the understanding of extreme events and the capacity of marine life to respond to them. The conditions exhibited by present extreme events may be a harbinger of what may become normal in the future. As a consequence, pursuing this research effort may also help us to better understand the responses of marine organisms and ecosystems to future climate change.

Widespread deoxygenation of temperate lakes.

The concentration of dissolved oxygen in aquatic systems helps to regulate biodiversity1,2, nutrient biogeochemistry3, greenhouse gas emissions4, and the quality of drinking water5. The long-term declines in dissolved oxygen concentrations in coastal and ocean waters have been linked to climate warming and human activity6,7, but little is known about the changes in dissolved oxygen concentrations in lakes. Although the solubility of dissolved oxygen decreases with increasing water temperatures, long-term lake trajectories are difficult to predict. Oxygen losses in warming lakes may be amplified by enhanced decomposition and stronger thermal stratification8,9 or oxygen may increase as a result of enhanced primary production10. Here we analyse a combined total of 45,148 dissolved oxygen and temperature profiles and calculate trends for 393 temperate lakes that span 1941 to 2017. We find that a decline in dissolved oxygen is widespread in surface and deep-water habitats. The decline in surface waters is primarily associated with reduced solubility under warmer water temperatures, although dissolved oxygen in surface waters increased in a subset of highly productive warming lakes, probably owing to increasing production of phytoplankton. By contrast, the decline in deep waters is associated with stronger thermal stratification and loss of water clarity, but not with changes in gas solubility. Our results suggest that climate change and declining water clarity have altered the physical and chemical environment of lakes. Declines in dissolved oxygen in freshwater are 2.75 to 9.3 times greater than observed in the world's oceans6,7 and could threaten essential lake ecosystem services2,3,5,11.

A 200-million-year delay in permanent atmospheric oxygenation.

The rise of atmospheric oxygen fundamentally changed the chemistry of surficial environments and the nature of Earth's habitability1. Early atmospheric oxygenation occurred over a protracted period of extreme climatic instability marked by multiple global glaciations2,3, with the initial rise of oxygen concentration to above 10-5 of the present atmospheric level constrained to about 2.43 billion years ago4,5. Subsequent fluctuations in atmospheric oxygen levels have, however, been reported to have occurred until about 2.32 billion years ago4, which represents the estimated timing of irreversible oxygenation of the atmosphere6,7. Here we report a high-resolution reconstruction of atmospheric and local oceanic redox conditions across the final two glaciations of the early Palaeoproterozoic era, as documented by marine sediments from the Transvaal Supergroup, South Africa. Using multiple sulfur isotope and iron-sulfur-carbon systematics, we demonstrate continued oscillations in atmospheric oxygen levels after about 2.32 billion years ago that are linked to major perturbations in ocean redox chemistry and climate. Oxygen levels thus fluctuated across the threshold of 10-5 of the present atmospheric level for about 200 million years, with permanent atmospheric oxygenation finally arriving with the Lomagundi carbon isotope excursion at about 2.22 billion years ago, some 100 million years later than currently estimated.

Dynamic in situ detection in iRhizo-Chip reveals diurnal fluctuations of Bacillus subtilis in the rhizosphere.

Effective colonization by microbe in the rhizosphere is critical for establishing a beneficial symbiotic relationship with the host plant. Bacillus subtilis, a soil-dwelling bacterium that is commonly found in association with plants and their rhizosphere, has garnered interest for its potential to enhance plant growth, suppress pathogens, and contribute to sustainable agricultural practices. However, research on the dynamic distribution of B. subtilis within the rhizosphere and its interaction mechanisms with plant roots remains insufficient due to limitations in existing in situ detection methodologies. To achieve dynamic in situ detection of the rhizosphere environment, we established iRhizo-Chip, a microfluidics-based platform. Using this device to investigate microbial behavior within the rhizosphere, we found obvious diurnal fluctuations in the growth of B. subtilis in the rhizosphere. Temporal dynamic analysis of rhizosphere dissolved oxygen (DO), pH, dissolved organic carbon, and reactive oxygen species showed that diurnal fluctuations in the growth of B. subtilis are potentially related to a variety of environmental factors. Spatial dynamic analysis also showed that the spatial distribution changes of B. subtilis and DO and pH were similar. Subsequently, through in vitro control experiments, we proved that rhizosphere DO and pH are the main driving forces for diurnal fluctuations in the growth of B. subtilis. Our results show that the growth of B. subtilis is driven by rhizosphere DO and pH, resulting in diurnal fluctuations, and iRhizo-Chip is a valuable tool for studying plant rhizosphere dynamics.

Nasal High-Flow Therapy during Neonatal Endotracheal Intubation.

BACKGROUND: Neonatal endotracheal intubation often involves more than one attempt, and oxygen desaturation is common. It is unclear whether nasal high-flow therapy, which extends the time to desaturation during elective intubation in children and adults receiving general anesthesia, can improve the likelihood of successful neonatal intubation on the first attempt. METHODS: We performed a randomized, controlled trial to compare nasal high-flow therapy with standard care (no nasal high-flow therapy or supplemental oxygen) in neonates undergoing oral endotracheal intubation at two Australian tertiary neonatal intensive care units. Randomization of intubations to the high-flow group or the standard-care group was stratified according to trial center, the use of premedication for intubation (yes or no), and postmenstrual age of the infant (≤28 or >28 weeks). The primary outcome was successful intubation on the first attempt without physiological instability (defined as an absolute decrease in the peripheral oxygen saturation of >20% from the preintubation baseline level or bradycardia with a heart rate of <100 beats per minute) in the infant. RESULTS: The primary intention-to-treat analysis included the outcomes of 251 intubations in 202 infants; 124 intubations were assigned to the high-flow group and 127 to the standard-care group. The infants had a median postmenstrual age of 27.9 weeks and a median weight of 920 g at the time of intubation. A successful intubation on the first attempt without physiological instability was achieved in 62 of 124 intubations (50.0%) in the high-flow group and in 40 of 127 intubations (31.5%) in the standard-care group (adjusted risk difference, 17.6 percentage points; 95% confidence interval [CI], 6.0 to 29.2), for a number needed to treat of 6 (95% CI, 4 to 17) for 1 infant to benefit. Successful intubation on the first attempt regardless of physiological stability was accomplished in 68.5% of the intubations in the high-flow group and in 54.3% of the intubations in the standard-care group (adjusted risk difference, 15.8 percentage points; 95% CI, 4.3 to 27.3). CONCLUSIONS: Among infants undergoing endotracheal intubation at two Australian tertiary neonatal intensive care units, nasal high-flow therapy during the procedure improved the likelihood of successful intubation on the first attempt without physiological instability in the infant. (Funded by the National Health and Medical Research Council; Australian New Zealand Clinical Trials Registry number, ACTRN12618001498280.).

Mineral protection regulates long-term global preservation of natural organic carbon.

The balance between photosynthetic organic carbon production and respiration controls atmospheric composition and climate1,2. The majority of organic carbon is respired back to carbon dioxide in the biosphere, but a small fraction escapes remineralization and is preserved over geological timescales3. By removing reduced carbon from Earth's surface, this sequestration process promotes atmospheric oxygen accumulation2 and carbon dioxide removal1. Two major mechanisms have been proposed to explain organic carbon preservation: selective preservation of biochemically unreactive compounds4,5 and protection resulting from interactions with a mineral matrix6,7. Although both mechanisms can operate across a range of environments and timescales, their global relative importance on 1,000-year to 100,000-year timescales remains uncertain4. Here we present a global dataset of the distributions of organic carbon activation energy and corresponding radiocarbon ages in soils, sediments and dissolved organic carbon. We find that activation energy distributions broaden over time in all mineral-containing samples. This result requires increasing bond-strength diversity, consistent with the formation of organo-mineral bonds8 but inconsistent with selective preservation. Radiocarbon ages further reveal that high-energy, mineral-bound organic carbon persists for millennia relative to low-energy, unbound organic carbon. Our results provide globally coherent evidence for the proposed7 importance of mineral protection in promoting organic carbon preservation. We suggest that similar studies of bond-strength diversity in ancient sediments may reveal how and why organic carbon preservation-and thus atmospheric composition and climate-has varied over geological time.

Oceanic nutrient rise and the late Miocene inception of Pacific oxygen-deficient zones.

The modern Pacific Ocean hosts the largest oxygen-deficient zones (ODZs), where oxygen concentrations are so low that nitrate is used to respire organic matter. The history of the ODZs may offer key insights into ocean deoxygenation under future global warming. In a 12-My record from the southeastern Pacific, we observe a >10‰ increase in foraminifera-bound nitrogen isotopes (15N/14N) since the late Miocene (8 to 9 Mya), indicating large ODZs expansion. Coinciding with this change, we find a major increase in the nutrient content of the ocean, reconstructed from phosphorus and iron measurements of hydrothermal sediments at the same site. Whereas global warming studies cast seawater oxygen concentrations as mainly dependent on climate and ocean circulation, our findings indicate that modern ODZs are underpinned by historically high concentrations of seawater phosphate.

Machine learning-based inverse design for electrochemically controlled microscopic gradients of O2 and H2O2.

A fundamental understanding of extracellular microenvironments of O2 and reactive oxygen species (ROS) such as H2O2, ubiquitous in microbiology, demands high-throughput methods of mimicking, controlling, and perturbing gradients of O2 and H2O2 at microscopic scale with high spatiotemporal precision. However, there is a paucity of high-throughput strategies of microenvironment design, and it remains challenging to achieve O2 and H2O2 heterogeneities with microbiologically desirable spatiotemporal resolutions. Here, we report the inverse design, based on machine learning (ML), of electrochemically generated microscopic O2 and H2O2 profiles relevant for microbiology. Microwire arrays with suitably designed electrochemical catalysts enable the independent control of O2 and H2O2 profiles with spatial resolution of ∼101 µm and temporal resolution of ∼10° s. Neural networks aided by data augmentation inversely design the experimental conditions needed for targeted O2 and H2O2 microenvironments while being two orders of magnitude faster than experimental explorations. Interfacing ML-based inverse design with electrochemically controlled concentration heterogeneity creates a viable fast-response platform toward better understanding the extracellular space with desirable spatiotemporal control.

A Neandertal dietary conundrum: Insights provided by tooth enamel Zn isotopes from Gabasa, Spain.

The characterization of Neandertals' diets has mostly relied on nitrogen isotope analyses of bone and tooth collagen. However, few nitrogen isotope data have been recovered from bones or teeth from Iberia due to poor collagen preservation at Paleolithic sites in the region. Zinc isotopes have been shown to be a reliable method for reconstructing trophic levels in the absence of organic matter preservation. Here, we present the results of zinc (Zn), strontium (Sr), carbon (C), and oxygen (O) isotope and trace element ratio analysis measured in dental enamel on a Pleistocene food web in Gabasa, Spain, to characterize the diet and ecology of a Middle Paleolithic Neandertal individual. Based on the extremely low δ66Zn value observed in the Neandertal's tooth enamel, our results support the interpretation of Neandertals as carnivores as already suggested by δ15N isotope values of specimens from other regions. Further work could help identify if such isotopic peculiarities (lowest δ66Zn and highest δ15N of the food web) are due to a metabolic and/or dietary specificity of the Neandertals.

Microbial diversity and abundance vary along salinity, oxygen, and particle size gradients in the Chesapeake Bay.

Marine snow and other particles are abundant in estuaries, where they drive biogeochemical transformations and elemental transport. Particles range in size, thereby providing a corresponding gradient of habitats for marine microorganisms. We used standard normalized amplicon sequencing, verified with microscopy, to characterize taxon-specific microbial abundances, (cells per litre of water and per milligrams of particles), across six particle size classes, ranging from 0.2 to 500 µm, along the main stem of the Chesapeake Bay estuary. Microbial communities varied in salinity, oxygen concentrations, and particle size. Many taxonomic groups were most densely packed on large particles (in cells/mg particles), yet were primarily associated with the smallest particle size class, because small particles made up a substantially larger portion of total particle mass. However, organisms potentially involved in methanotrophy, nitrite oxidation, and sulphate reduction were found primarily on intermediately sized (5-180 µm) particles, where species richness was also highest. All abundant ostensibly free-living organisms, including SAR11 and Synecococcus, appeared on particles, albeit at lower abundance than in the free-living fraction, suggesting that aggregation processes may incorporate them into particles. Our approach opens the door to a more quantitative understanding of the microscale and macroscale biogeography of marine microorganisms.