diff --git a/workspace/NHD_navigate.Rmd b/workspace/NHD_navigate.Rmd
index fe2a879cdb9dd2fab60a4fc50ff4a70151c0065a..e905aa86d5903e02d51faf069046114b0ee8d1c9 100644
--- a/workspace/NHD_navigate.Rmd
+++ b/workspace/NHD_navigate.Rmd
@@ -352,12 +352,14 @@ if(needs_layer(out_gpkg, nsegment_raw)) {
   # Write out nsegments composed of nhdflowlines
   write_sf(ncombined, out_gpkg, nsegment_raw)
 
-  # Aggregate flowlines per POI_ID to a single segment, carrying over useful information
-  nsegments<-ncombined %>% group_by(POI_ID) %>%  mutate(VA_MA = na_if(VA_MA, -9998)) %>%
+  # Dissolve flowlines to aggregated segments
+  nsegments <- ncombined %>% group_by(POI_ID) %>% arrange(desc(LevelPathI), desc(Hydroseq)) %>%
     mutate(TT_Hours = (LENGTHKM * 3280.84) / (VA_MA * 3600)) %>% 
-    summarize(TotalLength = sum(LENGTHKM), TotalDA = max(TotDASqKM), HW = max(StartFlag), TT = sum(TT_Hours)) %>%
+    summarize(TotalLength = sum(LENGTHKM),TotalDA = max(TotDASqKM), HW = max(StartFlag), TT = sum(PathTimeMA),
+            do_union=FALSE) %>% st_cast("MULTILINESTRING")  %>%
     inner_join(st_drop_geometry(hucgagessegs) %>% select(COMID, Hydroseq, DnHydroseq),   by = c("POI_ID" = "COMID"))
-
+  
+  
   # produce a short data frame for populating TO_POI for downstream segment
   to_from<-segPOIs %>% inner_join(st_drop_geometry(ncombined), by = c("DnHydroseq" = "Hydroseq")) %>%
     select(COMID.x, Hydroseq, DnHydroseq, POI_ID)