Gaussian process for estimating parameters of partial differential equations and its application to the Richards equation