1. A computer-implemented method for estimating uncertainty in a bathymetry database comprising the steps of: obtaining conditional probability densities of bathymetric uncertainty due to horizontal error and bottom slope using a Monte Carlo technique on representative sets of soundings data from the bathymetry database, the representative sets being associated with a training set area; training a Bayesian Network to provide causal relationships between the horizontal error and the bottom slope and the bathymetric uncertainty based on the conditional probability densities; creating a first histogram of the bathymetric uncertainty from the Bayesian Network for an area similar to the training set area; and estimating an uncertainty based on the first histogram.