Accurately estimating aboveground biomass (AGB) in forest ecosystems facilitates efficient resource management, carbon accounting, and conservation efforts. This study examines the relationship between predictors from Landsat-9 remote sensing data and several topographical features. While Landsat-9 provides reliable data crucial for long-term monitoring, it is part of a broader suite of available remote sensing technologies. We employ machine learning algorithms such as Extreme Gradient Boosting (XGBoost), Support Vector Regression (SVR), and Random Forest (RF), alongside linear regression techniques like Multiple Linear Regression (MLR). The primary objectives of this study encompass two key aspects. Firstly, the research methodically selects optimal predictor combinations from four distinct variable groups: Landsat-9 (L1) data, a fusion of Landsat-9 data and Vegetation-based indices (L2), and the integration of Landsat-9 data with the Shuttle Radar Topography Mission Digital Elevation Model (SRTM DEM) (L3) and the combination of best predictors (L4) derived from L1, L2, and L3. Secondly, the research systematically assesses the effectiveness of different algorithms to identify the most precise method for establishing any potential relationship between field-measured AGB and predictor variables. Our study revealed that the Random Forest (RF) model was the most efficient method utilizing Landsat-9 OLI and SRTM DEM (L3) predictors, achieving remarkable accuracy. This conclusion was reached by assessing its outstanding performance when compared to an independent validation dataset. The RF model exhibited remarkable accuracy, presenting relative mean absolute error (RMAE), relative root mean square error (RRMSE), and R2 values of 14.33%, 22.23%, and 0.81, respectively. The XGBoost model is the subsequent choice with RMAE, RRMSE, and R2 values of 15.54%, 23.85%, and 0.77, respectively. The study further highlights the significance of specific spectral bands, notably B4 and B5 from Landsat 9 OLI data, in capturing spatial AGB distribution patterns. Integration of vegetation-based indices, including TNDVI, NDVI, RVI, and GNDVI, further refines AGB mapping precision. Elevation, slope, and the Topographic Wetness Index (TWI) are crucial proxies for representing biophysical and biological mechanisms impacting AGB. Through the utilization of openly accessible fine-resolution data and employing the RF algorithm, the research demonstrated promising outcomes in the identification of optimal predictor-algorithm combinations for forest AGB mapping. This comprehensive approach offers a valuable avenue for informed decision-making in forest management, carbon assessment, and ecological monitoring initiatives.