The decadal variability of the North Atlantic thermohaline circulation (THC) is investigated within a three-dimensional ocean circulation model using the conditional nonlinear optimal perturbation method. The results show that the optimal initial perturbations of temperature and salinity exciting the strongest decadal THC variations have similar structures: the perturbations are mainly in the northwestern basin at a depth ranging from ?1500 to ?3000 m. These temperature and salinity perturbations act as the optimal precursors for future modifications of the THC, highlighting the importance of observations in the northwestern basin to monitor the variations of temperature and salinity at depth. The decadal THC variation in the nonlinear model initialized by the optimal salinity perturbations is much stronger than that caused by the optimal temperature perturbations, indicating that salinity variations might play a relatively important role in exciting the decadal THC variability. Moreover, the decadal THC variations in the tangent linear and nonlinear models show remarkably different characteristics, suggesting the importance of nonlinear processes in the decadal variability of the THC.