Knowledge of natural climate variability is essential to better constrain the uncertainties in projections of twenty-first-century climate change 1–5. The past 2,000 years (2 kyr) have emerged as a critical interval in this endeavour, with sufficient length to characterize natural decadal-to-centennial scale change, known external climate forcings 6 and with distinctive patterns of spatiotemporal temperature variations 7. However, reconstructions for the full 2 kyr interval are not available for the global ocean, a primary heat reservoir 8 and an important regulator of global climate on longer timescales 9–11. Here we present a global ocean sea surface temperature (SST) synthesis (Ocean2k SST synthesis) spanning the Common Era, which shows a cooling trend that is similar, within uncertainty, to that simulated by realistically forced climate models for the past millennium. We use the simulations to identify the climate forcing(s) consistent with reconstructed SST variations during the past millennium. The oceans mediate the response of global climate to natural and anthropogenic forcings. Yet for the past 2,000 years — a key interval for understanding the present and future climate response to these forcings — global sea surface temperature changes and the underlying driving mechanisms are poorly constrained. Here we present a global synthesis of sea surface temperatures for the Common Era (ce) derived from 57 individual marine reconstructions that meet strict quality control criteria. We observe a cooling trend from 1 to 1800 ce that is robust against explicit tests for potential biases in the reconstructions. Between 801 and 1800 ce, the surface cooling trend is qualitatively consistent with an independent synthesis of terrestrial temperature reconstructions, and with a sea surface temperature composite derived from an ensemble of climate model simulations using best estimates of past external radiative forcings. Climate simulations using single and cumulative forcings suggest that the ocean surface cooling trend from 801 to 1800 ce is not primarily a response to orbital forcing but arises from a high frequency of explosive volcanism. Our results show that repeated clusters of volcanic eruptions can induce a net negative radiative forcing that results in a centennial and global scale cooling trend via a decline in mixed-layer oceanic heat content.