Los estudios de diversidad acuática plantean una declinación del número de taxones en relación al incremento altitudinal, debido a la baja presión de oxígeno o a cierta homogeneidad del hábitat, entre otros. Las investigaciones se realizaron generalmente a una escala local o subestimando la gran heterogeneidad ambiental que caracteriza a estas zonas a pesar de su aspecto desértico e inhóspito. El objetivo del trabajo fue determinar el patrón espacial de la diversidad de macroinvertebrados acuáticos en la ecorregión de Puna y altos Andes de la provincia de Jujuy, analizando los componentes alfa y beta de la diversidad de ensambles de macroinvertebrados acuáticos entre distintas cuencas (6). Se realizó un diseño a escala regional, en 19 sitios (3500 – 4600 msnm). Se realizó un análisis de cluster para clasificar los sitios en grupos con similitudes ambientales (temperatura media anual, precipitación media anual y altitud fueron significativos). Los sitios se agruparon en 3 áreas: a) Olaroz – Tuzgle, b) Río Grande – Yavi – Miraflores y c) Vilama, esta clasificación explica el 99 % de la variabilidad. La riqueza específica (S) y los valores de diversidad verdadera (DV), señalan al grupo a) como el de valores significativamente menores (4000 msnm), a diferencia de los grupos b) y c) con modelos que se ajustan significativamente mostrando un incremento de S y DV en relación al incremento de la altitud. Comparando la composición de especies entre los grupos, se encontraron diferencias en la diversidad beta dadas principalmente por un recambio de especies.