Supersonic turbulence is believed to decay rapidly within a flow crossing time irrespective of the degree of magnetization. However, this consensus of decaying magnetohydrodynamic (MHD) turbulence relies on local isothermal simulations, which are unable to investigate the role of global magnetic fields and structures. Utilizing three-dimensional MHD simulations including interstellar cooling and heating, we investigate decaying MHD turbulence within cold neutral medium sheets embedded in warm neutral medium. Early evolution is consistent with previous studies characterized rapid decay of turbulence with the decaying time shorter than a flow crossing time and power-law temporal decay of turbulent kinetic energy with slope of -1. If initial magnetic fields are strong and perpendicular to the sheet, however long term evolutions of kinetic energy shows that a significant amount of turbulent energy still remains even after ten flow crossing times, and decaying rate is reduced as field strengths increase. We analyse power spectra of remaining turbulence to show that incompressible, in-plane motions dominate.