Basemap options within common R spatial packages are not entirely satisfactory when mapping along estuarine coastlines. USGS provides a Hydrography layer that is very useful, but is not built into said packages. This function adds the USGS Hydrography layer to a tmap object that has already been created. In view mode, the output is turned into an interactive leaflet object. In plot mode, it remains a tmap object. This function works with vector objects but is currently untested on rasters.