[1] This work presents the results of a local empirical model that describes the behavior of the ionospheric F 2 region peak. The model was developed using nearly 25 years of incoherent scatter radar (ISR) measurements made at the Arecibo Observatory (AO) between 1985 and 2009. The model describes the variability of the F 2 peak frequency ( f o F 2 ) and F 2 peak height (h m F 2 ) as a function of local time, season, and solar activity for quiet-to-moderate geomagnetic activity conditions (Kp < 4+). Our results show that the solar activity control of h m F 2 and f o F 2 over Arecibo can be better described by a new proxy of the solar flux (F 107P ), which is presented here. The variation of h m F 2 parameter with F 107P is virtually linear, and only a small saturation of the f o F 2 parameter is observed at the highest levels of solar flux. The winter anomaly and asymmetries in the variation of the modeled parameters between equinoxes were detected during the analyses and have been taken into account by the AO model. Comparisons of ISR data with international reference ionosphere (IRI) model predictions indicate that both CCIR and URSI modes overestimate f o F 2 during the daytime and underestimate it at night. As expected, this underestimation is not observed in the AO model. Our analyses also show that the h m F 2 parameter predicted by the IRI modes shows a saturation point, which causes h m F 2 to be underestimated at high solar activity. The underestimation increases with higher levels of solar activity. Finally, we also found that IRI predictions of the seasonal variability of f o F 2 and h m F 2 over Arecibo can be improved by using a small correction that varies with solar activity and local time.Citation: Brum, C. G. M., F. S. Rodrigues, P. T. dos Santos, A. C. Matta, N. Aponte, S. A. Gonzalez, and E. Robles (2011), A modeling study of f o F 2 and h m F 2 parameters measured by the Arecibo incoherent scatter radar and comparison with IRI model predictions for solar cycles 21, 22, and 23,