# load libraries
library(sf)

# load and crop data
data <- list(NULL)
sf_use_s2(FALSE)
for (layer in layers){
  data[[layer]] <- st_crop(st_read(paste0(nhdplusTools_data_dir(),"/data.gpkg"), layer = layer), y = extent_bbox)
}