Improved modeling of sea level patterns by incorporating self‐attraction and loading

[1] We implement the effects of gravitational self-attraction and loading (SAL) into a global baroclinic ocean circulation model and investigate effects on sea level patterns, ocean circulation, and density distributions. We compute SAL modifications as an additional force on the water masses at every time step by decomposing the field of ocean bottom pressure anomalies into spherical harmonic functions and then applying Love numbers to account for the elastic properties of the solid Earth. Considering SAL in the postprocessing turns out to be insufficient, especially in coastal waters and on subweekly time scales, where SAL modifies local sea level by around 0.6–0.8 cm on average; in the open ocean, changes mostly remain around 0.3 cm. Modifications of water velocities as well as of heat and salt distributions are modeled, yet they are small. Simple parameterizations of SAL effects currently used in a number of ocean circulation models suffer from the process's inhomogeneity in space and time. These parameterizations improve the modeled sea level patterns but fail to reproduce SAL impacts on circulation and density distributions. We therefore suggest to explicitly consider the full SAL effect in ocean circulation models, especially when investigating sea level variations faster than around 4 days.