Methods and apparatus which characterize the ionospheric error across a network of GNSS reference stations are presented. The method relies on dual-frequency phase measurements in a geometry-free linear combination. The data are filtered for ambiguities and the characteristic parameters of the ionosphere. In combination with filter results from other combinations of phase measurements (ionosphere free combination), the physically-based model provides rapid and reliable ambiguity resolution.