Hessian-Based Response Surface Approximations For Uncertainty Quantification In Large-Scale Statistical Inverse Problems, With Applications To Groundwater Flow