Here is another way to see this. As noticed by Theo in the comments to the OP, the center of $\mathcal{E}$ is the endomorphism category of $\mathcal{E}$ as and $\mathcal{E}$-$\mathcal{E}$-bimodule category. That it is, its objects consist of the functors $\mathcal{E} \to \mathcal{E}$ equipped with coherence data with respect to the left and right $\mathcal{E}$-actions.
When you take the tensor product of tensor categories $\mathcal{C} \boxtimes \mathcal{D}$, then the left and right actions break apart in the usual way.
For example the $\mathcal{C}$ actions "only act on the $\mathcal{C}$ components". What this means is that we have a formula
$$\mathcal{Z}(\mathcal{C} \boxtimes \mathcal D) = End_{C \boxtimes D, C \boxtimes D} (\mathcal{C} \boxtimes \mathcal{D}) $$
$$ \cong Fun_{C,C}(\mathcal{C}, Fun_{D,D}(\mathcal{D}, \mathcal{C} \boxtimes \mathcal{D}))$$
$$\cong Fun_{C,C}(\mathcal{C}, \mathcal{C} \boxtimes Fun_{D,D}(\mathcal{D}, \mathcal{D}))$$
$$\cong Fun_{C,C}(\mathcal{C}, \mathcal{C}) \boxtimes Fun_{D,D}(\mathcal{D}, \mathcal{D})) $$
$$ = \mathcal{Z}(\mathcal{C}) \boxtimes \mathcal{Z}(D)$$
I think the only one of these isomorphisms which, perhaps, isn't straigtforward is the second one, where we identify $Fun_{D,D}(\mathcal{D}, \mathcal{C} \boxtimes \mathcal{D}) \cong \mathcal{C} \boxtimes Fun_{D,D}(\mathcal{D}, \mathcal{D}))$. Note that this doesn't use anything about the tensor category structure of $\mathcal{C}$. So for fusion categories this is totally obvious since as underlying categories then are simply finite sums of $Vect$. It also holds for finite linear categories, though this takes more work. I think it is likely to hold for some larger classes of categories as well, but I am not sure - you might have to be careful about what tensor product you are taking.