Ocean Model Intercomparison Project Research Articles

Experimental design for the Marine Ice Sheet–Ocean Model Intercomparison Project – phase 2 (MISOMIP2)

Abstract. The Marine Ice Sheet–Ocean Model Intercomparison Project – phase 2 (MISOMIP2) is a natural progression of previous and ongoing model intercomparison exercises that have focused on the simulation of ice-sheet and ocean processes in Antarctica. The previous exercises motivate the move towards realistic configurations, as well as more diverse model parameters and resolutions. The main objective of MISOMIP2 is to investigate the performance of existing ocean and coupled ice-sheet–ocean models in a range of Antarctic environments through comparisons to observational data. We will assess the status of ice-sheet–ocean modelling as a community and identify common characteristics of models that are best able to capture observed features. As models are highly tuned based on present-day data, we will also compare their sensitivity to prescribed abrupt atmospheric perturbations leading to either very warm or slightly warmer ocean conditions compared to the present day. The approach of MISOMIP2 is to welcome contributions of models as they are, including global and regional configurations, but we request standardized variables and common grids for the outputs. We target the analysis at two specific regions, the Amundsen Sea and the Weddell Sea, since they describe two different ocean environments and have been relatively well observed compared to other areas of Antarctica. An observational “MIPkit” synthesizing existing ocean and ice-sheet observations for a common period is provided to evaluate ocean and ice-sheet models in these two regions.

An assessment of equatorial Atlantic interannual variability in Ocean Model Intercomparison Project (OMIP) simulations

An assessment of equatorial Atlantic interannual variability in Ocean Model Intercomparison Project (OMIP) simulations

Impact of increased resolution on Arctic Ocean simulations in Ocean Model Intercomparison Project phase 2 (OMIP-2)

Abstract. This study evaluates the impact of increasing resolution on Arctic Ocean simulations using five pairs of matched low- and high-resolution models within the OMIP-2 (Ocean Model Intercomparison Project phase 2) framework. The primary objective is to assess whether a higher resolution can mitigate typical biases in low-resolution models and improve the representation of key climate-relevant variables. We reveal that increasing the horizontal resolution contributes to a reduction in biases in mean temperature and salinity and improves the simulation of the Atlantic water layer and its decadal warming events. A higher resolution also leads to better agreement with observed surface mixed-layer depth, cold halocline base depth and Arctic gateway transports in the Fram and Davis straits. However, the simulation of the mean state and temporal changes in Arctic freshwater content does not show improvement with increased resolution. Not all models achieve improvements for all analyzed ocean variables when spatial resolution is increased so it is crucial to recognize that model numerics and parameterizations also play an important role in faithful simulations. Overall, a higher resolution shows promise in improving the simulation of key Arctic Ocean features and processes, but efforts in model development are required to achieve more accurate representations across all climate-relevant variables.

Simulations of ocean deoxygenation in the historical era: insights from forced and coupled models

Ocean deoxygenation due to anthropogenic warming represents a major threat to marine ecosystems and fisheries. Challenges remain in simulating the modern observed changes in the dissolved oxygen (O2). Here, we present an analysis of upper ocean (0-700m) deoxygenation in recent decades from a suite of the Coupled Model Intercomparison Project phase 6 (CMIP6) ocean biogeochemical simulations. The physics and biogeochemical simulations include both ocean-only (the Ocean Model Intercomparison Project Phase 1 and 2, OMIP1 and OMIP2) and coupled Earth system (CMIP6 Historical) configurations. We examine simulated changes in the O2 inventory and ocean heat content (OHC) over the past 5 decades across models. The models simulate spatially divergent evolution of O2 trends over the past 5 decades. The trend (multi-model mean and spread) for upper ocean global O2 inventory for each of the MIP simulations over the past 5 decades is 0.03 ± 0.39×1014 [mol/decade] for OMIP1, −0.37 ± 0.15×1014 [mol/decade] for OMIP2, and −1.06 ± 0.68×1014 [mol/decade] for CMIP6 Historical, respectively. The trend in the upper ocean global O2 inventory for the latest observations based on the World Ocean Database 2018 is −0.98×1014 [mol/decade], in line with the CMIP6 Historical multi-model mean, though this recent observations-based trend estimate is weaker than previously reported trends. A comparison across ocean-only simulations from OMIP1 and OMIP2 suggests that differences in atmospheric forcing such as surface wind explain the simulated divergence across configurations in O2 inventory changes. Additionally, a comparison of coupled model simulations from the CMIP6 Historical configuration indicates that differences in background mean states due to differences in spin-up duration and equilibrium states result in substantial differences in the climate change response of O2. Finally, we discuss gaps and uncertainties in both ocean biogeochemical simulations and observations and explore possible future coordinated ocean biogeochemistry simulations to fill in gaps and unravel the mechanisms controlling the O2 changes.

