Given the strong recent interest in the decadal timescale variability and the potential for its predictability, it is critical to identify dynamics that carry inherent decadal-scale predictability. This work enhances our understanding and prediction capability of the subsurface signature of the decadal variability in the eastern North Pacific upwelling systems using reanalysis products and a set of eddy-resolving ocean model simulations.