Plotting sf file ontop of ggmap

Question

I've been trying to plot a rivers system over a ggmap and hitting a bunch of walls. Hopefully there is a good solution.

Here is where I'm getting the rivers data: https://data.review.fao.org/map/catalog/srv/api/records/6a53d768-1e20-46ea-92a8-c4040286057d

Loading in:

basemap <- get_stamenmap(bbox = c(left = 149.5, bottom = -35.9, right = 151.5, top = -32.5),
 zoom = 3, maptype = 'terrain-background')

rr <- st_read("rivers_australia_37252/rivers_australia_37252.shp")
 box = c(xmin = 145, ymin = -37, xmax = 155, ymax = -30)
rivers <- st_crop(rr, box)

class(rivers)
[1] "sf"         "data.frame"

Plotting and problem code

ggmap(basemap) +
 geom_sf(data=rivers, inherit.aes = FALSE)
# Error in st_cast.POINT(x[[1]], to, ...) : cannot create MULTILINESTRING from POINT


ggmap(basemap) +
 geom_sf(data=rivers, aes(geometry), inherit.aes = FALSE)
# Error in is.finite(x) : default method not implemented for type 'list'

I then tried unlist() and it came up with a fortify error. Any suggestions of how to transform the data or what to add in the geom_sf() code would be appreciated. Thanks!