The Lake Natron region is one of the potential geothermal resource areas in Tanzania, this is based on previous studies that focused on geology and the characteristics of the hot springs found in this region. However, the exact locations and sizes of the geothermal system(s)in this region are still unknown and require further exploration. This study used aeromagnetic data, Digital Elevation Model (DEM) data, multispectral remote sensing Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data that is characterized with spatial resolution of up to 15 m and extended wavelength range to Shortwave Infrared Region (SWIR) of the electromagnetic spectrum, to identify potential areas for geothermal systems. The datasets were processed to produce geological structures, indicator mineral alterations for geothermal areas and land surface temperature (LST) anomalies, which were used as proxies to the geothermal systems. Remote sensing data were processed using various techniques to improve the signals from geological features of interest. In the processing of ASTER data, the processing techniques included band ratios, false composite images, and supervised classification. Field observation was conducted to verify the results along with geological mapping and identification of geological features around geothermal manifestations. Aeromagnetic and SRTM data analysis allowed to interpret surface and subsurface structures, and compare them with the field observations. Results have revealed that the main structural trend is NE-SW, which agrees with the field evidence. This trend reflects the regional tectonic framework and the influence of major faults. Field mapping realized the presence of geothermal manifestations such as hot springs with temperatures range between 32℃ and 51℃, and pH range between 8 and 11. All hot springs are located along the mapped structures and major faults of the East African Rift System. Some parts of the study area are characterized by LST anomalies and indicator mineral alterations such as clays, sulfates, carbonates and ferric iron, these are characteristic geothermal manifestations in geothermal systems. Overall observation has also revealed a good match between remote sensing mapped mineral alterations, LST anomalies and other geothermal surface manifestations, this includes also geothermal favorable structure settings such as bend in normal faults, displacement transfer zones and horse tail structures. Despite the prominent results from this study, there is still a need for detailed investigation using high spectral and spatial resolution remote sensing data, time series thermal infrared (TIR) data, the best fit LST inversion method, along with other geological information (e.g. MT/TM geophysical surveys and geochemistry) to further validate the identified geothermal system.