Abstract This paper focuses on implementing a decentralized rainwater harvesting (DRH) program for rural cities under stakeholder interests and purchasing power restrictions. The study is conducted by considering: (1) a no tax incentive scheme (NTIS); (2) a flat rate tax incentive scheme (FTIS); and (3) a progressive tax incentive scheme (PTIS). A three-stage decision framework is proposed. In the first stage, a water balance simulation model is developed to estimate the amount of rainwater that can be collected at individual households using different capacity rain barrels and to determine the threshold rain barrel capacity. In the second stage, the water balance simulation model is run multiple times with varying input parameters to obtain the probability distributions for rainwater used and water demand for different households using different capacity rain barrels, which are then used to generate scenarios for the third stage optimization model. In the third stage, a multi-objective stochastic mixed integer linear programming (Mo-SMILP) model is proposed that aims to determine the household participation rate and cumulative city-wide DRH capacity considering tax incentive schemes, stakeholder interests, and purchasing power restrictions. The proposed Mo-SMILP model is solved by using an improved augmented epsilon constraint method (AUGMECON 2) in order to understand the trade-off between city-wide household reliability, city-wide household savings, and city government savings. A case study of a small rural city in the State of Texas in the US is used to illustrate the effectiveness of the proposed decision framework. The results shows a significant decrease in city-wide household reliability when stakeholder interests and purchasing power are considered. In addition, the PTIS is best at improving city-wide household reliability and savings. NTIS and FTIS are best at improving city government savings depending on the household preference for reliability and savings.