Ensemble average and ensemble variance behavior of unsteady, one-dimensional groundwater flow in unconfined, heterogeneous aquifers: an exact second-order model