The spatial sparsity and temporal discontinuity of station-based SAT data do not allow to fully understand Antarctic surface air temperature (SAT) variations over the last decades. Generating spatiotemporally continuous SAT fields using spatial interpolation represents an approach to address this problem. This study proposed a backpropagation artificial neural network (BPANN) optimized by a genetic algorithm (GA) to estimate the monthly SAT fields of the Antarctic continent for the period 1960–2019. Cross-validations demonstrate that the interpolation accuracy of GA-BPANN is higher than that of two benchmark methods, i.e., BPANN and multiple linear regression (MLR). The errors of the three interpolation methods feature month-dependent variations and tend to be lower (larger) in warm (cold) months. Moreover, the annual SAT had a significant cooling trend during 1960–1989 (trend = −0.07°C/year;
p
=
0.04
) and a significant warming trend during 1990–2019 (trend = 0.06°C/year;
p
=
0.05
). The monthly SAT did not show consistent cooling or warming trends in all months, e.g., SAT did not show a significant cooling trend in January and December during 1960–1989 and a significant warming trend in January, June, July, and December during 1990–2019. Furthermore, the Antarctic SAT decreases with latitude and the distance away from the coastline, but the eastern Antarctic is overall colder than the western Antarctic. Spatiotemporal inconsistencies on SAT trends are apparent over the Antarctic continent, e.g., most of the Antarctic continent showed a cooling trend during 1960–1989 (trend = −0.20∼0°C/year;
p
=
0.01
∼
0.27
) with a peak over the central part of the eastern Antarctic continent, while the entire Antarctic continent showed a warming trend during 1990–2019 (trend = 0∼0.10°C/year;
p
=
0.04
∼
0.42
) with a peak over the higher latitudes.