A network of six NPK long-term field trials was carried out on different soils of citrus-producing regions of São Paulo state, Brazil, in order to estimate quantitative relations of fruit yield to NPK fertilization and to determine parameters for fertilizer recommendation based on soil testing and leaf analysis. The experiments were set up in an incomplete factorial design 1/2 43 with 32 treatments, with four yearly rates of N (30, 100, 170 and 240 kg N/ha), P (9, 27, 45 and 63 kg P/ha) and K (25, 91, 157 and 223 kg K/ha). Four to seven harvests were recorded for the six experiments. Response surfaces of the type y = bo + b1N + b11N2 + b2P + b22P2 + b3K + b33K2 + b12NP + b13NK + b23PK were adjusted to the average yields of each trial. Correlation were established for yield increases, expressed as relative yields, and results of soil analysis of P and K, and leaf analysis of N. Soil samples taken at 0-20 cm depth in the beginning of each experiment were analyzed for resin extractable P and exchangeable K using an ion-exchange resin procedure. Yield responses for phosphorus and potassium applications were observed respectively in soils with less than 20 mg dm-3 of P and 20 mmolc dm-3 of K+. Yield responses to nitrogen were related to the total content of nitrogen in leaves, being largest for N values of 23 g kg-1 and smallest for N of 28 kg-1. With these field information, a practical approach for fertilizer recommendation for citrus, based on soil analysis for P and K and leaf analysis for N, was developed.