• <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)流抑制策略
    日韩一卡2卡3卡4卡2021年| 欧美日韩视频精品一区| 看免费成人av毛片| 亚洲视频免费观看视频| 国产精品久久久人人做人人爽| 日本vs欧美在线观看视频| a级毛片在线看网站| 精品一区二区三区四区五区乱码 | 日韩不卡一区二区三区视频在线| 色视频在线一区二区三区| 99久国产av精品国产电影| www.自偷自拍.com| 啦啦啦 在线观看视频| 最新的欧美精品一区二区| 国产97色在线日韩免费| 国产在视频线精品| 老熟女久久久| 久久久国产精品麻豆| 91aial.com中文字幕在线观看| 制服诱惑二区| 老鸭窝网址在线观看| 免费高清在线观看日韩| 国产欧美亚洲国产| 亚洲av日韩精品久久久久久密 | 美女国产高潮福利片在线看| 美女福利国产在线| 伦理电影免费视频| 亚洲国产中文字幕在线视频| √禁漫天堂资源中文www| 丰满乱子伦码专区| av在线播放精品| 日日爽夜夜爽网站| 免费黄色在线免费观看| 精品人妻在线不人妻| 狠狠婷婷综合久久久久久88av| 少妇人妻 视频| 少妇的丰满在线观看| 一级毛片电影观看| 国产精品久久久久久精品电影小说| 亚洲少妇的诱惑av| 久久精品aⅴ一区二区三区四区| 亚洲精品美女久久av网站| 国产在线免费精品| 欧美亚洲 丝袜 人妻 在线| 国产激情久久老熟女| av不卡在线播放| 人人妻,人人澡人人爽秒播 | 亚洲一码二码三码区别大吗| 国产亚洲欧美精品永久| www日本在线高清视频| 女人高潮潮喷娇喘18禁视频| 亚洲av电影在线进入| 69精品国产乱码久久久| 国产1区2区3区精品| 妹子高潮喷水视频| 在线天堂中文资源库| 久久精品亚洲av国产电影网| 观看美女的网站| 一级黄片播放器| 午夜福利乱码中文字幕| 亚洲一区二区三区欧美精品| 国产一区二区 视频在线| 亚洲精品成人av观看孕妇| 成年动漫av网址| 欧美日韩亚洲综合一区二区三区_| 久久免费观看电影| 亚洲av日韩在线播放| 黄频高清免费视频| 美女大奶头黄色视频| 在线观看免费日韩欧美大片| 欧美日韩一级在线毛片| 精品一区二区三区四区五区乱码 | 亚洲国产最新在线播放| 伊人久久国产一区二区| 黄频高清免费视频| 亚洲av男天堂| 亚洲国产精品一区三区| 亚洲国产欧美日韩在线播放| 亚洲成色77777| 人人妻人人澡人人看| 国产av国产精品国产| 国产精品久久久久久精品古装| 尾随美女入室| av福利片在线| netflix在线观看网站| 精品亚洲乱码少妇综合久久| 男女之事视频高清在线观看 | 午夜激情久久久久久久| av不卡在线播放| 搡老岳熟女国产| 波野结衣二区三区在线| 亚洲精品aⅴ在线观看| 咕卡用的链子| 国产激情久久老熟女| 精品久久久久久电影网| 最近最新中文字幕大全免费视频 | a级毛片黄视频| 99久久综合免费| 国产av国产精品国产| 日韩av免费高清视频| 国产欧美日韩综合在线一区二区| 国产高清不卡午夜福利| 少妇被粗大的猛进出69影院| 飞空精品影院首页| 国产成人a∨麻豆精品| 男女高潮啪啪啪动态图| 狂野欧美激情性xxxx| 久久久国产一区二区| 久久综合国产亚洲精品| 久久人妻熟女aⅴ| 免费观看av网站的网址| 国产 一区精品| 国产日韩一区二区三区精品不卡| 中文字幕制服av| 美女高潮到喷水免费观看| 久久久久精品久久久久真实原创| 色播在线永久视频| 日本欧美视频一区| 十八禁人妻一区二区| 亚洲国产精品一区三区| 99久久人妻综合| 久久国产精品大桥未久av| 亚洲国产av新网站| 国产精品嫩草影院av在线观看| 午夜影院在线不卡| 久久久久久久久久久久大奶| 欧美日韩视频高清一区二区三区二| 热re99久久精品国产66热6| 亚洲久久久国产精品| 亚洲色图 男人天堂 中文字幕| 免费观看性生交大片5| 18禁动态无遮挡网站| 亚洲精品美女久久av网站| 操出白浆在线播放| 国产精品熟女久久久久浪| h视频一区二区三区| 老鸭窝网址在线观看| 亚洲欧洲精品一区二区精品久久久 | 久久久久精品国产欧美久久久 | 多毛熟女@视频| 一边摸一边抽搐一进一出视频| 成年动漫av网址| 亚洲一区二区三区欧美精品| 久久韩国三级中文字幕| 国产亚洲精品第一综合不卡| 交换朋友夫妻互换小说| 午夜免费鲁丝| 久久久久国产精品人妻一区二区| 国产在线一区二区三区精| 久久精品aⅴ一区二区三区四区| 亚洲久久久国产精品| 亚洲中文av在线| 麻豆av在线久日| 国产精品秋霞免费鲁丝片| 多毛熟女@视频| 王馨瑶露胸无遮挡在线观看| 婷婷色av中文字幕| 国产成人欧美| 久久国产精品男人的天堂亚洲| 在线天堂中文资源库| 国产成人精品久久久久久| 大香蕉久久成人网| 精品一区在线观看国产| 国产精品蜜桃在线观看| 久久人妻熟女aⅴ| 成人三级做爰电影| 亚洲欧美一区二区三区国产| 久久av网站| 中文字幕人妻熟女乱码| 精品少妇内射三级| av在线老鸭窝| 91aial.com中文字幕在线观看| 丁香六月欧美| 美女中出高潮动态图| 亚洲人成电影观看| 久久国产亚洲av麻豆专区| 无遮挡黄片免费观看| 中文天堂在线官网| 精品第一国产精品| tube8黄色片| 精品免费久久久久久久清纯 | 丝袜人妻中文字幕| 国产精品 欧美亚洲| 精品少妇内射三级| 日韩制服丝袜自拍偷拍| 国产一卡二卡三卡精品 | 桃花免费在线播放| 热re99久久国产66热| 久久人人爽人人片av| 91成人精品电影| 欧美黄色片欧美黄色片| 亚洲七黄色美女视频| 欧美日韩福利视频一区二区| 人人澡人人妻人| 亚洲精品日本国产第一区| 国产在线一区二区三区精| 高清视频免费观看一区二区| 亚洲av欧美aⅴ国产| 操美女的视频在线观看| 亚洲精品av麻豆狂野| 亚洲熟女毛片儿| 久久精品人人爽人人爽视色| 中文字幕最新亚洲高清| 午夜福利视频在线观看免费| 亚洲精品视频女| 我的亚洲天堂| 日韩精品免费视频一区二区三区| 精品国产一区二区久久| 亚洲欧美激情在线| 国产福利在线免费观看视频| 午夜福利视频精品| 久久亚洲国产成人精品v| 国产麻豆69| 国产精品久久久久久精品电影小说| 十八禁网站网址无遮挡| 亚洲国产欧美在线一区| 中文字幕亚洲精品专区| 亚洲一区二区三区欧美精品| 黄片播放在线免费| 国产成人欧美| 美女高潮到喷水免费观看| 午夜av观看不卡| 欧美精品av麻豆av| 99久国产av精品国产电影| 久久精品国产a三级三级三级| 成人免费观看视频高清| 国产成人啪精品午夜网站| 欧美 日韩 精品 国产| 高清黄色对白视频在线免费看| 久久人人爽av亚洲精品天堂| 秋霞伦理黄片| 亚洲精品乱久久久久久| 成年av动漫网址| 欧美在线一区亚洲| netflix在线观看网站| 免费看不卡的av| 免费不卡黄色视频| 如何舔出高潮| 1024视频免费在线观看| 久久精品亚洲熟妇少妇任你| 大香蕉久久网| 中国国产av一级| 亚洲av在线观看美女高潮| 亚洲av中文av极速乱| 9191精品国产免费久久| 免费少妇av软件| 国产精品久久久久久人妻精品电影 | 这个男人来自地球电影免费观看 | 少妇精品久久久久久久| videosex国产| 国产精品秋霞免费鲁丝片| 男男h啪啪无遮挡| 中文字幕亚洲精品专区| 国产精品香港三级国产av潘金莲 | 亚洲 欧美一区二区三区| 国产淫语在线视频| 操美女的视频在线观看| 狠狠精品人妻久久久久久综合| 欧美在线一区亚洲| 一本大道久久a久久精品| 人妻一区二区av| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产欧美日韩在线播放| 亚洲人成电影观看| 九色亚洲精品在线播放| 亚洲av电影在线进入| 制服人妻中文乱码| 男人添女人高潮全过程视频| 最黄视频免费看| 精品少妇内射三级| 国产爽快片一区二区三区| 女性生殖器流出的白浆| 成人国产麻豆网| 毛片一级片免费看久久久久| www.自偷自拍.com| 最近中文字幕2019免费版| 亚洲国产中文字幕在线视频| 午夜影院在线不卡| 女人精品久久久久毛片| 99久久人妻综合| 国产黄频视频在线观看| 国产精品一区二区精品视频观看| 天堂中文最新版在线下载| 亚洲色图综合在线观看| 久久av网站| 久久久国产一区二区| 亚洲av成人不卡在线观看播放网 | 制服丝袜香蕉在线| netflix在线观看网站| 国产精品二区激情视频| 伦理电影免费视频| 两个人免费观看高清视频| 少妇人妻久久综合中文| 中文字幕制服av| 十八禁人妻一区二区| 欧美国产精品va在线观看不卡| 国产成人精品福利久久| 日本av手机在线免费观看| 在线观看免费视频网站a站| 菩萨蛮人人尽说江南好唐韦庄| 天堂俺去俺来也www色官网| 满18在线观看网站| 99久国产av精品国产电影| 亚洲成色77777| 国产成人啪精品午夜网站| 热re99久久精品国产66热6| 色94色欧美一区二区| 性少妇av在线| 亚洲国产精品成人久久小说| 国产精品国产av在线观看| 69精品国产乱码久久久| 一级毛片我不卡| 美女中出高潮动态图| 少妇人妻久久综合中文| 国产深夜福利视频在线观看| 观看av在线不卡| 麻豆乱淫一区二区| 日韩av不卡免费在线播放| 我要看黄色一级片免费的| 国产在线免费精品| 高清不卡的av网站| 久久久久精品人妻al黑| 天美传媒精品一区二区| 欧美精品高潮呻吟av久久| 国产精品偷伦视频观看了| av在线播放精品| 高清在线视频一区二区三区| 国产福利在线免费观看视频| 亚洲av电影在线观看一区二区三区| 黄色 视频免费看| 欧美最新免费一区二区三区| 亚洲精品国产一区二区精华液| 精品人妻一区二区三区麻豆| 国产成人91sexporn| 女人被躁到高潮嗷嗷叫费观| 99香蕉大伊视频| 黄网站色视频无遮挡免费观看| 老鸭窝网址在线观看| 人妻人人澡人人爽人人| 国产熟女欧美一区二区| 99久久99久久久精品蜜桃| 国产xxxxx性猛交| 别揉我奶头~嗯~啊~动态视频 | 国产精品久久久久成人av| 成人免费观看视频高清| 啦啦啦在线免费观看视频4| 亚洲成av片中文字幕在线观看| 久久午夜综合久久蜜桃| 午夜福利,免费看| 国产人伦9x9x在线观看| 免费看不卡的av| 黑人猛操日本美女一级片| 国产精品香港三级国产av潘金莲 | 国产精品嫩草影院av在线观看| 国产有黄有色有爽视频| 亚洲自偷自拍图片 自拍| 亚洲av电影在线进入| 久久 成人 亚洲| 在线天堂中文资源库| 黄色怎么调成土黄色| 人妻一区二区av| 久久久久人妻精品一区果冻| kizo精华| 一级片'在线观看视频| 黄网站色视频无遮挡免费观看| 国产成人精品福利久久| 大码成人一级视频| 97精品久久久久久久久久精品| 中文字幕人妻丝袜一区二区 | 新久久久久国产一级毛片| 这个男人来自地球电影免费观看 | 久久精品人人爽人人爽视色| 亚洲av福利一区| 亚洲精品在线美女| 欧美日韩综合久久久久久| 搡老乐熟女国产| 日韩制服丝袜自拍偷拍| 国产精品 国内视频| 久久 成人 亚洲| 国产日韩欧美在线精品| 菩萨蛮人人尽说江南好唐韦庄| 日韩制服丝袜自拍偷拍| 亚洲免费av在线视频| 久久人人97超碰香蕉20202| 男女边摸边吃奶| 午夜激情久久久久久久| 交换朋友夫妻互换小说| 一本—道久久a久久精品蜜桃钙片| 亚洲精品国产av成人精品| 久久毛片免费看一区二区三区| 一级毛片电影观看| 韩国av在线不卡| 人人妻人人澡人人爽人人夜夜| 又粗又硬又长又爽又黄的视频| 亚洲欧美激情在线| 免费高清在线观看日韩| 国产高清国产精品国产三级| 亚洲国产精品国产精品| 午夜日本视频在线| 日韩欧美精品免费久久| a级片在线免费高清观看视频| 最新在线观看一区二区三区 | 国产乱来视频区| 国产极品天堂在线| 亚洲精品国产区一区二| 欧美日韩一级在线毛片| 亚洲少妇的诱惑av| 国产男人的电影天堂91| 热99久久久久精品小说推荐| 两个人看的免费小视频| 国产在线视频一区二区| 精品国产国语对白av| 国产熟女欧美一区二区| 老司机影院成人| 美女大奶头黄色视频| 亚洲五月色婷婷综合| 十八禁人妻一区二区| 亚洲av综合色区一区| 免费观看a级毛片全部| 中文字幕亚洲精品专区| 亚洲av男天堂| 久久精品熟女亚洲av麻豆精品| 欧美国产精品一级二级三级| 午夜福利,免费看| 亚洲av电影在线观看一区二区三区| 成年动漫av网址| 国产精品亚洲av一区麻豆 | 久久久亚洲精品成人影院| 人人妻人人添人人爽欧美一区卜| 亚洲国产毛片av蜜桃av| 日本色播在线视频| 亚洲精品日本国产第一区| 国产成人精品福利久久| 搡老乐熟女国产| 男人添女人高潮全过程视频| 日韩精品有码人妻一区| 天堂俺去俺来也www色官网| 在线看a的网站| 99久久99久久久精品蜜桃| 亚洲在久久综合| 成人影院久久| 国产又爽黄色视频| 女人精品久久久久毛片| 中文字幕亚洲精品专区| 又大又爽又粗| 咕卡用的链子| 欧美激情极品国产一区二区三区| 久久ye,这里只有精品| 天天影视国产精品| 少妇人妻久久综合中文| 如何舔出高潮| 久久性视频一级片| 亚洲精品国产区一区二| 菩萨蛮人人尽说江南好唐韦庄| 少妇精品久久久久久久| 少妇人妻久久综合中文| 黄色 视频免费看| 久久免费观看电影| 免费观看a级毛片全部| 侵犯人妻中文字幕一二三四区| 这个男人来自地球电影免费观看 | 日韩一卡2卡3卡4卡2021年| 天天影视国产精品| 中国国产av一级| www.av在线官网国产| 天天躁日日躁夜夜躁夜夜| 永久免费av网站大全| 亚洲欧美日韩另类电影网站| 波多野结衣一区麻豆| 亚洲美女搞黄在线观看| 国产精品人妻久久久影院| 国产精品久久久人人做人人爽| 校园人妻丝袜中文字幕| 99精国产麻豆久久婷婷| 国产精品免费视频内射| 看非洲黑人一级黄片| 在线观看人妻少妇| 国产精品 欧美亚洲| 亚洲精品美女久久久久99蜜臀 | 如何舔出高潮| 各种免费的搞黄视频| 国产午夜精品一二区理论片| 丁香六月欧美| 久久精品国产亚洲av高清一级| 曰老女人黄片| 男女下面插进去视频免费观看| 无限看片的www在线观看| 看免费成人av毛片| 亚洲国产中文字幕在线视频| 少妇精品久久久久久久| 桃花免费在线播放| 国产欧美日韩一区二区三区在线| 久久99一区二区三区| 少妇猛男粗大的猛烈进出视频| 99re6热这里在线精品视频| 制服人妻中文乱码| 国产午夜精品一二区理论片| 亚洲国产中文字幕在线视频| 精品国产一区二区三区四区第35| 亚洲成国产人片在线观看| 国产在线一区二区三区精| 午夜激情久久久久久久| 亚洲精品美女久久久久99蜜臀 | 乱人伦中国视频| 熟女少妇亚洲综合色aaa.| 亚洲成av片中文字幕在线观看| 99九九在线精品视频| 电影成人av| 欧美日韩亚洲高清精品| 亚洲国产欧美在线一区| 99国产精品免费福利视频| 国产免费一区二区三区四区乱码| 老鸭窝网址在线观看| 飞空精品影院首页| 精品久久蜜臀av无| 欧美黄色片欧美黄色片| 一区二区三区乱码不卡18| 1024视频免费在线观看| 韩国精品一区二区三区| 欧美另类一区| 90打野战视频偷拍视频| 亚洲欧美精品综合一区二区三区| 国产极品粉嫩免费观看在线| 成年动漫av网址| 久久99热这里只频精品6学生| 狠狠婷婷综合久久久久久88av| 亚洲欧美激情在线| 操美女的视频在线观看| kizo精华| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品视频女| 丝瓜视频免费看黄片| 国产精品欧美亚洲77777| www.av在线官网国产| 国产片特级美女逼逼视频| 哪个播放器可以免费观看大片| 一级爰片在线观看| 日本爱情动作片www.在线观看| 精品免费久久久久久久清纯 | 午夜福利,免费看| 一区二区三区乱码不卡18| 日韩大码丰满熟妇| 极品少妇高潮喷水抽搐| 久久久久精品国产欧美久久久 | 人人妻,人人澡人人爽秒播 | 91精品伊人久久大香线蕉| 天天操日日干夜夜撸| 国产一级毛片在线| 色综合欧美亚洲国产小说| 欧美日本中文国产一区发布| a级毛片黄视频| 看非洲黑人一级黄片| 成年动漫av网址| 黄色毛片三级朝国网站| 天天躁夜夜躁狠狠久久av| 999久久久国产精品视频| 国产亚洲午夜精品一区二区久久| videosex国产| 精品少妇一区二区三区视频日本电影 | 久久影院123| 欧美日韩亚洲综合一区二区三区_| 麻豆精品久久久久久蜜桃| 桃花免费在线播放| 男女床上黄色一级片免费看| 亚洲欧美清纯卡通| 黄片无遮挡物在线观看| 日日啪夜夜爽| 青青草视频在线视频观看| 啦啦啦 在线观看视频| 在现免费观看毛片| 欧美精品一区二区免费开放| 亚洲欧美色中文字幕在线| 久久青草综合色| 欧美精品人与动牲交sv欧美| av线在线观看网站| av天堂久久9| 狠狠婷婷综合久久久久久88av| 蜜桃国产av成人99| 美女扒开内裤让男人捅视频| 午夜福利一区二区在线看| 91老司机精品| 不卡视频在线观看欧美| 精品一区二区三区av网在线观看 | 国产精品一二三区在线看| 国产精品偷伦视频观看了| 精品国产超薄肉色丝袜足j| 在线观看免费日韩欧美大片| 久久久久久人人人人人| 国产成人91sexporn| 精品一区二区三卡| 国产亚洲最大av| 日本黄色日本黄色录像| 少妇猛男粗大的猛烈进出视频| 在线天堂最新版资源| 国产老妇伦熟女老妇高清| 日韩av在线免费看完整版不卡| 久久精品国产a三级三级三级| 亚洲av电影在线观看一区二区三区| 亚洲欧美精品自产自拍| 国产爽快片一区二区三区| 久久ye,这里只有精品| 欧美xxⅹ黑人| 伊人久久大香线蕉亚洲五| 美女扒开内裤让男人捅视频| 色94色欧美一区二区| 久久久精品区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 精品一区在线观看国产| 亚洲精品美女久久久久99蜜臀 | 在线 av 中文字幕| 在线观看免费日韩欧美大片| 国产精品久久久久久久久免| 一区二区av电影网| 国产乱人偷精品视频| 国产男人的电影天堂91|