1) Si $abcd = 0$ es decir, si al menos uno entre $a,b,c,d$ es cero, todos son claramente cero. Por lo tanto, hay una solución: $$a=b=c=d=0$$
2) Si $abcd \neq 0$ , dejemos que $A=\ln(|a|), B=\ln(|b|), C=\ln(|c]), D=\ln(|d|).$
Tomando los valores absolutos y luego los logaritmos de las 4 ecuaciones, obtenemos el siguiente sistema lineal homogéneo $$\begin{bmatrix}1&1&1&-1\\-1&1&1&1\\1&-1&1&1\\1&1&-1&1 \end{bmatrix}\begin{bmatrix}A\\B\\C\\D \end{bmatrix}=\begin{bmatrix}0\\0\\0\\0 \end{bmatrix}$$
El determinante de la matriz $M$ del sistema es $16 \neq 0$ por lo que el núcleo de $M$ se reduce a $$(A,B,C,D)=(0,0,0,0)=(\ln(1),\ln(1),\ln(1),\ln(1))$$ Por lo tanto
$$|a|=1, |b|=1, |c|=1, |d|=1 \ \ \ (1)$$
Puede que no todas las soluciones de (1) sean soluciones del sistema inicial (porque, al tomar valores absolutos, posiblemente hayamos ampliado el conjunto de soluciones).
Por lo tanto, tenemos que comprobar las 16 combinaciones de signos posibles para $a,b,c$ y $d$ . Haciendo esto, quedan 8 soluciones: $$(a,b,c,d)=(-1, -1, -1, -1), (-1, 1, -1, 1), (-1, 1, 1, -1), (-1, -1, 1, 1), (1, 1, -1, -1), (1, -1, -1, 1) , (1, -1, 1, -1), (1, 1, 1, 1).$$