The tropical thermocline plays an important role in regulating equatorial sea surface temperature (SST); at present, it is still poorly simulated in the state-of-the-art climate models. In this paper, thermocline biases in the tropical North Pacific are investigated using the newly released CMIP6 historical simulations. It is found that CMIP6 models tend to produce an overly shallow thermocline in the northwestern tropics, accompanied by a deep thermocline in the northeastern tropics. A pronounced thermocline strength bias arises in the tropical northeastern Pacific, demonstrating a dipole structure with a sign change at about 8°N. These thermocline biases are accompanied with biases in the simulations of oceanic circulations, including a too weak North Equatorial Countercurrent (NECC), a reduction in water exchanges between the subtropics and the equatorial regions, and an eastward extension of the equatorward interior water transport. The causes of these thermocline biases are further analyzed. The thermocline bias is primarily caused by the model deficiency in simulating the surface wind stress curl, which can be further attributed to the longstanding double-ITCZ bias in the tropical North Pacific. Besides, thermocline strength bias can be partly attributed to the poor prescription of oceanic background diffusivity. By constraining the diffusivity to match observations, the thermocline strength in the tropical northeastern Pacific is greatly increased.