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

    基于DSAS和FA的1959—2002年黃河三角洲海岸線演變規(guī)律及影響因素研究*

    2015-03-08 08:48:11杜國云李雪艷
    海洋與湖沼 2015年3期

    劉 鵬 王 慶 戰(zhàn) 超 王 昕 杜國云 李雪艷

    (魯東大學(xué)海岸研究所 煙臺 264025)

    海岸線是海岸帶環(huán)境的重要組成部分, 受全球氣候變化、海平面上升和流域—海域人類活動的影響,其變化具有復(fù)雜、敏感和多變的特性(莊振業(yè)等, 2003;諶艷珍等, 2010)。從海岸地貌學(xué)角度看, 海岸線變化是海岸帶沖淤演變的直接體現(xiàn), 是泥沙、動力等因素在特定相對海平面變化背景和海陸邊界條件下共同作用的結(jié)果(Fenster et al, 2001; 豐愛平等, 2006; 薛春汀, 2009)。在三角洲地區(qū), 因受入海河流和海洋的雙重動力作用, 其海岸線變化的速率和方向較一般海岸又有明顯差異(李安龍等, 2004; Cui et al, 2011;栗云召等, 2012)。另一方面, 海岸帶和流域日益強(qiáng)烈的人類活動, 特別是最近數(shù)十年來大規(guī)模的流域土地利用、水庫調(diào)蓄、灘涂(淺海)圍墾導(dǎo)致的入海水沙、海岸動力變化, 對海岸地貌及海岸線變化的影響也日趨劇烈(王穎等, 2011; 李蒙蒙等, 2013; 戰(zhàn)超等,2013)。

    作為全球著名的陸相淤積性弱潮三角洲, 最近數(shù)十年來黃河三角洲地區(qū)的河口、岸線及地貌變遷非常頻繁, 而且隨著黃河三角洲高效生態(tài)經(jīng)濟(jì)區(qū)被列入我國國家戰(zhàn)略規(guī)劃區(qū), 其城市規(guī)劃、產(chǎn)業(yè)結(jié)構(gòu)與布局、環(huán)境保護(hù)和資源利用, 對未來岸線變遷及海岸地貌沖淤趨勢預(yù)測的要求變得更高。鑒于前人對黃河口沉積動力過程、黃河入海懸沙擴(kuò)散、廢棄河口—亞三角洲及若干重要岸段的動力、地貌演變已有較多研究(Hu et al, 1998; Li et al, 2000; Li et al, 2001; Chu et al, 2006;Wang et al, 2006; Qiao et al, 2010; Cui et al, 2011), 為揭示最近50年尺度上三角洲岸線變化與黃河河口河道變遷、河流入海輸沙數(shù)量變化及海岸沉積動力過程的相互關(guān)系, 將黃河三角洲海岸作為一個統(tǒng)一的地貌體,研究其在數(shù)十年尺度上的海岸地貌沖淤演變機(jī)制與趨勢。本文基于黃河三角洲1959年至2002年之間的多期實(shí)測海岸線數(shù)據(jù), 利用數(shù)字岸線分析系統(tǒng)(DSAS)和岸線分形分析(FA), 對黃河三角洲岸線演變規(guī)律及影響因素進(jìn)行研究。

    圖1 黃河三角洲位置、范圍、不同時期海岸線及DSAS分析所用的基線和切線分布Fig.1 Location, range, and the coastlines in different periods of the Yellow River delta and the distributions of baseline and transects for DSAS (Digital Shoreline Analysis System) analysis

    1 研究區(qū)概況

    黃河三角洲瀕臨渤海灣南部和萊州灣西部, 系由古代、近代和現(xiàn)代三角洲組成的亞三角洲聯(lián)合體。其中, 近、現(xiàn)代三角洲形成于 1855年黃河奪大清河流路入渤海以來。本文所指黃河三角洲海岸線范圍,西北端起于灣灣溝附近, 東南端止于淄脈溝附近。最近 50多年來, 由于黃河下游河道變遷, 黃河口位置發(fā)生了多次遷移(圖1)。1960年黃河口自神仙溝口向東遷至老神仙溝口, 1964年向西遷到釣口, 1976年向東南遷到清水溝口, 1996年又向北遷到清水溝河口沙嘴北側(cè)的清8口。每次新的河口形成后, 新的河道迅速向海延伸, 形成大規(guī)模的河口沙嘴及亞三角洲, 而廢棄河口—亞三角洲及水下三角洲因泥沙消失而由淤進(jìn)轉(zhuǎn)為蝕退。與此同時, 河口附近的海岸線也隨之快速地變化。

    黃河三角洲為典型的弱潮海岸。沿岸潮汐主要受M2分潮控制, 在五號樁東北方向近岸海域存在一無潮點(diǎn), 其西為瀕臨渤海灣的三角洲北岸, 其南為瀕臨萊州灣的三角洲東岸。三角洲沿岸平均潮差約1m(0.6—1.3m), 大部分岸段為不正規(guī)半日潮, 僅神仙溝附近岸段為不正規(guī)全日潮。行水期河口潮流以平行于河道的往復(fù)流為主, 口門兩側(cè)海岸則為大致平行岸線的往復(fù)流(薛鴻超等, 1996)。根據(jù)1987年在三角洲東北岸外的黃河海港碼頭北端處(水深 7m)的觀測結(jié)果, 波高大于1.5m的波浪出現(xiàn)頻率為11.8%, 常浪向SSE-S、強(qiáng)浪向NE-NNE, 最大波高5.3m、波向NE、周期8.0s。

    黃河口為著名的多沙河口。但是, 河口輸沙的年內(nèi)分配、年際差異極大。據(jù)利津水文站 1950—2007年的觀測數(shù)據(jù)統(tǒng)計, 黃河口年輸沙量最大高達(dá) 20.98億t, 最小年份僅0.16億t。其中, 多年汛期平均的入海沙量占全年 84.4%, 年汛期輸沙量最大高達(dá) 19.26億t, 最小年份僅0.07億t。黃河泥沙來源為黃土高原,主要由粉砂、粘土構(gòu)成, 粒徑小于0.063mm的顆粒占94.2%。這些細(xì)粒泥沙在波浪作用下極易啟動但難以沉積, 其入海后主要被潮流搬運(yùn)、遷移, 為黃河三角洲海岸地貌形成演變提供了豐富的物質(zhì)基礎(chǔ)。

    2 數(shù)據(jù)及方法

    2.1 數(shù)據(jù)源及預(yù)處理

    由于最近10多年來黃河三角洲人工岸線迅速增加, 再考慮到歷史時期實(shí)測海岸地形、水深資料的限制, 本文選取測量于1959年、1968年、1974年、1980年、1984年、1992年、1996年、2000年、2002年的9期黃河三角洲附近海域水深圖和海圖為數(shù)據(jù)源。其中, 對各年份水深圖或海圖上的海岸線位置, 均利用相同年份成像的衛(wèi)星遙感影像進(jìn)行校正, 并補(bǔ)齊缺失海岸線數(shù)據(jù)的局部岸段岸線。文中所指海岸線均為大潮平均高潮線, 所用岸線數(shù)據(jù)詳情如表1所示。

    在定量分析前, 先對岸線數(shù)據(jù)進(jìn)行預(yù)處理, 其處理流程如下。首先, 將紙質(zhì)水深圖或海圖掃描數(shù)字化并得到 RGB格式圖像, 將其橢球體和基準(zhǔn)面設(shè)置為WGS84。其次, 根據(jù)圖上標(biāo)注的經(jīng)緯度選取約 30個控制點(diǎn)進(jìn)行配準(zhǔn)。然后, 提取矢量化并配準(zhǔn)后的水深圖或海圖上的岸線信息。最后, 根據(jù)海岸線變化速率計算需要, 對所得不同時期的岸線數(shù)據(jù)進(jìn)行空間疊加。在空間疊加前, 在ArcGIS10.1軟件的支持下, 將經(jīng)過預(yù)處理的水深圖或海圖統(tǒng)一轉(zhuǎn)換到墨卡托投影(UTM, Zone 50), 橢球體和基準(zhǔn)面均為WGS 84。

    表1 黃河三角洲歷史海岸線數(shù)據(jù)情況一覽表Tab.1 Historical coastline data source of the Yellow River delta

    2.2 數(shù)字岸線計算

    采用美國地質(zhì)調(diào)查局(USGS)研發(fā)的數(shù)字海岸線分析系統(tǒng)(DSAS, Digital Shoreline Analysis System),對黃河三角洲岸線進(jìn)行定量化分析。DSAS分析借助于 ArcGIS平臺, 通過回歸分析方法計算出海岸線的變化速率, 從而定量揭示大范圍、長時間序列的海岸線變遷歷史。目前, DSAS系統(tǒng)已在國內(nèi)外很多有關(guān)岸線變遷的研究中得到成功應(yīng)用, 并取得了有價值的研究成果(Jabaloy-Sánchez et al, 2010; Kuleli et al,2011; 張曉祥等, 2014; Jabaloy-Sánchez et al, 2014;Kaliraj et al, 2014)。

    利用 DSAS進(jìn)行數(shù)字海岸線分析的主要計算流程如下。首先, 結(jié)合已有的岸線數(shù)據(jù), 通過多次擬合、修改和調(diào)整, 選定了作為岸線數(shù)字計算基準(zhǔn)的基線(圖1b)。其次, 統(tǒng)一設(shè)置1500m的采樣間距, 在進(jìn)行多次擬合、修正和調(diào)整后, 生成了109條垂直于上述基線并與所有岸線相交的等間距切線(圖 1b)。然后,分別選取相鄰的兩期岸線, 利用DSAS系統(tǒng)對109個斷面進(jìn)行計算。

    根據(jù)岸線變遷及地貌沖淤分析需要, 計算了不同年份岸線長度、不同時段岸線年平均變化速率、終點(diǎn)變化速率以及前進(jìn)、后退、進(jìn)退平衡岸線所占比例等參數(shù), 用于對黃河三角洲海岸線變遷進(jìn)行定量分析研究。其中, 岸線長度、岸線年平均變化速率反映三角洲岸線的總體進(jìn)退情況, 而終點(diǎn)變化速率以及前進(jìn)、后退和平衡岸線所占比例則反映岸線變化的空間差異。

    岸線年平均變化速率的計算公式如下:

    式中, Ri,i+1是相鄰年份間的岸線年平均變化率(i=1—8);Si+1為第 i+1期海岸線與研究區(qū)陸域縱、豎邊界所圍區(qū)域的面積(i=1—8);Si為第i期海岸線與研究區(qū)陸域縱、豎邊界所圍區(qū)域的面積(i=1—8); Li為第i期海岸線的長度為第 i+1期與第i期海岸線年份數(shù)的差值(i=1—8)。

    終點(diǎn)變化速率(End Point Rate)的計算公式如下:

    式中, Ei,j是相鄰年份間沿某條切線m的岸線終點(diǎn)變化速率; dj為沿切線m第j期海岸線到基線的距離;di為沿切線m第i期海岸線到基線的距離; ΔYj,i為沿切線m第j期與第i期海岸線年份數(shù)的差值。

    前進(jìn)岸線、后退岸線和進(jìn)退平衡岸線所占比例的計算公式如下:

    式中, Pi,i+1是相鄰年份之間前進(jìn)、后退和進(jìn)退平衡岸線所占比例(i=1—8);Ti,i+1為參與計算的切線總數(shù)(i=1—8);ti,i+1為岸線處于前進(jìn)、后退和平衡的切線數(shù)量(i=1—8)。考慮到傳統(tǒng)紙質(zhì)海圖的極限精度為0.1mm、作業(yè)分辨率為0.2mm, 而本文所用紙質(zhì)海圖(水深圖)的最小比例尺為1 : 150000, 故以±30m/a作為劃分標(biāo)準(zhǔn), 若岸線變化速率大于 30m/a, 則為前進(jìn)岸線; 若岸線變化速率小于–30m/a, 則為后退岸線;若岸線變化速率介于–30—30m/a之間, 則確定為進(jìn)退平衡岸線。

    2.3 岸線分維數(shù)計算

    岸線分維數(shù)表示岸線平面輪廓形態(tài)的復(fù)雜程度,而岸線形態(tài)復(fù)雜程度又與地貌沖淤有關(guān), 即岸線分維數(shù)與海岸侵蝕強(qiáng)度呈正相關(guān)。本文用 Matlab基于盒子法通過編程計算不同時期的岸線分維數(shù)。計算流程如下。首先, 用 ArcGIS將岸線矢量數(shù)據(jù)轉(zhuǎn)換為單一圖層?xùn)鸥駭?shù)據(jù), 并進(jìn)行二值化運(yùn)算。其次, 生成覆蓋全部研究區(qū)的正方形格網(wǎng), 計算子格網(wǎng)邊長統(tǒng)一為r0時覆蓋全部岸線所需最少網(wǎng)格數(shù)。然后, 逐次減小r0并使之趨于0, 得到一系列離散的r0值。最后, 在對數(shù)坐標(biāo)系中對進(jìn)行線性回歸, 所得直線段斜率即為岸線分維數(shù)。如果某時期的岸線分維數(shù)變大, 則岸線變得曲折、復(fù)雜; 分維數(shù)變小, 則岸線變得平滑、簡單。分維數(shù)(D)的計算公式如下:

    圖2 1959年到2002年黃河三角洲岸線長度及其年平均變化速率Fig.2 Coastline length and annual average change rate of the Yellow River delta from 1959 to 2002

    3 研究結(jié)果

    3.1 岸線長度及其年平均速率變化

    從1959年到2002年, 以1984年為轉(zhuǎn)折點(diǎn), 黃河三角洲海岸線長度變化趨勢相反(圖2)。在1984年以前, 岸線長度總體上呈減小趨勢, 從 1959年的264703.21m減小到1984年的219959.11m, 減少幅度為16.9%。其中, 從1968年到1974年發(fā)生了一次短暫逆轉(zhuǎn), 由1968年的249696.32m增加到1974年的271577.65m, 增加了 8.8%。在 1984年以后, 岸線長度則總體上呈增加趨勢, 從1984年的219959.11m增加到 2002年的302928.23m, 增加幅度高達(dá) 37.7%。其中, 從1992年到2000年又存在短暫逆轉(zhuǎn), 從1992年的268630.06m減少到2000年的262574.1m, 減少了2.3%。

    從1959年到2002年, 以1984年為轉(zhuǎn)折點(diǎn), 黃河三角洲岸線的年平均變化速度變化過程包括總體上分別由小變大的兩個階段(圖2)。在1984年以前, 岸線年平均變化速率數(shù)值較大, 除了從1968年到1974年有短期的銳減外, 其變化趨勢是總體上由小增大,從1959—1968年的119.15m/a增大到1980—1984年的360.58m/a。在1984年以后, 岸線變化速率數(shù)值較小, 且前期速率為負(fù)值、后期為正值, 速率變化的總趨勢是由小變大, 但變化幅度不顯著, 從1984—1992年的–43.67m/a增大到2000—2002年的31.33m/a。

    從1959年到2002年, 黃河三角洲岸線長度變化和進(jìn)退速率變化均具有顯著的階段性, 即以 1984年為界存在突變, 岸線長度由總體變短轉(zhuǎn)為總體增長,岸線進(jìn)退速率變化則以此為界由前、后兩個顯著不同的階段組成。對比發(fā)現(xiàn), 1984年以前的岸線長度總體減小對應(yīng)于同期的岸線變化速率總體增大, 1984年以后的岸線長度總體增加對應(yīng)于同期的岸線變化速率總體減小。其中, 1968—1974年的岸線長度短暫增加對應(yīng)于同期的岸線速率減小, 而 1992—2000年的岸線長度短暫減小對應(yīng)于同期的岸線速率增加。

    3.2 岸線分維數(shù)變化

    從1959年到2002年, 以1984年為轉(zhuǎn)折點(diǎn), 黃河三角洲岸線分維數(shù)的變化趨勢相反(圖3)。在1984年以前, 三角洲岸線分維數(shù)變化總體呈減小趨勢, 從1959年的 1.5780減小到 1984年的 1.2997, 減少了17.6%。其中, 從1968年到1974年岸線分維數(shù)變化趨勢存在短暫的逆轉(zhuǎn), 從 1968年的 1.4157增加到1974年的1.5299, 增加了8.1%。在1984年以后, 三角洲岸線分維數(shù)變化則呈總體增加趨勢, 從 1984年的1.2997增加到2002年的1.4163, 增加了9.0%。其中, 從1992年到2000年岸線分維數(shù)變化趨勢存在短暫的逆轉(zhuǎn), 從 1992年的 1.3951減少到 2000年的1.3335m, 減少了4.4%。

    如前所述, 從1959年到2002年, 黃河三角洲岸線的長度和進(jìn)退速率變化具有顯著的階段性, 1984年以前岸線的長度減小、進(jìn)退速率增加, 1984年以后岸線長度增加、進(jìn)退速率減小。因此, 比較而言, 1959年到1984年的三角洲岸線進(jìn)退速率較大、岸線長度減小, 同期岸線的平面展布形態(tài)變得平滑、簡單;1984年到2002年的三角洲岸線進(jìn)退速率較小、岸線長度增加, 同期岸線的平面展布形態(tài)變得曲折、復(fù)雜。

    圖3 1959—2002年黃河三角洲岸線分維數(shù)變化Fig.3 Coastline fractal dimensions change of the Yellow River delta from 1959 to 2002

    圖4 1959—2002年黃河三角洲前進(jìn)、后退和平衡岸線所占比例Fig.4 Percentage of forward、backward and balanced coastline of the Yellow River delta from 1959 to 2002

    3.3 前進(jìn)、后退和平衡岸線比例構(gòu)成變化

    從1959年到2002年, 除1974—1980年(43.16%)、1984—1992 年(35.29%)、1992—1996 年(17.76%)外,黃河三角洲前進(jìn)岸線占全部岸線比例均在 50%以上(圖4)。其中, 1980—1984年、1996—2000年均超過70%, 分別為 76.53%、73.58%(圖 4)。其變化過程由三個自高而低的階段構(gòu)成, 三個階段之間存在兩次由低而高的突變(圖4)。其中, 第一個階段從1959年到1980年, 前進(jìn)岸線比例由56.48%(1959—1968年)降至 43.16%(1974—1980年), 后退及平衡岸線比例之和則相應(yīng)由 43.52%升至 56.85%; 第二個階段從1980年到 1996年, 前進(jìn)岸線比例由 75.53%(1980—1984年)降至17.76%(1992—1996年), 后退及平衡岸線比例之和則相應(yīng)由 24.47%升至 82.24%; 第三個階段由 1996年到 2002年, 前進(jìn)岸線比例由 73.58%(1996—2000年)降至56.48%(2000—2002年), 后退及平衡岸線比例之和則相應(yīng)由26.42%升至43.52%。

    在后退岸線和平衡岸線中, 后退岸線比例超過50%以上者, 僅出現(xiàn)在 1974—1980年(50.53%)和1984—1992年(55.88%), 其余時段均在 31%以下(圖4)。其中, 1974—1980年后退岸線比例超過了前進(jìn)岸線。平衡岸線所占比例超過 50%以上者, 僅出現(xiàn)在1992—1996年(56.07%), 其余時段則均在30%以下。其中, 1968—1974年平衡岸線(27.37%)超過了后退岸線(18.95%), 而 1992—1996年平衡岸線則超過了前進(jìn)岸線(17.76%)和后退岸線(26.17%)之和。

    3.4 岸線變化速率的時空分布

    根據(jù)沿各切線計算的終點(diǎn)變化速率, 在 1959年到2002年的40多年中, 黃河三角洲岸線變化具有顯著的階段性和區(qū)域性(圖 5)。從時間上看, 岸線速率變化過程包括三個不同的階段。在從1959年到1974年的15年中, 1959—1968年清水溝以西岸段變化速率較大、以東變化不明顯, 而到 1968—1974年幾乎全部岸段均轉(zhuǎn)為低速率階段。在從1974年到1996年的22年中, 1974—1984年總體上岸線變化速率最大,到1992—1996年轉(zhuǎn)為低速率階段, 1984—1992年岸線變化速率則介于其間。在從1996年到2002年的6年中, 岸線變化速率數(shù)值總體上較大, 而且 1996—2000年和2000—2002年兩個階段的變化速率數(shù)值差別不大, 但是同一岸段的變化方向卻大致相反。

    從1959年到1968年, 黃河三角洲清水溝以東岸線進(jìn)退速率均很小; 以西則岸線進(jìn)退速率很大, 且前進(jìn)段與后退段在空間上交替出現(xiàn)。其中, 前進(jìn)段位于釣河口、老神仙溝口附近, 在兩河口之間及釣河口以西為后退段。從1968年到1974年, 黃河三角洲岸線總體穩(wěn)定, 岸線進(jìn)退速率均很小, 僅神仙溝口與釣河口之間局部岸段進(jìn)退速率較大。從1974年到1980年,除釣河口以西、清水溝口附近岸線以較大速率前進(jìn)外,黃河三角洲大部分岸段岸線發(fā)生低速后退, 且速率沿平行海岸方向波動。

    從1980年到1984年, 黃河三角洲岸線總體上高速前進(jìn), 但不同岸段速度差別很大, 以神仙溝口附近最大。從1984年到1992年, 除行水河口附近岸段以較大速率前進(jìn)外, 大部分岸線發(fā)生低速后退和前進(jìn)。從1992年到1996年, 三角洲岸線整體穩(wěn)定, 僅河口附近有小范圍的低速前進(jìn)。從1996年到2002年, 位于中段的清水溝口到神仙溝口之間岸線穩(wěn)定, 但三角洲兩翼部位的岸線進(jìn)退卻較顯著, 而且前進(jìn)段與后退段在空間上發(fā)生頻繁交替。其中, 同一岸段1996—2000年和2000—2002年的岸線遷移方向大致相反。

    從岸線進(jìn)退的空間分布看, 在從1959年到2002年的 40多年中, 黃河三角洲岸線不同階段的前進(jìn)岸段均位于該時期的河口附近, 而岸線后退和基本岸段均位于前進(jìn)岸段之間及三角洲北、東兩翼(圖 5)。其中, 前期河口均位于黃河三角洲北部渤海灣沿岸,此時岸線進(jìn)退主要發(fā)生在河口附近, 東部的萊州灣沿岸則進(jìn)退不顯著。后期河口位于三角洲東部萊州灣沿岸, 此時除河口附近岸線發(fā)生劇烈進(jìn)退外, 北部的渤海灣沿岸進(jìn)退變化也非常劇烈。

    4 討論

    塑造黃河三角洲海岸的物質(zhì)基礎(chǔ)是黃河入海泥沙, 而黃河入海泥沙最近 60年來總體上呈顯著的減少趨勢。根據(jù)利津水文站觀測數(shù)據(jù), 1954—1964年黃河入海泥沙量11.7億t/a, 1964—1976年11.2億t/a, 而1976—1996年只有5.3億t/a, 1996年以后則不到2.6億t/a。進(jìn)一步的分析表明, 黃河入海泥沙的減少過程具有顯著的突變性, 年入海泥沙量于 1968、1985、1996年突變, 汛期泥沙量于1959、1985、1996年突變, 非汛期泥沙量則突變于1975年(樊輝等, 2009)。其中, 從1959年到1985年, 黃河入海泥沙數(shù)總量較多(達(dá) 5億 t以上), 故三角洲平原以較大的速率向海淤進(jìn), 海岸地貌在縱向上向夷平方向演替, 海岸線變得平滑、簡單, 整個三角洲岸線長度減小; 從 1985年到 2002年, 入海泥沙數(shù)量大幅度減少, 三角洲海岸的淤積作用減弱、侵蝕作用增強(qiáng), 海岸線變得曲折、復(fù)雜, 整個三角洲岸線長度增大(圖 2、3、4)。據(jù)此推測, 黃河入海泥沙總量5億t可能是三角洲岸線變化趨勢的轉(zhuǎn)折點(diǎn), 即利津水文站年泥沙總量在 5億t以上時, 三角洲岸線變短、曲折性降低, 而5億t以下則岸線變長、曲折性增加。

    受黃河入海泥沙數(shù)量減少的影響, 從 1959年到1980年雖然黃河三角洲海岸總體上為淤積前進(jìn)、岸線總長度減小, 但前進(jìn)岸線的比例卻發(fā)生了顯著減少(圖 4)。導(dǎo)致這種現(xiàn)象的因素, 除入海泥沙數(shù)量的影響以外, 還受到河口河道地貌演變的影響。其中,從1959年到1960年為神仙溝流路的單股河道末期,從1960年起在四號樁剖面開始出汊, 從1968年起釣口流路單股成河, 從1972年起在羅13斷面上游附近出汊, 從1976年起又在清水溝出汊, 直至1980年清水溝流路單股河道形成。因此, 在1980年前的21年中絕大部分時間不存在單股河道, 有大量河流輸沙在入海前即堆積在沿河低洼部位。但是, 1980年清水溝流路河口單股成河后, 因河道突然束狹、河床縱比降大, 河流輸沙中實(shí)際入海的比例顯著增加。因此, 雖然位于河口段上游的利津站泥沙量減少, 從1980年到1984年前進(jìn)岸線比例卻呈短暫的突然增加(圖4)。此后, 隨著清水溝單股成河后河道向海延長, 河床縱比降變小, 河道輸沙能力減小, 有大量泥沙淤積在河灘上, 其結(jié)果是前進(jìn)岸線的比例再次減小(圖 4)。

    圖5 1959—2002年黃河三角洲不同時期的岸線終點(diǎn)速率及其分布Fig.5 The End Point Rate and distribution of coastline during different periods of the Yellow River delta from 1959 to 2002

    從海岸線進(jìn)退速率變化及其空間分布看, 在從1959年到 2002年的 40多年中, 1968—1974年、1992—1996年兩個階段中的黃河三角洲岸線總體上進(jìn)退變化不顯著(圖 5)。這說明在黃河入海泥沙總體減少背景下, 不管河口處于三角洲北部還是東部, 黃河三角洲岸線在40年尺度上曾兩度達(dá)到或接近平衡狀態(tài), 其原因應(yīng)與黃河三角洲特有的海岸沉積動力體系有關(guān)。如前所述, 黃河口潮流為平行河道的往復(fù)潮流, 而河口兩側(cè)海岸為大致平行岸線、垂直河口河道的往復(fù)潮流。根據(jù) 1984年觀測, 在三角洲東北部沿岸常年存在弧形高速流帶, 在三角洲東部沿岸則存在平行海岸的季節(jié)性高流速帶, 其中心區(qū)域的潮流流速均在120cm/s以上(薛鴻超等, 1996)。另據(jù)研究表明, 三角洲東北部岸外海域波高 1.5m波浪的擾動水深6.6—10.8m, 與上述潮流高流速帶所處位置的水深相當(dāng), 即波浪擾動區(qū)和強(qiáng)潮流區(qū)在空間上一致(陳沈良等, 2004)。河口往復(fù)流帶、近岸高流速帶、波浪破碎帶的這種空間配置關(guān)系, 有利于河流入海泥沙沿海岸的縱向輸運(yùn)、堆積, 進(jìn)而起到海岸夷平的地貌效應(yīng)。簡單而言, 在新河口形成初期, 河口均較兩側(cè)海岸顯著凹進(jìn)、距口外高流速帶較遠(yuǎn), 大量泥沙在河口附近堆積, 河口附近岸線淤進(jìn)迅速, 而原廢棄河口附近則蝕退迅速。此后, 隨著河口沙嘴形成并向海突出、延伸, 河口往復(fù)流帶與高流速帶距離減小甚至重疊, 大量泥沙輸送到河口以外并沿岸輸運(yùn)、堆積, 因而海岸淤積、侵蝕均變緩, 淤進(jìn)岸線和蝕退岸線比例均減小, 而進(jìn)退平衡岸線比例增大, 其中又尤其以1992—1996年最為顯著(圖5)。

    5 結(jié)論

    綜上所述, 基于岸線數(shù)字分析和分形分析結(jié)果,從1959年到1984年黃河三角洲岸線進(jìn)退速率增大、岸線長度減小, 岸線的平面輪廓展布形態(tài)變得平滑、簡單; 從1984年到2002年三角洲岸線進(jìn)退速率變小、岸線長度增加, 岸線的平面輪廓展布形態(tài)變得曲折、復(fù)雜。另一方面, 前進(jìn)、后退及沖淤進(jìn)退平衡岸線的比例構(gòu)成變化具有顯著的階段性??傮w上, 前進(jìn)岸線所占比例呈不斷減小趨勢, 但在 1980—1984年、1996—2000年又有短暫的突然增大。從岸線變化速率及其空間分布看, 在40年中黃河三角洲岸線變化曾經(jīng)于 1968—1974年、1992—1996年兩度達(dá)到或接近平衡狀態(tài)。根據(jù)初步分析, 1959—2002年黃河三角洲岸線演變是在黃河入海泥沙總體減少背景下發(fā)生的, 是黃河口地貌變遷、海岸沉積動力系統(tǒng)共同作用的結(jié)果。

    豐愛平, 夏東興, 谷東起等, 2006. 萊州灣南岸海岸侵蝕過程與原因研究. 海洋科學(xué)進(jìn)展, 24(1): 83—90

    王 穎, 季小梅, 2011. 中國海陸過渡帶——海岸海洋環(huán)境特征與變化研究. 地理科學(xué), 31(2): 129—135

    莊振業(yè), 林振宏, 劉志杰等, 2003. 海平面變化及其海岸響應(yīng).海洋地質(zhì)動態(tài), 19(7): 1—12

    李安龍, 李廣雪, 曹立華等, 2004. 黃河三角洲廢棄葉瓣海岸侵蝕與岸線演化. 地理學(xué)報, 59(5): 731—737

    李蒙蒙, 王 慶, 張安定等, 2013. 最近 50年來萊州灣西——南部淤泥質(zhì)海岸地貌演變研究. 海洋通報, 32(2): 141—151

    張曉祥, 王偉瑋, 嚴(yán)長清等, 2014. 南宋以來江蘇海岸帶歷史海岸線時空演變研究. 地理科學(xué), 34(3): 344—351

    陳沈良, 張國安, 谷國傳, 2004. 黃河三角洲海岸強(qiáng)侵蝕機(jī)理及治理對策. 水利學(xué)報, (7): 1—6, 13

    戰(zhàn) 超, 王 慶, 夏艷玲等, 2013. 膠東半島南部典型海灣地貌過程對灘涂養(yǎng)殖的響應(yīng). 海洋與湖沼, 44(2): 283—291

    栗云召, 于君寶, 韓廣軒等, 2012. 基于遙感的黃河三角洲海岸線變化研究. 海洋科學(xué), 36(4): 99—106

    諶艷珍, 方國智, 倪 金等, 2010. 遼河口海岸線近百年來的變遷. 海洋學(xué)研究, 28(2): 14—21

    樊 輝, 劉艷霞, 黃海軍, 2009. 1950—2007年黃河入海水沙通量變化趨勢及突變特征. 泥沙研究, (5): 9—16

    薛春汀, 2009. 7000年來渤海西岸、南岸海岸線變遷. 地理科學(xué), 29(2): 217—222

    薛鴻超, 謝金贊, 1996. 中國海岸帶水文. 北京: 海洋出版社,1—265

    Chu Z X, Sun X G, Zhai S K et al, 2006. Changing pattern of accretion/erosion of the modern Yellow River (Huanghe)subaerial delta, China: Based on remote sensing images.Marine Geology, 227(1—2): 13—30

    Cui B L, Li X Y, 2011. Coastline change of the Yellow River estuary and its response to the sediment and runoff (1976—2005). Geomorphology, 127(1—2): 32—40

    Fenster M S, Dolan R, Morton R A, 2001. Coastal storms and shoreline change: signal or noise? Journal of Coastal Research, 17(3): 714—720

    Hu C H, Ji Z W, Wang T, 1998. Dynamic characteristics of sea currents and sediment dispersion in the Yellow River Estuary.International Journal of Sediment Research, 13(2): 16—26

    Jabaloy-Sánchez A, Lobo F J, Azor A et al, 2010. Human-driven coastline changes in the Adra River deltaic system, southeast Spain. Geomorphology, 119(1—2): 9—22

    Jabaloy-Sánchez A, Lobo F J, Azor A et al, 2014. Six thousand years of coastline evolution in the Guadalfeo deltaic system (southern Iberian Peninsula). Geomorphology, 206: 374—391

    Kaliraj S, Chandrasekar N, Magesh N S, 2014. Impacts of wave energy and littoral currents on shoreline erosion/accretion along the south-west coast of Kanyakumari, Tamil Nadu using DSAS and geospatial technology. Environmental Earth Sciences, 71(10): 4523—4542

    Kuleli T, Guneroglu A, Karsli F et al, 2011. Automatic detection of shoreline change on coastal Ramsar wetlands of Turkey.Ocean Engineering, 38(10): 1141—1149

    Li G X, Tang Z S, Yue S H et al, 2001. Sedimentation in the shear front off the Yellow River mouth. Continental Shelf Research, 21(6—7): 607—625

    Li G X, Zhuang K L, Wei H L, 2000. Sedimentation in the Yellow River delta. Part III. Seabed erosion and diapirism in the abandoned subaqueous delta lobe. Marine Geology,168(1—4): 129—144

    Qiao S Q, Shi X F, Zhu A M et al, 2010. Distribution and transport of suspended sediments off the Yellow River(Huanghe) mouth and the nearby Bohai Sea. Estuarine,Coastal and Shelf Science, 86(3): 337—344

    Wang H J, Yang Z S, Li G X et al, 2006. Wave climate modeling on the abandoned Huanghe (Yellow River) delta lobe and related deltaic erosion. Journal of Coastal Research, 22(4):906—918

    天美传媒精品一区二区| www.av在线官网国产| 真实男女啪啪啪动态图| 亚洲三级黄色毛片| 国产v大片淫在线免费观看| 高清毛片免费看| 国产高清三级在线| 日产精品乱码卡一卡2卡三| 日韩一本色道免费dvd| 在线免费观看不下载黄p国产| 美女xxoo啪啪120秒动态图| 精品久久久久久久久av| 777米奇影视久久| 看非洲黑人一级黄片| 天天躁夜夜躁狠狠久久av| 中文字幕av成人在线电影| 麻豆久久精品国产亚洲av| 免费播放大片免费观看视频在线观看| 国产高清不卡午夜福利| 国产男女超爽视频在线观看| 成人毛片60女人毛片免费| 日韩欧美 国产精品| 日韩欧美精品v在线| 日韩一本色道免费dvd| 少妇人妻一区二区三区视频| 免费看不卡的av| 亚洲精品成人av观看孕妇| videossex国产| 18禁裸乳无遮挡免费网站照片| 在线天堂最新版资源| 人人妻人人看人人澡| 天天一区二区日本电影三级| 国产精品一区二区性色av| 亚洲欧美一区二区三区黑人 | 精品久久久久久久人妻蜜臀av| 热re99久久精品国产66热6| 免费看不卡的av| a级一级毛片免费在线观看| 日韩一区二区三区影片| 丰满少妇做爰视频| 我的老师免费观看完整版| 国产白丝娇喘喷水9色精品| 久热这里只有精品99| 三级经典国产精品| 99视频精品全部免费 在线| 我的老师免费观看完整版| 日韩伦理黄色片| 中国三级夫妇交换| 成人特级av手机在线观看| 我要看日韩黄色一级片| 欧美一区二区亚洲| 97人妻精品一区二区三区麻豆| 白带黄色成豆腐渣| 在线看a的网站| 欧美成人午夜免费资源| 99久国产av精品国产电影| 亚洲精品日本国产第一区| 亚洲成人一二三区av| 国产成人精品一,二区| 麻豆成人av视频| 男人添女人高潮全过程视频| 久热这里只有精品99| 在线播放无遮挡| 成人午夜精彩视频在线观看| 国产午夜精品一二区理论片| 在线观看国产h片| 丰满乱子伦码专区| 欧美一区二区亚洲| 国产精品久久久久久精品电影小说 | 亚洲国产精品国产精品| 赤兔流量卡办理| 欧美最新免费一区二区三区| 蜜臀久久99精品久久宅男| 国产精品国产三级国产av玫瑰| 国产欧美另类精品又又久久亚洲欧美| 国产av不卡久久| 欧美潮喷喷水| 亚洲图色成人| 日韩一区二区视频免费看| 99re6热这里在线精品视频| 国产精品一区二区在线观看99| 日本熟妇午夜| 久久鲁丝午夜福利片| 夜夜看夜夜爽夜夜摸| 男女边摸边吃奶| 下体分泌物呈黄色| 51国产日韩欧美| 在线 av 中文字幕| 晚上一个人看的免费电影| 日韩在线高清观看一区二区三区| av天堂中文字幕网| 亚洲精品亚洲一区二区| 亚洲av福利一区| av国产久精品久网站免费入址| 国产亚洲一区二区精品| 日韩在线高清观看一区二区三区| 肉色欧美久久久久久久蜜桃 | 老司机影院成人| 99热这里只有精品一区| kizo精华| 天堂中文最新版在线下载 | 日本与韩国留学比较| 国内揄拍国产精品人妻在线| av.在线天堂| 久久人人爽av亚洲精品天堂 | 免费不卡的大黄色大毛片视频在线观看| 男女边摸边吃奶| 国产亚洲5aaaaa淫片| 国产女主播在线喷水免费视频网站| 九九在线视频观看精品| 久久99蜜桃精品久久| 亚洲av中文av极速乱| 22中文网久久字幕| 亚洲真实伦在线观看| 五月伊人婷婷丁香| 99热网站在线观看| 欧美日韩在线观看h| 亚洲av.av天堂| 国产真实伦视频高清在线观看| 啦啦啦在线观看免费高清www| 久久精品综合一区二区三区| 国模一区二区三区四区视频| 在线免费观看不下载黄p国产| 国产色爽女视频免费观看| 狠狠精品人妻久久久久久综合| 亚洲欧美一区二区三区国产| 国产成人精品一,二区| www.av在线官网国产| 国产白丝娇喘喷水9色精品| 少妇猛男粗大的猛烈进出视频 | 日韩一区二区三区影片| 天堂俺去俺来也www色官网| 日本三级黄在线观看| 国产在线男女| 人人妻人人澡人人爽人人夜夜| 日韩av不卡免费在线播放| 菩萨蛮人人尽说江南好唐韦庄| 国产精品女同一区二区软件| 下体分泌物呈黄色| av在线app专区| 18+在线观看网站| 欧美人与善性xxx| 中国美白少妇内射xxxbb| 日韩成人伦理影院| 亚洲欧美一区二区三区黑人 | 女的被弄到高潮叫床怎么办| www.av在线官网国产| 亚洲精品乱码久久久久久按摩| 亚洲自拍偷在线| 午夜福利高清视频| 肉色欧美久久久久久久蜜桃 | 免费看日本二区| 另类亚洲欧美激情| 大话2 男鬼变身卡| 少妇的逼好多水| 亚洲精品成人av观看孕妇| 亚洲,欧美,日韩| 国产 一区精品| 视频中文字幕在线观看| 亚洲av一区综合| 男人和女人高潮做爰伦理| 欧美丝袜亚洲另类| 久久午夜福利片| 少妇人妻一区二区三区视频| 久久久国产一区二区| 久久精品熟女亚洲av麻豆精品| 亚洲三级黄色毛片| 欧美区成人在线视频| 91久久精品电影网| 久久精品国产自在天天线| 婷婷色综合大香蕉| 亚洲精品成人久久久久久| 成人国产av品久久久| 欧美成人一区二区免费高清观看| 欧美zozozo另类| 午夜福利视频1000在线观看| 亚洲精品第二区| 亚洲精品456在线播放app| 人人妻人人看人人澡| 我的老师免费观看完整版| 高清日韩中文字幕在线| 亚洲精品影视一区二区三区av| 欧美3d第一页| 亚洲在久久综合| 婷婷色麻豆天堂久久| 日韩欧美 国产精品| 美女被艹到高潮喷水动态| 黄色一级大片看看| 你懂的网址亚洲精品在线观看| 国产成人91sexporn| 欧美成人一区二区免费高清观看| 又爽又黄a免费视频| 国产毛片a区久久久久| 日韩av不卡免费在线播放| 免费电影在线观看免费观看| 亚洲aⅴ乱码一区二区在线播放| 午夜亚洲福利在线播放| 亚洲欧美中文字幕日韩二区| 狠狠精品人妻久久久久久综合| 狠狠精品人妻久久久久久综合| 小蜜桃在线观看免费完整版高清| 网址你懂的国产日韩在线| 午夜福利高清视频| 日韩不卡一区二区三区视频在线| 国产男女内射视频| 亚洲,一卡二卡三卡| 亚洲欧美中文字幕日韩二区| 日日啪夜夜爽| 十八禁网站网址无遮挡 | 久久韩国三级中文字幕| 久久韩国三级中文字幕| 国产成人aa在线观看| 成年免费大片在线观看| 蜜桃亚洲精品一区二区三区| 日本猛色少妇xxxxx猛交久久| 亚洲最大成人av| 嫩草影院新地址| 男人爽女人下面视频在线观看| 成年av动漫网址| 日韩国内少妇激情av| 黄片无遮挡物在线观看| 黄色怎么调成土黄色| 久久99热6这里只有精品| 国产成人一区二区在线| av播播在线观看一区| 欧美国产精品一级二级三级 | 成人美女网站在线观看视频| 伊人久久精品亚洲午夜| 91aial.com中文字幕在线观看| 国产综合精华液| 久久久久久久大尺度免费视频| 人妻 亚洲 视频| 国产精品福利在线免费观看| 午夜爱爱视频在线播放| 精品国产露脸久久av麻豆| 日本猛色少妇xxxxx猛交久久| 在线 av 中文字幕| 热re99久久精品国产66热6| 亚洲天堂国产精品一区在线| 人人妻人人看人人澡| 男人舔奶头视频| 色哟哟·www| 51国产日韩欧美| 91精品一卡2卡3卡4卡| 国产真实伦视频高清在线观看| 日本爱情动作片www.在线观看| 成人美女网站在线观看视频| 国产精品无大码| 国产av国产精品国产| 丝袜美腿在线中文| 97热精品久久久久久| 日韩一区二区三区影片| 精品一区二区三卡| 欧美激情国产日韩精品一区| 午夜福利网站1000一区二区三区| 国产成人91sexporn| 国国产精品蜜臀av免费| 亚洲人与动物交配视频| 色5月婷婷丁香| 亚洲成人中文字幕在线播放| 日日摸夜夜添夜夜爱| 日本av手机在线免费观看| 99久久精品国产国产毛片| 精品人妻一区二区三区麻豆| 高清视频免费观看一区二区| 中文欧美无线码| 搡女人真爽免费视频火全软件| 男女国产视频网站| 黄色视频在线播放观看不卡| 狠狠精品人妻久久久久久综合| 国产在线男女| 亚洲成人av在线免费| 国产一区二区在线观看日韩| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人美女网站在线观看视频| 亚洲欧美成人精品一区二区| 免费av观看视频| 人人妻人人看人人澡| 欧美日韩精品成人综合77777| 国产免费一区二区三区四区乱码| 国产精品嫩草影院av在线观看| 中国美白少妇内射xxxbb| 亚洲精品久久久久久婷婷小说| 少妇猛男粗大的猛烈进出视频 | 五月伊人婷婷丁香| 男人和女人高潮做爰伦理| 中文资源天堂在线| 精品一区二区免费观看| 国产老妇伦熟女老妇高清| 色5月婷婷丁香| 久久精品久久久久久噜噜老黄| 老司机影院成人| 日韩,欧美,国产一区二区三区| 在线观看三级黄色| 亚洲精品456在线播放app| 亚洲av日韩在线播放| 久久ye,这里只有精品| 又爽又黄无遮挡网站| 精品国产一区二区三区久久久樱花 | 欧美老熟妇乱子伦牲交| 久久99精品国语久久久| 成人无遮挡网站| 成人午夜精彩视频在线观看| av免费观看日本| 亚洲最大成人手机在线| 亚洲欧美成人精品一区二区| 草草在线视频免费看| 精品一区二区三区视频在线| 欧美+日韩+精品| 久久久久久久久大av| 国产黄色视频一区二区在线观看| 特级一级黄色大片| 美女国产视频在线观看| 乱系列少妇在线播放| 国产精品久久久久久久久免| 国产亚洲av嫩草精品影院| 免费人成在线观看视频色| 欧美日韩精品成人综合77777| 亚洲欧美清纯卡通| 午夜免费观看性视频| 久久99热这里只有精品18| 亚洲成人久久爱视频| 亚洲国产日韩一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 网址你懂的国产日韩在线| 国产女主播在线喷水免费视频网站| 美女国产视频在线观看| 成人无遮挡网站| 成人漫画全彩无遮挡| 在现免费观看毛片| 亚洲人成网站高清观看| 免费av观看视频| 99久久九九国产精品国产免费| 国产成人免费观看mmmm| 久久国内精品自在自线图片| 国产91av在线免费观看| 最近的中文字幕免费完整| 亚洲av中文字字幕乱码综合| 日韩欧美 国产精品| 免费电影在线观看免费观看| 精品国产露脸久久av麻豆| 亚洲美女视频黄频| 久久久色成人| 成人高潮视频无遮挡免费网站| 精品午夜福利在线看| 中文字幕亚洲精品专区| 国产av码专区亚洲av| 有码 亚洲区| 特大巨黑吊av在线直播| 久久久久久久久久成人| av在线老鸭窝| 夫妻性生交免费视频一级片| 国产精品偷伦视频观看了| 草草在线视频免费看| 亚洲av欧美aⅴ国产| 97在线视频观看| 亚洲成人一二三区av| 日韩国内少妇激情av| 国产精品一及| 中文字幕免费在线视频6| 免费高清在线观看视频在线观看| 日本爱情动作片www.在线观看| 又黄又爽又刺激的免费视频.| av在线观看视频网站免费| 久久午夜福利片| 国产色爽女视频免费观看| 乱系列少妇在线播放| 别揉我奶头 嗯啊视频| 亚洲欧美一区二区三区黑人 | 一级毛片我不卡| 亚洲av不卡在线观看| 欧美日韩亚洲高清精品| 夜夜看夜夜爽夜夜摸| 丰满人妻一区二区三区视频av| 久久久久久伊人网av| 秋霞在线观看毛片| 91精品一卡2卡3卡4卡| 午夜日本视频在线| 久久人人爽av亚洲精品天堂 | 91久久精品电影网| 国产免费视频播放在线视频| 小蜜桃在线观看免费完整版高清| 亚洲综合色惰| 秋霞在线观看毛片| 亚洲欧美精品自产自拍| 国产毛片在线视频| 国产午夜精品一二区理论片| 色5月婷婷丁香| 亚洲成人中文字幕在线播放| av国产精品久久久久影院| 亚洲成人一二三区av| 少妇人妻精品综合一区二区| 欧美高清成人免费视频www| 五月伊人婷婷丁香| 中文乱码字字幕精品一区二区三区| 国产又色又爽无遮挡免| 国产成年人精品一区二区| 新久久久久国产一级毛片| 寂寞人妻少妇视频99o| 嫩草影院新地址| 久久久精品免费免费高清| 亚洲av.av天堂| 亚洲av不卡在线观看| 国产精品偷伦视频观看了| 亚洲av免费高清在线观看| 天天一区二区日本电影三级| 在线观看国产h片| 午夜精品一区二区三区免费看| 黄片无遮挡物在线观看| 亚洲精品日本国产第一区| 美女脱内裤让男人舔精品视频| 欧美zozozo另类| 国产午夜精品一二区理论片| 亚洲精品色激情综合| 亚洲国产成人一精品久久久| 自拍偷自拍亚洲精品老妇| 国产 精品1| 永久网站在线| 97热精品久久久久久| 黄色一级大片看看| 欧美日韩视频高清一区二区三区二| 国产久久久一区二区三区| 神马国产精品三级电影在线观看| 国产在视频线精品| 国产一区二区三区综合在线观看 | 成人免费观看视频高清| 久久久国产一区二区| 国产综合懂色| 精品久久久久久电影网| 国产91av在线免费观看| 欧美日韩精品成人综合77777| 在线免费十八禁| 97在线人人人人妻| 国产黄频视频在线观看| 国产免费视频播放在线视频| 亚洲三级黄色毛片| 男女下面进入的视频免费午夜| 亚洲精品aⅴ在线观看| 2018国产大陆天天弄谢| 久久久久久久大尺度免费视频| 国产极品天堂在线| 好男人视频免费观看在线| 精品视频人人做人人爽| 亚洲自偷自拍三级| 国产精品麻豆人妻色哟哟久久| 亚洲精品成人久久久久久| 一级黄片播放器| 国产永久视频网站| 国产老妇伦熟女老妇高清| 亚洲国产高清在线一区二区三| 亚洲国产最新在线播放| 成人亚洲精品一区在线观看 | 免费不卡的大黄色大毛片视频在线观看| 久久99热6这里只有精品| 国产亚洲5aaaaa淫片| 狂野欧美白嫩少妇大欣赏| 91久久精品电影网| 能在线免费看毛片的网站| 国产黄片美女视频| 在线观看免费高清a一片| 不卡视频在线观看欧美| 免费大片18禁| 久久亚洲国产成人精品v| 国国产精品蜜臀av免费| 欧美一级a爱片免费观看看| 国产精品国产三级国产av玫瑰| 亚洲va在线va天堂va国产| 国产黄片美女视频| av国产精品久久久久影院| 欧美人与善性xxx| 丰满乱子伦码专区| 别揉我奶头 嗯啊视频| 免费人成在线观看视频色| 美女主播在线视频| 久久久久国产精品人妻一区二区| 久久女婷五月综合色啪小说 | 日日摸夜夜添夜夜爱| 国产成人午夜福利电影在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜免费男女啪啪视频观看| 大片电影免费在线观看免费| 乱码一卡2卡4卡精品| 成人国产麻豆网| 亚洲精品456在线播放app| 国产精品成人在线| 狂野欧美白嫩少妇大欣赏| 91久久精品电影网| av女优亚洲男人天堂| 日本-黄色视频高清免费观看| 欧美日韩国产mv在线观看视频 | 国产精品99久久久久久久久| 国产国拍精品亚洲av在线观看| 精品99又大又爽又粗少妇毛片| 国产精品成人在线| 99热这里只有是精品在线观看| 婷婷色麻豆天堂久久| 在线免费观看不下载黄p国产| 亚洲av欧美aⅴ国产| 久久久久网色| 国产午夜福利久久久久久| 99久久精品国产国产毛片| 肉色欧美久久久久久久蜜桃 | .国产精品久久| 69人妻影院| 一级黄片播放器| 99re6热这里在线精品视频| 美女xxoo啪啪120秒动态图| a级一级毛片免费在线观看| 亚洲av福利一区| av黄色大香蕉| 国产黄色免费在线视频| 少妇人妻 视频| 久久精品国产a三级三级三级| 插逼视频在线观看| 91久久精品国产一区二区成人| 欧美+日韩+精品| 激情五月婷婷亚洲| 国产成人a区在线观看| 精品一区二区免费观看| 亚洲精品乱码久久久v下载方式| eeuss影院久久| 老司机影院毛片| 国产午夜精品一二区理论片| 丰满人妻一区二区三区视频av| 日本三级黄在线观看| 国产精品久久久久久精品电影| 国产精品久久久久久精品古装| 亚洲欧美一区二区三区黑人 | 亚洲精品亚洲一区二区| av卡一久久| 国产一区二区三区综合在线观看 | freevideosex欧美| 亚洲精品第二区| 亚洲精品,欧美精品| 中文字幕免费在线视频6| 日韩大片免费观看网站| 日韩一区二区视频免费看| 青青草视频在线视频观看| 国产探花极品一区二区| 一区二区三区四区激情视频| 日本熟妇午夜| 亚洲aⅴ乱码一区二区在线播放| 欧美激情国产日韩精品一区| 欧美一区二区亚洲| 日韩,欧美,国产一区二区三区| 少妇高潮的动态图| 少妇裸体淫交视频免费看高清| 久久久久国产网址| 丰满少妇做爰视频| 永久免费av网站大全| 建设人人有责人人尽责人人享有的 | 毛片一级片免费看久久久久| 插逼视频在线观看| 国产片特级美女逼逼视频| 国产精品偷伦视频观看了| 九九爱精品视频在线观看| 亚洲精品色激情综合| 久久久欧美国产精品| 国产高清有码在线观看视频| h日本视频在线播放| 最后的刺客免费高清国语| av.在线天堂| av女优亚洲男人天堂| 国产黄色免费在线视频| 国产精品蜜桃在线观看| 国产精品一及| 亚洲丝袜综合中文字幕| 中国国产av一级| 三级经典国产精品| 男插女下体视频免费在线播放| 国产爱豆传媒在线观看| 热re99久久精品国产66热6| 美女国产视频在线观看| 久热这里只有精品99| 校园人妻丝袜中文字幕| 九九在线视频观看精品| 欧美xxxx性猛交bbbb| 亚洲天堂av无毛| 伊人久久国产一区二区| 精品久久国产蜜桃| 亚洲国产日韩一区二区| 国产亚洲精品久久久com| 亚洲精品日韩在线中文字幕| 极品教师在线视频| 草草在线视频免费看| 国产伦精品一区二区三区四那| 嫩草影院新地址| 高清午夜精品一区二区三区| 国产老妇伦熟女老妇高清| 欧美亚洲 丝袜 人妻 在线| 免费看日本二区| 美女国产视频在线观看| 成人综合一区亚洲| 五月玫瑰六月丁香| 亚洲丝袜综合中文字幕| 久久综合国产亚洲精品| 精品久久久久久久末码| 国产男女超爽视频在线观看| 国模一区二区三区四区视频| 日本免费在线观看一区| 亚洲一区中文字幕在线| 亚洲激情五月婷婷啪啪| 国产无遮挡羞羞视频在线观看| 婷婷色综合大香蕉| 亚洲美女搞黄在线观看| 精品一区在线观看国产| 国产熟女欧美一区二区| 欧美 日韩 精品 国产| 男女高潮啪啪啪动态图| 色网站视频免费| 黄色 视频免费看| www.精华液| 香蕉丝袜av| 18禁观看日本| 成人18禁高潮啪啪吃奶动态图|