I've came across the Jefimenko's equations, which are the general solutions of Maxwell's Equations and are compatible with Special Relativity.
They are formulated in terms of the retarded potential: $$ t_r = t - \frac{|\mathbf r - \mathbf{r'}|}{c} \\ \varphi(\mathbf r, t) = \frac{1}{4\pi\varepsilon_0} \int{\frac{\rho(\mathbf r', t_r)}{|\mathbf r - \mathbf{r'}|} \; \mathrm d^3 \mathbf r'} $$ where $t_r$ is the time taken for an observer at $\mathbf{r'}$, to see a change in the potential (or in the field).
It is known that Newton's Theory of Gravity, is not compatible with special relativity and, also, the field propagates throughout space immediately. Would proposing the following be a folly: $$ \phi(\mathbf r, t) = G \int{\frac{\rho(\mathbf r', t_r)}{|\mathbf r - \mathbf{r'}|} \; \mathrm d^3 \mathbf r'} $$ Is there something that doesn't allow you to do this? I know proposing this equation, is pure phenomenology and intuition (because if it describes EM-waves, why not adjust it to "predict" gravity ones?). There it goes my question: would this predict waves of the $\mathbf g$ field? For the Maxwell Equations to predict EM-waves, you need the $\mathbf E$ and $\mathbf B$ fields (two fields), but now, you have only the $\mathbf g$ field.
I did a simulation with only one field that satisfies the potential I proposed and we get the following. It reminds me of grvitational waves in GR. Here is the simulation of a normal $\mathbf g$ field without retarded potential. We might conclude that a wave equation is hiding in here, no?
My question is: Is it possible to describe relativistic gravity without GR using this potential? Why it haven't been done before? Are there any inconsistencies in the reasoning?