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

    基于環(huán)境一號(hào)衛(wèi)星的大洋河河口海域營養(yǎng)鹽遙感反演

    2012-12-26 02:34:34趙冬至楊建洪國家海洋環(huán)境監(jiān)測(cè)中心海洋遙感室遼寧大連116023
    中國環(huán)境科學(xué) 2012年1期
    關(guān)鍵詞:營養(yǎng)鹽河口反射率

    王 林,趙冬至,楊建洪 (國家海洋環(huán)境監(jiān)測(cè)中心海洋遙感室,遼寧 大連 116023)

    基于環(huán)境一號(hào)衛(wèi)星的大洋河河口海域營養(yǎng)鹽遙感反演

    王 林,趙冬至*,楊建洪 (國家海洋環(huán)境監(jiān)測(cè)中心海洋遙感室,遼寧 大連 116023)

    綜合利用大洋河河口海域水體調(diào)查的實(shí)測(cè)數(shù)據(jù)及環(huán)境一號(hào) CCD影像數(shù)據(jù),進(jìn)行了該區(qū)域營養(yǎng)鹽反演研究.結(jié)果表明,對(duì)數(shù)函數(shù)擬合NO2--N與Rrs(Band2)效果最佳,而線性函數(shù)擬合NH4+-N與Rrs(Band2)、NO3--N與Rrs(Band1)、TIN與Rrs(Band1)效果最佳,相關(guān)性R2分別可達(dá)0.939、0.935、0.945、0.970,未發(fā)現(xiàn)較佳的PO43--P反演模型;營養(yǎng)鹽在可見光波段并非光敏感性物質(zhì),其與等效遙感反射率存在相關(guān)性的主要原因在于河口海域陸源輸入物質(zhì)典型的擴(kuò)散特征,使得營養(yǎng)鹽與某些水色參數(shù)存在較高的相關(guān)性,NO2--N、NH4+-N、NO3--N 以及TIN與有色溶解性有機(jī)物吸收系數(shù)ag(400)的相關(guān)性R2分別可達(dá)0.836、0.808、0.786、0.854,這決定了在河口海域進(jìn)行營養(yǎng)鹽遙感反演的可行性;從衛(wèi)星產(chǎn)品專題圖可看出,大洋河河口海域NO2--N、NH4+-N、NO3--N以及TIN都存在一個(gè)明顯的高值區(qū),該高值區(qū)營養(yǎng)鹽濃度由河口向外逐漸降低,直至外海最低.

    環(huán)境一號(hào)衛(wèi)星;大洋河河口;營養(yǎng)鹽;遙感反演

    河口海域是陸海相互作用耦合帶和生產(chǎn)力最高的區(qū)域,生源元素來源豐富,各種因素(包括水動(dòng)力、生物地球化學(xué)過程及人類活動(dòng)等)十分復(fù)雜,其生態(tài)環(huán)境系統(tǒng)和生物資源的變化與人類活動(dòng)和經(jīng)濟(jì)發(fā)展關(guān)系密切[1].近岸海域濃度偏高的營養(yǎng)鹽所帶來的富營養(yǎng)化問題是我國沿岸水域最突出的環(huán)境問題之一[2].傳統(tǒng)的水質(zhì)監(jiān)測(cè)方法需要耗費(fèi)較多的人力、物力與財(cái)力,而且采樣點(diǎn)布設(shè)的數(shù)量有限,監(jiān)測(cè)周期較長,難以有效的對(duì)整體水域進(jìn)行大尺度、動(dòng)態(tài)的監(jiān)測(cè).遙感技術(shù)以獨(dú)特的優(yōu)勢(shì)為水質(zhì)監(jiān)測(cè)開辟了新的途徑,它可以實(shí)現(xiàn)對(duì)水體快速、大范圍、低成本、周期性的動(dòng)態(tài)監(jiān)測(cè),具有不可替代的優(yōu)越性[3-4].然而,營養(yǎng)鹽并非光敏感性物質(zhì),在自然水體中吸收、散射光能微弱,對(duì)其進(jìn)行遙感反演具有一定難度,國內(nèi)外相關(guān)研究成果較少.目前,營養(yǎng)鹽遙感反演方法主要分為間接法和直接法.間接法主要是基于特定海域營養(yǎng)鹽與某些光敏感性物質(zhì)(葉綠素、海溫等)的關(guān)系,建立營養(yǎng)鹽遙感反演模型[5].直接法則主要利用偏最小二乘法、微分光譜技術(shù)等獲得營養(yǎng)鹽與光譜數(shù)據(jù)的關(guān)系,實(shí)現(xiàn)營養(yǎng)鹽濃度的遙感反演[6-12].

    MODIS、HY-1等海洋水色衛(wèi)星空間分辨率較低,對(duì)于河口海域水體環(huán)境要素監(jiān)測(cè)存在一定難度. 我國于2008年9月6日成功發(fā)射了第一顆專用于環(huán)境監(jiān)測(cè)的環(huán)境一號(hào)衛(wèi)星,其攜帶的CCD相機(jī)較高的空間分辨率可有效監(jiān)測(cè)面積相對(duì)較小的河口海域,而較短的重返周期則有利于實(shí)現(xiàn)對(duì)水體的大規(guī)模實(shí)時(shí)監(jiān)測(cè).目前為止利用該衛(wèi)星數(shù)據(jù)對(duì)河口海域營養(yǎng)鹽遙感監(jiān)測(cè)的相關(guān)研究尚未見報(bào)道.本研究利用2009年5月6日大洋河河口海域的實(shí)測(cè)數(shù)據(jù),建立該區(qū)域營養(yǎng)鹽遙感反演模型,并將其應(yīng)用于環(huán)境一號(hào) CCD影像上,從而為分析大洋河河口海域營養(yǎng)鹽含量時(shí)空變化特征、富營養(yǎng)化狀況以及其它污染問題研究提供依據(jù).

    1 材料與方法

    1.1 研究區(qū)域與站位布設(shè)

    圖1 大洋河河口海域水體調(diào)查站位Fig.1 Stations of the investigation cruise in the Dayang River Estuary

    研究區(qū)域選擇在我國北方大洋河河口海域.大洋河河口海域位于大洋河流域南部,具有明顯的由陸地向海洋過渡性特點(diǎn),自然條件復(fù)雜,生態(tài)環(huán)境脆弱.而大洋河流域處于遼東半島的東北部,遼寧省營口市以西,丹東市以東,流域范圍主要包括岫巖滿族自治縣大部分區(qū)域,鳳城市西南部和東港市西部區(qū)域,流域面積 6515.25km2.調(diào)查于2009年5月6日進(jìn)行,依據(jù)河口海域水體擴(kuò)散特征,可體現(xiàn)徑流輸入物質(zhì)的濃度梯度變化,共布設(shè)站點(diǎn)13個(gè),如圖1所示.

    1.2 數(shù)據(jù)獲取與分析

    1.2.1 光譜數(shù)據(jù)獲取 現(xiàn)場(chǎng)光譜測(cè)量采用水面以上測(cè)量法[13],實(shí)驗(yàn)儀器為美國 ASD公司生產(chǎn)的便攜式雙通道可見光-近紅外光譜輻射計(jì).儀器觀測(cè)平面與太陽入射平面的夾角為135°,儀器與海面法線方向的夾角為 40°,以避免絕大部分的太陽直射反射,并減少船舶陰影的影響.在每個(gè)站位進(jìn)行2次獨(dú)立的測(cè)量.數(shù)據(jù)處理后得到遙感反射率Rrs,按如下公式計(jì)算:

    式中, Lsw為測(cè)量水體的輻亮度;r為氣-水界面對(duì)天空漫散射光的反射率,一般情況可取 0.028,平靜海面時(shí)取0.026[14];Lsky為天空漫散射光的輻亮度;ρp為標(biāo)準(zhǔn)灰板的反射率;Lp為測(cè)量標(biāo)準(zhǔn)板的輻亮度.

    光譜數(shù)據(jù)采集于09:00~15:00之間完成,到達(dá)1、2號(hào)站點(diǎn)時(shí),由于時(shí)間較早,不滿足光譜測(cè)定所需太陽高度角要求,僅進(jìn)行了水色和化學(xué)參數(shù)測(cè)量,未進(jìn)行現(xiàn)場(chǎng)光譜測(cè)量,共獲取光譜數(shù)據(jù)11組.

    獲取現(xiàn)場(chǎng)水體ASD光譜遙感反射率后,根據(jù)環(huán)境一號(hào)衛(wèi)星CCD相機(jī)的光譜響應(yīng)函數(shù)進(jìn)行實(shí)測(cè)Rrs(λ)的波段等效處理,具體公式如下:

    式中, Rrs(Bandx)為HJ星CCD傳感器Bandx波段的等效遙感反射率;Rrs(λ)為ASD光譜儀現(xiàn)場(chǎng)獲取遙感反射率;Fs(λ)為日地平均距離處大氣層外太陽輻照度;Sx(λ)為 Bandx波段的光譜響應(yīng)函數(shù).

    1.2.2 水色參數(shù)獲取 色素濃度測(cè)量采用熒光法[15],測(cè)量時(shí)使用美國 Turner公司生產(chǎn)的TD-700實(shí)驗(yàn)室熒光儀.數(shù)據(jù)處理后可得葉綠素a濃度(Chl-a)和脫鎂葉綠素濃度(Pheo),將兩濃度之和作為總?cè)~綠素濃度(TChl).熒光儀定標(biāo)采用Sigma公司生產(chǎn)的Chl-a標(biāo)準(zhǔn)樣品.

    總懸浮物濃度(TSM)測(cè)量采用重量法,稱量時(shí)使用有效數(shù)字至少4位的電子天平稱量.

    有色溶解性有機(jī)物(CDOM)吸收系數(shù)測(cè)量采用實(shí)驗(yàn)室分光光度計(jì)法[16],測(cè)量時(shí)使用Hitachi公司生產(chǎn)的 U-3010型紫外-可見分光光度計(jì).CDOM濃度高低用400nm波段處的吸收系數(shù)ag(400)表示.

    共獲取水色參數(shù)數(shù)據(jù)13組.

    1.2.3 營養(yǎng)鹽等參數(shù)獲取 硝酸鹽(NO3--N)、亞硝酸鹽(NO2--N)以及銨鹽(NH4+-N)測(cè)量采用鹽酸萘乙二胺分光光度法,樣品測(cè)量使用7230分光光度計(jì),將NO3--N、NO2--N以及NH4+-N濃度之和作為總?cè)芙庑詿o機(jī)氮(TIN).無機(jī)磷(PO43--P)測(cè)量采用磷鉬蘭法[17],樣品測(cè)量使用7230分光光度計(jì).需要說明的是,營養(yǎng)鹽測(cè)量后發(fā)現(xiàn)13號(hào)站點(diǎn)NO3--N數(shù)值異常,建模及相關(guān)性分析時(shí)未采用.共獲取 NO2--N、NH4+-N以及PO43--P數(shù)據(jù)13組,NO3-N數(shù)據(jù)12組.

    鹽度(Sal)數(shù)據(jù)采用 Hydrolab多參數(shù)水質(zhì)儀現(xiàn)場(chǎng)測(cè)量獲得.到達(dá) 13號(hào)站點(diǎn)時(shí),測(cè)量鹽度所用的多參數(shù)水質(zhì)儀供電蓄電池電量不足,未進(jìn)行鹽度測(cè)量.共獲取鹽度數(shù)據(jù)12組.

    1.2.4 衛(wèi)星數(shù)據(jù)獲取 在HJ-1A衛(wèi)星和HJ-1B衛(wèi)星上均裝載的2臺(tái)CCD相機(jī)設(shè)計(jì)原理完全相同,以星下點(diǎn)對(duì)稱放置,平分視場(chǎng)、并行觀測(cè),聯(lián)合完成對(duì)地刈幅寬度為700km、地面像元分辨率為30m、4個(gè)譜段的推掃成像.本研究未能獲取5月6日外業(yè)調(diào)查當(dāng)日衛(wèi)星數(shù)據(jù),而于5月24日獲取一景HJ-1A CCD2影像.

    CCD影像的預(yù)處理主要包括:輻射校正、幾何校正和大氣校正等.

    輻射校正:利用絕對(duì)定標(biāo)系數(shù)將 CCD圖像DN值轉(zhuǎn)換為輻亮度圖像的公式為:

    式中, A為絕對(duì)定標(biāo)系數(shù)增益;L0為絕對(duì)定標(biāo)系數(shù)偏移量,轉(zhuǎn)換后輻亮度單位為 W/(m2?sr?μm).其中A和L0的具體參數(shù)由《環(huán)境減災(zāi)星座A/B星各載荷在軌絕對(duì)輻射定標(biāo)系數(shù)》[18]提供(中國資源衛(wèi)星應(yīng)用中心,2009).

    幾何校正:選用經(jīng)過空間配準(zhǔn)的 1:5萬地形圖為參考圖像,在ENVI 4.7支持下對(duì)環(huán)境一號(hào)衛(wèi)星CCD影像進(jìn)行幾何校正,圖像重采樣采用最近鄰點(diǎn)法,總誤差控制在0.5個(gè)像元內(nèi).

    大氣校正:采用ENVI 4.7軟件下的FLAASH大氣校正模塊對(duì)環(huán)境一號(hào)衛(wèi)星CCD數(shù)據(jù)進(jìn)行大氣校正,輸入相關(guān)參數(shù),計(jì)算得到大氣校正后的遙感反射率圖像.

    1.3 算法精度評(píng)價(jià)指標(biāo)

    本研究營養(yǎng)鹽反演算法精度評(píng)價(jià)指標(biāo)采用均方根誤差(erms)和相對(duì)誤差(er).

    2 結(jié)果與討論

    2.1 營養(yǎng)鹽遙感反演模型建立

    將經(jīng)HJ-1A CCD2光譜響應(yīng)函數(shù)處理后的4個(gè)波段的等效遙感反射率分別與營養(yǎng)鹽濃度進(jìn)行相關(guān)性分析,結(jié)果如表 1所示.可發(fā)現(xiàn), Band1 (430~520nm)、Band2(520~600nm) 以及 Band3 (630~690nm)的等效遙感反射率均與 NO2--N、NH4+-N、NO3--N以及TIN存在較高的相關(guān)性,除Rrs(Band3)與NO2--N相關(guān)性相對(duì)較低外,其余R2均高于 0.85,且大部分在 0.90以上;而 Band4 (760~900nm)的等效遙感反射率與所有營養(yǎng)鹽參數(shù)的相關(guān)性均非常低,可見該波段在營養(yǎng)鹽反演中并不適用;值得注意的是,未發(fā)現(xiàn)任何波段與PO43--P存在明顯的相關(guān)性.總體看來,HJ-1A CCD2 Band1、Band2以及Band3可用于NO2-N、NH4+-N、NO3--N以及TIN的遙感反演,結(jié)合表1確定將Band2作為反演NO2--N、NH4+-N的敏感波段,Band1作為反演NO3--N、TIN的敏感波段,而PO43--P不存在任何敏感波段.

    表1 河口海域營養(yǎng)鹽參數(shù)與4個(gè)波段的等效遙感反射率間相關(guān)性分析結(jié)果Table 1 Correlation analysis results between nutrients and equivalent remote sensing reflectance

    分別將Rrs(Band2)與NO2--N和NH4+-N、Rrs(Band1)與NO3--N和TIN進(jìn)行線性、對(duì)數(shù)以及e指數(shù)模型回歸分析,算法精度見表2.NO2--N反演時(shí),對(duì)數(shù)模型表現(xiàn)最佳,相關(guān)性R2為0.939,平均相對(duì)誤差為 26.9%,均方根誤差為0.09148μmol/L,線性模型次之,e指數(shù)模型最差; NH4+-N反演時(shí),線性模型與對(duì)數(shù)模型表現(xiàn)近似,線性模型相對(duì)較佳,e指數(shù)模型最差;NO3--N 和TIN反演時(shí),線性模型表現(xiàn)最佳,對(duì)數(shù)模型次之,e指數(shù)模型最差.此外,NO3--N反演模型的平均相對(duì)誤差普遍偏高,分析原因在于NO3--N 2個(gè)極低值點(diǎn)(0.31429,0.25714μmol/L)反演相對(duì)誤差較大造成.因此,本研究選擇對(duì)數(shù)模型反演 NO2--N,而選擇線性模型反演NH4+-N、NO3--N以及TIN.

    表2 NO2--N、NH4+-N、NO3--N以及TIN反演精度統(tǒng)計(jì)Table 2 The statistics of various nutrients retrieval accuracy

    2.2 營養(yǎng)鹽遙感反演機(jī)制分析

    大量研究表明,NO2--N、NH4+-N以及NO3--N在可見光波段吸收、散射光能微弱,但其與遙感反射率存在的相關(guān)性該如何解釋成為本研究亟需解決的關(guān)鍵問題.將河口海域理、化以及光學(xué)參數(shù)等作相關(guān)性分析,主要包括:代表徑流擴(kuò)散狀況的指數(shù)Sal;水色參數(shù)TChl、TSM以及ag(400);營養(yǎng)鹽參數(shù)NO2--N、NH4+-N、NO3--N、TIN以及PO43--P,分析結(jié)果如表 3所示.Sal與其他因子(除PO43--P外)均存在一定的相關(guān)性,其中,與顆粒態(tài)物質(zhì)的相關(guān)性明顯低于溶解態(tài)物質(zhì)(除 PO43--P外),說明河口海域的污染物質(zhì)主要來自陸源徑流輸入,而且徑流輸入的溶解態(tài)物質(zhì)均勻擴(kuò)散到海水中,且就本研究而言,線性擴(kuò)散特征明顯,而顆粒態(tài)物質(zhì)的線性擴(kuò)散特征相對(duì)較弱;此外,NO2--N、NH4+-N、NO3--N以及TIN與水色參數(shù)中溶解態(tài)ag(400)的相關(guān)性R2分別可達(dá)0.836、0.808、0.786、0.854,普遍高于顆粒態(tài)水色參數(shù)TChl、TSM;值得注意的是,本研究中PO43--P與其他參數(shù)均不存在明顯的相關(guān)性,可見其來源和擴(kuò)散特征與其他參數(shù)不同,在今后研究中需引起重視.基于以上分析,可知河口海域陸源徑流輸入物質(zhì)典型的擴(kuò)散特征,及在此基礎(chǔ)上使得營養(yǎng)鹽參數(shù)與某些水色參數(shù)存在較高的相關(guān)性是進(jìn)行河口海域營養(yǎng)鹽遙感反演的重要依據(jù).

    2.3 衛(wèi)星遙感產(chǎn)品的實(shí)現(xiàn)

    由于受天氣、海況以及衛(wèi)星過境時(shí)間等多種因素影響,本研究所能獲取距5月6日外業(yè)調(diào)查時(shí)間最近且圖像質(zhì)量最好的衛(wèi)星數(shù)據(jù)為5月24日HJ-1A CCD2影像.本文營養(yǎng)鹽反演基礎(chǔ)在于河口海域陸源徑流輸入物質(zhì)典型的擴(kuò)散特征,而北方地區(qū)5月份為枯水季節(jié),河口徑流量較少且相對(duì)穩(wěn)定,因此基于5月6日實(shí)測(cè)數(shù)據(jù)所建模型可用于5月24日衛(wèi)星數(shù)據(jù).將NO2--N、NH4+-N、NO3--N以及TIN反演模型應(yīng)用于環(huán)境一號(hào)衛(wèi)星CCD數(shù)據(jù),得到營養(yǎng)鹽濃度衛(wèi)星產(chǎn)品專題圖(圖 2).NO2--N、NH4+-N、NO3--N以及TIN都表現(xiàn)為在大洋河河口有一個(gè)明顯的高值區(qū),該高值區(qū)營養(yǎng)鹽濃度由河口向外逐漸降低,直至外海最低.因此,大洋河河口海域營養(yǎng)鹽遙感反演具有可行性,從而為河口海域這一特殊區(qū)域水體中營養(yǎng)鹽的遙感反演研究提供參考.

    表3 河口海域理、化以及光學(xué)參數(shù)間相關(guān)性分析結(jié)果Table 3 Correlation analysis results among physical, chemical and optical parameters

    圖2 NO2--N、NH4+-N、NO3--N以及TIN衛(wèi)星產(chǎn)品專題圖Fig.2 NO2--N, NH4+-N, NO3--N and TIN satellite products

    3 結(jié)論

    3.1 HJ-1A CCD2 Band2波段為 NO2--N 和NH4+-N反演的敏感波段,Band1為 NO3--N 和TIN反演的敏感波段,由對(duì)數(shù)函數(shù)擬合NO2--N與Rrs(Band2)效果最佳,相關(guān)性R2可達(dá)0.939,而由線性函數(shù)擬合NH4+-N與Rrs(Band2)、NO3--N與Rrs(Band1)、TIN與 Rrs(Band1)效果最佳,相關(guān)性R2分別可達(dá)0.935、0.945、0.970,未發(fā)現(xiàn)較佳的PO43--P遙感反演模型.

    3.2 Sal與 NO2--N、NH4+-N、NO3--N 以及ag(400)存在較高的相關(guān)性,NO2--N、NH4+-N、 NO3--N與 ag(400)的相關(guān)性分別可達(dá) 0.836、 0.808、0.786,因此,河口海域陸源徑流輸入物質(zhì)典型的擴(kuò)散特征,及在此基礎(chǔ)上使得營養(yǎng)鹽參數(shù)與某些水色參數(shù)存在較高的相關(guān)性是進(jìn)行河口海域營養(yǎng)鹽遙感反演的重要依據(jù).

    3.3 從衛(wèi)星產(chǎn)品專題圖可看出,大洋河河口海域NO2--N、NH4+-N、NO3--N以及TIN都存在一個(gè)明顯的高值區(qū),該高值區(qū)營養(yǎng)鹽濃度由河口向外逐漸降低,直至外海最低.

    [1] 胡瑩瑩,王菊英,張志鋒,等.遼河口近岸海域水體營養(yǎng)物推薦基準(zhǔn)值的制定方法 [J]. 中國環(huán)境科學(xué), 2011,31(6):996-1000.

    [2] 李小斌,陳楚群,施 平,等.珠江口海域總無機(jī)氮的遙感提取研究 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2007,27(2):313-318.

    [3] Pekka H, Jenni V, Tuula H, et al. Detection of water quality using simulated satellite data and semi2emp irical algorithms in Finland [J]. The Science of the total Environment, 2001,268:107-121.

    [4] Dmitry P, Robert S, Anton K, et al. Operational algorithm for the retrieval of water quality in the GreatLakes [J]. Remote Sensing of Environment, 2005,97:352-370.

    [5] Ana Silió-Calzada, Annick Bricaud, Bernard Gentili. Estimates of sea surface nitrate concentrations from sea surface temperature and chlorophyll concentration in upwelling areas a case study for the Benguela system [J]. Remote Sensing of Environment, 2008, 112:3173-3180.

    [6] 許 勇.江蘇沿海無機(jī)氮磷營養(yǎng)鹽遙感監(jiān)測(cè)研究 [D]. 南京:南京師范大學(xué), 2009.

    [7] 劉 征,賀軍亮,彭 林,等.黃壁莊水庫總氮、總磷含量與反射光譜特征的關(guān)系 [J]. 石家莊學(xué)院學(xué)報(bào), 2009,11(3):45-55.

    [8] 趙旭陽,劉 征,賀軍亮,等.黃壁莊水庫水質(zhì)參數(shù)遙感反演研究[J]. 地理與地理信息科學(xué), 2007,23(6):46-49.

    [9] 鞏彩蘭,尹 球,匡定波.黃浦江水質(zhì)指標(biāo)與反射光譜特征的關(guān)系分析 [J]. 遙感學(xué)報(bào), 2006,10(6):69-75.

    [10] 雷 坤,鄭丙輝,王 橋.基于中巴地球資源1號(hào)衛(wèi)星的太湖表層水體水質(zhì)遙感 [J]. 環(huán)境科學(xué)學(xué)報(bào), 2004,24(3):376-380.

    [11] 王學(xué)軍,馬 廷.應(yīng)用遙感技術(shù)監(jiān)測(cè)和評(píng)價(jià)太湖水質(zhì)狀況 [J].環(huán)境科學(xué), 2000,21(6):65-68.

    [12] 王建平,程聲通,賈海峰,等.用 TM 影像進(jìn)行湖泊水色反演研究的人工神經(jīng)網(wǎng)絡(luò)模型 [J]. 環(huán)境科學(xué), 2003,24(2):18-24.

    [13] NASA/TM-2003-211621/Rev4 Ocean optics protocols for satellite ocean color sensor validation above-water radiance and remote sensing reflectance measurement and analysis protocols [S].

    [14] Mobley C D. Estimation of the remote-sensing reflectance from above-surface measurements [J]. Applied Optics, 1999,38(36): 7442-7455.

    [15] NASA/TM-2002-210004/Rev3 Ocean optics protocols for satellite ocean color sensor validation fluorometric chlorophyll a: sampling, laboratory methods, and data analysis protocols [S].

    [16] NASA/TM-2000-209966/Rev2 Ocean optics protocols for satellite ocean color sensor validation determination of spectral absorption coefficients of particles, dissolved material and phytoplankton for discrete water samples [S].

    [17] GB/T 17378.4-2007 海洋監(jiān)測(cè)規(guī)范 [S].

    [18] 環(huán)境減災(zāi)星座A/B星各載荷在軌絕對(duì)輻射定標(biāo)系數(shù) [R]. 北京:中國資源衛(wèi)星應(yīng)用中心, 2009.

    Nutrients retrieval in the Dayang River Estuary based on HJ-1 satellite remote-sensed imagery.

    WANG Lin, ZHAO Dong-zhi*, YANG Jian-hong (Department of Ocean Remote Sensing,National Marine Environmental Monitoring Center, Dalian 116023, China) . China Environmental Science, 2012,32(1):136~141

    On the basis of in situ measurement data in the Dayang River Estuary and HJ-1 satellite CCD imagery, the--N) concentration had the best correlation with the remote sensing reflectance (Rrs) in Band 2 (520~600nm), following a natural logarithmic function (R2=0.939). The ammonia (NH4+-N) concentration was best correlated with Rrsin Band 2 as well, following a linear function (R2=0.935). Besides, the nitrate (NO3--N) and total inorganic nitrogen (TIN) were best correlated with Rrsin Band 1 (430~520nm), likewise following linear functions (R2=0.945 and 0.970 respectively). However, good relationship between the phosphate (PO43--P) and Rrsin any band was not found. Theoretically, nutrients were not optically sensitive substance in the visible light band, thus these empirical relationships between nutrients and equivalent remote sensing reflectance were largely due to the typical diffusion characteristics of river input materials, which indirectly led to high correlativity between nutrients and some optical coefficients, like the absorption of yellow substance (ag(400)). The indirect correlation proposes the potential feasibility of nutrients retrieval in the estuary waters through satellite remote-sensed imagery. Applying the retrieval algorithm to the HJ-1 satellite data, a euphotic zone with extremely high concentration of nutrients was detected in the Dayang River Estuary, and various nutrients gradually declined from the estuary to open seas.

    HJ-1 satellite;Dayang River Estuary;nutrients;remote sensing retrieval nutrients retrieval algorithms were studied. The nitrite (NO2

    2011-04-16

    國家自然科學(xué)基金資助項(xiàng)目(40876091);國家海洋公益性行業(yè)科研專項(xiàng)(201005030);國家海洋公益性重點(diǎn)項(xiàng)目(200805064);國家“863”項(xiàng)目(2007AA092003);上海市科學(xué)技術(shù)委員會(huì)科研計(jì)劃(08DZ1206304)

    * 責(zé)任作者, 研究員, dzzhao@nmemc.gov.cn

    X145

    A

    1000-6923(2012)01-0136-06

    致謝:感謝國家海洋環(huán)境監(jiān)測(cè)中心海洋遙感室全體工作人員的協(xié)助.

    王 林(1981-),男,河北滄州人,研習(xí)員,碩士,主要從事海洋光學(xué)與水質(zhì)遙感研究.發(fā)表論文5篇.

    猜你喜歡
    營養(yǎng)鹽河口反射率
    影響Mini LED板油墨層反射率的因素
    近岸水體異源遙感反射率產(chǎn)品的融合方法研究
    具有顏色恒常性的光譜反射率重建
    涼水河子河營養(yǎng)鹽滯留能力評(píng)估
    化學(xué)腐蝕硅表面結(jié)構(gòu)反射率影響因素的研究*
    電子器件(2017年2期)2017-04-25 08:58:37
    瓊東海域冬季、夏季營養(yǎng)鹽結(jié)構(gòu)特征及其對(duì)浮游植物生長的影響
    他們?yōu)槭裁催x擇河口
    河口,我們的家
    特殊的河口水
    河口
    天堂动漫精品| 麻豆一二三区av精品| 欧美日韩黄片免| 色播亚洲综合网| 美女国产高潮福利片在线看| 久久精品国产亚洲av香蕉五月| 真人一进一出gif抽搐免费| 精品国产一区二区三区四区第35| 欧美最黄视频在线播放免费| 国产成人精品久久二区二区91| 免费在线观看视频国产中文字幕亚洲| 日本欧美视频一区| 日本vs欧美在线观看视频| 国产精品乱码一区二三区的特点 | 99精品欧美一区二区三区四区| 两性午夜刺激爽爽歪歪视频在线观看 | 性欧美人与动物交配| 国产亚洲精品久久久久5区| 精品人妻1区二区| 制服人妻中文乱码| 国产av又大| 久久香蕉精品热| 一二三四在线观看免费中文在| 亚洲精品在线观看二区| 午夜免费观看网址| 男男h啪啪无遮挡| 黄片大片在线免费观看| 黑人巨大精品欧美一区二区mp4| 精品熟女少妇八av免费久了| 女性被躁到高潮视频| 久久这里只有精品19| 国产主播在线观看一区二区| 丁香欧美五月| 国产精品一区二区精品视频观看| 亚洲国产精品久久男人天堂| 久久久久久亚洲精品国产蜜桃av| 久久精品91无色码中文字幕| 91精品三级在线观看| 国产xxxxx性猛交| 又紧又爽又黄一区二区| 看免费av毛片| www国产在线视频色| 久久精品人人爽人人爽视色| 一区福利在线观看| 亚洲av成人一区二区三| 中文字幕av电影在线播放| 可以免费在线观看a视频的电影网站| 母亲3免费完整高清在线观看| 亚洲激情在线av| 亚洲va日本ⅴa欧美va伊人久久| 人成视频在线观看免费观看| 在线永久观看黄色视频| 国产精品久久久av美女十八| 国产三级在线视频| 国产在线观看jvid| 在线观看免费视频网站a站| 亚洲自拍偷在线| 亚洲国产欧美日韩在线播放| 久久婷婷成人综合色麻豆| 亚洲熟女毛片儿| 国产亚洲精品一区二区www| 美女扒开内裤让男人捅视频| 51午夜福利影视在线观看| 正在播放国产对白刺激| 香蕉丝袜av| 色综合亚洲欧美另类图片| 精品国产超薄肉色丝袜足j| 欧美乱色亚洲激情| 国产精品秋霞免费鲁丝片| 大型av网站在线播放| 黄片小视频在线播放| 天天躁夜夜躁狠狠躁躁| 日韩欧美在线二视频| 制服诱惑二区| 欧美成人午夜精品| 亚洲人成伊人成综合网2020| 久久亚洲真实| 午夜影院日韩av| 一级黄色大片毛片| 久久人人爽av亚洲精品天堂| 美女午夜性视频免费| 欧美日本亚洲视频在线播放| 国产av在哪里看| 色哟哟哟哟哟哟| 国产av在哪里看| 久久狼人影院| 久久亚洲精品不卡| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久水蜜桃国产精品网| 黄片小视频在线播放| 女性生殖器流出的白浆| 午夜a级毛片| 亚洲最大成人中文| 欧美日韩亚洲国产一区二区在线观看| 波多野结衣一区麻豆| АⅤ资源中文在线天堂| 国产在线观看jvid| 熟妇人妻久久中文字幕3abv| 成人亚洲精品av一区二区| 波多野结衣高清无吗| 日韩国内少妇激情av| 男人的好看免费观看在线视频 | 午夜免费激情av| 欧美黑人欧美精品刺激| АⅤ资源中文在线天堂| av有码第一页| 亚洲男人天堂网一区| videosex国产| 久久久国产精品麻豆| 免费高清在线观看日韩| 亚洲av日韩精品久久久久久密| 国产精品久久久久久亚洲av鲁大| 18禁观看日本| 91精品国产国语对白视频| 亚洲成av片中文字幕在线观看| 欧美在线黄色| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品综合久久久久久久免费 | 精品欧美一区二区三区在线| 国产国语露脸激情在线看| 国产亚洲欧美98| 国产成人av教育| 日韩欧美免费精品| 精品久久久久久久人妻蜜臀av | 亚洲精品一区av在线观看| 久久精品91蜜桃| 亚洲一区二区三区不卡视频| 久久精品亚洲精品国产色婷小说| 亚洲国产日韩欧美精品在线观看 | 亚洲中文av在线| 亚洲av电影在线进入| 亚洲av日韩精品久久久久久密| 女性生殖器流出的白浆| 久久久久久国产a免费观看| 99国产综合亚洲精品| 精品乱码久久久久久99久播| 日本五十路高清| 妹子高潮喷水视频| 亚洲成人久久性| 久久久久久免费高清国产稀缺| 亚洲成人免费电影在线观看| 欧美一区二区精品小视频在线| 亚洲国产毛片av蜜桃av| 亚洲精品在线美女| 欧美一区二区精品小视频在线| 免费在线观看日本一区| 国产aⅴ精品一区二区三区波| 国产人伦9x9x在线观看| 午夜福利成人在线免费观看| 免费在线观看黄色视频的| 一级,二级,三级黄色视频| 精品国产乱码久久久久久男人| 亚洲av第一区精品v没综合| 51午夜福利影视在线观看| 色播亚洲综合网| 国产精华一区二区三区| 18禁黄网站禁片午夜丰满| 国产成人精品久久二区二区91| 国产精品av久久久久免费| 日韩成人在线观看一区二区三区| 午夜福利一区二区在线看| 国产蜜桃级精品一区二区三区| 亚洲第一青青草原| 亚洲男人天堂网一区| 巨乳人妻的诱惑在线观看| 精品久久久久久,| 亚洲,欧美精品.| 久久久久亚洲av毛片大全| 久久精品亚洲熟妇少妇任你| 一区福利在线观看| 午夜免费成人在线视频| 亚洲午夜理论影院| 中文字幕精品免费在线观看视频| 日韩国内少妇激情av| 一边摸一边抽搐一进一出视频| 亚洲色图综合在线观看| 久久午夜综合久久蜜桃| 久久国产精品影院| 精品人妻在线不人妻| 高清黄色对白视频在线免费看| 国产单亲对白刺激| 成人国产综合亚洲| 午夜福利高清视频| 午夜福利,免费看| a在线观看视频网站| 午夜福利成人在线免费观看| 成年版毛片免费区| 国产不卡一卡二| 18禁美女被吸乳视频| 少妇的丰满在线观看| 精品乱码久久久久久99久播| 一级,二级,三级黄色视频| 成人18禁在线播放| 国产精品98久久久久久宅男小说| 日韩免费av在线播放| 美女大奶头视频| 久久国产精品男人的天堂亚洲| 久久久久久久久免费视频了| 国产一区二区三区在线臀色熟女| 一二三四社区在线视频社区8| 99精品欧美一区二区三区四区| 在线观看免费视频网站a站| 美女免费视频网站| 欧美日本亚洲视频在线播放| 日本 av在线| 国产成人精品无人区| av网站免费在线观看视频| 亚洲国产中文字幕在线视频| 一级a爱视频在线免费观看| av视频免费观看在线观看| av网站免费在线观看视频| 一进一出抽搐动态| 欧美性长视频在线观看| 一区福利在线观看| 亚洲五月婷婷丁香| 日韩国内少妇激情av| 十八禁网站免费在线| 老司机在亚洲福利影院| 午夜福利18| 在线观看免费日韩欧美大片| 欧美日韩瑟瑟在线播放| 淫秽高清视频在线观看| 久久热在线av| 国产av一区二区精品久久| 非洲黑人性xxxx精品又粗又长| 少妇 在线观看| videosex国产| 中文字幕久久专区| 久久精品国产清高在天天线| 精品午夜福利视频在线观看一区| 少妇粗大呻吟视频| 亚洲五月婷婷丁香| 黄色成人免费大全| 99国产综合亚洲精品| 欧美久久黑人一区二区| 国产av一区二区精品久久| 很黄的视频免费| 一区二区三区激情视频| 午夜亚洲福利在线播放| 日日夜夜操网爽| 两个人免费观看高清视频| 免费不卡黄色视频| 999久久久精品免费观看国产| 国产高清有码在线观看视频 | 久久香蕉国产精品| 黄色a级毛片大全视频| 午夜a级毛片| 老司机午夜十八禁免费视频| 久久久久久免费高清国产稀缺| 少妇裸体淫交视频免费看高清 | 欧美日韩一级在线毛片| 久久精品91蜜桃| 两个人视频免费观看高清| 国产1区2区3区精品| 在线永久观看黄色视频| 国产精品久久电影中文字幕| 色在线成人网| 亚洲人成伊人成综合网2020| 日本免费一区二区三区高清不卡 | 纯流量卡能插随身wifi吗| a级毛片在线看网站| av电影中文网址| 91精品三级在线观看| 午夜久久久在线观看| 亚洲少妇的诱惑av| 热99re8久久精品国产| 少妇 在线观看| 性色av乱码一区二区三区2| 免费在线观看日本一区| 国产成年人精品一区二区| 欧美日韩亚洲国产一区二区在线观看| 亚洲国产精品999在线| 麻豆国产av国片精品| 国产精品免费视频内射| 啪啪无遮挡十八禁网站| 亚洲一区中文字幕在线| 日日干狠狠操夜夜爽| 91精品国产国语对白视频| 欧美激情 高清一区二区三区| 精品日产1卡2卡| 精品人妻1区二区| 好男人在线观看高清免费视频 | 91麻豆精品激情在线观看国产| 亚洲一区高清亚洲精品| 男女午夜视频在线观看| 日本vs欧美在线观看视频| 麻豆久久精品国产亚洲av| 麻豆av在线久日| 搞女人的毛片| 美女高潮到喷水免费观看| 国产av一区在线观看免费| 色在线成人网| 最新在线观看一区二区三区| 欧美日本亚洲视频在线播放| 97超级碰碰碰精品色视频在线观看| 久久久久久免费高清国产稀缺| 一级片免费观看大全| 欧美绝顶高潮抽搐喷水| or卡值多少钱| 亚洲国产欧美网| 99国产精品免费福利视频| 亚洲精品国产一区二区精华液| 国产精品久久久久久亚洲av鲁大| 国产成人av教育| 国产成人欧美| 大香蕉久久成人网| 亚洲电影在线观看av| 欧美成人性av电影在线观看| 少妇熟女aⅴ在线视频| 成在线人永久免费视频| 看黄色毛片网站| 亚洲久久久国产精品| 黄色丝袜av网址大全| 国产野战对白在线观看| 久久久久久久久久久久大奶| 夜夜躁狠狠躁天天躁| 乱人伦中国视频| 90打野战视频偷拍视频| 可以在线观看的亚洲视频| 精品一区二区三区av网在线观看| 国产精品99久久99久久久不卡| 免费看a级黄色片| 久久精品成人免费网站| 人人澡人人妻人| 亚洲美女黄片视频| 女人被狂操c到高潮| 亚洲精品在线观看二区| 97超级碰碰碰精品色视频在线观看| 91字幕亚洲| 午夜精品在线福利| 动漫黄色视频在线观看| 午夜精品在线福利| 12—13女人毛片做爰片一| 两性夫妻黄色片| 亚洲三区欧美一区| 亚洲成av人片免费观看| 久久久久久免费高清国产稀缺| 搞女人的毛片| 嫩草影院精品99| 热re99久久国产66热| 亚洲中文av在线| 欧美最黄视频在线播放免费| 亚洲男人天堂网一区| 久久精品国产99精品国产亚洲性色 | 久久久久精品国产欧美久久久| 国产欧美日韩一区二区三区在线| 久久精品国产亚洲av高清一级| 欧美最黄视频在线播放免费| 村上凉子中文字幕在线| 一级,二级,三级黄色视频| 久久影院123| ponron亚洲| 两个人视频免费观看高清| 午夜成年电影在线免费观看| 夜夜看夜夜爽夜夜摸| АⅤ资源中文在线天堂| 精品久久久久久成人av| 成人国产一区最新在线观看| 午夜亚洲福利在线播放| 午夜成年电影在线免费观看| 一进一出抽搐gif免费好疼| 午夜亚洲福利在线播放| 成人18禁在线播放| 久久久久久人人人人人| 成人免费观看视频高清| 非洲黑人性xxxx精品又粗又长| 国产熟女xx| 久久精品国产亚洲av香蕉五月| 可以在线观看的亚洲视频| 欧美不卡视频在线免费观看 | √禁漫天堂资源中文www| 成人av一区二区三区在线看| 国产精品av久久久久免费| 国产精品一区二区在线不卡| 97碰自拍视频| 日韩大码丰满熟妇| 看片在线看免费视频| 婷婷六月久久综合丁香| bbb黄色大片| 欧美不卡视频在线免费观看 | 亚洲成人免费电影在线观看| 亚洲九九香蕉| 色综合婷婷激情| 黑丝袜美女国产一区| 欧美乱码精品一区二区三区| а√天堂www在线а√下载| 岛国在线观看网站| 一级片免费观看大全| 久热爱精品视频在线9| 亚洲五月天丁香| 精品少妇一区二区三区视频日本电影| 成年人黄色毛片网站| 每晚都被弄得嗷嗷叫到高潮| 18禁黄网站禁片午夜丰满| 国产精品久久久久久亚洲av鲁大| 国产单亲对白刺激| 欧美黑人精品巨大| 欧美久久黑人一区二区| 国产一区二区激情短视频| 在线国产一区二区在线| 18禁美女被吸乳视频| 岛国在线观看网站| 亚洲精品久久国产高清桃花| 麻豆久久精品国产亚洲av| 国产欧美日韩综合在线一区二区| 亚洲伊人色综图| 97超级碰碰碰精品色视频在线观看| 亚洲自拍偷在线| 精品日产1卡2卡| 男女做爰动态图高潮gif福利片 | 久久精品国产亚洲av香蕉五月| 正在播放国产对白刺激| 久久精品影院6| 看片在线看免费视频| 成在线人永久免费视频| 一边摸一边抽搐一进一出视频| a级毛片在线看网站| 夜夜看夜夜爽夜夜摸| 国产熟女午夜一区二区三区| av中文乱码字幕在线| 久久久久久亚洲精品国产蜜桃av| 欧美在线黄色| 久久国产精品人妻蜜桃| 变态另类丝袜制服| av电影中文网址| 国产成人精品久久二区二区91| 国产精品99久久99久久久不卡| 高清毛片免费观看视频网站| 午夜福利高清视频| 成人国语在线视频| 国产精品久久久人人做人人爽| 久久热在线av| 成人亚洲精品一区在线观看| 欧美精品啪啪一区二区三区| 国产成人影院久久av| 婷婷六月久久综合丁香| 午夜久久久久精精品| 十八禁网站免费在线| 男男h啪啪无遮挡| 欧美激情高清一区二区三区| 亚洲av五月六月丁香网| 黄片大片在线免费观看| 成人国产一区最新在线观看| 大码成人一级视频| 国产一卡二卡三卡精品| 国产激情欧美一区二区| 丝袜美腿诱惑在线| 国产精品免费视频内射| 午夜精品久久久久久毛片777| 久久人人97超碰香蕉20202| 国产精品亚洲av一区麻豆| 国产主播在线观看一区二区| 国内精品久久久久精免费| 91成年电影在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久视频播放| 老汉色∧v一级毛片| 国产成人av教育| 国产成人系列免费观看| 最新美女视频免费是黄的| 国产黄a三级三级三级人| av电影中文网址| 久久久久久久久久久久大奶| 免费在线观看影片大全网站| 首页视频小说图片口味搜索| 亚洲激情在线av| 中文字幕av电影在线播放| 日韩视频一区二区在线观看| 国产精品一区二区精品视频观看| 精品欧美国产一区二区三| 午夜免费成人在线视频| 精品久久久久久成人av| 亚洲熟妇中文字幕五十中出| 亚洲一区二区三区不卡视频| 久久精品成人免费网站| 成人国产一区最新在线观看| 精品欧美一区二区三区在线| 午夜福利,免费看| 午夜日韩欧美国产| 69av精品久久久久久| 亚洲熟女毛片儿| 美女 人体艺术 gogo| 国产精品久久久av美女十八| 欧美久久黑人一区二区| 极品教师在线免费播放| 女同久久另类99精品国产91| 日本一区二区免费在线视频| 国产在线观看jvid| 黄片播放在线免费| 精品国产乱码久久久久久男人| 亚洲欧美精品综合一区二区三区| 国产熟女午夜一区二区三区| 亚洲自拍偷在线| 日韩欧美国产在线观看| 成人免费观看视频高清| 欧美日韩精品网址| 欧美日韩福利视频一区二区| 国产精品久久久av美女十八| 黄片播放在线免费| 免费观看人在逋| 一二三四在线观看免费中文在| 国产欧美日韩一区二区精品| 国产野战对白在线观看| 最新美女视频免费是黄的| or卡值多少钱| 一级黄色大片毛片| 国产麻豆成人av免费视频| 久久人人97超碰香蕉20202| 国产97色在线日韩免费| 日韩欧美一区二区三区在线观看| 日日爽夜夜爽网站| 日日干狠狠操夜夜爽| 国产一区二区三区视频了| 欧美丝袜亚洲另类 | 亚洲国产精品久久男人天堂| 国产精品九九99| 国产成+人综合+亚洲专区| 一个人观看的视频www高清免费观看 | 日本欧美视频一区| 高潮久久久久久久久久久不卡| 欧美激情高清一区二区三区| 精品国产亚洲在线| 国产97色在线日韩免费| 免费在线观看完整版高清| 精品卡一卡二卡四卡免费| 中亚洲国语对白在线视频| 激情视频va一区二区三区| 无人区码免费观看不卡| 国产免费av片在线观看野外av| 亚洲色图综合在线观看| 免费人成视频x8x8入口观看| 欧美一区二区精品小视频在线| 午夜福利成人在线免费观看| 国产精品九九99| 欧美日韩精品网址| 热99re8久久精品国产| 久久精品国产亚洲av香蕉五月| 成人精品一区二区免费| 久久精品91无色码中文字幕| 99国产极品粉嫩在线观看| 亚洲国产精品久久男人天堂| 午夜福利免费观看在线| 国产精品一区二区精品视频观看| 美女 人体艺术 gogo| av视频免费观看在线观看| 午夜日韩欧美国产| 十八禁人妻一区二区| 久久人人精品亚洲av| 黄片播放在线免费| 国产高清有码在线观看视频 | 波多野结衣巨乳人妻| 中亚洲国语对白在线视频| 久久久久亚洲av毛片大全| av天堂久久9| 香蕉丝袜av| 亚洲中文字幕一区二区三区有码在线看 | 悠悠久久av| 亚洲欧美精品综合久久99| 91精品三级在线观看| 又黄又爽又免费观看的视频| 别揉我奶头~嗯~啊~动态视频| 一区福利在线观看| 在线国产一区二区在线| 桃红色精品国产亚洲av| 成人三级黄色视频| 中文字幕av电影在线播放| 亚洲精华国产精华精| 韩国精品一区二区三区| www.熟女人妻精品国产| 精品国产乱码久久久久久男人| 国产一区二区三区视频了| 午夜福利成人在线免费观看| 久久久国产精品麻豆| 亚洲人成网站在线播放欧美日韩| 欧美中文日本在线观看视频| 日韩三级视频一区二区三区| 亚洲欧美精品综合一区二区三区| aaaaa片日本免费| 人人妻人人爽人人添夜夜欢视频| 久久久久国产精品人妻aⅴ院| bbb黄色大片| 18禁国产床啪视频网站| 亚洲专区字幕在线| 一本久久中文字幕| 国产真人三级小视频在线观看| 日日夜夜操网爽| 1024视频免费在线观看| 午夜福利18| 亚洲精品av麻豆狂野| 这个男人来自地球电影免费观看| 欧美av亚洲av综合av国产av| 一级毛片精品| 十八禁人妻一区二区| 啦啦啦免费观看视频1| 精品高清国产在线一区| 午夜福利视频1000在线观看 | 国内精品久久久久精免费| 久久久久久久久久久久大奶| 涩涩av久久男人的天堂| 欧美最黄视频在线播放免费| 欧美成人一区二区免费高清观看 | 久久国产精品人妻蜜桃| 满18在线观看网站| 国产精品影院久久| 亚洲欧美日韩高清在线视频| cao死你这个sao货| 女人高潮潮喷娇喘18禁视频| 免费在线观看日本一区| 免费看十八禁软件| 精品久久久久久,| 看片在线看免费视频| 禁无遮挡网站| 99re在线观看精品视频| 久久精品国产综合久久久| 麻豆国产av国片精品|