Evaluation of the CMCC global eddying ocean model for the Ocean Model Intercomparison Project (OMIP2)

Abstract. This paper describes the global eddying ocean–sea ice simulation produced at the Euro-Mediterranean Center on Climate Change (CMCC) obtained following the experimental design of the Ocean Model Intercomparison Project phase 2 (OMIP2). The eddy-rich model (GLOB16) is based on the NEMOv3.6 framework, with a global horizontal resolution of 1/16∘ and 98 vertical levels and was originally designed for an operational short-term ocean forecasting system. Here, it is driven by one multi-decadal cycle of the prescribed JRA55-do atmospheric reanalysis and runoff dataset in order to perform a long-term benchmarking experiment. To assess the accuracy of simulated 3D ocean fields and highlight the relative benefits of resolving mesoscale processes, the GLOB16 performances are evaluated via a selection of key climate metrics against observational datasets and two other NEMO configurations at lower resolutions: an eddy-permitting resolution (ORCA025) and a non-eddying resolution (ORCA1) designed to form the ocean–sea ice component of the fully coupled CMCC climate model. The well-known biases in the low-resolution simulations are significantly improved in the high-resolution model. The evolution and spatial pattern of large-scale features (such as sea surface temperature biases and winter mixed-layer structure) in GLOB16 are generally better reproduced, and the large-scale circulation is remarkably improved compared to the low-resolution oceans. We find that eddying resolution is an advantage in resolving the structure of western boundary currents, the overturning cells, and flow through key passages. GLOB16 might be an appropriate tool for ocean climate modeling efforts, even though the benefit of eddying resolution does not provide unambiguous advances for all ocean variables in all regions.

Decadal variability of oxygen uptake, export, and storage in the Labrador Sea from observations and CMIP6 models

The uptake of dissolved oxygen from the atmosphere via air-sea gas exchange and its physical transport away from the region of uptake are crucial for supplying oxygen to the deep ocean. This process takes place in a few key regions that feature intense oxygen uptake, deep water formation, and physical oxygen export. In this study we analyze one such region, the Labrador Sea, utilizing the World Ocean Database (WOD) to construct a 65–year oxygen content time series in the Labrador Sea Water (LSW) layer (0–2200 m). The data reveal decadal variability associated with the strength of deep convection, with a maximum anomaly of 27 mol m–2 in 1992. There is no long-term trend in the time series, suggesting that the mean oxygen uptake is balanced by oxygen export out of the region. We compared the time series with output from nine models of the Ocean Model Intercomparison Project phase 1 in the Climate Model Intercomparison Project phase 6, (CMIP6-OMIP1), and constructed a “model score” to evaluate how well they match oxygen observations. Most CMIP6-OMIP1 models score around 50/100 points and the highest score is 57/100 for the ensemble mean, suggesting that improvements are needed. All of the models underestimate the maximum oxygen content anomaly in the 1990s. One possible cause for this is the representation of air-sea gas exchange for oxygen, with all models underestimating the mean uptake by a factor of two or more. Unrealistically deep convection and biased mean oxygen profiles may also contribute to the mismatch. Refining the representation of these processes in climate models could be vital for enhanced predictions of deoxygenation. In the CMIP6-OMIP1 multi-model mean, oxygen uptake has its maximum in 1980–1992, followed by a decrease in 1994–2006. There is a concurrent decrease in export, but oxygen storage also changes between the two periods, with oxygen accumulated in the first period and drained out in the second. Consequently, the change in oxygen export (5%) is much less than that in uptake (28%), suggesting that newly ventilated LSW which remains in the formation region acts to buffer the linkage between air-sea gas exchange and oxygen export.

