diff --git a/mapping_draft.Rmd b/mapping_draft.Rmd
index 4b8d0d2..4d4ecba 100644
--- a/mapping_draft.Rmd
+++ b/mapping_draft.Rmd
@@ -887,6 +887,15 @@ unzip("data/National_Forest_Inventory_Woodland_Scotland_2017.zip", exdir = "data
 forest_inventory <- st_read("data/National_Forest_Inventory_Woodland_Scotland_2017.shp")
 forest_inventory_sp <- readOGR("./data", "National_Forest_Inventory_Woodland_Scotland_2017") 
 
+# Download data for scenic areas
+# https://opendata.arcgis.com/datasets/d7f6b987c7224a72a185ce012258d500_23.zip
+# International Union for Conservation of Nature (IUCN), Category V Protected Landscapes
+
+download.file("https://opendata.arcgis.com/datasets/d7f6b987c7224a72a185ce012258d500_23.zip", destfile = "data/ScenicAreas.zip")
+unzip("data/ScenicAreas.zip", exdir = "data")
+scenicareas <- st_read("data/ScenicAreas.shp")
+scenicareas_sp <- readOGR("./data", "ScenicAreas") 
+
 # Set symmetrical CRS for analysis below
 st_crs(sssi) <- 27700
 st_crs(ecs_sf) <- 27700