polartoolkit.fetch

class EarthDataDownloader

Bases: object

Adapted from IcePack: icepack/icepack Either pulls login details from pre-set environment variables, or prompts user to input username and password.

basal_melt(variable='w_b')

Antarctic ice shelf basal melt rates for 1994-2018 from satellite radar altimetry. from Adusumilli et al.[1].

accessed from http://library.ucsd.edu/dc/object/bb0448974g

reading files and preprocessing from supplied jupyternotebooks: sioglaciology/ice_shelf_change

Units are in m/yr

Parameters:

variable (str) – choose which variable to load, either ‘w_b’ for basal melt rate, ‘w_b_interp’, for basal melt rate with interpolated values, and ‘w_b_uncert’ for uncertainty

Returns:

Returns a dataarray of basal melt rate values

Return type:

xr.DataArray

References

[1]

Susheel Adusumilli, Helen Amanda Fricker, Brooke Medley, Laurie Padman, and Matthew R. Siegfried. Interannual variations in meltwater input to the Southern Ocean from Antarctic ice shelves. Nature Geoscience, 13(9):616–620, September 2020. doi:10.1038/s41561-020-0616-z.

basement(region=None, spacing=None, registration=None)

Load a grid of basement topography. Offshore and sub-Ross Ice Shelf basement topography. from Tankersley et al.[2], Tankersley et al.[3] and Lindeque et al.[4].

Elevations are referenced to WGS84 ellipsoid.

Parameters:

• region (tuple[float, float, float, float], optional) – bounding region to clip the loaded grid to, by default doesn’t clip

• spacing (float, optional) – grid spacing to resample the loaded grid to, by default 5e3

• registration (str, optional) – change registration with either ‘p’ for pixel or ‘g’ for gridline registration, by default is None.

Returns:

dataarray of basement depths

Return type:

xr.DataArray

References

[2] (1,2)

M. Tankersley, H. J. Horgan, C. S. Siddoway, F. Caratori Tontini, and K. J. Tinto. Basement Topography and Sediment Thickness Beneath Antarctica’s Ross Ice Shelf. Geophysical Research Letters, May 2022. doi:10.1029/2021GL097371.

[3] (1,2)

Matthew Tankersley, Huw J. Horgan, Christine S. Smith Siddoway, Fabio Caratori Tontini, and Kirsty Tinto. Basement topography and sediment thickness beneath Antarctica’s Ross Ice Shelf imaged with airborne magnetic data. 2022. doi:10.1594/PANGAEA.941238.

[4] (1,2)

Ansa Lindeque, Karsten Gohl, Florian Wobbe, and Gabriele Uenzelmann-Neben. Pre-glacial to glacial sediment thickness grids for the Southern Pacific Margin of West Antarctica, NetCDF files. 2016. doi:10.1594/PANGAEA.864906.

bedmachine(layer, reference='eigen-6c4', region=None, spacing=None, registration=None, **kwargs)

Load BedMachine v3 data, from Morlighem[5].

Accessed from NSIDC via https://nsidc.org/data/nsidc-0756/versions/3. Also available from ldeo-glaciology/pangeo-bedmachine

Referenced to the EIGEN-6C4 geoid. To convert to be ellipsoid-referenced, we add the geoid grid. use reference=’ellipsoid’ to include this conversion in the fetch call.

Surface and ice thickness are in ice equivalents. Actual snow surface is from REMA [6], and has had firn thickness added(?) to it to get Bedmachine Surface.

To get snow surface: surface+firn To get firn and ice thickness: thickness+firn

Here, icebase will return a grid of surface-thickness This should be the same as snow-surface - (firn and ice thickness)

Parameters:

• layer (str) – choose which layer to fetch: ‘bed’, ‘dataid’, ‘errbed’, ‘firn’, ‘geoid’, ‘mapping’, ‘mask’, ‘source’, ‘surface’, ‘thickness’; ‘icebase’ will give results of surface-thickness

• reference (str) – choose whether heights are referenced to ‘eigen-6c4’ geoid or the ‘ellipsoid’ (WGS84), by default is eigen-6c4’

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (str or int, optional) – grid spacing to resample the loaded grid to, by default 500m

• registration (str, optional) – change registration with either ‘p’ for pixel or ‘g’ for gridline registration, by default is None.

• kwargs (Any) – additional kwargs to pass to resample_grid

Returns:

Returns a loaded, and optional clip/resampled grid of Bedmachine.

Return type:

xr.DataArray

References

[5]

M. Morlighem. MEaSUREs BedMachine Antarctica, Version 3. 2022.

[6] (1,2)

Ian M. Howat, Claire Porter, Benjamin E. Smith, Myoung-Jong Noh, and Paul Morin. The Reference Elevation Model of Antarctica. The Cryosphere, 13(2):665–674, February 2019. doi:10.5194/tc-13-665-2019.

bedmap2(layer, reference='eigen-gl04c', region=None, spacing=None, registration=None, fill_nans=False, **kwargs)

Load bedmap2 data as xarray.DataArrays from Fretwell et al.[7]. accessed from https://ramadda.data.bas.ac.uk/repository/entry/show?entryid=fa5d606c-dc95-47ee-9016-7a82e446f2f2.

All grids are by default referenced to the EIGEN-GL04C geoid. Use the reference=’ellipsoid’ to convert to the WGS-84 ellipsoid or reference=’eigen-6c4’ to convert to the EIGEN-6c4 geoid.

Unlike Bedmachine data, Bedmap2 surface and icethickness contain NaN’s over the ocean, instead of 0’s. To fill these NaN’s with 0’s, set fill_nans=True. Note, this only makes since if the reference is the geoid, therefore, if reference=’ellipsoid and fill_nans=True, the nan’s will be filled before converting the results to the geoid (just for surface, since thickness isn’t relative to anything).

Parameters:

• layer (str) – choose which layer to fetch: “bed”, “coverage”, “grounded_bed_uncertainty”, “icemask_grounded_and_shelves”, “lakemask_vostok”, “rockmask”, “surface”, “thickness”, “thickness_uncertainty_5km”, “gl04c_geiod_to_WGS84”, “icebase”

