Abstract. We propose a method to examine the current status of the ozone recovery attributed to changes of ozone-depleting substances (ODS) in the stratosphere. The total column ozone (TCO3) datasets used are based on the ground-based (by the Dobson and/or Brewer spectrophotometer) measurements, satellite observations (from the Solar Backscatter Ultraviolet (SBUV) and Ozone Mapping and Profiler Suite (OMPS) instruments), and output of reanalyses (Multi-Sensor Reanalysis version 2 (MSR2) and Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA2)). The TCO3 time series are calculated for selected sites in the mid-latitudes of the Northern Hemisphere (NH, 35–60∘ N), which are station locations with long-term TCO3 observations archived at the World Ozone and Ultraviolet Radiation Data Centre (WOUDC). The TCO3 monthly means (1980–2020) are averaged over the April–September period to obtain TCO3 time series for the warm sub-period of the year. Two types of the averaged TCO3 time series are considered: the original one and non-proxy time series with removed natural variability by a standard multiple regression model. The TCO3 time series were smoothed by the locally weighted scatterplot smoother (LOWESS) and the super smoother (SS). The smoothed TCO3 values in 1980, 1988, 1997, and 2020 were used to build ozone recovery indices (ORIs) in 2020. These are key years in the equivalent effective stratospheric chlorine (EESC) time series for the period 1980–2020, i.e., the stratosphere was only slightly contaminated by ODS in 1980, 1988 is the year in which the EESC value is equal to its value at the end (2020), and in 1997, the EESC maximum was in mid-latitude stratosphere. The first proposed ORI, ORI1, is the normalized difference between the TCO3 values in 2020 and 1988. The second one, ORI2, is the percentage of the recovered TCO3 in 2020 since the ODS maximum. Following these definitions, the corresponding reference ranges (from −0.5 % to 1 % for ORI1 and from 40 % to 60 % for ORI2) are obtained by analyzing a set of possible EESC time series simulated via the Goddard automailer. The ozone recovery phases are classified comparing the current ORI values and their uncertainty ranges (by the bootstrapping) with these reference ranges. In the analyzed TCO3 time series, for specific combinations of datasets, data types, and the smoother used, we find faster (for ORI1 or ORI2 above the reference range) and slower (for ORI1 or ORI2 below the reference range) recovery in 2020 than that inferred from the EESC change, and a continuation of the TCO3 decline after the EESC peak (ORI2<0 %). Strong signal of the slower TCO3 recovery is found in Toronto, Hohenpeissenberg, Hradec Kralove, and Belsk. A continuation of ozone decline after the turnaround in ODS concentration is found in both the original and non-proxy time series from WOUDC (Toronto), SBUV and OMPS (Toronto, Arosa, Hohenpeissenberg, Uccle, Hradec Kralove, and Belsk), and MERRA2 data (Arosa, Hohenpeissenberg, Hradec Kralove, and Belsk).