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

    基于MODIS衛(wèi)星遙感數(shù)據(jù)的西北太平洋初級(jí)生產(chǎn)力與環(huán)境參數(shù)的相關(guān)性

    2020-09-10 07:22:44李璇陳文忠
    海洋開(kāi)發(fā)與管理 2020年4期

    李璇 陳文忠

    摘要:為掌握西北太平洋初級(jí)生產(chǎn)力的變化規(guī)律,科學(xué)預(yù)防赤潮災(zāi)害頻發(fā)和海洋環(huán)境污染,保障海洋漁業(yè)等產(chǎn)業(yè)的發(fā)展,文章利用2007—2018年的MODIS衛(wèi)星遙感數(shù)據(jù),分別分析西北太平洋2個(gè)典型區(qū)域的葉綠素濃度、光合有效輻射、海表溫度和凈初級(jí)生產(chǎn)力的時(shí)空變化特征及其相關(guān)性,通過(guò)相關(guān)性對(duì)比分析不同區(qū)域凈初級(jí)生產(chǎn)力的主要影響因素及其差異性,并根據(jù)異常氣候條件(“厄爾尼諾”現(xiàn)象)年份的情況驗(yàn)證上述相關(guān)性結(jié)論。研究結(jié)果表明:區(qū)域A位于中高緯度的親潮和黑潮寒暖流交界處,區(qū)域B受副熱帶環(huán)流和黑潮共同作用且臺(tái)風(fēng)較活躍,由于二者營(yíng)養(yǎng)鹽水平不同,凈初級(jí)生產(chǎn)力的主要影響因素具有差異性;由于酶的活性受海表溫度影響較大,營(yíng)養(yǎng)鹽濃度較高海域(區(qū)域A)的凈初級(jí)生產(chǎn)力主要考慮海表溫度的影響;營(yíng)養(yǎng)鹽濃度較低海域(區(qū)域B)的凈初級(jí)生產(chǎn)力受葉綠素濃度和海表溫度影響較大,其中與海表溫度相關(guān)性較高的原因主要是葉綠素濃度(浮游植物生物量)受海表溫度影響較大;2個(gè)區(qū)域的凈初級(jí)生產(chǎn)力受光合有效輻射影響較小;異常氣候條件年份(2015年)中,2個(gè)區(qū)域凈初級(jí)生產(chǎn)力和相關(guān)環(huán)境參數(shù)的變化特征可驗(yàn)證上述相關(guān)性結(jié)論。

    關(guān)鍵詞:凈初級(jí)生產(chǎn)力;葉綠素濃度;海表溫度;光合有效輻射;西北太平洋

    Abstract: In order to grasp the changes of primary productivity in the Northwest Pacific,scientific prevention of frequent red tide disasters and pollution of the ocean environment,and safeguard the development of industries such as ocean fisheries.This paper used MODIS satellite remote sensing data (2007—2018) to analyze the spatiotemporal characteristics and correlations of chlorophyll concentration,photosynthetically active radiation,sea surface temperature and net primary productivity in two typical regions of the Northwest Pacific.Through the correlation comparison,the main influencing factors and differences of net primary productivity in different areas were analyzed,and the above correlation conclusions based on the years of abnormal climate (“El Nio” phenomenon) were verified.The research results showed that: Area A was located at the junction of pro-tide and Kuroshio cold and warm currents at mid-high latitudes.Area B was affected by subtropical circulation and Kuroshio,and typhoons were active.Due to the different nutrient levels,the main factors affecting net primary productivity were different.Because enzyme activity is greatly affected by sea surface temperature,the net primary productivity in areas with high nutrient concentrations (area A) is mainly affected by sea surface temperature.Net primary productivity in areas with low nutrient concentrations (area B) is mainly affected by chlorophyll concentration and sea surface temperature,and the higher correlation with sea surface temperature is mainly due to the chlorophyll concentration (phytoplankton biomass) being greatly affected by sea surface temperature.The net primary productivity of the two areas was less affected by photosynthetically active radiation.In the year of abnormal climate (2015),the change characteristics of the net primary productivity and related environmental parameters of the two areas could verify the above-mentioned correlation conclusions.

    Key words:Net primary productivity,Chlorophyll concentration,Sea surface temperature,Photosynthetically active radiation,Northwest Pacific Ocean

    0 引言

    海洋初級(jí)生產(chǎn)力(ocean primary productivity,OPP)是指海洋底棲植物和浮游植物等自養(yǎng)生物在單位面積上光合作用的能力,一般表示為每天或每年在單位面積上所產(chǎn)生的有機(jī)物質(zhì)數(shù)量。其中,凈初級(jí)生產(chǎn)力(net primary productivity,NPP)是指自養(yǎng)生物在初級(jí)生產(chǎn)過(guò)程中有部分能量被自身消耗,剩余的用于自身生長(zhǎng)的能量。由于海洋生物的所有過(guò)程都與初級(jí)生產(chǎn)力有一定的關(guān)系,海洋初級(jí)生產(chǎn)力一直是海洋科學(xué)研究的核心和熱點(diǎn)之一。海洋初級(jí)生產(chǎn)力為海洋生態(tài)系統(tǒng)的運(yùn)轉(zhuǎn)提供能量來(lái)源,也可用于估算漁業(yè)產(chǎn)量,同時(shí)在全球碳循環(huán)中發(fā)揮十分重要的作用。隨著近年來(lái)我國(guó)赤潮災(zāi)害頻發(fā)和二氧化碳含量升高等海洋環(huán)境問(wèn)題日趨嚴(yán)重,初級(jí)生產(chǎn)力逐漸作為海洋環(huán)境評(píng)估和赤潮監(jiān)測(cè)的重要參數(shù)。

    海洋面積較大且復(fù)雜多變,獲取大范圍和長(zhǎng)時(shí)間序列數(shù)據(jù)較困難。衛(wèi)星遙感技術(shù)的發(fā)展可實(shí)現(xiàn)對(duì)高分辨率和長(zhǎng)時(shí)間序列海洋數(shù)據(jù)的獲取,其中最多應(yīng)用于水色研究的為Sea WiFS和MODIS數(shù)據(jù)。衛(wèi)星遙感所獲取的海洋信息一般基于海洋表面,對(duì)于初級(jí)生產(chǎn)力的研究應(yīng)針對(duì)多種因素進(jìn)行探討和分析。Ryther等[1]在光飽和的情況下研究初級(jí)生產(chǎn)力與葉綠素濃度之間的關(guān)系,為反演初級(jí)生產(chǎn)力模型的產(chǎn)生奠定基礎(chǔ)。Behrenfeld等[2]通過(guò)分析大量的實(shí)測(cè)數(shù)據(jù),發(fā)現(xiàn)可利用葉綠素濃度、海表溫度和光照條件等衛(wèi)星數(shù)據(jù)進(jìn)行反演,從而得到初級(jí)生產(chǎn)力產(chǎn)品,并建立VGPM模型;該模型中的重要參數(shù)都可通過(guò)衛(wèi)星數(shù)據(jù)獲得,可幫助研究者較便捷地獲取初級(jí)生產(chǎn)力的分布特點(diǎn),此后對(duì)于反演方法的修改和優(yōu)化大部分基于該模型;由該模型的參數(shù)可以看出,影響海洋浮游植物初級(jí)生產(chǎn)力的因素有很多,包括海水營(yíng)養(yǎng)鹽濃度、光照條件、溫度、葉綠素濃度、海水透明度和生物攝食等。

    國(guó)內(nèi)外學(xué)者對(duì)于初級(jí)生產(chǎn)力的影響因素進(jìn)行大量研究。陳興群等[3]對(duì)太平洋部分海域的初級(jí)生產(chǎn)力和葉綠素濃度進(jìn)行調(diào)查分析,發(fā)現(xiàn)葉綠素濃度與溫度和鹽度都存在一定的相關(guān)性,且不同區(qū)域的相關(guān)性不同;潘友聯(lián)[4]對(duì)葉綠素a濃度和初級(jí)生產(chǎn)力進(jìn)行研究,發(fā)現(xiàn)不同區(qū)域和不同季度的初級(jí)生產(chǎn)力變化很大,并不能單純地用葉綠素濃度表征海域的初級(jí)生產(chǎn)力,且由于影響初級(jí)生產(chǎn)力的環(huán)境因子有很多,葉綠素濃度與初級(jí)生產(chǎn)力之間雖有很強(qiáng)的相關(guān)性但有時(shí)截然不同;林智濤[5]研究我國(guó)南海風(fēng)場(chǎng)和溫度對(duì)于初級(jí)生產(chǎn)力的影響,發(fā)現(xiàn)風(fēng)場(chǎng)、海表溫度和葉綠素濃度與浮游植物的初級(jí)生產(chǎn)力具有很強(qiáng)的相關(guān)性;叢丕福[6]對(duì)我國(guó)陸架海整體進(jìn)行海洋葉綠素和初級(jí)生產(chǎn)力的研究,發(fā)現(xiàn)二者都具有一定的季節(jié)變化特征,且葉綠素濃度與海溫呈顯著負(fù)相關(guān);Ryther等[1]最早提出浮游植物進(jìn)行光合作用的速率與葉綠素濃度之間在一定條件下具有一定的相關(guān)性,并提出表征二者關(guān)系的簡(jiǎn)單經(jīng)驗(yàn)?zāi)J剑嘘P(guān)“從生物量到生產(chǎn)力”的研究模式也由此展開(kāi);Carr等[7]研究表明控制NPP模型之間差異性和相似性的2個(gè)主要因素是輸入葉綠素?cái)?shù)據(jù)的選擇以及與光合速率的相關(guān)度,初級(jí)生產(chǎn)力與葉綠素濃度和光照條件之間具有一定的相關(guān)性;Gong等[8]對(duì)我國(guó)東海的初級(jí)生產(chǎn)力和葉綠素濃度進(jìn)行分析,發(fā)現(xiàn)西北部海域具有明顯的季節(jié)變化特征,但東南部海域的季節(jié)變化不明顯,因此初級(jí)生產(chǎn)力和葉綠素濃度之間的關(guān)系對(duì)于不同海域來(lái)說(shuō)是不固定的。

    光合作用是浮游植物在葉綠素、酶和一定的光照條件下,將自然界中存在的無(wú)機(jī)碳合成為自身有機(jī)碳的過(guò)程。因此,光照條件、葉綠素濃度和酶的活性對(duì)于光合速率有很大的影響,而溫度是影響酶的活性的重要因子。目前國(guó)內(nèi)外對(duì)于初級(jí)生產(chǎn)力影響因素的研究主要集中于葉綠素濃度,且研究?jī)?nèi)容較單一,主要表明上述3種因子對(duì)于初級(jí)生產(chǎn)力具有一定的影響性,但并未將3種影響因子進(jìn)行對(duì)比以及系統(tǒng)研究和總結(jié),也很少對(duì)每種情況下所得出的結(jié)果進(jìn)行驗(yàn)證并分析主要影響因子,且大部分研究對(duì)象是整個(gè)海域。本研究利用MODIS衛(wèi)星遙感產(chǎn)品,研究葉綠素濃度、光合有效輻射和海表溫度3種環(huán)境參數(shù)與西北太平洋不同海域凈初級(jí)生產(chǎn)力的相關(guān)性,通過(guò)對(duì)比得出主要影響因素并分析其差異性,同時(shí)選取異常氣候條件(“厄爾尼諾”現(xiàn)象)的年份對(duì)上述結(jié)論進(jìn)行驗(yàn)證。

    1 數(shù)據(jù)和研究區(qū)域

    1.1 數(shù)據(jù)

    中分辨率成像光譜儀(MODIS)是搭載在TERRA和AQUA衛(wèi)星上的用于監(jiān)測(cè)全球環(huán)境變化的重要儀器,其在1~2 d內(nèi)可監(jiān)測(cè)整個(gè)地球,光譜分布為0.405~14.385 μm,共有36個(gè)光譜波段,分別應(yīng)用于大氣、陸地和海洋水色遙感。MODIS屬于高信噪比儀器,具有數(shù)據(jù)存儲(chǔ)高效、信息豐富、數(shù)據(jù)獲取快和覆蓋范圍廣等特點(diǎn),有非常大的應(yīng)用價(jià)值,在海洋環(huán)境監(jiān)測(cè)和研究中得到廣泛應(yīng)用。按照處理級(jí)別,MODIS數(shù)據(jù)產(chǎn)品分為0級(jí)、1級(jí)、2級(jí)、3級(jí)、4級(jí)和5級(jí)(以上)。

    本研究選取的數(shù)據(jù)為葉綠素a濃度(CHL)、光合有效輻射(PAR)、海表溫度(SST)和凈初級(jí)生產(chǎn)力(NPP),其中凈初級(jí)生產(chǎn)力的MODIS衛(wèi)星數(shù)據(jù)主要來(lái)源于美國(guó)俄勒岡大學(xué)的官方網(wǎng)站。數(shù)據(jù)模型為AQUA/MODIS反演的Level-3月均產(chǎn)品,其中凈初級(jí)生產(chǎn)力的數(shù)據(jù)是由Standard VGPM模型反演得到的產(chǎn)品。該模型的計(jì)算公式為:

    1.2 研究區(qū)域

    西北太平洋具有較復(fù)雜的海洋環(huán)境,對(duì)于全球氣候變化、碳循環(huán)和海洋漁業(yè)等都具有一定的影響[9]。我國(guó)位于太平洋西岸,是受西北太平洋影響的國(guó)家之一,因此對(duì)于西北太平洋凈初級(jí)生產(chǎn)力的研究十分重要[10]。本研究選取西北太平洋部分海域?yàn)檠芯繀^(qū)域,分別為區(qū)域A(40°~50°N,150°~170°E)和區(qū)域B(20°~30°N,130°~150°E)(圖1)。

    區(qū)域A位于西北太平洋中高緯度的親潮區(qū)域,親潮為北部寒流,北赤道流的黑潮暖流和親潮寒流交匯于此。親潮具有高含氧量和高含鹽量的特點(diǎn),大量營(yíng)養(yǎng)鹽隨黑潮涌來(lái)形成漁場(chǎng),因此該區(qū)域的海表溫度變化較顯著,且營(yíng)養(yǎng)鹽濃度較高,有利于研究海表溫度和凈初級(jí)生產(chǎn)力的相關(guān)性。楊海燕等[11]研究發(fā)現(xiàn)該區(qū)域的營(yíng)養(yǎng)鹽濃度一般在5~20 μmol/L,在西北太平洋屬于營(yíng)養(yǎng)鹽較高的區(qū)域。

    區(qū)域B為副熱帶環(huán)流區(qū)的黑潮區(qū)域,沿途海溫通常受黑潮的影響,因此溫度對(duì)其凈初級(jí)生產(chǎn)力有一定的影響。該區(qū)域緯度較低、溫度較高以及營(yíng)養(yǎng)鹽濃度較低,營(yíng)養(yǎng)鹽濃度常年低于0.5 μmol/L,冬季對(duì)流作用加強(qiáng),營(yíng)養(yǎng)鹽濃度略有增大[11]。葉綠素濃度是營(yíng)養(yǎng)鹽的表征之一,因此有利于研究葉綠素濃度和海表溫度與凈初級(jí)生產(chǎn)力的相關(guān)性。

    區(qū)域A和區(qū)域B處于不同的經(jīng)緯度帶,分別位于親潮區(qū)域和黑潮區(qū)域,且營(yíng)養(yǎng)鹽濃度差距較大,因此可很好地研究?jī)舫跫?jí)生產(chǎn)力的空間分布和季節(jié)變化特點(diǎn)。此外,2個(gè)區(qū)域均位于寬闊大洋,可去除陸源因素的影響,有利于研究?jī)舫跫?jí)生產(chǎn)力與葉綠素濃度、海表溫度和光合有效輻射3個(gè)環(huán)境因素的相關(guān)性。

    2 結(jié)果與分析

    2.1 各區(qū)域環(huán)境參數(shù)的變化特征

    本研究根據(jù)2007-2018年MODIS衛(wèi)星數(shù)據(jù)產(chǎn)品得到區(qū)域A和區(qū)域B的環(huán)境參數(shù),經(jīng)進(jìn)一步處理,分析各區(qū)域環(huán)境參數(shù)的年平均分布、月平均變化和年際月平均變化。

    2.1.1 年平均分布

    2007—2018年區(qū)域A和區(qū)域B各參數(shù)的年平均分布如圖2和圖3所示。

    由圖2和圖3可以看出:①區(qū)域A 的凈初級(jí)生產(chǎn)力大部分表現(xiàn)為低緯度區(qū)域較高、高緯度區(qū)域較低的特征,西北小部分區(qū)域受陸地因素影響存在一定的異常現(xiàn)象;葉綠素濃度僅在西北部分靠近陸地區(qū)域較高,且低緯度區(qū)域較低、高緯度區(qū)域較高;海表溫度和光合有效輻射的時(shí)空分布表現(xiàn)出很好的緯度分化,即高緯度區(qū)域較低、低緯度區(qū)域較高。②區(qū)域B的凈初級(jí)生產(chǎn)力也表現(xiàn)出一定的緯度分化,即高緯度區(qū)域較高、低緯度區(qū)域較低;葉綠素濃度與凈初級(jí)生產(chǎn)力的分布情況大體相同,且較區(qū)域A偏低,大部分區(qū)域低于0.2 mg/m3;海表溫度和光合有效輻射的分布情況與凈初級(jí)生產(chǎn)力剛好相反,均表現(xiàn)為高緯度區(qū)域較低、低緯度區(qū)域較高。

    2.1.2 月平均變化

    2007—2018年區(qū)域A和區(qū)域B各參數(shù)的月平均變化如圖4所示。

    由圖4可以看出:①區(qū)域A的凈初級(jí)生產(chǎn)力在7—9月較高,呈單峰現(xiàn)象,表現(xiàn)為夏季最高和冬季最低,變化趨勢(shì)與李一凡等[12]的研究結(jié)果具有一致性;葉綠素濃度有2個(gè)峰值,分別出現(xiàn)在5月和10月;海表溫度在7—9月較高,也呈單峰現(xiàn)象,與凈初級(jí)生產(chǎn)力的季節(jié)變化基本相同;光合有效輻射在6月呈單峰現(xiàn)象。②區(qū)域B的凈初級(jí)生產(chǎn)力在1—3月較高,在9月最低,且自10月起有上升趨勢(shì),變化特征與陳小燕[13]的研究結(jié)果相同;葉綠素濃度與凈初級(jí)生產(chǎn)力的季節(jié)變化基本相同,均為冬季最高和夏季最低;海表溫度在7—9月最高,與凈初級(jí)生產(chǎn)力的季節(jié)變化相反;光合有效輻射在7月最高。

    2.1.3 年際月平均變化

    2007—2018年區(qū)域A和區(qū)域B各參數(shù)的年際月平均變化如圖5所示。

    由圖5可以看出,2007—2018年區(qū)域A和區(qū)域B的年際月平均變化趨勢(shì)基本相同。①區(qū)域A的凈初級(jí)生產(chǎn)力略有上升,葉綠素濃度和海表溫度也呈較小的上升趨勢(shì),光合有效輻射的變化較小;②區(qū)域B除海表溫度呈上升趨勢(shì)外,其余參數(shù)的變化較平緩。

    總體來(lái)說(shuō),冬季區(qū)域B的凈初級(jí)生產(chǎn)力較區(qū)域A高,而夏、秋季偏低;但區(qū)域B的葉綠素濃度常年處于較低水平,即使冬季較高仍小于區(qū)域A的最低值;區(qū)域B的海表溫度和光合有效輻射均高于區(qū)域A??梢?jiàn),不同區(qū)域環(huán)境參數(shù)的季節(jié)變化具有一定的差異性,雖然葉綠素濃度是反演估算凈初級(jí)生產(chǎn)力的重要參數(shù),但凈初級(jí)生產(chǎn)力與葉綠素濃度的變化并不完全相同,即影響凈初級(jí)生產(chǎn)力的因素有很多?;诖耍狙芯繉?duì)不同區(qū)域凈初級(jí)生產(chǎn)力的主要影響因素進(jìn)行分析,探討其差異性以及產(chǎn)生差異的原因。

    2.2 各區(qū)域環(huán)境參數(shù)的線性相關(guān)性

    2.2.1 區(qū)域A

    2007—2018年區(qū)域A各參數(shù)逐月平均的相關(guān)性分析如圖6所示。

    由圖6可以看出:區(qū)域A的凈初級(jí)生產(chǎn)力與海表溫度的線性關(guān)系最好,即相關(guān)性最高,相關(guān)系數(shù)約為0.74,具體來(lái)說(shuō)即凈初級(jí)生產(chǎn)力隨著海表溫度的升高而升高;凈初級(jí)生產(chǎn)力與葉綠素濃度和光合有效輻射的相關(guān)性較低,相關(guān)系數(shù)分別約為0.45和0.50;葉綠素濃度與海表溫度的相關(guān)性較低,相關(guān)系數(shù)約為0.19,即葉綠素濃度受海表溫度的影響較小。

    由于區(qū)域A在葉綠素濃度較高的季節(jié)其凈初級(jí)生產(chǎn)力并未在較高水平,葉綠素濃度不是其凈初級(jí)生產(chǎn)力的主要影響因素。區(qū)域A具有高營(yíng)養(yǎng)鹽的特點(diǎn),在這種情況下主要考慮光照條件和海表溫度的影響[14]。在非極地區(qū)域,光照條件對(duì)于凈初級(jí)生產(chǎn)力的影響主要體現(xiàn)在水體的垂直分布上,即光照隨水深的增加而減弱[14]。區(qū)域A位于中高緯度區(qū)域,海表溫度常年處于20℃以下,平均值僅為8.33℃;同時(shí),浮游植物須在酶的作用下進(jìn)行光合作用,而酶的活性直接受溫度的影響,大部分浮游植物進(jìn)行光合作用的最適宜溫度約為20℃[2]。因此,該區(qū)域較低的海表溫度直接限制浮游植物的光合作用速率,凈初級(jí)生產(chǎn)力受海表溫度的影響較大。

    2.2.2 區(qū)域B

    2007—2018年區(qū)域B各參數(shù)逐月平均的相關(guān)性分析如圖7所示。

    由圖7可以看出:區(qū)域B的凈初級(jí)生產(chǎn)力與葉綠素濃度的相關(guān)性較高且呈正相關(guān)關(guān)系,相關(guān)系數(shù)約為0.87;凈初級(jí)生產(chǎn)力與光合有效輻射的相關(guān)性很低,相關(guān)系數(shù)約為0.10;凈初級(jí)生產(chǎn)力與海表溫度的相關(guān)性很高,且呈負(fù)相關(guān)關(guān)系,相關(guān)系數(shù)約為0.92;葉綠素濃度與海表溫度的相關(guān)性較高,且呈負(fù)相關(guān)關(guān)系,相關(guān)系數(shù)約為0.89。

    區(qū)域B屬于低營(yíng)養(yǎng)鹽區(qū)域,浮游植物的初級(jí)生產(chǎn)力受營(yíng)養(yǎng)鹽控制,即使在光照和溫度都非常適宜的情況下,浮游植物仍然無(wú)法暴發(fā)。營(yíng)養(yǎng)鹽對(duì)于浮游植物凈初級(jí)生產(chǎn)力的影響主要體現(xiàn)在生物量上[14],由于所有浮游植物都含有葉綠素a,葉綠素a濃度是浮游植物生物量的重要表征[15]。因此,區(qū)域B的葉綠素濃度較低,且其凈初級(jí)生產(chǎn)力受葉綠素濃度的影響較大。此外,海表溫度可體現(xiàn)海水的層化程度和上升流的強(qiáng)度,而海水的層化程度對(duì)于表層營(yíng)養(yǎng)鹽濃度有一定的影響。區(qū)域B冬季的海表溫度較低,混合層較其他季節(jié)最淺,因此其冬季的表層營(yíng)養(yǎng)鹽濃度最高,葉綠素濃度最高,凈初級(jí)生產(chǎn)力最高。這與陳小燕[13]的研究結(jié)論一致。

    綜上所述:①區(qū)域A的凈初級(jí)生產(chǎn)力與海表溫度的相關(guān)性最高,較低的溫度嚴(yán)重限制酶的活性,因此雖然其葉綠素濃度較區(qū)域B一直處于較高水平,但較低的海表溫度限制其光合作用速率,其凈初級(jí)生產(chǎn)力在冬季最低;②區(qū)域B的凈初級(jí)生產(chǎn)力與葉綠素濃度和海表溫度的相關(guān)性較高,其中與海表溫度相關(guān)性較高的原因主要是葉綠素濃度受海表溫度的影響,因此其凈初級(jí)生產(chǎn)力在冬季最高,且高于區(qū)域A;③與區(qū)域B相比,區(qū)域A各環(huán)境參數(shù)與凈初級(jí)生產(chǎn)力的相關(guān)性均較低;④由于光照對(duì)于海域凈初級(jí)生產(chǎn)力的影響主要體現(xiàn)在極地區(qū)域,2個(gè)區(qū)域的凈初級(jí)生產(chǎn)力與光合有效輻射的相關(guān)性都較低。

    2.3 異常氣候條件下相關(guān)性的驗(yàn)證

    根據(jù)前文分析,2個(gè)區(qū)域的凈初級(jí)生產(chǎn)力與海表溫度都有較高的相關(guān)性,其中區(qū)域B的葉綠素濃度與海表溫度的相關(guān)性也很高。由于研究期間的2015年為強(qiáng)“厄爾尼諾”年份,本研究選取該年份對(duì)2個(gè)區(qū)域的凈初級(jí)生產(chǎn)力、葉綠素濃度和海表溫度與2007—2018年整體進(jìn)行逐月平均對(duì)比,從而驗(yàn)證各環(huán)境參數(shù)相關(guān)性的結(jié)論。

    2.3.1 區(qū)域A

    2015年和2007—2018年區(qū)域A各參數(shù)的逐月平均對(duì)比如圖8所示。

    由圖8可以看出,與2007—2018年相比,2015年:①海表溫度在2-5月表現(xiàn)為正異常,6—10月表現(xiàn)為負(fù)異常,此后負(fù)異常逐漸減弱;②凈初級(jí)生產(chǎn)力在3—5月表現(xiàn)為正異常,6-10月基本表現(xiàn)為負(fù)異常(除7月外),此后負(fù)異常逐漸減弱;③葉綠素濃度有一定的變化,但與凈初級(jí)生產(chǎn)力的變化不存在一致性。

    區(qū)域A處于中高緯度區(qū)域,海洋的熱傳遞占比約為70%,因此該區(qū)域的熱傳遞主要集中于黑潮延伸體和親潮區(qū)域,屬于暖寒交界處。黑潮控制著太平洋中緯度區(qū)域的經(jīng)向輸送,而黑潮異常與“厄爾尼諾”現(xiàn)象具有很高的正相關(guān)性,因此當(dāng)西北太平洋的海表溫度變化時(shí),黑潮區(qū)域的海表溫度隨之變化[16]。

    除海表溫度等基本環(huán)境因素,該區(qū)域的凈初級(jí)生產(chǎn)力還受近年來(lái)氣溶膠沉降等多種因素的影響,例如:2015年7月凈初級(jí)生產(chǎn)力與歷年平均相比較高,可能是當(dāng)年3月底至4月初我國(guó)華北地區(qū)有大范圍沙塵天氣[17],沙塵傳輸?shù)轿鞅碧窖笠话阌?~2個(gè)月的時(shí)滯,導(dǎo)致5—6月西北太平洋的氣溶膠傳輸量變大[12]。

    2.3.2 區(qū)域B

    2015年和2007—2018年區(qū)域B各參數(shù)的逐月平均對(duì)比如圖9所示。

    由圖9可以看出,與2007—2018年相比,2015年:①海表溫度在1-3月表現(xiàn)為負(fù)異常,4—6月表現(xiàn)為正異常,7—11月基本表現(xiàn)為負(fù)異常,12月表現(xiàn)為正異常;②葉綠素濃度與海表溫度的情況相反,即1—3月表現(xiàn)為正異常,4—5月表現(xiàn)為負(fù)異常,6—11月表現(xiàn)為正異常;③凈初級(jí)生產(chǎn)力與葉綠素濃度的情況基本相同。

    區(qū)域B位于西北太平洋副熱帶區(qū)域和黑潮區(qū)域,受“厄爾尼諾”現(xiàn)象的影響較大,即海表溫度波動(dòng)較大,且常年處于較高水平。2015年西北太平洋的海表溫度因“厄爾尼諾”現(xiàn)象有一定的波動(dòng),夏、秋季已達(dá)到超強(qiáng)“厄爾尼諾”標(biāo)準(zhǔn),因此大部分海域的海表溫度在夏、秋季表現(xiàn)為負(fù)異常。區(qū)域A和區(qū)域B海表溫度的波動(dòng)情況與雷星宇[18]關(guān)于2015年“厄爾尼諾”現(xiàn)象對(duì)西太平洋海表溫度影響的研究結(jié)果相同,但區(qū)域A受臺(tái)風(fēng)影響較小,其海表溫度的波動(dòng)主要是由“厄爾尼諾”現(xiàn)象、黑潮作用和緯度變化共同導(dǎo)致的正常波動(dòng)。

    2015年“厄爾尼諾”現(xiàn)象對(duì)于西北太平洋的影響主要體現(xiàn)在臺(tái)風(fēng)季提前,因此區(qū)域B自2015年年初即有臺(tái)風(fēng)頻繁經(jīng)過(guò)[19],導(dǎo)致其凈初級(jí)生產(chǎn)力即固碳量在短時(shí)間內(nèi)出現(xiàn)正異常。小尺度的氣候變化對(duì)于凈初級(jí)生產(chǎn)力的影響主要體現(xiàn)在低營(yíng)養(yǎng)鹽海域,且臺(tái)風(fēng)對(duì)于西北太平洋副熱帶區(qū)域的碳循環(huán)影響較大[20]。區(qū)域B位于西北太平洋的寬闊海域,營(yíng)養(yǎng)躍層較深[21],臺(tái)風(fēng)過(guò)境時(shí)營(yíng)養(yǎng)鹽濃度僅存在微小的波動(dòng)[13]。總體來(lái)說(shuō),區(qū)域B在臺(tái)風(fēng)過(guò)境時(shí)凈初級(jí)生產(chǎn)力明顯升高。

    綜上所述,異常氣候條件下區(qū)域A和區(qū)域B的環(huán)境參數(shù)與凈初級(jí)生產(chǎn)力的相關(guān)性與前述結(jié)論具有一致性。

    3 結(jié)語(yǔ)

    葉綠素濃度可表征浮游植物的生物量,是估算海洋初級(jí)生產(chǎn)力的重要參數(shù)之一[22]。海表溫度和光照條件也是影響浮游植物光合作用的重要因素[23-24],對(duì)于海洋初級(jí)生產(chǎn)力都有一定的影響。而對(duì)于不同海域來(lái)說(shuō),海洋環(huán)境的差異使初級(jí)生產(chǎn)力的主要影響因素有所不同。

    本研究選取海表溫度、葉綠素濃度和光合有效輻射3個(gè)環(huán)境參數(shù),分析其對(duì)不同海域凈初級(jí)生產(chǎn)力的影響。①對(duì)于非極地區(qū)域來(lái)說(shuō),凈初級(jí)生產(chǎn)力受光合有效輻射的影響較小[14]。②區(qū)域A在營(yíng)養(yǎng)鹽濃度較高的情況下主要考慮海表溫度的影響,即酶的活性受海表溫度影響較大。③區(qū)域B在營(yíng)養(yǎng)鹽濃度較低的情況下主要考慮浮游植物生物量的影響,即受葉綠素濃度影響較大[14];由于該區(qū)域位于黑潮區(qū)域,溫度對(duì)于海水混合層深度具有一定的影響,從而影響海水表層營(yíng)養(yǎng)鹽濃度,在低營(yíng)養(yǎng)鹽情況下其影響性很大,最終影響浮游植物生物量,導(dǎo)致葉綠素濃度偏低[13]。④2個(gè)區(qū)域的環(huán)境參數(shù)都有一定的緯度分化特征。⑤對(duì)于異常氣候條件(“厄爾尼諾”現(xiàn)象)下的2015年來(lái)說(shuō),2個(gè)區(qū)域凈初級(jí)生產(chǎn)力的變化特征與本研究的相關(guān)性結(jié)論基本一致。

    本研究選取的2個(gè)區(qū)域具有一定的特殊性和代表性,其中區(qū)域A位于中高緯度的親潮和黑潮寒暖流交界處,區(qū)域B受副熱帶環(huán)流和黑潮共同作用且臺(tái)風(fēng)較活躍。由于2個(gè)區(qū)域的營(yíng)養(yǎng)鹽水平不同,凈初級(jí)生產(chǎn)力的主要影響因素具有差異性,即葉綠素濃度不能完全表征凈初級(jí)生產(chǎn)力。面對(duì)當(dāng)今全球氣候變暖和異常天氣增多的形勢(shì),我國(guó)應(yīng)深入研究西北太平洋初級(jí)生產(chǎn)力的變化,從而預(yù)防赤潮災(zāi)害頻發(fā)和海洋環(huán)境污染等情況,保障海洋漁業(yè)等產(chǎn)業(yè)的發(fā)展。凈初級(jí)生產(chǎn)力的變化受多種因素的影響,很難將每種因素逐一歸結(jié),本研究?jī)H分析少數(shù)基本海洋環(huán)境參數(shù)與凈初級(jí)生產(chǎn)力的相關(guān)性,未來(lái)將收集更多的海洋環(huán)境資料,進(jìn)一步研究?jī)舫跫?jí)生產(chǎn)力的影響因素。

    參考文獻(xiàn)

    [1] RYTHER J H,YENTSCH C S.The estimation of phytoplankton production in the ocean from chlorophyll and Light Data1[J].Limnology and Oceanography,1957,2:281-286.

    [2] BEHRENFELD M J,F(xiàn)ALKOWSKI P G.Photosynthetic rates derived from satellite-based chlorophyll concentration[J].Limnology and Oceanography,1997,42:1-20.

    [3] 陳興群,張明,陳其煥.熱帶太平洋西部及赤道暖水區(qū)的初級(jí)生產(chǎn)力[J].海洋學(xué)報(bào)(中文版),2002(1):86-95.

    [4] 潘友聯(lián).葉綠素與海洋初級(jí)生產(chǎn)力[J].海洋科學(xué),1987(1):63-65.

    [5] 林智濤.南海北部與西部初級(jí)生產(chǎn)力影響因子探究[D].湛江:廣東海洋大學(xué),2017.

    [6] 叢丕福.海洋葉綠素遙感反演及海洋初級(jí)生產(chǎn)力估算研究[D].北京:中國(guó)科學(xué)院研究生院(遙感應(yīng)用研究所),2006.

    [7] CARR M E,F(xiàn)RIEDRICHS M A M,SCHMELTZ M,et al.A comparison of global estimates of marine primary production from ocean color[J].Deep-Sea Research,Part Ⅱ:Topical Studies in Oceanography,2006,53(5/7):741-770.

    [8] GONG G,WEN Y,WANG B,et al.Seasonal variation of chlorophyll a concentration,primary production and environmental conditions in the subtropical East China Sea[J].Deep-Sea Research,Part Ⅱ:Topical Studies in Oceanography,2003,50(6/7):1219-1236.

    [9] LONGHURST A R.Ecological geography of the sea.Second edition[M].Sandiego:Academic Press,2007.

    [10] 寧修仁.西太平洋初級(jí)生產(chǎn)力和對(duì)國(guó)內(nèi)研究工作的幾點(diǎn)建議[J].東海海洋,1984(3):78-83,72.

    [11] 楊海燕,毛新燕,郭新宇.基于WOD數(shù)據(jù)集的西北太平洋混合層內(nèi)營(yíng)養(yǎng)鹽濃度初步研究[J].中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版),2018,48(8):1-9.

    [12] 李一凡,陳文忠.基于MODIS和CALIOP衛(wèi)星遙感數(shù)據(jù)的氣溶膠光學(xué)厚度與海洋初級(jí)生產(chǎn)力相關(guān)性[J].中國(guó)環(huán)境科學(xué),2017,37(1):76-86.

    [13] 陳小燕.基于遙感的長(zhǎng)時(shí)間序列浮游植物的多尺度變化研究[D].杭州:浙江大學(xué),2013.

    [14] 王玲.富營(yíng)養(yǎng)化條件下太湖梅梁灣藻類初級(jí)生產(chǎn)力的光、溫影響研究[D].南京:南京信息工程大學(xué),2016.

    [15] MORALES J,STUART V,PLATT T,et al.Handbook of satellite remote sensing image interpretation:applications for marine living resources conservation and management[R].EU PRESPO and IOCCG,2011.

    [16] 王閃閃.黑潮、厄爾尼諾:南方濤動(dòng)和太平洋年代際濤動(dòng)的相互聯(lián)系及對(duì)氣候影響的研究[D].蘭州:蘭州大學(xué),2015.

    [17] 李珊珊,潘濤,閆靜,等.2015年春季北京市一次沙塵天氣過(guò)程分析[J].環(huán)境科學(xué)與技術(shù),2016,39(4):137-143,160.

    [18] 雷星宇.2015-2016年超強(qiáng)厄爾尼諾事件對(duì)西太副高的影響研究[A].中國(guó)氣象學(xué)會(huì).第35屆中國(guó)氣象學(xué)會(huì)年會(huì) S7 東亞氣候、極端氣候事件變異機(jī)理及氣候預(yù)測(cè)[C].北京:中國(guó)氣象學(xué)會(huì),2018:138-155.

    [19] 錢伊恬,徐邦琪.2015年西北太平洋臺(tái)風(fēng)季提早展開(kāi):2015/2016超級(jí)厄爾尼諾的影響[J].大氣科學(xué)學(xué)報(bào),2016,39(6):788-800.

    [20] LIU F,CHEN C,ZHAN H.Decadal variability of chlorophyll a in the South China Sea:a possible mechanism[J].Chinese Journal of Oceanology and Limnology,2012,30(6):1054-1062.

    [21] GARCIA H E,LOCARNINI R A,BOYER T P,et al.World Ocean Atlas 2009[Z].2010.

    [22] 劉廣鵬,高磊,朱禮鑫,等.長(zhǎng)江口冬季葉綠素活性熒光及遙感分析[J].海洋環(huán)境科學(xué),2015,34(2):245-250.

    [23] FURUYA K,HASEGAWA O,YOSHIKAWA T,et al.Photosynthesis-irradiance relationship of phytoplankton and primary production in the vicinity of Kuroshio warm core ring in spring[J].Journal of Oceanography,1998,54(5):545-552.

    [24] HIROFUMI A,KOZO T.A 9-year time-series of planktonic foraminifer fluxes and environmental change in the Bering sea and the central subarctic Pacific Ocean,1990—1999[J].Progress in Oceanography,2007,72(4):343-363.

    91久久精品电影网| 成人漫画全彩无遮挡| 一区二区三区乱码不卡18| 亚洲精品中文字幕在线视频 | 亚洲怡红院男人天堂| 国产综合精华液| 老司机影院成人| av黄色大香蕉| 国产精品.久久久| 久久99精品国语久久久| 免费大片18禁| 国产一级毛片在线| 大香蕉久久网| 蜜臀久久99精品久久宅男| 精品久久国产蜜桃| 免费黄色在线免费观看| 91精品一卡2卡3卡4卡| 国产免费福利视频在线观看| 午夜激情久久久久久久| 日本欧美国产在线视频| 国产淫语在线视频| 亚洲精品乱码久久久久久按摩| 中文精品一卡2卡3卡4更新| 99九九在线精品视频 | 国产男人的电影天堂91| 乱系列少妇在线播放| 欧美一级a爱片免费观看看| 久久久久久久久久久免费av| 少妇人妻久久综合中文| 国产成人a∨麻豆精品| 嘟嘟电影网在线观看| 中国美白少妇内射xxxbb| 男人舔奶头视频| 夫妻午夜视频| 妹子高潮喷水视频| 五月开心婷婷网| 观看av在线不卡| 欧美日韩视频高清一区二区三区二| 久久久久国产精品人妻一区二区| 男人添女人高潮全过程视频| 男女边摸边吃奶| 免费久久久久久久精品成人欧美视频 | 午夜免费观看性视频| 三级国产精品欧美在线观看| 黑人猛操日本美女一级片| 亚洲欧美中文字幕日韩二区| 亚洲欧美成人精品一区二区| 日本色播在线视频| 最近中文字幕高清免费大全6| 天堂俺去俺来也www色官网| 另类亚洲欧美激情| 亚洲美女黄色视频免费看| 青春草亚洲视频在线观看| 一本一本综合久久| 精品一品国产午夜福利视频| 国产亚洲5aaaaa淫片| 成人影院久久| 欧美日韩亚洲高清精品| 国产精品一区二区三区四区免费观看| 少妇人妻久久综合中文| 国产片特级美女逼逼视频| 一级爰片在线观看| 一级,二级,三级黄色视频| 夫妻午夜视频| 香蕉精品网在线| 免费看光身美女| 亚洲天堂av无毛| 成人18禁高潮啪啪吃奶动态图 | 高清不卡的av网站| 乱人伦中国视频| 国产精品一区二区三区四区免费观看| 草草在线视频免费看| 久久久久久久精品精品| 丝袜在线中文字幕| 中文资源天堂在线| 亚洲天堂av无毛| 天堂8中文在线网| 久久人人爽人人片av| 一级毛片aaaaaa免费看小| 成年美女黄网站色视频大全免费 | 亚洲精品国产成人久久av| 亚洲欧美中文字幕日韩二区| 精品视频人人做人人爽| 男女边摸边吃奶| 欧美国产精品一级二级三级 | 纵有疾风起免费观看全集完整版| 午夜久久久在线观看| 在现免费观看毛片| 日韩欧美一区视频在线观看 | 五月天丁香电影| 老女人水多毛片| av福利片在线观看| 爱豆传媒免费全集在线观看| 麻豆成人午夜福利视频| a级片在线免费高清观看视频| 美女大奶头黄色视频| 国产成人精品婷婷| 国产视频内射| 亚洲国产日韩一区二区| 少妇人妻 视频| 日日摸夜夜添夜夜爱| 香蕉精品网在线| 尾随美女入室| 免费人妻精品一区二区三区视频| 亚洲精品久久午夜乱码| 伊人久久国产一区二区| 色网站视频免费| 一本色道久久久久久精品综合| 欧美 亚洲 国产 日韩一| 黄色欧美视频在线观看| 一级,二级,三级黄色视频| 偷拍熟女少妇极品色| 嘟嘟电影网在线观看| 黄色日韩在线| 一级毛片电影观看| 亚洲国产成人一精品久久久| 天堂中文最新版在线下载| 国国产精品蜜臀av免费| 中文字幕亚洲精品专区| 在线天堂最新版资源| 日韩av免费高清视频| 精品久久久精品久久久| 国产精品国产av在线观看| 国产永久视频网站| 精品久久国产蜜桃| 91精品国产九色| h日本视频在线播放| 亚洲av.av天堂| 日韩亚洲欧美综合| 中文资源天堂在线| 97超视频在线观看视频| 99热这里只有是精品在线观看| 免费看光身美女| 菩萨蛮人人尽说江南好唐韦庄| 一本大道久久a久久精品| 一级毛片aaaaaa免费看小| 在线观看免费日韩欧美大片 | 中文字幕免费在线视频6| 人人妻人人澡人人爽人人夜夜| 免费不卡的大黄色大毛片视频在线观看| 日韩亚洲欧美综合| 岛国毛片在线播放| 人人妻人人澡人人看| 一级毛片我不卡| 男女无遮挡免费网站观看| 亚洲精品久久久久久婷婷小说| 国产黄片视频在线免费观看| 日韩视频在线欧美| 69精品国产乱码久久久| 成年人午夜在线观看视频| 人人妻人人看人人澡| 国产真实伦视频高清在线观看| 黄色毛片三级朝国网站 | 9色porny在线观看| 精品亚洲乱码少妇综合久久| 三级国产精品欧美在线观看| 久久6这里有精品| av在线app专区| 欧美精品一区二区大全| 日日爽夜夜爽网站| 国产精品国产av在线观看| 看非洲黑人一级黄片| 国产在视频线精品| 不卡视频在线观看欧美| 九色成人免费人妻av| 日本黄色片子视频| 午夜福利在线观看免费完整高清在| 多毛熟女@视频| 国产一级毛片在线| 在线观看av片永久免费下载| 国产成人免费无遮挡视频| 又黄又爽又刺激的免费视频.| 天天躁夜夜躁狠狠久久av| 综合色丁香网| 十八禁高潮呻吟视频 | 18禁动态无遮挡网站| 亚洲第一区二区三区不卡| 狂野欧美激情性xxxx在线观看| 特大巨黑吊av在线直播| 91午夜精品亚洲一区二区三区| 久久99热这里只频精品6学生| 中文在线观看免费www的网站| 22中文网久久字幕| 18禁在线无遮挡免费观看视频| 夜夜爽夜夜爽视频| 3wmmmm亚洲av在线观看| 国产黄色视频一区二区在线观看| 人人妻人人澡人人爽人人夜夜| 国产男女内射视频| 内地一区二区视频在线| 日本与韩国留学比较| 一区二区三区乱码不卡18| 18+在线观看网站| 国产极品粉嫩免费观看在线 | 成人二区视频| 精品人妻偷拍中文字幕| 久久午夜福利片| 国产欧美日韩精品一区二区| 亚洲国产av新网站| 久久精品国产亚洲网站| 久久人妻熟女aⅴ| 视频中文字幕在线观看| 久久狼人影院| 精品一区二区免费观看| 性色av一级| 人体艺术视频欧美日本| 久久久久久久精品精品| 午夜免费男女啪啪视频观看| 18禁裸乳无遮挡动漫免费视频| 免费在线观看成人毛片| 日本-黄色视频高清免费观看| 18禁在线播放成人免费| videos熟女内射| 蜜臀久久99精品久久宅男| 亚洲va在线va天堂va国产| a级片在线免费高清观看视频| 日本vs欧美在线观看视频 | 亚洲伊人久久精品综合| 建设人人有责人人尽责人人享有的| av福利片在线观看| 国产高清有码在线观看视频| 80岁老熟妇乱子伦牲交| 各种免费的搞黄视频| 婷婷色av中文字幕| 日韩一本色道免费dvd| av卡一久久| 成人午夜精彩视频在线观看| 久久久久视频综合| 久久久久久久久久久丰满| 这个男人来自地球电影免费观看 | tube8黄色片| 久久国产乱子免费精品| 韩国av在线不卡| 最近中文字幕2019免费版| av天堂久久9| 国产一区有黄有色的免费视频| 一区在线观看完整版| 女性被躁到高潮视频| 精品酒店卫生间| 国产精品欧美亚洲77777| 亚洲精品乱码久久久v下载方式| 男女啪啪激烈高潮av片| 亚洲四区av| 国产成人freesex在线| 美女国产视频在线观看| 色视频www国产| 美女脱内裤让男人舔精品视频| 观看美女的网站| 日本欧美视频一区| 亚洲伊人久久精品综合| 久久久久久久国产电影| 欧美日韩视频精品一区| av女优亚洲男人天堂| 久久久久久久久久久免费av| 久久青草综合色| 亚洲丝袜综合中文字幕| 爱豆传媒免费全集在线观看| 国产乱人偷精品视频| 街头女战士在线观看网站| 日韩亚洲欧美综合| 肉色欧美久久久久久久蜜桃| 热re99久久国产66热| 蜜桃在线观看..| 亚洲国产欧美在线一区| av在线app专区| 国产熟女午夜一区二区三区 | 国产精品国产av在线观看| 成人国产av品久久久| 亚洲伊人久久精品综合| 在线亚洲精品国产二区图片欧美 | 国产高清国产精品国产三级| 免费人成在线观看视频色| 亚洲av成人精品一二三区| 丰满饥渴人妻一区二区三| 寂寞人妻少妇视频99o| 美女福利国产在线| 日韩,欧美,国产一区二区三区| 一级毛片电影观看| a级一级毛片免费在线观看| 18禁裸乳无遮挡动漫免费视频| 中文字幕制服av| 国产女主播在线喷水免费视频网站| 国产免费视频播放在线视频| 亚洲国产日韩一区二区| 嫩草影院新地址| 另类亚洲欧美激情| 啦啦啦在线观看免费高清www| 免费黄网站久久成人精品| 亚洲欧美精品专区久久| 久久99热6这里只有精品| 高清黄色对白视频在线免费看 | 亚洲av欧美aⅴ国产| 99九九在线精品视频 | 在现免费观看毛片| av天堂中文字幕网| 国产成人a∨麻豆精品| 国产精品99久久久久久久久| 自拍欧美九色日韩亚洲蝌蚪91 | 2018国产大陆天天弄谢| 亚洲国产毛片av蜜桃av| 亚洲综合色惰| 自线自在国产av| 欧美日韩在线观看h| 欧美精品亚洲一区二区| 亚洲av欧美aⅴ国产| 久久国内精品自在自线图片| 国产熟女欧美一区二区| 夜夜骑夜夜射夜夜干| 少妇高潮的动态图| 国精品久久久久久国模美| 国产成人一区二区在线| 久久人人爽av亚洲精品天堂| 日韩成人伦理影院| 在线观看国产h片| 久久精品熟女亚洲av麻豆精品| 桃花免费在线播放| 观看免费一级毛片| 我的女老师完整版在线观看| 十八禁高潮呻吟视频 | 成年女人在线观看亚洲视频| 大片电影免费在线观看免费| 欧美一级a爱片免费观看看| 我的女老师完整版在线观看| 少妇精品久久久久久久| 久久久久视频综合| 国产色婷婷99| 中国三级夫妇交换| 亚洲av成人精品一区久久| 少妇人妻久久综合中文| 涩涩av久久男人的天堂| 免费不卡的大黄色大毛片视频在线观看| av线在线观看网站| 亚洲国产精品专区欧美| 久久久久久久久久久丰满| 久久久久精品久久久久真实原创| 成人二区视频| 日韩中字成人| 丝瓜视频免费看黄片| 国产极品粉嫩免费观看在线 | 亚洲第一av免费看| 大片电影免费在线观看免费| 秋霞伦理黄片| 国产国拍精品亚洲av在线观看| 日本91视频免费播放| 日本wwww免费看| 久久久久久久国产电影| 国产精品秋霞免费鲁丝片| 中文字幕亚洲精品专区| 五月伊人婷婷丁香| 亚洲美女搞黄在线观看| 少妇被粗大猛烈的视频| 亚洲美女搞黄在线观看| 日日撸夜夜添| 观看免费一级毛片| 中文精品一卡2卡3卡4更新| 午夜影院在线不卡| 精品午夜福利在线看| 老女人水多毛片| 七月丁香在线播放| 十八禁高潮呻吟视频 | 另类亚洲欧美激情| 18禁在线播放成人免费| 日本黄色日本黄色录像| 日日爽夜夜爽网站| 国产男女超爽视频在线观看| 精品久久久精品久久久| 日本黄大片高清| 亚洲四区av| 三级国产精品欧美在线观看| 乱系列少妇在线播放| 久久热精品热| 久久99热6这里只有精品| 一区二区三区精品91| 韩国高清视频一区二区三区| 观看av在线不卡| 中文字幕制服av| 最新中文字幕久久久久| 国产乱人偷精品视频| 国产成人91sexporn| 国产在线一区二区三区精| 草草在线视频免费看| 国产在线一区二区三区精| 美女国产视频在线观看| 夜夜爽夜夜爽视频| 成人特级av手机在线观看| 国产成人午夜福利电影在线观看| 观看美女的网站| 亚洲国产日韩一区二区| 好男人视频免费观看在线| 亚洲成人av在线免费| 最后的刺客免费高清国语| 欧美日韩一区二区视频在线观看视频在线| 成年人午夜在线观看视频| 男的添女的下面高潮视频| 亚洲欧美一区二区三区国产| 亚洲欧美精品专区久久| 亚洲人成网站在线播| 少妇被粗大猛烈的视频| 新久久久久国产一级毛片| 18禁裸乳无遮挡动漫免费视频| 亚洲第一av免费看| 91久久精品国产一区二区三区| kizo精华| 日本av手机在线免费观看| 黄色配什么色好看| 精品国产国语对白av| 最近中文字幕2019免费版| 免费观看的影片在线观看| 永久网站在线| 久久精品夜色国产| 日本色播在线视频| 精品人妻一区二区三区麻豆| 国产亚洲欧美精品永久| av视频免费观看在线观看| 国产高清国产精品国产三级| 成年人免费黄色播放视频 | 精品少妇黑人巨大在线播放| 精品少妇内射三级| 制服丝袜香蕉在线| 一级黄片播放器| 天天躁夜夜躁狠狠久久av| 99久国产av精品国产电影| 亚洲av二区三区四区| 久久99热这里只频精品6学生| 日韩一区二区视频免费看| 日韩制服骚丝袜av| 久久人人爽人人片av| 亚洲国产精品一区二区三区在线| 美女内射精品一级片tv| 国产av一区二区精品久久| 777米奇影视久久| av黄色大香蕉| 日本午夜av视频| 亚洲av免费高清在线观看| 女人久久www免费人成看片| 黄色怎么调成土黄色| 中文资源天堂在线| 精品一区在线观看国产| 免费人成在线观看视频色| 亚洲av福利一区| 国产深夜福利视频在线观看| 日韩一本色道免费dvd| 最近手机中文字幕大全| 午夜福利,免费看| 美女xxoo啪啪120秒动态图| 久久精品熟女亚洲av麻豆精品| 各种免费的搞黄视频| 乱码一卡2卡4卡精品| 精品国产乱码久久久久久小说| 国产黄色视频一区二区在线观看| 制服丝袜香蕉在线| 成人无遮挡网站| 免费av不卡在线播放| 国精品久久久久久国模美| 人妻一区二区av| 国产国拍精品亚洲av在线观看| 亚洲av不卡在线观看| 自拍偷自拍亚洲精品老妇| 欧美精品亚洲一区二区| 寂寞人妻少妇视频99o| 夜夜骑夜夜射夜夜干| 欧美日韩国产mv在线观看视频| 中文乱码字字幕精品一区二区三区| 国产精品一区二区性色av| av国产久精品久网站免费入址| 欧美日本中文国产一区发布| xxx大片免费视频| 99九九线精品视频在线观看视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧洲国产日韩| 亚洲av综合色区一区| 最新的欧美精品一区二区| 国产免费一区二区三区四区乱码| 99热这里只有是精品50| 在线观看免费视频网站a站| 九色成人免费人妻av| 日韩欧美一区视频在线观看 | 免费黄网站久久成人精品| 国产欧美日韩综合在线一区二区 | 黑丝袜美女国产一区| 最近手机中文字幕大全| 免费观看在线日韩| 久久国内精品自在自线图片| 肉色欧美久久久久久久蜜桃| 最近手机中文字幕大全| 国产精品久久久久久久电影| 精品久久久精品久久久| 国产熟女欧美一区二区| 亚洲va在线va天堂va国产| 全区人妻精品视频| 男男h啪啪无遮挡| 亚洲国产精品999| 一级毛片久久久久久久久女| 亚洲性久久影院| 色视频在线一区二区三区| 内射极品少妇av片p| 亚洲精品久久午夜乱码| 亚洲欧美日韩东京热| 免费观看性生交大片5| 久久久久国产网址| 久久国产亚洲av麻豆专区| 人妻制服诱惑在线中文字幕| 亚洲色图综合在线观看| 久久鲁丝午夜福利片| 综合色丁香网| 人妻制服诱惑在线中文字幕| 久久6这里有精品| 国产黄片美女视频| 精品少妇内射三级| 伦理电影大哥的女人| 丰满少妇做爰视频| av.在线天堂| 老司机亚洲免费影院| 永久网站在线| 久久久久精品性色| 欧美精品人与动牲交sv欧美| 99热全是精品| 观看av在线不卡| 日日啪夜夜爽| 久久久久视频综合| 美女中出高潮动态图| 成人特级av手机在线观看| 久久av网站| 日日爽夜夜爽网站| 大香蕉久久网| 9色porny在线观看| 欧美激情国产日韩精品一区| 亚洲精品乱久久久久久| 亚洲三级黄色毛片| 午夜久久久在线观看| 80岁老熟妇乱子伦牲交| a级毛片免费高清观看在线播放| 欧美高清成人免费视频www| 成人漫画全彩无遮挡| 久久精品国产亚洲av涩爱| 国产av国产精品国产| 久久久久久久久久久免费av| 日韩欧美精品免费久久| 亚洲第一区二区三区不卡| 麻豆成人午夜福利视频| 2022亚洲国产成人精品| 国产一级毛片在线| 亚洲美女视频黄频| 亚洲国产毛片av蜜桃av| 亚洲欧洲日产国产| 人人妻人人看人人澡| 一级毛片久久久久久久久女| 成年人免费黄色播放视频 | 中文字幕人妻熟人妻熟丝袜美| 97超视频在线观看视频| 秋霞伦理黄片| 丰满迷人的少妇在线观看| 久久鲁丝午夜福利片| 大片免费播放器 马上看| 免费观看av网站的网址| 简卡轻食公司| 欧美亚洲 丝袜 人妻 在线| 亚洲精品456在线播放app| 国产欧美日韩综合在线一区二区 | 纵有疾风起免费观看全集完整版| 少妇 在线观看| 91aial.com中文字幕在线观看| 高清毛片免费看| 日本vs欧美在线观看视频 | 国国产精品蜜臀av免费| 亚洲国产精品999| 91精品国产国语对白视频| 色网站视频免费| 五月天丁香电影| 精品人妻偷拍中文字幕| 久久99热6这里只有精品| 自拍偷自拍亚洲精品老妇| 免费看日本二区| 女人精品久久久久毛片| 日韩精品免费视频一区二区三区 | 新久久久久国产一级毛片| 久久久久视频综合| 亚洲精品一二三| 六月丁香七月| 亚洲电影在线观看av| 女性被躁到高潮视频| 久久久午夜欧美精品| 亚洲精品第二区| 免费人成在线观看视频色| 亚州av有码| 国产精品不卡视频一区二区| 国产免费一区二区三区四区乱码| 午夜福利网站1000一区二区三区| 熟女人妻精品中文字幕| 国产成人精品福利久久| 男男h啪啪无遮挡| 伊人亚洲综合成人网| 久久6这里有精品| 高清午夜精品一区二区三区| 午夜激情福利司机影院| av福利片在线| 日日啪夜夜爽| 激情五月婷婷亚洲| 亚洲,欧美,日韩| 午夜福利视频精品| 成年美女黄网站色视频大全免费 | 男人和女人高潮做爰伦理| 欧美精品一区二区免费开放| 观看免费一级毛片| 久久综合国产亚洲精品| 全区人妻精品视频| a 毛片基地| 最黄视频免费看| 黄色欧美视频在线观看| 成人18禁高潮啪啪吃奶动态图 | 香蕉精品网在线| 七月丁香在线播放| 午夜免费观看性视频| 午夜免费鲁丝| 精品国产露脸久久av麻豆|