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

    西太平洋暖池海洋熱浪在2020~2022 三年拉尼娜事件爆發(fā)背景下的演變特征、爆發(fā)機(jī)制及其影響研究

    2024-05-06 06:51:00鄭飛張小娟曹庭偉
    大氣科學(xué) 2024年1期
    關(guān)鍵詞:天數(shù)通量葉綠素

    鄭飛 張小娟 曹庭偉 , 3

    1 中國科學(xué)院大氣物理研究所國際氣候與環(huán)境研究中心(ICCES), 北京 100029

    2 甘肅省水利水電勘測設(shè)計研究院有限責(zé)任公司, 蘭州 730000

    3 中國科學(xué)院大學(xué), 北京 100049

    1 引言

    全球變暖已導(dǎo)致地球氣候系統(tǒng)能量失衡,尤其是海洋上層溫度的極度升高(Cheng et al., 2019)。最近在海洋中發(fā)現(xiàn)了極端高溫現(xiàn)象,即海洋熱浪(Marine Heat wave,MHW),其是海表溫度(Sea Surface Temperature,SST)的持續(xù)偏暖事件,可持續(xù)數(shù)天到數(shù)月甚至更長時間,范圍亦能跨越數(shù)千公里的海洋(Pearce and Feng, 2013; Hobday et al.,2016)。MHW 事件在過去幾十年中已經(jīng)變得更頻繁、更廣泛和更強(qiáng)烈(Fr?licher et al., 2018; Oliver et al., 2018; Darmaraki et al., 2019; Laufk?tter et al.,2020),不斷創(chuàng)記錄的MHW 事件已經(jīng)對海洋生態(tài)系統(tǒng)和社會經(jīng)濟(jì)造成了廣泛且嚴(yán)重的損害(Garrabou et al., 2009; Oliver et al., 2017; Fr?licher and Laufk?tter, 2018; Cheung and Fr?licher, 2020; Dalton et al., 2020)。

    由 于 西 太 平 洋 暖 池(Western Pacific Warm Pool, WPWP)具有全球最高的SST 而存在強(qiáng)烈的海氣相互作用過程,其在局地乃至全球氣候系統(tǒng)中至關(guān)重要(Yan et al., 1992; 柏穎, 2016)。最新研究表明,在過去幾十年里,WPWP 區(qū)域表層和次表層的MHW 爆發(fā)頻次和持續(xù)時間明顯增加,且次表層的MHW 在La Ni?a 事件期間比El Ni?o 事件期間較強(qiáng)且更頻繁(Hu et al., 2021; Spillman et al.,2021; Holbrook et al., 2022; Zhang et al., 2022)。WPWP 區(qū)域持續(xù)增加的MHW 已經(jīng)對嚴(yán)重依賴海洋資源來維持經(jīng)濟(jì)的沿海國家及其周邊地區(qū)造成了重大影響(Spillman et al., 2021; Holbrook et al., 2022;Noh et al., 2022; Hamdeno et al., 2022)。而WPWP區(qū)域MHW 主要受與ENSO 等相關(guān)的海洋—大氣動力過程的驅(qū)動(Holbrook et al., 2019; Huang et al., 2021; Lee et al., 2022; 張 小 娟 和 鄭 飛, 2022),但目前對其的演變特征和爆發(fā)機(jī)制等的研究還相對較少。此外,考慮到在高能耗溫室氣體排放情景下,MHW 頻次的最大變化將發(fā)生在熱帶西太平洋和北冰洋(Fr?licher et al., 2018; Darmaraki et al., 2019;Qiu et al., 2021),故需要進(jìn)一步探討西太平洋暖池MHW 的演變特征、物理機(jī)制及生態(tài)影響等。

    本文首先分析了1982~2022 年WPWP 區(qū)MHW各屬性的時間序列變化、空間分布特征和季節(jié)變化,接著重點探討了在2020~2022 三年拉尼娜事件爆發(fā)背景下該區(qū)域MHW 的時空演變特征,并利用混合層熱收支方程分析了近三年MHW 的爆發(fā)機(jī)制,最后探討了其對海洋表層和次表層浮游生物濃度的影響。

    2 資料和方法

    2.1 數(shù)據(jù)說明

    本文采用的資料主要包括中國科學(xué)院大氣物理研究所研發(fā)的MHW 高分辨率數(shù)據(jù)產(chǎn)品,大氣海洋再分析資料和葉綠素濃度-a,其具體信息可參考表1。MHW 高分辨率數(shù)據(jù)產(chǎn)品是基于1982~2022 年美國國家海洋和大氣管理局(NOAA)的高分辨率(1/4°)逐日最優(yōu)插值海表溫度(Daily Optimum Interpolation Sea Surface Temperature V2,OISST)數(shù)據(jù)統(tǒng)計的MHW 日數(shù)據(jù)。產(chǎn)品在確定閾值時參照了前人的方法(Hobday et al., 2016),其對MHW的定義為:對于每個格點,日平均SST 連續(xù)至少五天高于30 年歷史基準(zhǔn)期(1983~2012)第90 百分位(閾值)的時段,定義為一次MHW 事件。這樣每個格點每天都存在一個閾值,其隨時空變化可以更真實地分析MHW 的變化特征。本文對MHW的定義與前人對MHW 的定義略有不同,主要采用了原始數(shù)據(jù)中計算出沒有任何平滑的更真實的閾值和氣候平均值,因為平滑的氣候可以人為地放大MHW 事件的強(qiáng)度并延長其持續(xù)時間,另外該研究也處理了多一天的閏年2 月 29 日的 MHW,可以用來準(zhǔn)確分析MHW 的日變化過程(Zhang et al.,2022)。葉綠素濃度-a 已被用作浮游生物量的代表(Behrenfeld et al., 2005; Blondeau-Patissier et al.,2014; Hu et al., 2022),故本文也利用其代表浮游生物量的變化。

    表1 研究所用數(shù)據(jù)Table 1 Datasets used in this study

    2.2 MHW 的定義和分析方法

    本文基于MHW 數(shù)據(jù),對WPWP 區(qū)域每個格點統(tǒng)計了MHW 特征屬性。包括MHW 每年的爆發(fā)頻次、累積強(qiáng)度(MHW 期間每日強(qiáng)度的總和)、最大強(qiáng)度(MHW 期間日強(qiáng)度的最大值)、平均強(qiáng)度(MHW 發(fā)生時期平均每天的海溫異常)、總天數(shù)(每年MHW 發(fā)生的天數(shù))和持續(xù)時間(每年MHW 發(fā)生的天數(shù)除以爆發(fā)頻次)。

    本文主要利用合成分析的方法對La Ni?a 與MHW 各屬性的關(guān)系進(jìn)行了分析(黃嘉佑, 2000),并利用奇異值分解(Singular value decomposition,SVD)方法(魏鳳英, 1999)來確定WPWP 區(qū)MHW的累積強(qiáng)度與葉綠素濃度-a 的時空相關(guān)性,本文選?。?0°S~35°N,100°~180°E)范圍內(nèi)的 MHW作為 SVD 右場,同范圍的葉綠素濃度-a 作為SVD左場進(jìn)行了 SVD 分解。

    進(jìn)一步地,為了評估2020~2022 年WPWP 區(qū)海氣熱收支和海洋動力過程對MHW 期間海溫變化的貢獻(xiàn),本文使用混合層熱收支方程進(jìn)行了海溫異常分析(Jin et al., 2006; Ren and Jin, 2013; Ren and Wang, 2020)。在WPWP 區(qū) 域 求 平 均 值(5°S~15°N, 120°E~160°E),線性化的混合層熱收支方程可寫為

    公式(1)中的平均上升流平流可以進(jìn)一步分解為以下近似(An et al., 1999; Ren and Jin, 2013):

    其中,MC 是平均環(huán)流的影響,ZA 是緯向平流反饋,EK 是Ekman 抽吸反饋,TH 是溫躍層反饋,NDH 是非線性動力加熱項, TD 是凈熱通量項,R是殘余項。

    3 西太平洋暖池MHW 變化特征

    本節(jié)首先分析了1982~2022 年西太平洋暖池平均日SST 異常和MHW 的覆蓋面積。圖1 是統(tǒng)計得到的1982~2022 年WPWP 區(qū)平均日SST 異常和MHW 覆蓋面積的時間序列。自1982 年衛(wèi)星觀測時期以來,WPWP 經(jīng)歷了顯著的表層變暖,尤其是在2000 年之后(圖1a),這可能是人為強(qiáng)迫和內(nèi)部變率的結(jié)果(Oliver et al., 2018; Laufk?tter et al., 2020)。相應(yīng)地,近年來在WPWP 區(qū)出現(xiàn)了更頻繁的、面積更大的MHW(圖1b)。在過去四十一年,WPWP 的MHW 覆蓋面積(某一天的綜合面積)呈顯著增加趨勢,其中1998 年、2010/2011 年、2013/2014 年、2016/2017 年 和2020~2022 年MHW 的覆蓋面積超過平均覆蓋面積的四倍。在這些MHW 覆蓋面積顯著的年份中,除了2014 年,其他年份都為連續(xù)La Ni?a 年。而且,最近三年的SST 正異常和MHW 覆蓋面積都達(dá)到了四十一年來的歷史最高,這可能與2020~2022 年三年La Ni?a 的發(fā)生有關(guān)。由上可知,西太平洋暖池MHW 覆蓋面積的變化與全球增暖的趨勢和連續(xù)La Ni?a 的發(fā)生密切相關(guān),本次研究主要討論連續(xù)La Ni?a 對西太平洋暖池MHW 各屬性的影響。前人研究表明全球變暖是強(qiáng)MHW 越來越頻發(fā)的重要 原 因 之 一(Oliver et al, 2018; Laufk?tter et al.,2020),故對于全球增暖趨勢和La Ni?a 的年際變率對MHW 影響的歸因分析,可在后續(xù)研究中進(jìn)一步深入探討。

    1982~2022 年WPWP 區(qū) 域 六 種MHW 屬 性(爆發(fā)頻次、累積強(qiáng)度、最大強(qiáng)度、平均強(qiáng)度、總天數(shù)和持續(xù)時間)的空間分布如圖2 所示。WPWP的MHW 爆發(fā)頻次、總天數(shù)和持續(xù)時間空間分布比較一致,其大值區(qū)主要集中在5°N 以北(圖2a,e和f)。此外,WPWP 區(qū)MHW 爆發(fā)頻次整體偏高,平均每年發(fā)生5 次左右(圖2a)。MHW 的累積強(qiáng)度存在一個明顯的大值區(qū),主要位于(2°N~15°N,125°E~140°E(圖2b)。MHW 的最大強(qiáng)度和平均強(qiáng)度也有一致的空間分布特征,大部分區(qū)域最大強(qiáng)度為1.65°C 左右,平均強(qiáng)度為1.1°C 左右。總體而言,WPWP 區(qū)年平均MHW 屬性表現(xiàn)出爆發(fā)頻次高,總天數(shù)多,平均強(qiáng)度較小的特征。

    基于以上分析,WPWP 區(qū)域在2020~2022 年MHW 覆蓋面積達(dá)到了41 年來最高。為了進(jìn)一步探究近三年WPWP 區(qū)域MHW 各屬性的異??臻g分布特征,本文分別對2020~2022 年平均MHW的爆發(fā)頻次、累積強(qiáng)度、最大強(qiáng)度、平均強(qiáng)度、持續(xù)時間和總天數(shù)計算了其基于41 年氣候態(tài)(1982~2022 年)的異常值(圖3)。如圖3 所示,WPWP的MHW 各屬性整體都有正異常,MHW 的爆發(fā)頻次、累積強(qiáng)度和總天數(shù)存在顯著的正異常大值區(qū),主要位于125°E~145°E 海域(圖3a,b,e)。其中,MHW 爆發(fā)頻次較往年增加最顯著,大部分區(qū)域異常值達(dá)到了6 次(較常年增加了一倍),表明在連續(xù)三年La Ni?a 事件背景下,WPWP 區(qū)的MHW爆發(fā)頻次顯著增加(圖3a)。同理,WPWP 大部分區(qū)域MHW 的累積強(qiáng)度異常值達(dá)到了96°C d,總天數(shù)的異常值更是達(dá)到了74 天(圖3b,e)。表明連續(xù)La Ni?a 事件的發(fā)生可促使此區(qū)域MHW變得更加持續(xù)和強(qiáng)度增強(qiáng)(圖3b)。對于MHW的最大強(qiáng)度和平均強(qiáng)度,其較往年的異常值較?。▓D3c,d)。由圖3f 所示,近三年MHW 的持續(xù)時間在其本身較大的區(qū)域有最強(qiáng)的正異常。綜上可意識到,2020~2022 年WPWP 區(qū)MHW 表現(xiàn)出較往年更頻繁和持續(xù),平均強(qiáng)度增強(qiáng)較小的特征,這可能主要與連續(xù)三年的La Ni?a 事件爆發(fā)相關(guān)。此外,前人研究表明,持續(xù)時間的MHW 的累積強(qiáng)度也會增強(qiáng),其會對海洋生物造成長期的持續(xù)的毀滅性影響(Holbrook et al., 2022)。

    圖3 2020~2022 年西太平洋暖池區(qū)年平均海洋熱浪(MHW)屬性相對于氣候態(tài)(1982~2022)的異常:(a)MHW 的爆發(fā)頻次(次數(shù));(b)累積強(qiáng)度(單位:°C d);(c)最大強(qiáng)度(單位:°C);(d)平均強(qiáng)度(單位:°C);(e)總天數(shù)(單位:d);(f)持續(xù)時間(單位:d)。打點區(qū)表示通過顯著性檢驗(異常值超過一倍標(biāo)準(zhǔn)差)。Fig.3 Anomalies of the annual mean MHW in the WPWP from 2020 to 2022 relative to the climatological mean from 1982 to 2022: (a) MHW frequency (counts); (b) cumulative intensity (units: °C d); (c) maximum intensity (units: °C); (d) intensity mean (units: °C); (e) total days (units: d);(f) durations (units: d).The dotted area indicates passing the significance test (i.e., the anomalies exceed the standard deviation).

    為了驗證WPWP 區(qū)域MHW 與多年La Ni?a的相關(guān)性,本文進(jìn)一步分析了研究區(qū)域MHW 各屬性區(qū)域平均的年變化與連續(xù)La Ni?a 事件的關(guān)系。圖4 即為MHW 各屬性區(qū)域平均的年變化,其中橙色條代表連續(xù)兩年及以上的La Ni?a 事件,藍(lán)色星號為三年連續(xù)事件。在過去41 年里,MHW 的累積強(qiáng)度、爆發(fā)頻次、總天數(shù)和持續(xù)時間在統(tǒng)計上顯著增加,且各屬性值較大的年份與連續(xù)La Ni?a 事件有較好的對應(yīng)關(guān)系,尤其是在2020~2022 年這四個屬性達(dá)到了歷史極值(圖4a,d-f)。除此之外,MHW 的累積強(qiáng)度、爆發(fā)頻次和總天數(shù)在2013/2014、2016/2017 和2020~2022 年超過了各屬性長期平均值長的兩倍(圖4a,d-e)。同時,MHW 的最大強(qiáng)度和累積強(qiáng)度與連續(xù)三年La Ni?a也有較強(qiáng)的關(guān)系,其都在1998 年達(dá)到了歷史最高(圖4b,c)。這表明近年來在連續(xù)La Ni?a 事件的背景下,WPWP 區(qū)域爆發(fā)了累積強(qiáng)度更強(qiáng)、更頻繁的、持續(xù)時間更長的、覆蓋面積更廣的MHW事件。

    圖4 1982~2022 年WPWP 區(qū)年平均MHW 屬性的時間序列:(a)MHW 的累積強(qiáng)度(單位:°C d);(b)最大強(qiáng)度(單位:°C);(c)平均強(qiáng)度(單位:°C);(d)爆發(fā)頻次(次數(shù));(e)總天數(shù)(單位:d);(f)持續(xù)時間(單位:d)。紅色柱狀條表示值超過MHW 各屬性整個時間序列平均值的兩倍,橙色條代表連續(xù)La Ni?a 年,帶菱形為三年連續(xù)事件。Fig.4 Time series of the annual mean MHW properties in the WPWP from 1982 to 2022: (a) Cumulative intensity (units: °C d); (b) maximum intensity (units: °C); (c) intensity mean (units: °C); (d) frequency; (e) total days (units: d); (f) durations of MHW (units: d).The red bars in (a)-(f) denote the values that are more than twice the averages of the full-time series, the orange bars denote the multiyear La Ni?a events, and the diamonds denote the three-year La Ni?a events.

    基于以上結(jié)果,本文研究揭示了WPWP 區(qū)域MHW 的累積強(qiáng)度、爆發(fā)頻次和總天數(shù)在近幾十年顯著增加,并且受連續(xù)La Ni?a 事件的影響更為顯著。本文進(jìn)一步分析了WPWP 區(qū)域平均月MHW累積強(qiáng)度(一個月內(nèi)MHW 總的強(qiáng)度)、爆發(fā)頻次(一個月內(nèi)發(fā)生MHW 的次數(shù))和總天數(shù)(一個月內(nèi)MHW 總?cè)諗?shù))的季節(jié)變化。可以看到,1982~2022 年WPWP 區(qū)域月MHW 累積強(qiáng)度和總天數(shù)存在季節(jié)變化,其在3 月份和9 月份達(dá)到峰值,而爆發(fā)頻次沒有明顯的季節(jié)變化(圖5a, c, e)。對于MHW 累積強(qiáng)度、爆發(fā)頻次和總天數(shù)在連續(xù)La Ni?a 年(圖4 選擇的連續(xù)La Ni?a 年)合成的季節(jié)變化,月MHW 累積強(qiáng)度和總天數(shù)顯著增加,其只在9 月份達(dá)到最大值(圖5b,d)。同理,連續(xù)La Ni?a 年月MHW 爆發(fā)頻次也有所增加,但是沒有明顯的季節(jié)變化(圖5f)。因此,WPWP 區(qū)MHW累積強(qiáng)度和總天數(shù)的季節(jié)變化受到連續(xù)La Ni?a 事件的影響。

    圖5 WPWP 區(qū)MHW 屬性的季節(jié)變化。1982 年1 月至2022 年12 月的月MHW(a)累積強(qiáng)度、(c)總天數(shù)和(e)爆發(fā)頻次的合成。(b) 連續(xù) La Ni?a 年的月MHW(b)累積強(qiáng)度、(d)總天數(shù)和(f)爆發(fā)頻次的合成。藍(lán)色誤差棒代表均值的95%置信區(qū)間。Fig.5 Seasonal variations of the MHW properties in the WPWP.A composite of the monthly MHW (a) cumulative intensity, (c) total days, and (e)frequency from January 1982 to December 2022.A composite of the monthly MHW (b) cumulative intensity, (d) total days, and (f) frequency for multiyear La Ni?a events.The blue error bars represent the confidence interval of the mean (95%).

    以上統(tǒng)計分析表明,在連續(xù)三年La Ni?a 事件的背景下,2020~2022 年WPWP 區(qū)MHW 的覆蓋面積、累積強(qiáng)度、爆發(fā)頻次和總天數(shù)都達(dá)到了41年以來的最高值,且MHW 的累積強(qiáng)度和總天數(shù)在9 月最顯著。由圖6a 可以看到,2020~2022 年在連續(xù)La Ni?a 的背景下,WPWP 區(qū)MHW 事件爆發(fā)頻次是全球最顯著的。且2020~2022 年WPWP區(qū)域MHW 事件在一年中任何季節(jié)都有發(fā)生,最顯著的事件發(fā)生在2020 年的夏季(圖6b)。

    圖6 2020~2022 年(a)MHW 爆發(fā)頻次空間分布和(b)WPWP 區(qū)域平均出現(xiàn)MHW 事件的季節(jié)變化。(a)中黑色框為WPWP 區(qū)(5°S~15°N,120°E~160°E),(b)中黑色線代表海表溫度(SST),綠色線為90%閾值,藍(lán)色線為SST 氣候態(tài)平均,橙色陰影為統(tǒng)計出現(xiàn)的MHW 事件。Fig.6 (a) Spatial distribution of the outbreaked frequency of MHW from 2020 to 2022 and (b) seasonal variations in the area mean MHW in the WPWP.The black frame in (a) represents the WPWP region (5°S-15°N, 120°E-160°E); the black line in (b) represents the sea surface temperature(SST), the green line represents the 90% threshold of SST, the blue line represents the climatic state of SST, and the orange shade represents the statistical MHW.

    4 2020~2022 年西太平洋暖池MHW的驅(qū)動機(jī)制分析

    上文分析表明,2020~2022 年在連續(xù)三年La Ni?a 事件背景下,WPWP 區(qū)域MHW 事件的部分屬性都達(dá)到了歷史極值,且爆發(fā)頻次為全球最顯著。前人研究表明全球變暖是西太暖池區(qū)MHW 頻發(fā)和變強(qiáng)的重要原因(Sen et al., 2020; 張小娟和鄭飛,2022)。三年La Ni?a 是一個年際氣候現(xiàn)象,前人研究也強(qiáng)調(diào)了其對暖池MHW 的調(diào)制作用(Huang et al., 2021; Lee et al., 2022),因此本文主要關(guān)注了這次超長La Ni?a 事件對暖池MHW 的影響。為了探討近三年WPWP 區(qū)域MHW 事件爆發(fā)的驅(qū)動機(jī)制,本小節(jié)接下來分析了WPWP 區(qū)域線性化的混合層熱收支方程(公式5.4),用于研究動力過程和熱力學(xué)過程對MHW 爆發(fā)的貢獻(xiàn)。由圖7i 可看出,2020~2022 年MHW 累積強(qiáng)度峰值出現(xiàn)在2020 年夏季(7~9 月),因此本文重點分析2020年7~9 月期間上層海洋的溫度收支過程。在最強(qiáng)MHW 事件發(fā)生前期(2020 年2~6 月),暖海溫異常開始逐步發(fā)展以儲存熱量,其主要由平均環(huán)流項異常、非線性動力加熱項異常主導(dǎo)(圖7a,e)。當(dāng)MHW 累積強(qiáng)度達(dá)到最強(qiáng)時(2020 年7 月之后),凈海表熱通量項異常、緯向平流異常和殘余項異常的貢獻(xiàn)變得更突出(圖7b,f,g)。此時海氣凈熱通量異常和緯向平流異常兩個過程共同推動了暖溫度異常的變化,而其它海洋動力項對WPWP 區(qū)域的海溫正異常貢獻(xiàn)不顯著。緯向平流項的正異常是與連續(xù)三年La Ni?a 事件背景下,東風(fēng)持續(xù)的正異常帶來的西太平洋暖池暖水堆積有關(guān)。對于Ekman 抽吸反饋、溫躍層反饋,其對WPWP區(qū)域海溫正異常貢獻(xiàn)不顯著(圖7c-e)。此外,值得注意的是,2020~2022 年凈海表熱通量項的時間序列波動與海溫傾向的序列波動基本一致,緯向平流在2020 年6 月之后也一直處于正異常狀態(tài)(圖7b, f),表明整個時期WPWP 區(qū)域的MHW事件的形成主要是由海氣凈熱通量和緯向平流異常過程所引起的。

    圖7 2020~2022 年WPWP 區(qū)域混合層熱收支和月MHW 累積強(qiáng)度。(a)平均環(huán)流影響(MC);(b)緯向平流反饋(ZA);(c)Ekman抽吸反饋(EK);(d)溫躍層反饋(TH);(e)非線性動力加熱項(NDH);(f)凈熱通量項(Qnet);(g)殘余項(Res)。(a)-(g)中黑線為混合層溫度異常收支。(h)各項反饋的占比,(i)WPWP 區(qū)月MHW 累積強(qiáng)度。Fig.7 Mixed layer heat budget and monthly MHW cumulative intensity in the WPWP from 2020 to 2022: (a) Effect of mean circulation (MC);(b) zonal advection feedback (ZA); (c) Ekman pumping feedback (EK); (d) thermocline feedback (TH); (e) nonlinear dynamical heating (NDH); (f) net heat flux (Qnet); (g) residual term (Res).The black line denotes the mixed layer temperature anomaly budget.(h) Proportion of seven feedbacks.(i) Cumulative intensity of the monthly MHW in the WPWP region.

    基于以上的分析可知,凈海表熱通量在2020~2022 年WPWP 區(qū)域MHW 事件的形成和維持中起到了關(guān)鍵作用,因此本文進(jìn)一步分析了凈海表熱通量的各項異常的貢獻(xiàn)(公式5.2)。在凈海表熱通量各項異常中,向下短波輻射和凈潛熱通量的異常值波動較大,而向下長波輻射、向上長波輻射、感熱通量和向上短波輻射異常值較?。▓D8)。其中,向下短波輻射正異常對上層海洋變暖起主導(dǎo)作用,其序列波動與海溫傾向波動較一致(圖8d)。這可能是由于大氣內(nèi)部變率導(dǎo)致降水減少,造成云量減少,使得更多的短波輻射加熱海洋(Wu et al.,2006; He et al., 2018)。2020~2022 年凈向下短波輻射正異常是凈海表熱通量中使得MHW 頻繁發(fā)生的主要驅(qū)動因素。

    5 2020~2022 年西太平洋暖池MHW對海洋生態(tài)的影響

    基于前面的研究分析,可明確在連續(xù)La Ni?a事件背景下,2020~2022 年WPWP 區(qū)域MHW 事件為1982 年以來最顯著,這期間最強(qiáng)MHW 的發(fā)生(2020 年7~9 月)主要由凈海表熱通量和海洋動力過程共同主導(dǎo)。前人研究表明,MHW 的發(fā)生通常會降低浮游植物生物量(Noh et al., 2022;Hamdeno et al., 2022),因此本文同時利用SVD 方法探討了WPWP 區(qū)域頻繁發(fā)生的MHW 對浮游生物量的影響。如圖9 所示,2020~2022 年MHW累積強(qiáng)度與海洋生態(tài)指標(biāo)海表面葉綠素-a 濃度兩者SVD 分解的第一對模態(tài)的空間分布和時間序列。可以看到,第一對SVD 模態(tài)的方差貢獻(xiàn)率為50%,且對應(yīng)時間序列的相關(guān)系數(shù)較高(0.84)(圖10c)。由圖10 可知,在拉尼娜事件的盛期,WPWP 區(qū)域MHW 事件累積強(qiáng)度呈現(xiàn)正異常分布,海表面葉綠素-a 濃度場主要呈現(xiàn)出負(fù)異常分布,二者為負(fù)相關(guān)協(xié)同變化。例如在2020 年5~11 月期間,MHW累積強(qiáng)度增加,而海表面葉綠素-a 濃度主要呈現(xiàn)降低趨勢(圖10c)。

    為了進(jìn)一步探討2020~2022 年WPWP 區(qū)域浮游生物量在次表層以及溫躍層的變化,本文選取了包含2020~2022 年WPWP 區(qū)域葉綠素-a 濃度數(shù)據(jù)的三個Argo 浮標(biāo),這三個浮標(biāo)的活動范圍為(11°N~14°N, 125°E~158°E),時 間 范 圍 是2020 年8 月至2022 年12 月。本文根據(jù)這3 個浮標(biāo)觀測的溫躍層葉綠素-a 濃度探討了WPWP 區(qū)域MHW 事件對海表層以下浮游生物量的影響。圖10呈現(xiàn)了研究區(qū)表層至220 m 溫躍層深度浮游生物的三維結(jié)構(gòu)(衛(wèi)星數(shù)據(jù)的海表面分布和3 個浮標(biāo)葉綠素-a 濃度剖面的垂直分布)及其在不同深度層的季節(jié)變化。由圖10 可以看出,近表層和混合層的浮游植物生物量較低,其葉綠素-a 濃度基本都小于0.1 mg m-3,尤其是在混合層最淺的時候。浮游生物量主要聚集在混合層下方到200 m 左右的溫躍層(0.25 mg m-3左右),且具有季節(jié)變化,其在每年的11 月左右有極值。此外,在頻繁的MHW 事件爆發(fā)背景下,近三年研究區(qū)域的浮游生物量在近表層、混合層、混合層以下以及最多層都有顯著的減少(圖10d)。2020~2022 年WPWP 區(qū)域浮游生物量整體都呈下降趨勢,表明頻繁發(fā)生的MHW事件不利于暖池區(qū)域海洋上層浮游生物的生存和生長。

    6 結(jié)論與討論

    本文分析了1982~2022 年WPWP 區(qū)域MHW事件的空間分布特征和季節(jié)—年際變化,并重點探討了MHW 在2020~2022 年連續(xù)三年La Ni?a 事件背景下頻繁爆發(fā)的驅(qū)動機(jī)制與生態(tài)影響,得到以下幾點結(jié)論:

    (1)近幾十年來,尤其是在2000 年之后,WPWP區(qū)域MHW 事件的覆蓋面積、爆發(fā)頻次、總天數(shù)、累積強(qiáng)度等屬性均呈現(xiàn)顯著增加的趨勢。與此同時,在季節(jié)—年際尺度上,WPWP 區(qū)域MHW 事件與連續(xù)La Ni?a 事件的爆發(fā)相關(guān),其在2020~2022年連續(xù)La Ni?a 事件背景下MHW 爆發(fā)頻次達(dá)到全球最高,且MHW 的覆蓋面積、爆發(fā)頻次、總天數(shù)、累積強(qiáng)度達(dá)到41 年來最高。

    (2)利用混合層熱收支方程探討了2020~2022 年WPWP 區(qū)域MHW 事件的驅(qū)動機(jī)制。最強(qiáng)MHW 事件的發(fā)生(2020 年夏)主要是由凈海表熱通量中向下短波輻射和海洋動力過程中的緯向平流導(dǎo)致。

    (3)揭示MHW 事件與海洋生態(tài)指標(biāo)葉綠素-a 濃度在時空尺度上呈現(xiàn)負(fù)相關(guān)的協(xié)同變化關(guān)系,當(dāng)MHW 為正異常分布時,海表面葉綠素-a 濃度呈現(xiàn)降低趨勢,其在2020 年7~9 月最為典型。進(jìn)一步利用浮標(biāo)表層至溫躍層葉綠素-a 濃度數(shù)據(jù)分析表明,近表層和混合層的浮游生物量較低,浮游植物生物量主要聚集在海洋混合層下方到200 m左右的溫躍層,且具有季節(jié)變化。同時,在本文研究區(qū)域,近三年浮游生物量在近表層、混合層、混合層以下都有顯著的減少,與頻繁發(fā)生的MHW 事件導(dǎo)致暖池區(qū)域海洋上層生物量整體呈下降趨勢有關(guān)。

    此外,本文初步利用混合層熱收支方程分析了MHW 期間海溫的變化,對于MHW 事件的天氣尺度成因分析需要進(jìn)一步細(xì)化。由于本文所用的海洋溫度和洋流數(shù)據(jù)是GODAS 提供的,而GODAS的海洋再分析數(shù)據(jù)是由NCEP2 提供的動量通量、熱通量等驅(qū)動出來的,其與NCEP2 相匹配,故本文對西太暖池MHW 的驅(qū)動機(jī)制分析采用的是分辨率較低的海面通量數(shù)據(jù)NCEP2,本研究計劃下一步研究利用其它較高分辨率的海面通量數(shù)據(jù)(比如ERA5 等)進(jìn)一步對比驗證MHW 期間海溫的變化。對于遙相關(guān)過程,例如PDO 對WPWP 區(qū)域MHW事件的影響也需要后續(xù)深入研究。與此同時,WPWP區(qū)域MHW 事件跟ENSO 本身關(guān)系密切,我們是否可以根據(jù)ENSO 預(yù)報系統(tǒng)開展暖池區(qū)域MHW的預(yù)測研究可做進(jìn)一步嘗試。

    猜你喜歡
    天數(shù)通量葉綠素
    本周連漲天數(shù)居前個股
    本周連跌天數(shù)居前個股
    冬小麥田N2O通量研究
    提取葉綠素
    桃樹葉綠素含量與SPAD值呈極顯著正相關(guān)
    葉綠素家族概述
    生日謎題
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    由松針制取三種葉綠素鈉鹽及其穩(wěn)定性的研究
    春、夏季長江口及鄰近海域溶解甲烷的分布與釋放通量
    亚洲综合色惰| 人妻少妇偷人精品九色| 国产成人aa在线观看| 18禁黄网站禁片午夜丰满| 欧美丝袜亚洲另类 | 国产精品1区2区在线观看.| 男人舔奶头视频| 久久亚洲精品不卡| 亚洲精品国产成人久久av| 在线观看午夜福利视频| 免费看日本二区| 成人一区二区视频在线观看| 变态另类成人亚洲欧美熟女| 午夜a级毛片| 亚洲狠狠婷婷综合久久图片| 亚洲成人久久爱视频| 精品人妻一区二区三区麻豆 | 99久久九九国产精品国产免费| 日韩人妻高清精品专区| 日日撸夜夜添| 亚洲黑人精品在线| 热99re8久久精品国产| 色尼玛亚洲综合影院| 国产午夜精品论理片| 亚洲中文字幕日韩| 久久精品久久久久久噜噜老黄 | 日韩,欧美,国产一区二区三区 | 久久精品国产亚洲网站| 久久久精品大字幕| 久久久久国产精品人妻aⅴ院| 亚洲专区中文字幕在线| 国产精品国产三级国产av玫瑰| 69人妻影院| 非洲黑人性xxxx精品又粗又长| 老女人水多毛片| 看免费成人av毛片| 人人妻人人澡欧美一区二区| 亚洲狠狠婷婷综合久久图片| 亚洲欧美激情综合另类| 大又大粗又爽又黄少妇毛片口| 婷婷精品国产亚洲av| 国产午夜精品久久久久久一区二区三区 | 日本在线视频免费播放| 亚洲va日本ⅴa欧美va伊人久久| 亚洲一区二区三区色噜噜| 精品一区二区免费观看| 久久亚洲真实| 国产av一区在线观看免费| 成人高潮视频无遮挡免费网站| 男人舔奶头视频| 麻豆国产97在线/欧美| 99久久精品国产国产毛片| 精华霜和精华液先用哪个| 亚洲精品国产成人久久av| 在线免费十八禁| 成人美女网站在线观看视频| 1000部很黄的大片| 国产成年人精品一区二区| 一边摸一边抽搐一进一小说| 最近最新中文字幕大全电影3| 欧美最新免费一区二区三区| 91麻豆精品激情在线观看国产| 嫩草影院入口| 国产综合懂色| 999久久久精品免费观看国产| 日韩欧美一区二区三区在线观看| www.www免费av| 亚洲精品乱码久久久v下载方式| 欧美绝顶高潮抽搐喷水| av中文乱码字幕在线| 欧美日韩中文字幕国产精品一区二区三区| 国内精品久久久久久久电影| 国产伦在线观看视频一区| 久久欧美精品欧美久久欧美| 久久人人精品亚洲av| 韩国av一区二区三区四区| 草草在线视频免费看| 婷婷精品国产亚洲av| h日本视频在线播放| 久久久精品大字幕| 国产综合懂色| 国产色爽女视频免费观看| 亚洲不卡免费看| 女的被弄到高潮叫床怎么办 | 悠悠久久av| 亚洲精品在线观看二区| 日韩国内少妇激情av| 看十八女毛片水多多多| 1000部很黄的大片| 日韩欧美国产一区二区入口| 亚洲第一区二区三区不卡| 波多野结衣高清无吗| 91午夜精品亚洲一区二区三区 | 热99re8久久精品国产| 国产精品一区www在线观看 | 免费看光身美女| 亚洲经典国产精华液单| 成人av一区二区三区在线看| 午夜免费成人在线视频| av在线蜜桃| 欧美绝顶高潮抽搐喷水| 国产精品国产高清国产av| 又黄又爽又刺激的免费视频.| 看十八女毛片水多多多| 精品福利观看| 村上凉子中文字幕在线| 亚洲内射少妇av| 成年人黄色毛片网站| 亚洲欧美精品综合久久99| av黄色大香蕉| 搞女人的毛片| 日日啪夜夜撸| 亚洲精品色激情综合| 免费观看人在逋| 97热精品久久久久久| 欧美3d第一页| 天堂网av新在线| 伦精品一区二区三区| 免费看光身美女| 久久精品国产亚洲av涩爱 | 久久久久久国产a免费观看| 色综合站精品国产| 99久久无色码亚洲精品果冻| 午夜老司机福利剧场| 51国产日韩欧美| 国产视频内射| 啦啦啦啦在线视频资源| 亚洲avbb在线观看| 日日夜夜操网爽| 黄色配什么色好看| 波多野结衣高清无吗| 国产高清视频在线观看网站| 啪啪无遮挡十八禁网站| 午夜福利成人在线免费观看| 麻豆久久精品国产亚洲av| 国产精品日韩av在线免费观看| 夜夜夜夜夜久久久久| 成年版毛片免费区| 五月玫瑰六月丁香| xxxwww97欧美| 全区人妻精品视频| 欧美日韩中文字幕国产精品一区二区三区| 欧美高清性xxxxhd video| 在线天堂最新版资源| 床上黄色一级片| 99热这里只有是精品在线观看| 人妻制服诱惑在线中文字幕| 哪里可以看免费的av片| 成年免费大片在线观看| 黄色日韩在线| 亚洲18禁久久av| 中文字幕精品亚洲无线码一区| 亚洲性久久影院| 搡女人真爽免费视频火全软件 | 久久99热6这里只有精品| 又黄又爽又免费观看的视频| 精品久久久久久久久av| 久久草成人影院| 国产精品久久久久久久电影| 少妇裸体淫交视频免费看高清| 欧美高清性xxxxhd video| 不卡一级毛片| 女人十人毛片免费观看3o分钟| 真实男女啪啪啪动态图| 久久久午夜欧美精品| 欧美zozozo另类| 春色校园在线视频观看| 在线播放无遮挡| 国产一区二区在线观看日韩| 51国产日韩欧美| 黄色一级大片看看| 深夜精品福利| 成熟少妇高潮喷水视频| 1000部很黄的大片| 日日干狠狠操夜夜爽| 久久人妻av系列| 99国产极品粉嫩在线观看| 日韩av在线大香蕉| 欧美一级a爱片免费观看看| 神马国产精品三级电影在线观看| 久久亚洲真实| av在线亚洲专区| 精品一区二区免费观看| 可以在线观看的亚洲视频| 日韩国内少妇激情av| 99热只有精品国产| 成年免费大片在线观看| 久久香蕉精品热| a级毛片免费高清观看在线播放| 精品99又大又爽又粗少妇毛片 | 干丝袜人妻中文字幕| 国产视频一区二区在线看| 久久久久久久午夜电影| 亚洲美女搞黄在线观看 | 亚洲内射少妇av| 久久草成人影院| 永久网站在线| 亚洲av熟女| 天堂网av新在线| 国产精品无大码| 午夜免费激情av| 韩国av一区二区三区四区| 日日撸夜夜添| 亚洲av不卡在线观看| 免费观看在线日韩| 国内毛片毛片毛片毛片毛片| av专区在线播放| 色噜噜av男人的天堂激情| 亚洲国产色片| 在线观看66精品国产| 久久久久久久午夜电影| 久久久久精品国产欧美久久久| 变态另类成人亚洲欧美熟女| 男女做爰动态图高潮gif福利片| 亚洲av第一区精品v没综合| 婷婷丁香在线五月| 欧美日韩综合久久久久久 | 午夜福利视频1000在线观看| 亚洲第一电影网av| 亚洲第一区二区三区不卡| 国产精品日韩av在线免费观看| 欧美中文日本在线观看视频| 少妇熟女aⅴ在线视频| 日韩欧美精品v在线| 深爱激情五月婷婷| 亚洲人与动物交配视频| 又黄又爽又刺激的免费视频.| 国产高清不卡午夜福利| 国产午夜精品久久久久久一区二区三区 | 国内精品久久久久精免费| 一区二区三区高清视频在线| 色吧在线观看| 免费搜索国产男女视频| 美女 人体艺术 gogo| 18+在线观看网站| 一区二区三区高清视频在线| 亚洲熟妇中文字幕五十中出| 一边摸一边抽搐一进一小说| 直男gayav资源| 婷婷亚洲欧美| www日本黄色视频网| 最近最新免费中文字幕在线| 中文字幕av在线有码专区| 草草在线视频免费看| 国产日本99.免费观看| 亚洲电影在线观看av| 校园人妻丝袜中文字幕| 成人综合一区亚洲| av国产免费在线观看| 国产欧美日韩一区二区精品| 在线天堂最新版资源| 男女之事视频高清在线观看| 人人妻人人看人人澡| 91麻豆精品激情在线观看国产| 啦啦啦啦在线视频资源| 成人毛片a级毛片在线播放| 毛片一级片免费看久久久久 | 亚洲18禁久久av| 人妻久久中文字幕网| 在线观看av片永久免费下载| 在现免费观看毛片| 看十八女毛片水多多多| 制服丝袜大香蕉在线| 国产成人a区在线观看| 3wmmmm亚洲av在线观看| 久久精品国产亚洲网站| 精品一区二区三区人妻视频| 俺也久久电影网| 免费看光身美女| 久久久久久久午夜电影| 成人美女网站在线观看视频| 国产91精品成人一区二区三区| 性欧美人与动物交配| 我的老师免费观看完整版| 夜夜爽天天搞| 搡老熟女国产l中国老女人| 天堂√8在线中文| 亚洲第一电影网av| 五月伊人婷婷丁香| 在线免费观看不下载黄p国产 | 97超级碰碰碰精品色视频在线观看| 亚洲国产高清在线一区二区三| 亚洲欧美日韩卡通动漫| 国产精品亚洲美女久久久| 国产一区二区在线av高清观看| 黄色丝袜av网址大全| 美女免费视频网站| 可以在线观看的亚洲视频| 中文亚洲av片在线观看爽| 日本a在线网址| 人人妻人人看人人澡| 亚洲18禁久久av| 久久精品国产鲁丝片午夜精品 | 国产在线精品亚洲第一网站| 99久久中文字幕三级久久日本| 久久草成人影院| 神马国产精品三级电影在线观看| 国产麻豆成人av免费视频| 一区二区三区激情视频| 中亚洲国语对白在线视频| 久久精品国产亚洲av涩爱 | 国产精品人妻久久久影院| 欧美高清性xxxxhd video| 无遮挡黄片免费观看| bbb黄色大片| 国产一级毛片七仙女欲春2| 一个人观看的视频www高清免费观看| 尾随美女入室| 极品教师在线免费播放| 精品久久久久久久人妻蜜臀av| 亚州av有码| 日本三级黄在线观看| 国产午夜精品久久久久久一区二区三区 | 欧美黑人巨大hd| 在线播放无遮挡| 嫩草影视91久久| 成年人黄色毛片网站| 国产乱人伦免费视频| 国产伦人伦偷精品视频| 乱码一卡2卡4卡精品| 嫩草影视91久久| 夜夜看夜夜爽夜夜摸| 观看免费一级毛片| 国产午夜精品久久久久久一区二区三区 | 日韩欧美在线乱码| 精品人妻1区二区| 国产伦人伦偷精品视频| 日韩大尺度精品在线看网址| 久久精品夜夜夜夜夜久久蜜豆| 国产成人aa在线观看| 在线观看66精品国产| 舔av片在线| 国产高清三级在线| 亚洲专区中文字幕在线| 成人午夜高清在线视频| 99热这里只有精品一区| 成人特级黄色片久久久久久久| 欧美一区二区国产精品久久精品| 国产亚洲欧美98| 两性午夜刺激爽爽歪歪视频在线观看| 久久婷婷人人爽人人干人人爱| www.www免费av| 国产精品电影一区二区三区| 色在线成人网| 一区二区三区免费毛片| 精品久久久噜噜| 中国美女看黄片| 国产在线男女| 国产av在哪里看| 我要搜黄色片| 欧美极品一区二区三区四区| 国产精品98久久久久久宅男小说| 91麻豆精品激情在线观看国产| 99riav亚洲国产免费| 91狼人影院| 久久久久久久久久黄片| 啦啦啦啦在线视频资源| 国产探花极品一区二区| 婷婷六月久久综合丁香| 男插女下体视频免费在线播放| 小蜜桃在线观看免费完整版高清| 国产三级在线视频| 亚洲精品影视一区二区三区av| 国产精品野战在线观看| 夜夜看夜夜爽夜夜摸| 亚洲av免费在线观看| 非洲黑人性xxxx精品又粗又长| 成人欧美大片| 自拍偷自拍亚洲精品老妇| 美女xxoo啪啪120秒动态图| 亚洲男人的天堂狠狠| 在线a可以看的网站| 亚州av有码| 精品日产1卡2卡| 人妻夜夜爽99麻豆av| 舔av片在线| 国产精品久久久久久久久免| 麻豆成人av在线观看| 波多野结衣高清作品| 国产视频内射| 欧美区成人在线视频| 国产一区二区亚洲精品在线观看| 国产精品伦人一区二区| 亚洲国产精品合色在线| 国产伦精品一区二区三区视频9| 非洲黑人性xxxx精品又粗又长| 一个人看的www免费观看视频| 国产精品,欧美在线| videossex国产| 色综合婷婷激情| 久久婷婷人人爽人人干人人爱| 99riav亚洲国产免费| 小蜜桃在线观看免费完整版高清| 他把我摸到了高潮在线观看| 成人av在线播放网站| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美三级三区| 亚洲成a人片在线一区二区| 动漫黄色视频在线观看| 亚洲无线在线观看| 日韩大尺度精品在线看网址| 欧美一区二区国产精品久久精品| 国产午夜精品论理片| 校园人妻丝袜中文字幕| 国产白丝娇喘喷水9色精品| 午夜福利在线观看吧| 女同久久另类99精品国产91| 成年人黄色毛片网站| 12—13女人毛片做爰片一| 国产成人aa在线观看| 午夜福利在线在线| 一级av片app| eeuss影院久久| 蜜桃久久精品国产亚洲av| 久久99热这里只有精品18| 国产在视频线在精品| 精品乱码久久久久久99久播| av专区在线播放| 国产精品爽爽va在线观看网站| av在线亚洲专区| 国产高清有码在线观看视频| 精品不卡国产一区二区三区| 真人一进一出gif抽搐免费| 成人毛片a级毛片在线播放| 非洲黑人性xxxx精品又粗又长| 国产精品福利在线免费观看| 特大巨黑吊av在线直播| 免费av毛片视频| 亚洲在线观看片| 亚洲精品粉嫩美女一区| 日日摸夜夜添夜夜添av毛片 | 久久久久久久久中文| 免费大片18禁| 国产黄片美女视频| 99久久精品国产国产毛片| 国产女主播在线喷水免费视频网站 | 熟女人妻精品中文字幕| 色在线成人网| 久久久久久久久久久丰满 | 精品久久久久久成人av| 99热只有精品国产| 久久久久久久久中文| 国产黄色小视频在线观看| 婷婷亚洲欧美| 日本欧美国产在线视频| 99国产极品粉嫩在线观看| 国产v大片淫在线免费观看| 99久久精品热视频| 淫妇啪啪啪对白视频| 欧美一区二区国产精品久久精品| 不卡视频在线观看欧美| 97热精品久久久久久| av在线蜜桃| 亚洲欧美激情综合另类| 成人国产麻豆网| 一本精品99久久精品77| 国产一区二区激情短视频| 人人妻人人看人人澡| 最近中文字幕高清免费大全6 | 日本免费a在线| 一个人看的www免费观看视频| 久久九九热精品免费| 大型黄色视频在线免费观看| 亚洲精品456在线播放app | 国内久久婷婷六月综合欲色啪| 悠悠久久av| 两人在一起打扑克的视频| 在线观看一区二区三区| 亚洲av成人精品一区久久| 欧美一区二区国产精品久久精品| 亚洲va日本ⅴa欧美va伊人久久| 久久久久国内视频| 午夜福利高清视频| 女人被狂操c到高潮| 国产蜜桃级精品一区二区三区| 免费在线观看日本一区| 亚洲av一区综合| 一级黄色大片毛片| 简卡轻食公司| 精华霜和精华液先用哪个| 精品人妻一区二区三区麻豆 | 一个人免费在线观看电影| 精品福利观看| 成人性生交大片免费视频hd| 一级a爱片免费观看的视频| 乱码一卡2卡4卡精品| 欧美日韩中文字幕国产精品一区二区三区| 国内精品一区二区在线观看| 国产亚洲av嫩草精品影院| 午夜福利在线在线| 国内久久婷婷六月综合欲色啪| 少妇的逼好多水| 中文资源天堂在线| 亚洲五月天丁香| 欧美一区二区精品小视频在线| 久久精品国产99精品国产亚洲性色| 丰满的人妻完整版| 日日摸夜夜添夜夜添小说| 日本免费a在线| 久久香蕉精品热| 非洲黑人性xxxx精品又粗又长| 一a级毛片在线观看| 亚洲在线自拍视频| 免费观看人在逋| 全区人妻精品视频| 久久99热6这里只有精品| 国产免费av片在线观看野外av| 免费观看精品视频网站| 日本a在线网址| 精品一区二区三区视频在线观看免费| 久久99热6这里只有精品| 他把我摸到了高潮在线观看| 色综合婷婷激情| 欧美成人a在线观看| 1000部很黄的大片| 婷婷丁香在线五月| 久久精品影院6| 天天一区二区日本电影三级| 亚洲最大成人手机在线| 深夜a级毛片| 亚洲人成网站高清观看| 免费看美女性在线毛片视频| 日本黄色视频三级网站网址| 久久99热6这里只有精品| 亚洲综合色惰| 欧美日本视频| videossex国产| 神马国产精品三级电影在线观看| 久久精品国产鲁丝片午夜精品 | 成年版毛片免费区| 国产成人aa在线观看| 日韩精品有码人妻一区| 天堂av国产一区二区熟女人妻| 色播亚洲综合网| 22中文网久久字幕| 欧美另类亚洲清纯唯美| 少妇的逼好多水| 国产精品一区二区三区四区免费观看 | 成人av在线播放网站| www日本黄色视频网| 成人综合一区亚洲| 亚洲欧美日韩卡通动漫| 午夜福利高清视频| 亚洲图色成人| 国产精品自产拍在线观看55亚洲| 九九在线视频观看精品| 欧美丝袜亚洲另类 | 在线国产一区二区在线| 欧美xxxx黑人xx丫x性爽| 久久精品人妻少妇| 久久精品国产鲁丝片午夜精品 | 亚洲精品一卡2卡三卡4卡5卡| bbb黄色大片| 99久久精品热视频| 国产亚洲精品av在线| 亚洲欧美清纯卡通| 亚洲国产精品sss在线观看| 午夜精品一区二区三区免费看| 狠狠狠狠99中文字幕| 两个人的视频大全免费| 免费黄网站久久成人精品| 成年女人永久免费观看视频| 免费av不卡在线播放| x7x7x7水蜜桃| 亚洲一级一片aⅴ在线观看| 欧美日韩国产亚洲二区| 一本久久中文字幕| 美女大奶头视频| 成人一区二区视频在线观看| 99久久无色码亚洲精品果冻| 久久欧美精品欧美久久欧美| 国产 一区 欧美 日韩| 天美传媒精品一区二区| 91久久精品电影网| 直男gayav资源| 欧美性感艳星| 欧美日韩黄片免| 人妻久久中文字幕网| 欧美日韩中文字幕国产精品一区二区三区| 亚洲专区国产一区二区| 99热这里只有是精品在线观看| 美女高潮喷水抽搐中文字幕| 天天一区二区日本电影三级| 熟女人妻精品中文字幕| 香蕉av资源在线| 搞女人的毛片| 熟女人妻精品中文字幕| 日日撸夜夜添| 伊人久久精品亚洲午夜| 国产成人aa在线观看| 久久久久久久久大av| 搞女人的毛片| 国产美女午夜福利| 色吧在线观看| 欧美日韩瑟瑟在线播放| 51国产日韩欧美| 国产精品,欧美在线| 少妇熟女aⅴ在线视频| 三级毛片av免费| 国产精品,欧美在线| 黄色配什么色好看| 非洲黑人性xxxx精品又粗又长| 99热只有精品国产| 成人国产综合亚洲| 非洲黑人性xxxx精品又粗又长| 日韩欧美三级三区| 色综合色国产| 麻豆精品久久久久久蜜桃| 亚洲18禁久久av| 一个人看的www免费观看视频| 久久99热6这里只有精品| 99在线人妻在线中文字幕| 麻豆成人av在线观看| 国产精品久久久久久久久免| 欧美日韩中文字幕国产精品一区二区三区|