We analyze the variability of mean age of air (AoA) and of the local effects of the stratospheric residual circulation and eddy mixing on AoA within the framework of the isentropic zonal mean continuity equation. AoA for the period 1988–2013 has been simulated with the Lagrangian chemistry transport model CLaMS driven by ERA‐Interim winds and diabatic heating rates. Model simulated AoA in the lower stratosphere shows good agreement with both in situ observations and satellite observations from Michelson Interferometer for Passive Atmospheric Sounding, even regarding interannual variability and changes during the last decade. The interannual variability throughout the lower stratosphere is largely affected by the quasi‐biennial‐oscillation‐induced circulation and mixing anomalies, with year‐to‐year AoA changes of about 0.5 years. The decadal 2002–2012 change shows decreasing AoA in the lowest stratosphere, below about 450 K. Above, AoA increases in the Northern Hemisphere and decreases in the Southern Hemisphere. Mixing appears to be crucial for understanding AoA variability, with local AoA changes resulting from a close balance between residual circulation and mixing effects. Locally, mixing increases AoA at low latitudes (40°S–40°N) and decreases AoA at higher latitudes. Strongest mixing occurs below about 500 K, consistent with the separation between shallow and deep circulation branches. The effect of mixing integrated along the air parcel path, however, significantly increases AoA globally, except in the polar lower stratosphere. Changes of local effects of residual circulation and mixing during the last decade are supportive of a strengthening shallow circulation branch in the lowest stratosphere and a southward shifting circulation pattern above.