alaska <- c(63.5888, -154.4931)
dat <- data.frame(lng = alaska[2], lat = alaska[1])
states<- sf::read_sf("~/Downloads/cb_2018_us_state_500k/cb_2018_us_state_500k.shp")%>%
  sf::st_transform('+proj=longlat +datum=WGS84')
Alaska<- states %>%
  dplyr::filter(NAME%in%c("Alaska"))

rivers<-sf::read_sf("~/Downloads/SOA_rivers_levees_AOIs 2/SOA_rivers_levees_AOIs.shp") %>%
  sf::st_transform('+proj=longlat +datum=WGS84')

PA<-sf::read_sf("~/Downloads/Alaska_DEC_Sensitive_Areas/Alaska_DEC_Sensitive_Areas.shp") %>%
  sf::st_transform('+proj=longlat +datum=WGS84')

        
leaflet() %>%
  addProviderTiles(providers$Esri.WorldTopoMap) %>%
  setView(lng = alaska[2], lat = alaska[1], zoom = 3.30) %>%
  addMarkers(data = dat, popup = "Alaska") %>%
  addPolygons(data=Alaska, fill=F, color="green", opacity=1) %>%
  addPolylines(data=rivers,fill=F, color="blue", opacity=1) %>%
  addPolygons(data=PA, fill=F, color="#e38c62", opacity=1) %>%
  addMarkers(lng=-154.693756,lat=58.597813,popup= "Katmai National Park and Preserve, Surface Area = 4,093,077 acres") %>%
  addMarkers(lng=-136.900208,lat=58.665806,popup= "Glacier Bay National Park and Preserve, Surface Area = 3,283,000 acres") %>%
  addMarkers(lng=-142.79843,lat=65.08979,popup= "Yukon-Charley Rivers National Preserve, Surface Area = 2,525,512 acres")%>%
  addMarkers(lng=-134.4895,lat=57.6853,popup= "Admiralty Island National Monument, Surface Area = 955,747 acres")%>%
  addMarkers(lng=-154.323502,lat=60.412697,,popup= "Lake Clark National Park and Preserve, Surface Area = 4,030,015 acres")%>%
  addMarkers(lng=-165.408418,lat=64.497922,popup= "Bering Land Bridge National Preservet, Surface Area = 2,697,391 acres")%>%
  addMarkers(lng=-142.985687,lat=61.710445,popup= "Wrangell-St. Elias National Park & Preserve, Surface Area = 13,175,799 acres")%>%
  addMarkers(lng=-130.686710,lat=55.668243,popup= "Misty Fjords National Monument, Surface Area = 2,294,343 acres")%>%
  addMarkers(lng=-151.197418,lat=63.129887,popup= "Denali National Park & Preserve, Surface Area = 6,045,153 acres")%>%
  addMarkers(lng=-153.3241,lat=67.6855,popup= "Gates Of The Arctic National Park & Preserve, Surface Area = 8,472,506 acres")%>%
  addMarkers(lng=-159.213175,lat=67.406386,popup= "Kobuk Valley National Park, Surface Area = 1,750,716 acres")%>%
  addMarkers(lng=-157.814709,lat=56.870864,popup= "Aniakchak National Monument & Preserve, Surface Area = 601,294 acres")%>%
  addMarkers(lng=-149.816360,lat=60.043777,popup= "Kenai Fjords National Park, Surface Area = 669,984 acres")%>%
  addMarkers(lng=-159.9070,lat=68.1353,popup= "Noatak National Preserve, Surface Area = 6,569,904 acres")%>%
  addMarkers(lng=-163.5054,lat=67.4181,popup= "Cape Krusenstern National Monument, Surface Area = 649,082 acres")

Resources:

Protected Areas: https://gis.data.alaska.gov/datasets/agio-hub::alaska-national-parks-preserves-monuments/about

Alaska outline: https://www.census.gov/geographies/mapping-files/time-series/geo/carto-boundary-file.html

Filter shapefile to only Alaska: https://stackoverflow.com/questions/76028013/outline-multiple-polygons-r-leaflet-to-show-territories

Popups and Markers Help: https://rstudio.github.io/leaflet/articles/popups.html#popups

Solved Error when trying to install package (sf): https://github.com/r-spatial/sf/issues/1839

Solved Error on sf failure to project: https://stackoverflow.com/questions/57223853/how-to-apply-spatial-polygons-to-leaflet-map-using-shp

Polygons and Polylines for leaflet: https://rstudio.github.io/leaflet/articles/shapes.html#rectangles

Shapefile for river: https://www.arcgis.com/home/item.html?id=46068a898cdb4d48b051e27a48c580de

Other information about Protected Areas surface area and coordinates were googled.