• <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| 午夜日本视频在线| 人妻一区二区av| 国产精品久久久久久久久免| 久久精品久久精品一区二区三区| 中文精品一卡2卡3卡4更新| 精品人妻偷拍中文字幕| 人体艺术视频欧美日本| 女性生殖器流出的白浆| 久久精品国产亚洲av涩爱| 亚洲精品久久成人aⅴ小说| 国产乱人偷精品视频| 侵犯人妻中文字幕一二三四区| 青青草视频在线视频观看| 欧美日韩一级在线毛片| 国产精品国产三级专区第一集| 人成视频在线观看免费观看| 你懂的网址亚洲精品在线观看| 国产男女超爽视频在线观看| 亚洲第一av免费看| 狠狠精品人妻久久久久久综合| 热re99久久精品国产66热6| 国产黄色视频一区二区在线观看| 日韩免费高清中文字幕av| 丝袜美足系列| 久久精品久久精品一区二区三区| 人人澡人人妻人| 精品一区在线观看国产| 亚洲欧洲日产国产| xxx大片免费视频| 一本—道久久a久久精品蜜桃钙片| 日韩精品有码人妻一区| 中文字幕人妻熟女乱码| 欧美日韩一级在线毛片| 汤姆久久久久久久影院中文字幕| 肉色欧美久久久久久久蜜桃| 老司机影院毛片| 90打野战视频偷拍视频| 岛国毛片在线播放| 国产成人精品无人区| 国产欧美日韩综合在线一区二区| 国产人伦9x9x在线观看 | 久久精品国产a三级三级三级| 18禁动态无遮挡网站| 熟妇人妻不卡中文字幕| 精品酒店卫生间| 亚洲精品视频女| 极品人妻少妇av视频| 两个人看的免费小视频| 少妇熟女欧美另类| 少妇的逼水好多| 美女脱内裤让男人舔精品视频| 成人手机av| 亚洲成国产人片在线观看| 日本爱情动作片www.在线观看| 久久久久网色| 国产精品嫩草影院av在线观看| 一本久久精品| 丰满少妇做爰视频| 亚洲国产最新在线播放| 两个人看的免费小视频| 亚洲色图综合在线观看| 成人影院久久| 久久影院123| 丰满迷人的少妇在线观看| 久久午夜综合久久蜜桃| 国产探花极品一区二区| videosex国产| 五月伊人婷婷丁香| 9色porny在线观看| 亚洲四区av| 韩国av在线不卡| 日日啪夜夜爽| videosex国产| 人人妻人人澡人人看| 国产日韩一区二区三区精品不卡| 日韩成人av中文字幕在线观看| 丝袜脚勾引网站| 极品少妇高潮喷水抽搐| 欧美激情极品国产一区二区三区| 成人手机av| 2018国产大陆天天弄谢| 丰满少妇做爰视频| 丁香六月天网| 熟女电影av网| 日韩熟女老妇一区二区性免费视频| 国产黄频视频在线观看| 久久久久人妻精品一区果冻| 国产精品久久久久久av不卡| 高清欧美精品videossex| 综合色丁香网| 在线观看免费高清a一片| 久久韩国三级中文字幕| 91精品国产国语对白视频| 久久精品国产综合久久久| 嫩草影院入口| 日韩免费高清中文字幕av| 99久国产av精品国产电影| 水蜜桃什么品种好| 天堂俺去俺来也www色官网| 一本久久精品| 嫩草影院入口| 春色校园在线视频观看| 夜夜骑夜夜射夜夜干| 国语对白做爰xxxⅹ性视频网站| 亚洲色图综合在线观看| 黑人欧美特级aaaaaa片| 人人澡人人妻人| 黑人猛操日本美女一级片| 中文字幕最新亚洲高清| 又黄又粗又硬又大视频| 成人二区视频| 成年动漫av网址| 亚洲成色77777| 久久精品亚洲av国产电影网| 国产成人91sexporn| 国产男人的电影天堂91| 国产男人的电影天堂91| 亚洲欧洲精品一区二区精品久久久 | 久久久久久久国产电影| 建设人人有责人人尽责人人享有的| 97人妻天天添夜夜摸| 97精品久久久久久久久久精品| 亚洲欧美清纯卡通| 国产精品99久久99久久久不卡 | 精品国产乱码久久久久久男人| 亚洲欧美精品自产自拍| 精品人妻偷拍中文字幕| 久久精品亚洲av国产电影网| 亚洲色图 男人天堂 中文字幕| 日本wwww免费看| 亚洲一级一片aⅴ在线观看| 色哟哟·www| www.精华液| 久久热在线av| 少妇猛男粗大的猛烈进出视频| 精品国产国语对白av| 美女国产高潮福利片在线看| 欧美精品亚洲一区二区| 免费黄色在线免费观看| 韩国高清视频一区二区三区| 亚洲精品国产av蜜桃| 免费黄频网站在线观看国产| 国产av一区二区精品久久| 国产精品香港三级国产av潘金莲 | 久久久久久久久久人人人人人人| 人体艺术视频欧美日本| 黑丝袜美女国产一区| 男女啪啪激烈高潮av片| 久久精品国产自在天天线| 一级毛片 在线播放| 在线观看三级黄色| 国产亚洲欧美精品永久| 亚洲av电影在线观看一区二区三区| 亚洲精品,欧美精品| 精品一区二区免费观看| 中文精品一卡2卡3卡4更新| 久久热在线av| 国产免费一区二区三区四区乱码| 亚洲欧美中文字幕日韩二区| 亚洲精品在线美女| 国产精品久久久久久精品电影小说| 国产精品免费视频内射| 亚洲,一卡二卡三卡| 啦啦啦中文免费视频观看日本| 好男人视频免费观看在线| 午夜91福利影院| 五月天丁香电影| av片东京热男人的天堂| av一本久久久久| 制服人妻中文乱码| 亚洲,一卡二卡三卡| 下体分泌物呈黄色| 下体分泌物呈黄色| 国产无遮挡羞羞视频在线观看| 麻豆精品久久久久久蜜桃| 在线亚洲精品国产二区图片欧美| 天天影视国产精品| 亚洲视频免费观看视频| 夫妻性生交免费视频一级片| 国产女主播在线喷水免费视频网站| 91国产中文字幕| 亚洲av.av天堂| 夜夜骑夜夜射夜夜干| 侵犯人妻中文字幕一二三四区| 欧美av亚洲av综合av国产av | 在线天堂中文资源库| 亚洲精品,欧美精品| 中文字幕最新亚洲高清| 国产精品免费视频内射| 久久精品aⅴ一区二区三区四区 | 男女国产视频网站| 国产成人aa在线观看| 国产精品香港三级国产av潘金莲 | 97精品久久久久久久久久精品| 亚洲国产看品久久| 国产成人精品久久二区二区91 | 69精品国产乱码久久久| 日本爱情动作片www.在线观看| 建设人人有责人人尽责人人享有的| 日韩精品免费视频一区二区三区| 午夜福利在线观看免费完整高清在| 99热网站在线观看| 午夜91福利影院| 91成人精品电影| 黄频高清免费视频| 久久久久久久久免费视频了| 天天躁狠狠躁夜夜躁狠狠躁| 五月伊人婷婷丁香| 欧美成人午夜精品| 交换朋友夫妻互换小说| 99久久综合免费| 国产一区有黄有色的免费视频| 欧美亚洲 丝袜 人妻 在线| 丝袜美足系列| 国产在线视频一区二区| 男女国产视频网站| 欧美97在线视频| 99香蕉大伊视频| 人妻少妇偷人精品九色| 国产欧美亚洲国产| 日本vs欧美在线观看视频| 9191精品国产免费久久| 精品99又大又爽又粗少妇毛片| 国产欧美亚洲国产| 99热国产这里只有精品6| 热re99久久国产66热| 国产深夜福利视频在线观看| 最近中文字幕2019免费版| 美女国产视频在线观看| 日韩免费高清中文字幕av| 亚洲国产毛片av蜜桃av| 久久99一区二区三区| 欧美精品国产亚洲| 老司机亚洲免费影院| 久久精品熟女亚洲av麻豆精品| 欧美精品一区二区大全| 日韩一卡2卡3卡4卡2021年| 久久青草综合色| 欧美成人精品欧美一级黄| 欧美黄色片欧美黄色片| 亚洲av国产av综合av卡| 国产人伦9x9x在线观看 | 婷婷色麻豆天堂久久| 搡老乐熟女国产| 国产精品国产三级国产专区5o| 秋霞伦理黄片| 亚洲精品视频女| 黄片播放在线免费| 国产精品蜜桃在线观看| 中文字幕人妻丝袜一区二区 | 久久这里只有精品19| 少妇猛男粗大的猛烈进出视频| 日韩制服丝袜自拍偷拍| 久久久久久久久久人人人人人人| 一本色道久久久久久精品综合| 免费日韩欧美在线观看| 成年女人在线观看亚洲视频| 国产精品 国内视频| 伦理电影免费视频| 亚洲成人一二三区av| 美女高潮到喷水免费观看| 精品久久久久久电影网| 免费在线观看视频国产中文字幕亚洲 | 多毛熟女@视频| 午夜免费鲁丝| 亚洲成人一二三区av| 五月开心婷婷网| 国产精品亚洲av一区麻豆 | 中文乱码字字幕精品一区二区三区| 欧美变态另类bdsm刘玥| 如日韩欧美国产精品一区二区三区| 日韩成人av中文字幕在线观看| 婷婷色av中文字幕| 精品酒店卫生间| 精品国产一区二区久久| 一区二区三区乱码不卡18| 国产色婷婷99| 国产黄色免费在线视频| 国产成人av激情在线播放| 纯流量卡能插随身wifi吗| 日本-黄色视频高清免费观看| 国产精品亚洲av一区麻豆 | 亚洲欧美精品综合一区二区三区 | 在线观看免费日韩欧美大片| 天堂8中文在线网| 国产爽快片一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 久久久久国产一级毛片高清牌| 水蜜桃什么品种好| 久久久精品免费免费高清| 高清黄色对白视频在线免费看| 国产熟女午夜一区二区三区| 有码 亚洲区| 国产一区二区 视频在线| 国产综合精华液| videossex国产| 亚洲av免费高清在线观看| 久久久久久久亚洲中文字幕| 成人亚洲精品一区在线观看| 哪个播放器可以免费观看大片| av天堂久久9| 欧美成人精品欧美一级黄| 亚洲男人天堂网一区| 亚洲av中文av极速乱| 欧美日韩亚洲国产一区二区在线观看 | 老司机影院毛片| 精品一区二区免费观看| 亚洲精品av麻豆狂野| 成年人午夜在线观看视频| 国产精品蜜桃在线观看| 大香蕉久久网| 五月伊人婷婷丁香| 久久精品久久精品一区二区三区| 国产无遮挡羞羞视频在线观看| 黄片小视频在线播放| 国产极品天堂在线| 国产日韩一区二区三区精品不卡| 亚洲精品第二区| 精品人妻在线不人妻| av网站免费在线观看视频| 一边亲一边摸免费视频| 欧美日韩精品网址| 久久av网站| 街头女战士在线观看网站| 久久久精品94久久精品| 9热在线视频观看99| 精品国产一区二区三区四区第35| 精品国产乱码久久久久久男人| 一本色道久久久久久精品综合| 欧美 亚洲 国产 日韩一| 免费黄频网站在线观看国产| 熟女电影av网| 国产精品亚洲av一区麻豆 | 美国免费a级毛片| 精品亚洲成国产av| 黑丝袜美女国产一区| 99久久中文字幕三级久久日本| 亚洲精品国产av成人精品| 香蕉精品网在线| 三级国产精品片| 午夜福利网站1000一区二区三区| 亚洲综合精品二区| 中文字幕人妻丝袜制服| 咕卡用的链子| 国产高清不卡午夜福利| 久久午夜福利片| www.av在线官网国产| 校园人妻丝袜中文字幕| 在线天堂中文资源库| av天堂久久9| 午夜免费男女啪啪视频观看| 国产男女内射视频| 午夜免费观看性视频| 又粗又硬又长又爽又黄的视频| 91成人精品电影| 亚洲伊人久久精品综合| 新久久久久国产一级毛片| 极品人妻少妇av视频| 精品99又大又爽又粗少妇毛片| 久久精品国产自在天天线| 亚洲av综合色区一区| 丝袜人妻中文字幕| 一级黄片播放器| 亚洲av国产av综合av卡| 亚洲国产色片| 精品人妻偷拍中文字幕| 一二三四在线观看免费中文在| 成人毛片60女人毛片免费| 欧美 亚洲 国产 日韩一| 肉色欧美久久久久久久蜜桃| 女性被躁到高潮视频| 最近最新中文字幕大全免费视频 | 日韩大片免费观看网站| 国产成人精品一,二区| 中文字幕色久视频| 日本欧美视频一区| 国产亚洲欧美精品永久| 一区二区日韩欧美中文字幕| 男女高潮啪啪啪动态图| 尾随美女入室| 亚洲精品av麻豆狂野| 岛国毛片在线播放| videossex国产| 精品国产一区二区三区久久久樱花| 欧美日韩精品网址| 亚洲精品自拍成人| 亚洲精品在线美女| 久久精品久久久久久噜噜老黄| 天天躁夜夜躁狠狠躁躁| 亚洲av男天堂| 亚洲国产看品久久| 国产片内射在线| 亚洲欧美成人精品一区二区| xxxhd国产人妻xxx| 久久女婷五月综合色啪小说| 最近手机中文字幕大全| 免费看不卡的av| 欧美日韩精品成人综合77777| 国产免费视频播放在线视频| 热re99久久精品国产66热6| 伦精品一区二区三区| 亚洲欧美成人综合另类久久久| 高清欧美精品videossex| 午夜免费男女啪啪视频观看| 免费在线观看黄色视频的| 在现免费观看毛片| 黄片小视频在线播放| 少妇人妻久久综合中文| 97人妻天天添夜夜摸| 精品酒店卫生间| 国产精品三级大全| 成年av动漫网址| 男的添女的下面高潮视频| 亚洲,一卡二卡三卡| 国产精品香港三级国产av潘金莲 | 日韩人妻精品一区2区三区| 国产成人aa在线观看| 日产精品乱码卡一卡2卡三| 久久久精品区二区三区| 中文字幕人妻熟女乱码| 自线自在国产av| 日韩制服骚丝袜av| 国产一区有黄有色的免费视频| 久久久国产一区二区| 成人午夜精彩视频在线观看| 成人漫画全彩无遮挡| 91精品三级在线观看| 亚洲精品久久午夜乱码| 黄色一级大片看看| 丰满少妇做爰视频| av网站免费在线观看视频| 日韩一区二区三区影片| 婷婷色综合大香蕉| 欧美成人午夜免费资源| 国产男人的电影天堂91| 宅男免费午夜| 中文字幕av电影在线播放| 亚洲成人av在线免费| 宅男免费午夜| 五月开心婷婷网| 毛片一级片免费看久久久久| 精品国产乱码久久久久久男人| 亚洲av免费高清在线观看| 精品福利永久在线观看| 99久国产av精品国产电影| 18禁国产床啪视频网站| 精品久久久久久电影网| 亚洲av综合色区一区| 国产极品天堂在线| 久久久久久人妻| 成人国语在线视频| 另类亚洲欧美激情| 黄色毛片三级朝国网站| 亚洲欧美一区二区三区黑人 | 狠狠婷婷综合久久久久久88av| a 毛片基地| 国产黄频视频在线观看| 亚洲图色成人| 大片电影免费在线观看免费| 亚洲久久久国产精品| av卡一久久| 国产福利在线免费观看视频| 另类精品久久| 亚洲综合色惰| 人妻 亚洲 视频| 国产综合精华液| 老司机影院毛片| a 毛片基地| 国产精品一国产av| av在线app专区| 一边亲一边摸免费视频| 国产伦理片在线播放av一区| 99久久综合免费| 大话2 男鬼变身卡| 婷婷色av中文字幕| 18禁裸乳无遮挡动漫免费视频| 亚洲人成网站在线观看播放| 卡戴珊不雅视频在线播放| 日本免费在线观看一区| 亚洲三区欧美一区| 亚洲综合色惰| 国产不卡av网站在线观看| 一二三四中文在线观看免费高清| 亚洲av日韩在线播放| 激情视频va一区二区三区| 久久久精品国产亚洲av高清涩受| 久久精品熟女亚洲av麻豆精品| 色视频在线一区二区三区| 熟女电影av网| 建设人人有责人人尽责人人享有的| 一区二区av电影网| 成人午夜精彩视频在线观看| 国产av国产精品国产| 大话2 男鬼变身卡| 丝袜在线中文字幕| 色吧在线观看| 两个人免费观看高清视频| 日本爱情动作片www.在线观看| 中文字幕另类日韩欧美亚洲嫩草| 久久久久国产精品人妻一区二区| 国产精品国产三级专区第一集| 97精品久久久久久久久久精品| 成年美女黄网站色视频大全免费| 精品卡一卡二卡四卡免费| 精品国产一区二区三区久久久樱花| 人妻人人澡人人爽人人| 亚洲国产精品一区二区三区在线| 亚洲精品在线美女| 各种免费的搞黄视频| 18+在线观看网站| 午夜日韩欧美国产| 日韩人妻精品一区2区三区| 亚洲成人一二三区av| 18禁裸乳无遮挡动漫免费视频| 在线观看免费视频网站a站| 久久精品久久精品一区二区三区| 大话2 男鬼变身卡| 中文乱码字字幕精品一区二区三区| 日本-黄色视频高清免费观看| 99久久综合免费| 亚洲久久久国产精品| 成人漫画全彩无遮挡| 欧美变态另类bdsm刘玥| www.熟女人妻精品国产| 久久精品aⅴ一区二区三区四区 | 毛片一级片免费看久久久久| 伦理电影大哥的女人| 啦啦啦在线免费观看视频4| 另类亚洲欧美激情| 国产精品三级大全| 日韩精品免费视频一区二区三区| 欧美日韩亚洲国产一区二区在线观看 | 色哟哟·www| 人妻 亚洲 视频| 亚洲欧洲日产国产| 久热这里只有精品99| 免费久久久久久久精品成人欧美视频| 91成人精品电影| 亚洲伊人久久精品综合| 国产在线一区二区三区精| 日韩精品有码人妻一区| 另类精品久久| 欧美精品一区二区免费开放| 啦啦啦啦在线视频资源| 18在线观看网站| 99热全是精品| 久久精品亚洲av国产电影网| 三级国产精品片| 亚洲av中文av极速乱| 欧美日韩av久久| 一区在线观看完整版| 亚洲国产看品久久| 国产av国产精品国产| 一级片免费观看大全| 亚洲国产成人一精品久久久| 男女边吃奶边做爰视频| 亚洲成色77777| 人妻系列 视频| 亚洲天堂av无毛| 国产日韩欧美在线精品| 啦啦啦视频在线资源免费观看| 国产免费福利视频在线观看| 女人被躁到高潮嗷嗷叫费观| 中文字幕人妻丝袜一区二区 | 蜜桃国产av成人99| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲第一av免费看| 国产成人免费无遮挡视频| 成人毛片60女人毛片免费| 自线自在国产av| √禁漫天堂资源中文www| 成人亚洲欧美一区二区av| 人妻系列 视频| 大码成人一级视频| 国产在线视频一区二区| 最近2019中文字幕mv第一页| 亚洲欧美成人综合另类久久久| 少妇的逼水好多| 春色校园在线视频观看| 亚洲天堂av无毛| 80岁老熟妇乱子伦牲交| 精品少妇内射三级| 国产精品国产三级专区第一集| 女人高潮潮喷娇喘18禁视频| 天美传媒精品一区二区| 大话2 男鬼变身卡| 韩国精品一区二区三区| 久久久久人妻精品一区果冻| 国产成人精品久久久久久| 熟妇人妻不卡中文字幕| 国产亚洲最大av| 国产熟女欧美一区二区| 国产日韩欧美视频二区| 欧美精品国产亚洲| 成人免费观看视频高清| 天堂俺去俺来也www色官网| 七月丁香在线播放| 看免费成人av毛片| 亚洲一区中文字幕在线| 精品一区二区免费观看| 高清视频免费观看一区二区| a 毛片基地| 成人国语在线视频| 在线精品无人区一区二区三| 九色亚洲精品在线播放| 精品国产一区二区久久| 亚洲成av片中文字幕在线观看 |