We assembled a dataset of 14C-based productivity measurements to understand the critical variables required for accurate assessment of daily depth-integrated phytoplankton carbon fixation (PP(PPeu)u) from measurements of sea surface pigment concentrations (Csat)(Csat). From this dataset, we developed a light-dependent, depth-resolved model for carbon fixation (VGPM) that partitions environmental factors affecting primary production into those that influence the relative vertical distribution of primary production (Pz)z) and those that control the optimal assimilation efficiency of the productivity profile (P(PBopt). The VGPM accounted for 79% of the observed variability in Pz and 86% of the variability in PPeu by using measured values of PBopt. Our results indicate that the accuracy of productivity algorithms in estimating PPeu is dependent primarily upon the ability to accurately represent variability in Pbopt. We developed a temperature-dependent Pbopt model that was used in conjunction with monthly climatological images of Csat sea surface temperature, and cloud-corrected estimates of surface irradiance to calculate a global annual phytoplankton carbon fixation (PPannu) rate of 43.5 Pg C yr‒1. The geographical distribution of PPannu was distinctly different than results from previous models. Our results illustrate the importance of focusing Pbopt model development on temporal and spatial, rather than the vertical, variability.

Species distribution models (Maxent) predicting the distribution of two Vulnerable Marine Ecosystems (VME): the reef-forming Scleractinian coral Desmophyllum pertusum and the aggregations forming Hexactinellid sponge Pheronema carpenteri. Both of these species are VME indicator taxa and form habitat that enhance deep-sea diversity (Ross and Howell, 2013). Maps of the likely distribution of the habitat formed by these two species will enable efficient Marine Spatial planning to facilitate their conservation. This work was performed at the University of Plymouth in 2021. A GIS layer is provided for each species.

This dataset comprises the global frequency, classification and distribution of marine heat waves (MHWs) from 1996-2020, in Chauhan et al. 2023 (https://doi.org/10.3389/fmars.2023.1177571). The classification was done based on their attributes and using different baselines. Daily SST values were extracted from the NOAA-OISST v2 high-resolution (0.25°) dataset from 1982-2020. MHWs were detected using the method presented by Hobday et al. 2016 (https://doi.org/10.1016/j.pocean.2015.12.014), and by using the 95th percentile of the accumulated temperature distribution to flag the extreme events. A shifting baseline of 8-year rolling period was selected between the years 1982-1996, since this period shows relatively stable maximum values of temperature across different ocean regions. The shifting baseline aims to account for the decadal changes of westerly winds, temperatures and ocean gyres circulations. The classification was done using the KMeans clustering algorithm to identify the relevant features of MHWs and classify them into separate groups based on feature similarities. This algorithm takes MHW features, namely duration, maximum intensity, rate onset and rate decline, as input vectors and applies clustering in the 4-dimensional feature space where each data point represents an MHW event. Note that all the MHWs features are standardized because unequal variances can put more weight on variables with smaller variances. This record comprehends the geospatial datasets of: Average number of MHW days per year (i.e., the sum of all MHW days divided by the total number of years, 1996-2020). Average cumulative intensity per year (i.e., the sum of cumulative intensity divided by the total number of years, 1996-2020). Total number of MHW events across the different periods averaged on the total number of years (1989-2020). The period 1982-1988 was only used as an initial baseline without calculating MHWs. Spatial distribution of three MHW categories: moderate MHWs, abrupt and Intense MHWs and extreme MHWs; displaying the total number of MHW days normalized by the number of years considered (i.e., 1989-2020). Distribution of Extreme MHWs across the different periods (A) 1989-1996, (B) 1997-2004, (C) 2005-2012, (D) 2013-2020. The relative frequency (γ) is a ratio of extreme MHWs in a specific period and all extreme MHWs in the same cluster for all periods.