Bayesian Nonlinear PDE Inference via Gaussian Process Collocation with Application to the Richards Equation