Include model structure in pairwise comparisons with multivariate data

Question

I\'m working with a multivariate soil community dataset collected from a randomized complete block design. Block is contributing to a lot of the variation in soil community