Polynomial reconstruction of staggered unstructured vector fields has been considered. Coefficients of such polynomials are determined by the least squares method. Reduction in the rank of the least squares systems caused by the over-specification of the divergence may lead to difficulties. This has been investigated. The rank of these systems may be further reduced depending on the mesh geometry, or they may become ill conditioned. Guidelines for solving such linear systems have been presented.