Emergent trade-offs among plasticity strategies in mixotrophs.

Marine mixotrophs combine phagotrophy and phototrophy to acquire the resources they need for growth. Metabolic plasticity, the ability for individuals to dynamically alter their relative investment between different metabolic processes, allows mixotrophs to efficiently exploit variable environmental conditions. Different mixotrophs may vary in how quickly they respond to environmental stimuli, with slow-responding mixotrophs exhibiting a significant lag between a change in the environment and the resulting change metabolic strategy. In this study, we develop a model of mixotroph metabolic strategy and explore how the rate of the plastic response affects the seasonality, competitive fitness, and biogeochemical role of mixotroph populations. Fast-responding mixotrophs are characterized by more efficient resource use and higher average growth rates than slow-responding mixotrophs because any lag in the plastic response following a change in environmental conditions creates a mismatch between the mixotroph's metabolic requirements and their resource acquisition. However, this mismatch also results in increased storage of unused resources that support growth under future nutrient-limited conditions. As a result of this trade-off, mixotroph biomass and productivity are maximized at intermediate plastic response rates. Furthermore, the trade-off represents a mechanism for coexistence between fast-responding and slow-responding mixotrophs. In mixed communities, fast-responding mixotrophs are numerically dominant, but slow-responding mixotrophs persist at low abundance due to the provisioning effect that emerges as a result of their less efficient resource acquisition strategy. In addition to increased competitive ability, fast-responding mixotrophs are, on average, more autotrophic than slow-responding mixotrophs. Notably, these trade-offs associated with mixotroph response rate arise without including an explicit physiological cost associated with plasticity, a conclusion that may provide insight into evolutionary constraints of metabolic plasticity in mixotrophic organisms. When an explicit cost is added to the model, it alters the competitive relationships between fast- and slow-responding mixotrophs. Faster plastic response rates are favored by lower physiological costs as well as higher amplitude seasonal cycles.

Phytoplankton thermal trait parameterization alters community structure and biogeochemical processes in a modeled ocean.

Phytoplankton exhibit diverse physiological responses to temperature which influence their fitness in the environment and consequently alter their community structure. Here, we explored the sensitivity of phytoplankton community structure to thermal response parameterization in a modelled marine phytoplankton community. Using published empirical data, we evaluated the maximum thermal growth rates (µmax ) and temperature coefficients (Q10 ; the rate at which growth scales with temperature) of six key Phytoplankton Functional Types (PFTs): coccolithophores, cyanobacteria, diatoms, diazotrophs, dinoflagellates, and green algae. Following three well-documented methods, PFTs were either assumed to have (1) the same µmax and the same Q10 (as in to Eppley, 1972), (2) a unique µmax but the same Q10 (similar to Kremer et al., 2017), or (3) a unique µmax and a unique Q10 (following Anderson et al., 2021). These trait values were then implemented within the Massachusetts Institute of Technology biogeochemistry and ecosystem model (called Darwin) for each PFT under a control and climate change scenario. Our results suggest that applying a µmax and Q10 universally across PFTs (as in Eppley, 1972) leads to unrealistic phytoplankton communities, which lack diatoms globally. Additionally, we find that accounting for differences in the Q10 between PFTs can significantly impact each PFT's competitive ability, especially at high latitudes, leading to altered modeled phytoplankton community structures in our control and climate change simulations. This then impacts estimates of biogeochemical processes, with, for example, estimates of export production varying by ~10% in the Southern Ocean depending on the parameterization. Our results indicate that the diversity of thermal response traits in phytoplankton not only shape community composition in the historical and future, warmer ocean, but that these traits have significant feedbacks on global biogeochemical cycles.

A Lagrangian model for drifting ecosystems reveals heterogeneity-driven enhancement of marine plankton blooms.

