В работе предложен оригинальный метод построения глобальных ионосферных карт полного электронного содержания (TEC), основанный на фазоразностном подходе к анализу измерений фаз навигационных радиосигналов GNSS на паре когерентных рабочих частот на распределенной сети наземных приёмников мировой сети IGS. В предложенном подходе используется представление ионосферы в виде тонкого слоя с распределением TEC, заданным усечённым разложением в ряд по сферическим гармоникам в геомагнитной системе координат локального времени. Коэффициенты разложения при этом определяются МНК с ограничением на положительность TEC, что реализовано за счёт решения задачи линейной дополнительности. Предложенный метод не требует оценки дифференциальных кодовых задержек аппаратуры спутников и приёмников, что позволяет комбинировать в рамках одного алгоритма данные различных GNSS, например GPS, ГЛОНАСС и Galileo. Представлены результаты тестирования предложенного метода на данных синтезированных наблюдений с использованием реальной геометрии спутников GNSS, приёмников IGS и ионосферы, заданной на основе модели NeQuick2. Приводятся результаты сопоставления реальных глобальных ионосферных карт, полученных предложенным методом, и карт центра CODE.