diff --git a/R/02-fogo-island-figure.R b/R/02-fogo-island-figure.R
index 367492c..dfee56c 100644
--- a/R/02-fogo-island-figure.R
+++ b/R/02-fogo-island-figure.R
@@ -12,8 +12,8 @@ lapply(libs, require, character.only = TRUE)
 
 
 ### Data ----
-fogo <- readRDS('output/fogo-island-polygons.Rds')
-roads <- readRDS('output/fogo-roads.Rds')
+fogo <- st_read('output/fogo-island-polygons.gpkg')
+roads <- st_read('output/fogo-roads.gpkg')
 
 # CRS
 utm <- st_crs('+proj=utm +zone=21 ellps=WGS84')
diff --git a/R/04-newfoundland-figure.R b/R/04-newfoundland-figure.R
index aaec6b4..03d3b8b 100644
--- a/R/04-newfoundland-figure.R
+++ b/R/04-newfoundland-figure.R
@@ -12,7 +12,7 @@ lapply(libs, require, character.only = TRUE)
 
 
 ### Data ----
-nl <- readRDS('output/newfoundland-polygons.Rds')
+nl <- st_read('output/newfoundland-polygons.gpkg')
 
 # CRS
 utm <- st_crs('+proj=utm +zone=21 ellps=WGS84')
diff --git a/R/05-fogo-inset-newfoundland-figure.R b/R/05-fogo-inset-newfoundland-figure.R
index 4b97b43..e695c9b 100644
--- a/R/05-fogo-inset-newfoundland-figure.R
+++ b/R/05-fogo-inset-newfoundland-figure.R
@@ -13,10 +13,10 @@ lapply(libs, require, character.only = TRUE)
 
 
 ### Data ----
-fogo <- readRDS('output/fogo-island-polygons.Rds')
-roads <- readRDS('output/fogo-roads.Rds')
+fogo <- st_read('output/fogo-island-polygons.gpkg')
+roads <- st_read('output/fogo-roads.gpkg')
 
-nl <- readRDS('output/newfoundland-polygons.Rds')
+nl <- st_read('output/newfoundland-polygons.gpkg')
 
 
 # CRS