Description: We present a new surface-atmospheric dataset for driving ocean–sea-ice models based on Japanese 55-year atmospheric reanalysis (JRA-55), referred to here as JRA55-do. The JRA55-do dataset aims to replace the CORE interannual forcing version 2 (hereafter called the CORE dataset), which is currently used in the framework of the Coordinated Ocean-ice Reference Experiments (COREs) and the Ocean Model Intercomparison Project (OMIP). A major improvement in JRA55-do is the refined horizontal grid spacing (∼ 55 km) and temporal interval (3 hr). The data production method for JRA55-do essentially follows that of the CORE dataset, whereby the surface fields from an atmospheric reanalysis are adjusted relative to reference datasets. To improve the adjustment method, we use high-quality products derived from satellites and from several other atmospheric reanalysis projects, as well as feedback on the CORE dataset from the ocean modelling community. Notably, the surface air temperature and specific humidity are adjusted using multi-reanalysis ensemble means. In JRA55-do, the downwelling radiative fluxes and precipitation, which are affected by an ambiguous cloud parameterisation employed in the atmospheric model used for the reanalysis, are based on the reanalysis products. This approach represents a notable change from the CORE dataset, which imported independent observational products. Consequently, the JRA55-do dataset is more self-contained than the CORE dataset, and thus can be continually updated in near real-time. The JRA55-do dataset extends from 1958 to the present, with updates expected at least annually. This paper details the adjustments to the original JRA-55 fields, the scientific rationale for these adjustments, and the evaluation of JRA55-do. The adjustments successfully corrected the biases in the original JRA-55 fields. The globally averaged features are similar between the JRA55-do and CORE datasets, implying that JRA55-do can suitably replace the CORE dataset for use in driving global ocean–sea-ice models.

Description: Large amounts of the greenhouse gas methane are released from the seabed to the water column1, where it may be consumed by aerobic methanotrophic bacteria2. The size and activity of methanotrophic communities, which determine the amount of methane consumed in the water column, are thought to be mainly controlled by nutrient and redox dynamics3–7. Here, we report repeated measurements of methanotrophic activity and community size at methane seeps west of Svalbard, and relate them to physical water mass properties and modelled ocean currents. We show that cold bottom water, which contained a large number of aerobic methanotrophs, was displaced by warmer water with a considerably smaller methanotrophic community within days. Ocean current simulations using a global ocean/sea-ice model suggest that this water mass exchange is consistent with short-term variations in the meandering West Spitsbergen Current. We conclude that the shift from an offshore to a nearshore position of the current can rapidly and severely reduce methanotrophic activity in the water column. Strong fluctuating currents are common at many methane seep systems globally, and we suggest that they affect methane oxidation in the water column at other sites, too.

Description: Highlights: • Inter-annual to decadal variability in AMOC from CORE-II simulations is presented. • AMOC variability shows three stages, with maximum transports in mid- to late-1990s. • North Atlantic temporal variability features are in good agreement among simulations. • Such agreements suggest variability is dictated by the atmospheric data sets. • Simulations differ in spatial structures of variability due to ocean dynamics. Simulated inter-annual to decadal variability and trends in the North Atlantic for the 1958–2007 period from twenty global ocean – sea-ice coupled models are presented. These simulations are performed as contributions to the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II). The study is Part II of our companion paper (Danabasoglu et al., 2014) which documented the mean states in the North Atlantic from the same models. A major focus of the present study is the representation of Atlantic meridional overturning circulation (AMOC) variability in the participating models. Relationships between AMOC variability and those of some other related variables, such as subpolar mixed layer depths, the North Atlantic Oscillation (NAO), and the Labrador Sea upper-ocean hydrographic properties, are also investigated. In general, AMOC variability shows three distinct stages. During the first stage that lasts until the mid- to late-1970s, AMOC is relatively steady, remaining lower than its long-term (1958–2007) mean. Thereafter, AMOC intensifies with maximum transports achieved in the mid- to late-1990s. This enhancement is then followed by a weakening trend until the end of our integration period. This sequence of low frequency AMOC variability is consistent with previous studies. Regarding strengthening of AMOC between about the mid-1970s and the mid-1990s, our results support a previously identified variability mechanism where AMOC intensification is connected to increased deep water formation in the subpolar North Atlantic, driven by NAO-related surface fluxes. The simulations tend to show general agreement in their temporal representations of, for example, AMOC, sea surface temperature (SST), and subpolar mixed layer depth variabilities. In particular, the observed variability of the North Atlantic SSTs is captured well by all models. These findings indicate that simulated variability and trends are primarily dictated by the atmospheric datasets which include the influence of ocean dynamics from nature superimposed onto anthropogenic effects. Despite these general agreements, there are many differences among the model solutions, particularly in the spatial structures of variability patterns. For example, the location of the maximum AMOC variability differs among the models between Northern and Southern Hemispheres.

