Коэффициенты, рассчитанные с помощью fit.variogram, не соответствуют графику с подбором модели (пакет gstat R)

Я подгоняю экспоненциальную модель к эмпирической вариограмме, используя пакет gstat. Мой набор данных (данные) имеет значения широты и долготы в десятичных градусах. Диапазон, оцененный с использованием fit.variogram(), не соответствует диапазону, который я ожидал бы, глядя на график модели. См. код и рисунки ниже.

Это вывод из fit.variogram():

   model    psill    range
1   Nug 26.66719  0.00000
2   Exp 26.26643 25.76976

Отсюда я ожидаю, что диапазон, соответствующий порогу ~ 53 (26,66719 + 26,26643), будет равен 25,76976. Однако, если посмотреть на график, диапазон, соответствующий порогу ~53, составляет >60 км.

library(gstat)  
library(sp)

coordinates(data)= ~ longitude+latitude
proj4string(data)=CRS("+longlat")
TheVariogram=variogram(Svmean~1, data=data, cutoff=90)
plot(TheVariogram) #to get seed parameters
lzn.fit <- fit.variogram(TheVariogram, model = vgm(psill=60, model="Exp", nugget=20, range=70),fit.method=6)
plot(TheVariogram, pl=F, model=lzn.fit, xlab="Distance (km)", ylab=expression(Semivariance (dB^2)))