Adjoint-based gradient estimation from gray-box solutions of unknown conservation laws