Combinations of various numerical models and datasets with diverse observation characteristics have been used to assess the mass evolution of ice sheets. As a consequence, a wide range of estimates have been produced using markedly different methodologies, data, approximation methods and model assumptions. Current attempts to reconcile these estimates using simple combination methods are unsatisfactory, as common sources of errors across different methodologies may not be accurately quantified (e.g. systematic biases in models). Here we provide a general approach which deals with this issue by considering all data sources simultaneously, and, crucially, by reducing the dependence on numerical models. The methodology is based on exploiting the different space-time characteristics of the relevant ice-sheet processes, and using statistical smoothing methods to establish the causes of the observed change. In omitting direct dependence on numerical models, the methodology provides a novel means for assessing glacio-isostatic adjustment and climate models alike, using remote-sensing datasets. This is particularly advantageous in Antarctica, where in situ measurements are difficult to obtain. We illustrate the methodology by using it to infer Antarctica's mass trend from 2003 to 2009 and produce surface mass-balance anomaly estimates to validate the RACMO2.1 regional climate model.