Estimated state of ocean for climate research Version:03b
This version of dataset is an updated version of dataset shown by Osafune et al. (2015) (see them for more details). Our four-dimensional variational (4D-VAR) data synthesis system, which consists of an ocean circulation model, its adjoint and an optimization system, is based on a system developed as a part of the Japan Agency for Marine-Earth Science and Technology (JAMSTEC)-Kyoto University collaborative program (known as "the K7 consortium"). The Ocean General Circulation Model (OGCM) is version 3 of the GFDL Modular Ocean Model (MOM3) , which is equipped with several sophisticated schemes; e.g., Noh scheme for mixed layer physics, the Gent and McWilliams (GM) scheme for isopycnal mixing , and quicker advection scheme. Major physical parameter values are determined through a variational optimization procedure . The horizontal resolution is 1o in both latitude and longitude, with 46 vertical levels spaced from 10 m near the sea surface to 400 m at the bottom. This model is better able to reproduce ocean circulation processes and is expected to form a platform suitable for the use of the 4D-VAR adjoint model.
4D-VAR adjoint data assimilation approach is applied [5,6]. The adjoint codes of the OGCM were obtained using the Tangent linear and Adjoint Model Compiler (TAMC)  and the Transformation of Algorithms in Fortran (TAF) . In the 4D-VAR approach, optimized 4-dimensional datasets are sought by minimizing a cost function [9,10].
The assimilated elements in this study are subsurface temperature and salinity, Sea-surface Temperature (SST), Sea-surface height anomaly (SSHA) and Global Mean Sea Level (GMSL). The subsurface data is from EN4 dataset which was quality controlled using a comprehensive set of objective checks developed at the Hadley Centre of the UK Meteorological Office . This dataset is largely composed of observations from the World Ocean Database 2009  and supplemented by data from the GTSPP (Global Temperature and Salinity Profile Program) and Argo autonomous profiling floats . The SST data is from Reynolds and OISST, and SSHA data is derived from high-precision multi-satellite altimetry products distributed by Copernicus Marine Service. The GMSL data is monthly data until 2013, which is reconstructed based on Church and White (2011) and published by Ocean and Atmosphere unit in the Commonwealth Scientific and Industrial Research Organisation (CSIRO) (http://www.cmar.csiro.au/sealevel/sl_data_cmar.html). All observational data except for GMSL were averaged onto 1o by 1o bins and then compiled as series of 10-day means for the SST and SSHA data and monthly means for the subsurface and GMSL data. The control variables are surface fluxes (for net-heat, fresh water, and momentum) and oceanic initial conditions. The assimilation window is 58 years during 1957-2014.
The synthesis of available observations and a pelagic ecosystem model based on nitrogen cycle produces a dynamically self-consistent dataset. Optimized 4-dimensional datasets are sought by minimizing a cost function on the basis of Green's function approach . The assimilated elements are the climatological monthly mean nitrate from WOA05, monthly mean ocean color data from SeaWiFS, and annual mean chlorophyll-a from WOA98 as detritus.
-  R. C. Pacanowski, S. M. Griffies, The MOM 3 Manual, Geophysical Fluid Dynamics Laboratory/NOAA, Princeton, USA, p.680 (1999).
-  P. R. Gent, J. C. McWilliams, J. Phys. Oceanogr., 20, 150 (1990).
-  Y. Noh, Geophys. Res. Lett., 31, L23305 (2004).
-  D. Menemenlis et al., Mon. Weath. Rev., 133, 1224 (2005).
-  Y. Sasaki, Mon. Weather Rev., 98, 875(1970)
-  C. Wunsch, The Ocean Circulation Inverse Problem, Cambridge Univ. Press, New York, 442 pp (1996).
-  R. Giering, T. Kaminski, Recipes for Adjoint Code Construction, ACM Trans. On Math. Software, 24 (4), 437 (1998).
-  R. Giering, T. Kaminski, Applying TAF to generate efficient derivative code of Fortran 77-95 programs, Proceedings in Applied Mathematics and Mechanics, 2 (1), 54 (2003).
-  J. Marotzke et al., J. Geophys. Res., 104, c12, 29529 (1999).
-  D. Stammer et al., J. Geophys. Res., 107, C9, 3118 (2002).
-  B. Ingleby, M. Huddleston, J. Mar. Sys., 65, 158 (2007).
-  Boyer et al., World Ocean Database 2005, NOAA Atlas NESDIS 60, US Gov. Print. Off., Washington DC, (2006).
-  Gould, Deep Sea Res., II, 52, 529 (2005)
|Variables, abbreviations||Potential temperature [°C],||tmp|
|Horizontal velocity u[m/s] v[m/s],||vel|
|Surface heat flux [cal/m2/s],||shf|
|Surface freshwater flux [m/s],||sff|
|Wind stress τx[N/m2] τy[N/m2],||tau|
|Dissolved inorganic carbon [μmol/kg],||dic|
|Dissolved Oxygen [μmol/L],||oxy|
|Resolution||Horizontal 1°x1°, Vertical 46 levels|
|Period||1957-2014 (Ver. 03b)|
|File format||Monthly data in netcdf format:
“k7oda_[XXX]_[YYYY][MM]00_[VVV].nc” where “[XXX]” is model variable, “[YYYY]” year, “[MM]” month, and “[VVV]”version, respectively.
Reference papers to be cited
Osafune, S., S. Masuda, N. Sugiura, T. Doi (2015) Evaluation of the applicability of the Estimated State of the Global Ocean for Climate Research (ESTOC) data set, Geophys. Res. Lett., 42, 12, 4903–4911.
Doi, T., S. Osafune, N. Sugiura, S. Kouketsu, A. Murata, S. Masuda, and T. Toyoda (2015) Multi-decadal change in the dissolved inorganic carbon in a long-term ocean state estimation, Journal of Advances in Modeling Earth Systems, 7, 4, 1885–1990