Intracellular calcium (Ca 2þ) cycling dynamics in cardiac myocytes are spatiotemporally generated by stochastic events arising from a spatially distributed network of coupled Ca 2þ release units that interact with an intertwined mitochondrial network. In this study, we developed a spatiotemporal ventricular myocyte model that integrates mitochondria-related Ca 2þ cycling components into our previously developed ventricular myocyte model consisting of a three-dimensional Ca 2þ release unit network. Mathematical formulations of mitochondrial membrane potential, mitochondrial Ca 2þ cycling, mitochondrial permeability transition pore stochastic opening and closing, intracellular reactive oxygen species signaling, and oxidized Ca 2þ /calmodulin-dependent protein kinase II signaling were incorporated into the model. We then used the model to simulate the effects of mitochondrial depolarization on mitochondrial Ca 2þ cycling, Ca 2þ spark frequency, and Ca 2þ amplitude, which agree well with experimental data. We also simulated the effects of the strength of mitochondrial Ca 2þ uniporters and their spatial localization on intracellular Ca 2þ cycling properties, which substantially affected diastolic and systolic Ca 2þ levels in the mitochondria but exhibited only a small effect on sarcoplasmic reticulum and cytosolic Ca 2þ levels under normal conditions. We show that mitochondrial depolarization can cause Ca 2þ waves and Ca 2þ alternans, which agrees with previous experimental observations. We propose that this new, to our knowledge, spatiotemporal ventricular myocyte model, incorporating properties of mitochondrial Ca 2þ cycling and reactive-oxygen-species-dependent signaling, will be useful for investigating the effects of mitochondria on intracellular Ca 2þ cycling and action potential dynamics in ventricular myocytes.