m2 <- lmer(flow~topo+period+topo*period+(1|patch)+(1|period)+(1|tod))

flow.lsm <- lsmeans(m2, "topo", by="period")
con      <- contrast(flow.lsm, method="eff", interaction=TRUE)