The mixed-layer depth in the Ocean Model Intercomparison Project (OMIP): impact of resolving mesoscale eddies

Abstract. The ocean mixed layer is the interface between the ocean interior and the atmosphere or sea ice and plays a key role in climate variability. It is thus critical that numerical models used in climate studies are capable of a good representation of the mixed layer, especially its depth. Here we evaluate the mixed-layer depth (MLD) in six pairs of non-eddying (1∘ grid spacing) and eddy-rich (up to 1/16∘) models from the Ocean Model Intercomparison Project (OMIP), forced by a common atmospheric state. For model evaluation, we use an updated MLD dataset computed from observations using the OMIP protocol (a constant density threshold). In winter, low-resolution models exhibit large biases in the deep-water formation regions. These biases are reduced in eddy-rich models but not uniformly across models and regions. The improvement is most noticeable in the mode-water formation regions of the Northern Hemisphere. Results in the Southern Ocean are more contrasted, with biases of either sign remaining at high resolution. In eddy-rich models, mesoscale eddies control the spatial variability in MLD in winter. Contrary to a hypothesis that the deepening of the mixed layer in anticyclones would make the MLD larger globally, eddy-rich models tend to have a shallower mixed layer at most latitudes than coarser models do. In addition, our study highlights the sensitivity of the MLD computation to the choice of a reference level and the spatio-temporal sampling, which motivates new recommendations for MLD computation in future model intercomparison projects.

Impact of atmospheric forcing uncertainties on Arctic and Antarctic sea ice simulations in CMIP6 OMIP models

Abstract. Atmospheric reanalyses are valuable datasets for driving ocean–sea ice general circulation models and for proposing multidecadal reconstructions of the ocean–sea ice system in polar regions. However, these reanalyses exhibit biases in these regions. It was previously found that the representation of Arctic and Antarctic sea ice in models participating in the Ocean Model Intercomparison Project Phase 2 (OMIP2, using the updated Japanese 55-year atmospheric reanalysis, JRA55-do) was significantly more realistic than in OMIP1 (forced by the atmospheric state from the Coordinated Ocean-ice Reference Experiments version 2, CORE-II). To understand why, we study the sea ice concentration budget and its relations to surface heat and momentum fluxes as well as the connections between the simulated ice drift and the ice concentration, the ice thickness and the wind stress in a subset of three models (CMCC-CM2-SR5, MRI-ESM2-0 and NorESM2-LM). These three models are representative of the ensemble and are the only ones to provide the surface fluxes and the tendencies of ice concentrations attributed to dynamic and thermodynamic processes required for the ice concentration budget analysis. The sea ice simulations of two other models (EC-Earth3 and MIROC6) forced by both CORE-II and JRA55-do reanalysis are also included in the analysis. It is found that negative summer biases in high-ice-concentration regions and positive biases in the Canadian Arctic Archipelago (CAA) and central Weddell Sea (CWS) regions are reduced from OMIP1 to OMIP2 due to surface heat flux changes. Net shortwave radiation fluxes provide key improvements in the Arctic interior, CAA and CWS regions. There is also an influence of improved surface wind stress in OMIP2 giving better winter Antarctic ice concentration and the Arctic ice drift magnitude simulations near the ice edge. The ice velocity direction simulations in the Beaufort Gyre and the Pacific and Atlantic sectors of the Southern Ocean in OMIP2 are also improved owing to surface wind stress changes. This study provides clues on how improved atmospheric reanalysis products influence sea ice simulations. Our findings suggest that attention should be paid to the radiation fluxes and winds in atmospheric reanalyses in polar regions.

FIO-ESM v2.0 CORE2-forced experiment for the CMIP6 Ocean Model Intercomparison Project

FIO-ESM v2.0 CORE2-forced experiment for the CMIP6 Ocean Model Intercomparison Project

TIMCOM model datasets for the CMIP6 Ocean Model Intercomparison Project

