Earth-observing satellites, such Landsat, provide many multitemporal images of earth, either water body or land. Using spectral water body characteristics and field measurement, bathymetry (water depth) of the study area can be calculated for each image recorded/acquired at different time. Four multitemporal spatial images were used for generating bathymetry images which were arranged as multilayer images. Bathymetric forecasting is needed for a dynamic (rapidly change) area, such as an estuary of a river with a lot of sediments. Bathymetric forecasting in the study area utilized linear and quadratic regressions. Spatial image layers representing standard errors, constants of linear and quadratic equations were generated from the cubic raster database and were arranged in a cubic raster database as well . The values of the layers were classified into several classes and showed with distinctive colors to ease in visual observation. The bathymetric forecasting (forward or backward) can ...