With a global GSSTF2 and NCEP/NCAR reanalysis database and observation data at the Yong Xing station of Xisha Island in the South China Sea, we simulated the turbulent sensible and latent heat flux at sea surface in Chinese and neighboring seas (hereafter termed as China seas) using a common bulk method with some improved parameters. Comparing the simulated results with the observed and reanalyzed data, the improvement yielded higher accuracy, a smaller mean square deviation within 10 W/m2, and a smaller average relative error at about 25%. In addition, spatial resolution was improved to 0.1°×0.1°. The simulation is able to replay the main features of regional and seasonal variation in turbulent heat fluxes, and also the general pattern of heat flux changes during the summer monsoon outbreak in the South China Sea. With a global GSSTF2 and NCEP / NCAR reanalysis database and observation data at the Yong Xing station of Xisha Island in the South China Sea, we simulated the turbulent sensible and latent heat flux at sea surface in Chinese and neighboring seas (hereafter termed as China seas Comparing the simulated results with the observed and reanalyzed data, the improvement yielded higher accuracy, a smaller mean square deviation within 10 W / m2, and a smaller average relative error at about 25%. In addition, spatial resolution was improved to 0.1 ° × 0.1 °. The simulation is able to replay the main features of regional and seasonal variation in turbulent heat fluxes, and also the general pattern of heat flux changes during the summer monsoon outbreak in the South. China Sea.