Marine plankton play a crucial role in carbon storage, global climate, and ecosystem function. Planktonic ecosystems are embedded in patches of water that are continuously moving, stretching, and diluting. These processes drive inhomegeneities on a range of scales, with implications for the integrated ecosystem properties, but are hard to characterize. We present a theoretical framework that accounts for all these aspects; tracking the water patch hosting a drifting ecosystem along with its physical, environmental, and biochemical features. The theory resolves patch dilution and internal physical mixing as a function of oceanic strain and diffusion. Ecological dynamics are parameterized by an idealized nutrient and phytoplankton population and we specifically capture the time evolution of the biochemical spatial variances to represent within-patch heterogeneity. We find that, depending only on the physical processes to which the water patch is subjected, the plankton biomass response to a resource perturbation can vary in size up to six times. This work indicates that we must account for these processes when interpreting and modeling marine ecosystems and provides a framework with which to do so.

Mixoplankton and mixotrophy: future research priorities.

Phago-mixotrophy, the combination of photoautotrophy and phagotrophy in mixoplankton, organisms that can combine both trophic strategies, have gained increasing attention over the past decade. It is now recognized that a substantial number of protistan plankton species engage in phago-mixotrophy to obtain nutrients for growth and reproduction under a range of environmental conditions. Unfortunately, our current understanding of mixoplankton in aquatic systems significantly lags behind our understanding of zooplankton and phytoplankton, limiting our ability to fully comprehend the role of mixoplankton (and phago-mixotrophy) in the plankton food web and biogeochemical cycling. Here, we put forward five research directions that we believe will lead to major advancement in the field: (i) evolution: understanding mixotrophy in the context of the evolutionary transition from phagotrophy to photoautotrophy; (ii) traits and trade-offs: identifying the key traits and trade-offs constraining mixotrophic metabolisms; (iii) biogeography: large-scale patterns of mixoplankton distribution; (iv) biogeochemistry and trophic transfer: understanding mixoplankton as conduits of nutrients and energy; and (v) in situ methods: improving the identification of in situ mixoplankton and their phago-mixotrophic activity.

Global climate-change trends detected in indicators of ocean ecology.

Strong natural variability has been thought to mask possible climate-change-driven trends in phytoplankton populations from Earth-observing satellites. More than 30 years of continuous data were thought to be needed to detect a trend driven by climate change1. Here we show that climate-change trends emerge more rapidly in ocean colour (remote-sensing reflectance, Rrs), because Rrs is multivariate and some wavebands have low interannual variability. We analyse a 20-year Rrs time series from the Moderate Resolution Imaging Spectroradiometer (MODIS) aboard the Aqua satellite, and find significant trends in Rrs for 56% of the global surface ocean, mainly equatorward of 40°. The climate-change signal in Rrs emerges after 20 years in similar regions covering a similar fraction of the ocean in a state-of-the-art ecosystem model2, which suggests that our observed trends indicate shifts in ocean colour-and, by extension, in surface-ocean ecosystems-that are driven by climate change. On the whole, low-latitude oceans have become greener in the past 20 years.

High Growth Rate of Diatoms Explained by Reduced Carbon Requirement and Low Energy Cost of Silica Deposition.

The rapid growth of diatoms makes them one of the most pervasive and productive types of plankton in the world's ocean, but the physiological basis for their high growth rates remains poorly understood. Here, we evaluate the factors that elevate diatom growth rates, relative to other plankton, using a steady-state metabolic flux model that computes the photosynthetic C source from intracellular light attenuation and the carbon cost of growth from empirical cell C quotas, across a wide range of cell sizes. For both diatoms and other phytoplankton, growth rates decline with increased cell volume, consistent with observations, because the C cost of division increases with size faster than photosynthesis. However, the model predicts overall higher growth rates for diatoms due to reduced C requirements and the low energetic cost of Si deposition. The C savings from the silica frustule are supported by metatranscriptomic data from Tara Oceans, which show that the abundance of transcripts for cytoskeleton components in diatoms is lower than in other phytoplankton. Our results highlight the importance of understanding the origins of phylogenetic differences in cellular C quotas, and suggest that the evolution of silica frustules may play a critical role in the global dominance of marine diatoms. IMPORTANCE This study addresses a longstanding issue regarding diatoms, namely, their fast growth. Diatoms, which broadly are phytoplankton with silica frustules, are the world's most productive microorganisms and dominate in polar and upwelling regions. Their dominance is largely supported by their high growth rate, but the physiological reasoning behind that characteristic has been obscure. In this study, we combine a quantitative model and metatranscriptomic approaches and show that diatoms' low carbon requirements and low energy costs for silica frustule production are the key factors supporting their fast growth. Our study suggests that the effective use of energy-efficient silica as a cellular structure, instead of carbon, enables diatoms to be the most productive organisms in the global ocean.

