Relationships between optical backscattering, particulate organic carbon, and phytoplankton carbon in the oligotrophic South China Sea basin

Biogeochemical Model Optimization by Using Satellite-Derived Phytoplankton Functional Type Data and BGC-Argo Observations in the Northern South China Sea

Marine biogeochemical models have been widely used to understand ecosystem dynamics and biogeochemical cycles. To resolve more processes, models typically increase in complexity, and require optimization of more parameters. Data assimilation is an essential tool for parameter optimization, which can reduce model uncertainty and improve model predictability. At present, model parameters are often adjusted using sporadic in-situ measurements or satellite-derived total chlorophyll-a concentration at sea surface. However, new ocean datasets and satellite products have become available, providing a unique opportunity to further constrain ecosystem models. Biogeochemical-Argo (BGC-Argo) floats are able to observe the ocean interior continuously and satellite phytoplankton functional type (PFT) data has the potential to optimize biogeochemical models with multiple phytoplankton species. In this study, we assess the value of assimilating BGC-Argo measurements and satellite-derived PFT data in a biogeochemical model in the northern South China Sea (SCS) by using a genetic algorithm. The assimilation of the satellite-derived PFT data was found to improve not only the modeled total chlorophyll-a concentration, but also the individual phytoplankton groups at surface. The improvement of simulated surface diatom provided a better representation of subsurface particulate organic carbon (POC). However, using satellite data alone did not improve vertical distributions of chlorophyll-a and POC. Instead, these distributions were improved by combining the satellite data with BGC-Argo data. As the dominant variability of phytoplankton in the northern SCS is at the seasonal timescale, we find that utilizing monthly-averaged BGC-Argo profiles provides an optimal fit between model outputs and measurements in the region, better than using high-frequency measurements.