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

    釣鉤深度和浸泡時間對東太平洋公海長鰭金槍魚延繩釣漁獲性能的影響研究?

    2018-12-13 08:00:12劉莉莉虞聰達臧迎亮
    關鍵詞:漁獲金槍魚水層

    劉莉莉, 周 成, 虞聰達, 鄭 基, 臧迎亮

    (1. 浙江海洋大學水產學院, 浙江 舟山 316022; 2. 浙江省海洋漁業(yè)裝備技術研究重點實驗室, 浙江 舟山 316022; 3. 中國海洋大學水產學院, 山東 青島 266003)

    釣鉤深度和浸泡時間對東太平洋公海長鰭金槍魚延繩釣漁獲性能的影響研究?

    劉莉莉1,2, 周 成3**, 虞聰達1,2, 鄭 基1,2, 臧迎亮1,2

    (1. 浙江海洋大學水產學院, 浙江 舟山 316022; 2. 浙江省海洋漁業(yè)裝備技術研究重點實驗室, 浙江 舟山 316022; 3. 中國海洋大學水產學院, 山東 青島 266003)

    為了探究金槍魚延繩釣的釣鉤深度和浸泡時間對釣捕性能和漁獲特征的影響,本文基于2015年9月—2016年2月收集的東太平洋公海長鰭金槍魚漁業(yè)調查數據,利用局部加權回歸和線性回歸分別擬合了釣鉤深度分布與捕獲的長鰭金槍魚叉長分布的關系,統(tǒng)計了各鉤號和叉長組的上鉤漁獲物的存活比例,并利用Logistic回歸模型對釣鉤浸泡時間與捕撈率的關系進行了統(tǒng)計分析。研究表明:漁獲物的平均叉長隨著釣鉤水層的加深而略微增大,192 m(8號鉤所處水層)為臨界水層,該水層以淺的長鰭金槍魚叉長明顯小于該水層以深的叉長。鉤號、叉長均與金槍魚的活體比例呈較明顯的相關性,這表明捕撈水層越深或魚體越大,漁獲存活率越高。以每11個連續(xù)浮球為單位計算浸泡時間,漁獲率隨著浸泡用時呈現波動式增長的動態(tài)響應,并在11和16 h呈現2個峰值。Logistic回歸模型顯示,釣鉤水深對捕撈效率的影響極顯著,167.57 m水深為捕撈效率密度最高的水層,若釣鉤水層分布于124~211 m將有更高的捕撈效率。

    金槍魚延繩釣; 釣鉤深度; 浸泡時間; 捕撈率; 漁獲特征

    長鰭金槍魚(Thunnusalalunga)是高度洄游性的上層營養(yǎng)級海洋捕食者,也是重要的商業(yè)捕撈物種,分布于全球 50°N 和40°S之間除赤道區(qū)域附近以外的溫熱帶海域[1-4]。2015年長鰭金槍魚總產量達到223 013 t,占全球金槍魚和類金槍魚總捕撈量的4.6%,其中,延繩釣捕獲產量占長鰭金槍魚總產量的72%[5](FAO統(tǒng)計數據),近50%的捕撈量來自于太平洋海域。據中西太平洋金槍魚委員會(WCPFC)統(tǒng)計,2015年太平洋長鰭金槍魚捕撈量為102 997 t,我國大陸金槍魚延繩釣捕撈產量約占10%。

    目前,捕撈壓力的持續(xù)增加正導致長鰭金槍魚生物量的減少[6],盡管該種群未處于過度捕撈以及正在過度捕撈的狀態(tài),但為了避免生物量的進一步減少,需要降低延繩釣捕撈死亡率[7]。種群資源評估通常假設捕撈率(CPUE)為種群資源豐度的指數,但由于可捕系數(Catchability)一般不是常數,名義努力量和捕撈產量之間并沒有必然的關系。捕撈率與種群生物量之間的關系也取決于可捕系數的大小。事實上,捕撈率的改變可能是由于漁具結構或者操作所導致的捕撈性能的改變,并非種群資源量發(fā)生變化所致[8-9]。因此,正確理解捕撈率和漁具效率之間的關系,對于基于漁業(yè)數據的種群資源評估和管理極其重要。對于被動且高選擇性的延繩釣捕撈,釣捕率除受時空和海洋環(huán)境的影響外,還與漁具的結構、垂直水層分布和捕撈策略等因素有關,如釣捕深度、時間、餌料及浸泡時間等[10-14],彼此關系仍未被全面認識和應用。

    為此,本研究本著探究金槍魚延繩釣漁具結構和捕撈方式對漁具性能和漁獲特征的內在關系出發(fā),利用2015年9月—2016年2月東太平洋公海長鰭金槍魚漁業(yè)生產船調查數據對不同鉤深下的捕撈率、叉長分布和釣捕死亡比例,以及不同釣鉤浸泡時間與捕撈率之間的關系進行統(tǒng)計分析,旨在探明延繩釣捕撈特性的規(guī)律,以期為提高漁獲效率和漁具選擇性提供重要的科學依據。

    1 材料與方法

    1.1 調查船、時間及海域

    調查作業(yè)漁船為浙江省大洋世家股份有限公司所屬的超低溫金槍魚延繩釣船“新世紀71號”,該船總長56.5 m、型寬8.5 m、型深3.65 m、總噸位634 t、主機功率735 kW。以長鰭金槍魚為主要目標種群,作業(yè)于法屬波利尼西亞塔希提島以東的東太平洋公海。2015年9月—2016年2月期間在14°S~23°S,115°W~124°W水域內共調查49個站點(見圖1)。

    (A(實心圓):從投繩結束處起繩;B(實心三角):從投繩開始處起繩。右上方小圖陰影部分表示該調查區(qū)域。A with filled circle indicates counter-retrieved-operation; B with filled triangle indicates forward-retrieved-operation.(In the zoom-out plot at upper right the survey area is shaded in grey.)

    圖1 延繩釣投放鉤位置

    Fig.1 Distribution of deployment locations in survey of tuna longline

    1.2 漁具結構和捕撈操作

    延繩釣主繩為8股尼龍單絲編織,直徑5.4 mm;支線結構為自動掛扣(帶附環(huán))連接2 m直徑4 mm的紅色聚酯繩(PES),通過一枚箱型轉環(huán)連接18 m直徑1.3 mm尼龍單絲,通過一枚S型轉環(huán)連接6 m直徑1 mm 尼龍單絲,最后懸掛一枚圓型鉤,總長26 m;浮子繩為80股編織的聚酯繩(PES),直徑5.5 mm,長25 m;浮子材料為ABS注塑,直徑360 mm,靜浮力24.5 kg。

    投繩過程中平均船速9.6 kn,投繩平均速度6 m/s,每籃配置26條支線,2支線之間的主繩長度為33 m,每次投放釣鉤3 770枚。投繩起止時間通常為5:00~11:00,持續(xù)6 h;起鉤起止時間通常為14:00~次日5:00,持續(xù)15 h。起鉤方式包括從投繩開始處起繩(起頭)和從投繩結束處起繩(起尾),起鉤平均船速4.5 kn。

    1.3 調查采樣方法

    調查記錄的捕撈技術參數包括總鉤數、兩浮子間的鉤數、出繩速度、船速、兩鉤間的時間間隔、投繩起始船位、投繩終止船位、航向、投繩開始時間、投繩結束時間、起繩開始時間、起繩結束時間和捕撈產量等;使用卷尺(GC-F5019,精度1.0 mm)和臺秤(TGT-500,精度200 g)隨機測定上鉤漁獲物的叉長和體重,并記錄釣捕的鉤號和生存狀態(tài);以每10籃為單位,記錄起鉤期間全部漁獲物的尾數。

    1.4 數據處理

    1.4.1 延繩釣形態(tài)估計 在無環(huán)境(海流)影響下2個浮球之間延繩釣呈懸鏈線狀態(tài),對2個浮球之間懸鏈線對稱一側的釣鉤從淺至深進行編號,假設每2個浮球之間的懸鏈線形狀相同且在操作過程中浮球之間的距離不變,按照Yoshihara[15]提出的延繩釣理論深度計算方法預測各釣鉤的布置水深和與最近浮球間的水平距離:

    Dj=F+B+

    (1)

    (2)

    式中:Dj為釣鉤理論深度(m);Hj為釣鉤與最近浮球之間的水平距離;F為浮子繩長度(m);B為支線長度(m);L為兩浮子間主繩長度;N為每籃支線數,共26枚;θ為延繩釣主繩在浮子繩末端接合點處切線與水平方向的夾角(°),與短縮率k(2個連續(xù)浮球間的水平距離與主繩長度的比值)具有以下關系:

    (3)

    L=(N+1)·S;

    (4)

    (5)

    式中:H為兩相鄰浮子之間的水平距離(m);S為兩相鄰支線間的主繩長度(m);Vs和Vb分別為投繩速度和船速(m/s)。

    實際生產時,由于水流的影響,漁具形狀和短縮率往往發(fā)生改變,釣鉤的實際深度一般小于懸鏈線理論深度。通常采用上浮率(上浮深度占理論深度的百分比)表示釣鉤上浮程度。在不同漁場和時間,上浮率可達11%~46%[10]。例如:太平洋中北部海域金槍魚延繩釣的平均上浮率為35%[16];法屬波利尼西亞專屬經濟區(qū)東北部區(qū)域的平均上浮率為19%[17];夏威夷外海區(qū)域2個不同時期的上浮率分別為46%和32%[10]?;谙嗨坪S虻墓烙嬛担狙芯考僭O平均上浮率19%。

    (6)

    (7)

    各單位內的釣獲率(CPUE,尾/千鉤)計算公式為:

    (8)

    式中:Ci為第i單位內的漁獲數量;Ei為第i單位內的釣鉤數。

    1.4.3 漁獲性能在時間和空間上的變化 考慮到漁獲性能(魚體大小和漁獲率)在空間和時間上的動態(tài)變化可能是非線性的,應用局部加權回歸擬合叉長和鉤號、漁獲率和浸泡時間的關系。局部加權法即取一定比例的局部數據,在這部分子集中擬合多項式回歸曲線,以便觀察到數據在局部展現出來的規(guī)律和趨勢。

    1.4.4 鉤深和捕撈效率的關系 延繩釣的作業(yè)過程可看作對兩浮子間不同水深釣鉤的漁獲進行抽樣,并重復M次試驗,漁獲結果以捕獲(表示為“1”)和未捕獲(表示為“0”)的二值型分類變量表示,屬于0-1響應的二項分布,因此對鉤深與釣捕效率之間的關系采用連接函數為二項分布族中l(wèi)ogit邏輯函數的Logistic回歸模型表達。假設魚群僅分布在釣鉤布置的水層范圍,響應變量為某深度下的捕撈效率,表示為該深度內的累積捕撈效率,如100 m深度的捕撈效率為100 m水深內所有釣鉤的漁獲數量占總漁獲量的比例。通過該模型(表示如下)可預測捕撈成功的概率:

    (9)

    式中:P表示漁獲效率,即某一水深內釣獲的漁獲數量(尾數)占總漁獲數量的比例;D為釣鉤深度(m);β0為截距;β1為模型回歸系數,顯著性水平設為0.05。

    2 結果

    2.1 不同釣鉤的漁獲物特征

    延繩釣相鄰兩浮子的水平距離為712 m,釣鉤的理論計算深度范圍為78~287 m,每籃內有2枚釣鉤理論上處于相同深度,如圖2所示。通過上浮率修正,可獲得延繩釣各釣鉤的估算深度(見表1)。

    表1 延繩釣各釣鉤的估算深度

    Note:①Hook number;②Estimated depth

    調查期間總共隨機抽樣532尾長鰭金槍魚,占總釣獲數量的24.5%。其中,7號鉤(釣鉤號中位數)的捕獲數量最多,達100尾, 1號鉤無漁獲,2號鉤和13號鉤分別僅釣獲7尾和5尾。捕獲的長鰭金槍魚叉長83~110 cm,其中優(yōu)勢叉長95~105 cm,占總數的77.4%。利用多重t檢驗分析不同鉤號上釣漁獲物的叉長(見表2),結果表明2號鉤和13號鉤釣獲的漁獲物叉長與其它鉤均沒有明顯差異(Pgt;0.1),3、4、5、6、7號鉤組之間無明顯差異, 8、9、10、11、12號鉤組之間同樣無明顯差異,但2組之間任意2鉤(號)釣獲的長鰭金槍魚叉長存在明顯差異(Plt;0.1)。不同鉤號釣獲的長鰭金槍魚的尾數(紅色)及叉長頻次分布(藍色)見圖3。分別使用局部加權回歸(紅實線)和線性回歸(綠實線)擬合叉長和鉤號的關系,結果表明2~7號鉤釣獲的長鰭金槍魚的叉長基本接近,隨著鉤深的增加,釣獲的長鰭金槍魚個體大小有略微的增大,8號鉤所處水層(192 m)為一個臨界水層,該水層以深釣獲的長鰭金槍魚的叉長明顯大于該水層以淺釣獲的個體叉長。

    (實線表示懸鏈線理論預測的主繩形態(tài);虛線表示經過修正后的主繩形態(tài)。釣鉤下方的數字表示鉤號。White solid line is the configuration of mainline predicted by catenary algorithm, and white dotted line is the configuration of mainline after correction of shoaling. Ringed numbers below the hooks indicate the hook numbers.)

    圖2 兩浮子間延繩釣釣鉤垂直分布

    注:顯著性水平為0.1。Note:0.1 for significance level.

    (分別使用局部加權回歸(紅實線)和線性回歸(綠實線)擬合叉長和鉤號的關系,紅虛線表示標準差。Locally weighted regression (red line) and linear regression (blue line) are taken for interpreting the relationship between the fitted fork length and hook numbers, dotted line indicate standard deviation.)

    圖3 不同鉤號漁獲物的叉長分布(藍色)及漁獲數量(紅色)

    Fig.3 Fork length distribution (blue) and number of individuals (red) caught by different hooks

    按鉤號和叉長組統(tǒng)計的上鉤漁獲物活體比例(活體數占總數的百分比)見圖4。其中,鉤號與活體比例存在明顯的相關性(Plt;0.05),這表明捕撈水層越深,魚體的存活率越高。不同叉長組的存活率在叉長91 cm及以上的組別中同樣表現出了較為明顯的相關關系,盡管87~88.9和89~90.9兩個叉長組的存活率分別為1和0,但樣本數極少(分別為1和3)。因此,結果表明漁獲物上鉤后的活體比例隨著魚體大小的增加而增加。

    圖4 按鉤號和叉長組統(tǒng)計的上鉤漁獲物活體比例

    2.2 浸泡時間對漁獲率的影響

    試驗期間共投放鉤32次,每次投放14個單位。統(tǒng)計顯示釣鉤最長浸泡時間達1 340 min,最短117 min。試驗期間共釣獲長鰭金槍魚2167尾,漁獲率最高19尾/千鉤,其中釣鉤浸泡時間14~17 h、20~22 h和9~11 h 3個時間段漁獲率較多。對單位浸泡時間(每小時)平均漁獲率進行局部加權回歸,結果顯示,漁獲率整體上隨著浸泡時間的增加呈現波動式的增長,并在11和16 h呈現2個峰值,浸泡時間為9~16 h的單位數量較多(見圖5)。

    (紅實線和虛線表示局部加權回歸及其標準差,右側色階表示各浸泡時間下單位的數量。Red lines stand for locally weighted regression (solid) and standard deviation (dotted). On the right is color gradation indicating the numbers of unit at certain period of soaking.)

    圖5 漁獲率在浸泡時間內的動態(tài)變化

    Fig.5 Dynamic response of catch rate over period of soaking.

    2.3 釣鉤分布深度和漁獲效率的Logistic模型

    Logistic回歸模型顯示釣鉤水深對捕撈效率的影響是極其顯著的(Plt;0.05,見表3),釣鉤分布水深的漁獲效率分布擬合模型(釣鉤分布水深的漁獲效率分布)表示如下:

    (10)

    模型表示167.57 m水層的捕撈概率密度最大。延繩釣鉤深作為有限且數量一定的離散變量,更加集中地布置在該水層上下水體范圍內是更為有效的捕撈方法,因為深于或淺于該水深的釣鉤,其釣捕效率都下降。如果以該最高漁獲物釣捕水深(167.57m)為基準,獲得80%的漁獲物的釣鉤深度置信區(qū)間為124 ~211 m(見圖6,A和B分別表示捕撈效率的分布函數圖和概率密度函數圖),即釣鉤布置在此水深區(qū)間內,將有更高的捕撈成功概率。

    表3 Logistic回歸模型參數估計Table 3 Estimates of modeled parameters by Logistic regression

    Note:①Coefficients;②Estimate;③Standard deviation;④Z statistics

    (紅色實點為鉤深對應漁獲效率的觀測值,藍線為Logistic擬合曲線,陰影部分表示80%漁獲率的水深范圍。圖A和B分別表示漁獲效率的分布函數和概率密度函數。Red dots are observed catch efficiency against hook depth; blue line shows Logistic-fitted curve; shading region indicates the depth range of 80% catch efficiency. Sub-graph A and B show the distribution and probability density, respectively.)

    圖6 漁獲效率和鉤深的Logistic回歸

    Fig.6 Logistic regression between catch efficiency and hook depth

    3 討論

    3.1 長鰭金槍魚個體大小的垂直分布特征

    金槍魚的垂直分布特征與個體大小存在相關性,個體較小的魚分布在較淺水層,這種現象在長鰭金槍魚[18-22]、大眼金槍魚(Thunnusobesus)[23-24]、黃鰭金槍魚(Thunnusalbacares)[25-27]等的檔案式標記追蹤試驗中均有報道。一種觀點認為,金槍魚幼魚無法向成魚一樣能充分發(fā)揮自身熱量保存的機制,因此不能下潛到更深的水層[28];而另一種觀點認為,幼魚在成長過程中食性發(fā)生了改變,成魚往往在更深的水層進行攝食[18]。延繩釣通常捕獲90~105 cm叉長范圍的成年長鰭金槍魚個體[7],本研究抽樣的叉長范圍為83~110 cm。魚體大小存在隨著釣鉤深度增加而增大的關系,但由于體長范圍較窄導致這種關系并不十分密切(見圖3,對比8~13號鉤與1~7號鉤的叉長分布),這可能是由于本研究釣鉤分布在60~230 m之間,體長更小的魚往往棲息在表層,被捕獲的概率較低。例如,Cosgrove等[20]的衛(wèi)星檔案標記數據顯示,叉長78 cm的長鰭金槍魚平均棲息深度為19 m,并在50 m以淺水體內垂直移動;Childers等[29]的標記放流試驗認為,絕大多數時間內長鰭金槍魚幼魚分布在50 m以淺的混合層。然而,魚體大小與棲息水深分布的正相關性應從漁獲物上鉤的晝夜時段以及捕撈區(qū)域進一步予以考證,因為金槍魚種群通常隨著晝夜交替具有“W”或“V”型的大范圍垂直移動模式,即夜晚主要棲息于混合層以淺的表層暖水域,黎明時刻下潛至更深的冷水域。這種行為模式也因不同水域有顯著差異,如溫帶水域的長鰭金槍魚并不表現出白天下潛的行為運動模式[21]。因此,今后的研究應考慮不同捕撈區(qū)域以及漁獲物釣獲時間等因素對于長鰭金槍魚個體大小垂直分布特征的影響。

    3.2 捕撈率在浸泡時間內的動態(tài)響應

    捕獲率在浸泡時間內的動態(tài)響應是極為復雜的,Ward等[30]在分析影響捕撈率的概率事件時,考慮的因素包括餌料入水后可能滑落、可能被海洋動物捕食而發(fā)生丟失、隨著時間推移引誘效果持續(xù)下降、魚體上鉤后可能咬斷釣鉤或者支線而逃逸、或者可能被捕食者掠食等。受各種因素的影響,先前關于浸泡時間對與捕撈率影響的研究結果差異較大。這些差異一方面來源于研究對象或研究區(qū)域的不同;另一方面則歸于研究方法的進化,早期的研究結果傾向于浸泡時間對于捕撈率的不確定性[31-33]。隨著更加復雜的模型和多因素分析方法的引入,能夠更加深入和詳細地解釋兩者關系[34]。本研究排除了因受鯨魚干擾造成的產量下降的所有鉤次,捕撈區(qū)域的鯊魚豐度較低(僅7次鯊魚捕獲的投放鉤次記錄)并且作業(yè)區(qū)域風浪影響較小,因此上鉤漁獲被捕食而丟失以及餌料和上鉤漁獲受海流作用而丟失的可能性降低,這解釋了漁獲率整體上與浸泡時間的正相關性。然而,漁獲率并非隨著浸泡時間單調遞增,而是呈現波動式地增長(11和16 h表現出兩個峰值,或者,13和18 h時捕撈率顯著降低)。其中,第一個峰值與宋利明等[35-36]的研究結論一致,第二個峰值則可能與捕撈時機有關,即某個自然時段目標種群因晝夜垂直移動的范圍與釣鉤分布范圍重疊,由此造成漁獲率增加,例如,有研究認為晝夜時期對漁獲率有著顯著的影響,如黃昏時期作業(yè)的釣鉤對于捕撈率具有積極作用,而黎明時期對于捕撈率影響不大[30]。

    3.3 最佳捕撈效率下的鉤深分布

    繩釣各釣鉤捕撈概率可簡單地表達為一個關系式,即單位時間的捕撈概率密度與釣鉤影響魚群的深度范圍的乘積,并在浸泡時間內進行積分。假設各釣鉤對魚群的影響范圍(釣鉤所在的水深周圍一定的垂直范圍內能夠吸引魚群)和浸泡時間相同,釣鉤的捕撈概率則取決于這個概率密度函數。由于金槍魚的棲息環(huán)境偏好,導致在某一深度內分布最密集。根據Logistic統(tǒng)計模型,本研究認為167.57 m水層的捕撈概率密度最大。盡管長鰭金槍魚的垂直分布范圍很大,但如能將釣鉤布置于捕撈對象概率密度高的水層范圍內,將會獲得最佳的釣捕性能。本研究選取80%捕撈效率作為一個候選捕撈策略,以最高漁獲物釣捕水深(167.57 m)為基準,作業(yè)水深向上和向下延伸至區(qū)間為124~211 m的范圍內無疑有更高的捕撈率。

    鉤位深度精確布置最常用的方法為調整短縮率、浮子繩長、支線長以及投繩的時間或距離間隔。例如,若將本研究中的釣鉤更加集中地布置于124~211 m水層,需要增加浮子繩或支線的長度,并通過降低出繩速度增大短縮率。盡管釣鉤深度的調整是極其方便的,然而投放鉤策略的調整需要更加全面的考慮,如延繩釣總長度、釣鉤分布密度、投鉤數量與勞動力的匹配、作業(yè)循環(huán)的連貫性等,任何因素的調整可能會打破整個策略的“平衡”。例如,降低出繩速度會導致空間內釣鉤密度的增加,釣鉤之間干擾強烈或吸引力重疊(當投繩采用時間固定模式時);或者導致總投鉤數的減少,努力量和勞動力沒有得到有效利用(當投繩采用距離固定模式時)。

    4 結論

    本研究利用東太平洋公海長鰭金槍魚漁業(yè)調查數據,擬合了釣鉤深度分布與捕獲的長鰭金槍魚叉長分布的關系,統(tǒng)計了各鉤號和叉長組的上鉤漁獲物的存活比例,并通過Logistic回歸模型對釣捕水深和捕撈效率之間的關系進行了統(tǒng)計分析,研究結論如下:

    (1)魚體大小與釣鉤的深度之間存在微弱的正相關,192 m(8號鉤所處水層)為臨界水層,該水層以淺的長鰭金槍魚叉長明顯小于該水層以深的個體;捕撈水層越深,被捕魚體的存活率越高。

    (2)釣鉤浸泡時間大多集中在9~16 h,漁獲率整體上隨著浸泡時間的增加呈現波動式的增長,并在11和16 h時表現出2個峰值。

    (3)釣鉤水深對捕撈效率的影響極其顯著,167.57 m水層的捕撈概率密度最大,以該水深為基準,80%的漁獲物在124~211 m水深之間捕獲。

    [2] Williams A J, Allain V, Nicol S J, et al. Vertical behavior and diet of albacore tuna(Thunnusalalunga) vary with latitude in the South Pacific Ocean[J]. Deep-Sea Research II, 2015, 113: 154-169.

    [3] Zárate V O, Babcock E A. Estimating individual growth variability in albacore(Thunnusalalunga) from the North Atlantic stock: Aging for assessment purposes[J]. Fisheries Research, 2015, 180: 54-66.

    [4] Wells D, Kohin S, Teo S, et al. Age and growth of North Pacific albacore(Thunnusalalunga): Implications for stock assessment[J]. Fisheries Research, 2013, 147: 55-62.

    [5] FAO. The State of World Fisheries and Aquaculture[M]. Rome: FAO, 2016.

    [6] Langley A, Hampton J. Stock assessment of albacore tuna in the south Pacific Ocean[R]. //Working Paper SA-WP-3, 1st Scientific Committee meeting of the Western and Central Pacific Fisheries Commission. New Caledonia: Noumea, 2005.

    [7] Harley S, Williams P, Nicol S, et al. The western and central Pacific tuna fishery: 2014 overview and status of stocks[R]. //Tuna Fisheries Assessment Report 15. New Caledonia: Noumea, 2015.

    [8] Campbell R A, Young J W. Monitoring the behaviour of longline gears and the depth and time of fish capture in the Australian Eastern Tuna and Billfish Fishery[J]. Fisheries Research, 2012, 119-120: 48-65.

    [9] Yokawa K, Takeuchi Y, Okazaki M, et al. Standardizations of CPUE of blue marlin and white marlin caught by Japanese longliners in the Atlantic Ocean[J]. Col. Vol. Sci. Pap. ICCAT, 2001, 53: 345-355.

    [10] Boggs C H. Depth, capture time, and hooked longevity of longline-caught pelagic fish: timing bites of fish with chips[J]. Fishery Bulletin, 1992, 90: 642-658.

    [11] Satoh K, Kusaga I, Takeda S, et al. Tuna long-line tests using monofilament Nylon[J]. Nippon Suisan Gakkaishi, 1990, 56(10): 1605-1609.

    [12] Bigelow K A, Hampton J, Miyabe N. Application of a habitat-based model to estimate effective longline fishing effort and relative abundance of Pacific bigeye tuna(Thunnusobesus)[J]. Fisheries Oceanography, 2002, 11: 143-155.

    [13] Ward P, Myers R A. Inferring the depth distribution of catchability for pelagic fishes and correcting for variations in the depth of longline fishing gear[J]. Canadian Journal of Fisheries and Aquatic Sciences, 2005, 62: 1130-1142.

    [14] Takeuchi Y. Is historically available hooks-per-basket information enough to standardize actual hooks-per-basket effects on CPUE? Preliminary simulation approach[J]. ICCAT Coll Vol Sci Pap, 2001, 53: 356-364.

    [15] Yoshihara T. Distribution of fishes caught by the longline -II. Vertical distribution[J]. Bulletin of the Japanese Society of Scientific Fisheries, 1951, 16: 370-374.

    [16] Bigelow K A, Michael K, Musyl M K, et al. Pelagic longline gear depth and shoaling[J]. Fisheries Research, 2006, 77: 173-183.

    [17] Bach P, Gaernter D, Menkes C, et al. Effects of the gear deployment strategy and current shear on pelagic longline shoaling[J]. Fisheries Research, 2009, 95: 55-64.

    [19] Cosgrove R, Arregui I, Arrizablaga H, et al. Predation of pop-up satellite archival tagged albacore(Thunnusalalunga)[ J]. Fisheries Research, 2015, 162: 48-52.

    [20] Cosgrove R, Arregui I, Arrizabalaga H, et al. New insights to behaviour of North Atlantic albacore tuna(Thunnusalalunga) observed with pop-up satellite archival tags[J]. Fisheries Research, 2014, 150: 89-99.

    [21] Williams A J, Allain V, Nicol S, et al. Vertical behavior and diet of albacore tuna(Thunnusalalunga) vary with latitude in the South Pacific Ocean[J]. Deep-Sea Research II, 2015, 113: 154-169.

    [22] 翟天晨, 戴小杰, 朱江峰, 等. 東太平洋長鰭金槍魚個體大小與釣獲深度的關系[J]. 海洋漁業(yè), 2015, 37(1): 10-16.

    Zhai Tianchen, Dai Xiaojie, Zhu Jiangfeng. Preliminary research on size difference and depth environment ofThunnusalalungain East Pacific[J]. Marine Fisheries, 2015, 37(1): 10-16.

    [23] Fuller D W, Schaefer K M, Hampton J, et al. Vertical movements, behavior, and habitat of bigeye tuna(Thunnusobesus) in the equatorial central Pacific Ocean[J]. Fisheries Research, 2015, 172: 57-70.

    [24] Musyl M K, Brill R W, Boggs C H, et al. Vertical movements of bigeye tuna(Thunnusobesus) associated with islands, buoys, and seamounts near the main Hawaiian Islands from archival tagging data[J]. Fisheries Oceanography, 2003, 12(3): 152-169.

    [25] Dagorn L, Josse E, Bach P. Individual differences in horizontal movements of yellowfin tuna(Thunnusalbacares) in nearshore areas in French Polynesia, determined using ultrasonic telemetry[J]. Aquatic Living Resources, 2000, 13: 193-202.

    [26] Schaefer K M, Fuller D W, Block B A. Movements, behavior, and habitat utilization of yellowfin tuna(Thunnusalbacares) in the Pacific Ocean off Baja California, Mexico, determined from archival tag data analyses, including unscented Kalman filtering[J]. Fisheries Research, 2011, 112: 22-37.

    [27] Schaefer K M, Fuller D W, Aldana G. Movements, behavior, and habitat utilization of yellowfin tuna(Thunnusalbacares) in waters surrounding the Revillagigedo Islands Archipelago Biosphere Reserve, Mexico[J]. Fisheries Oceanography, 2014, 23(1): 65-82.

    [28] Dickson K A, Johnson N M, Donley J M, et al. Ontogenetic changes in characteristics required for endothermy in juvenile black skipjack tuna(Euthynnuslineatus)[J]. Journal of Experimental Biology, 2000, 203(20): 3077-3087.

    [29] Childers J, Snyser S, Kohin S. Migration and behavior of juvenile North Pacific albacore(Thunnusalalunga)[J]. Fisheries Oceanography, 2011, 20(3): 157-173.

    [30] Ward P, Myers R A, Blanchard W. Fish lost at sea: the effect of soak time on pelagic longline catches[J]. Fishery Bulletin, 2005, 102(1): 179-195.

    [31] L?kkeborg S, Pina T. Effects of setting time, setting direction and soak time on longline catch rates[J]. Fisheries Research, 1997, 32: 213-222.

    [32] Sivasubramaniam K. Relation between soaking time and catch of tunas in longline fisheries[J]. Bulletin of the Japanese Society of Scientific Fisheries, 1961, 27: 835-845.

    [33] Skud B E. Factors Affecting Longline Catch and Effort: Ⅲ. Bait Loss and Competition[R]. Seattle: The International Pacific Halibut Commission Scientific Report, 1978: 64-66.

    [34] Carruthers E H, Neilson J D, Smith S C. Overlooked bycatch mitigation opportunities in pelagic longline fisheries: Soak time and temperature effects on swordfish(Xiphiasgladius) and blue shark(Prionace glauca) catch[J]. Fisheries Research, 2011, 108: 112-120.

    [35] 宋利明, 徐偉云, 曹道梅, 等. 金槍魚延繩釣釣具的最適浸泡時間[J]. 中國水產科學, 2013(2): 346-350.

    Song Lingming, Xu Weiyun, Cao Daomei, et al. Optimum soak time of tuna longline gear[J]. Journal of Fishery Sciences of China, 2013(2): 346-350.

    [36] 宋利明, 李冬靜, 劉海陽, 等. 主捕長鰭金槍魚延繩釣釣具的最適浸泡時間[J]. 上海海洋大學學報, 2014,23(2): 290-296.

    Song Lingming, Li Dongjing, Liu Haiyang, et al. Optimum soak time of pelagic longline gear targeting albacore tuna[J]. Journal of Shanghai Ocean University, 2014, 23(2): 290-296.

    責任編輯 朱寶象

    IdentifyingtheCatchCharacteristicsofAlbacoreTuna(Thunnusalalunga)LonglineAssociatedwiththeDepthandSoakingTimeofHooksintheHighSeasofEasternPacificOcean

    LIU Li-Li1, 2, ZHOU Cheng3, YU Cong-Da1, 2, ZHENG Ji1, 2, ZANG Ying-Liang1, 2

    (1. College of Fisheries, Zhejiang Ocean University, Zhoushan 316022, China; 2. Key Laboratory of Marine Fishery Equipment and Technology of Zhejiang, Zhoushan 316022, China; 3. College of Fisheries, Ocean University of China, Qingdao 266003, China)

    For the intention of further understanding the fishing capability and catch characteristics associated with gear and fishing operation, this study, on the basis of tuna longline fishery-dependent data from a survey in the high seas of eastern Pacific Ocean spanning from September 2015 to February 2016, performed the locally weighted regression and linear regression of hook depths to the fork length distribution, and statistics of at-vessel survival rate grouped by hooks number and fork length, as well as the dynamic response of catch rate over the period of soaking based on a logistic regression model. The results showed that, despite a less powerful relationship, the fish size remained to increase with the hooking depth. Hook No. 8 located at a critical depth that the fork length of albacore tuna caught above this depth was significantly smaller than that bellow this depth. At-vessel survival rates showed evidently positive dependence with the hook depth and fork length. Catch rate responding to the soaking time based on eleven consecutive floats as unit to calculate exhibited a growth in fluctuation with bimodal characteristics at period of 11 and 16 h. Logistic regression model suggested a significant effect of hook depth on the catch efficiency. The water layer for the highest density of catch efficiency located at depth of 167.57 m. An alternative strategy with the hooks deploying at depth ranging from 124 to 211 m will result in a more considerable fishing efficiency.

    tuna longline; hook depth; soaking time; catch rate; catch characteristics

    S973.3

    A

    1672-5174(2018)01-040-09

    10.16441/j.cnki.hdxb. 20170076

    劉莉莉, 周成, 虞聰達, 等. 釣鉤深度和浸泡時間對東太平洋公海長鰭金槍魚延繩釣漁獲性能的影響研究[J]. 中國海洋大學學報(自然科學版), 2018, 48(1): 40-48.

    LIU Li-Li, ZHOU Cheng, YU Cong-Da, et al. Identifying the catch characteristics of albacore tuna (Thunnusalalunga) longline associated with the depth and soaking time of hooks in the high seas of eastern Pacific Ocean[J].Periodical of Ocean University of China, 2018, 48(1): 40-48.

    浙江省公益技術研究社會發(fā)展項目(2016C33083);國家自然科學基金青年基金項目(41606110;41506151);浙江省自然科學基金青年基金項目(LQ14C190002);2015年東太平洋公海金槍魚資源探捕項目資助

    Supported by Social Development Project of Zhejiang Public Welfare Technology Research (2016C33083); National Natural Science Foundation (41606110; 41506151); Natural Science Foundation of Zhejiang Province (LQ14C190002); Project for High Seas of East Pacific Ocean Resources Exploration in 2015

    2017-02-23;

    2017-04-29

    劉莉莉(1991-),女,博士,講師。E-mail:lililiuxing@hotmail.com

    ?? 通訊作者:E-mail:zhoucheng286@126.com

    猜你喜歡
    漁獲金槍魚水層
    黃渤海不同水層中浮游植物對灰霾添加的響應
    金槍魚淹死事件
    趣味(數學)(2022年3期)2022-06-02 02:32:46
    可以加速的金槍魚機器人
    軍事文摘(2021年22期)2022-01-18 06:22:00
    金槍魚與海豚
    首批國家級海洋捕撈漁獲物定點上岸漁港廣西唯一入選名額落戶南漁港
    農家之友(2020年9期)2020-04-21 02:18:42
    東海北部桁桿拖蝦漁獲組成分析
    國內新型遠洋金槍魚圍網漁船首航
    商周刊(2018年19期)2018-12-06 09:49:41
    南海北部深海區(qū)燈光罩網漁獲物組成及漁獲率的時空分布
    高壓電纜大截面分割導體焊接后的機械性能及緩沖阻水層設計
    電線電纜(2018年4期)2018-08-31 05:57:30
    水稻水層管理田間試驗總結
    国产成人av激情在线播放| 欧美日韩一区二区视频在线观看视频在线| 熟女av电影| videosex国产| bbb黄色大片| 国产亚洲午夜精品一区二区久久| 婷婷色综合www| 乱人伦中国视频| 亚洲国产日韩一区二区| 如日韩欧美国产精品一区二区三区| 多毛熟女@视频| 超碰成人久久| 丝瓜视频免费看黄片| 精品人妻熟女毛片av久久网站| 亚洲精品视频女| 成人三级做爰电影| 美女国产高潮福利片在线看| 亚洲精品久久成人aⅴ小说| 免费女性裸体啪啪无遮挡网站| 久久久精品国产亚洲av高清涩受| 精品卡一卡二卡四卡免费| 国产一级毛片在线| 五月开心婷婷网| 国产 一区精品| 制服人妻中文乱码| 妹子高潮喷水视频| 免费在线观看完整版高清| 亚洲国产欧美一区二区综合| 丰满乱子伦码专区| 麻豆精品久久久久久蜜桃| 丰满饥渴人妻一区二区三| 成年动漫av网址| 国产精品一区二区在线观看99| 99香蕉大伊视频| 国产成人av激情在线播放| 日韩中文字幕欧美一区二区 | 国产xxxxx性猛交| 黑丝袜美女国产一区| 国产欧美日韩综合在线一区二区| 好男人视频免费观看在线| 999久久久国产精品视频| 成年人免费黄色播放视频| 91国产中文字幕| 又粗又硬又长又爽又黄的视频| 午夜91福利影院| 黑人巨大精品欧美一区二区蜜桃| 亚洲男人天堂网一区| 亚洲人成电影观看| 欧美日韩视频高清一区二区三区二| 亚洲精品第二区| 久热这里只有精品99| 永久免费av网站大全| 色综合欧美亚洲国产小说| 我要看黄色一级片免费的| 国产日韩欧美亚洲二区| 国产福利在线免费观看视频| 亚洲欧洲精品一区二区精品久久久 | 校园人妻丝袜中文字幕| 狠狠精品人妻久久久久久综合| 老汉色av国产亚洲站长工具| 国产一区二区在线观看av| 丝袜美足系列| 亚洲第一av免费看| 久久久久久久久久久久大奶| 成人漫画全彩无遮挡| 亚洲熟女毛片儿| 亚洲欧美成人精品一区二区| 69精品国产乱码久久久| 亚洲成人手机| 免费少妇av软件| 欧美xxⅹ黑人| 亚洲av电影在线观看一区二区三区| 秋霞在线观看毛片| 91精品三级在线观看| 亚洲熟女精品中文字幕| av国产久精品久网站免费入址| 亚洲av成人不卡在线观看播放网 | 人人妻人人爽人人添夜夜欢视频| 欧美日韩一级在线毛片| av卡一久久| 伊人久久大香线蕉亚洲五| 五月开心婷婷网| 一区二区av电影网| 老司机深夜福利视频在线观看 | 日韩制服丝袜自拍偷拍| 国产伦人伦偷精品视频| 国产精品欧美亚洲77777| 97人妻天天添夜夜摸| 日韩一区二区视频免费看| 国产精品 国内视频| 国产人伦9x9x在线观看| 免费看av在线观看网站| 久久女婷五月综合色啪小说| 日日撸夜夜添| 欧美日韩精品网址| 人妻一区二区av| 欧美日韩视频精品一区| 国精品久久久久久国模美| www.自偷自拍.com| 男女免费视频国产| 精品第一国产精品| 久久精品国产a三级三级三级| 久久国产精品大桥未久av| 国产97色在线日韩免费| 在线 av 中文字幕| av视频免费观看在线观看| 这个男人来自地球电影免费观看 | 欧美日韩亚洲综合一区二区三区_| 国产97色在线日韩免费| 亚洲精品成人av观看孕妇| 国产伦理片在线播放av一区| 国产国语露脸激情在线看| 99久久99久久久精品蜜桃| 精品亚洲成a人片在线观看| 欧美日韩一区二区视频在线观看视频在线| 国产视频首页在线观看| 亚洲伊人色综图| 在线天堂最新版资源| 在线观看三级黄色| 啦啦啦啦在线视频资源| 亚洲中文av在线| 最新在线观看一区二区三区 | 国产精品一二三区在线看| 亚洲国产日韩一区二区| 在线观看三级黄色| 婷婷色av中文字幕| 美女脱内裤让男人舔精品视频| 成年女人毛片免费观看观看9 | 午夜精品国产一区二区电影| 叶爱在线成人免费视频播放| 亚洲精品国产色婷婷电影| 在线观看www视频免费| 啦啦啦 在线观看视频| 亚洲av成人精品一二三区| 一级毛片我不卡| 久久国产精品男人的天堂亚洲| 天天操日日干夜夜撸| 国产人伦9x9x在线观看| 亚洲国产精品国产精品| 色精品久久人妻99蜜桃| 最近中文字幕高清免费大全6| 国产爽快片一区二区三区| 叶爱在线成人免费视频播放| 国产亚洲最大av| 老汉色∧v一级毛片| 夫妻性生交免费视频一级片| 国产片内射在线| 亚洲精品中文字幕在线视频| 亚洲七黄色美女视频| 欧美另类一区| 亚洲国产精品国产精品| 成人毛片60女人毛片免费| 亚洲精品国产一区二区精华液| 秋霞在线观看毛片| 精品国产露脸久久av麻豆| 高清av免费在线| 亚洲三区欧美一区| 一级毛片黄色毛片免费观看视频| 欧美中文综合在线视频| 18禁观看日本| 在线看a的网站| 丰满乱子伦码专区| 女人高潮潮喷娇喘18禁视频| 大香蕉久久成人网| 亚洲自偷自拍图片 自拍| 久久久精品免费免费高清| 男女边摸边吃奶| 国产成人av激情在线播放| 人人妻,人人澡人人爽秒播 | 欧美人与性动交α欧美精品济南到| 男男h啪啪无遮挡| 久久狼人影院| 久久久久精品久久久久真实原创| 侵犯人妻中文字幕一二三四区| 99国产精品免费福利视频| 99热网站在线观看| 九色亚洲精品在线播放| 亚洲国产日韩一区二区| 免费av中文字幕在线| 婷婷色av中文字幕| 欧美精品一区二区大全| 久久ye,这里只有精品| 久久人人97超碰香蕉20202| 免费女性裸体啪啪无遮挡网站| 精品一区在线观看国产| 亚洲少妇的诱惑av| 亚洲美女黄色视频免费看| 女人爽到高潮嗷嗷叫在线视频| 在线观看免费视频网站a站| 如日韩欧美国产精品一区二区三区| 日韩欧美精品免费久久| 一本一本久久a久久精品综合妖精| 亚洲国产av影院在线观看| 黄色视频在线播放观看不卡| 老汉色av国产亚洲站长工具| 国产亚洲午夜精品一区二区久久| 精品亚洲成a人片在线观看| 婷婷色麻豆天堂久久| 啦啦啦视频在线资源免费观看| 青春草视频在线免费观看| 哪个播放器可以免费观看大片| av卡一久久| 精品人妻在线不人妻| 精品福利永久在线观看| 免费久久久久久久精品成人欧美视频| 韩国精品一区二区三区| 国产成人精品福利久久| 色吧在线观看| 亚洲专区中文字幕在线 | 999精品在线视频| 国产伦人伦偷精品视频| 久久99热这里只频精品6学生| tube8黄色片| 国产一区二区三区av在线| 中文字幕精品免费在线观看视频| 大香蕉久久网| 青春草亚洲视频在线观看| 人人妻,人人澡人人爽秒播 | 日韩中文字幕视频在线看片| 男人舔女人的私密视频| 久久精品亚洲av国产电影网| a 毛片基地| 久久鲁丝午夜福利片| 亚洲av福利一区| 国产精品久久久久久久久免| 十八禁人妻一区二区| 日韩欧美精品免费久久| 亚洲中文av在线| 丰满乱子伦码专区| 国产色婷婷99| 成人国产av品久久久| 交换朋友夫妻互换小说| 9191精品国产免费久久| 精品少妇黑人巨大在线播放| 黑人猛操日本美女一级片| 69精品国产乱码久久久| 国产乱来视频区| 国产有黄有色有爽视频| 男男h啪啪无遮挡| 1024香蕉在线观看| av网站在线播放免费| 免费黄网站久久成人精品| 亚洲av中文av极速乱| 亚洲欧美一区二区三区久久| 美女高潮到喷水免费观看| 国产精品亚洲av一区麻豆 | 国产伦人伦偷精品视频| 又粗又硬又长又爽又黄的视频| 亚洲av综合色区一区| 久久女婷五月综合色啪小说| 在线观看三级黄色| 亚洲成人一二三区av| 叶爱在线成人免费视频播放| 亚洲精品国产色婷婷电影| 中文字幕制服av| av.在线天堂| 亚洲av福利一区| 亚洲熟女毛片儿| 国产一卡二卡三卡精品 | 青青草视频在线视频观看| 亚洲美女搞黄在线观看| 免费高清在线观看日韩| e午夜精品久久久久久久| netflix在线观看网站| 国产一卡二卡三卡精品 | 成人免费观看视频高清| 精品一区二区三区四区五区乱码 | 亚洲av电影在线观看一区二区三区| 一区二区三区激情视频| 天天操日日干夜夜撸| 99九九在线精品视频| 国产伦理片在线播放av一区| 亚洲欧洲精品一区二区精品久久久 | 色94色欧美一区二区| 黑人欧美特级aaaaaa片| 久久婷婷青草| 日韩中文字幕视频在线看片| 免费黄频网站在线观看国产| 性少妇av在线| 熟妇人妻不卡中文字幕| 日本午夜av视频| 最近2019中文字幕mv第一页| 精品免费久久久久久久清纯 | 亚洲av在线观看美女高潮| 免费人妻精品一区二区三区视频| 亚洲欧美精品综合一区二区三区| 欧美在线一区亚洲| 在现免费观看毛片| 免费高清在线观看日韩| 嫩草影院入口| 日韩熟女老妇一区二区性免费视频| 国产人伦9x9x在线观看| 久久国产乱子伦精品免费另类| 在线观看日韩欧美| 不卡一级毛片| 欧美成狂野欧美在线观看| 国产精品免费视频内射| 大陆偷拍与自拍| 夜夜躁狠狠躁天天躁| 午夜a级毛片| 欧美黄色淫秽网站| 久久人人精品亚洲av| 午夜精品久久久久久毛片777| 一级a爱视频在线免费观看| 一夜夜www| 一区二区三区国产精品乱码| 国产成年人精品一区二区| 成人18禁在线播放| 午夜福利一区二区在线看| 色综合欧美亚洲国产小说| 男女午夜视频在线观看| 桃色一区二区三区在线观看| 亚洲伊人色综图| 免费在线观看黄色视频的| 麻豆成人av在线观看| 两个人视频免费观看高清| 香蕉久久夜色| 黑丝袜美女国产一区| 久久 成人 亚洲| 亚洲成av片中文字幕在线观看| www.熟女人妻精品国产| 婷婷六月久久综合丁香| 99久久久亚洲精品蜜臀av| 叶爱在线成人免费视频播放| 人人澡人人妻人| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品国产高清国产av| 淫秽高清视频在线观看| 成人18禁高潮啪啪吃奶动态图| 精品久久久久久,| 99精品久久久久人妻精品| av电影中文网址| 国产91精品成人一区二区三区| 桃色一区二区三区在线观看| 欧美日韩黄片免| 久久亚洲真实| 女人爽到高潮嗷嗷叫在线视频| 国产私拍福利视频在线观看| aaaaa片日本免费| 亚洲欧美激情在线| 香蕉丝袜av| 黑人欧美特级aaaaaa片| 美女午夜性视频免费| 欧美色视频一区免费| 国产成+人综合+亚洲专区| 日本精品一区二区三区蜜桃| www国产在线视频色| 久久久久久人人人人人| 69av精品久久久久久| 日韩大尺度精品在线看网址 | 最近最新中文字幕大全免费视频| 亚洲情色 制服丝袜| 国产男靠女视频免费网站| 麻豆av在线久日| 在线观看舔阴道视频| 少妇被粗大的猛进出69影院| 精品国产一区二区久久| 女同久久另类99精品国产91| 久久精品人人爽人人爽视色| 欧美国产日韩亚洲一区| 一区二区三区高清视频在线| 香蕉丝袜av| 人人澡人人妻人| 国产成人系列免费观看| 制服诱惑二区| 岛国在线观看网站| 欧美激情 高清一区二区三区| 亚洲欧美精品综合久久99| 美女免费视频网站| 欧美日韩亚洲国产一区二区在线观看| 两个人免费观看高清视频| 欧美另类亚洲清纯唯美| 国产精品久久久久久精品电影 | 黄网站色视频无遮挡免费观看| 亚洲精品久久国产高清桃花| 女同久久另类99精品国产91| 久久午夜亚洲精品久久| 天天躁夜夜躁狠狠躁躁| 非洲黑人性xxxx精品又粗又长| 此物有八面人人有两片| 国产单亲对白刺激| 男女下面进入的视频免费午夜 | 精品国内亚洲2022精品成人| 国产一区二区激情短视频| 99国产精品99久久久久| 亚洲性夜色夜夜综合| 国产精品亚洲美女久久久| 国产亚洲精品久久久久5区| 欧美性长视频在线观看| 一级,二级,三级黄色视频| 黄色成人免费大全| 午夜免费激情av| 亚洲精品国产区一区二| 在线观看免费午夜福利视频| 男女午夜视频在线观看| 国产亚洲欧美精品永久| 日韩欧美一区二区三区在线观看| 国产日韩一区二区三区精品不卡| АⅤ资源中文在线天堂| 色播在线永久视频| 黄频高清免费视频| 老司机在亚洲福利影院| 成人手机av| 身体一侧抽搐| 国产精品秋霞免费鲁丝片| 涩涩av久久男人的天堂| 色av中文字幕| 久久中文看片网| 看免费av毛片| 精品久久久久久久毛片微露脸| 咕卡用的链子| 一个人观看的视频www高清免费观看 | 日韩欧美一区二区三区在线观看| 淫妇啪啪啪对白视频| 黄色视频不卡| 国产亚洲精品一区二区www| 国产伦人伦偷精品视频| 久热爱精品视频在线9| 黄片大片在线免费观看| ponron亚洲| 欧美色欧美亚洲另类二区 | 亚洲情色 制服丝袜| 亚洲熟妇熟女久久| 女人爽到高潮嗷嗷叫在线视频| 亚洲七黄色美女视频| 黄色成人免费大全| 一区在线观看完整版| 亚洲全国av大片| 女人精品久久久久毛片| 亚洲精品一区av在线观看| 欧美人与性动交α欧美精品济南到| 侵犯人妻中文字幕一二三四区| 精品免费久久久久久久清纯| 色综合婷婷激情| 日韩三级视频一区二区三区| 在线播放国产精品三级| 美女大奶头视频| 丝袜美足系列| 久久久国产欧美日韩av| 欧美日韩一级在线毛片| 欧美日韩精品网址| www.999成人在线观看| 亚洲欧洲精品一区二区精品久久久| 国产又爽黄色视频| 麻豆成人av在线观看| 51午夜福利影视在线观看| 日日摸夜夜添夜夜添小说| 岛国视频午夜一区免费看| 亚洲第一电影网av| 国产又色又爽无遮挡免费看| 日本撒尿小便嘘嘘汇集6| 黄色 视频免费看| 一级,二级,三级黄色视频| 一本综合久久免费| 久久久国产精品麻豆| 亚洲精品av麻豆狂野| 日韩欧美三级三区| 大码成人一级视频| 亚洲精品一区av在线观看| 美女 人体艺术 gogo| 成人av一区二区三区在线看| 久久精品91无色码中文字幕| 曰老女人黄片| 日日干狠狠操夜夜爽| 在线观看舔阴道视频| 久久精品人人爽人人爽视色| 亚洲 欧美一区二区三区| 日韩欧美三级三区| 多毛熟女@视频| 美女午夜性视频免费| 免费不卡黄色视频| 老熟妇仑乱视频hdxx| 十八禁人妻一区二区| 国产精品98久久久久久宅男小说| 狠狠狠狠99中文字幕| av电影中文网址| 精品久久久久久久毛片微露脸| 久久 成人 亚洲| 19禁男女啪啪无遮挡网站| 怎么达到女性高潮| 国产成年人精品一区二区| 午夜久久久久精精品| 国产免费av片在线观看野外av| 亚洲avbb在线观看| 美女午夜性视频免费| 91麻豆av在线| 国产精品久久久av美女十八| 欧美中文日本在线观看视频| 久久国产乱子伦精品免费另类| 亚洲成a人片在线一区二区| 亚洲全国av大片| 99久久国产精品久久久| 99久久久亚洲精品蜜臀av| 别揉我奶头~嗯~啊~动态视频| 亚洲一区二区三区色噜噜| 国产熟女午夜一区二区三区| 悠悠久久av| 真人一进一出gif抽搐免费| 身体一侧抽搐| av福利片在线| 波多野结衣高清无吗| 女人精品久久久久毛片| 久久天堂一区二区三区四区| 中文字幕高清在线视频| 人成视频在线观看免费观看| 久久久久久人人人人人| 欧美日韩亚洲国产一区二区在线观看| 久久精品影院6| 久久精品亚洲熟妇少妇任你| 两性夫妻黄色片| 村上凉子中文字幕在线| 日韩欧美在线二视频| 亚洲精品粉嫩美女一区| 女性被躁到高潮视频| 真人一进一出gif抽搐免费| 欧美av亚洲av综合av国产av| 亚洲在线自拍视频| 性色av乱码一区二区三区2| 午夜两性在线视频| 两人在一起打扑克的视频| 久久天堂一区二区三区四区| 久久精品国产综合久久久| 久久精品成人免费网站| 国产av在哪里看| 黄色丝袜av网址大全| 国产成人精品在线电影| 麻豆成人av在线观看| x7x7x7水蜜桃| 91麻豆精品激情在线观看国产| 成人精品一区二区免费| 亚洲精华国产精华精| 天天躁夜夜躁狠狠躁躁| svipshipincom国产片| 亚洲无线在线观看| 麻豆成人av在线观看| 老司机午夜十八禁免费视频| 亚洲激情在线av| 搡老熟女国产l中国老女人| 久久久久九九精品影院| 亚洲自拍偷在线| 精品第一国产精品| 免费女性裸体啪啪无遮挡网站| 国产午夜福利久久久久久| 黄色片一级片一级黄色片| 婷婷精品国产亚洲av在线| 可以免费在线观看a视频的电影网站| avwww免费| 9色porny在线观看| 精品久久久久久久毛片微露脸| 侵犯人妻中文字幕一二三四区| 性少妇av在线| 欧美黄色片欧美黄色片| 两性午夜刺激爽爽歪歪视频在线观看 | 神马国产精品三级电影在线观看 | 午夜福利成人在线免费观看| 国产一区二区在线av高清观看| 久久久久国产精品人妻aⅴ院| 午夜精品久久久久久毛片777| 国产精品香港三级国产av潘金莲| 国产1区2区3区精品| 国产精品久久视频播放| 国产男靠女视频免费网站| 欧美久久黑人一区二区| 亚洲欧美一区二区三区黑人| 欧美一区二区精品小视频在线| 法律面前人人平等表现在哪些方面| 视频区欧美日本亚洲| 在线观看免费视频日本深夜| 最新在线观看一区二区三区| 免费高清视频大片| 777久久人妻少妇嫩草av网站| 很黄的视频免费| 黄片小视频在线播放| 一本大道久久a久久精品| 国产熟女午夜一区二区三区| 一区二区日韩欧美中文字幕| 精品卡一卡二卡四卡免费| 亚洲五月婷婷丁香| 搡老岳熟女国产| 99精品在免费线老司机午夜| 一级,二级,三级黄色视频| 最好的美女福利视频网| 怎么达到女性高潮| av超薄肉色丝袜交足视频| 涩涩av久久男人的天堂| 一级毛片女人18水好多| 国产精品久久久久久人妻精品电影| www.自偷自拍.com| 九色国产91popny在线| 最近最新中文字幕大全免费视频| 伊人久久大香线蕉亚洲五| 成人18禁高潮啪啪吃奶动态图| 18禁美女被吸乳视频| 天天躁夜夜躁狠狠躁躁| 最近最新中文字幕大全电影3 | 国产欧美日韩一区二区精品| 国产高清激情床上av| 在线观看舔阴道视频| 亚洲色图 男人天堂 中文字幕| 国产又色又爽无遮挡免费看| 99久久99久久久精品蜜桃| 如日韩欧美国产精品一区二区三区| 亚洲片人在线观看| 女人爽到高潮嗷嗷叫在线视频| 99在线视频只有这里精品首页| 黄色女人牲交| 亚洲专区字幕在线| 可以在线观看的亚洲视频| 1024香蕉在线观看| 99国产精品一区二区三区| 日本一区二区免费在线视频| 精品乱码久久久久久99久播| 女性被躁到高潮视频| 日本撒尿小便嘘嘘汇集6| 国产一区二区三区综合在线观看|