Modeling Water-Mass Distributions in the Modern and LGM Ocean: Circulation Change and Isopycnal and Diapycnal Mixing

Paleoproxy observations suggest that deep-ocean water-mass distributions were different at the Last Glacial Maximum than they are today. However, even modern deep-ocean water-mass distributions are not completely explained by observations of the modern ocean circulation. This paper investigates two processes that influence deep-ocean water-mass distributions: 1) interior downwelling caused by vertical mixing that increases in the downward direction and 2) isopycnal mixing. Passive tracers are used to assess how changes in the circulation and in the isopycnal-mixing coefficient impact deep-ocean water-mass distributions in an idealized two-basin model. We compare two circulations, one in which the upper cell of the overturning reaches to 4000-m depth and one in which it shoals to 2500-m depth. Previous work suggests that in the latter case the upper cell and the abyssal cell of the overturning are separate structures. Nonetheless, high concentrations of North Atlantic Water (NAW) are found in our model’s abyssal cell: these tracers are advected into the abyssal cell by interior downwelling caused by our vertical mixing profile, which increases in the downward direction. Further experiments suggest that the NAW concentration in the deep South Atlantic Ocean and in the deep Pacific Ocean is influenced by the isopycnal-mixing coefficient in the top 2000 m of the Southern Ocean. Both the strength and the vertical profile of isopycnal mixing are important for setting deep-ocean tracer concentrations. A 1D advection–diffusion model elucidates how NAW concentration depends on advective and diffusive processes.

Internal salt content: a useful framework for understanding the oceanic branch of the water cycle

The global water cycle is dominated by an atmospheric branch which transfers fresh water away from subtropical regions and an oceanic branch which returns that fresh water from subpolar and tropical regions. Salt content is commonly used to understand the oceanic branch because surface freshwater fluxes leave an imprint on ocean salinity. However, freshwater fluxes do not actually change the amount of salt in the ocean and – in the mean – no salt is transported meridionally by ocean circulation. To study the processes which determine ocean salinity we introduce a new variable: “internal salt” and its counterpart “internal fresh water”. Precise budgets for internal salt in salinity coordinates relate meridional and diahaline transport to surface freshwater forcing, ocean circulation and mixing, and reveal the pathway of fresh water in the ocean. We apply this framework to a 1° global ocean model. We find that in order for fresh water to be exported from the ocean’s tropical and subpolar regions to the subtropics, salt must be mixed across the salinity surfaces that bound those regions. In the tropics, this mixing is achieved by parameterized vertical mixing, along-isopycnal mixing, and numerical mixing associated with truncation errors in the model’s advection scheme, while along-isopycnal mixing dominates at high latitudes. We analyze the internal freshwater budgets of the Indo-Pacific and Atlantic Ocean basins and identify the transport pathways between them which redistribute fresh water added through precipitation, balancing asymmetries in freshwater forcing between the basins.

Mixing Inferred from an Ocean Climatology and Surface Fluxes

This study provides observation-based estimates, determined by inverse methods, of horizontal and isopycnal eddy diffusion coefficients KH and KI, respectively, the small-scale mixing coefficient D, and the diathermohaline streamfunction Ψ. The inverse solution of Ψ represents the ocean circulation in Absolute Salinity SA and Conservative Temperature Θ coordinates. The authors suggest that the observation-based estimate of Ψ will be useful for comparison with equivalent diagnostics from numerical climate models. The estimates of KH and KI represent horizontal eddy mixing in the mixed layer and isopycnal eddy mixing in the ocean interior, respectively. This study finds that the solution for D and KH are comparable to existing estimates. The solution for KI is one of the first observation-based global and full-depth constrained estimates of isopycnal mixing and indicates that KI is an order of magnitude smaller than KH. This suggests that there is a large vertical variation in the eddy mixing coefficient, which is generally not included in ocean models. With ocean models being very sensitive to the choice of isopycnal mixing, this result suggests that further investigation of the spatial structure of isopycnal eddy mixing from observations is required.

Exploring the isopycnal mixing and helium–heat paradoxes in a suite of Earth system models

This paper uses a suite of Earth system models which simulate the distribution of He isotopes and radiocarbon to examine two paradoxes in Earth science, each of which results from an inconsistency between theoretically motivated global energy balances and direct observations. The helium–heat paradox refers to the fact that helium emissions to the deep ocean are far lower than would be expected given the rate of geothermal heating, since both are thought to be the result of radioactive decay in Earth's interior. The isopycnal mixing paradox comes from the fact that many theoretical parameterizations of the isopycnal mixing coefficient ARedi that link it to baroclinic instability project it to be small (of order a few hundred m2 s−1) in the ocean interior away from boundary currents. However, direct observations using tracers and floats (largely in the upper ocean) suggest that values of this coefficient are an order of magnitude higher. Helium isotopes equilibrate rapidly with the atmosphere and thus exhibit large gradients along isopycnals while radiocarbon equilibrates slowly and thus exhibits smaller gradients along isopycnals. Thus it might be thought that resolving the isopycnal mixing paradox in favor of the higher observational estimates of ARedi might also solve the helium paradox, by increasing the transport of mantle helium to the surface more than it would radiocarbon. In this paper we show that this is not the case. In a suite of models with different spatially constant and spatially varying values of ARedi the distribution of radiocarbon and helium isotopes is sensitive to the value of ARedi. However, away from strong helium sources in the southeastern Pacific, the relationship between the two is not sensitive, indicating that large-scale advection is the limiting process for removing helium and radiocarbon from the deep ocean. The helium isotopes, in turn, suggest a higher value of ARedi below the thermocline than is seen in theoretical parameterizations based on baroclinic growth rates. We argue that a key part of resolving the isopycnal mixing paradox is to abandon the idea that ARedi has a direct relationship to local baroclinic instability and to the so-called "thickness" mixing coefficient AGM.