• reference (str) – choose whether heights are referenced to the ‘eigen-6c4’ geoid, the WGS84 ellipsoid, ‘ellipsoid’, or by default the ‘eigen-gl04c’ geoid.

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (str or int, optional) – grid spacing to resample the loaded grid to, by default 10e3

• registration (str, optional,) – choose between ‘g’ (gridline) or ‘p’ (pixel) registration types, by default is the original type of the grid

• fill_nans (bool, optional,) – choose whether to fill nans in ‘surface’ and ‘thickness’ with 0. If converting to reference to the geoid, will fill nan’s before conversion, by default is False

• kwargs (Any) – additional kwargs to pass to resample_grid

Returns:

Returns a loaded, and optional clip/resampled grid of Bedmap2.

Return type:

xr.DataArray

References

[7] (1,2)

P. Fretwell, H. D. Pritchard, D. G. Vaughan, J. L. Bamber, N. E. Barrand, R. E. Bell, C. Bianchi, R. G. Bingham, D. D. Blankenship, G. Casassa, G. Catania, D. Callens, H. Conway, A. J. Cook, H. F. J. Corr, D. Damaske, V. Damm, F. Ferraccioli, R. Forsberg, S. Fujita, Y. Gim, P. Gogineni, J. A. Griggs, R. C. A. Hindmarsh, P. Holmlund, J. W. Holt, R. W. Jacobel, A. Jenkins, W. Jokat, T. A. Jordan, E. C. King, J. Kohler, W. Krabill, M. Riger-Kusk, K. A. Langley, G. Leitchenkov, C. Leuschen, B. Luyendyk, K. Matsuoka, J. Mouginot, F. O. Nitsche, Y. Nogi, O. A. Nost, S. V. Popov, E. Rignot, D. M. Rippin, A. Rivera, J. Roberts, N. Ross, M. J. Siegert, A. M. Smith, D. Steinhage, M. Studinger, B. Sun, B. K. Tinto, B. C. Welch, D. Wilson, D. A. Young, C. Xiangbin, and A. Zirizzotti. Bedmap2: improved ice bed, surface and thickness datasets for Antarctica. The Cryosphere, 7(1):375–393, 2013. doi:10.5194/tc-7-375-2013.

bedmap_points(version, region=None)

Load bedmap point data, choose from Bedmap 1, 2 or 3

version == ‘bedmap1’

from Lythe and Vaughan[8]. accessed from https://data.bas.ac.uk/full-record.php?id=GB/NERC/BAS/PDC/01619

version == ‘bedmap2’

from Fretwell et al.[7]. accessed from https://data.bas.ac.uk/full-record.php?id=GB/NERC/BAS/PDC/01616

version == ‘bedmap3’

from Fremand et al.[9]. accessed from https://data.bas.ac.uk/full-record.php?id=GB/NERC/BAS/PDC/01614

Parameters:

• version (str) – choose between ‘bedmap1’, ‘bedmap2’, or ‘bedmap3’ point data

• region (tuple[float, float, float, float], optional) – add a GMT region to subset the data by, by default None

Returns:

Return a dataframe, optionally subset by a region

Return type:

pd.DataFrame

References

[8]

Matthew B. Lythe and David G. Vaughan. BEDMAP: A new ice thickness and subglacial topographic model of Antarctica. Journal of Geophysical Research: Solid Earth, 106(B6):11335–11351, 2001. doi:10.1029/2000JB900449.

[9]

Alice Fremand, Peter Fretwell, Julien Bodart, Hamish Pritchard, A. Aitken, Jonathan Bamber, R.E. Bell, C Bianchi, Robert Bingham, Donald Blankenship, G. Casassa, Ginny Catania, K Christianson, H Conway, Hugh Corr, X. Cui, D Damaske, V. Damn, R. Drews, G. Eagles, O. Eisen, H. Eisermann, Fausto Ferraccioli, Rene Forsberg, S Franke, S. Fujita, Y. Gim, V. Goel, P. Gogineni, J.S. Greenbaum, B. Hills, Richard Hindmarsh, P. Holmlund, N. Holschuh, J.W. Holt, A. Humbert, R. Jacobel, D. Jansen, Adrian Jenkins, Wilfried Jokat, Tom Jordan, Edward King, Jack Kohler, W. Krabill, K. Langley, J. Lee, G. Leitchenkov, C. Leuschen, B. Luyendyk, J. MacGregor, E. MacKie, Kenichi Matsuoka, Mathieu Morlighem, Jeremie Mouginot, Frank Nitsche, Y. Nogi, O. Nost, John Paden, F. Pattyn, S. Popov, M. Riger-Kusk, Eric Rignot, David Rippin, Andres Rivera, Jason Roberts, Neil Ross, A. Ruppel, Dustin Schroeder, Martin Siegert, Andrew Smith, D. Steinhage, M. Studinger, B. Sun, I. Tabacco, Kirsty Tinto, S. Urbini, David Vaughan, B. Welch, D.S. Wilson, Duncan Young, and A. Zirizzotti. BEDMAP3 - Ice thickness, bed and surface elevation for Antarctica - standardised data points. October 2022. doi:10.5285/91523FF9-D621-46B3-87F7-FFB6EFCD1847.

crustal_thickness(version, region=None, spacing=None, registration=None)

Load 1 of x ‘versions’ of Antarctic crustal thickness grids.

version=’shen-2018’ Crustal thickness (excluding ice layer) from Shen et al.[10]. Accessed from https://sites.google.com/view/weisen/research-products?authuser=0

version=’an-2015’ Crustal thickness (distance from solid (ice and rock) top to Moho discontinuity) from An et al.[11]. Accessed from http://www.seismolab.org/model/antarctica/lithosphere/index.html#an1s File is the AN1-CRUST model, paper states “Moho depths and crustal thicknesses referred to below are the distance from the solid surface to the Moho. We note that this definition of Moho depth is different from that in the compilation of AN-MOHO”. Unclear, but seems moho depth is just negative of crustal thickness. Not sure if its to the ice surface or ice base.

Parameters:

• version (str) – Either ‘shen-2018’, will add later: ‘lamb-2020’, ‘an-2015’, ‘baranov’, ‘chaput’, ‘crust1’, ‘szwillus’, ‘llubes’, ‘pappa’, ‘stal’

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (int, optional) – grid spacing to resample the loaded grid to, by default spacing is read from downloaded files

