matt-long / cmip6hack-ocean-o2 Goto Github PK

Examine upper ocean O2 change as a function of temperature.

Here are results from CMIP5, including observations of Ito et al. (2017).

Started in
notebooks/o2_v_temp.ipynb

Oxygen change results from the combination of direct temperature effects (O2sat) and indirect effects of circulation and respiration on AOU. CMIP5 model should consistent behavior in the extratropics, but divergent behavior in the tropics, leading a complicated picture at the global scale.

Let's make this plot for CMIP6.

Hi everyone,

I would like to focus on the equatorial circulation, and how it changes between models and under forcing.
I think this could be helpful for the context of several bgc projects and also the tropical SST gradient project.

Ill start of with a basic evaluation. Please join me if you are interested, or tell me what aspect of this could be important for other projects.

Comparisons of simulated O2 fields to WOA:

• maps of 200-600 m O2 with RMSE, bias, pattern corr
• volumes at 80, 50, 5 mmol/m^3
• etc.

What models have oxygen data? Well, some have it annually and others monthly. In notebooks/o2_v_temp.ipynb I wrote the following code to select the appropriate union of models.

Should this be moved to util.py to be shared across notebooks?

experiments = ['historical', 'ssp585']

def get_models(table_id, variable_id):
    # all models
    models = set(uni_dict['source_id']['values'])

    for experiment_id in experiments:
        query = dict(experiment_id=experiment_id, variable_id=variable_id, 
                     table_id=table_id, grid_label='gn')  
        cat = col.search(**query)
        models = models.intersection({model for model in cat.df.source_id.unique().tolist()})

    # ensure the CESM2 models are not included (oxygen was erroneously submitted to the archive)
    return models - {'CESM2-WACCM', 'CESM2'}
    
models = {}    

# look for models with O2
# prioritize getting data from annual table, if present
models['Oyr.o2'] = get_models('Oyr', 'o2')
models['Omon.o2'] = get_models('Omon', 'o2') - models['Oyr.o2']
have_oxygen = models['Oyr.o2'].union(models['Omon.o2'])

# find models with temperature that also have O2
models['Omon.theta'] = get_models('Omon', 'thetao').intersection(have_oxygen)

models_all = list(models['Omon.theta'])
models

How does the volume of hypoxic waters change? How does the volume of anoxic waters change?

Hypoxic threshold = 60 mmol/m^3
Anoxic threshold = 5 mmol/m^3

Copy content from discourse
https://discourse.pangeo.io/t/ocean-oxygen-in-a-warming-world/68