Understanding variations in atmospheric ozone in the Arctic is difficult because there are only a few long-term records of vertical ozone profiles in this region. We present 12 years of ozone profiles from February 2005 to February 2017 at four sites: Summit Station, Greenland; Ny-Ålesund, Svalbard, Norway; and Alert and Eureka, Nunavut, Canada. These profiles are created by combining ozonesonde measurements with ozone profile retrievals using data from the Microwave Limb Sounder (MLS). This combination creates a high-quality dataset with low uncertainty values by relying on in situ measurements of the maximum altitude of the ozonesondes (∼30 km) and satellite retrievals in the upper atmosphere (up to 60 km). For each station, the total column ozone (TCO) and the partial column ozone (PCO) in four atmospheric layers (troposphere to upper stratosphere) are analyzed. Overall, the seasonal cycles are similar at these sites. However, the TCO over Ny-Ålesund starts to decline 2 months later than at the other sites. In summer, the PCO in the upper stratosphere over Summit Station is slightly higher than at the other sites and exhibits a higher standard deviation. The decrease in PCO in the middle and upper stratosphere during fall is also lower over Summit Station. The maximum value of the lower- and middle-stratospheric PCO is reached earlier in the year over Eureka. Trend analysis over the 12-year period shows significant trends in most of the layers over Summit and Ny-Ålesund during summer and fall. To understand deseasonalized ozone variations, we identify the most important dynamical drivers of Arctic ozone at each level. These drivers are chosen based on mutual selected proxies at the four sites using stepwise multiple regression (SMR) analysis of various dynamical parameters with deseasonalized data. The final regression model is able to explain more than 80 % of the TCO and more than 70 % of the PCO in almost all of the layers. The regression model provides the greatest explanatory value in the middle stratosphere. The important proxies of the deseasonalized ozone time series at the four sites are tropopause pressure (TP) and equivalent latitude (EQL) at 370 K in the troposphere, the quasi-biennial oscillation (QBO) in the troposphere and lower stratosphere, the equivalent latitude at 550 K in the middle and upper stratosphere, and the eddy heat flux (EHF) and volume of polar stratospheric clouds throughout the stratosphere.