熱帶季節(jié)內(nèi)振蕩的年際變化與增暖的耦合

熱帶季節(jié)內(nèi)振蕩的年際變化與增暖的耦合

1氣候變化對帶頭氣候的影響由于mad-ulan的隨機(jī)性，它通常也被稱為mad-ulan風(fēng)場和壓縮場，以便首先在熱帶大氣中發(fā)現(xiàn)風(fēng)場和壓縮場。ISO是熱帶大氣在7~90天時間尺度的最強(qiáng)振蕩信號,有明顯的區(qū)域分布特征,熱帶印度洋和熱帶西太平洋最強(qiáng);時域上,ISO強(qiáng)度有明顯的季節(jié)變化,冬春季強(qiáng),夏季最弱。由于其時間尺度介于月、季之間,因而與長期天氣變化和短期氣候異常均有密切關(guān)聯(lián),一直受到氣象學(xué)界的廣泛重視。一些研究還表明,季節(jié)內(nèi)振蕩既有明顯的季節(jié)變化也存在明顯的年際變化,許多研究將ISO的年際變化與ENSO現(xiàn)象建立了聯(lián)系,最近有研究還注意到了熱帶ISO的年代際變化。全球氣候正在并且將繼續(xù)變暖已逐漸成為世界上大多數(shù)科學(xué)家的共識,令人感興趣的問題是這種氣候變暖趨勢會對海洋和大氣環(huán)流系統(tǒng)產(chǎn)生什么影響?許多研究發(fā)現(xiàn),20世紀(jì)70年代中期氣候系統(tǒng)出現(xiàn)了一次躍變(這與全球地面氣溫的再度增暖是對應(yīng)的),導(dǎo)致ENSO事件演變、青藏高原熱量源匯及夏季風(fēng)活動均產(chǎn)生了大的變化。Slingo等在研究熱帶ISO與熱帶海洋SST之間的關(guān)聯(lián)時也發(fā)現(xiàn),在年際變化上二者顯示了一個弱的并不顯著的關(guān)系,但在年代際變化上二者似乎存在更好的關(guān)聯(lián)。這給予我們一個提示,是否熱帶ISO的較長尺度的變化與增暖有關(guān)?考慮到中國科學(xué)院大氣物理研究所發(fā)展的耦合模式已經(jīng)成功地用于氣候研究的許多方面,在本文我們將利用該耦合模式進(jìn)行數(shù)值模擬并從以下二方面進(jìn)行研究:一是自然變率(無人類活動影響)下熱帶ISO的變化特征;二是氣候增暖背景(有人類活動影響)下熱帶ISO的可能變化。這是因為任何實際的氣候變化都是由氣候自然變化與人類活動、其他外部強(qiáng)迫因子引起的變化疊加而成的,而目前普遍認(rèn)為,過去100年的全球增暖可能主要是由人類活動向大氣排放溫室氣體造成的,理論研究和模式模擬也證實了大氣溫室氣體增加能引起全球平均地面氣溫升高。通過對上述兩個方面的對比分析研究,就可以區(qū)分出熱帶ISO的哪些變化是由人類活動(溫室氣體增加)引起?哪些變化可能是自然的變化?由于只有用模式才能產(chǎn)生自然變率,因此,數(shù)值模擬研究是開展上述兩方面研究的唯一手段。Slingo等利用大氣模式比較計劃(AMIP)中15個AGCM模擬ISO的結(jié)果顯示,大氣環(huán)流模式模擬ISO普遍存在以下不足:低估了ISO振蕩強(qiáng)度,未能體現(xiàn)出ISO的季節(jié)性特點,模擬的季節(jié)內(nèi)振蕩周期偏短等。根據(jù)簡單海-氣耦合模式能較好地模擬ISO這一特點,李薇等利用完全耦合模式對ISO模擬效果開展了一些敏感性試驗和評估,研究結(jié)果表明,利用完全耦合模式再現(xiàn)ISO的季節(jié)性特征更真實,反映了海氣耦合對ISO的調(diào)制作用,這也表明耦合模式較單純大氣環(huán)流模式優(yōu)越,再考慮到研究的時間尺度較長,因此,采用耦合氣候模式進(jìn)行熱帶大氣季節(jié)內(nèi)振蕩的模擬將是較好的選擇。本文的研究目的是利用FGOALS1.0_g耦合模式控制試驗?zāi)M結(jié)果探討在工業(yè)革命以前(忽略人類活動影響)熱帶季節(jié)內(nèi)振蕩年際以上尺度的變化特征,也就是在自然變化下熱帶季節(jié)內(nèi)振蕩年際以上時間尺度的變化特征,一方面了解無人類活動影響情形下熱帶季節(jié)內(nèi)振蕩在較長時間尺度上的變化特性,另一方面也為進(jìn)一步研究人類活動對季節(jié)內(nèi)振蕩的影響(將另文討論)奠定基礎(chǔ)。2模型總結(jié)、數(shù)據(jù)和方法2.1海洋分量模式本文中所使用的耦合氣候系統(tǒng)模式是中國科學(xué)院大氣物理研究所大氣科學(xué)和地球流體力學(xué)數(shù)值模擬國家重點實驗室(LASG)最新發(fā)展的第四代耦合氣候系統(tǒng)模式FGOALS1.0_g(FlexibleGlobalOcean-Atmosphere-LandSystemModel,Version1.0_g),該耦合模式利用一個通量耦合器,將海洋、大氣、陸面和海冰四個分量模式耦合在一起,在模式界面上考慮了熱量、動量和淡水通量的耦合,并且沒有使用任何形式的通量訂正技術(shù)。其中的海洋分量模式是在LICOM1.0的基礎(chǔ)上改動而來的,主要變化將模式的水平范圍從65°N擴(kuò)展到北極,同時因受計算機(jī)能力的限制,還將模式的水平分辨率由0.5°×0.5°降為1°×1°。改變以后的海洋模式計算速度大為提高,而模式的基本性能與原來差別不大。關(guān)于該耦合模式的進(jìn)一步說明可以參見Yu等的研究工作。該耦合模式的大氣分量模式是GAMIL1.0,該模式是Wang等將NCAR的大氣環(huán)流模式CAM2中的物理過程引進(jìn)到一個新發(fā)展的格點大氣動力框架中發(fā)展而成的。GAMIL1.0的水平分辨率大約為2.8°×2.8°,模式可以比較成功地模擬大氣的基本型式。FGOALS1.0_g中的海冰分量模式是一個熱力-動力學(xué)海冰模式。該海冰模式在垂直方向有5層,其動力學(xué)過程是基于彈性-塑性流變學(xué)原理,而熱力過程則采用了一種能量守恒的處理辦法。鑒于耦合器要求海洋和海冰模式的格點必須完全一致,在本文中作者將海冰模式的水平分辨率提高到與海洋環(huán)流模式LICOM完全一致的1°×1°的經(jīng)緯網(wǎng)格,而不是原先的將北極點旋轉(zhuǎn)到大陸上的正交曲線坐標(biāo),同時海冰模式也采用LICOM完全一樣的海陸分布。關(guān)于該耦合模式的具體描述可以參見Yu等的研究工作。利用該耦合模式,我們完成了多個長期數(shù)值模擬試驗,其中包括工業(yè)革命前控制試驗以及大氣中二氧化碳濃度從工業(yè)革命前的280×10-6以1%等比增加到560×10-6的加倍試驗。本文的第I部分主要利用工業(yè)革命前控制試驗(簡稱控制試驗,即模擬的自然變率)模擬結(jié)果研究模擬的自然變率下大氣ISO的一些變化特征;第II部分則利用對比控制試驗和加倍試驗?zāi)M的ISO現(xiàn)象,探討全球變暖對ISO的可能影響。2.2與實測資料的對比本文所用資料主要有二類:一類資料是作為實測資料的56年NCEP/NCAR逐日再分析資料,資料的時段為1948~2003年,該類資料主要用于對模式模擬結(jié)果評估對比時使用。另一類資料就是利用上述耦合氣候模式模擬的模式結(jié)果資料,我們選取了與二氧化碳濃度從工業(yè)革命前的280×10-6以1%等比增加到560×10-6的加倍試驗(共70年)所對應(yīng)的控制試驗中60年逐日模擬結(jié)果資料,包括高層(200hPa)和低層(850hPa)經(jīng)、緯向風(fēng),下墊面海表水溫、降水等要素資料,取60年主要是與實測資料的56年NCEP/NCAR資料近似等長,以方便對比。本文所使用的診斷方法主要有ISO指標(biāo)場方法、時空譜方法、小波分析方法等。3對比各階段的實際應(yīng)用結(jié)果—模擬的自然變率下大氣季節(jié)內(nèi)振蕩的基本特征在進(jìn)行自然變率下熱帶ISO年際與年代際變化特征的研究之前,首先分析用FGOALS1.0_g耦合模式控制試驗?zāi)M結(jié)果模擬的季節(jié)內(nèi)振蕩的一些基本特征,如活躍區(qū)位置、季節(jié)變化及ISO周期等。Slingo等在利用大氣模式比較計劃(AMIP)中15個模式檢驗AGCM模擬ISO的能力后指出,對氣候基本態(tài)的準(zhǔn)確描述是模式能夠刻畫出實際的大氣季節(jié)內(nèi)振蕩的先決條件。我們將FGOALS1.0_g耦合模式模擬的氣候基本態(tài)與實際的氣候基本態(tài)(由NCEP/NCAR資料計算)進(jìn)行了比較(圖略),發(fā)現(xiàn)FGOALS1.0_g耦合模式能很好再現(xiàn)實際的氣候基本態(tài),這為FGOALS1.0_g耦合模式能較好地模擬季節(jié)內(nèi)振蕩奠定了-定基礎(chǔ)。為較好地描述ISO的區(qū)域特征、季節(jié)變化以及年際、年代際變化特征,本文通過創(chuàng)建一個ISO強(qiáng)度指標(biāo)場來實現(xiàn),并借助這個ISO強(qiáng)度指標(biāo)場,通過對比模擬結(jié)果與實際結(jié)果(NCEP/NCAR資料得到)來檢驗在季節(jié)內(nèi)尺度上模式的性能。ISO強(qiáng)度指標(biāo)場的計算過程大致如下:首先,對選定的區(qū)域(既可以是全球也可以是整個熱帶區(qū)域)的每一點逐日時間序列進(jìn)行20~100天帶通濾波;然后,利用滑動窗口技術(shù),取滑動步長為1個日歷月份,并以100天為窗口寬度來計算其方差以代表該日歷月份的方差即ISO強(qiáng)度;最后,為消除局地方差的影響,同時計算該點的總方差(用濾波前的時間序列),并取季節(jié)內(nèi)方差與總方差的比值代表所選點的ISO強(qiáng)度。利用上述方法得到的是一個隨時間、空間變化的場,由于它描述的是ISO強(qiáng)度(相對值),故稱之為ISO指標(biāo)場。圖1給出了利用NCEP/NCAR再分析資料由上述指標(biāo)場(要素取200hPa速度勢)計算的1月、4月、7月、10月(分別代表冬、春、夏、秋四季)1949~2002年多年平均的全球ISO強(qiáng)度空間分布。此圖代表了全球ISO強(qiáng)度季節(jié)變化的平均特征。因20~100天濾波因素,丟掉了兩頭部分信息,為保持完整性,我們采用去掉頭尾兩年不完整結(jié)果(以下處理類似)。由圖1可以看出,ISO指標(biāo)較好地描述了季節(jié)內(nèi)振蕩的空間分布和季節(jié)變化特征:ISO活躍區(qū)(強(qiáng)度高值區(qū))主要位于熱帶印度洋和熱帶西太平洋;有明顯季節(jié)變化,冬春季強(qiáng),夏秋季弱;冬季,ISO活躍區(qū)位置偏南,夏季,ISO活躍區(qū)位置偏北。圖2是由FGOALS1.0_g耦合模式控制試驗?zāi)M結(jié)果分析得到的不同季節(jié)200hPa速度勢的全球ISO強(qiáng)度空間分布圖(58年平均)。由圖2可以發(fā)現(xiàn),首先,模擬的自然變率下的ISO活躍區(qū)主要位于印度洋與西太平洋,這與近六十年來的觀測分析結(jié)果基本接近,雖然自然變率下的分析結(jié)果由于不包括人類活動影響,因而與實測結(jié)果不完全具有可比性,但此結(jié)果表明,大氣ISO活躍區(qū)的位置主要位于印度洋與西太平洋是其基本特征之一;其次,模擬的自然變率下的ISO冬季主要活躍區(qū)的位置偏南,夏季偏北的特征也很明顯,這與觀測分析結(jié)果的變動趨勢完全一致;第三,模擬的自然變率下的ISO冬春強(qiáng)、夏季弱的季節(jié)變化也較明顯;這與實測分析結(jié)果一致。但也可以看出,模擬的ISO強(qiáng)度要比觀測結(jié)果弱,這可能有兩方面的原因,其一可能是模式本身的問題,因為模擬的ISO強(qiáng)度偏弱是單純大氣環(huán)流模式的一個較為普遍的現(xiàn)象,耦合模式可能同樣不能完全避免,這與大氣模式中所取的對流參數(shù)化方案有一定關(guān)系;其二可能就是自然變率狀況下ISO的特征,關(guān)于這一問題可以通過對比考慮人類活動影響的加倍試驗(模擬結(jié)果將在第II部分進(jìn)行討論)。圖3進(jìn)一步給出了用NCEP/NCAR再分析資料以及海氣耦合控制試驗?zāi)M結(jié)果資料對赤道海區(qū)(10°S~10°N)200hPa速度勢進(jìn)行時空譜分析得到的結(jié)果,時空譜計算是對每一年資料逐年計算,最后計算多年平均。其中,圖3a為NCEP/NCAR再分析資料計算結(jié)果(56年平均),圖3b為海氣耦合控制試驗?zāi)M結(jié)果(60年平均)。與NCEP/NCAR結(jié)果(圖3a)相比,耦合模式模擬的ISO周期(圖3b)明顯偏弱,且僅在高頻段有一小峰值,表明模擬的自然變率下的ISO周期偏弱。這同樣可能有兩方面原因,其一是數(shù)值模式模擬ISO結(jié)果的一個較為普遍現(xiàn)象,Slingo等在AMIP比較計劃中所使用的15個模式大多都有這樣的特點,FGOALS1.0_g耦合模式同樣也未能避免;其二可能與自然變率本身有關(guān)。是哪方面的原因?qū)⒃诒狙芯康牡贗I部分作進(jìn)一步的探討。4模擬的自然變率帶為了探討熱帶ISO強(qiáng)度的年際與年代際變化特征,我們定義一個用以描述熱帶ISO強(qiáng)度隨時間變化的指標(biāo)。對上面獲得的ISO強(qiáng)度指標(biāo)場,我們?nèi)≌麄€熱帶區(qū)域(30°S~30°N,0°~360°E),對該區(qū)域的ISO強(qiáng)度指標(biāo)場進(jìn)行EOF展開,將所獲得的空間第一特征向量(占總方差的26.54%)對應(yīng)的時間系數(shù)作為熱帶ISO強(qiáng)度指數(shù),由于熱帶ISO強(qiáng)度的季節(jié)變化明顯,計算前進(jìn)行了消除季節(jié)變化的處理以便突出年際與年代際變化?？臻g第一特征向量的空間分布顯示出(圖略),熱帶印度洋與熱帶西太平洋是ISO強(qiáng)度高值區(qū),表明空間第一特征向量代表了熱帶ISO的主要空間分布特征。圖4給出了對應(yīng)空間第一特征向量的時間系數(shù),即熱帶ISO強(qiáng)度指數(shù)的年際與年代際變化曲線,圖中的年際與年代際變化曲線系利用小波的濾波特性由小波分解重構(gòu)方法計算得出,此圖反映了模擬的自然變率下熱帶ISO強(qiáng)度的年際與年代際變化。由圖4可以看出,模擬的自然變率下的熱帶ISO既存在明顯的年際變化也存在年代際變化。對年際變化的周期特征,我們用Morlet小波進(jìn)行了時頻分析,結(jié)果發(fā)現(xiàn)(圖略),模擬的自然變率下熱帶ISO強(qiáng)度存在2~10年之間的準(zhǔn)周期變化,其中,前40年10年周期較明顯,但在有些時段也存在2~3年左右的周期變化;后20年以5~7年左右準(zhǔn)周期為主。圖4b表明,模擬的58年自然變率下熱帶ISO強(qiáng)度的年代際變化也較明顯,就所研究的時段而言,總體變化趨勢是下降的。最初的15年熱帶ISO強(qiáng)度較強(qiáng),然后逐漸減弱,在35~53年這個時段最弱,之后又有所增強(qiáng)。為了進(jìn)一步說明年代際變化的差異,我們選取了兩個時段(如圖4b),選取第2年至第19年這18年代表熱帶ISO強(qiáng)度高指數(shù)時段,而35~52年這個時段代表熱帶ISO強(qiáng)度低指數(shù)時段。分別計算了這兩個時段模擬的全球ISO強(qiáng)度在不同季節(jié)的空間分布圖。在熱帶ISO強(qiáng)度較強(qiáng)時段(圖5),與平均情況相比(圖2),ISO強(qiáng)度與范圍有增強(qiáng)與擴(kuò)大趨勢,但主要是春季和秋季變化稍明顯,冬季并不明顯。在熱帶ISO強(qiáng)度較弱時段(圖6),ISO強(qiáng)度與范圍則明顯減弱與縮小,但同樣也存在季節(jié)性差異,冬季變化不大。如果就高低強(qiáng)度指數(shù)時段之間進(jìn)行對比,自然變率下的熱帶ISO強(qiáng)度的年代際變化是相當(dāng)明顯的,但是存在季節(jié)性差異,冬季變化不大,春秋季最明顯。給出了用NCEP/NCAR資料計算得出的偏冷階段(圖7)和偏暖階段(圖8)不同季節(jié)多年平均全球大氣ISO強(qiáng)度的空間分布,階段的劃分是依據(jù)全球地面平均氣溫變化特征(圖略)確定的,1948~1976為偏冷階段,1977~2003為偏暖階段。通過圖7、8的對比可以發(fā)現(xiàn),偏暖階段大氣ISO無論是強(qiáng)度還是活躍區(qū)范圍均比偏冷階段大,但存在季節(jié)性差異,夏秋季十分明顯,冬季不明顯。實測結(jié)果在兩個不同階段的變化特征與模擬的自然變率下的高/低指數(shù)階段的變化特征基本一致,但實測分析結(jié)果似乎表明增暖對應(yīng)高指數(shù)階段,偏冷對應(yīng)低指數(shù)階段,即增暖可能是ISO強(qiáng)度增強(qiáng)、范圍擴(kuò)大的原因。5自然變率下帶能量的年際變化研究熱帶ISO的傳播特征,通常都是采用時空譜方法,本文也采用此方法。另外,為了方便研究年際與年代際變化,時空譜計算中也采用了滑動窗口技術(shù),并取窗口寬度為365天(1年),滑動步長為1個日歷月份,這樣得到的是逐月的時空譜值。圖9給出了用FGOALS/IAP海氣耦合模式控制試驗?zāi)M的赤道地區(qū)(10°S~10°N)大氣200hPa速度勢東傳(圖9a)和西傳(圖9b)1~6波ISO時空功率譜比值(1~6波ISO時空功率譜之和與ISO時空功率譜總和之比)的年代際變化曲線,用以描述自然變率下熱帶ISO主要的東傳和西傳波能量隨時間的變化。其中,1~6波ISO時空功率譜比值是前6個東傳(或西傳)波時空功率譜之和與季節(jié)內(nèi)振蕩總時空功率譜之比,采用比值的目的是比較相對大小。由圖9可以看出,東傳1~6波季節(jié)內(nèi)振蕩能量在總季節(jié)內(nèi)振蕩能量中所占比例平均約為60%,比西傳(38%左右)明顯高,這與觀測中東傳占優(yōu)的實際情況相似,前6個波平均季節(jié)內(nèi)振蕩能量巳占季節(jié)內(nèi)振蕩總能量的98%(東傳與西傳之和),表明熱帶ISO以前6個波為主,這也與實際觀測相似。由圖9還可以發(fā)現(xiàn),無論是東傳波還是西傳波,雖然1~6波ISO時空功率譜比值存在明顯的年代際變化,但這種變化不存在明顯的上升或下降趨勢,它僅圍繞著平均值波動,用低層850hPa速度勢所得結(jié)果相同(圖略)。這似乎表明,自然變率下的熱帶大氣東傳或西傳的ISO能量雖然也存在年際與年代際變化,但從長期來看基本維持一種平衡狀態(tài)。作為對比,我們給出的近六十年來實際狀況下的對應(yīng)結(jié)果,圖10是用1948~2003年共56年NCEP/NCAR逐日再分析資料分析得到的近六十年來赤道地區(qū)(10°S~10°N)緯向東傳或西傳1~6波之和季節(jié)內(nèi)振蕩能量相對比值的年代際變化曲線。由圖1