NASA’s ICESat-2 has been providing sea ice freeboard measurements across the polar regions since October 2018. In spite of the outstanding spatial resolution and precision of ICESat-2, the spatial sparsity of the data can be a critical issue for sea ice monitoring. This study employs a geostatistical approach (i.e., ordinary kriging) to characterize the spatial autocorrelation of the ICESat-2 freeboard measurements (ATL10) to estimate weekly freeboard variations in 2019 for the entire Ross Sea area, including where ICESat-2 tracks are not directly available. Three variogram models (exponential, Gaussian, and spherical) are compared in this study. According to the cross-validation results, the kriging-estimated freeboards show correlation coefficients of 0.56–0.57, root mean square error (RMSE) of ~0.12 m, and mean absolute error (MAE) of ~0.07 m with the actual ATL10 freeboard measurements. In addition, the estimated errors of the kriging interpolation are low in autumn and high in winter to spring, and low in southern regions and high in northern regions of the Ross Sea. The effective ranges of the variograms are 5–10 km and the results from the three variogram models do not show significant differences with each other. The southwest (SW) sector of the Ross Sea shows low and consistent freeboard over the entire year because of the frequent opening of wide polynya areas generating new ice in this sector. However, the southeast (SE) sector shows large variations in freeboard, which demonstrates the advection of thick multiyear ice from the Amundsen Sea into the Ross Sea. Thus, this kriging-based interpolation of ICESat-2 freeboard can be used in the future to estimate accurate sea ice production over the Ross Sea by incorporating other remote sensing data.