Computer models have been widely used to evaluate the impact of agronomic management on nitrogen (N) dynamics in subsurface drained fields. However, they have not been evaluated as to their ability to capture the variability of nitrate-nitrogen (NO₃–N) concentration in subsurface drainage at a wide range of N application rates due to possible errors in the simulation of other system components. The objective of this study was to evaluate the performance of Root Zone Water Quality Model2 (RZWQM2) in simulating the response of NO₃–N concentration in subsurface drainage to N application rate. A 16-yr field study conducted in Iowa at nine N rates (0–252 kg N ha⁻¹) from 1989 to 2004 was used to evaluate the model, based on a previous calibration with data from 2005 to 2009 at this site. The results showed that the RZWQM2 model performed “satisfactorily” in simulating the response of NO₃–N concentration in subsurface drainage to N fertilizer rate with 0.76, 0.49, and −3% for the Nash-Sutcliffe efficiency, the ratio of the root mean square error to the standard deviation, and percent bias, respectively. The simulation also identified that the N application rate required to achieve the maximum contaminant level for the annual average NO₃–N concentration was similar to field-observed data. This study supports the use of RZWQM2 to predict NO₃–N concentration in subsurface drainage at various N application rates once it is calibrated for the local condition.