AbstractWe study the time fractional diffusion equation , , in a bounded domain with an elliptic operator and a locally Lipschitz nonlinearity on fractional Sobolev spaces, subjected to the homogeneous Dirichlet boundary condition. Data have not been measured at the initial time , but at a final time , that is, is given instead of . The problem is, therefore, called an inverse initial‐value problem. We first establish the well‐posedness of this problem on fractional Sobolev spaces and the regularity of the solution by assuming only the local Lipschitz continuity of . Second, an susceptible‐infected (shortly, SI) model with heterogeneity and a Navier–Stokes equation have been exemplified. Finally, a spatial ‐estimate for the solution and its gradient has been provided. The essential tools are asymptotic behaviours of Mittag–Leffler functions, fractional power spaces, fractional Sobolev spaces and embedding, weighted functional spaces, and estimates for heat semigroup.