Since the discovery of z ∼ 6 quasars two decades ago, studies of their Lyα-transparent proximity zones have largely focused on their utility as a probe of cosmic reionization. But even when in a highly ionized intergalactic medium, these zones provide a rich laboratory for determining the timescales that govern quasar activity and the concomitant growth of their supermassive black holes. In this work, we use a suite of 1D radiative transfer simulations of quasar proximity zones to explore their timedependent behaviour for activity timescales from ∼ 10 3 -10 8 years. The sizes of the simulated proximity zones, as quantified by the distance at which the smoothed Lyα transmission drops below 10% (denoted R p ), are in excellent agreement with observations, with the exception of a handful of particularly small zones that have been attributed to extremely short 10 4 lifetimes. We develop a physically motivated semi-analytic model of proximity zones which captures the bulk of their equilibrium and non-equilibrium behaviour, and use this model to investigate how quasar variability on < ∼ 10 5 year timescales is imprinted on the distribution of observed proximity zone sizes. We show that large variations in the ionizing luminosity of quasars on timescales of < ∼ 10 4 years are disfavored based on the good agreement between the observed distribution of R p and our model prediction based on "lightbulb" (i.e. steady constant emission) light curves.