Optimality-based approach for computationally efficient modeling of phytoplankton growth, chlorophyll-to-carbon, and nitrogen-to-carbon ratios

Chen, Bingzhang and Smith, Sherwood Lan (2018) Optimality-based approach for computationally efficient modeling of phytoplankton growth, chlorophyll-to-carbon, and nitrogen-to-carbon ratios. Ecological Modelling, 385. pp. 197-212. ISSN 0304-3800 (https://doi.org/10.1016/j.ecolmodel.2018.08.001)

[thumbnail of Chen-Smith-EM-2018-Optimality-based-approach-for-computationally-efficient-modeling-of-phytoplankton-growth]
Preview
Text. Filename: Chen_Smith_EM_2018_Optimality_based_approach_for_computationally_efficient_modeling_of_phytoplankton_growth.pdf
Accepted Author Manuscript
License: Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 logo

Download (804kB)| Preview

Abstract

To increase the efficiency of computing phytoplankton growth rate (μ), chlorophyll-to-carbon (θ) and nitrogen- to-carbon ratios (QN) in three-dimensional ocean circulation models, it is preferable to directly calculate θ and QN from ambient environmental factors instead of treating them as independent tracers. Optimality-based modeling has emerged as a novel and efficient approach to fulfill this task. However, it is still unclear precisely how the response of optimality-based models differs from conventional models. We compare a recent optimality- based phytoplankton model (PAHLOW model), based on which the familiar Droop function can be derived, to a commonly used Monod-type (MONOD) model. The two models generate similar patterns of μ with some im- portant differences. Compared to the MONOD model, the PAHLOW model predicts higher μ under light lim- itation. The PAHLOW model also predicts that θ decreases with decreasing light under dim light and predicts decreasing QN with increasing light even at constant nutrient levels. Compared to the MONOD model, these features of the PAHLOW model qualitatively agree better with laboratory data. The PAHLOW model also suffers from a few shortcomings including the underestimation of θ under very low light and two times of computation time compared to the MONOD model. The two models generate striking differences of QN and θ in a one- dimensional implementation. Validation of such patterns will require more direct in situ measurements of μ, θ and QN.