Using the observational sea surface temperature(SST) data,the interannual variation of winter North Pacific Oscillation(NPO) is simulated by the National Center for Atmospheric Research(NCAR) Community Atmosphere Model Version 5(CAM5) with or without the inclusion of wind gustiness effect in the parameterization scheme of air-sea turbulent heat fluxes.Results show that,with the inclusion of wind gustiness effect,the interannaul variation of winter NPO can be reasonably reproduced by the model,and the correlation coefficient between 1979—1999 observed and simulated NPO is increased from 0.09 to 0.57.The SST anomaly in the eastern tropical Pacific plays an important role on the interannual variation of winter NPO,by influencing the air-sea turbulent heat fluxes.By considering the wind gustiness effect in the parameterization scheme of air-sea turbulent heat fluxes in CAM5,the simulating capability of interannual variation of sea surface turbulent heat fluxes in the eastern tropical Pacific is significantly improved,leading to the simulations of interannual variation of winter NPO be improved significantly.