PDE penalization for spatial fields smoothing