Chlorophyll a fluorescence induction (FI) kinetics, in the microseconds to the second range, reflects the overall performance of the photosynthetic apparatus. In this paper, we have developed a novel FI model, using a rule-based kinetic Monte Carlo method, which incorporates not only structural and kinetic information on PSII, but also a simplified photosystem I. This model has allowed us to successfully simulate the FI under normal or different treatment conditions, i.e., with different levels of measuring light, under 3-(3',4'-dichlorophenyl)-1,1-dimethylurea treatment, under 2,5-dibromo-3-methyl-6-isopropyl-p-benzoquinone treatment, and under methyl viologen treatment. Further, using this model, we have systematically studied the mechanistic basis and factors influencing the FI kinetics. The results of our simulations suggest that (1) the J step is caused by the two-electron gate at the Q B site; (2) the I step is caused by the rate limitation of the plastoquinol re-oxidation in the plastoquinone pool. This new model provides a framework for exploring impacts of modifying not only kinetic but also structural parameters on the FI kinetics.