• <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)行為研究——以黃河三角洲為例
    搡老乐熟女国产| 午夜福利网站1000一区二区三区| 3wmmmm亚洲av在线观看| 久久ye,这里只有精品| 亚洲国产最新在线播放| 新久久久久国产一级毛片| 伦理电影大哥的女人| 国产精品国产三级专区第一集| 中文精品一卡2卡3卡4更新| av国产精品久久久久影院| 韩国av在线不卡| videos熟女内射| 色94色欧美一区二区| 精华霜和精华液先用哪个| 丰满乱子伦码专区| 51国产日韩欧美| 男女无遮挡免费网站观看| 九九在线视频观看精品| 久久久久久伊人网av| 国产av国产精品国产| 超碰97精品在线观看| 亚洲精华国产精华液的使用体验| 日韩av不卡免费在线播放| 一级片'在线观看视频| 日韩欧美精品免费久久| 观看免费一级毛片| 男女免费视频国产| 夫妻午夜视频| 久久精品国产亚洲网站| 欧美激情国产日韩精品一区| 国产 一区精品| 9色porny在线观看| 男人狂女人下面高潮的视频| 五月开心婷婷网| 欧美精品一区二区免费开放| 黄色怎么调成土黄色| 国产乱来视频区| 少妇被粗大的猛进出69影院 | 中文字幕免费在线视频6| 久久av网站| 这个男人来自地球电影免费观看 | 麻豆成人午夜福利视频| 在线观看免费视频网站a站| 高清黄色对白视频在线免费看 | 18+在线观看网站| 精品人妻偷拍中文字幕| 美女xxoo啪啪120秒动态图| 99热6这里只有精品| 人妻 亚洲 视频| 欧美精品一区二区大全| 噜噜噜噜噜久久久久久91| 少妇人妻精品综合一区二区| 色5月婷婷丁香| 日本与韩国留学比较| 亚洲av电影在线观看一区二区三区| 最近2019中文字幕mv第一页| 午夜视频国产福利| 在线观看免费日韩欧美大片 | 青春草国产在线视频| 国产又色又爽无遮挡免| 精华霜和精华液先用哪个| av国产精品久久久久影院| 亚洲av不卡在线观看| 久久鲁丝午夜福利片| 日本色播在线视频| av国产久精品久网站免费入址| av线在线观看网站| 亚洲第一区二区三区不卡| 免费不卡的大黄色大毛片视频在线观看| 插阴视频在线观看视频| 日本欧美视频一区| 免费高清在线观看视频在线观看| 久久6这里有精品| 中文字幕亚洲精品专区| 女人久久www免费人成看片| 欧美日韩精品成人综合77777| 啦啦啦中文免费视频观看日本| 18禁动态无遮挡网站| 日韩,欧美,国产一区二区三区| 亚洲美女黄色视频免费看| 少妇裸体淫交视频免费看高清| 成人二区视频| 精品国产一区二区三区久久久樱花| 99热网站在线观看| 人妻少妇偷人精品九色| 精品少妇黑人巨大在线播放| 极品人妻少妇av视频| 国产在视频线精品| 少妇被粗大的猛进出69影院 | 久久av网站| 国产一区二区三区av在线| 视频中文字幕在线观看| 又爽又黄a免费视频| 一级毛片黄色毛片免费观看视频| 国产国拍精品亚洲av在线观看| 麻豆成人午夜福利视频| 内射极品少妇av片p| 最黄视频免费看| 久久国产亚洲av麻豆专区| 亚洲精品乱码久久久v下载方式| 成人亚洲精品一区在线观看| 亚洲av.av天堂| 国产精品熟女久久久久浪| 国产欧美日韩精品一区二区| 欧美日韩在线观看h| 最近中文字幕高清免费大全6| 三上悠亚av全集在线观看 | 国产老妇伦熟女老妇高清| 久热这里只有精品99| 亚洲一区二区三区欧美精品| 日日爽夜夜爽网站| 少妇精品久久久久久久| 性高湖久久久久久久久免费观看| 最新的欧美精品一区二区| 日韩中字成人| 欧美日本中文国产一区发布| 国产欧美日韩一区二区三区在线 | 黄色欧美视频在线观看| 免费大片黄手机在线观看| 亚洲精品亚洲一区二区| 久久久久精品久久久久真实原创| 亚洲欧美日韩卡通动漫| 国产乱来视频区| 成人美女网站在线观看视频| 色94色欧美一区二区| 日韩精品有码人妻一区| 午夜福利视频精品| 性色av一级| 国产精品不卡视频一区二区| 男人舔奶头视频| 精品久久国产蜜桃| av女优亚洲男人天堂| 又黄又爽又刺激的免费视频.| 久久久久久久久久人人人人人人| 国产一区二区三区av在线| 少妇裸体淫交视频免费看高清| 亚洲欧美一区二区三区国产| 国产高清三级在线| h日本视频在线播放| 国产黄频视频在线观看| 国产精品女同一区二区软件| 亚洲人与动物交配视频| 日本欧美视频一区| 亚洲经典国产精华液单| 欧美日韩视频精品一区| 亚洲精品国产色婷婷电影| 国产精品伦人一区二区| 大话2 男鬼变身卡| 亚洲,一卡二卡三卡| 精品亚洲成a人片在线观看| 纵有疾风起免费观看全集完整版| 亚洲欧美日韩卡通动漫| 欧美一级a爱片免费观看看| 久久热精品热| 毛片一级片免费看久久久久| 久久久精品免费免费高清| 亚洲真实伦在线观看| 久久久久人妻精品一区果冻| 啦啦啦中文免费视频观看日本| 五月伊人婷婷丁香| 久久影院123| 久久午夜福利片| 香蕉精品网在线| 99re6热这里在线精品视频| 色视频在线一区二区三区| 少妇被粗大猛烈的视频| 大片电影免费在线观看免费| 欧美日韩国产mv在线观看视频| 啦啦啦啦在线视频资源| 成人午夜精彩视频在线观看| 成年美女黄网站色视频大全免费 | 国产亚洲5aaaaa淫片| 夫妻午夜视频| 亚洲国产色片| 一个人免费看片子| 人体艺术视频欧美日本| 成人国产av品久久久| 精品少妇黑人巨大在线播放| a级毛色黄片| 精品亚洲乱码少妇综合久久| 欧美日本中文国产一区发布| 久久99蜜桃精品久久| 大码成人一级视频| 能在线免费看毛片的网站| 在线亚洲精品国产二区图片欧美 | 国产 精品1| 色哟哟·www| 人妻人人澡人人爽人人| 亚洲精品第二区| 色视频在线一区二区三区| 男人爽女人下面视频在线观看| 三上悠亚av全集在线观看 | 99久久中文字幕三级久久日本| 久久韩国三级中文字幕| 纯流量卡能插随身wifi吗| 国产欧美亚洲国产| 欧美日韩精品成人综合77777| 精品少妇久久久久久888优播| 国产精品久久久久久精品古装| 欧美+日韩+精品| 中文欧美无线码| 青春草亚洲视频在线观看| 97精品久久久久久久久久精品| 日韩精品免费视频一区二区三区 | 在现免费观看毛片| 久久国产精品男人的天堂亚洲 | 久久精品国产a三级三级三级| 三级经典国产精品| 久久久久久伊人网av| 成年人午夜在线观看视频| 两个人免费观看高清视频 | 国产精品99久久99久久久不卡 | 大话2 男鬼变身卡| 国产一区二区在线观看av| 久久国产亚洲av麻豆专区| 亚洲av在线观看美女高潮| 黄色日韩在线| 十分钟在线观看高清视频www | av黄色大香蕉| 边亲边吃奶的免费视频| 国产成人精品婷婷| 久久精品久久久久久噜噜老黄| 91在线精品国自产拍蜜月| 日韩一区二区三区影片| 97超碰精品成人国产| 亚洲av国产av综合av卡| 国产亚洲最大av| 国产精品99久久99久久久不卡 | 在线观看www视频免费| 精品少妇黑人巨大在线播放| 亚洲美女视频黄频| 亚洲精品日本国产第一区| 亚洲欧美成人精品一区二区| 在线观看av片永久免费下载| 激情五月婷婷亚洲| 狂野欧美白嫩少妇大欣赏| 内射极品少妇av片p| 亚洲国产精品专区欧美| 2021少妇久久久久久久久久久| a级一级毛片免费在线观看| 久久热精品热| 91久久精品国产一区二区三区| 人妻少妇偷人精品九色| 国产又色又爽无遮挡免| 久久99精品国语久久久| 又粗又硬又长又爽又黄的视频| 日韩 亚洲 欧美在线| 尾随美女入室| 亚洲天堂av无毛| 丰满迷人的少妇在线观看| 国产精品久久久久久久久免| 精品熟女少妇av免费看| 国产精品女同一区二区软件| 一个人免费看片子| 三级国产精品欧美在线观看| 少妇裸体淫交视频免费看高清| 91成人精品电影| 国产精品一区二区性色av| 桃花免费在线播放| 国产片特级美女逼逼视频| h视频一区二区三区| 久久久午夜欧美精品| av在线播放精品| 亚洲欧美日韩另类电影网站| 五月天丁香电影| 国产成人精品一,二区| 国产女主播在线喷水免费视频网站| √禁漫天堂资源中文www| 少妇人妻精品综合一区二区| 97超碰精品成人国产| 日本免费在线观看一区| 黄色配什么色好看| 国产精品国产三级专区第一集| 亚洲欧美成人综合另类久久久| 久久久久久久精品精品| 另类精品久久| 在线精品无人区一区二区三| 日韩伦理黄色片| 嫩草影院入口| 99热这里只有精品一区| 国产精品久久久久久久电影| 又粗又硬又长又爽又黄的视频| 免费黄色在线免费观看| 久久久久久久久久成人| 丰满人妻一区二区三区视频av| 亚洲精品日本国产第一区| 日日摸夜夜添夜夜添av毛片| 免费观看a级毛片全部| 国产免费一区二区三区四区乱码| 国产色婷婷99| 久久99精品国语久久久| 91久久精品国产一区二区成人| 久久精品国产鲁丝片午夜精品| 男女免费视频国产| av天堂久久9| 自线自在国产av| 三级国产精品欧美在线观看| 丰满迷人的少妇在线观看| 亚洲精品国产av蜜桃| freevideosex欧美| 亚洲精品第二区| 两个人的视频大全免费| 国产精品嫩草影院av在线观看| 国产白丝娇喘喷水9色精品| 一级毛片久久久久久久久女| 人人妻人人爽人人添夜夜欢视频 | 欧美 亚洲 国产 日韩一| 亚洲中文av在线| 免费大片18禁| 国产伦理片在线播放av一区| 国产日韩一区二区三区精品不卡 | 亚洲欧美日韩东京热| 噜噜噜噜噜久久久久久91| 成年女人在线观看亚洲视频| 有码 亚洲区| 一级爰片在线观看| 在线观看三级黄色| 久久国产精品男人的天堂亚洲 | 日本色播在线视频| 香蕉精品网在线| 街头女战士在线观看网站| 亚洲欧美日韩另类电影网站| 久久久午夜欧美精品| 国产精品99久久99久久久不卡 | 高清av免费在线| 国产色爽女视频免费观看| 亚洲伊人久久精品综合| 老熟女久久久| 亚洲av成人精品一二三区| 欧美日韩亚洲高清精品| 99热国产这里只有精品6| 亚洲av综合色区一区| 婷婷色麻豆天堂久久| 国产精品99久久99久久久不卡 | 2018国产大陆天天弄谢| 亚洲三级黄色毛片| 国产精品三级大全| 午夜免费男女啪啪视频观看| 日本猛色少妇xxxxx猛交久久| 精品久久久噜噜| 97超视频在线观看视频| 精品人妻一区二区三区麻豆| 18禁裸乳无遮挡动漫免费视频| av在线app专区| 亚洲av成人精品一二三区| 亚洲av不卡在线观看| 一区在线观看完整版| 在线观看www视频免费| 免费高清在线观看视频在线观看| 丰满饥渴人妻一区二区三| 男人舔奶头视频| 少妇 在线观看| 久久久亚洲精品成人影院| 丝袜脚勾引网站| 亚洲欧美精品自产自拍| 男人爽女人下面视频在线观看| 日本-黄色视频高清免费观看| 亚洲国产最新在线播放| 精品熟女少妇av免费看| 亚洲国产色片| 亚洲精品亚洲一区二区| 日本黄色片子视频| 亚洲熟女精品中文字幕| 嘟嘟电影网在线观看| 如日韩欧美国产精品一区二区三区 | 国产黄色免费在线视频| 国语对白做爰xxxⅹ性视频网站| 亚洲精品一二三| 丰满乱子伦码专区| 亚洲国产精品999| 成人18禁高潮啪啪吃奶动态图 | 欧美激情国产日韩精品一区| 亚洲色图综合在线观看| 亚洲久久久国产精品| 高清av免费在线| 中文字幕亚洲精品专区| 交换朋友夫妻互换小说| 能在线免费看毛片的网站| 久久久久久久久久成人| 国产亚洲欧美精品永久| 国产深夜福利视频在线观看| 男女啪啪激烈高潮av片| 亚洲国产最新在线播放| 麻豆成人av视频| 一级a做视频免费观看| 晚上一个人看的免费电影| 亚洲精品aⅴ在线观看| 一区二区三区免费毛片| 少妇被粗大的猛进出69影院 | 啦啦啦啦在线视频资源| 午夜福利在线观看免费完整高清在| 成人18禁高潮啪啪吃奶动态图 | 99热全是精品| 免费观看av网站的网址| 久久6这里有精品| 欧美高清成人免费视频www| 色吧在线观看| 欧美成人精品欧美一级黄| 边亲边吃奶的免费视频| 亚洲高清免费不卡视频| 欧美激情国产日韩精品一区| 免费久久久久久久精品成人欧美视频 | 国产乱人偷精品视频| 亚洲欧美精品专区久久| 在现免费观看毛片| 成人免费观看视频高清| 丰满迷人的少妇在线观看| 精品视频人人做人人爽| 桃花免费在线播放| 在线亚洲精品国产二区图片欧美 | 观看美女的网站| 国产在线视频一区二区| 一区二区三区免费毛片| 日本黄大片高清| 乱系列少妇在线播放| 国产永久视频网站| 国产美女午夜福利| 黄片无遮挡物在线观看| 高清毛片免费看| 成人国产av品久久久| 成人综合一区亚洲| 夫妻午夜视频| 亚洲精华国产精华液的使用体验| 久久综合国产亚洲精品| 欧美激情极品国产一区二区三区 | 亚洲精品国产成人久久av| 在线播放无遮挡| 精品少妇久久久久久888优播| 性高湖久久久久久久久免费观看| 高清欧美精品videossex| 日韩欧美 国产精品| 国产黄色视频一区二区在线观看| 2022亚洲国产成人精品| 亚洲精品456在线播放app| 亚洲欧美成人综合另类久久久| 精品人妻一区二区三区麻豆| 国产男女内射视频| 国产成人一区二区在线| 精品少妇内射三级| 高清黄色对白视频在线免费看 | 少妇人妻 视频| 中文字幕人妻熟人妻熟丝袜美| av又黄又爽大尺度在线免费看| 国产欧美日韩精品一区二区| 亚洲国产精品一区三区| 性色av一级| 精品人妻偷拍中文字幕| 热re99久久精品国产66热6| 五月伊人婷婷丁香| 在线观看人妻少妇| 成人特级av手机在线观看| 在线精品无人区一区二区三| 国模一区二区三区四区视频| 在线看a的网站| 久久精品国产亚洲网站| 亚洲欧洲国产日韩| 大又大粗又爽又黄少妇毛片口| 国产亚洲最大av| 成年人免费黄色播放视频 | 日韩成人伦理影院| 久久青草综合色| 成人亚洲欧美一区二区av| 国产亚洲精品久久久com| 男女啪啪激烈高潮av片| 女的被弄到高潮叫床怎么办| 菩萨蛮人人尽说江南好唐韦庄| 午夜影院在线不卡| 成人漫画全彩无遮挡| 日本黄色日本黄色录像| 高清欧美精品videossex| 亚洲av日韩在线播放| 51国产日韩欧美| 日本午夜av视频| 精品亚洲成国产av| 欧美国产精品一级二级三级 | 亚洲婷婷狠狠爱综合网| 纵有疾风起免费观看全集完整版| 成人二区视频| 一本色道久久久久久精品综合| 国产伦精品一区二区三区四那| 观看美女的网站| 最近中文字幕2019免费版| 老司机亚洲免费影院| 久久久久久久久久久久大奶| 天堂俺去俺来也www色官网| 又大又黄又爽视频免费| 国产淫语在线视频| 免费大片18禁| 高清午夜精品一区二区三区| 日韩一区二区三区影片| 日本vs欧美在线观看视频 | 欧美性感艳星| 国产男女超爽视频在线观看| 亚洲国产色片| 欧美日韩精品成人综合77777| 乱人伦中国视频| 亚洲欧洲日产国产| 黄色配什么色好看| 久久国内精品自在自线图片| 日日啪夜夜撸| 十八禁高潮呻吟视频 | 一个人看视频在线观看www免费| 欧美日韩国产mv在线观看视频| 成人毛片a级毛片在线播放| 国产69精品久久久久777片| 视频区图区小说| 亚洲三级黄色毛片| av国产精品久久久久影院| 日本-黄色视频高清免费观看| 建设人人有责人人尽责人人享有的| 亚洲欧美成人综合另类久久久| 一级毛片aaaaaa免费看小| 少妇裸体淫交视频免费看高清| 精品酒店卫生间| 热99国产精品久久久久久7| 久久99蜜桃精品久久| 精品亚洲成a人片在线观看| 久久ye,这里只有精品| 日韩熟女老妇一区二区性免费视频| 麻豆成人午夜福利视频| 午夜91福利影院| 亚洲激情五月婷婷啪啪| av线在线观看网站| 91久久精品国产一区二区三区| 三级经典国产精品| 国产精品熟女久久久久浪| 自拍偷自拍亚洲精品老妇| 高清在线视频一区二区三区| 国产高清国产精品国产三级| 国产男女超爽视频在线观看| 少妇熟女欧美另类| 少妇人妻久久综合中文| 国产精品国产av在线观看| 91久久精品国产一区二区成人| 亚洲精品aⅴ在线观看| 国产视频内射| 国产精品国产三级国产av玫瑰| 日日摸夜夜添夜夜添av毛片| 日韩免费高清中文字幕av| 亚洲精品,欧美精品| 91成人精品电影| 9色porny在线观看| 天堂中文最新版在线下载| 男人狂女人下面高潮的视频| 大码成人一级视频| 久久久久国产网址| 水蜜桃什么品种好| 国产一区二区在线观看av| 午夜免费男女啪啪视频观看| 成人无遮挡网站| 这个男人来自地球电影免费观看 | 波野结衣二区三区在线| 少妇 在线观看| 日本vs欧美在线观看视频 | 好男人视频免费观看在线| 男女啪啪激烈高潮av片| 亚洲国产精品国产精品| 成人漫画全彩无遮挡| 精品熟女少妇av免费看| 极品教师在线视频| 日韩av免费高清视频| 久久久久久久久久久免费av| 一级黄片播放器| 国产成人免费观看mmmm| 午夜福利,免费看| 日本欧美国产在线视频| 日本午夜av视频| 精品一品国产午夜福利视频| 男男h啪啪无遮挡| 久久久午夜欧美精品| 老女人水多毛片| 校园人妻丝袜中文字幕| 男的添女的下面高潮视频| 丝瓜视频免费看黄片| 久久久久人妻精品一区果冻| 欧美三级亚洲精品| 中文字幕制服av| 亚洲精品国产成人久久av| 精品国产露脸久久av麻豆| 国产精品久久久久久精品古装| 大片电影免费在线观看免费| h日本视频在线播放| 麻豆精品久久久久久蜜桃| 国产精品人妻久久久影院| 免费看光身美女| 国产精品99久久99久久久不卡 | 精品久久久精品久久久| 久久99蜜桃精品久久| 欧美97在线视频| 亚洲成人手机| 国产午夜精品久久久久久一区二区三区| 少妇精品久久久久久久| 亚洲欧美日韩卡通动漫| freevideosex欧美| 久久久久久久亚洲中文字幕| 啦啦啦视频在线资源免费观看| 婷婷色麻豆天堂久久| 午夜影院在线不卡| 在线观看免费高清a一片| 久久久午夜欧美精品| 日本91视频免费播放| 一区二区三区四区激情视频| 亚洲精品亚洲一区二区| 又黄又爽又刺激的免费视频.| 免费观看性生交大片5| 亚洲精品乱码久久久久久按摩| 一个人免费看片子| 啦啦啦视频在线资源免费观看| 2022亚洲国产成人精品| 国产色婷婷99| 久久午夜福利片| 免费观看a级毛片全部| 久久毛片免费看一区二区三区| 亚洲综合精品二区| 欧美bdsm另类|