I lack time to give you the details, but the generalThe formula for the specific case is $$\det K=\det(A+B+C+D)\det(A-B+C-D)\det(A+iB-C-iD)\det(A-iB-C+iD).$$ More generally, for a block-circulant matrix with $n$ square blocks $A_0,\ldots,A_{n-1}$, the formula is $$\det K=\prod_{\omega^n=1}\det(A_0+\omega A_1+\ldots+\omega^{n-1}A_{n-1}).$$$$\det K=\prod_{\omega^n=1}\det(A_0+\omega A_1+\cdots+\omega^{n-1}A_{n-1}).$$ At least, youTo see easily by linear combinationsthis, observe that if some of the factors above vanishes$K$ is block-diagonalisable, then $$K=U^*{\rm diag}(A_0+\alpha A_1+\cdots+\alpha^{n-1}A_{n-1},A_0+\alpha^2 A_1+\cdots+\alpha^{2(n-1)}A_{n-1},\ldots)U$$ where $\det K$ vanishes too$\alpha=\exp\frac{2i\pi}n$ and $$U=\frac1{\sqrt n}((\alpha^{(i-1)(j-1)}I_d))_{1\le i,j\le n}.$$ Hereabove, the blocks $A_j$ are $d\times d$.
Perhaps I'll write details tomorrow morning This shows the formula, up to the factor $|\det U|^2$, which is easily seen to be equal to $1$.