this is not an answer to the question "what matrix model produces the eigenvalue distribution (2)", but it does explain why (2) cannot be an extrapolation between the $\beta$-Ginibre ensembles.
The eigenvalue distribution in the Ginibre ensemble is remarkably complicated, it only has the simple form (2) for $\beta=2$, see equations 1,2,3 of Eigenvalue statistics of the real Ginibre ensemble. For that reasonIf $\beta=4$, instead of a factor $|\lambda_i-\lambda_j|^4$ the repulsion is of the form $|\lambda_i-\lambda_j|^2 |\lambda_i-\bar{\lambda}_j|^2$, and for $\beta=1$ one has three repulsion factors: one factor $|\lambda_i-\lambda_j|^2 |\lambda_i-\bar{\lambda}_j|^2$ between complex eigenvalues $\lambda_i$, one factor $(\mu_i-\mu_j)$ between real eigenvalues $\mu_i$, and one factor $|\lambda_i-\mu_j|^2$ between a real and a complex eigenvalue.
So you see, the "interpolating distribution" (2) is not a natural object for complex eigenvalues --- unlikeeigenvalue statistics in the case of real eigenvalues, where it has the same form (1) for $\beta=1,2,4$complex plane.