Site-Specific Multiple Stressor Assessments Based on High Frequency Surface Observations and an Earth System Model

Global Earth system models are often enlisted to assess the impacts of climate variability and change on marine ecosystems. In this study, we compare high frequency (daily) outputs of potential ecosystem stressors, such as sea surface temperature and surface pH, and associated variables from an Earth system model (GFDL ESM4.1) with high frequency time series from a global network of moorings to directly assess the capacity of the model to resolve local biogeochemical variability on time scales from daily to interannual. Our analysis indicates variability in surface temperature is most consistent between ESM4.1 and observations, with a Pearson correlation coefficient of 0.93 and bias of 0.40°C, followed by variability in surface salinity. Physical variability is reproduced with greater accuracy than biogeochemical variability, and variability on seasonal and longer time scales is more consistent between the model and observations than higher frequency variability. At the same time, the well-resolved seasonal and longer timescale variability is a reasonably good predictor, in many cases, of the likelihood of extreme events. Despite limited model representation of high frequency variability, model and observation-based assessments of the fraction of days experiencing surface T-pH and T-Ω_arag multistressor conditions show reasonable agreement, depending on the stressor combination and threshold definition. We also identify circumstances in which some errors could be reduced by accounting for model biases.