In this paper, we study strong solutions and decay estimates for an incompressible Vlasov-magnetohydrodynamic (Vlasov-MHD) model arising in magnetized plasmas. This model consists of the Vlasov equation and the incompressible MHD equations, which interact mutually through the Lorentz force. It can be readily verified that the system has two equilibria: \((\bar{f},\bar{u},\bar{B})=(0,0,0)\) and \(( \tilde{f},\tilde{u},\tilde{B})=((2\pi )^{- \frac{3}{2}}e^{- \frac{|v|^2}{2}},0,0)\) . For any prescribed time interval, we establish the existence and uniqueness of strong solutions under appropriately small initial perturbations that depend on the length of the interval. In both the whole-space and periodic settings, we further derive the corresponding decay estimates for the fluid velocity and magnetic field over the interval of existence. The primary difficulty stems from the lack of drag-type or Fokker–Planck dissipation in the kinetic equation, along with the trilinear coupling generated by the Lorentz force. To address this issue, we combine the characteristic method, support propagation for the kinetic component, refined energy estimates, and Fourier analysis.
