The time series of ozone columns measured with the SBUV satellite instrument over three subarctic stations (Saint Petersburg, Harestua, and Kiruna) are analyzed. The daily and monthly mean ozone values in the layers of 0-25, 25-60, and 0-60 km are compared with the results of simulations with RSHU and EMAC numerical models for the period of 2000-2015. Model data are in good agreement with satellite data both in general and in the cases of rapid short-term ozone loss. However, there are some differences between the models and measurements as well as between the two considered models. These differences require the more detailed analysis in order to modify model parameters. Experimental data demonstrate the increase in ozone columns in the layer of 25-60 km which amounts to 2.1 +/- 0.7, 2.4 +/- 0.7, and 1.5 +/- 0.8% per decade for Saint Petersburg, Harestua, and Kiruna stations, respectively. The results of numerical simulations do not contradict these estimates.