Model-based variance estimation in two-dimensional systematic sampling