Bibliography - Bonita L Samuels

1. Griffies, Stephen, Robert W Hallberg, A Pirani, Bonita L Samuels, and Michael Winton, et al., January 2009: Coordinated ocean-ice reference experiments (COREs). Ocean Modelling, 26(1-2), doi:10.1016/j.ocemod.2008.08.007.
   [ Abstract ]

   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.

2. Anderson, Whit G., Anand Gnanadesikan, Robert W Hallberg, John Dunne, and Bonita L Samuels, June 2007: Impact of ocean color on the maintenance of the Pacific Cold Tongue. Geophysical Research Letters, 34, L11609, doi:10.1029/2007GL030100.
   [ Abstract ]

   The impact of the penetration length scale of shortwave radiation into the surface ocean is investigated with a fully coupled ocean, atmosphere, land and ice model. Oceanic shortwave radiation penetration is assumed to depend on the chlorophyll concentration. As chlorophyll concentrations increase the distribution of shortwave heating becomes shallower. This change in heat distribution impacts mixed-layer depth. This study shows that removing all chlorophyll from the ocean results in a system that tends strongly towards an El Niño state—suggesting that chlorophyll is implicated in maintenance of the Pacific cold tongue. The regions most responsible for this response are located off-equator and correspond to the oligotrophic gyres. Results from a suite of surface chlorophyll perturbation experiments suggest a potential positive feedback between chlorophyll concentration and a non-local coupled response in the fully coupled ocean-atmosphere system.

3. Gnanadesikan, Anand, Stephen Griffies, and Bonita L Samuels, 2007: Effects in a climate model of slope tapering in neutral physics schemes. Ocean Modelling, 16(1-2), doi:10.1016/j.ocemod.2006.06.004.
   [ Abstract ]

   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 S_max which in combination with the thickness diffusivity A_gm defines a maximum eddy-induced advective transport A_gm*S_max. In this paper, we examine the impact of changing S_max 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.

4. Gnanadesikan, Anand, Keith W Dixon, Stephen Griffies, Ventakramani Balaji, M Barreiro, J A Beesley, W F Cooke, Thomas L Delworth, R Gerdes, Matthew J Harrison, Isaac Held, William J Hurlin, H C Lee, Z Liang, G Nong, Ronald C Pacanowski, Anthony Rosati, J L Russell, Bonita L Samuels, Qian Song, Michael J Spelman, Ronald J Stouffer, C Sweeney, G A Vecchi, Michael Winton, Andrew T Wittenberg, Fanrong Zeng, Rong Zhang, and John Dunne, 2006: GFDL's CM2 Global Coupled Climate Models. Part II: The baseline ocean simulation. Journal of Climate, 19(5), doi:10.1175/JCLI3630.1.
   [ Abstract ]

   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.

5. Griffies, Stephen, Anand Gnanadesikan, Keith W Dixon, John Dunne, R Gerdes, Matthew J Harrison, Anthony Rosati, J L Russell, Bonita L Samuels, Michael J Spelman, Michael Winton, and Rong Zhang, 2005: Formulation of an ocean model for global climate simulations. Ocean Science, 1, 45-79.
   [ Abstract PDF ]

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

6. Sweeney, C, Anand Gnanadesikan, Stephen Griffies, Matthew J Harrison, Anthony Rosati, and Bonita L Samuels, 2005: Impacts of shortwave penetration depth on large-scale ocean circulation and heat transport. Journal of Physical Oceanography, 35(6), 1103-1119.
   [ Abstract PDF ]

   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.

7. Gnanadesikan, Anand, Richard D Slater, and Bonita L Samuels, 2003: Sensitivity of water mass transformation and heat transport to subgridscale mixing in coarse-resolution ocean models. Geophysical Research Abstracts, 30(18), 1967, doi:10.1029/2003GL018036.
   [ Abstract PDF ]

   This paper considers the impact of the parameterization of subgridscale mixing on ocean heat transport in coarse-resolution ocean models of the type used in coupled climate models. Increasing the vertical diffusion increases poleward heat transport in both hemispheres. Increasing lateral diffusion associated with transient eddies increases poleward heat transport in the southern hemisphere while decreasing it in the northern hemisphere. The results are interpreted in the context of a simple analytical model.