The participation of the Taiwan Multi-scale Community Ocean Model (TIMCOM) in the Ocean Model Intercomparison Project (OMIP) experiments is introduced here, as part of phase 6 of the Coupled Model Intercomparison Project (CMIP6). Two ocean–sea ice model experiments are compared: (a) OMIP1, forced by the Coordinated Ocean-Ice Reference Experiments Phase II data (1948–2009), and (b) OMIP2, forced by JRA55-do data (1958–2018). The observed annual means and the interannual variability of physical states are reasonably captured in both experiments, but improved mean temperatures and salinities are found in OMIP2. The weaker winds and stronger freshwater discharge in the OMIP2 forcing contribute to some simulated differences between OMIP1 and OMIP2. Many patterns and biases are similar to those found in other modeling efforts, confirming the common systematic biases. However, a few unique features are found in this study, including the recent increase of the Atlantic Meridional Overturning Circulation (AMOC) that has been observed in the last decade and a generally higher Drake Passage transport. The enhanced AMOC can be explained by the recent cooling event over the North Atlantic, which thermally increased the surface density flux. The higher Drake Passage transport compared to observations is possibly linked to a stronger bottom cell of meridional circulation and a smaller Antarctic sea-ice extent.

Evaluation of an emergent feature of sub-shelf melt oscillations from an idealized coupled ice sheet–ocean model using FISOC (v1.1) – ROMSIceShelf (v1.0) – Elmer/Ice (v9.0)

Abstract. Changes in ocean-driven basal melting have a key influence on the stability of ice shelves, the mass loss from the ice sheet, ocean circulation, and global sea level rise. Coupled ice sheet–ocean models play a critical role in understanding future ice sheet evolution and examining the processes governing ice sheet responses to basal melting. However, as a new approach, coupled ice sheet–ocean systems come with new challenges, and the impacts of solutions implemented to date have not been investigated. An emergent feature in several contributing coupled models to the 1st Marine Ice Sheet–Ocean Model Intercomparison Project (MISOMIP1) was a time-varying oscillation in basal melt rates. Here, we use a recently developed coupling framework, FISOC (v1.1), to connect the modified ocean model ROMSIceShelf (v1.0) and ice sheet model Elmer/Ice (v9.0), to investigate the origin and implications of the feature and, more generally, the impact of coupled modeling strategies on the simulated basal melt in an idealized ice shelf cavity based on the MISOMIP setup. We found the spatial-averaged basal melt rates (3.56 m yr−1) oscillated with an amplitude ∼0.7 m yr−1 and approximate period of ∼6 years between year 30 and 100 depending on the experimental design. The melt oscillations emerged in the coupled system and the standalone ocean model using a prescribed change of cavity geometry. We found that the oscillation feature is closely related to the discretized ungrounding of the ice sheet, exposing new ocean, and is likely strengthened by a combination of positive buoyancy–melt feedback and/or melt–geometry feedback near the grounding line, and the frequent coupling of ice geometry and ocean evolution. Sensitivity tests demonstrate that the oscillation feature is always present, regardless of the choice of coupling interval, vertical resolution in the ocean model, tracer properties of cells ungrounded by the retreating ice sheet, or the dependency of friction velocities to the vertical resolution. However, the amplitude, phase, and sub-cycle variability of the oscillation varied significantly across the different configurations. We were unable to ultimately determine whether the feature arises purely due to numerical issues (related to discretization) or a compounding of multiple physical processes amplifying a numerical artifact. We suggest a pathway and choices of physical parameters to help other efforts understand the coupled ice sheet–ocean system using numerical models.

Subsurface Warm Biases in the Tropical Atlantic and Their Attributions to the Role of Wind Forcing and Ocean Vertical Mixing

