It is shown how to non-perturbatively define a random matrix model that captures key physics of ${\cal N}{=}2$ Jackiw-Teitelboim (JT) supergravity, going well beyond the perturbative topological expansion defined recently by Turiaci and Witten. A decomposition into an infinite family of certain multicritical models is derived, leading to the definition of a non-linear differential equation from which the physics may be computed. BPS states are naturally described by the model. The non-perturb...