• registration (str, optional) – change registration with either ‘p’ for pixel or ‘g’ for gridline registration, by default is None.

Returns:

Returns a loaded, and optional clip/resampled grid of crustal thickness.

Return type:

xr.DataArray

References

[10] (1,2)

Weisen Shen, Douglas A. Wiens, Sridhar Anandakrishnan, Richard C. Aster, Peter Gerstoft, Peter D. Bromirski, Samantha E. Hansen, Ian W. D. Dalziel, David S. Heeszel, Audrey D. Huerta, Andrew A. Nyblade, Ralph Stephen, Terry J. Wilson, and J. Paul Winberry. The crust and upper mantle structure of Central and West Antarctica from bayesian inversion of Rayleigh wave and receiver functions. Journal of Geophysical Research: Solid Earth, 123(9):7824–7849, 2018. doi:10.1029/2017JB015346.

[11] (1,2)

Meijian An, Douglas A. Wiens, Yue Zhao, Mei Feng, Andrew A. Nyblade, Masaki Kanao, Yuansheng Li, Alessia Maggi, and Jean-Jacques Lévêque. S-velocity model and inferred Moho topography beneath the Antarctic Plate from Rayleigh waves: Antarctic S-velocities and Moho. Journal of Geophysical Research: Solid Earth, 120(1):359–383, 2015. doi:10.1002/2014JB011332.

deepbedmap(region=None, spacing=None, registration=None)

Load DeepBedMap data, from Leong and Horgan[12] and Leong and Horgan[13].

Parameters:

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (str or int, optional) – grid spacing to resample the loaded grid to, by default 10e3

• registration (str, optional) – change registration with either ‘p’ for pixel or ‘g’ for gridline registration, by default is None.

Returns:

Returns the filepath of DeepBedMap.

Return type:

str

References

[12]

Wei Ji Leong and Huw Joseph Horgan. DeepBedMap: a deep neural network for resolving the bed topography of Antarctica. The Cryosphere, 14(11):3687–3705, November 2020. doi:10.5194/tc-14-3687-2020.

[13]

Wei Ji Leong and Huw Joseph Horgan. DeepBedMap: A super-resolution neural network created bed topography of Antarctica. November 2020. doi:10.5281/zenodo.4054246.

etopo(region=None, spacing=None, registration=None)

Loads a grid of Antarctic topography from ETOPO1 from ETOPO1 1 Arc-Minute Global Relief Model: Procedures, Data Sources and Analysis[14]. Originally at 10 arc-min resolution, reference to mean sea-level.

originally from https://www.ncei.noaa.gov/access/metadata/landing-page/bin/iso?id=gov.noaa.ngdc.mgg.dem:316 Accessed via the Fatiando data repository fatiando-data/earth-topography-10arcmin

Parameters:

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (str or int, optional) – grid spacing to resample the loaded grid to, by default 10e3

• registration (str, optional) – change registration with either ‘p’ for pixel or ‘g’ for gridline registration, by default is None.

Returns:

Returns a loaded, and optional clip/resampled grid of topography.

Return type:

xr.DataArray

References

[14]

ETOPO1 1 Arc-Minute Global Relief Model: Procedures, Data Sources and Analysis. 2009. doi:10.7289/V5C8276M.

geoid(region=None, spacing=None, registration=None, **kwargs)

Loads a grid of Antarctic geoid heights derived from the EIGEN-6C4 from Förste et al.[15] spherical harmonic model of Earth’s gravity field. Originally at 10 arc-min resolution. Negative values indicate the geoid is below the ellipsoid surface and vice-versa. To convert a topographic grid which is referenced to the ellipsoid to be referenced to the geoid, add this grid. To convert a topographic grid which is referenced to the geoid to be reference to the ellipsoid, subtract this grid.

originally from https://dataservices.gfz-potsdam.de/icgem/showshort.php?id=escidoc:1119897 Accessed via the Fatiando data repository fatiando-data/earth-geoid-10arcmin

Parameters:

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (str or int, optional) – grid spacing to resample the loaded grid to, by default 10e3

• registration (str, optional) – change registration with either ‘p’ for pixel or ‘g’ for gridline registration, by default is None.

• kwargs (Any) – additional kwargs to pass to resample_grid

Returns:

Returns a loaded, and optional clip/resampled grid of geoid height.

Return type:

xr.DataArray

References

[15] (1,2)

Christoph Förste, \relax Sean.L. Bruinsma, Oleg Abrikosov, Jean-Michel Lemoine, Jean Charles Marty, Frank Flechtner, G. Balmino, F. Barthelmes, and R. Biancale. EIGEN-6C4 The latest combined global gravity field model including GOCE data up to degree and order 2190 of GFZ Potsdam and GRGS Toulouse. 2014. doi:10.5880/ICGEM.2015.1.

geomap(version='faults', region=None)

Data from GeoMAP accessed from https://doi.pangaea.de/10.1594/PANGAEA.951482?format=html#download

from Cox et al.[16] and Cox et al.[17].

Parameters:

• version (str, optional) – choose which version to retrieve, “faults”, “units”, “sources”, or “quality”, by default “faults”

• region (tuple[float, float, float, float], optional) – return only data within this region, by default None

Returns:

Returns a geodataframe

Return type:

gpd.GeoDataFrame

References

[16]

Simon C. Cox, Belinda Smith Lyttle, Samuel Elkind, Christine Smith Siddoway, Paul Morin, Giovanni Capponi, Tamer Abu-Alam, Matilda Ballinger, Lauren Bamber, Brett Kitchener, Luigi Lelli, Jasmine Mawson, Alexie Millikin, Nicola Dal Seno, Louis Whitburn, Tristan White, Alex Burton-Johnson, Laura Crispini, David Elliot, Synnøve Elvevold, John Goodge, Jacqueline Halpin, Joachim Jacobs, Adam P. Martin, Eugene Mikhalsky, Fraser Morgan, Phil Scadden, John Smellie, and Gary Wilson. A continent-wide detailed geological map dataset of Antarctica. Scientific Data, 10(1):250, May 2023. doi:10.1038/s41597-023-02152-9.

[17]

