Phaeocystis antarctica is an important phytoplankton species in the Southern Ocean. We incorporated P. antarctica into the biogeochemical elemental cycling ocean model to study Southern Ocean ecosystem dynamics and biogeochemistry. The optimum values of ecological parameters for Phaeocystis were sought through synthesizing laboratory and field observations, and the model output was evaluated with observed chlorophyll a, carbon biomass, and nutrient distributions. Several factors have been proposed to control Southern Ocean ecosystem structure, including light adaptation, iron uptake capability, and loss processes. Optimum simulation results were obtained when P. antarctica had a relatively high α (P-I curve initial slope) value and a higher half-saturation constant for iron uptake than other phytoplankton. Simulation results suggested that P. antarctica had a competitive advantage under low irradiance levels, especially in the Ross Sea and Weddell Sea. However, the distributions of P. antarctica and diatoms were also strongly influenced by iron availability. Although grazing rates had an influence on total biomass, our simulations did not show a strong influence of grazing pressure in the competition between P. antarctica and diatoms. However, limited observations and the relative simplicity of zooplankton in our model suggest further research is needed. Overall, P. antarctica contributed ∼13% of annual primary production and ∼19% of sinking carbon export in the Southern Ocean (>40°S) in our best case simulation. At higher latitudes (>60°S) P. antarctica accounts for ∼23% of annual primary production and ∼30% of sinking carbon export.

Biogeochemical cycling in the Southern Ocean (SO) plays a key role in the global sea-air CO₂balance and in the ocean anthropogenic carbon inventory (Ito et al., 2010; Khatiwala et al., 2009; Sarmiento et al., 2004). Some previous studies suggest a decreasing trend in the Southern Ocean carbon sink (Le Quéré et al., 2007; Lovenduski et al., 2007; Wetzel et al., 2005). We investigate the interannual and decadal variations in sea-air CO₂flux and phytoplankton production in the SO with hindcast simulations by an ocean biogeochemical model. Decreasing trends in sinking POC and primary production are found from 1979 to 2003, concurrent with a decreasing trend in carbon uptake from the atmosphere. Simulations show substantial interannual and decadal variability in productivity. The sea-air CO₂flux is significantly correlated with sinking POC, especially in high productivity regions of the Southern Ocean. Both mixed layer depths and iron concentrations are important to the long-term trends in production and phytoplankton community structure. Sea ice cover also plays an important role at high latitudes. Variability in dust deposition in recent decades has little influence on total SO productivity and carbon uptake, however, there are regional impacts near dust source regions. Accurately representing mixed layer depths and their impacts on phytoplankton light stress are critical for understanding how climate change impacts SO ecosystems and biogeochemistry.

Organic macromolecules constitute a high percentage of remote sea spray. They enter the atmosphere through adsorption onto bubbles followed by bursting at the ocean surface, and go on to influence the chemistry of the fine mode aerosol. We present a global estimate of mixed-layer macromolecular distributions, driven by offline marine systems model output. The approach permits estimation of oceanic concentrations and bubble film surface coverages for several classes of organic compound. Mixed layer levels are computed from the output of a global ocean ecodynamics model by relating the macromolecules to standard biogeochemical tracers. Steady state is assumed for labile forms, and for longer-lived components we rely on ratios to existing transported variables. Adsorption is then represented through conventional Langmuir isotherms, with equilibria deduced from laboratory analogs. Open water concentrations locally exceed one micromolar carbon for the total of proteins, polysaccharides and refractory heteropolycondensates. The shorter-lived lipids remain confined to regions of strong biological activity. Results are evaluated against available measurements for all compound types, and agreement is generally well within an order of magnitude. Global distributions are further estimated for both fractional coverage of bubble films at the air–water interface and the two-dimensional concentration excess. Overall, we show that macromolecular mapping provides a novel tool for the comprehension of oceanic surfactant patterns. These results may prove useful in planning field experiments and assessing the potential response of surface chemical behaviors to global change.

This work documents the first version of the U.S. Department of Energy (DOE) new Energy Exascale Earth System Model (E3SMv1). We focus on the standard resolution of the fully coupled physical model designed to address DOE mission-relevant water cycle questions. Its components include atmosphere and land (110-km grid spacing), ocean and sea ice (60 km in the midlatitudes and 30 km at the equator and poles), and river transport (55 km) models. This base configuration will also serve as a foundation for additional configurations exploring higher horizontal resolution as well as augmented capabilities in the form of biogeochemistry and cryosphere configurations. The performance of E3SMv1 is evaluated by means of a standard set of Coupled Model Intercomparison Project Phase 6 (CMIP6) Diagnosis, Evaluation, and Characterization of Klima simulations consisting of a long preindustrial control, historical simulations (ensembles of fully coupled and prescribed SSTs) as well as idealized CO2 forcing simulations. The model performs well overall with biases typical of other CMIP-class models, although the simulated Atlantic Meridional Overturning Circulation is weaker than many CMIP-class models. While the E3SMv1 historical ensemble captures the bulk of the observed warming between preindustrial (1850) and present day, the trajectory of the warming diverges from observations in the second half of the twentieth century with a period of delayed warming followed by an excessive warming trend. Using a two-layer energy balance model, we attribute this divergence to the model's strong aerosol-related effective radiative forcing (ERFari+aci = −1.65 W/m2) and high equilibrium climate sensitivity (ECS = 5.3 K).