Take a look at my explorations of multivariate normal and an alternative here. I wonder if it could make sense for your application to abandon modelling the correlation matrix as a whole hierarchically, but nonetheless model each individual correlation hierarchically (presumably on the Fisher’s-Z scale).