Physics-based Penalization for Hyperparameter Estimation in Gaussian Process Regression