diff --git a/NAMESPACE b/NAMESPACE
index 70471a8..2cd8db1 100644
--- a/NAMESPACE
+++ b/NAMESPACE
@@ -65,6 +65,7 @@ importFrom(dplyr,everything)
 importFrom(dplyr,filter)
 importFrom(dplyr,if_any)
 importFrom(dplyr,mutate)
+importFrom(dplyr,pull)
 importFrom(dplyr,rename)
 importFrom(dplyr,select)
 importFrom(dplyr,tbl)
diff --git a/R/get_subset.R b/R/get_subset.R
index c835e9c..74683af 100644
--- a/R/get_subset.R
+++ b/R/get_subset.R
@@ -19,6 +19,7 @@
 #' @importFrom arrow open_dataset
 #' @importFrom dplyr select filter collect `%>%` everything if_any any_of distinct rename
 #' @importFrom sf st_set_crs write_sf st_sfc st_point st_bbox
+
 get_subset <- function(
   id = NULL, 
   comid = NULL,  
@@ -55,7 +56,7 @@ get_subset <- function(
     .new <- query_set_sink(.new, sink = outfile, overwrite = overwrite)
   }
 
-  query_subset(.new)
+  query_subset(query = .new)
 }
 
 
diff --git a/R/query_subset.R b/R/query_subset.R
index 1888cac..b081425 100644
--- a/R/query_subset.R
+++ b/R/query_subset.R
@@ -14,17 +14,17 @@ query_subset <- function(query) {
     type = class(identifier)
   )
   
-  network <- query_source_layer(query$source, "network")
+  net <- query_source_layer(query$source, "network")
   
   if(!is.null(suppressWarnings(origin$vpuid))){
-   network = dplyr::filter(network, vpuid == !!origin$vpuid)
+   net = dplyr::filter(net, vpuid == !!origin$vpuid)
   } 
   
-  network <- network |>
+  network <- net |>
     dplyr::select(dplyr::any_of(c("id", "toid", "divide_id", "poi_id"))) |>
     dplyr::distinct() |>
     dplyr::collect()
-
+  
   topology <- suppressWarnings(nhdplusTools::get_sorted(network, outlets = origin$toid))
   
   topology$toid[nrow(topology)] <- NA
@@ -36,9 +36,23 @@ query_subset <- function(query) {
 
   all_identifiers <-
     all_identifiers[!is.na(all_identifiers)]
+  
+  if('lakes' %in% query$layers){
+    lake_id <- net |> 
+      dplyr::select(id, hl_uri, poi_id) |> 
+      dplyr::filter(!is.na(hl_uri)) |> 
+      dplyr::filter(id %in% !!unique(all_identifiers)) |> 
+      dplyr::distinct() |>
+      dplyr::collect() |> 
+      dplyr::filter(grepl("LAKE", hl_uri)) |> 
+      dplyr::pull(poi_id)
+  
+  } else {
+    lake_id <- NULL
+  }
 
   query$vpuid <- suppressWarnings({ origin$vpuid })
-  query$requested <- all_identifiers
+  query$requested <- c(all_identifiers, lake_id)
 
   query_extract(query)
 }
diff --git a/hfsubsetR.Rproj b/hfsubsetR.Rproj
index 21a4da0..913c085 100644
--- a/hfsubsetR.Rproj
+++ b/hfsubsetR.Rproj
@@ -1,4 +1,5 @@
 Version: 1.0
+ProjectId: 0a2ccedc-d9ac-44e5-be01-229a60452314
 
 RestoreWorkspace: Default
 SaveWorkspace: Default
diff --git a/man/hfsubsetR-package.Rd b/man/hfsubsetR-package.Rd
index 5cc27fd..97e6cca 100644
--- a/man/hfsubsetR-package.Rd
+++ b/man/hfsubsetR-package.Rd
@@ -6,6 +6,8 @@
 \alias{hfsubsetR-package}
 \title{hfsubsetR: Hydrofabric Subsetter}
 \description{
+\if{html}{\figure{logo.png}{options: style='float: right' alt='logo' width='120'}}
+
 Subset Hydrofabric Data in R.
 }
 \seealso{