8. Toggweiler, J R., and Bonita L Samuels, 1998: On the ocean's large-scale circulation near the limit of no vertical mixing. Journal of Physical Oceanography, 28(9), 1832-1852.
   [ Abstract PDF ]

   By convention, the ocean's large-scale circulation is assumed to be a thermohaline overturning driven by the addition and extraction of buoyancy at the surface and vertical mixing in the interior. Previous work suggests that the overturning should die out as vertical mixing rates are reduced to zero. In this paper, a formal energy analysis is applied to a series of ocean general circulation models to evaluate changes in the large-scale circulation over a range of vertical mixing rates. Two different model configurations are used. One has an open zonal channel and an Antarctic Circumpolar Current (ACC). The other configuration does not. The authors find that a vigorous large-scale circulation persists at the limit of no mixing in the model with a wind-driven ACC. A wind-powered overturning circulation linked to the ACC can exist without vertical mixing and without much energy input from surface buoyancy forces.

9. Toggweiler, J R., E Tziperman, Y Feliks, Kirk Bryan, Stephen Griffies, and Bonita L Samuels, 1996: Reply. Journal of Physical Oceanography, 26(6), 1106-1110.
   [ Abstract PDF ]

   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

10. Toggweiler, J R., and Bonita L Samuels, 1995: Effect of Drake Passage on the global thermohaline circulation. Deep-Sea Research, Part I, 42(4), 477-500.
    [ Abstract PDF ]

    The Ekman divergence around Antarctica raises a large amount of deep water to the ocean's surface. The regional Ekman transport moves the upwelled deep water northward out of the circumpolar zone. The divergence and northward surface drift combine, in effect, to remove deep water from the interior of the ocean. This wind-driven removal process is facilitated by a unique dynamic constraint operating in the latitude band containing Drake Passage. Through a simple model sensitivity experiment we show that the upwelling and removal of deep water in the circumpolar belt may be quantitatively related to the formation of new deep water in the northern North Atlantic. These results show that stronger winds in the south can induce more deep water formation in the north and more deep outflow through the South Atlantic. The fact that winds in the southern hemisphere might influence the formation of deep water in the North Atlantic brings into question long-standing notions about the forces that drive the ocean's thermohaline circulation.

11. Toggweiler, J R., and Bonita L Samuels, 1995: Effect of sea ice on the salinity of Antarctic bottom waters. Journal of Physical Oceanography, 25(9), 1980-1997.
    [ Abstract PDF ]

    Brine rejection during the formation of Antarctic sea ice is known to enhance the salinity of dense shelf waters in the Weddell and Ross Seas. As these shelf waters flow off the shelves and descend to the bottom, they entrain ambient deep water to create new bottom water. It is not uncommon for ocean modelers to modify salinity boundary conditions around Antarctica in an attempt to include a "sea ice effect" in their models. However, the degree to which Antarctic salinities are enhanced is usually not quantified or defended. In this paper, studies of shelf hydrography and delta¹⁸O are reviewed to assess the level of salinity enhancement appropriate for ocean general circulation models. The relevant quantities are 1) the salinity difference between the water masses modified on the shelves and the final offshelf flow and 2) the flux of salt (or freshwater) that gives rise to this salinity difference. Onshelf/offshelf salinity changes in the Weddell and Ross Seas appear to be fairly small, 0.15-0.20 salinity units. The quantity of brine needed to produce this salinification is equivalent to the salt drained from <0.50 m of new sea ice every year. Salt fluxes and salinity distributions from three GCM simulations are then compared. The first model has its surface salinities simply restored to the Levitus observations. Levitus restoring produces a slight freshening in the area of the Weddell and Ross Sea shelves. The global-mean bottom-water salinity in this model is 34.57 psu, which 0.16 units less than observed. The second model includes a very modest salinity enhancement in the area of the Weddell and Ross Sea shelves. This produces a salt flux equivalent to the formation of ~ 0.50 m yr^-1 of new sea ice. Even though this amount of salt input is close to the amount observed, global-average deep salinities in the second model are only 0.02 units greater than the deep salinities in the first model. The third model includes a large salinity enrichment, which is applied throughout the Weddell and Ross embayments without regard to water depth. Its deep salinities are 0.18 units higher than the deep salinities in the first model, but the amount of salt pumped into the model greatly exceeds the salt flux in nature. The authors conclude that salt from sea ice is probably not a major influence on the salinity of Antarctic bottom waters. Predicted salinities in ocean GCMs are too fresh because of circulation deficiencies, not because of inadequate boundary conditions. Models that employ large salinity modifications near Antarctica run the risk of grossly distorting the processes of deep-water formation.