Global patterns in marine organic matter stoichiometry driven by phytoplankton ecophysiology.

The proportion of major elements in marine organic matter links cellular processes to global nutrient, oxygen and carbon cycles. Differences in the C:N:P ratios of organic matter have been observed between ocean biomes, but these patterns have yet to be quantified from the underlying small-scale physiological and ecological processes. Here we use an ecosystem model that includes adaptive resource allocation within and between ecologically distinct plankton size classes to attribute the causes of global patterns in the C:N:P ratios. We find that patterns of N:C variation are largely driven by common physiological adjustment strategies across all phytoplankton, while patterns of N:P are driven by ecological selection for taxonomic groups with different phosphorus storage capacities. Although N:C varies widely due to cellular adjustment to light and nutrients, its latitudinal gradient is modest because of depth-dependent trade-offs between nutrient and light availability. Strong latitudinal variation in N:P reflects an ecological balance favouring small plankton with lower P storage capacity in the subtropics, and larger eukaryotes with a higher cellular P storage capacity in nutrient-rich high latitudes. A weaker N:P difference between southern and northern hemispheres, and between the Atlantic and Pacific oceans, reflects differences in phosphate available for cellular storage. Despite simulating only two phytoplankton size classes, the emergent global variability of elemental ratios resembles that of all measured species, suggesting that the range of growth conditions and ecological selection sustain the observed diversity of stoichiometry among phytoplankton.

A nutrient relay sustains subtropical ocean productivity.

The expansive gyres of the subtropical ocean account for a significant fraction of global organic carbon export from the upper ocean. In the gyre interior, vertical mixing and the heaving of nutrient-rich waters into the euphotic layer sustain local productivity, in turn depleting the layers below. However, the nutrient pathways by which these subeuphotic layers are themselves replenished remain unclear. Using a global, eddy-permitting simulation of ocean physics and biogeochemistry, we quantify nutrient resupply mechanisms along and across density surfaces, including the contribution of eddy-scale motions that are challenging to observe. We find that mesoscale eddies (10 to 100 km) flux nutrients from the shallow flanks of the gyre into the recirculating interior, through time-varying motions along density surfaces. The subeuphotic layers are ultimately replenished in approximately equal contributions by this mesoscale eddy transport and the remineralization of sinking particles. The mesoscale eddy resupply is most important in the lower thermocline for the whole subtropical region but is dominant at all depths within the gyre interior. Subtropical gyre productivity may therefore be sustained by a nutrient relay, where the lateral transport resupplies nutrients to the thermocline and allows vertical exchanges to maintain surface biological production and carbon export.

Attribution of Space-Time Variability in Global-Ocean Dissolved Inorganic Carbon.

