Vietnam is one of the most affected countries by the global sea-level rise, in which the Mekong Delta is the most heavily impacted area due to the low-lying position. The study aims to simulate and delineate the changing of elevation by combining spatial geostatistics and Kriging interpolation techniques. The assumption of inland water level increases at different sea-level rise scenarios by interpolating 967 elevations geo-rectified data from the topographic map at 1/250,000 scale. The result showed that the exponential is the best model for interpolation of elevation data, distance A0 = 585,900 m; coefficient R2 = 0.985; and residue sum of square = 0.0011. Inland inundated water levels maps with sea-level rise scenarios from 0.2 to 2.8 m were delineated. The submersion started when the water level rose to 0.6 m and the total submerged area gradually increasing; until the water levels rose to 2.8 m, the whole region completely flooded. The submerged areas increase, resulting in the loss of cultivated land and food production. However, further research needs to be considered a more detailed study on population, food production, and tidal due to land loss. This study proved that spatial analysis is an effective tool for simulating surface water levels. It is of great significance for controlling water levels changes in the Delta