1. A method comprising: obtaining an inverse model of a spatial distribution of a physical property of subsurface formations made using time-varying measurements acquired by sensors deployed proximate the subsurface formations, wherein the spatial distribution of the physical property varies in at least in at least two dimensions; obtaining a set of possible models of the spatial distribution based on the time-varying measurements; obtaining a set of model parameters; reducing a number of model parameters in the set by a covariance-free compression transform; mapping upper and lower limits of a physical bound value to an orthogonal space defining a model polytope including a geometric region of feasible models; performing at least one of random and geometric sampling of the model polytope in a reduced-dimensional space to generate an equi-feasible ensemble of models, the reduced-dimensional space consisting of a hypercube approximation; evaluating probable model samples based on misfits from among the equi-feasible ensemble of models, the misfits being determined by forward simulation for time-varying measurements and at least in two-dimensions, and data misfit rejection; computing final uncertainties based on the equi-feasible model ensemble; and displaying the final uncertainties in at least one map.