|| We present an inversion method, called "column inversion", to retrieve two-dimensional variation of the depth to the Moho discontinuity from Bouguer gravity anomaly data, assuming that a density jump between the lower crust and upper mantle is regionally invariant. The inversion method is applied to reveal the apatial undulation of the Moho discontinuity in the Chugoku and Shikoku districts, southwest Japan. The result shows that the Moho is deep in the land area and shallow in the sea area. The Moho is anomalously deep in the west regions of Shikoku and Chugoku districts, where the prominent negative gravity anomaly is observed. The two-dimensional variation of the Moho depth is quite similar to that determined by inversion of the first P-arrival time data. The column inversion of the Bouguer gravity anomaly data is proved to be useful for determining the fine structure of the Moho discontinuity undulation.