This study investigates the hypsometric integral (HI) correlations with lithology and neotectonics from northern Romanian Carpathians, Rodna Mountains. Two types of DEM were used in the analysis (a 30 m DEM resolution with specific grid size and a non-gridded 30 m DEM resolution). Several methods have been applied in order to calculate HI: the CalHypso GIS extension was used to automatically extract multiple hypsometric curves from digital elevation model (DEM) and to calculate the main statistics related to the HI by applying polynomial fits; local indices of spatial autocorrelation were applied using Anselin Local Moran’s I and Getis-Ord Gi* in order to see if HI distribution has spatial patterns values and selection of valid squares (500-m2 grid using 30-m DEM using ArcMap 10.1) for the lithological analysis. The values of the hypsometric kurtosis density, kurtosis, skew and density skew are increasing eastwards. Values of hypsometric skewness for the northern slope are in range between 0.290-0.816 and southern part between 0.414-0.507. The results show that HI values are higher than 0.5 in the areas with recent tectonic influences especially in northern part near Dragoş Vodă fault system. The Oligocene, Miocene, Pleistocene and Holocene formations have low values (0.5). By comparing the HI values from the two analyses (500 m grid size DEM and the DEM without grid size) of the DEM we find slight differences in the HI values, especially in the northern part of the range in Repedea, Negoiescu Mare, Fântâna and Lala basins.