We present a new framework for global ocean–sea-ice model simulations based on phase 2 of the Ocean Model Intercomparison Project (OMIP-2), making use of the surface dataset based on the Japanese 55-year atmospheric reanalysis for driving ocean–sea-ice models (JRA55-do). We motivate the use of OMIP-2 over the framework for the first phase of OMIP (OMIP-1), previously referred to as the Coordinated Ocean–ice Reference Experiments (COREs), via the evaluation of OMIP-1 and OMIP-2 simulations from 11 state-of-the-science global ocean–sea-ice models. In the present evaluation, multi-model ensemble means and spreads are calculated separately for the OMIP-1 and OMIP-2 simulations and overall performance is assessed considering metrics commonly used by ocean modelers. Both OMIP-1 and OMIP-2 multi-model ensemble ranges capture observations in more than 80 % of the time and region for most metrics, with the multi-model ensemble spread greatly exceeding the difference between the means of the two datasets. Many features, including some climatologically relevant ocean circulation indices, are very similar between OMIP-1 and OMIP-2 simulations, and yet we could also identify key qualitative improvements in transitioning from OMIP-1 to OMIP-2. For example, the sea surface temperatures of the OMIP-2 simulations reproduce the observed global warming during the 1980s and 1990s, as well as the warming slowdown in the 2000s and the more recent accelerated warming, which were absent in OMIP-1, noting that the last feature is part of the design of OMIP-2 because OMIP-1 forcing stopped in 2009. A negative bias in the sea-ice concentration in summer of both hemispheres in OMIP-1 is significantly reduced in OMIP-2. The overall reproducibility of both seasonal and interannual variations in sea surface temperature and sea surface height (dynamic sea level) is improved in OMIP-2. These improvements represent a new capability of the OMIP-2 framework for evaluating process-level responses using simulation results. Regarding the sensitivity of individual models to the change in forcing, the models show well-ordered responses for the metrics that are directly forced, while they show less organized responses for those that require complex model adjustments. Many of the remaining common model biases may be attributed either to errors in representing important processes in ocean–sea-ice models, some of which are expected to be reduced by using finer horizontal and/or vertical resolutions, or to shared biases and limitations in the atmospheric forcing. In particular, further efforts are warranted to resolve remaining issues in OMIP-2 such as the warm bias in the upper layer, the mismatch between the observed and simulated variability of heat content and thermosteric sea level before 1990s, and the erroneous representation of deep and bottom water formations and circulations. We suggest that such problems can be resolved through collaboration between those developing models (including parameterizations) and forcing datasets. Overall, the present assessment justifies our recommendation that future model development and analysis studies use the OMIP-2 framework.

The impact of horizontal resolution (1/12° to 1/50°; 6 to 1.5 km at midlatitudes) on Gulf Stream separation, penetration, and variability is quantified in a series of identical North Atlantic experiments. The questions the authors seek to address are twofold: 1) Is the realism of the modeled solution increased as resolution is increased? 2) How robust is the modeled mesoscale and submesoscale eddy activity as a function of grid spacing and how representative is it of interior quasigeostrophic (QG) or surface quasigeostrophic (SQG) turbulence? This study shows that (i) the representation of Gulf Stream penetration and associated recirculating gyres shifts from unrealistic to realistic when the resolution is increased to 1/50° and when the nonlinear effects of the submesoscale eddies intensifies the midlatitude jet and increases its penetration eastward, (ii) the penetration into the deep ocean drastically increases with resolution and closely resembles the observations, and (iii) surface power spectra in the 70–250-km mesoscale range are independent of the horizontal resolution and of the latitude and are representative of 2D QG and SQG turbulence.

