范伶俐,張福穎
① 廣東海洋大學,廣東 湛江 524088;② 南京信息工程大學,江蘇 南京 210044
El Ni?o事件次年中國東部夏季氣溫分布特征研究
范伶俐①*,張福穎②
① 廣東海洋大學,廣東 湛江 524088;② 南京信息工程大學,江蘇 南京 210044
2016-08-15收稿,2016-10-13接受
廣東海洋大學“創(chuàng)新強校工程”建設項目(GDOU2016050245)
利用HadISST和NCEP/NCAR再分析資料,采用統(tǒng)計學方法分析了1949年以來El Ni?o事件次年夏季中國東部氣溫分布特征。結果表明,混合型El Ni?o次年夏季,中國東部氣溫距平表現(xiàn)為山東—淮河流域及華南沿海出現(xiàn)負距平,其余地區(qū)為正距平;東部型,則表現(xiàn)為中國東部地區(qū)整體的負距平氣溫,且北方降溫幅度明顯大于長江以南地區(qū);中部型,主要表現(xiàn)為長江以北氣溫負距平遠大于長江以南。結合環(huán)流形勢分析發(fā)現(xiàn),三類El Ni?o事件次年夏季的赤道地區(qū)(5°S~5°N)異常緯向垂直環(huán)流均表現(xiàn)為三圈型,但3類異常Walker環(huán)流不僅強度有差異,緯向范圍也有所不同。500 hPa上西太平洋副高的強度和范圍,是造成三類El Ni?o事件次年夏季中國東部地區(qū)氣溫距平分布差異的主要原因之一,850 hPa上西南暖濕氣流的輸送對其也存在明顯影響。
中國東部
El Ni?o
夏季氣溫
大氣環(huán)流
近幾十年來,由于人類活動及自然因素等綜合影響,全球變暖已成為不爭的事實。在全球變暖的背景下,極端天氣氣候事件趨強趨多,尤其是區(qū)域性極端高溫事件頻繁發(fā)生,給社會、經(jīng)濟和人民生活造成了嚴重的影響和損失(徐金芳等,2009)。2013和2016年夏季中國東部高溫熱浪事件對人民的生活及生產(chǎn)造成了重要影響。研究夏季氣溫的變化規(guī)律及影響機理有助于提高高溫事件的預測準確率,同時為政府部門應對提供依據(jù),在減少生命財產(chǎn)損失方面具有積極的意義。
近年來我國夏季高溫熱浪的頻次、日數(shù)和強度總體呈增多、增強趨勢(葉殿秀等,2013)。1961—2003年中國極端高溫日數(shù)的兩個大值區(qū)分別位于中國的西北和東南地區(qū),大部分站點極端高溫日數(shù)是增多的(Zhang et al.,2008)。黃丹青和錢永甫(2008)指出,全國極端高溫日數(shù)的線性趨勢為正趨勢,內蒙古和新疆北部為大趨勢區(qū),西南則為小趨勢區(qū)。周長艷等(2013)指出,四川高溫天數(shù)50 a來總體呈顯著增多趨勢,20世紀90年代以來增多趨勢尤為顯著。張志薇等(2011)研究發(fā)現(xiàn)華中地區(qū)近年來年平均高溫日數(shù)和高溫日平均最高氣溫在20世紀90年代之后呈現(xiàn)上升的趨勢,尤其在2000年之后,表現(xiàn)得更加明顯。張茜和李棟梁(2011)分析東北夏季氣溫的新特征時發(fā)現(xiàn),近60 a來東北夏季氣溫總體呈上升趨勢,尤其是1990年以后增溫明顯,區(qū)域差異主要表現(xiàn)在東北部升溫早而西南部遲。以上的研究均表明近年來中國整體的夏季平均氣溫和極端高溫日數(shù)的增多趨勢是確定的,但由于不同地區(qū)受到的環(huán)流形勢及各地地形的影響不同,中國夏季高溫事件的空間分布存在地區(qū)差異。
對高溫天氣的成因也做了諸多研究。楊群等(2011)、陳麗華等(2010)和嚴文蓮等(2012)均發(fā)現(xiàn)西太平洋副高的持續(xù)異常是導致南方異常高溫的主要原因。方宇凌和簡茂球(2011)發(fā)現(xiàn),在2003年6月29日—8月15日的3段持續(xù)性高溫期間,大氣升溫主要是由非絕熱加熱造成,而水平溫度平流對升溫起負貢獻。鄒海波等(2015)發(fā)現(xiàn)2013年盛夏中國中東部地區(qū)溫度異常偏高,主要是由中低層異常偏強的下沉運動造成的。許多學者(楊義文等,2004;劉春玲等,2005;雷楊娜等,2009;史軍等,2009;朱周平和錢炳強,2010)還發(fā)現(xiàn),高溫天氣與海溫密切相關。
這些研究豐富了對全球變暖背景下中國夏季高溫的認識。但對不同類型El Ni?o次年熱帶太平洋—東亞的海氣耦合過程變化及其對中國東部夏季氣溫的影響,關注較少。那么,在不同類型El Ni?o事件次年,中國東部夏季氣溫的空間分布特征是否相同?對應的熱帶太平洋—東亞海氣特征如何?為了解決上述問題,本文根據(jù)中國氣象局國家氣候中心確定的El Ni?o事件分類,通過分析合成的三類El Ni?o事件次年的海氣要素分布特征,細致對比在東部型、中部型和混合型El Ni?o次年,中國東部夏季氣溫空間分布特征的差異,為進一步探討中國東部夏季氣溫空間分布規(guī)律,深入理解三類 El Ni?o事件對其影響的機理提供科學基礎,為提升高溫天氣的預報、預測水平服務。
資料包括:1)海表溫度資料集HadISST(Hadley Centre Sea Ice and Sea Surface Temperature data set)中的全球月平均海表溫度資料,空間分辨率為1°×1°(Rayner et al.,2003)。2)逐日和逐月的美國國家環(huán)境預測中心/國家大氣環(huán)境中心(NCEP/NCAR)的再分析資料,包括氣溫逐月資料、多層逐月風場、垂直速度等再分析資料,空間分辨率為2.5°×2.5°(Kalnay et al.,1996)。
夏季的定義為6—8月(JJA)。中國東部是指105°E以東的中國區(qū)域。研究時段選取為1949—2016年。
海洋作為大氣的下邊界,對氣候的形成和變化具有重要的影響。本節(jié)探討三類El Nio事件次年,夏季海表溫度異常(Sea Surface Temperature Anomaly,SSTA)分布特征。
表1 1949年以來El Nio事件發(fā)生時間、類型劃分及事件影響年(新指數(shù))
Table 1 El Nio events lists according to spatial categories(new index)
序號起止年月長度/月SSTA峰值/℃峰值月份∑SSTA/℃海溫強度爆發(fā)類型影響年11957年4月—1958年3月121.119.1中等E+C195821963年7月—1963年12月60.682.9極弱E196431965年6月—1966年2月91.1127.5中等E196641968年12月—1969年6月70.954.5弱C196951972年4月—1973年2月111.61212.1中等E197361982年5月—1983年8月162.212月—次年1月20.4極強C198471986年10月—1988年1月161.39月15.4強C198881991年5月—1992年5月131.2411.0中等C199291994年10月—1995年2月51.0124.1弱C1995101997年4月—1998年5月142.51223.5極強E+C1998112002年6月—2003年2月91.3116.7弱C2003122006年8月—2007年1月61.1125.2弱E+C2007132009年6月—2010年4月111.3129.2中等E2010142014年9月—2016年5月212.92015年11月30.2極強E+C2016
注:資料來自中國氣象局國家氣候中心(ENSO_Characteristic_Table.pdf).
圖1 合成的El Nio事件次年夏季的海溫距平分布(單位:℃) a.混合型;b.東部型;c.中部型Fig.1 Composite SSTA of (a)EP+CP El Nio events (b)EP El Nio events and (c)CP El Nio events during the following summer(units:℃)
圖2 合成的三類El Nio事件次年夏季中國東部氣溫距平分布(單位:℃) a.混合型;b.東部型;c.中部型Fig.2 Composite air temperature anomaly over East China of (a)EP+CP El Nio events (b)EP El Nio events and (c)CP El Nio events during the following summer(units:℃)
氣溫作為一種氣象要素,其變化必然受到大氣環(huán)流形勢的影響和制約。什么原因造成了三類 El Nio事件次年夏季中國東部氣溫分布差異?為了弄清楚這一問題,首先從Walker環(huán)流出發(fā)診斷三類El Nio事件次年夏季中國東部氣溫分布差異(圖3)。
3.1 Walker環(huán)流特征分析
圖3 合成的El Nio次年夏季5°S~5°N平均Walker環(huán)流距平(相對于的1981—2010年氣候平均場)的分布(陰影表示擴大至1 000倍的垂直速度,單位:10 hPa/s) a.混合型;b.東部型;c.中部型Fig.3 Composite anomalous(relative to the 1981—2010 climate mean field) Walker circulation between 5°S—5°N of (a)EP+CP El Nio events (b)EP El Nio events and (c)CP El Nio events during the following summer(vertical velocity multiplied by a factor of 1 000;units:10 hPa·s-1)
3.2 500 hPa位勢高度場環(huán)流特征分析
其次,分析三類事件次年夏季其對應的500 hPa高度場異常情況。從氣候平均的夏季500 hPa位勢高度場(圖4)可以看出,我國東部不同地區(qū)夏季受到各天氣系統(tǒng)影響的程度是不同的,東北地處高緯,夏季受到西太副高的直接影響較小,主要受鄂霍次克海阻高影響;華北受低緯的副高及中高緯的阻高的直接影響都不大,但中高緯該地區(qū)淺槽的深淺、冷空氣南下頻次的變化會受到中高緯經(jīng)向環(huán)流調整的影響,從而影響華北氣溫;江淮地區(qū)受副高的直接影響比華北大,再加上印度洋和孟加拉灣的水汽在青藏高原的影響下沿副高的西北側北上,影響江淮的氣溫;華南地區(qū)的高溫日數(shù)主要受到西太平洋副高的影響。
圖4 1981—2010年夏季氣候平均50 hPa位勢高度場(單位:gpm)Fig.4 Cimate mean 500 hPa potential height field during the period of 1981—2010(unit:gpm)
可見,中國東部各地區(qū)夏季受到主要環(huán)流系統(tǒng)的影響不同,直接影響系統(tǒng)的強弱對各地區(qū)氣溫分布特征起到?jīng)Q定性作用,當然,不同的環(huán)流系統(tǒng)之間也會通過大氣環(huán)流的調整而互相作用,從而對各地區(qū)氣溫產(chǎn)生影響。
圖5 合成的El Nio事件次年夏季的平均500 hPa位勢高度距平場(相對于的1981—2010年氣候平均場;單位:gpm) a.混合型;b.東部型;c.中部型Fig.5 Composite anomalous(relative to the 1981—2010 climate mean field) 500 hPa potential height of (a)EP+CP El Nio events (b)EP El Nio events and (c)CP El Nio events during the following summer(units:gpm)
圖5b顯示東部型,洋面上500 hPa位勢高度都有不同程度的增高,中國東部陸面上500 hPa位勢高度卻有不同程度的降低,中國東部地區(qū)經(jīng)向環(huán)流較混合型有所加強,利于冷空氣南下,氣溫較常年顯著偏低,此時西太平洋副熱帶高壓強度減弱,覆蓋范圍也減小了,華南也出現(xiàn)了一定程度的氣溫負距平。
圖5c顯示中部型,中國東部,經(jīng)向環(huán)流相對較強,利于冷空氣南下影響長江以北地區(qū),這種分布形態(tài)有利于北方氣溫的偏低,西太平洋副高位置偏南,形狀更加扁平,副高影響范圍增大且強度稍有增加,表明副髙對長江以南的影響加強,使得長江以南降溫幅度較小,而對長江及以北的影響有所減弱,從而促使氣溫下降。
3.3 低空和高空風場特征
高、低空風場的異常也會對氣溫產(chǎn)生影響,分析三類El Nio事件次年夏季高空、低空風場異常的合成特征,并結合位勢高度場的分析,探討風場對中國東部夏季氣溫的影響。
圖6 合成的混合型El Nio事件次年夏季風場距平場(相對于的1981—2010年氣候平均場;單位:m/s)a.850 hPa;b.200 hPaFig.6 Composite anomalous wind:(a)850 hPa;(b) 200 hPa(relative to the 1981—2010 climate mean field) of EP+CP El Nio events during the following summer(units:m·s-1)
圖7 合成的東部型El Nio事件次年夏季風場距平場(相對于的1981—2010年氣候平均場;單位:m/s)a.850hPa;b.200 hPaFig.7 Composite anomalous wind (a)850 hPa (b)200 hPa(relative to the 1981—2010 climate mean field) of EP El Nio events during the following summer(units:m·s-1)
圖8 合成的中部型El Nio事件次年夏季風場距平場(相對于的1981—2010年氣候平均場;單位:m/s)a.850 hPa;b.200 hPaFig.8 Composite anomalous wind:(a)850 hPa;(b)200 hPa(relative to the 1981—2010 climate mean field) of CP El Nio events during the following summer(units:m·s-1)
3)500 hPa上西太平洋副高的強度和范圍,是造成三類El Nio事件次年的夏季中國東部地區(qū)氣溫距平分布差異的主要原因之一?;旌闲虴l Nio事件次年夏季中國東部陸面上500 hPa位勢高度距平呈“三明治”型分布,東北東部—山東半島、華南地區(qū)為正異常,其余地區(qū)則有不同程度的降低,且負異常值幅度由北向南遞減;東部型,次年夏季中國東
部陸地上500 hPa位勢高度都有不同程度的降低,經(jīng)向環(huán)流較混合型有所加強,利于冷空氣南下,氣溫較常年顯著偏低,此時西太平洋副熱帶高壓強度減弱,覆蓋范圍減小,華南出現(xiàn)了一定程度的氣溫負距平。中部型,中國東部,經(jīng)向環(huán)流相對較強,利于冷空氣南下影響長江以北地區(qū),北方氣溫的偏低,西太平洋副高位置偏南,形狀更加扁平,副高影響范圍增大且強度稍有增加,副高對長江以南的影響加強,長江以南降溫幅度較小。
4)距平風場分析表明,850 hPa上西南暖濕氣流的輸送對El Nio次年的夏季中國東部氣溫距平空間分布有明顯影響。
致謝:NASA、ECMWF提供了AURA/MLS、CALIPSO/CALIOP、Interim資料的在線下載服務。
References)
Ashok K,Behera S K,Rao S A,et al.,2007.El Nio Modoki and its possible teleconnection[J].J Geophys Res,112,C11007.doi:10.1029/2006JC003798.
Chen D,Lian T,Fu C B,et al.,2015.Strong influence of westerly wind bursts on El Nio diversity[J].Nature Geoscience,8:339-345.doi:10.1038/ngeo2399.
陳麗華,周率,黨建濤,等,2010.2006年盛夏川渝地區(qū)高溫干旱氣候形成的物理機制研究[J].氣象,36(5):85-91. Chen L H,Zhou S,Dang J T,et al.,2010.Research on physical mechanism of summer severe climatic high temperature and drought disasters in Sichuan and Chongqing area in 2006[J].Meteor Mon,36(5):85-91.(in Chinese).
方宇凌,簡茂球,2011.2003年夏季華南持續(xù)高溫天氣過程及熱力診斷[J].熱帶海洋學報,30(3):30-37. Fang Y L,Jian M Q,2011.Diagnosis study of persistent heat waves in South China during summer 2003[J].J Trop Oceanography,30(3):30-37.(in Chinese).
黃丹青,錢永甫,2008.我國極端溫度事件的定義和趨勢分析[J].中山大學學報(自然科學版),47(3):112-116. Huang D Q,Qian Y F,2008.The definition of daily mean temperature extremes over china and its trend[J].Acta Scientiarum Naturalium Universitatis Sunyatseni,47(3):112-116.(in Chinese).
Kalnay E,Kanamitsu M,Kistler R,et al.,1996.The NCEP/NCAR 40-year reanalysis project[J].Bull Amer Meteor Soc,77:437-470.
雷楊娜,龔道溢,張自銀,等,2009.中國夏季高溫日數(shù)時空變化及其環(huán)流背景[J].地理研究,28(3):653-662. Lei Y N,Gong D Y,Zhang Z Y,et al.,2009.Spatial-temporal characteristics of high-temperature events in summer in eastern China and the associated atmospheric circulation[J].Geographical Research,28(3):653-662.(in Chinese).
李曉燕,翟盤茂,2000.ENSO事件指數(shù)與指標研究[J].氣象學報,58(1):102-109. Li X Y,Zhai P M,2000.On indices and indicators of ENSO episodes[J].Acta Meteorology Sinica,58(1):102-109.(in Chinese).
劉春玲,張強,許有鵬,等,2005.近半個世紀ENSO事件對長江三角洲地區(qū)氣候的影響[J].氣象,31(3):12-16. Liu C L,Zhang Q,Xu Y P,et al.,2005.Impacts of ENSO events on climatic changes during last 50 years in the Yangtze Delta[J].Meteor Mon,31(3):12-16.(in Chinese).
Rasmusson E M,Carpenter T H,1982.Variations in tropical sea surface temperature and surface wind fields associated with the Southern Oscillation/El Nio[J].Mon Wea Rev,110(5):354-384.
Rayner N,Parker D,Horton E,et al.,2003.Global analyses of sea sur-face temperature,sea ice,and night marine air temperature since thelate nineteenth century[J].J Geophys Res,108(D14):4407.
史軍,丁一匯,崔林麗,2009.華東極端高溫氣候特征及成因分析[J].大氣科學,33(2):347-358. Shi J,Ding Y H,Cui L L,2009.Climatic characteristics of extreme maximum temperature in East China and its causes[J].Chin J Atmos Sci,33(2):347-358.(in Chinese).
Weng H Y,Ashok K,Behera S K,et al.,2007.Impacts of recent El Nio Modoki on dry/wet conditions in the Pacific rim during boreal summer [J].Climate Dyn,29:113-129.doi:10.1007/s00382-007-0234-0.
徐金芳,鄧振錦,陳敏,2009.中國髙溫熱浪危害特征的研究綜述[J].干旱氣象,27(2):163-167. Xu J F,Deng Z J,Chen M,2009.Summary of studying on characteristics of high temperature and heat wave damage in China[J].Journal of Arid Meteorology,27(2):163-167.(in Chinese).
Yeh S W,Kug J S,Dewitte B,et al.,2009.El Nio in a changing climate[J].Nature,461(7263):511-514.
嚴文蓮,韓桂榮,魏建蘇,等,2012.2010年江蘇省連續(xù)性高溫特點及成因分析[J].氣象科學,32(3):293-297. Yan W L,Han G R,Wei J S,et al.,2012.Analysis on characteristic and cause of continuous high temperature of Jiangsu Province in 2010[J].J Meteor Sci,32(3):293-297.(in Chinese).
楊群,晏理華,周長志,等,2011.2009年銅仁地區(qū)高溫干旱特征及成因診斷分析[J].高原氣象,30(4):1018-1026. Yang Q,Yan L H,Zhou C Z,et al.,2011.Characteristic and diagnostic analysis of high temperature and drought in Tongren area in 2009[J].Plateau Meteor,30(4):1018-1026.(in Chinese).
楊義文,許力,龔振淞,2004.2003年北半球大氣環(huán)流及中國氣候異常特征[J].氣象,30(4):20-25. Yang Y W,Xu L,Gong Z S,2004.General circulation over the Northern Hemispherein 2003 and its impact on the climate in China[J].Meteor Mon,30(4):20-25.(in Chinese).
葉殿秀,尹繼福,陳正洪,等,2013.1961—2010年我國夏季高溫熱浪時空變化特征[J].氣候變化研究進展,9(1):15-20. Ye D X,Yin J F,Chen Z H,et al.,2013.Spatiotemporal change characteristics of summer heat waves in China in 1961—2010[J].Adv Climate Change Res,9(1):15-20.(in Chinese).
Zhang F,Gao H,Cui X,2008.Frequency of extreme high temperature days in China 1961—2003[J].Weather,63(2):46-49.
張茜,李棟梁,2011.東北及鄰近地區(qū)夏季氣溫異常的新特征及對大氣環(huán)流的響應[J].高原氣象,30(6):1604-1614. Zhang Q,Li D L,2011.New anomaly feature of summer air temperature in Northeast China and its responding to atmospheric general circulation[J].Plateau Meteor,30(6):1604-1614.(in Chinese).
張志薇,王式功,尚可政,等,2011.華中地區(qū)近年高溫事件及大氣環(huán)流成因分析[J].蘭州大學學報(自然科學版),47(2):50-55. Zhang Z W,Wang S G,Shang K Z,et al.,2011.Analysis of the atmospheric circulation anomalies of large-scale extreme high temperature events over central China in recent 50 years[J].Journal of Lanzhou University(Natural Sciences),47(2):50-55.(in Chinese).
周長艷,張順謙,齊冬梅,等,2013.近50年四川高溫變化特征及其影響[J].髙原氣象,32(6):1720-1728. Zhou C Y,Zhang S Q,Qi D M,et al.,2013.Variation features and impact of high temperature in Sichuan during last 50 years[J].Plateau Meteor,32(6):1720-1728.(in Chinese).
朱周平,錢炳強,2010.中國東部地區(qū)夏季高溫南北不均勻分布特征分析[J].氣象,36(11):26-31. Zhu Z P,Qian B Q,2010.Analysis of the meridional heterogeneous distribution of summer heat wave over Eastern China[J].Meteor Mon,36(11):26-31.(in Chinese).
鄒海波,吳珊珊,單九生,等,2015.2013年盛夏中國中東部高溫天氣的成因分析[J].氣象學報,73(3):481-495. Zhou H B,Wu S S,Shan J S,et al.,2015.Diagnostic study of the severe high temperature event over Mid-East China in 2013 summer[J].Acta Meteorologica Sinica,73(3):481-495.(in Chinese).
Study on the characteristics of air temperature in summer following El Nio events over East China
FAN Lingli1,ZHANG Fuying2
1GuangdongOceanUniversity,Zhanjiang524088,China;
2NanjingUniversityofInformationScience&Technology,Nanjing210044,China
Using the HadISST and NCEP/NCAR reanalysis data,a diagnosis and analysis was carried out for the abnormal temperature over eastern China(EC) during El Nio events in the following summer beginning from 1949.The results show that,in the summer after EP+CP El Nio,the negative temperature anomaly appeared in the Shandong-Huaihe Basin and along the coast of southern China,while the other areas showed positive anomaly.The EP El Nio,EC showed negative anomaly throughout,and the range in the north was greater than that to the south of the Yangtze River.For the CP El Nio,the temperature throughout the north was down more than that in the south of the Yangtze River.The corresponding atmospheric circulation analyses indicated that,in the summer after the three types of El Nio events,the anomalous zonal vertical circulation of the equatorial region(5°S—5°N) was shown to appear as three circles,but their anomolies not only had differences in strength,the latitude ranges were also different.In the summer after EP+CP El Nio events or EP El Nio events,there were two abnormal Walker circulations in the equatorial Pacific,while in the summer following CP El Nio events there were two Walker circulation anomalies in the equatorial India Ocean.The intensity and range of the subtropical high in the Western Pacific Ocean at 500 hPa was one of the main reasons for the difference of temperature distribution in EC during the summer after the three types of El Nio events.In the summer following EP+CP El Nio events,the 500 hPa geopotential height anomaly showed “sandwich” type distribution over eastern China.The eastern part of northeast China,the Shandong Peninsula and southern China area were positive anomalies,while the other regions were reduced to some extent,and the negative anomaly amplitude decreased from north to south.In the summer after EP El Nio events,the 500 hPa geopotential height had different degrees of reduction over EC,the meridional circulation was strengthened,which was conducive to the cold air travelling southward,thus the temperature was significantly lower than normal,the intensity of the subtropical high in the Western Pacific was weakened,its coverage decreased,and there was a certain degree of negative temperature anomaly in southern China.In the summer after CP El Nio events,the meridional circulation was relatively strong at the 500 hPa geopotential height over EC,which was conducive to the cold air travelling southward to the north of the Yangtze River region,the temperature was low in the north,the West Pacific subtropical high travelled southward,with a rather flat shape,the effect range of the subtropical high increased and its strength increased slightly,and the effect of the subtropical high on southern China strengthened,so that the temperature in southern China increased.However,the effect of the subtropical high to the north of the Yangtze River weakened and caused the temperature there to drop;the transportation of the southwest warm and humid air at 850 hPa also had a significant effect on this.
(責任編輯:劉菲)
范伶俐,張福穎,2016.El Nio事件次年中國東部夏季氣溫分布特征研究[J].大氣科學學報,39(6):845-854. Fan L L,Zhang F Y,2016.Study on the characteristics of air temperature in summer following El Nio events over East China[J].Trans Atmos Sci,39(6):845-854.
10.13878/j.cnki.dqkxxb.20130427001.(in Chinese).
10.13878/j.cnki.dqkxxb.20160815003
*聯(lián)系人,E-mail:fanlingli@126.com