Over the last few decades, the development of new GIS-based computer algorithms, automated relief analysis tools, and free access to datasets have attracted the attention of researchers to study the morphology of the Earth's surface for various purposes. For instance, the automated extraction of geomorphological settings using geoinformatics techniques has become a crucial aspect of present environmental analysis and paleo landscape modelling efforts. Moreover, morphological data obtained from highly accurate digital elevation models (DEMs), such as LiDAR-derived DEMs, can provide valuable information related to landscape modelling and landform classification processes. In geoarchaeological investigations, the different types of landforms ranging from large-scale features such as plains and mountains to local topography such as small hills and valleys, play a significant role in human behavior, particularly in terms of habitation practices. Therefore, in this work, we aimed to apply an automated GIS-based landform classification method to characterize the relief of over 350 Late Bronze Age (LBA) settlements belonging to the Noua Culture (NC) (cal. 1500/1450-1100 BCE) located in northeast Romania. For this purpose, we used an adapted version of the Topographic Position Index (TPI) methodology developed by Weiss A.D. and implemented as an extension for ESRI ArcView 3.x. by Jenness J., which we abbreviated as DEV. This approach involved two steps: (i) applying the standard deviation of TPI for the mean elevation (DEV) around each LBA site analyzed (1000 m buffer zone) and (ii) classifying the archaeological site's location using six slope position classes (first method) or ten morphological classes by combining the parameters from two small-DEV and large-DEV neighborhood sizes (second method). Our findings suggest that the populations belonging to the Noua Culture preferred to establish their settlements on hilltops close to steep slopes and/or on the small hills located in large valleys, indicating a strong connection between the landform patterns and habitation practices during the Late Bronze Age. From the perspective of the GIS-based methodology applied in this study, the automated relief analysis using TPI-based DEV combined with LiDAR-derived DEMs and other geomorphological variables (e.g., terrain slope) can be easily integrated into various geoarchaeological studies due to its great applicability. Furthermore, the TPI-based DEV techniques bring many improvements to conventional geoarchaeological surveys: the fast and low-cost performance of relief analysis at both small and large scales, the ability to divide the landscape into ten landform classes and replicate the analysis for various archaeological contexts, and the ability to describe prehistoric human behavior based on certain geographical datasets. Additionally, the outcomes contribute to archaeological predictive modelling in the heterogeneous landscape of Northeast Romania.