  We study the complete time evolution of scalar fields propagating in
space-times of higher dimensional Lifshitz Black Holes with dynamical critical
exponent $z=2$, obtained from a theory including the most general quadratic
curvature corrections to Einstein-Hilbert gravity in $D$ dimensions. We also
computed the quasinormal spectrum after performing a numerical integration and
solving exactly the Klein-Gordon equation obeyed by the massive scalar field.
We found that quasinormal modes are purely imaginary for all dimensions.
