AbstractGlacial‐isostatic adjustment (GIA) is the key process controlling relative sea‐level (RSL) and paleo‐topography. The viscoelastic response of the solid Earth is controlled by its viscosity structure. Therefore, the appropriate choice of Earth structure for GIA models is still an important area of research in geodynamics. We construct 18 3D Earth structures that are derived from seismic tomography models and are geodynamically constrained. We consider uncertainties in 3D viscosity structures that arise from variations in the conversion from seismic velocity to temperature variations (factor r) and radial viscosity profiles (RVP). We apply these Earth models to a 3D GIA model, VILMA, to investigate the influence of such structure on RSL predictions. The variabilities in 3D Earth structures and RSL predictions are investigated for globally distributed sites and applied for comparisons with regional 1D models for ice center (North America, Antarctica) and peripheral regions (Central Oregon Coast, San Jorge Gulf). The results from 1D and 3D models reveal substantial influence of lateral viscosity variations on RSL. Depending on time and location, the influence of factor r and/or RVP can be reverse, for example, the same RVP causes lowest RSL in Churchill and largest RSL in Oregon. Regional 1D models representing the structure beneath the ice and 3D models show similar influence of factor r and RVP on RSL prediction. This is not the case for regional 1D models representing the structure beneath peripheral regions indicating the dependence on the 3D Earth structure. The 3D Earth structures of this study are made available.