The inventory and variability of oceanic dissolved inorganic carbon (DIC) is driven by the interplay of physical, chemical, and biological processes. Quantifying the spatiotemporal variability of these drivers is crucial for a mechanistic understanding of the ocean carbon sink and its future trajectory. Here, we use the Estimating the Circulation and Climate of the Ocean-Darwin ocean biogeochemistry state estimate to generate a global-ocean, data-constrained DIC budget and investigate how spatial and seasonal-to-interannual variability in three-dimensional circulation, air-sea CO2 flux, and biological processes have modulated the ocean sink for 1995-2018. Our results demonstrate substantial compensation between budget terms, resulting in distinct upper-ocean carbon regimes. For example, boundary current regions have strong contributions from vertical diffusion while equatorial regions exhibit compensation between upwelling and biological processes. When integrated across the full ocean depth, the 24-year DIC mass increase of 64 Pg C (2.7 Pg C year-1) primarily tracks the anthropogenic CO2 growth rate, with biological processes providing a small contribution of 2% (1.4 Pg C). In the upper 100 m, which stores roughly 13% (8.1 Pg C) of the global increase, we find that circulation provides the largest DIC gain (6.3 Pg C year-1) and biological processes are the largest loss (8.6 Pg C year-1). Interannual variability is dominated by vertical advection in equatorial regions, with the 1997-1998 El Niño-Southern Oscillation causing the largest year-to-year change in upper-ocean DIC (2.1 Pg C). Our results provide a novel, data-constrained framework for an improved mechanistic understanding of natural and anthropogenic perturbations to the ocean sink.

Ocean mover's distance: using optimal transport for analysing oceanographic data.

Remote sensing observations from satellites and global biogeochemical models have combined to revolutionize the study of ocean biogeochemical cycling, but comparing the two data streams to each other and across time remains challenging due to the strong spatial-temporal structuring of the ocean. Here, we show that the Wasserstein distance provides a powerful metric for harnessing these structured datasets for better marine ecosystem and climate predictions. The Wasserstein distance complements commonly used point-wise difference methods such as the root-mean-squared error, by quantifying differences in terms of spatial displacement in addition to magnitude. As a test case, we consider chlorophyll (a key indicator of phytoplankton biomass) in the northeast Pacific Ocean, obtained from model simulations, in situ measurements, and satellite observations. We focus on two main applications: (i) comparing model predictions with satellite observations, and (ii) temporal evolution of chlorophyll both seasonally and over longer time frames. The Wasserstein distance successfully isolates temporal and depth variability and quantifies shifts in biogeochemical province boundaries. It also exposes relevant temporal trends in satellite chlorophyll consistent with climate change predictions. Our study shows that optimal transport vectors underlying the Wasserstein distance provide a novel visualization tool for testing models and better understanding temporal dynamics in the ocean.

Trophic interactions with heterotrophic bacteria limit the range of Prochlorococcus.

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.

Bottom-Heavy Trophic Pyramids Impair Methylmercury Biomagnification in the Marine Plankton Ecosystems.

Methylmercury (CH3Hg+, MMHg) in the phytoplankton and zooplankton, which form the bottom of marine food webs, is a good predictor of MMHg in top predators, including humans. Therefore, evaluating the potential exposure of MMHg to higher trophic levels (TLs) requires a better understanding of relationships between MMHg biomagnification and plankton dynamics. In this study, a coupled ecological/physical model with 366 plankton types of different sizes, biogeochemical functions, and temperature tolerance is used to simulate the relationships between MMHg biomagnification and the ecosystem structure. The study shows that the MMHg biomagnification becomes more significant with increasing TLs. Trophic magnification factors (TMFs) in the lowest two TLs show the opposite spatial pattern to TMFs in higher TLs. The low TMFs are usually associated with a short food-chain length. The less bottom-heavy trophic pyramids in the oligotrophic oceans enhance the MMHg trophic transfer. The global average TMF is increased from 2.3 to 2.8 in the warmer future with a medium climate sensitivity of 2.5 °C. Our study suggests that if there are no mitigation measures for Hg emission, MMHg in the high-trophic-level plankton is increased more dramatically in the warming future, indicating greater MMHg exposure for top predators such as humans.

Abrupt shifts in 21st-century plankton communities.

