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

    基于GIS的桑溝灣及周圍海域海帶養(yǎng)殖適宜性評價*

    2020-02-13 07:01:32孫倩雯尚偉濤于良巨姜曉鵬蔡碧瑩常麗榮肖露陽
    漁業(yè)科學進展 2020年1期
    關鍵詞:海帶鹽度圖層

    孫倩雯 劉 慧 尚偉濤 于良巨 姜曉鵬 蔡碧瑩 常麗榮 肖露陽

    基于GIS的桑溝灣及周圍海域海帶養(yǎng)殖適宜性評價*

    孫倩雯1,2劉 慧2①尚偉濤3,4于良巨4姜曉鵬3,4蔡碧瑩1,2常麗榮5肖露陽5

    (1. 上海海洋大學水產與生命學院 上海 201306; 2. 中國水產科學研究院黃海水產研究所 青島 266071; 3. 中國科學院大學 北京 100049; 4. 中國科學院海岸帶環(huán)境過程與生態(tài)修復重點實驗室 中國科學院煙臺海岸帶研究所 煙臺 264003; 5. 威海長青海洋科技股份有限公司 榮成 264316)

    目前水產養(yǎng)殖面臨用海沖突、單位面積和人均生產率偏低、養(yǎng)殖業(yè)的生態(tài)和經濟成本較高等問題。為充分利用海區(qū)的自然生產力,提高養(yǎng)殖效率和效益,需要對海區(qū)進行適宜性評價,選擇最適宜的水域進行養(yǎng)殖。以桑溝灣及周圍海域主要的養(yǎng)殖品種海帶()作為評價對象,選取光照、溫度、流速、無機氮、鹽度、深度作為適宜性評價指標,利用遙感技術確定海區(qū)養(yǎng)殖布局,通過野外調查和數(shù)據模擬獲取養(yǎng)殖海區(qū)環(huán)境參數(shù)數(shù)據,根據動態(tài)能量學(Dynamic energy budget, DEB)模型——STELLA模型敏感性分析結果,并結合層次分析法計算評價指標權重,海帶生長相關環(huán)境因子的強制函數(shù)擬合得到單因子評分曲線進行評分。最后,采用線性加權疊加分析方法得到海帶養(yǎng)殖適宜性評價結果。適宜性評價過程以地理信息系統(tǒng)(Geographical information systems, GIS)作為技術支持,利用GIS空間插值功能生成光照、溫度、流速、無機氮、鹽度、深度專題圖層,采用柵格計算功能對各專題圖層進行疊加,得到綜合多因素的海帶養(yǎng)殖適宜性評分和適宜性等級。評價結果顯示,在不考慮用海沖突的情況下,桑溝灣及周圍海域海帶養(yǎng)殖適宜性分數(shù)分布在0~6.7范圍內,分值主要集中在4~6分,屬于中等適宜,占研究區(qū)總面積的67%,最適宜和不適宜分別占研究區(qū)總面積的23%和10%,無一般適宜區(qū)。評分較高的區(qū)域主要分布在北部愛蓮灣和楮島東部海域周圍,適宜性評分較低的區(qū)域主要位于近岸水深較淺的海域。根據適宜性評價結果可進行水域空間規(guī)劃和安排養(yǎng)殖生產計劃,為管理部門開展養(yǎng)殖分區(qū)和海洋功能區(qū)劃提供參考。

    水產養(yǎng)殖;適宜性評價;養(yǎng)殖區(qū)選址;GIS

    養(yǎng)殖容量的基礎是水域承載力的計算,傳統(tǒng)的養(yǎng)殖容量評估一般是基于營養(yǎng)鹽的供應或初級生產力水平,來估算水域可以養(yǎng)殖的藻類或者貝類的生物量(方建光等, 1996)。但養(yǎng)殖容量并未考慮養(yǎng)殖生物的其他生理需求,例如光照、溫度、鹽度等,而這些條件對于養(yǎng)殖生物是否能夠健康快速生長同樣重要。養(yǎng)殖適宜性評價通過全面評估生物的環(huán)境適應性,來選擇最適合養(yǎng)殖的水域,與養(yǎng)殖容量評估互為補充,可以更好地指導養(yǎng)殖規(guī)劃。

    未來水產養(yǎng)殖規(guī)劃管理將基于生態(tài)學方法,平衡社會、經濟、環(huán)境關系,盡可能縮小各因素間的不利影響(FAO, 2017)。Gentry等(2017)利用生理、異速生長和生長理論的創(chuàng)新方法,對全球海水養(yǎng)殖潛力進行了量化評估,結果顯示,中國的養(yǎng)殖產量超過預計潛力,反映出中國已開發(fā)大面積的海域發(fā)展水產養(yǎng)殖,存在超容量養(yǎng)殖。超高密度的水產養(yǎng)殖會帶來不利的環(huán)境影響,造成海域使用沖突,并產生高昂的社會經濟損失(Ottinger, 2016)。因此,為充分利用水域空間和自然生產力,增加養(yǎng)殖產量和生產效率,管理部門和養(yǎng)殖企業(yè)都需要對養(yǎng)殖區(qū)開展適宜性評價,從而了解水域適合養(yǎng)殖的品種和養(yǎng)殖規(guī)模,以便進行合理的空間規(guī)劃與布局。

    水產養(yǎng)殖品種的存活和生長在很大程度上依靠水域自然環(huán)境,水文、氣候、水化學要素以及初級生產力等自然環(huán)境因子是影響?zhàn)B殖品種生長的主要因素。因此,水產養(yǎng)殖適宜性評價以國家政策法規(guī)、海洋功能區(qū)劃和水環(huán)境標準為依據,通過養(yǎng)殖生物的生理生態(tài)特性、生長所需的環(huán)境條件和養(yǎng)殖水域的環(huán)境要素進行分析比較,選擇適合養(yǎng)殖生物生長的水域,為進行科學的養(yǎng)殖布局提供參考。由于養(yǎng)殖適宜性評價涉及的數(shù)據量大,故普遍采用地理信息系統(tǒng)(Geographical information systems, 簡稱GIS) (Foresman, 1998)的邏輯判斷、評價分析和可視化展示功能,進行空間插值和專題圖層疊加,對水產養(yǎng)殖區(qū)進行適宜性評價分析。在適宜性評價的基礎上,通過進一步整合模型運算和環(huán)境參數(shù),還可以對特定水域的水產養(yǎng)殖現(xiàn)狀及發(fā)展前景做出科學合理的評價與預測(劉慧等, 2018)。

    1 研究方法

    1.1 研究區(qū)和數(shù)據來源

    1.1.1 研究區(qū)概況 桑溝灣是位于山東半島東端的半封閉海灣(37°01¢~37°09¢N,122°24¢~122°35¢E),面積為133.3 km2,桑溝灣及周圍愛蓮灣和楮島部分海域總面積約為448 km2(36°58¢~37°12¢N, 122°23¢E~ 122°42¢) (圖1)。其地理環(huán)境優(yōu)越,海底平坦,水質較好,灣內平均水深為7.5 m,有不規(guī)則半日潮(張朝暉等, 2007)。年平均溫度為13℃,2月水溫最低,月平均溫度為1.8℃;8月水溫最高,月平均溫度為24℃。環(huán)境調查發(fā)現(xiàn),桑溝灣內全年海水鹽度范圍為29~32,太陽輻射強度為75~313 w/m2,平均日照時長為12 h (史潔, 2009),適宜發(fā)展水產養(yǎng)殖。海帶()作為當?shù)刂饕B(yǎng)殖品種,年產量達10萬t (干重)(蔡碧瑩等, 2019),且商品化程度高,因此選取海帶作為適宜性評價對象。

    圖1 海帶養(yǎng)殖適宜性評價研究區(qū)地理位置及主要等深線

    (數(shù)據來源: Copernicus Open Access Hub. https://scihub.copernicus.eu/)

    Fig.1 The location and depth contours of the studied area for kelp aquaculture suitability assessment

    (Source: Copernicus Open Access Hub. https://scihub.copernicus.eu/)

    1.1.2 數(shù)據來源 根據海帶生長發(fā)育的生理需求,選取光照、溫度、流速、無機氮、鹽度、深度6個環(huán)境參數(shù)用于養(yǎng)殖適宜性評價。其中,桑溝灣及周圍海域地理信息數(shù)據由遙感獲取(Source: Copernicus Open Access Hub. https://scihub.copernicus.eu/);流速和水溫數(shù)據來源于FVCOM模型(Finite volume coastal ocean model)基于動力方程對桑溝灣模擬得到(Xuan, 2016);光照采用養(yǎng)殖海區(qū)平均海表光照強度,通過晴雨天數(shù)以及每月日照時長計算獲得,數(shù)據來源于2011年中國氣象局榮成市天氣預報;無機氮和鹽度數(shù)據來源于對桑溝灣及周圍海域2011年4個季度(4、8、10月、2012年1月)的大面調查。

    1.2 評價指標的選擇

    適宜性評價指標根據海帶生理生態(tài)學特點,并結合桑溝灣養(yǎng)殖區(qū)地理和水環(huán)境特點進行篩選。影響海帶生長的主要因素有光照、溫度、流速、無機氮、鹽度、深度。光照對海帶生長有重要的影響,光照強弱影響海帶光合作用,且光能夠刺激大型海藻對營養(yǎng)鹽的吸收(Lobban, 1994)。深度通過影響海帶接受光照,進而影響海帶生長,過淺的水深,會限制海帶長度的生長。海帶生長受溫度影響明顯,溫度超過17.5℃時,海帶葉片末梢枯爛率超過其生長率,當溫度降低時,長海帶()對NO3?的吸收速率降低(Suzuki, 2008; Harlin, 1978);水動力是大型藻類生產的一個關鍵因素,在其他環(huán)境因素不受限時,若大型藻類生長的主流速度持續(xù)低于飽和水平,其生產率降低,且較弱的水動力影響營養(yǎng)鹽的補充(Leigh, 1987; Wheeler, 1980)。無機氮是海帶生長所需的主要營養(yǎng)鹽,藻類利用氮元素合成自身生長的蛋白質(Lobban, 1996)。鹽度變化影響海水藻類光合作用同化速率(王憲等, 1991)。這些環(huán)境因子的變化對海帶生長產生重要的影響,因此,將其作為環(huán)境因素的評價指標,并匯總各評價指標對于海帶生長的參數(shù)范圍(表1)。

    1.3 權重計算

    不同的環(huán)境因子對海帶生長的影響程度不同,需要根據各評價指標的重要程度賦予相應的權重。深度作為限制因素,不賦予權重;光照、溫度、流速、無機氮、鹽度的重要性等級排序參考海帶生長動態(tài)能量學(Dynamic energy budget, DEB)模型——STELLA模型敏感性分析結果(表2) (蔡碧瑩等, 2019)。DEB模型可描述海帶生長與環(huán)境參數(shù)的關系,模擬結果與實測值擬合度高,可較好地反演海帶真實的生長過程 (蔡碧瑩等, 2019)。通過敏感性分析得出各環(huán)境因素對海帶生長的影響程度排序,結合層次分析法(AHP) (Saaty, 1977)構造對比矩陣,計算各環(huán)境因子的權重(表4)。

    表1 海帶適宜性評價指標

    Tab.1 The criteria for assessment of culture suitability of kelp S. japonica

    1.4 評分

    1.4.1 單因子評分 影響海帶生長的環(huán)境因子評分采用8分制,1~8分表示環(huán)境條件對海帶養(yǎng)殖適宜程度由低到高。采用海帶生長相關環(huán)境參數(shù)的強制函數(shù)來擬合評分曲線。強制函數(shù)表示不同環(huán)境因子對海帶生長速率的影響,每個環(huán)境參數(shù)的強制函數(shù)均來自下列經驗公式:

    溫度評分曲線采用經驗公式溫度方程(EPA, 1985):

    式中,x為溫度生態(tài)幅,opt為最適宜溫度,當opt時,x=max(溫度生態(tài)幅上限)。

    鹽度評分曲線采用鹽度限制函數(shù)(Martins, 2002):

    式中,opt為最適鹽度,當

    流速評分方程采用藻類生長模型和水動力模型相耦合的流速曲線(Barr, 2008)

    式中,為流速。

    海帶對氮的吸收特征符合飽和吸收動力學,因此,無機氮方程采用模擬Monod或Michaelis-Menten動力學方程(Kitadai, 2003):

    式中,N為總氮濃度,c為N的半飽和同化系數(shù)。

    光照評分采用光照限制函數(shù)(Steele, 1962):

    式中,為海帶表面光強,opt為海帶光合作用最適光強。

    將表1中的參數(shù)值代入到對應的經驗公式中得到單因子評分曲線(圖2),根據曲線得到1~8分中每個分數(shù)段對應的參數(shù)范圍,以此作為海帶養(yǎng)殖適宜性評分的依據。深度作為限制因素,深度<5 m,分值設置為0;深度>5 m,分值設置為1。

    表2 海帶對各環(huán)境參數(shù)變化的敏感度(蔡碧瑩等, 2019)

    Tab.2 Sensitivity of kelp to the variation of environmental parameters (Cai et al, 2019)

    將桑溝灣及周圍海域光照、溫度、流速、無機氮、鹽度、深度的實測和模擬數(shù)據導入ArcGIS,采用空間插值方法生成各環(huán)境參數(shù)對應的專題圖層。按照評分依據,對各個環(huán)境參數(shù)的專題圖層進行重分類處理,每幅專題圖層中,按照評分標準賦予各參數(shù)范圍對應的分值,生成光照、溫度、流速、無機氮、鹽度、深度6幅單因子評分圖層。

    1.4.2 多指標綜合評分 水產養(yǎng)殖適宜性評價涉及到多種因素,由于這些因素的重要程度不同,因而不能將指標圖層簡單地進行疊加。采用線性加權疊加分析(Malczewski, 2000)計算綜合適宜性評分,計算公式如下:

    式中,Ai為評價對象像元的適宜性分數(shù),為指標的權重,為像元指標的分數(shù)。

    將光照、溫度、流速、無機氮、鹽度、深度單因子評分圖層進行重采樣,根據疊加分析,利用ArcGIS柵格計算功能將各單因子評分圖層進行疊加,輸出 4個季度的適宜性評分圖層,圖層顯示,綜合各因素后每個季度的適宜性評分(1~8)。將每個季度的適宜性評分重新分類,劃分為最適宜、中等適宜、一般適宜和不適宜4個等級(表3),不同的適宜性等級表示海帶養(yǎng)殖不同的適宜程度。最后,綜合4個季度的適宜性評分,進行柵格圖層疊加,生成最終適宜性評分圖層和適宜性等級圖層。

    表3 海帶養(yǎng)殖適宜性等級劃分

    Tab.3 The classification for culture suitability level of the kelp

    表4 海帶養(yǎng)殖適宜性評價指標權重

    Tab.4 The weighted indices of kelp culture suitability

    2 結果

    海帶養(yǎng)殖適宜性評價結果顯示,桑溝灣及其周圍海域適宜性分數(shù)分布在0~6.7范圍內,由桑溝灣近岸到離岸較遠海域,適宜性分值逐漸增大后減小,評分較高的區(qū)域主要分布在北部愛蓮灣和楮島東部海域,適宜性評分較低的區(qū)域主要位于近岸海域。適宜性分值主要集中在4~6分,屬于中等適宜(298.9 km2),占研究區(qū)總面積的67%,最適宜(103.7 km2)和不適宜(45.4 km2)分別占研究區(qū)總面積的23%和10%,沒有分值分布在一般適宜(圖3和圖5)。

    2011年4個季度適宜性評價結果顯示,春季適宜性分值較高的區(qū)域呈塊狀分布,夏季適宜性評分呈階梯狀,評分分值由桑溝灣灣內向灣外方向遞增,且分值差距最大,秋季適宜性高分區(qū)主要分布在愛蓮灣和楮島東部海域附近,冬季較春、夏、秋季,高分值區(qū)向桑溝灣灣內遷移(圖4a)。春、秋和冬季中等適宜海區(qū)面積占比重最大,分別占研究區(qū)總面積的56%、86%和78%,夏季最適宜海區(qū)面積最大,約占研究區(qū)總面積的65%,海帶養(yǎng)殖期間,4個季節(jié)不適宜養(yǎng)殖區(qū)面積均為45.4 km2,占研究區(qū)總面積的10% (圖4b和圖5)。

    適宜性評價結果與現(xiàn)有養(yǎng)殖布局對比顯示,養(yǎng)殖區(qū)分布密集,養(yǎng)殖筏架主要分布在中等適宜區(qū)域,約6 km2的筏架處于不適宜區(qū),19.8 km2的筏架處于最適區(qū)。從遙感影像顯示,約62.7 km2的最適宜區(qū)域并未安置養(yǎng)殖筏架(圖6)。實際養(yǎng)殖中,楮島東部海域最適宜區(qū)域海帶長度約為4 m,而楮島近岸不適宜養(yǎng)殖區(qū)的海帶長度最多長到3 m。

    圖3 2011年桑溝灣及周圍海域海帶養(yǎng)殖適宜性評分和適宜性等級

    圖4 2011年4個季度桑溝灣及周圍海帶養(yǎng)殖適宜性評分(a)和適宜性等級(b)

    3 討論

    現(xiàn)有養(yǎng)殖布局(圖6)顯示,桑溝灣及其周圍海域的養(yǎng)殖已經達到一定的規(guī)模,生產面臨挑戰(zhàn)。通過對海帶適宜性評價,確定海帶最適宜的養(yǎng)殖區(qū)域,以期獲得最佳的養(yǎng)殖結果,同時,可緩解現(xiàn)有的養(yǎng)殖壓力。適宜性評價關鍵在于環(huán)境因素評分的設定,已有的評分方法多采用經驗打分和等間距打分(Radiarta, 2008),使得評分結果主觀性較大。Gentry等(2017)利用生理、異速生長和生長理論的方法,以多物種GPI值(Growth performance index)作為一個養(yǎng)殖生長適宜性的評價標準,GPI值由von Bertalanffy生長方程(Froehlich, 2016)和具體養(yǎng)殖品種參數(shù)求得,高的GPI值視為有更好的生長條件,對水產養(yǎng)殖潛力進行了宏觀評估。本研究針對桑溝灣及周圍海域具體的養(yǎng)殖品種進行適宜性評價,采用影響海帶生長的經驗公式作為評分曲線,能較好地反映出海帶生長對各環(huán)境因素的需求范圍,且通過評分曲線,各環(huán)境因素實測值都有與之對應的評分,評分結果較為科學。

    海帶養(yǎng)殖對季節(jié)依賴性較強,4個季節(jié)適宜性有較大的差異。春季高分區(qū)呈塊狀分布可能由于春季海帶生長較快,海區(qū)無機氮濃度處于一年當中最低值,因此,高分值區(qū)分散在無機氮較為充足的區(qū)域。夏季最適宜面積大于春季、秋季和冬季,可能由于夏季降水多,陸源輸送使得海區(qū)無機氮得到補充,且此時部分海帶已收割,海帶自身的阻流作用減弱,水交換較好,處于比較適宜的環(huán)境狀態(tài)。秋季適宜性評分較高的區(qū)域主要分布在桑溝灣灣外、愛蓮灣和楮島東部海域附近,此時,海表面光強比較平均,相對于灣內,灣外水交換較快,無機氮可得到及時補充,且水溫比較適宜,所以更適宜海帶生長。冬季適宜性評分較高的區(qū)域逐漸向近岸處遷移,數(shù)據模擬顯示(Xuan, 2016),由于冬季水溫較低,靠近灣內水溫相對較高,適宜海帶生長,且冬季灣外風浪大,風浪攪動水底泥沙,實測數(shù)據顯示,冬季水體透明度降低,影響海帶受光,不適宜海帶生長(平仲良, 1993)。

    適宜性評價中,采用海帶生長DEB模型——STELLA模型進行敏感性分析,結果顯示,鹽度的敏感度大于溫度和無機氮(蔡碧瑩等, 2019)。雖然,在養(yǎng)殖環(huán)境中,水體鹽度發(fā)生改變會對海帶生長產生較大的影響,但桑溝灣沒有大型河流匯入,多年平均降雨量處于中等水平,鹽度的年變化較小,基本處于適宜海帶生長的鹽度范圍,且桑溝灣及周圍海域鹽度強制函數(shù)()=0.91~0.96 (蔡碧瑩等, 2019),相較溫度和無機氮,鹽度對海帶生長限制較小。敏感性分析雖未涉及流速參數(shù),但由于水動力影響灣內外水交換和無機氮的輸送,水交換帶動營養(yǎng)鹽補充(史潔等, 2010),所以,流速較無機氮對海帶生長影響更大,此外,溫度影響海水流動(馮士筰等, 1999),因而影響海帶生長的環(huán)境因子重要性依次為光照>溫度>流速>無機氮>鹽度。

    根據實際生產結果,海帶養(yǎng)殖適宜性評價結果得到了驗證,證明該評價結果較為科學。在北部愛蓮灣海域,淺水區(qū)水溫較高,水溫升高導致海帶腐爛,海帶長度最多長到3 m;位于最適宜海區(qū)海帶長度可達4 m左右。楮島近岸海域周圍(水深<5 m的紅色區(qū)域),由于受到水深限制,海帶長度<3m,達不到一般正常尺寸(4 m),不適宜安排海帶養(yǎng)殖。對于灣外水深超過30 m的區(qū)域,雖然,從適宜性評價結果上屬于最適宜或中等適宜,但實際海帶養(yǎng)殖方面,由于風浪較大,養(yǎng)殖設施損失較大,從經濟利益角度考慮并不適宜安排海帶養(yǎng)殖。楮島東部海域附近有大面積的海域屬于最適宜養(yǎng)殖,但該區(qū)域養(yǎng)殖筏架主要集中在近岸,離岸較近海域養(yǎng)殖易受沿岸人類活動影響,因此,可適當將海帶養(yǎng)殖筏架后移。

    4 結論

    本研究通過適宜性評價對桑溝灣及周圍海域海帶養(yǎng)殖的適宜性進行了評分,劃分出不同的適宜性等級。結果表明,適宜性評分較高的區(qū)域位于愛蓮灣及楮島東部海域周圍,評分較低的區(qū)域分布在近岸水深較淺的區(qū)域(水深<5 m)。評價結果顯示了桑溝灣及周圍海域海帶養(yǎng)殖區(qū)適宜性分布,可根據適宜性評價結果調整現(xiàn)有的養(yǎng)殖布局,將海帶養(yǎng)殖區(qū)優(yōu)先安排在海帶生長最適宜區(qū)或適宜性評分較高的海區(qū)。

    GIS作為有效的適宜性評價工具,可將海域時間和空間上非連續(xù)的環(huán)境參數(shù),可視化的反映其空間分布特征,并簡化數(shù)據計算和多圖層疊加過程,增加了評價結果的精確性和科學性。對于一片養(yǎng)殖水域進行適宜性評價,需要綜合考慮養(yǎng)殖過程涉及到的各個方面,例如養(yǎng)殖品種的生長狀況、養(yǎng)殖活動對生態(tài)環(huán)境的影響、養(yǎng)殖產量及利潤的評估,養(yǎng)殖風險的預測分析等。這些模型的結合使水產養(yǎng)殖適宜性綜合評價更為全面,所呈現(xiàn)出的評價結果更科學合理,可對養(yǎng)殖場選址提出建議,為決策者和利益相關者提供參考,有助于我國開展基于生態(tài)系統(tǒng)的水產養(yǎng)殖空間規(guī)劃以及可持續(xù)的水產養(yǎng)殖管理工作。

    Barr NG, Kloeppel A, Rees TAV,. Wave surge increases rates of growth and nutrient uptake in the green seaweedmaintained at low bulk flow velocities. Aquatic Biology, 2008, 3(2): 179–186

    Cai BY, Zhu CB, Liu H,. Model simulated growth of kelpin Sanggou Bay. Progress in Fishery Sciences, 2019, 40(3): 31–41 [蔡碧瑩, 朱長波, 劉慧, 等. 桑溝灣海帶生長的模型預測. 漁業(yè)科學進展, 2019,40(3): 31–41]

    Chen GL, Wang DS. Experience in pilot production management of kelp farming. China Fisheries, 1958(4): 10 [陳根祿, 王東室. 海帶養(yǎng)殖試點生產管理中的幾點體會. 中國水產, 1958(4): 10]

    Chen DY, Wang JX. Observation of the relationship between the growth and water temperature of kelp in southern Zhejiang. Zhejiang Agricultural Science, 1964(2): 89–93 [陳達義, 汪進興. 海帶在浙南沿海生長發(fā)育與水溫關系的觀察. 浙江農業(yè)科學, 1964(2): 89–93]

    Duarte P, Meneses R, Hawkins AJS,. Mathematical modelling to assess the carrying capacity for multi-species culture within coastal waters. Ecological Modelling, 2003, 168(1–2): 109–143

    EPA (Environmental Protection Agency, USA). Rates, constants, and kinetics. Formulations in surface water quality modeling: 2nd Edition, 1985, 455

    Fang JG, Sun HL, Kuang SH,. Assessing the carrying capacity of Sanggou Bay for culture of kelp. Marine Fisheries Research, 1996(2): 7–17 [方建光, 孫慧玲, 匡世煥, 等. 桑溝灣海帶養(yǎng)殖容量的研究. 海洋水產研究, 1996(2): 7–17]

    FAO (Food and Agriculture Organization of the United Nations). Ninth session FAO COFI sub-committee on aquaculture, Rome, 2017, http://www.fao.org/fishery/aquaculture/en

    Feng SZ, Li FQ, Li SJ. Introduction to ocean science. Beijing: High Education Press, 1999 [馮士筰, 李鳳歧, 李少菁. 海洋科學導論. 北京: 高等教育出版社, 1999]

    Foresman TW. The history of geographic information system perspectives from the pioneers. Englewood: Prentice Hall PTR, 1998

    Froehlich HE, Gentry RR, Halpern BS. Synthesis and comparative analysis of physiological tolerance and life-history growth traits in marine aquaculture species. Aquaculture, 2016, 460: 75–82

    Gentry RR, Halley EF, Grimm D,. Mapping the global potential for marine aquaculture. Nature Ecology and Evolution, 2017, 1(9): 1317–1324

    Harlin MM. Nitrate uptake byspp. (Chlorophyceae): Application to aquaculture systems. Aquaculture, 1978, 15(4): 373–376

    Kitadai Y, Kadowaki S. The growth process and N, P uptake rates ofcultured in coastal fish farms. Aquaculture Science, 2003, 51(1): 15–23

    Leigh EG, Paine RT, Quinn JF,. Wave energy and intertidal productivity. Proceedings of the National Academy of Sciences of the United States of America, 1987, 84(5): 1314–1318

    Liu H, Cai BY. Advance in research and application on aquaculture carrying capacity. Progress in Fishery Sciences, 2018, 39(3): 158–166 [劉慧, 蔡碧瑩. 水產養(yǎng)殖容量研究進展及應用. 漁業(yè)科學進展, 2018, 39(3): 158–166]

    Lobban CS, Harrison PJ. Seaweed ecology and physiology. Cambridge: Cambridge University Press, 1996

    Lobban CS, Harrison PJ. Seaweed ecology and physiology. Cambridge: Cambridge University Press, 1994, 123–162

    Malczewski J. On the use of weighted linear combination method in GIS: Common and best practice approach. Transaction in GIS, 2000, 4(1): 5–22

    Martins I, Marques JC. A model for the growth of opportunistic macroalgae (sp) in tidal estuaries. Estuarine, Coastal and Shelf Science, 2002, 55(2): 247–257

    Ottinger M, Clauss K, Kuenzer C. Aquaculture: Relevance, distribution, impacts and spatial assessments – A review. Ocean and Coastal Management, 2016, 119: 244–266

    Ping ZL. Calculating the suspended matter content in the Yellow Sea with the seawater transparency data and NOAA satellite data. Oceanologia et Limnologia Sinica, 1993, 24(1): 24–30, 118 [平仲良. 用實測海水透明度數(shù)據和NOAA衛(wèi)星數(shù)據計算黃海懸浮體含量. 海洋與湖沼, 1993, 24(1): 24–30, 118]

    Radiarta IN, Saitoh SI, Miyazono A. GIS-based multi-criteria evaluation models for identifying suitable sites for Japanese scallop () aquaculture in Funka Bay, southwestern Hokkaido, Japan. Aquaculture, 2008, 284(1–4): 127–135

    Saaty TL. A scaling method for priorities in hierarchical structures. Journal of Mathematical Psychology, 1977, 15(3): 234–281

    Shi J, Wei H, Zhao L,. Study on ecosystem model of multi-species culture in Sanggou Bay: Ⅱ Simulation of ecosystem and the circulation of nutrients. Progress in Fishery Sciences, 2010, 31(4): 36–42 [史潔, 魏皓, 趙亮, 等. 桑溝灣多元養(yǎng)殖生態(tài)模型研究:Ⅱ生態(tài)環(huán)境模擬與生源要素循環(huán). 漁業(yè)科學進展, 2010, 31(4): 36–42]

    Shi J. Numerical study on the influences of physical processes on the aquaculture carrying capacity in a semi-enclosed bay. Doctoral Dissertation of Ocean University of China, 2009 [史潔. 物理過程對半封閉海灣養(yǎng)殖容量影響的數(shù)值研究. 中國海洋大學博士研究生學位論文, 2009]

    Steele JH. Environmental control of photosynthesis in the sea. Limnology and Oceanography, 1962, 7(2): 137–150

    Suzuki S, Furuya K, Kawai T,. Effect of seawater temperature on the productivity ofin the Uwa Sea, southern Japan. Journal of Applied Phycology, 2008, 20(5): 833–844

    Wang X, Li WQ. Effects of salinity and pH on photosynthesis rate of marine algae. Marine Environmental Science, 1991, 10(1): 37–40 [王憲, 李文權. 鹽度、pH對海洋藻類光合作用速率的影響. 海洋環(huán)境科學, 1991, 10(1): 37–40]

    Wheeler WN. Effect of boundary layer transport on the fixation of carbon by the giant kelp. Marine Biology, 1980, 56(2): 103–110

    Wu RJ, Zhang XL, Zhu MY,. A model for the growth of Haidai () in aquaculture. Marine Science Bullentin, 2009, 28(2): 34–40 [吳榮軍, 張學雷, 朱明遠, 等. 養(yǎng)殖海帶的生長模型研究. 海洋通報, 2009, 28(2): 34–40]

    Xuan JL, Yang ZQ, Huang DJ,. Tidal residual current and its role in the mean flow on the Changjiang bank. Journal of Marine Systems 2016, 154: 66–81

    Zhang DM, Liao GR, Yang QM. Research on the relationship between the coastal current andraising Ⅱ. An experiment of comparsion between the current speed andgrowth. Journal of Shandong College of Oceanology, 1982, 12(3): 73–79 [張定民, 繆國榮, 楊清明. 沿岸流與海帶養(yǎng)殖關系的研究Ⅱ. 流速對海帶生長的影響. 山東海洋學院學報, 1982, 12(3): 73–79]

    Zhang ZH, Lv JB, Ye SF,. Values of marine ecosystem services in Sanggou Bay. Chinese Journal of Applied Ecology, 2007, 18(11): 2540?2547 [張朝暉, 呂吉斌, 葉屬峰, 等. 桑溝灣海洋生態(tài)系統(tǒng)的服務價值. 應用生態(tài)學報, 2007, 18(11): 2540–2547]

    GIS-Based Assessment for Culture Suitability of Kelp () in Sanggou Bay and the Surrounding Waters

    SUN Qianwen1,2, LIU Hui2①, SHANG Weitao3,4, YU Liangju4, JIANG Xiaopeng3,4, CAI Biying1,2, CHANG Lirong5, XIAO Luyang5

    (1.201306; 2.266071; 3.100049; 4.264003; 5.264316)

    With the expansion of aquaculture, there is an increasing concern regarding issues such as conflicts with multi-sea use, low productivity, and higher economic costs. The selection of suitable sites for aquaculture is essential for utilization of sea space and improvement of the efficiency and profits of aquaculture industry. Kelp () is the main culture species in Sanggou Bay and the surrounding waters. The assessment of suitability for kelp culture has been conducted with factors, including light intensity, water temperature, velocity, dissolved inorganic nitrogen, salinity, and depth. The remote sensing is employed to obtain the information of the existing culture areas. The environment variables were collected from field surveys along with model simulations. The Dynamic Energy Budget (DEB) model–STELLA model coupled with analytic hierarchy process was applied to calculate weighted-sums of scores for suitability evaluation. The result showed that the suitability scores were 0~6.7 with the main distribution range of 4~6. 23% of area was estimated to be the most suitable for kelp culture, whereas 10% of total study area was not suitable. The high score areas are mainly located in Ailian Bay and Chudao eastern surrounding waters, whereas the low score areas are mainly located in the shallow waters near the shore. The results have provided a guideline to help in site selection and in the management of kelp culture.

    Aquaculture; Suitability assessment; Site selection; GIS

    S931.3

    A

    2095-9869(2020)01-0041-10

    10.19663/j.issn2095-9869.20181024001

    * 科技部國際創(chuàng)新合作專項“基于生態(tài)系統(tǒng)的水產養(yǎng)殖空間規(guī)劃研究”(2016YFE0112600)和歐盟地平線2020項目(633476-H2020-SFS-2014-2015)共同資助[This work was supported by the Key Programme for International Cooperation on Scientific and Technological Innovation, Ministry of Science and Technology (2016YFE0112600), and Optimizing Space Available for European Aquaculture (AquaSpace) (633476-H2020-SFS-2014-2015)]. 孫倩雯, E-mail: gwen_sun@126.com

    劉 慧,研究員,E-mail: liuhui@ysfri.ac.cn

    2018-10-24,

    2018-11-15

    http://www.yykxjz.cn/

    孫倩雯, 劉慧, 尚偉濤, 于良巨, 姜曉鵬, 蔡碧瑩, 常麗榮, 肖露陽. 基于GIS的桑溝灣及周圍海域海帶養(yǎng)殖適宜性評價. 漁業(yè)科學進展, 2020, 41(1): 41–50

    Sun QW, Liu H, Shang WT, Yu LJ, Jiang XP, Cai BY, Chang LR, Xiao LY. GIS-Based assessment for culture suitability of kelp () in Sanggou Bay and the surrounding waters. Progress in Fishery Sciences, 2020, 41(1): 41–50

    LIU Hui, E-mail: liuhui@ysfri.ac.cn

    (編輯 馬璀艷)

    猜你喜歡
    海帶鹽度圖層
    即興
    飛碟探索(2023年3期)2023-11-23 16:08:36
    如果沒有你
    飛碟探索(2023年1期)2023-05-30 22:59:34
    夢的解析
    第一財經(2020年10期)2020-10-15 03:17:55
    第一財經(2019年7期)2019-07-25 17:48:18
    巧用混合圖層 制作抽象動感森林
    圖層法在地理區(qū)域圖讀圖中的應用
    地理教學(2015年18期)2016-01-06 12:00:38
    鹽度和pH對細角螺耗氧率和排氨率的影響
    鹽度脅迫對入侵生物福壽螺的急性毒性效應
    適用于高鹽度和致密巖層驅油的表面活性劑
    淺析“遞層優(yōu)化法”在礦井制圖中的應用
    河南科技(2014年6期)2014-04-04 08:00:42
    伊人久久精品亚洲午夜| 一区二区三区免费毛片| 十八禁国产超污无遮挡网站| 色综合站精品国产| 亚洲成色77777| 男女那种视频在线观看| 人妻少妇偷人精品九色| 搡老乐熟女国产| 2022亚洲国产成人精品| 国产在线一区二区三区精| 男人和女人高潮做爰伦理| 爱豆传媒免费全集在线观看| 人人妻人人澡欧美一区二区| 久久亚洲国产成人精品v| 国产精品爽爽va在线观看网站| 街头女战士在线观看网站| 午夜精品一区二区三区免费看| 一区二区三区免费毛片| 午夜激情久久久久久久| 欧美日韩在线观看h| 国产成人freesex在线| 高清毛片免费看| 欧美zozozo另类| 赤兔流量卡办理| 大片免费播放器 马上看| 美女黄网站色视频| 中文字幕av成人在线电影| 久久久久免费精品人妻一区二区| 女人被狂操c到高潮| 一个人免费在线观看电影| 精品不卡国产一区二区三区| 欧美激情国产日韩精品一区| 国产高清三级在线| av国产免费在线观看| 在线a可以看的网站| 国产在线男女| 欧美高清性xxxxhd video| 麻豆精品久久久久久蜜桃| 五月天丁香电影| 久久精品综合一区二区三区| 最新中文字幕久久久久| 久久久久久久亚洲中文字幕| 国产淫片久久久久久久久| 一二三四中文在线观看免费高清| 久久久精品94久久精品| 久久久久久久久久黄片| 婷婷色麻豆天堂久久| 国产一区二区三区综合在线观看 | 日韩 亚洲 欧美在线| 综合色av麻豆| 国产成人福利小说| 99久久精品热视频| 自拍偷自拍亚洲精品老妇| 日韩av不卡免费在线播放| 欧美日本视频| 国产单亲对白刺激| 亚洲最大成人av| 高清视频免费观看一区二区 | 青春草国产在线视频| 精品久久国产蜜桃| 91aial.com中文字幕在线观看| 亚洲精品日韩av片在线观看| 亚洲精品亚洲一区二区| 一区二区三区高清视频在线| 国产午夜精品论理片| 欧美丝袜亚洲另类| 亚洲精品久久午夜乱码| 久久精品夜夜夜夜夜久久蜜豆| 日韩av在线免费看完整版不卡| 白带黄色成豆腐渣| 99re6热这里在线精品视频| 一级毛片aaaaaa免费看小| 成人亚洲精品av一区二区| 三级毛片av免费| 国产精品久久久久久av不卡| 91久久精品国产一区二区三区| 亚洲av成人av| 成人毛片a级毛片在线播放| 成年女人看的毛片在线观看| 国产黄片美女视频| 亚洲美女搞黄在线观看| 搡老妇女老女人老熟妇| 国产黄a三级三级三级人| 久久久久久久午夜电影| freevideosex欧美| 天天躁夜夜躁狠狠久久av| 波野结衣二区三区在线| 国产在视频线在精品| 三级国产精品片| 精品欧美国产一区二区三| 国产精品久久久久久精品电影| 成人毛片a级毛片在线播放| 一级毛片aaaaaa免费看小| 亚洲精品456在线播放app| 中文乱码字字幕精品一区二区三区 | 国产精品1区2区在线观看.| 欧美成人午夜免费资源| h日本视频在线播放| 国产午夜福利久久久久久| 联通29元200g的流量卡| 色哟哟·www| 欧美bdsm另类| 精品一区二区三区人妻视频| 人妻夜夜爽99麻豆av| 精品久久久久久成人av| 99热全是精品| 夫妻午夜视频| 成人漫画全彩无遮挡| 国产精品一及| 三级男女做爰猛烈吃奶摸视频| 99视频精品全部免费 在线| 亚洲经典国产精华液单| 中国美白少妇内射xxxbb| 国产探花在线观看一区二区| 亚洲av在线观看美女高潮| 国产综合懂色| 中文字幕久久专区| 国产精品麻豆人妻色哟哟久久 | 午夜视频国产福利| 日韩一区二区视频免费看| 国产亚洲精品av在线| 九草在线视频观看| 亚洲人成网站在线观看播放| av在线老鸭窝| 中文在线观看免费www的网站| 在线观看人妻少妇| 国产欧美日韩精品一区二区| 久久久欧美国产精品| 久久精品国产亚洲av涩爱| 精品国产露脸久久av麻豆 | 网址你懂的国产日韩在线| 国产精品.久久久| 夫妻性生交免费视频一级片| 成人综合一区亚洲| 天堂网av新在线| 国产激情偷乱视频一区二区| 国产成人午夜福利电影在线观看| 国产淫语在线视频| 大又大粗又爽又黄少妇毛片口| 国产精品麻豆人妻色哟哟久久 | 91狼人影院| 国产精品不卡视频一区二区| 男女下面进入的视频免费午夜| 亚洲成人中文字幕在线播放| 国产不卡一卡二| 久久精品国产自在天天线| 中文欧美无线码| 2022亚洲国产成人精品| 神马国产精品三级电影在线观看| 三级经典国产精品| 人体艺术视频欧美日本| 九草在线视频观看| 人人妻人人澡人人爽人人夜夜 | av在线老鸭窝| 中文在线观看免费www的网站| 联通29元200g的流量卡| 国产淫片久久久久久久久| 成人av在线播放网站| 国产人妻一区二区三区在| 欧美另类一区| 久久精品综合一区二区三区| 床上黄色一级片| 亚洲精品国产成人久久av| 久久久精品欧美日韩精品| 国产亚洲午夜精品一区二区久久 | 日韩欧美 国产精品| 亚洲欧美一区二区三区黑人 | 亚洲国产精品sss在线观看| 国产欧美另类精品又又久久亚洲欧美| 1000部很黄的大片| 国产高清不卡午夜福利| 国产av码专区亚洲av| 国产老妇女一区| 两个人的视频大全免费| 精品人妻偷拍中文字幕| 只有这里有精品99| 神马国产精品三级电影在线观看| 亚洲乱码一区二区免费版| 国产精品综合久久久久久久免费| av在线观看视频网站免费| 色综合亚洲欧美另类图片| www.色视频.com| 国产精品av视频在线免费观看| 舔av片在线| 一个人观看的视频www高清免费观看| 老师上课跳d突然被开到最大视频| 日本三级黄在线观看| 欧美变态另类bdsm刘玥| 美女脱内裤让男人舔精品视频| 日韩人妻高清精品专区| 人妻夜夜爽99麻豆av| 免费看a级黄色片| 国产三级在线视频| 国产精品福利在线免费观看| 久久精品国产自在天天线| 日韩欧美三级三区| 国产精品爽爽va在线观看网站| 一个人免费在线观看电影| 国模一区二区三区四区视频| 日韩亚洲欧美综合| 69av精品久久久久久| 舔av片在线| 嫩草影院入口| 国模一区二区三区四区视频| 午夜精品在线福利| 97精品久久久久久久久久精品| 蜜桃亚洲精品一区二区三区| 亚洲av在线观看美女高潮| 国产成人精品一,二区| 美女被艹到高潮喷水动态| 91精品一卡2卡3卡4卡| 成人无遮挡网站| 六月丁香七月| 亚洲欧美中文字幕日韩二区| 秋霞在线观看毛片| 99热全是精品| 晚上一个人看的免费电影| 午夜久久久久精精品| 日韩欧美一区视频在线观看 | 国产精品人妻久久久影院| 国产欧美另类精品又又久久亚洲欧美| 在线免费观看的www视频| 搡老乐熟女国产| av一本久久久久| 乱人视频在线观看| 久久午夜福利片| 97精品久久久久久久久久精品| 色尼玛亚洲综合影院| 国产成人午夜福利电影在线观看| 免费av不卡在线播放| 国产永久视频网站| 黄片无遮挡物在线观看| 精品国产三级普通话版| 国产成年人精品一区二区| 狂野欧美激情性xxxx在线观看| 69av精品久久久久久| 亚洲最大成人手机在线| 麻豆乱淫一区二区| 久久久久久久国产电影| 老女人水多毛片| 午夜免费男女啪啪视频观看| 精品不卡国产一区二区三区| 国产午夜精品久久久久久一区二区三区| 最新中文字幕久久久久| 国产精品三级大全| 精品久久久久久久末码| 久久久久九九精品影院| 国产精品嫩草影院av在线观看| 啦啦啦啦在线视频资源| 在线观看人妻少妇| 亚洲精品一二三| 十八禁网站网址无遮挡 | 高清午夜精品一区二区三区| 国产成人精品福利久久| av.在线天堂| 看免费成人av毛片| 久久久久久久久久成人| 少妇高潮的动态图| 久久精品国产亚洲av涩爱| 欧美日韩综合久久久久久| 在线观看一区二区三区| 国产成人福利小说| 性色avwww在线观看| 免费在线观看成人毛片| 丝袜喷水一区| 波野结衣二区三区在线| 黄片wwwwww| 国产综合精华液| 水蜜桃什么品种好| 日本一本二区三区精品| 夜夜看夜夜爽夜夜摸| 99久国产av精品| 晚上一个人看的免费电影| 成人综合一区亚洲| 一个人看的www免费观看视频| 蜜臀久久99精品久久宅男| 国产精品人妻久久久影院| 午夜福利在线观看吧| 欧美 日韩 精品 国产| 91狼人影院| 午夜福利成人在线免费观看| kizo精华| 午夜福利网站1000一区二区三区| 天堂√8在线中文| 人妻少妇偷人精品九色| 国产激情偷乱视频一区二区| 亚洲精品乱码久久久久久按摩| 国产男人的电影天堂91| 国产精品熟女久久久久浪| 夜夜爽夜夜爽视频| 国产av码专区亚洲av| 亚洲国产av新网站| 插阴视频在线观看视频| 亚洲久久久久久中文字幕| 搞女人的毛片| 亚洲精品aⅴ在线观看| 一级片'在线观看视频| 久久久色成人| 青春草视频在线免费观看| 午夜福利高清视频| 色哟哟·www| 99热6这里只有精品| 欧美日本视频| 亚洲成人av在线免费| 三级男女做爰猛烈吃奶摸视频| 少妇熟女aⅴ在线视频| 亚洲av不卡在线观看| 亚洲aⅴ乱码一区二区在线播放| freevideosex欧美| 日韩欧美一区视频在线观看 | videos熟女内射| 少妇高潮的动态图| 国产极品天堂在线| 综合色av麻豆| 天堂俺去俺来也www色官网 | 国产精品av视频在线免费观看| 午夜福利在线观看吧| 亚洲久久久久久中文字幕| 国产精品一区二区在线观看99 | 天美传媒精品一区二区| 水蜜桃什么品种好| 91狼人影院| 国产 一区精品| 女人十人毛片免费观看3o分钟| 中文字幕制服av| 欧美区成人在线视频| 久久久国产一区二区| 在线免费十八禁| 欧美变态另类bdsm刘玥| 日韩av在线免费看完整版不卡| 一级毛片电影观看| 建设人人有责人人尽责人人享有的 | 如何舔出高潮| 成年版毛片免费区| 亚洲自偷自拍三级| 99热6这里只有精品| 欧美激情国产日韩精品一区| 亚洲美女视频黄频| 最近的中文字幕免费完整| 国产伦在线观看视频一区| 欧美日韩在线观看h| 天天躁日日操中文字幕| 国产色爽女视频免费观看| 日韩人妻高清精品专区| 18禁在线播放成人免费| 精品亚洲乱码少妇综合久久| 国产精品一区二区性色av| 国产精品精品国产色婷婷| 亚洲精品中文字幕在线视频 | 成人美女网站在线观看视频| 成人毛片a级毛片在线播放| 18禁在线播放成人免费| 精品亚洲乱码少妇综合久久| 精品国产露脸久久av麻豆 | 国产精品熟女久久久久浪| 80岁老熟妇乱子伦牲交| av国产免费在线观看| 美女内射精品一级片tv| 国产综合懂色| 中文天堂在线官网| 深爱激情五月婷婷| 午夜福利成人在线免费观看| 老司机影院毛片| 日本猛色少妇xxxxx猛交久久| 久久精品久久久久久久性| 人人妻人人澡人人爽人人夜夜 | 日本一二三区视频观看| 亚洲精品一二三| 99久久精品热视频| 乱系列少妇在线播放| 国产精品不卡视频一区二区| 国产三级在线视频| 最后的刺客免费高清国语| 欧美一级a爱片免费观看看| 永久免费av网站大全| 国产精品女同一区二区软件| 精品久久久久久成人av| 国产女主播在线喷水免费视频网站 | 欧美日韩视频高清一区二区三区二| 欧美成人精品欧美一级黄| 纵有疾风起免费观看全集完整版 | 2021少妇久久久久久久久久久| 亚洲在线观看片| 婷婷色麻豆天堂久久| 国产免费视频播放在线视频 | 亚洲精品乱久久久久久| 最近2019中文字幕mv第一页| 国内精品一区二区在线观看| 2018国产大陆天天弄谢| 国产精品伦人一区二区| 欧美区成人在线视频| 国产亚洲午夜精品一区二区久久 | 成人亚洲精品av一区二区| 日本-黄色视频高清免费观看| 久久精品综合一区二区三区| 91精品国产九色| 男人爽女人下面视频在线观看| 看黄色毛片网站| 亚洲av男天堂| 日韩亚洲欧美综合| 国产午夜精品论理片| 国产精品一及| 国产精品国产三级国产专区5o| av天堂中文字幕网| 亚洲精品国产成人久久av| 精品久久久久久成人av| 中文天堂在线官网| 欧美成人a在线观看| 波多野结衣巨乳人妻| 日韩精品有码人妻一区| 五月天丁香电影| 欧美激情久久久久久爽电影| 中文精品一卡2卡3卡4更新| 男人狂女人下面高潮的视频| 夜夜看夜夜爽夜夜摸| 久久综合国产亚洲精品| 男女那种视频在线观看| 午夜日本视频在线| 18禁在线无遮挡免费观看视频| 日韩在线高清观看一区二区三区| 在线观看一区二区三区| 亚洲欧美一区二区三区国产| 国产白丝娇喘喷水9色精品| 国产黄a三级三级三级人| 国产久久久一区二区三区| av播播在线观看一区| 2021少妇久久久久久久久久久| 国内精品一区二区在线观看| 高清av免费在线| 亚洲av不卡在线观看| 国内少妇人妻偷人精品xxx网站| 嫩草影院新地址| 99久久精品一区二区三区| 国产伦精品一区二区三区视频9| 国产美女午夜福利| 国产av码专区亚洲av| 国产精品女同一区二区软件| 国产免费又黄又爽又色| 草草在线视频免费看| 99久久精品国产国产毛片| 日韩,欧美,国产一区二区三区| 久久久精品欧美日韩精品| 国产探花在线观看一区二区| 国产一区二区三区av在线| 亚洲最大成人手机在线| 亚洲欧美日韩东京热| 床上黄色一级片| av线在线观看网站| 午夜福利在线在线| 国产亚洲av嫩草精品影院| 99久久九九国产精品国产免费| 欧美高清性xxxxhd video| a级一级毛片免费在线观看| 精品久久久久久电影网| 女人被狂操c到高潮| 男女边摸边吃奶| 国内精品宾馆在线| 人妻制服诱惑在线中文字幕| 在线免费观看的www视频| 午夜福利成人在线免费观看| 日日啪夜夜撸| 精品一区二区三区视频在线| 久久久久久久久久人人人人人人| 搡老妇女老女人老熟妇| 亚洲av中文字字幕乱码综合| 国产精品国产三级国产av玫瑰| 狠狠精品人妻久久久久久综合| 久久99精品国语久久久| 亚洲四区av| 色网站视频免费| 久久精品久久精品一区二区三区| 午夜免费男女啪啪视频观看| 高清av免费在线| 日韩电影二区| 久久久久久久久久黄片| 欧美高清性xxxxhd video| 街头女战士在线观看网站| 国产精品美女特级片免费视频播放器| 夜夜看夜夜爽夜夜摸| 亚洲aⅴ乱码一区二区在线播放| 午夜免费男女啪啪视频观看| 亚洲av成人av| 九九爱精品视频在线观看| av免费观看日本| 国产精品一二三区在线看| 91久久精品电影网| 午夜福利视频1000在线观看| 99re6热这里在线精品视频| 有码 亚洲区| 天堂网av新在线| 免费观看a级毛片全部| 久久精品久久久久久久性| 亚洲精品自拍成人| 国产亚洲精品久久久com| 夫妻性生交免费视频一级片| 一二三四中文在线观看免费高清| 特级一级黄色大片| 18禁在线播放成人免费| 小蜜桃在线观看免费完整版高清| 国产男人的电影天堂91| 亚洲人成网站高清观看| 国产成人91sexporn| 欧美成人午夜免费资源| 街头女战士在线观看网站| 秋霞在线观看毛片| av网站免费在线观看视频 | 色视频www国产| 国产免费一级a男人的天堂| 一区二区三区免费毛片| 国产成人精品福利久久| 国产综合精华液| 欧美高清成人免费视频www| 又粗又硬又长又爽又黄的视频| 中文字幕人妻熟人妻熟丝袜美| 久久99热这里只有精品18| 超碰97精品在线观看| 亚洲精品视频女| 久久国内精品自在自线图片| 国产精品1区2区在线观看.| 天天躁日日操中文字幕| 欧美不卡视频在线免费观看| 校园人妻丝袜中文字幕| 亚洲欧美成人精品一区二区| 色综合站精品国产| 一级毛片aaaaaa免费看小| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲最大成人av| 可以在线观看毛片的网站| 亚洲av福利一区| 成人美女网站在线观看视频| 国产精品一及| 乱码一卡2卡4卡精品| 久久久久久久久大av| 中文字幕亚洲精品专区| 亚洲精品日韩av片在线观看| 青春草国产在线视频| 国产永久视频网站| 国产精品久久久久久精品电影小说 | av免费在线看不卡| 色尼玛亚洲综合影院| 成人国产麻豆网| 色尼玛亚洲综合影院| 一级毛片 在线播放| 免费黄频网站在线观看国产| 亚洲精品乱久久久久久| 最近最新中文字幕大全电影3| 亚洲18禁久久av| 可以在线观看毛片的网站| 成人特级av手机在线观看| 午夜福利网站1000一区二区三区| 一区二区三区免费毛片| 美女脱内裤让男人舔精品视频| 插阴视频在线观看视频| 美女脱内裤让男人舔精品视频| 97热精品久久久久久| .国产精品久久| 国模一区二区三区四区视频| 九九在线视频观看精品| 日本av手机在线免费观看| 伦精品一区二区三区| 偷拍熟女少妇极品色| 一个人看视频在线观看www免费| 69人妻影院| 国产毛片a区久久久久| 毛片一级片免费看久久久久| 成年版毛片免费区| 国产精品久久久久久av不卡| 波野结衣二区三区在线| 一本久久精品| 国产免费福利视频在线观看| 在现免费观看毛片| 国产精品综合久久久久久久免费| 国产人妻一区二区三区在| 日本-黄色视频高清免费观看| 好男人在线观看高清免费视频| 亚洲精品aⅴ在线观看| 国产精品久久久久久精品电影小说 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费观看的影片在线观看| 人妻夜夜爽99麻豆av| 亚洲欧美一区二区三区国产| 亚洲精品,欧美精品| 99久国产av精品| 菩萨蛮人人尽说江南好唐韦庄| 又黄又爽又刺激的免费视频.| 亚洲欧洲日产国产| 熟妇人妻久久中文字幕3abv| 亚洲精品久久午夜乱码| 网址你懂的国产日韩在线| 日本一二三区视频观看| 一区二区三区高清视频在线| 99热全是精品| 国产一区二区在线观看日韩| 老女人水多毛片| 精品久久国产蜜桃| 色5月婷婷丁香| 日本爱情动作片www.在线观看| 精品酒店卫生间| 亚洲欧美日韩卡通动漫| 亚洲av二区三区四区| 午夜福利视频1000在线观看| 久久久久久九九精品二区国产| a级毛色黄片| 日韩伦理黄色片| 最近最新中文字幕大全电影3| 亚洲欧洲日产国产| 欧美人与善性xxx| 在线观看免费高清a一片| 国产熟女欧美一区二区| 国产片特级美女逼逼视频| 免费高清在线观看视频在线观看| 亚洲av不卡在线观看| 国产精品国产三级国产av玫瑰| av专区在线播放| 欧美一区二区亚洲|