Simon Christopher Cox, Belinda Smith Lyttle, Samuel Elkind, Christine S. Smith Siddoway, Paul Morin, Giovanni Capponi, Tamer Abu-Alam, Matilda Ballinger, Lauren Bamber, Brett Kitchener, Luigi Lelli, Jasmine F. Mawson, Alexie Millikin, Nicola Dal Seno, Louis Whitburn, Tristan White, Alex Burton-Johnson, Laura Crispini, David Elliot, Synnøve Elvevold, John W. Goodge, Jacqueline A. Halpin, Joachim Jacobs, Eugene Mikhalsky, Adam P. Martin, Fraser Morgan, John Smellie, Phil Scadden, and Gary S. Wilson. The GeoMAP (v.2022-08) continent-wide detailed geological dataset of Antarctica. 2023. doi:10.1594/PANGAEA.951482.

get_fetches()

get all the fetch functions defined in this module.

Returns:

names of each fetch function

Return type:

list[str, tuple[float, float, float, float] ]

ghf(version, region=None, spacing=None, registration=None, **kwargs)

Load 1 of 6 ‘versions’ of Antarctic geothermal heat flux data.

version=’an-2015’ From An et al.[18]. Accessed from http://www.seismolab.org/model/antarctica/lithosphere/index.html

version=’martos-2017’ From Martos et al.[19] and Martos[20].

version=’shen-2020’: From Shen et al.[21]. Accessed from https://sites.google.com/view/weisen/research-products?authuser=0 Used https://paperform.co/templates/apps/direct-download-link-google-drive/ to generate a direct download link from google drive page. https://drive.google.com/uc?export=download&id=1Fz7dAHTzPnlytuyRNctk6tAugCAjiqzR

version=’burton-johnson-2020’ From Burton-Johnson et al.[22]. Accessed from supplementary material Choose for either of grid, or the point measurements

version=’losing-ebbing-2021’ From Lösing and Ebbing[23] and Lösing and Ebbing[24].

version=’aq1’ From Stål et al.[25] and Stål et al.[26].

Parameters:

• version (str) – Either ‘burton-johnson-2020’, ‘losing-ebbing-2021’, ‘aq1’,

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (int, optional) – grid spacing to resample the loaded grid to, by default spacing is read from downloaded files

• registration (str, optional) – change registration with either ‘p’ for pixel or ‘g’ for gridline registration, by default is None.

• kwargs (Any) – if version=’burton-johnson-2020’, then kwargs are passed to return point measurements instead of the grid.

Returns:

Returns a loaded, and optional clip/resampled grid of GHF data.

Return type:

xr.DataArray

References

[18]

Meijian An, Douglas A. Wiens, Yue Zhao, Mei Feng, Andrew Nyblade, Masaki Kanao, Yuansheng Li, Alessia Maggi, and Jean-Jacques Lévêque. Temperature, lithosphere-asthenosphere boundary, and heat flux beneath the Antarctic Plate inferred from seismic velocities. Journal of Geophysical Research: Solid Earth, 120(12):8720–8742, December 2015. doi:10.1002/2015JB011917.

[19]

Yasmina M. Martos, Manuel Catalán, T. A. Jordan, Alexander Golynsky, Dmitry Golynsky, Graeme Eagles, and David G. Vaughan. Heat Flux Distribution of Antarctica Unveiled. Geophysical Research Letters, 44:1–10, November 2017. doi:10.1002/2017GL075609.

[20]

Yasmina M. Martos. Antarctic geothermal heat flux distribution and estimated Curie Depths, links to gridded files. 2017. doi:10.1594/PANGAEA.882503.

[21]

Weisen Shen, Douglas A. Wiens, Andrew J. Lloyd, and Andrew A. Nyblade. A geothermal heat flux map of Antarctica empirically constrained by seismic structure. Geophysical Research Letters, 2020. doi:10.1029/2020GL086955.

[22]

Alex Burton-Johnson, Ricarda Dziadek, and Carlos Martin. Geothermal heat flow in Antarctica: current and future directions. The Cryosphere Discussions, pages 1–45, 2020. doi:10.5194/tc-2020-59.

[23]

M. Lösing and J. Ebbing. Predicting Geothermal Heat Flow in Antarctica With a Machine Learning Approach. Journal of Geophysical Research: Solid Earth, June 2021. doi:10.1029/2020JB021499.

[24]

Mareen Lösing and Jörg Ebbing. Predicted Antarctic Heat Flow and Uncertainties using Machine Learning. 2021. doi:10.1594/PANGAEA.930237.

[25]

Tobias Stål, Anya M. Reading, Jacqueline A. Halpin, and Joanne M. Whittaker. Antarctic Geothermal Heat Flow Model: Aq1. Geochemistry, Geophysics, Geosystems, February 2021. doi:10.1029/2020GC009428.

[26]

Tobias Stål, Anya M. Reading, Jacqueline A. Halpin, and Joanne Whittaker. Antarctic geothermal heat flow model: Aq1. 2020. doi:10.1594/PANGAEA.924857.

gia(version, region=None, spacing=None, registration=None)

Load 1 of 1 ‘versions’ of Antarctic glacial isostatic adjustment grids.

version=’stal-2020’ From Stål et al.[27] and Stål et al.[28].

Parameters:

• version (str) – For now the only option is ‘stal-2020’,

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (int, optional) – grid spacing to resample the loaded grid to, by default spacing is read from downloaded files

• registration (str, optional) – change registration with either ‘p’ for pixel or ‘g’ for gridline registration, by default is None.

Returns:

Returns a loaded, and optional clip/resampled grid of GIA data.

Return type:

xr.DataArray

References

[27]

Tobias Stål, Anya M. Reading, Jacqueline A. Halpin, Steven J. Phipps, and Joanne M. Whittaker. The Antarctic Crust and Upper Mantle: A Flexible 3D Model and Software Framework for Interdisciplinary Research. Frontiers in Earth Science, 8:577502, November 2020. doi:10.3389/feart.2020.577502.

[28]

Tobias Stål, Anya M. Reading, Jacqueline A. Halpin, Phipps Steven J, and Joanne M. Whittaker. The Antarctic crust and upper mantle: a flexible 3D model and framework for interdisciplinary research. Zenodo, August 2020. doi:10.5281/zenodo.4003423.

