Один из вариантов - использовать stat_summary
.
. Это добавит еще один эстет c к графику и, таким образом, еще одну легенду, далеко от другой.
airquality %>%
mutate(Month=factor(Month)) %>%
ggplot(aes(x=Day, y=Temp, col=Month)) +
geom_line() +
stat_summary(aes(lwd="Nederland totaal"), fun=mean, geom="line", col="black") +
theme_classic() +
theme(legend.title = element_blank()) +
guides(lwd = guide_legend(order = 1))

Но преимущество этого метода заключается в том, что вам не нужно рассчитывать средние значения за год и вручную добавлять его в переменную провинции.