Description: A global ocean model with 1/12∘ horizontal resolution is used to assess the seasonal cycle of surface Eddy Kinetic Energy (EKE). The model reproduces the salient features of the observed mean surface EKE, including amplitude and phase of its seasonal cycle in most parts of the ocean. In all subtropical gyres of the Pacific and Atlantic, EKE peaks in summer down to a depth of ∼350 m, below which the seasonal cycle is weak. Investigation of the possible driving mechanisms reveals the seasonal changes in the thermal interactions with the atmosphere to be the most likely cause of the summer maximum of EKE. The development of the seasonal thermocline in spring and summer is accompanied by stronger mesoscale variations in the horizontal temperature gradients near the surface which corresponds, by thermal wind balance, to an intensification of mesoscale velocity anomalies towards the surface.

Description: Highlights: • Mean circulation patterns are assessed and Kuroshio transport is underestimated. • Water mass distribution is compared and analyzed within COREII models. • Main biases of deep MLDs result from the inaccurate Kuroshio separation. • Reasonable modeled tropical dynamics but a discrepancy from the surface wind. Abstract: We evaluate the mean circulation patterns, water mass distributions, and tropical dynamics of the North and Equatorial Pacific Ocean based on a suite of global ocean-sea ice simulations driven by the CORE-II atmospheric forcing from 1963-2007. The first three moments (mean, standard deviation and skewness) of sea surface height and surface temperature variability are assessed against observations. Large discrepancies are found in the variance and skewness of sea surface height and in the skewness of sea surface temperature. Comparing with the observation, most models underestimate the Kuroshio transport in the Asian Marginal seas due to the missing influence of the unresolved western boundary current and meso-scale eddies. In terms of the Mixed Layer Depths (MLDs) in the North Pacific, the two observed maxima associated with Subtropical Mode Water and Central Mode Water formation coalesce into a large pool of deep MLDs in all participating models, but another local maximum associated with the formation of Eastern Subtropical Mode Water can be found in all models with different magnitudes. The main model bias of deep MLDs results from excessive Subtropical Mode Water formation due to inaccurate representation of the Kuroshio separation and of the associated excessively warm and salty Kuroshio water. Further water mass analysis shows that the North Pacific Intermediate Water can penetrate southward in most models, but its distribution greatly varies among models depending not only on grid resolution and vertical coordinate but also on the model dynamics. All simulations show overall similar large scale tropical current system, but with differences in the structures of the Equatorial Undercurrent. We also confirm the key role of the meridional gradient of the wind stress curl in driving the equatorial transport, leading to a generally weak North Equatorial Counter Current in all models due to inaccurate CORE-II equatorial wind fields. Most models show a larger interior transport of Pacific subtropical cells than the observation due to the overestimated transport in the Northern Hemisphere likely resulting from the deep pycnocline

Description: Version 1.0.2 (2016)

Description: The Agulhas Current, the western boundary current of the South Indian Ocean, has been shown to play an important role in the connectivity between the Indian and Atlantic oceans. The greater Agulhas Current system is highly dominated by mesoscale dynamics. To investigate their influence on the regional and global circulations, a family of high-resolution ocean general circulation model configurations based on the NEMO code has been developed. Horizontal resolution refinement is achieved by embedding “nests” covering the South Atlantic and the western Indian oceans at 1/10∘ (INALT10) and 1/20∘ (INALT20) within global hosts with coarser resolutions. Nests and hosts are connected through two-way interaction, allowing the nests not only to receive boundary conditions from their respective host but also to feed back the impact of regional dynamics onto the global ocean. A double-nested configuration at 1/60∘ resolution (INALT60) has been developed to gain insights into submesoscale processes within the Agulhas Current system. Large-scale measures such as the Drake Passage transport and the strength of the Atlantic meridional overturning circulation are rather robust among the different configurations, indicating the important role of the hosts in providing a consistent embedment of the regionally refined grids into the global circulation. The dynamics of the Agulhas Current system strongly depend on the representation of mesoscale processes. Both the southward-flowing Agulhas Current and the northward-flowing Agulhas Undercurrent increase in strength with increasing resolution towards more realistic values, which suggests the importance of improving mesoscale dynamics as well as bathymetric slopes along this narrow western boundary current regime. The exploration of numerical choices such as lateral boundary conditions and details of the implementation of surface wind stress forcing demonstrates the range of solutions within any given configuration.

Description: Sea surface height (SSH) variability in the tropical Pacific is determined by surface fluxes of which momentum flux is the most important one. The new atmospheric forcing data set JRA-55-do (Tsujino, H. et al. (submitted)) offers the chance of an improved representation of tropical SSH variability. We present a comparison of SSH variability of the tropical Pacific in two global ORCA025 (Nemo3.6) configurations, forced with COREv2 and JRA-55-do atmospheric datasets. While the amplitude of SSH variability in COREv2-runs is in close agreement with altimeter observations, a reduced interannual variability of wind stress in JRA-55-do leads to weaker SSH-variability in the tropical Pacific. A lagged correlation analyses of SSH with climate indices such as ENSO or Southern Oscillation shows a higher agreement between JRA-55-do and altimeter observations than between COREv2 and observations, suggesting an improved representation of the processes that determine SSH variability