Drake, Henri F., Shanice Bailey, Raphael Dussin, Stephen M Griffies, John P Krasting, Graeme A MacGilchrist, Geoffrey Stanley, Jan-Erik Tesdal, and Jan D Zika, March 2025: Water mass transformation budgets in finite-volume generalized vertical coordinate ocean models. Journal of Advances in Modeling Earth Systems, 17(3), DOI:10.1029/2024MS004383. Abstract
Water Mass Transformation (WMT) theory provides conceptual tools that in principle enable innovative analyses of numerical ocean models; in practice, however, these methods can be challenging to implement and interpret, and therefore remain under-utilized. Our aim is to demonstrate the feasibility of diagnosing all terms in the water mass budget and to exemplify their usefulness for scientific inquiry and model development by quantitatively relating water mass changes, overturning circulations, boundary fluxes, and interior mixing. We begin with a pedagogical derivation of key results of classical WMT theory. We then describe best practices for diagnosing each of the water mass budget terms from the output of Finite-Volume Generalized Vertical Coordinate (FV-GVC) ocean models, including the identification of a non-negligible remainder term as the spurious numerical mixing due to advection scheme discretization errors. We illustrate key aspects of the methodology through the analysis of a polygonal region of the Greater Baltic Sea in a regional demonstration simulation using the Modular Ocean Model v6 (MOM6). We verify the convergence of our WMT diagnostics by brute-force, comparing time-averaged (“offline”) diagnostics on various vertical grids to timestep-averaged (“online”) diagnostics on the native model grid. Finally, we briefly describe a stack of xarray-enabled Python packages for evaluating WMT budgets in FV-GVC models (culminating in the new xwmb package), which is intended to be model-agnostic and available for community use and development.
Haine, Thomas W., Stephen M Griffies, Geoffrey Gebbie, and Wenrui Jiang, July 2025: A review of Green's function methods for tracer timescales and pathways in ocean models. Journal of Advances in Modeling Earth Systems, 17(7), DOI:10.1029/2024MS004637. Abstract
Understanding advective-diffusive dispersal of trace substances in environmental fluids like the global ocean is a ubiquitous challenge in geophysics. Since the turn of the millennium, substantial progress has been made in the theory, implementation in models, and application of such tracers in oceanography. For the first time, this progress is reviewed here in a synthetic way. We focus on tracer techniques in ocean models, including real and virtual tracers that diagnose timescale information, and we emphasize the connection to the Green's function that solves the advection-diffusion equation. Implementation of these techniques in ocean models is explained in an accessible way. We present example applications of these techniques to questions concerning ocean circulation, transport of biogeochemicals, and paleoceanography, including future opportunities.
Lobo, Matthew, Stephen M Griffies, and Wenda Zhang, April 2025: Vertical structure of baroclinic instability in a three-layer quasi-geostrophic model over a sloping bottom. Journal of Physical Oceanography, 55(4), DOI:10.1175/JPO-D-24-0130.1341-359. Abstract
We use a three-layer quasigeostrophic model to study the effects of sloping bottom topography on baroclinic instability in an ocean-like setting. Three layers allow for an interior potential vorticity gradient, a feature that is precluded by the use of a two-layer model. Material conservation of quasigeostrophic potential vorticity (QG PV) is expressed in terms of sloping background density interfaces, perturbations to the sloping interfaces, and the bottom slope. Linear stability analysis and numerical simulations of initial linear growth are used to demonstrate the dependence of topographic modifications to baroclinically unstable modes on background shear and stratification. Instabilities are classified according to their vertical structure, and the nature of topographic modification to baroclinically unstable modes is shown to correspond to instability type. Surface-intensified vertical shear of the background flow supports surface-intensified baroclinic instability that is less sensitive to bottom topography, compared to uniformly sheared background flow. When background shear is constant or close-to-constant, topography has a leading order effect on the vertical structure of the most baroclinically unstable mode. Upper-ocean stratification changes the lateral and vertical scales of the most unstable mode, with weaker upper-ocean stratification supporting high-wavenumber surface-intensified baroclinically unstable modes that are less sensitive to topography. We also present and assess a novel criterion that predicts the degree of surface intensification of the fastest-growing baroclinically unstable mode using only prescribed background properties, highlighting the role of interior potential vorticity in setting the most unstable mode’s vertical structure.
Gregory, Jonathan M., Jonah Bloch-Johnson, Matthew P Couldrey, Eleftheria Exarchou, Stephen M Griffies, Till Kuhlbrodt, Emily R Newsom, Oleg A Saenko, Tatsuo Suzuki, Quran Wu, Shogo Urakawa, and Laure Zanna, March 2024: A new conceptual model of global ocean heat uptake. Climate Dynamics, 62, DOI:10.1007/s00382-023-06989-z1669-1713. Abstract
We formulate a new conceptual model, named “MT2”, to describe global ocean heat uptake, as simulated by atmosphere–ocean general circulation models (AOGCMs) forced by increasing atmospheric CO
, as a function of global-mean surface temperature change T and the strength of the Atlantic meridional overturning circulation (AMOC, M). MT2 has two routes whereby heat reaches the deep ocean. On the basis of circumstantial evidence, we hypothetically identify these routes as low- and high-latitude. In low latitudes, which dominate the global-mean energy balance, heat uptake is temperature-driven and described by the two-layer model, with global-mean T as the temperature change of the upper layer. In high latitudes, a proportion p (about 14%) of the forcing is taken up along isopycnals, mostly in the Southern Ocean, nearly like a passive tracer, and unrelated to T. Because the proportion p depends linearly on the AMOC strength in the unperturbed climate, we hypothesise that high-latitude heat uptake and the AMOC are both affected by some characteristic of the unperturbed global ocean state, possibly related to stratification. MT2 can explain several relationships among AOGCM projections, some found in this work, others previously reported:
Ocean heat uptake efficiency correlates strongly with the AMOC.
Global ocean heat uptake is not correlated with the AMOC.
Transient climate response (TCR) is anticorrelated with the AMOC.
T projected for the late twenty-first century under high-forcing scenarios correlates more strongly with the effective climate sensitivity than with the TCR.
Jansen, Malte, Alistair Adcroft, Stephen M Griffies, and Ian Grooms, December 2024: The averaged hydrostatic Boussinesq ocean equations in generalized vertical coordinates. Journal of Advances in Modeling Earth Systems, 16(12), DOI:10.1029/2024MS004506. Abstract
Due to their limited resolution, numerical ocean models need to be interpreted as representing filtered or averaged equations. How to interpret models in terms of formally averaged equations, however, is not always clear, particularly in the case of hybrid or generalized vertical coordinate models, which limits our ability to interpret the model results and to develop parameterizations for the unresolved eddy contributions. We here derive the averaged hydrostatic Boussinesq equations in generalized vertical coordinates for an arbitrary thickness-weighted average. We then consider various special cases and discuss the extent to which the averaged equations are consistent with existing ocean model formulations. As previously discussed, the momentum equations in existing depth-coordinate models are best interpreted as representing Eulerian averages (i.e., averages taken at fixed depth), while the tracer equations can be interpreted as either Eulerian or thickness-weighted isopycnal averages. Instead we find that no averaging is fully consistent with existing formulations of the parameterizations in semi-Lagrangian discretizations of generalized vertical coordinate ocean models such as MOM6. A coordinate-following average would require “coordinate-aware” parameterizations that can account for the changing nature of the eddy terms as the coordinate changes. Alternatively, the model variables can be interpreted as representing either Eulerian or (thickness-weighted) isopycnal averages, independent of the model coordinate that is being used for the numerical discretization. Existing parameterizations in generalized vertical coordinate models, however, are not always consistent with either of these interpretations, which, respectively, would require a three-dimensional divergence-free eddy tracer advection or a form-stress parameterization in the momentum equations.
Khatri, Hermant, Stephen M Griffies, Benjamin A Storer, Michele Buzzicotti, Hussein Aluie, Maike Sonnewald, Raphael Dussin, and Andrew Shao, June 2024: A scale-dependent analysis of the barotropic vorticity budget in a global ocean simulation. Journal of Advances in Modeling Earth Systems, 16(6), DOI:10.1029/2023MS003813. Abstract
The climatological mean barotropic vorticity budget is analyzed to investigate the relative importance of surface wind stress, topography, planetary vorticity advection, and nonlinear advection in dynamical balances in a global ocean simulation. In addition to a pronounced regional variability in vorticity balances, the relative magnitudes of vorticity budget terms strongly depend on the length-scale of interest. To carry out a length-scale dependent vorticity analysis in different ocean basins, vorticity budget terms are spatially coarse-grained. At length-scales greater than 1,000 km, the dynamics closely follow the Topographic-Sverdrup balance in which bottom pressure torque, surface wind stress curl and planetary vorticity advection terms are in balance. In contrast, when including all length-scales resolved by the model, bottom pressure torque and nonlinear advection terms dominate the vorticity budget (Topographic-Nonlinear balance), which suggests a prominent role of oceanic eddies, which are of km in size, and the associated bottom pressure anomalies in local vorticity balances at length-scales smaller than 1,000 km. Overall, there is a transition from the Topographic-Nonlinear regime at scales smaller than 1,000 km to the Topographic-Sverdrup regime at length-scales greater than 1,000 km. These dynamical balances hold across all ocean basins; however, interpretations of the dominant vorticity balances depend on the level of spatial filtering or the effective model resolution. On the other hand, the contribution of bottom and lateral friction terms in the barotropic vorticity budget remains small and is significant only near sea-land boundaries, where bottom stress and horizontal viscous friction generally peak.
Krasting, John P., Stephen M Griffies, Jan-Erik Tesdal, Graeme A MacGilchrist, Rebecca L Beadling, and Christopher M Little, December 2024: Steric sea level rise and relationships with model drift and water mass representation in GFDL CM4 and ESM4. Journal of Climate, 37(24), DOI:10.1175/JCLI-D-23-0591.1. Abstract
Density-driven steric seawater changes are a leading-order contributor to global mean sea level rise. However, intermodel differences in the magnitude and spatial patterns of steric sea level rise exist at regional scales and often emerge during the spinup and preindustrial control integrations of climate models. Steric sea level results from an eddy-permitting climate model, GFDL CM4, are compared with a lower-resolution counterpart, GFDL-ESM4. The results from both models are examined through basin-scale heat budgets and watermass analysis, and we compare the patterns of ocean heat uptake, redistribution, and sea level differ in ocean-only [i.e., Ocean Model Intercomparison Project (OMIP)] and coupled climate configurations. After correcting for model drift, both GFDL CM4 and GFDL-ESM4 simulate nearly equivalent ocean heat content change and global sea level rise during the historical period. However, the GFDL CM4 model exhibits as much as a 40% increase in surface ocean heat uptake in the Southern Ocean and subsequent increases in horizontal export to other ocean basins after bias correction. The results suggest regional differences in the processes governing Southern Ocean heat export, such as the formation of Antarctic Intermediate Water (AAIW), Subpolar Mode Water (SPMW), and gyre transport between the two models, and that sea level changes in these models cannot be fully bias-corrected. Since the process-level differences between the two models are evident in the preindustrial control simulations of both models, these results suggest that the control simulations are important for identifying and correcting sea level–related model biases.
Neme, Julia, Matthew H England, Andrew McC Hogg, Hermant Khatri, and Stephen M Griffies, January 2024: The role of bottom friction in mediating the response of the Weddell Gyre circulation to changes in surface stress and buoyancy fluxes. Journal of Physical Oceanography, 54(1), DOI:10.1175/JPO-D-23-0165.1217-236. Abstract
The Weddell Gyre is one of the dominant features of the Southern Ocean circulation and its dynamics have been linked to processes of climatic relevance. Variability in the strength of the gyre’s horizontal transport has been linked to heat transport toward the Antarctic margins and changes in the properties and rates of export of bottom waters from the Weddell Sea region to the abyssal global ocean. However, the precise physical mechanisms that force variability in the Weddell’s lateral circulation across different time scales remain unknown. In this study, we use a barotropic vorticity budget from a mesoscale eddy active model simulation to attribute changes in gyre strength to variability in possible driving processes. We find that the Weddell Gyre’s circulation is sensitive to bottom friction associated with the overflowing dense waters at its western boundary. In particular, an increase in the production of dense waters at the southwestern continental shelf strengthens the bottom flow at the gyre’s western boundary, yet this drives a weakening of the depth-integrated barotropic circulation via increased bottom friction. Strengthening surface winds initially accelerate the gyre, but within a few years the response reverses once dense water production and export increases. These results reveal that the gyre can weaken in response to stronger surface winds, putting into question the traditional assumption of a direct relationship between surface stress and gyre strength in regions where overflowing dense water forms part of the depth-integrated flow.
Prend, Channing J., Graeme A MacGilchrist, Georgy E Manucharyan, Rachel Q Pang, Ruth Moorman, Andrew F Thompson, Stephen M Griffies, Matthew R Mazloff, Lynne D Talley, and Sarah T Gille, January 2024: Ross Gyre variability modulates oceanic heat supply toward the West Antarctic continental shelf. Communications Earth and Environment, 5, 47, DOI:10.1038/s43247-024-01207-y. Abstract
West Antarctic Ice Sheet mass loss is a major source of uncertainty in sea level projections. The primary driver of this melting is oceanic heat from Circumpolar Deep Water originating offshore in the Antarctic Circumpolar Current. Yet, in assessing melt variability, open ocean processes have received considerably less attention than those governing cross-shelf exchange. Here, we use Lagrangian particle release experiments in an ocean model to investigate the pathways by which Circumpolar Deep Water moves toward the continental shelf across the Pacific sector of the Southern Ocean. We show that Ross Gyre expansion, linked to wind and sea ice variability, increases poleward heat transport along the gyre’s eastern limb and the relative fraction of transport toward the Amundsen Sea. Ross Gyre variability, therefore, influences oceanic heat supply toward the West Antarctic continental slope. Understanding remote controls on basal melt is necessary to predict the ice sheet response to anthropogenic forcing.
Deficiencies in upper ocean vertical mixing parameterizations contribute to tropical upper ocean biases in global coupled general circulation models, affecting their simulated ocean heat uptake and ENSO variability. To better understand these deficiencies, we develop a suite of ocean model experiments including both idealized single column models and realistic global simulations. The vertical mixing parameterizations are first evaluated using large eddy simulations as a baseline to assess uncertainties and evaluate their implied turbulent mixing. Global models are then developed following NOAA/GFDL's 0.25° nominal horizontal grid spacing OM4 (uncoupled) configuration of the MOM6 ocean model, with various modifications that target biases in the original model. We test several enhancements to the existing mixing schemes and evaluate them against observational constraints from Tropical Atmosphere Ocean moorings and Argo floats. In particular, we find that we can improve the diurnal variability of mixing in OM4 via modifications to its surface boundary layer mixing scheme, and can improve the net mixing in the upper thermocline by reducing the background vertical viscosity, allowing for more realistic, less diffuse currents. The improved OM4 model better represents the mixing, leading to improved diurnal deep-cycle variability, a more realistic time-mean tropical thermocline structure, and a better Pacific Equatorial Undercurrent.
Using a recently developed 1/12th degree regional ocean model, we establish a link between U.S. East Coast sea level variability and offshore upper ocean heat content change. This link manifests as a cross-shore mass redistribution driven by an offshore thermosteric sea level response to subsurface warming or cooling. Approximately 50% of simulated monthly to interannual coastal sea level variance south of Cape Hatteras can be statistically accounted for by this mechanism, realized as a function of regional ocean hypsometry, gyre scale warming, and the depth dependence of density change. This response to offshore warming explains the nonstationarity of U.S. East Coast sea level covariance, a specifically observed and modeled behavior after ~ 2010. Since approximately 2010, elevated rates of sea level rise south of Cape Hatteras can be partly explained as the result of shoreward mass redistribution due to offshore subsurface warming within the North Atlantic subtropical gyre. These results reveal a mechanism that connects local coastal sea level to a broader region and identifies the influence of regional heat content changes on coastal sea level. This analysis presents a framework for identifying new regions that may be susceptible to enhanced sea level rise due to ocean warming and helps bridge the gap between quantifying large scale change and anticipating local coastal impacts that can make flooding and storm surge more acutely damaging.
Zhang, Wenda, Stephen M Griffies, Robert Hallberg, Yi-Hung Kuo, and Christopher L P Wolfe, June 2024: The role of surface potential vorticity in the vertical structure of mesoscale eddies in wind-driven ocean circulations. Journal of Physical Oceanography, 54(6), DOI:10.1175/JPO-D-23-0203.11243-1266. Abstract
The vertical structure of ocean eddies is generally surface-intensified, commonly attributed to the dominant baroclinic modes arising from the boundary conditions (BCs). Conventional BC considerations mostly focus on either flat- or rough-bottom conditions. The impact of surface buoyancy anomalies—often represented by surface potential vorticity (PV) anomalies—has not been fully explored. Here, we study the role of the surface PV in setting the vertical distribution of eddy kinetic energy (EKE) in an idealized adiabatic ocean model driven by wind stress. The simulated EKE profile in the extratropical ocean tends to peak at the surface and have an e-folding depth typically smaller than half of the ocean depth. This vertical structure can be reasonably represented by a single surface quasigeostrophic (SQG) mode at the energy-containing scale resulting from the large-scale PV structure. Due to isopycnal outcropping and interior PV homogenization, the surface meridional PV gradient is substantially stronger than the interior PV gradient, yielding surface-trapped baroclinically unstable modes with horizontal scales comparable to or smaller than the deformation radius. These surface-trapped eddies then grow in size both horizontally and vertically through an inverse energy cascade up to the energy-containing scale, which dominates the vertical distribution of EKE. As for smaller horizontal scales, the EKE distribution decays faster with depth. Guided by this interpretation, an SQG-based scale-aware parameterization of the EKE profile is proposed. Preliminary offline diagnosis of a high-resolution simulation shows the proposed scheme successfully reproducing the dependence of the vertical structure of EKE on the horizontal grid resolution.
Buzzicotti, Michele, Benjamin A Storer, Hermant Khatri, Stephen M Griffies, and Hussein Aluie, June 2023: Spatio-temporal coarse-graining decomposition of the global ocean geostrophic kinetic energy. Journal of Advances in Modeling Earth Systems, 15(6), DOI:10.1029/2023MS003693. Abstract
We expand on a recent determination of the first global energy spectrum of the ocean's surface geostrophic circulation (Storer et al., 2022, https://doi.org/10.1038/s41467-022-33031-3) using a coarse-graining (CG) method. We compare spectra from CG to those from spherical harmonics by treating land in a manner consistent with the boundary conditions. While the two methods yield qualitatively consistent domain-averaged results, spherical harmonics spectra are too noisy at gyre-scales (>1,000 km). More importantly, spherical harmonics are inherently global and cannot provide local information connecting scales with currents geographically. CG shows that the extra-tropics mesoscales (100–500 km) have a root-mean-square (rms) velocity of ∼15 cm/s, which increases to ∼30–40 cm/s locally in the Gulf Stream and Kuroshio and to ∼16–28 cm/s in the ACC. There is notable hemispheric asymmetry in mesoscale energy-per-area, which is higher in the north due to continental boundaries. We estimate that ≈25%–50% of total geostrophic energy is at scales smaller than 100 km, and is un(der)-resolved by pre-SWOT satellite products. Spectra of the time-mean circulation show that most of its energy (up to 70%) resides in stationary eddies with characteristic scales smaller than (<500 km). This highlights the preponderance of “standing” small-scale structures in the global ocean due to the temporally coherent forcing by boundaries. By coarse-graining in space and time, we compute the first spatio-temporal global spectrum of geostrophic circulation from AVISO and NEMO. These spectra show that every length-scale evolves over a wide range of time-scales with a consistent peak at ≈200 km and ≈2–3 weeks.
The use of coarse resolution and strong grid-scale dissipation has prevented global ocean models from simulating the correct kinetic energy level. Recently parameterizing energy backscatter has been proposed to energize the model simulations. Parameterizing backscatter reduces long-standing North Atlantic sea surface temperature (SST) and associated surface current biases, but the underlying mechanism remains unclear. Here, we apply backscatter in different geographic regions to distinguish the different physical processes at play. We show that an improved Gulf Stream path is due to backscatter acting north of the Grand Banks to maintain a strong deep western boundary current. An improved North Atlantic Current path is due to backscatter acting around the Flemish Cap, with likely an improved nearby topography-flow interactions. These results suggest that the SST improvement with backscatter is partly due to the resulted strengthening of resolved currents, whereas the role of improved eddy physics requires further research.
Falasca, Fabrizio, Andrew Brettin, Laure Zanna, Stephen M Griffies, Jianjun Yin, and Ming Zhao, June 2023: Exploring the nonstationarity of coastal sea level probability distributions. Environmental Data Science, 2, e16, DOI:10.1017/eds.2023.10. Abstract
Studies agree on a significant global mean sea level rise in the 20th century and its recent 21st century acceleration in the satellite record. At regional scale, the evolution of sea level probability distributions is often assumed to be dominated by changes in the mean. However, a quantification of changes in distributional shapes in a changing climate is currently missing. To this end, we propose a novel framework quantifying significant changes in probability distributions from time series data. The framework first quantifies linear trends in quantiles through quantile regression. Quantile slopes are then projected onto a set of four orthogonal polynomials quantifying how such changes can be explained by independent shifts in the first four statistical moments. The framework proposed is theoretically founded, general, and can be applied to any climate observable with close-to-linear changes in distributions. We focus on observations and a coupled climate model (GFDL-CM4). In the historical period, trends in coastal daily sea level have been driven mainly by changes in the mean and can therefore be explained by a shift of the distribution with no change in shape. In the modeled world, robust changes in higher order moments emerge with increasing concentration. Such changes are driven in part by ocean circulation alone and get amplified by sea level pressure fluctuations, with possible consequences for sea level extremes attribution studies.
Loose, Nora, Gustavo Marques, Alistair Adcroft, Scott D Bachman, Stephen M Griffies, Ian Grooms, Robert Hallberg, and Malte Jansen, December 2023: Comparing two parameterizations for the restratification effect of mesoscale eddies in an isopycnal ocean model. Journal of Advances in Modeling Earth Systems, 15(12), DOI:10.1029/2022MS003518. Abstract
There are two distinct parameterizations for the restratification effect of mesoscale eddies: the Greatbatch and Lamb (1990, GL90, https://journals.ametsoc.org/view/journals/phoc/20/10/1520-0485_1990_020_1634_opvmom_2_0_co_2.xml?tab_body=abstract-display) parameterization, which mixes horizontal momentum in the vertical, and the Gent and McWilliams (1990, GM90, https://journals.ametsoc.org/view/journals/phoc/20/1/1520-0485_1990_020_0150_imiocm_2_0_co_2.xml) parameterization, which flattens isopycnals adiabatically. Even though these two parameterizations are effectively equivalent under the assumption of quasi-geostrophy, GL90 has been used much less than GM90, and exclusively in z-coordinate models. In this paper, we compare the GL90 and GM90 parameterizations in an idealized isopycnal coordinate model, both from a theoretical and practical perspective. From a theoretical perspective, GL90 is more attractive than GM90 for isopycnal coordinate models because GL90 provides an interpretation that is fully consistent with thickness-weighted isopycnal averaging, while GM90 cannot be entirely reconciled with any fully isopycnal averaging framework. From a practical perspective, the GL90 and GM90 parameterizations lead to extremely similar energy levels, flow and vertical structure, even though their energetic pathways are very different. The striking resemblance between the GL90 and GM90 simulations persists from non-eddying through eddy-permitting resolution. We conclude that GL90 is a promising alternative to GM90 for isopycnal coordinate models, where it is more consistent with theory, computationally more efficient, easier to implement, and numerically more stable. Assessing the applicability of GL90 in realistic global ocean simulations with hybrid coordinate schemes should be a priority for future work.
Morrison, Adele K., Wilma G C Huneke, Julia Neme, Paul Spence, Andrew McC Hogg, Matthew H England, and Stephen M Griffies, September 2023: Sensitivity of Antarctic shelf waters and abyssal overturning to local winds. Journal of Climate, 36(18), DOI:10.1175/JCLI-D-22-0858.16465-6479. Abstract
Winds around the Antarctic continental margin are known to exert a strong control on the local ocean stratification and circulation. However, past work has largely focused on the ocean response to changing winds in limited regional sectors and the circumpolar dynamical response to polar wind change remains uncertain. In this work, we use a high-resolution global ocean–sea ice model to investigate how dense shelf water formation and the temperature of continental shelf waters respond to changes in the zonal and meridional components of the polar surface winds. Increasing the zonal easterly wind component drives an enhanced southward Ekman transport in the surface layer, raising sea level over the continental shelf and deepening coastal isopycnals. The downward isopycnal movement cools the continental shelf, as colder surface waters replace warmer waters below. However, in this model the zonal easterly winds do not impact the strength of the abyssal overturning circulation, in contrast to past idealized model studies. Instead, increasing the meridional wind speed strengthens the abyssal overturning circulation via a sea ice advection mechanism. Enhanced offshore meridional wind speed increases the northward export of sea ice, resulting in decreased sea ice thickness over the continental shelf. The reduction in sea ice coverage leads to increased air–sea heat loss, sea ice formation, brine rejection, dense shelf water formation, and abyssal overturning circulation. Increasing the meridional winds causes warming at depth over most of the continental shelf, due to a heat advection feedback associated with the enhanced overturning circulation.
Storer, Benjamin A., Michele Buzzicotti, Hermant Khatri, Stephen M Griffies, and Hussein Aluie, December 2023: Global cascade of kinetic energy in the ocean and the atmospheric imprint. Science Advances, 9(51), DOI:10.1126/sciadv.adi7420. Abstract
Here, we present an estimate for the ocean's global scale transfer of kinetic energy (KE), across scales from 10 to 40,000 km. Oceanic KE transfer between gyre scales and mesoscales is induced by the atmosphere’s Hadley, Ferrel, and polar cells, and the intertropical convergence zone induces an intense downscale KE transfer. Upscale transfer peaks at 300 gigawatts across mesoscales of 120 km in size, roughly one-third the energy input by winds into the oceanic general circulation. Nearly three quarters of this “cascade” occurs south of 15°S and penetrates almost the entire water column. The mesoscale cascade has a self-similar seasonal cycle with characteristic lag time of ≈27 days per octave of length scales; transfer across 50 km peaks in spring, while transfer across 500 km peaks in summer. KE of those mesoscales follows the same cycle but peaks ≈40 days after the peak cascade, suggesting that energy transferred across a scale is primarily deposited at a scale four times larger.
Swart, Neil C., Torge Martin, Rebecca L Beadling, Jia-Jia Chen, Christopher Danek, Matthew H England, Riccardo Farneti, Stephen M Griffies, Tore Hattermann, Judith Hauck, F Alexander Haumann, André Jüling, Qian Li, John Marshall, Morven Muilwijk, Andrew G Pauling, Ariaan Purich, Inga J Smith, and Max Thomas, December 2023: The Southern Ocean Freshwater Input from Antarctica (SOFIA) Initiative: Scientific objectives and experimental design. Geoscientific Model Development, 16(24), DOI:10.5194/gmd-16-7289-20237289–7309. Abstract
As the climate warms, the grounded ice sheet and floating ice shelves surrounding Antarctica are melting and releasing additional freshwater into the Southern Ocean. Nonetheless, almost all existing coupled climate models have fixed ice sheets and lack the physics required to represent the dominant sources of Antarctic melt. These missing ice dynamics represent a key uncertainty that is typically unaccounted for in current global climate change projections. Previous modelling studies that have imposed additional Antarctic meltwater have demonstrated regional impacts on Southern Ocean stratification, circulation, and sea ice, as well as remote changes in atmospheric circulation, tropical precipitation, and global temperature. However, these previous studies have used widely varying rates of freshwater forcing, have been conducted using different climate models and configurations, and have reached differing conclusions on the magnitude of meltwater–climate feedbacks. The Southern Ocean Freshwater Input from Antarctica (SOFIA) initiative brings together a team of scientists to quantify the climate system response to Antarctic meltwater input along with key aspects of the uncertainty. In this paper, we summarize the state of knowledge on meltwater discharge from the Antarctic ice sheet and ice shelves to the Southern Ocean and explain the scientific objectives of our initiative. We propose a series of coupled and ocean–sea ice model experiments, including idealized meltwater experiments, historical experiments with observationally consistent meltwater input, and future scenarios driven by meltwater inputs derived from stand-alone ice sheet models. Through coordinating a multi-model ensemble of simulations using a common experimental design, open data archiving, and facilitating scientific collaboration, SOFIA aims to move the community toward better constraining our understanding of the climate system response to Antarctic melt.
Tesdal, Jan-Erik, Graeme A MacGilchrist, Rebecca L Beadling, Stephen M Griffies, John P Krasting, and Paul J Durack, March 2023: Revisiting Interior Water Mass Responses to Surface Forcing Changes and the Subsequent Effects on Overturning in the Southern Ocean. Journal of Geophysical Research: Oceans, 128(3), DOI:10.1029/2022JC019105. Abstract
Two coupled climate models, differing primarily in horizontal resolution and treatment of mesoscale eddies, were used to assess the impact of perturbations in wind stress and Antarctic ice sheet (AIS) melting on the Southern Ocean meridional overturning circulation (SO MOC), which plays an important role in global climate regulation. The largest impact is found in the SO MOC lower limb, associated with the formation of Antarctic Bottom Water (AABW), which in both models is enhanced by wind and weakened by AIS meltwater perturbations. Even though both models under the AIS melting perturbation show similar AABW transport reductions of 4–5 Sv (50%–60%), the volume deflation of AABW south of 30°S is four times greater in the higher resolution simulation (−20 vs. −5 Sv). Water mass transformation (WMT) analysis reveals that surface-forced dense water formation on the Antarctic shelf is absent in the higher resolution and reduced by half in the lower resolution model in response to the increased AIS melting. However, the decline of the AABW volume (and its inter-model difference) far exceeds the surface-forced WMT changes alone, which indicates that the divergent model responses arise from interactions between changes in surface forcing and interior mixing processes. This model divergence demonstrates an important source of uncertainty in climate modeling, and indicates that accurate shelf processes together with scenarios accounting for AIS melting are necessary for robust projections of the deep ocean's response to anthropogenic forcing.
We highlight the differing roles of vorticity and strain in the transport of coarse-grained scalars at length scales larger than ℓ by smaller scale (subscale) turbulence. We use the first term in a multiscale gradient expansion due to Eyink [J. Fluid Mech. 549, 159 (2006)], which exhibits excellent correlation with the exact subscale physics when the partitioning length ℓ is any scale smaller than that of the spectral peak. We show that unlike subscale strain, which acts as an anisotropic diffusion/antidiffusion tensor, subscale vorticity's contribution is solely a conservative advection of coarse-grained quantities by an eddy-induced nondivergent velocity, v∗, that is proportional to the curl of vorticity. Therefore, material (Lagrangian) advection of coarse-grained quantities is accomplished not by the coarse-grained flow velocity,¯¯¯uℓ, but by the effective velocity,
¯¯¯uℓ+v∗, the physics of which may improve commonly used LES models.
We use two coupled climate models, GFDL-CM4 and GFDL-ESM4, to investigate the physical response of the Southern Ocean to changes in surface wind stress, Antarctic meltwater, and the combined forcing of the two in a pre-industrial control simulation. The meltwater cools the ocean surface in all regions except the Weddell Sea, where the wind stress warms the near-surface layer. The limited sensitivity of the Weddell Sea surface layer to the meltwater is due to the spatial distribution of the meltwater fluxes, regional bathymetry, and large-scale circulation patterns. The meltwater forcing dominates the Antarctic shelf response and the models yield strikingly different responses along West Antarctica. The disagreement is attributable to the mean-state representation and meltwater-driven acceleration of the Antarctic Slope Current (ASC). In CM4, the meltwater is efficiently trapped on the shelf by a well resolved, strong, and accelerating ASC which isolates the West Antarctic shelf from warm offshore waters, leading to strong subsurface cooling. In ESM4, a weaker and diffuse ASC allows more meltwater to escape to the open ocean, the West Antarctic shelf does not become isolated, and instead strong subsurface warming occurs. The CM4 results suggest a possible negative feedback mechanism that acts to limit future melting, while the ESM4 results suggest a possible positive feedback mechanism that acts to accelerate melt. Our results demonstrate the strong influence the ASC has on governing changes along the shelf, highlighting the importance of coupling interactive ice sheet models to ocean models that can resolve these dynamical processes.
Light, Charles X., Brian K Arbic, Paige E Martin, Laurent Brodeau, J Thomas Farrar, and Stephen M Griffies, et al., March 2022: Effects of grid spacing on high-frequency precipitation variance in coupled high-resolution global ocean–atmosphere models. Climate Dynamics, DOI:10.1007/s00382-022-06257-6. Abstract
High-frequency precipitation variance is calculated in 12 different free-running (non-data-assimilative) coupled high resolution atmosphere–ocean model simulations, an assimilative coupled atmosphere–ocean weather forecast model, and an assimilative reanalysis. The results are compared with results from satellite estimates of precipitation and rain gauge observations. An analysis of irregular sub-daily fluctuations, which was applied by Covey et al. (Geophys Res Lett 45:12514–12522, 2018. https://doi.org/10.1029/2018GL078926) to satellite products and low-resolution climate models, is applied here to rain gauges and higher-resolution models. In contrast to lower-resolution climate simulations, which Covey et al. (2018) found to be lacking with respect to variance in irregular sub-daily fluctuations, the highest-resolution simulations examined here display an irregular sub-daily fluctuation variance that lies closer to that found in satellite products. Most of the simulations used here cannot be analyzed via the Covey et al. (2018) technique, because they do not output precipitation at sub-daily intervals. Thus the remainder of the paper focuses on frequency power spectral density of precipitation and on cumulative distribution functions over time scales (2–100 days) that are still relatively “high-frequency” in the context of climate modeling. Refined atmospheric or oceanic model grid spacing is generally found to increase high-frequency precipitation variance in simulations, approaching the values derived from observations. Mesoscale-eddy-rich ocean simulations significantly increase precipitation variance only when the atmosphere grid spacing is sufficiently fine (< 0.5°). Despite the improvements noted above, all of the simulations examined here suffer from the “drizzle effect”, in which precipitation is not temporally intermittent to the extent found in observations.
We describe an idealized primitive-equation model for studying mesoscale turbulence and leverage a hierarchy of grid resolutions to make eddy-resolving calculations on the finest grids more affordable. The model has intermediate complexity, incorporating basin-scale geometry with idealized Atlantic and Southern oceans and with non-uniform ocean depth to allow for mesoscale eddy interactions with topography. The model is perfectly adiabatic and spans the Equator and thus fills a gap between quasi-geostrophic models, which cannot span two hemispheres, and idealized general circulation models, which generally include diabatic processes and buoyancy forcing. We show that the model solution is approaching convergence in mean kinetic energy for the ocean mesoscale processes of interest and has a rich range of dynamics with circulation features that emerge only due to resolving mesoscale turbulence.
Naveira Garabato, Alberto C., Xiaolong Yu, Jörn Callies, Roy Barkan, Kurt L Polzin, Eleanor E Frajka-Williams, Christian E Buckingham, and Stephen M Griffies, January 2022: Kinetic energy transfers between mesoscale and submesoscale motions in the open ocean's upper layers. Journal of Physical Oceanography, 52(1), DOI:10.1175/JPO-D-21-0099.175-97. Abstract
Mesoscale eddies contain the bulk of the ocean’s kinetic energy (KE), but fundamental questions remain on the cross-scale KE transfers linking eddy generation and dissipation. The role of submesoscale flows represents the key point of discussion, with contrasting views of submesoscales as either a source or a sink of mesoscale KE. Here, the first observational assessment of the annual cycle of the KE transfer between mesoscale and submesoscale motions is performed in the upper layers of a typical open-ocean region. Although these diagnostics have marginal statistical significance and should be regarded cautiously, they are physically plausible and can provide a valuable benchmark for model evaluation. The cross-scale KE transfer exhibits two distinct stages, whereby submesoscales energize mesoscales in winter and drain mesoscales in spring. Despite this seasonal reversal, an inverse KE cascade operates throughout the year across much of the mesoscale range. Our results are not incompatible with recent modeling investigations that place the headwaters of the inverse KE cascade at the submesoscale, and that rationalize the seasonality of mesoscale KE as an inverse cascade-mediated response to the generation of submesoscales in winter. However, our findings may challenge those investigations by suggesting that, in spring, a downscale KE transfer could dampen the inverse KE cascade. An exploratory appraisal of the dynamics governing mesoscale–submesoscale KE exchanges suggests that the upscale KE transfer in winter is underpinned by mixed layer baroclinic instabilities, and that the downscale KE transfer in spring is associated with frontogenesis. Current submesoscale-permitting ocean models may substantially understate this downscale KE transfer, due to the models’ muted representation of frontogenesis.
Turbulent mixing in the ocean surface boundary layer leads to the presence of a surface mixed layer. This mixed layer is important for many phenomena including large-scale ocean dynamics, ocean-atmosphere coupling, and biological and biogeochemical processes. Analysis of the ocean mixed layer requires one to estimate its vertical extent, for which there are various definitions. Correspondingly, there are uncertainties on how to best identify an ocean surface mixed layer for a given application. We propose defining the mixed layer depth (MLD) from energetic principles through the potential energy (PE). The PE based MLD is based on the concept of PE anomaly, which measures the stratification of a layer of seawater by estimating its energetic distance from a well-mixed state. We apply the PE anomaly to diagnose the MLD as the depth to which a given energy could homogenize a layer of seawater. We evaluate the MLD defined by common existing methods and demonstrate that they contain a wide range of PE anomalies for the same MLD, particularly evident for deep winter mixed layers. The MLD defined from the PE anomaly ensures a more consistent MLD identified for a large range of stratifications. Furthermore, the PE method relates to the turbulent kinetic energy budget of the ocean surface boundary layer, which is fundamental to upper ocean mixing processes and parameterizations. The resulting MLD is more representative of active boundary layer turbulence, and is more robust to small anomalies in seawater properties.
Storer, Benjamin A., Michele Buzzicotti, Hermant Khatri, Stephen M Griffies, and Hussein Aluie, September 2022: Global energy spectrum of the general oceanic circulation. Nature Communications, 13, 5314, DOI:10.1038/s41467-022-33031-3. Abstract
Advent of satellite altimetry brought into focus the pervasiveness of mesoscale eddies O(100) km in size, which are the ocean’s analogue of weather systems and are often regarded as the spectral peak of kinetic energy (KE). Yet, understanding of the ocean’s spatial scales has been derived mostly from Fourier analysis in small "representative” regions that cannot capture the vast dynamic range at planetary scales. Here, we use a coarse-graining method to analyze scales much larger than what had been possible before. Spectra spanning over three decades of length-scales reveal the Antarctic Circumpolar Current as the spectral peak of the global extra-tropical circulation, at ≈ 104 km, and a previously unobserved power-law scaling over scales larger than 103 km. A smaller spectral peak exists at ≈ 300 km associated with mesoscales, which, due to their wider spread in wavenumber space, account for more than 50% of resolved surface KE globally. Seasonal cycles of length-scales exhibit a characteristic lag-time of ≈ 40 days per octave of length-scales such that in both hemispheres, KE at 102 km peaks in spring while KE at 103 km peaks in late summer. These results provide a new window for understanding the multiscale oceanic circulation within Earth’s climate system, including the largest planetary scales.
Vogt, Linus, Friedrich A Burger, Stephen M Griffies, and Thomas L Frölicher, May 2022: Local drivers of marine heatwaves: A global analysis with an Earth system model. Frontiers in Climate, DOI:10.3389/fclim.2022.847995. Abstract
Marine heatwaves (MHWs) are periods of extreme warm ocean temperatures that can have devastating impacts on marine organisms and socio-economic systems. Despite recent advances in understanding the underlying processes of individual events, a global view of the local oceanic and atmospheric drivers of MHWs is currently missing. Here, we use daily-mean output of temperature tendency terms from a comprehensive fully coupled coarse-resolution Earth system model to quantify the main local processes leading to the onset and decline of surface MHWs in different seasons. The onset of MHWs in the subtropics and mid-to-high latitudes is primarily driven by net ocean heat uptake associated with a reduction of latent heat loss in all seasons, increased shortwave heat absorption in summer and reduced sensible heat loss in winter, dampened by reduced vertical mixing from the non-local portion of the K-Profile Parameterization boundary layer scheme (KPP) especially in summer. In the tropics, ocean heat uptake is reduced and lowered vertical local mixing and diffusion cause the warming. In the subsequent decline phase, increased ocean heat loss to the atmosphere due to enhanced latent heat loss in all seasons together with enhanced vertical local mixing and diffusion in the high latitudes during summer dominate the temperature decrease globally. The processes leading to the onset and decline of MHWs are similar for short and long MHWs, but there are differences in the drivers between summer and winter. Different types of MHWs with distinct driver combinations are identified within the large variability among events. Our analysis contributes to a better understanding of MHW drivers and processes and may therefore help to improve the prediction of high-impact marine heatwaves.
The discrete baroclinic modes of quasigeostrophic theory are incomplete, and the incompleteness manifests as a loss of information in the projection process. The incompleteness of the baroclinic modes is related to the presence of two previously unnoticed stationary step-wave solutions of the Rossby wave problem with flat boundaries. These step waves are the limit of surface quasigeostrophic waves as boundary buoyancy gradients vanish. A complete normal-mode basis for quasigeostrophic theory is obtained by considering the traditional Rossby wave problem with prescribed buoyancy gradients at the lower and upper boundaries. The presence of these boundary buoyancy gradients activates the previously inert boundary degrees of freedom. These Rossby waves have several novel properties such as the presence of multiple modes with no internal zeros, a finite number of modes with negative norms, and the fact that their vertical structures form a basis capable of representing any quasigeostrophic state with a differentiable series expansion. These properties are a consequence of the Pontryagin-space setting of the Rossby wave problem in the presence of boundary buoyancy gradients (as opposed to the usual Hilbert-space setting). We also examine the quasigeostrophic vertical velocity modes and derive a complete basis for such modes as well. A natural application of these modes is the development of a weakly nonlinear wave-interaction theory of geostrophic turbulence that takes topography into account.
Yassin, Houssam, and Stephen M Griffies, December 2022: Surface quasigeostrophic turbulence in variable stratification. Journal of Physical Oceanography, 52(12), DOI:10.1175/JPO-D-22-0040.12995-3013. Abstract
Numerical and observational evidence indicates that, in regions where mixed layer instability is active, the surface geostrophic velocity is largely induced by surface buoyancy anomalies. Yet, in these regions, the observed surface kinetic energy spectrum is steeper than predicted by uniformly stratified surface quasigeostrophic theory. By generalizing surface quasigeostrophic theory to account for variable stratification, we show that surface buoyancy anomalies can generate a variety of dynamical regimes depending on the stratification’s vertical structure. Buoyancy anomalies generate longer-range velocity fields over decreasing stratification and shorter-range velocity fields over increasing stratification. As a result, the surface kinetic energy spectrum is steeper over decreasing stratification than over increasing stratification. An exception occurs if the near-surface stratification is much larger than the deep-ocean stratification. In this case, we find an extremely local turbulent regime with surface buoyancy homogenization and a steep surface kinetic energy spectrum, similar to equivalent barotropic turbulence. By applying the variable stratification theory to the wintertime North Atlantic, and assuming that mixed layer instability acts as a narrowband small-scale surface buoyancy forcing, we obtain a predicted surface kinetic energy spectrum between k−4/3 and k−7/3, which is consistent with the observed wintertime k−2 spectrum. We conclude by suggesting a method of measuring the buoyancy frequency’s vertical structure using satellite observations.
Couldrey, Matthew P., Jonathan M Gregory, Fabio Boeira Dias, Peter Dobrohotoff, Catia M Domingues, Oluwayemi Garuba, and Stephen M Griffies, et al., January 2021: What causes the spread of model projections of ocean dynamic sea-level change in response to greenhouse gas forcing?Climate Dynamics, 56, DOI:10.1007/s00382-020-05471-4155-187. Abstract
Sea levels of different atmosphere–ocean general circulation models (AOGCMs) respond to climate change forcing in different ways, representing a crucial uncertainty in climate change research. We isolate the role of the ocean dynamics in setting the spatial pattern of dynamic sea-level (ζ) change by forcing several AOGCMs with prescribed identical heat, momentum (wind) and freshwater flux perturbations. This method produces a ζ projection spread comparable in magnitude to the spread that results from greenhouse gas forcing, indicating that the differences in ocean model formulation are the cause, rather than diversity in surface flux change. The heat flux change drives most of the global pattern of ζ change, while the momentum and water flux changes cause locally confined features. North Atlantic heat uptake causes large temperature and salinity driven density changes, altering local ocean transport and ζ. The spread between AOGCMs here is caused largely by differences in their regional transport adjustment, which redistributes heat that was already in the ocean prior to perturbation. The geographic details of the ζ change in the North Atlantic are diverse across models, but the underlying dynamic change is similar. In contrast, the heat absorbed by the Southern Ocean does not strongly alter the vertically coherent circulation. The Arctic ζ change is dissimilar across models, owing to differences in passive heat uptake and circulation change. Only the Arctic is strongly affected by nonlinear interactions between the three air-sea flux changes, and these are model specific.
The mission for AGU's Journal of Advances in Modeling Earth Systems (JAMES) is to publish original research papers that advance the science that underlies Earth system models and emerges from their use. JAMES′ scope concerns the outer envelope of the Earth system including the atmosphere, oceans, land surface, and cryosphere, and with research focusing on the attendant physics, chemistry, and biology. JAMES considers models as instantiations of scientific theories and concepts, and the journal publishes papers that expand our capabilities to model, understand, and predict the Earth system. In this editorial, we present general principles as well as specific notions that guide the strategy of JAMES′ editors in realizing the journal's mission.
Holmes, Ryan M., Jan D Zika, and Stephen M Griffies, et al., July 2021: The geography of numerical mixing in a suite of global ocean models. Journal of Advances in Modeling Earth Systems, 13(7), DOI:10.1029/2020MS002333. Abstract
Numerical mixing, defined here as the physically spurious tracer diffusion due to the numerical discretization of advection, is known to contribute to biases in ocean models. However, quantifying numerical mixing is nontrivial, with most studies utilizing targeted experiments in idealized settings. Here, we present a water mass transformation-based method for quantifying numerical mixing that can be applied to any conserved variable in general circulation models. Furthermore, the method can be applied within individual fluid columns to provide spatial information. We apply the method to a suite of global ocean model simulations with differing grid spacings and subgrid-scale parameterizations. In all configurations numerical mixing drives diathermal heat transport of comparable magnitude to that associated with explicit parameterizations. Numerical mixing is prominent in the tropical thermocline, where it is sensitive to the vertical diffusivity and resolution. At colder temperatures numerical mixing is sensitive to the presence of explicit neutral diffusion, suggesting that it may act as a proxy for neutral diffusion when it is explicitly absent. Comparison of otherwise equivalent 1/4° and 1/10° configurations with grid-scale dependent horizontal viscosity shows only a modest enhancement in numerical mixing at 1/4°. However, if the lateral viscosity is maintained while resolution is increased then numerical mixing is reduced by almost 35%. This result suggests that the common practice of reducing viscosity in order to maximize permitted variability must be considered carefully. Our results provide a detailed view of numerical mixing in ocean models and pave the way for improvements in parameter choices and numerical methods.
Hsu, Chia-Wei, Jianjun Yin, Stephen M Griffies, and Raphael Dussin, May 2021: A mechanistic analysis of tropical Pacific dynamic sea level in GFDL-OM4 under OMIP-I and OMIP-II forcings. Geoscientific Model Development, 14(5), DOI:10.5194/gmd-14-2471-20212471-2502. 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.
Khatri, Hermant, and Stephen M Griffies, et al., September 2021: Role of mixed-layer instabilities in the seasonal evolution of eddy kinetic energy spectra in a global submesoscale permitting simulation. Geophysical Research Letters, 48(18), DOI:10.1029/2021GL094777. Abstract
A submesoscale-permitting global ocean simulation is used to study the upper ocean turbulence in high kinetic energy (KE) regions. Submesoscale processes peak in winter so that the geostrophic KE spectra tend to be relatively shallow in winter (~k-2) with steeper spectra in summer (~k-3). This transition in KE spectral scaling has two phases. In the first phase (late autumn), KE spectra show the presence of two spectral regimes:~k-3 power-law in mesoscales and ~k-2 power-law in submesoscales. The first phase appears with the onset of mixed-layer instabilities, which convert available potential energy into KE, and this process results in a flattening of KE spectra at submesoscales. However, KE spectra at longer wavelengths follow ~k-3 scaling associated with a forward enstrophy transfer. In the second phase (late winter), KE produced through mixed-layer instabilities is transferred to larger scales, and k-2 power-law also develops in mesoscales.
Lockwood, Joseph W., Carolina O Dufour, Stephen M Griffies, and Michael Winton, April 2021: On the role of the Antarctic Slope Front on the occurrence of the Weddell Sea polynya under climate change. Journal of Climate, 34(7), DOI:10.1175/JCLI-D-20-0069.12529-2548. Abstract
This study investigates the occurrence of the Weddell Sea polynya under an idealized climate change scenario by evaluating simulations from climate models of different ocean resolutions. The GFDL-CM2.6 climate model, with roughly 3.8-km horizontal ocean grid spacing in the high latitudes, forms a Weddell Sea polynya at similar time and duration under idealized climate change forcing as under preindustrial forcing. In contrast, all convective models forming phase 5 of the Coupled Model Intercomparison Project (CMIP5) show either a cessation or a slowdown of Weddell Sea polynya events under climate warming. The representation of the Antarctic Slope Current and related Antarctic Slope Front is found to be key in explaining the differences between the two categories of models, with these features being more realistic in CM2.6 than in CMIP5. In CM2.6, the freshwater input driven by sea ice melt and enhanced runoff found under climate warming largely remains on the shelf region since the slope front restricts the lateral spread of the freshwater. In contrast, for most CMIP5 models, open-ocean stratification is enhanced by freshening since the absence of a slope front allows coastal freshwater anomalies to spread into the open ocean. This enhanced freshening contributes to the slowdown the occurrence of Weddell Sea polynyas. Hence, an improved representation of Weddell Sea shelf processes in current climate models is desirable to increase our ability to predict the fate of the Weddell Sea polynyas under climate change.
McDougall, Trevor J., Paul M Barker, Ryan M Holmes, Rich Pawlowicz, Stephen M Griffies, and Paul J Durack, October 2021: The interpretation of temperature and salinity variables in numerical ocean model output and the calculation of heat fluxes and heat content. Geoscientific Model Development, 14(10), DOI:10.5194/gmd-14-6445-2021. Abstract
The international Thermodynamic Equation of Seawater 2010 (TEOS-10) defined the enthalpy and entropy of seawater, thus enabling the global ocean heat content to be calculated as the volume integral of the product of in situ density, ρ, and potential enthalpy, h0 (with reference sea pressure of 0 dbar). In terms of Conservative Temperature, Θ, ocean heat content is the volume integral of ρc0pΘ, where c0p is a constant “isobaric heat capacity”.
However, many ocean models in the Coupled Model Intercomparison Project Phase 6 (CMIP6) as well as all models that contributed to earlier phases, such as CMIP5, CMIP3, CMIP2, and CMIP1, used EOS-80 (Equation of State – 1980) rather than the updated TEOS-10, so the question arises of how the salinity and temperature variables in these models should be physically interpreted, with a particular focus on comparison to TEOS-10-compliant observations. In this article we address how heat content, surface heat fluxes, and the meridional heat transport are best calculated using output from these models and how these quantities should be compared with those calculated from corresponding observations. We conclude that even though a model uses the EOS-80, which expects potential temperature as its input temperature, the most appropriate interpretation of the model's temperature variable is actually Conservative Temperature. This perhaps unexpected interpretation is needed to ensure that the air–sea heat flux that leaves and arrives in atmosphere and sea ice models is the same as that which arrives in and leaves the ocean model.
We also show that the salinity variable carried by present TEOS-10-based models is Preformed Salinity, while the salinity variable of EOS-80-based models is also proportional to Preformed Salinity. These interpretations of the salinity and temperature variables in ocean models are an update on the comprehensive Griffies et al. (2016) paper that discusses the interpretation of many aspects of coupled Earth system models.
Saenko, Oleg A., Jonathan M Gregory, Stephen M Griffies, Matthew P Couldrey, and Fabio Boeira Dias, March 2021: Contribution of ocean physics and dynamics at different scales to heat uptake in low-resolution AOGCMs. Journal of Climate, 34(6), DOI:10.1175/JCLI-D-20-0652.12017-2035. Abstract
Using an ensemble of atmosphere–ocean general circulation models (AOGCMs) in an idealized climate change experiment, this study quantifies the contributions to ocean heat uptake (OHU) from ocean physical parameterizations and resolved dynamical processes operating at different scales. Analysis of heat budget diagnostics reveals a leading-order global heat balance in the subsurface upper ocean in a steady state between the large-scale circulation warming it and mesoscale processes cooling it, and shows that there are positive contributions from processes on all scales to the subsurface OHU during climate change. There is better agreement among the AOGCMs in the net OHU than in the individual scales/processes contributing to it. In the upper ocean and at high latitudes, OHU is dominated by small-scale diapycnal processes. Below 400 m, OHU is dominated by the superresidual transport, representing large-scale ocean dynamics combined with all parameterized mesoscale and submesoscale eddy effects. Weakening of the AMOC leads to less heat convergence in the subpolar North Atlantic and less heat divergence at lower latitudes, with a small overall effect on the net Atlantic heat content. At low latitudes, the dominance of advective heat redistribution is contrary to the diffusive OHU mechanism assumed by the commonly used upwelling-diffusion model. Using a density water-mass framework, it is found that most of the OHU occurs along isopycnal directions. This feature of OHU is used to accurately reconstruct the global vertical ocean warming profile from the surface heat flux anomalies, supporting advective (rather than diffusive) models of OHU and sea level rise.
Boeira Dias, Fabio, Catia M Domingues, S J Marsland, and Stephen M Griffies, et al., February 2020: On the superposition of mean advective and eddy-induced transports in global ocean heat and salt budgets. Journal of Climate, 33(3), DOI:10.1175/JCLI-D-19-0418.1. Abstract
Ocean thermal expansion is a large contributor to observed sea-level rise which is expected to continue into the future. However, large uncertainties exist in sea-level projections among climate models, partially due to inter-model differences in ocean heat uptake and redistribution of buoyancy. Here, the mechanisms of vertical ocean heat and salt transport are investigated in quasi-steady state model simulations using the Australian Community Climate and Earth System Simulator Ocean Model (ACCESS-OM2). New insights into the net effect of key physical processes are gained within the super-residual transport (SRT) framework. In this framework, vertical tracer transport is dominated by downward fluxes associated with the large-scale ocean circulation and upward fluxes induced by mesoscale eddies, with two distinct physical regimes. In the upper-ocean, where high-latitude watermasses are formed by mixed-layer processes, through cooling or salinification, the SRT counter-acts those processes by transporting heat and salt downward. In contrast, in the ocean interior, the SRT opposes dianeutral diffusion via upward fluxes of heat and salt, with about 60% of the vertical heat transport occurring in the Southern Ocean. Overall, the SRT is largely responsible for removing newly-formed watermasses from the mixed-layer into the ocean interior, where they are eroded by dianeutral diffusion. Unlike the classical advective-diffusive balance, dianeutral diffusion is bottom-intensified above rough bottom topography, allowing an overturning cell to develop in alignment with recent theories. Implications are discussed for understanding the role of vertical tracer transport on the simulation of ocean climate and sea level.
We describe the baseline coupled model configuration and simulation characteristics of GFDL's Earth System Model Version 4.1 (ESM4.1), which builds on component and coupled model developments at GFDL over 2013–2018 for coupled carbon‐chemistry‐climate simulation contributing to the sixth phase of the Coupled Model Intercomparison Project. In contrast with GFDL's CM4.0 development effort that focuses on ocean resolution for physical climate, ESM4.1 focuses on comprehensiveness of Earth system interactions. ESM4.1 features doubled horizontal resolution of both atmosphere (2° to 1°) and ocean (1° to 0.5°) relative to GFDL's previous‐generation coupled ESM2‐carbon and CM3‐chemistry models. ESM4.1 brings together key representational advances in CM4.0 dynamics and physics along with those in aerosols and their precursor emissions, land ecosystem vegetation and canopy competition, and multiday fire; ocean ecological and biogeochemical interactions, comprehensive land‐atmosphere‐ocean cycling of CO2, dust and iron, and interactive ocean‐atmosphere nitrogen cycling are described in detail across this volume of JAMES and presented here in terms of the overall coupling and resulting fidelity. ESM4.1 provides much improved fidelity in CO2 and chemistry over ESM2 and CM3, captures most of CM4.0's baseline simulations characteristics, and notably improves on CM4.0 in (1) Southern Ocean mode and intermediate water ventilation, (2) Southern Ocean aerosols, and (3) reduced spurious ocean heat uptake. ESM4.1 has reduced transient and equilibrium climate sensitivity compared to CM4.0. Fidelity concerns include (1) moderate degradation in sea surface temperature biases, (2) degradation in aerosols in some regions, and (3) strong centennial scale climate modulation by Southern Ocean convection.
This paper provides a primer on the mathematical, physical, and numerical foundations of ocean models that are formulated using finite volume generalized vertical coordinate equations and that use the vertical Lagrangian‐remap method to evolve the ocean state. We consider the mathematical structure of the governing ocean equations in both their strong formulation (partial differential equations) and weak formulation (finite volume integral equations), thus enabling an understanding of their physical content and providing a physical‐mathematical framework to develop numerical algorithms. A connection is made between the Lagrangian‐remap method and the ocean equations as written using finite volume generalized vertical budgets. Thought experiments are offered to exemplify the mechanics of the vertical Lagrangian‐remap method and to compare with other methods used for ocean model algorithms.
Hirschi, J, B Barnier, C Böning, A Biastoch, A T Blaker, A C Coward, S Danilov, Sybren Drijfhout, K Getzlaff, and Stephen M Griffies, et al., April 2020: The Atlantic meridional overturning circulation in high resolution models. Journal of Geophysical Research: Oceans, 125(4), DOI:10.1029/2019JC015522. Abstract
The Atlantic meridional overturning circulation (AMOC) represents the zonally integrated stream function of meridional volume transport in the Atlantic Basin. The AMOC plays an important role in transporting heat meridionally in the climate system. Observations suggest a heat transport by the AMOC of 1.3 PW at 26°N ‐ a latitude which is close to where the Atlantic northward heat transport is thought to reach its maximum. This shapes the climate of the North Atlantic region as we know it today. In recent years there has been significant progress both in our ability to observe the AMOC in nature and to simulate it in numerical models. Most previous modeling investigations of the AMOC and its impact on climate have relied on models with horizontal resolution that does not resolve ocean mesoscale eddies and the dynamics of the Gulf Stream/North Atlantic Current system. As a result of recent increases in computing power, models are now being run that are able to represent mesoscale ocean dynamics and the circulation features that rely on them. The aim of this review is to describe new insights into the AMOC provided by high‐resolution models. Furthermore, we will describe how high‐resolution model simulations can help resolve outstanding challenges in our understanding of the AMOC.
Kiss, A E., Andrew McC Hogg, N Hannah, Fabio Boeira Dias, G Brassington, Matthew A Chamberlain, C Chapman, Peter Dobrohotoff, Catia M Domingues, E R Duran, Matthew H England, R Fiedler, Stephen M Griffies, A Heerdegen, P Heil, Ryan M Holmes, A Klocker, S J Marsland, Adele K Morrison, J Munroe, P Oke, M Nikurashin, G S Pilo, O Richet, A Savita, P Spence, K D Stewart, and Marshall L Ward, et al., February 2020: ACCESS-OM2: A Global Ocean-Sea Ice Model at Three Resolutions. Geoscientific Model Development, 13(2), DOI:10.5194/gmd-13-401-2020. Abstract
We introduce a new version of the ocean-sea ice implementation of the Australian Community Climate and Earth System Simulator, ACCESS-OM2. The model has been developed with the aim of being aligned as closely as possible with the fully coupled (atmosphere-land-ocean-sea ice) ACCESS-CM2. Importantly, the model is available at three different horizontal resolutions: a coarse resolution (nominally 1° horizontal grid spacing), an eddy-permitting resolution (nominally 0.25°) and an eddy-rich resolution (0.1° with 75 vertical levels), where the eddy-rich model is designed to be incorporated into the Bluelink operational ocean prediction and reanalysis system. The different resolutions have been developed simultaneously, both to allow testing at lower resolutions and to permit comparison across resolutions. In this manuscript, the model is introduced and the individual components are documented. The model performance is evaluated across the three different resolutions, highlighting the relative advantages and disadvantages of running ocean-sea ice models at higher resolution. We find that higher resolution is an advantage in resolving flow through small straits, the structure of western boundary currents and the abyssal overturning cell, but that there is scope for improvements in sub-grid scale parameterisations at the highest resolution.
Rahaman, H, U Srinivasu, Swapna Panickal, J V Durgadoo, and Stephen M Griffies, et al., January 2020: An assessment of the Indian Ocean mean state and seasonal cycle in a suite of interannual CORE-II simulations. Ocean Modelling, 145, 101503, DOI:10.1016/j.ocemod.2019.101503. Abstract
We present an analysis of annual and seasonal mean characteristics of the Indian Ocean circulation and water masses from 16 global ocean-sea-ice model simulations that follow the Coordinated Ocean-ice Reference Experiments (CORE) interannual protocol (CORE-II). All simulations show a similar large-scale tropical current system, but with differences in the Equatorial Undercurrent. Most CORE-II models simulate the structure of the Cross Equatorial Cell (CEC) in the Indian Ocean. We uncover a previously unidentified secondary pathway of northward cross-equatorial transport along 75 degrees E, thus complementing the pathway near the Somali Coast. This secondary pathway is most prominent in the models which represent topography realistically, thus suggesting a need for realistic bathymetry in climate models. When probing the water mass structure in the upper ocean, we find that the salinity profiles are closer to observations in geopotential (level) models than in isopycnal models. More generally, we find that biases are model dependent, thus suggesting a grouping into model lineage, formulation of the surface boundary, vertical coordinate and surface salinity restoring. Refinement in model horizontal resolution (one degree versus 1/4 degree) does not significantly improve simulations, though there are some marginal improvements in the salinity and barrier layer results. The results in turn suggest that a focus on improving physical parameterizations (e.g. boundary layer processes) may offer more near-term advances in Indian Ocean simulations than refined grid resolution.
We present a neutral diffusion operator appropriate for an ocean model making use of general vertical coordinates. The diffusion scheme uses polynomial reconstructions in the vertical, along with a horizontally local but vertically nonlocal stencil for estimates of tracer fluxes. These fluxes are calculated on a vertical grid that is the superset of model columns in a neutral density space. Using flux-limiters, the algorithm dissipates tracer extrema locally, and no new extrema are created. A demonstration using a linear equation of state in an idealized configuration shows that the algorithm is perfectly neutral. When using the nonlinear TEOS-10 equation of state with a constant reference pressure, the algorithm compares nearly exactly to a case discretized onto isopycnal surfaces and using along-layer diffusion. The algorithm's cost is comparable to that of tracer advection and can be readily implemented into ocean general circulation models.
Tsujino, Hiroyuki, Shogo Urakawa, Stephen M Griffies, Gokhan Danabasoglu, Alistair Adcroft, A E Amaral, T Arsouze, M Bentsen, Raffaele Bernardello, C Böning, A Bozec, Eric P Chassignet, S Danilov, and Raphael Dussin, et al., August 2020: Evaluation of global ocean–sea-ice model simulations based on the experimental protocols of the Ocean Model Intercomparison Project phase 2 (OMIP-2). Geoscientific Model Development, 13(8), DOI:10.5194/gmd-13-3643-2020. Abstract
We present a new framework for global ocean–sea-ice model simulations based on phase 2 of the Ocean Model Intercomparison Project (OMIP-2), making use of the surface dataset based on the Japanese 55-year atmospheric reanalysis for driving ocean–sea-ice models (JRA55-do). We motivate the use of OMIP-2 over the framework for the first phase of OMIP (OMIP-1), previously referred to as the Coordinated Ocean–ice Reference Experiments (COREs), via the evaluation of OMIP-1 and OMIP-2 simulations from 11 state-of-the-science global ocean–sea-ice models. In the present evaluation, multi-model ensemble means and spreads are calculated separately for the OMIP-1 and OMIP-2 simulations and overall performance is assessed considering metrics commonly used by ocean modelers. Both OMIP-1 and OMIP-2 multi-model ensemble ranges capture observations in more than 80 % of the time and region for most metrics, with the multi-model ensemble spread greatly exceeding the difference between the means of the two datasets. Many features, including some climatologically relevant ocean circulation indices, are very similar between OMIP-1 and OMIP-2 simulations, and yet we could also identify key qualitative improvements in transitioning from OMIP-1 to OMIP-2. For example, the sea surface temperatures of the OMIP-2 simulations reproduce the observed global warming during the 1980s and 1990s, as well as the warming slowdown in the 2000s and the more recent accelerated warming, which were absent in OMIP-1, noting that the last feature is part of the design of OMIP-2 because OMIP-1 forcing stopped in 2009. A negative bias in the sea-ice concentration in summer of both hemispheres in OMIP-1 is significantly reduced in OMIP-2. The overall reproducibility of both seasonal and interannual variations in sea surface temperature and sea surface height (dynamic sea level) is improved in OMIP-2. These improvements represent a new capability of the OMIP-2 framework for evaluating process-level responses using simulation results. Regarding the sensitivity of individual models to the change in forcing, the models show well-ordered responses for the metrics that are directly forced, while they show less organized responses for those that require complex model adjustments. Many of the remaining common model biases may be attributed either to errors in representing important processes in ocean–sea-ice models, some of which are expected to be reduced by using finer horizontal and/or vertical resolutions, or to shared biases and limitations in the atmospheric forcing. In particular, further efforts are warranted to resolve remaining issues in OMIP-2 such as the warm bias in the upper layer, the mismatch between the observed and simulated variability of heat content and thermosteric sea level before 1990s, and the erroneous representation of deep and bottom water formations and circulations. We suggest that such problems can be resolved through collaboration between those developing models (including parameterizations) and forcing datasets. Overall, the present assessment justifies our recommendation that future model development and analysis studies use the OMIP-2 framework.
Storm surge and coastal flooding caused by tropical cyclones (hurricanes) and extratropical cyclones (nor'easters) pose a threat to communities along the Atlantic coast of the United States. Climate change and sea level rise are altering the statistics of these extreme events in a rather complex fashion. Here we use a fully-coupled global weather/climate modeling system (GFDL CM4) to study characteristics of extreme daily sea level (ESL) along the US Atlantic coast and their response to global warming. We find that under natural weather processes, the Gulf of Mexico coast is most vulnerable to storm surge and related ESL. New Orleans is a striking hotspot with the highest surge efficiency in response to storm winds. Under a 1% per year atmospheric CO2 increase on centennial time scales, the anthropogenic signal in ESL is robust along the US East Coast. It can emerge from the background variability as soon as in twenty years, or even before global sea level rise is taken into account. The regional dynamic sea level rise induced by the weakening of the Atlantic meridional overturning circulation facilitates this early emergence, especially during wintertime coastal flooding associated with nor’easters. Along the Gulf Coast, ESL is sensitive to the modification of hurricane characteristics under the CO2 forcing.
We document the configuration and emergent simulation features from the Geophysical Fluid Dynamics Laboratory (GFDL) OM4.0 ocean/sea‐ice model. OM4 serves as the ocean/sea‐ice component for the GFDL climate and Earth system models. It is also used for climate science research and is contributing to the Coupled Model Intercomparison Project version 6 Ocean Model Intercomparison Project (CMIP6/OMIP). The ocean component of OM4 uses version 6 of the Modular Ocean Model (MOM6) and the sea‐ice component uses version 2 of the Sea Ice Simulator (SIS2), which have identical horizontal grid layouts (Arakawa C‐grid). We follow the Coordinated Ocean‐sea ice Reference Experiments (CORE) protocol to assess simulation quality across a broad suite of climate relevant features. We present results from two versions differing by horizontal grid spacing and physical parameterizations: OM4p5 has nominal 0.5° spacing and includes mesoscale eddy parameterizations and OM4p25 has nominal 0.25° spacing with no mesoscale eddy parameterization.
MOM6 makes use of a vertical Lagrangian‐remap algorithm that enables general vertical coordinates. We show that use of a hybrid depth‐isopycnal coordinate reduces the mid‐depth ocean warming drift commonly found in pure z* vertical coordinate ocean models. To test the need for the mesoscale eddy parameterization used in OM4p5, we examine the results from a simulation that removes the eddy parameterization. The water mass structure and model drift are physically degraded relative to OM4p5, thus supporting the key role for a mesoscale closure at this resolution.
We revisit the challenges and prospects for ocean circulation models following Griffies et al. (2010). Over the past decade, ocean circulation models evolved through improved understanding, numerics, spatial discretization, grid configurations, parameterizations, data assimilation, environmental monitoring, and process-level observations and modeling. Important large scale applications over the last decade are simulations of the Southern Ocean, the Meridional Overturning Circulation and its variability, and regional sea level change. Submesoscale variability is now routinely resolved in process models and permitted in a few global models, and submesoscale effects are parameterized in most global models. The scales where nonhydrostatic effects become important are beginning to be resolved in regional and process models. Coupling to sea ice, ice shelves, and high-resolution atmospheric models has stimulated new ideas and driven improvements in numerics. Observations have provided insight into turbulence and mixing around the globe and its consequences are assessed through perturbed physics models. Relatedly, parameterizations of the mixing and overturning processes in boundary layers and the ocean interior have improved. New diagnostics being used for evaluating models alongside present and novel observations are briefly referenced. The overall goal is summarizing new developments in ocean modeling, including: how new and existing observations can be used, what modeling challenges remain, and how simulations can be used to support observations.
Gregory, Jonathan M., and Stephen M Griffies, et al., November 2019: Concepts and terminology for sea level--mean, variability and change, both local and global. Surveys in Geophysics, 40(6), DOI:10.1007/s10712-019-09525-z. Abstract
Changes in sea level lead to some of the most severe impacts of anthropogenic climate change. Consequently, they are a subject of great interest in both scientific research and public policy. This paper defines concepts and terminology associated with sea level and sea-level changes in order to facilitate progress in sea-level science, in which communication is sometimes hindered by inconsistent and unclear language. We identify key terms and clarify their physical and mathematical meanings, make links between concepts and across disciplines, draw distinctions where there is ambiguity, and propose new terminology where it is lacking or where existing terminology is confusing. We include formulae and diagrams to support the definitions.
The water mass transformation (WMT) framework weaves together circulation, thermodynamics, and biogeochemistry into a description of the ocean that complements traditional Eulerian and Lagrangian methods. In so doing, a WMT analysis renders novel insights and predictive capabilities for studies of ocean physics and biogeochemistry. In this review, we describe fundamentals of the WMT framework and illustrate its practical analysis capabilities. We show how it provides a robust methodology to characterize and quantify the impact of physical processes on buoyancy and other thermodynamic fields. We also detail how to extend WMT to insightful analysis of biogeochemical cycles.
Groeskamp, S, Paul M Barker, Trevor J McDougall, Ryan Abernathey, and Stephen M Griffies, August 2019: VENM An algorithm to accurately calculate neutral slopes and gradients. Journal of Advances in Modeling Earth Systems, 11(7), DOI:10.1029/2019MS001613. Abstract
Mesoscale eddies stir along the neutral plane, and the resulting neutral diffusion is a fundamental aspect of subgrid scale tracer transport in ocean models. Calculating neutral diffusion traditionally involves calculating neutral slopes and three‐dimensional tracer gradients. The calculation of the neutral slope traditionally occurs by computing the ratio of the horizontal to vertical locally referenced potential density derivative. However, this approach is problematic in regions of weak vertical stratification, prompting the use of a variety of ad hoc regularization methods that can lead to rather nonphysical dependencies for the resulting neutral tracer gradients. Here we use a VErtical Non‐local Method “VENM”, a search algorithm that requires no ad hoc regularization and significantly improves the numerical accuracy of calculating neutral slopes, neutral tracer gradients and associated neutral diffusive fluxes. We compare and contrast VENM against a more traditional method, using an independent objective neutrality condition combined with estimates of spurious diffusion, heat transport and water mass transformation rates. VENM is more accurate, both physically and numerically and should form the basis for future efforts involving neutral diffusion calculations from observations and possibly numerical model simulations.
We describe GFDL's CM4.0 physical climate model, with emphasis on those aspects that may be of particular importance to users of this model and its simulations. The model is built with the AM4.0/LM4.0 atmosphere/land model and OM4.0 ocean model. Topics include the rationale for key choices made in the model formulation, the stability as well as drift of the pre‐industrial control simulation, and comparison of key aspects of the historical simulations with observations from recent decades. Notable achievements include the relatively small biases in seasonal spatial patterns of top‐of‐atmosphere fluxes, surface temperature, and precipitation; reduced double Intertropical Convergence Zone bias; dramatically improved representation of ocean boundary currents; a high quality simulation of climatological Arctic sea ice extent and its recent decline; and excellent simulation of the El Niño‐Southern Oscillation spectrum and structure. Areas of concern include inadequate deep convection in the Nordic Seas; an inaccurate Antarctic sea ice simulation; precipitation and wind composites still affected by the equatorial cold tongue bias; muted variability in the Atlantic Meridional Overturning Circulation; strong 100 year quasi‐periodicity in Southern Ocean ventilation; and a lack of historical warming before 1990 and too rapid warming thereafter due to high climate sensitivity and strong aerosol forcing, in contrast to the observational record. Overall, CM4.0 scores very well in its fidelity against observations compared to the Coupled Model Intercomparison Project Phase 5 generation in terms of both mean state and modes of variability and should prove a valuable new addition for analysis across a broad array of applications.
Six recent Langmuir turbulence parameterization schemes and five traditional schemes are implemented in a common single column modeling framework and consistently compared. These schemes are tested in scenarios versus matched large eddy simulations (LES), across the globe with realistic forcing (JRA55‐do, WAVEWATCH‐III simulated waves) and initial conditions (Argo), and under realistic conditions as observed at ocean moorings. Traditional non‐Langmuir schemes systematically under‐predict LES vertical mixing under weak convective forcing, while Langmuir schemes vary in accuracy. Under global, realistic forcing Langmuir schemes produce 6% (‐1% to 14% for 90% confidence) or 5.2 m (‐0.2 m to 17.4 m for 90% confidence) deeper monthly mean mixed layer depths (MLD) than their non‐Langmuir counterparts, with the greatest differences in extratropical regions, especially the Southern Ocean in austral summer. Discrepancies among Langmuir schemes are large (15% in MLD standard deviation over the mean): largest under wave‐driven turbulence with stabilizing buoyancy forcing, next largest under strongly wave‐driven conditions with weak buoyancy forcing, and agreeing during strong convective forcing. Non‐Langmuir schemes disagree with each other to a lesser extent, with a similar ordering. Langmuir discrepancies obscure a cross‐scheme estimate of the Langmuir effect magnitude under realistic forcing, highlighting limited understanding and numerical deficiencies. Maps of the regions and seasons where the greatest discrepancies occur are provided to guide further studies and observations.
Naveira Garabato, Alberto C., Eleanor E Frajka-Williams, Carl P Spingys, Sonya Legg, Kurt L Polzin, A Forryan, E Povl Abrahamsen, Christian E Buckingham, and Stephen M Griffies, July 2019: Rapid mixing and exchange of deep-ocean waters in an abyssal boundary current. Proceedings of the National Academy of Sciences, 116(27), DOI:10.1073/pnas.1904087116. Abstract
The overturning circulation of the global ocean is critically shaped by deep-ocean mixing, which transforms cold waters sinking at high latitudes into warmer, shallower waters. The effectiveness of mixing in driving this transformation is jointly set by two factors: the intensity of turbulence near topography and the rate at which well-mixed boundary waters are exchanged with the stratified ocean interior. Here, we use innovative observations of a major branch of the overturning circulation—an abyssal boundary current in the Southern Ocean—to identify a previously undocumented mixing mechanism, by which deep-ocean waters are efficiently laundered through intensified near-boundary turbulence and boundary–interior exchange. The linchpin of the mechanism is the generation of submesoscale dynamical instabilities by the flow of deep-ocean waters along a steep topographic boundary. As the conditions conducive to this mode of mixing are common to many abyssal boundary currents, our findings highlight an imperative for its representation in models of oceanic overturning.
Nurser, A J G., and Stephen M Griffies, September 2019: Relating the diffusive salt flux just below the ocean surface to boundary freshwater and salt fluxes. Journal of Physical Oceanography, 49(9), DOI:10.1175/JPO-D-19-0037.1. Abstract
We detail the physical means whereby boundary transfers of freshwater and salt induce diffusive fluxes of salinity. Our considerations focus on the kinematic balance between the diffusive fluxes of salt and freshwater, with this balance imposed by mass conservation for an element of seawater. The flux balance leads to a specific form for the diffusive salt flux immediately below the ocean surface and, in the Boussinesq approximation, to a specific form for the salinity flux. This note clarifies conceptual and formulational ambiguities in the literature concerning the surface boundary condition for the salinity equation and for the contribution of freshwater to the buoyancy budget.
Ponte, R M., M Carson, M Cirano, Catia M Domingues, S Jevrejeva, M Marcos, G Mitchum, R S W van der Wal, P L Woodworth, M Ablain, F Ardhuin, V Ballu, M Becker, J Benveniste, F Birol, E Bradshaw, A Cazenave, P De Mey-Frémaux, F Durand, Tal Ezer, L-L Fu, I Fukumori, K Gordon, M Gravelle, and Stephen M Griffies, et al., July 2019: Towards Comprehensive Observing and Modeling Systems for Monitoring and Predicting Regional to Coastal Sea Level. Frontiers in Marine Science, 6, DOI:10.3389/fmars.2019.00437. Abstract
A major challenge for managing impacts and implementing effective mitigation measures and adaptation strategies for coastal zones affected by future sea level (SL) rise is our limited capacity to predict SL change at the coast on relevant spatial and temporal scales. Predicting coastal SL requires the ability to monitor and simulate a multitude of physical processes affecting SL, from local effects of wind waves and river runoff to remote influences of the large-scale ocean circulation on the coast. Here we assess our current understanding of the causes of coastal SL variability on monthly to multi-decadal timescales, including geodetic, oceanographic and atmospheric aspects of the problem, and review available observing systems informing on coastal SL. We also review the ability of existing models and data assimilation systems to estimate coastal SL variations and of atmosphere-ocean global coupled models and related regional downscaling efforts to project future SL changes. We discuss (1) observational gaps and uncertainties, and priorities for the development of an optimal and integrated coastal SL observing system, (2) strategies for advancing model capabilities in forecasting short-term processes and projecting long-term changes affecting coastal SL, and (3) possible future developments of sea level services enabling better connection of scientists and user communities and facilitating assessment and decision making for adaptation to future coastal SL change.
Randall, David A., C M Bitz, Gokhan Danabasoglu, A S Denning, P R Gent, Andrew Gettelman, and Stephen M Griffies, et al., July 2019: 100 Years of Earth System Model Development In A Century of Progress in Atmospheric and Related Sciences: Celebrating the American Meteorological Society Centennial, Boston, MA, Meteorological Monographs, American Meteorological Society, 59, DOI:10.1175/AMSMONOGRAPHS-D-18-0018.112.1-12.66. Abstract
Today’s global Earth system models began as simple regional models of tropospheric weather systems. Over the past century, the physical realism of the models has steadily increased, while the scope of the models has broadened to include the global troposphere and stratosphere, the ocean, the vegetated land surface, and terrestrial ice sheets. This chapter gives an approximately chronological account of the many and profound conceptual and technological advances that made today’s models possible. For brevity, we omit any discussion of the roles of chemistry and biogeochemistry, and terrestrial ice sheets.
Stammer, Detlef, A Bracco, Krishna M Achuta Rao, L Beal, Nathaniel L Bindoff, P Braconnot, Wenju Cai, D Chen, Matthew Collins, Gokhan Danabasoglu, Boris Dewitte, Riccardo Farneti, Baylor Fox-Kemper, John C Fyfe, and Stephen M Griffies, et al., July 2019: Ocean climate observing requirements in support of Climate Research and Climate Information. Frontiers in Marine Science, 6, DOI:10.3389/fmars.2019.00444. Abstract
Natural variability and change of the Earth’s climate have significant global societal impacts. With its large heat and carbon capacity and relatively slow dynamics, the ocean plays an integral role in climate, and provides an important source of predictability at seasonal and longer timescales. In addition, the ocean provides the slowly evolving lower boundary to the atmosphere, both driving and modifying atmospheric weather and climate. Understanding and monitoring ocean climate variability and change, to constrain and initialize models as well as identify model biases for improved climate hindcasting and prediction requires a more scale-sensitive, long-term observing system. A climate observing system has requirements that significantly differ from, and sometimes are orthogonal to, those of other applications. In general terms, they can be summarized by the simultaneous need for both large spatial and long temporal coverage, and by the accuracy and stability required for detecting the local climate signals.
This paper reviews the requirements of a climate observing system in terms of space and time scales, and revisits the question of which parameters such a system should encompass to meet future strategic goals of the World Climate Research Program (WCRP), with emphasis on ocean and sea-ice covered areas. It considers global as well as regional aspects that should be accounted for in designing observing systems in individual basins. Furthermore, the paper discusses which data-driven products are required to meet WCRP research and modeling needs, and ways to obtain them through data synthesis and assimilation approaches. Finally, it addresses the need for scientific capacity building and international collaboration in support of the collection of high-quality measurements over the large spatial scales and long time-scales required for climate research, bridging the scientific rational to the required resources for implementation.
Ocean surface winds determine energy, material and momentum fluxes through the air-sea interface. Accounting for wind variability in time and space is thus essential to reliably analyze and simulate ocean circulation and the dynamics of marine ecosystems. Here, we present an assessment of surface winds from three widely used atmospheric reanalysis products (NCEP/NCAR, ERA-Interim and JRA-55) and their corresponding ocean forcing data sets (CORE v2.1, DFS v5.2 and JRA55-do), which include corrections for use in ocean simulations. We compared wind patterns most relevant to ocean circulation (surface wind stress, its curl and estimates of induced vertical upwelling velocity) across global and regional scales, with added emphasis on the main Eastern Boundary Upwelling Ecosystems (EBUEs). All products provided consistent large-scale patterns in surface winds and wind stress, although agreement was reduced for indices involving the calculation of spatial derivatives, like wind stress curl and Ekman pumping. Fidelity with respect to a reference reanalysis based on blended satellite and buoy observations (CCMP v2.0) improved in more recent, higher resolution products like JRA-55 and ERA-Interim. Adjustments applied when deriving ocean forcing data sets from atmospheric reanalysis robustly improved wind speed and wind stress vectors, but degraded wind stress curl (and implied Ekman upwelling) in two of the three ocean forcing products considered (DFS v5.2 and CORE v2.1).
At regional scales, we found significant inconsistencies in equatorial and polar regions, as well as in coastal areas. In EBUEs, upwelling favorable winds were weaker in atmospheric reanalysis products and ocean forcing data sets than estimates based on CCMP v2.0 and QuikSCAT. All reanalysis products featured lower amplitude seasonal cycles and contrasting patterns of low frequency variability within each EBUE, including the presence of sudden changes in mean upwelling only for some products.
Taken together, our results highlight the importance of incorporating uncertainties in wind forcing into ocean simulation experiments and retrospective analysis, and of correcting reanalysis products for ocean forcing data sets. Despite the continued improvement in the quality of wind data sets, prevailing limitations in reanalysis models demonstrate the need to confirm global products against regional measurements whenever possible and improve correction strategies across multiple ocean-relevant wind properties.
Meltwater from the Antarctic Ice Sheet is projected to cause up to one metre of sea-level rise by 2100 under the highest greenhouse gas concentration trajectory (RCP8.5) considered by the Intergovernmental Panel on Climate Change (IPCC). However, the effects of meltwater from the ice sheets and ice shelves of Antarctica are not included in the widely used CMIP5 climate models, which introduces bias into IPCC climate projections. Here we assess a large ensemble simulation of the CMIP5 model ‘GFDL ESM2M’ that accounts for RCP8.5-projected Antarctic Ice Sheet meltwater. We find that, relative to the standard RCP8.5 scenario, accounting for meltwater delays the exceedance of the maximum global-mean atmospheric warming targets of 1.5 and 2 degrees Celsius by more than a decade, enhances drying of the Southern Hemisphere and reduces drying of the Northern Hemisphere, increases the formation of Antarctic sea ice (consistent with recent observations of increasing Antarctic sea-ice area) and warms the subsurface ocean around the Antarctic coast. Moreover, the meltwater-induced subsurface ocean warming could lead to further ice-sheet and ice-shelf melting through a positive feedback mechanism, highlighting the importance of including meltwater effects in simulations of future climate.
In this paper we study upwelling pathways and timescales of Circumpolar Deep Water (CDW) in a hierarchy of models using a Lagrangian particle tracking method. Lagrangian timescales of CDW upwelling decrease from 87 years to 31 years to 17 years as the ocean resolution is refined from 1° to 0.25° to 0.1°. We attribute some of the differences in timescale to the strength of the eddy fields, as demonstrated by temporally degrading high resolution model velocity fields. Consistent with the timescale dependence, we find that an average Lagrangian particle completes 3.2 circumpolar loops in the 1° model in comparison to 0.9 loops in the 0.1° model. These differences suggest that advective timescales and thus inter-basin merging of upwelling CDW may be overestimated by coarse resolution models, potentially affecting the skill of centennial scale climate change projections.
Oceanic heat uptake (OHU) is a significant source of uncertainty in both the transient and equilibrium responses to increasing the planetary radiative forcing. OHU differs among climate models and is related in part to their representation of vertical and lateral mixing. This study examines the role of ocean model formulation – specifically the choice of vertical coordinate and strength of background diapycnal diffusivity (Kd) – in the millennial-scale near-equilibrium climate response to a quadrupling of atmospheric CO2. Using two fully-coupled Earth System Models (ESMs) with nearly identical atmosphere, land, sea ice, and biogeochemical components, it is possible to independently configure their ocean model components with different formulations and produce similar near-equilibrium climate responses. The SST responses are similar between the two models (r2 = 0.75, global average ∼ 4.3 °C) despite their initial pre-industrial climate mean states differing by 0.4 °C globally. The surface and interior responses of temperature and salinity are also similar between the two models. However, the Atlantic Meridional Overturning Circulation (AMOC) responses are different between the two models, and the associated differences in ventilation and deep water formation have an impact on the accumulation of dissolved inorganic carbon in the ocean interior. A parameter sensitivity analysis demonstrates that increasing the amount of Kd produces very different near-equilibrium climate responses within a given model. These results suggest that the impact of the ocean vertical coordinate on the climate response is small relative to the representation of sub-gridscale mixing.
O'Rourke, Amanda K., Brian K Arbic, and Stephen M Griffies, March 2018: Frequency-domain analysis of atmospherically forced versus intrinsic ocean surface kinetic energy variability in GFDL’s CM2-O model hierarchy. Journal of Climate, 31(5), DOI:10.1175/JCLI-D-17-0024.1. Abstract
Low-frequency variability at the ocean surface can be excited by both atmospheric forcing, such as in exchanges of heat and momentum, and by the intrinsic nonlinear transfer of energy between mesoscale ocean eddies. Recent studies have shown that nonlinear eddy interactions can excite an energy transfer from high to low frequencies analogous to the transfer of energy from high to low wavenumbers (small to large spatial scales) in quasi-twodimensional turbulence. As the spatial inverse cascade is driven by oceanic eddies, the process of energy exchange across frequencies may be sensitive to ocean model resolution. Here a cross-spectrum diagnostic is applied to GFDLs CM2-O hierarchy of fully coupled ocean-atmosphere models to address the transfer of ocean surface kinetic energy between high and low frequencies. The cross-spectral diagnostic allows for a comparison of the relative contributions of coupled atmospheric forcing through wind stress and the intrinsic advection to low-frequency ocean surface kinetic energy. Diagnostics of energy flux and transfer within the frequency domain are compared between three coupled models with horizontal ocean resolutions of 1, 1/4th, and 1/10th degree to address the importance of resolving eddies in the driving of energy to low frequencies in coupled models.
The Pacific equatorial cold tongue plays a leading role in Earth’s strongest and most predictable climate signals. To illuminate the processes governing cold tongue temperatures, the upper-ocean heat budget is explored using the GFDL FLOR coupled GCM. Starting from the exact temperature budget for layers of time-varying thickness, the layer temperature tendency terms are studied using hourly-, daily-, and monthly-mean output from a 30-year simulation driven by present-day radiative forcings. The budget is then applied to (1) a surface mixed layer whose temperature is highly correlated with SST, in which the air-sea heat flux is balanced mainly by downward diffusion of heat across the layer base; and (2) a thicker advective layer that subsumes most of the vertical mixing, in which the air-sea heat flux is balanced mainly by monthly-scale advection. The surface warming from shortwave fluxes and submonthly meridional advection, and the subsurface cooling from monthly vertical advection, are both shown to be essential to maintain the cold tongue thermal stratification against the destratifying effects of vertical mixing. Although layer undulations strongly mediate the tendency terms on diurnal-to-interannual scales, the 30-year-mean tendencies are found to be well summarized by analogous budgets developed for stationary but spatially-varying layers. The results are used to derive practical simplifications of the exact budget, to support the analyses in Part II of this paper and to facilitate broader application of heat budget analyses when evaluating and comparing climate simulations.
The heat budget of the Pacific equatorial cold tongue (ECT) is explored using the GFDL FLOR coupled GCM and ocean reanalyses, leveraging the two-layer framework developed in Part I. Despite FLOR’s relatively weak meridional stirring by tropical instability waves (TIWs), the model maintains a reasonable SST and thermocline depth in the ECT via two compensating biases: (1) enhanced monthly-scale vertical advective cooling below the surface mixed layer (SML), due to overly cyclonic off-equatorial wind stress which acts to cool the equatorial source waters; and (2) an excessive SST contrast between the ECT and off-equator, which boosts the equatorward heat transport by TIWs. FLOR’s strong advective cooling at the SML base is compensated by strong downward diffusion of heat out of the SML, which then allows FLOR’s ECT to take up a realistic heat flux from the atmosphere. Correcting FLOR’s climatological SST and wind stress biases via flux adjustment (FA) leads to weaker deep advective cooling of the ECT, which then erodes the upper-ocean thermal stratification, enhances vertical mixing, and excessively deepens the thermocline. FA does strengthen FLOR’s meridional shear of the zonal currents in the east Pacific, but this does not amplify the simulated TIWs nor their equatorward heat transport, likely due to FLOR’s coarse zonal ocean resolution. The analysis suggests that to advance coupled simulations of the ECT, improved winds and surface heat fluxes must go hand in hand with improved subseasonal and parameterized ocean processes. Implications for model development and the tropical Pacific observing system are discussed.
Roberts, Malcolm J., Pier Luigi Vidale, Catherine A Senior, Helene T Hewitt, C Bates, S Berthou, Ping Chang, J H Christensen, S Danilov, M-E Demory, and Stephen M Griffies, et al., November 2018: The benefits of global high-resolution for climate simulation: process-understanding and the enabling of stakeholder decisions at the regional scale. Bulletin of the American Meteorological Society, 99(11), DOI:10.1175/BAMS-D-15-00320.1. Abstract
A perspective on current and future capabilities in global high-resolution climate simulation for assessing climate risks over next few decades, including advances in process representation and analysis, justifying the emergence of dedicated, coordinated experimental protocols.
The timescales of the Paris Climate Agreement indicate urgent action is required on climate policies over the next few decades, in order to avoid the worst risks posed by climate change. On these relatively short timescales the combined effect of climate variability and change are both key drivers of extreme events, with decadal timescales also important for infrastructure planning. Hence, in order to assess climate risk on such timescales, we require climate models to be able to represent key aspects of both internally driven climate variability, as well as the response to changing forcings.
In this paper we argue that we now have the modelling capability to address these requirements - specifically with global models having horizontal resolutions considerably enhanced from those typically used in previous IPCC and CMIP exercises. The improved representation of weather and climate processes in such models underpins our enhanced confidence in predictions and projections, as well as providing improved forcing to regional models, which are better able to represent local-scale extremes (such as convective precipitation). We choose the global water cycle as an illustrative example, because it is governed by a chain of processes for which there is growing evidence of the benefits of higher resolution. At the same time it comprises key processes involved in many of the expected future climate extremes (e.g. flooding, drought, tropical and mid-latitude storms).
Stammer, Detlef, A Bracco, P Braconnot, G P Brasseur, Stephen M Griffies, and Ed Hawkins, November 2018: Science Directions in a Post COP21 World of Transient Climate Change: Enabling Regional to Local Predictions in Support of Reliable Climate Information. Earth's Future, 6(11), DOI:10.1029/2018EF000979. Abstract
During recent decades, through theoretical considerations and analyses of observations and model simulations, the scientific community has fundamentally advanced our understanding of the coupled climate system, thereby establishing that humans affect the Earth's climate. Resulting from this remarkable accomplishment, the COP21 agreement marks a historic turning point for climate research by calling for actionable regional climate change information on time scales from seasonal to centuries for the benefit of humanity, as well as living and nonliving elements of the Earth environment. Out of the underlying United National Framework Convention on climate Change process, improving seamless regional climate forecast capabilities emerges as a key challenge for the international research community. Addressing it requires a multiscale approach to climate predictions. Here we offer a vision that emphasizes enhanced scientific understanding of regional to local climate processes as the foundation for progress. The scientific challenge is extreme due to the rich complexity of interactions and feedbacks between regional and global processes, each of which affects the global climate trajectory. To gain the necessary scientific insight and to turn it into actionable climate information require technical development, international coordination, and a close interaction between the science and stakeholder communities.
We identify Lagrangian coherent vortices in a global mesoscale eddy-permitting ocean model using the rotation-based method of Haller et al. (2016). We present an analysis of the acute sensitivity of the identification results to varying the method’s free parameters, and develop physically justified parameter choices that allow for systematic vortex identification. In contrast to prior vortex studies, we probe the broad spectrum of coherency in the ocean by determining free parameter choices that partition the spectrum into distinct coherency classes, allowing for the identification of strictly coherent, moderately coherent, and leaky vortices. Our tuning methodology is grounded in a combination of sensitivity analysis, convergence tests, and consideration of the ocean model’s physics. To aid in this process, we introduce the Coherency Index, a novel Lagrangian diagnostic for mathematically quantifying the degree of material coherency of a Lagrangian vortex. We aim for this manuscript and the accompanying open-access code to serve as a manual and toolset for the oceanographer interested in harnessing a rigorous Lagrangian method to uncover coherent structures in ocean models and observations.
Tsujino, Hiroyuki, Shogo Urakawa, Hideyuki Nakano, J Small, W M Kim, Stephen G Yeager, Gokhan Danabasoglu, Tatsuo Suzuki, J L Bamber, M Bentsen, C Böning, A Bozec, Eric P Chassignet, Enrique N Curchitser, Fabio Boeira Dias, Paul J Durack, and Stephen M Griffies, et al., October 2018: JRA-55 based surface dataset for driving ocean–sea-ice models (JRA55-do). Ocean Modelling, 130, DOI:10.1016/j.ocemod.2018.07.002. Abstract
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.
Van Roekel, L, Alistair Adcroft, Gokhan Danabasoglu, Stephen M Griffies, B Kauffman, William G Large, Michael Levy, and Brandon G Reichl, et al., November 2018: The KPP boundary layer scheme for the ocean: revisiting its formulation and benchmarking one‐dimensional simulations relative to LES. Journal of Advances in Modeling Earth Systems, 10(11), DOI:10.1029/2018MS001336. Abstract
We evaluate the Community ocean Vertical Mixing (CVMix) project version of the K‐profile parameterization (KPP) for modeling upper ocean turbulent mixing. For this purpose, one‐dimensional KPP simulations are compared across a suite of oceanographically relevant regimes against horizontally averaged large eddy simulations (LES). We find the standard configuration of KPP consistent with LES across many forcing regimes, supporting its physical basis. Our evaluation also motivates recommendations for KPP “best practices” within ocean circulation models, and identifies areas where further research is warranted.
The original treatment of KPP recommends the matching of interior diffusivities and their gradients to the KPP predicted values computed in the ocean surface boundary layer (OSBL). However, we find that difficulties in representing derivatives of rapidly changing diffusivities near the base of the OSBL can lead to loss of simulation fidelity. To mitigate this difficulty, we propose and evaluate two computationally simpler approaches: (1) match to the internal predicted diffusivity alone, (2) set the KPP diffusivity to zero at the OSBL base.
We find the KPP entrainment buoyancy flux to be sensitive to vertical grid resolution and details of how to diagnose the KPP boundary layer depth. We modify the KPP turbulent shear velocity parameterization to reduce resolution dependence. Additionally, an examination of LES vertical turbulent scalar flux budgets shows that the KPP parameterized non‐local tracer flux is incomplete due to the assumption that it solely redistributes the surface tracer flux. This result motivates further studies of the non‐local flux parameterization.
van Sebille, E, Stephen M Griffies, Ryan Abernathey, T P Adams, P Berloff, A Biastoch, B Blanke, Eric P Chassignet, Yu Cheng, C J Cotter, E Deleersnijder, K Döös, and Henri F Drake, et al., January 2018: Lagrangian ocean analysis: fundamentals and practices. Ocean Modelling, 121, DOI:10.1016/j.ocemod.2017.11.008. Abstract
Lagrangian analysis is a powerful way to analyse the output of ocean circulation models and other ocean velocity data such as from altimetry. In the Lagrangian approach, large sets of virtual particles are integrated within the three-dimensional, time-evolving velocity fields. Over several decades, a variety of tools and methods for this purpose have emerged. Here, we review the state of the art in the field of Lagrangian analysis of ocean velocity data, starting from a fundamental kinematic framework and with a focus on large-scale open ocean applications. Beyond the use of explicit velocity fields, we consider the influence of unresolved physics and dynamics on particle trajectories. We comprehensively list and discuss the tools currently available for tracking virtual particles. We then showcase some of the innovative applications of trajectory data, and conclude with some open questions and an outlook. The overall goal of this review paper is to reconcile some of the different techniques and methods in Lagrangian ocean analysis, while recognising the rich diversity of codes that have and continue to emerge, and the challenges of the coming age of petascale computing.
The Geophysical Fluid Dynamics Laboratory (GFDL) has recently developed two global coupled GCMs, FLOR and HiFLOR, which are now being utilized for climate research and seasonal predictions. Compared to their predecessor CM2.1, the new versions have improved ocean/atmosphere physics and numerics, and refinement of the atmospheric horizontal grid from 220 km (CM2.1) to 55 km (FLOR) and 26 km (HiFLOR). Both FLOR and HiFLOR demonstrate greatly improved simulations of the tropical Pacific annual‐mean climatology, with FLOR practically eliminating any equatorial cold bias in sea surface temperature. An additional model experiment (LOAR1) using FLOR's ocean/atmosphere physics, but with the atmospheric grid coarsened toward that of CM2.1, is used to further isolate the impacts of the refined atmospheric grid versus the improved physics and numerics. The improved ocean/atmosphere formulations are found to produce more realistic tropical Pacific patterns of sea surface temperature and rainfall, surface heat fluxes, ocean mixed layer depths, surface currents, and tropical instability wave (TIW) activity; enhance the near‐surface equatorial upwelling; and reduce the inter‐centennial warm drift of the tropical Pacific upper ocean. The atmospheric grid refinement further improves these features, and also improves the tropical Pacific surface wind stress, implied Ekman and Sverdrup transports, subsurface temperature and salinity structure, and heat advection in the equatorial upper ocean. The results highlight the importance of nonlocal air‐sea interactions in the tropical Pacific climate system, including the influence of off‐equatorial surface fluxes on the equatorial annual‐mean state. Implications are discussed for improving future simulations, observations, and predictions of tropical Pacific climate.
Yamamoto, A, J B Palter, Carolina O Dufour, Stephen M Griffies, Daniele Bianchi, M Claret, John P Dunne, Ivy Frenger, and Eric D Galbraith, October 2018: Roles of the ocean mesoscale in the horizontal supply of mass, heat, carbon and nutrients to the Northern Hemisphere subtropical gyres. Journal of Geophysical Research: Oceans, 123(10), DOI:10.1029/2018JC013969. Abstract
Horizontal transport at the boundaries of the subtropical gyres plays a crucial role in providing the nutrients that fuel gyre primary productivity, the heat that helps restratify the surface mixed layer, and the dissolved inorganic carbon (DIC) that influences air‐sea carbon exchange. Mesoscale eddies may be an important component of these horizontal transports; however, previous studies have not quantified the horizontal tracer transport due to eddies across the subtropical gyre boundaries. Here we assess the physical mechanisms that control the horizontal transport of mass, heat, nutrients and carbon across the North Pacific and North Atlantic subtropical gyre boundaries using the eddy‐rich ocean component of a climate model (GFDL's CM2.6) coupled to a simple biogeochemical model (mini‐BLING). Our results suggest that horizontal transport across the gyre boundaries supplies a substantial amount of mass and tracers to the ventilated layer of both Northern Hemisphere subtropical gyres, with the Kuroshio and Gulf Stream acting as main exchange gateways. Mass, heat, and DIC supply is principally driven by the time‐mean circulation, whereas nutrient transport differs markedly from the other tracers, as nutrients are mainly supplied to both subtropical gyres by down‐gradient eddy mixing across gyre boundaries. A budget analysis further reveals that the horizontal nutrient transport, combining the roles of both mean and eddy components, is responsible for more than three quarters of the total nutrient supply into the subtropical gyres, surpassing a recent estimate based on a coarse resolution model and thus further highlighting the importance of horizontal nutrient transport.
Dufour, Carolina O., Adele K Morrison, Stephen M Griffies, Ivy Frenger, Hannah Zanowski, and Michael Winton, October 2017: Preconditioning of the Weddell Sea polynya by the ocean mesoscale and dense water overflows. Journal of Climate, 30(19), DOI:10.1175/JCLI-D-16-0586.1. Abstract
TheWeddell Sea polynya is a large opening in the open-ocean sea ice cover associated with intense deep convection in the ocean. A necessary condition to form and maintain a polynya is the presence of a strong subsurface heat reservoir. This study investigates the processes that control the stratification and hence the build-up of the subsurface heat reservoir in theWeddell Sea. To do so, a climate model run for 200 years under preindustrial forcing with two eddying resolutions in the ocean (0.25° CM2.5 and 0.10° CM2.6) is investigated. Over the course of the simulation, CM2.6 develops two polynyas in the Weddell Sea, while CM2.5 exhibits quasi-continuous deep convection, but no polynyas, exemplifying that deep convection is not a sufficient condition for a polynya to occur. CM2.5 features a weaker subsurface heat reservoir than CM2.6 due to weak stratification associated with episodes of gravitational instability and enhanced vertical mixing of heat, resulting in an erosion of the reservoir. In contrast, in CM2.6, the water column is more stably stratified, allowing the subsurface heat reservoir to build up. The enhanced stratification in CM2.6 arises from its refined horizontal grid spacing and resolution of topography which allows, in particular, a better representation of the restratifying effect by transient mesoscale eddies and of the overflows of dense waters along the continental slope.
Goddard, P, Carolina O Dufour, Jianjun Yin, Stephen M Griffies, and Michael Winton, October 2017: CO2-Induced Ocean Warming of the Antarctic Continental Shelf in an Eddying Global Climate Model. Journal of Geophysical Research: Oceans, 122(10), DOI:10.1002/2017JC012849. Abstract
Ocean warming near the Antarctic ice shelves has critical implications for future ice sheet mass loss and global sea level rise. A global climate model with an eddying ocean is used to quantify the mechanisms contributing to ocean warming on the Antarctic continental shelf in an idealized 2xCO2 experiment. The results indicate that relatively large warm anomalies occur both in the upper 100 m and at depths above the shelf floor, which are controlled by different mechanisms. The near-surface ocean warming is primarily a response to enhanced onshore advective heat transport across the shelf break. The deep shelf warming is initiated by onshore intrusions of relatively warm Circumpolar Deep Water (CDW), in density classes that access the shelf, as well as the reduction of the vertical mixing of heat. CO2-induced shelf freshening influences both warming mechanisms. The shelf freshening slows vertical mixing by limiting gravitational instabilities and the upward diffusion of heat associated with CDW, resulting in the build-up of heat at depth. Meanwhile, freshening near the shelf break enhances the lateral density gradient of the Antarctic Slope Front (ASF) and disconnect isopycnals between the shelf and CDW, making cross-ASF heat exchange more difficult. However, at several locations along the ASF, the cross-ASF heat transport is less inhibited and heat can move onshore. Once onshore, lateral and vertical heat advection work to disperse the heat anomalies across the shelf region. Understanding the inhomogeneous Antarctic shelf warming will lead to better projections of future ice sheet mass loss.
Hewitt, Helene T., M J Bell, Eric P Chassignet, A Czaja, D Ferreira, and Stephen M Griffies, et al., December 2017: Will high-resolution global ocean models benefit coupled predictions on short-range to climate timescales?Ocean Modelling, 120, DOI:10.1016/j.ocemod.2017.11.002. Abstract
As the importance of the ocean in the weather and climate system is increasingly recognised, operational systems are now moving towards coupled prediction not only for seasonal to climate timescales but also for short-range forecasts. A three-way tension exists between the allocation of computing resources to refine model resolution, the expansion of model complexity/capability, and the increase of ensemble size. Here we review evidence for the benefits of increased ocean resolution in global coupled models, where the ocean component explicitly represents transient mesoscale eddies and narrow boundary currents. We consider lessons learned from forced ocean/sea-ice simulations; from studies concerning the SST resolution required to impact atmospheric simulations; and from coupled predictions. Impacts of the mesoscale ocean in western boundary current regions on the large-scale atmospheric state have been identified. Understanding of air-sea feedback in western boundary currents is modifying our view of the dynamics in these key regions. It remains unclear whether variability associated with open ocean mesoscale eddies is equally important to the large-scale atmospheric state. We include a discussion of what processes can presently be parameterised in coupled models with coarse resolution non-eddying ocean models, and where parameterizations may fall short. We discuss the benefits of resolution and identify gaps in the current literature that leave important questions unanswered.
MacKinnon, J A., M H Alford, Joseph K Ansong, Brian K Arbic, A Barna, B P Briegleb, F O Bryan, Maarten C Buijsman, Eric P Chassignet, Gokhan Danabasoglu, S Diggs, Stephen M Griffies, Robert Hallberg, S R Jayne, M Jochum, J Klymak, E Kunze, William G Large, Sonya Legg, B Mater, and Angelique Melet, et al., November 2017: Climate Process Team on Internal-Wave Driven Ocean Mixing. Bulletin of the American Meteorological Society, 98(11), DOI:10.1175/BAMS-D-16-0030.1. Abstract
Recent advances in our understanding of internal-wave driven turbulent mixing in the ocean interior are summarized. New parameterizations for global climate ocean models, and their climate impacts, are introduced.
Diapycnal mixing plays a primary role in the thermodynamic balance of the ocean and, consequently, in oceanic heat and carbon uptake and storage. Though observed mixing rates are on average consistent with values required by inverse models, recent attention has focused on the dramatic spatial variability, spanning several orders of magnitude, of mixing rates in both the upper and deep ocean. Away from ocean boundaries, the spatio-temporal patterns of mixing are largely driven by the geography of generation, propagation and dissipation of internal waves, which supply much of the power for turbulent mixing. Over the last five years and under the auspices of US CLIVAR, a NSF- and NOAA-supported Climate Process Team has been engaged in developing, implementing and testing dynamics-based parameterizations for internal-wave driven turbulent mixing in global ocean models. The work has primarily focused on turbulence 1) near sites of internal tide generation, 2) in the upper ocean related to wind-generated near inertial motions, 3) due to internal lee waves generated by low-frequency mesoscale flows over topography, and 4) at ocean margins. Here we review recent progress, describe the tools developed, and discuss future directions.
Orr, James C., R G Najjar, Olivier Aumont, Laurent Bopp, J L Bullister, Gokhan Danabasoglu, Scott C Doney, John P Dunne, J-C Dutay, H D Graven, Stephen M Griffies, and Jasmin G John, et al., June 2017: Biogeochemical protocols and diagnostics for the CMIP6 Ocean Model Intercomparison Project (OMIP). Geoscientific Model Development, 10(6), DOI:10.5194/gmd-10-2169-2017. Abstract
The Ocean Model Intercomparison Project (OMIP) focuses on the physics and biogeochemistry of the ocean component of Earth System Models participating in the sixth phase of the Coupled Model Intercomparison Project (CMIP6). OMIP aims to provide standard protocols and diagnostics for ocean models, while offering a forum to promote their common assessment and improvement. It also offers to compare solutions of the same ocean models when forced with reanalysis data (OMIP simulations) versus when integrated within fully coupled Earth System Models (CMIP6). Here we detail simulation protocols and diagnostics for OMIP's biogeochemical and inert chemical tracers. These passive-tracer simulations will be coupled online to ocean circulation models, initialized with observational data or output from a model spin up, and forced by repeating the 1948–2009 surface fluxes of heat, fresh water, and momentum. These so-called OMIP-BGC simulations include three inert chemical tracers (CFC-11, CFC-12, SF6 and biogeochemical tracers (e.g., dissolved inorganic carbon, carbon isotopes, alkalinity, nutrients, and oxygen). Modelers will use their preferred prognostic BGC model but should follow common guidelines for gas exchange and carbonate chemistry. Simulations include both natural and total carbon tracers. The required forced simulation (omip1) will be initialized with gridded observational climatologies. An optional forced simulation (omip1-spunup) will be initialized instead with BGC fields from a long model spin up, preferably for 2000 years or more and forced by repeating the same 62-year meteorological forcing. That optional run will also include abiotic tracers of total dissolved inorganic carbon and radiocarbon, CTabio and 14CTabio, to assess deep-ocean ventilation and distinguish the role of physics vs. biology. These simulations will be forced by observed atmospheric histories of the three inert gases and CO2 as well as carbon isotope ratios of CO2. OMIP-BGC simulation protocols are founded on those from previous phases of the Ocean Carbon-Cycle Model Intercomparison Project. They have been merged and updated to reflect improvements concerning gas exchange, carbonate chemistry, and new data for initial conditions and atmospheric gas histories. Code is provided to facililtate their implementation.
Spence, P, Ryan M Holmes, Andrew McC Hogg, Stephen M Griffies, K D Stewart, and Matthew H England, August 2017: Localized rapid warming of West Antarctic subsurface waters by remote winds. Nature Climate Change, 7(8), DOI:10.1038/nclimate3335. Abstract
The highest rates of Antarctic glacial ice mass loss are occurring to the west of the Antarctica Peninsula in regions where warming of subsurface continental shelf waters is also largest. However, the physical mechanisms responsible for this warming remain unknown. Here we show how localized changes in coastal winds off East Antarctica can produce significant subsurface temperature anomalies (>2 °C) around much of the continent. We demonstrate how coastal-trapped barotropic Kelvin waves communicate the wind disturbance around the Antarctic coastline. The warming is focused on the western flank of the Antarctic Peninsula because the circulation induced by the coastal-trapped waves is intensified by the steep continental slope there, and because of the presence of pre-existing warm subsurface water offshore. The adjustment to the coastal-trapped waves shoals the subsurface isotherms and brings warm deep water upwards onto the continental shelf and closer to the coast. This result demonstrates the vulnerability of the West Antarctic region to a changing climate.
Improvements in the horizontal resolution of global ocean models, motivated by the horizontal resolution requirements for specific flow features, has advanced modelling capabilities into the dynamical regime dominated by mesoscale variability. In contrast, the choice of the vertical grid remains a subjective choice, and it is not clear that efforts to improve vertical resolution adequately support their horizontal counterparts. Indeed, considering that the bulk of the vertical ocean dynamics (including convection) are parameterized, it is not immediately obvious what the vertical grid is supposed to resolve. Here, we propose that the primary purpose of the vertical grid in a hydrostatic ocean model is to resolve the vertical structure of horizontal flows, rather than to resolve vertical motion. With this principle we construct vertical grids based on their abilities to represent baroclinic modal structures commensurate with the theoretical capabilities of a given horizontal grid. This approach is designed to ensure that the vertical grids of global ocean models complement (and, importantly, to not undermine) the resolution capabilities of the horizontal grid. We find that for z-coordinate global ocean models, at least 50 well-positioned vertical levels are required to resolve the first baroclinic mode, with an additional 25 levels per subsequent mode. High-resolution ocean-sea ice simulations are used to illustrate some of the dynamical enhancements gained by improving the vertical resolution of a 1/10° global ocean model. These enhancements include substantial increases in the sea surface height variance (∼ 30% increase south of 40°S), the barotropic and baroclinic eddy kinetic energies (up to 200% increase on and surrounding the Antarctic continental shelf and slopes), and the overturning streamfunction in potential density space (near-tripling of the Antarctic Bottom Water cell at 65°S).
Swapna, P, J Jyoti, R Krishnan, N Sandeep, and Stephen M Griffies, October 2017: Multi-decadal weakening of Indian summer monsoon circulation induces an increasing northern Indian Ocean sea level. Geophysical Research Letters, 44(20), DOI:10.1002/2017GL074706. Abstract
North Indian Ocean sea level has shown significant increase during last 3-4 decades. Analyses of long-term climate datasets and ocean model sensitivity experiments identify a mechanism for multi-decadal sea level variability relative to global mean. Our results indicate that north Indian Ocean sea level rise is accompanied by a weakening summer monsoon circulation. Given that Indian Ocean meridional heat transport is primarily regulated by the annual cycle of monsoon winds, weakening of summer monsoon circulation has resulted in reduced upwelling off Arabia and Somalia and decreased southward heat-transport, and corresponding increase of heat storage in the north Indian Ocean. These changes in-turn lead to increased retention of heat and increased thermosteric sea level rise in the north Indian Ocean, especially in the Arabian Sea. These findings imply that rising north Indian Ocean sea level due to weakening of monsoon circulation demand adaptive strategies to enable a resilient South Asian population.
Upwelling of global deep waters to the sea surface in the Southern Ocean closes the global overturning circulation and is fundamentally important for oceanic uptake of carbon and heat, nutrient resupply for sustaining oceanic biological production, and the melt rate of ice shelves. However, the exact pathways and role of topography in Southern Ocean upwelling remain largely unknown. Here we show detailed upwelling pathways in three dimensions, using hydrographic observations and particle tracking in high-resolution models. The analysis reveals that the northern-sourced deep waters enter the Antarctic Circumpolar Current via southward flow along the boundaries of the three ocean basins, before spiraling southeastward and upward through the Antarctic Circumpolar Current. Upwelling is greatly enhanced at five major topographic features, associated with vigorous mesoscale eddy activity. Deep water reaches the upper ocean predominantly south of the Antarctic Circumpolar Current, with a spatially nonuniform distribution. The timescale for half of the deep water to upwell from 30° S to the mixed layer is ~60–90 years.
Relatively rapid re-emergence of anthropogenic carbon (Cant) in the Equatorial Pacific is of potential importance for its impact on the carbonate buffering capacity of surface seawater, and thereby impeding the ocean's ability to further absorb Cant from the atmosphere. We explore the mechanisms sustaining Cant re-emergence (upwelling) from the thermocline to surface layers by applying water mass transformation diagnostics to a global ocean/sea-ice/biogeochemistry model. We find that the upwelling rate of Cant (0.4 PgC yr-1) from the thermocline to the surface layer is almost twice as large as air-sea Cant fluxes (0.203 PgC yr-1). The upwelling of Cant from the thermocline to the surface layer can be understood as a two-step process: the first being due to diapycnal diffusive transformation fluxes and the second due to surface buoyancy fluxes. We also find that this re-emergence of Cant decreases dramatically during the 1982/1983 and 1997/1998 El Niño events.
Danabasoglu, Gokhan, Stephen G Yeager, W M Kim, Erik Behrens, M Bentsen, D Bi, A Biastoch, R Bleck, C Böning, A Bozec, V M Canuto, Christophe Cassou, Eric P Chassignet, A C Coward, S Danilov, N Diansky, H Drange, Riccardo Farneti, E Fernandez, P G Fogli, G Forget, Yosuke Fujii, Stephen M Griffies, A Gusev, Patrick Heimbach, A Howard, M Ilicak, T Jung, Alicia R Karspeck, M Kelley, William G Large, A Leboissetier, J Lu, G Madec, S J Marsland, S Masina, A Navarra, A J George Nurser, Anna Pirani, Anastasia Romanou, D Salas y Mélia, and Bonita L Samuels, et al., January 2016: North Atlantic Simulations in Coordinated Ocean-ice Reference Experiments phase II (CORE-II). Part II: Inter-Annual to Decadal Variability. Ocean Modelling, 97, DOI:10.1016/j.ocemod.2015.11.007. Abstract
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 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.
Model and observational studies have concluded that geothermal heating significantly alters the global overturning circulation and the properties of the widely-distributed Antarctic Bottom Waters. Here we test two distinct geothermal heat flux datasets under different experimental designs in a fully coupled model that mimics the control run of a typical Coupled Model Intercomparison Project (CMIP) climate model. Our regional analysis reveals that bottom temperature and transport changes, due to the inclusion of geothermal heating, are propagated throughout the water column, most prominently in the Southern Ocean, with the background density structure and major circulation pathways acting as drivers of these changes. Whilst geothermal heating enhances Southern Ocean abyssal overturning circulation by 20-50%, upwelling of warmer deep waters and cooling of upper ocean waters within the Antarctic Circumpolar Current (ACC) region decrease its transport by 3 to 5 Sv. The transient responses in regional bottom temperature increases exceed 0.1°C. The combination of large scale features that we show act to transport anomalies far from their geothermal source all exist in the Southern Ocean. Such features include steeply sloping isopycnals, weak abyssal stratification, voluminous southward flowing deep waters and exported bottom waters, the ACC, and the polar gyres. Recently the Southern Ocean has been identified as a prime region for deep ocean warming; geothermal heating should be included in climate models to ensure accurate representation of these abyssal temperature changes.
Gregory, Jonathan M., N Bouttes-Mauhourat, Stephen M Griffies, Helmuth Haak, William J Hurlin, J H Jungclaus, M Kelley, W G Lee, J Marshall, Anastasia Romanou, Oleg A Saenko, Detlef Stammer, and Michael Winton, November 2016: The Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP) contribution to CMIP6: Investigation of sea-level and ocean climate change in response to CO2 forcing. Geoscientific Model Development, 9(11), DOI:10.5194/gmd-9-3993-2016. Abstract
The Flux-Anomaly-Forced Model Intercomparison Project (FAFMIP) aims to investigate the spread in simulations of sea-level and ocean climate change in response to CO2 forcing by atmosphere-ocean general circulation models (AOGCMs). It is particularly motivated by the uncertainties in projections of ocean heat uptake, global-mean sea-level rise due to thermal expansion and the geographical patterns of sea-level change due to ocean density and circulation change. FAFMIP has three tier-1 experiments, in which prescribed surface flux perturbations of momentum, heat and freshwater respectively are applied to the ocean in separate AOGCM simulations. All other conditions are as in the pre-industrial control. The prescribed fields are typical of pattern and magnitude of changes in these fluxes projected by AOGCMs for doubled CO2 concentration. Five groups have tested the experimental design with existing AOGCMs. Their results show diversity in the pattern and magnitude of changes, with some common qualitative features. Heat and water flux perturbation cause the dipole in sea-level change in the North Atlantic, while momentum and heat flux perturbation cause the gradient across the Antarctic Circumpolar Current. The Atlantic Meridional Overturning Circulation (AMOC) declines in response to the heat flux perturbation, and there is a strong positive feedback on this effect due to the consequent cooling of sea surface temperature in the North Atlantic, which enhances the local heat input to the ocean. The momentum and water flux perturbations do not substantially affect the AMOC. Heat is taken up largely as a passive tracer in the Southern Ocean, which is the region of greatest heat input, but elsewhere heat is actively redistributed towards lower latitude. Future analysis of these and other phenomena with the wider range of CMIP6 FAFMIP AOGCMs will benefit from new diagnostics of temperature and salinity tendencies, which will enable investigation of the model spread in behaviour in terms of physical processes as formulated in the models.
Griffies, Stephen M., Gokhan Danabasoglu, Paul J Durack, Alistair Adcroft, V Balaji, C Böning, Eric P Chassignet, Enrique N Curchitser, Julie Deshayes, H Drange, Baylor Fox-Kemper, Peter J Gleckler, Jonathan M Gregory, Helmuth Haak, Robert Hallberg, Helene T Hewitt, David M Holland, Tatiana Ilyina, J H Jungclaus, Y Komuro, John P Krasting, William G Large, S J Marsland, S Masina, Trevor J McDougall, A J George Nurser, James C Orr, Anna Pirani, Fangli Qiao, Ronald J Stouffer, Karl E Taylor, Anne Marie Treguier, Hiroyuki Tsujino, P Uotila, M Valdivieso, Michael Winton, and Stephen G Yeager, September 2016: OMIP contribution to CMIP6: experimental and diagnostic protocol for the physical component of the Ocean Model Intercomparison Project. Geoscientific Model Development, 9(9), DOI:10.5194/gmd-9-3231-2016. Abstract
The Ocean Model Intercomparison Project (OMIP) aims to provide a framework for evaluating, understanding, and improving the ocean and sea-ice components of global climate and earth system models contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses these aims in two complementary manners: (A) by providing an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing, (B) by providing a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) offering details for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows that of the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II have become the standard method to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP (Scenario MIP), as well as the ocean-sea ice OMIP simulations. The bulk of this paper offers scientific rationale for saving these diagnostics.
Ilicak, M, H Drange, Q Wang, R Gerdes, Y Aksenov, David A Bailey, M Bentsen, A Biastoch, A Bozec, C Böning, Christophe Cassou, Eric P Chassignet, A C Coward, B Curry, Gokhan Danabasoglu, S Danilov, E Fernandez, P G Fogli, Yosuke Fujii, Stephen M Griffies, Doroteaciro Iovino, Alexandra Jahn, T Jung, William G Large, Craig Lee, C Lique, J Lu, S Masina, A J George Nurser, C Roth, D Salas y Mélia, and Bonita L Samuels, et al., April 2016: An assessment of the Arctic Ocean in a suite of interannual CORE-II simulations. Part III: Hydrography and fluxes. Ocean Modelling, 100, DOI:10.1016/j.ocemod.2016.02.004. Abstract
In this paper we compare the simulated Arctic Ocean in fifteen global ocean-sea ice models in the framework of the Coordinated Ocean-ice Reference Experiments, phase II (CORE-II). Most of these models are the ocean and sea-ice components of the coupled climate models used in the Coupled Model Intercomparison Project Phase 5 (CMIP5) experiments. We mainly focus on the hydrography of the Arctic interior, the state of Atlantic Water layer and heat and volume transports at the gateways of the Davis Strait, the Bering Strait, the Fram Strait and the Barents Sea Opening. We found that there is a large spread in temperature in the Arctic Ocean between the models, and generally large differences compared to the observed temperature at intermediate depths. Warm bias models have a strong temperature anomaly of inflow of the Atlantic Water entering the Arctic Ocean through the Fram Strait. Another process that is not represented accurately in the CORE-II models is the formation of cold and dense water, originating on the eastern shelves. In the cold bias models, excessive cold water forms in the Barents Sea and spreads into the Arctic Ocean through the St. Anna Through. There is a large spread in the simulated mean heat and volume transports through the Fram Strait and the Barents Sea Opening. The models agree more on the decadal variability, to a large degree dictated by the common atmospheric forcing. We conclude that the CORE-II model study helps us to understand the crucial biases in the Arctic Ocean. The current coarse resolution state-of-the-art ocean models need to be improved in accurate representation of the Atlantic Water inflow into the Arctic and density currents coming from the shelves.
The Southern Ocean plays a dominant role in anthropogenic oceanic heat uptake. Strong northward transport of the heat content anomaly limits warming of the sea surface temperature in the uptake region and allows the heat uptake to be sustained. Using an eddy-rich global climate model, the processes controlling the northward transport and convergence of the heat anomaly in the mid-latitude Southern Ocean are investigated in an idealized 1% yr−1 increasing CO2 simulation. Heat budget analyses reveal that different processes dominate to the north and south of the main convergence region. The heat transport northward from the uptake region in the south is driven primarily by passive advection of the heat content anomaly by the existing time mean circulation, with a smaller 20% contribution from enhanced upwelling. The heat anomaly converges in the mid-latitude deep mixed layers, because there is not a corresponding increase in the mean heat transport out of the deep mixed layers northward into the mode waters. To the north of the deep mixed layers, eddy processes drive the warming and account for nearly 80% of the northward heat transport anomaly. The eddy transport mechanism results from a reduction in both the diffusive and advective southward eddy heat transports, driven by decreasing isopycnal slopes and decreasing along-isopycnal temperature gradients on the northern edge of the peak warming.
The Intergovernmental Panel on Climate Change (IPCC) fifth assessment of projected global and regional ocean temperature change is based on global climate models that have coarse (∼100-km) ocean and atmosphere resolutions. In the Northwest Atlantic, the ensemble of global climate models has a warm bias in sea surface temperature due to a misrepresentation of the Gulf Stream position; thus, existing climate change projections are based on unrealistic regional ocean circulation. Here we compare simulations and an atmospheric CO2 doubling response from four global climate models of varying ocean and atmosphere resolution. We find that the highest resolution climate model (∼10-km ocean, ∼50-km atmosphere) resolves Northwest Atlantic circulation and water mass distribution most accurately. The CO2 doubling response from this model shows that upper-ocean (0-300 m) temperature in the Northwest Atlantic Shelf warms at a rate nearly twice as fast as the coarser models and nearly three times faster than the global average. This enhanced warming is accompanied by an increase in salinity due to a change in water mass distribution that is related to a retreat of the Labrador Current and a northerly shift of the Gulf Stream. Both observations and the climate model demonstrate a robust relationship between a weakening Atlantic Meridional Overturning Circulation (AMOC) and an increase in the proportion of Warm-Temperate Slope Water entering the Northwest Atlantic Shelf. Therefore, prior climate change projections for the Northwest Atlantic may be far too conservative. These results point to the need to improve simulations of basin and regional-scale ocean circulation.
Tseng, Y-H, Hongyang Lin, Han-Ching Chen, K Thompson, M Bentsen, C Böning, A Bozec, Christophe Cassou, Eric P Chassignet, C H Chow, Gokhan Danabasoglu, S Danilov, Riccardo Farneti, P G Fogli, Yosuke Fujii, Stephen M Griffies, M Ilicak, T Jung, S Masina, A Navarra, L Patara, and Bonita L Samuels, et al., August 2016: North and Equatorial Pacific Ocean Circulation in the CORE-II Hindcast Simulations. Ocean Modelling, 104, DOI:10.1016/j.ocemod.2016.06.003. 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.
Wang, Q, M Ilicak, R Gerdes, H Drange, Y Aksenov, David A Bailey, M Bentsen, A Biastoch, A Bozec, C Böning, Christophe Cassou, Eric P Chassignet, A C Coward, B Curry, Gokhan Danabasoglu, S Danilov, E Fernandez, P G Fogli, Yosuke Fujii, Stephen M Griffies, Doroteaciro Iovino, Alexandra Jahn, T Jung, William G Large, Craig Lee, C Lique, J Lu, S Masina, A J George Nurser, B Rabe, C Roth, D Salas y Mélia, and Bonita L Samuels, et al., March 2016: An assessment of the Arctic Ocean in a suite of interannual CORE-II simulations. Part I: Sea ice and solid freshwater. Ocean Modelling, 99, DOI:10.1016/j.ocemod.2015.12.008. Abstract
The Arctic Ocean simulated in fourteen global ocean-sea ice models in the framework of the Coordinated Ocean-ice Reference Experiments, phase II (CORE II) is analyzed. The focus is on the Arctic sea ice extent, the solid freshwater (FW) sources and solid freshwater content (FWC). Available observations are used for model evaluation. The variability of sea ice extent and solid FW budget is more consistently reproduced than their mean state in the models. The descending trend of September sea ice extent is well simulated in terms of the model ensemble mean. Models overestimating sea ice thickness tend to underestimate the descending trend of September sea ice extent. The models underestimate the observed sea ice thinning trend by a factor of two. When averaged on decadal time scales, the variation of Arctic solid FWC is contributed by those of both sea ice production and sea ice transport, which are out of phase in time. The solid FWC decreased in the recent decades, caused mainly by the reduction in sea ice thickness. The models did not simulate the acceleration of sea ice thickness decline, leading to an underestimation of solid FWC trend after 2000. The common model behaviour, including the tendency to underestimate the trend of sea ice thickness and March sea ice extent, remains to be improved.
Wang, Q, M Ilicak, R Gerdes, H Drange, Y Aksenov, David A Bailey, M Bentsen, A Biastoch, A Bozec, C Böning, Christophe Cassou, Eric P Chassignet, A C Coward, B Curry, Gokhan Danabasoglu, S Danilov, E Fernandez, P G Fogli, Yosuke Fujii, Stephen M Griffies, Doroteaciro Iovino, Alexandra Jahn, T Jung, William G Large, Craig Lee, C Lique, J Lu, S Masina, A J George Nurser, B Rabe, C Roth, D Salas y Mélia, and Bonita L Samuels, et al., March 2016: An assessment of the Arctic Ocean in a suite of interannual CORE-II simulations. Part II: Liquid freshwater. Ocean Modelling, 99, DOI:10.1016/j.ocemod.2015.12.009. Abstract
The Arctic Ocean simulated in fourteen global ocean-sea ice models in the framework of the Coordinated Ocean-ice Reference Experiments, phase II (CORE-II) is analyzed in this study. The focus is on the Arctic liquid freshwater (FW) sources and freshwater content (FWC). The models agree on the interannual variability of liquid FW transport at the gateways where the ocean volume transport determines the FW transport variability. The variation of liquid FWC is induced by both the surface FW flux (associated with sea ice production) and lateral liquid FW transport, which are in phase when averaged on decadal time scales. The liquid FWC shows an increase starting from the mid-1990s, caused by the reduction of both sea ice formation and liquid FW export, with the former being more significant in most of the models. The mean state of the FW budget is less consistently simulated than the temporal variability. The model ensemble means of liquid FW transport through the Arctic gateways compare well with observations. On average, the models have too high mean FWC, weaker upward trends of FWC in the recent decade than the observation, and low consistency in the temporal variation of FWC spatial distribution, which needs to be further explored for the purpose of model development.
Biastoch, A, J V Durgadoo, Adele K Morrison, E van Sebille, W Weijer, and Stephen M Griffies, December 2015: Atlantic multi-decadal oscillation covaries with Agulhas leakage. Nature Communications, 6, 10082, DOI:10.1038/ncomms10082. Abstract
The interoceanic transfer of seawater between the Indian Ocean and the Atlantic, ‘Agulhas leakage’, forms a choke point for the overturning circulation in the global ocean. Here, by combining output from a series of high-resolution ocean and climate models with in situ and satellite observations, we construct a time series of Agulhas leakage for the period 1870–2014. The time series demonstrates the impact of Southern Hemisphere westerlies on decadal timescales. Agulhas leakage shows a correlation with the Atlantic Multi-decadal Oscillation on multi-decadal timescales; the former leading by 15 years. This is relevant for climate in the North Atlantic.
Downes, S M., Riccardo Farneti, P Uotila, Stephen M Griffies, S J Marsland, David A Bailey, Erik Behrens, M Bentsen, D Bi, A Biastoch, C Böning, A Bozec, V M Canuto, Eric P Chassignet, Gokhan Danabasoglu, S Danilov, N Diansky, H Drange, P G Fogli, A Gusev, A Howard, M Ilicak, T Jung, M Kelley, William G Large, A Leboissetier, Matthew C Long, J Lu, S Masina, A Mishra, A Navarra, A J George Nurser, L Patara, and Bonita L Samuels, et al., October 2015: An assessment of Southern Ocean water masses and sea ice during 1988-2007 in a suite of inter-annual CORE-II simulations. Ocean Modelling, 94, DOI:10.1016/j.ocemod.2015.07.022. Abstract
We characterize the representation of the Southern Ocean water mass structure and sea ice within a suite of 15 global ocean-ice models run with the Coordinated Ocean-ice Reference Experiment Phase II (CORE-II) protocol. The main focus is the representation of the present (1988–2007) mode and intermediate waters, thus framing an analysis of winter and summer mixed layer depths; temperature, salinity, and potential vorticity structure; and temporal variability of sea ice distributions. We also consider the inter-annual variability over the same 20 year period. Comparisons are made between models as well as to observation-based analyses where available.
The CORE-II models exhibit several biases relative to Southern Ocean observations, including an underestimation of the model mean mixed layer depths of mode and intermediate water masses in March (associated with greater ocean surface heat gain), and an overestimation in September (associated with greater high latitude ocean heat loss and a more northward winter sea-ice extent). In addition, the models have cold and fresh/ warm and salty water column biases centered near 50°S. Over the 1988–2007 period, the CORE-II models consistently simulate spatially variable trends in sea-ice concentration, surface freshwater fluxes, mixed layer depths, and 200–700 m ocean heat content. In particular, sea-ice coverage around most of the Antarctic continental shelf is reduced, leading to a cooling and freshening of the near surface waters. The shoaling of the mixed layer is associated with increased surface buoyancy gain, except in the Pacific where sea ice is also influential. The models are in disagreement, despite the common CORE-II atmospheric state, in their spatial pattern of the 20-year trends in the mixed layer depth and sea-ice.
This study examines the role of processes transporting tracers across the Polar Front (PF) in the depth interval between the surface and major topographic sills, which we refer to as the “PF core”. A preindustrial control simulation of an eddying climate model coupled to a biogeochemical model (CM2.6-miniBLING, 0.1° ocean model) is used to investigate the transport of heat, carbon, oxygen and phosphate across the PF core, with a particular focus on the role of mesoscale eddies. We find that the total transport across the PF core results from an ubiquitous Ekman transport that drives the upwelled tracers to the north, and a localized opposing eddy transport that induces tracer leakages to the south at major topographic obstacles. In the Ekman layer, the southward eddy transport only partially compensates the northward Ekman transport, while below the Ekman layer, the southward eddy transport dominates the total transport but remains much smaller in magnitude than the near-surface northward transport. Most of the southward branch of the total transport is achieved below the PF core, mainly through geostrophic currents. We find that the eddy diffusive transport reinforces the southward eddy advective transport for carbon and heat, and opposes it for oxygen and phosphate. Eddy advective transport is likely to be the leading-order component of eddy-induced transport for all four tracers. However, eddy diffusive transport may provide a significant contribution to the southward eddy heat transport due to strong along-isopycnal temperature gradients.
Farneti, Riccardo, S M Downes, Stephen M Griffies, S J Marsland, Erik Behrens, M Bentsen, D Bi, A Biastoch, C Böning, A Bozec, V M Canuto, Eric P Chassignet, Gokhan Danabasoglu, S Danilov, N Diansky, H Drange, P G Fogli, A Gusev, Robert Hallberg, A Howard, M Ilicak, T Jung, M Kelley, William G Large, A Leboissetier, Matthew C Long, J Lu, S Masina, A Mishra, A Navarra, A J George Nurser, L Patara, and Bonita L Samuels, et al., September 2015: An assessment of Antarctic Circumpolar Current and Southern Ocean Meridional Overturning Circulation during 1958–2007 in a suite of interannual CORE-II simulations. Ocean Modelling, 93, DOI:10.1016/j.ocemod.2015.07.009. Abstract
In the framework of the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II), we present an analysis of the representation of the Antarctic Circumpolar Current (ACC) and Southern Ocean Meridional Overturning Circulation (MOC) in a suite of seventeen global ocean-sea ice models. We focus on the mean, variability and trends of both the ACC and MOC over the 1958–2007 period, and discuss their relationship with the surface forcing. We aim to quantify the degree of eddy saturation and eddy compensation in the models participating in CORE-II, and compare our results with available observations, previous fine-resolution numerical studies and theoretical constraints. Most models show weak ACC transport sensitivity to changes in forcing during the past five decades, and they can be considered to be in an eddy saturated regime. Larger contrasts arise when considering MOC trends, with a majority of models exhibiting significant strengthening of the MOC during the late 20th and early 21st century. Only a few models show a relatively small sensitivity to forcing changes, responding with an intensified eddy-induced circulation that provides some degree of eddy compensation, while still showing considerable decadal trends. Both ACC and MOC interannual variability are largely controlled by the Southern Annular Mode (SAM). Based on these results, models are clustered into two groups. Models with constant or two-dimensional (horizontal) specification of the eddy-induced advection coefficient κ show larger ocean interior decadal trends, larger ACC transport decadal trends and no eddy compensation in the MOC. Eddy-permitting models or models with a three-dimensional time varying κ show smaller changes in isopycnal slopes and associated ACC trends, and partial eddy compensation. As previously argued, a constant in time or space κ is responsible for a poor representation of mesoscale eddy effects and cannot properly simulate the sensitivity of the ACC and MOC to changing surface forcing. Evidence is given for a larger sensitivity of the MOC as compared to the ACC transport, even when approaching eddy saturation. Future process studies designed for disentangling the role of momentum and buoyancy forcing in driving the ACC and MOC are proposed.
Goddard, P, Jianjun Yin, Stephen M Griffies, and Shaoqing Zhang, February 2015: An extreme event of sea-level rise along the Northeast coast of North America in 2009–2010. Nature Communications, 6, 6346, DOI:10.1038/ncomms7346. Abstract
The coastal sea levels along the Northeast Coast of North America show significant year-to-year fluctuations in a general upward trend. The analysis of long-term tide gauge records identified an extreme sea-level rise (SLR) event during 2009–10. Within this 2-year period, the coastal sea level north of New York City jumped by 128 mm. This magnitude of interannual SLR is unprecedented (a 1-in-850 year event) during the entire history of the tide gauge records. Here we show that this extreme SLR event is a combined effect of two factors: an observed 30% downturn of the Atlantic meridional overturning circulation during 2009–10, and a significant negative North Atlantic Oscillation index. The extreme nature of the 2009–10 SLR event suggests that such a significant downturn of the Atlantic overturning circulation is very unusual. During the twenty-first century, climate models project an increase in magnitude and frequency of extreme interannual SLR events along this densely populated coast.
We characterize impacts on heat in the ocean climate system from transient ocean mesoscale eddies. Our tool is a suite of centennial-scale 1990 radiatively forced numerical climate simulations from three GFDL coupled models comprising the CM2-O model suite. CM2-O models differ in their ocean resolution: CM2.6 uses a 0.1° ocean grid, CM2.5 uses an intermediate grid with 0.25° spacing, and CM2-1deg uses a nominally 1.0° grid.
Analysis of the ocean heat budget reveals that mesoscale eddies act to transport heat upward in a manner that partially compensates (or offsets) for the downward heat transport from the time mean currents. Stronger vertical eddy heat transport in CM2.6 relative to CM2.5 accounts for the significantly smaller temperature drift in CM2.6. The mesoscale eddy parameterization used in CM2-1deg also imparts an upward heat transport, yet it differs systematically from that found in CM2.6. This analysis points to the fundamental role that ocean mesoscale features play in transient ocean heat uptake. In general, the more accurate simulation found in CM2.6 provides an argument for either including a rich representation of the ocean mesoscale in model simulations of the mean and transient climate, or for employing parameterizations that faithfully reflect the role of eddies in both lateral and vertical heat transport.
South Atlantic transports, as simulated by a global ocean-sea ice model forced with the Coordinated Ocean-ice Reference Experiments version 2 (CORE-II) interannually varying air-sea reanalysis data sets, are analyzed for the period 1958–2007. The ocean-sea ice model is configured at three different resolutions: from eddy-permitting to coarsened grid spacing. A particular focus is given to the effect of eddy fluxes and inter-ocean exchanges on the South Atlantic Meridional Overturning Circulation (SAMOC), as well as on the main factors contributing to the interannual variability during the integration period. Differences between refined and coarsened grid spacing models are more evident in coastal areas and in regions of high eddy activities. Major discrepancies are associated to both the parametrization of eddy fluxes and the coarse representation of the bathymetry. The refined grid spacing model produces higher values of both SAMOC index, defined as the maximum of the zonally-integrated northward cumulative volume transport (CVT) from surface to bottom across ∼ 34° S, and meridional heat transport (MHT). All models show high correlations between SAMOC index and MHT, as well as a strengthening of the transports in the 1980-2007 period. The strengthening of the SAMOC index is mainly dominated by surface and mode waters in all models. In surface and intermediate layers, the regions contributing to this trend are located east of 40° W. These changes are compensated by the strengthening of the poleward transport in deeper layers, mostly in the western part of the basin. The MHT trend is connected with the combined effect of a heat transport increase through the Drake Passage and a reduction of the heat loss through the eastern section between Africa and Antarctica, mainly associated with a strengthening in heat entering into the basin through the Agulhas system.
Snow, K, Andrew McC Hogg, S M Downes, B M Sloyan, M L Bates, and Stephen M Griffies, May 2015: Sensitivity of abyssal water masses to overflow parameterisations. Ocean Modelling, 89, DOI:10.1016/j.ocemod.2015.03.004. Abstract
Antarctic Bottom Water (AABW) and North Atlantic Deep Water (NADW) control the abyssal limb of the global overturning circulation and play a major role in oceanic heat uptake and carbon storage. However, current general circulation models are unable to resolve the observed AABW and NADW formation and transport processes. One key process, that of overflows, motivates the application of overflow parameterisations. We present a sensitivity study of both AABW and NADW properties to three current parameterisations using a z∗z∗-coordinate ocean-sea ice model within a realistic-topography sector of the Atlantic Ocean.
Overflow parameterisations that affect only tracer equations are compared to a fully dynamical Lagrangian point particle method. An overflow parameterisation involving partial convective mixing of tracers is most efficient at transporting dense NADW water downslope. This parameterisation leads to a maximum mean increase in density in the north of 0.027 kg m−3 and a decrease in age of 525 years (53%). The relative change in density and age in the south is less than 30% of that in the north for all overflow parameterisations. The reduced response in the south may result from the differing dense water formation and overflow characteristics of AABW compared to NADW. Alternative approaches may be necessary to improve AABW representation in z∗z∗-coordinate ocean climate models.
Bueti, M R., Isaac Ginis, L Rothstein, and Stephen M Griffies, September 2014: Tropical Cyclone-Induced Thermocline Warming and its Regional and Global Impacts. Journal of Climate, 27(18), DOI:10.1175/JCLI-D-14-00152.1. Abstract
Strong surface winds of a hurricane locally cool the surface and warm the subsurface waters via turbulent mixing processes. While the surface cool anomalies generally decay in roughly a month, the warm subsurface anomalies can persist over a seasonal cycle. We examine questions related to the magnitude and cumulative footprint of subsurface warm anomalies forced by tropical cyclones during the combined global tropical cyclone seasons, making use of a global ocean model forced by tropical cyclones.
Simulations of the 2004-2005 tropical cyclone season are conducted with and without tropical cyclone wind forcing, blended with the daily Coordinated Ocean-ice Reference Experiments (COREs) atmospheric state. Physical characteristics of cyclone-forced surface and subsurface anomalies are elucidated. In particular, we examine the spatial extent and magnitude of tropical cyclone forced subsurface warm anomalies over the course of an entire season. This analysis allows us to estimate the contribution of cyclone-induced anomalies to the ocean heat content and sea surface temperature, and to understand anomalous meridional heat transport.
Globally, there is a maximum accumulated heat uptake 4.1·1021J, with the greatest regional contributions in the North Atlantic (1.7·1021J), West Pacific (1.5·1021J), and East Pacific (1.7·1021J). We find an export of heat from the subtropics to the tropics via rapid advective pathways, most notably in the West Pacific. These warm anomalies tend to remain in the equatorial band, with potential implications for the tropical climate system.
Simulation characteristics from eighteen global ocean–sea-ice coupled models are presented with a focus on the mean Atlantic meridional overturning circulation (AMOC) and other related fields in the North Atlantic. These experiments use inter-annually varying atmospheric forcing data sets for the 60-year period from 1948 to 2007 and are performed as contributions to the second phase of the Coordinated Ocean-ice Reference Experiments (CORE-II). The protocol for conducting such CORE-II experiments is summarized. Despite using the same atmospheric forcing, the solutions show significant differences. As most models also differ from available observations, biases in the Labrador Sea region in upper-ocean potential temperature and salinity distributions, mixed layer depths, and sea-ice cover are identified as contributors to differences in AMOC. These differences in the solutions do not suggest an obvious grouping of the models based on their ocean model lineage, their vertical coordinate representations, or surface salinity restoring strengths. Thus, the solution differences among the models are attributed primarily to use of different subgrid scale parameterizations and parameter choices as well as to differences in vertical and horizontal grid resolutions in the ocean models. Use of a wide variety of sea-ice models with diverse snow and sea-ice albedo treatments also contributes to these differences. Based on the diagnostics considered, the majority of the models appear suitable for use in studies involving the North Atlantic, but some models require dedicated development effort.
Fan, Yalin, Shian-Jiann Lin, Stephen M Griffies, and Mark A Hemer, May 2014: Simulated Global Swell and Wind Sea Climate and Their responses to Anthropogenic Climate Change at the End of the 21st Century. Journal of Climate, 27(10), DOI:10.1175/JCLI-D-13-00198.1.
Fan, Yalin, and Stephen M Griffies, June 2014: Impacts of parameterized Langmuir turbulence and non-breaking wave mixing in global climate simulations. Journal of Climate, 27(12), DOI:10.1175/JCLI-D-13-00583.1. Abstract
We assess the impacts of parameterized upper ocean wave mixing on global climate simulations through modification to the K-profile ocean boundary layer parameterization (KPP; Large et al 1994) in a coupled atmosphere-ocean-wave global climate model. We consider three parameterizations and focus on impacts to high latitude ocean mixed layer depths and related ocean diagnostics. The McWilliams and Sullivan (2000) parameterization (MS2000) adds a Langmuir turbulence enhancement to the non-local component of KPP. We find that the Langmuir turbulence induced mixing provided by this parameterization is too strong in winter, producing overly deep mixed layers, and of minimal impact in summer. The Smyth et al (2002) parameterization modifies MS2000 by adding a stratification effect to restrain the turbulence enhancement under weak stratification conditions (e.g., winter) and to magnify the enhancement under strong stratification conditions. The Smyth et al (2002) scheme improves the simulated winter Mixed Layer Depth in our simulations, with mixed layer deepening in the Labrador Sea and shoaling in the Weddell and Ross Seas. Enhanced vertical mixing through parameterized Langmuir turbulence, coupled with enhanced lateral transport associated with parameterized mesoscale/submesoscale eddies, are found to be key elements for improving mixed layer simulations. Secondary impacts include strengthening the Atlantic Meridional Overturning Circulation and reducing the Antarctic Circumpolar Current. The Qiao et al (2004) non-breaking wave parameterization is the third scheme we assess. It adds a wave orbital velocity to the Reynolds stress calculation, and provides the strongest summer mixed layer deepening in the Southern Ocean among the three experiments, but with weak impacts during winter.
Farneti, Riccardo, S Dwivedi, Fred Kucharski, F Molteni, and Stephen M Griffies, September 2014: On Pacific Subtropical Cell Variability over the Second Half of the Twentieth Century. Journal of Climate, 27(18), DOI:10.1175/JCLI-D-13-00707.1. Abstract
The evolution of the Pacific subtropical cells (STC) is presented for the period 1948-2007. Using ocean models of different resolutions forced with interannually varying atmospheric forcing data sets, the mechanisms responsible for the observed STC weakening and late recovery during the period of study are analyzed. As a result of the STC weakening (strengthening), warming (cooling) trends are found in the equatorial Pacific sea surface temperatures (SSTs). Model results agree well with observed estimates of STC transport, convergence, and equatorial SST anomalies. It is shown that subtropical atmospheric variability is the primary driver of the STC and equatorial SST low-frequency evolution, and responsible for both the slowdown during the second half of the twentieth century and the rebound at the end of the century. Subtropically-forced STC variability is identified as a major player in the generation of equatorial Pacific decadal SST anomalies, pacing tropical Pacific natural climate variability on inter-decadal time scales, as observed in historical records. The natural mode of variability has implications for the evolution of equatorial SST in the coming decades under the concomitant effects of climate change.
We provide an assessment of sea level simulated in a suite of global ocean-sea ice models using the interannual CORE atmospheric state to determine surface ocean boundary buoyancy and momentum fluxes. These CORE-II simulations are compared amongst themselves as well as to observation-based estimates. We focus on the final 15 years of the simulations (1993-2007), as this is a period where the CORE-II atmospheric state is well sampled, and it allows us to compare sea level related fields to both satellite and in situ analyses. The ensemble mean of the CORE-II simulations broadly agree with various global and regional observation-based analyses during this period, though with the global mean thermosteric sea level rise biased low relative to observation-based analyses. The simulations reveal a positive trend in dynamic sea level in the west Pacific and negative trend in the east, with this trend arising from wind shifts and regional changes in upper 700 m ocean heat content. The models also exhibit a thermosteric sea level rise in the subpolar North Atlantic associated with a transition around 1995/1996 of the North Atlantic Oscillation to its negative phase, and the advection of warm subtropical waters into the subpolar gyre. Sea level trends are predominantly associated with steric trends, with thermosteric effects generally far larger than halosteric effects, except in the Arctic and North Atlantic. There is a general anti-correlation between thermosteric and halosteric effects for much of the World Ocean, associated with density compensated changes.
The small-slope approximation to the full three-dimensional diffusion tensor of epineutral diffusion gives exactly the same tracer flux as the commonly-used projected non-orthogonal diffusive flux of layered ocean models and of theoretical studies. The epineutral diffusion achieved by this small-slope approximation is not exactly in the direction of the correct epineutral tracer gradient. That is, the use of the small-slope approximation leads to a very small flux of tracer in a direction in which there is no epineutral gradient of tracer. For (the tracer) temperature or salinity, the difference between the correct epineutral gradient and the small-slope approximation to it is proportional to neutral helicity.
We also make the point that small-scale turbulent mixing processes act to diffuse tracer isotropically (i.e. the same in each spatial direction) and hence it is strictly a misnomer to call this process “dianeutral diffusion” or “vertical diffusion”. This realization also has implications for the diffusion tensor.
Despite slow rates of ocean mixing, observational and modeling studies suggest that buoyancy is redistributed to all depths of the ocean on surprisingly short interannual to decadal time scales. The mechanisms responsible for this redistribution remain poorly understood. This work uses an Earth System Model to evaluate the global steady state ocean buoyancy (and related steric sea level) budget, its interannual variability, and its transient response to a doubling of CO2 over 70 years, with a focus on the deep ocean. At steady state, the simple view of vertical advective-diffusive balance for the deep ocean holds at low- to mid-latitudes. At higher latitudes, the balance depends on a myriad of additional terms, namely mesoscale and submesoscale advection, convection and overflows from marginal seas, and terms related to the nonlinear equation of state. These high-latitude processes rapidly communicate anomalies in surface buoyancy forcing to the deep ocean locally; the deep, high-latitude changes then influence the large-scale advection of buoyancy to create transient deep buoyancy anomalies at lower latitudes. Following a doubling of atmospheric carbon dioxide concentrations, the high latitude buoyancy sinks are suppressed by a slowdown in convection and reduced dense water formation. This change is accompanied by a slowing of both upper and lower cells of the global meridional overturning circulation, reducing the supply of dense water to low latitudes beneath the pycnocline and the commensurate flow of light waters to high latitudes above the pycnocline. By this mechanism, changes in high latitude buoyancy are communicated to the global deep ocean on relatively fast advective timescales.
We have developed a one-way nested Indian Ocean regional model. The model combines the National Oceanic and Atmospheric Administration (NOAA) Geophysical Fluid Dynamics Laboratory’s (GFDL) Modular Ocean Model (MOM4p1) at global climate model resolution (nominally one degree), and a regional Indian Ocean MOM4p1 configuration with 25 km horizontal resolution and 1 meter vertical resolution near the surface. Inter-annual global simulations with Coordinated Ocean-Ice Reference Experiments (CORE-II) surface forcing over years 1992-2005 provide surface boundary conditions. We show that relative to the global simulation, (i) biases in upper ocean temperature, salinity and mixed layer depth are reduced, (ii) sea surface height and upper ocean circulation are closer to observations, and (iii) improvements in model simulation can be attributed to refined resolution, more realistic topography and inclusion of seasonal river runoff. Notably, the surface salinity bias is reduced to less than 0.1 psu over the Bay of Bengal using relatively weak restoring to observations, and the model simulates the strong, shallow halocline often observed in the north Bay of Bengal. There is marked improvement in subsurface salinity and temperature, as well as mixed layer depth in the Bay of Bengal. Major seasonal signatures in observed sea surface height anomaly in the tropical Indian Ocean, including the coastal waveguide around the Indian peninsula, are simulated with great fidelity. The use of realistic topography and seasonal river runoff brings the three dimensional structure of the East India Coastal Current and West India Coastal Current much closer to observations. As a result, the incursion of low salinity Bay of Bengal water into the south-eastern Arabian Sea is more realistic.
Spence, P, and Stephen M Griffies, et al., July 2014: Rapid subsurface warming and circulation changes of Antarctic coastal waters by poleward shifting winds. Geophysical Research Letters, 41(13), DOI:10.1002/2014GL060613. Abstract
The southern hemisphere westerly winds have been strengthening and shifting poleward since the 1950s. This wind trend is projected to persist under continued anthropogenic forcing, but the impact of the changing winds on Antarctic coastal heat distribution remains poorly understood. Here we show that a poleward wind shift at the latitudes of the Antarctic Peninsula can produce an intense warming of subsurface coastal waters that exceeds 2 °C at 200-700 m depth. The model simulated warming results from a rapid advective heat flux induced by weakened near-shore Ekman pumping, and is associated with weakened coastal currents. This analysis shows that anthropogenically induced wind changes can dramatically increase the temperature of ocean water at ice sheet grounding lines and at the base of floating ice shelves around Antarctica, with potentially significant ramifications for global sea level rise.
We investigate the influence of ocean component resolution on simulation of climate sensitivity using variants of the GFDL CM2.5 climate model incorporating eddy-resolving (1/10o) and eddy-parameterizing (1o) ocean resolutions. Two parameterization configurations of the coarse-resolution model are used yielding a three-model suite with significant variation in the transient climate response (TCR). The variation of TCR in this suite and in an enhanced group of 10 GFDL models is found to be strongly associated with the control climate Atlantic meridional overturning circulation (AMOC) magnitude and its decline under forcing. We find it is the AMOC behavior rather than resolution per se that accounts for most of the TCR differences. A smaller difference in TCR stems from the eddy-resolving model having more Southern Ocean surface warming than the coarse models.
Bi, D, and Stephen M Griffies, et al., March 2013: The ACCESS coupled model: description, control climate and evaluation. Australian Meteorological and Oceanographic Journal, 63(1), 41-64. Abstract
The Australian Community Climate and Earth System Simulator coupled model
(ACCESS-CM) has been developed at the Centre for Australian Weather and Climate Research (CAWCR), a partnership between CSIRO1
and the Bureau of Meteorology. It is built by coupling the UK Met Office atmospheric unified model (UM),
and other sub-models as required, to the ACCESS ocean model, which consists of
the NOAA/GFDL2 ocean model MOM4p1 and the LANL3 sea-ice model CICE4.1,
under the CERFACS4 OASIS3.2–5 coupling framework. The primary goal of the
ACCESS-CM development is to provide the Australian climate community with
a new generation fully coupled climate model for climate research, and to participate in phase five of the Coupled Model Inter-comparison Project (CMIP5).
This paper describes the ACCESS-CM framework and components, and presents
the control climates from two versions of the ACCESS-CM, ACCESS1.0 and ACCESS1.3, together with some fields from the 20th
century historical experiments,
as part of model evaluation. While sharing the same ocean sea-ice model (except
different setups for a few parameters), ACCESS1.0 and ACCESS1.3 differ from
each other in their atmospheric and land surface components: the former is configured with the UK Met Office HadGEM2 (r1.1) atmospheric physics and the Met
Office Surface Exchange Scheme land surface model version 2, and the latter with
atmospheric physics similar to the UK Met Office Global Atmosphere 1.0 including modifications performed at CAWCR and the CSIRO Community Atmosphere
Biosphere Land Exchange land surface model version 1.8. The global average
annual mean surface air temperature across the 500-year preindustrial control
integrations show a warming drift of 0.35 °C in ACCESS1.0 and 0.04 °C in ACCESS1.3. The overall skills of ACCESS-CM in simulating a set of key climatic fields
both globally and over Australia significantly surpass those from the preceding
CSIRO Mk3.5 model delivered to the previous coupled model inter-comparison.
However, ACCESS-CM, like other CMIP5 models, has deficiencies in various aspects, and these are also discussed.
Bi, D, and Stephen M Griffies, et al., March 2013: ACCESS-OM: the ocean and sea-ice core of the ACCESS coupled model. Australian Meteorological and Oceanographic Journal, 63(1), 213-232. Abstract
The Australian Community Climate and Earth System Simulator Ocean Model
(ACCESS-OM), a global coupled ocean and sea-ice model, has been developed
at the Centre for Australian Weather and Climate Research 1. It is aimed to serve
the Australian climate sciences community, including the Bureau of Meteorology,
CSIRO2 and Australian universities, for ocean climate research. ACCESS-OM
comprises the NOAA/GFDL3
Modular Ocean Model version 4p1; the LANL4
Sea-ice Model version 4.1, a data atmospheric model; and the CERFACS5 OASIS3.25
coupler, which constrains data exchange between the sub-models. ACCESS-OM
has been functioning as the ocean and sea-ice coupling core of the ACCESS coupled model, one of the Australian models participating in the Coupled Model Inter-comparison Project phase 5. This paper describes the ACCESS-OM sub-models, coupler, coupling strategy and framework. A selection of key metrics from
an ACCESS-OM benchmark simulation, which has run for 500 years using the
Coordinated Ocean-ice Reference Experiments normal year forcing, is presented
and compared with observations to evaluate the model performance. It shows
ACCESS-OM simulates the global ocean and sea-ice climate generally comparably to the results from other ocean sea-ice models of the same class (Griffies et al.
2009). For example, the global ocean volume-averaged temperature undergoes
minor evolution. The maximum transport of North Atlantic overturning circulation is 18.5 Sv and the Antarctic Circumpolar Current transport through Drake
Passage is 150 Sv, both in fair agreement with the observations; and the sea-ice
coverage has reasonable distribution and annual cycle. Measured against other
ocean sea-ice models and observations, ACCESS-OM is an appropriate tool for
Australia’s future ocean climate modelling efforts.
Bony, Sandrine, Isaac M Held, and Stephen M Griffies, et al., September 2013: Carbon Dioxide and Climate: Perspectives on a Scientific Assessment In Climate Science for Serving Society Research, Modeling and Prediction Priorities, DOI:10.1007/978-94-007-6692-1_14. Abstract
Many of the findings of the Charney Report on CO2-induced climate change published in 1979 are still valid, even after 30 additional years of climate research and observations. This paper considers the reasons why the report was so prescient, and assesses the progress achieved since its publication. We suggest that emphasis on the importance of physical understanding gained through the use of theory and simple models, both in isolation and as an aid in the interpretation of the results of General Circulation Models, provided much of the authors’ insight at the time. Increased emphasis on these aspects of research is likely to continue to be productive in the future, and even to constitute one of the most efficient routes towards improved climate change assessments.
We describe carbon system formulation and simulation characteristics of two new global coupled carbon-climate Earth System Models, ESM2M and ESM2G. These models demonstrate good climate fidelity as described in Part I while incorporating explicit and consistent carbon dynamics. The two models differ almost exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4.1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. On land, both ESMs include a revised land model to simulate competitive vegetation distributions and functioning, including carbon cycling among vegetation, soil and atmosphere. In the ocean, both models include new biogeochemical algorithms including phytoplankton functional group dynamics with flexible stoichiometry. Preindustrial simulations are spun up to give stable, realistic carbon cycle means and variability. Significant differences in simulation characteristics of these two models are described. Due to differences in oceanic ventilation rates (Part I) ESM2M has a stronger biological carbon pump but weaker northward implied atmospheric CO2 transport than ESM2G. The major advantages of ESM2G over ESM2M are: improved representation of surface chlorophyll in the Atlantic and Indian Oceans and thermocline nutrients and oxygen in the North Pacific. Improved tree mortality parameters in ESM2G produced more realistic carbon accumulation in vegetation pools. The major advantages of ESM2M over ESM2G are reduced nutrient and oxygen biases in the Southern and Tropical Oceans.
Frankcombe, L M., P Spence, Andrew McC Hogg, Matthew H England, and Stephen M Griffies, November 2013: Sea level changes forced by Southern Ocean winds. Geophysical Research Letters, 40(21), DOI:10.1002/2013GL058104. Abstract
On regional scales, changes in sea level are significantly affected by local dynamical changes. Westerly winds over the Southern Ocean have been strengthening and shifting southward in recent decades, and this change is projected to continue in the future. This study applies wind forcing anomalies to an eddy-permitting ocean model to study the dynamical response to a Southern Hemisphere westerly wind increase and/or southward shift. It is shown that the applied wind anomalies result in a change in sea surface slope across the Antarctic Circumpolar Current such that a fall in sea level occurs around the Antarctic continental margin. The Antarctic Circumpolar Current transport and regional sea level are particularly sensitive to latitudinal shifts in the wind, with a much more muted response found when only wind strengthening is applied. In addition to the local sea level changes, Southern Ocean winds also have a global effect through changing ocean heat content and the global overturning circulation.
Griffies, Stephen M., and Anne Marie Treguier, 2013: Ocean Circulation Models and Modeling In Ocean Circulation and Climate: A 21st Century Perspective, Second Edition, International Geophysics Vol 103, Elsevier, 521-551. Abstract
This chapter focuses on numerical models used to understand and predict large-scale circulation, such as the circulation comprising basin and global scales. It is organized according to two themes. The first addresses physical and numerical topics forming a foundation for ocean models. We focus here on the science of ocean models, in which we ask questions about fundamental processes and develop the mathematical equations for ocean thermo-hydrodynamics. We also touch upon various methods used to represent the continuum ocean fluid with a discrete computer model, raising such topics as the finite volume formulation of the ocean equations; the choice for vertical coordinate; the complementary issues related to horizontal gridding; and the pervasive questions of subgrid scale parameterizations. The second theme of this chapter concerns the applications of ocean models, in particular how to design an experiment and how to analyze results. This material forms the basis for ocean modelling, with the aim being to mechanistically describe, interpret, understand, and predict emergent features of the simulated, and ultimately the observed, ocean.
The impact of climate warming on the upper layer of the Bering Sea is investigated by using a high-resolution coupled global climate model. The model is forced by increasing atmospheric CO2 at a rate of 1% per year until CO2 reaches double its initial value (after 70 years), after which it is held constant. In response to this forcing, the upper layer of the Bering Sea warms by about 2�C in the southeastern shelf and by a little more than 1�C in the western basin. The wintertime ventilation to the permanent thermocline weakens in the western Bering Sea. After CO2 doubling, the southeastern shelf of the Bering Sea becomes almost ice-free in March, and the stratification of the upper layer strengthens in May and June. Changes of physical condition due to the climate warming would impact the pre-condition of spring bio-productivity in the southeastern shelf.
Marsland, S J., and Stephen M Griffies, et al., March 2013: Evaluation of ACCESS climate model ocean metrics in CMIP5 simulations. Australian Meteorological and Oceanographic Journal, 63(1), 101-119. Abstract
Global and regional diagnostics are used to evaluate the ocean performance of
the Australian Community Climate and Earth System Simulator coupled model
(ACCESS-CM) contributions to the Climate Model Intercomparison Project phase
5 (CMIP5). Two versions of ACCESS-CM have been submitted to CMIP; namely
CSIRO-BOM ACCESS1.0 and CSIRO-BOM ACCESS1.3. Results from six of the
core CMIP5 experiments (piControl, historical, rcp45, rcp85, 1pctCO2, and abrupt 4xCO2) are evaluated for each of the two ACCESS-CM model versions. Overall,
both model versions exhibit a reasonable and stable representation of key diagnostics of ocean climate performance in the pre-industrial control simulations,
including a meridional overturning circulation with North Atlantic Deep Water maxima in the range 22–24 Sv, and a poleward heat transport maximum of around
1.5 PW. For the projected climate change scenarios considered the ACCESS-CM
results are in reasonable agreement with responses found in other CMIP models,
with the familiar ocean warming, and reduction in strength of meridional overturning and poleward heat transport. Drifts in the control simulations of both
global ocean salinity and global sea-level are opposite in sign for ACCESS1.0 and
ACCESS1.3, suggesting problems exist in the closure of the hydrological cycle.
The simulation of ocean climate change over the historical period shows a weak
response compared to observations, which manifests as a late response of ocean
warming and sea level rise starting around 1990 in the model, compared to the
mid 1960s in observations. Further historical simulations are underway to ascertain if this late response in ACCESS is a robust model feature, or just low frequency variability. If the weak response over the historical period proves robust,
the likely cause is a too strong cooling from atmospheric aerosols. Broadening
the set of experiments to further investigate the relative warming response of the
ACCESS-CM to greenhouse gases compared to the cooling response to aerosols
is underway, and preliminary results do suggest that the cooling due to aerosols is
strong in the historical simulations.
Climate models simulate a wide range of climate changes at high northern latitudes in response to increased CO2. They also have substantial disagreement on projected changes of the Atlantic meridional overturning circulation (AMOC). Here we use two pairs of closely related climate models - each containing members with large and small AMOC declines - to explore the influence of AMOC decline on the high latitude response to increased CO2. The models with larger AMOC decline have less high latitude warming and sea ice decline than their small AMOC decline counterpart. By examining differences in the perturbation heat budget of the 40�90�N region, it is shown that AMOC decline diminishes the warming by weakening poleward ocean heat transport and increasing the ocean heat uptake. The cooling impact of this AMOC forced surface heat flux perturbation difference is enhanced by shortwave feedback and diminished by longwave feedback and atmospheric heat transport differences. The magnitude of the AMOC decline within model pairs is positively related to the magnitudes of control climate AMOC and Labrador Sea convection. Because the 40degree 90degree N region accounts for up to 40% of the simulated global ocean heat uptake over one hundred years, the process described here influences the global heat uptake efficiency.
Siedler, G A., Stephen M Griffies, J Gould, and J A Church, December 2013: Ocean Circulation and Climate: A 21st Century Perspective , Second Edition, International Geophysics V. 103: Academic Press, 904pp.
Straneo, F, Patrick Heimbach, Olga V Sergienko, G Hamilton, G Catania, Stephen M Griffies, and Robert Hallberg, et al., August 2013: Challenges to Understand the Dynamic Response of Greenland's Marine Terminating Glaciers to Oceanic and Atmospheric Forcing. Bulletin of the American Meteorological Society, 94(8), DOI:10.1175/BAMS-D-12-00100.1. Abstract
The recent retreat and speedup of outlet glaciers, as well as enhanced surface melting around the ice sheet margin, have increased Greenland's contribution to sea level rise to 0.6±0.1 mm/yr and its discharge of freshwater into the North Atlantic. The widespread, near-synchronous glacier retreat, and its coincidence with a period of oceanic and atmospheric warming, suggest a common climate driver. Evidence points to the marine margins of these glaciers as the region from which changes propagated inland. Yet the forcings and mechanisms behind these dynamic responses are poorly understood and either missing or crudely parameterized in climate and ice sheet models. Resulting projected sea level rise contributions from Greenland by 2100 remain highly uncertain.
This paper summarizes current state of knowledge and highlights key physical aspects of Greenland's coupled ice-sheet/ocean/atmosphere system. Three research thrusts are identified to yield fundamental insights into ice sheet, ocean, sea ice and atmosphere interactions, their role in Earth's climate system, and probable trajectories of future changes: (1) focused process studies addressing critical glacier, ocean, atmosphere and coupled dynamics; (2) sustained observations at key sites; and (3) inclusion of relevant dynamics in Earth System Models.
Understanding the dynamic response of Greenland's glaciers to climate forcing constitutes both a scientific and technological frontier given the challenges of obtaining the appropriate measurements from the glaciers' marine termini and the complexity of the dynamics involved, including the coupling of the ocean, atmosphere, glacier and sea ice systems. Interdisciplinary and international cooperation are crucial to making progress on this novel and complex problem.
Capsule: An interdisciplinary and multi-faceted approach is needed to understand the forcings and mechanisms behind the recent retreat and acceleration of Greenland's glaciers and its implications for future sea level rise
We examine the influence of alternative ocean and atmosphere subcomponents on climate model simulation of transient sensitivities by comparing three GFDL climate models used for the CMIP5. The base model ESM2M is closely related to GFDL's CMIP3 climate model CM2.1, and makes use of a depth coordinate ocean component. The second model, ESM2G, is identical to ESM2M but makes use of an isopycnal coordinate ocean model. We compare the impact of this "ocean swap" with an "atmosphere swap" that produces the CM3 climate model by replacing the AM2 atmosphere with AM3 while retaining a depth coordinate ocean model. The atmosphere swap is found to have much larger influence on sensitivities of global surface temperature and Northern Hemisphere sea ice cover. The atmosphere swap also introduces a multi-decadal response timescale through its indirect influence on heat uptake. Despite significant differences in their interior ocean mean states, the ESM2M and ESM2G simulations of these metrics of climate change are very similar, except for an enhanced high latitude salinity response accompanied by temporarily advancing sea ice in ESM2G. In the ESM2G historical simulation this behavior results in the establishment of a strong halocline in the subpolar North Atlantic during the early 20th century and an associated cooling which are counter to observations in that region. The Atlantic meridional overturning declines comparably in all three models.
The influence of changing ocean currents on climate change is evaluated by comparing an earth system model’s response to increased CO2 with and without an ocean circulation response. Inhibiting the ocean circulation response, by specifying a seasonally-varying preindustrial climatology of currents, has a much larger influence on the heat storage pattern than on the carbon storage pattern. The heat storage pattern without circulation changes resembles carbon storage (either with or without circulation changes) more than it resembles the heat storage when currents are allowed to respond. This is shown to be due to the larger magnitude of the redistribution transport – the change in transport due to circulation anomalies acting on control climate gradients – for heat than for carbon. The net ocean heat and carbon uptake are slightly reduced when currents are allowed to respond. Hence, ocean circulation changes potentially act to warm the surface climate. However, the impact of the reduced carbon uptake on radiative forcing is estimated to be small while the redistribution heat transport shifts ocean heat uptake from low to high latitudes increasing its cooling power. Consequently, global surface warming is significantly reduced by circulation changes. Circulation changes also shift the pattern of warming from broad northern hemisphere amplification to a more structured pattern with reduced warming at subpolar latitudes in both hemispheres and enhanced warming near the equator.
Bates, M L., Stephen M Griffies, and Matthew H England, December 2012: A dynamic, embedded Lagrangian model for ocean climate models, Part I: Theory and implementation. Ocean Modelling, 59-60, DOI:10.1016/j.ocemod.2012.05.004. Abstract
A framework for embedding a Lagrangian model within ocean climate models that employ horizontal Eulerian grids is presented. The embedded Lagrangian model can be used to explicitly represent processes that are at the subgrid scale to the Eulerian model. The framework is applied to open ocean deep convection and gravity driven downslope flows, both of which are subgridscale in the present generation of level coordinate ocean climate models. In order to apply the embedded Lagrangian framework to these processes, it is necessary to partition the mass and momentum of the model into an Eulerian system and a Lagrangian system. This partitioning allows the Lagrangian model to transport seawater using a more appropriate set of dynamics.
A number of schemes suitable for implementation in the embedded Lagrangian model are derived. Two dynamically passive schemes are derived that emulate existing parameterisations and two dynamically active schemes are also derived that evolve Lagrangian parcels of water (termed ”blobs”) according to a set of physical equations. Some details of the implementation into the Geophysical Fluid Dynamics Modular Ocean Model are also given. Finally, results are presented that show that the dynamically passive schemes are able to emulate their Eulerian counterparts to within roundoff error in idealised test cases.
Bates, M L., Stephen M Griffies, and Matthew H England, December 2012: A dynamic, embedded Lagrangian model for ocean climate models, Part II: Idealised overflow tests. Ocean Modelling, 59-60, DOI:10.1016/j.ocemod.2012.08.003. Abstract
Dense gravity current overflows occur in several regions throughout the world and are an important process in the meridional overturning circulation. Overflows are poorly represented in coarse resolution level coordinate ocean climate models. Here, the embedded Lagrangian model formulated in the companion paper of Bates et al., 2012 is used in two idealised test cases to examine the effect on the representation of dense gravity driven plumes, as well as the effect on the circulation of the bulk ocean in the Eulerian model. The results are compared with simulations with no parameterisation for overflows, as well as simulations that use traditional hydrostatic overflow schemes.
The use of Lagrangian “blobs” is shown to improve three key characteristics that are poorly represented in coarse resolution level coordinate models: (1) the depth of the plume, (2) the along slope velocity of the plume, and (3) the response of the bulk ocean to the bottom boundary layer. These improvements are associated with the more appropriate set of dynamics satisfied by the blobs, leading to a more physically sound representation. Experiments are also conducted to examine sensitivity to blob parameters. The blob parameters are examined over a large parameter space.
We present results for simulated climate and climate change from a newly developed high-resolution global climate model (GFDL CM2.5). The GFDL CM2.5 model has an atmospheric resolution of approximately 50 Km in the horizontal, with 32 vertical levels. The horizontal resolution in the ocean ranges from 28 Km in the tropics to 8 Km at high latitudes, with 50 vertical levels. This resolution allows the explicit simulation of some mesoscale eddies in the ocean, particularly at lower latitudes.
We present analyses based on the output of a 280 year control simulation; we also present results based on a 140 year simulation in which atmospheric CO2 increases at 1% per year until doubling after 70 years.
Results are compared to the GFDL CM2.1 climate model, which has somewhat similar physics but coarser resolution. The simulated climate in CM2.5 shows marked improvement over many regions, especially the tropics, including a reduction in the double ITCZ and an improved simulation of ENSO. Regional precipitation features are much improved. The Indian monsoon and Amazonian rainfall are also substantially more realistic in CM2.5.
The response of CM2.5 to a doubling of atmospheric CO2 has many features in common with CM2.1, with some notable differences. For example, rainfall changes over the Mediterranean appear to be tightly linked to topography in CM2.5, in contrast to CM2.1 where the response is more spatially homogeneous. In addition, in CM2.5 the near-surface ocean warms substantially in the high latitudes of the Southern Ocean, in contrast to simulations using CM2.1.
We describe the physical climate formulation and simulation characteristics of two new global coupled carbon-climate Earth System Models, ESM2M and ESM2G. These models demonstrate similar climate fidelity as the Geophysical Fluid Dynamics Laboratory’s previous CM2.1 climate model while incorporating explicit and consistent carbon dynamics. The two models differ exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4.1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. Differences in the ocean mean state include the thermocline depth being relatively deep in ESM2M and relatively shallow in ESM2G compared to observations. The crucial role of ocean dynamics on climate variability is highlighted in the El Niño-Southern Oscillation being overly strong in ESM2M and overly weak ESM2G relative to observations. Thus, while ESM2G might better represent climate changes relating to: total heat content variability given its lack of long term drift, gyre circulation and ventilation in the North Pacific, tropical Atlantic and Indian Oceans, and depth structure in the overturning and abyssal flows, ESM2M might better represent climate changes relating to: surface circulation given its superior surface temperature, salinity and height patterns, tropical Pacific circulation and variability, and Southern Ocean dynamics. Our overall assessment is that neither model is fundamentally superior to the other, and that both models achieve sufficient fidelity to allow meaningful climate and earth system modeling applications. This affords us the ability to assess the role of ocean configuration on earth system interactions in the context of two state-of-the-art coupled carbon-climate models.
This paper develops an analysis framework to identify how physical processes, as represented in ocean climate models, impact the evolution of global mean sea level. The formulation utilizes the coarse grained equations appropriate for an ocean model, and starts from the vertically integrated mass conservation equation in its Lagrangian form. Global integration of this kinematic equation results in an evolution equation for global mean sea level that depends on two physical processes: boundary fluxes of mass and the non-Boussinesq steric effect. The non-Boussinesq steric effect itself contains contributions from boundary fluxes of buoyancy; interior buoyancy changes associated with parameterized subgrid scale processes; and motion across pressure surfaces. The non-Boussinesq steric effect can be diagnosed in either volume conserving Boussinesq or mass conserving non-Boussinesq ocean circulation models, with differences found to be negligible. We find that surface heating is the dominant term affecting sea level arising from buoyancy fluxes, contributing to a net positive tendency to global mean sea level, largely due to low latitude heating and because the thermal expansion coefficient is much larger in the tropics than high latitudes. Subgrid scale effects from parameterized quasi-Stokes transport, vertical diffusion, cabbeling, and thermobaricity are also found to be significant, each resulting in a reduction of global mean sea level. Sea level rise through low latitude heating is largely compensated by a sea level drop from poleward eddy heat transport and ocean mixing. Spatial variations in the thermal expansion coefficient provide an essential modulation of how physical effects from mixing and eddy induced advective transport impact global mean sea level.
This paper examines spurious dianeutral transport within a suite of ocean models (GOLD, MITgcm, MOM, and ROMS). We quantify such transport through a global diagnostic that computes the reference potential energy, whose evolution arises solely through transport between density classes. Previous studies have focused on the importance of accurate tracer advection schemes in reducing the spurious transport and closure. The present study highlights complementary issues associated with momentum transport. Spurious dianeutral transport is shown to be directly proportional to the lateral grid Reynolds number (ReΔ), with such transport significantly reduced when ReΔ<10.
Simulations with the isopycnal model GOLD provide a benchmark for the smallest level of spurious dianeutral transport realizable in our model suite. For idealized simulations with a linear equation of state, GOLD exhibits identically zero spurious dianeutral mixing, and thus maintains a constant reference potential energy when all physical mixing processes are omitted. Amongst the non-isopycnal models tested in idealized simulations, ROMS generally produces smaller spurious dianeutral mixing than MITgcm or MOM, since ROMS makes use of a higher order upwind-biased scheme for momentum transport that enforces a small ReΔ. In contrast, MITgcm and MOM both employ unbiased (centered) discretizations of momentum transport, and therefore rely on lateral friction operators to control the grid Reynolds number. We find that a lateral shear-dependent Smagorinsky viscosity provides an effective means to locally reduce ReΔ, and thus to reduce spurious dianeutral transport in MITgcm and MOM.
In addition to four idealized simulations, we quantify spurious dianeutral transport in realistic global ocean climate simulations using GOLD and MOM with a realistic equation of state for seawater, both with and without mesoscale eddies in the resolved flow field. The GOLD simulations have detectable levels of spurious cabbeling from along isopycnal advective truncation errors. Significantly larger spurious dianeutral transport arises in a non-eddying MOM simulation. In an eddying MOM simulation, spurious dianeutral transport is larger still but is reduced by increasing momentum friction.
Lorbacher, K, S J Marsland, J A Church, Stephen M Griffies, and Detlef Stammer, June 2012: Rapid barotropic sea-level rise from ice-sheet melting. Journal of Geophysical Research: Oceans, 117, C06003, DOI:10.1029/2011JC007733. Abstract
Sea-level rise associated with idealized Greenland and Antarctic ice-sheet melting events is examined using a global coupled ocean sea-ice model that has a free surface formulation and thus can simulate fast barotropic motions. The perturbation experiments follow the Coordinated Ocean-ice Reference Experiment (CORE) version III. All regions of the global ocean experience a sea-level rise within 7-8 days of the initialization of a polar meltwater input of 0.1 Sv (1 Sv = 106 m3 s−1). The fast adjustment contrasts sharply with the slower adjustment associated with the smaller steric sea-level evolution that is also connected with melt events. The global mean sea-level rises by 9 mm yr−1 when this forcing is applied either from Greenland or Antarctica. Nevertheless, horizontal inter-basin gradients in sea level remain. For climate adaption in low-lying coastal and island regions, it is critical that the barotropic sea-level signal associated with melt events is taken into consideration, as it leads to a fast sea-level rise from melting ice-sheets for the bulk of the global ocean. A linear relation between sea-level rise and global meltwater input is further supported by experiments in which idealized melting occurs only in a region east or west of the Antarctic Peninsula, and when melting rates are varied between 0.01 Sv and 1.0 Sv. The results indicate that in ocean models that do not explicitly represent the barotropic signal, the barotropic component of sea-level rise can be added off-line to the simulated steric signal.
The Geophysical Fluid Dynamics Laboratory (GFDL) has developed a coupled general circulation model (CM3) for atmosphere, oceans, land, and sea ice. The goal of CM3 is to address emerging issues in climate change, including aerosol-cloud interactions, chemistry-climate interactions, and coupling between the troposphere and stratosphere. The model is also designed to serve as the physical-system component of earth-system models and models for decadal prediction in the near-term future, for example, through improved simulations in tropical land precipitation relative to earlier-generation GFDL models. This paper describes the dynamical core, physical parameterizations, and basic simulation characteristics of the atmospheric component (AM3) of this model.
Relative to GFDL AM2, AM3 includes new treatments of deep and shallow cumulus convection, cloud-droplet activation by aerosols, sub-grid variability of stratiform vertical velocities for droplet activation, and atmospheric chemistry driven by emissions with advective, convective, and turbulent transport. AM3 employs a cubed-sphere implementation of a finite-volume dynamical core and is coupled to LM3, a new land model with eco-system dynamics and hydrology.
Most basic circulation features in AM3 are simulated as realistically, or more so, than in AM2. In particular, dry biases have been reduced over South America. In coupled mode, the simulation of Arctic sea ice concentration has improved. AM3 aerosol optical depths, scattering properties, and surface clear-sky downward shortwave radiation are more realistic than in AM2. The simulation of marine stratocumulus decks and the intensity distributions of precipitation remain problematic, as in AM2.
The last two decades of the 20th century warm in CM3 by .32°C relative to 1881-1920. The Climate Research Unit (CRU) and Goddard Institute for Space Studies analyses of observations show warming of .56°C and .52°C, respectively, over this period. CM3 includes anthropogenic cooling by aerosol cloud interactions, and its warming by late 20th century is somewhat less realistic than in CM2.1, which warmed .66°C but did not include aerosol cloud interactions. The improved simulation of the direct aerosol effect (apparent in surface clear-sky downward radiation) in CM3 evidently acts in concert with its simulation of cloud-aerosol interactions to limit greenhouse gas warming in a way that is consistent with observed global temperature changes.
We estimate water mass transformation rates resulting from surface buoyancy fluxes and interior diapycnal fluxes in the region south of 30°S in the ECCO model based state estimation and three free-running coupled climate models. The meridional transport of deep and intermediate waters across 30°S agrees well between models and observationally based estimates in the Atlantic Ocean, but not in the Indian and Pacific where the model based estimates are much smaller. Associated with this, in the models about half the southward flowing deep water is converted into lighter waters and half to denser bottom waters, whereas the observationally-based estimates convert most of the inflowing deep water to bottom waters. In the models, both Antarctic Intermediate Water (AAIW) and Antarctic Bottom Water (AABW) are formed primarily via an interior diapycnal transformation rather than being transformed at the surface via heat or freshwater fluxes. Given the small vertical diffusivity specified in the models in this region, we conclude that other processes such as cabbeling and thermobaricity must be playing an important role in water mass transformation. Finally, in the models, the largest contribution of the surface buoyancy fluxes in the Southern Ocean is to convert Upper Circumpolar Deep Water (UCDW) and Antarctic Intermediate Water (AAIW) into lighter Sub-Antarctic Mode Water (SAMW) and Antarctic Intermediate Water (AAIW).
A parameterization for the restratification by finite-amplitude, submesoscale, mixed layer eddies, formulated as an overturning streamfunction, has been recently proposed to approximate eddy fluxes of density and other tracers. Here, the technicalities of implementing the parameterization in the coarse-resolution ocean component of global climate models are made explicit, and the primary impacts on model solutions of implementing the parameterization are discussed. Three global ocean general circulation models including this parameterization are contrasted with control simulations lacking the parameterization. The MLE parameterization behaves as expected and fairly consistently in models differing in discretization, boundary layer mixing, resolution, and other parameterizations. The primary impact of the parameterization is a shoaling of the mixed layer, with the largest effect in polar winter regions. Secondary impacts include strengthening the Atlantic meridional overturning while reducing its variability, reducing CFC and tracer ventilation, modest changes to sea surface temperature and air-sea fluxes, and an apparent reduction of sea ice basal melting.
The distribution of radiocarbon (14C) in the ocean and atmosphere has fluctuated on timescales ranging from seasons to millennia. It is thought that these fluctuations partly reflect variability in the climate system, offering a rich potential source of information to help understand mechanisms of past climate change. Here, a long simulation with a new, coupled model is used to explore the mechanisms that redistribute 14C within the Earth system on inter-annual to centennial timescales. The model, CM2Mc, is a lower-resolution version of the Geophysical Fluid Dynamics Laboratory's CM2M model, uses no flux adjustments, and incorporates a simple prognostic ocean biogeochemistry model including 14C. The atmospheric 14C and radiative boundary conditions are held constant, so that the oceanic distribution of 14C is only a function of internal climate variability. The simulation displays previously-described relationships between tropical sea surface 14C and the model-equivalents of the El Niño Southern Oscillation and Indonesian Throughflow. Sea surface 14C variability also arises from fluctuations in the circulations of the subarctic Pacific and Southern Ocean, including North Pacific decadal variability, and episodic ventilation events in the Weddell Sea that are reminiscent of the Weddell Polynya of 1974–1976. Interannual variability in the air-sea balance of 14C is dominated by exchange within the belt of intense Southern Westerly winds, rather than at the convective locations where the surface 14C is most variable. Despite significant interannual variability, the simulated impact on air-sea exchange is an order of magnitude smaller than the recorded atmospheric 14C variability of the past millennium. This result partly reflects the importance of variability in the production rate of 14C in determining atmospheric 14C, but may also reflect an underestimate of natural climate variability, particularly in the Southern Westerly winds.
This paper documents time mean simulation characteristics from the ocean and sea ice components in a new coupled climate model developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). The climate model, known as CM3, is formulated with effectively the same ocean and sea ice components as the earlier GFDL climate model, CM2.1, yet with extensive developments made to the atmosphere and land model components. Both CM2.1 and CM3 show stable mean climate indices, such as large scale circulation and sea surface temperatures (SSTs). There are notable improvements in the CM3 climate simulation relative to CM2.1, including a modified SST bias pattern and reduced biases in the Arctic sea ice cover. We anticipate SST differences between CM2.1 and CM3 in lower latitudes through analysis of the atmospheric fluxes at the ocean surface in corresponding Atmospheric Model Intercomparison Project (AMIP) simulations. In contrast, SST changes in the high latitudes are dominated by ocean and sea ice effects absent in AMIP simulations. The ocean interior simulation in CM3 is generally warmer than CM2.1, which adversely impacts the interior biases.
Griffies, Stephen M., and Gokhan Danabasoglu, May 2011: Physical ocean fields in CMIP5. Clivar Exchanges, 16(2), 32-34. PDF
Griffies, Stephen M., July 2011: Preface to the Ocean Modelling special issue on ocean eddies. Ocean Modelling, 39(1-2), 1. PDF
The applicability of Modular Ocean Model version 4 (MOM4p1) as a code base to study regional physical oceanographic phenomena is presented, highlighting features recently implemented for use in limited area domains. Central to the successful operation of limited area model applications are the inclusion of a comprehensive suite of open boundary conditions, turbulence closure and vertical discretization. The open boundary problem, in particular, is considered and we present the open boundary condition implementation and performance in limited area model configurations corresponding to three realistic test cases. These tests represent typical configurations the physical oceanographer may encounter, and consist of (1) a coastal shelf model application where a two-tiered model configuration is used for down-scaling from a coarse grid model to supply sufficiently resolved boundary values for active cross-shelf open boundaries of a regional model; (2) tidal response of a gulf with one open boundary across the mouth of the gulf; (3) response of a coastal region to the passage of a tropical cyclone, where the open boundaries behave in primarily a passive capacity. Although the code base used in the test cases is MOM4p1, emphasis is placed on general features of the tests that are necessary for the scientific and operational use of any limited area model, hence key findings may be applied to limited area models in general.
The study of climate impacts on Living Marine Resources (LMRs) has increased rapidly in recent years with the availability of climate model simulations contributed to the assessment reports of the Intergovernmental Panel on Climate Change (IPCC). Collaboration between climate and LMR scientists and shared understanding of critical challenges for such applications are essential for developing robust projections of climate impacts on LMRs. This paper assesses present approaches for generating projections of climate impacts on LMRs using IPCC-class climate models, recommends practices that should be followed for these applications, and identifies priority developments that could improve current projections. Understanding of the climate system and its representation within climate models has progressed to a point where many climate model outputs can now be used effectively to make LMR projections. However, uncertainty in climate model projections (particularly biases and inter-model spread at regional to local scales), coarse climate model resolution, and the uncertainty and potential complexity of the mechanisms underlying the response of LMRs to climate limit the robustness and precision of LMR projections. A variety of techniques including the analysis of multi-model ensembles, bias corrections, and statistical and dynamical downscaling can ameliorate some limitations, though the assumptions underlying these approaches and the sensitivity of results to their application must be assessed for each application. Developments in LMR science that could improve current projections of climate impacts on LMRs include improved understanding of the multi-scale mechanisms that link climate and LMRs and better representations of these mechanisms within more holistic LMR models. These developments require a strong baseline of field and laboratory observations including long time-series and measurements over the broad range of spatial and temporal scales over which LMRs and climate interact. Priority developments for IPCC-class climate models include improved model accuracy (particularly at regional and local scales), inter-annual to decadal-scale predictions, and the continued development of earth system models capable of simulating the evolution of both the physical climate system and biosphere. Efforts to address these issues should occur in parallel and be informed by the continued application of existing climate and LMR models.
Yin, Jianjun, J E Overland, Stephen M Griffies, Aixue Hu, Joellen L Russell, and Ronald J Stouffer, August 2011: Different magnitudes of projected subsurface ocean warming around Greenland and Antarctica. Nature Geoscience, 4(8), DOI:10.1038/ngeo1189. Abstract
The observed acceleration of outlet glaciers and ice flows in Greenland and Antarctica is closely linked to ocean warming, especially in the subsurface layer. Accurate projections of ice-sheet dynamics and global sea-level rise therefore require information of future ocean warming in the vicinity of the large ice sheets. Here we use a set of 19 state-of-the-art climate models to quantify this ocean warming in the next two centuries. We find that in response to a mid-range increase in atmospheric greenhouse-gas concentrations, the subsurface oceans surrounding the two polar ice sheets at depths of 200–500 m warm substantially compared with the observed changes thus far6, 7, 8. Model projections suggest that over the course of the twenty-first century, the maximum ocean warming around Greenland will be almost double the global mean, with a magnitude of 1.7–2.0 °C. By contrast, ocean warming around Antarctica will be only about half as large as global mean warming, with a magnitude of 0.5–0.6 °C. A more detailed evaluation indicates that ocean warming is controlled by different mechanisms around Greenland and Antarctica. We conclude that projected subsurface ocean warming could drive significant increases in ice-mass loss, and heighten the risk of future large sea-level rise.
Simulations from a fine-resolution global coupled model, the Geophysical Fluid Dynamics Laboratory
Climate Model, version 2.4 (CM2.4), are presented, and the results are compared with a coarse version of the
same coupled model, CM2.1, under idealized climate change scenarios. A particular focus is given to the
dynamical response of the Southern Ocean and the role played by the eddies—parameterized or permitted—
in setting the residual circulation and meridional density structure. Compared to the case in which eddies are
parameterized and consistent with recent observational and idealized modeling studies, the eddy-permitting
integrations of CM2.4 show that eddy activity is greatly energized with increasing mechanical and buoyancy
forcings, buffering the ocean to atmospheric changes, and the magnitude of the residual oceanic circulation
response is thus greatly reduced. Although compensation is far from being perfect, changes in poleward eddy
fluxes partially compensate for the enhanced equatorward Ekman transport, leading to weak modifications in
local isopycnal slopes, transport by the Antarctic Circumpolar Current, and overturning circulation. Since the
presence of active ocean eddy dynamics buffers the oceanic response to atmospheric changes, the associated
atmospheric response to those reduced ocean changes is also weakened. Further, it is hypothesized that
present numerical approaches for the parameterization of eddy-induced transports could be too restrictive
and prevent coarse-resolution models from faithfully representing the eddy response to variability and change
in the forcing fields.
We present a physically and numerically motivated boundary-value problem for each vertical ocean column, whose solution yields a parameterized mesoscale eddy-induced transport streamfunction. The new streamfunction is a nonlocal function of the properties of the fluid column. It is constructed to have a low baroclinic mode vertical structure and to smoothly transition through regions of weak stratification such as boundary layers or mode waters. It requires no matching conditions or regularization in unstratified regions; it satisfies boundary conditions of zero transport at the ocean surface and bottom; and it provides a sink of available potential energy for each vertical seawater column, but not necessarily at each location within the column. Numerical implementation of the methodology requires the solution of a one-dimensional tridiagonal problem for each vertical column. To illustrate the approach, we present an analytical example based on the nonlinear Eady problem and two numerical simulations.
We overview problems and prospects in ocean circulation models, with emphasis on certain developments aiming to
enhance the physical integrity and flexibility of large-scale models used to study global climate. We also consider elements
of observational measures rendering information to help evaluate simulations and to guide development priorities.
http://www.oceanobs09.net/blog/?p=88
Kopp, Robert E., J X Mitrovica, Stephen M Griffies, Jianjun Yin, C C Hay, and Ronald J Stouffer, et al., December 2010: The impact of Greenland melt on local sea levels: a partially coupled analysis of dynamic and static equilibrium effects in idealized water-hosing experiments. Climatic Change, 103(3-4), DOI:10.1007/s10584-010-9935-1. Abstract
Local sea level can deviate from mean global sea level because of both dynamic sea level (DSL) effects, resulting from oceanic and atmospheric circulation and temperature and salinity distributions, and changes in the static equilibrium (SE) sea level configuration, produced by the gravitational, elastic, and rotational effects of mass redistribution. Both effects will contribute to future sea level change. To compare their magnitude, we simulated the effects of Greenland Ice Sheet (GIS) melt by conducting idealized North Atlantic “water-hosing” experiments in a climate model unidirectionally coupled to a SE sea level model. At current rates of GIS melt, we find that geographic SE patterns should be challenging but possible to detect above dynamic variability. At higher melt rates, we find that DSL trends are strongest in the western North Atlantic, while SE effects will dominate in most of the ocean when melt exceeds ~20 cm equivalent sea level.
Rienecker, M M., Stephen M Griffies, and Anthony Rosati, et al., September 2010: Synthesis and Assimilation Systems: Essential Adjuncts to the Global Ocean Observing System In OceanObs’09: Sustained Ocean Observations and Information for Society, Vol. 2, ESA Publication, DOI:doi:10.5270/OceanObs09.pp.31.
The unphysical virtual salt flux (VSF) formulation widely used in the ocean component of climate models has the potential to cause systematic and significant biases in modeling the climate system and projecting its future evolution. Here a freshwater flux (FWF) and a virtual salt flux version of the Geophysical Fluid Dynamics Laboratory Climate Model version 2.1 (GFDL CM2.1) are used to evaluate and quantify the uncertainties induced by the VSF formulation. Both unforced and forced runs with the two model versions are performed and compared in detail. It is found that the differences between the two versions are generally small or statistically insignificant in the unforced control runs and in the runs with a small external forcing. In response to a large external forcing, however, some biases in the VSF version become significant, especially the responses of regional salinity and global sea level. However, many fundamental aspects of the responses differ only quantitatively between the two versions. An unexpected result is the distinctly different ENSO responses. Under a strong external freshwater forcing, the great enhancement of the ENSO variability simulated by the FWF version does not occur in the VSF version and is caused by the overexpansion of the top model layer. In summary, the principle assumption behind using virtual salt flux is not seriously violated and the VSF model has the ability to simulate the current climate and project near-term climate evolution. For some special studies such as a large hosing experiment, however, both the VSF formulation and the use of the FWF in the geopotential coordinate ocean model could have some deficiencies and one should be cautious to avoid them.
A set of state-of-the-science climate models are used to investigate global sea level rise (SLR) patterns
induced by ocean dynamics in twenty-first-century climate projections. The identified robust features include
bipolar and bihemisphere seesaws in the basin-wide SLR, dipole patterns in the North Atlantic and North
Pacific, and a beltlike pattern in the Southern Ocean. The physical and dynamical mechanisms that cause
these patterns are investigated in detail using version 2.1 of the Geophysical Fluid Dynamics Laboratory
(GFDL) Coupled Model (CM2.1). Under the Intergovernmental Panel on Climate Change’s (IPCC) Special
Report on Emissions Scenarios (SRES) A1B scenario, the steric sea level changes relative to the global mean
(the local part) in different ocean basins are attributed to differential heating and salinity changes of various
ocean layers and associated physical processes. As a result of these changes, water tends to move from the
ocean interior to continental shelves. In the North Atlantic, sea level rises north of the Gulf Stream but falls to
the south. The dipole pattern is induced by a weakening of the meridional overturning circulation. This
weakening leads to a local steric SLR east of North America, which drives more waters toward the shelf,
directly impacting northeastern North America. An opposite dipole occurs in the North Pacific. The dynamic
SLR east of Japan is linked to a strong steric effect in the upper ocean and a poleward expansion of the
subtropical gyre. In the Southern Ocean, the beltlike pattern is dominated by the baroclinic process during
the twenty-first century, while the barotropic response of sea level to wind stress anomalies is significantly
delayed.
Coordinated Ocean-ice Reference Experiments (COREs) are presented as a tool to explore the behaviour of global ocean-ice models under forcing from a common atmospheric dataset. We highlight issues arising when designing coupled global ocean and sea ice experiments, such as difficulties formulating a consistent forcing methodology and experimental protocol. Particular focus is given to the hydrological forcing, the details of which are key to realizing simulations with stable meridional overturning circulations.
The atmospheric forcing from [Large, W., Yeager, S., 2004. Diurnal to decadal global forcing for ocean and sea-ice models: the data sets and flux climatologies. NCAR Technical Note: NCAR/TN-460+STR. CGD Division of the National Center for Atmospheric Research] was developed for coupled-ocean and sea ice models. We found it to be suitable for our purposes, even though its evaluation originally focussed more on the ocean than on the sea-ice. Simulations with this atmospheric forcing are presented from seven global ocean-ice models using the CORE-I design (repeating annual cycle of atmospheric forcing for 500 years). These simulations test the hypothesis that global ocean-ice models run under the same atmospheric state produce qualitatively similar simulations. The validity of this hypothesis is shown to depend on the chosen diagnostic. The CORE simulations provide feedback to the fidelity of the atmospheric forcing and model configuration, with identification of biases promoting avenues for forcing dataset and/or model development.
Griffies, Stephen M., Alistair Adcroft, V Balaji, Robert Hallberg, Sonya Legg, Torge Martin, and Anna Pirani, et al., February 2009: Sampling Physical Ocean Field in WCRP CMIP5 Simulations: CLIVAR Working Group on Ocean Model Development (WGOMD) Committee on CMIP5 Ocean Model Output, International CLIVAR Project Office, CLIVAR Publication Series No. 137, 56pp. PDF
Legg, Sonya, Tal Ezer, Stephen M Griffies, Robert Hallberg, and L Jackson, et al., May 2009: Improving oceanic overflow representation in climate models: The gravity current entrainment climate process team. Bulletin of the American Meteorological Society, 90(5), DOI:10.1175/2008BAMS2667.1. Abstract
Oceanic overflows are bottom-trapped density currents originating in semienclosed basins, such as the Nordic seas, or on continental shelves, such as the Antarctic shelf. Overflows are the source of most of the abyssal waters, and therefore play an important role in the large-scale ocean circulation, forming a component of the sinking branch of the thermohaline circulation. As they descend the continental slope, overflows mix vigorously with the surrounding oceanic waters, changing their density and transport significantly. These mixing processes occur on spatial scales well below the resolution of ocean climate models, with the result that deep waters and deep western boundary currents are simulated poorly. The Gravity Current Entrainment Climate Process Team was established by the U.S. Climate Variability and Prediction (CLIVAR) Program to accelerate the development and implementation of improved representations of overflows within large-scale climate models, bringing together climate model developers with those conducting observational, numerical, and laboratory process studies of overflows. Here, the organization of the Climate Process Team is described, and a few of the successes and lessons learned during this collaboration are highlighted, with some emphasis on the well-observed Mediterranean overflow. The Climate Process Team has developed several different overflow parameterizations, which are examined in a hierarchy of ocean models, from comparatively well-resolved regional models to the largest-scale global climate models.
Griffies, Stephen M., and Alistair Adcroft, 2008: Formulating the equations of ocean models In Ocean Modeling in an Eddying Regime, Geophysical Monograph 177, M. W. Hecht, and H. Hasumi, eds., Washington, DC, American Geophysical Union, 281-318. Abstract PDF
We formulate mathematical equations describing the thermo-hydrodynamics of the ocean and introduce certain numerical methods employed by models used for ocean simulations.
Griffies, Stephen M., H Banks, and Anna Pirani, 2008: Furthering the science of ocean climate modelling. Clivar Exchanges, 13(1), 281-318. PDF
Pirani, Anna, Stephen M Griffies, and H Banks, 2008: Report from the CLIVAR Working Group on ocean model development (WGOMD). Clivar Exchanges, 13(1), 30-32. PDF
In
many global ocean climate models, mesoscale eddies are parameterized as
along isopycnal diffusion and eddy-induced advection (or equivalently
skew-diffusion). The eddy-induced advection flattens isopycnals and acts as
a sink of available potential energy, whereas the isopycnal diffusion mixes
tracers along neutral directions. While much effort has gone into estimating
diffusivities associated with this closure, less attention has been paid to
the details of how this closure (which tries to flatten isopycnals)
interacts with the mixed layer (in which vertical mixing tries to drive the
isopycnals vertical). In order to maintain numerical stability, models often
stipulate a maximum slope Smax which in combination with
the thickness diffusivity Agm defines a maximum
eddy-induced advective transport Agm*Smax.
In this paper, we examine the impact of changing Smax
within the GFDL global coupled climate model. We show that this parameter
produces significant changes in wintertime mixed layer depth, with
implications for wintertime temperatures in key regions, the distribution of
precipitation, and the vertical structure of heat uptake. Smaller changes
are seen in details of ventilation and currents, and even smaller changes as
regards the overall hydrography. The results suggest that not only the value
of the coefficient, but the details of the tapering scheme, need to be
considered when comparing isopycnal mixing schemes in models.
Griffies, Stephen M., C Böning, and Anne Marie Treguier, 2007: Design considerations for coordinated ocean-ice reference experiments. Flux News, 3, 3-5. PDF
The formulation and simulation characteristics of two new global coupled climate models developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) are described. The models were designed to simulate atmospheric and oceanic climate and variability from the diurnal time scale through multicentury climate change, given our computational constraints. In particular, an important goal was to use the same model for both experimental seasonal to interannual forecasting and the study of multicentury global climate change, and this goal has been achieved.
Two versions of the coupled model are described, called CM2.0 and CM2.1. The versions differ primarily in the dynamical core used in the atmospheric component, along with the cloud tuning and some details of the land and ocean components. For both coupled models, the resolution of the land and atmospheric components is 2° latitude × 2.5° longitude; the atmospheric model has 24 vertical levels. The ocean resolution is 1° in latitude and longitude, with meridional resolution equatorward of 30° becoming progressively finer, such that the meridional resolution is 1/3° at the equator. There are 50 vertical levels in the ocean, with 22 evenly spaced levels within the top 220 m. The ocean component has poles over North America and Eurasia to avoid polar filtering. Neither coupled model employs flux adjustments.
The control simulations have stable, realistic climates when integrated over multiple centuries. Both models have simulations of ENSO that are substantially improved relative to previous GFDL coupled models. The CM2.0 model has been further evaluated as an ENSO forecast model and has good skill (CM2.1 has not been evaluated as an ENSO forecast model). Generally reduced temperature and salinity biases exist in CM2.1 relative to CM2.0. These reductions are associated with 1) improved simulations of surface wind stress in CM2.1 and associated changes in oceanic gyre circulations; 2) changes in cloud tuning and the land model, both of which act to increase the net surface shortwave radiation in CM2.1, thereby reducing an overall cold bias present in CM2.0; and 3) a reduction of ocean lateral viscosity in the extratropics in CM2.1, which reduces sea ice biases in the North Atlantic.
Both models have been used to conduct a suite of climate change simulations for the 2007 Intergovernmental Panel on Climate Change (IPCC) assessment report and are able to simulate the main features of the observed warming of the twentieth century. The climate sensitivities of the CM2.0 and CM2.1 models are 2.9 and 3.4 K, respectively. These sensitivities are defined by coupling the atmospheric components of CM2.0 and CM2.1 to a slab ocean model and allowing the model to come into equilibrium with a doubling of atmospheric CO2. The output from a suite of integrations conducted with these models is freely available online (see http://nomads.gfdl.noaa.gov/).
Manuscript received 8 December 2004, in final form 18 March 2005
We study the reaction of a global ocean–sea ice model to an increase of fresh water input into the northern North Atlantic under different surface boundary conditions, ranging from simple restoring of surface salinity to the use of an energy balance model (EBM) for the atmosphere. The anomalous fresh water flux is distributed around Greenland, reflecting increased melting of the Greenland ice sheet and increasing fresh water export from the Arctic Ocean. Depending on the type of surface boundary condition, the large circulation reacts with a slow-down of overturning and gyre circulations. Restoring of the total or mean surface salinity prevents a large scale redistribution of the salinity field that is apparent under mixed boundary conditions and with the EBM. The control run under mixed boundary conditions exhibits large and unrealistic oscillations of the meridional overturning. Although the reaction to the fresh water flux anomaly is similar to the response with the EBM, mixed boundary conditions must thus be considered unreliable. With the EBM, the waters in the deep western boundary current initially become saltier and a new fresh water mass forms in the north-eastern North Atlantic in response to the fresh water flux anomaly around Greenland. After an accumulation period of several decades duration, this new North East Atlantic Intermediate Water spreads towards the western boundary and opens a new southward pathway at intermediate depths along the western boundary for the fresh waters of high northern latitudes.
The current generation of coupled climate models run at the Geophysical Fluid Dynamics Laboratory (GFDL) as part of the Climate Change Science Program contains ocean components that differ in almost every respect from those contained in previous generations of GFDL climate models. This paper summarizes the new physical features of the models and examines the simulations that they produce. Of the two new coupled climate model versions 2.1 (CM2.1) and 2.0 (CM2.0), the CM2.1 model represents a major improvement over CM2.0 in most of the major oceanic features examined, with strikingly lower drifts in hydrographic fields such as temperature and salinity, more realistic ventilation of the deep ocean, and currents that are closer to their observed values. Regional analysis of the differences between the models highlights the importance of wind stress in determining the circulation, particularly in the Southern Ocean. At present, major errors in both models are associated with Northern Hemisphere Mode Waters and outflows from overflows, particularly the Mediterranean Sea and Red Sea.
Jackett, D R., Trevor J McDougall, R Feistel, D Wright, and Stephen M Griffies, 2006: Algorithms for Density, Potential Temperature, Conservative Temperature, and the Freezing Temperature of Seawater. Journal of Atmospheric and Oceanic Technology, 23(12), DOI:10.1175/JTECH1946.1. Abstract
Algorithms are presented for density, potential temperature, conservative temperature, and the freezing temperature of seawater. The algorithms for potential temperature and density (in terms of potential temperature) are updates to routines recently published by McDougall et al., while the algorithms involving conservative temperature and the freezing temperatures of seawater are new. The McDougall et al. algorithms were based on the thermodynamic potential of Feistel and Hagen; the algorithms in this study are all based on the “new extended Gibbs thermodynamic potential of seawater” of Feistel. The algorithm for the computation of density in terms of salinity, pressure, and conservative temperature produces errors in density and in the corresponding thermal expansion coefficient of the same order as errors for the density equation using potential temperature, both being twice as accurate as the International Equation of State when compared with Feistel’s new equation of state. An inverse function relating potential temperature to conservative temperature is also provided. The difference between practical salinity and absolute salinity is discussed, and it is shown that the present practice of essentially ignoring the difference between these two different salinities is unlikely to cause significant errors in ocean models.
Griffies, Stephen M., 2005: Some ocean model fundamentals In Ocean Weather Forecasting: An Integrated View of Oceanography, Berlin, Germany, Springer, 19-74. Abstract
The purpose of these lectures is to present elements of the equations and algorithms used in numerical models of the large-scale ocean circulation. Such models generally integrate the ocean's primitive equations, which are based on Newton's Laws applied to a continuum fluid under hydrostatic balance in a spherical geometry, along with linear irreversible thermodynamics and subgrid scale (SGS) parameterizations. During formulations of both the kinematics and dynamics, we highlight issues related to the use of a generalized vertical coordinate. The vertical coordinate is arguably the most critical element determining how a model is designed and applications to which a model is of use.
This paper summarizes the formulation of the ocean component to the Geophysical Fluid Dynamics Laboratory's (GFDL) climate model used for the 4th IPCC Assessment (AR4) of global climate change. In particular, it reviews the numerical schemes and physical parameterizations that make up an ocean climate model and how these schemes are pieced together for use in a state-of-the-art climate model. Features of the model described here include the following: (1) tripolar grid to resolve the Arctic Ocean without polar filtering, (2) partial bottom step representation of topography to better represent topographically influenced advective and wave processes, (3) more accurate equation of state, (4) three-dimensional flux limited tracer advection to reduce overshoots and undershoots, (5) incorporation of regional climatological variability in shortwave penetration, (6) neutral physics parameterization for representation of the pathways of tracer transport, (7) staggered time stepping for tracer conservation and numerical efficiency, (8) anisotropic horizontal viscosities for representation of equatorial currents, (9) parameterization of exchange with marginal seas, (10) incorporation of a free surface that accommodates a dynamic ice model and wave propagation, (11) transport of water across the ocean free surface to eliminate unphysical "virtual tracer flux" methods, (12) parameterization of tidal mixing on continental shelves. We also present preliminary analyses of two particularly important sensitivities isolated during the development process, namely the details of how parameterized subgridscale eddies transport momentum and tracers.
The impact of changes in shortwave radiation penetration depth on the global ocean circulation and heat transport is studied using the GFDL Modular Ocean Model (MOM4) with two independent parameterizations that use ocean color to estimate the penetration depth of shortwave radiation. Ten to eighteen percent increases in the depth of 1% downwelling surface irradiance levels results in an increase in mixed layer depths of 3-20 m in the subtropical and tropical regions with no change at higher latitudes. While 1D models have predicted that sea surface temperatures at the equator would decrease with deeper penetration of solar irradiance, this study shows a warming, resulting in a 10% decrease in the required restoring heat flux needed to maintain climatological sea surface temperatures in the eastern equatorial Atlantic and Pacific Oceans. The decrease in the restoring heat flux is attributed to a slowdown in heat transport (5%) from the Tropics and an increase in the temperature of submixed layer waters being transported into the equatorial regions. Calculations were made using a simple relationship between mixed layer depth and meridional mass transport. When compared with model diagnostics, these calculations suggest that the slowdown in heat transport is primarily due to off-equatorial increases in mixed layer depths. At higher latitudes (5°-40°), higher restoring heat fluxes are needed to maintain sea surface temperatures because of deeper mixed layers and an increase in storage of heat below the mixed layer. This study offers a way to evaluate the changes in irradiance penetration depths in coupled ocean-atmosphere GCMs and the potential effect that large-scale changes in chlorophyll a concentrations will have on ocean circulation.
Griffies, Stephen M., 2004: Fundamentals of Ocean Climate Models, Princeton, NJ: Princeton University Press, 518 pp. Abstract
This book sets forth the physical, mathematical, and numerical foundations of computer models used to understand and predict the global ocean climate system. Aimed at students and researchers of ocean and climate science who seek to understand the physical content of ocean model equations and numerical methods for their solution, it is largely general in formulation and employs modern mathematical techniques. It also highlights certain areas of cutting-edge research.
Stephen Griffies presents material that spans a broad spectrum of issues critical for modern ocean climate models. Topics are organized into parts consisting of related chapters, with each part largely self-contained. Early chapters focus on the basic equations arising from classical mechanics and thermodynamics used to rationalize ocean fluid dynamics. These equations are then cast into a form appropriate for numerical models of finite grid resolution. Basic discretization methods are described for commonly used classes of ocean climate models. The book proceeds to focus on the parameterization of phenomena occurring at scales unresolved by the ocean model, which represents a large part of modern oceanographic research. The final part provides a tutorial on the tensor methods that are used throughout the book, in a general and elegant fashion, to formulate the equations.
This manual provides a detailed description of the analytical, numerical, and computational aspects of the MOM4 ocean model.
Griffies, Stephen M., 2003: An introduction to linear predictability analysis In Global Climate, Rodó, X., and F. A. Comín, eds., Berlin, Springer-Verlag, 80-101.
Griffies, Stephen M., 2003: An introduction to ocean climate modeling In Global Climate, Rodó, X., and F. A. Comín, eds., Berlin, Springer-Verlag, 55-79.
This paper details a free surface method using an explicit time stepping scheme for use in z-coordinate ocean models. One key property that makes the method especially suitable for climate simulations is its very stable numerical time stepping scheme, which allows for the use of a long density time step, as commonly employed with coarse-resolution rigid-lid models. Additionally, the effects of the undulating free surface height are directly incorporated into the baroclinic momentum and tracer equations. The novel issues related to local and global tracer conservation when allowing for the top cell to undulate are the focus of this work. The method presented here is quasi-conservative locally and globally of tracer when the baroclinic and tracer time steps are equal. Important issues relevant for using this method in regional as well as large-scale climate models are discussed and illustrated, and examples of scaling achieved on parallel computers provided.
Stocker, T F., Thomas L Delworth, Stephen M Griffies, Isaac M Held, V Ramaswamy, and Brian J Soden, et al., 2001: Physical climate processes and feedbacks In Climate Change 2001: The Scientific Basis. Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge, UK, Cambridge University Press, 418-470.
Griffies, Stephen M., 2000: Review of "Ocean Modeling and Parameterization," E. P. Chassignet and J. Verron, eds., 1998, Kluwer Academic Publishers. Bulletin of the American Meteorological Society, 81(3), 591-593.
Griffies, Stephen M., C Böning, F O Bryan, Eric P Chassignet, R Gerdes, H Hasumi, A C Hirst, Anne Marie Treguier, and D Webb, 2000: Developments in ocean climate modelling. Ocean Modelling, 2, 123-192. Abstract PDF
This paper presents some research developments in primitive equation ocean models which could impact the ocean component of realistic global climate models aimed at large-scale, low frequency climate simulations and predictions. It is written primarily to an audience of modellers concerned with the ocean component of climate models, although not necessarily experts in the design and implementation of ocean model algorithms.
Griffies, Stephen M., and Robert Hallberg, 2000: Biharmonic friction with a Smagorinsky-like viscosity for use in large-scale eddy-permitting ocean models. Monthly Weather Review, 128(8), 2935-2946. Abstract PDF
This paper discusses a numerical closure, motivated from the ideas of Smagorinsky, for use with a biharmonic operator. The result is a highly scale-selective, state-dependent friction operator for use in eddy-permitting geophysical fluid models. This friction should prove most useful for large-scale ocean models in which there are multiple regimes of geostrophic turbulence. Examples are provided from primitive equation geopotential and isopycnal-coordinate ocean models.
This paper discusses spurious diapycnal mixing associated with the transport of density in a z-coordinate ocean model. A general method, based on the work of Winters and collaborators, is employed for empirically diagnosing an effective diapycnal diffusivity corresponding to any numerical transport process. This method is then used to quantify the spurious mixing engendered by various numerical representations of advection. Both coarse and fine resolution examples are provided that illustrate the importance of adequately resolving the admitted scales of motion in order to maintain a small amount of mixing consistent with that measured within the ocean's pycnocline. Such resolution depends on details of the advection scheme, momentum and tracer dissipation, and grid resolution. Vertical transport processes, such as convective adjustment, act as yet another means to increase the spurious mixing introduced by dispersive errors from numerical advective fluxes.
Pacanowski, Ronald C., and Stephen M Griffies, 1999: The MOM3 Manual, GFDL Ocean Group Technical Report No. 4, Princeton, NJ: NOAA/Geophysical Fluid Dynamics Laboratory, 680 pp.
Schneider, Tapio, and Stephen M Griffies, 1999: A conceptual framework for predictability studies. Journal of Climate, 12(10), 3133-3155. Abstract PDF
A conceptual framework is presented for a unified treatment of issues arising in a variety of predictability studies. The predictive power (PP), a predictability measure based on information-theoretical principles, lies at the center of the framework. The PP is invarient under linear coordinate transformations and applies to multivariate predictions irrespective of assumptions about the probability distribution of prediction errors. For univariate Gaussian predictions, the PP reduces to conventional predictability measures that are based upon the ratio of the rms error of a model prediction over the rms error of the climatological mean prediction.
Since climatic variability on intraseasonal to interdecadal timescales follows an approximately Gaussian distribution, the emphasis of this paper is on multivariate Gaussian random variables. Predictable and unpredictable components of multivariate Gaussian systems can be distinguished by predictable component analysis, a procedure derived from discriminant analysis: seeking components with large PP leads to an eigenvalue problem, whose solution yields uncorrelated components that are ordered by PP from largest to smallest.
In a discussion of the application of the PP and the predictable component analysis in different types of predictability studies, studies are considered that use either ensemble integrations of numerical models or autoregressive models fitted to observed or simulated data.
An investigation of simulated multidecadal variability of the North Atlantic illustrates the proposed methodology. Reanalyzing an ensemble of integrations of the Geophysical Fluid Dynamics Laboratory coupled general circulation model confirms and refines earlier findings. With an autoregressive model fitted to a single integration of the same model, it is demonstrated that similar conclusions can be reached without resorting to computationally costly ensemble integrations.
This paper formulates tracer stirring arising from the Gent-McWilliams (GM) eddy-induced transport in terms of a skew-diffusive flux. A skew-diffusive tracer flux is directed normal to the tracer gradient, which is in contrast to a diffusive tracer flux directed down the tracer gradient. Analysis of the GM skew flux provides an understanding of the physical mechanisms prescribed by GM stirring, which is complementary to the more familiar advective flux perspective. Additionally, it unifies the tracer mixing operators arising from Redi isoneutral diffusion and GM stirring. This perspective allows for a computationally efficient and simple manner in which to implement the GM closure in z-coordinate models. With this approach, no more computation is necessary than when using isoneutral diffusion alone. Additionally, the numerical realization of the skew flux is significantly smoother than the advective flux. The reason is that to compute the skew flux, no gradient of the diffusivity or isoneutral slope is taken, whereas such a gradient is needed for computing the advective flux. The skew-flux formulation also exposes a striking cancellation of terms that results when the GM diffusion coefficient is identical to the Redi isoneutral diffusion coefficient. For this case, the horizontal components to the tracer flux are aligned down the horizontal tracer gradient, and the resulting computational cost of Redi diffusion plus GM skew diffusion is roughly half that needed for Redi diffusion alone.
Griffies, Stephen M., Anand Gnanadesikan, Ronald C Pacanowski, V D Larichev, J K Dukowicz, and R D Smith, 1998: Isoneutral diffusion in a z-coordinate ocean model. Journal of Physical Oceanography, 28(5), 805-830. Abstract PDF
This paper considers the requirements that must be satisfied in order to provide a stable and physically based isoneutral tracer diffusion scheme in a z-coordinate ocean model. Two properties are emphasized: 1) downgradient orientation of the diffusive fluxes along the neutral directions and 2) zero isoneutral diffusive flux of locally referenced potential density. It is shown that the Cox diffusion scheme does not respect either of these properties, which provides an explanation for the necessity to add a nontrivial background horizontal diffusion to that scheme. A new isoneutral diffusion scheme is proposed that aims to satisfy the stated properties and is found to require no horizontal background diffusion.
Bryan, Kirk, and Stephen M Griffies, 1997: Predictability of North Atlantic climate on decadal times scales estimated using a coupled ocean-atmosphere model. International WOCE Newsletter, 26, 5-9.
Atmospheric weather systems become unpredictable beyond a few weeks, but climate variations can be predictable over much longer periods because of the coupling of the ocean and atmosphere. With the use of a global coupled ocean-atmosphere model, it is shown that the North Atlantic may have climatic predictability on the order of a decade or longer. These results suggest that variations of the dominant multidecadal sea surface temperature patterns in the North Atlantic, which have been associated with changes in climate over Eurasia, can be predicted if an adequate and sustainable system for monitoring the Atlantic Ocean exists.
The North Atlantic is one of the few places on the globe where the atmosphere is linked to the deep ocean through air-sea interaction. While the internal variability of the atmosphere by itself is only predictable over a period of one to two weeks, climate variations are potentially predictable for much longer periods of months or even years because of coupling with the ocean. This work presents details from the first study to quantify the predictability for simulated multidecadal climate variability over the North Atlantic. The model used for this purpose is the GFDL coupled ocean-atmosphere climate model used extensively for studies of global warming and natural climate variability. This model contains fluctuations of the North Atlantic and high-latitude oceanic circulation with variability concentrated in the 40-60 year range. Oceanic predictability is quantified through analysis of the time-dependent behavior of large-scale empirical orthogonal function (EOF) patterns for the meridional stream function, dynamic topography, 170 m temperature, surface temperature and surface salinity. The results indicate that predictability in the North Atlantic depends on three main physical mechanisms. The first involves the oceanic deep convection in the subpolar region which acts to integrate atmospheric fluctuations, thus providing for a red noise oceanic response as elaborated by Hasselmann. The second involves the large-scale dynamics of the thermohaline circulation, which can cause the oceanic variations to have an oscillatory character on the multidecadal time scale. The third involves non-local effects on the North Atlantic arising from periodic anomalous fresh water transport advecting southward from the polar regions in the East Greenland Current. When the multidecadal oscillatory variations of the thermohaline circulation are active, the first and second EOF patterns for the North Atlantic dynamic topography have predictability time scales on the order of 10-20 y, whereas EOF-1 of SST has predictability time scales of 5-7 y. When the thermohaline variability has weak multidecadal power, the Hasselmann mechanism is dominant and the predictability is reduced by at least a factor of two. When the third mechanism is in an extreme phase, the North Atlantic dynamic topography patterns realize a 10-20 year predictability time scale. Additional analysis of SST in the Greenland Sea, in a region associated with the southward propagating fresh water anomalies, indicates the potential for decadal scale predictability for this high latitude region as well. The model calculations also allow insight into regional variations of predictability, which might be useful information for the design of a monitoring system for the North Atlantic. Predictability appears to break down most rapidly in regions of active convection in the high-latitude regions of the North Atlantic.
The comment by Rahmstorf suggests that a numerical problem in Tziperman et al. (1994, TTFB) leads to a noisy E - P field that invalidates TTFB's conclusions. The authors eliminate the noise, caused by the Fourier filtering used in the model, and show that TTFB's conclusions are still valid. Rahmstorf questions whether a critical value in the freshwater forcing separates TTFB's stable and unstable runs. By TTFB's original definition, the unstable runs in both TTFB and in Rahmstorf's comment have most definitely crossed a stability transition point upon switching to mixed boundary conditions. Rahmstorf finally suggests that the instability mechanism active in TTFB is a fast convective mechanism, not the slow advective mechanism proposed in TTFB. The authors show that the timescale of the instability is, in fact, consistent with the advective mechanism
Griffies, Stephen M., and E Tziperman, 1995: A linear thermohaline oscillator driven by stochastic atmospheric forcing. Journal of Climate, 8(10), 2440-2453. Abstract PDF
The interdecadal variability of a stochastically forced four-box model of the oceanic meridional thermohaline circulation (TMC) is described and compared to the THC variability in the coupled ocean-atmosphere GCM of Delworth, Manabe, and Stouffer. The box model is placed in linearly stable thermally dominant mean state under mixed boundary conditions. A linear stability analysis of this state reveals one damped oscillatory THC mode in addition to purely damped modes. The variability of the model under a moderate amount of stochastic forcing, meant to emulate the random variability of the atmosphere affecting the coupled model's interdecadal THC variability, is studied. A linear interpretation, in which the damped oscillatory mode is of primary importance, is sufficient for understanding the mechanism accounting for the stochastically forced variability. Direct comparison of the variability in the box model and coupled GCM reveals common qualitative aspects. Such a comparison supports, although does not verify, the hypothesis that the coupled model's THC variability can be interpreted as the result of atmospheric weather exciting a linear damped oscillatory THC mode.
Griffies, Stephen M., and Kirk Bryan, 1994: Predictability of North Atlantic climate variability on multidecadal time scales In The Atlantic Climate Change Program, Proceedings from the principal investigators meeting, NOAA, University Corporation for Atmospheric Research, 77-80. Abstract
A major goal of the ACCP program is to gain the understanding of North Atlantic climate variability required for making predictions. An essential first step in this direction is to assess the predictability of Atlantic climate variability from models. A methodology for doing this was first proposed by Lorenz (1965) for atmospheric models. Recently, predictability studies have been extended to coupled atmosphere-ocean models in connection with the El Niño/Southern Oscillation phenomenon (e.g., Cane and Zebiak, 1987; Goswami and Shukla, 1991). At present, no operational monitoring system exists to provide proper initial conditions for the ocean on a global basis or even for the North Atlantic. The goal of this study is to use the GFDL climate model to determine the value, in terms of practical prediction of multi-decadal climate variability, of an operational, deep-sea observing system. We present here preliminary results toward this goal.
