Biochemical models of bone remodeling consider the population dynamics of bone cells, which are regulated by numerous chemical signaling factors. Temporal bone cell dynamics are modeled using ordinary differential equations (ODEs), while spatiotemporal dynamics are modeled using partial differential ODEs are the most common technique for mathematical modeling of bone remodeling but cannot explicitly include geometric and transport effects. Spatiotemporal models that incorporate these effects more accurately depict the bone remodeling process. For example, continuous PDEs can model the migration of osteoclasts and osteoblasts to specific locations within the remodeling site. These are important steps in bone remodeling that ODEs cannot resolve. However, PDEs are more computationally expensive than ODEs because they include spatial and temporal effects. ABMs are less widely adopted for spatiotemporal modeling of bone remodeling . Like PDEs, ABMs can model cell movement and how the spatial positioning influences the bone remodeling cycle. However, ABMs are discrete rather than continuous, so their computational intensity depends on the number of agents and the algorithms used to execute their interaction rules. In ABMs, cells are represented as agents that follow rules to move, proliferate, transform, die, and/or secrete signaling factors. The rules governing these cell actions consider the surrounding cell environment and probabilities for introducing stochasticity into the rules.