Assimilation methods, meant to constrain divergence of model trajectory from reality using observations, do not exactly satisfy the physical laws governing the model state variables. This allows mismatches in the analysis in the vicinity of observation locations where the effect of assimilation is most prominent. These mismatches are usually mitigated either by the model dynamics in between the analysis cycles and/or by assimilation at the next analysis cycle. However, if the observations coverage is limited in space, as it was in the ocean before the Argo era, these mechanisms may be insufficient to dampen the mismatches, which we call shocks, and they may remain and grow. Here we show through controlled experiments, using real and simulated observations in two different ocean models and assimilation systems, that such shocks are generated in the ocean at the lateral boundaries of the moored buoy network. They thrive and propagate westward as Rossby waves along these boundaries. However, these shocks are essentially eliminated by the assimilation of near-homogenous global Argo distribution. These findings question the fidelity of ocean reanalysis products in the pre-Argo era. For example, a reanalysis that ignores Argo floats and assimilates only moored buoys, wrongly represents 2008 as a negative Indian Ocean Dipole year.