Abstract Realistic ocean subsurface simulations of thermal structure and variation are critically important to success in climate prediction and projection; currently, substantial systematic subsurface biases still exist in the state-of-the-art ocean and climate models. In this paper, subsurface biases in the tropical Atlantic Ocean (TA) are investigated by analyzing simulations from the Ocean Model Intercomparison Project (OMIP) and conducting ocean-only experiments that are based on the Parallel Ocean Program, version 2 (POP2). The subsurface biases are prominent in almost all OMIP simulations, characterized by two warm-bias patches off the equator. By conducting two groups of POP2-based ocean-only experiments, two potential origins of the biases are explored, including uncertainties in wind forcing and vertical mixing parameterization, respectively. It is illustrated that the warm bias near 10°N can be slightly reduced by modulating the prescribed wind field, and the warm biases over the entire basin are significantly reduced by reducing background diffusivity in the ocean interior in ways to match observations. By conducting a heat-budget analysis, it is found that the improved subsurface simulations are attributed to the enhanced cooling effect by constraining the vertical mixing diffusivity in terms of the observational estimate, implying that overestimation of vertical mixing is primarily responsible for the subsurface warm biases in the TA. Since the climate simulation is very sensitive to the vertical mixing parameterization, more accurate representations of ocean vertical mixing are clearly needed in ocean and climate models. Significance Statement The purpose of our study is to analyze the characteristics of subsurface temperature biases in the tropical Atlantic Ocean and to investigate the causes for the biases. This is important because subsurface biases greatly reduce the reliability of models in climate prediction and projection. It is found that significant subsurface warm biases arise in 100–150 m over the entire tropical Atlantic basin and the biases are mainly attributed to overestimated ocean vertical mixing. Our work highlights that subsurface ocean simulations are highly sensitive to vertical mixing parameterization, and further research is necessary for more accurate representations of ocean vertical mixing in ocean and climate modeling.

Ice Algae Model Intercomparison Project phase 2 (IAMIP2)

Abstract. Ice algae play a fundamental role in shaping sea-ice-associated ecosystems and biogeochemistry. This role can be investigated by field observations; however the influence of ice algae at the regional and global scales remains unclear due to limited spatial and temporal coverage of observations and because ice algae are typically not included in current Earth system models. To address this knowledge gap, we introduce a new model intercomparison project (MIP), referred to here as the Ice Algae Model Intercomparison Project phase 2 (IAMIP2). IAMIP2 is built upon the experience from its previous phase and expands its scope to global coverage (both Arctic and Antarctic) and centennial timescales (spanning the mid-20th century to the end of the 21st century). Participating models are three-dimensional regional and global coupled sea-ice–ocean models that incorporate sea-ice ecosystem components. These models are driven by the same initial conditions and atmospheric forcing datasets by incorporating and expanding the protocols of the Ocean Model Intercomparison Project, an endorsed MIP of the Coupled Model Intercomparison Project phase 6 (CMIP6). Doing so provides more robust estimates of model bias and uncertainty and consequently advances the science of polar marine ecosystems and biogeochemistry. A diagnostic protocol is designed to enhance the reusability of the model data products of IAMIP2. Lastly, the limitations and strengths of IAMIP2 are discussed in the context of prospective research outcomes.

SITool (v1.0) – a new evaluation tool for large-scale sea ice simulations: application to CMIP6 OMIP

Abstract. The Sea Ice Evaluation Tool (SITool) described in this paper is a performance metrics and diagnostics tool developed to evaluate the skill of Arctic and Antarctic model reconstructions of sea ice concentration, extent, edge location, drift, thickness, and snow depth. It is a Python-based software and consists of well-documented functions used to derive various sea ice metrics and diagnostics. Here, SITool version 1.0 (v1.0) is introduced and documented, and is then used to evaluate the performance of global sea ice reconstructions from nine models that provided sea ice output under the experimental protocols of the Coupled Model Intercomparison Project phase 6 (CMIP6) Ocean Model Intercomparison Project with two different atmospheric forcing datasets: the Coordinated Ocean-ice Reference Experiments version 2 (CORE-II) and the updated Japanese 55-year atmospheric reanalysis (JRA55-do). Two sets of observational references for the sea ice concentration, thickness, snow depth, and ice drift are systematically used to reflect the impact of observational uncertainty on model performance. Based on available model outputs and observational references, the ice concentration, extent, and edge location during 1980–2007, as well as the ice thickness, snow depth, and ice drift during 2003–2007 are evaluated. In general, model biases are larger than observational uncertainties, and model performance is primarily consistent compared to different observational references. By changing the atmospheric forcing from CORE-II to JRA55-do reanalysis data, the overall performance (mean state, interannual variability, and trend) of the simulated sea ice areal properties in both hemispheres, as well as the mean ice thickness simulation in the Antarctic, the mean snow depth, and ice drift simulations in both hemispheres are improved. The simulated sea ice areal properties are also improved in the model with higher spatial resolution. For the cross-metric analysis, there is no link between the performance in one variable and the performance in another. SITool is an open-access version-controlled software that can run on a wide range of CMIP6-compliant sea ice outputs. The current version of SITool (v1.0) is primarily developed to evaluate atmosphere-forced simulations and it could be eventually extended to fully coupled models.

