By using the Lusztig's a-function, we give a combinatorial algorithm for Gelfand-Kirillov dimensions of simple highest weight modules of Lie algebras sp 2n , so2n and so2n+1 in terms of their highest weights. Then we determine the associated varieties of highest weight Harish-Chandra modules of Lie groups Sp(2n, R), SO * (2n), SO(2, 2n − 1) and SO(2, 2n − 2) by computing their Gelfand-Kirillov dimensions.