Remineralization of organic matter by heterotrophic organisms regulates the biological sequestration of carbon, thereby mediating atmospheric CO2. While surface nutrient supply impacts the elemental ratios of primary production, stoichiometric control by remineralization remains unclear. Here we develop a mechanistic description of remineralization and its stoichiometry in a marine microbial ecosystem model. The model simulates the observed elemental plasticity of phytoplankton and the relatively constant, lower C:N of heterotrophic biomass. In addition, the model captures the observed decreases in the C:N of more labile dissolved organic matter (DOM) and the C:N of its remineralization with depth, which are driven by a switch in the dominant source of DOM from phytoplankton to heterotrophic biomass. Only a model version with targeted remineralization of N-rich components is able to simulate the observed profiles of preferential remineralization of N relative to C and the elevated C:N of bulk DOM. The model suggests that more labile substrates are associated with C-limited heterotrophic growth and not with preferential remineralization, while less labile substrates are associated with growth limited by processing rates and with preferential remineralization. The resulting patterns of variable remineralization stoichiometry mediate the extent to which a proportional increase in carbon production resulting from changes in phytoplankton stoichiometry can increase the efficiency of the biological pump. Results emphasize the importance of understanding the physiology of both phytoplankton and heterotrophs for anticipating changes in biologically driven ocean carbon storage.

Ammonia oxidation to nitrite and its subsequent oxidation to nitrate provides energy to the two populations of nitrifying chemo-autotrophs in the energy-starved dark ocean, driving a coupling between reduced inorganic nitrogen (N) pools and production of new organic carbon (C) in the dark ocean. However, the relationship between the flux of new C production and the fluxes of N of the two steps of oxidation remains unclear. Here, we show that, despite orders-of-magnitude difference in cell abundances between ammonia oxidizers and nitrite oxidizers, the two populations sustain similar bulk N-oxidation rates throughout the deep waters with similarly high affinities for ammonia and nitrite under increasing substrate limitation, thus maintaining overall homeostasis in the oceanic nitrification pathway. Our observations confirm the theoretical predictions of a redox-informed ecosystem model. Using balances from this model, we suggest that consistently low ammonia and nitrite concentrations are maintained when the two populations have similarly high substrate affinities and their loss rates are proportional to their maximum growth rates. The stoichiometric relations between the fluxes of C and N indicate a threefold to fourfold higher C-fixation efficiency per mole of N oxidized by ammonia oxidizers compared to nitrite oxidizers due to nearly identical apparent energetic requirements for C fixation of the two populations. We estimate that the rate of chemoautotrophic C fixation amounts to similar to 1 x 10(13) to similar to 2 x 10(13) mol of C per year globally through the flux of similar to 1 x 10(14) to similar to 2 x 10(14) mol of N per year of the two steps of oxidation throughout the dark ocean.

Marine microbes form the base of ocean food webs and drive ocean biogeochemical cycling. Yet little is known about the ability of microbial populations to adapt as they are advected through changing conditions. Here, we investigated the interplay between physical and biological timescales using a model of adaptation and an eddy-resolving ocean circulation climate model. Two criteria were identified that relate the timing and nature of adaptation to the ratio of physical to biological timescales. Genetic adaptation was impeded in highly variable regimes by nongenetic modifications but was promoted in more stable environments. An evolutionary trade-off emerged where greater short-term nongenetic transgenerational effects (low-gamma strategy) enabled rapid responses to environmental fluctuations but delayed genetic adaptation, while fewer short-term transgenerational effects (high-gamma strategy) allowed faster genetic adaptation but inhibited short-term responses. Our results demonstrate that the selective pressures for organisms within a single water mass vary based on differences in generation timescales resulting in different evolutionary strategies being favored. Organisms that experience more variable environments should favor a low-gamma strategy. Furthermore, faster cell division rates should be a key factor in genetic adaptation in a changing ocean. Understanding and quantifying the relationship between evolutionary and physical timescales is critical for robust predictions of future microbial dynamics.

Methylmercury is greatly bioconcentrated and biomagnified in marine plankton ecosystems, and these communities form the basis of marine food webs. Therefore, the evaluation of the potential exposure of methylmercury to higher trophic levels, including humans, requires a better understanding of its distribution in the ocean and the factors that control its biomagnification. In this study, a coupled physical/ecological model is used to simulate the trophic transfer of monomethylmercury (MMHg) in a marine plankton ecosystem. The model includes phytoplankton, a microbial community, herbivorous zooplankton (HZ), and carnivorous zooplankton (CZ). The model captures both shorter food chains in oligotrophic regions, with small HZ feeding on small phytoplankton, and longer chains in higher nutrient conditions, with larger HZ feeding on larger phytoplankton and larger CZ feeding on larger HZ. In the model, trophic dilution occurs in the food webs that involve small zooplankton, as the grazing fluxes of small zooplankton are insufficient to accumulate more MMHg in themselves than in their prey. The model suggests that biomagnification is more prominent in large zooplankton and that the microbial community plays an important role in the trophic transfer of MMHg. Sensitivity analyses show that with increasing body size, the sensitivity of the trophic magnification ratio to grazing, mortality rates, and food assimilation efficiency (AE(C)) increases, while the sensitivity to excretion rates decreases. More predation or a longer zooplankton lifespan may lead to more prominent biomagnification, especially for large species. Because lower AE(C) results in more predation, modeled ratios of MMHg concentrations between large plankton are doubled or even tripled when the AE(C) decreases from 50% to 10%. This suggests that the biomagnification of large zooplankton is particularly sensitive to food assimilation efficiency.

