Krause, J. W., Stukel, M. R., Taylor, A. G., Taniguchi, D. A. A., De Verneil, A., & Landry, M. R. (2016). Net biogenic silica production and the contribution of diatoms to new production and organic matter export in the Costa Rica Dome ecosystem. J Plankton Res , 38 (2), 216–229.
Abstract: We determined the net rate of biogenic silica (bSiO2) production and estimated the diatom contribution to new production and organic matter export in the Costa Rica Dome during summer 2010. The shallow thermocline significantly reduces bSiO2 dissolution rates below the mixed layer, leading to significant enhancement of bSiO2 relative to organic matter (silicate-pump condition). This may explain why deep export of bSiO2 in this region is elevated by an order of magnitude relative to comparable systems. Diatom carbon, relative to autotrophic carbon, was low (<3%); however, the contribution of diatoms to new production averaged 3 and 13% using independent approaches. The 4-old discrepancy between methods may be explained by a low average C:Si ratio ( approximately 1.4) for the net produced diatom C relative to the net produced bSiO2. We speculate that this low production ratio is not the result of reduced C, but may arise from a significant contribution of non-diatom silicifying organisms to bSiO2 production. The contribution of diatoms to organic matter export was minor (5.7%). These results, and those of the broader project, suggest substantial food-web transformation of diatom organic matter in the euphotic zone, which creates enriched bSiO2 relative to organic matter within the exported material.
Selph, K. E., Landry, M. R., Taylor, A. G., Gutierrez-Rodriguez, A., Stukel, M. R., Wokuluk, J., et al. (2016). Phytoplankton production and taxon-specific growth rates in the Costa Rica Dome. J Plankton Res , 38 (2), 199–215.
Abstract: During summer 2010, we investigated phytoplankton production and growth rates at 19 stations in the eastern tropical Pacific, where winds and strong opposing currents generate the Costa Rica Dome (CRD), an open-ocean upwelling feature. Primary production (14C-incorporation) and group-specific growth and net growth rates (two-treatment seawater dilution method) were estimated from samples incubated in situ at eight depths. Our cruise coincided with a mild El Nino event, and only weak upwelling was observed in the CRD. Nevertheless, the highest phytoplankton abundances were found near the dome center. However, mixed-layer growth rates were lowest in the dome center ( approximately 0.5-0.9 day-1), but higher on the edge of the dome ( approximately 0.9-1.0 day-1) and in adjacent coastal waters (0.9-1.3 day-1). We found good agreement between independent methods to estimate growth rates. Mixed-layer growth rates of Prochlorococcus and Synechococcus were largely balanced by mortality, whereas eukaryotic phytoplankton showed positive net growth ( approximately 0.5-0.6 day-1), that is, growth available to support larger (mesozooplankton) consumer biomass. These are the first group-specific phytoplankton rate estimates in this region, and they demonstrate that integrated primary production is high, exceeding 1 g C m-2 day-1 on average, even during a period of reduced upwelling.
Conroy, B. J., Steinberg, D. K., Stukel, M. R., Goes, J. I., & Coles, V. J. (2016). Meso- and microzooplankton grazing in the Amazon River plume and western tropical North Atlantic. Limnol. Oceanogr. , 61 (3), 825–840.
Stukel, M. R., & Ducklow, H. W. (2017). Stirring Up the Biological Pump: Vertical Mixing and Carbon Export in the Southern Ocean. Global Biogeochem. Cycles , 31 (9), 1420–1434.
Stukel, M. R., Kahru, M., Benitez-Nelson, C. R., Décima, M., Goericke, R., Landry, M. R., et al. (2015). Using Lagrangian-based process studies to test satellite algorithms of vertical carbon flux in the eastern North Pacific Ocean. J. Geophys. Res. Oceans , 120 (11), 7208–7222.
Brzezinski, M. A., Krause, J. W., Bundy, R. M., Barbeau, K. A., Franks, P., Goericke, R., et al. (2015). Enhanced silica ballasting from iron stress sustains carbon export in a frontal zone within the California Current. J. Geophys. Res. Oceans , 120 (7), 4654–4669.
Stukel, M. R., Decima, M., & Kelly, T. B. (2018). A new approach for incorporating 15N isotopic data into linear inverse ecosystem models with Markov Chain Monte Carlo sampling. PLoS One , 13 (6), e0199123.
Abstract: Oceanographic field programs often use delta15N biogeochemical measurements and in situ rate measurements to investigate nitrogen cycling and planktonic ecosystem structure. However, integrative modeling approaches capable of synthesizing these distinct measurement types are lacking. We develop a novel approach for incorporating delta15N isotopic data into existing Markov Chain Monte Carlo (MCMC) random walk methods for solving linear inverse ecosystem models. We test the ability of this approach to recover food web indices (nitrate uptake, nitrogen fixation, zooplankton trophic level, and secondary production) derived from forward models simulating the planktonic ecosystems of the California Current and Amazon River Plume. We show that the MCMC with delta15N approach typically does a better job of recovering ecosystem structure than the standard MCMC or L2 minimum norm (L2MN) approaches, and also outperforms an L2MN with delta15N approach. Furthermore, we find that the MCMC with delta15N approach is robust to the removal of input equations and hence is well suited to typical pelagic ecosystem studies for which the system is usually vastly under-constrained. Our approach is easily extendable for use with delta13C isotopic measurements or variable carbon:nitrogen stoichiometry.
Coles, V. J., Stukel, M. R., Brooks, M. T., Burd, A., Crump, B. C., Moran, M. A., et al. (2017). Ocean biogeochemistry modeled with emergent trait-based genomics. Science , 358 (6367), 1149–1154.
Abstract: Marine ecosystem models have advanced to incorporate metabolic pathways discovered with genomic sequencing, but direct comparisons between models and “omics” data are lacking. We developed a model that directly simulates metagenomes and metatranscriptomes for comparison with observations. Model microbes were randomly assigned genes for specialized functions, and communities of 68 species were simulated in the Atlantic Ocean. Unfit organisms were replaced, and the model self-organized to develop community genomes and transcriptomes. Emergent communities from simulations that were initialized with different cohorts of randomly generated microbes all produced realistic vertical and horizontal ocean nutrient, genome, and transcriptome gradients. Thus, the library of gene functions available to the community, rather than the distribution of functions among specific organisms, drove community assembly and biogeochemical gradients in the model ocean.
