I think so. It's not obvious to me how.
gsplot(frame.plot=FALSE) %>%
background_color(col="grey90") %>%
grid(col = "white", lty = 1) %>%
par(tcl=-0.2, cex=0.75, xaxs= "r", yaxs = "r") %>%
axis(side = c(1,2), lwd = 0, lwd.ticks = 1, n.minor=1)
@jiwalker-usgs maybe I can pick your brain to see if we can use some lazy_eval to get those locations.