Geological studies in the Persian Gulf (PG) have revealed the existence of subseabed saltdomes. With suitable filtering of a high-resolution PG seabed topography, it is seen that the domes leave their signature in the seabed, i.e., numerous hills and valleys with amplitudes of several tens of meters and radii from a few up to tens of kilometers. It was suspected that the ``shark skin{''} of the PG seabed may affect the tidal residual flow. The interaction of tidal dynamics and these obstacles was investigated in a nonlinear hydrodynamic numerical tidal model of the PG. The model was first used to characterize flow patterns of residual currents generated by a tidal wave passing over symmetric, elongated and tilted obstacles. Thereafter it was applied to the entire PG. The model was forced at its open boundary by the four dominant tidal constituents residing in the PG. Each tidal constituent was simulated separately. Results, i.e., tidal residual currents in the PG, as depicted by Lagrangian trajectories reveal a stationary flow that is very rich in eddies. Each eddy can be identified with a topographic obstacle. This confirms that the tidal residual flow field is strongly influenced by the nonlinear interaction of the tidal wave with the bottom relief which, in turn, is deformed by salt-domes beneath the seabed. Different areas of maximum residual current velocities are identified for major tidal constituents. The pattern of trajectories indicates the presence of two main cyclonic gyres and several adjacent gyres rotating in opposite directions and a strong coastal current in the northern PG.