The GPU version of LASG/IAP Climate System Ocean Model version 3 (LICOM3) under the heterogeneous-compute interface for portability (HIP) framework and its large-scale application

Abstract. A high-resolution (1/20∘) global ocean general circulation model with graphics processing unit (GPU) code implementations is developed based on the LASG/IAP Climate System Ocean Model version 3 (LICOM3) under a heterogeneous-compute interface for portability (HIP) framework. The dynamic core and physics package of LICOM3 are both ported to the GPU, and three-dimensional parallelization (also partitioned in the vertical direction) is applied. The HIP version of LICOM3 (LICOM3-HIP) is 42 times faster than the same number of CPU cores when 384 AMD GPUs and CPU cores are used. LICOM3-HIP has excellent scalability; it can still obtain a speedup of more than 4 on 9216 GPUs compared to 384 GPUs. In this phase, we successfully performed a test of 1/20∘ LICOM3-HIP using 6550 nodes and 26 200 GPUs, and on a large scale, the model's speed was increased to approximately 2.72 simulated years per day (SYPD). By putting almost all the computation processes inside GPUs, the time cost of data transfer between CPUs and GPUs was reduced, resulting in high performance. Simultaneously, a 14-year spin-up integration following phase 2 of the Ocean Model Intercomparison Project (OMIP-2) protocol of surface forcing was performed, and preliminary results were evaluated. We found that the model results had little difference from the CPU version. Further comparison with observations and lower-resolution LICOM3 results suggests that the 1/20∘ LICOM3-HIP can reproduce the observations and produce many smaller-scale activities, such as submesoscale eddies and frontal-scale structures.

A mechanistic analysis of tropical Pacific dynamic sea level in GFDL-OM4 under OMIP-I and OMIP-II forcings

Abstract. The sea level over the tropical Pacific is a key indicator reflecting vertically integrated heat distribution over the ocean. Here, we use the Geophysical Fluid Dynamics Laboratory global ocean–sea ice model (GFDL-OM4) forced by both the Coordinated Ocean-Ice Reference Experiment (CORE) and Japanese 55-year Reanalysis (JRA-55)-based surface dataset for driving ocean–sea ice models (JRA55-do) atmospheric states (Ocean Model Intercomparison Project (OMIP) versions I and II) to evaluate the model performance and biases compared against available observations. We find persisting mean state dynamic sea level (DSL) bias along 9∘ N even with updated wind forcing in JRA55-do relative to CORE. The mean state bias is related to biases in wind stress forcing and geostrophic currents in the 4 to 9∘ N latitudinal band. The simulation forced by JRA55-do significantly reduces the bias in DSL trend over the northern tropical Pacific relative to CORE. In the CORE forcing, the anomalous westerly wind trend in the eastern tropical Pacific causes an underestimated DSL trend across the entire Pacific basin along 10∘ N. The simulation forced by JRA55-do significantly reduces the bias in DSL trend over the northern tropical Pacific relative to CORE. We also identify a bias in the easterly wind trend along 20∘ N in both JRA55-do and CORE, thus motivating future improvement. In JRA55-do, an accurate Rossby wave initiated in the eastern tropical Pacific at seasonal timescale corrects a biased seasonal variability of the northern equatorial countercurrent in the CORE simulation. Both CORE and JRA55-do generate realistic DSL variation during El Niño. We find an asymmetry in the DSL pattern on two sides of the Equator is strongly related to wind stress curl that follows the sea level pressure evolution during El Niño.

CAS-ESM2.0 Model Datasets for the CMIP6 Ocean Model Intercomparison Project Phase 1 (OMIP1)

