of release from the spark, τ and Po,ss will change. If Po(t) < Po,ss, an activation τ based on the opening rate of BK channels is used; otherwise, a deactivation τ based on the closing rate is used (see Figs. 7 and 8). The simulation does not contain any explicit spatial component; distance from a spark to a BK cluster is incorporated by the choice of [Ca(t)] to which each cluster is exposed.