Greater variability in renal function is associated with mortality in patients with chronic kidney disease (CKD). However, few studies have demonstrated the predictive value of renal function variability in relation to renal outcomes. This study investigates the predictive ability of different methods of determining estimated glomerular filtration rate (eGFR) variability for progression to renal replacement therapy (RRT) in CKD patients. This was a prospective observational study, which enrolled 1,862 CKD patients. The renal end point was defined as commencement of RRT. The variability in eGFR was measured by the area under the eGFR curve (AUC)%. A significant improvement in model prediction was based on the -2 log likelihood ratio statistic. During a median 28.7-month follow-up, there were 564 (30.3%) patients receiving RRT. In an adjusted Cox model, a smaller initial eGFR AUC%_12M (P < 0.001), a smaller peak eGFR AUC%_12M (P < 0.001), and a larger negative eGFR slope_12M (P < 0.001) were associated with a higher risk of renal end point. Two calculated formulas: initial eGFR AUC%_12M and eGFR slope_12M were the best predictors. Our results demonstrate that the greater eGFR variability by AUC% is associated with the higher risk of progression to RRT.