Los resultados de una simulación de hindcast que permite submesoescala desde 1990 hasta 2016 se utilizan para investigar las variaciones interanuales a decenales de los movimientos submesoescala. La región en la que nos enfocamos es el Pacífico Noroeste subtropical, incluyendo la contracorriente subtropical. La energía cinética (EC) submesoescala se caracteriza por una fuerte variabilidad interanual y decenal, mostrando magnitudes más grandes en 1996, 2003 y 2015, y magnitudes más pequeñas en 1999, 2009, 2010 y 2016. Estas variaciones se explican parcialmente por las de la liberación de energía potencial disponible (EPD) a submesoescala impulsada por la inestabilidad de la capa mixta en invierno. De hecho, esta liberación de EPD depende de la profundidad de la capa mixta y del gradiente de flotabilidad horizontal, ambos modulados por la Oscilación Decadal del Pacífico (ODP). Como resultado de la cascada inversa de EC, la variabilidad de la EC submesoescala posiblemente conduce a variaciones interanuales a decenales de la EC a mesoescala (EC de remolino (ECR)). Estos resultados muestran que los movimientos submesoescala son una posible vía para explicar el impacto asociado con la ODP en la variabilidad decenal de la ECR. La liberación de EPD en invierno estimada a partir de las observaciones de flotadores Argo varía de manera sincrónica con la de la simulación en las escalas de tiempo interanuales, lo que sugiere la capacidad de observación para diagnosticar la variabilidad de la EC submesoescala.