Search

The effect of biomass dynamics on the estimation of watercolumn primary production is analyzed, by coupling a primary production model to a simple growth equation for phytoplankton. The production model is formulated with depth- and time-resolved biomass, and placed in the context of earlier models, with emphasis on the canonical solution for watercolumn production. A relation between the canonical solution and the general solution for the case of an arbitrary depth-dependent biomass profile was derived, together with an analytical solution for watercolumn production in case of a depth dependent biomass profile described with the shifted Gaussian function. The analysis was further extended to the case of a time-dependent, mixed-layer biomass, and two additional analytical solutions to this problem were derived, the first in case of increasing mixed-layer biomass and the second in case of declining biomass. The solutions were tested with Hawaii Ocean Time-series data. The canonical solution for mixed-layer production has proven to be a good model for this data set. The shifted Gaussian function was demonstrated to be an accurate model for the measured biomass profiles and the shifted Gaussian parameters extracted from the measured profiles were further used in the analytical solution for watercolumn production and results compared with data. The influence of time-dependent biomass on mixed-layer production was studied through analytical solutions. Re-examining the Critical Depth Hypothesis we derived an expression for the daily increase in mixed-layer biomass. Finally, the work was placed in a remote sensing context and the time-dependent model for biomass related to the remotely sensed-biomass.

The rate of carbon assimilation in phytoplankton primary production models is mathematically prescribed with photosynthesis irradiance functions, which convert a light flux (energy) into a material flux (carbon). Information on this rate is contained in photosynthesis parameters: the initial slope and the assimilation number. The exactness of parameter values is crucial for precise calculation of primary production. Here we use a model of the daily production profile based on a suite of photosynthesis irradiance functions and extract photosynthesis parameters from in situ measured daily production profiles at the Hawaii Ocean Time-series station Aloha. For each function we recover parameter values, establish parameter distributions and quantify model skill. We observe that the choice of the photosynthesis irradiance function to estimate the photosynthesis parameters affects the magnitudes of parameter values as recovered from in situ profiles. We also tackle the problem of parameter exchange amongst the models and the effect it has on model performance. All models displayed little or no bias prior to parameter exchange, but significant bias following parameter exchange. The best model performance resulted from using optimal parameter values. Model formulation was extended further by accounting for spectral effects and deriving a spectral analytical solution for the daily production profile. The daily production profile was also formulated with time dependent growing biomass governed by a growth equation. The work on parameter recovery was further extended by exploring how to extract photosynthesis parameters from information on watercolumn production. It was demonstrated how to estimate parameter values based on a linearization of the full analytical solution for normalized watercolumn production and from the solution itself, without linearization. The paper complements previous works on photosynthesis irradiance models by analysing the skill and consistency of photosynthesis irradiance functions and parameters for modeling in situ production profiles. In light of the results obtained in this work we argue that the choice of the primary production model should reflect the available data and these models should be data driven regarding parameter estimation.