Please Numerical groundwater quality models (GQMs) often run at high computational cost resulting in long simulation times and complex parameter calibration that limit their practical applications. In this study, a novel reduced-order model (ROM) was developed for nitrate simulation in groundwater including a simple structure and with similar accuracy as more extensive GQMs. The proposed methodology for the development of ROM presents a solution for the problem in ROMs developed with eigenvectors, to make predictions into the future. The model performance was investigated by simulation of nitrate in the Karaj Aquifer, Iran. The dominant modes of spatiotemporal variation of nitrate during a five-year period was calculated by the model. The results revealed an excellent agreement between nitrate simulated by the ROM and the well-known Modular Transport 3D Multi Species (MT3DMS). The absolute error between the ROM and the MT3DMS was less than 0.5 mg/l in the most parts of the aquifer. Thus, results confirm that the use of ROM has advantages through a much simpler structure and shorter calculation times. Observed spatiotemporal variation of nitrate in the aquifer was well represented by the ROM simulations. The simplicity of the model makes it highly interesting also to other water resources problems.