In the last decade, convolution operators of matrix functions have received unusual attention due to their diverse applications. The contractivity properties of matrix convolution semigroups are studied and applications to harmonic functions on Lie groups and Riemannian symmetric spaces are discussed.