Here I present a novel framework for the solution of the inverse problem associated with coupled flow and geomechanics models of CO2 injection and subsurface storage, and this solution builds from the application of Gaussian Processes, MCMC and posterior predictive analysis to effectively incorporate relevant earthquake attributes (earthquake time, location and magnitude).