Seepage plays an important role in soil erosion in contour ridge systems. Seepage generated from subsurface flow causes hillslope instability by reducing the soil shear strength and mobilizing soil particles and can produce cave-like features called seepage undercuts that can lead to contour failure. As the main threat to contour plowing, seepage results in large amounts of soil erosion due to row grade. Models used to predict seepage over a time series will provide a basis for modeling soil erosion resulting from seepage in contour ridging systems. Understanding seepage and its effects will advance our knowledge regarding seepage erosion mechanisms in contour ridge systems. In this study, 23 treatments were arranged using an orthogonal rotatable central composite design to model a seepage time series, build a simple seepage prediction model and investigate the effects of row grade, field slope and ridge height on seepage discharge. Most of the seepage discharge time series followed an S-shaped curve. The seepage discharge processes were fit by an exponential model with a determination coefficients (R2) greater than 0.995. Furthermore, the physical meaning of the exponential model was consistent with the experimental results. The seepage discharge continuously increased before the inflection point and then decreased. Finally, the seepage discharge approached a steady value. The maximum seepage discharge growth rate was achieved within 14 min, and the seepage discharge became steady within 106 min. Second-order polynomial regression models were used to determine the total and predicted steady seepage discharge using independent variables of row grade, field slope and ridge height, which produced R2 values of 0.66 (p