Исследовано влияние флуктуаций показателя преломления на филаментацию субтераваттного фемтосекундного импульса на километровой трассе в турбулентной атмосфере. Предложен оригинальный метод численного моделирования филаментации импульса на протяженной атмосферной трассе, в котором используются аналитические решения для частотно-модулированного импульса в воздухе. Приведена детальная картина филаментации, происходящая в фазово-модулированном импульсе на расстояниях свыше 1 км.