本文由中国气象局气候变化专项CCSF201313,CCSF201410,CCSF201606、黑龙江省气象局自立项目HQ2015018和中国气象局核心业务发展专项CMAHX20160205资助;This work was supported by the Special Project in Climate Change of China Meteorological Administration CCSF201313,CCSF201410,CCSF201606,Self-reliance Project of Heilongjiang Province Meteorological Bureau HQ2015018 and Core Business Development Special Program of China Meteorological Administration CMAHX20160205.