With the goal to support effective water resource management, water quality models have gained popularity as tools for evaluating the distributions of pollutants and sediments. This work focuses on the application of the numerical solution of an advection-dispersion-reaction (ADR) water quality model for rivers and streams to a major Philippine waterway, the Pasig River. The water quality constituent is described by a system of reaction and advection-dispersion-reaction equations. The model and method are based on a previously used strategy where Guass-Jordan decomposition is applied to the matrix system and the resulting conservative form of the model is solved numerically using the fully implicit scheme and finite element method. The methodology is demonstrated by a case study in Pasig River involving the concentrations of total dissolved solids (TDS) obtained from the Department of Environment and Natural Resources (DENR) through the Pasig River Unified Monitoring Stations (PRUMS) report. Sensitivity analysis and parameter estimation are also applied to the model to assess which parameters influence the model output the most.