gravity(version, region=None, spacing=None, registration=None, **kwargs)

Loads 1 of 3 ‘versions’ of Antarctic gravity grids.

version=’antgg’ Antarctic-wide gravity data compilation of ground-based, airborne, and shipborne data, from Scheinert et al.[29]. Accessed from https://doi.pangaea.de/10.1594/PANGAEA.848168

version=’antgg-update’ Preliminary compilation of Antarctica gravity and gravity gradient data. Updates on 2016 AntGG compilation. Accessed from https://ftp.space.dtu.dk/pub/RF/4D-ANTARCTICA/

version=’eigen’ Earth gravity grid (eigen-6c4) at 10 arc-min resolution at 10km geometric (ellipsoidal) height from Förste et al.[15]. originally from https://dataservices.gfz-potsdam.de/icgem/showshort.php?id=escidoc:1119897 Accessed via the Fatiando data repository fatiando-data/earth-gravity-10arcmin

Parameters:

• version (str) – choose which version of gravity data to fetch.

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (str or int, optional) – grid spacing to resample the loaded grid to, by default 10e3

• registration (str, optional) – change registration with either ‘p’ for pixel or ‘g’ for gridline registration, by default is None.

• kwargs (Any) – additional kwargs to pass to resample_grid and set the anomaly_type.

Keyword Arguments:

anomaly_type (str) – either ‘FA’ or ‘BA’, for free-air and bouguer anomalies, respectively. For antgg-update can also be ‘DG’ for gravity disturbance, or ‘Err’ for error estimates.

Returns:

Returns a loaded, and optional clip/resampled grid of either observed, free-air or Bouguer gravity anomalies.

Return type:

xr.DataArray

References

[29]

M. Scheinert, F. Ferraccioli, J. Schwabe, R. Bell, M. Studinger, D. Damaske, W. Jokat, N. Aleshkova, T. Jordan, G. Leitchenkov, D. D. Blankenship, T. M. Damiani, D. Young, J. R. Cochran, and T. D. Richter. New Antarctic gravity anomaly grid for enhanced geodetic and geophysical studies in Antarctica. Geophysical Research Letters, 43(2):600–610, 2016. doi:10.1002/2015GL067439.

groundingline(version='depoorter-2013')

Load the file path of two versions of groundingline shapefiles

version = “depoorter-2013” from Depoorter et al.[30]. Supplement to Depoorter et al.[31].

version = “measures-v2” from Mouginot et al.[32]. accessed at https://nsidc.org/data/nsidc-0709/versions/2

Parameters:

version (str, optional) – choose which version to retrieve, by default “depoorter-2013”

Returns:

file path

Return type:

str

References

[30]

Mathieu A. Depoorter, Jonathan L. Bamber, Jennifer Griggs, Jan T. M. Lenaerts, Stefan R. M. Ligtenberg, Michiel R. van den Broeke, and Geir Moholdt. Antarctic masks (ice-shelves, ice-sheet, and islands), link to shape file. 2013. doi:10.1594/PANGAEA.819147.

[31]

M. A. Depoorter, J. L. Bamber, J. A. Griggs, J. T. M. Lenaerts, S. R. M. Ligtenberg, M. R. van den Broeke, and G. Moholdt. Calving fluxes and basal melt rates of Antarctic ice shelves. Nature, 502(7469):89–92, 2013. doi:10.1038/nature12567.

[32] (1,2)

J. Mouginot, B. Scheuchl, and E. Rignot. MEaSUREs Antarctic Boundaries for IPY 2007-2009 from Satellite Radar, Version 2. https://nsidc.org/data/nsidc-0709/versions/2, 2017.

ibcso(layer, region=None, spacing=None, registration=None)

Load IBCSO v2 data, from Dorschel et al.[33] and Dorschel et al.[34].

Parameters:

• layer (str) – choose which layer to fetch: ‘surface’, ‘bed’

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (str or int, optional) – grid spacing to resample the loaded grid to, by default

• registration (str, optional) – change registration with either ‘p’ for pixel or ‘g’ for gridline registration, by default is None.

Returns:

Returns a loaded, and optional clip/resampled grid of IBCSO data.

Return type:

xr.DataArray

References

[33] (1,2)

Boris Dorschel, Laura Hehemann, Sacha Viquerat, Fynn Warnke, Simon Dreutter, Yvonne Schulze Tenberge, Daniela Accettella, Lu An, Felipe Barrios, Evgenia Bazhenova, Jenny Black, Fernando Bohoyo, Craig Davey, Laura De Santis, Carlota Escutia Dotti, Alice C. Fremand, Peter T. Fretwell, Jenny A. Gales, Jinyao Gao, Luca Gasperini, Jamin S. Greenbaum, Jennifer Henderson Jencks, Kelly Hogan, Jong Kuk Hong, Martin Jakobsson, Laura Jensen, Johnathan Kool, Sergei Larin, Robert D. Larter, German Leitchenkov, Benoît Loubrieu, Kevin Mackay, Larry Mayer, Romain Millan, Mathieu Morlighem, Francisco Navidad, Frank O. Nitsche, Yoshifumi Nogi, Cécile Pertuisot, Alexandra L. Post, Hamish D. Pritchard, Autun Purser, Michele Rebesco, Eric Rignot, Jason L. Roberts, Marzia Rovere, Ivan Ryzhov, Chiara Sauli, Thierry Schmitt, Alessandro Silvano, Jodie Smith, Helen Snaith, Alex J. Tate, Kirsty Tinto, Philippe Vandenbossche, Pauline Weatherall, Paul Wintersteller, Chunguo Yang, Tao Zhang, and Jan Erik Arndt. The International Bathymetric Chart of the Southern Ocean Version 2. Scientific Data, 9(1):275, December 2022. doi:10.1038/s41597-022-01366-7.

[34] (1,2)

