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

    1880—2020年安徽蕪湖氣溫長序列構(gòu)建及年代際特征

    2022-11-16 04:40:38孫大兵張曉憶
    干旱氣象 2022年5期

    劉 蕾,李 鸞,張 麗,孫大兵,張曉憶

    (安徽省蕪湖市氣象局,安徽 蕪湖 241000)

    引言

    近1000 a中,20世紀(jì)以來氣候變暖最為顯著,尤其近50 a的變暖趨勢更為明顯[1]。研究發(fā)現(xiàn),中國百年氣溫的暖期分別出現(xiàn)在20世紀(jì)30—40年代和80—90年代,但存在明顯的區(qū)域性差異[2],早期增暖過程中氣候系統(tǒng)內(nèi)部過程影響較大,而近幾十年的增溫則更多受外強迫的影響[3]。其中,中國西北、東北、華北地區(qū)增暖最顯著,而西南地區(qū)增溫趨勢最弱[4]。由于資料來源、觀測時段不同,中國近百年增溫速率仍存在爭議[5-7],總體上百年增溫范圍在0.90~1.52℃之間[8]。完整的百年氣溫長序列是氣候變化分析的基礎(chǔ),鑒于1951年前全國氣象觀測臺站稀少,氣溫序列缺測較多,最初有學(xué)者利用氣溫等級對氣溫序列進行延長插補[9],或利用史料、樹木年輪等代用資料來完善,而近年來利用重建的再分析資料研究氣溫的百年變化成為新態(tài)勢[5-6,10],但這些代用資料、再分析資料普遍存在變化幅度小、極值不顯著的特征,無法準(zhǔn)確反映極端年份的氣溫變化,因此利用百年以上尺度的實測氣溫數(shù)據(jù)進行氣候變化分析至關(guān)重要。研究結(jié)果顯示,不同地區(qū)局地百年氣溫變化既有一定的共性,也存在一定的差異。如,廣州在20世紀(jì)40年代開始增溫,80年代至20世紀(jì)末進入快速增溫時期[11];大連在1930年和1982年存在2個增暖突變點,20世紀(jì)90年代之后增暖趨勢最顯著[12];上海在20世紀(jì)30—40年代、90年代至21世紀(jì)初升溫明顯[13];浙江20世紀(jì)30—40年代增溫不顯著,但在20世紀(jì)70年代之后增溫迅速[14]??梢姡瑥V州、大連、上海等城市的暖期與全國基本一致,而浙江與全國存在一定差異。

    長江中下游地區(qū)僅有南京、武漢、蕪湖、上海4個WMO(世界氣象組織)百年氣象站,但這些站的百年氣溫變化研究時間偏早,且在氣溫序列重建中對觀測時次導(dǎo)致的誤差、臺站搬遷、數(shù)據(jù)均一性檢驗等考慮較少[15-16],無法得到真實反映。蕪湖站作為我國14個世界級百年氣象站之一,同時也是安徽省唯一的世界級百年氣象站,其意義對于安徽省乃至中國非常重大。然而,目前蕪湖站近140 a的氣溫數(shù)據(jù)仍未數(shù)字化整理和質(zhì)量控制。為此,本文參考前人研究方法,對蕪湖站2個時段(1880—1937年、1952—2020年)月平均氣溫進行數(shù)字化整理及質(zhì)量控制,并結(jié)合英國東英吉利大學(xué)氣候研究中心(Climatic Research Unit,CRU)格點氣溫資料,采用逐步回歸分析方法,對觀測氣溫進行擬合并對缺測值(1938—1951年)插補,進而構(gòu)建1880—2020年百年月平均氣溫序列,在此基礎(chǔ)上探討氣溫的年代際變化特征。

    1 資料及來源

    1880年3月1日起,蕪湖海關(guān)氣象觀測站正式開始觀測,1937年11月蕪湖淪陷,氣象觀測被迫中止,建國后1952年1月11日開始恢復(fù)觀測。除1938—1951年、1952年1月部分數(shù)據(jù)缺測外,其他時段數(shù)據(jù)完整、可用性高。其中,1880—1937年日平均氣溫來源于手抄的“海關(guān)氣象觀測月總薄”,1952—2020年氣溫資料來源于安徽省氣象信息中心。另外,還使用了英國東英吉利大學(xué)氣候研究中心(CRU)發(fā)布的1901—2020年TS(time series)4.04格點氣溫逐月資料,空間分辨率為0.5°×0.5°,該數(shù)據(jù)在中國長江流域適用性強[17]。另外,安徽省行政邊界是基于安徽省自然資源廳審核批準(zhǔn)的審圖號為皖S(2021)3號的標(biāo)準(zhǔn)地圖制作,底圖無修改。

    首先,對蕪湖站氣溫資料進行質(zhì)量控制,即設(shè)定各月日平均氣溫3倍標(biāo)準(zhǔn)差為閾值,對1880—1937年、1952—2020年日平均氣溫進行質(zhì)量檢驗,若日平均氣溫超出閾值則記為缺測,氣溫缺測5 d以上的月份記為月平均氣溫缺測??傮w來看,蕪湖站月平均氣溫質(zhì)量較高,僅1923年10、11月和1937年12月、1952年1月平均氣溫缺測。將缺測的月平均氣溫暫用距離最近的CRU格點氣溫代替。

    2 蕪湖站月平均氣溫訂正和序列構(gòu)建

    2.1 1880—1937年和1952—2020年氣溫訂正及插補

    2.1.1 差值訂正

    1880—1937年氣溫觀測時次與建國后不同,由于觀測時次、計算方法不同,平均氣溫往往存在細微差別,需要將其訂正到統(tǒng)一標(biāo)準(zhǔn)上[14,18]。根據(jù)《地面氣象觀測規(guī)范》要求[19],蕪湖站1952—2020年平均氣溫計算采用02:00、08:00、14:00、20:00(北京時,下同)4次定時觀測值,本文稱之標(biāo)準(zhǔn)時次。然而,1880—1937年蕪湖站氣溫觀測時次并非如此(表1),1880—1885年、1904—1937年每日8次觀測,觀測時次依次為03:00、06:00、09:00、12:00、15:00、18:00、21:00、24:00,平均氣溫為8個時次的均值;1886—1903年每日4次觀測,觀測時次依次為03:00、09:00、15:00、21:00,平均氣溫為4個時次的均值??梢?,觀測時次、計算方法不同,必然導(dǎo)致誤差產(chǎn)生。因此,需利用差值訂正法將1880—1937年兩種不同觀測時次的氣溫月值統(tǒng)一訂正到標(biāo)準(zhǔn)時次。

    表1 蕪湖站不同時期氣溫觀測時次及日平均氣溫的計算時次Tab.1 Observation time of temperature and calculation time of daily mean temperature in different periods at Wuhu station

    為了使訂正數(shù)據(jù)更加準(zhǔn)確,考慮到蕪湖站自2007年起開始逐小時氣溫觀測,故利用2007—2020年數(shù)據(jù)計算各月不同觀測時次與標(biāo)準(zhǔn)時次平均氣溫差值進行訂正。具體方法是:對2007—2020年各月8個觀測 時 次(03:00、06:00、09:00、12:00、15:00、18:00、21:00、24:00)和4個標(biāo)準(zhǔn)時次的氣溫求平均后相減,用該差值訂正1880—1885年及1904—1937年的月平均氣溫,即非標(biāo)準(zhǔn)時次平均氣溫減去對應(yīng)差值,從而實現(xiàn)氣溫訂正。同樣,計算2007—2020年各月4個觀測時次(03:00、09:00、15:00、21:00)與4個標(biāo)準(zhǔn)時次的月平均氣溫差值,用該差值訂正1886—1903年的月平均氣溫。

    圖1是2007—2020年蕪湖站不同觀測時次與標(biāo)準(zhǔn)時次氣溫平均誤差的月際變化??梢钥闯?,除9月外,8次和4次觀測的月平均氣溫均較標(biāo)準(zhǔn)時次偏高,誤差為0.0~0.2℃,但4次觀測的誤差更大,說明觀測時次越密越有利于誤差減小。因此,1880—1885年和1904—1937年月平均氣溫數(shù)據(jù)質(zhì)量優(yōu)于1886—1903年。利用各月平均誤差進行相應(yīng)的氣溫訂正,可進一步消除因觀測時次不同而導(dǎo)致的誤差。

    圖1 2007—2020年蕪湖站不同觀測時次與標(biāo)準(zhǔn)時次氣溫平均誤差的月際變化Fig.1 Monthly variation of mean error of temperature at different observation time and standard time at Wuhu station during 2007-2020

    伴隨氣候變暖的加快,全球氣溫日較差呈現(xiàn)不同程度的下降[20]。那么,僅利用全球變暖背景下2007—2020年數(shù)據(jù)對1880—1937年月平均氣溫進行差值訂正是否可靠?為驗證這一問題,對2007—2020年全年及不同季節(jié)8次和4次平均氣溫誤差的線性傾向率進行統(tǒng)計(表2),發(fā)現(xiàn)8次和4次平均氣溫的誤差在2007—2020年變化趨勢非常緩慢,線性傾向率均未通過α=0.05的顯著性檢驗,其中春、夏季呈現(xiàn)緩慢增大趨勢,秋、冬季則呈現(xiàn)緩慢減小趨勢,而年平均氣溫誤差的變化基本可以忽略不計。雖然2007—2020年蕪湖站總體處于增溫階段,且日較差有所減小,但其變化趨勢對8次和4次平均氣溫誤差影響較小,故推測氣溫的年代際變化對不同觀測時次平均氣溫誤差影響有限,利用2007—2020年小時氣溫數(shù)據(jù)進行不同觀測時次的差值訂正是有效的。

    表2 2007—2020年全年、不同季節(jié)8次和4次觀測氣溫的平均誤差線性變化趨勢Tab.2 Linear variation tendency of average error of temperature at eight time and four time in the whole year and different seasons during 2007-2020單位:℃·(10 a)-1

    2.1.2 均一性訂正

    采用標(biāo)準(zhǔn)正態(tài)均一性檢驗(standard normal homogeneity test,SNHT)方法[21-22],對1952—2020年蕪湖站月平均氣溫資料進行均一性檢驗。在均一性檢驗過程中,建立數(shù)據(jù)質(zhì)量較高的參考序列非常必要。因此,選取蕪湖站周邊200 km內(nèi)溧陽、長豐、黟縣3站作為參考站(圖2),這3站自建站以來未遷過站、探測環(huán)境較穩(wěn)定,序列資料完整性好且質(zhì)量高。長豐站建站較晚,資料始于1967年,故選取1967—2020年3站月平均氣溫作為參考。

    圖2 蕪湖站及周邊相鄰的CRU格點與參考站空間分布Fig.2 The spatial distribution of Wuhu station and its adjacent CRU grid points,reference stations

    經(jīng)計算,3個參考站月平均氣溫與蕪湖站的相關(guān)系數(shù)高達0.99以上。SNHT均一性檢驗(0.01的顯著性水平)發(fā)現(xiàn),建國后1988年11月、1990年7月、1992年10月、1994年7月的平均氣溫出現(xiàn)非均一性突變點(表3)。經(jīng)核查蕪湖臺站沿革資料,發(fā)現(xiàn)1988年突變是因觀測場海拔整體升高所致,1992—1994年突變則是因觀測場北部建設(shè)辦公樓所致,而1990年突變原因不明。另外,自建站以來,蕪湖站分別在1955年4月1日(118°23′E,31°20′N;海拔15.2 m)、2006年1月1日(118°22′E,31°23′N;9.5 m)、2016年1月1日(118°21′E,31°20′N;海拔12.8 m)經(jīng)歷3次遷站(表略),遷站前后地理條件相似,并未造成氣溫的非均一性突變。

    1901—1937年蕪湖站周邊無觀測資料連續(xù)且質(zhì)量較好的參考站,故選取該站周邊相鄰且相關(guān)性最高的6個CRU格點(圖2)資料作為參考氣溫。統(tǒng)計顯示(表4),與距離蕪湖站最近的CRU格點氣溫相比,6個格點月平均氣溫與蕪湖站的均方根誤差明顯偏小,偏差也相對較小,相關(guān)系數(shù)更高,更接近觀測值。因此,用6個CRU格點的月平均氣溫作為參考氣溫更為合理。同樣,采用SNHT方法對1901—1937年蕪湖站實測月平均氣溫進行均一性檢驗,發(fā)現(xiàn)有6個數(shù)據(jù)出現(xiàn)突變(表3),且數(shù)值偏大。由于建國前無蕪湖臺站沿革信息,故無法推斷出現(xiàn)非均一性的原因。1880—1900年沒有CRU數(shù)據(jù)作為參考,故該時段內(nèi)蕪湖站月平均氣溫采用距平累積法和極值檢查,未發(fā)現(xiàn)明顯的突變點。

    表3 蕪湖站建國前后月平均氣溫的非均一性突變點及建國后突變原因Tab.3 Non-uniformity mutation points of monthly mean temperature at Wuhu station before and after the founding of new China and the reasons after the founding of new China

    表4 1901—1937年蕪湖站月平均氣溫與周邊CRU格點氣溫的統(tǒng)計關(guān)系Tab.4 Statistical relations between monthly observational temperature at Wuhu station and its adjacent grid point temperature of CRU during 1901-1937

    2.1.3 異常值插補

    由于1952—2020年、1901—1937年蕪湖站月平均氣溫與參考氣溫存在極顯著正相關(guān)關(guān)系,故分別建立這2個時段兩者的線性回歸方程,并利用回歸方程擬合數(shù)據(jù)對非均一性突變點進行插補。

    鑒于1952—2020年突變點較少,插補后的序列和原始序列差異不大(圖略)。1880—1937年蕪湖站經(jīng)差值、均一化訂正及插補后的氣溫序列略低于原始序列,特別在1891—1911年訂正效果明顯(圖3)??梢?,通過上述訂正及插補,可以消除觀測時次不同及非均一性導(dǎo)致的誤差,明顯提高了1880—1937年數(shù)據(jù)的可用性。

    圖3 1880—1937年蕪湖站原始和訂正插補的平均氣溫時間序列Fig.3 Time series of original and adjusted and interpolated mean temperature at Wuhu station during 1880-1937

    2.2 1938—1951年氣溫序列構(gòu)建

    1938—1951年蕪湖站觀測中斷,需利用CRU格點氣溫數(shù)據(jù)進行構(gòu)建。由于CRU格點氣溫與蕪湖站觀測氣溫存在一定偏差,直接代替勢必造成數(shù)據(jù)的非均一性。經(jīng)統(tǒng)計,蕪湖站鄰近的CRU格點氣溫比觀測值偏低,1901—1937年、1953—2020年(考慮到1952年1月氣溫缺測,故從1953年開始計算)分別偏低0.9、0.4℃。故基于這2個時段蕪湖站月觀測氣溫與上述周邊相鄰的6個CRU格點氣溫,采用多元逐步回歸分析方法,分別建立蕪湖站月平均氣溫的最優(yōu)擬合方程,對1938—1951年氣溫進行擬合,取平均值進行插補。多元逐步回歸計算公式[23]如下:

    式中:(℃)為站點擬合氣溫;yj(℃)為站點周邊相鄰的第j個格點氣溫;β0是站點氣溫擬合常數(shù);βj是站點周邊第j個格點的氣溫擬合系數(shù);ε(℃)是站點氣溫殘差;n為格點數(shù),本文取6。

    以構(gòu)建1901—1937年1月最優(yōu)擬合方程為例,首先計算該時段內(nèi)蕪湖站觀測氣溫與周邊多個CRU格點氣溫的相關(guān)系數(shù),選出6個相關(guān)系數(shù)最高的格點氣溫作為自變量,蕪湖站觀測氣溫作為因變量,建立多元逐步回歸方程;然后,基于1901—2020年上述CRU格點氣溫,利用多元逐步回歸擬合方程,擬合同期蕪湖站氣溫。同樣,基于1953—2020年蕪湖站觀測氣溫和CRU格點氣溫建立擬合方程,并擬合1901—2020年蕪湖站氣溫。以此類推,建立兩套蕪湖站月擬合氣溫序列。

    從圖4看出,蕪湖站觀測氣溫與基于1901—1937年、1953—2020年數(shù)據(jù)的擬合氣溫存在一定偏差,總體上前者較觀測平均偏高0.25℃,后者平均偏低0.19℃,但兩套擬合氣溫與觀測氣溫的變化趨勢基本一致,且2個偏差序列的標(biāo)準(zhǔn)差接近,說明2個時段數(shù)據(jù)擬合的氣溫偏差變化穩(wěn)定。1938年以前,兩套擬合氣溫絕大部分年份低于觀測氣溫,但基于1901—1937年數(shù)據(jù)擬合的氣溫更接近觀測;1953年以后,基于1901—1937年數(shù)據(jù)擬合的氣溫明顯高于觀測,而基于1953—2020年數(shù)據(jù)擬合的氣溫更接近觀測,其中20世紀(jì)90年代以前擬合氣溫略偏高,之后擬合氣溫略偏低。可見,1938—1951年月平均氣溫用兩套擬合氣溫序列的平均值進行插補,比CRU格點序列或單一的擬合氣溫序列更加準(zhǔn)確。

    圖4 蕪湖站1901—2020年訂正插補的觀測氣溫及擬合氣溫偏差時間序列Fig.4 Time series of adjusted and interpolated observation temperature and deviation of fitting temperature at Wuhu station during 1901-2020

    3 蕪湖平均氣溫的年代際特征

    3.1 變化特征

    圖5是1880—2020年蕪湖站全年及各季平均氣溫的變化曲線??梢钥闯?,近140 a蕪湖站年平均氣溫整體呈顯著上升趨勢(通過α=0.01的顯著性檢驗),氣候傾向率為0.06℃·(10 a)-1,低于中國百年氣溫增速[5-9],且重建的1938—1951年序列沒有明顯突變[圖5(a)]。與周邊上海站對比發(fā)現(xiàn),蕪湖、上海年平均氣溫在20世紀(jì)40年代均略有上升[16],這間接表明了重構(gòu)的1938—1951年數(shù)據(jù)能夠較好地反映缺測時段蕪湖站氣溫的年際變化。另外,蕪湖氣溫變化存在明顯的階段性,19世紀(jì)末為顯著偏冷時段,累積距平始終為負值;20世紀(jì)初至20年代為偏暖時段,隨后30年代偏冷,40年代偏暖,50—80年代又轉(zhuǎn)為明顯偏冷時段,90年代之后氣溫快速回升,累積距平持續(xù)為正值,增溫持續(xù)時間最長。其中,21世紀(jì)以來平均氣溫較19世紀(jì)末明顯偏高1.3℃。

    與年平均氣溫變化趨勢一致,近140 a蕪湖站四季氣溫也呈現(xiàn)增加趨勢,春季增溫最顯著,氣候傾向率為0.12℃·(10 a)-1,遠超出其他季節(jié),對年平均氣溫增溫的貢獻最大。春季氣溫在19世紀(jì)末至20世紀(jì)20年代為持續(xù)穩(wěn)定增溫階段,之后與年平均氣溫的變化基本一致[圖5(b)]。研究表明,1885年以來長江中下游地區(qū)四季氣溫百年增加0.59~0.92℃,春季增溫最大[24]。在同時段內(nèi),蕪湖春季增溫趨勢略高于長江中下游區(qū)域平均。夏、冬季氣溫的氣候傾向率與年平均氣溫接近,分別為0.06、0.07℃·(10 a)-1(通過α=0.01的顯著性檢驗)。夏季,除19世紀(jì)末氣溫異常偏低外,20世紀(jì)初至60年代氣溫的年代際差異較小,冷暖交替變化沒有春季顯著,而70—80年代為偏冷期,90年代之后轉(zhuǎn)為偏暖期,長期維持正距平[圖5(c)],這一點與年平均氣溫變化一致;冬季氣溫也呈現(xiàn)冷暖交替的年代際變化特征,19世紀(jì)末異常偏冷,20世紀(jì)初至20年代偏暖,30—80年代持續(xù)偏冷,90年代以后增暖趨勢顯著[圖5(e)]。秋季增溫最弱,氣候傾向率僅為0.02℃·(10 a)-1[圖5(d)],未通過α=0.01的顯著性檢驗,與年平均氣溫變化類似,在20世紀(jì)初至20年代、90年代之后存在2個暖期。綜上所述,近140 a蕪湖站年及季節(jié)(秋季除外)氣溫整體顯著上升,且存在冷暖交替的年代際變化特征,但不同季節(jié)冷暖變化時期略有不同。

    圖5 1880—2020年蕪湖站年(a)和春季(b)、夏季(c)、秋季(d)、冬季(e)平均氣溫時間序列及年代際變化Fig.5 The time series and inter-decadal variation of mean temperature in the whole year(a),spring(b),summer(c),autumn(d)and winter(e)at Wuhu station from 1880 to 2020

    鑒于四季氣溫的氣候傾向率與所選時段有很大關(guān)系,故計算四季氣溫在不同時段的氣候傾向率(圖6)??梢钥闯?,春季氣溫在各時段的增溫率均最大,遠高于其他季節(jié);冬季增溫率次之,與前面分析一致。雖然1880—2020年蕪湖增溫趨勢弱于全國平均[5-6],但隨著時間的推移和統(tǒng)計時段的縮短,各季節(jié)增溫趨勢逐漸增大,1980—2020年四季增溫趨勢最強,增溫率在0.33℃·(10 a)-1以上,其中春季增溫率高達0.65℃·(10 a)-1。然而,進入21世紀(jì)以來,四季增溫率明顯減弱,夏、秋季氣溫甚至出現(xiàn)下降趨勢,這一特征與東亞增溫停滯[25-26]相吻合,說明近20 a來蕪湖增溫減緩。KNIGHT等[25]研究指出,1998年全球氣溫達到最高點后,增溫率不再繼續(xù)上升,這主要是CRU資料的問題。事實上,蕪湖站氣溫觀測資料顯示21世紀(jì)以來增溫停滯現(xiàn)象客觀存在。盡管增溫有所放緩,但氣候變化是一個長期過程,需要從幾十年甚至上百年的尺度進行研究。

    圖6 不同時段蕪湖站季節(jié)平均氣溫的氣候傾向率Fig.6 Climatic tendency rates of seasonal mean temperature in different periods at Wuhu station

    3.2 周期特征

    對近140 a蕪湖站年和四季平均氣溫進行小波分析(圖7)??梢钥闯觯徍昙八募酒骄鶜鉁卮嬖?4 a以上時間尺度的周期振蕩,但該周期多處于邊緣效應(yīng)影響區(qū)域,故真實性有待商榷。另外,年及秋季(圖略)、冬季平均氣溫均存在40~50 a的周期振蕩,對應(yīng)該時段冷暖交替頻繁。此外,年、四季平均氣溫在不同時段還存在20~30 a的周期振蕩,其中春、冬季該周期貫穿整個時域,而年及夏、秋季該周期僅集中出現(xiàn)在19世紀(jì)末至20世紀(jì)60年代。

    圖7 1880—2020年蕪湖年(a)及春(b)、夏(c)、冬(d)季平均氣溫的小波功率譜Fig.7 The wavelet power spectrum of annual(a),spring(b),summer(c)and winter(d)mean temperature during 1880-2020 at Wuhu station

    4 結(jié)論

    本文利用蕪湖站百年氣溫觀測數(shù)據(jù)和周邊CRU格點再分析資料,在差值和均一化訂正基礎(chǔ)上,采用多元逐步回歸分析方法,構(gòu)建1880—2020年蕪湖月平均氣溫序列,并對年及四季平均氣溫的年代際變化特征進行分析??傮w上,蕪湖站氣溫原始序列個別存在非均一性突變,且1880—1937年月平均氣溫偏高,經(jīng)差值和均一化訂正后數(shù)據(jù)質(zhì)量得到明顯改善,年及季節(jié)平均氣溫均呈現(xiàn)冷暖交替的年代際變化,不同季節(jié)變化特征略有差異。具體結(jié)論如下:

    (1)不同觀測時次計算的月平均氣溫有所差異。與目前標(biāo)準(zhǔn)時次(02:00、08:00、14:00、20:00)相比,建國前8次(03:00、06:00、09:00、12:00、15:00、18:00、21:00、24:00)和4次(03:00、09:00、15:00、21:00)觀測的平均氣溫均存在誤差,且4次觀測的氣溫誤差更大。均一性檢驗發(fā)現(xiàn),局部環(huán)境改變會引起氣溫突變。經(jīng)差值和均一化訂正后,蕪湖1880—1937年月平均氣溫數(shù)據(jù)的可用性明顯提高。

    (2)基于蕪湖站1901—1937年、1953—2020年月觀測氣溫和1901—2020年CRU格點氣溫,得到的2個擬合氣溫序列的均值比CRU格點氣溫序列或單一擬合氣溫序列更接近觀測值,故利用2個擬合氣溫的均值對1938—1951年數(shù)據(jù)進行插補。

    (3)1880—2020年蕪湖年、季節(jié)平均氣溫均呈現(xiàn)冷暖交替的年代際變化,其中20世紀(jì)80年代末由冷到暖的交替最顯著。近140 a蕪湖四季一致升溫,春季升溫最明顯,夏季、冬季次之,秋季增溫趨勢不顯著。2000年之前蕪湖增溫速率不斷加快,但之后四季增溫有所減緩,特別是夏、秋季甚至出現(xiàn)弱的降溫趨勢。

    (4)近140 a蕪湖年和四季平均氣溫在不同時段存在20~30 a的周期振蕩,且年及秋、冬季平均氣溫還存在40~50 a的周期變化。

    观看美女的网站| 久久久精品免费免费高清| 亚洲色图av天堂| 亚洲精品乱码久久久久久按摩| 免费电影在线观看免费观看| 欧美一区二区亚洲| 在线免费观看不下载黄p国产| 直男gayav资源| 欧美变态另类bdsm刘玥| 1000部很黄的大片| 亚洲国产欧美在线一区| 亚洲av免费高清在线观看| 精品99又大又爽又粗少妇毛片| 国产精品久久久久久久久免| 欧美人与善性xxx| 日本熟妇午夜| 日韩欧美精品v在线| 日韩强制内射视频| 久久久久性生活片| 特级一级黄色大片| 日韩伦理黄色片| 淫秽高清视频在线观看| 亚洲av免费高清在线观看| 欧美成人午夜免费资源| 久久久久久久久久久免费av| 精品久久久久久电影网| 一级毛片我不卡| 亚洲av不卡在线观看| 国产日韩欧美在线精品| 国内精品美女久久久久久| 亚洲综合色惰| 亚洲欧美清纯卡通| 国产淫语在线视频| 91午夜精品亚洲一区二区三区| 国产一区亚洲一区在线观看| 国产精品嫩草影院av在线观看| 国产色婷婷99| 国产欧美日韩精品一区二区| 男女那种视频在线观看| 夜夜爽夜夜爽视频| 久久精品熟女亚洲av麻豆精品 | 国产综合懂色| 国产成人a区在线观看| 伦理电影大哥的女人| 3wmmmm亚洲av在线观看| videossex国产| 极品教师在线视频| 黄色一级大片看看| 99热网站在线观看| 国产淫片久久久久久久久| 久久久久久久亚洲中文字幕| 亚洲第一区二区三区不卡| 国产成人91sexporn| 国产精品久久视频播放| 99热这里只有是精品50| 日韩欧美 国产精品| 777米奇影视久久| 亚洲av免费高清在线观看| 成人美女网站在线观看视频| 边亲边吃奶的免费视频| 日本猛色少妇xxxxx猛交久久| 成人性生交大片免费视频hd| 成人午夜精彩视频在线观看| 街头女战士在线观看网站| 青春草亚洲视频在线观看| 国产黄色视频一区二区在线观看| 日韩在线高清观看一区二区三区| av免费观看日本| 欧美激情在线99| 国产91av在线免费观看| av天堂中文字幕网| 一级毛片 在线播放| 亚洲综合精品二区| av在线老鸭窝| 夜夜爽夜夜爽视频| 综合色丁香网| 亚洲国产精品成人综合色| 成年免费大片在线观看| 秋霞伦理黄片| 欧美+日韩+精品| 在线天堂最新版资源| 日韩欧美精品v在线| 九色成人免费人妻av| 97在线视频观看| 一级毛片电影观看| 中文精品一卡2卡3卡4更新| 噜噜噜噜噜久久久久久91| 嫩草影院新地址| 久久韩国三级中文字幕| 国产老妇伦熟女老妇高清| 干丝袜人妻中文字幕| 国产精品久久久久久久久免| 精品久久久久久久人妻蜜臀av| 黄色配什么色好看| 欧美日韩在线观看h| 国精品久久久久久国模美| 欧美高清性xxxxhd video| 黄色日韩在线| 少妇的逼好多水| 午夜爱爱视频在线播放| 看免费成人av毛片| 亚洲自偷自拍三级| 午夜激情福利司机影院| 国产毛片a区久久久久| 搡女人真爽免费视频火全软件| 国产高清有码在线观看视频| 天天躁日日操中文字幕| 九九在线视频观看精品| 亚洲乱码一区二区免费版| 日产精品乱码卡一卡2卡三| 丰满人妻一区二区三区视频av| 一级二级三级毛片免费看| 一区二区三区乱码不卡18| 国内精品宾馆在线| 天堂俺去俺来也www色官网 | av国产免费在线观看| 久久韩国三级中文字幕| 亚洲丝袜综合中文字幕| 91av网一区二区| 日日干狠狠操夜夜爽| 久久精品国产鲁丝片午夜精品| 精品国产露脸久久av麻豆 | 免费黄网站久久成人精品| 三级男女做爰猛烈吃奶摸视频| 免费av毛片视频| 亚洲最大成人手机在线| 午夜免费男女啪啪视频观看| 日韩精品有码人妻一区| 一本一本综合久久| 色5月婷婷丁香| 精品久久久久久久久亚洲| av国产久精品久网站免费入址| 国产成人a区在线观看| 久久久亚洲精品成人影院| 日本免费a在线| 亚洲精品乱码久久久v下载方式| 国产有黄有色有爽视频| 国产精品精品国产色婷婷| 久久久久久久大尺度免费视频| 欧美最新免费一区二区三区| 精品久久久久久久久av| 特大巨黑吊av在线直播| 国产免费一级a男人的天堂| 亚洲精品国产成人久久av| 又爽又黄无遮挡网站| 看十八女毛片水多多多| 丰满乱子伦码专区| 丝瓜视频免费看黄片| 永久网站在线| 天美传媒精品一区二区| 成人综合一区亚洲| 日韩伦理黄色片| 免费观看的影片在线观看| 久久久精品免费免费高清| 91av网一区二区| 精品久久久精品久久久| 91精品一卡2卡3卡4卡| 免费播放大片免费观看视频在线观看| 最近2019中文字幕mv第一页| 大香蕉久久网| 狂野欧美激情性xxxx在线观看| 特大巨黑吊av在线直播| 22中文网久久字幕| 亚洲av一区综合| 插逼视频在线观看| 美女黄网站色视频| 女人被狂操c到高潮| 国产精品一区二区在线观看99 | 久久精品国产鲁丝片午夜精品| av在线蜜桃| 美女xxoo啪啪120秒动态图| 99久久人妻综合| 日本黄大片高清| 男女国产视频网站| 午夜久久久久精精品| 午夜福利高清视频| 亚洲精品影视一区二区三区av| 观看免费一级毛片| 精品人妻视频免费看| videossex国产| 国产在线男女| 免费电影在线观看免费观看| 久久草成人影院| 亚洲国产高清在线一区二区三| or卡值多少钱| 亚洲av成人精品一区久久| 蜜臀久久99精品久久宅男| 午夜福利视频精品| 亚洲无线观看免费| 九色成人免费人妻av| 如何舔出高潮| 日韩精品青青久久久久久| 一个人观看的视频www高清免费观看| 国产色婷婷99| 婷婷色综合大香蕉| 99热6这里只有精品| 欧美不卡视频在线免费观看| 成人特级av手机在线观看| 免费黄网站久久成人精品| 久久久久久九九精品二区国产| 一夜夜www| 在线观看人妻少妇| eeuss影院久久| 国产精品嫩草影院av在线观看| 欧美激情国产日韩精品一区| 国产成人精品婷婷| 日本三级黄在线观看| 亚洲精品乱码久久久久久按摩| 六月丁香七月| 国产精品99久久久久久久久| 亚洲av福利一区| 久热久热在线精品观看| 国产免费一级a男人的天堂| 啦啦啦中文免费视频观看日本| 国产色婷婷99| 亚洲国产av新网站| 男人狂女人下面高潮的视频| 熟女电影av网| 青青草视频在线视频观看| 日日啪夜夜爽| 少妇被粗大猛烈的视频| 韩国高清视频一区二区三区| 国产成人精品久久久久久| 综合色丁香网| 国产成人午夜福利电影在线观看| 五月玫瑰六月丁香| 高清在线视频一区二区三区| 日本色播在线视频| 麻豆乱淫一区二区| 成人毛片a级毛片在线播放| 免费观看无遮挡的男女| 18+在线观看网站| av免费在线看不卡| 精品99又大又爽又粗少妇毛片| 日韩一本色道免费dvd| 六月丁香七月| 日日摸夜夜添夜夜添av毛片| 国产成人福利小说| 99热这里只有是精品50| 黄片wwwwww| 午夜爱爱视频在线播放| 欧美三级亚洲精品| 六月丁香七月| 欧美激情久久久久久爽电影| 中国国产av一级| 在线观看人妻少妇| 搡女人真爽免费视频火全软件| 秋霞在线观看毛片| av.在线天堂| 午夜亚洲福利在线播放| 亚洲精华国产精华液的使用体验| 亚洲精品乱久久久久久| 亚洲图色成人| 精华霜和精华液先用哪个| 精品人妻熟女av久视频| 高清午夜精品一区二区三区| 欧美激情在线99| 国产成人午夜福利电影在线观看| 日本猛色少妇xxxxx猛交久久| 别揉我奶头 嗯啊视频| 国产精品福利在线免费观看| 午夜爱爱视频在线播放| av在线蜜桃| 美女内射精品一级片tv| www.av在线官网国产| 久久精品综合一区二区三区| 亚洲在线观看片| 亚洲欧美成人综合另类久久久| 免费观看的影片在线观看| 一级毛片aaaaaa免费看小| 夜夜看夜夜爽夜夜摸| 97在线视频观看| 亚洲国产成人一精品久久久| 美女被艹到高潮喷水动态| 欧美区成人在线视频| 日韩三级伦理在线观看| 99热全是精品| 男女啪啪激烈高潮av片| 久久久久网色| 国内精品美女久久久久久| 欧美成人a在线观看| 日日摸夜夜添夜夜添av毛片| 国产亚洲午夜精品一区二区久久 | 精品不卡国产一区二区三区| 久久精品久久精品一区二区三区| 中国美白少妇内射xxxbb| 久久久国产一区二区| 别揉我奶头 嗯啊视频| 国产一区二区三区av在线| 最近2019中文字幕mv第一页| 午夜福利高清视频| 成人性生交大片免费视频hd| 国产精品久久久久久久电影| 91精品伊人久久大香线蕉| 国产伦一二天堂av在线观看| a级毛片免费高清观看在线播放| 久久久久久久亚洲中文字幕| 丰满人妻一区二区三区视频av| 中文字幕免费在线视频6| 2021天堂中文幕一二区在线观| 国产精品美女特级片免费视频播放器| 日韩欧美 国产精品| 插阴视频在线观看视频| 人体艺术视频欧美日本| 亚洲国产最新在线播放| 亚洲经典国产精华液单| 亚洲高清免费不卡视频| 免费av不卡在线播放| 久久精品久久精品一区二区三区| 不卡视频在线观看欧美| 一夜夜www| 中文资源天堂在线| 深爱激情五月婷婷| 欧美+日韩+精品| 激情五月婷婷亚洲| 免费人成在线观看视频色| 日韩欧美三级三区| 国产精品麻豆人妻色哟哟久久 | 欧美潮喷喷水| 中文精品一卡2卡3卡4更新| 老师上课跳d突然被开到最大视频| 久久精品久久久久久噜噜老黄| 男人狂女人下面高潮的视频| av免费观看日本| 国产69精品久久久久777片| 亚洲精品日韩在线中文字幕| 亚洲欧美中文字幕日韩二区| 1000部很黄的大片| 国产精品一区www在线观看| 高清午夜精品一区二区三区| 春色校园在线视频观看| av.在线天堂| 丰满少妇做爰视频| 成人亚洲欧美一区二区av| 欧美激情国产日韩精品一区| 欧美性感艳星| 久久久久久久久久成人| 最近中文字幕高清免费大全6| 亚洲精品一二三| a级毛片免费高清观看在线播放| 成人一区二区视频在线观看| 国产精品一区二区三区四区久久| 成人综合一区亚洲| 亚洲欧洲日产国产| 久久久欧美国产精品| 男女边吃奶边做爰视频| 精品久久久久久久人妻蜜臀av| 床上黄色一级片| 欧美一级a爱片免费观看看| 免费不卡的大黄色大毛片视频在线观看 | 成人毛片60女人毛片免费| 久久国内精品自在自线图片| 麻豆久久精品国产亚洲av| 久久国内精品自在自线图片| 三级国产精品欧美在线观看| 精品一区二区三区视频在线| 国产永久视频网站| 51国产日韩欧美| 大话2 男鬼变身卡| 能在线免费观看的黄片| 亚洲精品自拍成人| 国产精品一二三区在线看| 午夜日本视频在线| 麻豆av噜噜一区二区三区| 国产亚洲5aaaaa淫片| 人妻系列 视频| 亚洲真实伦在线观看| 免费在线观看成人毛片| 2021天堂中文幕一二区在线观| 麻豆精品久久久久久蜜桃| 国产av在哪里看| 99热6这里只有精品| 婷婷色麻豆天堂久久| 久久久久性生活片| 免费看a级黄色片| 青春草国产在线视频| 国产黄频视频在线观看| 久热久热在线精品观看| 搡女人真爽免费视频火全软件| 亚洲不卡免费看| 亚洲国产av新网站| 国产精品麻豆人妻色哟哟久久 | 伦精品一区二区三区| 在线观看美女被高潮喷水网站| 中国美白少妇内射xxxbb| 一级毛片久久久久久久久女| www.av在线官网国产| 一二三四中文在线观看免费高清| 国产男女超爽视频在线观看| 春色校园在线视频观看| 成人毛片a级毛片在线播放| 国产精品美女特级片免费视频播放器| 蜜臀久久99精品久久宅男| 韩国高清视频一区二区三区| 少妇被粗大猛烈的视频| 又爽又黄a免费视频| 最近最新中文字幕免费大全7| 在线观看人妻少妇| 国产黄a三级三级三级人| 成人欧美大片| 中国国产av一级| 亚洲,欧美,日韩| 69人妻影院| 免费大片黄手机在线观看| 久久韩国三级中文字幕| 你懂的网址亚洲精品在线观看| 国内精品宾馆在线| 亚洲精品自拍成人| 嘟嘟电影网在线观看| 男女下面进入的视频免费午夜| 99热这里只有精品一区| 色吧在线观看| 亚洲丝袜综合中文字幕| 深爱激情五月婷婷| 网址你懂的国产日韩在线| 午夜老司机福利剧场| 亚洲欧美精品自产自拍| 国产成人aa在线观看| 国产亚洲5aaaaa淫片| 国产一级毛片在线| 国产精品久久久久久精品电影| 乱码一卡2卡4卡精品| 一级片'在线观看视频| 99久久精品一区二区三区| 中国国产av一级| 精品少妇黑人巨大在线播放| 好男人视频免费观看在线| 免费观看在线日韩| 午夜激情欧美在线| 亚洲真实伦在线观看| 极品教师在线视频| 国产成人免费观看mmmm| 秋霞伦理黄片| 欧美日韩视频高清一区二区三区二| 九草在线视频观看| 最近手机中文字幕大全| 日韩一本色道免费dvd| 色视频www国产| 3wmmmm亚洲av在线观看| 丝袜喷水一区| 色综合站精品国产| 成人高潮视频无遮挡免费网站| 一二三四中文在线观看免费高清| 免费观看在线日韩| 99热全是精品| 最新中文字幕久久久久| 国产高清三级在线| 夫妻午夜视频| 亚洲最大成人中文| 婷婷色综合www| 22中文网久久字幕| 久热久热在线精品观看| av免费在线看不卡| 赤兔流量卡办理| 校园人妻丝袜中文字幕| 插逼视频在线观看| 午夜激情福利司机影院| 少妇丰满av| 亚洲天堂国产精品一区在线| 国产亚洲午夜精品一区二区久久 | 国产三级在线视频| 熟妇人妻不卡中文字幕| av又黄又爽大尺度在线免费看| 午夜视频国产福利| 免费观看无遮挡的男女| 永久网站在线| 色视频www国产| 搡女人真爽免费视频火全软件| 亚洲在线观看片| 国产高清国产精品国产三级 | 欧美日本视频| 成年人午夜在线观看视频 | av国产免费在线观看| 成年免费大片在线观看| 午夜日本视频在线| 国产黄色视频一区二区在线观看| 人妻夜夜爽99麻豆av| 久久国内精品自在自线图片| 国产男女超爽视频在线观看| 18禁在线无遮挡免费观看视频| 久久久久国产网址| 搡老乐熟女国产| a级一级毛片免费在线观看| 亚洲伊人久久精品综合| 中文在线观看免费www的网站| 国产三级在线视频| 99热6这里只有精品| 人妻夜夜爽99麻豆av| 日韩 亚洲 欧美在线| 男女啪啪激烈高潮av片| 青春草国产在线视频| 日本色播在线视频| 国产单亲对白刺激| 久久久久精品久久久久真实原创| 亚洲精品中文字幕在线视频 | 国产精品人妻久久久影院| 丝瓜视频免费看黄片| 国产成人91sexporn| 九九爱精品视频在线观看| 91aial.com中文字幕在线观看| 一级av片app| 亚洲熟女精品中文字幕| 特大巨黑吊av在线直播| 男女下面进入的视频免费午夜| 91精品国产九色| 亚洲欧美成人综合另类久久久| 亚洲图色成人| 美女国产视频在线观看| 别揉我奶头 嗯啊视频| 国产伦精品一区二区三区视频9| 免费看av在线观看网站| 婷婷色麻豆天堂久久| 成人鲁丝片一二三区免费| 日产精品乱码卡一卡2卡三| 水蜜桃什么品种好| 永久免费av网站大全| 久久久久网色| 亚洲精品亚洲一区二区| 最近2019中文字幕mv第一页| 亚洲av中文av极速乱| 国产中年淑女户外野战色| 汤姆久久久久久久影院中文字幕 | 亚洲,欧美,日韩| 成人二区视频| 中文欧美无线码| 韩国av在线不卡| 日产精品乱码卡一卡2卡三| 狂野欧美白嫩少妇大欣赏| 只有这里有精品99| 亚洲av电影不卡..在线观看| 日韩不卡一区二区三区视频在线| 国产亚洲av片在线观看秒播厂 | 一级a做视频免费观看| 久久这里有精品视频免费| 欧美性猛交╳xxx乱大交人| 免费看光身美女| 国精品久久久久久国模美| 中国美白少妇内射xxxbb| 噜噜噜噜噜久久久久久91| 国产精品99久久久久久久久| 校园人妻丝袜中文字幕| 三级经典国产精品| 亚洲欧美清纯卡通| 永久免费av网站大全| 国产精品嫩草影院av在线观看| 精品国产露脸久久av麻豆 | 亚洲久久久久久中文字幕| 久久久久久伊人网av| 国产精品麻豆人妻色哟哟久久 | 欧美日韩亚洲高清精品| 美女大奶头视频| 精品久久久久久成人av| 日韩制服骚丝袜av| 午夜久久久久精精品| 亚洲精品色激情综合| 在线观看人妻少妇| 亚洲18禁久久av| 成人午夜精彩视频在线观看| 91久久精品国产一区二区成人| 国产精品综合久久久久久久免费| 丝瓜视频免费看黄片| 狂野欧美激情性xxxx在线观看| 欧美区成人在线视频| 国产69精品久久久久777片| 国产一区亚洲一区在线观看| 久久久久国产网址| 国产麻豆成人av免费视频| 色综合亚洲欧美另类图片| 中国美白少妇内射xxxbb| 国产亚洲精品久久久com| 亚洲精品成人久久久久久| av专区在线播放| 精品久久久久久久末码| 国产黄频视频在线观看| 黄色欧美视频在线观看| 91久久精品国产一区二区成人| 亚洲激情五月婷婷啪啪| 中文字幕人妻熟人妻熟丝袜美| 综合色av麻豆| 免费看光身美女| 最后的刺客免费高清国语| 真实男女啪啪啪动态图| 国产精品久久久久久久电影| 男人爽女人下面视频在线观看| 我的老师免费观看完整版| 秋霞在线观看毛片| 少妇被粗大猛烈的视频| 99热这里只有是精品50| 麻豆国产97在线/欧美| 久久久精品免费免费高清| 日本av手机在线免费观看| 色网站视频免费| 国产淫片久久久久久久久| 久久精品人妻少妇| 三级毛片av免费| 精品人妻视频免费看| 亚洲天堂国产精品一区在线| 国产黄a三级三级三级人| 不卡视频在线观看欧美| 99久久精品国产国产毛片| 亚洲精品自拍成人| 蜜桃久久精品国产亚洲av| 97超视频在线观看视频| 天堂网av新在线| a级毛片免费高清观看在线播放| 国内精品一区二区在线观看| 性色avwww在线观看| 日韩欧美 国产精品| 天堂俺去俺来也www色官网 | 激情五月婷婷亚洲| 夜夜看夜夜爽夜夜摸| 精品99又大又爽又粗少妇毛片| 青春草亚洲视频在线观看| 两个人的视频大全免费| 又爽又黄a免费视频|