We explore the chaotic dynamics of the mass deformed ABJM model. To do so, we first perform a dimensional reduction of this model from 2 + 1-to 0 + 1dimensions, considering that the fields are spatially uniform. Working in the 't Hooft limit and tracing over ansatz configurations involving fuzzy two spheres, which are described in terms of the GRVV matrices with collective time dependence, we obtain a family of reduced effective Lagrangians and demonstrate that they have chaotic dynamics by computing the associated Lyapunov spectrum. In particular, we analyze in detail, how the largest Lyapunov exponent, λ L , changes as a function of E/N 2 . Depending on the structure of the effective potentials, we find either, where γ N (k, µ) is a constant determined in terms of the Chern-Simons coupling k, the mass µ, and the matrix level N . Using our results, we investigate the temperature dependence of the largest Lyapunov exponents and give upper bounds on the temperature above which λ L values comply with the MSS bound, λ L ≤ 2πT , and below which it will eventually be violated.