Results obtained as a solution to the geodetic boundary-value problem, where two different sources of input boundary data are combined together, are presented. On the upper part of the boundary defined by the Earth surface and separated from the Earth masses, the gravity anomalies are derived from GOCE gravity gradiometry data using a transformation and downward continuation procedure. On the bottom part of the boundary, the gravity anomalies derived from the terrestrial or marine gravity measurements are used. The boundary-value problem is then solved numerically by the finite element method. In some cases an independency of our solution from the highdegree global gravity models (e.g. EGM2008) can be considered as the advantage. The standard deviation of our model tested by GNSS/leveling method in Auvergne area is 10 cm.