Microbial activity mediates the fluxes of greenhouse gases. However, in the global models of the marine and terrestrial biospheres used for climate change projections, typically only photosynthetic microbial activity is resolved mechanistically. To move forward, we argue that global biogeochemical models need a theoretically grounded framework with which to constrain parameterizations of diverse microbial metabolisms. Here, we explain how the key redox chemistry underlying metabolisms provides a path towards this goal. Using this first-principles approach, the presence or absence of metabolic functional types emerges dynamically from ecological interactions, expanding model applicability to unobserved environments."Nothing is less real than realism. It is only by selection, by elimination, by emphasis, that we get at the real meaning of things." -Georgia O'Keefe Marine microbial activities fuel biogeochemical cycles that impact the climate, but global models do not account for the myriad physiological processes that microbes perform. Here the authors argue for a model framework that reinterprets the ocean as physics coupled to biologically-driven redox chemistry.

Anaerobic microbial activity in the ocean causes losses of bioavailable nitrogen and emission of nitrous oxide to the atmosphere, but its predictability at global scales remains limited. Resource ratio theory suggests that anaerobic activity becomes sustainable when the ratio of oxygen to organic matter supply is below the ratio required by aerobic metabolisms. Here, we demonstrate the relevance of this framework at the global scale using three-dimensional ocean datasets, providing a new interpretation of existing observations. Evaluations of the location and extent of anoxic zones and a diagnostic rate of pelagic nitrogen loss are consistent with previous estimates. However, we demonstrate that a threshold based on substrate-supply fluxes is qualitatively different from a threshold based solely on the ambient oxygen concentration. This implies that use of the flux-based threshold in global biogeochemical models can result in different predictions of anaerobic activity and nitrogen loss.

Organic matter constitutes a key reservoir in global elemental cycles. However, our understanding of the dynamics of organic matter and its accumulation remains incomplete. Seemingly disparate hypotheses have been proposed to explain organic matter accumulation: the slow degradation of intrinsically recalcitrant substrates, the depletion to concentrations that inhibit microbial consumption, and a dependency on the consumption capabilities of nearby microbial populations. Here, using a mechanistic model, we develop a theoretical framework that explains how organic matter predictably accumulates in natural environments due to biochemical, ecological, and environmental factors. Our framework subsumes the previous hypotheses. Changes in the microbial community or the environment can move a class of organic matter from a state of functional recalcitrance to a state of depletion by microbial consumers. The model explains the vertical profile of dissolved organic carbon in the ocean and connects microbial activity at subannual timescales to organic matter turnover at millennial timescales. The threshold behavior of the model implies that organic matter accumulation may respond nonlinearly to changes in temperature and other factors, providing hypotheses for the observed correlations between organic carbon reservoirs and temperature in past earth climates.

Prochlorococcus is both the smallest and numerically most abundant photosynthesizing organism on the planet. While thriving in the warm oligotrophic gyres, Prochlorococcus concentrations drop rapidly in higher-latitude regions. Transect data from the North Pacific show the collapse occurring at a wide range of temperatures and latitudes (temperature is often hypothesized to cause this shift), suggesting an ecological mechanism may be at play. An often used size-based theory of phytoplankton community structure that has been incorporated into computational models correctly predicts the dominance of Prochlorococcus in the gyres, and the relative dominance of larger cells at high latitudes. However, both theory and computational models fail to explain the poleward collapse. When heterotrophic bacteria and predators that prey nonspecifically on both Prochlorococcus and bacteria are included in the theoretical framework, the collapse of Prochlorococcus occurs with increasing nutrient supplies. The poleward collapse of Prochlorococcus populations then naturally emerges when this mechanism of "shared predation" is implemented in a complex global ecosystem model. Additionally, the theory correctly predicts trends in both the abundance and mean size of the heterotrophic bacteria. These results suggest that ecological controls need to be considered to understand the biogeography of Prochlorococcus and predict its changes under future ocean conditions. Indirect interactions within a microbial network can be essential in setting community structure.

Micro-scale microbial community dynamics can substantially alter the fate of sinking particulates in the ocean thus playing a key role in setting the vertical flux of particulate carbon in the ocean.

Anthropogenic habitat loss and climate change are reducing species' geographic ranges, increasing extinction risk and losses of species' genetic diversity. Although preserving genetic diversity is key to maintaining species' adaptability, we lack predictive tools and global estimates of genetic diversity loss across ecosystems. We introduce a mathematical framework that bridges biodiversity theory and population genetics to understand the loss of naturally occurring DNA mutations with decreasing habitat. By analyzing genomic variation of 10,095 georeferenced individuals from 20 plant and animal species, we show that genome-wide diversity follows a mutations-area relationship power law with geographic area, which can predict genetic diversity loss from local population extinctions. We estimate that more than 10% of genetic diversity may already be lost for many threatened and nonthreatened species, surpassing the United Nations' post-2020 targets for genetic preservation.