12. Toggweiler, J R., and Bonita L Samuels, 1993: New radiocarbon constraints on the upwelling of abyssal water to the ocean's surface In The Global Carbon Cycle, NATO ASI Series - Vol. 115. Heidelberg, Germany, Springer-Verlag, 333-366.
    [ Abstract ]

    The output from seven different ocean model simulations is compared on the basis of the Δ ¹⁴C difference between North Pacific deep water and Antarctic surface water. This set of models produces a range of North Pacific-Antarctic Δ ¹⁴C differences between -173% and -108%, all but the smallest of which are substantially larger than the actual pre-bomb difference, -80 to -110%. Predicted values are highly correlated with the quantity of mid-depth water which flows out of the Pacific to the south. A circulation in which most of the Antarctic bottom water flows back out of the basin at mid-depth produces the smallest North Pacific-Antarctic Δ ¹⁴C differences, whereas a circulation in which all the inflow of bottom water upwells through the thermocline produces tha largest and least realistic differences. According to the models, upwelled abyssal water becomes entrained into the wind-driven convergence of thermocline water toward the equator. When it reaches the surface it spreads to the north and south, producing a Δ ¹⁴C minimum along the equator. A detailed analysis of both pre-bomb and post-bomb Δ ¹⁴C data indicates that the oldest water in the tropical Pacific is actually found south of the equator and is associated with the upwelling off Peru, not the upwelling along the equator. Toggweiler et al. (1991) trace the low-Δ ¹⁴C signal in the Peru upwelling to deep water raised to the surface around Antarctica which is pushed northward into the thermocline by circumpolar winds. According to the models, even a small amount of abyssal water upwelling through the thermocline (~3x10⁶ m3 s^-1) leaves a characteristic signal in the surface Δ ¹⁴C distribution which is not observed. One is left with the general conclusion that there is very little upwelling associated with a top-to-bottom thermohaline circulation in the world ocean. Virtually all the upwelling of abyssal water to the ocean's surface occurs around Antarctica where it is mainly wind-forced. The implications of this conclusion for the carbon cycle are discussed.

13. Toggweiler, J R., and Bonita L Samuels, 1993: Is the magnitude of the deep outflow from the Atlantic Ocean actually governed by Southern Hemisphere winds? In The Global Carbon Cycle, NATO ASI Series - Vol. 115. Heidelberg, Germany, Springer-Verlag, 303-331.
    [ Abstract ]

    The large-scale overturning in the Atlantic Ocean and its export of salty water to the other basins of the ocean is usually thought of as a thermohaline process driven by the formation of dense bottom water in the isolated basins of the North Atlantic. In this paper the output from several different runs of a global ocean GCM is used to show that the inflow of upper kilometer water in the South Atlantic and the outflow of deep water varies in direct proportion to the westerly wind stress in the circumpolar region of the southern hemisphere. According to the results presented here, the production of dense bottom water in the North Atlantic makes it possible for an Atlantic overturning to exist, but southern hemisphere winds appear to determine the magnitude of the inflow and outflow. The connection between southern hemisphere winds and the Atlantic overturning is due to a unique dynamic constraint which operates in the latitude band of Drake Passage. This constraint suggests the possibility of a very simple relationship between the magnitude of the northward wind drift at the latitude of the tip of South America and the magnitude of the inflow and outflow from the Atlantic basin.

Direct link to page: http://www.gfdl.noaa.gov/bibliography/resultstest.php?author=1088