Three vectors $v_1,v_2,v_3$ lie in a hyperplane $H:\alpha x+\beta y+\gamma z+\delta t=0$, in this plane we have $Q(v,v):=(\alpha x+\beta y)^2-(\gamma z+\delta t)^2=0,\forall v\in H$. Thus by polarization $Q(v,w)=\frac 14 Q(v+w,v+w)-Q(v-w,v-w)=0$ for all $v,w\in H$ that yields a relation between columns of your matrix: if $v=(a,b,c,d), w=(A,B,C,D)$, then $Q(v,w)=\alpha^2 aA+\beta^2 bB+\alpha\beta(aB+bA)-\gamma^2 cC-\delta^2 dD-\gamma\delta(cC+dD)$.

Three vectors $v_1,v_2,v_3$ lie in a hyperplane $H:\alpha x+\beta y+\gamma z+\delta t=0$, in this plane we have $Q(v,v):=(\alpha x+\beta y)^2-(\gamma z+\delta t)^2=0,\forall v\in H$. Thus by polarization $Q(v,w)=\frac 14 (Q(v+w,v+w)-Q(v-w,v-w))=0$ for all $v,w\in H$ that yields a relation between columns of your matrix: if $v=(a,b,c,d), w=(A,B,C,D)$, then $Q(v,w)=\alpha^2 aA+\beta^2 bB+\alpha\beta(aB+bA)-\gamma^2 cC-\delta^2 dD-\gamma\delta(cD+dC)$.