This repository provides a set of essential variables to support research on forest loss driven by mining. All variables have been resampled to 30 arcsec spatial resolution (approximately 1 by 1 km at the equator) and are encoded in Geographic Tagged Image File Format (GeoTIFF). The grid extends from the longitude β180 to 180 degrees and from the latitude β90 to 90 degrees in the geographical reference system WGS84. Cells over water have no-data values. Below we describe the list of variables, sources, and processing
steps.area_of_mines_circa_2018.tif: mining area in square metres. This layer was derived from a global-scale data set of mining polygons [Maus et al., 202a,b0] available from [doi:10.1594/PANGAEA.910894] under CC BY-SA 4.0 license. The mining area for each 30 arcsec grid was calculated intersecting cells and mining
polygons.distance_to_mine_circa_2018.tif: distance to the nearest mine in metres. This layer was derived by calculating the Euclidean distance between each grid cell's centroid to the centroid of the closest grid cell with mine presence, i.e. cells where
area_of_mines_circa_2018.tif > 0.area_of_forest_cover_circa_2000.tif: area of forest cover in square metres. This layer was derived from the Global Forest Change (GFC) dataset [Hansen et al., 2013] version 1.7 available from [https://earthenginepartners.appspot.com/science-2013-global-forest/download_v1.7.html] under CC BY 4.0 license. We aggregated the GFC data from 1 arcsec to our 30 arcsec grid cells by summing the area of forest cover pixels weighted by their surface intersection with the 30 arcsec
cells.area_of_forest_cover_within_mines_circa_2000.tif: area of forest cover in square metres. This layer was derived using the same methods as
area_of_forest_cover_circa_2000.tif; however, it only includes forest area intersecting mining polygons, i.e. the on-site forest cover circa
2000.area_of_forest_cover_loss_yearly_from_2001_to_2019.tif: area of forest cover loss in square metres. This GeoTIFF file has 19 layers (one layer per year) starting from 2000. We aggregated the GFC data from 1 arcsec to our 30 arcsec grid cells by summing the area of forest loss pixels weighted by their surface intersection with the 30 arcsec
cells.ecoregions2017_code.tif: an integer with the ecoregions code (ECO_ID) rasterized from the Ecoregion 2017 polygons [Dinerstein et al., 2017; Resolve, 2017], which is available from [https://ecoregions2017.appspot.com/] under CC BY 4.0 license. The polygons were rasterized to a 30 arcsec grid by the major class present. The ecoregion class names corresponding to the GeoTIFF file values are available in the auxiliary file
ecoregions_2017_concordance_tbl.csv, which contains the following variables ECO_ID, ECO_NAME, BIOME_NUM, BIOME_NAME, where ECO_ID is a unique identifier.
The layers available from this repo can be stacked together with other variables essential for land-use modelling. Some of these variables are openly available at the same spatial extent and resolution, for example, grided population [NASA, 2018], elevation and slope [Amatulli et al., 2018a,b]