Marine diatoms are one of the marine phytoplankton functional groups, with high species diversity, playing important roles in the marine food web and carbon sequestration. In order to evaluate the species-specific responses of coastal diatoms to the combined effects of future ocean acidification (OA) and warming on the coastal diatoms, we conducted a semi-continuous incubation on the large centric diatom Thalassiosira sp. (~30 μm) and small pennate diatom Nitzschia closterium f.minutissima (~15 μm). A full factorial combination of two temperature levels (15 and 20°C) and pCO2 (400 and 1,000 ppm) was examined. The results suggest that changes in temperature played a more important role in regulating the physiology of Thalassiosira sp. and N. closterium f.minutissima than CO2. For Thalassiosira sp., elevated temperature significantly reduced the cellular particulate organic carbon (POC), particulate organic nitrogen (PON), particulate organic phosphate (POP), biogenic silica (BSi), chlorophyll a (Chl a), and protein contents, and the C:N ratio. CO2 only had significant effects on the growth rate and the protein content. However, for the smaller pennate diatom N. closterium f.minutissima, the growth rate, POC production rate, and the C:P ratio significantly increased with an elevated temperature, whereas the cellular POP and BSi contents significantly decreased. CO2 had significant effects on the POC production rate, cellular BSi, POC, and PON contents, the C:P, Si:C, N:P, and Si:P ratios, and sinking rate. The interaction between OA and warming showed mostly antagonistic effects on the physiology of both species. Overall, by comparison between the two species, CO2 played a more significant role in regulating the growth rate and sinking rate of the large centric diatom Thalassiosira sp., whereas had more significant effects on the elemental compositions of the smaller pennate diatom N. closterium f.minutissima. These results suggest differential sensitivities of different diatom species with different sizes and morphology to the changes in CO2/temperature regimes and their interactions.

In order to allow full comparability with other ocean acidification data sets, the R package seacarb (Gattuso et al, 2021) was used to compute a complete and consistent set of carbonate system variables, as described by Nisumaa et al. (2010). In this dataset the original values were archived in addition with the recalculated parameters (see related PI). The date of carbonate chemistry calculation by seacarb is 2022-07-22.