Boris Dorschel, Laura Hehemann, Sacha Viquerat, Fynn Warnke, Simon Dreutter, Yvonne Schulze Tenberge, Daniela Accettella, Lu An, Felipe Barrios, Evgenia A. Bazhenova, Jenny Black, Fernando Bohoyo, Craig Davey, Laura de Santis, Carlota Escutia Dotti, Alice C. Fremand, Peter T. Fretwell, Jenny A. Gales, Jinyao Gao, Luca Gasperini, Jamin S. Greenbaum, Jennifer Henderson Jencks, Kelly A. Hogan, Jong Kuk Hong, Martin Jakobsson, Laura Jensen, Johnathan Kool, Sergei Larin, Robert D. Larter, German L. Leitchenkov, Benoît Loubrieu, Kevin Mackay, Larry Mayer, Romain Millan, Mathieu Morlighem, Francisco Navidad, Frank-Oliver Nitsche, Yoshifumi Nogi, Cécile Pertuisot, Alexandra L. Post, Hamish D. Pritchard, Autun Purser, Michele Rebesco, Eric Rignot, Jason L. Roberts, Marzia Rovere, Ivan Ryzhov, Chiara Sauli, Thierry Schmitt, Alessandro Silvano, Jodie E. Smith, Helen Snaith, Alex J. Tate, Kirsty Tinto, Philippe Vandenbossche, Pauline Weatherall, Paul Wintersteller, Chunguo Yang, Tao Zhang, and Jan Erik Arndt. The International Bathymetric Chart of the Southern Ocean Version 2 (IBCSO v2). 2022. doi:10.1594/PANGAEA.937574.

ibcso_coverage(region)

Load IBCSO v2 data, from Dorschel et al.[33] and Dorschel et al.[34].

Parameters:

region (tuple[float, float, float, float]) – GMT-format region to subset the data from.

Returns:

Returns a geodataframe of a subset of IBCSO v2 point measurement locations

Return type:

gpd.GeoDataFrame

References

ice_vel(region=None, spacing=5000.0, registration=None, **kwargs)

MEaSUREs Phase-Based Antarctica Ice Velocity Map, version 1 from Mouginot et al.[35] and Mouginot et al.[36].

accessed from https://nsidc.org/data/nsidc-0754/versions/1#anchor-1 Data part of https://doi.org/10.1029/2019GL083826

Units are in m/yr

Parameters:

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (float,) – grid spacing to resample the loaded grid to, by default 5e3, original spacing is 450m

• registration (str, optional) – change registration with either ‘p’ for pixel or ‘g’ for gridline registration, by default is None.

• kwargs (Any) – additional keyword arguments to pass to resample_grid

Returns:

Returns a calculated grid of ice velocity in meters/year.

Return type:

xr.DataArray

References

[35]

J. Mouginot, E. Rignot, and B. Scheuchl. Continent-wide, interferometric SAR phase, mapping of Antarctic ice velocity. Geophysical Research Letters, 46(16):9710–9718, 2019. doi:10.1029/2019GL083826.

[36]

J. Mouginot, E. Rignot, and B. Scheuchl. MEaSUREs Phase Map of Antarctic Ice Velocity, Version 1. 2019. doi:10.5067/PZ3NJ5RXRH10.

imagery()

Load the file path of Antarctic imagery geotiff from LIMA from Bindschadler et al.[37]. accessed from https://lima.usgs.gov/

will replace with below once figured out login issue with pooch MODIS Mosaic of Antarctica: https://doi.org/10.5067/68TBT0CGJSOJ Assessed from https://daacdata.apps.nsidc.org/pub/DATASETS/nsidc0730_MEASURES_MOA2014_v01/geotiff/

Returns:

file path

Return type:

str

References

[37]

Robert Bindschadler, Patricia Vornberger, Andrew Fleming, Adrian Fox, Jerry Mullins, Douglas Binnie, Sara Jean Paulsen, Brian Granneman, and David Gorodetzky. The Landsat Image Mosaic of Antarctica. Remote Sensing of Environment, 112(12):4214–4226, December 2008. doi:10.1016/j.rse.2008.07.006.

magnetics(version, region=None, spacing=None, registration=None, **kwargs)

Load 1 of 3 ‘versions’ of Antarctic magnetic anomaly grid. from and Golynsky et al.[38].

version=’admap1’ ADMAP-2001 magnetic anomaly compilation of Antarctica. Accessed from http://admap.kopri.re.kr/databases.html

version=’admap2’ ADMAP2 magnetic anomaly compilation of Antarctica from Golynsky et al.[39]. Accessed from https://doi.pangaea.de/10.1594/PANGAEA.892723?format=html#download

version=’admap2_gdb’ Geosoft-specific .gdb abridged files Golynsky et al.[40]. Accessed from https://doi.pangaea.de/10.1594/PANGAEA.892722?format=html#download

Parameters:

• version (str) – Either ‘admap1’, ‘admap2’, or ‘admap2_gdb’

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (str or int, optional) – grid spacing to resample the loaded grid to, by default 10e3

• registration (str, optional,) – choose between ‘g’ (gridline) or ‘p’ (pixel) registration types, by default is the original type of the grid

• kwargs (Any) – key word arguments to pass to resample_grid.

Returns:

Returns a loaded, and optional clip/resampled grid of magnetic anomalies.

Return type:

xr.DataArray

References

[38]

Alexander Golynsky, Massimo Chiappini, Detlef Damaske, Fausto Ferraccioli, Carol A. Finn, Takemi Ishihara, Hyung Rae Kim, Luis Kovacs, Valery N. Masolov, Peter Morris, and Ralph von Frese. ADMAP — A Digital Magnetic Anomaly Map of the Antarctic. In Dieter Karl Fütterer, Detlef Damaske, Georg Kleinschmidt, Hubert Miller, and Franz Tessensohn, editors, Antarctica: Contributions to Global Earth Sciences, pages 109–116. Springer, Berlin, Heidelberg, 2006. doi:10.1007/3-540-32934-X_12.

[39]

Alexander V. Golynsky, Fausto Ferraccioli, Jong Kuk Hong, Dmitry A. Golynsky, Ralph R. B. von Frese, Duncan A. Young, Donald D. Blankenship, John W. Holt, Sergey V. Ivanov, Alexander V. Kiselev, Valery N. Masolov, Graeme Eagles, Karsten Gohl, Wilfried Jokat, Detlef Damaske, Carol A. Finn, Alan Aitken, Robin E. Bell, Egidio Armadillo, Tom A. Jordan, Jamin S. Greenbaum, Emanuele Bozzo, Giorgio Caneva, René Forsberg, Marta E. Ghidella, Jesús Galíndo-Zaldívar, Fernando Bohoyo, Yasmina M. Martos, Yoshifumi Nogi, Enrica Quartini, Hyung Rae Kim, and Jason L. Roberts. ADMAP2 Magnetic anomaly map of the Antarctic - links to grid (grd) files. 2018. doi:10.1594/PANGAEA.892723.

