Multivariate dependent variable with random factor

I've a large data set of ponds of different ages and different clusters/regions.

I want to test the effect of fixed factor age and random factor cluster on the environmental factors. The environmental factors are multivariate and consist of 14 different numerical variables.

I'm struggling to find a good method to execute this model in R