Antbounds Contents and data masking tips

The functions in this plugin for Antarctic Mapping Tools for Matlab are intended to simplify the process of Antarctic data masking and mapmaking. These functions are built on the MEaSURES Antarctic Boundaries for IPY 2007-2009 from Satellite Radar dataset Version 2, which has been kindly provided by Mouginot et al. (see citation below) and is described in full on the NSIDC site here. All necessary data are included in this File Exchange submission, so you don't need to download anything from the NSIDC.


Data Access

If you have downloaded this toolbox, you have all the data you need. These functions make it easy to import the grounding lines, coast lines, ice shelf outlines, and ice basins into Matlab:

Interpolation and masking functions

This toolbox contains five functions which are optimized to load Mouginot et al.'s Mask_Antarctica data and and interpolate to determine whether locations correspond to grounded ice, ice shelves, open ocean, etc. The inbasin function works the same way, and is convenient for limiting your dataset to contain only values within a given ice basin. In my work, I've found these functions and dist2mask have become absolutely invaluable. I hope they're just as helpful for you.

Data plotting

These functions make plotting pretty easy:

Tips for data masking

If you're working with a dataset that comes with its own mask (e.g., Bedmap2, RTopo-2, etc.), it might make sense to use the masks associated with those datasets. Otherwise, the logical is* functions listed above are in many cases the easiest to work with. For example, consider this 50 km resolution grid over Dronning Maud Land, which we create with psgrid:

[lat,lon] = psgrid('dronning maud land',3000,50);

Plot the grid, a grounding line, and a coast line for context:

plotps(lat,lon,'.','color',0.8*[1 1 1])
hold on
axis tight

To determine which grid points correspond to open ocean, use isopenocean:

ocean = isopenocean(lat,lon);

Perhaps your work is on Baudouin Ice Shelf. Outline it with a thick red line:


How many grid cells correspond to Baudouin Ice Shelf?

baudouin = isiceshelf(lat,lon,'Baudouin');
ans =

Okay, so 13 grid cells correspond to Baudouin Ice Shelf. Plot them in as large red x marks:


And we can add some junk to the map if we so desire:


Never heard of Baudouin Ice Shelf? Well guess how big it is:

[x,y] = antbounds_data('baudouin','xy');

pgon = polyshape(x/1000,y/1000,'Simplify',false); % req's R2017b+
ans =

So 33,000 km^2, or about the size of Moldova. Similarly, we can see how large of an area each major glacier drains using the basin_data function. For example, Totten:

[x_tot,y_tot] = basin_data('imbie refined','totten','xy');

axis image % removes whitespace, sets equal aspect ratio

And Totten drains an area of

pgon = polyshape(x_tot/1000,y_tot/1000,'Simplify',false); % req's R2017b+
ans =

About 550,000 km^2, or about the size of Spain. What's the average ice thickness in the Totten basin? We can load the bedmap2 data surrounding the Totten basin like this:

[lat,lon,thickness] = bedmap2_data('thickness',x_tot,y_tot);

Determine which Bedmap2 grid cells are drained by Totten (takes a few seconds):

totten = inbasin(lat,lon,'imbie refined','totten');

The Bedmap2 dataset is 1 km by 1 km, so the total number of grid cells should be about the same as the drainage basin area we calculated above:

ans =

And that looks like we're on the right track. The average ice thickness drained by Totten Glacier is simply

ans =

About 3.2 km thick.

We can mask out everything except the ice drained by Totten and plot:

thickness(~totten) = NaN;
cb = colorbar;
ylabel(cb,'ice thickness (m)')

Citing this dataset

If you use this dataset, please cite the following:

Author Info

This toolbox and supporting documentation were written by Chad A. Greene of the University of Texas Institute for Geophysics (UTIG), November 2016. Updated May 2017 for version 2 of the dataset.