Technical Abstract: Many complex agricultural system models have been developed, but a detailed description of expectations in model response and how they may be applied are rarely provided. We apply a sensitivity and uncertainty [first order error analysis (FOA) and Monte Carlo Latin Hypercube Sampling (LHS) simulation] analysis framework for evaluation of the WEPP and RZWQM erosion and water quality models. Sensitivity analysis results show the WEPP model response for runoff to be most sensitive to effective hydraulic conductivity and soil parameters used in the crusting factor adjustment. The soil loss response was most sensitive to erodibility factors and soil and management parameters influencing infiltration. RZWQM plant uptake, yield, and nitrate leaching output responses were most sensitive to plant growth input parameters and manure application rates. RZWQM output responses were also more sensitive to the average saturated conductivity of the entire soil profile than to individual soil layers. The FOA did not approximate the WEPP model responses for runoff and soil loss well due to model nonlinearity. Recommendations based on sensitivity analysis in conjunction with Monte Carlo LHS simulation are presented for evaluating agricultural system models. Derivation of output response error variances through Monte Carlo LHS simulation is also suggested for agricultural system models.