Atmospheric propagational phase variations are the dominant source of error for InSAR timeseries analysis, generally exceeding uncertainties from poor SNR or signal correlation. The spatial properties of these errors have been well studied, but their temporal dependence and correction have received much less attention to date. We present here an evaluation of the magnitude of tropospheric artifacts in derived time series after compensation using an algorithm that requires only the InSAR data themselves. The level of artifact reduction equals or exceeds that from many weather model based methods, while avoiding the need to access fine-scale atmosphere parameters globally at all times. Our method consists of identifying all points in an InSAR stack with consistently high correlation, and computing, then removing, a fit of the phase at each of these points with respect to elevation. Comparison with GPS truth yields a reduction of 3, from an rms misfit of 5-6 cm to ~2 cm over time. This algorithm can be readily incorporated into InSAR processing flows without need for outside information.