Un diagnóstico de atribución y causalidad en dos pasos está diseñado empleando análisis de espectro singular para descomponer la serie temporal climática atribuida en una matriz de trayectoria y luego se somete a un análisis de funciones ortogonales empíricas para identificar las fuerzas impulsoras en evolución, que finalmente pueden relacionarse con los principales modos climáticos a través de sus frecuencias independientes mediante análisis de wavelet. Los resultados de la aplicación de la región árida y propensa a sequías del sur de la región Intermontañosa de América del Norte se comparan con la forzamiento climático o de mayor escala diagnosticado a partir del análisis de características lentas utilizando las fuentes del balance de flujo de agua y energía. Se observan los siguientes resultados: (i) Los cambios entre los cuatro períodos de 20 años subsiguientes de 1930 a 2010 sugieren predominantemente un forzamiento inducido por el clima debido a la Oscilación Decadal del Pacífico y la Oscilación Multidecadales del Atlántico. (ii) Las influencias de la cobertura terrestre sobre el cambio de la cobertura terrestre son de magnitud considerablemente menor (en términos de porcentaje de área cubierta) cuya evolución temporal está bien documentada a partir de documentos de forestación. (iii) Los impulsores de los forzamientos inducidos por el clima en los últimos 20 años se identifican como la oscilación cuasi-bienal y la Oscilación del El Niño-Oscilación del Sur tanto por la atribución en dos pasos interanual como por los diagnósticos de causalidad con análisis de características lentas basados en escala mensual.