Linearized cosimulation of hydraulic conductivity, pressure head, and flux in saturated and unsaturated, heterogeneous porous media