Two levelling-based vertical datums have been used in North America, namely CGVD28 in Canada and NAVD88 in the USA and Mexico. Although the two datums will be replaced by a common and continent-wide vertical datum in a few years, their connection and unification are of great interest to the scientific and user communities. In this paper, the geodetic boundary value problem (GBVP) approach is studied as a rigorous method for connecting two or more vertical datums through computed datum offsets from a global equipotential surface defined by a GOCEbased geoid. The so-called indirect bias term, the effect of the GOCE geoid omission error, the effect of the systematic levelling datum errors and distortions, and the effect of the geodetic data errors on the datum unification are four important factors affecting the practical implementation of this approach. These factors are investigated numerically using the GNSS-levelling and tide gauge stations in Canada, the USA, Alaska, and Mexico. The results show that the indirect bias term can be omitted if a GOCE-based global geopotential model is used in gravimetric geoid computations. The omission of the indirect bias term simplifies the linear system of equations for the estimation of the datum offset(s). Because of the existing systematic levelling errors and distortions in the Canadian and US levelling networks, the datum offsets are investigated in eight smaller regions along the Canadian and US coastal areas. Using GNSSlevelling stations in the US coastal regions, the mean datum offset can be estimated with a 1 cm standard deviation if the GOCE geoid omission error is taken into account by means of the local gravity and topographic information. In the Cana-B B. Amjadiparvar