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

    黃河三角洲淺層地下水位年內(nèi)變化特征及影響因素

    2021-11-17 09:09:54常茂祥史經(jīng)昊葉思源王祥東李廣雪
    海洋科學(xué) 2021年10期
    關(guān)鍵詞:黃河三角洲蒸發(fā)量井水

    常茂祥, 史經(jīng)昊, 葉思源, 王祥東, 李廣雪

    黃河三角洲淺層地下水位年內(nèi)變化特征及影響因素

    常茂祥1, 2, 史經(jīng)昊3, 葉思源4, 王祥東3, 李廣雪1, 2

    (1.海底科學(xué)與探測(cè)技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室(中國(guó)海洋大學(xué)), 山東 青島 266100; 2.中國(guó)海洋大學(xué)海洋地球科學(xué)學(xué)院, 山東 青島 266100; 3.青島藍(lán)色地球大數(shù)據(jù)科技有限公司, 山東 青島 266400; 4.自然資源部青島海洋地質(zhì)研究所, 山東 青島 266071)

    基于黃河三角洲地區(qū)3個(gè)站位的1小時(shí)一次的淺層地下水位連續(xù)觀測(cè)記錄, 將淺層地下水位與同時(shí)期ERA5高分辨率再分析實(shí)際蒸發(fā)量數(shù)據(jù)結(jié)合, 并聯(lián)合潮位、降水量和徑流量等數(shù)據(jù), 通過(guò)快速傅里葉變換、時(shí)序分析等方法, 分析并闡明了黃河三角洲地區(qū)淺層地下水位變化的特征及機(jī)制。結(jié)果表明: (1) 不同的沉積環(huán)境導(dǎo)致了淺層地下水位整體變化的差異。表層沉積物滲透性相對(duì)較強(qiáng)區(qū)域, 淺層地下水位波動(dòng)劇烈; 而表層沉積物透水性較差地區(qū), 淺層地下水位在6月中上旬存在低谷, 但總體相對(duì)穩(wěn)定; (2) 潮汐對(duì)黃河三角洲淺層地下水位在水平方向的影響范圍至少可達(dá)7 km, 但不超過(guò)15 km。在其影響區(qū)域內(nèi), 淺層地下水位波動(dòng)滯后于潮汐的時(shí)間存在年內(nèi)變化, 分為2個(gè)時(shí)間上持續(xù)各6個(gè)月的區(qū)間, 二者數(shù)值相差約12 h; (3) 降水量與實(shí)際蒸發(fā)量是黃河三角洲淺層地下水位升降最主要的影響因素。此外, 農(nóng)業(yè)活動(dòng)也對(duì)淺層地下水位的變化有一定影響。對(duì)黃河三角洲淺層地下水位變化規(guī)律研究, 能夠?yàn)楸镜貐^(qū)土地鹽堿化、海水入侵災(zāi)害防治與生態(tài)保護(hù)提供科學(xué)依據(jù)。

    黃河三角洲; 淺層地下水位; 潮汐; 滯后性; 降水量; 實(shí)際蒸發(fā)量

    黃河三角洲擁有豐富的石油資源與生態(tài)資源[1-2], 對(duì)當(dāng)?shù)厣鐣?huì)經(jīng)濟(jì)進(jìn)步起到了重要推動(dòng)作用。然而黃河三角洲也是我國(guó)水資源最為短缺的地區(qū)之一, 該地區(qū)90%以上的區(qū)域受到不同程度的土壤鹽堿化影響, 淺層地下水僅小部分為全淡水[3-4], 長(zhǎng)期以來(lái)該區(qū)域的發(fā)展一直受到水資源短缺的阻礙與制約[4]。對(duì)淺層地下水變化規(guī)律及其成因的充分認(rèn)識(shí)在黃河三角洲區(qū)域尤為重要。黃河三角洲區(qū)域僅有少量的淺層地下水長(zhǎng)期觀測(cè)站, 對(duì)于該地區(qū)淺層地下水動(dòng)態(tài)及成因的研究還較為欠缺[5]。亟需結(jié)合可靠的數(shù)據(jù)對(duì)當(dāng)?shù)販\層地下水動(dòng)態(tài)和影響淺層地下水位的因素進(jìn)行全面闡述分析, 為準(zhǔn)確評(píng)估和利用地下水資源服務(wù)。

    影響大河三角洲地下水位變化的機(jī)制繁多, 主要包括自然因素和人為因素。自然因素主要包括: 降水、蒸發(fā)、附近河流徑流量、地質(zhì)條件等, 濱海濕地地區(qū)還受海潮作用; 人為因素則主要為地下水開(kāi)采、工程建設(shè)、農(nóng)業(yè)生產(chǎn)活動(dòng)等。國(guó)內(nèi)外諸多學(xué)者對(duì)大河三角洲地下水動(dòng)態(tài)及其影響因素進(jìn)行了廣泛研究。由于黃河三角洲特殊的濱海地理位置, 潮汐對(duì)該地淺層地下水位具有明顯影響。高茂生等[6]發(fā)現(xiàn)位于黃河三角洲濱海濕地的DY122井地下水位受潮汐波動(dòng)影響明顯, 并滯后潮汐波動(dòng)3 h。在萊州灣南岸, 6個(gè)水文觀測(cè)井淺層地下水位頻率振幅與潮汐規(guī)律表現(xiàn)一致, 每個(gè)觀測(cè)井淺層地下水位波動(dòng)滯后潮汐的時(shí)間唯一[7-8]。以上研究中水位滯后潮位波動(dòng)的時(shí)間只取決于地理位置, 并不隨時(shí)間變化。這些成果所用數(shù)據(jù)的時(shí)間跨度在1 d到3個(gè)月之間, 考慮到黃河三角洲復(fù)雜的地質(zhì)和氣象條件, 短時(shí)間尺度下得出的結(jié)論能否適用于更長(zhǎng)時(shí)間, 還有待探討。在珠江三角洲, Zhang等[9]發(fā)現(xiàn)天文大潮潮位與水位相關(guān)性要高于天文小潮, 說(shuō)明發(fā)生天文大潮時(shí)潮汐對(duì)地下水影響更大。除了潮汐, 降水對(duì)淺層地下水位的影響, 也被廣泛研究。袁瑞強(qiáng)等[10]結(jié)合黃河農(nóng)場(chǎng)地區(qū)3個(gè)淺層地下水觀測(cè)孔與降水資料, 研究結(jié)果表明降水是引起該地淺層地下水位上升的主要原因。孟加拉國(guó)濱海濕地3口監(jiān)測(cè)井的結(jié)果顯示, 淺層地下水位對(duì)于強(qiáng)降水有非??斓奶憫?yīng)[11]。安樂(lè)生等[5]研究了黃河三角洲小清河以北多口觀測(cè)井水位與降水的關(guān)系, 將區(qū)域內(nèi)水位動(dòng)態(tài)分為基本穩(wěn)定型和震蕩起伏型, 地下水動(dòng)態(tài)成因類型分別為“滲入?徑流或水文”型與“滲入?蒸發(fā)”型。在前人有關(guān)蒸發(fā)對(duì)淺層地下水動(dòng)態(tài)影響的研究中, 蒸發(fā)量多采用蒸發(fā)皿觀測(cè)數(shù)據(jù)或公式計(jì)算的潛在蒸發(fā)量數(shù)據(jù)。蒸發(fā)皿測(cè)得數(shù)據(jù)是有限水面在充分供水情況下的蒸發(fā)量, 不能代表真實(shí)的地表蒸發(fā)情況; 潛在蒸發(fā)量是陸面實(shí)際蒸發(fā)量的理論上限, 同樣不能體現(xiàn)自然界的實(shí)際蒸發(fā)情況[12]。因此利用蒸發(fā)皿蒸發(fā)量和潛在蒸發(fā)量分析淺層地下水的研究?jī)H能為了解蒸發(fā)對(duì)淺層地下水位的影響提供了參考, 想要得到更加可靠的定量結(jié)果, 就必須利用更能反映真實(shí)情況的蒸發(fā)數(shù)據(jù)。自然過(guò)程對(duì)大河三角洲淺層地下水動(dòng)態(tài)有重要影響, 但人類活動(dòng)的作用也不可忽視。劉勇等[4]發(fā)現(xiàn), 黃河三角洲地區(qū)地下水超采造成的地下水位劇烈下降會(huì)生成地下水降落漏斗并產(chǎn)生咸水入侵。Minderhoud等[13]發(fā)現(xiàn), 湄公河三角洲地下水過(guò)量開(kāi)采致使蓄水層壓實(shí), 地下水位迅速下降導(dǎo)致該區(qū)域在1991—2016年期間平均沉降量達(dá)約18 cm。在美國(guó)的密西西比河三角洲, Killian等[14]發(fā)現(xiàn)農(nóng)業(yè)灌溉致使的地下水位下降與地表水流量減少密切相關(guān)。對(duì)科羅拉多河三角洲地下水位和微重力的研究結(jié)果顯示: 上游排水導(dǎo)致的當(dāng)?shù)氐叵滤簧攥F(xiàn)象持續(xù)時(shí)間較短; 因抽取地下水產(chǎn)生的低地下水位導(dǎo)致大量河道水滲透散失[15]。

    本文分析了黃河三角洲淺層地下水位動(dòng)態(tài)特征及潮汐、降水量、實(shí)際蒸發(fā)量、黃河徑流量以及農(nóng)業(yè)活動(dòng)與淺層地下水位的關(guān)系, 首次將ERA5高分辨率再分析實(shí)際蒸發(fā)量數(shù)據(jù)與淺層地下水位結(jié)合分析, 揭示了淺層地下水位對(duì)于潮汐波動(dòng)的滯后響應(yīng)存在變化這一現(xiàn)象。

    1 研究區(qū)概況

    研究區(qū)位于黃河三角洲濱海濕地, 緊鄰渤海灣與萊州灣(圖1), 行政區(qū)劃上屬于山東省東營(yíng)市。海拔高程低, 約0至1 m, 東西向自然坡降約為1/10 000[16]。黃河輸沙量大, 歷史上河道變遷頻繁, 因此黃河三角洲沉積環(huán)境不穩(wěn)定, 沉積相復(fù)雜。黃河三角洲低潮線以上區(qū)域可分為兩部分, 即上三角洲平原和下三角洲平原。上三角洲平原包括泛濫平原、決口扇和天然堤等亞環(huán)境。泛濫平原沉積物較薄, 分布面積小, 多是黏土質(zhì)粉砂夾黃色粉砂紋層, 平均粒徑7.7~8.4 φ, 分選差。決口扇分布廣泛, 其沉積物和天然堤沉積物均以粗粉砂為主[17-19]。下三角洲平原中潮灘發(fā)育, 寬度可以超過(guò)10 km。潮灘潮上帶沼澤斑狀分布, 沉積物為黃色粉砂和褐色黏土質(zhì)粉砂互層?,F(xiàn)代河口兩側(cè)的潮間帶沉積物粒徑小, 主要是黃褐色黏土質(zhì)粉砂夾黃色粉砂透鏡體以及黃色粉砂夾黏土質(zhì)粉砂紋層[17]。在水文地質(zhì)分布上, 研究區(qū)處在全咸水區(qū)[4]。地質(zhì)構(gòu)造上位于濟(jì)陽(yáng)盆地的湖相斷陷東營(yíng)坳陷內(nèi), 是典型的古近系半斷陷盆地[18]。溫帶季風(fēng)氣候盛行于研究區(qū), 夏季炎熱多雨, 冬季寒冷干燥, 主要糧食作物為冬小麥。年平均降水量約600 mm, 實(shí)際蒸發(fā)量大于降水量, 年平均氣溫適宜, 約12 ℃[20]。

    2 數(shù)據(jù)來(lái)源與研究方法

    2.1 數(shù)據(jù)來(lái)源

    122-3井、123-2井和124井淺層地下水位使用荷蘭生產(chǎn)的Diver自動(dòng)監(jiān)測(cè)儀進(jìn)行監(jiān)測(cè), 年際水位變化大于1 m, 遠(yuǎn)大于大氣壓強(qiáng)值的變化。在進(jìn)行最終水位值的計(jì)算和分析過(guò)程中, 采用了當(dāng)?shù)仄骄髿鈮簭?qiáng)值1.04×105Pa來(lái)補(bǔ)償校正, 其波動(dòng)趨勢(shì)和范圍能夠近似代表實(shí)際水位動(dòng)態(tài)特征。124井記錄時(shí)間為2009年5月至12月, 其余2井記錄時(shí)間為2009年全年, 水位記錄間隔均為1 h。

    2009年潮位數(shù)據(jù)來(lái)自于孤東驗(yàn)潮站, 數(shù)據(jù)記錄頻率為1次/h。孤東驗(yàn)潮站采用XZY3-1型自動(dòng)水文氣象站連續(xù)監(jiān)測(cè)潮位, 所得數(shù)據(jù)由孤東海堤管理站定期收集校驗(yàn), 在海洋領(lǐng)域得到了較多應(yīng)用[21-24]。

    降水量和實(shí)際蒸發(fā)量數(shù)據(jù)采用2009年ERA5逐小時(shí)數(shù)據(jù), 數(shù)據(jù)范圍為118.5°E~119°E, 37.7°N~37.9°N。文中時(shí)間均采用北京時(shí)間。ERA5是歐洲中期天氣預(yù)報(bào)中心最新一代全球氣象再分析數(shù)據(jù)。相比前一代, ERA5具有更小的延遲性和更高的時(shí)空分辨率(1 h, 0.25°×0.25°), 并能夠更加精確的反映裸地蒸發(fā)量[25]。自發(fā)布以來(lái)ERA5數(shù)據(jù)得到了廣泛應(yīng)用, 并且許多研究證明ERA5數(shù)據(jù)實(shí)際應(yīng)用可靠。利用ERA5降水和溫度等數(shù)據(jù), Xue等[26]研究了一帶一路區(qū)域內(nèi)的干旱時(shí)空分布。黃建平等[27]利用ERA5實(shí)際蒸發(fā)量數(shù)據(jù), 分析了黃河流域近40 a的氣候變化。ERA5數(shù)據(jù)與其他現(xiàn)代再分析和實(shí)測(cè)的降水?dāng)?shù)據(jù)[28-29]以及海洋表層及低空風(fēng)數(shù)據(jù)[30-32]比對(duì)的結(jié)果良好。ERA5與中國(guó)746個(gè)氣象站的相對(duì)濕度數(shù)據(jù)匹配良好[33]。在多種再分析氣象數(shù)據(jù)中, ERA5能夠最好的反映東亞的云特性[34]。在第三極的西風(fēng)和季風(fēng)主導(dǎo)流域源區(qū), ERA5與256個(gè)地面觀測(cè)站降水?dāng)?shù)據(jù)隨海拔變化總體特征一致[35]。在山東及周邊地區(qū)ERA5總體適用性優(yōu)于ERA-Interim再分析數(shù)據(jù)[36]。

    圖1 研究區(qū)衛(wèi)星光學(xué)影像與淺層地下水位觀測(cè)井位分布(Landsat衛(wèi)星, 2018年12月31日)

    黃河徑流量與輸沙量數(shù)據(jù)來(lái)源于水利部黃河水利委員會(huì)官網(wǎng), 站點(diǎn)為東營(yíng)市利津縣的利津站。利津站是黃河干流最下游也是離研究區(qū)最近的重要水文控制站, 能有效反映2009年內(nèi)研究區(qū)黃河的徑流量和輸沙量的月變化。

    2.2 研究方法

    利用Python語(yǔ)言編程過(guò)程中為減少計(jì)算量, 本文采用DFT的快速算法, 即numpy庫(kù)的快速傅里葉變換(FFT)函數(shù): fft.rfft。

    時(shí)間滯后互相關(guān)(time lag cross-correlation)用來(lái)揭示兩個(gè)序列之間超前滯后關(guān)系。潮汐引起的地下水位波動(dòng)信號(hào)受多種長(zhǎng)周期因素干擾, 對(duì)水位數(shù)據(jù)進(jìn)行濾波將有利于更加準(zhǔn)確分析二者關(guān)系。為了探究水位波動(dòng)的滯后性, 本文對(duì)122-3井、123-2井的2009年水位數(shù)據(jù)以及潮位數(shù)據(jù)進(jìn)行滑動(dòng)平均, 再用原始數(shù)據(jù)減去滑動(dòng)平均, 得到去趨勢(shì)以后的數(shù)據(jù)用于時(shí)間滯后互相關(guān)分析。本文采取對(duì)處理后的水位數(shù)據(jù)逐時(shí)向前平移, 然后計(jì)算平移后水位和潮位皮爾森相關(guān)系數(shù)的時(shí)間滯后互相關(guān)算法。若水位數(shù)據(jù)向前平移小時(shí)后, 水位與潮位數(shù)據(jù)具有強(qiáng)相關(guān)性, 則可認(rèn)為水位滯后潮位h。

    在研究降水蒸發(fā)對(duì)水位的影響中, 主要運(yùn)用相關(guān)分析方法, 分析水位累計(jì)升幅降幅與累積降水量蒸發(fā)量的相關(guān)性。相關(guān)系數(shù)的計(jì)算仍采用皮爾森算法。皮爾森相關(guān)系數(shù)自20世紀(jì)由英國(guó)學(xué)者皮爾森提出以來(lái), 在衡量?jī)蓚€(gè)變量的線性相關(guān)程度中得到了廣泛應(yīng)用, 其公式為:

    3 結(jié)果與討論

    3.1 結(jié)果

    3.1.1 年內(nèi)淺層地下水位動(dòng)態(tài)變化

    2009年5—12月份, 124井水位明顯高于其他兩井, 平均水位分別高出同時(shí)期122-3井和123-2井3.25 m和2.31 m(表1)。雨季中124井水位波動(dòng)頻繁劇烈, 雨季過(guò)后程度減弱。1—2月份122-3井與123-2井水位平穩(wěn)。3—4月份122-3井水位緩慢上升, 123-2井則呈下降趨勢(shì)。5—6月份兩井水位波動(dòng)明顯, 在6月中上旬兩井均存在顯著的水位低谷。7月中旬后至9月份水位維持在較高水平, 后較平穩(wěn)下降(圖2)。整體而言, 124井水位較高, 水位不穩(wěn)定; 122-3井與123-2井水位較低, 除5—6月份外, 水位較為平穩(wěn)。

    表1 2009年內(nèi)研究區(qū)各觀測(cè)井水位動(dòng)態(tài)信息統(tǒng)計(jì)

    圖2 2009年內(nèi)研究區(qū)淺層地下水位與降水蒸發(fā)動(dòng)態(tài)

    3.1.2 潮汐與淺層地下水位

    潮汐效應(yīng)[37]指出, 海洋的潮汐作用會(huì)使濱海含水層中的地下水位隨海平面起落有規(guī)律地波動(dòng)。為研究黃河三角洲地區(qū)淺層地下水位與潮汐的關(guān)系, 利用快速傅里葉變換將孤東驗(yàn)潮站潮位數(shù)據(jù)從時(shí)域轉(zhuǎn)為頻域, 繪制頻率振幅譜??梢钥闯? 孤東潮汐數(shù)據(jù)頻譜在0.042 h–1和0.083 h–1處存在兩個(gè)主要的波峰, 且前者遠(yuǎn)大于后者(圖3a)。

    同樣基于快速傅里葉變換, 對(duì)研究區(qū)三口地下水觀測(cè)井水位進(jìn)行頻譜分析。122-3井與123-2井與潮汐數(shù)據(jù)頻譜特征一致, 在約0.042 h–1和0.083 h–1頻率處均有較明顯波峰, 且0.042 h–1左右振幅遠(yuǎn)大于0.083 h–1。122-3井在兩個(gè)頻率的振幅上遠(yuǎn)超123-2井(圖3a, b, c)。然而124井表現(xiàn)出了與潮汐、122-3井和123-2井水位明顯不同的頻譜, 124井水位頻譜未顯示明顯的振幅高值(圖3d)。

    除了頻域, 時(shí)域上的分析也可以揭示潮汐與地下水位的關(guān)系。由于124井在頻域分析中并未顯示海潮對(duì)其有明顯影響, 所以接下來(lái)只分析其余兩井。海潮波動(dòng)信號(hào)傳遞給地下水時(shí), 由于陸地的阻尼作用會(huì)產(chǎn)生延遲。潮位與122-3井、123-2井的水位時(shí)間滯后互相關(guān)分析結(jié)果如圖4所示。圖4橫坐標(biāo)指示向前平移水位數(shù)據(jù)的時(shí)間, 縱坐標(biāo)代表平移后水位數(shù)據(jù)與潮位的相關(guān)系數(shù)?;ハ嚓P(guān)分析顯示, 122-3井水位在向前平移0 h時(shí), 與潮汐相關(guān)系數(shù)最高, 但僅為0.01; 123-2井水位與潮汐相關(guān)系數(shù)在平移15 h時(shí)達(dá)到最大, 也僅為0.02(圖4)。

    圖3 潮位數(shù)據(jù)與淺層地下水位數(shù)據(jù)頻譜分析圖

    注: (a) 潮汐, (b) 122-3井, (c) 123-2井, (d) 124井

    圖4 淺層地下水位與潮汐時(shí)間滯后互相關(guān)結(jié)果

    122-3井與123-2井水位變化具有潮汐規(guī)律(圖3), 滯后性應(yīng)當(dāng)存在, 但這與潮位和水位時(shí)間滯后互相關(guān)分析得到的極低的相關(guān)系數(shù)(圖4)似乎矛盾。天文大潮日時(shí)潮汐振幅最大, 對(duì)地下水位的影響也最為明顯。所以為進(jìn)一步研究滯后性, 并考慮潮汐數(shù)據(jù)是否殘缺和相關(guān)時(shí)間段前后有無(wú)強(qiáng)降水過(guò)程等因素, 分別在2009年的每個(gè)月選取最接近天文大潮日的兩整日(共24 d)進(jìn)行分析(圖5)。

    圖5的12個(gè)子圖中, 潮汐與水位波動(dòng)均表現(xiàn)了明顯的全日潮特征, 這與頻域分析結(jié)果一致。不同地下水觀測(cè)井水位對(duì)潮汐影響的響應(yīng)存在差異, 122-3井波動(dòng)恒早于123-2井, 提前響應(yīng)時(shí)間穩(wěn)定, 約2 h; 122-3井1日內(nèi)的水位振幅是123-2井的2到3倍。在12個(gè)天文大潮期中, 潮汐波動(dòng)呈現(xiàn)兩種態(tài)勢(shì)。3月到8月潮汐曲線呈英文字母“W”形, 其余半年潮位呈現(xiàn)“M”形狀。在“W”形潮汐高潮的時(shí)刻, “M”形潮汐處于低潮, 反之亦然。以122-3井水位為例, “W”形潮汐下122-3井水位變化滯后潮汐時(shí)間為15~18 h; “M”形潮汐下122-3井水位波動(dòng)滯后潮汐2到6 h。

    3.1.3 降水量、實(shí)際蒸發(fā)量與淺層地下水位

    降水是對(duì)地下水水量的直接補(bǔ)給, 對(duì)地下水影響最為直接。黃河三角洲地區(qū)位于東亞季風(fēng)區(qū)內(nèi), 降水呈現(xiàn)明顯的季節(jié)性分布。2009年, 研究區(qū)年總降水量平均為650 mm, 5到9月份降水量為504 mm, 占全年降水量的78%。相比于122-3井和123-2井, 124井水位對(duì)于降水更加敏感。在降水量較多的7—9月份, 124井水位高于其他月份。124井水位在5月9日、6月19日、7月8日、8月22日、9月5日與10月30日左右均有較大幅度的迅速抬升, 這數(shù)次突增均精確對(duì)應(yīng)著高強(qiáng)度降水過(guò)程(圖2), 上述幾日的日降水量均遠(yuǎn)高于黃河三角洲地區(qū)日平均1.78 mm的降水量, 分別為40.04 mm、23.90 mm、26.51 mm、14.33 mm、23.50 mm與23.14 mm。這表明, 無(wú)論從季節(jié)尺度, 還是日際尺度, 124井水位均在很大程度上受降水調(diào)控。122-3井與123-2井水位也在降水較多的7~9月份表現(xiàn)出總體較高水位。對(duì)于較強(qiáng)降水事件, 兩井也有響應(yīng), 但程度較弱。例如4月15日、5月9日、6月19日和10月30日左右的水位抬升(圖2)。

    為更直觀地研究黃河三角洲地區(qū)降水對(duì)淺層地下水位影響, 對(duì)3個(gè)井位水位數(shù)據(jù)對(duì)應(yīng)研究時(shí)段的累積降水量和累計(jì)水位升幅的關(guān)系進(jìn)行分析(圖6)。3個(gè)井各自的累積降水量與累計(jì)水位升幅的時(shí)間序列曲線呈現(xiàn)極高相關(guān)性, 具體表現(xiàn)在各自井位的2條曲線拐點(diǎn)幾乎同步, 枯水期各自井位2條曲線基本一致。在2009年122-3井與123-2井累積降水量與累計(jì)水位升幅散點(diǎn)圖上(圖6b, d), 可以看到, 除了在累積降水量200~400 mm區(qū)間(占總數(shù)據(jù)量的0.08%)的散點(diǎn)與擬合直線有較明顯偏離外, 絕大部分點(diǎn)與直線擬合較好, 2井累積降水量與累計(jì)水位升幅相關(guān)系數(shù)分別為0.98和0.95(表2)。5—12月124井與同時(shí)期其他兩井相比, 累積降水量與累計(jì)水位升幅隨時(shí)間變化曲線更為一致, 散點(diǎn)圖最接近直線(圖6), 相關(guān)系數(shù)最高, 達(dá)0.99, 比相關(guān)系數(shù)最低的123-2井高出0.1(表2)。

    圖5 12個(gè)天文大潮期中潮汐與水位變化曲線

    注: 子圖中左上角字母a—i分別對(duì)應(yīng)1月8—9日、2月10—11日、3月28—29日、4月26—27日、5月25—26日、6月24—25日、7月23—24日、8月6—7日、9月22—23日、10月4—5日、11月18—19日與12月17—18日。

    圖6 2009年內(nèi)各觀測(cè)井累計(jì)水位升幅與累積降水量關(guān)系

    蒸發(fā)是自然環(huán)境下大氣回收水分的主要過(guò)程, 實(shí)際蒸發(fā)量與淺層地下水位的降低密不可分。2009年研究區(qū)實(shí)際蒸發(fā)量為687 mm, 略大于總降水量, 且明顯集中于夏季, 5—9月份實(shí)際蒸發(fā)量為468 mm, 占全年總蒸發(fā)量的68%。通過(guò)122-3井與123-2井各自累計(jì)水位降幅與累積蒸發(fā)量隨時(shí)間變化圖與相關(guān)關(guān)系圖(圖7a, b, c, d)可以發(fā)現(xiàn), 在2009年全年, 2條曲線比較貼合, 大部分散點(diǎn)與直線擬合較好, 相關(guān)系數(shù)均大于等于0.95。與降水對(duì)三口觀測(cè)井水位的影響規(guī)律類似, 在5—12月三井對(duì)比中, 124井水位下降與蒸發(fā)曲線貼合最緊密, 散點(diǎn)圖最接近擬合直線, 相關(guān)系數(shù)最高, 為0.99, 123-2井最低, 為0.90 (圖7, 表3)。122-3井與123-2井在5月至6月中旬水位下降速率高于蒸發(fā)速率, 之后至8月底低于蒸發(fā)速率(圖7a, c)。

    表2 2009年內(nèi)各觀測(cè)井累計(jì)水位升幅與累積降水量線性擬合結(jié)果

    圖7 2009年內(nèi)各觀測(cè)井累計(jì)水位降幅與累積蒸發(fā)量關(guān)系

    表3 2009年內(nèi)各觀測(cè)井累計(jì)水位降幅與累積蒸發(fā)量線性擬合結(jié)果

    3.2 討論

    3.2.1 沉積環(huán)境對(duì)淺層地下水位的影響

    124井與122-3井/123-2井水位整體動(dòng)態(tài)不同的主要原因是地理位置與沉積環(huán)境的差異。根據(jù)成國(guó)棟[17]的黃河三角洲現(xiàn)代沉積環(huán)境分區(qū)圖, 結(jié)合3個(gè)觀測(cè)井位的地理位置, 分析了3個(gè)井位所處的沉積環(huán)境。124井位于上三角洲平原, 周圍存在泛濫平原、決口扇和天然堤沉積相。122-3井處在上三角洲平原和下三角洲平原交界處, 周圍存在決口扇、河道以及潮灘等沉積相。123-2井則位于下三角洲平原的潮灘中, 離黃河口最近。124井周圍由于有防洪堤的保護(hù), 減弱了黃河泛濫的影響, 泛濫平原沉積物厚度較薄。該地沉積物主要以粗粉砂為主, 平均粒徑較大, 透水性較強(qiáng), 降水能有效入滲, 所以124井水位顯示出強(qiáng)烈的波動(dòng)特征(圖2)。下三角洲平原的潮灘表層沉積主要由粉砂與黏土質(zhì)粉砂組成, 粒徑小, 保水能力強(qiáng), 其下部沉積物以透水性差的黏土質(zhì)粉砂組成[6],致使潮灘區(qū)水分垂直運(yùn)動(dòng)差, 對(duì)降水等自然因素的響應(yīng)弱, 因此123-2井整體水位最為平穩(wěn)。位于上三角洲平原和下三角洲平原交界的122-3井周圍沉積物粒度總體介于124井和123-2井之間, 因此122-3井整體水位波動(dòng)程度明顯弱于124井水位, 但比123-2井水位稍明顯(圖2)。

    3.2.2 潮汐對(duì)淺層地下水位的影響

    黃河三角洲濱海濕地淺層地下水含水層之間連通性差, 淺層地下水水平運(yùn)動(dòng)能力差, 與潮汐沒(méi)有直接水力聯(lián)系[6]。但是潮汐作用于海岸時(shí)產(chǎn)生的沉積地層應(yīng)力效應(yīng)將通過(guò)含水層傳遞給淺層地下水, 使得淺層地下水位隨潮汐出現(xiàn)有規(guī)律的波動(dòng)[6, 8, 38-39]。Kim等[40]認(rèn)為潮汐對(duì)濟(jì)州島地下水位的水平影響距離可達(dá)3 km。利用黃河三角洲濱海濕地24 h的連續(xù)淺層地下水位和潮位數(shù)據(jù), 高茂生等[6]發(fā)現(xiàn), 距離驗(yàn)潮站7 km的水井淺層地下水位波動(dòng)與潮汐存在較強(qiáng)相關(guān)關(guān)系。楊再興[41]利用實(shí)測(cè)數(shù)據(jù)和軟件模擬, 認(rèn)為在北部灣地區(qū), 距離海岸3.7 km時(shí)地下水才不再明顯受潮汐影響。

    0.042 h–1和0.083 h–1的頻率分別對(duì)應(yīng)著24 h和12 h的周期。同時(shí)0.042 h–1頻率下的明顯高振幅, 表明研究區(qū)海潮24 h周期盛行(圖3a), 即研究區(qū)潮位具有近似全日潮的規(guī)律。122-3井與123-2井與潮汐數(shù)據(jù)高度相似的頻譜特征, 表明兩井水位均受潮汐波動(dòng)影響。二者相比, 122-3井頻率振幅更高(圖3b, c), 且12個(gè)天文大潮期中, 122-3井相對(duì)于123-2井, 對(duì)潮汐的響應(yīng)更早、水位波動(dòng)更大(圖5), 說(shuō)明122-3井更容易受到潮汐影響。頻率0.042 h–1和0.083 h–1在圖3d頻譜中未顯示明顯振幅表明124井水位波動(dòng)不存在某一特定周期, 受潮汐的調(diào)控微弱。122-3井、123-2井與124井距離海岸的最近垂直距離分別約為7 km、4 km和15 km。所以海潮對(duì)黃河三角洲濱海濕地淺層地下水位的水平調(diào)控范圍至少可達(dá)7 km。至15 km時(shí), 淺層地下水位幾乎不受潮汐影響。

    在前人有關(guān)滯后性的研究中, 某一井位淺層地下水位滯后潮汐的時(shí)間均為固定數(shù)值[6-8]。然而本文研究結(jié)果與前人不同。從時(shí)間滯后互相關(guān)分析結(jié)果中(圖4)可以看到, 潮位與平移后的水位相關(guān)系數(shù)恒小于0.03, 表現(xiàn)為不相關(guān), 這表明在全年時(shí)間尺度上, 相對(duì)于潮汐變化淺層地下水位波動(dòng)不存在特定滯后時(shí)間。圖5中可以明顯看出12個(gè)天文大潮期中潮汐超前水位的時(shí)間并不統(tǒng)一, 水位滯后“W”形和“M”形潮汐的時(shí)間相差約12 h。這種不統(tǒng)一的滯后響應(yīng), 可能受黃河三角洲年內(nèi)沉積環(huán)境、氣象條件、淺層地下水位以及潮汐變化的綜合影響。

    3.2.3 降水與蒸發(fā)對(duì)淺層地下水位的影響

    122-3井、123-2井和124井均能反映出水位上升與降水的密切聯(lián)系。在三個(gè)井位各自觀測(cè)時(shí)段內(nèi), 累計(jì)水位升幅與累積降水量相關(guān)系數(shù)均接近于1(表2),為極強(qiáng)相關(guān), 指征降水是三井水位上升的最重要的因素。但在圖6散點(diǎn)圖中, 122-3井與123-2井均有一部分散點(diǎn)與一元線性擬合直線偏離較大, 這種偏離在時(shí)間序列圖(圖6a, c)中也有體現(xiàn), 具體表現(xiàn)為: 5月9日與6月19日左右水位升幅相比于降水量的突增, 7月上旬平緩的水位增幅與快速累積的降水。這種偏離主要原因在于兩口井所在地沉積物的弱透水性。在粒徑小, 水力聯(lián)系弱的沉積環(huán)境下, 伴隨著長(zhǎng)期較低的降水量, 隨著蒸發(fā)量加大, 表層沉積物水分逐漸散失, 并且很難被側(cè)向徑流補(bǔ)給, 產(chǎn)生可以蓄水的大量干燥空間。5月9日與6月19日對(duì)應(yīng)著兩個(gè)較強(qiáng)的降水過(guò)程, 122-3井和123-2井所處的地層在此期間充分吸收水分, 產(chǎn)生水位升幅較大現(xiàn)象(圖2)。而在7月上旬, 雨季后, 沉積物孔隙中水分飽和, 水位也已恢復(fù)到正常水平(圖2), 此時(shí)降水已經(jīng)不能再使淺層地下水位產(chǎn)生較大的升幅??傮w而言, 降水對(duì)研究區(qū)三口觀測(cè)井水位上升起到了決定性作用。由于三口井所在區(qū)域沉積環(huán)境不同, 沉積物滲透性更好的124井淺層地下水位升高和降水關(guān)系最為密切, 位于潮灘地區(qū)的123-2井附近沉積物滲透性最差, 淺層地下水位抬升受降水影響相對(duì)最弱, 這也解釋了123-2井比122-3井離海洋更近, 但受潮汐影響卻更弱的現(xiàn)象。

    122-3井與123-2井在2009年全年中表現(xiàn)出了累計(jì)水位降幅與累積蒸發(fā)量之間極強(qiáng)的相關(guān)性(圖7a, b, c, d)。在5—12月三井的對(duì)比中: 124井相關(guān)性最高, 表現(xiàn)為累計(jì)水位降幅與累積蒸發(fā)量曲線高度擬合, 相關(guān)系數(shù)為0.99; 123-2井相關(guān)性最低, 表現(xiàn)為散點(diǎn)圖與直線偏離程度最大, 累計(jì)水位降幅與累積蒸發(fā)量偏離最大, 相關(guān)系數(shù)最低(圖7, 表3)。3井在各自數(shù)據(jù)時(shí)段內(nèi), 累積蒸發(fā)量與累計(jì)水位降幅接近于1的相關(guān)系數(shù)(表3), 體現(xiàn)了蒸發(fā)對(duì)水位下降的主導(dǎo)作用。122-3井與123-2井累計(jì)水位降幅與累積蒸發(fā)量曲線的偏離, 集中在5月至9月(圖7a, c)。這主要是由于5月至6月上旬黃河三角洲地區(qū)雨季尚未來(lái)臨, 降水仍保持較低水平, 但此期間氣溫快速回升, 蒸發(fā)量加大(圖2), 打破了之前降水與蒸發(fā)的平衡, 致使水位下降速率較之前增加, 這也是122-3井與123-2井在6月中上旬存在明顯水位低值(圖2)的原因。6月中旬以后, 研究區(qū)開(kāi)始受到東亞夏季風(fēng)影響, 降水突增, 地下水得到補(bǔ)充, 因此出現(xiàn)水位累計(jì)下降速率低于蒸發(fā)速率的情況。

    3.2.4 其他因素對(duì)淺層地下水位的影響

    黃河年內(nèi)徑流量與輸沙量變化明顯, 趨勢(shì)相同, 6—7月份最高, 8—12月份徑流輸沙量是1—5月份的數(shù)倍(圖8)。與黃河徑流量規(guī)律相反, 三口井水位在6—7月份處在最低水平。122-3井和123-2井8—12月份比1—5月份水位高, 但這主要是由7—9月份的高強(qiáng)度集中降水以及后續(xù)時(shí)間段實(shí)際蒸發(fā)量減少導(dǎo)致(圖2), 且水位在抬高程度上不及黃河徑流量在兩個(gè)時(shí)段上的差異顯著。同時(shí)研究區(qū)沉積物以粉砂和黏土質(zhì)粉砂為主, 淺層地下水水平運(yùn)動(dòng)差, 黃河和研究站位淺層地下水之間難以建立直接水力聯(lián)系。以上因素均表明黃河對(duì)于研究區(qū)淺層地下水位的影響有限。

    圖8 2009年內(nèi)黃河利津水文站月徑流量與輸沙量

    除了黃河徑流量, 農(nóng)業(yè)灌溉也是應(yīng)當(dāng)考慮的因素。研究區(qū)主要糧食作物是冬小麥。冬小麥在3、4月份處于拔節(jié)期, 這期間小麥生長(zhǎng)迅速, 需進(jìn)行引黃灌溉, 灌溉水逐漸下滲將會(huì)補(bǔ)充地下水位。所以在3—4月份, 盡管研究區(qū)蒸發(fā)量加大且降水量保持較低水平, 靠近農(nóng)田的122-3井水位仍呈現(xiàn)緩慢上升態(tài)勢(shì)。123-2井位于山東黃河三角洲國(guó)家級(jí)自然保護(hù)區(qū)內(nèi), 附近幾乎無(wú)農(nóng)業(yè)活動(dòng), 主要受自然因素調(diào)控, 因此3—4月保持著波動(dòng)下降狀態(tài)。

    4 結(jié)論

    本文以黃河三角洲2009年3口觀測(cè)井水位數(shù)據(jù)為基礎(chǔ), 結(jié)合該區(qū)域同時(shí)期潮位、降水量和實(shí)際蒸發(fā)量數(shù)據(jù), 深入研究了黃河三角洲淺層地下水位動(dòng)態(tài)變化特征及影響因素, 得出以下結(jié)論:

    1) 在黃河三角洲地區(qū), 不同的沉積環(huán)境造成了淺層地下水位動(dòng)態(tài)的差異。沉積物松散粒徑較大的區(qū)域淺層地下水位對(duì)降水響應(yīng)幅度大, 水位波動(dòng)明顯頻繁。臨近海洋的潮灘及附近區(qū)域, 沉積物透水性差, 淺層地下水位對(duì)外界環(huán)境改變響應(yīng)程度較弱, 但由于降水和蒸發(fā)的相互作用, 6月中上旬淺層地下水位急劇下降, 其余月份淺層地下水位呈現(xiàn)平穩(wěn)態(tài)勢(shì)。

    2) 研究區(qū)潮汐對(duì)淺層地下水位的影響平面上至少可達(dá)7 km, 離海岸15 km時(shí)潮汐對(duì)水位影響十分微弱。影響范圍內(nèi), 淺層地下水位與潮汐存在高度一致的頻譜, 均存在24 h與12 h的周期, 且24 h周期更為顯著。在滯后性研究中, 發(fā)現(xiàn)淺層地下水位對(duì)于潮汐的滯后時(shí)間存在變化這一新現(xiàn)象, 即滯后時(shí)間存在兩個(gè)區(qū)間, 每一個(gè)區(qū)間持續(xù)時(shí)間約0.5 a, 二者相差約12 h。

    3) 通過(guò)ERA5高分辨率降水量和實(shí)際蒸發(fā)量數(shù)據(jù), 定量分析了降水和蒸發(fā)對(duì)黃河三角洲淺層地下水位動(dòng)態(tài)的作用。淺層地下水位的升降分別與降水量、實(shí)際蒸發(fā)量存在極強(qiáng)相關(guān)關(guān)系, 降水和蒸發(fā)是導(dǎo)致黃河三角洲淺層地下水位升降最根本的因素。此外農(nóng)業(yè)灌溉也有一定影響。

    [1] Kuenzer C, Ottinger M, Liu G, et al.Earth observation-based coastal zone monitoring of the Yellow River Delta: Dynamics in China’s second largest oil producing region over four decades[J].Applied Geo-graphy, 2014, 55: 92-107.

    [2] Zhang J, Huang H, Bi H.Land subsidence in the modern Yellow River Delta based on InSAR time series analysis[J].Natural Hazards, 2015, 75(3): 2385-2397.

    [3] Fan X, Pedroli B, Liu G, et al.Soil salinity develop-ment in the Yellow River Delta in relation to ground-water dynamics[J].Land Degradation & Development, 2012, 23(2): 175-189.

    [4] 劉勇, 徐興永, 李培英, 等.黃河三角洲地下水動(dòng)態(tài)變化及其環(huán)境災(zāi)害效應(yīng)研究[C]//中國(guó)環(huán)境科學(xué)學(xué)會(huì).2016中國(guó)環(huán)境科學(xué)學(xué)會(huì)學(xué)術(shù)年會(huì)論文集(第三卷).北京: 中國(guó)環(huán)境科學(xué)出版社, 2016: 1240-1248.

    Liu Yong, Xu Xingyong, Li Peiying, et al.Dynamics of groundwater and its environmental disasters influen-ces research in the Yellow River Delta[C]//Chinese Society for Environmental Sciences.2016 Annual Con-ference of Chinese Society for Environmental Sciences (Volume 3).Beijing: China Environmental Science Press, 2016: 1240-1248.

    [5] 安樂(lè)生, 趙全升, 許穎.黃河三角洲淺層地下水位動(dòng)態(tài)特征及其成因[J].環(huán)境科學(xué)與技術(shù), 2013, 36(9): 51-56.

    An Lesheng, Zhao Quansheng, Xu Ying.Dynamic characteristics of the shallow groundwater table and its genesis in the Yellow River Delta[J].Environment Science & Technology, 2013, 36(9): 51-56.

    [6] 高茂生, 葉思源, 史貴軍, 等.潮汐作用下的濱海濕地淺層地下水動(dòng)態(tài)變化[J].水文地質(zhì)工程地質(zhì), 2010, 37(4): 24-27.

    Gao Maosheng, Ye Siyuan, Shi Guijun, et al.Oceanic tide-induced shallow groundwater regime fluctuations in coastal wetland[J].Hydrology & Engineering Geo-logy, 2010, 37(4): 24-27.

    [7] 蘇喬, 徐興永, 陳廣泉, 等.潮汐作用對(duì)萊州灣南岸冬季地下水位的影響研究[J].海洋科學(xué)進(jìn)展, 2017, 35(4): 568-578.

    Su Qiao, Xu Xingyong, Chen Guangquan, et al.Study on the impact of tides on groundwater table fluctuation in coastal aquifer[J].Advances in Marine Science, 2017, 35(4): 568-578.

    [8] 蘇喬, 徐興永, 陳廣泉, 等.潮汐作用影響下濱海地區(qū)地下水位的變化頻率和滯后性[J].海洋開(kāi)發(fā)與管理, 2018, 35(10): 81-85.

    Su Qiao, Xu Xingyong, Chen Guangquan, et al.Fre-quency and hysteresis of groundwater level influenced by tides[J].Ocean Development and Management, 2018, 35(10): 81-85.

    [9] Zhang X, Dong F, Dai H, et al.Influence of lunar semidiurnal tides on groundwater dynamics in estuarine aquifers[J].Hydrogeology Journal, 2020, 28: 1419-1429.

    [10] 袁瑞強(qiáng), 劉貫群, 宋獻(xiàn)方.現(xiàn)代黃河三角洲淺層地下水對(duì)降水的響應(yīng)[J].資源科學(xué), 2009, 31(9): 1514-1521.

    Yuan Ruiqiang, Liu Guanqun, Song Xianfang.Shallow groundwater response to rainfall infiltration in modern Yellow River Delta[J].Resources Science, 2009, 31(9): 1514-1521.

    [11] Schulthess U, Ahmed Z U, Aravindakshan S, et al.Farming on the fringe: Shallow groundwater dyna-mics and irrigation scheduling for maize and wheat in Bangladesh’s coastal delta[J].Field Crops Research, 2019, 239: 135-148.

    [12] 張國(guó)華, 郭磊, 劉健, 等.山東省陸面實(shí)際蒸發(fā)量估算及變化特征分析[J].人民黃河, 2014, 36(10): 26-29.

    Zhang Guohua, Guo Lei, Liu Jian, et al.Estimation and variation characteristics of land actual evapo-trans-piration in Shandong Province[J].Yellow River, 2014, 36(10): 26-29.

    [13] Minderhoud P S J, Erkens G, Pham V H, et al.Impacts of 25 years of groundwater extraction on subsidence in the Mekong delta, Vietnam[J].Environ-mental Research Letters, 2017, 12(6): 064006.

    [14] Killian C D, Asquith W H, Barlow J R B, et al.Characterizing groundwater and surface-water interac-tion using hydrograph-separation techniques and groundwater- level data throughout the Mississippi Delta[J].Hydrogeo-logy Journal, 2019, 27(6): 2167-2179.

    [15] Kennedy J, Rodríguez-Burgue?o J E, Ramírez-Hernández J.Groundwater response to the 2014 pulse flow in the Colorado River Delta[J].Ecological Engineering, 2017, 106: 715-724.

    [16] 劉一霖.黃河三角洲地面沉降時(shí)序InSAR技術(shù)監(jiān)測(cè)與地下流體開(kāi)采相關(guān)性分析[D].青島: 中國(guó)科學(xué)院海洋研究所, 2016.

    Liu Yilin.Land subsidence monitoring and correlation analysis with ground fluid exploitation using time series INSAR technique over the Yellow River Delta (China)[D].Qingdao: Institute of Oceanology, Chinese Academy of Sciences, 2016.

    [17] 成國(guó)棟.黃河三角洲現(xiàn)代沉積作用及模式[M].北京: 地質(zhì)出版社, 1991: 39-49.

    CHENG Guodong.Modern sedimentation and patterns in the Yellow River Delta[M].Beijing: Geology Press, 1991: 39-49.

    [18] 成國(guó)棟.黃河三角洲沉積地質(zhì)學(xué)[M].北京: 地質(zhì)出版社, 1997: 3-72.

    CHENG Guodong.Sedimentary geology in the Yellow River Delta[M].Beijing: Geology Press, 1997: 3-72.

    [19] 喬淑卿, 石學(xué)法.黃河三角洲沉積特征和演化研究現(xiàn)狀及展望[J].海洋科學(xué)進(jìn)展, 2010, 28(3): 408-416.

    Qiao Shuqing, Shi Xuefa.Status and prospect of studies on sedimentary characteristics and evolution of the Yellow River Delta[J].Advances in Marine Science, 2010, 28(3): 408-416.

    [20] 李高偉.黃河三角洲氣候變化及其對(duì)濕地影響研究[D].濟(jì)南: 山東師范大學(xué), 2015.

    Li Gaowei.Research of climate change in the Yellow River Delta and its impacts on wetland[D].Jinan: Shandong Normal University, 2015.

    [21] 李向陽(yáng), 陳沈良, 李為華, 等.黃河三角洲孤東海域沖淤塑造營(yíng)力研究[J].人民黃河, 2008, 30(10): 35-37.

    Li Xiangyang, Chen Shenliang, Li Weihua, et al.Study on scour and fill shape of Gudong sea area of the Yellow River Delta[J].Yellow River, 2008, 30(10): 35-37.

    [22] 陳一強(qiáng), 陳沈良, 馬宗文, 等.黃河三角洲孤東近岸懸沙質(zhì)量濃度變化和底沙再懸浮[J].海洋科學(xué)進(jìn)展, 2011, 29(2): 157-162.

    Chen Yiqiang, Chen Shenliang, Ma Zongwen, et al.Variations of suspended sediment concentration and resuspension of bottom sediments in the Gudong near-shore area of the Huanghe River Delta[J].Advances in Marine Science, 2011, 29(2): 157-162.

    [23] 程慧.近40年來(lái)黃河三角洲孤東近岸的沖淤演變及其影響因素[D].上海: 華東師范大學(xué), 2019.

    Cheng Hui.Coastal evolution and its influencing factors in Gudong nearshore area of the Yellow River Delta in recent 40 years[D].Shanghai: East China Normal University, 2019.

    [24] Wang G, Li P, Li Z, et al.Coastal dam inundation assessment for the Yellow River Delta: Measurements, analysis and scenario[J].Remote Sensing, 2020, 12(21): 3658.

    [25] Hersbach H, Bell B, Berrisford P, et al.The ERA5 global reanalysis[J].Quarterly Journal of the Royal Meteorological Society, 2020, 146(730): 1999- 2049.

    [26] Xue C, Wu H, Jiang X.Temporal and spatial change monitoring of drought grade based on ERA5 analysis data and BFAST method in the belt and road area during 1989-2017[J].Advances in Meteorology, 2019, 2019(428): 1-10.

    [27] 黃建平, 張國(guó)龍, 于海鵬, 等.黃河流域近40年氣候變化的時(shí)空特征[J].水利學(xué)報(bào), 2020, 51(9): 1048-1058.

    Huang Jianping, Zhang Guolong, Yu Haipeng, et al.Characteristics of climate change in the Yellow River basin during recent 40 years[J].Shuili Xuebao, 2020, 51(9): 1048-1058.

    [28] Xu X, Frey S K, Boluwade A, et al.Evaluation of variability among different precipitation products in the Northern Great Plains[J].Journal of Hydrology: Regio-nal Studies, 2019, 24: 100608.

    [29] Tarek M, Brissette F P, Arsenault R.Evalua-tion of the ERA5 reanalysis as a potential reference dataset for hydrological modeling over North Ame-rica[J].Hydrology and Earth System Sciences, 2020, 24(5): 2527-2544.

    [30] Olauson J.ERA5: The new champion of wind power modelling?[J].Renewable Energy, 2018, 126: 322-331.

    [31] Belmonte Rivas M, Stoffelen A.Charac-teri-zing ERA-Interim and ERA5 surface wind biases using ASCAT[J].Ocean Science, 2019, 15(3): 831-852.

    [32] Kalverla P C, Duncan Jr J B, Steeneveld G J, et al.Low-level jets over the North Sea based on ERA5 and observations: together they do better[J].Wind Energy Science, 2019, 4(2): 193-209.

    [33] Li Z, Yan Z, Zhu Y, et al.Homogenized daily relative humidity series in China during 1960–2017[J].Advan-ces in Atmospheric Sciences, 2020, 37(4): 318-327.

    [34] Yao B, Liu C, Yin Y, et al.Evaluation of cloud pro-perties from reanalyses over East Asia with a radiance- based approach[J].Atmospheric Measurement Techni-ques, 2020, 13(3): 1033-1049.

    [35] 孫赫, 蘇鳳閣, 黃敬恒, 等.第三極西風(fēng)和季風(fēng)主導(dǎo)流域源區(qū)降水呈現(xiàn)不同梯度特征[J].科學(xué)通報(bào), 2020, 65(1): 91-104.

    Sun He, Su Fengge, Huang Jingheng, et al.Con-trasting precipitation gradient characteristics between westerlies and monsoon dominated upstream river basins in the Third Pole[J].Chinese Science Bulletin, 2020, 65(1): 91-104.

    [36] 孟憲貴, 郭俊建, 韓永清.ERA5再分析數(shù)據(jù)適用性初步評(píng)估[J].海洋氣象學(xué)報(bào), 2018, 38(1): 91-99.

    Meng Xiangui, Guo Junjian, Han Yongqing.Preli-minarily assessment of ERA5 reanalysis data[J].Jour-nal of Marine Meteorology, 2018, 38(1): 91-99.

    [37] 陳崇希.地下水不穩(wěn)定井流計(jì)算方法[M].北京: 地質(zhì)出版社, 1990: 18-19.

    CHEN Chongxi.The calculation method of unstable well flow of groundwater[M].Beijing: Geology Press, 1990: 18-19.

    [38] Kim J H, Lee J, Cheong T J, et al.Use of time series analysis for the identification of tidal effect on groundwater in the coastal area of Kimje, Korea[J].Journal of Hydrology, 2005, 300(1/4): 188-198.

    [39] 劉艷.海潮引起的濱海地區(qū)地下水位波動(dòng)的研究[D].北京: 中國(guó)地質(zhì)大學(xué)(北京), 2016.

    Liu Yan.Tide-induced groundwater level fluctuations in coastal areas[D].Beijing: China University of Geos-ciences (Beijing), 2016.

    [40] Kim K Y, Seong H, Kim T, et al.Tidal effects on variations of fresh-saltwater interface and groundwater flow in a multilayered coastal aquifer on a volcanic island (Jeju Island, Korea)[J].Journal of Hydrology, 2006, 330(3/4): 525-542.

    [41] 楊再興.海潮波動(dòng)對(duì)濱海地區(qū)地下水位的影響[D].北京: 中國(guó)地質(zhì)大學(xué)(北京), 2011.

    Yang Zaixing.Effect of tidal fluctuations on ground water levels in coastal areas[D].Beijing: China Uni-versity of Geosciences (Beijing), 2011.

    Dynamic characteristics and causes of shallow groundwater level intra-annual changes in the Yellow River Delta

    CHANG Mao-xiang1, 2, SHI Jing-hao3, YE Si-yuan4, WANG Xiang-dong3, LI Guang-xue1, 2

    (1.Key Lab of Submarine Geosciences and Prospecting Techniques (Ocean University of China), Ministry of Education, Qingdao 266100, China; 2.College of Marine Geosciences, Ocean University of China, Qingdao 266100, China; 3.Qingdao Blue Earth Big Data Technology Co., Ltd, Qingdao 266400, China; 4.Qingdao Institute of Marine Geology, Ministry of Natural Resources of the People’s Republic of China, Qingdao 266071, China)

    the Yellow River Delta; shallow groundwater level; tide; hysteresis; precipitation; actual evaporation

    In this work, the characteristics and mechanics of shallow groundwater level fluctuations from continuous monitoring records of three monitoring wells in the Yellow River Delta and the influence of tide, precipitation, and runoff are analyzed using hourly shallow groundwater level data and contemporarily high-resolution ERA5 reanalysis actual evaporation data, adopting methods including fast Fourier transform and time series analysis.Observations showed that (1) the type of sedimentary environment leads to the difference in general conditions of the groundwater level.Surface sediment with relatively good permeability displays shallow groundwater level fluc-tuation that is more violent and frequent; whereas, surface sediment with relatively weak permeability displays a shallow groundwater level that reaches a striking minimum in early June but displays a stable condition in general.(2) Horizontal scope of the tidal influence on the shallow groundwater level in the Yellow River Delta could reach 7 km but not exceed 15 km.Within that scope, the time lag between the shallow groundwater level fluctuation and tide in 2009 is not constant—it can be divided into two 6-month intervals with a 12 h difference in between.(3) Precipitation and actual evaporation are the fundamental factors that affect the rise and fall of shallow groundwater in the Yellow River Delta.Agricultural activities also have a certain effect on the shallow groundwater level.Research on shallow groundwater level dynamics could provide a scientific basis for ecological protection, prevention, and control of local land salinization and seawater intrusion.

    Nov.10, 2020

    [The National Key Research and Development Program of China, No.2017YFE0133500; The Taishan Scholar Project grant to GX Li]

    P736.5

    A

    1000-3096(2021)10-0020-12

    10.11759/hykx20201110007

    2020-11-10;

    2021-03-02

    國(guó)家重點(diǎn)研發(fā)計(jì)劃政府間重點(diǎn)專項(xiàng)項(xiàng)目(2017YFE0133500); 泰山學(xué)者建設(shè)工程專項(xiàng)

    常茂祥(1996—), 山東煙臺(tái)人, 碩士研究生, 主要從事海洋地質(zhì)與遙感研究, E-mail: 915995011@qq.com; 史經(jīng)昊(1981—), 通信作者, 博士, 主要從事近海綜合環(huán)境調(diào)查和監(jiān)測(cè)工作, 以及沉積動(dòng)力學(xué)和海洋工程方面的研究, E-mail: 34098302@qq.com

    (本文編輯: 趙衛(wèi)紅)

    猜你喜歡
    黃河三角洲蒸發(fā)量井水
    井水的奧秘
    黃河三角洲保護(hù)區(qū)自然資源的開(kāi)發(fā)與保護(hù)
    1958—2013年沽源縣蒸發(fā)量變化特征分析
    1981—2010年菏澤市定陶區(qū)蒸發(fā)量變化特征分析
    新疆民豐縣地表水面蒸發(fā)量分析
    達(dá)孜縣夏秋季大小型蒸發(fā)量特征、影響因子與差異分析
    地球(2016年7期)2016-08-23 03:01:35
    盤一井水氡異常落實(shí)與分析
    黃驊埕古1井水位異常研究
    黃河三角洲濱海濕地維管植物多樣性現(xiàn)狀及保護(hù)策略
    城鎮(zhèn)化與生態(tài)環(huán)境交互協(xié)調(diào)行為研究——以黃河三角洲為例
    中文天堂在线官网| 寂寞人妻少妇视频99o| 99视频精品全部免费 在线| 国产精品.久久久| 老司机影院毛片| 国产探花极品一区二区| 波野结衣二区三区在线| 精品国产一区二区三区久久久樱花 | 欧美97在线视频| 天堂av国产一区二区熟女人妻| 日韩成人伦理影院| 少妇高潮的动态图| h日本视频在线播放| 大话2 男鬼变身卡| 纵有疾风起免费观看全集完整版 | 午夜亚洲福利在线播放| 欧美成人一区二区免费高清观看| 好男人在线观看高清免费视频| 久久人人爽人人片av| 亚洲av在线观看美女高潮| 精品一区二区三卡| 高清视频免费观看一区二区 | 亚洲国产精品成人综合色| 最后的刺客免费高清国语| 国产精品美女特级片免费视频播放器| 欧美3d第一页| 中文在线观看免费www的网站| 可以在线观看毛片的网站| av黄色大香蕉| 十八禁国产超污无遮挡网站| 一级毛片aaaaaa免费看小| 在线播放无遮挡| 高清视频免费观看一区二区 | 午夜老司机福利剧场| 国产黄色视频一区二区在线观看| 国产单亲对白刺激| 丝袜喷水一区| 亚洲精华国产精华液的使用体验| 精品少妇黑人巨大在线播放| 国产高清不卡午夜福利| 在线观看av片永久免费下载| 久久99精品国语久久久| 草草在线视频免费看| 国产成人精品婷婷| 深爱激情五月婷婷| 亚洲国产最新在线播放| 国产老妇伦熟女老妇高清| 久久久久性生活片| 亚洲av电影不卡..在线观看| 婷婷色av中文字幕| 天堂影院成人在线观看| 久久久精品免费免费高清| 国产伦理片在线播放av一区| 国产精品一区二区在线观看99 | 777米奇影视久久| 欧美性猛交╳xxx乱大交人| 国产精品久久视频播放| 国产视频首页在线观看| 国产精品日韩av在线免费观看| 日日干狠狠操夜夜爽| 国产亚洲av嫩草精品影院| 久久久久性生活片| 尤物成人国产欧美一区二区三区| 午夜免费观看性视频| 亚洲欧美日韩东京热| 欧美三级亚洲精品| 日本与韩国留学比较| 亚洲av男天堂| 国产高清国产精品国产三级 | 久久99蜜桃精品久久| 久久久久久伊人网av| 我要看日韩黄色一级片| 能在线免费看毛片的网站| 男人和女人高潮做爰伦理| 成人二区视频| 亚洲在线自拍视频| 十八禁网站网址无遮挡 | 亚洲综合色惰| 欧美日韩视频高清一区二区三区二| 亚洲综合精品二区| 久久久久九九精品影院| 国产精品久久久久久精品电影小说 | 国产探花在线观看一区二区| 网址你懂的国产日韩在线| 亚洲精品久久久久久婷婷小说| av在线天堂中文字幕| 又爽又黄无遮挡网站| 国产一区二区三区综合在线观看 | 只有这里有精品99| 天天躁日日操中文字幕| 夜夜爽夜夜爽视频| 只有这里有精品99| 日日干狠狠操夜夜爽| 欧美成人午夜免费资源| 欧美三级亚洲精品| 青春草亚洲视频在线观看| 国产亚洲一区二区精品| 在线a可以看的网站| 七月丁香在线播放| 精品一区二区三区视频在线| 免费av观看视频| 亚洲自拍偷在线| 最近中文字幕高清免费大全6| 午夜福利成人在线免费观看| 免费观看av网站的网址| 丰满乱子伦码专区| 久久精品久久久久久久性| 免费看光身美女| 男人舔女人下体高潮全视频| 亚洲婷婷狠狠爱综合网| 午夜免费观看性视频| 免费高清在线观看视频在线观看| 少妇猛男粗大的猛烈进出视频 | 国产精品久久久久久久电影| 乱人视频在线观看| 久久久久久久国产电影| 国产精品1区2区在线观看.| 久久久久久久久久黄片| 日本黄色片子视频| 三级经典国产精品| 一级片'在线观看视频| 国产精品美女特级片免费视频播放器| 黄色一级大片看看| 国产国拍精品亚洲av在线观看| 亚洲欧美一区二区三区国产| 内射极品少妇av片p| 成人国产麻豆网| 最近手机中文字幕大全| 亚洲国产日韩欧美精品在线观看| 亚洲av成人精品一二三区| 久久久久九九精品影院| 18+在线观看网站| 亚洲精品日韩av片在线观看| 亚洲最大成人手机在线| 国产午夜精品论理片| 精品久久久久久久久亚洲| 亚洲高清免费不卡视频| 国产在线男女| 日本色播在线视频| 亚洲精品国产av成人精品| 美女被艹到高潮喷水动态| 国产三级在线视频| 亚洲综合精品二区| 免费av观看视频| 26uuu在线亚洲综合色| 性色avwww在线观看| 成年女人在线观看亚洲视频 | 久久综合国产亚洲精品| 老女人水多毛片| 日韩欧美 国产精品| 青春草亚洲视频在线观看| 成年av动漫网址| 网址你懂的国产日韩在线| av免费观看日本| 淫秽高清视频在线观看| 国产精品久久视频播放| 欧美 日韩 精品 国产| 久久久精品免费免费高清| 成年女人看的毛片在线观看| 高清日韩中文字幕在线| 欧美zozozo另类| 久久99热6这里只有精品| 欧美潮喷喷水| 大话2 男鬼变身卡| 国产成人免费观看mmmm| 天美传媒精品一区二区| 国产中年淑女户外野战色| 国产一区有黄有色的免费视频 | 美女国产视频在线观看| 色播亚洲综合网| 91精品国产九色| 国产午夜精品论理片| 国产黄色免费在线视频| 2018国产大陆天天弄谢| 春色校园在线视频观看| 亚洲精品乱码久久久v下载方式| 熟妇人妻久久中文字幕3abv| 日日啪夜夜撸| 久久这里只有精品中国| 日韩欧美三级三区| 91精品伊人久久大香线蕉| 婷婷色av中文字幕| 亚洲国产日韩欧美精品在线观看| 男女边摸边吃奶| 五月玫瑰六月丁香| 午夜老司机福利剧场| 别揉我奶头 嗯啊视频| 国产一区二区三区综合在线观看 | 成人av在线播放网站| 国产精品麻豆人妻色哟哟久久 | 最近手机中文字幕大全| 精品亚洲乱码少妇综合久久| 内射极品少妇av片p| av在线观看视频网站免费| 激情五月婷婷亚洲| 能在线免费观看的黄片| 国产高清三级在线| 成年版毛片免费区| 婷婷色综合大香蕉| 天美传媒精品一区二区| 麻豆乱淫一区二区| 少妇的逼水好多| 91久久精品电影网| 亚洲av成人av| 在线 av 中文字幕| 夜夜爽夜夜爽视频| 亚洲精品aⅴ在线观看| 一边亲一边摸免费视频| 日韩亚洲欧美综合| 亚洲在线自拍视频| 三级国产精品欧美在线观看| 亚洲一级一片aⅴ在线观看| 丰满少妇做爰视频| 中文字幕人妻熟人妻熟丝袜美| 免费av不卡在线播放| 日韩精品有码人妻一区| 卡戴珊不雅视频在线播放| 天堂俺去俺来也www色官网 | 欧美日韩精品成人综合77777| 在线观看免费高清a一片| 免费看日本二区| 精品少妇黑人巨大在线播放| 国产 一区精品| 一夜夜www| 亚洲国产成人一精品久久久| 不卡视频在线观看欧美| 久久久久精品久久久久真实原创| 久久97久久精品| 久热久热在线精品观看| 国内精品一区二区在线观看| 3wmmmm亚洲av在线观看| 免费观看无遮挡的男女| 九九爱精品视频在线观看| .国产精品久久| 国产精品一及| 午夜福利视频1000在线观看| 亚洲av成人精品一区久久| 大话2 男鬼变身卡| 日韩av免费高清视频| 能在线免费观看的黄片| 国产午夜精品一二区理论片| 亚洲精品视频女| or卡值多少钱| 成年免费大片在线观看| 一个人看视频在线观看www免费| 日韩欧美精品v在线| 日韩成人av中文字幕在线观看| 亚洲精品456在线播放app| 一个人看视频在线观看www免费| 黄片wwwwww| 一区二区三区免费毛片| 一级爰片在线观看| 久久久久久九九精品二区国产| 91久久精品国产一区二区三区| 在线免费观看不下载黄p国产| 久久亚洲国产成人精品v| 深夜a级毛片| 国精品久久久久久国模美| 丝瓜视频免费看黄片| 99视频精品全部免费 在线| av又黄又爽大尺度在线免费看| 亚洲国产色片| 91精品一卡2卡3卡4卡| 国产一区二区三区av在线| 欧美激情在线99| 久久久久久久久久久丰满| 精品国产露脸久久av麻豆 | 国产麻豆成人av免费视频| 秋霞伦理黄片| 国产女主播在线喷水免费视频网站 | 国产精品熟女久久久久浪| 欧美变态另类bdsm刘玥| 国产真实伦视频高清在线观看| 一区二区三区四区激情视频| 欧美3d第一页| 久久鲁丝午夜福利片| 欧美bdsm另类| 久久6这里有精品| 免费看日本二区| 看免费成人av毛片| 一级av片app| 精品久久久久久久末码| 欧美xxxx黑人xx丫x性爽| 国产高清三级在线| 亚洲精品亚洲一区二区| 亚洲aⅴ乱码一区二区在线播放| 嘟嘟电影网在线观看| 成人亚洲精品一区在线观看 | 午夜激情福利司机影院| 老司机影院毛片| 美女xxoo啪啪120秒动态图| 极品少妇高潮喷水抽搐| 国产黄色免费在线视频| 国产精品久久视频播放| 亚洲丝袜综合中文字幕| 国产黄色免费在线视频| 99久久九九国产精品国产免费| 欧美bdsm另类| 最近手机中文字幕大全| 少妇熟女aⅴ在线视频| 九九久久精品国产亚洲av麻豆| 欧美成人a在线观看| 午夜久久久久精精品| freevideosex欧美| 国产亚洲精品久久久com| 尾随美女入室| 日韩伦理黄色片| 国产免费视频播放在线视频 | 日韩不卡一区二区三区视频在线| 2021少妇久久久久久久久久久| 少妇的逼水好多| 男女下面进入的视频免费午夜| 一本久久精品| 午夜视频国产福利| 欧美bdsm另类| 久久久久网色| 97超碰精品成人国产| 一级片'在线观看视频| 国产 一区 欧美 日韩| 搡女人真爽免费视频火全软件| 婷婷色av中文字幕| 国产三级在线视频| 超碰97精品在线观看| 婷婷色麻豆天堂久久| 久久久久精品性色| 国产乱人视频| 麻豆成人午夜福利视频| 一区二区三区高清视频在线| 乱系列少妇在线播放| 人妻夜夜爽99麻豆av| 久久国内精品自在自线图片| 你懂的网址亚洲精品在线观看| 国产成人a∨麻豆精品| 伊人久久精品亚洲午夜| 日本爱情动作片www.在线观看| 在线播放无遮挡| av一本久久久久| 日韩欧美一区视频在线观看 | 一个人看视频在线观看www免费| 青春草亚洲视频在线观看| 亚洲人成网站在线观看播放| 欧美另类一区| 国产又色又爽无遮挡免| 美女xxoo啪啪120秒动态图| 免费看不卡的av| 美女黄网站色视频| 超碰97精品在线观看| 成年版毛片免费区| 九色成人免费人妻av| 亚洲av不卡在线观看| 纵有疾风起免费观看全集完整版 | 高清在线视频一区二区三区| 久久99蜜桃精品久久| 丰满少妇做爰视频| 日韩欧美一区视频在线观看 | 在线观看美女被高潮喷水网站| 高清在线视频一区二区三区| 国产av在哪里看| av卡一久久| 在线播放无遮挡| 久久久久久久久久黄片| 一级黄片播放器| 国产伦在线观看视频一区| 日本午夜av视频| 少妇裸体淫交视频免费看高清| 99热这里只有是精品50| 人人妻人人澡欧美一区二区| 18禁在线无遮挡免费观看视频| 国产一区有黄有色的免费视频 | 中文精品一卡2卡3卡4更新| 免费看日本二区| 在线观看av片永久免费下载| 久久人人爽人人片av| 汤姆久久久久久久影院中文字幕 | 国产午夜福利久久久久久| 婷婷色综合www| 日本免费在线观看一区| 日韩亚洲欧美综合| 久久久亚洲精品成人影院| 51国产日韩欧美| 亚洲精品乱码久久久久久按摩| 少妇人妻精品综合一区二区| 日本与韩国留学比较| 女人被狂操c到高潮| 男的添女的下面高潮视频| 国产成人a∨麻豆精品| 人妻一区二区av| 亚州av有码| 亚洲精品456在线播放app| 亚洲精品aⅴ在线观看| 18禁在线无遮挡免费观看视频| 久久亚洲国产成人精品v| 国产成年人精品一区二区| 亚洲av电影不卡..在线观看| 男人舔奶头视频| 国产在视频线精品| kizo精华| 爱豆传媒免费全集在线观看| 女的被弄到高潮叫床怎么办| 欧美精品一区二区大全| 一区二区三区免费毛片| 国产成人freesex在线| 淫秽高清视频在线观看| 欧美激情国产日韩精品一区| 欧美zozozo另类| .国产精品久久| 亚洲va在线va天堂va国产| 欧美xxⅹ黑人| 狂野欧美白嫩少妇大欣赏| 神马国产精品三级电影在线观看| 又爽又黄无遮挡网站| 少妇被粗大猛烈的视频| 又爽又黄a免费视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲av成人精品一区久久| 欧美bdsm另类| 亚洲成人一二三区av| 日日啪夜夜撸| 国产精品一区www在线观看| 国产精品1区2区在线观看.| 久久久精品欧美日韩精品| 国产免费又黄又爽又色| 热99在线观看视频| 亚洲成人一二三区av| 欧美精品一区二区大全| 亚洲欧美清纯卡通| 久久久久久久久久久免费av| 18禁在线播放成人免费| 少妇的逼水好多| 午夜激情久久久久久久| 亚洲人成网站在线播| 免费高清在线观看视频在线观看| 麻豆乱淫一区二区| 22中文网久久字幕| 人妻少妇偷人精品九色| 大话2 男鬼变身卡| 春色校园在线视频观看| 成人亚洲精品一区在线观看 | 国产探花在线观看一区二区| 亚洲av成人精品一二三区| 国产一级毛片在线| 亚洲成人中文字幕在线播放| 精品人妻一区二区三区麻豆| 国产成人精品婷婷| 蜜臀久久99精品久久宅男| av在线天堂中文字幕| 国产伦精品一区二区三区四那| 国产精品一区二区三区四区久久| 少妇的逼水好多| 一个人免费在线观看电影| 黄色日韩在线| 好男人在线观看高清免费视频| 亚洲精品日韩在线中文字幕| 久久久久久久大尺度免费视频| 久久久久久久久大av| 日韩电影二区| 国产综合精华液| 国产精品一区二区在线观看99 | 成年版毛片免费区| 亚洲精品日本国产第一区| 亚洲精品久久久久久婷婷小说| 18禁动态无遮挡网站| 一级毛片 在线播放| 国产伦精品一区二区三区四那| 青青草视频在线视频观看| 国产精品1区2区在线观看.| 在线 av 中文字幕| 国产av在哪里看| 免费在线观看成人毛片| 亚洲精品日韩在线中文字幕| 久久精品夜色国产| 国产精品综合久久久久久久免费| 一本一本综合久久| 国内精品美女久久久久久| 69av精品久久久久久| 亚洲欧美日韩东京热| 国产精品女同一区二区软件| 久久亚洲国产成人精品v| 亚洲最大成人av| 秋霞伦理黄片| 少妇被粗大猛烈的视频| 99热这里只有精品一区| 国内揄拍国产精品人妻在线| 蜜桃亚洲精品一区二区三区| 日日啪夜夜撸| 内地一区二区视频在线| 国产国拍精品亚洲av在线观看| 国产午夜精品论理片| 69人妻影院| 国内精品一区二区在线观看| 国产高清有码在线观看视频| 熟妇人妻久久中文字幕3abv| 日韩强制内射视频| 国产又色又爽无遮挡免| 国产精品爽爽va在线观看网站| 秋霞伦理黄片| 国产亚洲午夜精品一区二区久久 | 精品一区二区三卡| 综合色丁香网| 日本免费在线观看一区| av女优亚洲男人天堂| 日韩av在线免费看完整版不卡| 三级毛片av免费| 免费看a级黄色片| 成人鲁丝片一二三区免费| 亚洲欧美精品专区久久| 国产单亲对白刺激| 看非洲黑人一级黄片| 亚洲av电影不卡..在线观看| 黄片无遮挡物在线观看| 成人二区视频| 亚洲av成人精品一二三区| 成人一区二区视频在线观看| 日本av手机在线免费观看| 夜夜看夜夜爽夜夜摸| 麻豆精品久久久久久蜜桃| 色综合亚洲欧美另类图片| 99久久精品一区二区三区| 女人十人毛片免费观看3o分钟| 精品一区二区三区人妻视频| 又爽又黄无遮挡网站| 免费不卡的大黄色大毛片视频在线观看 | 免费播放大片免费观看视频在线观看| av专区在线播放| 一二三四中文在线观看免费高清| 日韩欧美 国产精品| 中文字幕制服av| 亚洲av福利一区| 免费看光身美女| 日韩亚洲欧美综合| 亚洲成人中文字幕在线播放| 亚洲av电影不卡..在线观看| www.色视频.com| 99热这里只有是精品50| 久久久国产一区二区| 成人综合一区亚洲| 激情五月婷婷亚洲| 老司机影院毛片| 久久久久性生活片| 亚洲一级一片aⅴ在线观看| 亚洲性久久影院| 大香蕉97超碰在线| av国产免费在线观看| 大又大粗又爽又黄少妇毛片口| 久久99热这里只频精品6学生| 国内精品美女久久久久久| 国产精品蜜桃在线观看| 久久人人爽人人爽人人片va| 精品酒店卫生间| 最近2019中文字幕mv第一页| 亚洲国产精品专区欧美| 人妻一区二区av| 亚洲四区av| av播播在线观看一区| 免费黄网站久久成人精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 免费黄色在线免费观看| 日韩一本色道免费dvd| 国产精品国产三级国产专区5o| 日本午夜av视频| 国产淫语在线视频| 国产精品久久久久久av不卡| 日日啪夜夜撸| 国产一区有黄有色的免费视频 | 中文字幕av在线有码专区| 丰满少妇做爰视频| 日韩欧美精品v在线| 亚洲三级黄色毛片| 国产成人a∨麻豆精品| 日日啪夜夜撸| 蜜桃久久精品国产亚洲av| 亚洲国产精品专区欧美| 美女被艹到高潮喷水动态| 一级二级三级毛片免费看| 女人久久www免费人成看片| 亚洲综合精品二区| 中文字幕av成人在线电影| 如何舔出高潮| 亚洲一级一片aⅴ在线观看| 久久这里有精品视频免费| 乱人视频在线观看| 国产精品一区二区三区四区免费观看| 日韩欧美精品免费久久| 国内少妇人妻偷人精品xxx网站| 色视频www国产| 91精品一卡2卡3卡4卡| 少妇人妻一区二区三区视频| 免费看日本二区| 日本午夜av视频| 国产 一区 欧美 日韩| 激情五月婷婷亚洲| 免费大片黄手机在线观看| 久久6这里有精品| 国产精品一区二区在线观看99 | 欧美3d第一页| 草草在线视频免费看| 丝袜喷水一区| 精品一区在线观看国产| 你懂的网址亚洲精品在线观看| 久久久久久久大尺度免费视频| av一本久久久久| 国产午夜福利久久久久久| 亚洲人成网站在线播| 欧美性猛交╳xxx乱大交人| 国产免费一级a男人的天堂| 亚洲精品一区蜜桃| 日韩成人av中文字幕在线观看| 日韩人妻高清精品专区| 老女人水多毛片| 日韩av在线大香蕉| 国产三级在线视频| 可以在线观看毛片的网站| 亚洲va在线va天堂va国产| 黄片wwwwww| 日韩成人伦理影院|