As a member of the Chinese modeling groups, the coupled ocean-ice component of the Chinese Academy of Sciences’ Earth System Model, version 2.0 (CAS-ESM2.0), is taking part in the Ocean Model Intercomparison Project Phase 1 (OMIP1) experiment of phase 6 of the Coupled Model Intercomparison Project (CMIP6). The simulation was conducted, and monthly outputs have been published on the ESGF (Earth System Grid Federation) data server. In this paper, the experimental dataset is introduced, and the preliminary performances of the ocean model in simulating the global ocean temperature, salinity, sea surface temperature, sea surface salinity, sea surface height, sea ice, and Atlantic Meridional Overturning Circulation (AMOC) are evaluated. The results show that the model is at quasi-equilibrium during the integration of 372 years, and performances of the model are reasonable compared with observations. This dataset is ready to be downloaded and used by the community in related research, e.g., multi-ocean-sea-ice model performance evaluation and interannual variation in oceans driven by prescribed atmospheric forcing.

Impact of horizontal resolution on global ocean–sea ice model simulations based on the experimental protocols of the Ocean Model Intercomparison Project phase 2 (OMIP-2)

Abstract. 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.

The sensitivity of a depth-coordinate model to diapycnal mixing induced by practical implementations of the isopycnal tracer diffusion scheme

The isopycnal tracer diffusion scheme is commonly used in a coarse-resolution global ocean model to parameterize the diffusive effect of mesoscale eddy stirring. This scheme implemented in a depth-coordinate model often accompanies a fail-safe system to prevent numerical instability caused by large isopycnal diffusion around steeply sloped isopycnal surfaces. Its practical implementation artificially induces spurious diapycnal mixing across an isopycnal surface whose slope exceeds a prescribed slope maximum. This slope maximum widely ranges from 1/1000 to 3/10 in a recent global ocean model intercomparison project of Coordinated Ocean-ice Reference Experiments phase II. This study quantitatively investigates the effects of isopycnal tracer diffusion on meridional overturning circulation by calculating water mass transformation rate due to isopycnal diffusion with the use of Walin’s (1982) methodology. We focus on the following three effects of isopycnal diffusion in a depth-coordinate model: densification induced by nonlinearity of the equation of state, spurious diapycnal mixing in the interior ocean induced by practical implementation of the isopycnal diffusion scheme, and diapycnal mixing in the surface diabatic layer. It is revealed that a slope maximum that is too small (1/1000) leads to large spurious diapycnal mixing in the interior ocean. This diapycnal mixing supplies much buoyancy to the bottom water around Antarctica and makes it locally upwell within the Southern Ocean. It results in less northward export of bottom waters from the Southern Ocean into the Atlantic and the Indo-Pacific oceans and weakening of the bottom cell of meridional overturning circulation. This spurious diapycnal mixing is significantly suppressed if we use the slope maximum of more than 1/100 or tapering isopycnal diffusivity around steep isopycnal slopes.This study also shows that these changes in the implementation of the isopycnal diffusion scheme also affect the various aspects of the simulated Southern Ocean. The experiment tapering isopycnal diffusivity suppresses unrealistic open ocean deep convections and polynyas in the Weddell Sea by both diminishing spurious diapycnal mixing in the interior ocean and applying the strong diapycnal mixing in the surface boundary layer associated with the isopycnal diffusion scheme. It results in the smallest biases of the potential temperature distribution, the winter sea ice concentration and mixed layer depth in the Weddell Sea and the ACC transport in our coarse resolution model.

Eddy-resolving Simulation of CAS-LICOM3 for Phase 2 of the Ocean Model Intercomparison Project

A 61-year (1958–2018) global eddy-resolving dataset for phase 2 of the Ocean Model Intercomparison Project has been produced by the version 3 of Chinese Academy of Science, the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics/Institute of Atmospheric Physics (LASG/IAP) Climate system Ocean Model (CAS-LICOM3). The monthly and a part of the surface daily data in this study can be accessed on the Earth System Grid Federation (ESGF) node. Besides the details of the model and experiments, the evolutions and spatial patterns of large-scale and mesoscale features are also presented. The mesoscale features are reproduced well in the high-resolution simulation, as the mesoscale activities can contribute up to 50% of the total SST variability in eddy-rich regions. Also, the large-scale circulations are remarkably improved compared with the low-resolution simulation, such as the climatological annual mean SST (the RMSE is reduced from 0.59°C to 0.47°C, globally) and the evolution of Atlantic Meridional Overturning Circulation. The preliminary evaluation also indicates that there are systematic biases in the salinity, the separation location of the western boundary currents, and the magnitude of eddy kinetic energy. All these biases are worthy of further investigation.