Registration of asynchronously coupled transient simulation of MIS3 with interactive ice sheets and solid earth dynamics

assigned to @k204227

@k204227 @k204212

I downloaded the data from the data archive to /work/bk1192/from_Mistral/bk1192/WG1/WP1.1/MPIM_transient_deglaciation_asynchronous

In this folder you find the atmospheric 2D Data (BOT files) as well as the output of the PISM model (pism_-${year}.nc).

The time axis is structured as follows: Simulation year 4100 of the atmospheric files corresponds to year 49000 of the ice-sheet files. The atmospheric simulation is run for 10 years. Then the ice sheet model is run for 100 years based on the 10-year mean of the atmosspheric/oceanic/land model.

If you need any more information please feel to contact me.

Note, that after discussion in the last meeting we would like to also include the iceberg variables as defined from the cmip project: number_of_icebergs_per_unit_area: The number of icebergs per unit area. Unit: m-2 water_flux_into_sea_water_from_icebergs: The water flux into sea water is the freshwater entering as a result of precipitation, evaporation, river inflow, sea ice effects and water flux correction (if applied). The water flux into sea water from icebergs is due to the melting of the iceberg. In accordance with common usage in geophysical disciplines, "flux" implies per unit area, called "flux density" in physics. Unit: kg m-2 s-1 We can provide details once needed.

@m300467 @b380454 @m300792 Please provide the missing information for your requested new variables.

The problem with the topography and glacier mask is, that this variables already exist with decadal frequency in the dec MIPtable. So one needs to define a new MIPtable for this variables. The MIPtable name can be influenced by the submodel type (in which submodel does this play a role? will the LR topography and glacier mask be provided as well? are these HR variables provided on a different grid than any submodel grid (echam, vilma, mpiom, pism)?).

In case a new standard_name has to be registered for relative sea level, this can take several months and requires a lot further input. Existing standard_names can be browsed here.

variable	standard_name	CMOR label	frequency / MIPtable	units	cell_methods
relative sea level RSL	?	?	?	?	area: mean where sea time: mean
thermosteric change in mean sea level	thermosteric_change_in_mean_sea_level	zosto	?	m	area: mean where sea time: mean
HR topography	surface_altitude	orog	?	m	area: time: mean
HR glacier mask	land_ice_area_fraction	sftgif	?	%	area: time: mean
Water Flux into Sea Water from Icebergs	water_flux_into_sea_water_from_icebergs	ficeberg	?	kg m-2 s-1	area: mean where sea time: mean
number of icebergs	number_of_icebergs_per_unit_area	niceberg	(?)dec	m-2	area: time: mean

@m300467 In this asynchronously coupled simulations, mPISM output has centennial resolution and not decadal resolution? I did not think is run for 100 years means that the mPISM frequency changes from 10a to 100a, I rather expected 10 timesteps for this 100a. This would mean, that we need to create new MIP tables for the mPISM variables with 100a frequency.

@k204212 In these asynchronous simulations the mPISM output is only written out every 100 years. Some output exists also every 10 years, but it is o.k. for this simulation to use the 100 year frequency and likely for the other simulations of the past. But just a note for the future: There will likely be some future experiments where the output frequency of mPISM will change to one year. Is this something that we already need to consider now (e.g. for registering mPISM variables with an annual frequency)?

Regarding the missing variable information of 'number of icebergs': frequency = same as mpiom; units = m-2; cell_methods = time:mean. CMOR label could be niceberg or anything else in this direction or the name that is currently used.

will the LR topography and glacier mask be provided as well? are these HR variables provided on a different grid than any submodel grid (echam, vilma, mpiom, pism)?). --> The HR variables are provided on a 1/6 degree grid. It is used in the HR model and for the fluxes. The CR variables will also have to be stored, such as topography and glacier mask on the echam grids. @m300792 - also on the mPISM grid?!

@b380454 is the sea level change comparable to the following variable in the mip table? -->

global_average_sea_level_change

Global average sea level change is due to change in volume of the water in the ocean, caused by mass and/or density change, or to change in the volume of the ocean basins, caused by tectonics etc. It is sometimes called "eustatic", which is a term that also has other definitions. It differs from the change in the global average sea surface height relative to the centre of the Earth by the global average vertical movement of the ocean floor. Zero sea level change is an arbitrary level. Because global average sea level change quantifies the change in volume of the world ocean, it is not calculated necessarily by considering local changes in mean sea level.

@m300467 I had a short phone call with Uwe:

the thermosteric change in mean sea level (zosto / zostoga) should not be published for simulations with active VILMA or prescribed time dependent topography changes
the HR topography and glacier masks: either they replace the echam/jsbach variants, or one has to come up with a new MIP table name (eg. SLmon / SLdec for "Solid Land monthly / decadal")
a point that I did not bring up yet: we should also submit a cell_area variable for VILMA (assuming the VILMA grid is the 1/6 degree grid the glacier mask and topography are provided on), just as we did with areacello and areacella for MPIOM and echam/jsbach. If this is not the VILMA grid we should add the 1/6 degree grid submodel to the source description.

mentioned in issue #6 (closed)

mentioned in issue #9

Registration of asynchronously coupled transient simulation of MIS3 with interactive ice sheets and solid earth dynamics

Designs

Child items ...

Activity