The poles' topography depends on the presence, extent, and thickness of ice sheets and their associated isostatic effects. If the reconstruction corresponds to a warmer period, one must isostatically compensate for the ice loss on Greenland or Antarctica. This tool compensates for the absence of ice mass in the poles. The DEM values will be modified in the areas where there is a noticeable difference between ice and bedrock topography. The percentage of ice to be removed from the polar regions must be defined. The higher it is, the stronger the isostatic rebound, and the higher the resulting topography will be.
This tool uses present-day ice and bedrock topography to calculate the thickness of present-day ice covering Antarctica and Greenland. The thickness is determined by subtracting the bedrock topography from ice topography. Then that topography is raised by ⅓ of the current ice thickness if all the ice is removed. If only part of the ice is to be removed, the following formula applies:
Hc = Hi+0.3*Hd*Ir
Hc - Compensated elevation (new elevation value)
Hi - Initial elevation
Hd - Elevation difference between the bedrock and ice topographies (ice thickness)
Ir - Amount of ice removed (in %)
Note that this tool does not generate any map of the ice topography or thickness. It only lifts the bedrock topography to compensate for the loss of ice mass.
Step by step
Click on the icon
In the field How would you like to process the input DEM, choose Isostatic Compensation.
In Raster to be modified, choose the bedrock topography DEM (rotated to the time of reconstruction).
In Ice topography raster, choose the ice topography DEM (rotated to the time of reconstruction).
Enter the percentage of ice that is removed from the poles to calculate the isostatic rebound. The higher this percentage, the stronger the rebound, and the higher the resulting bedrock topography.
If you don’t want this tool to work in the entire map, you have the option to constrain the isostatic compensation to particular regions. You need to select a Mask layer. If you don’t want to use all polygons in the Mask layer, there are two ways to choose which ones to use.
You may select the polygons you want the compensation to be constrained to on the map canvas. In this case, check the box Selected features only.
If the regions you want to limit the compensation to are Greenland and Antarctica, you don’t need to select manually the tectonic plates. You can check the box Get polar regions automatically. In this case, the tool will search for the names of polar continental plates in the mask layer automatically and use only these for the isostatic compensation. Make sure your shapefile has an attribute field called NAME with the names of the tectonic plates.
This second option is designed to be used with the continental blocks’ shapefile. It is important that the NAME attribute field is present.
The table below shows the names of polar tectonic plates the algorithm recognizes. Every plate whose name contains any of the words listed in the left column will be used for the isostatic compensation (case insensitive). The column on the right shows the most common spelling for those plate names.
- Choose the Output file path, and click Run.