By utilizing spectral analysis methods and properties of fractional order operators, some new estimates have been established for the nonlinear terms appearing from the two-dimensional non-stationary magnetohydrodynamic flows in a general domain. Based on these, a series of energy decay properties are given. Furthermore, the large time decay behavior of the solution and that of its gradient have also been addressed in the L∞ and L2 spaces, respectively.