In this paper, we carry out in-depth research centering around the Harnack inequality for positive solutions to nonlinear heat equation on Finsler metric measure manifolds with weighted Ricci curvature \(\textrm{Ric}_{\infty }\) bounded below. Aim on this topic, we first give a volume comparison theorem of Bishop-Gromov type. Then we prove a weighted Poincaré inequality by using Whitney-type coverings technique and give a local uniform Sobolev inequality. Further, we obtain two mean value inequalities for positive subsolutions and supersolutions of a class of parabolic differential equations. From the mean value inequality, we also derive a local gradient estimate for positive solutions to heat equation. Finally, as the application of the mean value inequalities and weighted Poincaré inequality, we get the desired Harnack inequality for positive solutions to heat equation.
