Skip to content

Latest commit

 

History

History
120 lines (94 loc) · 3.28 KB

Fig1_Map_CCAMLR-MPA-mgt_MEASO.md

File metadata and controls

120 lines (94 loc) · 3.28 KB
Raw

Map of CCAMLR MPA Planning Domains and MEASO regions

Denisse Fierro Arcos

CCAMLR MPA Planning Domains and MEASO

In this notebook, we will create a map with the limits of the nine marine protected areas (MPA) planning domains from the Commission for the Conservation of Antarctic and Marine Living Resources (CCAMLR). These nine domains are used for planning and reporting on the development of MPAs. The limits of the Marine Ecosystem Assessment for the Southern Ocean (MEASO) are also included in the map. MEASO offers standardised regions that researchers can use to assess and report changes in Southern Ocean ecosystems.

Loading libraries

library(sf)
Linking to GEOS 3.10.2, GDAL 3.4.1, PROJ 8.2.1; sf_use_s2() is TRUE

library(dplyr)
Attaching package: 'dplyr'

The following objects are masked from 'package:stats':

    filter, lag

The following objects are masked from 'package:base':

    intersect, setdiff, setequal, union

library(stringr)
library(rnaturalearth)
library(ggplot2)
library(measoshapes)

Loading relevant shapefiles

#Location of CCAMLR shapefile
ccamlr_path <- "../../SO_shapefiles/CCAMLR_MPA_Planning_Domains/CCAMLR_MPA_Planning_Domains.shp"
#Loading shapefile
ccamlr <- read_sf(ccamlr_path) 

#Loading MEASO regions
measo <- measo_regions05 |> 
  #Remove temperate regions
  filter(str_detect(name, "T$", negate = T)) |> 
  #Matching CCAMLR CRS
  st_transform(st_crs(ccamlr))
old-style crs object detected; please recreate object with a recent sf::st_crs()

#Getting Antarctica boundaries
world <- ne_countries(returnclass = "sf") |> 
  st_transform(crs = st_crs(ccamlr))

Plotting Southern Ocean models

p1 <- ggplot()+
  #Plot CCAMLR
  geom_sf(data = ccamlr, aes(fill = Location))+
  #Choosing colourbling friendly palette
  scale_fill_brewer(type = "qual", palette = "Paired")+
  #Overlaying MEASO boundaries as lines
  geom_sf(data = measo, color = "#999999", fill = NA, linewidth = 0.5, 
          linetype = "dashed")+
  #Changing legend title
  guides(fill = guide_legend(title = "CCAMLR MPA domains"))+
  #Add world
  geom_sf(data = world)+
  #Setting map limits
  lims(x = c(-5774572.727594968, 5774572.727594968), 
       y = c(-5774572.727594968, 5774572.727594968))+
  #Applying B&W theme
  theme_bw()

p1

Saving plot

# Saving plot as image file
ggsave("../outputs/ccamlr_measo_map.pdf", p1, device = "pdf", width = 9,
       height = 7)

#Saving plot as R object for further processing
saveRDS(p1, "../outputs/ccamlr_measo_map.rds")