• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    中國冬季氣溫不同年代際的季節(jié)內(nèi)變化特征及成因分析

    2022-06-01 07:15:46馬鋒敏陳麗娟李想李維京許彬
    大氣科學(xué) 2022年3期
    關(guān)鍵詞:置信水平環(huán)流氣溫

    馬鋒敏 陳麗娟 李想 李維京 許彬

    1 江西省氣候中心, 南昌330096

    2 國家氣候中心中國氣象局氣候研究開放實(shí)驗(yàn)室/中國氣象局—南京大學(xué)氣候預(yù)測研究聯(lián)合實(shí)驗(yàn)室, 北京100081

    3 南京信息工程大學(xué)氣象災(zāi)害預(yù)報(bào)預(yù)警與評估協(xié)同創(chuàng)新中心, 南京210044

    1 引言

    氣候變化是長期大氣狀態(tài)變化的一種反映,而氣溫是氣候變化研究的重要基礎(chǔ)性問題(Overpeck et al., 2011; 任國玉等, 2014)。近百年來地球氣候正經(jīng)歷一次以變暖為主要特征的顯著變化,其中冬季變暖最為明顯(陳隆勛等, 1991; 王紹武等, 1998;梁蘇潔等, 2014)。我國處于東亞季風(fēng)區(qū),東亞冬季風(fēng)異常會(huì)導(dǎo)致我國冬季出現(xiàn)大范圍冷暖異常。近年來的研究顯示,不僅東亞冬季風(fēng)具有較大的年際變率和年代際變化(Wang, 2001; 左瑞亭等, 2005;施曉暉等, 2007; Wang et al., 2009; Wang and He,2012; Wang and Chen, 2014; Li et al., 2019),影響我國冬季氣溫的重要影響因子,如西伯利亞高壓、北極濤動(dòng)、北極海冰和海表溫度等也具有年代際變化(朱亞芬和錢維宏, 1998; Watanabe and Nitta,1999; 武炳義等, 2000; Kug et al., 2009; Jeong et al.,2011)。此外,影響因子之間的關(guān)系同樣發(fā)生了顯著的年代際變化(李維京等, 2013)。這種新的變化,增加了對冬季氣候預(yù)測的難度。因此,研究不同氣候背景下冬季氣溫異常成為氣候變化和氣候預(yù)測研究的重點(diǎn)和難點(diǎn)。在全球氣候變暖的背景下,我國冬季氣溫異常呈現(xiàn)出一些新的特征。季節(jié)內(nèi)極端天氣氣候事件頻發(fā),強(qiáng)寒潮、雨雪冰凍等災(zāi)害性天氣對國民生產(chǎn)和人們生活造成很大影響,社會(huì)各級部門和公眾對冬季氣候預(yù)測精細(xì)化服務(wù)的需求增加。如何科學(xué)認(rèn)識氣候異常的特征和演變規(guī)律,準(zhǔn)確預(yù)測冬季氣候異常程度已經(jīng)成為迫切需要解決的科學(xué)問題,對提高防災(zāi)減災(zāi)服務(wù)能力有重要的應(yīng)用價(jià)值。

    針對我國冬季季節(jié)平均氣溫主要模態(tài)的研究取得了許多有意義的成果(Wu and Wang, 2002; Jhun and Lee, 2004; 康麗華等, 2006, 2009; 朱艷峰等,2007; Wang et al., 2003, 2010; 范可和劉輝, 2013; 孫淑清等, 2014; 汪子琪等, 2017),多數(shù)研究認(rèn)為中國冬季氣溫主要模態(tài)是一致偏暖或偏冷。然而基于冬季平均氣溫進(jìn)行氣候預(yù)測及服務(wù)存在較大的局限性,尤其是近年來我國冬季經(jīng)常出現(xiàn)季節(jié)內(nèi)氣溫變化大,前冬和后冬氣候異常相反的觀測事實(shí)。例如,2007 年12 月,全國氣溫異常偏高,而2008 年1~2 月出現(xiàn)大范圍低溫雨雪冰凍災(zāi)害(王凌等,2008; 高輝等, 2008; 張慶云等, 2008; Zhou et al.,2009, 2011)。2020 年12 月,全國大部地區(qū)氣溫異常偏低,而2021 年1 月中旬以后,氣溫大幅度回升,2021 年2 月氣溫異常偏高。許多研究關(guān)注到我國冬季氣溫具有明顯的季節(jié)內(nèi)變化特征(Hu et al., 2005; 黃嘉佑和胡永云, 2006; 韋瑋等, 2014,2020 ; Zuo et al., 2015; 賈丹和簡茂球, 2015; 孫健等,2019),尤其是前冬和后冬氣溫的變化可能存在很大的差異(黃嘉佑和胡永云, 2006)。韋瑋等(2014)指出,我國前冬(11~12 月)和后冬(1~3 月)氣溫的變化具有前、后冬趨勢一致的同相演變和前、后冬趨勢相反的反相演變兩個(gè)主要模態(tài),且中高緯度地區(qū)大氣遙相關(guān)型位相的持續(xù)或轉(zhuǎn)換是造成不同演變模態(tài)的重要原因。孫健等(2019)指出中國冬季氣溫在月尺度上常出現(xiàn)前、后冬相反甚至冷暖交替的現(xiàn)象。劉艷等(2020)研究表明北疆冬季氣溫季節(jié)內(nèi)主要有前后距平一致和前后距平相反兩種特征。孟鑫等(2017)指出東北地區(qū)前冬(12 月)和后冬(1~2 月)氣溫在空間上具有全區(qū)變化一致型和東北—西南反相變化的蹺蹺板型。以上研究揭示了過去幾十年來我國冬季季節(jié)內(nèi)變化特征,沒有區(qū)分不同的年代際背景。而已有研究顯示在不同的年代際時(shí)間尺度下,冬季氣溫和影響因子的關(guān)系可能發(fā)生了變化,同時(shí)季節(jié)內(nèi)氣候異常特征也可能發(fā)生了變化,因此重新認(rèn)識不同的年代際時(shí)間尺度下冬季季節(jié)內(nèi)氣溫變化的新規(guī)律及其環(huán)流特征,對理解冬季氣候異常特征和季節(jié)內(nèi)變化的復(fù)雜性以及提高短期氣候預(yù)測準(zhǔn)確率都有重要價(jià)值。

    2 資料和方法

    文中所用的資料包括中國氣象局國家氣候中心整編的全國160 站1951~2020 年逐月氣溫。大氣環(huán)流資料為美國國家環(huán)境預(yù)報(bào)中心和國家大氣研究中心(NCEP/NCAR)逐月再分析資料中的近地面溫度場、海平面氣壓場、500 hPa 位勢高度場、200 hPa 緯向風(fēng)場,資料水平分辨率為2.5°×2.5°(Kalnay et al., 1996; Kistler et al., 2001)。本文針對1951~2019 年冬季進(jìn)行研究,文中冬季指當(dāng)年12 月、次年1 月、2 月的平均。

    本文使用的主要統(tǒng)計(jì)方法有:相關(guān)分析、Mann-Kendall 突變法、滑動(dòng)相關(guān)、線性回歸、合成分析、Student-t檢驗(yàn)、小波分析等(魏鳳英, 2007)以及基于季節(jié)的經(jīng)驗(yàn)正交函數(shù)分解(S-EOF)(Wang and An, 2005)。S-EOF 是Wang and An(2005)在研究太平洋海溫氣候變率時(shí)對經(jīng)驗(yàn)正交函數(shù)分解(EOF)進(jìn)行改進(jìn)發(fā)展而成的一種新方法,能夠很好地識別隨季節(jié)演變的變量場在時(shí)間上連續(xù)的空間特征和演變規(guī)律,從而得到隨季節(jié)演變的主導(dǎo)模態(tài)。Wu et al.(2009)、孫健等(2019)、韋瑋等(2014, 2020)利用S-EOF 開展了相關(guān)研究工作。本文利用S-EOF 方法分析不同年代際時(shí)間尺度下中國冬季氣溫季節(jié)內(nèi)變化主導(dǎo)模態(tài)和標(biāo)準(zhǔn)化時(shí)間序列,進(jìn)一步采用合成和回歸的方法分析主導(dǎo)模態(tài)對應(yīng)的大氣環(huán)流場異常特征。另外,采用滑動(dòng)相關(guān)(林學(xué)椿, 1978)的方法確定中國冬季各月氣溫之間關(guān)系的穩(wěn)定性,滑動(dòng)窗口長度取為21 a,滑動(dòng)相關(guān)值記在窗口的第11 年,例如1951~1971 年相關(guān)值記在1961 年。

    3 冬季氣溫的季節(jié)內(nèi)變化特征

    3.1 前冬與后冬的確定

    利用1951~2020 年月平均氣溫資料,計(jì)算全國160 站平均氣溫在冬季不同月份間的相關(guān)系數(shù)(表1)??梢钥闯?,無論是否去除線性趨勢,12月與1 月、12 月與2 月的氣溫相關(guān)系數(shù)均不顯著,而1 月與2 月氣溫的聯(lián)系十分緊密,去除線性趨勢前后的相關(guān)系數(shù)分別為0.49 和0.34,均超過了99%的置信水平。計(jì)算去趨勢后的1951~2019 年冬季12 月與1 月(圖1a)、1 月與2 月(圖1b)各站點(diǎn)的相關(guān)系數(shù),空間分布圖也顯示我國冬季存在明顯的季節(jié)內(nèi)差異。12 月與1 月氣溫僅在東北部分地區(qū)、西南地區(qū)存在顯著相關(guān);而1 月與2 月的氣溫相關(guān)達(dá)到顯著性檢驗(yàn)的范圍明顯擴(kuò)大,除我國東南沿海相關(guān)不顯著之外,全國其他地區(qū)相關(guān)系數(shù)均通過顯著性檢驗(yàn),尤其是長江以北的大部地區(qū)相關(guān)系數(shù)均通過99%的置信水平。

    圖1 1951~2019 年冬季(a)12 月與1 月、(b)1 月與2 月氣溫(去趨勢后)相關(guān)系數(shù)空間分布。深、淺陰影分別表示通過99%和95%的置信水平Fig. 1 Spatial distributions of the correlation coefficients between (a) December and January, (b) January and February with detrended air surface temperature in winter over China from 1951 to 2019. Dark and light shadings indicate the 99% and 95% confidence levels, respectively

    表1 1951~2019 年冬季不同月份間全國160 站平均氣溫的相關(guān)系數(shù)Table 1 Correlation coefficients of surface air temperature averaged in 160 stations over China between different months from December to February for 1951-2019

    前文分析的是整個(gè)研究時(shí)段內(nèi)冬季各月氣溫間的關(guān)系,為了認(rèn)識中國冬季季節(jié)內(nèi)各月氣溫之間的年際關(guān)系是否穩(wěn)定,進(jìn)一步計(jì)算冬季各月間氣溫的21 a 滑動(dòng)相關(guān)(圖2)。1 月與2 月的滑動(dòng)相關(guān)顯示,二者在20 世紀(jì)70 年代以后的相關(guān)達(dá)到95%的置信水平,此后的相關(guān)基本穩(wěn)定,且通過95%的置信水平。而12 月與1 月的滑動(dòng)相關(guān)沒有通過信度檢驗(yàn),僅12 月與2 月在1990 年前后達(dá)到95%的置信水平,但是其他時(shí)段的相關(guān)均不顯著。這表明在年代際尺度上,1 月和2 月的關(guān)系更密切。再次計(jì)算前冬(12 月)與后冬(1~2 月)的21 a 滑動(dòng)相關(guān)(圖2 藍(lán)線),二者僅在1990 年前后達(dá)到95%的置信水平,其他時(shí)段相關(guān)均較弱。

    圖2 1951~2019 年冬季不同月份間中國160 站平均氣溫21 a 滑動(dòng)相關(guān)。黑色虛線:12 月與1 月;綠色實(shí)線:12 月與2 月;紅色實(shí)線:1月與2 月;藍(lán)色實(shí)線:前冬與后冬;灰色點(diǎn)線:95%的置信水平臨界線Fig. 2 Sliding correlations of surface air temperature averaged in 160 stations over China between different months from December to February with a 21-year moving window for 1951-2019. The black dashed line: between December and January; green solid line: between December and February;red solid line: between January and February; blue solid line: between early winter (EW) and late winter (LW); the gray straight dotted lines represent 95% confidence levels

    為了考查前冬和后冬氣溫空間相關(guān)關(guān)系隨時(shí)間的變化,進(jìn)一步計(jì)算了1951~2019 年69 個(gè)冬季中國160 站前冬與后冬氣溫21 a 滑動(dòng)相關(guān)的空間分布,滑動(dòng)相關(guān)值記在窗口的第11 年,得出1961~2009 年時(shí)間序列的逐年空間分布,圖3 僅給出每間隔10 a 的21 a 滑動(dòng)相關(guān)分布??梢钥闯?,全國大部分地區(qū)的相關(guān)很弱。相關(guān)顯著的區(qū)域有20 世紀(jì)70 年代的內(nèi)蒙古中部和東北部、西南地區(qū)南部;20 世紀(jì)80 年代的華北東北部;20 世紀(jì)90年代的東南沿海和東北南部;21 世紀(jì)的東北中北部地區(qū)。其他時(shí)段的顯著相關(guān)區(qū)域范圍很小(圖略)。

    圖3 1951~2019 年中國前冬與后冬氣溫21 a 滑動(dòng)相關(guān)空間分布?;瑒?dòng)窗口時(shí)間分別為(a)1970、(b)1980、(c)1990、(d)2000。深、淺陰影分別表示通過99%和95%的置信水平Fig. 3 Spatial distributions of sliding correlation coefficients between early and late winter temperature over China with a 21-year moving window.The years of moving window are (a) 1970, (b) 1980, (c) 1990, (d) 2000, respectively. Dark and light shadings indicate the 99% and 95% confidence levels, respectively

    上述對前冬、后冬的時(shí)間和空間相關(guān)分析表明,1 月與2 月的顯著正相關(guān)特征總體是比較穩(wěn)定的,而前冬(12 月)與后冬(1~2 月)的相關(guān)僅在20世紀(jì)90 年代前后相關(guān)系數(shù)較大,因此本文將12 月定為前冬,1~2 月定為后冬。

    3.2 年代際變化階段的確定

    使用Mann-Kendall 非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法(魏鳳英, 2007),檢驗(yàn)1951~2019 年前冬(12 月)、后冬(1~2 月)氣溫序列的突變特征(圖4),給定置信水平95%的臨界線??梢钥闯觯岸?、后冬氣溫具有相似的氣候變化特征:20 世紀(jì)80 年代中期開始,前冬、后冬氣溫均存在明顯的上升趨勢;前冬增暖趨勢在20 世紀(jì)90 年代后期超過95%的置信水平,后冬增暖趨勢較前冬超前且明顯,在20 世紀(jì)90 年代初超過95%的置信水平,甚至超過99.9%的置信水平,表明后冬氣溫的上升趨勢十分顯著。根據(jù)UF(Forward)和UB(Backward)曲線交點(diǎn)的位置,前冬、后冬氣溫增暖趨勢的突變點(diǎn)是一致的,均在1986 年開始。季節(jié)內(nèi)氣溫在年代際尺度上的突變特征與我國冬季氣溫在1985 年前后發(fā)生轉(zhuǎn)折的特征(丁一匯和張莉, 2008; 李維京等, 2013)一致。

    圖4 1951~2019 年中國(a)前冬、(b)后冬氣溫Mann-Kendall 檢驗(yàn)曲線,點(diǎn)線為95%的置信水平臨界線Fig. 4 Statistical curves of the Mann-Kendall method test of temperature for (a) early and (b) late winter over China for 1951-2019. The dotted lines are the critical value of 95% confidence level

    進(jìn)一步計(jì)算氣溫在1986 年突變前后冬季不同月份間全國160 站平均氣溫的相關(guān)系數(shù)(表2),用以驗(yàn)證不同年代際尺度背景下前冬與后冬劃分的合理性。從表中看出,無論是否去除線性趨勢,1986 年前后12 月與1 月、2 月氣溫間的相關(guān)系數(shù)均不顯著,而1 月與2 月氣溫間的聯(lián)系更緊密。尤其是1986 年之后,1 月與2 月氣溫的相關(guān)性增強(qiáng),超過了95%的置信水平。從相關(guān)系數(shù)空間分布(圖5)看,1986 年之前,12 月與1 月氣溫相關(guān)僅在西南部分地區(qū)、內(nèi)蒙古局部地區(qū)存在顯著性(圖5a),而1 月與2 月的顯著正相關(guān)區(qū)域明顯增加,主要分布在我國西北大部和內(nèi)蒙古地區(qū)(圖5b)。1986 年之后,12 月與1 月氣溫僅在東北中部地區(qū)相關(guān)顯著(圖5c),1 月與2 月的正相關(guān)顯著區(qū)域明顯擴(kuò)大(圖5d),顯著性區(qū)域東移,主要分布在東北、華北、西北地區(qū)東部、西南地區(qū)東北部等地區(qū)。保留線性趨勢的結(jié)果分析差別不大,只是1986 年之后的1 月與2 月在長江中游的正相關(guān)顯著區(qū)域略有擴(kuò)大(圖略)。上述分析顯示了前冬與后冬的關(guān)系在兩個(gè)不同年代際階段的差異。

    表2 1986 年 之 前(1951~1985 年)和 之 后(1986~2019 年)冬季不同月份間全國160 站平均氣溫的相關(guān)系數(shù)Table 2 Correlation coefficients of surface air temperature averaged in 160 stations over China between different months in winter before 1986 (1951-1985) and after 1986(1986-2019)

    圖5 中國前、后冬氣溫在1986 年(a、b)之前(1951~1985 年)和(c、d)之后(1986~2019 年)階段的(a、c)12 月與1 月、(b、d)1 月與2 月相關(guān)系數(shù)空間分布。深、淺陰影分別表示通過99%和95%的置信水平Fig. 5 Spatial distributions of correlation coefficients for temperature between (a, c) December and January, (b, d) January and February over China before 1986 (1951-1985) and after 1986 (1986-2019). Dark and light shadings indicate 99% and 95% confidence levels, respectively

    基于以上分析結(jié)果,同時(shí)考慮到前冬、后冬的年代際變化,去掉1990 年前后各5 a,本文將前冬、后冬的年代際尺度變化時(shí)段劃分為1951~1984 年(P1 時(shí)段)和1996~2019 年(P2 時(shí)段)。以下分別針對這兩個(gè)時(shí)段及各自的氣候態(tài)進(jìn)行分析。

    4 冬季氣溫在不同年代際時(shí)間尺度季節(jié)內(nèi)變率的優(yōu)勢模態(tài)及成因

    基于上一節(jié)對年代際時(shí)間尺度的劃分,參考國家標(biāo)準(zhǔn)《GB/T 21983-2020 暖冬等級》(國家市場監(jiān)督管理總局和國家標(biāo)準(zhǔn)化管理委員會(huì), 2020)的閾值標(biāo)準(zhǔn),以1986 年前/后全國160 站平均的前冬、后冬氣溫標(biāo)準(zhǔn)化時(shí)間序列的0.33 個(gè)標(biāo)準(zhǔn)差為異常標(biāo)準(zhǔn),統(tǒng)計(jì)冷/暖期全國160 站平均的前冬、后冬氣溫明顯異常年的組合(表3)發(fā)現(xiàn),1986 年之前前冬、后冬一致偏暖或偏冷為6 年,比例為17.7%;而前冬和后冬冷暖反位相年為8 年,比例為23.6%。1986 年之后前冬、后冬一致偏暖或偏冷年為6 年,比例為25.0%;前冬和后冬冷暖反位相年為7 年,為29.2%。1986 年之后較冷期前冬、后冬持續(xù)偏暖的比例增加,持續(xù)偏冷的比例下降;1986 年之后前冬暖后冬冷的比例上升,前冬冷后冬暖的比例也上升,即1986 年后前冬和后冬冷暖反位相比例較1986 年之前升高,季節(jié)內(nèi)變率加大。

    表3 1986 年前(1951~1984 年)、后(1996~2019 年)中國160 站平均的前冬、后冬氣溫顯著異常年數(shù)統(tǒng)計(jì)Table 3 Statistics of the number of years for significantly abnormal temperature based on 160 stations over China in early and late winter before 1986 (1951-1984)/after 1986(1996~2019)

    為認(rèn)識不同年代際時(shí)間尺度下前、后冬氣溫演變的空間分布特征,分別以P1 和P2 時(shí)段全國160 站前冬、后冬氣溫距平場作為2 個(gè)變量構(gòu)成矩陣進(jìn)行S-EOF 分解(Wang and An, 2005),得到代表前冬和后冬的平均氣溫距平隨時(shí)間演變的空間模態(tài)及時(shí)間系數(shù)序列。為直觀分析實(shí)際的氣溫變化量級,文中用標(biāo)準(zhǔn)化的時(shí)間系數(shù)序列對氣溫場進(jìn)行回歸來表示S-EOF 的空間模態(tài)(韋瑋等, 2014,2020)。

    東亞冬季風(fēng)是影響我國冬季氣溫變化的關(guān)鍵環(huán)流系統(tǒng),冬季風(fēng)系統(tǒng)成員主要包括對流層低層的西伯利亞高壓和阿留申低壓,對流層中層的東亞大槽以及對流層高層的西風(fēng)急流等,環(huán)流系統(tǒng)的異常對我國冬季氣溫異常具有直接作用。為了研究不同氣候背景下我國前、后冬氣溫演變優(yōu)勢模態(tài)對應(yīng)的環(huán)流系統(tǒng)特征,文中采用合成分析的方法對海平面氣壓場、500 hPa 高度場、200 hPa 緯向風(fēng)場進(jìn)行分析,正負(fù)典型異常年的選取如表4 所示,以S-EOF 第一模態(tài)標(biāo)準(zhǔn)化時(shí)間系數(shù)正負(fù)0.5 倍標(biāo)準(zhǔn)差為標(biāo)準(zhǔn)選取典型年,分別選取P1 階段典型正異常年10 年,典型負(fù)異常年7 年;P2 階段典型正、負(fù)異常年各6 年,并用t檢驗(yàn)方法分析合成結(jié)果的顯著性。同時(shí)利用S-EOF1 標(biāo)準(zhǔn)化時(shí)間系數(shù)回歸前冬和后冬的環(huán)流場(圖略),獲得的結(jié)果與正、負(fù)異常年差值場合成信息一致。

    表4 P1 階段(1951~1984 年)和P2 階段(1996~2019 年)中國前冬、后冬氣溫S-EOF 第一模態(tài)標(biāo)準(zhǔn)化時(shí)間序列典型正、負(fù)異常年份Table 4 Typical positive and negative anomaly years for surface air temperature as identified in terms of the first S-EOF(season-reliant empirical orthogonal function) mode between early and late winter over China during P1 (1951-1984) and P2(1996-2019)

    4.1 1986 年之前(1951~1984 年)冬季氣溫季節(jié)內(nèi)變率的優(yōu)勢模態(tài)及環(huán)流特征

    圖6 給出1986 年之前(P1 階段)前冬、后冬氣溫距平S-EOF 分解前兩個(gè)模態(tài),其方差貢獻(xiàn)分別為36.1%和20.9%。S-EOF1 空間分布的前后冬演變(圖6a、b)顯示,前冬的氣溫異常變化幅度更大更顯著,除青藏高原東部地區(qū)外,均通過99%的置信水平,且全國氣溫基本為一致負(fù)異常變化型,氣溫變化的異常幅度從南向北遞增,最大值位于內(nèi)蒙古中部和新疆北部地區(qū)。后冬異常變化的幅度雖然不顯著,但是除內(nèi)蒙古和西北部分地區(qū)為弱的負(fù)異常外,全國其余地區(qū)均為正異常,與前冬的氣溫變化趨勢相反。S-EOF2 空間分布同樣存在區(qū)域性和季節(jié)內(nèi)演變差異(圖6c、d),前冬為弱暖位相,后冬為顯著冷位相??梢钥闯?,SEOF1 和S-EOF2 均表現(xiàn)出前后冬冷暖反位相的特征。限于篇幅,本文主要對S-EOF1 進(jìn)行分析。

    圖6 P1 階段(1951~1984 年)中國(a、c)前冬、(b、d)后冬氣溫距平S-EOF 分解(a、b)第一、(c、d)第二模態(tài)標(biāo)準(zhǔn)化時(shí)間序列回歸到氣溫場的空間分布及S-EOF1(e)標(biāo)準(zhǔn)化時(shí)間系數(shù)(SPC1)和(f)小波功率譜。圖a-d 中,實(shí)、虛等值線表示正、負(fù)值;深、淺陰影分別表示通過99%和95%的置信水平。圖f 中,黑線表示小波變換的影響錐曲線,紅色實(shí)線表示通過95%的置信水平,紅色虛線為95%置信水平的紅噪聲檢驗(yàn),藍(lán)色實(shí)線表示小波全譜Fig. 6 Spatial distributions of (a, b) the first and (c, d) second S-EOF mode regression on surface air temperature for (a, c) early and (b, d) late winter,(e) normalized principal component (SPC1) and (f) wave power spectrum of the first S-EOF mode over China during P1 period (1951-1984). In Figs.a-d, solid (dashed) contours represent positive (negative) values; dark and light shadings indicate the 99% and 95% confidence levels, respectively. In Fig. f, the black thick lines are the influence cone curve of wavelet transform; the red solid line represents 95% confidence level; the red dashed line represents 95% confidence level of the red noise test; the blue solid line is global wavelet power spectrum

    S-EOF1 模態(tài)對應(yīng)的標(biāo)準(zhǔn)化時(shí)間系數(shù)(圖6e)顯示比較清楚的年際變化特征。對時(shí)間系數(shù)進(jìn)行小波功率譜分析,主要有2~3 a 和7~10 a 周期,其中2~3 a 周期最為顯著,通過95%置信水平的紅噪聲檢驗(yàn),主要發(fā)生在20 世紀(jì)50 年代和60 年代中后期;7~10 a 周期主要發(fā)生在1970 年后,沒有通過顯著性檢驗(yàn)(圖6f)。

    圖7 為冬季氣溫前冬和后冬S-EOF1 模態(tài)典型正、負(fù)異常年的大氣環(huán)流合成差值場。當(dāng)冬季氣溫出現(xiàn)前冬冷后冬暖的特征時(shí),前冬,海平面氣壓場上歐亞大陸中高緯度地區(qū)為顯著正異常,北極地區(qū)為負(fù)異常,西伯利亞高壓明顯偏強(qiáng);北大西洋北部地區(qū)為顯著正異常,冰島低壓偏弱;北太平洋濤動(dòng)(NPO)處于正位相,即阿留申低壓和夏威夷高壓同時(shí)增強(qiáng),均通過95%的置信水平,導(dǎo)致東亞—北太平洋地區(qū)海陸氣壓梯度增大,東亞冬季風(fēng)增強(qiáng),有利于我國冬季氣溫大范圍偏低(圖7a)。后冬,歐亞大陸中高緯度地區(qū)環(huán)流異常較前冬減弱,北極地區(qū)維持負(fù)異常,強(qiáng)度較前冬減弱;烏拉爾山地區(qū)轉(zhuǎn)為顯著負(fù)距平中心,負(fù)值區(qū)域向東南延伸至貝加爾湖附近;NPO 維持正位相,阿留申群島仍為負(fù)距平,但強(qiáng)度減弱;北大西洋北部地區(qū)轉(zhuǎn)為負(fù)距平中心,冰島低壓加強(qiáng),通過95%的置信水平。該特征顯示后冬北方地區(qū)維持弱冷的主要環(huán)流信號來自NPO 特征的維持,而北大西洋到歐亞地區(qū)的環(huán)流經(jīng)向度減弱,影響我國冷暖的信號出現(xiàn)不一致,從而使得我國后冬偏冷的可能性減弱,出現(xiàn)有利于偏暖的轉(zhuǎn)折特征(圖7d)。

    對流層中層500 hPa 高度場顯示,前冬,從大西洋北部—西歐—烏拉爾山至東亞北部地區(qū)表現(xiàn)為“+ - + -”的典型歐亞型(EU)大氣遙相關(guān)正位相(Wallace and Gutzler, 1981; 汪寧等, 2017)特征。歐亞中高緯地區(qū)的環(huán)流異常特征顯著,呈“西北高東南低”的空間分布,烏拉爾山附近的廣大地區(qū)為正距平,而貝加爾湖以東和以南到我國30°N 附近的大部分地區(qū)為負(fù)距平中心,均通過99%的置信水平,副高偏弱偏東。烏拉爾山高壓偏強(qiáng),東亞大槽加深且位置偏西,有利于引導(dǎo)冷空氣影響我國大范圍地區(qū)(圖7b)。后冬,歐亞中高緯度地區(qū)異常信號較前冬減弱,為弱的歐亞型(EU)大氣遙相關(guān)型負(fù)位相。烏拉爾山以西轉(zhuǎn)為顯著負(fù)距平中心,高壓脊減弱,但是仍通過95%的置信水平,負(fù)距平區(qū)域向東南延伸至我國中西部地區(qū);貝加爾湖以東至鄂霍茨克海地區(qū)為正距平,東亞大槽較前冬變淺東移,副高偏強(qiáng)偏西,我國東部地區(qū)為高度場正異常(圖7e)。環(huán)流特征顯示冷空氣活動(dòng)明顯減弱,導(dǎo)致氣溫偏低程度較前冬減弱,我國西北部地區(qū)上空受高度場負(fù)異常的影響,而南方地區(qū)受高度場正異常的控制,有利于北方地區(qū)氣溫維持弱冷狀態(tài),而南方地區(qū)氣溫轉(zhuǎn)暖(圖6b)。

    許多研究指出,東亞高空急流是影響中國冬季氣溫的重要環(huán)流系統(tǒng)(高輝, 2007; 楊輝和李崇銀,2008; 況雪源等, 2009; 彭京備等, 2020),東亞副熱帶西風(fēng)急流偏強(qiáng)時(shí),東亞冬季風(fēng)增強(qiáng)(Yang et al.,2002; Jhun and Lee, 2004; 況雪源等, 2008; Liu et al.,2012)。對流層上層200 hPa 緯向風(fēng)場顯示,前冬,歐亞地區(qū)自北向南呈現(xiàn)“+ - + -”的距平分布(圖7c),大部地區(qū)可通過99%和95%的置信水平;東亞上空40°N 南北兩側(cè)緯向風(fēng)異常符號相反,即東亞溫帶急流區(qū)為顯著東風(fēng)異常,副熱帶西風(fēng)急流區(qū)為顯著西風(fēng)異常,緯向風(fēng)大值中心位于日本群島及其以西地區(qū),東亞溫帶急流偏弱,副熱帶西風(fēng)急流偏強(qiáng),有利于產(chǎn)生氣旋式環(huán)流異常進(jìn)而導(dǎo)致東亞大槽加深,同時(shí)在20°N 以南的南亞和東南亞地區(qū)上空為東風(fēng)距平,即副熱帶急流區(qū)南側(cè)緯向風(fēng)偏弱,產(chǎn)生反氣旋式環(huán)流異常,進(jìn)一步導(dǎo)致東亞冬季風(fēng)增強(qiáng),有利于冷空氣南侵而使我國氣溫偏低(賀圣平和王會(huì)軍, 2012; 姚慧茹和李棟梁, 2013a, 2013b)。后冬,東亞副熱帶西風(fēng)急流區(qū)的西風(fēng)異常強(qiáng)度減弱,顯著性也相應(yīng)地減弱,西風(fēng)距平范圍向北擴(kuò)大至60°N附近,從而不利于冷空氣大規(guī)模南侵;熱帶地區(qū)的東風(fēng)距平范圍隨著副熱帶西風(fēng)急流的北推而向北擴(kuò)大至30°N 附近,有利于我國南方氣溫的回升(圖7f)。

    圖7 P1 階段(1951~1984 年)中國前冬S-EOF 優(yōu)勢模態(tài)典型正、負(fù)異常年的合成差值場:(a)海平面氣壓場(間隔:2 hPa);(b)500 hPa 位勢高度場(間隔:20 gpm);(c)200 hPa 緯向風(fēng)場(間隔:2 m/s)。深、淺陰影分別表示通過99%和95%的置信水平。圖c 中,紅色輪廓線表示200 hPa 緯向風(fēng)30 m/s 的氣候平均值。(d-f)同(a-c),但為后冬的合成差值場Fig. 7 Composite differences between the typical positive and negative anomaly years of S-EOF dominant mode in EW (early winter) over China during P1 period (1951-1984): (a) Sea level pressure (contours interval: 2 hPa); (b) 500-hPa geopotential height (contours interval: 20 gpm); (c) 200-hPa zonal wind (contours interval: 2 m/s). Dark and light shadings indicate 99% and 95% confidence levels, respectively. In Fig. c, the red contour denotes the climatic mean of the 200-hPa zonal wind with the value of 30 m/s. (d-f) As in (a-c), but for the composite differences in LW (late winter)

    上述分析顯示P1 階段S-EOF1 模態(tài)對應(yīng)的正異常年,對流層低中高層環(huán)流異常的顯著性在前冬明顯高于后冬。前冬和后冬對流層高中低層均表現(xiàn)出一致的季節(jié)內(nèi)變化特征。前冬,冬季風(fēng)環(huán)流整體偏強(qiáng),有利于全國大部地區(qū)氣溫偏低;后冬,歐亞地區(qū)對流層環(huán)流發(fā)生較大調(diào)整,冬季風(fēng)環(huán)流整體減弱,其中西北太平洋地區(qū)對流層中高層的環(huán)流調(diào)整更明顯,副熱帶高度場增強(qiáng),熱帶東風(fēng)急流北擴(kuò),有利于我國南部地區(qū)氣溫升高,造成季節(jié)內(nèi)反位相變率增大。負(fù)異常年對應(yīng)的特征亦成立。

    4.2 1986 年之后(1996~2019 年)冬季氣溫季節(jié)內(nèi)變率的優(yōu)勢模態(tài)及環(huán)流特征

    圖8 給出1986 年之后(P2 階段,1996~2019 年)前冬、后冬氣溫距平S-EOF 分解前兩個(gè)模態(tài),其方差貢獻(xiàn)分別為36.1%和20.6%。SEOF1 空間分布的前后冬演變(圖8a、b)顯示,后冬的氣溫異常變化幅度更大更顯著,除青藏高原東部地區(qū)外,均通過99%的置信水平,且全國氣溫基本為一致正異常變化型,異常變化幅度從西南至東北遞增,最大值位于東北地區(qū)。前冬異常變化的幅度不顯著,空間負(fù)位相和正位相差異的零分界線從青藏高原東部沿長江上游向南轉(zhuǎn)到東南沿海地區(qū)。同P1 階段類似,P2 階段前、后冬第二模態(tài)空間分布同樣存在區(qū)域性和季節(jié)內(nèi)演變差異,前冬為顯著冷位相,后冬為弱暖位相。限于篇幅,本文主要對S-EO1 進(jìn)行分析。

    從S-EOF1 模態(tài)對應(yīng)的標(biāo)準(zhǔn)化時(shí)間系數(shù)(圖8e)看出,該模態(tài)年際變化特征更明顯,對時(shí)間系數(shù)進(jìn)行小波功率譜分析,主要有2~3 a、4a 和6a 的周期,均通過95%置信水平的紅噪聲檢驗(yàn)(圖8f),其中4 a 周期主要發(fā)生在20 世紀(jì)90 年代后期至21 世紀(jì)前期,2~3 a 周期主要發(fā)生在2005~2010年間和2015 年前后,6 a 周期主要發(fā)生在2005~2015 年間。

    圖8 同圖6,但為P2 階段(1996~2019 年)的空間分布Fig. 8 As in Fig. 6, but for spatial distributions during P2 period (1996-2019)

    SEOF1 模態(tài)顯示,我國北方地區(qū)(除東北外)大部到華中地區(qū)在前冬為負(fù)位相,后冬為顯著的正位相,即前后冬呈反位相變化;而南方地區(qū)中,西南大部前冬為正位相,后冬為負(fù)位相,即前后冬亦呈反位相變化;江南東部到華南東部前冬為正位相,后冬為顯著的正位相,基本呈同位相變化。即北方和西南大部地區(qū)前后冬氣溫變化不一致較東南部地區(qū)更明顯,前后冬氣溫變率更大。優(yōu)勢模態(tài)顯示前冬北方冷(暖)轉(zhuǎn)為后冬全國明顯暖(冷)的可能性大,即北方的季節(jié)內(nèi)反相更顯著。Ding et al.(2021)的最新研究顯示全球變暖后,華北地區(qū)出現(xiàn)強(qiáng)降溫事件增多,從極端事件頻發(fā)的角度也反映了北方地區(qū)冬季季節(jié)內(nèi)氣溫變率增大的現(xiàn)象。

    圖9 為P2 時(shí)段前、后冬S-EOF1 典型正、負(fù)異常年的大氣環(huán)流合成差值場。S-EOF1 模態(tài)對應(yīng)的正異常年,我國冬季氣溫出現(xiàn)前冬冷后冬明顯轉(zhuǎn)暖的變化特征時(shí),前冬海平面氣壓場顯示極區(qū)為顯著負(fù)距平,極渦偏向東北亞地區(qū),阿留申群島為負(fù)距平中心,通過了95%的置信水平;西伯利亞高壓略偏強(qiáng),但信號不顯著;中南半島、南海和西太平洋地區(qū)為正距平(圖9a)。環(huán)流特征顯示極地冷空氣擴(kuò)散,阿留申低壓較強(qiáng),冬季風(fēng)略偏強(qiáng),有利于我國北方大部氣溫略偏低,而西南和東南沿海受海平面氣壓場偏高影響,有利于氣溫略高,但信號不顯著。后冬,歐亞大陸海平面氣壓為顯著負(fù)異常,西伯利亞高壓顯著偏弱;北太平洋地區(qū)為正距平,阿留申低壓填塞,南亞、中南半島、南海到熱帶西太平洋地區(qū)均為正距平(圖9d)。環(huán)流特征顯示中高緯度冷空氣勢力弱,冬季風(fēng)偏弱;低緯度地區(qū)為寬廣的正距平,海平面氣壓偏強(qiáng),多數(shù)地區(qū)通過95%或99%的置信水平,歐亞大陸和南部海域的“北低南高”的海平面氣壓異常分布顯示冷空氣非常弱,從而有利于我國氣溫異常偏高。

    圖9 同圖7,但為P2 時(shí)段的合成差值場Fig. 9 As in Fig.7, but for composite differences during P2 period (1996-2019)

    500 hPa 高度場顯示,前冬極區(qū)主體為負(fù)距平,極渦中心伸展至鄂霍茨克海地區(qū),60°N 以北為負(fù)距平區(qū),負(fù)距平中心通過95%的置信水平;環(huán)繞極區(qū)的60°N 以南中緯度地區(qū)為正距平區(qū),呈“北低南高”的分布,烏拉爾山及以西高度場略偏強(qiáng),但異常不顯著;印度洋到西太平洋的低緯度地區(qū)高度場異常偏強(qiáng)(圖9b)。這種分布有利于弱冷空氣活動(dòng),從而造成北方弱冷,南方弱暖的特征。后冬,歐亞大陸高中低緯度地區(qū)環(huán)流異常顯著,歐亞大陸中高緯度地區(qū)表現(xiàn)為類似極地—?dú)W亞型(POL 型)遙相關(guān)(Barnston and Livezey, 1987)的南北向偶極子分布形態(tài),50°N 以北為負(fù)距平區(qū),50°N 以南為正距平區(qū),異常信號顯著,大部分地區(qū)通過99%的置信水平,烏拉爾山地區(qū)為顯著負(fù)距平區(qū),正距平中心從巴爾喀什湖延伸到貝加爾湖及日本群島上空,東亞大槽填塞,歐亞以緯向型環(huán)流為主,同時(shí)低緯度地區(qū)高度場異常偏強(qiáng),冷空氣活動(dòng)非常弱,有利于我國氣溫顯著偏高(圖9e)。

    對流層上層200 hPa 緯向風(fēng)場顯示,前冬東亞副熱帶高空急流的異常不顯著,西段偏強(qiáng),東段偏弱(圖9c),溫帶急流略偏強(qiáng),二者的鋒區(qū)主要影響我國北方地區(qū)。后冬東亞副熱帶西風(fēng)急流總體偏弱,其中東段通過95%的置信水平,溫帶急流異常偏強(qiáng),通過99%的置信水平,抑制冷空氣南下,有利于我國氣溫顯著偏高(圖9f)。前冬與后冬在低緯度地區(qū)沒有顯著性差異。

    上述分析顯示,從前冬到后冬,對流層高中低層表現(xiàn)出比較一致的季節(jié)內(nèi)變化特征,且對流層中低層環(huán)流季節(jié)內(nèi)異常的顯著性更強(qiáng),后冬氣溫和環(huán)流異常的顯著性明顯高于前冬。前冬到后冬的歐亞中高緯度環(huán)流發(fā)生較大調(diào)整,而低緯度的環(huán)流變化不大,因此對應(yīng)模態(tài)的正異常年,我國北方地區(qū)前冬冷到后冬明顯轉(zhuǎn)暖,季節(jié)內(nèi)反位相變率大,而南方地區(qū)為弱暖到暖,季節(jié)內(nèi)變率較小。

    以上兩節(jié)的分析表明,不同年代際背景下冬季氣溫季節(jié)內(nèi)變化前后不一致的特征在增強(qiáng),季節(jié)內(nèi)大部分地區(qū)的反位相變化成為優(yōu)勢模態(tài),這一結(jié)論與表3 分析的1986 年前后中國前冬、后冬氣溫顯著異常年數(shù)統(tǒng)計(jì)結(jié)果基本一致。該結(jié)果顯示在不同氣候變化背景下考慮冬季季節(jié)內(nèi)變化特征與以往的研究(韋瑋等, 2014, 2020; 孫健等, 2019)有所不同,分別在P1 和P2 時(shí)段內(nèi)進(jìn)行S-EOF 分析,去除了年代際差異的影響,更多地體現(xiàn)出年際變化分量,季節(jié)內(nèi)變化特征也更加突出。祁莉和泮琬楠(2021)將東亞氣溫歷史數(shù)據(jù)進(jìn)行去趨勢和去氣候態(tài)處理,再進(jìn)行S-EOF 分解,也獲得S-EOF1 為冬季季節(jié)內(nèi)氣溫前后反相的結(jié)果。

    P1、P2 時(shí)段,冬季氣溫的季節(jié)內(nèi)變化受到冬季風(fēng)環(huán)流系統(tǒng)關(guān)鍵成員的季節(jié)內(nèi)變率的影響。P1時(shí)段的前冬東亞冬季風(fēng)環(huán)流異常幅度大,而P2 時(shí)段后冬環(huán)流異常顯著。根據(jù)朱艷峰(2008)定義的冬季風(fēng)指數(shù),分別計(jì)算兩個(gè)時(shí)段優(yōu)勢模態(tài)標(biāo)準(zhǔn)化時(shí)間系數(shù)(PC1)與冬季風(fēng)指數(shù)的相關(guān),二者在P1時(shí)段前冬的相關(guān)系數(shù)為0.85,P2 時(shí)段后冬的相關(guān)系數(shù)為-0.77,均通過99.9%的置信水平;但二者在P1 時(shí)段后冬與P2 時(shí)段前冬的相關(guān)性均不顯著,與上述對大氣環(huán)流系統(tǒng)的分析一致。

    文中重點(diǎn)分析了兩個(gè)時(shí)段優(yōu)勢模態(tài)對應(yīng)的正異常年氣溫變化(P1 時(shí)段前冬強(qiáng)冷、后冬弱暖和P2時(shí)段前冬弱冷、后冬強(qiáng)暖)及冬季風(fēng)環(huán)流特征,對于相反的演變模態(tài),冬季風(fēng)環(huán)流異??臻g分布分別與上述分析相似,但距平符號相反。

    5 總結(jié)與討論

    本文利用中國160 站月平均氣溫和再分析數(shù)據(jù),通過分析冬季各月之間氣溫的相關(guān)及穩(wěn)定性,確定12 月為前冬,1 月、2 月為后冬。根據(jù)前、后冬氣溫增暖趨勢的突變特征確定1986 年之前和之后為兩個(gè)年代際變化階段。進(jìn)一步利用S-EOF 方法研究不同年代際時(shí)間尺度下冬季季節(jié)內(nèi)氣溫變化的優(yōu)勢模態(tài)以及相應(yīng)的關(guān)鍵環(huán)流系統(tǒng)異常特征。1986 年前/后S-EOF1 的方差貢獻(xiàn)均為36.1%,對應(yīng)的標(biāo)準(zhǔn)化時(shí)間系數(shù)顯示清楚的年際變化特征,但空間分布特征及時(shí)間演變規(guī)律有很大差異。冬季氣溫的季節(jié)內(nèi)變化受冬季風(fēng)系統(tǒng)的季節(jié)內(nèi)調(diào)整影響顯著,1986 年之前,前冬東亞冬季風(fēng)環(huán)流異常幅度顯著;1986 年之后,后冬的環(huán)流異常幅度顯著。主要結(jié)論如下:

    (1)1986 年之前,S-EOF1 空間分布顯示前冬的氣溫異常變化更顯著,北方地區(qū)前冬為顯著的負(fù)位相,后冬為弱的負(fù)位相,前后冬呈同位相變化;而南方地區(qū)前冬為顯著的負(fù)位相,后冬為弱的正位相,前后冬呈反位相變化,即南方地區(qū)前后冬氣溫變率更大。優(yōu)勢模態(tài)對應(yīng)的正異常年,北太平洋濤動(dòng)(NPO)從前冬至后冬維持正位相,而歐亞型中高緯度大氣遙相關(guān)(EU)從正位相轉(zhuǎn)為弱的負(fù)位相,西伯利亞高壓由強(qiáng)變?nèi)?,西風(fēng)急流由偏強(qiáng)偏南轉(zhuǎn)為偏弱,對流層高中低層的環(huán)流異常均為一致的季節(jié)內(nèi)變化特征。后冬,歐亞地區(qū)對流層環(huán)流發(fā)生較大調(diào)整,冬季風(fēng)環(huán)流整體減弱,其中西北太平洋地區(qū)對流層中高層的環(huán)流調(diào)整更明顯,副熱帶高度場增強(qiáng),熱帶東風(fēng)急流北擴(kuò),有利于我國南部地區(qū)氣溫升高,造成南方地區(qū)季節(jié)內(nèi)反位相變率增大。

    (2)1986 年之后,S-EOF1 空間分布顯示后冬的氣溫異常變化更顯著,我國北方地區(qū)(除東北外)大部到華中地區(qū)在前冬為弱的負(fù)位相,后冬為顯著的正位相,即前后冬呈反位相變化;而江南東部到華南東部前冬為弱的正位相,后冬為顯著的正位相,基本呈同位相變化;即北方大部地區(qū)前后冬氣溫變率更大。優(yōu)勢模態(tài)對應(yīng)的正異常年,歐亞大陸中高緯度地區(qū)南北向偶極型遙相關(guān)增強(qiáng)南移,低緯度熱帶海洋地區(qū)的環(huán)流異常維持顯著性,即前冬和后冬的環(huán)流差異集中在歐亞大陸中高緯度區(qū),其中對流層中低層的環(huán)流異常更顯著;導(dǎo)致我國北方從前冬維持弱冷轉(zhuǎn)為后冬的明顯偏暖,造成北方地區(qū)季節(jié)內(nèi)反位相變率增大。

    本文分析了不同年代際時(shí)間尺度我國冬季氣溫季節(jié)內(nèi)變化的新特征,并對相應(yīng)的冬季風(fēng)環(huán)流系統(tǒng)季節(jié)內(nèi)變化異常進(jìn)行了分析。冬季風(fēng)環(huán)流系統(tǒng)還將受到其他環(huán)流異常和外強(qiáng)迫因子異常的影響。有研究顯示海溫和北極海冰作為全球氣候系統(tǒng)的重要組成部分,是影響大氣環(huán)流異常的重要因素,對我國冬季氣溫變化有重要的調(diào)制作用(陶詩言和張慶云, 1998; 張慶云和王媛, 2006; Wu et al., 2011, 2015;Li and Zhang, 2013, 2015),北大西洋濤動(dòng)和ENSO 事件也可能是造成東亞冬季氣溫前、后冬反相的影響因子(祁莉和泮琬楠, 2021),有必要深入探討不同年代際時(shí)間尺度下,影響冬季氣溫季節(jié)內(nèi)演變的物理機(jī)制和外強(qiáng)迫信號,以期為冬季氣候預(yù)測業(yè)務(wù)提供前兆信號和可預(yù)報(bào)性信息。

    猜你喜歡
    置信水平環(huán)流氣溫
    氣溫的神奇“魔法”
    基于FY-3D和FY-4A的氣溫時(shí)空融合
    內(nèi)環(huán)流控溫技術(shù)應(yīng)用實(shí)踐與發(fā)展前景
    產(chǎn)品控制與市場風(fēng)險(xiǎn)之間的相互作用研究
    單因子方差分析法在卷煙均勻性檢驗(yàn)中的研究與應(yīng)用
    熱鹽環(huán)流方程全局弱解的存在性
    用VaR方法分析中國A股市場的風(fēng)險(xiǎn)
    謎底大揭秘
    與氣溫成反比的東西
    兩相坐標(biāo)系下MMC環(huán)流抑制策略
    色播亚洲综合网| 亚洲精品av麻豆狂野| 久久香蕉国产精品| av中文乱码字幕在线| 久热爱精品视频在线9| 久久久久国产精品人妻aⅴ院| 精品少妇一区二区三区视频日本电影| 黑人操中国人逼视频| 亚洲欧美日韩高清在线视频| 人妻丰满熟妇av一区二区三区| 成人欧美大片| 欧美性猛交╳xxx乱大交人| 99久久国产精品久久久| 久久久久久久久久黄片| 国产又色又爽无遮挡免费看| 亚洲av成人不卡在线观看播放网| 日韩欧美国产一区二区入口| 日本在线视频免费播放| 欧美乱妇无乱码| 欧美成人性av电影在线观看| 男女床上黄色一级片免费看| 一本大道久久a久久精品| 精品久久久久久久久久免费视频| 欧美成人免费av一区二区三区| 一本久久中文字幕| 成年人黄色毛片网站| 久久九九热精品免费| 欧美成人性av电影在线观看| 最近最新免费中文字幕在线| 床上黄色一级片| www日本黄色视频网| 男女午夜视频在线观看| 小说图片视频综合网站| 日日夜夜操网爽| 国产熟女xx| 最近在线观看免费完整版| 久久香蕉国产精品| 亚洲国产精品sss在线观看| 女人爽到高潮嗷嗷叫在线视频| 中文在线观看免费www的网站 | 亚洲国产精品合色在线| 日韩欧美在线二视频| 国产成人aa在线观看| 99riav亚洲国产免费| 亚洲精品美女久久久久99蜜臀| 国产91精品成人一区二区三区| АⅤ资源中文在线天堂| 丰满的人妻完整版| 十八禁网站免费在线| 老汉色av国产亚洲站长工具| www.999成人在线观看| 99热只有精品国产| 88av欧美| 少妇被粗大的猛进出69影院| 国产精品久久久人人做人人爽| 亚洲国产精品sss在线观看| 2021天堂中文幕一二区在线观| 国产爱豆传媒在线观看 | 久久久久久久精品吃奶| 这个男人来自地球电影免费观看| 99国产综合亚洲精品| 久久久久久大精品| or卡值多少钱| 美女免费视频网站| 欧美日韩黄片免| 曰老女人黄片| 婷婷亚洲欧美| 在线国产一区二区在线| 成人一区二区视频在线观看| 精品久久久久久久久久久久久| 亚洲成a人片在线一区二区| 99在线视频只有这里精品首页| av超薄肉色丝袜交足视频| 免费看a级黄色片| 99久久综合精品五月天人人| 极品教师在线免费播放| 国产精品免费视频内射| 高清在线国产一区| 毛片女人毛片| 黄色丝袜av网址大全| 亚洲欧美日韩东京热| 久久久久久免费高清国产稀缺| 国产av在哪里看| 国产一区二区在线观看日韩 | 精品午夜福利视频在线观看一区| 午夜免费激情av| 亚洲精品久久成人aⅴ小说| 麻豆一二三区av精品| 欧美久久黑人一区二区| 精品久久久久久成人av| 美女扒开内裤让男人捅视频| 国产精品 欧美亚洲| 国产aⅴ精品一区二区三区波| 精品国产乱子伦一区二区三区| 国产1区2区3区精品| 在线观看舔阴道视频| 国产精品乱码一区二三区的特点| 国产精品1区2区在线观看.| 国产午夜福利久久久久久| 一区二区三区高清视频在线| 99国产极品粉嫩在线观看| 久久久久久大精品| 美女扒开内裤让男人捅视频| 熟女电影av网| 国产欧美日韩精品亚洲av| 欧美日韩福利视频一区二区| 欧美乱色亚洲激情| 国产欧美日韩一区二区三| a在线观看视频网站| 亚洲人成网站高清观看| 91字幕亚洲| 国产精品自产拍在线观看55亚洲| 人妻久久中文字幕网| 亚洲欧美精品综合一区二区三区| 国产高清视频在线观看网站| 欧美色欧美亚洲另类二区| 国产成+人综合+亚洲专区| 亚洲成人中文字幕在线播放| 18禁裸乳无遮挡免费网站照片| 最近最新中文字幕大全电影3| 亚洲中文字幕日韩| 欧美日韩中文字幕国产精品一区二区三区| 国产单亲对白刺激| 一区二区三区激情视频| or卡值多少钱| 亚洲avbb在线观看| 999精品在线视频| 床上黄色一级片| 很黄的视频免费| 搞女人的毛片| 国产精品久久久人人做人人爽| 男人的好看免费观看在线视频 | 在线永久观看黄色视频| av福利片在线| 欧美日韩国产亚洲二区| 香蕉av资源在线| 亚洲色图 男人天堂 中文字幕| 亚洲 欧美一区二区三区| 神马国产精品三级电影在线观看 | av福利片在线观看| 琪琪午夜伦伦电影理论片6080| 无人区码免费观看不卡| 欧美日韩精品网址| 国产午夜精品久久久久久| 老司机午夜十八禁免费视频| 久久国产精品影院| 国产精品爽爽va在线观看网站| 久久精品国产99精品国产亚洲性色| 欧美不卡视频在线免费观看 | 亚洲中文日韩欧美视频| 婷婷六月久久综合丁香| 在线观看免费午夜福利视频| 婷婷亚洲欧美| 日本精品一区二区三区蜜桃| 国产单亲对白刺激| 三级国产精品欧美在线观看 | 日本撒尿小便嘘嘘汇集6| 亚洲国产精品合色在线| 麻豆国产97在线/欧美 | 婷婷精品国产亚洲av在线| 18禁美女被吸乳视频| 国产片内射在线| 免费在线观看影片大全网站| 两个人视频免费观看高清| 一本精品99久久精品77| 亚洲精品一区av在线观看| 久久久久久久精品吃奶| 真人做人爱边吃奶动态| 国产免费av片在线观看野外av| 久久九九热精品免费| 在线国产一区二区在线| 亚洲成人国产一区在线观看| 最新在线观看一区二区三区| 欧美成人午夜精品| 国产伦人伦偷精品视频| 看免费av毛片| 五月伊人婷婷丁香| 久99久视频精品免费| 欧美 亚洲 国产 日韩一| 精品国产美女av久久久久小说| 久9热在线精品视频| 欧美乱码精品一区二区三区| 亚洲,欧美精品.| 非洲黑人性xxxx精品又粗又长| 亚洲最大成人中文| 变态另类丝袜制服| 国产高清videossex| 男人舔女人的私密视频| 久久天躁狠狠躁夜夜2o2o| 日日爽夜夜爽网站| 一a级毛片在线观看| tocl精华| 在线十欧美十亚洲十日本专区| 国产私拍福利视频在线观看| 99热6这里只有精品| 黄色成人免费大全| 嫩草影院精品99| 国内久久婷婷六月综合欲色啪| 天天躁夜夜躁狠狠躁躁| 久久久久九九精品影院| 天堂√8在线中文| 国产精品亚洲av一区麻豆| 亚洲成人免费电影在线观看| 狠狠狠狠99中文字幕| 99热这里只有是精品50| 在线观看一区二区三区| 99在线视频只有这里精品首页| 欧美性长视频在线观看| 亚洲av电影在线进入| 日本三级黄在线观看| 大型黄色视频在线免费观看| 一本久久中文字幕| 老司机在亚洲福利影院| 国产亚洲精品久久久久5区| 久久国产精品人妻蜜桃| 日本 av在线| 99久久国产精品久久久| 在线观看66精品国产| 最近视频中文字幕2019在线8| 日韩欧美国产在线观看| 巨乳人妻的诱惑在线观看| 久久香蕉激情| 在线a可以看的网站| 极品教师在线免费播放| 淫妇啪啪啪对白视频| 国产在线观看jvid| 免费在线观看亚洲国产| 日韩欧美 国产精品| 国产精品久久久久久人妻精品电影| 精品一区二区三区av网在线观看| 久久亚洲精品不卡| 国产日本99.免费观看| 麻豆成人av在线观看| 亚洲九九香蕉| 怎么达到女性高潮| 精品国产美女av久久久久小说| xxxwww97欧美| 在线观看免费视频日本深夜| а√天堂www在线а√下载| 怎么达到女性高潮| 亚洲,欧美精品.| 久久久久国产一级毛片高清牌| 免费在线观看影片大全网站| 黄片大片在线免费观看| 午夜激情av网站| 欧美性猛交黑人性爽| www.熟女人妻精品国产| 啦啦啦免费观看视频1| 欧美成狂野欧美在线观看| 欧美绝顶高潮抽搐喷水| 欧美性猛交黑人性爽| 中文字幕最新亚洲高清| 日本 av在线| 国产日本99.免费观看| 精品国产亚洲在线| 这个男人来自地球电影免费观看| 白带黄色成豆腐渣| 欧美日韩乱码在线| 亚洲成av人片在线播放无| 国产成年人精品一区二区| 99久久综合精品五月天人人| 国产三级在线视频| 18禁黄网站禁片午夜丰满| 久久香蕉精品热| 制服丝袜大香蕉在线| 看免费av毛片| 18禁国产床啪视频网站| 亚洲第一欧美日韩一区二区三区| 国产av在哪里看| 后天国语完整版免费观看| 亚洲一区中文字幕在线| 精品国产乱码久久久久久男人| 久久久精品大字幕| 1024香蕉在线观看| 日本在线视频免费播放| 一a级毛片在线观看| x7x7x7水蜜桃| 国内少妇人妻偷人精品xxx网站 | 最新美女视频免费是黄的| 在线视频色国产色| 欧美丝袜亚洲另类 | 国产私拍福利视频在线观看| 欧美国产日韩亚洲一区| 男女下面进入的视频免费午夜| 国产日本99.免费观看| 免费在线观看影片大全网站| 久久九九热精品免费| 欧美日韩精品网址| 人成视频在线观看免费观看| 亚洲一区中文字幕在线| 成人午夜高清在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 99riav亚洲国产免费| 少妇熟女aⅴ在线视频| 露出奶头的视频| 国产成人av激情在线播放| 久久亚洲精品不卡| 国产激情久久老熟女| 又粗又爽又猛毛片免费看| 国产精华一区二区三区| 欧美成人午夜精品| 这个男人来自地球电影免费观看| 一级黄色大片毛片| 啪啪无遮挡十八禁网站| 免费无遮挡裸体视频| 国内毛片毛片毛片毛片毛片| netflix在线观看网站| 五月玫瑰六月丁香| 50天的宝宝边吃奶边哭怎么回事| 深夜精品福利| 两个人看的免费小视频| 亚洲成人国产一区在线观看| 国产又黄又爽又无遮挡在线| 丰满人妻一区二区三区视频av | 日韩av在线大香蕉| 亚洲av五月六月丁香网| 亚洲一码二码三码区别大吗| 亚洲午夜精品一区,二区,三区| 叶爱在线成人免费视频播放| 国产成+人综合+亚洲专区| 90打野战视频偷拍视频| 色综合婷婷激情| 床上黄色一级片| 久久热在线av| 免费电影在线观看免费观看| 九色国产91popny在线| 国产又色又爽无遮挡免费看| 国产激情久久老熟女| www日本黄色视频网| 国产97色在线日韩免费| 亚洲国产精品sss在线观看| av片东京热男人的天堂| 男人舔女人的私密视频| 久久久精品国产亚洲av高清涩受| 一区二区三区高清视频在线| av福利片在线| 91字幕亚洲| 在线观看日韩欧美| ponron亚洲| 麻豆国产av国片精品| 久久精品91蜜桃| 中国美女看黄片| 国产精品99久久99久久久不卡| 亚洲七黄色美女视频| 一本综合久久免费| 搡老妇女老女人老熟妇| 国内精品一区二区在线观看| 一个人观看的视频www高清免费观看 | 看片在线看免费视频| 国产精品,欧美在线| 久久久国产成人精品二区| 动漫黄色视频在线观看| 巨乳人妻的诱惑在线观看| 日韩欧美国产一区二区入口| 天天添夜夜摸| 成年女人毛片免费观看观看9| 色老头精品视频在线观看| 欧美黄色片欧美黄色片| 国产精品 国内视频| 日本三级黄在线观看| 久久久精品国产亚洲av高清涩受| 高潮久久久久久久久久久不卡| 国产精品影院久久| 免费人成视频x8x8入口观看| 亚洲 欧美一区二区三区| av国产免费在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲国产欧美网| 亚洲片人在线观看| 日本一区二区免费在线视频| 国产精品一区二区三区四区免费观看 | 久久久国产欧美日韩av| 国产成人影院久久av| 欧美高清成人免费视频www| 国产亚洲欧美在线一区二区| 99久久99久久久精品蜜桃| 亚洲精品美女久久av网站| 一级毛片女人18水好多| 亚洲精品美女久久久久99蜜臀| 夜夜躁狠狠躁天天躁| 天堂√8在线中文| 99精品在免费线老司机午夜| 亚洲av电影在线进入| 在线观看舔阴道视频| 丝袜美腿诱惑在线| 天天一区二区日本电影三级| 亚洲精华国产精华精| 久久国产精品影院| 午夜精品在线福利| 国产精品综合久久久久久久免费| 无人区码免费观看不卡| 日韩欧美国产一区二区入口| 成人国语在线视频| 亚洲片人在线观看| 男女那种视频在线观看| 国产精品九九99| 亚洲av中文字字幕乱码综合| 中文字幕熟女人妻在线| 久久久久免费精品人妻一区二区| 日本免费a在线| 不卡av一区二区三区| cao死你这个sao货| 一边摸一边做爽爽视频免费| 欧美成狂野欧美在线观看| 久久婷婷成人综合色麻豆| 久久久国产成人免费| 日本免费一区二区三区高清不卡| 夜夜夜夜夜久久久久| 在线视频色国产色| 国产精品永久免费网站| www.999成人在线观看| 日日爽夜夜爽网站| 亚洲无线在线观看| 成人永久免费在线观看视频| cao死你这个sao货| 国产精品久久久久久久电影 | 狠狠狠狠99中文字幕| 老司机在亚洲福利影院| 禁无遮挡网站| 伊人久久大香线蕉亚洲五| 97超级碰碰碰精品色视频在线观看| 精品久久久久久,| 久久久久国产一级毛片高清牌| 午夜a级毛片| 正在播放国产对白刺激| 国产黄a三级三级三级人| 午夜日韩欧美国产| 成人午夜高清在线视频| 欧美成人午夜精品| 欧美人与性动交α欧美精品济南到| 日本免费一区二区三区高清不卡| 国产精品亚洲av一区麻豆| 国产精品香港三级国产av潘金莲| xxx96com| 国产精品日韩av在线免费观看| 国产成人精品久久二区二区免费| 狠狠狠狠99中文字幕| 亚洲专区国产一区二区| 1024手机看黄色片| 日韩欧美国产一区二区入口| 人妻久久中文字幕网| 免费av毛片视频| 欧美久久黑人一区二区| 露出奶头的视频| 日本熟妇午夜| 高清毛片免费观看视频网站| 日韩免费av在线播放| 欧美色欧美亚洲另类二区| 巨乳人妻的诱惑在线观看| 在线播放国产精品三级| 国产三级黄色录像| 日本熟妇午夜| 日本一二三区视频观看| 久久精品人妻少妇| 久久久久国产一级毛片高清牌| 精品少妇一区二区三区视频日本电影| 午夜免费激情av| ponron亚洲| 制服诱惑二区| 一区二区三区激情视频| 亚洲国产看品久久| 亚洲欧美一区二区三区黑人| 中文资源天堂在线| 91成年电影在线观看| 国产高清videossex| 18禁黄网站禁片免费观看直播| 久久久国产成人免费| 国产精品一区二区精品视频观看| 51午夜福利影视在线观看| 伦理电影免费视频| 亚洲av成人不卡在线观看播放网| 香蕉丝袜av| 精品一区二区三区av网在线观看| 亚洲自拍偷在线| 午夜福利视频1000在线观看| 精品久久久久久久人妻蜜臀av| 亚洲国产精品999在线| 香蕉国产在线看| 国产精品免费视频内射| 一个人观看的视频www高清免费观看 | 婷婷精品国产亚洲av在线| 欧美zozozo另类| 成人永久免费在线观看视频| 国产精品爽爽va在线观看网站| 99国产精品99久久久久| 黄色视频不卡| 亚洲av成人精品一区久久| 久久精品国产清高在天天线| 亚洲国产精品久久男人天堂| av在线播放免费不卡| 人妻久久中文字幕网| 欧美在线黄色| 国产v大片淫在线免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 首页视频小说图片口味搜索| 亚洲真实伦在线观看| xxxwww97欧美| 国产av在哪里看| av片东京热男人的天堂| 国产男靠女视频免费网站| 国产成人av教育| 国产成人系列免费观看| 观看免费一级毛片| 俄罗斯特黄特色一大片| 国产精品亚洲一级av第二区| 99热6这里只有精品| 欧美黑人精品巨大| 高潮久久久久久久久久久不卡| 香蕉丝袜av| 国产亚洲精品久久久久久毛片| 女警被强在线播放| 男女下面进入的视频免费午夜| 日韩大尺度精品在线看网址| 免费观看人在逋| 国产精品久久久久久久电影 | 成年人黄色毛片网站| 男女午夜视频在线观看| 成人18禁高潮啪啪吃奶动态图| a级毛片a级免费在线| 婷婷精品国产亚洲av| 久久精品91无色码中文字幕| 国产乱人伦免费视频| 少妇被粗大的猛进出69影院| 成在线人永久免费视频| 黑人操中国人逼视频| 国产精品日韩av在线免费观看| 91老司机精品| 草草在线视频免费看| www.熟女人妻精品国产| 国产欧美日韩一区二区三| 国产一区二区在线观看日韩 | 中文亚洲av片在线观看爽| 一本久久中文字幕| 亚洲国产高清在线一区二区三| 久久性视频一级片| 日本免费一区二区三区高清不卡| 免费在线观看黄色视频的| 在线观看66精品国产| 伊人久久大香线蕉亚洲五| 禁无遮挡网站| 变态另类成人亚洲欧美熟女| 色精品久久人妻99蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 久久精品91无色码中文字幕| 国内精品一区二区在线观看| 波多野结衣高清作品| 高潮久久久久久久久久久不卡| 岛国视频午夜一区免费看| 一级毛片女人18水好多| 19禁男女啪啪无遮挡网站| 天堂av国产一区二区熟女人妻 | 91av网站免费观看| 天堂av国产一区二区熟女人妻 | 日韩大码丰满熟妇| 国产午夜福利久久久久久| 午夜福利视频1000在线观看| 亚洲欧美精品综合一区二区三区| 一个人观看的视频www高清免费观看 | 床上黄色一级片| 波多野结衣巨乳人妻| 国产一区二区在线av高清观看| 脱女人内裤的视频| 亚洲av成人av| 欧美最黄视频在线播放免费| 亚洲国产欧美人成| 成人精品一区二区免费| 一个人免费在线观看的高清视频| 日日摸夜夜添夜夜添小说| 欧美黑人欧美精品刺激| 正在播放国产对白刺激| 国产高清激情床上av| 欧美另类亚洲清纯唯美| 午夜福利18| 久久这里只有精品19| 日韩免费av在线播放| 给我免费播放毛片高清在线观看| 国产精品1区2区在线观看.| 亚洲五月天丁香| 日本黄大片高清| 午夜老司机福利片| 无遮挡黄片免费观看| 中出人妻视频一区二区| 男女做爰动态图高潮gif福利片| 日本在线视频免费播放| 18禁国产床啪视频网站| 国产探花在线观看一区二区| 琪琪午夜伦伦电影理论片6080| 国产伦在线观看视频一区| 亚洲黑人精品在线| 亚洲电影在线观看av| 日韩有码中文字幕| 亚洲自拍偷在线| 法律面前人人平等表现在哪些方面| 黄片小视频在线播放| 看片在线看免费视频| 男女之事视频高清在线观看| 久热爱精品视频在线9| 99在线视频只有这里精品首页| 精品人妻1区二区| 午夜福利成人在线免费观看| 亚洲精品色激情综合| 亚洲av成人一区二区三| 岛国在线免费视频观看| 全区人妻精品视频| 18禁黄网站禁片午夜丰满| 美女免费视频网站| 九色国产91popny在线| av视频在线观看入口| 神马国产精品三级电影在线观看 | 国产黄a三级三级三级人| 日本成人三级电影网站| 成人av在线播放网站| 国产男靠女视频免费网站| 校园春色视频在线观看| 99热只有精品国产|