A simplified quasigeostrophic (QG) analytical model together with an idealized numerical model are used to study the effect of uneven ice-ocean stress on the temporal evolution of the geostrophic current under sea ice. The tendency of the geostrophic velocity in the QG model is given as a function of the lateral gradient of vertical velocity and is further related to the ice-ocean stress with consideration of a surface boundary layer. Combining the analytical and numerical solutions, we demonstrate that the uneven stress between the ice and an initially surface-intensified, laterally sheared geostrophic current can drive an overturning circulation to trigger the displacement of isopycnals and modify the vertical structure of the geostrophic velocity. When the near-surface isopycnals become tilted in the opposite direction to the deeper ones, a subsurface velocity core is generated (via geostrophic set up). This mechanism should help understand the formation of subsurface currents in the edge of Chukchi and Beaufort Seas seen in observations. Furthermore, our solutions reveal a reversed flow extending from the bottom to the mid-depth, suggesting that the ice-induced overturning circulation potentially influences the currents in the deep layers of the Arctic Ocean, such as the Atlantic Water boundary current.