Increased temperature is one of the signals of global warming. Trends in land surface temperature can be used to measure climate change. This research aimed to investigate the variation of land surface temperature in Borneo island using a cubic spline method and a multivariate regression model. The island was divided into 8 regions each comprising 9 subregions. Land surface temperatures for each subregion from 2000 to 2019 was obtained from the National Aeronautics and Space Administration moderate resolution imaging spectroradiometer database. The average increase in temperature was 0.2oC per decade with a 95% confidence interval of (0.14, 0.27). The changes differed by region; a significant increase was seen in Sarawak, North Kalimantan, West Kalimantan, West-central Kalimantan, Central-east Kalimantan region, a slight decrease in Sabah and Brunei Darussalam (Sabah & Brunei) region, a slight increase in East Kalimantan and a stable trend in South Kalimantan.