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

    氣候變化對(duì)長(zhǎng)江中下游地區(qū)水稻產(chǎn)量的影響

    2014-03-22 10:49王保黃思先孫衛(wèi)國(guó)
    湖北農(nóng)業(yè)科學(xué) 2014年1期
    關(guān)鍵詞:水稻產(chǎn)量長(zhǎng)江中下游地區(qū)小波變換

    王保+黃思先+孫衛(wèi)國(guó)

    摘要:為了揭示區(qū)域氣候變化對(duì)長(zhǎng)江中下游地區(qū)水稻產(chǎn)量的影響及水稻相對(duì)氣象產(chǎn)量波動(dòng)的原因,采用統(tǒng)計(jì)分析方法和小波變換方法分析了近60年來(lái)長(zhǎng)江中下游地區(qū)水稻產(chǎn)量、水稻生長(zhǎng)季內(nèi)平均氣溫、降水量、氣溫日較差以及≥10 ℃活動(dòng)積溫的年際變化以及水稻產(chǎn)量與區(qū)域氣候變化之間的時(shí)頻結(jié)構(gòu)特征。結(jié)果表明,近60年來(lái)長(zhǎng)江中下游地區(qū)水稻產(chǎn)量呈波動(dòng)增長(zhǎng)趨勢(shì);20世紀(jì)80年代以后水稻生長(zhǎng)季內(nèi)平均氣溫升高、≥10 ℃活動(dòng)積溫增大,降水量減少趨勢(shì)不明顯,氣溫日較差呈顯著下降趨勢(shì)。水稻產(chǎn)量和區(qū)域氣候變化均存在年際和年代際尺度變化周期,且兩者之間的時(shí)頻特征存在一定的相似性;交叉小波分析表明,水稻產(chǎn)量與區(qū)域氣候變化之間存在2~4a、6a、8a、12a和14a尺度的共振頻率,且多以正相關(guān)關(guān)系為主,但在局部時(shí)域中也存在不同頻率尺度的負(fù)相關(guān),兩者之間的相關(guān)關(guān)系和相關(guān)程度隨振蕩頻率尺度的不同而不同。長(zhǎng)江中下游地區(qū)水稻產(chǎn)量波動(dòng)與區(qū)域氣候變化密切相關(guān),氣候變暖導(dǎo)致熱量條件改善對(duì)水稻產(chǎn)量的影響大于降水量減少對(duì)水稻產(chǎn)量的影響。

    關(guān)鍵詞:長(zhǎng)江中下游地區(qū);水稻產(chǎn)量;氣候變化;小波變換

    中圖分類號(hào):S162.5+3 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):0439-8114(2014)01-0043-09

    Effects of Climate Change on Rice Yield of the Middle and Lower Reaches Region of the Yangtze River

    WANG Bao1,HUANG Si-xian2,SUN Wei-guo3

    (1.Xianning Bureau of Meterology, Xianning 437100,Hubei,China;2. Ezhou Bureau of Meterology,Ezhou 436000,Hubei,China;

    3.Nanjing University of Information Science & Technology, Naning 210044,China)

    Abstract: To reveal the impact of the regional climate change on rice yield of the middle and lower reaches of the Yangtze River, and to understand the reasons for rice relative meteorological yield fluctuations, statistical analysis and wavelet transform were used to analyze time and the frequency variation features of rice production, the average temperature, precipitation, diurnal temperature, and≥10℃ active accumulated temperature of the Yangtze River region in the past 60 years, and the time-frequency structural features between rice production and regional climate change. The results showed that rice yield of the middle and lower reaches of the Yangtze River region in the past 60 years had fluctuations in the growth trend. The average temperature and the active accumulated temperature increased during the growing season. The reduction of precipitation was not obvious. Diurnal temperature appeared to decrease significantly. Between rice production and climate change, there were periodic characteristics of interannual and decadal changes. The time-frequency characteristics was certainly similar between the two. Cross-wavelet transform results showed that between rice production and climate change, there were the resonance frequency of the different scales, like 2-4a, 6a, 8a, 12a and 14a. There were more positive correlation between the main, but in the local time domain there was a negative correlation of different frequency scale. The two relationship and relevance were different, varying from oscillation periods scale. It was believed that rice yield fluctuations of middle and lower reaches of the Yangtze River was closely related with the climate change. The periodic feature was obvious, and the effect of the improvement of heat condition on increasing production was much more than the impact of reducing precipitation.

    Key words: the middle and lower reaches of the Yangtze River; rice yield; climate change; wavelet transform

    收稿日期:2013-07-16

    作者簡(jiǎn)介:王 保(1989-),女,湖北黃岡人,助理工程師,主要從事短期天氣預(yù)報(bào)及氣候變化和小波變換等方面的研究,(電話)

    15272681090(電子信箱)wangbao1129@163.com;通訊作者,黃思先,助理工程師,主要從事短期天氣預(yù)報(bào)及重要天氣過(guò)程診斷分析

    和氣候變化等方面的研究,(電話)15926018229(電子信箱)ilyyrhsx@163.com。

    氣候變化已成為不爭(zhēng)的事實(shí)。政府間氣候變化專業(yè)委員會(huì)(IPCC)的氣候變化評(píng)估報(bào)告[1-4]指出,全球氣候正發(fā)生著以氣候變暖為主要特征的顯著性變化,全球平均地表溫度從1861年以來(lái)一直在升高。關(guān)于氣候變化對(duì)作物產(chǎn)量的影響,我國(guó)學(xué)者已經(jīng)進(jìn)行了多方面的探究[5,6],尤其是氣候變暖背景下作物產(chǎn)量的變化規(guī)律以及農(nóng)業(yè)應(yīng)對(duì)氣候變化的措施等已成為目前研究的熱點(diǎn)問(wèn)題。20世紀(jì)以來(lái),全球氣候變暖越來(lái)越快,氣候變化對(duì)作物生產(chǎn)的影響也越來(lái)越大,特別是氣候變暖背景下,極端氣候事件發(fā)生頻率、持續(xù)時(shí)間和分布規(guī)律的變化對(duì)農(nóng)田生態(tài)系統(tǒng)的影響往往超過(guò)了氣候平均變率所帶來(lái)的影響[7]。

    研究作物產(chǎn)量與氣候因子之間的關(guān)系,以往大多采用常規(guī)統(tǒng)計(jì)學(xué)方法[8-11],但氣候變化具有不同時(shí)間尺度,區(qū)域平均氣溫和降水量的周期性變化以及極端氣候事件的發(fā)生,必然引起作物產(chǎn)量的變化和波動(dòng),因此,有必要采用新的時(shí)頻分析方法對(duì)氣候變化及作物產(chǎn)量的周期性進(jìn)行研究。孫衛(wèi)國(guó)等[12]、張明等[13]和蘇占勝等[14]用功率譜和交叉小波變換方法分析了作物產(chǎn)量與區(qū)域氣候變化之間的關(guān)系,發(fā)現(xiàn)兩者關(guān)系密切,周期性變化特征顯著。

    長(zhǎng)江中下游地區(qū)是我國(guó)水稻的主要種植區(qū),氣候變化將會(huì)帶來(lái)水稻產(chǎn)量的變化,但目前氣候變化對(duì)該區(qū)水稻產(chǎn)量影響的研究較少[15-19],對(duì)該區(qū)產(chǎn)量周期性的研究更不多見(jiàn),所以,研究該區(qū)水稻產(chǎn)量與區(qū)域氣候變化的關(guān)系具有重要現(xiàn)實(shí)意義。水稻生長(zhǎng)發(fā)育要求的最低溫度在10 ℃以上,南方三季稻要求≥10 ℃的活動(dòng)積溫達(dá)到7 000 ℃,雙季稻要求達(dá)到5 300 ℃,且氣溫日較差對(duì)作物產(chǎn)量形成影響較大,因此,需要分析水稻生長(zhǎng)季內(nèi)平均氣溫、降水量、氣溫日較差和≥10 ℃的活動(dòng)積溫的影響。研究根據(jù)1951~2010年長(zhǎng)江中下游6?。ò不铡⒑?、湖南、江蘇、江西和浙江)實(shí)際水稻產(chǎn)量和同期107個(gè)氣象站的水稻生長(zhǎng)季內(nèi)平均氣溫、降水量、氣溫日較差和≥10℃的活動(dòng)積溫資料,采用交叉小波變換方法,以相對(duì)氣象產(chǎn)量為研究對(duì)象,分析長(zhǎng)江中下游地區(qū)60年水稻產(chǎn)量與生育期內(nèi)這些氣候變量的時(shí)頻變化特征及耦合振蕩的關(guān)系,討論水稻產(chǎn)量波動(dòng)與氣候因子變化周期的相關(guān)性,通過(guò)區(qū)域氣候變化對(duì)水稻產(chǎn)量的影響以及水稻產(chǎn)量波動(dòng)的原因分析,為區(qū)域性作物產(chǎn)量評(píng)估和預(yù)測(cè)提供參考。

    1 材料與方法

    1.1 材料的處理

    水稻產(chǎn)量資料來(lái)源于國(guó)家統(tǒng)計(jì)局,根據(jù)1951~2010年長(zhǎng)江中下游地區(qū)6省水稻產(chǎn)量(t/hm2),經(jīng)一致性檢驗(yàn),建立區(qū)域平均水稻產(chǎn)量序列,樣本數(shù)n=60。同期的氣候資料來(lái)源于中國(guó)氣象局逐日氣候資料,包括平均氣溫(℃)、降水量(mm)、最高氣溫(℃)和最低氣溫(℃)記錄,根據(jù)這6省107個(gè)氣象站的資料,算得區(qū)域水稻生長(zhǎng)季(4~10月)內(nèi)平均氣溫、總降水量、平均氣溫日較差以及≥10 ℃的活動(dòng)積溫。計(jì)算時(shí),為了保證資料序列的平穩(wěn)隨機(jī)過(guò)程,對(duì)資料進(jìn)行了標(biāo)準(zhǔn)化處理。

    1.2 研究方法

    1.2.1 趨勢(shì)產(chǎn)量估算 影響作物產(chǎn)量的因素比較多,且各影響因素的變化特征及其時(shí)間尺度不同,主要的影響因素有人為因素、氣象因素和隨機(jī)“噪音”三方面,分別構(gòu)成3個(gè)產(chǎn)量分量[20],作物產(chǎn)量=趨勢(shì)產(chǎn)量+氣象產(chǎn)量+隨機(jī)“噪音”。隨機(jī)“噪音”所占比例很小,一般可忽略不計(jì);趨勢(shì)產(chǎn)量是實(shí)際產(chǎn)量的平穩(wěn)變化項(xiàng),表示產(chǎn)量的演變趨勢(shì),其大小決定于歷年的生產(chǎn)技術(shù)水平,其模擬方法很多,本研究采用直線滑動(dòng)平均法[12],采用15年滑動(dòng)步長(zhǎng)來(lái)消除短周期波動(dòng)的影響,算出趨勢(shì)產(chǎn)量。氣象產(chǎn)量為作物產(chǎn)量扣除趨勢(shì)產(chǎn)量得到,反映了氣象條件對(duì)產(chǎn)量的影響。氣象產(chǎn)量與趨勢(shì)產(chǎn)量的比值可以表示為相對(duì)氣象產(chǎn)量。為了消除生產(chǎn)力水平對(duì)水稻產(chǎn)量的影響,真實(shí)地反映氣象因子對(duì)其影響規(guī)律,本研究在進(jìn)行產(chǎn)量分析時(shí)以相對(duì)氣象產(chǎn)量作為研究對(duì)象。

    1.2.2 小波變換方法 小波變換方法[21]起源于Fourier分析方法,F(xiàn)ourier級(jí)數(shù)主要用來(lái)表征信號(hào)的周期性,但在大多數(shù)情況下信號(hào)都是非周期性的,于是引入Fourier變換,將非周期性信號(hào)分解為多個(gè)周期性信號(hào)的疊加來(lái)處理,但是Fourier分析方法只是將信號(hào)在頻率域內(nèi)展開(kāi),沒(méi)有包含時(shí)域的信息,丟棄的這些時(shí)域信息可能對(duì)某些應(yīng)用非常重要,所以需要引進(jìn)小波變換方法,將信號(hào)在時(shí)頻兩域的信息提取出來(lái),充分利用信號(hào)信息[12]。

    交叉小波變換[22]是一種新型信號(hào)分析技術(shù),與傳統(tǒng)的交叉譜方法相比更具優(yōu)越性,特別是用于兩個(gè)時(shí)間序列耦合振蕩行為的相關(guān)分析上,除了可以彌補(bǔ)傳統(tǒng)交叉譜分析方法不可避免的缺陷外,還能夠?qū)⑿〔ㄗ儞Q在時(shí)頻兩域都具有表征信號(hào)局部化特征的能力發(fā)揮到極致;交叉小波變換方法具有比較強(qiáng)的耦合信號(hào)分辨的能力,具有便于描述耦合信號(hào)在時(shí)頻域中分布狀況等優(yōu)點(diǎn)[23],經(jīng)過(guò)多次實(shí)踐應(yīng)用驗(yàn)證表明[24-26],該方法確實(shí)能夠充分反映兩時(shí)間序列相關(guān)振蕩的周期顯著性和兩時(shí)間序列的相互依賴關(guān)系,還可以分析出兩時(shí)間序列信號(hào)精細(xì)的時(shí)頻變化特征,效果很顯著。

    2 結(jié)果與分析

    2.1 水稻產(chǎn)量與氣候變化特征的關(guān)系

    2.1.1 產(chǎn)量波動(dòng) 水稻產(chǎn)量的波動(dòng)與氣候條件密不可分,圖1為1951~2010年長(zhǎng)江中下游地區(qū)水稻的實(shí)際產(chǎn)量(實(shí)線)、氣象產(chǎn)量(虛線)以及趨勢(shì)產(chǎn)量(點(diǎn)劃線)的時(shí)間變化圖,實(shí)際產(chǎn)量和趨勢(shì)產(chǎn)量均為左側(cè)縱坐標(biāo)軸,氣象產(chǎn)量為右側(cè)縱坐標(biāo)軸。從圖1可以看出,長(zhǎng)江中下游地區(qū)水稻的實(shí)際產(chǎn)量隨時(shí)間呈波動(dòng)增長(zhǎng)的趨勢(shì),年代際變化主要在1985年之前氣象產(chǎn)量波動(dòng)比較大,1985年之后波動(dòng)較平穩(wěn),且1950年代和1980年代的波動(dòng)相對(duì)于1960和1970年代的波動(dòng)要大得多,這與長(zhǎng)江中下游地區(qū)水稻生長(zhǎng)季內(nèi)平均氣溫、平均氣溫日較差以及≥10 ℃活動(dòng)積溫的波動(dòng)特點(diǎn)相對(duì)應(yīng)。

    長(zhǎng)江中下游地區(qū)水稻生長(zhǎng)季內(nèi)的氣象災(zāi)害主要有高溫?zé)岷Α⒌蜏乩浜?、雨澇、干旱、春季低溫陰雨、寒露風(fēng)、臺(tái)風(fēng)以及風(fēng)雹等。1959年長(zhǎng)江中下游地區(qū)水稻減產(chǎn)主要受嚴(yán)重的伏秋連旱影響,1960年水稻低產(chǎn)主要與該區(qū)春季冷暖空氣交替頻繁,早稻育秧期出現(xiàn)低溫陰雨天氣造成大量爛秧,秧苗質(zhì)量差等有關(guān);1961年的低產(chǎn)主要由干旱造成,長(zhǎng)江中下游出現(xiàn)了大范圍、程度較重的干旱,秧田龜裂,稻禾枯黃,最終導(dǎo)致大幅減產(chǎn);1962年該區(qū)降水偏多,出現(xiàn)了雨澇,但對(duì)水稻產(chǎn)量影響顯著的還是春季的低溫連陰雨天氣;1975年的產(chǎn)量主要因“寒露風(fēng)”的影響而降低;1977年的減產(chǎn)主要受部分地區(qū)(江西等地)高溫干旱的影響;1980年長(zhǎng)江中下游地區(qū)遇到了雨澇、春季低溫陰雨以及“寒露風(fēng)”等災(zāi)害,影響了水稻的生長(zhǎng),造成大幅減產(chǎn);2003年主要是高溫?zé)岷υ斐傻臏p產(chǎn)。

    2.1.2 氣候變化 圖2為長(zhǎng)江中下游地區(qū)1951~2010年水稻生長(zhǎng)季內(nèi)平均氣溫(a)、降水量(b)、平均氣溫日較差(c)和≥10 ℃活動(dòng)積溫(d)的時(shí)間變化圖,圖中直線(或曲線)為其變化趨勢(shì)線。由圖2可以看出水稻生長(zhǎng)季內(nèi)的氣候變化特征主要表現(xiàn)為平均氣溫和≥10 ℃活動(dòng)積溫的變化呈現(xiàn)二次曲線形式,且平均氣溫轉(zhuǎn)折點(diǎn)在1975年左右,≥10 ℃活動(dòng)積溫轉(zhuǎn)折點(diǎn)在1980年左右。水稻生長(zhǎng)季內(nèi)平均氣溫、降水量、平均氣溫日較差以及≥10 ℃活動(dòng)積溫與時(shí)間的相關(guān)系數(shù)分別為0.354、0.000、-0.548和0.360,其中,水稻生長(zhǎng)季內(nèi)平均氣溫、≥10 ℃活動(dòng)積溫以及平均氣溫日較差與時(shí)間相關(guān)關(guān)系通過(guò)了0.01水平的顯著性檢驗(yàn),降水量與時(shí)間的相關(guān)性沒(méi)有通過(guò)0.05水平的檢驗(yàn)。

    2.1.3 水稻產(chǎn)量與氣候變化相關(guān)分析 由相關(guān)分析可知,水稻產(chǎn)量與水稻生長(zhǎng)季內(nèi)平均氣溫、平均氣溫日較差以及≥10℃活動(dòng)積溫的相關(guān)系數(shù)都通過(guò)了0.01水平的顯著性檢驗(yàn),與降水量的相關(guān)系數(shù)沒(méi)有通過(guò)0.05水平的顯著性檢驗(yàn)。水稻產(chǎn)量與生長(zhǎng)季內(nèi)的平均氣溫日較差呈負(fù)相關(guān),因?yàn)殚L(zhǎng)江中下游地區(qū)容易出現(xiàn)高溫天氣,雖然氣溫日較差大有利于水稻有機(jī)物的積累,但是高溫會(huì)影響水稻的生長(zhǎng),已有研究表明,水稻在孕穗至抽穗揚(yáng)花期對(duì)溫度極其敏感[27],如果日均溫度高于32 ℃,日最高溫度高于35 ℃,水稻抽穗揚(yáng)花就受到影響,從而造成產(chǎn)量損失和品質(zhì)下降;水稻在灌漿期也最易受高溫危害[28],造成結(jié)實(shí)率和千粒重的降低,從而導(dǎo)致減產(chǎn)。水稻產(chǎn)量與降水量的相關(guān)性不強(qiáng),因?yàn)樵陂L(zhǎng)江中下游地區(qū),水源較充足,農(nóng)業(yè)用水一般不受自然降水量的影響,但當(dāng)降水量過(guò)大或持續(xù)時(shí)間較長(zhǎng),仍然會(huì)影響水稻的開(kāi)花和授粉,特別是在開(kāi)花期,連陰雨天氣將會(huì)使水稻花粉的傳播受精受阻,從而導(dǎo)致減產(chǎn)。水稻產(chǎn)量與平均氣溫和≥10 ℃活動(dòng)積溫呈正相關(guān)關(guān)系,主要反映的是水稻生長(zhǎng)對(duì)熱量條件的需求,水稻是喜熱作物,熱量條件對(duì)其產(chǎn)量有重要影響,有效熱量越多,水稻產(chǎn)量越高。

    2.2 水稻產(chǎn)量與氣候變化的時(shí)頻結(jié)構(gòu)

    相關(guān)分析只能初步分析氣候變化與水稻產(chǎn)量的簡(jiǎn)單相關(guān)關(guān)系,不能充分反映區(qū)域氣候變化對(duì)水稻產(chǎn)量的影響,而小波變換則能從時(shí)間域和頻率域上充分利用時(shí)間序列信號(hào),分析出水稻產(chǎn)量和氣候變化的局部時(shí)頻變化規(guī)律,通過(guò)交叉小波分析出兩者之間在時(shí)頻結(jié)構(gòu)上的相關(guān)關(guān)系。

    小波變換系數(shù)的實(shí)數(shù)部分包含給定時(shí)間和尺度下相對(duì)于其他時(shí)間和尺度信號(hào)的強(qiáng)度和位相兩方面的信息[26]。小波系數(shù)實(shí)數(shù)部分為正時(shí),表示水稻相對(duì)氣象產(chǎn)量距平、降水量距平、平均氣溫距平、氣溫日較差距平和≥10 ℃活動(dòng)積溫距平相對(duì)偏多,圖3中用實(shí)線表示;小波系數(shù)實(shí)數(shù)部分為負(fù)時(shí),表示相對(duì)偏少,圖3中用虛線表示;小波系數(shù)為0的地方則為突變點(diǎn)。圖3中的符號(hào)反映了振蕩的位相,等值線中心反映了不同尺度振蕩的振幅最大值。采用小波功率譜檢驗(yàn)顯著性,小波功率譜圖為小波功率譜與置信水平為95%的紅噪聲總體譜的比值,比值大于1.0表示通過(guò)95%置信水平紅噪聲檢驗(yàn)的顯著周期振蕩(實(shí)線),比值小于1.0表示未通過(guò)95%置信水平的紅噪聲檢驗(yàn)(虛線)。

    2.2.1 水稻產(chǎn)量波動(dòng)的時(shí)頻變化特征 圖3為長(zhǎng)江中下游地區(qū)水稻相對(duì)氣象產(chǎn)量距平的小波變換系數(shù)實(shí)數(shù)部分等值線圖和小波功率譜圖。由圖3a可以看出,水稻相對(duì)氣象產(chǎn)量波動(dòng)的年際尺度主要表現(xiàn)在2~4年和6年尺度的周期振蕩,且時(shí)域分布不均勻,局部化特征明顯;年代際尺度主要表現(xiàn)在10年和14~16年尺度的周期振蕩。

    年際尺度特征為2~4年左右的尺度周期信號(hào)主要表現(xiàn)在1954、1955年以及1976~1982年,同時(shí)也反映了在此周期尺度上,水稻相對(duì)氣象產(chǎn)量距平在1954、1976年及1980年處于偏少期,但1955、1979年和1982年是偏多的;6年左右的周期信號(hào)主要表現(xiàn)在1976~1981年,且呈現(xiàn)偏少—偏多—偏少的振蕩規(guī)律,具體為1976年水稻相對(duì)氣象產(chǎn)量距平偏少,1977~1980年偏多,1981年又呈偏少特征。

    年代際尺度特征,10年左右的振蕩信號(hào)主要表現(xiàn)在1970~1985年,水稻相對(duì)氣象產(chǎn)量距平呈現(xiàn)偏少—偏多的兩次周期振蕩;14~16年左右的振蕩信號(hào)比較強(qiáng),等值線比較密集,在全域都有表現(xiàn),但是1985年之前的振蕩比1985年之后的強(qiáng),1985年之后的振蕩呈逐漸增強(qiáng)趨勢(shì),水稻相對(duì)氣象產(chǎn)量存在偏多—偏少的4次振蕩,具體表現(xiàn)為1950~1957年偏多,1958~1965年偏少,1966~1973年偏多,1974~1981年偏少,1982~1988年偏多,1989~1996年偏少,1997~2002年偏多,2003~2010年偏少,且等值線閉合,預(yù)計(jì)2010年之后會(huì)出現(xiàn)水稻相對(duì)氣象產(chǎn)量偏多的情況。

    由圖3b可以看出,2~4年的周期振蕩在1957年之前以及1975~1983年比較強(qiáng),結(jié)合圖3a中的結(jié)果分析,2~4年左右的周期信號(hào)在1954年、1955年和1979年的周期比較明顯,周期特征顯著;6年尺度的周期振蕩在1957~1963年以及1976~1982年比較明顯,結(jié)合圖3a可以得出,1976~1982年的周期振蕩是通過(guò)檢驗(yàn)的,周期性顯著;10年左右的周期振蕩在1970年之前比較強(qiáng),結(jié)合圖3a可知,10年左右的周期振蕩信號(hào)沒(méi)有通過(guò)顯著性檢驗(yàn);14~16年左右的周期振蕩在1980年之前比較明顯,通過(guò)了95%置信水平的紅噪聲檢驗(yàn)。

    綜上所述,長(zhǎng)江中下游地區(qū)水稻相對(duì)氣象產(chǎn)量的周期性特征主要表現(xiàn)出2~4年、6年、14~16年的周期振蕩,時(shí)域分布不均勻。

    2.2.2 氣候變化波動(dòng)的時(shí)頻變化特征 圖4為長(zhǎng)江中下游地區(qū)水稻生長(zhǎng)季內(nèi)月平均氣溫距平的小波變換系數(shù)實(shí)數(shù)部分等值線圖和小波功率譜圖。由圖4a可以看出,水稻生長(zhǎng)季內(nèi)月平均氣溫距平波動(dòng)的年際尺度主要表現(xiàn)在3~4年、6年和8年尺度的周期振蕩,且時(shí)域分布不均勻;年代際尺度主要表現(xiàn)在10和14~16年尺度的周期振蕩。

    由圖4b可以看出,3~4年尺度的周期信號(hào)在1962~1970年、1989~1997年振蕩比較強(qiáng),結(jié)合圖4a可知,3~4年尺度的周期信號(hào)在1962~1970年和1989~1997年周期性顯著;6年尺度的周期振蕩不強(qiáng),在整個(gè)時(shí)間域上都沒(méi)有通過(guò)顯著性檢驗(yàn);8年尺度周期振蕩信號(hào)在1954~1960年比較強(qiáng),結(jié)合圖3a可知,8年尺度的周期振蕩在1954~1960年通過(guò)了顯著性檢驗(yàn);10年左右的周期振蕩在2000~2008年比較強(qiáng),結(jié)合圖3a可知,10年左右的周期振蕩在2000~2008年周期性比較明顯,通過(guò)了顯著性檢驗(yàn);14~16年尺度的周期振蕩比較弱,在整個(gè)時(shí)間域上都沒(méi)有通過(guò)顯著性檢驗(yàn)(等值線的值均<1.0)。說(shuō)明長(zhǎng)江中下游地區(qū)生長(zhǎng)季內(nèi)的平均氣溫距平的周期性特征主要表現(xiàn)在3~4年、8年和10年尺度的周期振蕩,時(shí)域分布不均勻。

    對(duì)長(zhǎng)江中下游地區(qū)水稻生長(zhǎng)季內(nèi)月降水量距平的小波變換系數(shù)實(shí)數(shù)部分等值線圖和小波功率譜圖分析可以看出,水稻生長(zhǎng)季內(nèi)月降水量距平波動(dòng)的年際尺度主要表現(xiàn)在1~2年、3~4年、6年和8年尺度的周期振蕩,時(shí)域分布不均勻;年代際尺度主要表現(xiàn)在10年、14年和18年尺度的周期振蕩。1~2年尺度的周期振蕩在1952~1957年、1966~1967年、1972年、1979~1982年以及1995~2002年比較強(qiáng),其中1~2年尺度的周期信號(hào)在1952~1957年、1979~1982年以及1995~2002年的周期性比較明顯,周期特征顯著;3~4年尺度的周期信號(hào)在1951~1957年、1971~1986年以及2000~2005年的周期性比較明顯,通過(guò)了顯著性檢驗(yàn);6年尺度的周期振蕩在1966~1982年比較強(qiáng),結(jié)合小波系數(shù)實(shí)數(shù)部分分析,此周期尺度上在時(shí)間的全域范圍內(nèi)都沒(méi)有通過(guò)顯著性檢驗(yàn);8年尺度周期振蕩在1958~1961年比較強(qiáng),8年尺度周期信號(hào)在1958~1961年顯著;10年左右的周期信號(hào)在1958~1978年的振蕩比較強(qiáng),總結(jié)得出10年左右的周期信號(hào)只在1958~1978年顯著;14年以上尺度的周期振蕩比較弱,均未通過(guò)95%置信水平的紅噪聲檢驗(yàn)。說(shuō)明長(zhǎng)江中下游地區(qū)水稻生長(zhǎng)季內(nèi)的降水量距平主要表現(xiàn)在1~2年、3~4年、8年和10年尺度的周期振蕩,時(shí)域分布不均勻。

    對(duì)長(zhǎng)江中下游地區(qū)水稻生長(zhǎng)季內(nèi)氣溫日較差距平的小波變換系數(shù)實(shí)數(shù)部分和小波功率譜分析可知,水稻生長(zhǎng)季內(nèi)氣溫日較差距平波動(dòng)的年際尺度主要表現(xiàn)在1~2年、3~4年、6年和8年尺度的周期振蕩,時(shí)域分布不均勻;年代際尺度主要表現(xiàn)在14~16年尺度的周期振蕩。1~2年的周期振蕩在1955年之前、1957~1961年、1966~1967年、1972~1973年、1976~1982以及1984~1994年比較強(qiáng),其中1~2年尺度的周期信號(hào)在1955年之前、1976~1982以及1984~1994年的周期特征顯著;3~4年尺度的周期信號(hào)只在1959年之前以及1971~1975年通過(guò)了檢驗(yàn), 3~4年尺度的周期振蕩在1959年之前周期效果明顯;6年尺度的周期振蕩在1961~1985年間比較強(qiáng),通過(guò)了95%置信水平的紅噪聲檢驗(yàn);8年尺度的周期振蕩在1989年之前都比較強(qiáng),結(jié)合小波系數(shù)實(shí)數(shù)部分分析結(jié)果可知,8年尺度周期信號(hào)主要在1963年之前以及1985~1989年顯著,周期特征明顯;14~16年尺度的周期振蕩在1956~1981年間比較強(qiáng), 14~16年尺度的周期信號(hào)只在1956~1981年通過(guò)了顯著性檢驗(yàn),周期性強(qiáng)。說(shuō)明長(zhǎng)江中下游地區(qū)水稻生長(zhǎng)季內(nèi)氣溫日較差距平主要存在1~2年、3~4年、6年、8年以及14~16年尺度的周期振蕩,高頻部分比較明顯。

    對(duì)長(zhǎng)江中下游地區(qū)水稻生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫距平的小波變換系數(shù)實(shí)數(shù)部分和小波功率譜分析可知,水稻生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫距平波動(dòng)的年際尺度主要表現(xiàn)在1~2年、3~4年、6年和8年尺度的周期振蕩,時(shí)域分布不均勻;年代際尺度主要表現(xiàn)在10年和14年左右尺度的周期振蕩。生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫在高頻部分周期性比低頻明顯,具體表現(xiàn)為1~2年尺度的周期振蕩在1951~1954年、1958~1962年、1969~1972年、1974~1976年、1979~1982年以及1993~2004年比較強(qiáng),其中1~2年尺度的周期信號(hào)在1951~1954年、1969~1972年及1993~1998年比較明顯,周期特征顯著;3~4年尺度的周期信號(hào)在1962~1970年、1975~1979年、1988~1995年通過(guò)了檢驗(yàn),周期性顯著;6年尺度的周期振蕩比較弱,均未通過(guò)95%置信水平的紅噪聲檢驗(yàn);8年尺度的周期振蕩在1954~1960年比較強(qiáng),結(jié)合小波系數(shù)實(shí)數(shù)部分分析可知,8年尺度的周期信號(hào)在1954~1960年通過(guò)了顯著性檢驗(yàn);10年尺度的周期信號(hào)在2000~2008年振蕩比較強(qiáng),10年尺度的周期信號(hào)在整個(gè)時(shí)域上沒(méi)有通過(guò)顯著性檢驗(yàn);14年尺度的周期信號(hào)在全時(shí)域的振蕩都不強(qiáng),說(shuō)明該周期尺度在整個(gè)時(shí)域中都不顯著。

    長(zhǎng)江中下游地區(qū)水稻生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫距平的周期特征主要反映在高頻部分的1~2年、3~4年以及8年尺度上,低頻部分均未通過(guò)顯著性檢驗(yàn),同時(shí)將≥10 ℃活動(dòng)積溫距平的小波變換與平均氣溫距平小波變換對(duì)比發(fā)現(xiàn),它們的變換特征非常相似。

    水稻相對(duì)氣象產(chǎn)量距平小波變換的周期特征與同期生長(zhǎng)季內(nèi)的平均氣溫、降水量、氣溫日較差以及≥10℃活動(dòng)積溫距平的小波變換有某種相似性,有共同的3~4年和8年的周期振蕩,說(shuō)明水稻產(chǎn)量變化與氣象變化可能存在某種周期頻率上的聯(lián)系,因此作水稻相對(duì)氣象產(chǎn)量與各氣候變量序列之間的交叉小波協(xié)譜圖和功率譜密度圖來(lái)進(jìn)行進(jìn)一步的分析。

    2.3 水稻產(chǎn)量與區(qū)域氣候變化的相關(guān)性

    2.3.1 水稻產(chǎn)量與水稻生長(zhǎng)季內(nèi)平均氣溫的時(shí)頻變化相關(guān)性 圖5為長(zhǎng)江中下游地區(qū)水稻相對(duì)氣象產(chǎn)量距平與水稻生長(zhǎng)季內(nèi)平均氣溫距平的交叉小波協(xié)譜圖和功率譜密度圖,反映了長(zhǎng)江中下游地區(qū)水稻產(chǎn)量與水稻生長(zhǎng)季內(nèi)平均氣溫之間的時(shí)頻相關(guān)分布特征。兩者的正相關(guān)關(guān)系表明水稻生長(zhǎng)季內(nèi)平均氣溫對(duì)水稻產(chǎn)量是促進(jìn)作用,負(fù)相關(guān)則表示水稻生長(zhǎng)季內(nèi)平均氣溫出現(xiàn)不適宜水稻生產(chǎn)的情況。

    圖5a表明,相對(duì)氣象產(chǎn)量與水稻生長(zhǎng)季內(nèi)平均氣溫之間的關(guān)系以正相關(guān)為主, 在局部時(shí)域中個(gè)別頻率尺度上兩者之間也存在負(fù)相關(guān)關(guān)系。兩者正相關(guān)振蕩主要表現(xiàn)在2~3年、4年和12年尺度的共振頻率上:2~3年尺度的正相關(guān)出現(xiàn)在1952~1955年和1962~1964年;4年尺度的正相關(guān)出現(xiàn)在1973~1983年;12年尺度正相關(guān)出現(xiàn)在1967~1974年和1990~1995年。兩者負(fù)相關(guān)振蕩主要表現(xiàn)在2~3年和8年尺度上:2~3年尺度的負(fù)相關(guān)表現(xiàn)在1956~1957年以及1982~1986年;8年尺度的負(fù)相關(guān)表現(xiàn)在1950~1974年。從圖5b可以看出,水稻相對(duì)氣象產(chǎn)量與水稻生長(zhǎng)季內(nèi)平均氣溫在頻率域中年際尺度周期上的相關(guān)程度明顯高于年代際尺度周期上的相關(guān)程度,2~3年、4年以及12年尺度是呈正相關(guān)的,其中以12年時(shí)間尺度上的正相關(guān)程度最大,4年時(shí)間尺度上的次之,在8年時(shí)間尺度上的負(fù)相關(guān)程度最大,兩者之間的相關(guān)關(guān)系也隨振蕩周期尺度的不同而不同。結(jié)合圖5a可以得出,水稻相對(duì)氣象產(chǎn)量與水稻生長(zhǎng)季內(nèi)平均氣溫之間存在2~3年、4年、8年和12年尺度的共振頻率,且相關(guān)程度最好的是8年尺度的負(fù)相關(guān)和12年尺度的正相關(guān),4年尺度的正相關(guān)程度也比較明顯,2~3年尺度的相關(guān)程度不高。

    2.3.2 水稻產(chǎn)量與水稻生長(zhǎng)季內(nèi)降水量的時(shí)頻變化相關(guān)性 長(zhǎng)江中下游地區(qū)水稻相對(duì)氣象產(chǎn)量距平與水稻生長(zhǎng)季內(nèi)降水量距平的交叉小波協(xié)譜和功率譜密度反映了長(zhǎng)江中下游地區(qū)水稻產(chǎn)量與水稻生長(zhǎng)季內(nèi)降水量之間的時(shí)頻相關(guān)分布特征。由分析可知,長(zhǎng)江中下游地區(qū)水稻生長(zhǎng)季內(nèi)的降水量對(duì)水稻的產(chǎn)量影響不大,但是過(guò)量和持續(xù)的降水則會(huì)造成產(chǎn)量的降低,表現(xiàn)為兩者的負(fù)相關(guān)關(guān)系,而正常降水則會(huì)對(duì)產(chǎn)量產(chǎn)生正面影響,表現(xiàn)為兩者的正相關(guān)關(guān)系。

    通過(guò)水稻相對(duì)氣象產(chǎn)量距平與水稻生長(zhǎng)季內(nèi)降水量距平的交叉小波協(xié)譜和功率譜密度分析可知,水稻相對(duì)氣象產(chǎn)量與水稻生長(zhǎng)季內(nèi)降水量之間的關(guān)系既有正相關(guān)關(guān)系,又有負(fù)相關(guān)關(guān)系。兩者正相關(guān)主要表現(xiàn)在2~4年、6年、8年和14~16年尺度的共振頻率上:2~4年尺度的正相關(guān)出現(xiàn)在1966~1975;6年尺度的正相關(guān)出現(xiàn)在1960~1968年;8年尺度的正相關(guān)出現(xiàn)在1966~1981年;14~16年尺度正相關(guān)出現(xiàn)在1974~2010年。兩者負(fù)相關(guān)主要表現(xiàn)在2~4年、6年和14年尺度上:2~4年尺度的負(fù)相關(guān)表現(xiàn)在1958年以前以及1975~1984年;6年尺度的負(fù)相關(guān)表現(xiàn)在1960~1983年;14年尺度的負(fù)相關(guān)主要表現(xiàn)在1957~1975年。通過(guò)交叉小波功率譜密度分析可知,水稻相對(duì)氣象產(chǎn)量與水稻生長(zhǎng)季內(nèi)降水量在頻率域中年際尺度周期上的相關(guān)程度明顯高于年代際尺度周期上的相關(guān),2~4年、6年尺度是呈負(fù)相關(guān)的,8年、14~16年尺度是呈正相關(guān)的,其中以14~16年尺度上的正相關(guān)程度最大,2~4年尺度上的負(fù)相關(guān)程度最大,兩者之間的相關(guān)關(guān)系也隨振蕩周期尺度的不同而不同。綜合分析可知,水稻相對(duì)氣象產(chǎn)量與水稻生長(zhǎng)季內(nèi)降水量之間存在2~4年、6年、8年和14~16年尺度的共振頻率。

    2.3.3 水稻產(chǎn)量與水稻生長(zhǎng)季內(nèi)氣溫日較差的時(shí)頻變化相關(guān)性 長(zhǎng)江中下游地區(qū)水稻相對(duì)氣象產(chǎn)量距平與水稻生長(zhǎng)季內(nèi)氣溫日較差距平的交叉小波協(xié)譜和交叉小波功率譜密度反映了長(zhǎng)江中下游地區(qū)水稻產(chǎn)量與水稻生長(zhǎng)季內(nèi)氣溫日較差之間的時(shí)頻相關(guān)分布特征。由分析結(jié)果可知,長(zhǎng)江中下游地區(qū)水稻生長(zhǎng)季內(nèi)氣溫日較差與水稻產(chǎn)量整體是呈負(fù)相關(guān)的,因此,兩者的正相關(guān)表示水稻生長(zhǎng)季內(nèi)氣溫日較差偏小,沒(méi)有出現(xiàn)極端氣溫對(duì)水稻的傷害,有利于水稻產(chǎn)量的積累,而負(fù)相關(guān)則表示水稻生長(zhǎng)季內(nèi)氣溫日較差偏大,日最高氣溫或最低氣溫超過(guò)了水稻的適宜生長(zhǎng)溫度,不利于水稻的生長(zhǎng)發(fā)育,導(dǎo)致水稻產(chǎn)量偏低。

    通過(guò)水稻相對(duì)氣象產(chǎn)量距平與水稻生長(zhǎng)季內(nèi)氣溫日較差距平的交叉小波協(xié)譜和功率譜密度分析可知,水稻相對(duì)氣象產(chǎn)量與水稻氣溫日較差之間的關(guān)系主要表現(xiàn)為正相關(guān)關(guān)系,局部也存在有負(fù)相關(guān)關(guān)系。兩者正相關(guān)振蕩主要表現(xiàn)在2~4年、6年和14~16年尺度的共振頻率上:2~4年尺度的正相關(guān)出現(xiàn)在除1969~1974年之外的全域范圍內(nèi),且1976~1982年表現(xiàn)比較明顯;6年尺度的正相關(guān)出現(xiàn)在1968~2010年,且在1978~1980年表現(xiàn)比較明顯;14~16年尺度正相關(guān)在全域都有表現(xiàn),在1980年之前表現(xiàn)得比較強(qiáng)烈,且在1973~1979年表現(xiàn)比較明顯,等值線很密集,波幅中心出現(xiàn)在1962年左右。兩者負(fù)相關(guān)振蕩主要表現(xiàn)在2~4年、6年、8年和14年尺度上:2~4年尺度的負(fù)相關(guān)表現(xiàn)在1969~1974年,表現(xiàn)比較弱;6年尺度的負(fù)相關(guān)表現(xiàn)在1956~1967年;8年尺度的負(fù)相關(guān)主要表現(xiàn)在1986年之前,且1977年之前表現(xiàn)比較明顯,波幅中心在1955~1960年;14年尺度的負(fù)相關(guān)出現(xiàn)在1983年之后,表現(xiàn)不強(qiáng)。交叉小波功率譜密度分析表明,兩者在頻率域中年際尺度周期上的相關(guān)程度明顯低于年代際尺度周期上的相關(guān),6年、8年尺度是呈負(fù)相關(guān)的,8年左右尺度的負(fù)相關(guān)最強(qiáng),2~4年、14~16年尺度是呈正相關(guān)的,其中以14~16年尺度上的正相關(guān)程度最大,25年以后,周期尺度越大,相關(guān)性越不強(qiáng),而且兩者之間的相關(guān)關(guān)系也隨振蕩周期尺度的不同而不同。綜合分析可知,兩者之間存在2~4年、6年、8年和14~16年尺度的共振頻率,總體的相關(guān)程度是最高的,其中14~16年尺度的正相關(guān)程度最高,8年左右尺度的負(fù)相關(guān)最大,2~4年和6年的相關(guān)程度不高。

    2.3.4 水稻產(chǎn)量與水稻生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫的時(shí)頻變化相關(guān)性 長(zhǎng)江中下游地區(qū)水稻相對(duì)氣象產(chǎn)量距平與水稻生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫距平的交叉小波協(xié)譜和功率譜密度反映了長(zhǎng)江中下游地區(qū)水稻產(chǎn)量與水稻生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫之間的時(shí)頻相關(guān)分布特征。水稻生長(zhǎng)要求最低氣溫在10 ℃以上,整個(gè)生長(zhǎng)季要求≥10 ℃的活動(dòng)積溫達(dá)到一定的值,當(dāng)水稻產(chǎn)量與水稻生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫成正相關(guān)關(guān)系時(shí)表明水稻生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫達(dá)到水稻生長(zhǎng)的要求,負(fù)相關(guān)則表示水稻生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫沒(méi)有滿足水稻生長(zhǎng)的要求。

    通過(guò)水稻相對(duì)氣象產(chǎn)量距平與水稻生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫距平的交叉小波協(xié)譜和功率譜密度分析可知,兩者之間的關(guān)系以正相關(guān)為主,局部時(shí)域存在不同頻率尺度的負(fù)相關(guān)。兩者正相關(guān)振蕩主要表現(xiàn)在2~3年、4年、6年和12年尺度的共振頻率上:2~3年尺度的正相關(guān)出現(xiàn)在1955年之前、1962~1964年、1966~1975以及1990年之后,表現(xiàn)不強(qiáng);4年尺度的正相關(guān)在全域都有表現(xiàn),但在1975~1982年表現(xiàn)比較強(qiáng),波幅中心在1978年左右;6年尺度的正相關(guān)出現(xiàn)在1995年之后,表現(xiàn)很不明顯;12年尺度正相關(guān)在全域都有弱表現(xiàn)。兩者負(fù)相關(guān)振蕩主要表現(xiàn)在2~3年和8年尺度上:2~3年尺度的負(fù)相關(guān)表現(xiàn)在1957年以及1983~1986年,表現(xiàn)很弱;8年尺度的負(fù)相關(guān)表現(xiàn)在1974年之前,表現(xiàn)比較明顯,且波幅中心出現(xiàn)在1960年左右。由交叉小波功率譜密度分析可知,在頻率域中年際尺度周期上的相關(guān)高于年代際尺度周期上的相關(guān),2~3年、4年、6年、12年尺度是正相關(guān),2~3年、8年尺度是呈負(fù)相關(guān)的,其中以12年時(shí)間尺度上的正相關(guān)程度最大,8年時(shí)間尺度上的負(fù)相關(guān)程度最大,兩者之間的相關(guān)關(guān)系也隨振蕩周期尺度的不同而不同。綜合分析可知,兩者之間存在2~3年、4年、6年、8年和12年尺度的共振頻率,12年尺度的正相關(guān)程度最大,8年尺度的負(fù)相關(guān)程度最大。

    3 結(jié)論與討論

    近60年來(lái)長(zhǎng)江中下游地區(qū)的氣候變化趨勢(shì)與全國(guó)氣候變化一致,其中水稻生長(zhǎng)季內(nèi)平均氣溫和≥10 ℃活動(dòng)積溫隨時(shí)間呈二次曲線變化;氣溫日較差呈減小趨勢(shì);降水量與時(shí)間的相關(guān)性較差。水稻產(chǎn)量呈波動(dòng)增長(zhǎng)趨勢(shì),波動(dòng)特點(diǎn)明顯,出現(xiàn)了大量的豐歉年。

    1)氣候變化對(duì)長(zhǎng)江中下游地區(qū)水稻產(chǎn)量有一定的影響,水稻生長(zhǎng)季內(nèi)平均氣溫的升高和≥10 ℃活動(dòng)積溫的增加對(duì)水稻產(chǎn)量是正面影響,氣溫日較差的增大則會(huì)導(dǎo)致產(chǎn)量下降,而降水量變化則與產(chǎn)量變化關(guān)系不大。

    2)長(zhǎng)江中下游地區(qū)水稻相對(duì)氣象產(chǎn)量與水稻生長(zhǎng)季內(nèi)的平均氣溫、降水量、氣溫日較差以及≥10 ℃活動(dòng)積溫的時(shí)頻結(jié)構(gòu)有一定的相似性,存在某些頻率尺度的顯著性變化周期,主要表現(xiàn)在高頻部分。

    3)長(zhǎng)江中下游地區(qū)水稻產(chǎn)量與水稻生長(zhǎng)季內(nèi)平均氣溫和≥10 ℃活動(dòng)積溫之間存在多尺度的共振頻率,且在低頻部分的相關(guān)程度比高頻部分稍好,相關(guān)程度最好的是8年的反位相的振蕩和12年尺度的正位相振蕩。

    4)長(zhǎng)江中下游地區(qū)水稻相對(duì)氣象產(chǎn)量與水稻生長(zhǎng)季內(nèi)降水量之間存在2~4年、6年、8年和14~16年尺度的共振頻率,相關(guān)程度隨振蕩周期尺度的不同而不同,且總體相關(guān)程度較低。

    5)長(zhǎng)江中下游地區(qū)水稻相對(duì)氣象產(chǎn)量與水稻生長(zhǎng)季內(nèi)氣溫日較差之間存在2~4年、6年、8年和14~16年尺度的共振頻率,總體相關(guān)程度在所有與水稻產(chǎn)量進(jìn)行交叉小波變換的變量中最高。

    長(zhǎng)江中下游地區(qū)水熱資源豐富,氣候變化對(duì)其影響比較明顯,影響機(jī)理比較復(fù)雜,長(zhǎng)江下游地區(qū)靠近海岸的區(qū)域,受海陸氣候的影響比較大,特別是夏季太平洋西岸的熱帶高壓對(duì)其直接影響更導(dǎo)致了氣候變化的復(fù)雜性,同時(shí)城市化發(fā)展較快,城市“熱島效應(yīng)”明顯,它們對(duì)水稻產(chǎn)量的影響還有待進(jìn)一步研究,這對(duì)于估算水稻產(chǎn)量和提高水稻種植管理技術(shù)從而提高產(chǎn)量有重要意義。

    參考文獻(xiàn):

    [1] HOUGHTON J T, JENKINS G J,EPHRAUMS J J. Climate Change: The IPCC Scientific Assessment(1990). Report Prepared for Intergovernmental Panel on Climate Change by Working Group I[M]. Cambridge,United Kingdom: Cambridge University Press,1990.

    [2] HOUGHTON J T, MEIRA FILHO L G, CALLANDER B A, et al. Climate Change 1995: The Science of Climate Change.Contribution of Working Group I to the Second Assessment Report of the Intergovernmental Panel on Climate Change[M]. Cambridge,United Kingdom:Cambridge University Press,1996.

    [3] MCCARTHY J J, CANZIANI O F, LEARY N A,et al. Climate Change 2001: Impacts, Adaptation and Vulnerability.Contribution of Working Group II to the Third Assessment Report of the Intergovernmental Panel on Climate Change[M]. Cambridge,United Kingdom:Cambridge University Press,2001.

    [4] METZ B, DAVIDSON O R, BOSCH P R, et al. Climate Change 2007: Mitigation of Climate Change. Contribution of Working Group III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change [M]. Cambridge,United Kingdom:Cambridge University Press,2007.

    [5] 《氣候變化與作物產(chǎn)量》編寫組.氣候變化與作物產(chǎn)量[M]. 北京:中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社,1992.

    [6] 林而達(dá).全球氣候變化對(duì)中國(guó)農(nóng)業(yè)影響的模擬[M].北京:中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社,1997.

    [7] LOBELL D B, FIELD C B, CAHILL K N, et al. Impacts of future climate change on California perennial crop yields: model projections with climate and crop uncertainties [J]. Agricultural and Forest Meteorology,2006,141:208-218.

    [8] 賈建英,郭建平,彭 妮.氣候變化對(duì)東北地區(qū)玉米產(chǎn)量的影響[J].安徽農(nóng)業(yè)科學(xué),2010,38(32):18309-18312,18342.

    [9] 千懷遂,石艷蕊,魏東嵐,等.氣候?qū)幽鲜∶藁óa(chǎn)量的影響及其變化研究[J].生態(tài)學(xué)報(bào),2000,20(6):1061-1068.

    [10] 何永坤,王裕文.重慶市三峽庫(kù)區(qū)氣候變化對(duì)糧食產(chǎn)量的影響分析[J].山區(qū)開(kāi)發(fā),2001(12):36-39.

    [11] 扶 定,王青林,祁玉良,等.同異分析法評(píng)價(jià)氣象因子與水稻產(chǎn)量關(guān)系的初步研究[J]. 中國(guó)農(nóng)學(xué)通報(bào),2009,25(18):140-145.

    [12] 孫衛(wèi)國(guó),程炳巖,楊沈斌,等.區(qū)域氣候變化對(duì)華東地區(qū)水稻產(chǎn)量的影響[J].中國(guó)農(nóng)業(yè)氣象,2011,32(2):227-234.

    [13] 張 明,劉仁進(jìn),唐大清,等.當(dāng)陽(yáng)市水稻產(chǎn)量與當(dāng)陽(yáng)氣候變化的相關(guān)分析[J].貴州氣象,2010,34(增刊):103-105.

    [14] 蘇占勝,陳曉光,黃 峰,等.寧夏山區(qū)小麥產(chǎn)量變化特征及其氣候變化的響應(yīng)[J].干旱地區(qū)農(nóng)業(yè)研究,2007,25(2):218-224.

    [15] 楊沈斌,申雙和,趙小艷,等.氣候變化對(duì)長(zhǎng)江中下游稻區(qū)水稻產(chǎn)量的影響[J].作物學(xué)報(bào),2010,36(9):1519-1528.

    [16] 石全紅,劉建剛,陳 阜,等.長(zhǎng)江中下游地區(qū)水稻產(chǎn)量差及分布特征研究[J].中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào),2012,17(1):33-39.

    [17] 葛道闊,金之慶.氣候及其變率變化對(duì)長(zhǎng)江中下游稻區(qū)水稻生產(chǎn)的影響[J].中國(guó)水稻科學(xué),2009,23(1):57-64.

    [18] 張 倩,趙艷霞,王春乙.長(zhǎng)江中下游地區(qū)高溫?zé)岷?duì)水稻的影響[J].災(zāi)害學(xué),2011,26(4):58-62.

    [19] 劉偉昌,張雪芬,余衛(wèi)東,等.長(zhǎng)江中下游水稻高溫?zé)岷r(shí)空分布規(guī)律研究[J].安徽農(nóng)業(yè)科學(xué),2009,37(14):6454-6457.

    [20] 中國(guó)農(nóng)業(yè)科學(xué)院.中國(guó)農(nóng)業(yè)氣象學(xué)[M]. 北京:中國(guó)農(nóng)業(yè)出版社,1999.

    [21] TORRENCE C, COMPO G P. A practical guide to wavelet analysis [J]. Bulletin of the American Meteorological Society,1998,79:61-78.

    [22] PROKOPH A, BILALIH E. Cross-wavelet analysis: a tool for detection of relationships between paleo climate proxy records [J]. Mathematical Geosciences,2008,40(5):575-586.

    [23] 孫衛(wèi)國(guó),程炳巖.交叉小波變換在區(qū)域氣候分析中的應(yīng)用[J].應(yīng)用氣象學(xué)報(bào),2008,19(4):479-487.

    [24] 孫衛(wèi)國(guó),程炳巖,李 榮.黃河源區(qū)徑流量與區(qū)域氣候變化的多時(shí)間尺度相關(guān)[J].地理學(xué)報(bào),2009,64(1):117-127.

    [25] 孫衛(wèi)國(guó),程炳巖,郭 渠.華北地區(qū)降水量蒸發(fā)差與西太平洋副熱帶高壓的關(guān)系[J].中國(guó)農(nóng)業(yè)氣象,2010,31(1):41-47.

    [26] 孫衛(wèi)國(guó),程炳巖,李 榮.黃河源區(qū)徑流量的季節(jié)變化及其與區(qū)域氣候的小波相關(guān)[J].中國(guó)沙漠,2010,30(3):712-721.

    [27] 周百方.夏秋高溫?zé)岷?duì)中稻結(jié)實(shí)率的影響及應(yīng)對(duì)措施[J].作物研究,2009,23(4):237-239.

    [28] 王才林,仲維功.高溫對(duì)水稻結(jié)實(shí)率的影響及其防御對(duì)策[J].江蘇農(nóng)業(yè)科學(xué),2004(1):15-18.

    [4] METZ B, DAVIDSON O R, BOSCH P R, et al. Climate Change 2007: Mitigation of Climate Change. Contribution of Working Group III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change [M]. Cambridge,United Kingdom:Cambridge University Press,2007.

    [5] 《氣候變化與作物產(chǎn)量》編寫組.氣候變化與作物產(chǎn)量[M]. 北京:中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社,1992.

    [6] 林而達(dá).全球氣候變化對(duì)中國(guó)農(nóng)業(yè)影響的模擬[M].北京:中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社,1997.

    [7] LOBELL D B, FIELD C B, CAHILL K N, et al. Impacts of future climate change on California perennial crop yields: model projections with climate and crop uncertainties [J]. Agricultural and Forest Meteorology,2006,141:208-218.

    [8] 賈建英,郭建平,彭 妮.氣候變化對(duì)東北地區(qū)玉米產(chǎn)量的影響[J].安徽農(nóng)業(yè)科學(xué),2010,38(32):18309-18312,18342.

    [9] 千懷遂,石艷蕊,魏東嵐,等.氣候?qū)幽鲜∶藁óa(chǎn)量的影響及其變化研究[J].生態(tài)學(xué)報(bào),2000,20(6):1061-1068.

    [10] 何永坤,王裕文.重慶市三峽庫(kù)區(qū)氣候變化對(duì)糧食產(chǎn)量的影響分析[J].山區(qū)開(kāi)發(fā),2001(12):36-39.

    [11] 扶 定,王青林,祁玉良,等.同異分析法評(píng)價(jià)氣象因子與水稻產(chǎn)量關(guān)系的初步研究[J]. 中國(guó)農(nóng)學(xué)通報(bào),2009,25(18):140-145.

    [12] 孫衛(wèi)國(guó),程炳巖,楊沈斌,等.區(qū)域氣候變化對(duì)華東地區(qū)水稻產(chǎn)量的影響[J].中國(guó)農(nóng)業(yè)氣象,2011,32(2):227-234.

    [13] 張 明,劉仁進(jìn),唐大清,等.當(dāng)陽(yáng)市水稻產(chǎn)量與當(dāng)陽(yáng)氣候變化的相關(guān)分析[J].貴州氣象,2010,34(增刊):103-105.

    [14] 蘇占勝,陳曉光,黃 峰,等.寧夏山區(qū)小麥產(chǎn)量變化特征及其氣候變化的響應(yīng)[J].干旱地區(qū)農(nóng)業(yè)研究,2007,25(2):218-224.

    [15] 楊沈斌,申雙和,趙小艷,等.氣候變化對(duì)長(zhǎng)江中下游稻區(qū)水稻產(chǎn)量的影響[J].作物學(xué)報(bào),2010,36(9):1519-1528.

    [16] 石全紅,劉建剛,陳 阜,等.長(zhǎng)江中下游地區(qū)水稻產(chǎn)量差及分布特征研究[J].中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào),2012,17(1):33-39.

    [17] 葛道闊,金之慶.氣候及其變率變化對(duì)長(zhǎng)江中下游稻區(qū)水稻生產(chǎn)的影響[J].中國(guó)水稻科學(xué),2009,23(1):57-64.

    [18] 張 倩,趙艷霞,王春乙.長(zhǎng)江中下游地區(qū)高溫?zé)岷?duì)水稻的影響[J].災(zāi)害學(xué),2011,26(4):58-62.

    [19] 劉偉昌,張雪芬,余衛(wèi)東,等.長(zhǎng)江中下游水稻高溫?zé)岷r(shí)空分布規(guī)律研究[J].安徽農(nóng)業(yè)科學(xué),2009,37(14):6454-6457.

    [20] 中國(guó)農(nóng)業(yè)科學(xué)院.中國(guó)農(nóng)業(yè)氣象學(xué)[M]. 北京:中國(guó)農(nóng)業(yè)出版社,1999.

    [21] TORRENCE C, COMPO G P. A practical guide to wavelet analysis [J]. Bulletin of the American Meteorological Society,1998,79:61-78.

    [22] PROKOPH A, BILALIH E. Cross-wavelet analysis: a tool for detection of relationships between paleo climate proxy records [J]. Mathematical Geosciences,2008,40(5):575-586.

    [23] 孫衛(wèi)國(guó),程炳巖.交叉小波變換在區(qū)域氣候分析中的應(yīng)用[J].應(yīng)用氣象學(xué)報(bào),2008,19(4):479-487.

    [24] 孫衛(wèi)國(guó),程炳巖,李 榮.黃河源區(qū)徑流量與區(qū)域氣候變化的多時(shí)間尺度相關(guān)[J].地理學(xué)報(bào),2009,64(1):117-127.

    [25] 孫衛(wèi)國(guó),程炳巖,郭 渠.華北地區(qū)降水量蒸發(fā)差與西太平洋副熱帶高壓的關(guān)系[J].中國(guó)農(nóng)業(yè)氣象,2010,31(1):41-47.

    [26] 孫衛(wèi)國(guó),程炳巖,李 榮.黃河源區(qū)徑流量的季節(jié)變化及其與區(qū)域氣候的小波相關(guān)[J].中國(guó)沙漠,2010,30(3):712-721.

    [27] 周百方.夏秋高溫?zé)岷?duì)中稻結(jié)實(shí)率的影響及應(yīng)對(duì)措施[J].作物研究,2009,23(4):237-239.

    [28] 王才林,仲維功.高溫對(duì)水稻結(jié)實(shí)率的影響及其防御對(duì)策[J].江蘇農(nóng)業(yè)科學(xué),2004(1):15-18.

    [4] METZ B, DAVIDSON O R, BOSCH P R, et al. Climate Change 2007: Mitigation of Climate Change. Contribution of Working Group III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change [M]. Cambridge,United Kingdom:Cambridge University Press,2007.

    [5] 《氣候變化與作物產(chǎn)量》編寫組.氣候變化與作物產(chǎn)量[M]. 北京:中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社,1992.

    [6] 林而達(dá).全球氣候變化對(duì)中國(guó)農(nóng)業(yè)影響的模擬[M].北京:中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社,1997.

    [7] LOBELL D B, FIELD C B, CAHILL K N, et al. Impacts of future climate change on California perennial crop yields: model projections with climate and crop uncertainties [J]. Agricultural and Forest Meteorology,2006,141:208-218.

    [8] 賈建英,郭建平,彭 妮.氣候變化對(duì)東北地區(qū)玉米產(chǎn)量的影響[J].安徽農(nóng)業(yè)科學(xué),2010,38(32):18309-18312,18342.

    [9] 千懷遂,石艷蕊,魏東嵐,等.氣候?qū)幽鲜∶藁óa(chǎn)量的影響及其變化研究[J].生態(tài)學(xué)報(bào),2000,20(6):1061-1068.

    [10] 何永坤,王裕文.重慶市三峽庫(kù)區(qū)氣候變化對(duì)糧食產(chǎn)量的影響分析[J].山區(qū)開(kāi)發(fā),2001(12):36-39.

    [11] 扶 定,王青林,祁玉良,等.同異分析法評(píng)價(jià)氣象因子與水稻產(chǎn)量關(guān)系的初步研究[J]. 中國(guó)農(nóng)學(xué)通報(bào),2009,25(18):140-145.

    [12] 孫衛(wèi)國(guó),程炳巖,楊沈斌,等.區(qū)域氣候變化對(duì)華東地區(qū)水稻產(chǎn)量的影響[J].中國(guó)農(nóng)業(yè)氣象,2011,32(2):227-234.

    [13] 張 明,劉仁進(jìn),唐大清,等.當(dāng)陽(yáng)市水稻產(chǎn)量與當(dāng)陽(yáng)氣候變化的相關(guān)分析[J].貴州氣象,2010,34(增刊):103-105.

    [14] 蘇占勝,陳曉光,黃 峰,等.寧夏山區(qū)小麥產(chǎn)量變化特征及其氣候變化的響應(yīng)[J].干旱地區(qū)農(nóng)業(yè)研究,2007,25(2):218-224.

    [15] 楊沈斌,申雙和,趙小艷,等.氣候變化對(duì)長(zhǎng)江中下游稻區(qū)水稻產(chǎn)量的影響[J].作物學(xué)報(bào),2010,36(9):1519-1528.

    [16] 石全紅,劉建剛,陳 阜,等.長(zhǎng)江中下游地區(qū)水稻產(chǎn)量差及分布特征研究[J].中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào),2012,17(1):33-39.

    [17] 葛道闊,金之慶.氣候及其變率變化對(duì)長(zhǎng)江中下游稻區(qū)水稻生產(chǎn)的影響[J].中國(guó)水稻科學(xué),2009,23(1):57-64.

    [18] 張 倩,趙艷霞,王春乙.長(zhǎng)江中下游地區(qū)高溫?zé)岷?duì)水稻的影響[J].災(zāi)害學(xué),2011,26(4):58-62.

    [19] 劉偉昌,張雪芬,余衛(wèi)東,等.長(zhǎng)江中下游水稻高溫?zé)岷r(shí)空分布規(guī)律研究[J].安徽農(nóng)業(yè)科學(xué),2009,37(14):6454-6457.

    [20] 中國(guó)農(nóng)業(yè)科學(xué)院.中國(guó)農(nóng)業(yè)氣象學(xué)[M]. 北京:中國(guó)農(nóng)業(yè)出版社,1999.

    [21] TORRENCE C, COMPO G P. A practical guide to wavelet analysis [J]. Bulletin of the American Meteorological Society,1998,79:61-78.

    [22] PROKOPH A, BILALIH E. Cross-wavelet analysis: a tool for detection of relationships between paleo climate proxy records [J]. Mathematical Geosciences,2008,40(5):575-586.

    [23] 孫衛(wèi)國(guó),程炳巖.交叉小波變換在區(qū)域氣候分析中的應(yīng)用[J].應(yīng)用氣象學(xué)報(bào),2008,19(4):479-487.

    [24] 孫衛(wèi)國(guó),程炳巖,李 榮.黃河源區(qū)徑流量與區(qū)域氣候變化的多時(shí)間尺度相關(guān)[J].地理學(xué)報(bào),2009,64(1):117-127.

    [25] 孫衛(wèi)國(guó),程炳巖,郭 渠.華北地區(qū)降水量蒸發(fā)差與西太平洋副熱帶高壓的關(guān)系[J].中國(guó)農(nóng)業(yè)氣象,2010,31(1):41-47.

    [26] 孫衛(wèi)國(guó),程炳巖,李 榮.黃河源區(qū)徑流量的季節(jié)變化及其與區(qū)域氣候的小波相關(guān)[J].中國(guó)沙漠,2010,30(3):712-721.

    [27] 周百方.夏秋高溫?zé)岷?duì)中稻結(jié)實(shí)率的影響及應(yīng)對(duì)措施[J].作物研究,2009,23(4):237-239.

    [28] 王才林,仲維功.高溫對(duì)水稻結(jié)實(shí)率的影響及其防御對(duì)策[J].江蘇農(nóng)業(yè)科學(xué),2004(1):15-18.

    猜你喜歡
    水稻產(chǎn)量長(zhǎng)江中下游地區(qū)小波變換
    長(zhǎng)江中下游地區(qū)梅雨期異常年降水及大氣熱源分布特征
    青藏高原春季土壤濕度對(duì)長(zhǎng)江中下游地區(qū)初夏短期氣候影響的數(shù)值模擬
    長(zhǎng)江中下游地區(qū)一次短時(shí)大暴雨天氣的中尺度分析
    論青藏高原溫度對(duì)長(zhǎng)江中下游地區(qū)降水的影響
    亚洲欧美日韩另类电影网站| 国产精品免费大片| 99国产精品免费福利视频| 男女边吃奶边做爰视频| 亚洲av.av天堂| 街头女战士在线观看网站| 一本色道久久久久久精品综合| 国产欧美日韩综合在线一区二区| 成人国语在线视频| 国国产精品蜜臀av免费| 国产成人av激情在线播放 | 寂寞人妻少妇视频99o| 国产亚洲av片在线观看秒播厂| 伊人久久国产一区二区| 国产亚洲av片在线观看秒播厂| 国产午夜精品一二区理论片| 男女无遮挡免费网站观看| 黄色视频在线播放观看不卡| 国产成人freesex在线| 国产亚洲一区二区精品| 91精品国产九色| 亚洲精品456在线播放app| 91精品国产国语对白视频| 成人毛片60女人毛片免费| 久久99蜜桃精品久久| 日本vs欧美在线观看视频| 国产精品欧美亚洲77777| 日韩成人av中文字幕在线观看| 能在线免费看毛片的网站| 欧美日韩av久久| 精品亚洲乱码少妇综合久久| 在线精品无人区一区二区三| 制服丝袜香蕉在线| 亚洲天堂av无毛| 美女大奶头黄色视频| 制服诱惑二区| 国产精品无大码| 亚洲精品亚洲一区二区| 欧美bdsm另类| 久久久久视频综合| 欧美日韩视频高清一区二区三区二| 免费久久久久久久精品成人欧美视频 | 国产成人免费观看mmmm| 色吧在线观看| 日韩成人av中文字幕在线观看| 色婷婷久久久亚洲欧美| 久久精品国产鲁丝片午夜精品| 少妇被粗大的猛进出69影院 | 青春草视频在线免费观看| 午夜精品国产一区二区电影| 精品人妻一区二区三区麻豆| 欧美成人午夜免费资源| 老女人水多毛片| 成人二区视频| 亚洲av男天堂| xxx大片免费视频| 在线观看三级黄色| 天美传媒精品一区二区| 在线观看美女被高潮喷水网站| 在线精品无人区一区二区三| 丰满饥渴人妻一区二区三| 国产精品一区www在线观看| 国产高清有码在线观看视频| 丝袜美足系列| 插逼视频在线观看| 熟女av电影| 国产伦理片在线播放av一区| 五月开心婷婷网| 亚洲av日韩在线播放| 国产极品天堂在线| 又大又黄又爽视频免费| 五月玫瑰六月丁香| 成人影院久久| 久久97久久精品| 亚洲人与动物交配视频| 精品99又大又爽又粗少妇毛片| 久久久久久久亚洲中文字幕| 熟女av电影| 免费观看在线日韩| 七月丁香在线播放| 男人爽女人下面视频在线观看| 精品视频人人做人人爽| 蜜臀久久99精品久久宅男| 国产成人av激情在线播放 | 国产不卡av网站在线观看| 日韩一本色道免费dvd| 中国三级夫妇交换| 水蜜桃什么品种好| 伦精品一区二区三区| 视频在线观看一区二区三区| 黑人欧美特级aaaaaa片| 三上悠亚av全集在线观看| 久久女婷五月综合色啪小说| 七月丁香在线播放| 亚洲欧美日韩另类电影网站| 99re6热这里在线精品视频| av福利片在线| 亚洲欧美清纯卡通| 国产午夜精品一二区理论片| 免费播放大片免费观看视频在线观看| 丰满少妇做爰视频| 国产成人午夜福利电影在线观看| 日韩av在线免费看完整版不卡| 伦理电影免费视频| 妹子高潮喷水视频| 欧美最新免费一区二区三区| 成人国语在线视频| 欧美精品亚洲一区二区| 考比视频在线观看| 色网站视频免费| 亚洲婷婷狠狠爱综合网| 九九久久精品国产亚洲av麻豆| 国产高清不卡午夜福利| 欧美日韩在线观看h| 9色porny在线观看| 久久精品国产亚洲网站| 少妇被粗大的猛进出69影院 | 18禁动态无遮挡网站| 日韩欧美一区视频在线观看| 两个人免费观看高清视频| 最后的刺客免费高清国语| 一级毛片电影观看| 一本久久精品| 自线自在国产av| 高清视频免费观看一区二区| 两个人的视频大全免费| 亚洲国产精品专区欧美| 在线观看美女被高潮喷水网站| 国产成人aa在线观看| 伊人久久国产一区二区| www.av在线官网国产| 久久热精品热| 极品人妻少妇av视频| 一区二区三区四区激情视频| 国产精品一二三区在线看| 免费少妇av软件| 亚洲av成人精品一区久久| 人人澡人人妻人| 99九九线精品视频在线观看视频| 精品人妻熟女av久视频| 久久精品国产鲁丝片午夜精品| 久久亚洲国产成人精品v| 日本黄大片高清| 日本wwww免费看| 91在线精品国自产拍蜜月| 99热这里只有是精品在线观看| 在线观看国产h片| 一区二区三区乱码不卡18| 欧美日韩亚洲高清精品| 蜜桃在线观看..| 一区二区日韩欧美中文字幕 | 日本色播在线视频| 久久久a久久爽久久v久久| 婷婷成人精品国产| 亚洲av不卡在线观看| 熟女av电影| 亚洲激情五月婷婷啪啪| 最近最新中文字幕免费大全7| 一二三四中文在线观看免费高清| 国产精品久久久久久精品古装| 人妻 亚洲 视频| 精品熟女少妇av免费看| 亚洲国产精品专区欧美| 热re99久久国产66热| 国产精品久久久久久精品古装| 亚洲精华国产精华液的使用体验| 两个人的视频大全免费| 日韩三级伦理在线观看| 久久午夜福利片| 亚洲图色成人| 成人亚洲精品一区在线观看| 国产熟女午夜一区二区三区 | 在线亚洲精品国产二区图片欧美 | 黑人欧美特级aaaaaa片| 人妻少妇偷人精品九色| 日韩大片免费观看网站| 精品人妻偷拍中文字幕| 国产av精品麻豆| 免费av不卡在线播放| 人成视频在线观看免费观看| 午夜福利在线观看免费完整高清在| 亚洲综合精品二区| 亚洲av电影在线观看一区二区三区| 嫩草影院入口| 亚洲精品456在线播放app| 午夜福利,免费看| 亚洲av国产av综合av卡| 丰满乱子伦码专区| 曰老女人黄片| 夜夜爽夜夜爽视频| 99国产综合亚洲精品| 亚洲欧美成人综合另类久久久| 一级二级三级毛片免费看| 狂野欧美白嫩少妇大欣赏| 人妻 亚洲 视频| 久久久亚洲精品成人影院| 自拍欧美九色日韩亚洲蝌蚪91| 最近手机中文字幕大全| av黄色大香蕉| 国产熟女欧美一区二区| 亚洲欧洲国产日韩| 天美传媒精品一区二区| 美女视频免费永久观看网站| 国产毛片在线视频| 日产精品乱码卡一卡2卡三| 久久久久人妻精品一区果冻| 99热全是精品| 97超碰精品成人国产| 欧美亚洲 丝袜 人妻 在线| 男女边吃奶边做爰视频| 狂野欧美激情性bbbbbb| 青春草视频在线免费观看| 22中文网久久字幕| 精品久久久久久电影网| 大片电影免费在线观看免费| 精品少妇黑人巨大在线播放| 人妻制服诱惑在线中文字幕| 看非洲黑人一级黄片| 久久ye,这里只有精品| 久久久久精品性色| 久久精品国产a三级三级三级| 新久久久久国产一级毛片| 精品少妇黑人巨大在线播放| 九草在线视频观看| 91成人精品电影| 中文字幕精品免费在线观看视频 | 综合色丁香网| 妹子高潮喷水视频| 日本黄色日本黄色录像| 欧美bdsm另类| 亚洲高清免费不卡视频| av在线app专区| 亚洲精品美女久久av网站| 亚洲精品日韩av片在线观看| 国精品久久久久久国模美| 久热久热在线精品观看| a级片在线免费高清观看视频| 亚洲精品久久久久久婷婷小说| 丰满迷人的少妇在线观看| 一级黄片播放器| 久久久久久久久久成人| 日韩av在线免费看完整版不卡| 亚洲欧洲日产国产| 热99国产精品久久久久久7| 欧美成人精品欧美一级黄| 啦啦啦视频在线资源免费观看| 纯流量卡能插随身wifi吗| 亚洲成人av在线免费| 日日撸夜夜添| 欧美xxxx性猛交bbbb| a级片在线免费高清观看视频| 蜜桃国产av成人99| 日韩视频一区二区在线观看| 成年人黄色毛片网站| 日韩 欧美 亚洲 中文字幕| 看免费av毛片| 极品人妻少妇av视频| 欧美黄色片欧美黄色片| 精品欧美一区二区三区在线| 18禁黄网站禁片午夜丰满| 精品人妻在线不人妻| 国产有黄有色有爽视频| 午夜福利一区二区在线看| a级毛片在线看网站| 欧美av亚洲av综合av国产av| 女人被躁到高潮嗷嗷叫费观| 成人av一区二区三区在线看| 久久人人97超碰香蕉20202| 曰老女人黄片| 丝瓜视频免费看黄片| 老熟妇乱子伦视频在线观看| 日本撒尿小便嘘嘘汇集6| 精品人妻1区二区| 中文字幕人妻熟女乱码| 少妇裸体淫交视频免费看高清 | 成人精品一区二区免费| 男女高潮啪啪啪动态图| 国产免费视频播放在线视频| 国产1区2区3区精品| 国产成人免费无遮挡视频| 另类精品久久| 日韩欧美一区视频在线观看| 热re99久久精品国产66热6| 91老司机精品| 啦啦啦中文免费视频观看日本| 啦啦啦视频在线资源免费观看| 久久国产亚洲av麻豆专区| 亚洲中文字幕日韩| 国产成人系列免费观看| 男女床上黄色一级片免费看| 成人亚洲精品一区在线观看| 欧美黑人欧美精品刺激| 美女高潮到喷水免费观看| 亚洲人成77777在线视频| 黄色视频,在线免费观看| 99riav亚洲国产免费| 高清黄色对白视频在线免费看| 18禁美女被吸乳视频| 久久香蕉激情| 激情视频va一区二区三区| 欧美日韩精品网址| 国产成人系列免费观看| 国内毛片毛片毛片毛片毛片| 久久久久久亚洲精品国产蜜桃av| 国产麻豆69| 国产片内射在线| 色婷婷久久久亚洲欧美| 两性午夜刺激爽爽歪歪视频在线观看 | 久久国产精品男人的天堂亚洲| 欧美老熟妇乱子伦牲交| 久久中文字幕一级| 18禁裸乳无遮挡动漫免费视频| 老司机深夜福利视频在线观看| 无限看片的www在线观看| 欧美大码av| 午夜福利视频精品| 男男h啪啪无遮挡| 欧美人与性动交α欧美软件| 在线十欧美十亚洲十日本专区| 在线观看www视频免费| 午夜福利欧美成人| 亚洲av日韩精品久久久久久密| 男女边摸边吃奶| 日本wwww免费看| 国产亚洲精品久久久久5区| 久久久精品免费免费高清| 深夜精品福利| 黑人巨大精品欧美一区二区mp4| 欧美另类亚洲清纯唯美| 一区二区三区精品91| 午夜福利免费观看在线| 深夜精品福利| 少妇被粗大的猛进出69影院| 一区二区日韩欧美中文字幕| 新久久久久国产一级毛片| 成人国产一区最新在线观看| 男人舔女人的私密视频| 亚洲中文字幕日韩| 精品一区二区三卡| 99精国产麻豆久久婷婷| 欧美精品一区二区免费开放| 夫妻午夜视频| 老司机深夜福利视频在线观看| xxxhd国产人妻xxx| 日韩免费高清中文字幕av| 91精品国产国语对白视频| 免费av中文字幕在线| 99精品在免费线老司机午夜| 久久久久久免费高清国产稀缺| 他把我摸到了高潮在线观看 | 人人妻人人澡人人爽人人夜夜| 国产精品久久久久久人妻精品电影 | 国产日韩一区二区三区精品不卡| 久热爱精品视频在线9| 80岁老熟妇乱子伦牲交| 两性夫妻黄色片| 国产aⅴ精品一区二区三区波| 不卡一级毛片| 91精品三级在线观看| 视频在线观看一区二区三区| 国产激情久久老熟女| 国产精品熟女久久久久浪| 91成人精品电影| 精品国产一区二区三区久久久樱花| 啪啪无遮挡十八禁网站| 99久久精品国产亚洲精品| 国产成人啪精品午夜网站| 一级,二级,三级黄色视频| 后天国语完整版免费观看| 欧美久久黑人一区二区| 久久精品亚洲熟妇少妇任你| 一区福利在线观看| 露出奶头的视频| 欧美成人免费av一区二区三区 | 精品福利永久在线观看| 国产男女内射视频| 一区在线观看完整版| 午夜精品久久久久久毛片777| 一区二区三区国产精品乱码| 国产欧美日韩一区二区三区在线| 亚洲国产精品一区二区三区在线| √禁漫天堂资源中文www| 精品国产国语对白av| 丁香欧美五月| 麻豆av在线久日| 久久影院123| 黄色成人免费大全| 久久亚洲真实| 日韩一卡2卡3卡4卡2021年| 成人影院久久| 岛国在线观看网站| 亚洲国产欧美日韩在线播放| 国产深夜福利视频在线观看| 18禁国产床啪视频网站| 黄色视频,在线免费观看| 少妇的丰满在线观看| 色综合婷婷激情| 久久精品国产a三级三级三级| 免费高清在线观看日韩| 精品午夜福利视频在线观看一区 | 久久国产精品影院| 久久中文看片网| 中文字幕人妻丝袜制服| 老司机午夜福利在线观看视频 | 999久久久精品免费观看国产| 午夜免费成人在线视频| 国产一区二区三区综合在线观看| 交换朋友夫妻互换小说| 美女扒开内裤让男人捅视频| 欧美成人午夜精品| 国产欧美日韩精品亚洲av| 欧美日韩成人在线一区二区| 精品国产乱子伦一区二区三区| 国产一区二区三区视频了| www.999成人在线观看| 女警被强在线播放| 欧美日韩国产mv在线观看视频| 亚洲第一av免费看| 久久精品成人免费网站| 伦理电影免费视频| 国产视频一区二区在线看| 亚洲精品粉嫩美女一区| 免费人妻精品一区二区三区视频| av国产精品久久久久影院| 韩国精品一区二区三区| 在线看a的网站| 中文字幕最新亚洲高清| 国产成人欧美| 国产午夜精品久久久久久| 日本撒尿小便嘘嘘汇集6| 丁香六月欧美| 久久婷婷成人综合色麻豆| 黄色视频不卡| 老司机午夜十八禁免费视频| 国产日韩一区二区三区精品不卡| 亚洲国产欧美网| 免费少妇av软件| 性色av乱码一区二区三区2| 窝窝影院91人妻| 精品视频人人做人人爽| 午夜老司机福利片| 97人妻天天添夜夜摸| 亚洲人成77777在线视频| 国产人伦9x9x在线观看| 久久久久久久国产电影| 十八禁高潮呻吟视频| 啦啦啦 在线观看视频| 亚洲人成伊人成综合网2020| 国产精品av久久久久免费| 精品一区二区三区四区五区乱码| 一边摸一边做爽爽视频免费| 欧美激情高清一区二区三区| 欧美精品av麻豆av| 久久久久久久大尺度免费视频| 国产成人精品无人区| 国产成人精品久久二区二区免费| 亚洲三区欧美一区| 一夜夜www| 亚洲色图综合在线观看| 搡老岳熟女国产| 亚洲熟女精品中文字幕| 变态另类成人亚洲欧美熟女 | 十八禁网站免费在线| 午夜成年电影在线免费观看| 极品少妇高潮喷水抽搐| 视频区图区小说| 午夜视频精品福利| 男人操女人黄网站| 在线观看舔阴道视频| 男女床上黄色一级片免费看| 成人三级做爰电影| 日本五十路高清| 亚洲色图 男人天堂 中文字幕| 亚洲欧美一区二区三区久久| www.999成人在线观看| 一级毛片电影观看| 999精品在线视频| videosex国产| 一级毛片电影观看| 亚洲精品中文字幕在线视频| cao死你这个sao货| 真人做人爱边吃奶动态| 大型黄色视频在线免费观看| 国产精品久久久久久人妻精品电影 | 国产99久久九九免费精品| 男女无遮挡免费网站观看| 一进一出抽搐动态| 日韩一区二区三区影片| 91av网站免费观看| 日日摸夜夜添夜夜添小说| xxxhd国产人妻xxx| 国产野战对白在线观看| 中文欧美无线码| 国产老妇伦熟女老妇高清| 啦啦啦 在线观看视频| 国产成人啪精品午夜网站| 亚洲少妇的诱惑av| 久久九九热精品免费| 国产精品99久久99久久久不卡| 视频区欧美日本亚洲| 正在播放国产对白刺激| 在线观看人妻少妇| 欧美人与性动交α欧美精品济南到| 亚洲午夜精品一区,二区,三区| 国产免费现黄频在线看| 最新美女视频免费是黄的| 天天躁狠狠躁夜夜躁狠狠躁| 欧美亚洲 丝袜 人妻 在线| 一区二区三区乱码不卡18| 高清视频免费观看一区二区| 午夜激情久久久久久久| 99久久人妻综合| 91国产中文字幕| 免费在线观看视频国产中文字幕亚洲| 十八禁网站免费在线| 成人18禁在线播放| 国产野战对白在线观看| 国产色视频综合| 美国免费a级毛片| 99国产精品免费福利视频| 欧美中文综合在线视频| 国产野战对白在线观看| 亚洲国产av新网站| 国产精品一区二区精品视频观看| 亚洲精品中文字幕一二三四区 | 午夜福利乱码中文字幕| 久久久久国产一级毛片高清牌| 日韩制服丝袜自拍偷拍| 久久天躁狠狠躁夜夜2o2o| 日韩免费高清中文字幕av| 亚洲av第一区精品v没综合| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av电影在线进入| 在线观看免费高清a一片| 视频区欧美日本亚洲| 精品国内亚洲2022精品成人 | 又黄又粗又硬又大视频| 80岁老熟妇乱子伦牲交| 欧美日韩视频精品一区| 亚洲专区国产一区二区| 日韩视频一区二区在线观看| 国产精品久久电影中文字幕 | 国产精品一区二区免费欧美| av视频免费观看在线观看| 在线观看免费午夜福利视频| 国产精品国产高清国产av | 亚洲性夜色夜夜综合| 亚洲国产欧美网| 国产精品一区二区在线观看99| 日日爽夜夜爽网站| 美女高潮到喷水免费观看| 中文字幕av电影在线播放| 蜜桃在线观看..| 久久午夜综合久久蜜桃| 精品视频人人做人人爽| 免费久久久久久久精品成人欧美视频| 国产野战对白在线观看| 亚洲成国产人片在线观看| 十八禁高潮呻吟视频| e午夜精品久久久久久久| 女人高潮潮喷娇喘18禁视频| 99热网站在线观看| tube8黄色片| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产成人一精品久久久| 丝袜人妻中文字幕| 老司机亚洲免费影院| 亚洲第一青青草原| 美女扒开内裤让男人捅视频| 99久久人妻综合| 欧美日韩福利视频一区二区| 久久久精品区二区三区| 久久久久国内视频| 老司机深夜福利视频在线观看| 一个人免费在线观看的高清视频| 真人做人爱边吃奶动态| www.999成人在线观看| 精品少妇久久久久久888优播| 日韩欧美一区视频在线观看| 一级片'在线观看视频| 国产亚洲精品第一综合不卡| 中文欧美无线码| 老熟妇仑乱视频hdxx| 久久这里只有精品19| 久久久国产欧美日韩av| 大片电影免费在线观看免费| 80岁老熟妇乱子伦牲交| 正在播放国产对白刺激| 在线天堂中文资源库| 久久精品亚洲精品国产色婷小说| 欧美精品一区二区大全| 麻豆乱淫一区二区| 久久九九热精品免费| 中文字幕制服av| 成人黄色视频免费在线看| 丝袜在线中文字幕| 国产成人av教育| 飞空精品影院首页| 制服诱惑二区| 99久久人妻综合| 91麻豆av在线| 久久久精品免费免费高清| 两个人看的免费小视频| 岛国毛片在线播放| 亚洲av片天天在线观看| 国产男女超爽视频在线观看| 色婷婷久久久亚洲欧美| 最近最新免费中文字幕在线| 在线观看www视频免费| 久久性视频一级片| 免费看十八禁软件| 一本一本久久a久久精品综合妖精| 亚洲国产中文字幕在线视频| h视频一区二区三区| 亚洲综合色网址|