[40]

Alexander V. Golynsky, Fausto Ferraccioli, Jong Kuk Hong, Dmitry A. Golynsky, Ralph R. B. von Frese, Duncan A. Young, Donald D. Blankenship, John W. Holt, Sergey V. Ivanov, Alexander V. Kiselev, Valery N. Masolov, Graeme Eagles, Karsten Gohl, Wilfried Jokat, Detlef Damaske, Carol A. Finn, Alan Aitken, Robin E. Bell, Egidio Armadillo, Tom A. Jordan, Jamin S. Greenbaum, Emanuele Bozzo, Giorgio Caneva, René Forsberg, Marta E. Ghidella, Jesús Galíndo-Zaldívar, Fernando Bohoyo, Yasmina M. Martos, Yoshifumi Nogi, Enrica Quartini, Hyung Rae Kim, and Jason L. Roberts. ADMAP2 Magnetic anomaly map of the Antarctic - links to ADMAP2 gdb files. 2018. doi:10.1594/PANGAEA.892721.

mass_change(version='ais_dhdt_floating')

Ice-sheet height and thickness changes from ICESat to ICESat-2. from Smith et al.[41].

Choose a version of the data to download with the format: “ais_VERSION_TYPE” where VERSION is “dhdt” for total thickness change or “dmdt” for corrected for firn-air content. TYPE is “floating” or “grounded”

add “_filt” to retrieve a filtered version of the data.

accessed from https://digital.lib.washington.edu/researchworks/handle/1773/45388

Units are in m/yr

Parameters:

version (str, optional) – choose which version to retrieve, by default “ais_dhdt_floating”

Returns:

Returns a calculated grid of Antarctic ice mass change in meters/year.

Return type:

xr.DataArray

References

[41]

B. Smith, Helen A. Fricker, Alex S. Gardner, Brooke Medley, Johan Nilsson, Fernando S. Paolo, Nicholas Holschuh, Susheel Adusumilli, Kelly Brunt, Bea Csatho, Kaitlin Harbeck, Thorsten Markus, Thomas Neumann, Matthew R. Siegfried, and H. Jay Zwally. Pervasive ice sheet mass loss reflects competing ocean and atmosphere processes. Science, pages eaaz5845, April 2020. doi:10.1126/science.aaz5845.

measures_boundaries(version)

Load various files from the MEaSUREs Antarctic Boundaries for IPY 2007-2009

from Mouginot et al.[32]. accessed at https://nsidc.org/data/nsidc-0709/versions/2

Parameters:

version (str,) – choose which file to retrieve from the following list: “Coastline”, “Basins_Antarctica”, “Basins_IMBIE”, “IceBoundaries”, “IceShelf”, “Mask”

Returns:

file path

Return type:

str

References

modis_moa(version='750m')

Load the MODIS MoA imagery in either “750m” or “125m” resolutions from MODIS Mosaic of Antarctica 2008-2009 (MOA2009) Image Map, Version 2[42] and Scambos et al.[43].

accessed from https://nsidc.org/data/nsidc-0593/versions/2

Parameters:

version (str, optional) – choose between “750m” or “125m” resolutions, by default “750m”

Returns:

filepath for either 750m or 125m MODIS MoA Imagery

Return type:

str

References

[42]

MODIS Mosaic of Antarctica 2008-2009 (MOA2009) Image Map, Version 2. https://nsidc.org/data/nsidc-0593/versions/2, 2021.

[43]

T.A. Scambos, T.M. Haran, M.A. Fahnestock, T.H. Painter, and J. Bohlander. MODIS-based Mosaic of Antarctica (MOA) data sets: Continent-wide surface morphology and snow grain size. Remote Sensing of Environment, 111(2-3):242–257, 2007. doi:10.1016/j.rse.2006.12.020.

moho(version, region=None, spacing=None, registration=None)

Load 1 of x ‘versions’ of Antarctic Moho depth grids.

version=’shen-2018’ Depth to the Moho relative to the surface of solid earth (bottom of ice/ocean) from Shen et al.[10]. Accessed from https://sites.google.com/view/weisen/research-products?authuser=0 Appears to be almost identical to crustal thickness from Shen et al. 2018

version=’an-2015’ This is fetch.crustal_thickness(version=’an-2015)* -1 Documentation is unclear whether the An crust model from An et al.[11] is crustal thickness or moho depths, or whether it makes a big enough difference to matter.

version=’pappa-2019’ from Pappa et al.[44]. Accessed from supplement material

Parameters:

• version (str) – Either ‘shen-2018’, ‘an-2015’, ‘pappa-2019’, will add later: ‘lamb-2020’, ‘baranov’, ‘chaput’, ‘crust1’, ‘szwillus’, ‘llubes’,

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (int, optional) – grid spacing to resample the loaded grid to, by default spacing is read from downloaded files

• registration (str, optional,) – choose between ‘g’ (gridline) or ‘p’ (pixel) registration types, by default is the original type of the grid

Returns:

Returns a loaded, and optional clip/resampled grid of crustal thickness.

Return type:

xr.DataArray

References

[44]

F. Pappa, J. Ebbing, and F. Ferraccioli. Moho Depths of Antarctica: Comparison of Seismic, Gravity, and Isostatic Results. Geochemistry, Geophysics, Geosystems, 20(3):1629–1645, March 2019. doi:10.1029/2018GC008111.

rema(version='1km', region=None, spacing=None, registration=None)

Load the REMA surface elevation data from Howat et al.[6]. The data are in EPSG3031 and reference to the WGS84 ellipsoid. To convert the data to be geoid-referenced, subtract a geoid model, which you can get from fetch.geoid().

Choose between “1km” or “500m” resolutions with parameter version.

accessed from https://www.pgc.umn.edu/data/rema/

Parameters:

