Water masses of the future

Water masses of the future

Scientific Motivation

Water masses have historically been defined using certain potential temperature and salinity bounds for each given one. In view of projected changes in heat content and freshwater fluxes, will the bounds defined in the mid-20th century still hold?
Other approaches, based on geometric consideration (see Juza et al. https://www.sciencedirect.com/science/article/pii/S0924796318302173), have been proposed. The goal here is to investigate if significant water mass properties changes are found between late 20th versus late 21th centuries) using T/S diagrams, as well as changes in water masses total volume. And if a significant change in temperature is detected, estimate the heat content change in some water masses of interest. There is room for discussion about what criterion can/should be used for water mass description.

Proposed Hacking

Participants can focus on one (or more) water mass of particular interest to him/her.

  1. Building T/S diagrams for multiple CMIP ocean models and compare 20th and 21th centuries, detect potential changes.

  2. Compute water mass volume evolution (traditional criterion and maybe new one(s))

  3. Estimate changes in heat content.

Anticipated Data Needs

CMIP (5/6) potential temperature and salinity

Anticipated Software Tools

  1. pangeo software stack (xarray/dask)
  2. computation of water mass volume and heat content is a 3 for loop, 4 if statements that can be optimized for performance with f2py or numba.

Desired Collaborators

Beginner/Intermediate python users are welcome, basic knowledge of observational oceanography (can read T/S diagram).

3 Likes

I’m interested in this project. It isn’t very related to my research (though I am an oceanographer) but I’m happy to collaborate with you on this

@raphaeldussin - this is a great idea. I’d be interested in following this project.

Raphael–this sort of application was exactly why we built xhistogram. Please consider trying it out.

Given the similarity in goals and tools, I’m wondering if it would make sense to combine this project with Plotting ocean variables in density coordinates (using xhistogram) (proposed by @cspencerjones). Over there, they are talking about plotting different variables in density space. Here it’s using T/S space. Do we really want two separate projects? Or can we collaborate?

This seems to be a sensible suggestion to me.

On the science side, I have been thinking about this problem quite a bit. How do you properly compare to watermass distributions? What is the proper distance metric.

I believe this can be answered through optimal transport theory:

The idea is to figure out how much work it takes to transform one distribution into another. This problem can be solved efficiently, and there is a cool python package for doing it

Particular to the ocean T/S problem is the question of the cost of moving water in T/S space. One could imagine an energy-based metric, which penalizes diapycnal transformations in favor of isopycnal ones.

@rabernat xhistogram seems like the kind of tool that will be very useful for the projects. I think it’s also a good idea to merge with Spencer’s project.

I guess the question here is how may people are gonna be working on this hack
and can we try to distribute the efforts in a way where everybody got a bit to chew on while keeping the idea of a working group.

I’m not super familiar with the format so we’ll need the organizers to help out with this. But this looks like some great fun coming up, looking forward to it!

I am interested in this project. We have been looking at LSW and finding it rather difficult to pin point the density ranges that define this water mass. I suspect this is the case for other water masses as well. Certainly a group exercise. Thanks for posting this.

I have put together a couple of quick examples that we can use as a starting point:

1 Like

@aromanou - I was reading a paper by JB Sallee on water masses of the Southern Ocean in CMIP5 historical runs (Sallee et al 2013). It said that they had to manually identify the density ranges corresponding to the water masses for each individual model, as they were not necessarily the same due to density biases between models and obs.