MINI REVIEW article
Sec. Ocean Observation
Putting It All Together: Adding Value to the Global Ocean and Climate Observing Systems With Complete Self-Consistent Ocean State and Parameter Estimates
- 1Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX, United States
- 2Institute for Geophysics, The University of Texas at Austin, Austin, TX, United States
- 3Jackson School of Geosciences, The University of Texas at Austin, Austin, TX, United States
- 4Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, United States
- 5Department of Earth, Atmospheric, and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA, United States
- 6Atmospheric and Environmental Research, Lexington, MA, United States
- 7Center für Erdsystemforschung und Nachhaltigkeit, Universität Hamburg, Hamburg, Germany
- 8Department of Earth and Planetary Sciences, Harvard University, Cambridge, MA, United States
- 9Scripps Institution of Oceanography, La Jolla, CA, United States
- 10Department of Physical Oceanography, Woods Hole Oceanographic Institution, Woods Hole, MA, United States
In 1999, the consortium on Estimating the Circulation and Climate of the Ocean (ECCO) set out to synthesize the hydrographic data collected by the World Ocean Circulation Experiment (WOCE) and the satellite sea surface height measurements into a complete and coherent description of the ocean, afforded by an ocean general circulation model. Twenty years later, the versatility of ECCO's estimation framework enables the production of global and regional ocean and sea-ice state estimates, that incorporate not only the initial suite of data and its successors, but nearly all data streams available today. New observations include measurements from Argo floats, marine mammal-based hydrography, satellite retrievals of ocean bottom pressure and sea surface salinity, as well as ice-tethered profiled data in polar regions. The framework also produces improved estimates of uncertain inputs, including initial conditions, surface atmospheric state variables, and mixing parameters. The freely available state estimates and related efforts are property-conserving, allowing closed budget calculations that are a requisite to detect, quantify, and understand the evolution of climate-relevant signals, as mandated by the Coupled Model Intercomparison Project Phase 6 (CMIP6) protocol. The solutions can be reproduced by users through provision of the underlying modeling and assimilation machinery. Regional efforts have spun off that offer increased spatial resolution to better resolve relevant processes. Emerging foci of ECCO are on a global sea level changes, in particular contributions from polar ice sheets, and the increased use of biogeochemical and ecosystem data to constrain global cycles of carbon, nitrogen and oxygen. Challenges in the coming decade include provision of uncertainties, informing observing system design, globally increased resolution, and moving toward a coupled Earth system estimation with consistent momentum, heat and freshwater fluxes between the ocean, atmosphere, cryosphere and land.
The central goal of the ECCO consortium is the production of global ocean state and parameter estimates in support of climate research. ECCO requires dynamical and kinematical consistency of its products, in particular, conservation of mass, heat, and salt throughout the estimation period. Avoiding shortcomings identified in atmospheric reanalysis (e.g., Bengtsson et al., 2004, 2007) and making optimal use of the sparse observational coverage calls for the use of smoothing methods from optimal estimation theory (Wunsch and Heimbach, 2007, 2013; Stammer et al., 2016). The ECCO method exploits information contained in observations both forward and backward in time, while avoiding unphysical perturbations of the time-evolving state that is being constrained. It is the only method that has been found to be practical and that avoids the shortcomings of reanalyses and combines the very diverse ocean data sets that we now have and will continue to collect. The underlying model serves as a “dynamical interpolator” between and beyond the often sparse and heterogeneously sampled observations (in space and time) of various types.
Among ECCO's early accomplishments was the production of the first generation of near-global ocean state estimates, covering the years 1992–1997 (Stammer et al., 2002, 2004; Stammer, 2003). The latest ECCO solution can be used to produce climatologies, based on most data available from the global observing system since the early 1990s, not only for temperature and salinity, but which also provides consistent three-dimensional flow fields and connected dynamical variables (e.g., sea level and bottom pressure), consistent surface forcing fields, and property budgets to explore the underlying dynamics (e.g., Ekman and Sverdrup transports, mixing, and vorticity fluxes) (Fukumori et al., 2018). Self-consistency among the range of state variables is invaluable for depictions of the global ocean, e.g., in terms of its overturning circulation (Cessi, 2019).
2. The Present
2.1. The ECCO Central Production
The ECCO estimation framework in production today has undergone a number of significant improvements and updates. Extending over the period 1992–2015 (an update to 2017 is currently under way), the latest product, ECCO version 4 release 3 (ECCOv4, Forget et al., 2015a; Fukumori et al., 2017), has increased horizontal and vertical resolution and covers the entire globe. The estimation framework has been extended to account for uncertain model parameters that are now routinely part of the inversion (Forget et al., 2015b). The production of the next-generation ECCO version 5 at higher spatial resolution is currently ongoing.
Observational data streams have vastly expanded (Fukumori et al., 2017), and the ways in which these are ingested into the estimation framework have been refined. The space-based backbone consists of daily along-track sea level anomalies from satellite altimetry (Forget and Ponte, 2015) relative to a mean dynamic topography (Andersen et al., 2016), monthly ocean bottom pressure anomalies from GRACE mascon solutions (Watkins et al., 2015), monthly sea surface temperature fields from passive microwave radiometry (Reynolds et al., 2002), monthly sea surface salinity fields from Aquarius (Vinogradova et al., 2014), and daily sea ice concentration fields (Peng et al., 2013; Meier et al., 2017). Major in-situ observing systems used in ECCO include the global array of Argo floats (Roemmich et al., 2009; Riser et al., 2016), ship-based CTD and XBT hydrographic profiles and gridded monthly climatological temperature and salinity fields from the World Ocean Atlas 2009 (WOA09, Antonov et al., 2010; Locarnini et al., 2010), tagged marine mammals (Roquet et al., 2013; Treasure et al., 2017), and ice-tethered profilers (ITPs) in the Arctic (Krishfield et al., 2008). The versatility of the estimation framework enables the inclusion of novel data sets, such as satellite and in-situ inferred electric conductivity as a measure of ocean heat content changes (Trossman and Tyler, 2019). Ocean mixing parameters have been inferred from a subset of Argo, ITP, and hydrographic observations (Cole et al., 2015; Whalen et al., 2015), and are starting to be included in the observational data streams [Trossman et al., in revision].
2.2. Selected Science Applications
Numerous scientific studies have been conducted with various ECCO solutions, leading to new insights into the ocean's role in climate. Partial summaries are in Wunsch et al. (2009), Wunsch and Heimbach (2013), and Fukumori et al. (2018). Here, we highlight two research areas and related studies that have been afforded by the latest ECCO solution. This review only allows for a compressed discussion.
2.2.1. Ocean Heat Content Changes During the Recent Surface Warming Slowdown (SWS) Period
Much attention has been given, both in the scientific literature and in public media to the apparent warming slowdown in global mean surface temperature (GMST) over the first decade of the twenty-first century compared to the 1990s (e.g., Medhaug et al., 2017). The focus on surface temperatures distracted from the fact that a volumetric index such as vertical integrals of heat content changes is a physically more complete climate indicator than (surface) area-based indices. In this context, Nieves et al. (2015) identified issues with several ocean reanalysis products in providing reliable vertical profiles of temperature changes. Figure 1 shows decadal trends in global mean ocean temperature as a function of depth over the periods 1993–2001, 2002–2010, and the difference between the two, from ECCOv4 and two ocean hydrographies. Decadal difference profiles are not available for Argo, which only reached its global coverage in about 2006. The figure is adapted from Nieves et al. (2015), which did not provide any uncertainty estimates for the hydrographies. Compared to the ocean reanalysis trends analyzed by Nieves et al. (2015), which exhibits large deviations from hydrography, ECCOv4 shows a more credible fit to hydrography trends over much of the depth range 0–1,500 m. ECCO's depicted uncertainty (gray shading) represents the formal standard error computed from a least squares linear trend fit to the monthly ECCO values and scaled to account for the effective degrees of freedom (i.e., residual autocorrelation) assuming the residuals of the fit behave as a first-order autoregressive (AR1) process. Note further that the two hydrographics are markedly different in the upper 800 m. ECCOv4 also reproduces the apparent slowdown in surface temperature trends as compared to an optimally interpolated blend of in-situ and satellite SST data. The analysis is set against the larger backdrop of full-depth ocean heat content changes over the last few decades. The latest ECCOv4 estimate produces a global mean heating rate of 0.48 ± 0.16 W m−2, which includes a 0.095 W m−2 geothermal flux (Wunsch, 2018). All uncertainties quoted are likely at lower bounds as they do not account for systematic errors. A full-depth analysis of vertical heat transport by Liang et al. (2017) shows the global mean heat flux imbalances to be small residuals of regionally large anomalies that underly contributions from multiple centers of action, that cooling layers at depths may result from adjustment to surface forcing centuries ago (Gebbie and Huybers, 2019), and the need for accurate budget closure. The use of Argo data since roughly 2006 and satellite altimetric data from 1993 onward in combination with dynamical consistency provides powerful constraints on the ECCO solution over the estimation period.
Figure 1. Decadal trends in global mean potential temperature as function of depth over the periods 1993–2001 (A), 2002–2010 (B), and their difference (C), inferred from two hydrographies, three ocean reanalysis and the ECCOv4 state estimate. Black: ECCOv4 (gray shading indicates formal standard error, see main text), dark blue: WOA (Levitus et al., 2012); red: Ishii (Ishii et al., 2005); purple: GODAS (Huang et al., 2010); green: SODA (Carton and Santorelli, 2008); light blue: ORAS4 (Balmaseda et al., 2012); yellow: SST (Reynolds et al., 2002). Adapted from National Academies of Sciences, Engineering, and Medicine (2016) (their Figures 4, 23).
2.2.2. Origins of North Atlantic Water Mass Volumetric Variability
Quantifying Atlantic water mass variability in terms of its volumetric composition over time provides a powerful diagnostic for determining the relative role of diabatic (locally forced) vs. adiabatic (induced via advection) processes (Forget et al., 2011; Speer and Forget, 2013). An approach is to consider the volume of water contained within temperature classes, following Walin (1982). Such a study has been conducted by Evans et al. (2017) over the period 2004–2011, which includes the marked reduction in the Atlantic Meridional Overturning Circulation (AMOC) inferred at 26 N from the RAPID mooring array (Roberts et al., 2013). Water mass volume anomalies in temperature classes, V(θ, t), between 26 and 45 N (Figure 2, top panels) were derived from a gridded Argo product (Roemmich-Gilson Argo Climatology, RGAC; Roemmich and Gilson, 2009) and ECCOv4. Both products reflect seasonal exchange of volume between the warmer surface waters (θ >18°C) and mode/central waters (θ between 10°C and 18°C), as well as interannual variability in volumetric contributions of subtropical mode water (θ ~ 18°C), among others. Determining water mass transformation rates between temperature classes, dV/dt, proves difficult for RGAC (Figure 2c), conceivably due to aliasing when sampling the mesoscale eddy field, but is feasible for ECCOv4 (Figure 2d). The analysis reveals negative volume anomalies during the winters of 2009/10 and 2010/11. For temperature classes larger than 15°C these anomalies are consistent with diabatic changes inferred from air-sea heat flux diagnostics (Figure 2e). However, for temperatures below 15°C, the adiabatic component as diagnosed from ECCOv4 (Figure 2f) explains the bulk of the volumetric census anomalies. The study provides compelling evidence that wind-driven transport anomalies led to a southward shift in the mean structure of the interior subtropical gyre circulation, weakening northward volume transport at 26 N. Evidence for the role of such advective signals has previously been gathered across an isolated line of latitude from the RAPID mooring array at 26 N (Cunningham et al., 2013).
Figure 2. Top panels: Volume anomaly in temperature classes, V(θ, t), with respect to the time mean in the North Atlantic between 26 and 45 N from (a) RGAC and (b) ECCOv4. Lower panels: Total monthly dV/dt between 26 and 45 N from (c) RGAC and (d) ECCOv4. Also shown in (e) is the diabatic contribution to (d) inferred from monthly diathermal transformation due to air-sea heat fluxes, and (f) the adiabatic transformation in (d) implied by the volume change per temperature class due transport divergence between 26 and 45 N. For details, see Evans et al. (2017). ©American Meteorological Society. Used with permission.
2.3. Regional and Extended-Period Efforts
The tremendous computational cost involved in conducting the nonlinear least-squares optimization problem as well as the occurrence of strong nonlinearities have so far prevented the production of global eddy-resolving decadal state estimates. Instead, regional eddy-permitting estimates of limited duration have spun-off. These include the Southern Ocean State Estimate (SOSE, Mazloff et al., 2010), the Arctic Subpolar gyre sTate Estimate (ASTE, Nguyen et al., 2017), as well as estimates of the California Current System (Verdy et al., 2014), the tropical Pacific (Hoteit et al., 2010; Verdy et al., 2017), and the Gulf of Mexico (Gopalakrishnan et al., 2013). The versatility of the underlying ECCO infrastructure has facilitated these spin-offs. In turn, experience gained in the regional efforts has benefited the global estimation. Other non-ECCO related regional estimation efforts are summarized by Edwards et al. (2015).
An emerging emphasis has been on coupled ocean-sea ice estimation to account for Arctic and Southern Ocean sea ice. Dedicated efforts to develop a dynamic/thermodynamic sea ice model that fits within the estimation framework (Menemenlis et al., 2005; Heimbach et al., 2010; Losch et al., 2010; Fenty and Heimbach, 2013) led to an initial attempt at the global-scale coupled problem (Fenty et al., 2015). A major focus of ASTE is the finding of data used in Arctic research that are not necessarily part of global data repositories and assessing their use in state estimation (Nguyen et al., 2017). Emerging challenges are the use of satellite observations of sea ice (and snow) thickness, as well as remotely sensed drift data to constrain sea ice velocities.
Restricting estimates to the period with available satellite altimetric data limits the applicability of ECCO products for studies of decadal variability. This issue led to a dedicated effort by the German ECCO (GECCO) partners to extend the estimation period back to 1952 (Köhl and Stammer, 2008). Now in its second generation, GECCO2 has extended the period to cover the entire span of the available NCEP reanalysis, but at the cost of sparse observational coverage to constrain the solution prior to the 1990s. On these long timescales, there are issues with convergence of the optimization, which requires splitting the period into a number of smaller assimilation windows (Köhl, 2014). The challenge of quantifying uncertainties in the estimates that go along with changes in the observing system is exacerbated in the long state estimates.
2.4. Estimation Infrastructure, Data Access and Analysis Tools
A key enabling technology of ECCO is the ability to generate an adjoint version of the Massachusetts Institute of Technology general circulation model (MITgcm) for various configurations by means of algorithmic differentiation (Marotzke et al., 1999; Heimbach et al., 2005). Adjoint code generation via the open-source tool OpenAD (Utke et al., 2008) is being pursued. All ECCO state estimates are free-running solutions to the MITgcm and, as such, can be independently reproduced by users interested in performing new experiments (e.g., ocean response to idealized atmospheric wind stress forcing), determining the impact of new data constraints, and generating problem-specific model output (e.g., tracer dispersion). Extending the existing set of model-data misfit constraints is facilitated by the “generic cost” code framework introduced in ECCO v4. Instructions for re-running the model and complete model configurations (including parameters, initial conditions, atmospheric boundary conditions) are provided alongside the solutions (see Table 1).
ECCO products can be accessed via the ECCO webpage (ecco.jpl.nasa.gov). We are currently working to host the standard output fields on NASA's Physical Oceanography Distributed Data Center (PO.DAAC, podaac.jpl.nasa.gov), which will allow users to access the state estimate using several different technologies, including a new secure FTP-like interface (PO.DAAC Drive), Open-source Project for a Network Data Access Protocol (OPeNDAP), Thematic Realtime Environmental Distributed Data Services (THREDDS), and so-called web services enabling access via API protocols. A list of links to data products, model configurations, analysis tools and documentation is summarized in Table 1 in the Data Availability Statement below.
3. The Future
With the increasing accuracy and skill of the ECCO state estimates, new scientific frontiers come into view. Most of these are related to capturing coupled variability, representing secular changes, and closing property budgets across different components of the Earth system (Buizza et al., 2018). In the following, we sketch several coupled problems that appear on the horizon.
3.1. Increased Horizontal Resolution
An increase in horizontal resolution in future ECCO products is targeted to begin resolving the geostrophic eddy field and its impact on the mean circulation. A drawback is the increased degree of nonlinearity of the underlying estimation problem, and the question over which time period the linearization was implied by the adjoint model remains valid. Possible limitations to long assimilation windows have been raised by Lea et al. (2000) and Köhl and Willebrand (2002), among others. A number of computational and practical solutions in the context of estimating statistical properties rather than nonlinear features (“eddy-fitting”) and stabilizing the adjoint to improve controllability at high resolution have been proposed, e.g., by Hoteit et al. (2005), Abarbanel et al. (2010), Wang et al. (2014), and Gebbie and Hsieh (2017). Given the desire within ECCO for maintaining long assimilation windows, i.e. dynamical consistency, these methods are actively being pursued.
3.2. Coupled Ocean-Atmosphere Estimation
A natural extension of ECCO consists in the coupled ocean-atmosphere estimation problem, an avenue pursued by many reanalysis groups today (see Penny et al. [this issue] for a detailed review). Reasons include (i) the ability to close property budgets across the coupled system, (ii) to enable dynamical feedbacks, (iii) to obtain adjusted air-sea fluxes that are consistent with both ocean and atmosphere dynamics, (iv) to infer a coupled state that is balanced with respect to the underlying modeling framework and thus potentially more suitable for initializing extended predictions A major challenge consists in the disparity between oceanic and atmospheric time scales, the time window of validity of the model linearization, which in the atmosphere amounts to synoptic time scales (and in the ocean to resolved eddy turnover time scales), and implications for adjoint model stability for long assimilation windows.
Initial efforts at extending the ECCO capabilities to a fully coupled Earth system model are being conducted using intermediate complexity atmosphere/land models, such as the PlaSim model (Fraedrich et al., 2005; Blessing et al., 2014). This coupled model, called CEN Adjoint Model (CENAM), was put together such that an algorithmic differentiation tool can be used to construct its adjoint for state and parameter estimation purposes. Stammer et al. (2018) present a pilot study for computing adjoint sensitivities of the coupled climate system. To overcome strong nonlinearities, synchronization with observations approaches from dynamical systems theory are being explored to stabilize the adjoint model (Abarbanel et al., 2010; Lyu et al., 2018).
Complementary efforts to understand sensitivities of the ocean to the atmospheric state from a high-end atmospheric model are being conducted in preparation for coupling (Strobach et al., 2018). Other avenues using weak or hybrid coupled assimilation as well as approximate adjoints are also being pursued.
3.3. Coupled Ocean-Ice Sheet Estimation
There is mounting evidence that the increased mass loss from the polar ice sheets, Greenland (WCRP Global Sea Level Budget Group, 2018) and Antarctica (The IMBIE team, 2018), observed over the last two decades is linked to ocean circulation changes that have brought about warmer waters to the grounding zones of marine-terminating glaciers and ice shelves. Concerns over the implications of rising sea levels call for the joint treatment of the coupled ocean-ice sheet system. Substantial progress is being made, both with asynchronous coupling between the MITgcm and the Ice Sheet System Model (ISSM, Seroussi et al., 2017) as well as with synchronous, property-conserving coupling between the MITgcm's ocean and ice stream/shelf model (Goldberg et al., 2018; Jordan et al., 2018). The availability of adjoint models of all of these components, along with at least annually resolved satellite observations at Antarctica's marine margins, offer the prospect of developing a tightly coupled, skillful estimation system.
3.4. Coupled Ocean-Biogeochemistry and Ecology Estimation
The advent of profiling floats equipped with biogeochemical (BGC) sensors presents a revolution in data density for constraining BGC and ecosystem models. The software exists to assimilate these measurements along with remote sensing of ocean color into models (e.g., Gregg et al., 2009; Song et al., 2016; Verdy and Mazloff, 2017). BGC ocean property observations constrain many aspects of the Earth system, such that coupling not only informs the carbon system and ocean health, but also improves many other components of the Earth system models. Another thrust is the development of the ECCO-Darwin project, which combines physical and biological observations with the coupled framework of the eddy-permitting ECCO and Darwin ecosystem models (Follows and Dutkiewicz, 2011), but with significant remaining obstacles (Dutkiewicz et al., 2018).
3.5. Uncertainty Quantification (UQ) and Optimal Observing Network Design
Although formally an integral part of state and parameter estimation, deriving formal uncertainties accompanying the optimal estimates adds another level of computational complexity (National Academies of Sciences, Engineering, and Medicine, 2012). This has so far prevented most ocean reanalysis (or estimation) projects from dealing comprehensively with UQ. In the context of derivative-based estimation, identification of key metrics (or quantities) of interest enables the development of a formal chain that propagates the uncertainties from observations and uncertain parameters (priors) through the inference (i.e., posterior uncertainties at the optimal estimate) to the derived metrics of interest (Kalmikov and Heimbach, 2014, 2018). This Hessian-based framework lends itself to conducting optimal observing system design studies (see Fujii et al., under review) that provide valuable information on the optimal placement of available observational assets to maximize their utility in constraining key oceanographic quantities of interest (Köhl and Stammer, 2004).
3.6. Synergistic Use of Products and Model
While the state estimates are the central product of ECCO estimation, the virtue of their physical consistency is best realized by their analysis in conjunction with the underlying ocean general circulation model. The state estimates provide descriptions of the ocean, whereas the model affords its explanation; e.g., why is the ocean state what it is and why does it change as it does? As experience is gained, application of state estimation is expanding from drawing inferences from sampling the estimates akin to observations to quantitatively analyzing processes by utilizing the complete physics embodied in the model. Examples of such include analyses of property budgets that are closed without unresolved components (e.g., Buckley et al., 2015; Piecuch et al., 2017; Ponte and Piecuch, 2018), tracing origins and fate of ocean water masses (e.g., Fukumori et al., 2004; Gao et al., 2011; Qu et al., 2013) and quantifying causal mechanisms controlling the ocean (e.g., Fukumori et al., 2015; Pillar et al., 2016, 2018; Jones et al., 2018; Smith and Heimbach, 2019). The model's adjoint offers a unique tool in such efforts by providing an efficient means to evaluate physical dependencies among different quantities of interest. While the fidelity of state estimation will continue to evolve, existing systems provide a means to understanding and explaining what they do already resolve of the ocean. The full exploitation of state estimation requires a holistic approach and is ripe for innovation.
4. Concluding Remarks
The past two decades have seen substantial progress in the development and production of rigorous global ocean state and parameter estimates in support of climate research. That development has relied in part on the availability of continuous climate-quality records of quasi-global coverage, beginning with satellite altimetry (since 1992), satellite gravimetry (since 2003), and hydrographic profiles from the Argo float program (globally since ca. 2006). Sustaining such observing systems over long periods of time to build a climate record is a key imperative of ocean and climate monitoring (National Academies of Sciences, Engineering, and Medicine, 2017). The underlying computational estimation approaches used in the model-data synthesis serve several purposes: (i) they extract optimal information from the sparse and heterogeneous observational streams that constitute the Global Ocean Observing System (GOOS), (ii) they provide a quantitative framework for hypothesis testing and model parameter calibration, and (iii) they enable a quantitative understanding of the underlying dynamical and physical processes that have been learned jointly from observations and models. Much of what these approaches offer, for rigorous climate model calibration and initialization, remains under-explored. Realizing their full potential faces substantial practical hurdles but is indispensable for tackling important issues in ocean climate science. Increasing horizontal resolution and moving toward a comprehensive coupled Earth system estimation framework are major thrusts for the decade ahead.
ECCO strives to make all of its products available online to the scientific community. This includes the state estimates, ancillary fields to perform accurate budget calculations, the complete model configuration to reproduce the state estimate, analysis tools, as well as documentation. Table 1 provides a comprehensive list of links to these resources.
All authors listed have made a substantial, direct and intellectual contribution to the work, and approved it for publication.
Major support for ECCO is provided by NASA's Physical Oceanography program via a contract to JPL/Caltech, with additional support through NASA's Modeling, Analysis and Prediction program, the Cryosphere Science program, and the Computational Modeling and Cyberinfrastructure program. Supplemental funding was obtained throughout the years via standard grants to individual team members from NSF, NOAA, and ONR.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We are grateful to NASA's Physical Oceanography Program for continued support of ECCO throughout the years, and to NASA's High-End Computing (HEC) Program for providing outstanding supercomputing resources and user support at the NASA Advanced Supercomputing (NAS) Division at Ames Research Center. We gratefully acknowledge all groups and programs who provided observational data sets used to constrain the ECCO estimate. We are grateful for the help of various support staff throughout the years, in particular Charmaine King and Diana Spiegel at MIT, and Sue Rodriguez at UT Austin.
Andersen, O., Knudsen, P., and Stenseng, L. (2016). “The DTU13 MSS (mean sea surface) and MDT (mean dynamic topography) from 20 years of satellite altimetry,” in IGFS 2014 (Cham: Springer International Publishing), 111–121.
Antonov, J. I., Seidov, D., Boyer, T. P., Locarnini, R. A., Mishonov, A. V., Garcia, H. E., et al. (2010). World Ocean Atlas 2009, Vol. 2: Salinity. Washington, DC: NOAA Atlas NESDIS 69, U.S. Government Printing Office.
Bengtsson, L., Haines, K., Hodges, K. I., Arkin, P., Berrisford, P., Bougeault, P., et al. (2007). The need for a dynamical climate reanalysis. Bull. Am. Meteorol. Soc. 88, 495–501. doi: 10.1175/BAMS-88-4-495
Blessing, S., Kaminski, T., Lunkeit, F., Matei, I., Giering, R., Köhl, A., et al. (2014). Testing variational estimation of process parameters and initial conditions of an earth system model. Tellus A 66, 769. doi: 10.3402/tellusa.v66.22606
Buckley, M. W., Ponte, R. M., Forget, G., and Heimbach, P. (2015). Determining the Origins of Advective Heat Transport Convergence Variability in the North Atlantic. J. Clim. 28, 3943–3956. doi: 10.1175/JCLI-D-14-00579.1
Buizza, R., Brönnimann, S., Haimberger, L., Laloyaux, P., Martin, M. J., Fuentes, M., et al. (2018). The EU-FP7 ERA-CLIM2 project contribution to advancing science and production of earth system climate reanalyses. Bull. Am. Meteorol. Soc. 99, 1003–1014. doi: 10.1175/BAMS-D-17-0199.1
Cole, S. T., Wortham, C., Kunze, E., and Owens, W. B. (2015). Eddy stirring and horizontal diffusivity from Argo float observations: geographic and depth variability. Geophys. Res. Lett. 42, 3989–3997. doi: 10.1002/2015GL063827
Cunningham, S. A., Roberts, C. D., Frajka-Williams, E., Johns, W. E., Hobbs, W., Palmer, M. D., et al. (2013). Atlantic meridional overturning circulation slowdown cooled the subtropical ocean. Geophys. Res. Lett. 40, 6202–6207. doi: 10.1002/2013GL058464
Evans, D. G., Toole, J., Forget, G., Zika, J. D., Naveira Garabato, A. C., Nurser, A. J. G., et al. (2017). Recent wind-driven variability in atlantic water mass distribution and meridional overturning circulation. J. Phys. Oceanogr. 47, 633–647. doi: 10.1175/JPO-D-16-0089.1
Forget, G., Campin, J. M., Heimbach, P., Hill, C. N., Ponte, R. M., and Wunsch, C. (2015a). ECCO version 4: an integrated framework for non-linear inverse modeling and global ocean state estimation. Geosci. Model Dev. 8, 3071–3104. doi: 10.5194/gmd-8-3071-2015
Forget, G., Ferreira, D., and Liang, X. (2015b). On the observability of turbulent transport rates by Argo: supporting evidence from an inversion experiment. Ocean Sci. 11, 839–853. doi: 10.5194/os-11-839-2015
Fukumori, I., Lee, T., Cheng, B., and Menemenlis, D. (2004). The origin, pathway, and destination of Niño-3 water estimated by a simulated passive tracer and its adjoint. J. Phys. Oceanogr. 34, 582–604. doi: 10.1175/2515.1
Fukumori, I., Wang, O., Llovel, W., Fenty, I., and Forget, G. (2015). A near-uniform fluctuation of ocean bottom pressure and sea level across the deep ocean basins of the Arctic Ocean and the Nordic Seas. Progr. Oceanogr. 134, 152–172. doi: 10.1016/j.pocean.2015.01.013
Gao, S., Qu, T., and Fukumori, I. (2011). Effects of mixing on the subduction of South Pacific waters identified by a simulated passive tracer and its adjoint. Dyn. Atmos. Oceans 51, 45–54. doi: 10.1016/j.dynatmoce.2010.10.002
Goldberg, D. N., Snow, K., Holland, P., Jordan, J. R., Campin, J. M., Heimbach, P., et al. (2018). Representing grounding line migration in synchronous coupling between a marine ice sheet model and a z -coordinate ocean model. Ocean Model. 125, 45–60. doi: 10.1016/j.ocemod.2018.03.005
Gopalakrishnan, G., Cornuelle, B. D., Hoteit, I., Rudnick, D. L., and Owens, W. B. (2013). State estimates and forecasts of the loop current in the Gulf of Mexico using the MITgcm and its adjoint. J. Geophys. Res. Oceans 118, 3292–3314. doi: 10.1002/jgrc.20239
Gregg, W. W., Friedrichs, M. A. M., Robinson, A. R., Rose, K. A., Schlitzer, R., Thompson, K. R., et al. (2009). Skill assessment in ocean biological data assimilation. J. Mar. Syst. 76, 16–33. doi: 10.1016/j.jmarsys.2008.05.006
Heimbach, P., Hill, C., and Giering, R. (2005). An efficient exact adjoint of the parallel MIT General Circulation Model, generated via automatic differentiation. Future Generat. Comput. Syst. 21, 1356–1371. doi: 10.1016/j.future.2004.11.010
Heimbach, P., Menemenlis, D., Losch, M., Campin, J.-M., and Hill, C. (2010). On the formulation of sea-ice models. Part 2: lessons from multi-year adjoint sea-ice export sensitivities through the Canadian Arctic Archipelago. Ocean Model. 33, 145–158. doi: 10.1016/j.ocemod.2010.02.002
Hoteit, I., Cornuelle, B., and Heimbach, P. (2010). An eddy-permitting, dynamically consistent adjoint-based assimilation system for the tropical Pacific: Hindcast experiments in 2000. J. Geophys. Res. 115, C03001–C03023. doi: 10.1029/2009JC005437
Hoteit, I., Cornuelle, B., Köhl, A., and Stammer, D. (2005). Treating strong adjoint sensitivities in tropical eddy-permitting variational data assimilation. Q. J. R. Meteorol. Soc. 131, 3659–3682. doi: 10.1256/qj.05.97
Huang, B., Xue, Y., Zhang, D., Kumar, A., and McPhaden, M. J. (2010). The NCEP GODAS ocean analysis of the tropical pacific mixed layer heat budget on seasonal to interannual time scales. J. Clim. 23, 4901–4925. doi: 10.1175/2010JCLI3373.1
Ishii, M., Shouji, A., Sugimoto, S., and Matsumoto, T. (2005). Objective analyses of sea-surface temperature and marine meteorological variables for the 20th century using ICOADS and the Kobe Collection. Int. J. Climatol. 25, 865–879. doi: 10.1002/joc.1169
Jones, D. C., Forget, G., Sinha, B., Josey, S. A., Boland, E. J. D., Meijers, A. J. S., et al. (2018). Local and remote influences on the heat content of the labrador sea: an adjoint sensitivity study. J. Geophys. Res. Oceans 123, 2646–2667. doi: 10.1002/2018JC013774
Jordan, J. R., Holland, P. R., Goldberg, D., Snow, K., Arthern, R., Campin, J.-M., et al. (2018). Ocean-forced ice-shelf thinning in a synchronously coupled ice-ocean model. J. Geophys. Res. Oceans 125, 864–882. doi: 10.1002/2017JC013251
Kalmikov, A. G., and Heimbach, P. (2018). On barotropic mechanisms of uncertainty propagation in estimation of drake passage transport. arXiv.org. arXiv:1804.06033v2. Available online at: https://arxiv.org/abs/1804.06033v2
Köhl, A., and Willebrand, J. (2002). An adjoint method for the assimilation of statistical characteristics into eddy-resolving ocean models. Tellus A 54, 406–425. doi: 10.1034/j.1600-0870.2002.01294.x
Krishfield, R., Toole, J., Proshutinsky, A., and Timmermans, M. L. (2008). Automated ice-tethered profilers for seawater observations under pack ice in all seasons. J. Atmos. Oceanic Technol. 25, 2091–2105. doi: 10.1175/2008JTECHO587.1
Levitus, S., Antonov, J. I., Boyer, T. P., Baranova, O. K., Garcia, H. E., Locarnini, R. A., et al. (2012). World ocean heat content and thermosteric sea level change (0–2000 m), 1955–2010. Geophys. Res. Lett. 39:L10603. doi: 10.1029/2012GL051106
Liang, X., Piecuch, C. G., Ponte, R. M., Forget, G., Wunsch, C., and Heimbach, P. (2017). Change of the global ocean vertical heat transport over 1993–2010. J. Clim. 30, 5319–5327. doi: 10.1002/2015GL064156
Locarnini, R. A., Mishonov, A. V., Antonov, J. I., Boyer, T. P., Garcia, H. E., Baranova, O. K., et al. (2010). World Ocean Atlas 2009, Vol. 1, Temperature. Washington, DC: NOAA Atlas NESDIS 68, U.S. Government Printing Office.
Losch, M., Menemenlis, D., Campin, J.-M., Heimbach, P., and Hill, C. (2010). On the formulation of sea-ice models. Part 1: effects of different solver implementations and parameterizations. Ocean Model. 33, 129–144. doi: 10.1016/j.ocemod.2009.12.008
Marotzke, J., Giering, R., Zhang, K. Q., Stammer, D., Hill, C., and Lee, T. (1999). Construction of the adjoint MIT ocean general circulation model and application to Atlantic heat transport sensitivity. J. Geophys. Res. 104, 529–548.
National Academies of Sciences Engineering, and Medicine. (2012). “Assessing the Reliability of Complex Models,” in Mathematical and Statistical Foundations of Verification, Validation, and Uncertainty Quantification (Washington, DC: National Academies Press).
National Academies of Sciences Engineering, and Medicine. (2016). “Frontiers in Decadal Climate Variability: Proceedings of a Workshop,” in Proceedings of a Workshop (Washington, DC: National Academies Press).
Nguyen, A., Ocaña, V., Garg, V., Heimbach, P., Toole, J., Krishfield, R., et al. (2017). On the benefit of current and future ALPS data for improving arctic coupled ocean-sea ice state estimation. Oceanography 30, 69–73. doi: 10.5670/oceanog.2017.223
Peng, G., Meier, W. N., Scott, D. J., and Savoie, M. H. (2013). A long-term and reproducible passive microwave sea ice concentration data record for climate studies and monitoring. Earth Syst. Sci. Data 5, 311–318. doi: 10.5194/essd-5-311-2013
Piecuch, C. G., Ponte, R. M., Little, C. M., Buckley, M. W., and Fukumori, I. (2017). Mechanisms underlying recent decadal changes in subpolar North Atlantic Ocean heat content. J. Geophys. Res. Oceans 122, 7181–7197. doi: 10.1002/2017JC012845
Pillar, H. R., Johnson, H. L., Marshall, D. P., Heimbach, P., and Takao, S. (2018). Impacts of atmospheric reanalysis uncertainty on Atlantic overturning estimates at 25°N. J. Clim, 31, 8719–8744. doi: 10.1175/JCLI-D-18-0241.1
Qu, T., Gao, S., and Fukumori, I. (2013). Formation of salinity maximum water and its contribution to the overturning circulation in the North Atlantic as revealed by a global general circulation model. J. Geophys. Res. Oceans 118, 1982–1994. doi: 10.1002/jgrc.20152
Reynolds, R. W., Rayner, N. A., Smith, T. M., Stokes, D. C., and Wang, W. (2002). An improved in situ and satellite SST analysis for climate. J. Clim. 15, 1609–1625. doi: 10.1175/1520-0442(2002)015<1609:AIISAS>2.0.CO;2
Riser, S. C., Freeland, H. J., Roemmich, D., Wijffels, S., Troisi, A., Belbéoch, M., et al. (2016). Fifteen years of ocean observations with the global Argo array. Nat. Clim. Change 6, 145–153. doi: 10.1038/nclimate2872
Roberts, C. D., Waters, J., Peterson, K. A., Palmer, M. D., McCarthy, G. D., Frajka-Williams, E., et al. (2013). Atmosphere drives recent interannual variability of the Atlantic meridional overturning circulation at 26.5N. Geophys. Res. Lett. 40, 5164–5170. doi: 10.1002/grl.50930
Roemmich, D., and Gilson, J. (2009). The 2004–2008 mean and annual cycle of temperature, salinity, and steric height in the global ocean from the Argo Program. Progr. Oceanogr. 82, 81–100. doi: 10.1016/j.pocean.2009.03.004
Roemmich, D., Johnson, G., Riser, S., Davis, R., Gilson, J., Owens, W. B., et al. (2009). The argo program: observing the global oceans with profiling floats. Oceanography 22, 34–43. doi: 10.5670/oceanog.2009.36
Roquet, F., Wunsch, C., Forget, G., Heimbach, P., Guinet, C., Reverdin, G., et al. (2013). Estimates of the Southern Ocean general circulation improved by animal-borne instruments. Geophys. Res. Lett. 40, 6176–6180. doi: 10.1002/2013GL058304
Seroussi, H., Nakayama, Y., Larour, E., Menemenlis, D., Morlighem, M., Rignot, E., et al. (2017). Continued retreat of Thwaites Glacier, West Antarctica, controlled by bed topography and ocean circulation. Geophys. Res. Lett. 44, 6191–6199. doi: 10.1002/2017GL072910
Song, H., Edwards, C. A., Moore, A. M., and Fiechter, J. (2016). Data assimilation in a coupled physical-biogeochemical model of the California current system using an incremental lognormal 4-dimensional variational approach: Part 3—Assimilation in a realistic context using satellite and in situ observations. Ocean Model. 106, 159–172. doi: 10.1016/j.ocemod.2016.06.005
Speer, K., and Forget, G. (2013). “Global distribution and formation of mode waters,” in Ocean Circulation and Climate: A 21st Century Perspective, eds G. Siedler, S. M. Griffies, J. Gould, and J. A. Church (Amsterdam: Elsevier Ltd).
Stammer, D. (2003). Volume, heat, and freshwater transports of the global ocean circulation 1993–2000, estimated from a general circulation model constrained by World Ocean Circulation Experiment (WOCE) data. J. Geophys. Res. 108, 3007–7–23. doi: 10.1029/2001JC001115
Stammer, D., Balmaseda, M., Heimbach, P., Köhl, A., and Weaver, A. (2016). Ocean data assimilation in support of climate applications: status and perspectives. Ann. Rev. Mar. Sci. 8, 491–518. doi: 10.1146/annurev-marine-122414-034113
Stammer, D., Köhl, A., Vlasenko, A., Matei, I., Lunkeit, F., and Schubert, S. (2018). A pilot climate sensitivity study using the CEN coupled adjoint model (CESAM). J. Clim. 31, 2031–2056. doi: 10.1175/JCLI-D-17-0183.1
Stammer, D., Ueyoshi, K., Köhl, A., Large, W. G., Josey, S. A., and Wunsch, C. (2004). Estimating air-sea fluxes of heat, freshwater, and momentum through global ocean data assimilation. J. Geophys. Res. 109:C05023. doi: 10.1029/2003JC002082
Stammer, D., Wunsch, C., Giering, R., Eckert, C., Heimbach, P., Marotzke, J., et al. (2002). Global ocean circulation during 1992–1997, estimated from ocean observations and a general circulation model. J. Geophys. Res. 107, 3118–1–27. doi: 10.1029/2001JC000888
Strobach, E., Molod, A., Forget, G., Campin, J.-M., Hill, C., Menemenlis, D., et al. (2018). Consequences of different air-sea feedbacks on ocean using MITgcm and MERRA-2 forcing: implications for coupled data assimilation systems. Ocean Model. 132, 91–111. doi: 10.1016/j.ocemod.2018.10.006
Treasure, A., Roquet, F., Ansorge, I., Bester, M., Boehme, L., Bornemann, H., et al. (2017). Marine mammals exploring the oceans pole to pole: a review of the MEOP consortium. Oceanography 30, 132–138. doi: 10.5670/oceanog.2017.234
Utke, J., Naumann, U., Fagan, M., Tallent, N., Strout, M., Heimbach, P., et al. (2008). OpenAD/F: a modular open-source tool for automatic differentiation of fortran codes. ACM Trans. Math. Softw. 34, 1–36. doi: 10.1145/1377596.1377598
Verdy, A., Mazloff, M. R., Cornuelle, B. D., and Kim, S. Y. (2014). Wind-driven sea level variability on the California coast: an adjoint sensitivity analysis. J. Phys. Oceanogr. 44, 297–318. doi: 10.1175/JPO-D-13-018.1
Vinogradova, N. T., Ponte, R. M., Fukumori, I., and Wang, O. (2014). Estimating satellite salinity errors for assimilation of Aquarius and SMOS data into climate models. J. Geophys. Res. Oceans 119, 4732–4744. doi: 10.1002/2014JC009906
Watkins, M. M., Wiese, D. N., Yuan, D.-N., Boening, C., and Landerer, F. W. (2015). Improved methods for observing Earth's time variable mass distribution with GRACE using spherical cap mascons. J. Geophys. Res. Solid Earth 120, 2648–2671. doi: 10.1002/2014JB011547
Whalen, C. B., MacKinnon, J. A., Talley, L. D., and Waterhouse, A. F. (2015). Estimating the mean diapycnal mixing using a finescale strain parameterization. J. Phys. Oceanogr. 45, 1174–1188. doi: 10.1175/JPO-D-14-0167.1
Wunsch, C., and Heimbach, P. (2013). “Dynamically and kinematically consistent global ocean circulation and ice state estimates,” in Ocean Circulation and Climate: A 21st Century Perspective, eds G. Siedler, S. M. Griffies, J. Gould, J. A. Church (Amsterdam: Elsevier Ltd.), 553–579.
Keywords: ECCO, global ocean inverse modeling, optimal state and parameter estimation, adjoint method, ocean observations, coupled Earth system data assimilation, ocean reanalysis, global ocean circulation
Citation: Heimbach P, Fukumori I, Hill CN, Ponte RM, Stammer D, Wunsch C, Campin J-M, Cornuelle B, Fenty I, Forget G, Köhl A, Mazloff M, Menemenlis D, Nguyen AT, Piecuch C, Trossman D, Verdy A, Wang O and Zhang H (2019) Putting It All Together: Adding Value to the Global Ocean and Climate Observing Systems With Complete Self-Consistent Ocean State and Parameter Estimates. Front. Mar. Sci. 6:55. doi: 10.3389/fmars.2019.00055
Received: 14 November 2018; Accepted: 01 February 2019;
Published: 04 March 2019.
Edited by:John Siddorn, Met Office, United Kingdom
Reviewed by:Keith Haines, University of Reading, United Kingdom
Marcin Chrust, European Centre for Medium-Range Weather Forecasts, United Kingdom
Copyright © 2019 Heimbach, Fukumori, Hill, Ponte, Stammer, Wunsch, Campin, Cornuelle, Fenty, Forget, Köhl, Mazloff, Menemenlis, Nguyen, Piecuch, Trossman, Verdy, Wang and Zhang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Patrick Heimbach, firstname.lastname@example.org