Marine microbial communities sustain ocean food webs and mediate global elemental cycles. These communities will change with climate; these changes can be gradual or foreseeable but likely have much more substantial consequences when sudden and unpredictable. In a complex virtual marine microbial ecosystem, we find that climate change­driven shifts over the 21st century are often abrupt, large in amplitude and extent, and unpredictable using standard early warning signals. Phytoplankton with unique resource needs, especially fast-growing species such as diatoms, are more prone to abrupt shifts. Abrupt shifts in biomass, productivity, and community structure are concentrated in Atlantic and Pacific subtropics. Abrupt changes in environmental variables such as temperature and nutrients rarely precede these ecosystem shifts, indicating that rapid community restructuring can occur in response to gradual environmental changes, particularly in nutrient supply rate ratios.

Future phytoplankton diversity in a changing climate.

The future response of marine ecosystem diversity to continued anthropogenic forcing is poorly constrained. Phytoplankton are a diverse set of organisms that form the base of the marine ecosystem. Currently, ocean biogeochemistry and ecosystem models used for climate change projections typically include only 2-3 phytoplankton types and are, therefore, too simple to adequately assess the potential for changes in plankton community structure. Here, we analyse a complex ecosystem model with 35 phytoplankton types to evaluate the changes in phytoplankton community composition, turnover and size structure over the 21st century. We find that the rate of turnover in the phytoplankton community becomes faster during this century, that is, the community structure becomes increasingly unstable in response to climate change. Combined with alterations to phytoplankton diversity, our results imply a loss of ecological resilience with likely knock-on effects on the productivity and functioning of the marine environment.

Moving ecological and biogeochemical transitions across the North Pacific.

In the North Pacific Ocean, nutrient rich surface waters flow south from the subpolar gyre through a transitional region and into the subtropics. Along the way, nutrients are used, recycled, and exported, leading to lower biomass and a commensurate change in ecosystem structure moving southward. We focus on the region between the two gyres (the Transition Zone) using a coupled biophysical ocean model, remote sensing, floats, and cruise data to explore the nature of the physical, biogeochemical, and ecological fields in this region. Nonlinear interactions between biological processes and the meridional gradient in nutrient supply lead to sharp shifts across this zone. These transitions between a southern region with more uniform biological and biogeochemical properties and steep meridional gradients to the north are diagnosed from extrema in the first derivative of the properties with latitude. Some transitions like that for chlorophyll a (the transition zone chlorophyll front [TZCF]) experience large seasonal excursions while the location of the transitions in other properties moves very little. The seasonal shifts are not caused by changes in the horizontal flow field, but rather by the interaction of seasonal, depth related, forcing with the mean latitudinal gradients. Focusing on the TZCF as a case study, we express its phase velocity in terms of vertical nutrient flux and internal ecosystem processes, demonstrating their nearly equal influence on its motion. This framework of propagating biogeochemical transitions can be systematically expanded to better understand the processes that structure ecosystems and biogeochemistry in the North Pacific and beyond.

Global health effects of future atmospheric mercury emissions.

Mercury is a potent neurotoxin that poses health risks to the global population. Anthropogenic mercury emissions to the atmosphere are projected to decrease in the future due to enhanced policy efforts such as the Minamata Convention, a legally-binding international treaty entered into force in 2017. Here, we report the development of a comprehensive climate-atmosphere-land-ocean-ecosystem and exposure-risk model framework for mercury and its application to project the health effects of future atmospheric emissions. Our results show that the accumulated health effects associated with mercury exposure during 2010-2050 are $19 (95% confidence interval: 4.7-54) trillion (2020 USD) realized to 2050 (3% discount rate) for the current policy scenario. Our results suggest a substantial increase in global human health cost if emission reduction actions are delayed. This comprehensive modeling approach provides a much-needed tool to help parties to evaluate the effectiveness of Hg emission controls as required by the Minamata Convention.

Exploring biogeochemical and ecological redundancy in phytoplankton communities in the global ocean.

