Bathymetric effect on geoid modeling over the Great Lakes area

Bathymetry data over lake areas are not included in the current and previous NGS (National Geodetic Survey) geoid models. Lake surfaces are simply treated as land surfaces during the modeling regardless of the apparent density difference between water and rock, resulting in artificial masses that distort the model from the actual gravity field and the corresponding geoid surface. In this study, compiled high-resolution bathymetry data provided by National Centers for Environmental Information are used to identify the real volume of water bodies. Under the mass conservation principle, two strategies are deployed to properly account the water body bounded by the mean lake surface and the bathymetry indicated lake floor into the current NGS geoid modeling scheme, where the residual terrain modeling method is used to account for topographic effects. The first strategy condenses water bodies into equivalent rock masses, with the cost of changing the geometrical shape of the water body. The second one keeps the shape of the water body unchanged but replaces the water and rock densities inside each topographical column bounded by the geoid surface and the mean lake surface by an averaged density. Both strategies show up to 1-cm geoid changes when compared with the previous geoid model that does not consider bathymetric information. All three geoid models are evaluated by local GNSS/Leveling benchmarks and multi-year-multi-mission altimetry indicated mean lake surface heights. The results show that both strategies can improve the geoid model precision. And the second strategy yields more realistic results. Graphical Abstract