This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
pmip3:design:21k:final [2010/03/25 12:56] jypeter |
pmip3:design:21k:final [2010/03/26 14:14] jypeter |
||
---|---|---|---|
Line 6: | Line 6: | ||
===== Summary table ===== | ===== Summary table ===== | ||
- | The changes listed in this table are with respect to the [[pmip3:design:pi:final|PI]] configuration. | + | The changes listed in this table are //with respect to the [[pmip3:design:pi:final|PI]] configuration//. |
^ ^ PMIP3 ^ Minimum solution ^ | ^ ^ PMIP3 ^ Minimum solution ^ | ||
| Orbital parameters | [ **ecc** = 0.018994 ] - [ **obl** = 22.949° ] - [ **peri-180°** = 114.42° ] | | | | Orbital parameters | [ **ecc** = 0.018994 ] - [ **obl** = 22.949° ] - [ **peri-180°** = 114.42° ] | | | ||
| Date of vernal equinox | March 21 at Noon | | | | Date of vernal equinox | March 21 at Noon | | | ||
- | | Trace gases | [ **CO<sub>2</sub>** = 185 ppm ] - [ **CH<sub>4</sub>** = 350 ppb ] - [ **N<sub>2</sub>O** = 200 ppb ] - [ **CFC** = 0 ] - [ **O<sub>3</sub>** = same as in PI ] | | | + | | Trace gases | [ **CO<sub>2</sub>** = 185 ppm ] - [ **CH<sub>4</sub>** = 350 ppb ] - [ **N<sub>2</sub>O** = 200 ppb ]\\ [ **CFC** = 0 ] - [ **O<sub>3</sub>** = same as in PI ] | | |
| Aerosols | Same as in CMIP5 PI (see //[[final#dust_forcing|Dust forcing]]// note below) | | | | Aerosols | Same as in CMIP5 PI (see //[[final#dust_forcing|Dust forcing]]// note below) | | | ||
| Solar constant | Same as in CMIP5 PI | | | | Solar constant | Same as in CMIP5 PI | | | ||
| Vegetation | Same as in CMIP5 PI | | | | Vegetation | Same as in CMIP5 PI | | | ||
- | | Ice sheets | Ice sheet extent and related changes in topography provided ([[final#ice-sheet|see below]]) | | | + | | Ice sheets | Ice sheet extent and related changes in topography [[final#ice-sheet_extent|provided]] | | |
- | | Topography and coastlines | Land-sea mask and topography changes provided ([[final#land-sea_mask|see below]]) | Minimum changes ([[final#land-sea_mask|see below]]) | | + | | Land surface elevation and coastlines | [[final#land-sea_mask|Land-sea mask]] and [[final#land_surface_elevation|elevation changes]] provided | Minimum changes ([[final#land-sea_mask|see below]]) | |
- | | Ocean bathymetry | Up to groups depending on the flexibility of the model | | | + | | Ocean bathymetry | Up to groups, depending on the flexibility of their ocean model (see [[final#ocean_bathymetry|note below]]) | | |
| River outflow | Modified according to a river pathway map (provided, [[final#river_outflow|see below]]) | Same as in CMIP5 PI | | | River outflow | Modified according to a river pathway map (provided, [[final#river_outflow|see below]]) | Same as in CMIP5 PI | | ||
| Ice sheet mass balance | Add excess LGM freshwater to ocean ([[final#ice-sheet_mass_balance|see below]]) | Same as in CMIP5 PI | | | Ice sheet mass balance | Add excess LGM freshwater to ocean ([[final#ice-sheet_mass_balance|see below]]) | Same as in CMIP5 PI | | ||
- | | Mean ocean salinity | +1PSU everywhere\\ (to be added once at the beginning of the simulation) | | | + | | Mean ocean salinity | +1 PSU everywhere\\ (to be added __once__ at the beginning of the simulation) | | |
| Mean atmospheric surface pressure | Global average equal to PI value ([[final#initial_conditions|see below]]) | | | | Mean atmospheric surface pressure | Global average equal to PI value ([[final#initial_conditions|see below]]) | | | ||
Line 28: | Line 28: | ||
==== Insolation ==== | ==== Insolation ==== | ||
- | Note that insolation should follows PMIP requirements. Please check it carefuly using the following tables ([[http://pmip2.lsce.ipsl.fr/design/tables/inso_tables_21k.shtml|LGM insolation tables]]) | + | Note that insolation should follows PMIP requirements. Please check it carefuly using the following tables ([[http://pmip2.lsce.ipsl.fr/design/tables/inso_tables_21k.shtml|LGM insolation tables]]). |
- | Check carefully the date of the Vernal equinox, because it has implications to compare the paleo and PI seasonal cycles | + | Check carefully the date of the Vernal equinox, because it has implications when comparing the paleo and PI seasonal cycles. |
- | ==== Ice-sheet ==== | + | ==== Ice-sheets and related changes ==== |
- | The ice sheet provided for PMIP3/CMIP5 LGM experiments is a blended product obtained by averaging three different ice sheets: ICE6G provided by Dick Peltier, MOCA provided by Lev Tarasov and ANU provided by Kurt Lambeck. More information about this is available in the [[pmip3:design:pi:final:icesheet|LGM icesheet description]] page. | + | The ice sheet provided for PMIP3/CMIP5 LGM experiments is a blended product obtained by averaging three different ice sheets: ICE-6G v2.0 provided by Dick Peltier, MOCA provided by Lev Tarasov and ANU provided by Kurt Lambeck. More information about this is available in the [[pmip3:design:pi:final:icesheet|LGM ice-sheets description]] page. |
- | + | ||
- | The ice sheet should be implemented following PMIP2 protocol by considering the CHANGE in surface elevation (over ice sheets and land due to sea level lowering) in order that all models will be perturbed (forced) in the same way. This change in elevation would be added by each modeling group to the topography used in their control run. | + | |
| | ||
- | \\ \\ Data files are available on the [[pmip3:design:download#pmip3-cmip5_files_for_boundary_conditions|Files for PMIP3 experiments page]]. | + | Data files are available on the [[pmip3:design:download#pmip3-cmip5_files_for_boundary_conditions|Files for PMIP3 experiments]] page. |
- | ==== Land-sea mask ==== | + | Implementing the LGM **ice-sheets** implies changing also the **land-sea mask**, the **land surface elevation** and the **ocean bathymetry** as described below. |
- | The land-sea mask is provided by the ICE6G model, except for Antarctica where you should deduce it from the averaged ice-sheet. | + | === Land-sea mask === |
- | \\ Data files are available on the [[pmip3:design:download#pmip3-cmip5_files_for_boundary_conditions|Files for PMIP3 experiments page]]. | + | This mask should be prescribed using the data provided in [[target.nc|file.nc]] where the ''sftlf'' variable is 0 (%) over the ocean and 100 (%) over land. |
+ | If this proves difficult for a modelling group, the **__minimum__** implementation is to: | ||
+ | * fill in the Hudson Bay and the Barents sea, which were covered by the main northern hemisphere ice-sheets. | ||
+ | * close the Bering Strait (by changing the land-sea mask or imposing zero flux on water and sea-ice across the Bering Strait). | ||
+ | * make sure that the ocean model configuration keeps the following straits open: | ||
+ | * Denmark Strait (between Greenland and Ireland). | ||
+ | * Gibraltar Strait. | ||
+ | * //Indonesia Strait// (//check exact name with Olivier Marti//?). | ||
- | Note: the **minimum change** is to close the Bering Strait by changing the land-sea mask or imposing zero flux on water and sea-ice across the Bering Strait. | + | === Land surface elevation === |
+ | |||
+ | The change in surface elevation specified by the ''orog_diff'' variable provided in [[target.nc|file.nc]] should be added to the present-day surface elevation, so that all the participating models will be perturbed in the same way. | ||
+ | |||
+ | Note: this specified change assumes that the surface elevation over present-day oceans is zero (m). | ||
+ | |||
+ | === Ice-sheet extent === | ||
+ | |||
+ | The ice-sheet extent (aka //land-sea mask//) should be prescribed using the data provided in [[target.nc|file.nc]], where the ''sftgif'' variable is 0 (%) where there is no ice-sheet, or 100 (%).\\ This mask includes **ice shelves**, which should be included in models that allow this feature. | ||
+ | * The ice shelves can be defined from the ''sftgif'' and ''sftlf'' variables as areas where:\\ ''sftgif'' == 100 **and** ''sftlf'' == 0 | ||
+ | * Groups implementing the provided ice shelves should keep the Greenland-Iceland ridge open in the ocean configuration. | ||
+ | |||
+ | === Ocean bathymetry === | ||
+ | |||
+ | The mean sea-level change corresponding to the ice-sheet reconstruction to be used in the PMIP3 LGM experiment is approximatively NNN m. A more complete description of the change in bathymetry can be computed from the Peltier [[https://pmip3.lsce.ipsl.fr/under_construction_ask_jy|ICE-6G v2.0 reconstruction]]. | ||
+ | |||
+ | The ocean bathymetry will be adapted by each group, depending on the flexibility of their ocean model. | ||
+ | |||
+ | Modelling groups should pay special attention to: | ||
+ | * The bathymetry should be zero under the grounded ice-sheets. | ||
+ | * ... (see also the other ocean model configuration checks in this section). | ||
==== River outflow ==== | ==== River outflow ==== | ||
Line 55: | Line 80: | ||
You can use the river routing provided by Lev Tarasov to change the routing scheme in your code, so that the river pathways is consistent with the presence of the ice-sheet. | You can use the river routing provided by Lev Tarasov to change the routing scheme in your code, so that the river pathways is consistent with the presence of the ice-sheet. | ||
- | \\ Data files are available on the [[pmip3:design:download#pmip3-cmip5_files_for_boundary_conditions|Files for PMIP3 experiments page]]. | + | \\ Data files are available on the [[pmip3:design:download#pmip3-cmip5_files_for_boundary_conditions|Files for PMIP3 experiments]] page. |
==== Ice-sheet mass balance ==== | ==== Ice-sheet mass balance ==== | ||
Line 72: | Line 97: | ||
==== For Earth System Models with interactive carbon cycle ==== | ==== For Earth System Models with interactive carbon cycle ==== | ||
- | The simulations should be forced by the prescribed LGM CO2 concentrations. Please use the same protocol as in CMIP5 to store the diagnostic carbon fluxes and the variables needed for PCMIP (see | + | The simulations should be forced by the prescribed LGM CO<sub>2</sub> concentrations. Please use the same protocol as in CMIP5 to store the diagnostic carbon fluxes and the variables needed for PCMIP (see |
- | [[http://www.bridge.bris.ac.uk/projects/pcmip/experiments.html|pcmip Project]] and [[http://cmip-pcmdi.llnl.gov/cmip5/|CMIP5]]) | + | [[http://www.bridge.bris.ac.uk/projects/pcmip/experiments.html|PCMIP Project]] and [[http://cmip-pcmdi.llnl.gov/cmip5/|CMIP5]]) |
==== Dust forcing ==== | ==== Dust forcing ==== | ||
- | Some ESM have interative aerosols. In that case compute dust and associated forcing online, as in PI. If this is not the case then the recommendation is to keep dust and aerosols as in PI (I.E no change for 21 ka). | + | Some ESM have interactive aerosols. In that case compute dust and associated forcing online, as in PI. If this is not the case then the recommendation is to keep dust and aerosols as in PI (i.e. no change for 21 ka). |
===== Initial conditions ===== | ===== Initial conditions ===== | ||
- | The surface pressure field must be adjusted to the change in surface elevation over the continents. | ||
- | This can be done: | ||
- | * either by gradually changing the surface elevation in order to avoid generating gravity waves, | ||
- | * or by adjusting the initial pressure field to the LGM surface elevation. | ||
- | If you choose the second option, you must be careful to conserve atmospheric mass. | ||
- | The spin up procedure is up to group, following CMIP5 approach or from a previous cold state | + | * Mean ocean salinity: add +1 PSU everywhere once, at the beginning of the simulation. |
+ | |||
+ | * The surface pressure field must be adjusted to the change in surface elevation over the continents. This can be done: | ||
+ | * either by gradually changing the surface elevation in order to avoid generating gravity waves, | ||
+ | * or by adjusting the initial pressure field to the LGM surface elevation.\\ In this case, you must be careful to conserve atmospheric mass. | ||
- | Note several groups share the same ocean models for which an initial state can be provided from PMIP2 experiments. | + | * The spin up procedure is up to each group, following CMIP5 approach or from a previous cold state. |
- | If you need an initial state from a other group, you can contact Olivier.marti@lsce.ipsl.fr | + | * Note several groups share the same ocean models for which an initial state can be provided from PMIP2 experiments.\\ If you need an initial state from another group, you can get in touch with [[Olivier.marti@lsce.ipsl.fr]]. |
- | Groups with high resolution model for which it is too difficult to run long simulations should contact Olivier.marti@lsce.ipsl.fr to find the best alternative solution. | + | * Groups with high resolution model for which it is too difficult to run long simulations should contact [[Olivier.marti@lsce.ipsl.fr]] to find the best alternative solution. |
\\ \\ \\ | \\ \\ \\ |