Temperature, Salinity, and Stratification¶

Now that we have learned the basic thermodynamics of seawater and understood how temperature and salinity combine to determine the stratifcation of water columns, we are ready to examine the global patterns of temperature, salinity and stratification. These patterns are fascinating and reveal much about the underlying ocean circulation. We will also learn how to identify distinct water masses both in geographic space and in temperature / salinity coordinates.

import xarray as xr
from matplotlib import pyplot as plt
import hvplot.xarray
import gsw
import numpy as np
import holoviews as hv
from IPython.display import Image
hv.extension('bokeh')
plt.rcParams['figure.figsize'] = (12,7)

3D Temperature¶

def map_section_plot(data, label, data_range, levels, cmap='viridis'):
    ds = hv.Dataset(data)
    img = ds.to(hv.Image, ['lon', 'lat'], label=label)
    img = img.redim.range(**{data.name: data_range})

    posx = hv.streams.PointerX()
    vline = hv.DynamicMap(lambda x: hv.VLine(x or 180.5), streams=[posx])

    lats = np.arange(-90, 91)
    def cross_section(x):
        mesh = hv.QuadMesh((lats, woa.depth, data.sel(lon=x if x else 180, method='nearest').data[:-1]),
                           kdims=['lat', 'depth'], label=(label + ' Section'))
        contours = hv.operation.contours(mesh, levels=levels,
                                         overlaid=False, filled=True)
        return contours

    
    crosssection = hv.DynamicMap(cross_section, streams=[posx])
    layout = (img * vline) + crosssection 
    
    # set view options
    dict_spec = {'Image': {'plot': dict(colorbar=True, toolbar='above', width=600),
                           'style': dict(cmap=cmap)},
                 'NdOverlay': {'plot': dict(bgcolor='grey', invert_yaxis=True, width=600, height=200)},
                 'Layout': {'plot': dict(show_title=False)},
                 'Polygons': {'plot': dict(colorbar=True, invert_yaxis=True, width=600),
                              'style': dict(cmap=cmap, line_width=0.1)}
                 }
    layout.cols(1)
    layout = layout.opts(dict_spec)   
    return layout

map_section_plot(woa.t_an, 'Temperature', (-2,30), np.arange(-2,31), cmap='RdBu_r')

Observations¶

The warmest water is near the surface in the tropical regions, forming a bowl-shaped lens in the upper ~800 meters. This is called the thermocline.
There are clear zonal asymmetries in surface temperature. The western boundaries of the basins tend to be warmer, with colder water near the eastern boundaries.
There is something special going on at the equator, especially in the Pacific. It’s relatively cold, and the thermocline is shallow.
Below 1000 m, the water is mostly very cold (< 5 C).
In the polar regions, the coldest water is actually found near the surface! There is a mid-depth temperature maximum.

3D Salinity¶

map_section_plot(woa.s_an, 'Salinity', (32,38), np.arange(32,38.1, 0.25))

Observations¶

Salinity has a more complex structure than temperature.
The Atlantic is overall saltier than the Pacific.
The saltiest parts of the oceans are the extra-tropical-latitudes.
The equator is relatively fresh near the surface.
Salinity sections across the Atlantic reveal the interleaving of several different water masses including
- Mediterranean Water
- North Atlantic Deep Water
- Antarctic Intermediate Water
- Antarctic Bottom Water

3D Potential Density¶

map_section_plot(woa.sigma_0, 'Potential Density (sigma0)', (22.,28.), np.arange(22,28.1, 0.25))

map_section_plot(woa.sigma_4, 'Potential Density (sigma4)', (40.,46.5), np.arange(40,46.5, 0.25))

Intro to Physical Oceanography

Temperature, Salinity, and Stratification¶

Data Source: World Ocean Atlas¶

3D Temperature¶

Observations¶

3D Salinity¶

Observations¶

3D Potential Density¶