• <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法對圖們江干流含沙量年際變化的分析
    国产精品日韩av在线免费观看| 搡老熟女国产l中国老女人| 99久久综合精品五月天人人| 国产在线精品亚洲第一网站| 色老头精品视频在线观看| 日韩精品青青久久久久久| 亚洲片人在线观看| 久久中文字幕人妻熟女| 国产精品久久视频播放| 亚洲avbb在线观看| 婷婷亚洲欧美| 亚洲自拍偷在线| 亚洲第一青青草原| 久久精品91无色码中文字幕| 校园春色视频在线观看| 免费在线观看视频国产中文字幕亚洲| 国产精品综合久久久久久久免费| 亚洲久久久国产精品| 夜夜爽天天搞| 亚洲人成伊人成综合网2020| 一级毛片高清免费大全| 少妇粗大呻吟视频| 波多野结衣高清作品| 一个人观看的视频www高清免费观看 | 亚洲人成伊人成综合网2020| 看片在线看免费视频| 精品一区二区三区视频在线观看免费| 精品国产美女av久久久久小说| 免费高清在线观看日韩| 日本一区二区免费在线视频| 国产在线精品亚洲第一网站| 999久久久精品免费观看国产| 亚洲精品久久成人aⅴ小说| 国产亚洲精品久久久久5区| 99久久国产精品久久久| 国产激情欧美一区二区| 色哟哟哟哟哟哟| 91麻豆精品激情在线观看国产| 啪啪无遮挡十八禁网站| 熟女少妇亚洲综合色aaa.| 欧美绝顶高潮抽搐喷水| 久久国产乱子伦精品免费另类| 久久久国产精品麻豆| 亚洲人成电影免费在线| 亚洲欧洲精品一区二区精品久久久| 欧美不卡视频在线免费观看 | 久久中文字幕一级| 97超级碰碰碰精品色视频在线观看| 麻豆一二三区av精品| 999精品在线视频| 在线免费观看的www视频| 麻豆av在线久日| √禁漫天堂资源中文www| 欧美激情久久久久久爽电影| 午夜福利欧美成人| 99久久无色码亚洲精品果冻| 日韩欧美在线二视频| 又黄又爽又免费观看的视频| 一a级毛片在线观看| 男女之事视频高清在线观看| 精品久久蜜臀av无| 亚洲国产精品合色在线| 日韩精品免费视频一区二区三区| 狠狠狠狠99中文字幕| 最新在线观看一区二区三区| 亚洲国产日韩欧美精品在线观看 | 亚洲精品久久成人aⅴ小说| 欧美性猛交黑人性爽| 可以在线观看的亚洲视频| 50天的宝宝边吃奶边哭怎么回事| 免费在线观看视频国产中文字幕亚洲| 国产精品1区2区在线观看.| 国产成人精品久久二区二区免费| 亚洲男人的天堂狠狠| 50天的宝宝边吃奶边哭怎么回事| 国产伦在线观看视频一区| 日本成人三级电影网站| 成人av一区二区三区在线看| 午夜福利在线在线| av在线播放免费不卡| 深夜精品福利| 日韩欧美免费精品| 黄片大片在线免费观看| 欧美色欧美亚洲另类二区| 亚洲男人的天堂狠狠| 美国免费a级毛片| 午夜福利在线在线| 免费看美女性在线毛片视频| 亚洲第一欧美日韩一区二区三区| 日韩精品中文字幕看吧| 黄色女人牲交| 制服诱惑二区| 99国产精品99久久久久| 在线视频色国产色| 69av精品久久久久久| 亚洲avbb在线观看| 亚洲成人免费电影在线观看| 波多野结衣高清作品| 日本一区二区免费在线视频| 精品国产国语对白av| 女人高潮潮喷娇喘18禁视频| 精品免费久久久久久久清纯| 法律面前人人平等表现在哪些方面| 777久久人妻少妇嫩草av网站| 久久精品aⅴ一区二区三区四区| 日本 欧美在线| 精品久久久久久久久久久久久 | 成在线人永久免费视频| 久久精品国产99精品国产亚洲性色| 熟妇人妻久久中文字幕3abv| 久久婷婷人人爽人人干人人爱| 这个男人来自地球电影免费观看| 91老司机精品| 91字幕亚洲| 亚洲狠狠婷婷综合久久图片| 久久草成人影院| 成人国产综合亚洲| 宅男免费午夜| 99精品久久久久人妻精品| 一区二区三区国产精品乱码| 亚洲自拍偷在线| 亚洲国产精品成人综合色| 国产一区在线观看成人免费| 在线观看午夜福利视频| 亚洲人成网站高清观看| 91字幕亚洲| 制服丝袜大香蕉在线| 成熟少妇高潮喷水视频| 亚洲人成网站在线播放欧美日韩| 国产一区二区在线av高清观看| 热99re8久久精品国产| 女人爽到高潮嗷嗷叫在线视频| 美女大奶头视频| 99久久无色码亚洲精品果冻| 一级毛片女人18水好多| 少妇粗大呻吟视频| 免费女性裸体啪啪无遮挡网站| 成人三级做爰电影| 一二三四在线观看免费中文在| 久久中文看片网| 日韩有码中文字幕| 成人永久免费在线观看视频| 午夜a级毛片| www.www免费av| 日日爽夜夜爽网站| 村上凉子中文字幕在线| 免费人成视频x8x8入口观看| 国产人伦9x9x在线观看| 国产爱豆传媒在线观看 | 久久久久久亚洲精品国产蜜桃av| 一二三四在线观看免费中文在| 亚洲熟妇中文字幕五十中出| 老司机午夜十八禁免费视频| 欧美激情久久久久久爽电影| 亚洲成人国产一区在线观看| 免费在线观看成人毛片| 曰老女人黄片| 黄片小视频在线播放| 日韩有码中文字幕| 亚洲一区中文字幕在线| 18禁国产床啪视频网站| 人人妻人人澡人人看| 在线十欧美十亚洲十日本专区| 草草在线视频免费看| 给我免费播放毛片高清在线观看| 国产亚洲精品综合一区在线观看 | 99国产精品一区二区蜜桃av| 最近最新中文字幕大全电影3 | 久久久久九九精品影院| 热re99久久国产66热| 波多野结衣av一区二区av| 亚洲国产欧美一区二区综合| 免费在线观看日本一区| 国产在线观看jvid| 青草久久国产| 欧美绝顶高潮抽搐喷水| 中文字幕人妻丝袜一区二区| 免费在线观看亚洲国产| 国产精品免费视频内射| 琪琪午夜伦伦电影理论片6080| 性欧美人与动物交配| 欧美精品亚洲一区二区| 变态另类丝袜制服| 亚洲欧美精品综合久久99| 国产亚洲av高清不卡| 亚洲精品粉嫩美女一区| 老司机福利观看| 亚洲三区欧美一区| 琪琪午夜伦伦电影理论片6080| 老司机福利观看| 麻豆久久精品国产亚洲av| 99国产综合亚洲精品| 国产成人啪精品午夜网站| 99久久久亚洲精品蜜臀av| 91av网站免费观看| 黑人欧美特级aaaaaa片| 村上凉子中文字幕在线| 男人的好看免费观看在线视频 | 久久精品人妻少妇| 夜夜爽天天搞| 午夜日韩欧美国产| 久久婷婷人人爽人人干人人爱| 国产精品美女特级片免费视频播放器 | 久99久视频精品免费| 国产伦一二天堂av在线观看| 一级a爱片免费观看的视频| 黄频高清免费视频| 一级作爱视频免费观看| 性欧美人与动物交配| 亚洲国产欧美网| 国产精品久久久久久亚洲av鲁大| 亚洲在线自拍视频| 色综合亚洲欧美另类图片| 午夜精品在线福利| 熟妇人妻久久中文字幕3abv| 男女床上黄色一级片免费看| 亚洲色图 男人天堂 中文字幕| 色综合婷婷激情| 国产av一区二区精品久久| 欧美色视频一区免费| 久久精品亚洲精品国产色婷小说| 免费电影在线观看免费观看| 操出白浆在线播放| 国产精品免费视频内射| 露出奶头的视频| 18禁黄网站禁片午夜丰满| 在线播放国产精品三级| 欧美性猛交黑人性爽| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩福利视频一区二区| 国产又黄又爽又无遮挡在线| 亚洲人成网站在线播放欧美日韩| 18禁国产床啪视频网站| 99riav亚洲国产免费| 国产精品永久免费网站| www日本在线高清视频| 久久国产精品男人的天堂亚洲| 亚洲无线在线观看| 成年免费大片在线观看| 国产伦在线观看视频一区| 欧美三级亚洲精品| 99热只有精品国产| 亚洲色图 男人天堂 中文字幕| 国产激情偷乱视频一区二区| 老司机靠b影院| 69av精品久久久久久| 国产1区2区3区精品| 亚洲片人在线观看| 女人被狂操c到高潮| 夜夜看夜夜爽夜夜摸| 99国产综合亚洲精品| 精品一区二区三区av网在线观看| 精品国产一区二区三区四区第35| bbb黄色大片| 一级黄色大片毛片| 日日干狠狠操夜夜爽| 欧美最黄视频在线播放免费| 国产一区在线观看成人免费| 啪啪无遮挡十八禁网站| 久久香蕉激情| 国产亚洲欧美在线一区二区| 久久久久久久精品吃奶| 最新在线观看一区二区三区| 国产精品一区二区免费欧美| 少妇 在线观看| 美女午夜性视频免费| 久久精品国产亚洲av高清一级| 婷婷丁香在线五月| 777久久人妻少妇嫩草av网站| 色综合亚洲欧美另类图片| 一本久久中文字幕| 午夜免费鲁丝| 亚洲国产欧美日韩在线播放| 久久中文看片网| 亚洲色图 男人天堂 中文字幕| 免费在线观看日本一区| 97碰自拍视频| 亚洲黑人精品在线| 美国免费a级毛片| 啦啦啦 在线观看视频| 日韩欧美在线二视频| 亚洲美女黄片视频| 国产精品美女特级片免费视频播放器 | 国产精品自产拍在线观看55亚洲| 午夜免费鲁丝| 亚洲第一青青草原| 久久久久精品国产欧美久久久| 国产野战对白在线观看| 中文字幕人妻熟女乱码| 亚洲国产高清在线一区二区三 | 热re99久久国产66热| 亚洲av美国av| 欧美大码av| 国内揄拍国产精品人妻在线 | 亚洲精品在线美女| 丝袜人妻中文字幕| 成年版毛片免费区| 精品免费久久久久久久清纯| 搡老熟女国产l中国老女人| 在线天堂中文资源库| 色尼玛亚洲综合影院| 波多野结衣av一区二区av| 三级毛片av免费| 亚洲一区二区三区色噜噜| 母亲3免费完整高清在线观看| 久久精品国产综合久久久| 欧美激情极品国产一区二区三区| 久久人妻av系列| 久热这里只有精品99| 精品久久蜜臀av无| 在线视频色国产色| 成人国语在线视频| 午夜福利在线在线| 精品无人区乱码1区二区| 老司机靠b影院| 黄片大片在线免费观看| 久久香蕉精品热| 国产成人精品久久二区二区91| 国产熟女xx| 成年版毛片免费区| 亚洲七黄色美女视频| 午夜日韩欧美国产| 91成年电影在线观看| 50天的宝宝边吃奶边哭怎么回事| 两个人免费观看高清视频| 国产成人欧美在线观看| 免费看a级黄色片| 国产成人av教育| 午夜视频精品福利| 最近在线观看免费完整版| 一区二区三区激情视频| 在线观看午夜福利视频| 精品不卡国产一区二区三区| 亚洲七黄色美女视频| 成人永久免费在线观看视频| 一个人免费在线观看的高清视频| 看片在线看免费视频| 亚洲激情在线av| 一卡2卡三卡四卡精品乱码亚洲| 在线视频色国产色| 精品欧美一区二区三区在线| 美女午夜性视频免费| www国产在线视频色| 欧美日本亚洲视频在线播放| 久久久精品国产亚洲av高清涩受| 国产精品综合久久久久久久免费| 精品欧美国产一区二区三| 亚洲人成网站在线播放欧美日韩| 亚洲精华国产精华精| 亚洲色图av天堂| 日日干狠狠操夜夜爽| 99国产精品一区二区三区| 欧美性猛交╳xxx乱大交人| 亚洲av片天天在线观看| 露出奶头的视频| 看黄色毛片网站| 色综合站精品国产| 中文字幕高清在线视频| 日韩视频一区二区在线观看| 在线国产一区二区在线| 俄罗斯特黄特色一大片| 欧美绝顶高潮抽搐喷水| 国产成人精品无人区| 亚洲国产精品合色在线| 亚洲国产欧美一区二区综合| 国产黄片美女视频| 18禁裸乳无遮挡免费网站照片 | 制服丝袜大香蕉在线| 亚洲午夜理论影院| 深夜精品福利| 日日爽夜夜爽网站| 免费高清在线观看日韩| 熟妇人妻久久中文字幕3abv| 制服丝袜大香蕉在线| 欧美绝顶高潮抽搐喷水| 亚洲国产精品成人综合色| 成人国产综合亚洲| 亚洲精华国产精华精| 亚洲专区国产一区二区| 99riav亚洲国产免费| 熟女电影av网| 欧美人与性动交α欧美精品济南到| 悠悠久久av| 欧美日韩一级在线毛片| 韩国精品一区二区三区| xxxwww97欧美| 一边摸一边做爽爽视频免费| 亚洲九九香蕉| 国产欧美日韩一区二区精品| 成在线人永久免费视频| 午夜老司机福利片| 欧美 亚洲 国产 日韩一| 亚洲成人国产一区在线观看| xxxwww97欧美| 亚洲欧美精品综合一区二区三区| 天堂动漫精品| 午夜福利成人在线免费观看| 黄色 视频免费看| 免费搜索国产男女视频| 免费一级毛片在线播放高清视频| 美女午夜性视频免费| 午夜激情福利司机影院| 一边摸一边做爽爽视频免费| 国产精品精品国产色婷婷| 亚洲人成网站高清观看| 丁香欧美五月| 亚洲人成电影免费在线| 亚洲专区国产一区二区| 少妇裸体淫交视频免费看高清 | 午夜免费激情av| 色播在线永久视频| 男女之事视频高清在线观看| 757午夜福利合集在线观看| 亚洲av熟女| 黄色成人免费大全| 国产精华一区二区三区| 亚洲 国产 在线| 国产成人一区二区三区免费视频网站| 最好的美女福利视频网| 91麻豆av在线| 国产一卡二卡三卡精品| 丁香六月欧美| 极品教师在线免费播放| 岛国在线观看网站| 国产精品久久久久久人妻精品电影| 亚洲精品一卡2卡三卡4卡5卡| 久久久久九九精品影院| 国产爱豆传媒在线观看 | 在线观看舔阴道视频| 日韩欧美免费精品| 国产精品98久久久久久宅男小说| 欧美性猛交╳xxx乱大交人| 亚洲国产欧美网| 搡老妇女老女人老熟妇| 女同久久另类99精品国产91| 亚洲五月婷婷丁香| 欧美一区二区精品小视频在线| 免费在线观看成人毛片| 久久狼人影院| 精品一区二区三区av网在线观看| 亚洲av电影在线进入| 美国免费a级毛片| 欧美在线一区亚洲| 精品国产超薄肉色丝袜足j| 真人一进一出gif抽搐免费| ponron亚洲| 最新在线观看一区二区三区| 午夜福利一区二区在线看| 久久国产精品人妻蜜桃| 国产精品日韩av在线免费观看| 国产伦在线观看视频一区| 国产精品久久久久久人妻精品电影| 一边摸一边做爽爽视频免费| 黄色片一级片一级黄色片| 美女大奶头视频| www.熟女人妻精品国产| 国产一区二区在线av高清观看| 国产亚洲精品一区二区www| 亚洲黑人精品在线| 亚洲精品中文字幕在线视频| 桃色一区二区三区在线观看| 久久99热这里只有精品18| 高潮久久久久久久久久久不卡| 一个人免费在线观看的高清视频| 亚洲午夜精品一区,二区,三区| 在线观看午夜福利视频| 亚洲全国av大片| 十分钟在线观看高清视频www| 日韩欧美一区视频在线观看| 精品一区二区三区av网在线观看| 听说在线观看完整版免费高清| 亚洲va日本ⅴa欧美va伊人久久| 最近最新免费中文字幕在线| 久久久久久九九精品二区国产 | 亚洲一区二区三区不卡视频| 啪啪无遮挡十八禁网站| 12—13女人毛片做爰片一| 欧美丝袜亚洲另类 | 国产高清视频在线播放一区| 国产单亲对白刺激| 黄片小视频在线播放| 热re99久久国产66热| 亚洲av熟女| 精品一区二区三区视频在线观看免费| 欧美色欧美亚洲另类二区| 久久精品夜夜夜夜夜久久蜜豆 | 91麻豆精品激情在线观看国产| 成人一区二区视频在线观看| 国产亚洲精品久久久久5区| 欧美国产精品va在线观看不卡| 国产高清有码在线观看视频 | 国产精品久久久久久亚洲av鲁大| 国产激情久久老熟女| 国产欧美日韩精品亚洲av| 国产伦一二天堂av在线观看| 99riav亚洲国产免费| 伊人久久大香线蕉亚洲五| 狂野欧美激情性xxxx| 少妇被粗大的猛进出69影院| 成人精品一区二区免费| 亚洲国产欧美一区二区综合| 亚洲专区中文字幕在线| 免费在线观看亚洲国产| 久久99热这里只有精品18| 国产精品野战在线观看| 免费人成视频x8x8入口观看| 午夜视频精品福利| 99久久综合精品五月天人人| 十八禁网站免费在线| 中文字幕精品免费在线观看视频| 禁无遮挡网站| 亚洲欧美日韩无卡精品| 精品久久久久久成人av| 19禁男女啪啪无遮挡网站| 欧美精品啪啪一区二区三区| 国产欧美日韩一区二区精品| 啦啦啦免费观看视频1| 我的亚洲天堂| 听说在线观看完整版免费高清| 国产精品1区2区在线观看.| 国产精品美女特级片免费视频播放器 | 人妻丰满熟妇av一区二区三区| 日韩欧美免费精品| 在线观看www视频免费| 十分钟在线观看高清视频www| 国产私拍福利视频在线观看| 日日爽夜夜爽网站| 男女视频在线观看网站免费 | videosex国产| 欧美乱色亚洲激情| 99riav亚洲国产免费| 精品人妻1区二区| 国产人伦9x9x在线观看| 国产精品精品国产色婷婷| 亚洲人成网站在线播放欧美日韩| 亚洲成a人片在线一区二区| 1024香蕉在线观看| 一本综合久久免费| 午夜免费激情av| 国产野战对白在线观看| 黄片播放在线免费| 亚洲欧美一区二区三区黑人| 99久久久亚洲精品蜜臀av| 成熟少妇高潮喷水视频| 国产精品九九99| 久久久久免费精品人妻一区二区 | 日本熟妇午夜| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩精品免费视频一区二区三区| 看片在线看免费视频| a级毛片在线看网站| 午夜福利在线在线| 九色国产91popny在线| 欧美日韩乱码在线| 国产成人欧美| 亚洲国产中文字幕在线视频| 成人亚洲精品一区在线观看| 国语自产精品视频在线第100页| 男人操女人黄网站| 久久国产精品男人的天堂亚洲| 成人一区二区视频在线观看| 亚洲国产欧美一区二区综合| 国产av在哪里看| 久久国产精品人妻蜜桃| 精品免费久久久久久久清纯| 少妇 在线观看| 午夜精品久久久久久毛片777| 97碰自拍视频| 欧美黄色淫秽网站| 久久久国产欧美日韩av| 成年人黄色毛片网站| 又紧又爽又黄一区二区| 成年人黄色毛片网站| 欧美一区二区精品小视频在线| 精品无人区乱码1区二区| 久久国产乱子伦精品免费另类| 熟女少妇亚洲综合色aaa.| 一a级毛片在线观看| 久久午夜综合久久蜜桃| 国内揄拍国产精品人妻在线 | 久久国产精品影院| tocl精华| 午夜福利在线在线| aaaaa片日本免费| 99re在线观看精品视频| 日韩欧美在线二视频| 国产精品二区激情视频| 国产激情偷乱视频一区二区| 精品国产乱码久久久久久男人| 成年版毛片免费区| 一本一本综合久久| 人人妻人人澡欧美一区二区| 99久久精品国产亚洲精品| 午夜a级毛片| 脱女人内裤的视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品一卡2卡三卡4卡5卡| 最近在线观看免费完整版| 久久久国产欧美日韩av| www.熟女人妻精品国产| 国产一区二区在线av高清观看| 国产伦在线观看视频一区| 这个男人来自地球电影免费观看| 久久久久久免费高清国产稀缺| 午夜激情福利司机影院| 久久精品aⅴ一区二区三区四区| 50天的宝宝边吃奶边哭怎么回事| 国产一区二区三区在线臀色熟女| 国产一级毛片七仙女欲春2 | 少妇被粗大的猛进出69影院| 韩国av一区二区三区四区| 亚洲色图 男人天堂 中文字幕| 两个人免费观看高清视频| 又黄又粗又硬又大视频| 一进一出抽搐动态|