Global warming will have an impact on nature in many ways, including rising sea levels and an increasing spread of infectious diseases. Land surface temperature is one of the many indicators that can be used to measure climate change on both a local and global scale. This study aims to analyze the change in land surface temperatures on New Guinea Island using a cubic spline method, autoregressive model, and multivariate regression. New Guinea Island was divided into 5 regions each consisting of 9 subregions. The data of each subregion was obtained from the National Aeronautics and Space Administration moderate resolution imaging spectroradiometer database from 2000 to 2019. The average change in temperature was +0.012o C per decade. However, the changes differed by region; significantly decreasing in the northwest at -0.107o C per decade (95% CI: -0.207, -0.007), significantly increasing in the south at 0.201o C per decade (95% CI: 0.069, 0.333), and remaining stable in the centralnorth, southeast and northeast.