In Nielsen and Chuang (page:379), it shows how to represent a 1 qubit depolarizing channel in operator-sum representation. $$ \mathcal{E}_1(\rho)=pI/2+(1-p)\rho =(1-3p/4)\rho+p/4(X\rho X+Y\rho Y+Z\rho Z) $$
How to write a 2 qubit depolarizing channel in operator-sum representation?
$$ \mathcal{E}_2(\rho)=pI/4 +(1-p)\rho $$