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

    基于223Ra和224Ra的桑溝灣海底地下水排放通量

    2017-04-14 09:16:42王希龍杜金洲張經(jīng)
    海洋學(xué)報(bào) 2017年4期
    關(guān)鍵詞:潮汐活度表層

    王希龍,杜金洲*,張經(jīng)

    (1.華東師范大學(xué) 河口海岸學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海 200062)

    基于223Ra和224Ra的桑溝灣海底地下水排放通量

    王希龍1,杜金洲1*,張經(jīng)1

    (1.華東師范大學(xué) 河口海岸學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室,上海 200062)

    海底地下水排放(SGD)是陸地向海洋輸送水量和營(yíng)養(yǎng)物質(zhì)的重要通道之一,對(duì)沿海物質(zhì)通量及其生物地球化學(xué)循環(huán)有重要的影響,對(duì)生態(tài)環(huán)境起著不可忽視的作用。本文運(yùn)用天然放射性同位素223Ra和224Ra示蹤估算了我國(guó)北方典型養(yǎng)殖基地桑溝灣的海底地下水排放通量。結(jié)果表明,海底地下水樣尤其是間隙水中Ra活度[224Ra=(968±31) dpm/(100 L),223Ra=(31.4±4.9) dpm/(100 L),n=9]遠(yuǎn)高于表層海水[224Ra=(38.7±2.0) dpm/(100 L),223Ra=(1.70±0.50) dpm/(100 L),n=21]。假設(shè)穩(wěn)態(tài)條件下,考慮Ra的各源、匯項(xiàng),利用Ra平衡模型,估算出桑溝灣SGD排放通量為(0.23~1.03)×107m3/d。潮周期內(nèi)的觀測(cè)結(jié)果顯示,漲潮時(shí),水力梯度較小,SGD排放變?nèi)?,落潮時(shí),水力梯度較大,導(dǎo)致了相對(duì)較多的SGD排放。在一個(gè)潮周期間,基于223Ra和224Ra得到的SGD排放通量平均為0.39×107m3/d。潮汐動(dòng)力下的SGD排放平均占總SGD排放的61%,因此桑溝灣沿岸的地下水排放主要受潮汐動(dòng)力的影響,并對(duì)海水組成及海陸間物質(zhì)交換有顯著貢獻(xiàn)。

    鐳同位素;海底地下水排放;鐳庫(kù)模型;潮汐動(dòng)力;桑溝灣

    1 引言

    人類活動(dòng)影響對(duì)近海海域生態(tài)系統(tǒng)的主要表現(xiàn)之一就是陸海相互作用過程,而海底地下水排放(submarine groundwater discharge, SGD)是陸海相互作用過程中重要而往往容易被忽視的重要一環(huán)。SGD是全球水循環(huán)系統(tǒng)中重要的過程之一,由于地下水中的許多成分(如硝酸鹽)濃度與海水相比高出許多倍,SGD可能成為沿岸海域營(yíng)養(yǎng)物質(zhì)的重要載體和輸運(yùn)途徑。同時(shí),通過SGD攜帶的物質(zhì)可以改變近海營(yíng)養(yǎng)物質(zhì)的組成和結(jié)構(gòu),從而改變傳統(tǒng)的沿岸海域的生物地球化學(xué)循環(huán)模式[1—3]。

    自然界中存在4種天然放射性鐳同位素,226Ra(T1/2=1 600 a)、228Ra(T1/2=5.7 a)、224Ra(T1/2=3.7 d)和223Ra(T1/2=11.4 d)。目前這4種天然放射性鐳同位素已經(jīng)被廣泛應(yīng)用于SGD的研究中[4—10]。雖然我國(guó)SGD的研究起步較晚,但是相關(guān)工作在我國(guó)的三大河口及典型海灣都有所涉及[11—20]。

    桑溝灣位于山東省東部沿海(37°01′~37°09′N, 122°24′~122°35′E),是我國(guó)北方典型的養(yǎng)殖基地,北、西、南三面為陸地環(huán)抱,灣口向東,為半封閉型海灣(圖1)。桑溝灣水域廣闊,海灣面積163.2 km2,海底地勢(shì)平坦,由西向東逐漸傾斜,坡度較小,平均水深7~8 m,最大水深15~17 m[21]。桑溝灣沿岸無大河注入,入灣的較大河流有桑溝河、崖頭河、沽河、小落河等,年平均徑流總量為(1.7~2.3)×108m3,年輸沙量為17.1×104t。

    圖1 桑溝灣地理位置(左)及采樣站位圖(右),大面站、地下水和河水樣品采集于2015年5月,連續(xù)觀測(cè)樣品采集于2012年6月Fig.1 Location of Sanggou Bay (left) and the station locations in the Sanggou Bay (right): surface water, submarine groundwater, and river water were collected in May 2015; samples for time series station were collected in June 2012

    桑溝灣的潮汐類型為不正規(guī)半日潮,平均大、小潮差分別為1.47 m和0.57 m,最大潮差可達(dá)2.61 m。同時(shí)桑溝灣氣候溫和濕潤(rùn),變化幅度較小,月平均降水量為68.02 mm。灣內(nèi)水體溫度變化范圍較小,年平均水溫13℃左右,鹽度變化不大,年平均鹽度為31.76,1月份最低,6月份最高,鹽度變化受黃海沿岸流及降水的影響較為顯著[22]。

    作為我國(guó)典型立體示范養(yǎng)殖區(qū),歷史上,桑溝灣的水產(chǎn)資源十分豐富,在漁業(yè)生產(chǎn)旺季時(shí),曾有網(wǎng)產(chǎn)青魚數(shù)萬斤的記錄。自1957年在該灣推廣海帶筏式養(yǎng)殖生產(chǎn)技術(shù)以來,桑溝灣的水產(chǎn)活動(dòng)走上了以養(yǎng)為主、養(yǎng)捕結(jié)合的道路。但是由于對(duì)該灣水產(chǎn)資源的管理還存在不合理的地方,使得灣內(nèi)的水動(dòng)力情況及陸源物的輸入受到較大的影響。此前對(duì)于2012年桑溝灣SGD通量及其對(duì)營(yíng)養(yǎng)鹽收支平衡的評(píng)估已有報(bào)道[23],本文則是在此基礎(chǔ)上,利用短半衰期核素223Ra和224Ra對(duì)桑溝灣的SGD進(jìn)行再評(píng)估,并強(qiáng)調(diào)潮汐動(dòng)力對(duì)SGD的影響。

    2 采樣及實(shí)驗(yàn)方法

    2015年5月20日至27日期間,利用當(dāng)?shù)卮桓鶕?jù)養(yǎng)殖生物的不同分區(qū),以潛水泵采集桑溝灣表層水體樣品,樣品采集體積約60 L;同時(shí)進(jìn)行了桑溝灣沿岸地下井水(民用井)、沙灘間隙水和河流樣品的采集,地下井水樣品的采集體積為20 L,沙灘間隙水的采集體積為10 L,河流樣品的采集體積為40 L。所有的樣品均用于測(cè)定其中的224Ra和223Ra比活度,采樣站位如圖1所示。對(duì)于表層海水樣品和周圍的河流水樣,采用孔徑0.5 μm的濾芯進(jìn)行過濾除去水樣中的懸浮顆粒物,地下井水和沙灘間隙水樣品采用孔徑0.45 μm的醋酸纖維濾膜進(jìn)行過濾。同時(shí),2012年6月6日至7日在桑溝灣的南部定點(diǎn)(圖1中SGTS站)進(jìn)行了27 h的連續(xù)觀測(cè),每3 h采一次樣品,采樣及處理方法同上。

    3 結(jié)果

    3.1 桑溝灣223Ra和224Ra的活度分布

    224Ra和223Ra的活度分布如圖2所示。此處224Ra的活度值是指在測(cè)定的總224Ra活度中減掉其母體228Th的貢獻(xiàn)獲得的224Ra活度。表層水體中224Ra的活度范圍為18.9~81.9 dpm/(100 L),平均為(38.7±16.7) dpm/(100 L) (n=21);223Ra的活度范圍為0~2.89 dpm/(100 L),平均為(1.62±0.76) dpm/(100 L) (n=21)。在表層水體的活度分布等值線圖中可以看出224Ra和223Ra的活度均表現(xiàn)出明顯的近岸高、離岸低的趨勢(shì),但在研究區(qū)域的東北部均存在一個(gè)小幅的活度增長(zhǎng)趨勢(shì),推測(cè)這部分相對(duì)高的活度值可能來自桑溝灣以北的愛蓮灣的貢獻(xiàn)。愛蓮灣位于山東半島最東端,與桑溝灣相鄰(圖1),同樣是養(yǎng)殖型海灣,在氣候條件和海洋環(huán)境上與桑溝灣相同,但是桑溝灣為半封閉型海灣,愛蓮灣則相對(duì)開闊,與外海的交換也更充分[26],因此,此處的高值可能來自愛蓮灣的陸源輸入。此外,在圖1中站位SG-3、SG-9、SG-10、SG-11、SG-15處分別采集了底層水體樣品,其表、底層水體中224Ra和223Ra的活度列于表1中。從表1中可以看出,桑溝灣的淺水區(qū)(水深低于8 m)水體在垂向上基本混合均勻——水溫、鹽度和Ra的活度均沒有明顯的變化;在水深較深的SG-3站位處,測(cè)得的表層水體中224Ra比底層稍高,但是水溫和鹽度的變化并不大。

    海底地下水樣品(井水和沙灘間隙水)的采樣結(jié)果同樣顯示在圖2中。從圖中可以看出海底地下水樣品中的Ra活度變化范圍較大。井水樣品中的224Ra活度范圍是38.3~269 dpm/(100 L),平均為(104±73.8) dpm/(100 L) (n=8);223Ra活度范圍是0.03~5.36 dpm/(100 L),平均為(2.06±1.85) dpm/(100 L) (n=8)。岸邊沙灘間隙水中224Ra和223Ra的活度[224Ra:116~2 003,(968±653) dpm/(100 L);223Ra:1.78~50.9,(26.2±17.8) dpm/(100 L),n=9]則比井水高近10倍。

    圖2 桑溝灣表層海水224Ra(a)和223Ra(b)活度等值線圖及海底地下水(實(shí)心圓)、河流中(空心圓)的活度圖Fig.2 Contour plots of 223Ra and 224Ra in the surface water of Sanggou Bay, and dots plots for submarine groundwater (filled circles) and river water (open circle) around Sanggou Bay

    表1 桑溝灣部分站位的表、底層水的溫、鹽和鐳活度

    Tab.1 Temperature, salinity and Ra activities in several stations of Sanggou Bay

    站位水深/m層次水溫/℃鹽度224Ra誤差223Ra誤差dpm·(100L)-1SG?318 2表13 530 522 41 50 980 28底11 431 516 00 90 840 26SG?99 7表1531 943 42 21 010 34底14 231 940 41 71 080 38SG?109 1表16 231 440 41 71 920 60底15 931 332 82 30 760 33SG?118 5表--42 62 11 550 54底--34 62 63 430 99SG?158 0表17 131 236 11 62 080 63底16 131 436 61 81 310 46

    圖3 2015年5月表層水體(a)、地下井水(b)、間隙水(c)和河流(d)中224Ra vs. 223Ra及其相關(guān)性數(shù)據(jù)Fig.3 Plots of 224Ra vs. 223Ra for all samples, surface water (a), groundwater (b), pore water (c) and river water (d). The results of correlation analysis were also shown in the plots

    3.2 同位素比值

    近岸水體中Ra同位素的主要來源有河流沉積物的解吸和含鹽地下水的排放[5,27]。從沉積物上釋放之后,224Ra和223Ra不斷通過吸附在沉積物上的Th同位素的衰變產(chǎn)生,但是這種輸入在一個(gè)較長(zhǎng)的時(shí)間必然達(dá)到穩(wěn)定。對(duì)于224Ra和223Ra來講,即使具有相同的來源也會(huì)因?yàn)檩斎脒^程中的衰變而產(chǎn)生不同的活度比。圖3顯示了224Ra與223Ra在不同水體中的比值。從圖3中可以看出表層水體、地下井水、間隙水和河水中224Ra與223Ra之間的線性擬合直線斜率分別為20.9 (r=0.94,P<0.01,n=21)、39.7 (r=0.88,P<0.01,n=8)、31.8 (r=0.93,P<0.01,n=9)和16.8 (r=0.91,P<0.01,n=5),并表現(xiàn)出顯著的線性相關(guān)性。與地下水(井水和沙灘間隙水)相比,桑溝灣表層水體中的活度比值較低,這可能是由于地下水中224Ra不斷地由母體產(chǎn)生,而當(dāng)其進(jìn)入海水后,母體的貢獻(xiàn)減小,224Ra不斷衰變?cè)斐傻?。而與表層水體相比,河水中的224Ra與223Ra的活度及其比值均較低,因此表層水體中的224Ra/223Ra活度比(AR).很可能是由于Ra同位素在通過地下水向海輸入過程中的衰變?cè)斐傻摹?/p>

    4 討論

    4.1 桑溝灣水體表觀年齡

    水體表觀年齡可以用來描述水體的遷移擴(kuò)散過程,是了解水體及其物質(zhì)在海洋中輸送的水動(dòng)力過程的基礎(chǔ)[28-29]。要估算桑溝灣的海底地下水排放通量,也必然要先得到桑溝灣水體的停留時(shí)間。天然存在的Ra同位素,基于一個(gè)短半衰期和一個(gè)長(zhǎng)半衰期核素的活度比可用于確定近岸水體的表觀年齡[30]。近岸水體中Ra的添加項(xiàng)主要來自河流、沉積物以及地下水的輸入,Ra的損失項(xiàng)則主要由與低Ra活度水體的混合及其自身放射性衰變?cè)斐伞<僭O(shè)系統(tǒng)處于穩(wěn)態(tài),那么基于Moore等[30]的方法,針對(duì)224Ra和223Ra可以寫出如下平衡方程:

    (1)

    (2)

    式中,F(xiàn)224Ra和F223Ra分別代表輸入系統(tǒng)的224Ra和223Ra的總通量;I224Ra和I223Ra分別代表系統(tǒng)中224Ra和223Ra的總量;λ224和λ223分別代表224Ra和223Ra的衰變常數(shù),分別為0.189 d-1和0.060 8 d-1;τ代表系統(tǒng)的水體表觀年齡。用公式(1)除以公式(2)可以得到:

    (3)

    將公式(3)整理并求解水體表觀年齡τ,可以得到:

    (4)

    (5)

    式中,V是桑溝灣的平均水體體積,12.2×108m3;Q=P/T,P是潮周期內(nèi)的納潮量,為3.74×108m3,T是1個(gè)潮周期的時(shí)間,為0.517 d;In在這里代表河流對(duì)桑溝灣的凈輸入量,為5.48×105m3/d;b則代表回流因子,即外海水在研究區(qū)域內(nèi)的貢獻(xiàn)比例[30],基于水體、鹽度和226Ra的三端元混合模型,可以得到b=0.743[23]。因此,可以估算出桑溝灣的水體沖刷時(shí)間為6.56 d,與本研究結(jié)果在誤差范圍內(nèi)基本一致。因此,本研究中采用t=6.1 d作為桑溝灣水體滯留時(shí)間。

    4.2 過剩Ra庫(kù)的估算

    桑溝灣表層水體中短半衰期核素224Ra和223Ra的主要源、匯項(xiàng)包括河流輸入、與外海水的混合、自身的衰變、沉積物的貢獻(xiàn)以及SGD的輸入等。我們先前的研究表明,黃海對(duì)桑溝灣的水體貢獻(xiàn)比例平均為0.743,并且黃海水體中224Ra的平均活度為4.99 dpm/(100 L),223Ra則為0.20 dpm/(100 L)[23],由此可以得出桑溝灣表層水體中每個(gè)站位的過剩223Ra和224Ra。由于桑溝灣水體水深較淺,平均為7.5 m,因此我們采用表層水體中的Ra活度代表整個(gè)水柱中的Ra活度,結(jié)合桑溝灣的水域面積和平均水深,可以得到過剩223Ra和224Ra庫(kù)分別為1.80×1010和4.28×1011dpm。結(jié)合上述桑溝灣的水體滯留時(shí)間,可以獲得其過剩223Ra和224Ra通量分別為2.95×109和7.02×1010dpm/d。類似地,對(duì)于河流輸入的Ra通量,采用先前獲得的河流的貢獻(xiàn)比(0.045)[23]乘以河流端元樣品的Ra活度[223Ra=0.64 dpm/(100 L),224Ra=8.95 dpm/(100 L)]和河流的徑流量[(1.7~2.3)×108m3/a[21],本文采用中間值2.0×108m3/a],即河流輸入的223Ra和224Ra通量分別為1.58×105和2.21×106dpm/d。由此可見,相對(duì)于前述桑溝灣過剩通量,河流的貢獻(xiàn)可以忽略。

    對(duì)于桑溝灣水柱中的Ra來講,底層沉積物孔隙中Ra的向上擴(kuò)散是其另一個(gè)可能的來源。為了估算桑溝灣底層沉積物的貢獻(xiàn),我們根據(jù)公式(6)對(duì)Ra的最大擴(kuò)散通量進(jìn)行了估算[33]:

    (6)

    式中,F(xiàn)dif表示底層沉積物的釋放通量,單位dpm/d;Ds表示鐳在沉積物中的擴(kuò)散系數(shù),一般根據(jù)鐳在水體中的擴(kuò)散系數(shù)與沉積物的孔隙度求得[34—35],此處取4.0×10-6cm2/s;Kd表示分布系數(shù)(80),A為研究區(qū)域的水域面積(1.63×108m2);對(duì)于224Ra考慮其母體228Th的貢獻(xiàn),但由于223Ra的母體227Th的半衰期較短(18.7 d),227Th的母體227Ac雖然半衰期較長(zhǎng)(21.8 a),但在海水表層沉積物中的活度未曾有報(bào)道,因此對(duì)于表層沉積物中228Th的活度采用黃海的表層沉積物報(bào)道值(4.418 Bq/g)[36],再根據(jù)沙灘間隙水中的224Ra/223Ra活度比(42.8)對(duì)沉積物貢獻(xiàn)的223Ra進(jìn)行估算,分別得到沉積物貢獻(xiàn)的224Ra通量為2.1×108dpm/d,223Ra通量為4.8×106dpm/d。與上述桑溝灣的過剩通量結(jié)果相比,沉積物貢獻(xiàn)的223Ra和224Ra小2~3個(gè)數(shù)量級(jí),此值更印證了Wang等[23]的報(bào)道中沉積物的貢獻(xiàn)可以忽略這一結(jié)論。值得一提的是盡管溫度等環(huán)境因素可能會(huì)影響鐳同位素在沉積物和水體中的擴(kuò)散,但變化范圍不可能達(dá)到2~3個(gè)數(shù)量級(jí),因此,此處并未考慮溫度等參數(shù)對(duì)擴(kuò)散系數(shù)的影響。由于223Ra和224Ra的半衰期較短,其自身衰變是不可忽略項(xiàng),衰變損失項(xiàng)可由實(shí)測(cè)Ra庫(kù)值和衰變系數(shù)及上述水體滯留時(shí)間算得,分別為1.01×109dpm/d和5.31×1010dpm/d。因此以上述過剩Ra通量減去河流、沉積物的貢獻(xiàn)和衰變量,即可得到由SGD貢獻(xiàn)的223Ra和224Ra通量,分別為1.94×109dpm/d和1.71×1010dpm/d。

    4.3 SGD的估算

    由于近岸水體中的Ra核素不斷被外海的低Ra活度海水稀釋,所以觀測(cè)到近岸高值必然來自陸源的貢獻(xiàn)[37]。為了計(jì)算SGD,先要知道由SGD貢獻(xiàn)的過剩223Ra或224Ra和地下水端元中的223Ra或224Ra活度。通過上述計(jì)算,SGD輸送的223Ra和224Ra的通量已經(jīng)確定,因此,僅需要確定地下水端元的223Ra和224Ra的活度。既然海底地下水貢獻(xiàn)了桑溝灣的過剩223Ra和224Ra,那么地下水端元中的Ra活度必然高于桑溝灣水體中的Ra活度,同時(shí)其活度比(AR)必然大于過剩224Ra與223Ra的比值,據(jù)此篩選出地下水端元的224Ra的活度范圍為127~2 003 dpm/(100 L),平均值為737 dpm/(100 L),223Ra的活度范圍為4.06~50.9 dpm/(100 L),平均值為18.9 dpm/(100 L)。將上述算得的SGD輸送的223Ra和224Ra通量分別除以地下水端元中的223Ra和224Ra的平均活度,可以得到桑溝灣地區(qū)的SGD入海通量為(0.23~1.03)×107m3/d,此值與文獻(xiàn)報(bào)道的桑溝灣2012年6月的SGD通量[(2.59~3.07)×107m3/d][23]相比略小,可能與示蹤劑、季節(jié)、降雨以及養(yǎng)殖等存在一定的關(guān)系。在先前的報(bào)道中采用長(zhǎng)半衰期的226Ra和228Ra作為示蹤劑估算桑溝灣的SGD排放量,桑溝灣的水體年齡尺度上為數(shù)天,對(duì)于226Ra和228Ra來講,其衰變可以忽略;本研究中采用半衰較短的223Ra和224Ra作為示蹤劑,盡管海底地下水的采樣都選在最接近海灣的位置甚至在退潮時(shí)的沙灘,但確實(shí)存在一定的衰變影響,因此本研究得到的SGD結(jié)果在一定程度上低估了桑溝灣的SGD排放??偨Y(jié)榮成地區(qū)的降雨量信息可以發(fā)現(xiàn),不管是在年代較遠(yuǎn)的1959—1980年[21]還是在1981—2010年(中國(guó)氣象數(shù)據(jù)網(wǎng)),一年當(dāng)中5月的降水量始終比6月少,有時(shí)甚至僅有6月的1/2,進(jìn)而通過降水造成的5月地下徑流比6月也少,從而SGD排放中淡水部分相對(duì)減少。此外,受到桑溝灣海水養(yǎng)殖活動(dòng)的影響,5月海帶開始收獲,6月(夏季)海帶基本收獲完畢,水體流動(dòng)性加強(qiáng)[38],風(fēng)浪及潮汐作用下SGD的排放也可能會(huì)相應(yīng)變快。

    4.4 潮汐動(dòng)力下的桑溝灣SGD的估算

    4.4.1 潮周期內(nèi)Ra同位素的釋放通量

    在近岸海灣、河口等地區(qū),潮汐作用通常會(huì)成為影響近岸水體與外海交換的一個(gè)重要因素。為了評(píng)價(jià)桑溝灣地區(qū)潮汐作用對(duì)SGD輸送的影響,本文于2012年6月6日至7日在桑溝灣南部(37°3′18″N, 122°31′42″E)的魚排上進(jìn)行了27 h潮周期內(nèi)連續(xù)的定點(diǎn)觀測(cè),每3 h采一次樣品,樣品的采集及處理方法與上述相同。觀測(cè)期間鹽度從低潮時(shí)的30.5變化到高潮時(shí)的31.1。同樣地,Ra同位素的活度也表現(xiàn)潮周期內(nèi)的浮動(dòng)變化,與鹽度變化相反,表現(xiàn)為在高潮時(shí)活度較低而在低潮時(shí)活度較高(圖4)。這與高潮時(shí)帶有較低Ra活度的外海水進(jìn)入桑溝灣,低潮退出時(shí)帶入沿岸的高活度Ra有關(guān)。

    圖4 桑溝灣水體鹽度(a)、224Ra (b)和223Ra (c)活度在潮周期內(nèi)的變化情況Fig.4 The tidal level condition and the measured salinity (a), 224Ra (b), 223Ra (c) activities of Sanggou Bay during time series sampling

    在潮汐作用下,Ra同位素在海灣內(nèi)的交換過程可以表示為:

    (7)

    (8)

    其中水體的交換量可以表示為:

    Wexch=A×H×n,

    (9)

    4.4.2 潮汐動(dòng)力下的SGD排放通量

    考慮到桑溝灣Ra同位素的源匯,以及潮周期定點(diǎn)采樣的弊端,為了更合理的評(píng)估桑溝灣潮汐動(dòng)力下的SGD排放通量,我們采用如下公式進(jìn)行計(jì)算[15]:

    (10)

    式中,Ratotal是潮周期內(nèi)的測(cè)得的每一個(gè)樣品的Ra活度,Rabkgd是外海Ra活度的背景值,為了規(guī)避潮汐外的其他影響因素,以潮周期內(nèi)Ra活度的最低值作為此背景值,同樣地為了避免偶然性的影響將第一個(gè)最高值去掉,h是水深,A是海灣的面積,τ是水體滯留時(shí)間,Ragw是地下水端元中的Ra活度。由于Ratotal是基于連續(xù)站觀測(cè)的每一個(gè)時(shí)間點(diǎn)的Ra的活度,因此就能夠產(chǎn)生一系列的SGD輸入通量(圖5)。對(duì)于224Ra和223Ra,可以分別得到潮汐動(dòng)力下SGD的排放通量分別為(0.07~0.18)×107和(0.14~1.23)×107m3/d,平均值分別為0.13×107m3/d和0.64×107m3/d。此結(jié)果與2012年6月基于表層水樣的計(jì)算結(jié)果相比略小,與上述2015年5月基于Ra庫(kù)的結(jié)果相差不大,這反映了桑溝灣在雨量較少的春季,潮汐動(dòng)力對(duì)于SGD的排放起著至關(guān)重要的作用。需要指出的是與長(zhǎng)半衰期核素相比,短半衰期核素對(duì)環(huán)境變化比較敏感,且測(cè)量誤差較大,會(huì)對(duì)結(jié)果造成一定的影響。因此,上述結(jié)果存在較大范圍的浮動(dòng)。

    圖5 基于潮周期內(nèi)224Ra (a)和223Ra (b)得到的SGD通量的周期變化Fig.5 SGD variation during the tidal cycle based on 224Ra (a) and 223Ra (b)

    由于不同的地理環(huán)境和沿岸沉積物的各項(xiàng)異性,SGD通量在空間上變化很大,這一點(diǎn)在僅就我國(guó)沿岸海域及小海灣的SGD研究就可以看出(表2)。從北到南,Ma等[41]和Wang等[42]對(duì)萊州灣的SGD和SFGD(地下淡水排放)進(jìn)行了研究發(fā)現(xiàn)萊州灣洪季的SFGD和SGD的通量分別相當(dāng)于黃河年徑流量的0.57~0.88以及7.35~8.57倍。郭占榮等[12-14]分別運(yùn)用Ra和Rn同位素先后評(píng)價(jià)了福建省隆教灣、九龍江河口和五緣灣的SGD通量,其中九龍江地區(qū)的地下水輸入量是河流徑流量的4倍多。Wang等[43]對(duì)海南三亞灣珊瑚礁系統(tǒng)的SGD進(jìn)行了研究,發(fā)現(xiàn)SGD是海南三亞灣近岸水體酸化的重要誘因。以上研究表明對(duì)于近岸海域尤其是封閉半封閉的小海灣的SGD的研究是十分必要的,SGD對(duì)于類似的小海灣的生態(tài)環(huán)境具有不可忽略的影響。

    4.5 SGD估算的不確定性分析

    5 結(jié)論

    通過我國(guó)典型立體養(yǎng)殖區(qū)——桑溝灣的海水、地下水和河流樣品中短半衰期核素224Ra和223Ra的分布規(guī)律對(duì)春季桑溝灣的水體滯留時(shí)間和SGD排放通量進(jìn)行了研究,發(fā)現(xiàn)桑溝灣的水體滯留時(shí)間平均為6.1 d;運(yùn)用Ra庫(kù)模型示蹤估算得到2015年5月桑溝灣的SGD排放通量為(0.23~1.03)×107m3/d,盡管這一SGD通量存在著40%~62%的不確定度,但是這一通量預(yù)示著通過SGD排放的重要物質(zhì)通量可能會(huì)對(duì)桑溝灣的生源要素的收支平衡具有重要影響。此外,潮周期內(nèi)的連續(xù)觀測(cè)表明,桑溝灣Ra同位素的交換量分別是2.57×1010(224Ra)和0.40×1010dpm/d (223Ra),相應(yīng)SGD排放通量的平均值分別為0.13×107和0.64×107m3/d。對(duì)比發(fā)現(xiàn),盡管多種驅(qū)動(dòng)力促使著SGD的排放,但是在桑溝灣這樣一個(gè)半封閉的養(yǎng)殖型海灣中,潮汐動(dòng)力發(fā)揮著重要的作用。

    [1] Johannes R E. The ecological significance of the submarine discharge of groundwater[J]. Marine Ecology Progress Series, 1980, 3(4): 365-373.

    [2] Maher D T, Santos I R, Golsby-Smith L, et al. Groundwater-derived dissolved inorganic and organic carbon exports from a mangrove tidal creek: the missing mangrove carbon sink[J]. Limnology and Oceanography, 2013, 58(2): 475-488.

    [3] Kwon E Y, Kim G, Primeau F, et al. Global estimate of submarine groundwater discharge based on an observationally constrained radium isotope model[J]. Geophysical Research Letters, 2014, 41(23): 8438-8444.

    [4] Webster I T, Hancock G J, Murray A S. Use of radium isotopes to examine pore-water exchange in an estuary[J]. Limnology and Oceanography, 1994, 39(8): 1917-1927.

    [5] Rama, Moore W S. Using the radium quartet for evaluating groundwater input and water exchange in salt marshes[J]. Geochimica et Cosmochimica Acta, 1996, 60(23): 4645-4652.

    [6] Hancock G J, Webster I T, Ford P W, et al. Using Ra isotopes to examine transport processes controlling benthic fluxes into a shallow estuarine lagoon[J]. Geochimica et Cosmochimica Acta, 2000, 64(21): 3685-3699.

    [7] Krest J M, Harvey J W. Using natural distributions of short-lived radium isotopes to quantify groundwater discharge and recharge[J]. Limnology and Oceanography, 2003, 48(1): 290-298.

    [8] Beck A J, Rapaglia J P, Cochran J K, et al. Radium mass-balance in Jamaica Bay, NY: evidence for a substantial flux of submarine groundwater[J]. Marine Chemistry, 2007, 106(3/4): 419-441.

    [9] Colbert S L, Hammond D E. Shoreline and seafloor fluxes of water and short-lived Ra isotopes to surface water of San Pedro Bay, CA[J]. Marine Chemistry, 2008, 108(1/2): 1-17.

    [10] Moore W S, Beck M, Riedel T, et al. Radium-based pore water fluxes of silica, alkalinity, manganese, DOC, and uranium: a decade of studies in the German Wadden Sea[J]. Geochimica et Cosmochimica Acta, 2011, 75(21): 6535-6555.

    [11] 黃磊. 九龍江河口區(qū)的地下水輸入研究[D]. 廈門: 廈門大學(xué), 2009.

    Huang Lei. Research on groundwater discharge into Jiulongjiang Estuary[D]. Xiamen: Xiamen University, 2009.

    [12] 郭占榮, 黃磊, 劉花臺(tái), 等. 鐳同位素示蹤隆教灣的海底地下水排泄[J]. 地球?qū)W報(bào), 2008, 29(5): 647-652.

    Guo Zhanrong, Huang Lei, Liu Huatai, et al. The estimation of submarine inputs of groundwater to a coastal bay using radium isotopes[J]. Acta Geoscientica Sinica, 2008, 29(5): 647-652.

    [13] 郭占榮, 黃磊, 袁曉婕, 等. 用鐳同位素評(píng)價(jià)九龍江河口區(qū)的地下水輸入[J]. 水科學(xué)進(jìn)展, 2011, 22(1): 118-125.

    Guo Zhanrong, Huang Lei, Yuan Xiaojie, et al. Estimating submarine groundwater discharge to the Jiulong River estuary using Ra isotopes[J]. Advances in Water Science, 2011, 22(1): 118-125.

    [14] 郭占榮, 李開培, 袁曉婕, 等. 用氡-222評(píng)價(jià)五緣灣的地下水輸入[J]. 水科學(xué)進(jìn)展, 2012, 23(2): 263-270.

    Guo Zhanrong, Li Kaipei, Yuan Xiaojie, et al. Assessment of submarine groundwater discharge into the Wuyuan Bay via continuous Radon-222 measurements[J]. Advances in Water Science, 2012, 23(2): 263-270.

    [15] Peterson R N, Burnett W C, Taniguchi M, et al. Radon and radium isotope assessment of submarine groundwater discharge in the Yellow River delta, China[J]. Journal of Geophysical Research, 2008, 113(C9): C09021.

    [16] Tse K C, Jiao J J. Estimation of submarine groundwater discharge in Plover Cove, Tolo Harbour, Hong Kong by222Rn[J]. Marine Chemistry, 2008, 111(3/4): 160-170.

    [17] Su Ni, Du Jinzhou, Moore W S, et al. An examination of groundwater discharge and the associated nutrient fluxes into the estuaries of eastern Hainan Island, China using226Ra[J]. Science of the Total Environment, 2011, 409(19): 3909-3918.

    [18] Gu Hequan, Moore W S, Zhang Lei, et al. Using radium isotopes to estimate the residence time and the contribution of submarine groundwater discharge (SGD) in the Changjiang effluent plume, East China Sea[J]. Continental Shelf Research, 2012, 35: 95-107.

    [19] Liu Q, Dai M, Chen W, et al. How significant is submarine groundwater discharge and its associated dissolved inorganic carbon in a river-dominated shelf system[J]. Biogeosciences, 2012, 9(5): 1777-1795.

    [20] Ji Tao, Du Jinzhou, Moore W S, et al. Nutrient inputs to a Lagoon through submarine groundwater discharge: the case of Laoye Lagoon, Hainan, China[J]. Journal of Marine Systems, 2013, 111-112: 253-262.

    [21] 夏東興, 王文海, 劉傳信, 等. 中國(guó)海灣志·第三分冊(cè)·山東半島北部和東部海灣·第九章[M]. 北京: 海洋出版社, 1991.

    Xia Dongxing, Wang Wenhai, Liu Chuanxin, et al. The Journal of Chinese Gulf-The Third Section: Bays of the Northern and Eastern Shandong Peninsula-Chapter 9[M]. Beijing: China Ocean Press, 1991.

    [22] 毛興華, 張為先, 張建中. 桑溝灣增養(yǎng)殖環(huán)境綜合調(diào)查研究[M]. 青島: 青島出版社, 1988.

    Mao Xinghua, Zhang Weixian, Zhang Jianzhong. The Comprehensive Investigation of Aquaculture in Sanggou Bay[M]. Qingdao: Qingdao Press, 1988.

    [23] Wang Xilong, Du Jinzhou, Ji Tao, et al. An estimation of nutrient fluxes via submarine groundwater discharge into the Sanggou Bay-a typical multi-species culture ecosystem in China[J]. Marine Chemistry, 2014, 167: 113-122.

    [24] Moore W S. Sampling228Ra in the deep ocean[J]. Deep-Sea Research and Oceanographic Abstracts, 1976, 23(7): 647-651.

    [25] 谷河泉, 杜金洲, 吳梅桂, 等. 鐳延遲符合計(jì)數(shù)器(RaDeCC)測(cè)量海水中的224Ra和223Ra[J]. 海洋環(huán)境科學(xué), 2015, 34(4): 570-577.

    Gu Hequan, Du Jinzhou, Wu Meigui, et al. Precise determination of224Ra and223Ra in seawaters By RaDeCC[J]. Marine Environmental Science, 2015, 34(4): 570-577.

    [26] 李磊, 任景玲, 劉素美, 等. 桑溝灣溶解態(tài)無機(jī)砷的分布、季節(jié)變化及影響因素[J]. 環(huán)境科學(xué), 2014, 35(7): 2705-2713.

    Li Lei, Ren Jingling, Liu Sumei, et al. Distribution, seasonal variation and influence factors of dissolved inorganic arsenic in the Sanggou Bay[J]. Environmental Science, 2014, 35(7): 2705-2713.

    [27] Moore W S. Radium-228 in the south Atlantic bight[J]. Journal of Geophysical Research, 1987, 92(C5): 5177-5190.

    [28] Dronkers J, Zimmerman J T F. Some principles of mixing in tidal lagoons //Lasserre P, et al. (Ed.) Coastal lagoons: Proceedings of the international symposium on coastal lagoons, Bordeaux, France, 8-14 September 1981[J]. Oceanologica Acta, 1982, 4(S): 107-117.

    [29] Monsen N E, Cloern J E, Lucas L V, et al. A comment on the use of flushing time, residence time, and age as transport time scales[J]. Limnology and Oceanography, 2002, 47(5): 1545-1553.

    [30] Moore W S, Blanton J O, Joye S B. Estimates of flushing times, submarine groundwater discharge, and nutrient fluxes to Okatee Estuary, South Carolina[J]. Journal of Geophysical Research, 2006, 111(C9): C09006, doi: 10.1029/2005JC003041.

    [31] Sanford L P, Boicourt W C, Rives S R. Model for estimating tidal flushing of small embayments[J]. Journal of Waterway, Port, Coastal, and Ocean Engineering, 1992, 118(6): 635-654.

    [32] Wang Guizhi, Wang Zhangyong, Zhai Weidong, et al. Net subterranean estuarine export fluxes of dissolved inorganic C, N, P, Si, and total alkalinity into the Jiulong River estuary, China[J]. Geochimica et Cosmochimica Acta, 2015, 149: 103-114.

    [33] Krest J M, Moore W S, Rama.226Ra and228Ra in the mixing zones of the Mississippi and Atchafalaya Rivers: indicators of groundwater input[J]. Marine Chemistry, 1999, 64(3): 129-152.

    [34] 蔡立勝, 方建光, 董雙林. 桑溝灣養(yǎng)殖海區(qū)沉積物-海水界面氮、磷營(yíng)養(yǎng)鹽的通量[J]. 海洋水產(chǎn)研究, 2004, 25(4): 57-64.

    Cai Lisheng, Fang Jianguang, Dong Shuanglin. Preliminary studies on nitrogen and phosphorus fluxes between seawater and sediment in Sungo Bay[J]. Marine Fisheries Research, 2004, 25(4): 57-64.

    [35] Li Yuanhui, Gregory S. Diffusion of ions in sea water and in deep-sea sediments[J]. Geochimica et Cosmochimica Acta, 1974, 38(5): 703-714.

    [36] 洪華生. 中國(guó)區(qū)域海洋學(xué)——化學(xué)海洋學(xué)[M]. 北京: 海洋出版社, 2012.

    Hong Huasheng. Regional Oceanography of China Seas-Chemical Oceanography[M]. Beijing: China Ocean Press, 2012.

    [37] Moore W S. Seasonal distribution and flux of radium isotopes on the southeastern U.S. continental shelf[J]. Journal of Geophysical Research, 2007, 112(C10): C10013.

    [38] 陳聚法, 趙俊, 孫耀, 等. 桑溝灣貝類養(yǎng)殖水域沉積物再懸浮的動(dòng)力機(jī)制及其對(duì)水體中營(yíng)養(yǎng)鹽的影響[J]. 海洋水產(chǎn)研究, 2007, 28(3): 105-111.

    Chen Jufa, Zhao Jun, Sun Yao, et al. Dynamic mechanism of sediment resuspension with its effects on content of nutrients in water in the shellfish culture area of Sanggou Bay[J]. Marine Fisheries Research, 2007, 28(3): 105-111.

    [39] Garcia-Orellana J, Cochran J K, Bokuniewicz H, et al. Time-series sampling of223Ra and224Ra at the inlet to Great South Bay (New York): a strategy for characterizing the dominant terms in the Ra budget of the bay[J]. Journal of Environmental Radioactivity, 2010, 101(7): 582-588.

    [40] 杜鵬. 夏季桑溝灣三維水質(zhì)數(shù)值模擬[D]. 青島: 中國(guó)海洋大學(xué), 2009.

    Du Peng. The three-dimensional water quality numerical simulation in Sanggou Bay in summer[D]. Qingdao: Ocean University of China, 2009.

    [41] Ma Qian, Li Hailong, Wang Xuejing, et al. Estimation of seawater-groundwater exchange rate: case study in a tidal flat with a large-scale seepage face (Laizhou Bay, China)[J]. Hydrogeology Journal, 2015, 23(2): 265-275.

    [42] Wang Xuejing, Li Hailong, Jiao J J, et al. Submarine fresh groundwater discharge into Laizhou Bay comparable to the Yellow River flux[J]. Scientific Reports, 2015, 5: 8814.

    [43] Wang Guizhi, Jing Wenping, Wang Shuling, et al. Coastal acidification induced by tidal-driven submarine groundwater discharge in a coastal coral reef system[J]. Environmental Science & Technology, 2014, 48(22): 13069-13075.

    [44] Kim G, Ryu J W, Yang H S, et al. Submarine groundwater discharge (SGD) into the Yellow Sea revealed by228Ra and226Ra isotopes: implications for global silicate fluxes[J]. Earth and Planetary Science Letters, 2005, 237(1/2): 156-166.

    [45] Wu Zijun, Zhou Huaiyang, Zhang Shuai, et al. Using222Rn to estimate submarine groundwater discharge (SGD) and the associated nutrient fluxes into Xiangshan Bay, East China Sea[J]. Marine Pollution Bulletin, 2013, 73(1): 183-191.

    [46] Xu Bochao, Burnett W, Dimova N, et al. Hydrodynamics in the Yellow River Estuary via radium isotopes: ecological perspectives[J]. Continental Shelf Research, 2013, 66: 19-28.

    [47] Wen Tingyu, Du Jinzhou, Ji Tao, et al. Use of222Rn to trace submarine groundwater discharge in a tidal period along the coast of Xiangshan, Zhejiang, China[J]. Journal of Radioanalytical and Nuclear Chemistry, 2014, 299(1): 53-60.

    [48] Luo Xin, Jiao J J, Moore W S, et al. Submarine groundwater discharge estimation in an urbanized embayment in Hong Kong via short-lived radium isotopes and its implication of nutrient loadings and primary production[J]. Marine Pollution Bulletin, 2014, 82(1/2): 144-154.

    [49] Xia Dong, Yu Zhigang, Xu Bochao, et al. Variations of hydrodynamics and submarine groundwater discharge in the Yellow River Estuary under the influence of the water-sediment regulation scheme[J]. Estuaries and Coasts, 2016, 39(2): 333-343.

    Submarine groundwater discharge into Sanggou Bay traced by223Ra and224Ra

    Wang Xilong1, Du Jinzhou1, Zhang Jing1

    (1.StateKeyLaboratoryofEstuarineandCoastalResearch,EastChinaNormalUniversity,Shanghai200062,China)

    Submarine groundwater discharge (SGD) with inputs of nutrients and water into the ocean in certain regions plays a significant role in coastal material fluxes and their biogeochemical cycle, which cannot be ignored for its impact on the eco-system. In this paper, SGD into the typical aquaculture base of northern China, Sanggou Bay, was estimated using naturally occurring radium isotopes,223Ra and224Ra. Results showed that the radium activities of submarine groundwater, especially in the pore water [224Ra=(968±31) dpm/(100 L),223Ra=(31.4±4.9) dpm/(100 L),n=9], were obviously greater than those in the surface water [224Ra=(38.7±2.0) dpm/(100 L),223Ra=(1.70±0.50) dpm/(100 L),n=21]. Assuming steady state and using a radium mass balance model with sources and sinks, the SGD rate was estimated to be among (0.23~1.03) ×107m3/d. SGD flux varied from spring tide to ebb tide. During spring tide, SGD flux was smaller because of weak hydraulic gradient while it was much larger during ebb tide because of strong hydraulic gradient. During a tidal cycle, the SGD flux was averaged to be 0.39×107m3/d based on223Ra and224Ra. Tidal-driven SGD can account for 61% of the total SGD. Therefore, the main driving force of SGD into Sanggou Bay was tidal pumping and SGD of Sanggou Bay in coastal zone contributes a lot to seawater composition and material exchange between land and sea.

    radium isotopes; submarine groundwater discharge; balance model; tidal pumping; Sanggou Bay

    10.3969/j.issn.0253-4193.2017.04.002

    2016-05-08;

    2016-11-28。

    科技部973項(xiàng)目“多重壓力下近海生態(tài)系統(tǒng)可持續(xù)產(chǎn)出與適應(yīng)性管理的科學(xué)基礎(chǔ)”(2011CB409801);國(guó)家自然科學(xué)基金項(xiàng)目“我國(guó)黃海海底地下水排泄及營(yíng)養(yǎng)鹽入海通量的研究”(41376089)。

    王希龍(1989—),女,山東省聊城市人,博士研究生,從事海洋同位素化學(xué)方向研究。E-mail:xuanfeng698547@126.com

    *通信作者:杜金洲,男,教授,從事同位素海洋學(xué)和環(huán)境放射化學(xué)的研究。E-mail:jzdu@sklec.ecnu.edu.cn

    P734.2

    A

    0253-4193(2017)04-0016-12

    王希龍,杜金洲,張經(jīng). 基于223Ra和224Ra的桑溝灣海底地下水排放通量[J].海洋學(xué)報(bào),2017,39(4):16—27,

    Wang Xilong, Du Jinzhou, Zhang Jing. Submarine groundwater discharge into Sanggou Bay traced by223Ra and224Ra[J]. Haiyang Xuebao,2017,39(4):16—27, doi:10.3969/j.issn.0253-4193.2017.04.002

    猜你喜歡
    潮汐活度表層
    潮汐與戰(zhàn)爭(zhēng)(上)
    半潛式平臺(tái)表層卡套管處理與認(rèn)識(shí)
    海洋石油(2021年3期)2021-11-05 07:43:10
    水體表層沉積物對(duì)磷的吸收及釋放研究進(jìn)展
    絕美海灘
    CaO-SiO2-FeO-P2O5-Al2O3脫磷渣系中組元活度的計(jì)算
    潮汐式灌溉控制系統(tǒng)的設(shè)計(jì)及應(yīng)用
    電子制作(2017年9期)2017-04-17 03:00:56
    干法紙的潮汐
    生活用紙(2016年6期)2017-01-19 07:36:25
    核電廠惰性氣體排放活度濃度的估算
    氬弧熔覆原位合成Ti(C,N)-WC增強(qiáng)鎳基表層復(fù)合材料的研究
    焊接(2015年6期)2015-07-18 11:02:25
    鉆井液處理劑溶液活度測(cè)量方法對(duì)比
    神马国产精品三级电影在线观看| 久久久久国产精品人妻aⅴ院| 国产伦一二天堂av在线观看| 亚洲中文日韩欧美视频| 国产探花极品一区二区| 国产成人freesex在线 | 精品国内亚洲2022精品成人| 国产精品人妻久久久久久| 久久亚洲国产成人精品v| 99久久九九国产精品国产免费| 欧美最新免费一区二区三区| 天天躁夜夜躁狠狠久久av| av.在线天堂| 丰满的人妻完整版| 国内精品久久久久精免费| 一个人免费在线观看电影| 亚洲综合色惰| 亚洲三级黄色毛片| 亚洲欧美日韩卡通动漫| 亚洲熟妇熟女久久| 日韩精品中文字幕看吧| 国内久久婷婷六月综合欲色啪| 韩国av在线不卡| 晚上一个人看的免费电影| 亚洲第一电影网av| 啦啦啦韩国在线观看视频| 日日摸夜夜添夜夜添av毛片| 亚洲色图av天堂| 一区二区三区高清视频在线| 久久人人精品亚洲av| 精品国产三级普通话版| 变态另类丝袜制服| 嫩草影院精品99| 99久久无色码亚洲精品果冻| 男女视频在线观看网站免费| 欧美色视频一区免费| 亚洲成a人片在线一区二区| 老司机福利观看| 久久午夜亚洲精品久久| 日本一二三区视频观看| 欧美性猛交╳xxx乱大交人| 亚洲精品粉嫩美女一区| 亚洲成人中文字幕在线播放| av在线观看视频网站免费| 国产精品亚洲一级av第二区| 真人做人爱边吃奶动态| 国产在线精品亚洲第一网站| 色视频www国产| 亚洲av免费高清在线观看| 国产麻豆成人av免费视频| 老师上课跳d突然被开到最大视频| 婷婷精品国产亚洲av| avwww免费| av视频在线观看入口| 国产单亲对白刺激| 免费无遮挡裸体视频| 看非洲黑人一级黄片| 最近中文字幕高清免费大全6| 免费看av在线观看网站| 久久久久免费精品人妻一区二区| 亚洲成人中文字幕在线播放| 久久久久国内视频| 啦啦啦观看免费观看视频高清| 永久网站在线| 狂野欧美激情性xxxx在线观看| 亚洲人与动物交配视频| 赤兔流量卡办理| 午夜福利高清视频| 舔av片在线| 欧美3d第一页| 国产白丝娇喘喷水9色精品| 欧美bdsm另类| 国产成人影院久久av| 不卡视频在线观看欧美| 久久国产乱子免费精品| 亚洲精品影视一区二区三区av| 亚洲国产精品国产精品| 国产亚洲精品久久久久久毛片| 国产精品一及| 能在线免费观看的黄片| 久久精品国产亚洲av涩爱 | 国产成人freesex在线 | 精品人妻一区二区三区麻豆 | 日本a在线网址| 蜜桃亚洲精品一区二区三区| 婷婷亚洲欧美| 久久中文看片网| 午夜福利视频1000在线观看| 又黄又爽又刺激的免费视频.| 久久精品91蜜桃| 久久人人精品亚洲av| 哪里可以看免费的av片| 深爱激情五月婷婷| 熟女人妻精品中文字幕| 男女那种视频在线观看| 亚洲人成网站在线播放欧美日韩| 99久久精品热视频| 丰满乱子伦码专区| 久久综合国产亚洲精品| 日韩欧美一区二区三区在线观看| av免费在线看不卡| 菩萨蛮人人尽说江南好唐韦庄 | 偷拍熟女少妇极品色| 亚洲专区国产一区二区| 国产免费男女视频| 中文资源天堂在线| 桃色一区二区三区在线观看| 国产av一区在线观看免费| 欧美高清成人免费视频www| 麻豆久久精品国产亚洲av| 看非洲黑人一级黄片| 日韩强制内射视频| 春色校园在线视频观看| 欧美成人免费av一区二区三区| 亚洲人成网站在线观看播放| 亚洲中文字幕日韩| www日本黄色视频网| 一级毛片电影观看 | 91在线观看av| 永久网站在线| 色综合亚洲欧美另类图片| 深爱激情五月婷婷| 国产成人aa在线观看| 乱码一卡2卡4卡精品| 国产精品一区二区三区四区免费观看 | 老熟妇仑乱视频hdxx| www.色视频.com| 两个人视频免费观看高清| 在线免费观看不下载黄p国产| 最近最新中文字幕大全电影3| 两性午夜刺激爽爽歪歪视频在线观看| 日韩精品青青久久久久久| 国产精品一及| 日本a在线网址| 久久精品国产亚洲网站| 午夜精品一区二区三区免费看| 人妻制服诱惑在线中文字幕| 我的老师免费观看完整版| 国产精品一及| 六月丁香七月| 老司机午夜福利在线观看视频| 国产 一区 欧美 日韩| 国产老妇女一区| 亚洲乱码一区二区免费版| 亚洲熟妇中文字幕五十中出| 综合色丁香网| 欧美zozozo另类| 久久久久国内视频| 美女内射精品一级片tv| 日本一二三区视频观看| 又爽又黄无遮挡网站| 亚洲人成网站在线播| 淫妇啪啪啪对白视频| 亚洲欧美日韩东京热| 一级黄色大片毛片| 22中文网久久字幕| 午夜福利18| 丰满人妻一区二区三区视频av| 日韩欧美一区二区三区在线观看| 成人午夜高清在线视频| 国产v大片淫在线免费观看| 国产亚洲精品综合一区在线观看| 国产精品一及| 成年女人毛片免费观看观看9| 亚洲av五月六月丁香网| 男女边吃奶边做爰视频| 国产久久久一区二区三区| 日本黄大片高清| 日韩强制内射视频| 男女那种视频在线观看| 日本与韩国留学比较| 国产精品久久视频播放| 日日摸夜夜添夜夜添av毛片| 免费不卡的大黄色大毛片视频在线观看 | 国产av麻豆久久久久久久| 日本与韩国留学比较| 欧美成人一区二区免费高清观看| 精品久久久久久久人妻蜜臀av| 2021天堂中文幕一二区在线观| 岛国在线免费视频观看| 99热这里只有是精品50| 97碰自拍视频| 免费观看人在逋| 国产aⅴ精品一区二区三区波| 欧美一区二区精品小视频在线| 免费大片18禁| 亚洲av成人精品一区久久| 人人妻,人人澡人人爽秒播| 亚洲不卡免费看| 日韩欧美 国产精品| 欧美性感艳星| 成人综合一区亚洲| 国产一区二区三区av在线 | 国产高清三级在线| 波多野结衣高清无吗| 美女xxoo啪啪120秒动态图| 亚洲不卡免费看| 国产爱豆传媒在线观看| 久久午夜福利片| 国产激情偷乱视频一区二区| 草草在线视频免费看| 简卡轻食公司| 小说图片视频综合网站| 天美传媒精品一区二区| av在线老鸭窝| 久久精品夜夜夜夜夜久久蜜豆| 好男人在线观看高清免费视频| 国产伦一二天堂av在线观看| 国国产精品蜜臀av免费| 中文在线观看免费www的网站| 日韩精品青青久久久久久| 亚洲高清免费不卡视频| 成人永久免费在线观看视频| 色综合站精品国产| 成人美女网站在线观看视频| 亚洲精品日韩av片在线观看| 一本一本综合久久| 在线观看66精品国产| 久久精品综合一区二区三区| 性欧美人与动物交配| 你懂的网址亚洲精品在线观看 | 国产黄a三级三级三级人| 91精品国产九色| 成年女人看的毛片在线观看| 亚洲婷婷狠狠爱综合网| 欧美激情国产日韩精品一区| 有码 亚洲区| 国产91av在线免费观看| 日本与韩国留学比较| 禁无遮挡网站| av福利片在线观看| 亚洲av一区综合| 成人美女网站在线观看视频| 日本 av在线| 一级黄片播放器| 高清午夜精品一区二区三区 | 成人亚洲精品av一区二区| 人妻夜夜爽99麻豆av| 免费搜索国产男女视频| 天天一区二区日本电影三级| 欧美性猛交╳xxx乱大交人| 亚洲av成人精品一区久久| av国产免费在线观看| 乱系列少妇在线播放| 欧美绝顶高潮抽搐喷水| 亚洲精品国产成人久久av| 亚洲图色成人| 成人三级黄色视频| 国产视频内射| 国产乱人视频| 一本一本综合久久| 国产亚洲精品久久久久久毛片| 亚洲av美国av| 一夜夜www| 国产免费一级a男人的天堂| 一级黄色大片毛片| 最好的美女福利视频网| 变态另类丝袜制服| 久久韩国三级中文字幕| 91在线精品国自产拍蜜月| 99riav亚洲国产免费| 欧美高清性xxxxhd video| 日本熟妇午夜| 婷婷六月久久综合丁香| 女人十人毛片免费观看3o分钟| 中文字幕精品亚洲无线码一区| 日韩成人av中文字幕在线观看 | 啦啦啦观看免费观看视频高清| 三级国产精品欧美在线观看| 久久久久久久久久黄片| 观看美女的网站| 国产真实乱freesex| 欧美高清成人免费视频www| 午夜福利在线在线| 国产精品人妻久久久影院| 亚洲最大成人手机在线| 成人亚洲精品av一区二区| 色5月婷婷丁香| 麻豆国产av国片精品| 99热这里只有是精品在线观看| 亚洲av中文字字幕乱码综合| 国产精品久久久久久亚洲av鲁大| 亚洲性夜色夜夜综合| 亚洲成a人片在线一区二区| 久久久久久久久久成人| 天天躁夜夜躁狠狠久久av| 精华霜和精华液先用哪个| 俄罗斯特黄特色一大片| 欧洲精品卡2卡3卡4卡5卡区| 亚洲自拍偷在线| 最近的中文字幕免费完整| 精品熟女少妇av免费看| 久久久久国内视频| 三级男女做爰猛烈吃奶摸视频| 亚洲av第一区精品v没综合| a级一级毛片免费在线观看| 日韩欧美三级三区| 十八禁国产超污无遮挡网站| 日韩欧美三级三区| av在线老鸭窝| 深夜精品福利| av在线老鸭窝| 深夜精品福利| 国产在视频线在精品| 少妇丰满av| 一级av片app| 又粗又爽又猛毛片免费看| 久久久色成人| 亚洲成人久久性| 久久久久国内视频| 看十八女毛片水多多多| 天天躁日日操中文字幕| 色哟哟哟哟哟哟| 国产成人影院久久av| 乱码一卡2卡4卡精品| 国产精品一区二区性色av| 免费看美女性在线毛片视频| 久久草成人影院| av.在线天堂| 亚洲av美国av| 少妇人妻精品综合一区二区 | a级毛色黄片| 变态另类成人亚洲欧美熟女| 少妇的逼好多水| 免费看av在线观看网站| 九九在线视频观看精品| 狂野欧美激情性xxxx在线观看| 国产精品久久电影中文字幕| 一区二区三区高清视频在线| 一级毛片久久久久久久久女| 国产大屁股一区二区在线视频| 日日啪夜夜撸| 六月丁香七月| 亚洲在线自拍视频| 少妇人妻一区二区三区视频| 国产成人影院久久av| 欧美日本视频| 男人舔女人下体高潮全视频| 精品欧美国产一区二区三| 久久精品国产亚洲av涩爱 | 免费看美女性在线毛片视频| 老师上课跳d突然被开到最大视频| 日韩人妻高清精品专区| av卡一久久| 在线国产一区二区在线| 国产高清激情床上av| 久久久国产成人免费| 在线国产一区二区在线| 色综合亚洲欧美另类图片| 精品无人区乱码1区二区| 精品福利观看| 综合色av麻豆| 别揉我奶头 嗯啊视频| 97碰自拍视频| 日韩欧美精品免费久久| 午夜免费男女啪啪视频观看 | 亚洲成a人片在线一区二区| 国产毛片a区久久久久| 国产av不卡久久| 热99re8久久精品国产| 国产精品1区2区在线观看.| 亚洲无线在线观看| 国产乱人偷精品视频| 亚洲精品亚洲一区二区| 国产精品久久视频播放| 日韩成人伦理影院| av免费在线看不卡| 女人被狂操c到高潮| 夜夜看夜夜爽夜夜摸| 一个人免费在线观看电影| 午夜精品一区二区三区免费看| 精品久久久久久久久亚洲| 亚洲国产精品成人久久小说 | 国产一级毛片七仙女欲春2| 校园人妻丝袜中文字幕| 成人美女网站在线观看视频| 3wmmmm亚洲av在线观看| 国产黄色视频一区二区在线观看 | 校园春色视频在线观看| 日韩一本色道免费dvd| 国产精品久久久久久av不卡| 免费av观看视频| 国产精品一及| 久久草成人影院| 精品人妻偷拍中文字幕| av卡一久久| 国产午夜精品论理片| 十八禁网站免费在线| 欧美高清性xxxxhd video| 人人妻,人人澡人人爽秒播| 男人的好看免费观看在线视频| 一个人看视频在线观看www免费| 日本精品一区二区三区蜜桃| 欧美成人精品欧美一级黄| 岛国在线免费视频观看| 99在线视频只有这里精品首页| 亚洲成人av在线免费| 国产欧美日韩精品亚洲av| 99久国产av精品| 人妻夜夜爽99麻豆av| 日韩av在线大香蕉| 久久人妻av系列| 精品人妻视频免费看| 国产片特级美女逼逼视频| 黄色日韩在线| 免费看美女性在线毛片视频| 国产伦精品一区二区三区四那| 国产不卡一卡二| 女人十人毛片免费观看3o分钟| 国产精品电影一区二区三区| 国产精品综合久久久久久久免费| 我的女老师完整版在线观看| 老司机午夜福利在线观看视频| 又粗又爽又猛毛片免费看| 老司机福利观看| 国语自产精品视频在线第100页| 色哟哟·www| 蜜桃久久精品国产亚洲av| 波多野结衣高清无吗| 国产黄色小视频在线观看| 人妻丰满熟妇av一区二区三区| 日本与韩国留学比较| 两个人视频免费观看高清| 亚洲国产日韩欧美精品在线观看| 欧美性感艳星| 国产黄片美女视频| 一本一本综合久久| 婷婷精品国产亚洲av在线| 亚洲精品456在线播放app| 国产高清激情床上av| 最近最新中文字幕大全电影3| 久久人人爽人人爽人人片va| 亚洲中文日韩欧美视频| 成人午夜高清在线视频| 成人特级黄色片久久久久久久| 在线免费十八禁| av天堂在线播放| 性欧美人与动物交配| 精品国产三级普通话版| 精品久久久久久久久久久久久| 丰满人妻一区二区三区视频av| 精品国内亚洲2022精品成人| 偷拍熟女少妇极品色| 色综合色国产| 亚洲精品国产av成人精品 | 亚洲精品456在线播放app| 搡女人真爽免费视频火全软件 | 搡女人真爽免费视频火全软件 | 久久99热这里只有精品18| 亚洲在线自拍视频| 国产精品精品国产色婷婷| 久久精品国产亚洲av涩爱 | 国产aⅴ精品一区二区三区波| 99在线人妻在线中文字幕| 亚洲国产日韩欧美精品在线观看| 日韩av在线大香蕉| 寂寞人妻少妇视频99o| 一级黄片播放器| 欧美日韩精品成人综合77777| 又爽又黄a免费视频| 哪里可以看免费的av片| 尾随美女入室| 一区二区三区免费毛片| 亚洲成人久久性| 欧美丝袜亚洲另类| 91av网一区二区| 99热这里只有精品一区| 香蕉av资源在线| 美女被艹到高潮喷水动态| 久久久精品欧美日韩精品| 婷婷亚洲欧美| 国内精品宾馆在线| 神马国产精品三级电影在线观看| 中文字幕免费在线视频6| 国产成年人精品一区二区| 亚洲精华国产精华液的使用体验 | 成熟少妇高潮喷水视频| 亚洲av免费在线观看| 免费一级毛片在线播放高清视频| 欧美一区二区亚洲| 小说图片视频综合网站| 久久精品夜色国产| 欧美日韩精品成人综合77777| 久久人妻av系列| or卡值多少钱| 免费无遮挡裸体视频| 精品欧美国产一区二区三| 欧美日本亚洲视频在线播放| 亚洲图色成人| 一级毛片电影观看 | 国产久久久一区二区三区| 天堂√8在线中文| 久久久久久久久久成人| 男插女下体视频免费在线播放| 国产爱豆传媒在线观看| 亚洲av电影不卡..在线观看| 长腿黑丝高跟| 国产免费男女视频| 免费观看的影片在线观看| 男人舔奶头视频| 蜜桃久久精品国产亚洲av| 亚洲内射少妇av| 日韩,欧美,国产一区二区三区 | 亚洲婷婷狠狠爱综合网| 亚洲国产高清在线一区二区三| 99热全是精品| 一进一出好大好爽视频| 国产精品人妻久久久久久| 男女视频在线观看网站免费| 精品人妻偷拍中文字幕| 亚洲av不卡在线观看| 免费不卡的大黄色大毛片视频在线观看 | 日韩精品中文字幕看吧| 老司机影院成人| 成人三级黄色视频| 精品久久久噜噜| 观看美女的网站| 亚洲第一区二区三区不卡| 最新在线观看一区二区三区| 精品乱码久久久久久99久播| 国产伦在线观看视频一区| 久久久久九九精品影院| 日产精品乱码卡一卡2卡三| 亚洲国产精品国产精品| 噜噜噜噜噜久久久久久91| 三级男女做爰猛烈吃奶摸视频| 亚洲人与动物交配视频| 九九久久精品国产亚洲av麻豆| 真人做人爱边吃奶动态| 欧美成人a在线观看| 男女啪啪激烈高潮av片| 国产精品永久免费网站| 国产精品日韩av在线免费观看| 网址你懂的国产日韩在线| 亚洲国产色片| 成人国产麻豆网| 美女被艹到高潮喷水动态| 日韩av在线大香蕉| 中文字幕av成人在线电影| 亚洲国产精品久久男人天堂| 露出奶头的视频| 一级黄色大片毛片| 亚洲成a人片在线一区二区| 国产精品,欧美在线| 欧美+日韩+精品| 一本一本综合久久| 精品久久久久久久久亚洲| 啦啦啦韩国在线观看视频| 欧美在线一区亚洲| 亚洲精品亚洲一区二区| 日韩欧美国产在线观看| 精品久久久久久久久av| 大香蕉久久网| 简卡轻食公司| 99久久精品热视频| 中文资源天堂在线| 99久国产av精品国产电影| 日韩欧美精品v在线| 亚洲欧美日韩高清在线视频| 69av精品久久久久久| 99久久久亚洲精品蜜臀av| 欧美成人免费av一区二区三区| 麻豆国产97在线/欧美| 人人妻人人澡人人爽人人夜夜 | 三级国产精品欧美在线观看| 成人亚洲精品av一区二区| 久久精品综合一区二区三区| 国产精品一区二区三区四区免费观看 | 国产精品1区2区在线观看.| 少妇熟女欧美另类| 亚洲av美国av| 久久久久国产精品人妻aⅴ院| 亚洲国产精品成人久久小说 | 国产成人影院久久av| 亚洲自拍偷在线| 青春草视频在线免费观看| 一个人观看的视频www高清免费观看| 嫩草影视91久久| 国产又黄又爽又无遮挡在线| 偷拍熟女少妇极品色| 国产免费一级a男人的天堂| 内地一区二区视频在线| 欧美性猛交黑人性爽| 国产精品久久久久久久电影| 久久99热6这里只有精品| 在线免费观看的www视频| 久久精品国产亚洲av涩爱 | 国产精品人妻久久久久久| 国产精华一区二区三区| 久久精品国产自在天天线| 午夜免费男女啪啪视频观看 | 国产三级在线视频| 久久久a久久爽久久v久久| 欧美日韩乱码在线| 1000部很黄的大片| 亚洲第一电影网av| av在线天堂中文字幕| 日本-黄色视频高清免费观看| 午夜精品国产一区二区电影 | 一级毛片我不卡| 日本色播在线视频| 一级毛片aaaaaa免费看小| av卡一久久| 国产午夜精品久久久久久一区二区三区 | 欧美一级a爱片免费观看看| 亚洲aⅴ乱码一区二区在线播放| 神马国产精品三级电影在线观看| 亚洲av成人精品一区久久| 久久韩国三级中文字幕| 2021天堂中文幕一二区在线观| 大又大粗又爽又黄少妇毛片口| videossex国产| 精华霜和精华液先用哪个| 午夜日韩欧美国产| 一本久久中文字幕|