Predicting the hotspots of antimicrobial resistance (AMR) in aquatics is crucial for managing associated risks. We developed an integrated modeling framework toward predicting the spatiotemporal abundance of antibiotics, indicator bacteria, and their corresponding antibiotic-resistant bacteria (ARB), as well as assessing the potential AMR risks to the aquatic ecosystem in a tropical reservoir. Our focus was on two antibiotics, sulfamethoxazole (SMX) and trimethoprim (TMP), and on Escherichia coli (E. coli) and its variant resistant to sulfamethoxazole-trimethoprim (EC_SXT). We validated the predictive model using withheld data, with all Nash-Sutcliffe efficiency (NSE) values above 0.79, absolute relative difference (ARD) less than 25%, and coefficient of determination (R 2 ) greater than 0.800 for the modeled targets. Predictions indicated concentrations of 1−15 ng/L for SMX, 0.5−5 ng/L for TMP, and 0 to 5 (log 10 MPN/100 mL) for E. coli and −1.1 to 3.5 (log 10 CFU/100 mL) for EC_SXT. Risk assessment suggested that the predicted TMP could pose a higher risk of AMR development than SMX, but SMX could possess a higher ecological risk. The study lays down a hybrid modeling framework for integrating a statistic model with a process-based model to predict AMR in a holistic manner, thus facilitating the development of a better risk management framework.