This paper presents global comparisons of fundamental global climate variables from a suite of four pairs of matched low- and high-resolution ocean and sea ice simulations that are obtained following the OMIP-2 protocol (Griffies et al., 2016) and integrated for one cycle (1958–2018) of the JRA55-do atmospheric state and runoff dataset (Tsujino et al., 2018). Our goal is to assess the robustness of climate-relevant improvements in ocean simulations (mean and variability) associated with moving from coarse (∼ 1^∘) to eddy-resolving (∼ 0.1^∘) horizontal resolutions. The models are diverse in their numerics and parameterizations, but each low-resolution and high-resolution pair of models is matched so as to isolate, to the extent possible, the effects of horizontal resolution. A variety of observational datasets are used to assess the fidelity of simulated temperature and salinity, sea surface height, kinetic energy, heat and volume transports, and sea ice distribution. This paper provides a crucial benchmark for future studies comparing and improving different schemes in any of the models used in this study or similar ones. The biases in the low-resolution simulations are familiar, and their gross features – position, strength, and variability of western boundary currents, equatorial currents, and the Antarctic Circumpolar Current – are significantly improved in the high-resolution models. However, despite the fact that the high-resolution models “resolve” most of these features, the improvements in temperature and salinity are inconsistent among the different model families, and some regions show increased bias over their low-resolution counterparts. Greatly enhanced horizontal resolution does not deliver unambiguous bias improvement in all regions for all models.

Plastic is the most abundant type of marine litter and it is found in all of the world’s oceans and seas, even in remote areas far from human activities. It is a major concern because plastics remain in the oceans for a long time. To address questions that are of great interest to the international community as it seeks to attend to the major sources of marine plastics in the ocean, we use particle tracking simulations to simulate the motions of mismanaged plastic waste and provide a quantitative global estimate of (1) where does the marine litter released into the ocean by a given country go and (2) where does the marine litter found on the coastline of a given country come from. The overall distribution of the modeled marine litter is in good agreement with the limited observations that we have at our disposal and our results illustrate how countries that are far apart are connected via a complex web of ocean pathways (see interactive websitehttps://marinelitter.coaps.fsu.edu). The tables summarizing the statistics for all world countries are accessible from the supplemental information in .pdf or .csv formats.

Marine plastic pollution poses a potential threat to the ecosystem, but the sources and their magnitudes remain largely unclear. Existing bottom-up emission inventories vary among studies for two to three orders of magnitudes (OMs). Here, we adopt a top-down approach that uses observed dataset of sea surface plastic concentrations and an ensemble of ocean transport models to reduce the uncertainty of global plastic discharge. The optimal estimation of plastic emissions in this study varies about 1.5 OMs: 0.70 (0.13–3.8 as a 95% confidence interval) million metric tons yr⁻¹ at the present day. We find that the variability of surface plastic abundance caused by different emission inventories is higher than that caused by model parameters. We suggest that more accurate emission inventories, more data for the abundance in the seawater and other compartments, and more accurate model parameters are required to further reduce the uncertainty of our estimate.

Following the 2010 Deepwater Horizon accident of a massive blow-out in the Gulf of Mexico, scientists from government, industry, and academia collaborated to advance oil spill modeling and share best practices in model algorithms, parameterizations, and application protocols. This synergy was greatly enhanced by research funded under the Gulf of Mexico Research Initiative (GoMRI), a 10-year enterprise that allowed unprecedented collection of observations and data products, novel experiments, and international collaborations that focused on the Gulf of Mexico, but resulted in the generation of scientific findings and tools of broader value. Operational oil spill modeling greatly benefited from research during the GoMRI decade. This paper provides a comprehensive synthesis of the related scientific advances, remaining challenges, and future outlook. Two main modeling components are discussed: Ocean circulation and oil spill models, to provide details on all attributes that contribute to the success and limitations of the integrated oil spill forecasts. These forecasts are discussed in tandem with uncertainty factors and methods to mitigate them. The paper focuses on operational aspects of oil spill modeling and forecasting, including examples of international operational center practices, observational needs, communication protocols, and promising new methodologies.