The class content has moved here:
https://github.com/obidam/ds2-2020/tree/ds2-2021
- Fork this github repository to your own account
- Add a folder under ds2-2021/projects, name it with your last name
- Work out your project codes/documentation and stage/push it to your folder
- Once finalized, create a pull request to the main branch
Description: Because of the human driven intensification of the greenhouse effect, the ocean is warming. You can compute ocean heat content (OHC) and its trend with a regression (linear or not) for the entire ocean time series and extrapolate to the future, e.g. what is the expected ocean warming for the horizon 2100.
Instead of working globally, you can study the ocean warming locally. In that case, you can plot the local slopes of the different OHC time series and deduce where the ocean warming is moderate and where it is strong.
Data:
By default, you can use the EN4 dataset that is an interpolation of all available ocean observations (of temperature and salinity) onto a regular space/time grid.
This dataset can be accessed this way:
from intake import open_catalog
catalog_url = 'https://raw.githubusercontent.com/obidam/ds2-2021/main/ds2_data_catalog.yml'
cat = open_catalog(catalog_url)
ds = cat["en4"].to_dask()
or:
fs = gcsfs.GCSFileSystem(project="ds2class-2021")
gcsmap = fs.get_mapper("ds2data/EN.4.2.1.f.analysis.g10.zarr")
ds = xr.open_zarr(gcsmap)
But, if you wish, you can try to use climate model simulations data that accessible here: CMIP6 data.
Description: Sea level increases because of changes in currents (dynamic effect) and because of ocean density changes (steric effect). Compute ocean density changes contribution to Sea level rises (thermosteric effect) and demonstrate that it is the driver of regional sea level change trends.
Data: You can use the EN4 dataset that is an interpolation of all available ocean observations (of temperature and salinity) onto a regular space/time grid.
This dataset can be accessed this way:
from intake import open_catalog
catalog_url = 'https://raw.githubusercontent.com/obidam/ds2-2021/main/ds2_data_catalog.yml'
cat = open_catalog(catalog_url)
ds = cat["en4"].to_dask()
Sea Level data accessible on the catalog as well here:
ds = cat["sea_surface_height"].to_dask()
Description: The Coupled Model Intercomparison Project is a framework within which a number of research centres produce predictions of the future evolution of the climate. Many climate centres use an ensemble approach, in which they produce several simulations for each scenario, using the same model configuration. By comparing these multiple simulations (“the ensemble"), we can estimate the probability that the climate will evolve in a certain way.
Data: Google host a part of these data via their Public Datasets programme (https://cloud.google.com/blog/products/data-analytics/new-climate-model-data-now-google-public-datasets).
This dataset can be accessed this way:
df = pd.read_csv('https://storage.googleapis.com/cmip6/cmip6-zarr-consolidated-stores.csv')
df_ssh = df.query("activity_id=='ScenarioMIP' & table_id == 'Omon' & variable_id == 'zos' & institution_id == 'IPSL'")
gcs = gcsfs.GCSFileSystem(token='anon')
zstore = df_ssh.zstore.values[-1]
mapper = gcs.get_mapper(zstore)
ds = xr.(mapper, consolidated=True)