No reproducible example, so only a verbose description:
1) extract the year
2a) use ddply(data, .(year), numcolwise(mean))
or alternatively:
2b) use reshape2:::melt in combination with
ddply(melteddata, .(year, variable), summarize, mean=mean(value))