obj = AOI::getAOI("Colorado Springs") %>% findNHD() %>% add_retro() names(obj) [1] "comid" "id" "ogc_fid" "fdate" "resolution" "gnis_id" [7] "gnis_name" "lengthkm" "reachcode" "flowdir" "wbareacomi" "ftype" [13] "fcode" "streamleve" "streamorde" "streamcalc" ...