• version (str, optional,) – choose which resolution to fetch, either “1km” or “500m”, by default is “1km”

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (str or int, optional) – grid spacing to resample the loaded grid to, by default 10e3

• registration (str, optional,) – choose between ‘g’ (gridline) or ‘p’ (pixel) registration types, by default is the original type of the grid

Returns:

Returns a loaded, and optional clip/resampled grid of the REMA DEM.

Return type:

xr.DataArray

References

resample_grid(grid, initial_spacing=None, initial_region=None, initial_registration=None, spacing=None, region=None, registration=None, **kwargs)

Resample a grid to a new spacing, region, and/or registration. Method of resampling depends on comparison with initial and supplied values for spacing, region, and registration. If initial values not supplied, will try and extract them from the grid.

Parameters:

• grid (str | xr.DataArray) – grid to resample

• initial_spacing (float | None, optional) – spacing of input grid, if known, by default None

• initial_region (tuple[float, float, float, float] | None, optional) – region of input grid, if known, by default None

• initial_registration (str | None, optional) – registration of input grid, if known, by default None

• spacing (float | None, optional) – new spacing for grid, by default None

• region (tuple[float, float, float, float] | None, optional) – new region for grid, by default None

• registration (str | None, optional) – new registration for grid, by default None

• kwargs (dict[str, str]) –

Returns:

grid, either resampled or same as original depending on inputs. If no resampling, and supplied grid is a filepath, returns filepath.

Return type:

str | xr.DataArray

rosetta_gravity(version='gravity')

Load either a shapefile of ROSETTA-ice flightlines, a dataframe of ROSETTA-Ice airborne gravity data over the Ross Ice Shelf, or a dataframe of ROSETTA-Ice density values from the density inversion. from Tinto et al.[45]. Accessed from https://www.usap-dc.org/view/project/p0010035 This is only data from the first 2 of the 3 field seasons. Columns: Line Number: The ROSETTA-Ice survey line number, >1000 are tie lines Latitude (degrees): Latitude decimal degrees WGS84 Longitude (degrees): Longitude decimal degrees WGS84 unixtime (seconds): The number of seconds that have elapsed since midnight (00:00:00 UTC) on January 1st, 1970 Height (meters): Height above WGS84 ellipsoid x (meters): Polar stereographic projected coordinates true to scale at 71° S y (meters): Polar stereographic projected coordinates true to scale at 71° S FAG_levelled (mGal): Levelled free air gravity

Parameters:

version (str, optional) –

Returns:

Returns a dataframe containing the gravity, density, or flightline data

Return type:

pd.DataFrame

References

[45] (1,2)

K. J. Tinto, L. Padman, C. S. Siddoway, S. R. Springer, H. A. Fricker, I. Das, F. Caratori Tontini, D. F. Porter, N. P. Frearson, S. L. Howard, M. R. Siegfried, C. Mosbeux, M. K. Becker, C. Bertinato, A. Boghosian, N. Brady, B. L. Burton, W. Chu, S. I. Cordero, T. Dhakal, L. Dong, C. D. Gustafson, S. Keeshin, C. Locke, A. Lockett, G. O’Brien, J. J. Spergel, S. E. Starke, M. Tankersley, M. G. Wearing, and R. E. Bell. Ross Ice Shelf response to climate driven by the tectonic imprint on seafloor bathymetry. Nature Geoscience, 12(6):441–449, 2019. doi:10.1038/s41561-019-0370-2.

rosetta_magnetics()

Load a dataframe of ROSETTA-Ice airborne magnetics data over the Ross Ice Shelf from Tinto et al.[45].

Accessed from https://www.usap-dc.org/view/project/p0010035 Columns: Line Number: The ROSETTA-Ice survey line number, >1000 are tie lines Latitude (degrees): Latitude decimal degrees WGS84 Longitude (degrees): Longitude decimal degrees WGS84 unixtime (seconds): The number of seconds that have elapsed since midnight (00:00:00 UTC) on January 1st, 1970 H_Ell (meters): Height above WGS84 ellipsoid x (meters): Polar stereographic projected coordinates true to scale at 71° S y (meters): Polar stereographic projected coordinates true to scale at 71° S Mag_anomaly (nT): magnetic anomaly

Returns:

Returns a dataframe containing the data

Return type:

pd.DataFrame

References

sample_shp(name)

Load the file path of sample shapefiles

Parameters:

name (str) – chose which sample shapefile to load, either ‘Disco_deep_transect’ or ‘Roosevelt_Island’

Returns:

file path as a string

Return type:

str

sediment_thickness(version, region=None, spacing=None, registration=None)

Load 1 of 4 ‘versions’ of sediment thickness data.

version=’ANTASed’ From Baranov et al.[46]. Accessed from https://www.itpz-ran.ru/en/activity/current-projects/antased-a-new-sediment-model-for-antarctica/

version=’tankersley-2022’ From Tankersley et al.[2], Tankersley et al.[3].

version=’lindeque-2016’ From Lindeque et al.[4] and .

version=’GlobSed’ From Straume et al.[47]. Accessed from https://ngdc.noaa.gov/mgg/sedthick/

Parameters:

• version (str,) – choose which version of data to fetch.

• region (tuple[float, float, float, float], optional) – GMT-format region to clip the loaded grid to, by default doesn’t clip

• spacing (str or int, optional) – grid spacing to resample the loaded grid to, by default 10e3

• registration (str, optional) – change registration with either ‘p’ for pixel or ‘g’ for gridline registration, by default is None.

Returns:

Returns a loaded, and optional clip/resampled grid of sediment thickness.

Return type:

xr.DataArray

References

[46]

A. Baranov, A. Morelli, and A. Chuvaev. ANTASed – An Updated Sediment Model for Antarctica. Frontiers in Earth Science, 9:722699, August 2021. doi:10.3389/feart.2021.722699.

[47]

E. O. Straume, C. Gaina, S. Medvedev, K. Hochmuth, K. Gohl, J. M. Whittaker, R. Abdul Fattah, J. C. Doornenbal, and J. R. Hopper. GlobSed: Updated Total Sediment Thickness in the World’s Oceans. Geochemistry, Geophysics, Geosystems, 20(4):1756–1772, April 2019. doi:10.1029/2018GC008115.