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

    長江口汛期鳳鱭捕撈量年際變化研究

    2020-06-10 01:45:14倪健夫郭弘藝唐文喬張旭光
    海洋漁業(yè) 2020年2期
    關鍵詞:輸沙量長江口年際

    倪健夫 ,郭弘藝,唐文喬,,張 亞,張旭光

    (1. 上海海洋大學海洋動物系統(tǒng)分類與進化上海高校重點實驗室,上海 201306; 2. 上海海洋大學科技部海洋生物學國際聯(lián)合研究中心,上海 201306; 3. 上海海洋大學水產(chǎn)科學國家級實驗教學示范中心,上海 201306)

    鳳鱭(Coiliamystus)俗稱鳳尾魚、烤籽魚,隸屬于鯡形目(Clupeiformes),鳀科(Engraulidae),在我國沿海均有分布,被劃分為長江、閩江和珠江3個生態(tài)型群體,是一種短距離溯河洄游魚類[1-2]。長江口是我國鳳鱭漁業(yè)的主要產(chǎn)區(qū),1960—1998年鳳鱭年均產(chǎn)量為1 192 t[3],在20世紀80年代,產(chǎn)量更是高達長江口漁獲物總量的48.6%[4]。自20世紀90年代開始,隨著長江中下游流域社會經(jīng)濟的快速發(fā)展、水域污染的加劇及對資源的長期過度捕撈,鳳鱭捕撈量迅速下滑,生物學參數(shù)和懷卵量等繁殖特征均發(fā)生了顯著改變[5-8]。而2003年后三峽工程的蓄水,改變了長江入海徑流和泥沙輸運原有的季節(jié)與年際變化格局[9]。河口環(huán)境因子的改變,勢必影響到生態(tài)系統(tǒng)中生物群落的分布格局和漁業(yè)資源的適應性,鳳鱭捕撈量呈現(xiàn)近一步大幅下降[10]。近年來,已有學者對鳳鱭的資源變動[6-8, 10-11]、攝食習性[12-13]、生活史[14]和種群鑒別[15-19]及脂肪含量等生化成分[20-21]進行了研究報道,但未見關于長江口鳳鱭捕撈量格局轉(zhuǎn)變的研究。格局轉(zhuǎn)變貫序t檢驗法(sequential t-test testing analysis of regime shifts,STARS)可探測穩(wěn)態(tài)改變的發(fā)生時間、次數(shù)、持續(xù)時間及各穩(wěn)態(tài)階段均值,較傳統(tǒng)的Mann-Kendall檢驗和滑動t檢驗,對時間序列較短的數(shù)據(jù)檢驗更為嚴謹和準確[22-23],目前已應用于海洋學、生態(tài)學、地質(zhì)與氣候科學、環(huán)境科學、漁業(yè)資源等領域[22-29]。本文利用1982—2014年長江口鳳鱭捕撈量數(shù)據(jù),結合捕撈船數(shù)及環(huán)境因子水溫、徑流和輸沙量進行年際格局轉(zhuǎn)變研究,找出格局轉(zhuǎn)變年份。并結合廣義可加模型(generalized additive model,GAM)定量分析了鳳鱭捕撈量與捕撈船數(shù)及環(huán)境因子的相關性,旨在弄清長江口鳳鱭捕撈量格局轉(zhuǎn)變的影響因子,為鳳鱭資源的保護及可持續(xù)利用提供基礎數(shù)據(jù)。

    1 材料與方法

    1.1 數(shù)據(jù)來源

    本研究中,1982—2014年長江口鳳鱭汛期(4—6月)的捕撈船數(shù)(艘·年-1)、汛期總捕撈量(t·年-1)、單船全汛捕撈量(kg·年-1)數(shù)據(jù)來源于上海市漁政監(jiān)督管理處。選取距離長江口最近的長江干流水文監(jiān)測站大通站的年輸沙量(億t·年-1)和徑流量(億m3·年-1)數(shù)據(jù)用于分析輸沙量和徑流量對鳳鱭汛期捕撈量的影響,大通站1982—2014年年輸沙量和徑流量數(shù)據(jù)來源于水利部長江水利委員會發(fā)布的長江泥沙公報[30]。長江口水溫(℃)數(shù)據(jù)來源于東北亞海洋觀測系統(tǒng)(NEAR-GOOS)海洋預報系統(tǒng)網(wǎng)[31],網(wǎng)站提供時間跨度為1982年1月—2014年12月長江口(30°40′~31°40′N、 122°00′~122°20′E)海表水溫,分辨率為月。本研究選取鳳鱭捕撈期4—6月水溫數(shù)據(jù)。

    1.2 數(shù)據(jù)處理和模型建立

    采用Origin 8.6軟件對數(shù)據(jù)進行處理和繪圖,并用互相關函數(shù)(cross-correlation function)對鳳鱭捕撈量和環(huán)境因子的水溫、徑流和輸沙量進行互相關分析(cross-correlation analysis)。

    采用格局轉(zhuǎn)變貫序t檢驗(STARS)的方法計算鳳鱭捕撈量、捕撈船數(shù)和環(huán)境因子水溫、徑流量、輸沙量的年際格局轉(zhuǎn)變指數(shù)(regime shift index,RSI)。t0時刻的RSI值公式為[32]:

    RSI=ttest(t0)=

    鳳鱭捕撈量與捕撈船數(shù)及環(huán)境因子間的關系進一步采用廣義可加模型(GAM)擬合,將鳳鱭總捕撈量和單船全汛捕撈量分別作為模型響應變量,以捕撈船數(shù)、水溫、徑流和輸沙量作為解釋變量建模。由于4—6月水溫之間互相顯著相關,會造成模型多重共線性,同時考慮到5月為主汛期,月捕撈量比例高達61.2%[10],因此,選用5月水溫數(shù)據(jù)作為解釋因子用于建模。本研究中,采用逐步回歸的方法篩選模型變量,即將解釋變量逐一輸入到模型中,根據(jù)F檢驗評估變量的顯著性,若P<0.05時,保留此變量。依次增加變量,最后根據(jù)模型赤池信息量準則(Akaike information criterion,AIC)越小即表明模型擬合效果越佳的原則[33-34],確定最終模型。GAM模型擬合和繪圖均采用R 3.1.1軟件中的mgcv包實現(xiàn)[33]。

    2 結果與分析

    2.1 捕撈船數(shù)和捕撈量的年際格局轉(zhuǎn)變

    研究結果顯示,1982—2014年長江口鳳鱭汛期捕撈船數(shù)變幅為59~330艘·年-1,表現(xiàn)為先增后減,總體呈下降趨勢。而總捕撈量和單船全汛捕撈量變幅分別為7.9~3 252.0 t·年-1和130.0~12 945.5 kg·年-1,均呈下降趨勢。利用STARS研究1982—2014年中長江口鳳鱭汛期捕撈船數(shù)、總捕撈量和單船全汛捕撈量格局轉(zhuǎn)變的年份,分析轉(zhuǎn)變前后特征變化(圖1~圖3)。

    結果顯示,捕撈船數(shù)分別在1990年、2000年和2006年發(fā)生3次格局轉(zhuǎn)變,對應RSI值分別為0.992、-0.412和-0.763。在1990年格局轉(zhuǎn)變前,捕撈船數(shù)保持在平均191艘·年-1,1990年后捕撈船數(shù)大幅增加,上升至平均298艘·年-1,增幅達55.5%;2000年之后捕撈船數(shù)開始下降,平均為207艘·年-1;至2006年捕撈船數(shù)進一步大幅減少,平均為86艘·年-1。

    汛期總捕撈量分別在1997年和2007年發(fā)生2次格局轉(zhuǎn)變,對應RSI值分別為-0.580和-0.405,即2次轉(zhuǎn)變均為下降態(tài)勢。1997年格局轉(zhuǎn)變之前總捕撈量較高,平均高達1 844.15t·年-1,1997年之后總捕撈量降低,平均為863.88 t·年-1,較格局轉(zhuǎn)變前下降了53.2%;至2007年第二次格局轉(zhuǎn)變后,總捕撈量銳減至75.6 t·年-1,幾不成汛。

    圖1 長江口1982—2014年鳳鱭汛期捕撈船數(shù)的年際格局轉(zhuǎn)變(a)和格局轉(zhuǎn)變指數(shù)(b)Fig.1 Annual variations, regime shift (a) and regime shift index (b)of Coilia mystus fishing boat numbers in the Yangtze River Estuary during 1982—2014

    圖2 長江口1982—2014年鳳鱭汛期總捕撈量的年際格局轉(zhuǎn)變(a)和格局轉(zhuǎn)變指數(shù)(b)Fig.2 Annual variations, regime shift (a) and regime shift index (b)of Coilia mystus total catch in the Yangtze River Estuary in fishing season during 1982—2014

    圖3 長江口1982—2014年鳳鱭汛期單船全汛捕撈量的年際格局轉(zhuǎn)變(a)和格局轉(zhuǎn)變指數(shù)(b)Fig.3 Annual variations, regime shift (a) and regime shift index (b)of Coilia mystus CPUE in the Yangtze River Estuary in fishing season during 1982—2014

    1982—2014年中,鳳鱭單船全汛捕撈量亦經(jīng)歷2次下降的年際格局轉(zhuǎn)變,轉(zhuǎn)變年份為1990年和2007年,對應RSI值分別為-0.806和-0.371。1990年之前,單船全汛捕撈量平均為10 460.1 kg·年-1;格局轉(zhuǎn)變后減少至4 711.6 kg·年-1,下降了55.0%;到2007年后進一步減少為720.8 kg·年-1,降幅達84.7%(圖3)。

    綜上可見,在1990年年際格局轉(zhuǎn)變后,單船全汛捕撈量下降了55.0%,相應捕撈船數(shù)提高了55.5%,延遲了總捕撈量的下降,甚至在1995年出現(xiàn)了曇花一現(xiàn)的總捕撈量峰值,但至1997年后總捕撈量大幅下滑。捕撈強度的提高,僅短期提升總捕撈量。至2006—2007年格局轉(zhuǎn)變后,汛期捕撈船數(shù)、總捕撈量和單船全汛捕撈量均大幅減少。

    2.2 環(huán)境因子影響

    2.2.1 長江口水溫

    由圖4可知,長江口4月、5月和6月水溫在1982—2014年33年內(nèi)總體均呈增高趨勢。水溫在1997格局轉(zhuǎn)變,RSI系數(shù)分別為0.465、0.498和0.404。4月、5月和6月平均水溫分別為11.3℃、13.7℃和17.5℃;1997年后,4月、5月和6月平均水溫分別增至13.3℃、16.0℃和19.7℃?;ハ嚓P分析結果顯示,4月、5月和6月長江口水溫均與鳳鱭總捕撈量及單船全汛捕撈量呈顯著負相關關系(P<0.05)(圖5,圖6)。

    圖4 長江口1982—2014年4—6月水溫的年際格局轉(zhuǎn)變(a)和格局轉(zhuǎn)變指數(shù)(b)Fig.4 Annual variations, regime shift (a) and regime shift index (b)of the surface water temperature in the Yangtze River Estuary from April to June during 1982—2014

    2.2.2 徑流量和輸沙量

    1982—2014年大通站徑流量和輸沙量變化范圍分別為6 671~12 490億m3·年-1和0.72~5.05億t·年-1,總體均呈下降趨勢。由圖7、圖8可見,徑流量和輸沙量均在2006年發(fā)生下降的年際格局轉(zhuǎn)變,對應RSI分別為-0.110和-0.862。格局轉(zhuǎn)變前,平均徑流量和輸沙量分別為9 217億m3·年-1和3.44億t·年-1,2006年后,降低為8 268億m3·年-1和1.24億t·年-1,較之前分別下降10.4%和63.9%。

    互相關分析結果顯示,徑流量與鳳鱭總捕撈量、單船全汛捕撈量均無顯著相關關系(P>0.05)(圖9)。但輸沙量與鳳鱭總捕撈量、單船全汛捕撈量均呈顯著正相關關系(P<0.05),且時滯年份為0時,正相關系數(shù)最大(圖10)。

    圖5 長江口1982—2014年鳳鱭汛期總捕撈量與長江口4—6月水溫的互相關分析Fig.5 Cross-correlation analysis between Coilia mystus total catch in fishing season and surface water temperature in the Yangtze River Estuary from April to June during 1982—2014注:虛線代表P=0.05顯著水平Note: The dashed line indicates the significance level (P=0.05)

    圖6 長江口1982—2014年鳳鱭單船全汛捕撈量與長江口4—6月水溫的互相關分析Fig.6 Cross-correlation analysis between Coilia mystus CPUE and surface water temperature in the Yangtze River Estuary from April to June during 1982—2014注:虛線代表P=0.05顯著水平Note: The dashed line indicates the significance level (P=0.05)

    圖7 長江大通水文站1982—2014年徑流量年際格局轉(zhuǎn)變(a)和格局轉(zhuǎn)變指數(shù)(b)Fig.7 Annual variations, regime shift(a) and regime shift index (b)of water discharge at Datong Hydrological Station during 1982—2014

    圖8 長江大通水文站1982—2014年輸沙量年際格局轉(zhuǎn)變(a)和格局轉(zhuǎn)變指數(shù)(b)Fig.8 Annual variations, regime shift (a) and regime shift index (b)of sediment discharge at Datong Hydrological Station during 1982—2014

    圖9 長江口1982—2014年鳳鱭汛期總捕撈量(a)和單船全汛捕撈量(b)與大通站徑流量的互相關分析Fig.9 Cross-correlation analysis between Coilia mystus total catch in fishing season and water discharge (a), and between Coilia mystus CPUE and water discharge (b) at Datong Hydrological Station during 1982—2014注:虛線代表P=0.05顯著水平Note: The dashed line indicates the significance level (P=0.05)

    2.3 GAM模型結果

    GAM模型在構建過程中需要逐步加入各個解釋變量及其交互項,根據(jù)F檢驗和AIC值來選擇最優(yōu)模型,GAM模型的最終選擇結果見表1和表2。模型對長江口鳳鱭汛期總捕撈量和單船全汛捕撈量的總偏差解釋率分別為97.4%和90.2%,表明2個模型均較好擬合各解釋變量與捕撈量之間的關系。

    F檢驗表明,徑流量對鳳鱭總捕撈量影響不顯著(P>0.05),而捕撈船數(shù)、水溫和輸沙量對鳳鱭總捕撈量具有顯著影響(P<0.05)(表1)。其中,捕撈船數(shù)對鳳鱭總捕撈量影響最大,呈非線性正相關關系(圖11),偏差解釋率高達89.2%,而環(huán)境因子中的水溫和輸沙量對鳳鱭總捕撈量亦存在影響,偏差解釋率分別為4.9%和3.3%(表1)。

    GAM模型還顯示,對鳳鱭單船全汛捕撈量影響最大的因子為輸沙量,偏差解釋率高達55.2%,隨著輸沙量的增加,單船全汛捕撈量隨之增大(圖12-a)。其次為捕撈船數(shù),偏差解釋率為33.5%,二者亦呈非線性的正相關關系。捕撈船數(shù)在50~150艘·年-1區(qū)間內(nèi),單船捕撈量隨著捕撈船數(shù)的增加而增大;而在捕撈船數(shù)超過150艘·年-1后,單船捕撈量隨之增加不顯著(圖12-b)。水溫的影響最弱,偏差解釋率為1.5%(表2)。水溫與鳳鱭單船全汛捕撈量呈線性負相關,隨著水溫升高,單船全汛捕撈量隨之降低(圖12-c)。

    圖10 長江口1982—2014年鳳鱭汛期總捕撈量(a)和單船全汛捕撈量(b)與大通站輸沙量的互相關分析Fig.10 Cross-correlation analysis between Coilia mystus total catch in fishing season and sediment discharge (a),and between Coilia mystus CPUE and sediment discharge (b) at Datong Hydrological Station during 1982—2014注:虛線代表P=0.05顯著水平Note: The dashed line indicates the significance level (P=0.05)

    圖11 廣義可加模型(GAM)所揭示的顯著影響長江口鳳鱭汛期總捕撈量變化的因子:(a)捕撈船數(shù), (b)水溫, (c)輸沙量 Fig.11 Generalized additive model (GAM) derived effects of fishing boat number(a), water temperature(b), sediment discharge(c), on total catch of Coilia mystus in the Yangtze River Estuary during 1982—2014注:S (因子)為光滑樣條函數(shù)的擬合值,表示其對鳳鱭總捕撈量的影響。實線表示鳳鱭總捕撈量的期望值,上下兩側(cè)虛線表示方程的95%置信區(qū)間Note: S(factor) is the fitted value of smoothing spline functions, representing their impacts on total catch.The solid lines represent the expected values of total catch; the dotted lines represent the 95% confidence intervals of equations

    模型因子Model factors殘差自由度df殘差偏差Residual deviance偏差變化量Deviance variation累計解釋偏差/%Accumulation of deviance explainedPAICNull32.0097.22+S(捕撈船數(shù))28.1810.5086.7289.2<0.00166.29+S(水溫)22.795.754.7594.10.002 956.47+S(年輸沙量)14.812.453.3097.40.012 642.23

    表2 長江口鳳鱭單船全汛捕撈量與相關因子GAM檢驗Tab.2 Significance of explanatory variables effects on Coilia mystus CPUE in the Yangtze River Estuary during 1982—2014 as determined by GAM

    圖12 廣義可加模型(GAM)所揭示的顯著影響長江口鳳鱭單船全汛捕撈量變化的因子:(a)年輸沙量, (b)捕撈船數(shù), (c)水溫Fig.12 Generalized additive model (GAM) derived effects of sediment discharge(a), fishing boat number(b),water temperature(c) on Coilia mystus CPUE in the Yangtze River Estuary during 1982—2014注:S (因子)為光滑樣條函數(shù)的擬合值,表示其對單船捕撈量的影響。實線表示單船全汛捕撈量的期望值,上下兩側(cè)虛線表示方程的95%置信區(qū)間Note: S(factor) is the fitted value of smoothing spline functions, representing their impacts on CPUE. The solid lines represent the expected values of CPUE; the dotted lines represent the 95% confidence intervals of equations

    3 討論

    1984—2014年長江口鳳鱭汛期總捕撈量和單船全汛捕撈量變幅分別為7.9~3 252.0 t·年-1和130.0~12 945.5 kg·年-1,均呈下降趨勢。張國祥和華家棟[4]研究指出,自20世紀70年代起長江口鳳鱭捕撈強度持續(xù)多年強化,使得鳳鱭的補充群體逐年減少,至80年代后鳳鱭資源量已明顯下降。在本研究中,1990年之前,捕撈船數(shù)保持在平均191艘·年-1,單船全汛捕撈量平均高達10 460.1 kg·年-1,捕撈量隨捕撈船數(shù)波動響應明顯;1990年后,捕撈船數(shù)增至207艘·年-1,單船全汛捕撈量反而減少至4 711.6 kg·年-1,期間僅在1995年出現(xiàn)“曇花一現(xiàn)”的高產(chǎn),捕撈量對捕撈船數(shù)的響應大幅減弱。STARS分析結果顯示,作為鳳鱭資源量波動指示因子之一的單船全汛捕撈量在1990年和2007年發(fā)生3次下降的年際格局轉(zhuǎn)變,其下降即預示鳳鱭資源的衰退。究其原因,1990年第一次格局轉(zhuǎn)變可能是由于過度捕撈最終導致鳳鱭種群調(diào)節(jié)機制被破壞,鳳鱭種群恢復能力已逐漸無法適應過高的捕撈強度,故而造成資源的衰退。

    STARS分析還顯示,捕撈船數(shù)、大通水文站徑流量和輸沙量都在2006年發(fā)生格局轉(zhuǎn)變,這與鳳鱭單船汛期捕撈量的第二次格局轉(zhuǎn)變時間節(jié)點基本吻合。長江口素有來水來沙量大的特點,1986年前通過大通水文站輸入長江口的平均年輸沙總量高達4.64億t·年-1,受溫度、鹽度、潮流、徑流等因素的影響,長江口門附近形成最大渾濁帶,長江口的高懸沙含量及復雜的行為動態(tài)在世界各河口也是非常少見的[35]。研究顯示,近三十年來,通過大通水文站的輸沙量呈波動下降態(tài)勢,三峽建庫后輸沙量降低更為顯著,2006—2014年長江大通站輸沙量年均僅為1.24億t,較之前下降了63.9%。其中,在不考慮沿程沖刷恢復的條件下,中等粒徑顆粒的泥沙減少更為明顯[9]。輸沙量的減少使得河口及其鄰近海岸動力作用和泥沙沉積過程發(fā)生變化,不可避免地對長江下游及河口的生態(tài)環(huán)境產(chǎn)生顯著影響,從而影響河口魚類及浮游生物的空間格局[36-38]。同時,由于入海泥沙量的減少所引起的河口水域海水透明度的增加,可能使得浮游植物高生產(chǎn)力區(qū)向海岸擴展,許多魚類的產(chǎn)卵索餌場位置亦會有相應的變化[39]。GAM模型亦顯示,輸沙量與鳳鱭單船全汛捕撈量存在顯著正相關關系,綜上推測,長江口輸沙量減少所造成的水文條件改變,可能對鳳鱭資源產(chǎn)生了負面影響。

    研究期間,通過大通站的徑流量在2006年后較之前下降10.4%。徑流量的降低可能主要受三峽工程蓄水的影響。DRINKWATER 和FRANK[38]曾提出河口徑流的變動可引起溫度、鹽度、濁度和溶解氧的變化,這些效應最終都影響到河口生物的豐度和分布類型。董方勇[40]認為水利工程造成河口區(qū)4—7月的徑流量減少,使得河口海水上溯的鹽度隨之發(fā)生變化,從而可能造成鳳鱭產(chǎn)卵場西移以及生殖洄游時間推遲。管衛(wèi)兵等[8]研究亦顯示,徑流量的減少對鳳鱭的生長和生殖具有較強的負面影響。劉凱等[10]報道,三峽蓄水后長江口鳳鱭生物學參數(shù)和汛期繁殖力均呈下降趨勢,絕對懷卵量和相對懷卵量較蓄水前分別下降17.38%和16.97%、全長和體質(zhì)量較蓄水前分別下降4.58%和3.55%。本研究中,GAM模型和互相關分析均顯示,徑流量對長江口鳳鱭捕撈量的影響未達顯著水平,但由于徑流時間格局的轉(zhuǎn)變至今尚不足10年,因此,仍然需要長期監(jiān)測和謹慎評估徑流量減少對鳳鱭乃至其他魚類的影響。

    此外,GAM模型還顯示長江口汛期水溫亦是影響長江口鳳鱭總捕撈量的重要環(huán)境因子之一。由于江海洄游型魚類在生殖洄游過程中一般并不攝食,依靠體內(nèi)積累的營養(yǎng)物質(zhì)完成溯河洄游的運動、代謝和生殖等過程,因此在生殖洄游開始前需要積累大量的脂肪[41]。鳳鱭這類短途溯河繁殖的洄游魚類,其卵巢粗脂肪含量達60%以上、肌肉粗脂肪含量達1.92%以上(鮮重),才開始生殖洄游[21]。管衛(wèi)兵等[8]的研究表明,水溫升高年份的鳳鱭生殖群體的豐滿度和性腺重量反而下降。可能的解釋是,水溫的升高使得魚類代謝增強、體內(nèi)營養(yǎng)物質(zhì)(粗脂肪等)的積累減緩[42-43],從而影響鳳鱭群體的成熟、集群和洄游過程,進而可能影響長江口鳳鱭資源量。

    除環(huán)境因子以外,長江口鰻苗網(wǎng)的兼捕對鳳鱭補充群體損害也極為嚴重[44],自2007年起,長江口鳳鱭的魚卵和仔稚魚豐度呈現(xiàn)大幅降低趨勢[45]。環(huán)境因子和捕撈因素的綜合影響可能是造成鳳鱭單船汛期捕撈量的第二次格局轉(zhuǎn)變的原因。但鳳鱭的捕撈量變化不能簡單地歸因于一種或幾種原因,在這一長期變化過程中更多的是生態(tài)因素與人為擾動相互作用的結果,關于更多影響因素和影響機制的研究也亟待展開。

    猜你喜歡
    輸沙量長江口年際
    北緯30°中層頂區(qū)域鈉與鐵原子層的結構和年際變化
    遼河干流遼中站河道沖淤變化趨勢研究
    氣候變化和人類活動對祖厲河輸沙量變化的影響分析
    頻率曲線在荊江三口輸沙量計算中的應用
    長江口鳳鱭繁殖群體的年齡結構和生長特性
    長江口及其鄰近水域仔稚魚種類組成及分布特征
    亞洲夏季風的年際和年代際變化及其未來預測
    與北大西洋接壤的北極海冰和年際氣候變化
    關于遼河干流河道沖淤量沿程變化規(guī)律探討
    地下水(2015年5期)2015-12-02 02:44:26
    基于M-K法對圖們江干流含沙量年際變化的分析
    新久久久久国产一级毛片| 一级毛片黄色毛片免费观看视频| 国产真实伦视频高清在线观看| 最近最新中文字幕免费大全7| 丰满迷人的少妇在线观看| 日本欧美国产在线视频| 久久久久久久久久久丰满| 久久av网站| 欧美精品人与动牲交sv欧美| 国产精品欧美亚洲77777| 午夜老司机福利剧场| 亚洲精品国产av蜜桃| 国产伦在线观看视频一区| 大话2 男鬼变身卡| 精品一区二区免费观看| 美女福利国产在线| 九九爱精品视频在线观看| 欧美+日韩+精品| 日韩人妻高清精品专区| 黑人猛操日本美女一级片| 久久午夜福利片| 久热久热在线精品观看| 人人妻人人看人人澡| 欧美精品人与动牲交sv欧美| 久久午夜综合久久蜜桃| 建设人人有责人人尽责人人享有的| 秋霞伦理黄片| 亚洲精品久久久久久婷婷小说| 免费久久久久久久精品成人欧美视频 | 久久久久久久大尺度免费视频| 成年女人在线观看亚洲视频| 国产免费福利视频在线观看| 少妇的逼好多水| 一本久久精品| 国产色婷婷99| 久久99一区二区三区| av福利片在线| 国产黄片视频在线免费观看| 内地一区二区视频在线| h日本视频在线播放| 国产精品无大码| 国产日韩一区二区三区精品不卡 | 韩国高清视频一区二区三区| 国产伦在线观看视频一区| a级毛片在线看网站| 一边亲一边摸免费视频| 国产中年淑女户外野战色| 国产伦精品一区二区三区四那| 我的老师免费观看完整版| 黑人高潮一二区| 欧美 亚洲 国产 日韩一| 国产日韩欧美亚洲二区| 久久国产精品大桥未久av | 日韩电影二区| 国产黄色免费在线视频| 久久午夜综合久久蜜桃| 韩国高清视频一区二区三区| 久久久久久久久久久免费av| 亚洲国产av新网站| 久久国内精品自在自线图片| 午夜免费观看性视频| 中文天堂在线官网| 视频中文字幕在线观看| 欧美成人午夜免费资源| 亚洲第一区二区三区不卡| 日韩大片免费观看网站| 国产精品久久久久久精品电影小说| 免费少妇av软件| 天天操日日干夜夜撸| 久久久国产一区二区| 国产精品嫩草影院av在线观看| 观看免费一级毛片| 男女免费视频国产| 久久久久人妻精品一区果冻| 亚洲欧洲国产日韩| 亚洲婷婷狠狠爱综合网| 国产色婷婷99| 国内少妇人妻偷人精品xxx网站| 一级a做视频免费观看| a 毛片基地| 久久久精品免费免费高清| 久久热精品热| 日本av手机在线免费观看| 一级av片app| 伊人久久精品亚洲午夜| 日日撸夜夜添| 国产精品99久久99久久久不卡 | 国产亚洲5aaaaa淫片| 成人18禁高潮啪啪吃奶动态图 | av在线老鸭窝| a级毛片在线看网站| 亚洲欧洲精品一区二区精品久久久 | 插阴视频在线观看视频| 精品亚洲成国产av| 国产极品天堂在线| 卡戴珊不雅视频在线播放| 蜜桃久久精品国产亚洲av| 久久久国产精品麻豆| 我的女老师完整版在线观看| 国产成人aa在线观看| 国产午夜精品一二区理论片| 观看av在线不卡| 一二三四中文在线观看免费高清| 自线自在国产av| 国产高清不卡午夜福利| 色吧在线观看| 少妇人妻精品综合一区二区| 九草在线视频观看| 久久久久久久久大av| 中文字幕精品免费在线观看视频 | 成年av动漫网址| 波野结衣二区三区在线| 久久久国产一区二区| 久久人人爽人人片av| 精品人妻熟女毛片av久久网站| 国产一区二区在线观看日韩| 亚洲va在线va天堂va国产| 国产一区有黄有色的免费视频| 日本黄大片高清| av一本久久久久| 在线观看www视频免费| 亚洲国产日韩一区二区| 亚洲av福利一区| 熟女人妻精品中文字幕| 伦理电影大哥的女人| 国产一区二区在线观看av| 日韩免费高清中文字幕av| 在线观看免费视频网站a站| 婷婷色麻豆天堂久久| 欧美日韩视频高清一区二区三区二| 国产精品久久久久久久久免| 制服丝袜香蕉在线| 精品少妇久久久久久888优播| 国产91av在线免费观看| 夫妻性生交免费视频一级片| 国产精品三级大全| 性色av一级| 国产在线免费精品| 少妇的逼好多水| 啦啦啦啦在线视频资源| av福利片在线观看| 激情五月婷婷亚洲| 国产精品伦人一区二区| 日本免费在线观看一区| 丰满迷人的少妇在线观看| 亚洲美女搞黄在线观看| 美女大奶头黄色视频| 青春草亚洲视频在线观看| 特大巨黑吊av在线直播| 三级国产精品欧美在线观看| 天美传媒精品一区二区| 一级爰片在线观看| av在线观看视频网站免费| 国精品久久久久久国模美| 新久久久久国产一级毛片| 22中文网久久字幕| 久久国产亚洲av麻豆专区| 丁香六月天网| 国产亚洲5aaaaa淫片| 国产成人a∨麻豆精品| 久久久久久久亚洲中文字幕| 欧美亚洲 丝袜 人妻 在线| 久久韩国三级中文字幕| 久热这里只有精品99| 日本av免费视频播放| av卡一久久| 国产淫语在线视频| 99热国产这里只有精品6| 成人漫画全彩无遮挡| 观看美女的网站| 99热国产这里只有精品6| 又粗又硬又长又爽又黄的视频| 18禁在线播放成人免费| 狂野欧美白嫩少妇大欣赏| a级一级毛片免费在线观看| 在线天堂最新版资源| 桃花免费在线播放| 久久久久久久精品精品| 国产精品99久久99久久久不卡 | 欧美日韩一区二区视频在线观看视频在线| 国产精品偷伦视频观看了| 男女国产视频网站| av女优亚洲男人天堂| 精品亚洲成a人片在线观看| 久久女婷五月综合色啪小说| 欧美日韩精品成人综合77777| 久久久久久久精品精品| 在线观看www视频免费| 免费看不卡的av| 日韩av不卡免费在线播放| 欧美另类一区| 乱系列少妇在线播放| 男人添女人高潮全过程视频| 久久精品国产自在天天线| 久久国内精品自在自线图片| 看十八女毛片水多多多| 搡老乐熟女国产| 高清在线视频一区二区三区| 在线观看免费日韩欧美大片 | 乱码一卡2卡4卡精品| 日韩一本色道免费dvd| 国精品久久久久久国模美| 黑人巨大精品欧美一区二区蜜桃 | 夜夜看夜夜爽夜夜摸| 18+在线观看网站| 亚洲人成网站在线观看播放| 尾随美女入室| 永久免费av网站大全| 美女国产视频在线观看| 夫妻午夜视频| 乱系列少妇在线播放| 欧美成人精品欧美一级黄| 18+在线观看网站| 在线观看av片永久免费下载| 日韩中文字幕视频在线看片| 午夜福利影视在线免费观看| 久久久国产欧美日韩av| 九九在线视频观看精品| 七月丁香在线播放| 美女cb高潮喷水在线观看| 一个人免费看片子| 99视频精品全部免费 在线| 亚洲av男天堂| 一级毛片我不卡| 精品人妻一区二区三区麻豆| 夜夜看夜夜爽夜夜摸| 欧美精品国产亚洲| 国产一级毛片在线| 国产极品天堂在线| 美女cb高潮喷水在线观看| 内射极品少妇av片p| 日韩中文字幕视频在线看片| 久久精品国产a三级三级三级| 汤姆久久久久久久影院中文字幕| 国产精品熟女久久久久浪| 国产探花极品一区二区| 少妇被粗大猛烈的视频| 最近2019中文字幕mv第一页| 国产有黄有色有爽视频| 久久99热这里只频精品6学生| 国产精品麻豆人妻色哟哟久久| 久久99蜜桃精品久久| 国产精品久久久久成人av| 欧美激情国产日韩精品一区| 多毛熟女@视频| 久久久久久久大尺度免费视频| 精品人妻偷拍中文字幕| 国产在视频线精品| 国产精品.久久久| 91久久精品国产一区二区三区| 老熟女久久久| 久久国产乱子免费精品| 国产亚洲一区二区精品| 婷婷色麻豆天堂久久| 久久国内精品自在自线图片| 大香蕉久久网| 国产高清国产精品国产三级| 久热这里只有精品99| 亚洲av.av天堂| 另类精品久久| 成人毛片60女人毛片免费| 91精品国产国语对白视频| 国产精品不卡视频一区二区| 久久久精品94久久精品| 日韩一区二区视频免费看| 亚洲国产av新网站| 熟妇人妻不卡中文字幕| 少妇被粗大猛烈的视频| 丝袜在线中文字幕| 国产黄频视频在线观看| 国产黄色视频一区二区在线观看| 三级国产精品欧美在线观看| 春色校园在线视频观看| 九九在线视频观看精品| 成人无遮挡网站| 欧美人与善性xxx| 国国产精品蜜臀av免费| 国产欧美日韩精品一区二区| 啦啦啦在线观看免费高清www| 亚洲精品一区蜜桃| 国产极品粉嫩免费观看在线 | 如何舔出高潮| 国产午夜精品久久久久久一区二区三区| 男人和女人高潮做爰伦理| 人人妻人人澡人人爽人人夜夜| 天堂8中文在线网| 国产欧美日韩精品一区二区| 国产男女内射视频| 精品视频人人做人人爽| 各种免费的搞黄视频| 欧美丝袜亚洲另类| 欧美日韩综合久久久久久| 免费不卡的大黄色大毛片视频在线观看| 如何舔出高潮| 久久久午夜欧美精品| 免费大片18禁| 久久久欧美国产精品| 久久久国产一区二区| 一个人免费看片子| 18禁裸乳无遮挡动漫免费视频| 亚洲国产成人一精品久久久| av福利片在线| 国产美女午夜福利| 永久免费av网站大全| 夜夜骑夜夜射夜夜干| 国产精品一区二区在线观看99| 91在线精品国自产拍蜜月| 国产伦理片在线播放av一区| 美女xxoo啪啪120秒动态图| 80岁老熟妇乱子伦牲交| 精品久久久久久久久av| 日韩欧美一区视频在线观看 | 日韩在线高清观看一区二区三区| 色哟哟·www| 黄色欧美视频在线观看| 国产国拍精品亚洲av在线观看| 精品人妻一区二区三区麻豆| 桃花免费在线播放| 免费黄频网站在线观看国产| 国产精品秋霞免费鲁丝片| 婷婷色av中文字幕| 91aial.com中文字幕在线观看| 国产精品99久久99久久久不卡 | 日韩av免费高清视频| 国产一区有黄有色的免费视频| 日本vs欧美在线观看视频 | 中文精品一卡2卡3卡4更新| 高清在线视频一区二区三区| 99久久中文字幕三级久久日本| 亚洲人成网站在线播| 如日韩欧美国产精品一区二区三区 | 久久 成人 亚洲| 精品少妇黑人巨大在线播放| av国产精品久久久久影院| 欧美激情国产日韩精品一区| 人人妻人人添人人爽欧美一区卜| 国产精品99久久久久久久久| 欧美日韩视频高清一区二区三区二| 国产精品人妻久久久久久| 男女边吃奶边做爰视频| 亚洲三级黄色毛片| 国产亚洲午夜精品一区二区久久| 国产欧美日韩精品一区二区| 国产伦精品一区二区三区视频9| 少妇人妻精品综合一区二区| 男女啪啪激烈高潮av片| 国产片特级美女逼逼视频| av天堂中文字幕网| 亚洲国产最新在线播放| 人体艺术视频欧美日本| 国产亚洲精品久久久com| 夫妻性生交免费视频一级片| 大片免费播放器 马上看| 亚洲久久久国产精品| 免费av中文字幕在线| 国产免费又黄又爽又色| 又粗又硬又长又爽又黄的视频| 日本欧美视频一区| 国国产精品蜜臀av免费| 久久亚洲国产成人精品v| 中国国产av一级| 亚洲人成网站在线观看播放| 亚洲国产精品一区三区| 不卡视频在线观看欧美| 乱码一卡2卡4卡精品| 精品久久久久久电影网| 99热这里只有精品一区| 日日撸夜夜添| 国产色婷婷99| 成人特级av手机在线观看| 国产午夜精品一二区理论片| 青春草亚洲视频在线观看| 日韩成人伦理影院| 欧美国产精品一级二级三级 | 国内揄拍国产精品人妻在线| av国产久精品久网站免费入址| 乱码一卡2卡4卡精品| 天堂中文最新版在线下载| 搡老乐熟女国产| 免费在线观看成人毛片| 国产高清不卡午夜福利| 97在线视频观看| 亚洲精品乱码久久久v下载方式| 中文天堂在线官网| 六月丁香七月| 中文资源天堂在线| 七月丁香在线播放| 一级av片app| 欧美精品人与动牲交sv欧美| 王馨瑶露胸无遮挡在线观看| 18+在线观看网站| 国产欧美另类精品又又久久亚洲欧美| 美女视频免费永久观看网站| 免费观看av网站的网址| 久久青草综合色| 国产欧美日韩精品一区二区| 精品国产一区二区久久| 最近的中文字幕免费完整| 国产有黄有色有爽视频| 成年人午夜在线观看视频| 男的添女的下面高潮视频| 久久久久国产精品人妻一区二区| 亚洲精品视频女| 丰满饥渴人妻一区二区三| 中文欧美无线码| 晚上一个人看的免费电影| 成人国产麻豆网| 精品亚洲成a人片在线观看| 日韩大片免费观看网站| av专区在线播放| 久久久久人妻精品一区果冻| 免费大片黄手机在线观看| 一二三四中文在线观看免费高清| 国产有黄有色有爽视频| 熟女人妻精品中文字幕| 欧美日韩在线观看h| 一级av片app| 亚洲欧洲国产日韩| 国产日韩欧美在线精品| 成人影院久久| 国产成人精品福利久久| 高清欧美精品videossex| 纵有疾风起免费观看全集完整版| 午夜激情福利司机影院| 婷婷色综合大香蕉| 久久精品国产a三级三级三级| 国产成人免费观看mmmm| 美女中出高潮动态图| 夫妻午夜视频| 精品人妻熟女毛片av久久网站| 91精品伊人久久大香线蕉| 免费av不卡在线播放| 亚洲四区av| 日韩av不卡免费在线播放| 色视频www国产| 精华霜和精华液先用哪个| 久久精品国产亚洲av天美| 桃花免费在线播放| 国产伦精品一区二区三区四那| 成人漫画全彩无遮挡| 岛国毛片在线播放| 久久国内精品自在自线图片| 婷婷色麻豆天堂久久| 久热久热在线精品观看| 一级毛片 在线播放| 免费人成在线观看视频色| 久久午夜综合久久蜜桃| 国产白丝娇喘喷水9色精品| 一本大道久久a久久精品| 久久热精品热| 国产一级毛片在线| 十分钟在线观看高清视频www | 美女国产视频在线观看| a级一级毛片免费在线观看| 久久人人爽av亚洲精品天堂| 久久久国产精品麻豆| 大香蕉97超碰在线| 只有这里有精品99| 午夜福利在线观看免费完整高清在| 国产亚洲精品久久久com| 男女啪啪激烈高潮av片| 亚洲av不卡在线观看| 少妇人妻 视频| 亚洲国产成人一精品久久久| 黑人高潮一二区| 午夜免费男女啪啪视频观看| av国产精品久久久久影院| 国产av精品麻豆| 精品人妻熟女av久视频| 美女主播在线视频| 国产亚洲av片在线观看秒播厂| 国产精品三级大全| 久久久久国产网址| 美女cb高潮喷水在线观看| 亚洲精品视频女| 免费看日本二区| 桃花免费在线播放| kizo精华| videossex国产| 天堂中文最新版在线下载| 国产精品久久久久成人av| 婷婷色麻豆天堂久久| 午夜精品国产一区二区电影| 精品国产一区二区久久| 99久久精品一区二区三区| 嫩草影院新地址| 人体艺术视频欧美日本| 亚洲精品成人av观看孕妇| 久久久久精品久久久久真实原创| 亚洲三级黄色毛片| 国产黄片视频在线免费观看| 天堂中文最新版在线下载| 国产精品国产三级国产专区5o| 少妇被粗大猛烈的视频| 中文字幕久久专区| 免费观看性生交大片5| 美女脱内裤让男人舔精品视频| 国产成人aa在线观看| 久久精品国产自在天天线| 在线观看一区二区三区激情| 青青草视频在线视频观看| 精品国产一区二区三区久久久樱花| 曰老女人黄片| 久久精品国产自在天天线| 夜夜爽夜夜爽视频| 精品少妇久久久久久888优播| 夜夜看夜夜爽夜夜摸| a级毛色黄片| 国产成人91sexporn| 少妇被粗大的猛进出69影院 | 亚洲国产av新网站| 少妇被粗大的猛进出69影院 | 桃花免费在线播放| 久久毛片免费看一区二区三区| 亚洲欧美日韩卡通动漫| 国产女主播在线喷水免费视频网站| 国产精品伦人一区二区| 中国三级夫妇交换| 麻豆成人av视频| 婷婷色综合www| 自拍偷自拍亚洲精品老妇| 涩涩av久久男人的天堂| 18禁在线播放成人免费| 十八禁网站网址无遮挡 | 国产毛片在线视频| 精品久久国产蜜桃| 国产精品国产av在线观看| 亚洲精品国产色婷婷电影| 久久婷婷青草| 亚洲综合色惰| 涩涩av久久男人的天堂| 一区二区av电影网| 国产乱人偷精品视频| 日日撸夜夜添| 丝袜喷水一区| 在线观看三级黄色| 精品一品国产午夜福利视频| 欧美少妇被猛烈插入视频| 18+在线观看网站| 人人澡人人妻人| 久久精品国产亚洲av涩爱| av在线播放精品| 成人免费观看视频高清| 精品国产国语对白av| 欧美日韩综合久久久久久| 免费观看的影片在线观看| 午夜福利网站1000一区二区三区| 麻豆成人午夜福利视频| 久久热精品热| 99热这里只有是精品在线观看| 国产一区二区三区av在线| 纯流量卡能插随身wifi吗| 少妇人妻久久综合中文| 乱码一卡2卡4卡精品| 一级片'在线观看视频| 亚洲丝袜综合中文字幕| 国产欧美另类精品又又久久亚洲欧美| 丝袜喷水一区| 国产精品一区二区三区四区免费观看| 搡老乐熟女国产| 色5月婷婷丁香| av在线播放精品| av天堂中文字幕网| 国产亚洲精品久久久com| 欧美精品人与动牲交sv欧美| 日韩一区二区三区影片| av免费观看日本| 久久精品国产a三级三级三级| 爱豆传媒免费全集在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲欧洲精品一区二区精品久久久 | 少妇丰满av| 观看免费一级毛片| 日本猛色少妇xxxxx猛交久久| 在线播放无遮挡| 最近中文字幕高清免费大全6| 最近中文字幕2019免费版| 国产无遮挡羞羞视频在线观看| 国产成人精品婷婷| 国产黄片美女视频| 国产视频内射| 国产男人的电影天堂91| 日韩熟女老妇一区二区性免费视频| 下体分泌物呈黄色| 亚洲精品日韩av片在线观看| 久久久久久人妻| 国产欧美另类精品又又久久亚洲欧美| 日韩视频在线欧美| 九色成人免费人妻av| 久久毛片免费看一区二区三区| 久久久久久久精品精品| 国产成人freesex在线| 免费观看无遮挡的男女| 亚洲成人av在线免费| 免费黄色在线免费观看| 女人精品久久久久毛片| 国产成人freesex在线| 最近的中文字幕免费完整| 丰满人妻一区二区三区视频av| 99re6热这里在线精品视频| 乱系列少妇在线播放| 一级毛片 在线播放| 久久久久国产网址| 韩国av在线不卡| 久久久久视频综合| 国产精品.久久久| 丰满饥渴人妻一区二区三| 亚洲内射少妇av| 欧美区成人在线视频| 亚洲成人一二三区av| 夫妻午夜视频| 肉色欧美久久久久久久蜜桃| 狠狠精品人妻久久久久久综合| 青春草亚洲视频在线观看| 纵有疾风起免费观看全集完整版|