Climate-change-induced alterations of oceanic conditions will lead to the ecological niches of some marine phytoplankton species disappearing, at least regionally. How will such losses affect the ecosystem and the coupled biogeochemical cycles? Here, we couch this question in terms of ecological redundancy (will other species be able to fill the ecological roles of the extinct species) and biogeochemical redundancy (can other species replace their biogeochemical roles). Prior laboratory and field studies point to a spectrum in the degree of redundancy. We use a global three-dimensional computer model with diverse planktonic communities to explore these questions further. The model includes 35 phytoplankton types that differ in size, biogeochemical function and trophic strategy. We run two series of experiments in which single phytoplankton types are either partially or fully eliminated. The niches of the targeted types were not completely reoccupied, often with a reduction in the transfer of matter from autotrophs to heterotrophs. Primary production was often decreased, but sometimes increased due to reduction in grazing pressure. Complex trophic interactions (such as a decrease in the stocks of a predator's grazer) led to unexpected reshuffling of the community structure. Alterations in resource utilization may cause impacts beyond the regions where the type went extinct. Our results suggest a lack of redundancy, especially in the 'knock on' effects on higher trophic levels. Redundancy appeared lowest for types on the edges of trait space (e.g. smallest) or with unique competitive strategies. Though highly idealized, our modelling findings suggest that the results from laboratory or field studies often do not adequately capture the ramifications of functional redundancy. The modelled, often counterintuitive, responses-via complex food web interactions and bottom-up versus top-down controls-indicate that changes in planktonic community will be key determinants of future ocean global change ecology and biogeochemistry.

Elucidating ecological complexity: Unsupervised learning determines global marine eco-provinces.

An unsupervised learning method is presented for determining global marine ecological provinces (eco-provinces) from plankton community structure and nutrient flux data. The systematic aggregated eco-province (SAGE) method identifies eco-provinces within a highly nonlinear ecosystem model. To accommodate the non-Gaussian covariance of the data, SAGE uses t-stochastic neighbor embedding (t-SNE) to reduce dimensionality. Over a hundred eco-provinces are identified with the density-based spatial clustering of applications with noise (DBSCAN) algorithm. Using a connectivity graph with ecological dissimilarity as the distance metric, robust aggregated eco-provinces (AEPs) are objectively defined by nesting the eco-provinces. Using the AEPs, the control of nutrient supply rates on community structure is explored. Eco-provinces and AEPs are unique and aid model interpretation. They could facilitate model intercomparison and potentially improve understanding and monitoring of marine ecosystems.

Biomagnification of Methylmercury in a Marine Plankton Ecosystem.

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 (AEC) 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 AEC results in more predation, modeled ratios of MMHg concentrations between large plankton are doubled or even tripled when the AEC decreases from 50% to 10%. This suggests that the biomagnification of large zooplankton is particularly sensitive to food assimilation efficiency.

Microbial feedbacks optimize ocean iron availability.

Iron is the limiting factor for biological production over a large fraction of the surface ocean because free iron is rapidly scavenged or precipitated under aerobic conditions. Standing stocks of dissolved iron are maintained by association with organic molecules (ligands) produced by biological processes. We hypothesize a positive feedback between iron cycling, microbial activity, and ligand abundance: External iron input fuels microbial production, creating organic ligands that support more iron in seawater, leading to further macronutrient consumption until other microbial requirements such as macronutrients or light become limiting, and additional iron no longer increases productivity. This feedback emerges in numerical simulations of the coupled marine cycles of macronutrients and iron that resolve the dynamic microbial production and loss of iron-chelating ligands. The model solutions resemble modern nutrient distributions only over a finite range of prescribed ligand source/sink ratios where the model ocean is driven to global-scale colimitation by micronutrients and macronutrients and global production is maximized. We hypothesize that a global-scale selection for microbial ligand cycling may have occurred to maintain "just enough" iron in the ocean.