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:icesheet:index [2009/09/08 16:41] jypeter Replaced the ICE-6G figures with cleaner ones |
pmip3:design:21k:icesheet:index [2009/10/09 08:25] jypeter |
||
---|---|---|---|
Line 2: | Line 2: | ||
====== Ice Sheet for PMIP3/CMIP5 simulations ====== | ====== Ice Sheet for PMIP3/CMIP5 simulations ====== | ||
- | This page synthetises the Ongoing work for the definition of a new ice-sheet reconstruction to be used in PMIP3/CMIP5. | + | |
+ | |||
+ | |||
+ | <note warning>This page synthesises the **ongoing work** for the definition of a new ice-sheet reconstruction to be used in LGM simulations for PMIP3/CMIP5.\\ \\ Use the [[index#page_bottom_tag|discussion]] panel at the end of the page for comments </note> | ||
===== Introduction ===== | ===== Introduction ===== | ||
- | Several PMIP participants are in favour of using a new ice-sheet reconstruction that better matches the geomorphologic and glaciological constraints than ICES5G used in PMIP2. | + | Several PMIP participants are in favour of using a new ice-sheet reconstruction that better matches the geomorphologic and glaciological constraints than ICE-5G used in PMIP2. |
We plan thus to have | We plan thus to have | ||
- a reference ice-sheet for the core simulations | - a reference ice-sheet for the core simulations | ||
- alternatives for groups interested to run sensitivity experiment in order to discuss the sensitivity of the model results to the uncertainties in boundary conditions. | - alternatives for groups interested to run sensitivity experiment in order to discuss the sensitivity of the model results to the uncertainties in boundary conditions. | ||
- | We are now considering the possibility to use a new ice-sheet reconstruction proposed by Lev Tarasov and/or Dick Peltier (see text below). The final decision on the ice-sheet will be taken in September. It will depend on what is available at that time, quality checking and climate model constraints. We will involve several PMIP participants in the process, and you are all invited to comment the proposition in the [[index#page_bottom_tag|discussion panel]] (at the bottom of this page). It is still possible that we have to find an alternative and less rigorous option by the end of September, which is the deadline to finalize the boundary condition for CMIP5. | + | We are now considering the possibility to use a new ice-sheet reconstruction proposed by Lev Tarasov and/or Dick Peltier (see text below). **The final decision on the ice-sheet will be taken in September**. It will depend on what is available at that time, quality checking and climate model constraints. We will involve several PMIP participants in the process, and you are all invited to comment on the proposition in the [[index#page_bottom_tag|discussion panel]] (at the bottom of this page). It is still possible that we have to find an alternative and less rigorous option by the end of September, which is the deadline to finalize the boundary condition for CMIP5. |
From what is ready today, you can find below | From what is ready today, you can find below | ||
- | * a description of what is going on for the new reconstruction using both ice5G constraints and a glaciological model (the Lev-MOCA proposition discussed for PMIP3 deglaciation sub-project) | + | * a description of what is going on for the new reconstruction using both ICE-5G constraints and a glaciological model (the Lev-MOCA proposition discussed for PMIP3 deglaciation sub-project) |
* a new reconstruction using a new version of Dick's model (ICE-6G) | * a new reconstruction using a new version of Dick's model (ICE-6G) | ||
- | Please have a look at those files and make all the comments that can help us to finalise the CMIP5 bounary conditions rapidly. | + | Please have a look at these files and make all the comments that can help us to finalise the CMIP5 boundary conditions rapidly. |
- | ===== Details on the proposed deglacial ice-reconstructions ===== | + | ===== Details on the proposed LGM and deglacial ice-reconstructions ===== |
==== Lev-MOCA ==== | ==== Lev-MOCA ==== | ||
- | === Northern Hemispheric components (GLAC-1) === | ||
- | The Greenland reconstruction is the GrB model described in Tarasov and Peltier( GJI 2002 and JGR 2003). It's a hand tuning of a 3D thermo-mechanically coupled glaciological shallow ice-sheet model and visco-elastic bedrock response model (using the same VM2 rheology used for ICE-5G) to fit the Relative Sea-Level (RSL) observations, the GRIP borehole temperature and age profiles, and approximate basal temperatures at CampCentury and Dye3 along with secondary constraints described in the references. | + | === Northern Hemisphere (GLAC-1) === |
- | The North American and Eurasian reconstructions are objective Bayesian calibrations of the MUN/UofT glacial systems model (above 3D model plus asynchronously coupled down-slope surface drainage/lake depth solver, and various other components such as a thermodynamic lake ice, sub-glacial till-deformation, ice-shelf represention,..., some of which are described in Tarasov and Peltier, QSR 2004, and Nature 2005, and a more complete description is currently being written up). The calibration involves approximately 30 ensemble parameters to capture uncertainties in deglacial uncertainties in deglacial climate and ice dynamics. Model runs cover a full glacial cycle. The climate forcing for the calibration is based on interpolation between a glacial state (based on a combination of PMIPI for North America and PMIP II for Eurasia) LGM model output temperature and precipitation climatologies along with further regional precipitation controls) and a present-day climatology. The interpolation between glacial and present-day states is weighted to the regional temperature chronology from the tuned Greenland reconstruction. Calibration targets include a large set of RSL observations, geologically-inferred deglacial ice-margin chronologies, and geodetic constraints. For the case of North America, the calibrated ensemble is further scored with respect to strand-lines (paleo lake level indicators) and Marine Limit (maximum level of marine inundation) observations. | + | The Greenland model is from Tarasov and Peltier (2002 and 2003), a glaciological model with hand-tuned climate adjustments to enforce fit to Relative Sea-Level (RSL) records and the GRIP borehole temperature record. It was also validated against observed rated of present day uplift for 3 sites and against GPS measurements for horizontal ice surface velocity. A variant of it is the Greenland component of ICE-5G and ICE-6G. |
- | The data sets for North America and Eurasia will include surface elevation, ice mask, lake depth, surface drainage routing pointer fields and ice flux components of the surface drainage. | + | The North American and Eurasian reconstructions are objective Bayesian calibrations of the MUN/UofT glacial systems model. The latter incorporates a 3D thermo-mechanically coupled (shallow) ice-sheet model, with permafrost resolving bed-thermal model, asynchronously coupled down-slope surface drainage/lake depth solver, and various other components such as a thermodynamic lake ice, sub-glacial till-deformation, bouyancy and temperature dependant calving law, ice-shelf represention, ..., some of which are described in Tarasov and Peltier, QSR 2004, and Nature 2005, and a more complete description is currently being written up). The visco-elastic bedrock response uses either the VM2 (as used in ICE-5G) or VM5a (used in ICE-6G) earth rheologies. Relative Sea Level is computed using a gravitationally self-consistent formalism similar to that of Peltier, except for an eustatic approximation for dealing with changing ocean masks and the lack of accounting for rotational effects (which are mostly significant for far-field RSL records, ie records that do not locally constrain ice load history). |
+ | Climate forcing involves an interpolation between present day observed climatologies and the set of highest resolution LGM fields from PMIP I and II data sets. The interpolation is weighted according to a glaciological inversion of the GRIP record for regional temperatures over the last glacial cycle. | ||
+ | |||
+ | The calibration involves approximately 30 (currently 36 for North Am, 29 for Eurasia) ensemble parameters to capture uncertainties in deglacial climate and ice dynamics. The majority of these parameters are used for the climate forcing, including weighting the inter-model (ie between PMIP models) EOFs for LGM monthly precip and temperature, regional desert elevation effects, and LGM atmospheric lapse rate. Other ensemble parameters adjust calving response, effective viscosity of subglacial till, strength of margin forcing, and flow parameters for ice-shelves. Model runs are forced to stay within uncertainties of the independently derived ice margin chronologies (Dyke, 2004 for North Am, Gyllencreutz et al, in preparation for Eurasia). | ||
+ | |||
+ | Model runs cover a full glacial cycle. North America and Eurasia are calibrated separately. Calibration targets include a large set of RSL observations, geologically-inferred deglacial ice-margin chronologies, and geodetic constraints. For the case of North America, the calibrated ensemble is further scored with respect to strand-lines (paleo lake level indicators) and Marine Limit (maximum level of marine inundation) observations. A key point is that model runs are penalized in proportion to the amount of margin forcing required. So the calibration is directed towards a climate forcing that is consistent with the margin chronology. | ||
+ | |||
+ | The model was calibrated using the ICE4G ice load reconstruction for Antarctica and the VM2 earth rheology because the ICE6G Antarctic chronology and VM5a earth model along with a much expanded geodetic dataset were provided by Dick Peltier only in early September, which left too little time to recalibrate the models. There is the added issue that the ICE6G Antarctic chronology lacks error bars. The expanded geodetic data-set for North America included significant revisions to the previous geodetic constraints. This along with the significant reduction in LGM ice volume in ICE6G Antarctic as compared to ICE4 and 5G rendered a significant misfit with the far-field Barbados RSL record. With the limited time, a somewhat blind and largely random 2000 member ensemble was generated along with a rerun of the best 300 previously calibrated parameter sets and some 200 attempts at hand-tuning. nn450 is the weighted distribution of 7 model runs that passed certain hard threshold constraints. nn9021 is the best (though "best" depends to a certain extent on the weighting between various constraints) single run from the previous calibration and nn445 is the weighted ensemble mean for that previous calibration. | ||
+ | |||
+ | The Eurasian calibration did converge, and aside from issues with the Norwegian fjords (the latter are also a problem for ICE6G), the calibration was generally successful. nn8234 is one of the best runs with the largest 26ka RSL contribution to the Barbados record. A single run was chosen to ensure consistency between drainage fields and the surface topography. The mean distribution for the calibration can be made available upon request. | ||
+ | |||
+ | In summary, the GLAC-1 submission provides a set of glaciological models that are derived from a plausible climate forcing based on PMIP1 and PMIP2 results for LGM and that fit independently derived ice margin chronologies. This provides strong constraints throughout deglaciation. For Eurasia, the smaller set of constraints (among other more speculative reasons) resulted in a successful calibration that reasonably well covered the available constraint set. North America has a much larger and much more diverse set of constraints (and I suspect a much more complicated ice/climate interaction history), so the calibration has never been able to fully satisfy the whole set of constraints (strandlines are for instance a challenge to fit given their high sensitivity to drainage choke point elevations). | ||
+ | |||
+ | Unfortunately, these glaciological models in combination with the ICE-6G chronology for Antarctica (and Patagonia) and Dick Peltier's VM5a earth rheology have at best a weak fit to the the LGM segment of the Barbados record. There is a significant tradeoff between Barbados fit and fit to other constraints. What is unclear at this stage is the extent to which this is due to deficiencies in the glaciological models, to problems with the ICE-6G Antarctic ice chronology, or possibly with inferred uncertainties in the Barbados record and with the VM5a earth rheology. | ||
+ | |||
+ | One possibility for resolving Barbados, is to take the 1.5 sigma upper limit of the previously calibrated ensemble for North America which almost reaches the inferred Barbados record for 26 to 21 ka. Dick Peltier and Rosemarie Drummond will cross-check this dataset. The problem with using ensemble bounds is that this is no longer a glaciologically self-consistent model and RSL fits have also deteriorated. | ||
+ | |||
+ | == References == | ||
+ | |||
+ | * [[http://www.atmosp.physics.utoronto.ca/people/lev/g6jgrpub.pdf|g6jgrpub.pdf]] | ||
+ | * [[http://www.atmosp.physics.utoronto.ca/people/lev/g5jpub.pdf|g5jpub.pdf]] | ||
+ | |||
+ | == Data == | ||
+ | |||
+ | <note warning>Note that the North Am model is calibrated with a stub Greenland by Ellesmere. When a final model is chosen, the North American and Greenland models can be easily melded.\\ \\ Also, please note, that none of the above have yet been published, so these files are provided for discussion and consideration. Any fields and data not selected as a PMIP boundary condition are not to be used for research without prior consent.</note> | ||
+ | |||
+ | * NetCDF of 21 to 10ka surface elevation, thickness, basal velocities for the nn9021 North Am model and nn8234 Eurasian model: [[http://pmip3.lsce.ipsl.fr/share/design/glac1/NAnn9021.nc|NAnn9021.nc]], [[http://pmip3.lsce.ipsl.fr/share/design/glac1/EAnn8234.nc|EAnn8234.nc]] | ||
+ | |||
+ | * NetCDF of 21 to 10ka for nn445 ensemble means and 1.5 sigma range: [[http://pmip3.lsce.ipsl.fr/share/design/glac1/NAnn445.nc|NAnn445.nc]] | ||
+ | |||
+ | * Other data available upon request from [[lev@mun.ca|Lev Tarasov]] | ||
+ | |||
+ | == Plots == | ||
+ | |||
+ | * RSL plots for North Am and Eurasia | ||
+ | * North Am: (with nn8234 Eurasia)\\ [[http://pmip3.lsce.ipsl.fr/share/design/glac1/rslGLAC1RSLcompAc.pdf|rslGLAC1RSLcompAc.pdf]] nn9021 (as calibrated, ICE4G Antarctica, VM2),nn9021 with ICE6G Antarctica/VM5a, nn7280, nn445\\ \\ [[http://pmip3.lsce.ipsl.fr/share/design/glac1/rslGLAC1RSLcompBc.pdf|rslGLAC1RSLcompBc.pdf]] nn9021, nn445, nn445+1sigma, nn445-1sigma\\ | ||
+ | * Eurasia:\\ [[http://pmip3.lsce.ipsl.fr/share/design/glac1/rslGLAC1RSLcompEc.pdf|rslGLAC1RSLcompEc.pdf]] nn8234 (with nn9021), nn8234(with nn7280), nn8191(with nn9021)\\ | ||
+ | * Barbados: [[http://pmip3.lsce.ipsl.fr/share/design/glac1/GLAC1RSLcompBbarb.pdf|GLAC1RSLcompBbarb.pdf]] (note will slightly over predict barbados RSL drop around LGM due to lack of rotational effects on RSL calculation) | ||
+ | |||
+ | * Surface elevation and basal ice velocity area maps for 21 ka\\ Note: lakes and drainage results are not included in these plots\\ [[http://pmip3.lsce.ipsl.fr/share/design/glac1/21kaHuvnn8234.pdf|21kaHuvnn8234.pdf]], [[http://pmip3.lsce.ipsl.fr/share/design/glac1/18kaHuvnn8234.pdf|18kaHuvnn8234.pdf]], [[http://pmip3.lsce.ipsl.fr/share/design/glac1/14kaHuvnn8234.pdf|14kaHuvnn8234.pdf]] | ||
=== Southern hemisphere === | === Southern hemisphere === | ||
- | Glaciological constraints have shown that the ice5G reconstruction was too high in East Antarctica. Several tests are now made using a revised version of ICE-5G (VM2=ICE-6G) to correct this aspect while keeping the global on see level and other geomorphological constraints. Dick Peltier’s group is working now on smoothing the refine model and taking the extension of the loading of the shelf out to a distance from the present coast consistent with a specific bathymetry contour. The final version should be ready by September. | + | Glaciological constraints have shown that the ICE-5G reconstruction was too high in East Antarctica. Several tests are now made using a revised version of ICE-5G (VM2=ICE-6G) to correct this aspect while keeping the global on sea level and other geomorphological constraints. Dick Peltier’s group is working now on smoothing the refined model and taking the extension of the loading of the shelf out to a distance from the present coast consistent with a specific bathymetry contour. The final version should be ready by September. |
Line 46: | Line 87: | ||
They are provided for evaluation and test, and do not correspond to the "official" boundary condition to be used for LGM simulations in CMIP5. | They are provided for evaluation and test, and do not correspond to the "official" boundary condition to be used for LGM simulations in CMIP5. | ||
- | Here are also two submitted papers where you can find information about the new ice-sheet that are provided by Dick. | + | |
+ | === References === | ||
+ | |||
+ | Note: get in touch with [[pmip2web@lsce.ipsl.fr|Jean-Yves Peterschmitt]] in you need a copy of the following papers. | ||
+ | |||
+ | * Peltier, W.R., 2009, //Closure of the budget of global sea level rise over the GRACE Era: The importance and magnitudes of the required corrections for global glacial isostatic adjustment//. Quat. Sci. Rev. 28, 1658-1674. | ||
+ | |||
+ | * Donald F. Argus and W. R. Peltier, 2009. //Constraining models of postglacial rebound using space geodesy: A detailed assessment of model ICE-5G (VM2) and its relatives//, Geophysical Journal International, submitted. | ||
+ | |||
+ | * W.R. Peltier and R. Drummond, 2008. //Rheological stratification of the lithosphere: A direct inference basd upon the geodetically observed pattern of the glacial isostatic adjustment of the North American continent//. geophys. res. Lett., vol. 35, L16314, doi:10.1029/2008GL034586. | ||
+ | |||
+ | |||
+ | === Data === | ||
+ | |||
+ | <note tip>You should download the {{:pmip3:design:21k:icesheet:ice-ng_compare_polar.pdf}} file if you want to look at detailed graphics comparing the different versions of ICE-nG (56 Mb, 77 pages)</note> | ||
Notes: | Notes: | ||
+ | * Get in touch with [[pmip2web@lsce.ipsl.fr|Jean-Yves Peterschmitt]] if you want to download the netcdf files used to generate the figures below | ||
* Click on a figure to get the full size version | * Click on a figure to get the full size version | ||
- | * You can compare the figures with the [[http://pmip2.lsce.ipsl.fr/design/ice5g /ice5g_v1_1.shtml|ICE-5G v1.1 (PMIP2)]] and [[http://pmip2.lsce.ipsl.fr/design/ice5g/index.shtml|ICE-5G v1.2]] versions avalaible on the PMIP2 web site | + | * You can compare the figures with the ICE-5G versions available on the PMIP2 web site |
+ | *[[http://pmip2.lsce.ipsl.fr/design/ice5g/ice5g_v1_1.shtml|ICE-5G v1.1 (PMIP2)]] | ||
+ | * [[http://pmip2.lsce.ipsl.fr/design/ice5g/index.shtml|ICE-5G v1.2]] | ||
{{:pmip3:design:21k:icesheet:ice6g_v1.0_00.0k_1deg.png?200| 0 K}} | {{:pmip3:design:21k:icesheet:ice6g_v1.0_00.0k_1deg.png?200| 0 K}} | ||
Line 57: | Line 115: | ||
{{:pmip3:design:21k:icesheet:ice6g_v1.0_21.0k_1deg.png?200|21 K}} | {{:pmip3:design:21k:icesheet:ice6g_v1.0_21.0k_1deg.png?200|21 K}} | ||
{{:pmip3:design:21k:icesheet:ice6g_v1.0_26.0k_1deg.png?200|26 K}} | {{:pmip3:design:21k:icesheet:ice6g_v1.0_26.0k_1deg.png?200|26 K}} | ||
+ | |||
+ | === Ice volume === | ||
+ | |||
+ | The following table has been supplied by | ||
+ | [[rmarie@atmosp.physics.utoronto.ca|Rosemarie Drummond]]. It compares | ||
+ | the ice amounts in the various regions in terms of eustatic sea level | ||
+ | impact inferred by assuming the ocean area to remain fixed to the | ||
+ | modern area (//ice amount relative to now in m of ocean assuming ocean | ||
+ | area = 360768576 km<sup>2</sup>//). | ||
+ | |||
+ | ^ ^ ICE-4G\\ //final// (1) ^^ ICE-5G\\ //v1.2// (2) ^^ ICE-6G\\ //v1.0// ^^ | ||
+ | ^ ^ 21k ^ 26k ^ 21k ^ 26k ^ 21k ^ 26k ^ | ||
+ | ^ North America\\ (incl Innuit area) | 64.24 | 54.92 | 81.47 | 83.71 | 79.82 | 88.14 | | ||
+ | ^ Greenland & Iceland | 6.38 | 5.43 | 2.49 | 2.45 | 2.36 | 2.34 | | ||
+ | ^ Fennoscandia | 10.39 | 8.91 | 11.19 | 11.79 | 10.22 | 12.31 | | ||
+ | ^ Barents / Kara Seas | 14.05 | 12.26 | 8.43 | 9.29 | 7.26 | 9.10 | | ||
+ | ^ UK | 0.42 | 0.35 | 1.48 | 1.65 | 0.58 | 0.82 | | ||
+ | ^ Patagonia | 0.55 | 0.47 | 0.55 | 0.55 | 0.83 | 0.87 | | ||
+ | ^ West Antarctica | 9.74 | 8.33 | 9.68 | 9.68 | 11.79 | 11.79 | | ||
+ | ^ East Antarctica | 8.35 | 7.12 | 8.36 | 8.36 | 1.45 | 1.44 | | ||
+ | ^ Total | **114.12** | **97.79** | **123.63** | **127.48** | **114.31** | **126.81** | | ||
+ | ^ ice-shelves (3) | | | | | 2.80 | 3.00 | | ||
+ | |||
+ | * (1) The version of the ICE-4G (VM2) model for which data are provided is not exactly the version that was employed in the original PMIP project but rather one from which a significant amount of Antarctic ice was simply eliminated (version employed for illustrative purposes in Peltier's paper for the Mt Hood volume of QSR). | ||
+ | |||
+ | * (2) The version of ICE-5G used for PMIP2/MOTIF was **v1.1** | ||
+ | |||
+ | * (3) The calculation only cares about grounded ice, so where it detects that the ice-model has ice that is below the predicted sea-level, it deems the ice to be melted and added to the ocean. The ice may be melted or be part of an ice-shelf eg the Ronne, Ross ice-shelves. This is a feature of the most recent version of the solver. | ||
+ | |||
+ | * The large change between ICE-4G (VM2) and ICE-5G (VM2) is due to the introduction of Art Dyke's margin chronology which had the ice being eliminated from North America earlier than in the original chronology obtained from John Andrews. This meant that more ice had to be removed in order to have the same impact on RSL during the ice-free Holocene interval. | ||
+ | |||
+ | * The large differences between ICE-5G (VM2) and ICE-6G (VM5a) are a combined result of the introduction of two further constraints: | ||
+ | * the analysis from Valerie Masson demonstrating that the elevation of East Antarctic ice was not significantly different than modern | ||
+ | * the results of the analyses based upon the application of the space geodetic constraints described in the Argus and Peltier paper. | ||
+ | |||
====== Conclusion ====== | ====== Conclusion ====== | ||
- | In September we should be able to compare the results of the revised version of ICE-45G with those of the Lev-Dick collaboration and take a decision based on comparisons with data and on the advice of several specialists | + | In September we should be able to compare the results of the revised version of ICE-6G with those of the Lev-Dick collaboration and take a decision based on comparisons with data and on the advice of the community |
~~DISCUSSION~~ | ~~DISCUSSION~~ | ||
/* standard page footer */ | /* standard page footer */ | ||
Line 69: | Line 163: | ||
[ [[pmip3:|PMIP3 Wiki Home]] ] - | [ [[pmip3:|PMIP3 Wiki Home]] ] - | ||
[ [[pmip3:wiki_help|Help!]] ] - | [ [[pmip3:wiki_help|Help!]] ] - | ||
- | [ [[wiki:syntax|Wiki syntax]] ] | + | [ [[wiki:syntax|Wiki syntax]] ] <html><a name="page_bottom_tag"> </a></html> - [ [[index#page_top_tag|Top]] ] |