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

    黃河流域?qū)幭亩蔚乇硭廴咎卣髋c風(fēng)險(xiǎn)評(píng)價(jià)

    2023-11-28 12:44:48趙增鋒付永亮邱小琮周瑞娟趙睿智楊強(qiáng)強(qiáng)
    中國(guó)環(huán)境科學(xué) 2023年11期
    關(guān)鍵詞:黃河流域寧夏濃度

    趙增鋒,付永亮,邱小琮,周瑞娟,趙睿智,楊強(qiáng)強(qiáng)

    黃河流域?qū)幭亩蔚乇硭廴咎卣髋c風(fēng)險(xiǎn)評(píng)價(jià)

    趙增鋒1,付永亮1,邱小琮2*,周瑞娟3,趙睿智1,楊強(qiáng)強(qiáng)1

    (1.寧夏大學(xué)土木與水利工程學(xué)院,寧夏 銀川 750021;2.寧夏大學(xué)生命科學(xué)學(xué)院,寧夏 銀川 750021;3.寧夏回族自治區(qū)生態(tài)環(huán)境監(jiān)測(cè)中心,寧夏 銀川 750000)

    為探明黃河流域?qū)幭亩蔚乇硭蟹臅r(shí)空變化特征及其影響因素,于2019年、2021年的春(4月)、夏(7月)、秋 (10月)對(duì)黃河流域?qū)幭亩?9個(gè)點(diǎn)位進(jìn)行6次采樣.將99個(gè)點(diǎn)位劃分為8個(gè)水系(黃河干流(S1)、清水河(S2)、苦水河(S3)、衛(wèi)寧段黃河北岸排水溝(S4)、衛(wèi)寧段黃河南岸排水溝(S5)、青石段黃河西岸排水溝(S6)、青石段黃河?xùn)|岸排水溝(S7)、城市景觀河湖(S8)),分析氟的時(shí)空分布特征,解析環(huán)境因子對(duì)氟富集的影響,并進(jìn)行生態(tài)和健康風(fēng)險(xiǎn)評(píng)價(jià).結(jié)果表明:黃河流域?qū)幭亩蔚乇硭鷿舛染禐?.46mg/L,高于中國(guó)飲用水標(biāo)準(zhǔn)值(1.0mg/L).受蒸發(fā)濃縮和降雨稀釋作用影響,氟濃度季節(jié)性變化表現(xiàn)為秋季>春季>夏季. 8個(gè)水系中只有黃河干流氟濃度未超標(biāo),苦水河、清水河和青石段黃河?xùn)|岸排水溝氟超標(biāo)最為嚴(yán)重. 黃河流域?qū)幭亩蔚乇硭鷳B(tài)風(fēng)險(xiǎn)均值達(dá)到低風(fēng)險(xiǎn)水平,中風(fēng)險(xiǎn)比率秋季為27.84%,春季為22.33%,夏季為11.11%,苦水河的氟生態(tài)風(fēng)險(xiǎn)最高,中風(fēng)險(xiǎn)比率高達(dá)80%.除黃河干流外,其他7個(gè)水系對(duì)兒童和成人均會(huì)產(chǎn)生健康風(fēng)險(xiǎn),且對(duì)兒童的危害作用遠(yuǎn)大于成人.氟的時(shí)空分布主要受自然因素影響,蒸發(fā)巖風(fēng)化溶解和蒸發(fā)濃縮是造成黃河流域?qū)幭亩畏患闹饕蛩?

    黃河流域;氟濃度;富集;生態(tài)風(fēng)險(xiǎn);健康風(fēng)險(xiǎn)

    氟(F)是自然界分布廣泛的鹵族元素,也是世界衛(wèi)生組織認(rèn)定的環(huán)境污染物,可以通過(guò)多種環(huán)境介質(zhì)(大氣、水體、生物體等)長(zhǎng)距離遷移或滯留[1].氟是動(dòng)物、植物和人類必不可少的微量營(yíng)養(yǎng)素[2],但過(guò)量的氟對(duì)幾乎所有的生物體都有毒性[3].據(jù)統(tǒng)計(jì),氟過(guò)量引起的人體中毒問(wèn)題已涉及全球50多個(gè)國(guó)家,有報(bào)道稱全球超過(guò)2.6億人受到氟中毒的影響,中國(guó)約有2200萬(wàn)~4500萬(wàn)人口受氟中毒影響[4-7].引起氟中毒的暴露途徑主要包括飲水[8-9]、飲茶[10]、皮膚接觸含氟水體或土壤[11-12]、食用含氟食物[9,13]等,其中飲水是造成氟中毒的最主要原因.因此,必須開展高氟水分布及成因的研究,以評(píng)估水安全,從而采取相應(yīng)措施降低氟的危害.由于地下水常作為飲用水源,與人類健康關(guān)系密切,因此關(guān)于水體中氟的研究以地下水居多,卻忽略了對(duì)地表水體中氟的研究.事實(shí)上,地表水體作為一種重要載體,地表水體中的生物在不同生長(zhǎng)階段都會(huì)受到氟毒性的影響,生態(tài)系統(tǒng)的結(jié)構(gòu)和功能也受到嚴(yán)重威脅.尤其作為飲用水源的地表水,氟污染甚至威脅到人類健康.因此,開展地表水體氟污染研究十分必要.

    黃河流域是我國(guó)重要的生態(tài)安全屏障和人口活動(dòng)的經(jīng)濟(jì)核心地帶,在我國(guó)經(jīng)濟(jì)社會(huì)發(fā)展和生態(tài)安全方面具有十分重要的地位[14].據(jù)寧夏回族自治區(qū)地表水環(huán)境質(zhì)量狀況月報(bào)顯示,氟化物是寧夏境內(nèi)黃河支流及入黃排水溝的主要污染指標(biāo)之一,其中以支流清水河、苦水河和都思兔河氟污染最為嚴(yán)重.寧夏也是地氟病多發(fā)之地,地方性氟中毒已經(jīng)成為一種嚴(yán)重危害當(dāng)?shù)鼐用裆眢w健康的地方病[15].此外,課題組研究還發(fā)現(xiàn)在寧夏高氟區(qū)水生生物的生物量和多樣性均較低.然而,對(duì)于黃河流域?qū)幭亩蔚乇硭w中氟化物的分布、影響因素及風(fēng)險(xiǎn)還缺乏系統(tǒng)的認(rèn)識(shí).本研究以黃河流域?qū)幭亩螢檠芯繀^(qū),分析該區(qū)域河流、湖泊及排水溝等地表水體中氟的污染特征,評(píng)價(jià)氟污染的生態(tài)和健康風(fēng)險(xiǎn).該研究是迄今為止黃河流域?qū)幭亩蔚乇硭廴疽?guī)模最大的研究,有助于更全面地了解高氟地表水的形成,為該區(qū)域氟污染管控提供依據(jù).

    1 材料與方法

    1.1 研究區(qū)域概況

    寧夏回族自治區(qū)35°14′~39°23′N, 104°17′~ 107°39′E,總面積6.64萬(wàn)km2,是唯一一個(gè)全境均屬于黃河流域的省區(qū).地形分為三大板塊:北部引黃灌區(qū)、中部干旱帶、南部山區(qū).地勢(shì)南高北低,呈階梯狀下降.最南端的六盤山區(qū)屬半濕潤(rùn)區(qū),衛(wèi)寧平原以北屬干旱區(qū),其他地區(qū)為半干旱區(qū)[16],為高氟水的形成創(chuàng)造了有利條件.

    1.2 樣品采集與預(yù)處理

    圖1 黃河流域?qū)幭亩尾蓸狱c(diǎn)布置

    黃河干流寧夏段以青銅峽大壩為界分為衛(wèi)寧段和青石段,根據(jù)海拔、地形、水系分布等情況,同時(shí)考慮寧夏渠道引水主要為黃河水,且大多進(jìn)行襯砌處理,渠道水水質(zhì)和黃河干流差異不大,將黃河流域?qū)幭亩蔚乇硭祫澐譃?個(gè)部分,共設(shè)置99個(gè)采樣點(diǎn)位(見圖1):分別為黃河干流(S1)7個(gè),清水河(S2)16個(gè),苦水河(S3)6個(gè),衛(wèi)寧段黃河北岸排水溝(S4)10個(gè),衛(wèi)寧段黃河南岸排水溝(S5)10個(gè),青石段黃河西岸排水溝(S6)29個(gè),青石段黃河?xùn)|岸排水溝(S7)13個(gè),城市景觀河湖(S8)8個(gè).冬季氣溫較低、冰層較厚,不利于野外大規(guī)模采樣,且冬季(枯水期)會(huì)有部分季節(jié)性水系斷流,因此選擇2019年和2021年春季(4月)、夏季(7月)和秋季(10月)進(jìn)行采樣調(diào)查,用1L采水器采集表層水下50cm深處的水樣,若水深<50cm,則采集表層的水樣,若水深>1.5m, 則采集表層下50cm與底層上50cm混合水樣.將采集的水樣帶回實(shí)驗(yàn)室,硝酸(pH<2)處理,用0.45 μm孔徑的微孔濾膜過(guò)濾,濾液于冰箱中4 ℃保存待測(cè).由于部分水系在采樣時(shí)斷流,共采集到580份樣品,其中2019年春季99份,夏季97份,秋季99份,2021年春季98份,夏季92份,秋季95份.

    1.3 環(huán)境因子測(cè)定

    采用YSI(Pro-plus)便攜式水質(zhì)分析儀現(xiàn)場(chǎng)測(cè)定水溫(WT)、電導(dǎo)率(EC)、總?cè)芙庑怨腆w(TDS)、氯離子(Cl-),水體中其他環(huán)境因子檢測(cè)方法為:pH值采用玻璃電極法(GB 6920-1986);氟化物(F-)采用離子選擇電極法(GB 7484-87);硫酸根離子(SO4-)采用重量法(GB 11899-89);化學(xué)需氧量(CODCr)采用重鉻酸鹽法(HJ 828-2017);總磷(TP)采用鉬酸銨分光光度法(GB/T 11893-1989),總堿度(TALK,以CaCO3計(jì))采用酸堿滴定法(HZ-HJ-SZ-0130),氨氮(NH3-N)采用納氏試劑分光光度法(HJ 535-2009).

    1.4 生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)模型

    生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)模型有熵值法、SSDs模型、PERA模型等,數(shù)據(jù)的獲得性是選擇評(píng)價(jià)模型的決定性因素之一[17].熵值法[18]對(duì)毒理數(shù)據(jù)要求較低,可操作性強(qiáng),因此采用熵值法對(duì)黃河流域?qū)幭亩畏廴緺顩r進(jìn)行量化表達(dá),計(jì)算模型見式(1):

    HQ=EC/WQC(1)

    式中:HQ表示危險(xiǎn)熵;EC表示氟暴露濃度;WQC表示水質(zhì)基準(zhǔn),取值為慢性基準(zhǔn)值,本研究參照鄭麗萍等[19]的研究,取2mg/L.由于缺乏氟的HQ標(biāo)準(zhǔn),根據(jù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)(GB 3838-2002)》[20]中不同水質(zhì)條件下污染物濃度的限值情況,選擇較為接近鋅的HQ作為本研究參照標(biāo)準(zhǔn).HQ 的具體范圍分別為 HQ<0.1,0.1£HQ<1,1£HQ<10,HQ310,分別代表評(píng)估區(qū)域無(wú)風(fēng)險(xiǎn)、低風(fēng)險(xiǎn)、中風(fēng)險(xiǎn)和高風(fēng)險(xiǎn)[21].

    1.5 健康風(fēng)險(xiǎn)評(píng)價(jià)模型

    重點(diǎn)考慮飲水?dāng)z入氟所導(dǎo)致的非致癌風(fēng)險(xiǎn)[22],采用美國(guó)EPA 推薦的非致癌評(píng)價(jià)風(fēng)險(xiǎn)模型[23],以兒童、成人為受體人群,對(duì)研究區(qū)氟的健康風(fēng)險(xiǎn)進(jìn)行評(píng)價(jià).計(jì)算模型見式(2)和(3):

    式中:CDI為單位體重日均暴露劑量,mg/kd;W為飲用水中非致癌物濃度,mg/L;IR為平均日攝水量,L/d; EF為暴露頻率,d/a;ED為暴露持續(xù)時(shí)間,a;ABW為平均體重,kg;AET為平均壽命,a;HI為個(gè)人平均非致癌風(fēng)險(xiǎn)指數(shù);RFD為非致癌物經(jīng)飲水途徑獲取的日攝入劑量,mg/(kg·d),各參數(shù)取值見表1.

    表1 健康風(fēng)險(xiǎn)評(píng)價(jià)參數(shù)取值[24]

    1.6 數(shù)據(jù)分析與繪圖

    采用Excel 2019計(jì)算氟的生態(tài)和健康風(fēng)險(xiǎn),采用SPSS 25進(jìn)行因子分析,采用Origin 2021繪制氟時(shí)空分布箱線圖和雙因子相關(guān)關(guān)系圖,采用ArcGIS 10.6繪制采樣點(diǎn)位圖和氟空間分布圖.

    2 結(jié)果與分析

    2.1 水體理化因子的時(shí)空分布

    水體理化指標(biāo)如表2所示.夏季W(wǎng)T最高,秋季最低.除2019年春季3個(gè)樣本呈弱酸性以外,其它樣本均為堿性水.EC值、TDS和TP濃度表現(xiàn)為春季>秋季>夏季,NH3-N濃度、TALK和F-濃度表現(xiàn)為秋季>春季>夏季,CODCr濃度和SO42-濃度表現(xiàn)為夏季>秋季>春季,Cl-濃度表現(xiàn)為夏季>春季>秋季.多數(shù)指標(biāo)秋季>夏季,可能是由于秋季干旱,受蒸發(fā)作用濃縮,夏季多雨,受降雨作用稀釋.CODCr、SO42-和Cl-在夏季濃度較高,表明夏季有大量可溶性硫酸鹽和氯化鹽溶解,同時(shí)表明夏季水體受有機(jī)物污染嚴(yán)重.TP在春季和夏季濃度較高,可能是由于春季和夏季的農(nóng)業(yè)生產(chǎn)活動(dòng)使磷肥通過(guò)農(nóng)田退水等形式進(jìn)入水體.除WT、pH值和TALK外,其它指標(biāo)的變異系數(shù)三個(gè)季節(jié)均屬于高等變異(CV>0.36),表明黃河流域?qū)幭亩尾煌c(diǎn)位間水質(zhì)差異較大,不同點(diǎn)位間降雨、蒸發(fā)等不平衡可能是造成差異的主要原因.

    表2 水體理化因子描述統(tǒng)計(jì)

    注:Min為最小值,Max為最大值,Mean為平均值,SD為標(biāo)準(zhǔn)差,CV表示變異系數(shù),WT單位為℃,pH值無(wú)量綱,EC單位為μS/cm,其余指標(biāo)單位為mg/L.

    2.2 氟時(shí)空分布特征

    黃河流域?qū)幭亩蔚乇硭鷿舛葧r(shí)空分布箱線圖見圖2,氟濃度空間分布見圖3.2019和2021兩年氟濃度均值為1.46mg/L,高于《生活飲用水衛(wèi)生標(biāo)準(zhǔn)(GB 5749-2022)》限值(1.0mg/L).圖2(a)可知,2019年氟濃度略高于2021年,可能是由于兩年的氣候存在差異,也可能是因?yàn)?021年因COVID-19人為活動(dòng)減少,表明人為活動(dòng)可能會(huì)影響地表水氟濃度.在同一年,氟濃度均值和超標(biāo)率季節(jié)上均表現(xiàn)為秋季>春季>夏季,且均值只有2021年夏季低于生活飲用水衛(wèi)生標(biāo)準(zhǔn)(GB 5749-2002)限值.秋季屬于旱季,F-受蒸發(fā)作用影響,濃縮富集,夏季為雨季,受降雨稀釋作用影響.氟超標(biāo)率、氟濃度和氟濃度均值的最大值均出現(xiàn)在2019年秋季,分別為68.69%、7.23mg/L和1.83mg/L,最小值均出現(xiàn)在2021年夏季,分別為29.35%、0.19mg/L和0.90mg/L.圖2(b)可知,8個(gè)水系中,2019年和2021兩年氟濃度均值表現(xiàn)為:S3> S7>S2>黃河流域?qū)幭亩?S4>S8>S5>S6>生活飲用水衛(wèi)生標(biāo)準(zhǔn)>S1.S3氟超標(biāo)率最大,為100%,其次為S2,為85.71%,只有S1未出現(xiàn)超標(biāo)現(xiàn)象.

    2.3 氟生態(tài)風(fēng)險(xiǎn)評(píng)估

    基于熵值法,以中國(guó)淡水生物慢性基準(zhǔn)值(2mg/L)為水質(zhì)基準(zhǔn),進(jìn)行了黃河流域?qū)幭亩蔚乇硭wF-的潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià),評(píng)價(jià)結(jié)果見圖4(a)和(b).黃河流域?qū)幭亩握w上面臨潛在的生態(tài)風(fēng)險(xiǎn),99%以上的樣本HQ30.1,HQ均值為0.73,屬于低風(fēng)險(xiǎn).風(fēng)險(xiǎn)級(jí)別在季節(jié)上表現(xiàn)為:秋季>春季>夏季,與F-濃度表現(xiàn)一致.秋季27.84%的樣本屬于中風(fēng)險(xiǎn),遠(yuǎn)大于夏季的11.11%,3個(gè)季節(jié)中只有夏季1.06%的樣本為無(wú)風(fēng)險(xiǎn).8個(gè)水系中,HQ均值表現(xiàn)為: S3>S7>S2>黃河流域?qū)幭亩尉?S4>S8>S5>S6>S1.S3的風(fēng)險(xiǎn)最大, 80%的樣本為中風(fēng)險(xiǎn),其次為S7和S2,中風(fēng)險(xiǎn)分別為37.18%和35.16%.S1的生態(tài)風(fēng)險(xiǎn)最小,沒(méi)有達(dá)到中風(fēng)險(xiǎn)的樣本,且無(wú)風(fēng)險(xiǎn)的樣本為4.76%.

    圖4 黃河流域?qū)幭亩蔚乇硭鷳B(tài)風(fēng)險(xiǎn)箱線圖

    2.4 氟健康風(fēng)險(xiǎn)評(píng)估

    通過(guò)口服攝入HI指數(shù)評(píng)估了黃河流域?qū)幭亩蔚乇硭w中F-對(duì)人類健康的潛在風(fēng)險(xiǎn),評(píng)價(jià)結(jié)果見圖5.盡管地表水體不直接作為飲用水源,但地表水和地下水高F-相伴而生,二者存在交換作用,雨季地表水補(bǔ)充地下水,旱季地下水補(bǔ)充地表水.在缺乏地下水F-濃度數(shù)據(jù)時(shí),對(duì)地表水F-健康風(fēng)險(xiǎn)評(píng)價(jià)也能很好的識(shí)別氟高風(fēng)險(xiǎn)區(qū).

    圖5 黃河流域?qū)幭亩蔚乇硭】碉L(fēng)險(xiǎn)

    由圖5(a)和(b)可知,2019和2021年對(duì)兒童的HI均值為1.38,高于人的健康風(fēng)險(xiǎn)閾值(1.0),成人的HI均值為0.78,低于人的健康風(fēng)險(xiǎn)閾值.HI指數(shù)的時(shí)空分布與氟濃度分布一致,季節(jié)上表現(xiàn)為秋季>春季>夏季.兒童HI指數(shù)超標(biāo)率為32.80%~63.40%,成人HI指數(shù)超標(biāo)率為12.70%~30.93%.由圖5(c)和(d)可知,8個(gè)水系中,兒童HI指數(shù)均值表現(xiàn)為:S3>S7>S2>黃河流域?qū)幭亩尉?S4>S8>S5>S6>HI標(biāo)準(zhǔn)>S1,只有S1對(duì)兒童不產(chǎn)生健康風(fēng)險(xiǎn).成人HI指數(shù)均值表現(xiàn)為:S3>S7>S2>HI標(biāo)準(zhǔn)>黃河流域?qū)幭亩尉?S4>S8>S5>S6>S1.成人和兒童HI指數(shù)超標(biāo)率均為S3最大,分別為80%和100%,S1最小,未出現(xiàn)超標(biāo)現(xiàn)象.

    3 討論

    3.1 F-生態(tài)與健康風(fēng)險(xiǎn)探析

    生態(tài)風(fēng)險(xiǎn)是由環(huán)境的自然變化或人類活動(dòng)引起的生態(tài)系統(tǒng)組成、結(jié)構(gòu)的改變而導(dǎo)致系統(tǒng)功能損失的可能性,而評(píng)價(jià)是定量預(yù)測(cè)各種風(fēng)險(xiǎn)源對(duì)生態(tài)系統(tǒng)產(chǎn)生風(fēng)險(xiǎn)的方法[17].健康風(fēng)險(xiǎn)評(píng)價(jià)是通過(guò)建立環(huán)境污染程度和人體健康的量化關(guān)系,評(píng)估影響人類健康疾病可能性的一種手段.生態(tài)風(fēng)險(xiǎn)和健康風(fēng)險(xiǎn)評(píng)價(jià)程序本質(zhì)一致,總體上分為4個(gè)步驟[12,17]:(1)危害識(shí)別;(2)暴露-效應(yīng)分析; (3)暴露評(píng)估; (4)風(fēng)險(xiǎn)表征.本研究中氟的生態(tài)風(fēng)險(xiǎn)是指氟暴露對(duì)本土水生生物所帶來(lái)的風(fēng)險(xiǎn),氟的健康風(fēng)險(xiǎn)是指飲水?dāng)z入氟對(duì)人體帶來(lái)的風(fēng)險(xiǎn),風(fēng)險(xiǎn)大小均與F-濃度分布表現(xiàn)一致.荊秀艷等[8]的研究認(rèn)為銀川平原地下水中的氟對(duì)成人和兒童的健康均會(huì)產(chǎn)生影響,且對(duì)兒童的危害作用遠(yuǎn)大于成人.寧夏的改水工程取得了較好的成效[25],但調(diào)查中發(fā)現(xiàn)仍有部分村民習(xí)慣性飲用窖水等水源,對(duì)氟中毒意識(shí)薄弱,因此應(yīng)加大宣傳力度,盡量避免直接飲用高氟水.

    高F-會(huì)對(duì)水生生物、微生物產(chǎn)生毒性,但水生生物及微生物在長(zhǎng)期的高F-環(huán)境中可以通過(guò)多種機(jī)制以適應(yīng)氟毒性,如微生物可通過(guò)吸附、礦化、吸收和積累、細(xì)胞外沉淀、酶轉(zhuǎn)化和細(xì)胞外排,以適應(yīng)氟毒性[26].本研究水質(zhì)基準(zhǔn)采用中國(guó)淡水生物慢性基準(zhǔn)值,HQ參照鋅的標(biāo)準(zhǔn),可能會(huì)造成評(píng)價(jià)結(jié)果寬大化,比如S1(黃河干流)氟濃度未超標(biāo),生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果卻達(dá)到低風(fēng)險(xiǎn)程度.其次,不同監(jiān)測(cè)點(diǎn)位間的水質(zhì)存在較大差異,如研究認(rèn)為鹽度會(huì)降低銅和鉛等重金屬的毒性[27-28],也有研究認(rèn)為鹽度過(guò)高或過(guò)低都會(huì)增大氨氮和亞硝酸鹽氮的毒性[29-30],然而還未見有關(guān)于鹽度等其他理化指標(biāo)對(duì)氟毒性起到協(xié)同或抑制的研究.因此今后應(yīng)采用該地區(qū)的優(yōu)勢(shì)淡水生物物種,同時(shí)考慮其他污染因子的協(xié)同或抑制作用進(jìn)行毒性試驗(yàn),以確定該地區(qū)的水生生物毒性基準(zhǔn)值,使評(píng)價(jià)結(jié)果更加客觀合理.同樣,目前很少人會(huì)直接飲用地表水或地下水,并且考慮飲水所攝入氟的健康風(fēng)險(xiǎn)在不同個(gè)體間也存在顯著差異.但是,在水源缺乏的旱區(qū),高氟水常被用于灌溉或牲畜飲用,氟可能會(huì)富集在糧食、蔬菜或肉類產(chǎn)品中,進(jìn)而對(duì)人類產(chǎn)生毒害作用.如邰蘇日嘎拉等[9]在固原市獲取的玉米和小麥存在氟含量大于國(guó)家規(guī)定的《食物中污染物的限量標(biāo)準(zhǔn)(GB 2762- 2005)》[31]的情況.事實(shí)上,兩種評(píng)價(jià)方法僅僅是區(qū)域氟污染程度的一種定量表征.通過(guò)分析結(jié)果識(shí)別區(qū)域氟污染狀況,識(shí)別氟高風(fēng)險(xiǎn)區(qū)及成因,為氟污染防控提供依據(jù).

    3.2 高F-水形成與富集機(jī)制

    地表水和地下水的化學(xué)成分基本相同,并由相同的地球化學(xué)過(guò)程控制[32].由Web of Science數(shù)據(jù)庫(kù)和CNKI數(shù)據(jù)庫(kù)可知,我國(guó)高氟地表水(>1mg/L)主要分布于北方干旱半干旱區(qū),如西北地區(qū)沿河西走廊的塔里木河流域、黑河流域、黃河流域的部分地區(qū),東部的海河流域、東北的西遼河流域、松花江流域的部分地區(qū).水體氟污染除少數(shù)地區(qū)是人為因素外,如采礦、冶煉、燃煤、陶瓷、玻璃、磚瓦、磷肥等生產(chǎn)過(guò)程產(chǎn)生的廢水、廢氣、廢渣通過(guò)多種形式進(jìn)入水體,造成水體氟污染[33],絕大多數(shù)水體氟污染是由于自然因素造成的[34].影響水體氟化物濃度的自然因素有巖石類型、氣候因素(蒸發(fā)、降水等)、海拔或起伏因素以及該地區(qū)的地質(zhì)和地貌等[35-36].水體中氟空間變異的主要水文地球化學(xué)過(guò)程包括溶解/沉淀、吸附/解吸和混合作用,受地質(zhì)條件及土壤淋溶較大的區(qū)域地表水氟濃度通常較高[37].以下從三個(gè)方面探討影響黃河流域?qū)幭亩蔚乇硭瓼-富集的因素.

    3.2.1 水系類型對(duì)地表水F-富集的影響 我國(guó)干旱內(nèi)陸水體均發(fā)源于平原周圍的山區(qū)[33].本研究8個(gè)水系中,S3(苦水河)、S2(清水河)和S7(青石段黃河?xùn)|岸排水溝)主要為自然水系,匯水面積大,水源主要為降雨和地下水補(bǔ)給,流量具有明顯的季節(jié)差異,并且干旱內(nèi)陸流域中下游地區(qū)廣泛分布著原生鹽漬化土壤、鹽土和風(fēng)沙土等,含有豐富的水溶性氟[5].降雨時(shí)對(duì)流域內(nèi)的土壤和巖石形成淋溶和沖刷作用,溶解在水體中的氟通過(guò)地表徑流的形式進(jìn)入水系,造成地表水氟濃度升高.高氟地表水和地下水往往是相伴而生的,在不同的水文條件下起到相互補(bǔ)充的作用,高氟地下水補(bǔ)充地表水時(shí),氟隨之遷移進(jìn)入地表水[38]. 地下水補(bǔ)給地表水的過(guò)程還會(huì)淋濾沿途含氟巖石,造成F-濃度進(jìn)一步升高.荊秀艷等[8]對(duì)銀川平原地下水氟的分布特征進(jìn)行了研究,結(jié)果表明S3和S7所處的區(qū)域氟濃度普遍大于1.5mg/L.邰蘇日嘎拉等[9]對(duì)固原市原州區(qū)高氟地區(qū)氟對(duì)人體健康風(fēng)險(xiǎn)進(jìn)行了研究,結(jié)果表明S2上游彭堡地區(qū)地下水氟濃度大于1.0mg/L.此外,謝殿榮[39]對(duì)寧夏地下水氟污染進(jìn)行了預(yù)測(cè),結(jié)果顯示S2、S3和S7所處的區(qū)域氟化物濃度大于1.5mg/L的概率達(dá)到80%以上,進(jìn)一步證實(shí)了自然水系中地表水氟濃度受高氟地下水補(bǔ)給的影響.盡管高氟水最終都匯入S1(黃河干流),但由于其水流量大,對(duì)氟起到稀釋作用,使得S1氟濃度最低.其它水系主要為人工水系,如城市景觀河湖、農(nóng)田退水排水溝等,匯水面積小,地表水和地下水交換弱,氟主要來(lái)源于沉積物的內(nèi)源釋放及農(nóng)田退水和生活污水中含的少量氟等,自然因素影響較小,氟濃度相對(duì)較低.

    3.2.2 地形地質(zhì)因素對(duì)地表水F-富集的影響 氟的富集是長(zhǎng)期地質(zhì)作用和地球化學(xué)演變的結(jié)果.我國(guó)地表水高氟區(qū)地形主要特征是地下水位埋藏淺、排泄不暢、交替緩慢,且地表水徑流流速慢,水巖作用時(shí)間長(zhǎng),有利于含氟物質(zhì)沉積和富集[38],氟濃度呈現(xiàn)地勢(shì)越低含量越高的分布規(guī)律[40].黃河流域?qū)幭亩胃叻w分布符合此規(guī)律,主要分布在海拔1100~ 1250m的平原地區(qū).氟主要來(lái)源是巖土中的含氟礦物[36],穆敏等[41]測(cè)定了寧夏中部干旱帶的粗顆粒礫石及黏土礦物成分,結(jié)果認(rèn)為粗顆粒礫石中石英、云母為主要成分,黏土礦物以綠泥石為主.這些礦物均具有較高的氟本底值,如天山、祁連山等構(gòu)造山系的云母類礦物中氟含量高達(dá)2.1%~5.2%[5].黏土中的氟含量較高[34],綠泥石作為寧夏中部干旱帶黏土的重要成分,它的釋放能力對(duì)水體氟含量具有較大影響.張曉燕等[42]對(duì)若爾蓋盆地表層沉積物的礦物組成及環(huán)境意義進(jìn)行了研究,結(jié)果認(rèn)為綠泥石受到基巖輸入和風(fēng)化的雙重影響,且容易發(fā)生水解,受流域物質(zhì)輸入的影響,河流沉積物中的氟含量略高于表層土壤.因此,對(duì)于黃河流域?qū)幭亩蝸?lái)說(shuō),巖石風(fēng)化后的產(chǎn)物,通過(guò)降雨匯流輸入河流,并通過(guò)水解作用大量釋放氟化物是該地區(qū)高氟水形成的最重要原因.地區(qū)的構(gòu)造對(duì)提高水體中氟化物的濃度也有重大影響,查寧夏回族自治區(qū)地質(zhì)圖可知,黃河流域?qū)幭亩尾蓸狱c(diǎn)位主要位于馬家灘褶斷束、走廊加里東褶皺帶和銀川地塹上,地層主要為第四系,其中氟化物濃度總體表現(xiàn)為馬家灘褶斷束>走廊加里東褶皺帶>銀川地塹.

    3.2.3 氣候及水的理化性質(zhì)對(duì)地表水F-富集的影響 對(duì)除pH值以外的10種理化因子進(jìn)行標(biāo)準(zhǔn)化,然后進(jìn)行主成分分析.經(jīng)KMO度量及Bartlett球型度檢驗(yàn),KMO檢驗(yàn)值為0.68,大于0.50,Bartlett球形度檢驗(yàn)顯著性為0.00,表明數(shù)據(jù)可以作PCA分析.選取特征值大于1的為主因子,共提取4個(gè)主因子(表2),可以解釋數(shù)據(jù)矩陣中總方差的62.85%,表明黃河流域?qū)幭亩嗡瘜W(xué)成分控制因素的復(fù)雜性.

    因子1解釋了總方差的28.25%,F-、EC、TDS、CODCr、Cl-和SO42-具有較高的載荷.對(duì)于氣候干旱、降水稀少的地區(qū),由大氣降水所輸入的可溶性離子可以忽略不計(jì),地表水化學(xué)組分主要來(lái)源于礦物的風(fēng)化溶解[43].離子關(guān)系是認(rèn)識(shí)水文地球化學(xué)過(guò)程的依據(jù),由表3可知,TDS與SO42-、Cl-呈極顯著正相關(guān)(<0.01),表明SO42-和Cl-是研究區(qū)地表水TDS的主要貢獻(xiàn)陰離子.TDS和Cl-的關(guān)系可以很好的指示蒸發(fā)作用[44].TDS與Cl-的極顯著正相關(guān)顯著正相關(guān)關(guān)系表明蒸發(fā)濃縮是該地區(qū)地表水成因的一個(gè)重要因素.蒸發(fā)和/或蒸發(fā)-蒸騰作用的過(guò)程會(huì)加劇地表水體和淺層地下水中F-的濃縮[45-47].研究區(qū)主要位于干旱半干旱區(qū),降雨量小于蒸發(fā)量.其中寧夏中部干旱帶年平均降水量180mm,蒸發(fā)量高達(dá)2100~2400mm[41],導(dǎo)致F-的富集.尤其在秋季,降雨少,蒸發(fā)大,加速F-的富集.自然條件下,地表水的Cl-和SO42-與巖石風(fēng)化有關(guān)[44]. Cl-主要來(lái)源于硅酸鹽巖和蒸發(fā)巖的風(fēng)化溶解,SO42-主要來(lái)源于蒸發(fā)巖礦物的溶解[43,48].三個(gè)季節(jié)中,F-與EC、TDS、Cl-、SO42-均呈極顯著正相關(guān)(<0.01)(表3),有研究也認(rèn)為TDS、EC與F-存在正相關(guān)關(guān)系[32,34],表明高鹽、高TDS對(duì)巖石中氟的釋放具有促進(jìn)作用.因此,因子1揭示了黃河流域?qū)幭亩嗡瘜W(xué)成分的主要控制因素為蒸發(fā)巖風(fēng)化和蒸發(fā)濃縮.

    因子2解釋了總方差的12.24%,WT、F-和CODCr具有較高的載荷,主要受季節(jié)變化的影響,季節(jié)的不同氣候不同,同時(shí)能反映氣候變化的影響. CODCr在因子2中呈正載荷,且在夏季與WT呈極顯著正相關(guān)(<0.01),表明夏季人類活動(dòng)強(qiáng)度大,導(dǎo)致大量有機(jī)污染物進(jìn)入地表水. F-在因子2中呈負(fù)載荷,且WT和F-在不同的季節(jié)相關(guān)性不強(qiáng),在年際上呈極顯著負(fù)相關(guān)(<0.01).本研究在季節(jié)上表現(xiàn)為平均WT越高的季節(jié)所對(duì)應(yīng)的流量越大,表明F-濃度受水體稀釋作用影響較大.不同的氣候條件下,原生礦物的分解速度、元素的生物地球化學(xué)循環(huán)強(qiáng)度和元素遷移的速度都具有明顯的差異[49].研究認(rèn)為溫差大,巖石風(fēng)化強(qiáng)烈,有利于含氟礦物的風(fēng)化裸露[50].研究區(qū)符合溫差大的特點(diǎn),且氣候差異明顯,呈南寒北暖、南濕北干的特點(diǎn).此外,凍融干濕交替循環(huán)作用會(huì)促進(jìn)粗顆粒礫石的裂化[41],干濕交替下隨著濕度的增加礫石元素釋放量升高[51].研究區(qū)多數(shù)水系為季節(jié)性水系,頻繁經(jīng)歷著干濕交替.凍融發(fā)生在枯水期,溫度低,也是WT和F-呈極顯著負(fù)相關(guān)的另一個(gè)重要原因.

    因子3解釋了總方差的11.62%,NH3-N、TP和TALK具有較高的載荷,主要表現(xiàn)為人類活動(dòng)的影響.如農(nóng)業(yè)生產(chǎn)活動(dòng)中使用的化肥農(nóng)藥通過(guò)農(nóng)業(yè)退水或土壤淋溶進(jìn)入地表水造成氟污染等[52],但普遍認(rèn)為人類活動(dòng)不是造成水體氟污染的主要因素[6,52].夏季F-與NH3-N呈顯著正相關(guān)(<0.05),可能是夏季頻繁的農(nóng)業(yè)生產(chǎn)活動(dòng)帶來(lái)了氟污染.淡水鹽漬化源于道路用鹽、人為加速風(fēng)化、污水排放、城市建設(shè)、肥料施用、礦山排水、資源開采、軟水劑使用、鹽水入侵以及水文變化和氣候變化等多種人為和地質(zhì)因素,Sujay等[53]將鹽離子與化學(xué)、生物和地質(zhì)參數(shù)之間的復(fù)雜相互關(guān)系以及對(duì)自然、社會(huì)和建筑環(huán)境的影響稱為“淡水鹽堿化綜合征(FSS)”.TALK在因子3中具有較高的正載荷,可能與人類活動(dòng)所引起的“FSS”有關(guān).

    因子4解釋了總方差的10.74%,pH值和TALK具有較高的載荷. TALK是CO32-、HCO3-、OH-及其他一些弱酸鹽類的總和.HCO3-主要來(lái)源于碳酸鹽巖和硅酸鹽巖的溶解[43,54],因此,因子4反映了碳酸鹽巖風(fēng)化的影響. pH值高的水中OH-可交換含氟礦物(黑云母、白云母、螢石)中的F-,從而增加水中F-濃度[32].較高的CO32-離子濃度起著催化劑的作用,可以增強(qiáng)F-從巖石中的溶解[34].大量研究表明, CO32-和HCO3-與F-會(huì)形成競(jìng)爭(zhēng)性解吸附,通過(guò)置換使F-釋出,直到達(dá)到解吸附平衡[52,55].春季TALK與F-呈極顯著正相關(guān)(<0.01),表明春季受競(jìng)爭(zhēng)性解吸附的影響較大.

    表2 主成分分析結(jié)果

    表3 黃河流域?qū)幭亩嗡w理化參數(shù)Pearson分析

    注:*表示顯著相關(guān)(<0.05),**表示極顯著相關(guān)(<0.01).

    4 結(jié)論

    4.1 黃河流域?qū)幭亩蔚乇硭鷿舛染禐?.46mg/L,高于《生活飲用水衛(wèi)生標(biāo)準(zhǔn)(GB 5749- 2022)》限值(1.0mg/L).氟濃度年際上表現(xiàn)為2019年>2021年,季節(jié)上表現(xiàn)為秋季(1.79mg/L)>春季(1.53mg/L)>夏季(1.04mg/L),水系類型上總體表現(xiàn)為自然水系>人工水系.苦水河、清水河和青石段黃河?xùn)|岸排水溝氟超標(biāo)最為嚴(yán)重,超標(biāo)率依次為100%、85.71%和58.97%,只有黃河干流未出現(xiàn)超標(biāo)現(xiàn)象.

    4.2 黃河流域?qū)幭亩蔚乇硭鷳B(tài)風(fēng)險(xiǎn)和健康風(fēng)險(xiǎn)的時(shí)空分布特征與氟濃度具有高度一致性.生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果表明: 8個(gè)水系均面臨潛在的生態(tài)風(fēng)險(xiǎn),99%以上的樣本達(dá)到低風(fēng)險(xiǎn)以上水平,其中苦水河的風(fēng)險(xiǎn)最大,80%的樣本為中風(fēng)險(xiǎn),黃河干流的生態(tài)風(fēng)險(xiǎn)最小.健康風(fēng)險(xiǎn)評(píng)價(jià)結(jié)果表明:除黃河干流外,其余7個(gè)水系的部分地區(qū)對(duì)兒童和成人均會(huì)產(chǎn)生健康風(fēng)險(xiǎn),且對(duì)兒童的危害作用遠(yuǎn)大于成人.

    4.3 造成黃河流域?qū)幭亩蔚乇硭廴镜闹饕蛩厥亲匀灰蛩?人為因素對(duì)氟污染貢獻(xiàn)不大.蒸發(fā)巖風(fēng)化溶解和蒸發(fā)濃縮是控制該地區(qū)氟富集的最主要因素.

    [1] 吳衛(wèi)紅.土-水-氣界面間氟的遷移機(jī)理及其生態(tài)效應(yīng)[D]. 杭州:浙江大學(xué), 2002.Wu W H. Migration mechanism and ecological effect of fluorine at soil-water-air interface [D]. Hangzhou: Zhejiang University, 2002.

    [2] Katiyar P, Pandey N, Sahu K K. Biological approaches of fluoride remediation: potential for environmental clean-up [J]. Environmental science and pollution research international, 2020,27(12):13044- 13055.

    [3] Calero P, Gurdo N, Nikel P I. Role of the CrcB transporter of Pseudomonas putida in the multi-level stress response elicited by mineral fluoride [J]. Environmental microbiology, 2022, 24(11):5082- 5104.

    [4] Mukherjee I, Singh U K. Exploring a variance decomposition approach integrated with the Monte Carlo method to evaluate groundwater fluoride exposure on the residents of a typical fluorosis endemic semi-arid tract of India [J]. Environmental Research, 2022, 203:111697.

    [5] 王根緒,程國(guó)棟.西北干旱區(qū)水中氟的分布規(guī)律及環(huán)境特征[J]. 地理科學(xué), 2000,(2):153-159.Wang G X, Cheng G D. Distribution and environmental characteristics of fluoride in water in arid area of Northwest China [J]. Scientia Geographica Sinica, 2000,(2):153-159.

    [6] 李鳳嫣,蔣天宇,余 濤,等.環(huán)境中氟的來(lái)源及健康風(fēng)險(xiǎn)評(píng)估研究進(jìn)展[J]. 巖礦測(cè)試, 2021,40(6):793-807. Li F Y, Jiang T Y, Xu T, et al. Review on sources of fluorine in the environment and health risk assessment [J]. Rock and Mineral Analysis, 2021,40(6):793-807.

    [7] Amini M, Mueller K, Abbaspour K C, et al. Statistical modeling of global geogenic fluoride contamination in groundwaters [J]. Environmental Science & Technology, 2008,42(10):3662-3668.

    [8] 荊秀艷,李小珍,王文姬,等.銀川平原地下水中氟分布特征及健康風(fēng)險(xiǎn)評(píng)價(jià)[J]. 環(huán)境科學(xué)與技術(shù), 2022,45(2):174-181. Jing X Y, Li X Z, Wang W J, et al. Distribution characteristics and health risk assessment of fluorine in groundwater in Yinchuan Plain [J]. Environmental Science and Technology, 2022,45(2):174-181.

    [9] 邰蘇日嘎拉,李永春,周文輝,等.寧夏固原市原州區(qū)高氟地區(qū)氟對(duì)人體健康的影響[J]. 巖礦測(cè)試, 2021,40(6):919-929. Tai S R G L, Li Y C, Zhou W H, et al. Effect of fluorine on human health in high-fluorine areas in Yuanzhou District, Guyuan City, Ningxia Autonomous Region [J]. Rock and Mineral Analysis, 2021, 40(6):919-929.

    [10] 陳 萍,魯 青,姜 泓,等.果洛州地方性氟骨癥流行現(xiàn)狀分析[J]. 中國(guó)地方病防治, 2023,38(2):110-112. Chen P, Lu Q, Jiang H, et al. Analysis of endemic skeletal fluorosis in Guoluo Prefecture [J]. Chinese Journal of Control of Endemic Diseases, 2023,38(2):110-112.

    [11] 劉靖宇,赫銀峰,耿建陽(yáng),等.典型干旱區(qū)綠洲土壤氟的遷移富集及健康風(fēng)險(xiǎn)評(píng)價(jià)[J]. 中國(guó)地質(zhì), 2023:1-16. Liu J Y, He Y F, Geng J Y, et al. Migration, enrichment and health risk assessment of soil fluoride in typical arid region of oasis [J]. Geology in China, 2023:1-16.

    [12] 李 淵.汾河流域飲用水源中氟和砷的分布特征及土地利用和植被變化的影響[D]. 太原:山西大學(xué), 2020. Li Y. Distribution characteristics of fluorine and arsenic in drinking water sources and their effects on land use and vegetation change in Fenhe River Basin. Taiyuan: Shanxi University, 2020.

    [13] 黃彤彤,王文臺(tái),楊俊林,等.食用燃煤烘烤后含氟辣椒中氟的人體吸收和健康風(fēng)險(xiǎn)閾值[J]. 生態(tài)毒理學(xué)報(bào), 2023:1-12. Huang T T, Wang W T, Yang J L, et al. Human absorption and health risk threshold of fluoride in fluorinated capsicum after roasting with edible coal [J]. Chinese Journal of Ecotoxicology, 2023:1-12.

    [14] 王淵釗.黃河流域?qū)幭亩紊鷳B(tài)環(huán)境保護(hù)與高質(zhì)量發(fā)展耦合協(xié)調(diào)度研究[D]. 蘭州:西北民族大學(xué), 2022. Wang Y Z. Study on coupling coordination degree of ecological environment protection and high quality development in Ningxia section of Yellow River Basin [D]. Lanzhou: Northwest University for Nationalities, 2022.

    [15] 于麗萍,孫 靜,陳阿麗,等.2014年寧夏飲水型地方性氟中毒情況分析[J]. 寧夏醫(yī)科大學(xué)學(xué)報(bào), 2017,39(3):290-292. Yu L P, Sun J, Chen A L, et al. Drinking water type endemic fluorosis in monitoring data of Ningxia in 2014 [J]. Journal of Ningxia Medical University, 2017,39(3):290-292.

    [16] 楊亞麗,王建英,緱曉輝,等.寧夏旅游氣候舒適度時(shí)空特征分析[J]. 寧夏大學(xué)學(xué)報(bào)(自然科學(xué)版), 2021,42(1):88-94. Yang Y L, Wang J Y, Gou X H, et al. Spatial and temporal characteristics of tourism climate comfort in Ningxia [J]. Journal of Ningxia University (Natural Science Edition), 2021,42(1):88-94.

    [17] 曾建軍,鄒明亮,郭建軍,等.生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)研究進(jìn)展綜述[J]. 環(huán)境監(jiān)測(cè)管理與技術(shù), 2017,29(1):1-5. Zeng J J, Zou M L, Guo J J, et al. Ecological risk assessment and its research progress [J]. Environmental Monitoring Management and Technology, 2017,29(1):1-5.

    [18] Liu X, Tu M, Wang S, et al. Research on freshwater water quality criteria, sediment quality criteria and ecological risk assessment of triclosan in China [J]. The Science of the Total Environment. 2021, 816:151616.

    [19] 鄭麗萍,龍 濤,鄧紹坡,等.我國(guó)氟化物的淡水水生生物基準(zhǔn)研究[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2015,31(6):923-927. Zheng L P, Long T, Deng S P, et al. Biological criteria of fluoride for freshwater hydrobiont in China [J]. Journal of Ecology and Rural Environment, 2015,31(6):923-927.

    [20] GB 3838-2002 地表水環(huán)境質(zhì)量標(biāo)準(zhǔn) [S]. 北京:中國(guó)標(biāo)準(zhǔn)出版社, 2002. GB 2762-2005 Environmental quality standards for surface water [S]. Beijing: Standards Press of China, 2002.

    [21] 劉宛昕,劉思彤,穆 迪,等.渤海灣鋅海水水質(zhì)基準(zhǔn)推導(dǎo)及潛在生態(tài)風(fēng)險(xiǎn)評(píng)價(jià)[J]. 環(huán)境化學(xué), 2023:1-12. Liu W X, Liu S T, Mu D, et al. Derivation and potential ecological risk assessment of zinc seawater in Bohai Bay [J]. Environmental Chemistry, 2023:1-12.

    [22] 宋曉光,蘆 巖,梁仕凱,等.張家口壩下地區(qū)高氟地下水成因分析與健康風(fēng)險(xiǎn)評(píng)價(jià)[J]. 地質(zhì)科技通報(bào), 2022,41(1):240-250. Song X G, Lu Y, Liang S K, et al. Analysis of high-fluoride groundwater formation mechanisms and assessment of health risk in Baxia region, Zhangjiakou [J]. Bulletin of Geological Science and Technology, 2022,41(1):240-250.

    [23] Usepa. Available information on assessment exposure from pesticides in food [R]. Washington D C:U.S.Programs, Environmental Protection Agency Office, 2000.

    [24] Wang W Z, Li Z, Su H, et al. Spatial and seasonal variability, control factors and health risk of fluoride in natural water in the Loess Plateau of China [J]. Journal of Hazardous Materials, 2022,434:128897.

    [25] 徐海明,劉明月,田 華,等.寧夏吳忠地區(qū)飲用水氟健康風(fēng)險(xiǎn)評(píng)價(jià)[J]. 環(huán)境與健康雜志, 2017,34(8):715-716 Xu H M, Liu M Y, Tian H, et al. Health risk assessment of drinking water fluoride in Wuzhong district of Ningxia [J]. Journal of Environment and Health, 2017,34(8):715-716.

    [26] Jianping Z, Hongmei W, Charles A C I, et al. Dissolution of fluorapatite byfluorescens P35 resulting in fluorine Release [J]. Geomicrobiology Journal, 2016,34(5):421-433.

    [27] 單陽(yáng)陽(yáng),李正炎.基于鹽度校正法的中國(guó)河口鉛水生生物水質(zhì)基準(zhǔn)制定[J]. 生態(tài)毒理學(xué)報(bào), 2020,15(5):272-287. Shan Y Y, Li Z Y. Derivation of water quality criteria of lead for the protection of chinese estuarine species based on salinity correction method [J]. Asian Journal of Ecotoxicology, 2020,15(5):272-287.

    [28] 馬超群,李正炎,胡 泓,等.基于鹽度校正的中國(guó)河口銅水生生物水質(zhì)基準(zhǔn)研究[J]. 中國(guó)海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2023,53(3): 82-97. Ma C Q, Li Z Y, Hu H, et al. Study of water quality criteria of copper for the protection of Chinese estuarine species based on salinity correction method [J]. Periodical of Ocean University of China, 2023,53(3):82-97.

    [29] 于建渤,黃學(xué)聃,陳 歡,等.不同鹽度條件下氨氮對(duì)斑節(jié)對(duì)蝦仔蝦的毒性研究[J]. 水產(chǎn)養(yǎng)殖, 2023,44(4):16-21. Yu J B, Huang X D, Chen H, et al. Toxicity of ammonia nitrogen onlarvae under different salinity conditions [J]. Journal of Aquaculture, 2023,44(4):16-21.

    [30] 崔宇航,黃學(xué)聃,陳 歡,等.不同鹽度條件下亞硝酸鹽氮對(duì)斑節(jié)對(duì)蝦仔蝦的毒性研究[J]. 水產(chǎn)養(yǎng)殖, 2022,43(11):1-5. Cui Y H, Huang X D, Chen H, et al. Toxicity of nitrite nitrogen onlarvae under different salinity conditions [J]. Journal of Aquaculture, 2022,43(11):1-5.

    [31] GB 2762-2005 食物中污染物的限量標(biāo)準(zhǔn) [S]. 北京:中國(guó)標(biāo)準(zhǔn)出版社, 2005. GB 2762-2005 Standard for limits of contaminants in food [S]. Beijing: Standards Press of China, 2005.

    [32] Jun X, Zhangdong J, Fei Z. Geochemical controls on fluoride concentrations in natural waters from the middle Loess Plateau, China [J]. Journal of Geochemical Exploration, 2015,159:252-261.

    [33] Meththika V, Prosun B. Fluoride in the environment: Sources, distribution and defluoridation [J]. Environmental Chemistry Letters, 2015,13(2):131-147.

    [34] Anirban C, Mrinal K A, Arnab M, et al. A critical review on geochemical and geological aspects of fluoride belts, fluorosis and natural materials and other sources for alternatives to fluoride exposure [J]. Journal of Hydrology, 2019,574:333-359.

    [35] Liu Y, Zhu W H. Environmental characteristics of regional groundwater in relation to fluoride poisoning in North China [J]. Environmental Geology and Water Sciences, 1991,18(1).

    [36] 何 錦,張福存,韓雙寶,等.中國(guó)北方高氟地下水分布特征和成因分析[J]. 中國(guó)地質(zhì), 2010,37(3):621-626. He J, Zhang F C, Han S B, et al. The distribution and genetic types of high-fluoride groundwater in northern China [J]. Geology in China, 2010,37(3):621-626.

    [37] 劉白薇.半干旱區(qū)水文地球化學(xué)演化規(guī)律及成因研究——以土默川平原為例[D]. 武漢:中國(guó)地質(zhì)大學(xué), 2019. Liu B W. Study on hydrogeochemical evolution and its genesis in semi-arid region: A case study of Tumochuan Plain [D]. Wuhan: China University of Geosciences, 2019.

    [38] 嵇曉燕,李 波,楊 凱,等.中國(guó)地表水氟化物時(shí)空分布特征初步研究[J]. 地球與環(huán)境, 2022,50(6):787-796. Ji X Y, Li B, Yang K, et al. Spatial and temporal distribution characteristics of fluoride in surface water of China [J]. Earth and Environment, 2022,50(6):787-796.

    [39] 謝殿榮.寧夏地下水原生氟污染空間分布預(yù)測(cè)[J]. 安全與環(huán)境工程, 2022,29(6):149-155. Xie D R. Spatial distribution prediction of primary fluoride pollution in groundwater in Ningxia Hui Autonomous Region [J]. Safety and Environmental Engineering, 2022,29(6):149-155.

    [40] 王連方,孫幸之,楊成忠,等.新疆奎屯一烏蘇山前傾斜平原地方性氟中毒的地理分布規(guī)律[J]. 中國(guó)地方病學(xué)雜志, 1982,(4):35-39. Wang L F, Sun X Z, Yang C Z, et al. Geographical distribution of endemic fluorosis in the front sloping plain of Kuitun-Wusu Mountain, Xinjiang [J]. Chinese Journal of Local Epidemiology, 1982,(4):35-39.

    [41] 穆 敏,李王成,王 潔,等.寧夏中部干旱帶粗顆粒土裂化過(guò)程研究[J]. 東北農(nóng)業(yè)大學(xué)學(xué)報(bào), 2022,53(6):88-96. Mu M, Li W C, Wang J, et al. Study on cracking process of coarse-grained soils in central dry zone of Ningxia [J]. Journal of Northeast Agricultural University, 2022,53(6):88-96.

    [42] 張曉燕,張春霞,李 佩.青藏高原東緣若爾蓋盆地表層沉積物粘土礦物組成及其環(huán)境意義探討[J]. 第四紀(jì)研究, 2022,42(2):435-448. Zhang X Y, Zhang C X, Li P. Composition and environmental significance of clay minerals in surface sediments of Ruoergai Basin, eastern margin of Qinghai-Tibet Plateau [J]. Quaternary Sciences. 2022,42(2):435-448.

    [43] 劉 敏,趙良元,李青云,等.長(zhǎng)江源區(qū)主要河流水化學(xué)特征、主要離子來(lái)源 [J]. 中國(guó)環(huán)境科學(xué), 2021,41(3):1243-1254.Liu M, Zhao L Y, LiQ Y, et al. Hydrochemical characteristics, main ion sources of main rivers in the source region of Yangtze River [J]. China Environmental Science, 2021,41(3):1243-1254.

    [44] 鄒嘉文,劉 飛,張靖坤.南水北調(diào)典型受水區(qū)淺層地下水水化學(xué)特征及成因 [J]. 中國(guó)環(huán)境科學(xué), 2022,42(5):2260-2268.Zou J W, Liu F, Zhang J K.Hydrochemical characteristics and formation mechanism of shallow groundwater in typical water- receiving areas of the South-to-North Water Diversion Project [J]. China Environmental Science, 2022,42(5):2260-2268

    [45] Mukherjee I, Singh U K. Groundwater fluoride contamination, probable release, and containment mechanisms: A review on Indian context [J]. Environmental Geochemistry and Health, 2018,40(6): 2259-2301.

    [46] Guo Q H, Wang Y X, Ma T, et al. Geochemical processes controlling the elevated fluoride concentrations in groundwaters of the Taiyuan basin, Northern China [J]. Journal of Geochemical Exploration, 2007, 93(1):1-12.

    [47] Rasool A, Farooqi A, Xiao T F, et al. A review of global outlook on fluoride contamination in groundwater with prominence on the Pakistan current situation [J]. Environmental Geochemistry and Health, 2018,40(4):1265-1281.

    [48] 王 建,韓海東,許君利,等.塔里木河流域出山徑流水化學(xué)特征研究 [J]. 中國(guó)環(huán)境科學(xué), 2021,41(4):1576-1587.Wang J, Han H D, Xu J L, et al. Hydrochemical characteristics of the mountain runoff in Tarim River Basin, China [J]. China Environmental Science, 2021,41(4):1576-1587.

    [49] 曹小虎.涑水盆地高氟地下水的分布及成因分析[J]. 中國(guó)水運(yùn)(下半月), 2011,11(8):190-191. Cao X H. Distribution and origin analysis of high-fluorine groundwater in the Sushui Basin [J]. China Water Transport (Second Half), 2011,11(8):190-191.

    [50] 邵琳琳,楊勝科,王文科,等.奎屯河流域水土中氟的分布規(guī)律[J]. 地球科學(xué)與環(huán)境學(xué)報(bào), 2006,(4):64-68. Shao L L, Yang S K, Wang W K, et al. Distribution regularity of fluorine in shallow groundwater in unsaturated soils of Kuitun River Basin, Xinjiang [J]. Journal of Earth Sciences and Environment, 2006,(4):64-68.

    [51] Liu Q, Li Y, Li W, et al. Effect of the release of gravel elements on soil nutrients and Jujube Fruit Yield under wet-and-dry cycles [J]. Agronomy, 2022,12(11):2881.

    [52] 羅文婷,張 鑫,高旭波.運(yùn)城盆地人類活動(dòng)氟污染的環(huán)境效應(yīng)研究[J]. 安全與環(huán)境工程, 2017,24(2):46-52. Luo W T, Zhang X, Gao X B. Impact of human activities on environmental fluoride pollution in Yuncheng Basin [J]. Safety and Environmental Engineering, 2017,24(2):46-52.

    [53] Sujay S K, Gene E L, Michael L P, et al. Freshwater salinization syndrome: from emerging global problem to managing risks [J]. Biogeochemistry, 2021,154(2):255-292.

    [54] 汪生斌,祁澤學(xué),王萬(wàn)平,等.格爾木河水化學(xué)特征及成因[J]. 水資源保護(hù), 2020,36(5):93-98. Wang S B, Qi Z X, Wang W W, et al. Hydrochemical characteristics and causes of formation of the Golmud River [J]. Water Resources Protection, 2020,36(5):93-98.

    [55] 韓占濤,張發(fā)旺,桂建業(yè),等.鹽池地區(qū)高氟地下水成因與低氟水找水方向研究[J]. 干旱區(qū)資源與環(huán)境, 2009,23(12):151-156. Han Z T, Zhang F W, Gui J Y, et al. Research on the genesis of high-fluoride groundwater in Yanchi area and direction for low- fluoride groundwater prospect [J]. Journal of Arid Land Resources and Environment, 2009,23(12):151-156.

    Characteristics and risk assessment of surface water fluorine pollution in Ningxia section of Yellow River Basin.

    ZHAO Zeng-feng1, FU Yong-liang1, QIU Xiao-cong2*, ZHOU Rui-juan3, ZHAO Rui-zhi1, YANG Qiang-qiang1

    (1.School of Civil and Hydraulic Engineering, Ningxia University, Yinchuan 750021, China;2.School of Life Science, Ningxia University, Yinchuan 750021, China;3.Ningxia Environment Monitoring Center, Yinchuan 750000, China).,2023,43(11):5800~5811

    In spring (April), summer (July) and autumn (October) of 2019 and 2021, six times were collected at 99 points in the Ningxia section of the Yellow River Basin. The 99 points were divided into 8 water systems which including the main stream of the Yellow River (S1), Qingshui River (S2), Kushui River (S3), Drainage ditch on the north bank of the Yellow River in Weining Section (S4), drainage ditch on the south bank of the Yellow River in Weining Section (S5), drainage ditch on the west bank of the Yellow River in Qingshi Section (S6), drainage ditch on the east bank of the Yellow River in Qingshi Section (S7), and urban landscape River and lake (S8), and the temporal and spatial distribution of fluorine was analyzed, the influence of environmental factors on fluorine enrichment was analyzed, and ecological and health risk assessment was carried out. The results showed that the average fluoride concentration of surface water in Ningxia section of the Yellow River Basin was 1.46mg/L, higher than the standard value of drinking water in China (1.0mg/L). Influenced by evaporation concentration and rainfall dilution, the fluorine concentration changed from high to low in autumn, spring and summer. The exceeding fluorine was the most serious in the Kushui River, Qingshui River and drainage ditch on the east bank of the Yellow River in Qingshi Section, and the main stream of the Yellow River was the only water system that was not polluted by fluoride. The average ecological risk of surface water fluoride in Ningxia section of the Yellow River Basin reached the low risk level, and the medium risk ratio was 27.84% in autumn, 22.33% in spring, and 11.11% in summer. The highest ecological risk of fluoride was found in the Kushui River, with the medium risk ratio reaching 80%. The results of the health risk assessment showed that except for the main stream of the Yellow River, other seven water systems posed health risks to both children and adults, and the harm to children was much greater than that to adults. The temporal and spatial distribution of fluorine was mainly affected by natural factors, and the dissolution of evaporite weathering and evaporation concentration were the main factors of fluorine enrichment in this area.

    yellow river basin;fluorine concentration;enrichment;ecological risk;health risk

    X522

    A

    1000-6923(2023)11-5800-12

    趙增鋒(1994-),男,河南商丘人,寧夏大學(xué)博士研究生,研究方向?yàn)樗料到y(tǒng)氟的環(huán)境行為及其與微生物互作機(jī)制.發(fā)表論文10余篇. 1402809124@qq.com.

    趙增鋒,付永亮,邱小琮,等.黃河流域?qū)幭亩蔚乇硭廴咎卣髋c風(fēng)險(xiǎn)評(píng)價(jià) [J]. 中國(guó)環(huán)境科學(xué), 2023,43(11):5800-5811.

    Zhao Z F, Fu Y L, Qiu X C, et al. Characteristics and risk assessment of surface water fluorine pollution in Ningxia section of Yellow River Basin [J]. China Environmental Science, 2023,43(11):5800-5811.

    2023-04-08

    寧夏高等學(xué)校一流學(xué)科建設(shè)(水利工程)項(xiàng)目(NXYLXK2021A03);寧夏回族自治區(qū)生態(tài)環(huán)境廳科研項(xiàng)目(2022015)

    * 責(zé)任作者, 教授, qiu_xc@nxu.edu.cn

    猜你喜歡
    黃河流域寧夏濃度
    生態(tài)環(huán)境部啟動(dòng)新一年度黃河流域“清廢行動(dòng)”
    拍照濃度有待降低
    神醫(yī)的煩惱——濃度與配比
    黃河流域燦爛的齊家文化——陶器
    收藏界(2019年2期)2019-10-12 08:26:10
    寧夏
    增強(qiáng)大局意識(shí) 提升黃河流域生態(tài)保護(hù)發(fā)展水平
    寧夏煤電博弈
    能源(2017年5期)2017-07-06 09:25:50
    三十六計(jì)之順手牽羊
    寧夏
    對(duì)改良的三種最小抑菌濃度測(cè)試方法的探討
    免费久久久久久久精品成人欧美视频| 精品国产超薄肉色丝袜足j| 日本五十路高清| 午夜福利,免费看| 麻豆av在线久日| 91麻豆av在线| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜福利在线观看吧| 狂野欧美激情性xxxx| 在线国产一区二区在线| 婷婷精品国产亚洲av在线| www.999成人在线观看| 这个男人来自地球电影免费观看| 国产高清videossex| 超碰成人久久| 啦啦啦免费观看视频1| 国内毛片毛片毛片毛片毛片| 日本撒尿小便嘘嘘汇集6| 国产精品野战在线观看| 神马国产精品三级电影在线观看 | av视频免费观看在线观看| av在线天堂中文字幕| 亚洲人成电影免费在线| 国产激情欧美一区二区| 日韩免费av在线播放| 欧美黑人欧美精品刺激| 黄片大片在线免费观看| 一级毛片女人18水好多| 欧美+亚洲+日韩+国产| 国产午夜精品久久久久久| 亚洲色图 男人天堂 中文字幕| av天堂久久9| 久热爱精品视频在线9| 精品一品国产午夜福利视频| 丝袜人妻中文字幕| 高潮久久久久久久久久久不卡| 精品少妇一区二区三区视频日本电影| 在线免费观看的www视频| 不卡av一区二区三区| www.自偷自拍.com| 一个人免费在线观看的高清视频| 欧美中文日本在线观看视频| 国产精品国产高清国产av| 国产在线观看jvid| 一个人观看的视频www高清免费观看 | 国产精品亚洲美女久久久| 此物有八面人人有两片| 日本免费a在线| 亚洲一区高清亚洲精品| 亚洲色图av天堂| 一区福利在线观看| 一区二区三区精品91| 午夜久久久久精精品| 少妇 在线观看| 亚洲成人精品中文字幕电影| 国产高清有码在线观看视频 | 男女之事视频高清在线观看| 国产精品亚洲av一区麻豆| 久久精品亚洲精品国产色婷小说| 国产成人系列免费观看| 满18在线观看网站| www.熟女人妻精品国产| 黄色 视频免费看| 天堂√8在线中文| 午夜精品在线福利| 久久久久久久久中文| 欧美+亚洲+日韩+国产| 国产一区二区三区在线臀色熟女| 亚洲一区高清亚洲精品| 成人国产综合亚洲| 欧美黑人欧美精品刺激| 99在线视频只有这里精品首页| 久久影院123| 亚洲熟女毛片儿| 国产精品,欧美在线| 丰满的人妻完整版| 一个人观看的视频www高清免费观看 | 黄色成人免费大全| 99热只有精品国产| 国产熟女午夜一区二区三区| 精品久久久精品久久久| av视频免费观看在线观看| 免费看十八禁软件| 精品久久久久久久人妻蜜臀av | 一区二区三区精品91| 91老司机精品| 亚洲伊人色综图| 久久人妻福利社区极品人妻图片| 久久久久精品国产欧美久久久| 亚洲人成网站在线播放欧美日韩| 欧美激情 高清一区二区三区| 成人手机av| 亚洲欧美日韩高清在线视频| 亚洲 国产 在线| 久久久久久大精品| 国产精品亚洲av一区麻豆| 女性被躁到高潮视频| 亚洲国产精品999在线| 18禁美女被吸乳视频| 一卡2卡三卡四卡精品乱码亚洲| 成年版毛片免费区| 久久久水蜜桃国产精品网| 久久午夜亚洲精品久久| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产精品999在线| 人成视频在线观看免费观看| 欧美久久黑人一区二区| 日韩精品青青久久久久久| 久久中文字幕人妻熟女| 亚洲第一欧美日韩一区二区三区| 亚洲成人免费电影在线观看| 日韩欧美三级三区| 女警被强在线播放| 久久久久国产精品人妻aⅴ院| 亚洲精品久久成人aⅴ小说| 亚洲精品久久国产高清桃花| 色播在线永久视频| 精品人妻在线不人妻| 这个男人来自地球电影免费观看| 亚洲av片天天在线观看| 久久精品aⅴ一区二区三区四区| 99riav亚洲国产免费| 午夜免费鲁丝| 亚洲av成人不卡在线观看播放网| 男女下面插进去视频免费观看| 色尼玛亚洲综合影院| 日韩中文字幕欧美一区二区| 国产精品自产拍在线观看55亚洲| 中出人妻视频一区二区| 999精品在线视频| 亚洲人成77777在线视频| 少妇的丰满在线观看| 精品福利观看| 国产成年人精品一区二区| 国产国语露脸激情在线看| 国产熟女xx| 国产精品香港三级国产av潘金莲| 精品欧美国产一区二区三| 真人做人爱边吃奶动态| 在线观看一区二区三区| 一区在线观看完整版| 国产亚洲精品综合一区在线观看 | www.自偷自拍.com| 免费女性裸体啪啪无遮挡网站| www.精华液| 精品国产超薄肉色丝袜足j| 亚洲无线在线观看| 亚洲第一电影网av| 日韩免费av在线播放| 午夜亚洲福利在线播放| 欧美一级a爱片免费观看看 | 国产欧美日韩一区二区三区在线| 大型av网站在线播放| 91大片在线观看| 午夜福利高清视频| 亚洲国产精品成人综合色| 成人18禁高潮啪啪吃奶动态图| 老汉色∧v一级毛片| or卡值多少钱| 亚洲精品国产一区二区精华液| 午夜两性在线视频| 黄色视频,在线免费观看| 欧美中文综合在线视频| 97碰自拍视频| 午夜免费激情av| 国产成人免费无遮挡视频| 久久久国产欧美日韩av| 久久九九热精品免费| 在线观看舔阴道视频| 国产av在哪里看| 香蕉丝袜av| 多毛熟女@视频| 久久香蕉精品热| 国产一区二区激情短视频| 狂野欧美激情性xxxx| 亚洲中文字幕日韩| 国产伦人伦偷精品视频| 欧美最黄视频在线播放免费| 亚洲国产精品999在线| 夜夜爽天天搞| 一级毛片高清免费大全| 99re在线观看精品视频| 99香蕉大伊视频| 超碰成人久久| 国产一区二区三区综合在线观看| 久久久国产成人精品二区| 国产亚洲av高清不卡| 国产精品日韩av在线免费观看 | 日韩欧美一区视频在线观看| 国产麻豆成人av免费视频| 欧美激情久久久久久爽电影 | 成人永久免费在线观看视频| 亚洲av成人av| 中文字幕色久视频| 国产欧美日韩一区二区三区在线| 这个男人来自地球电影免费观看| 免费在线观看视频国产中文字幕亚洲| 老汉色av国产亚洲站长工具| 午夜久久久在线观看| 一进一出抽搐gif免费好疼| 18禁美女被吸乳视频| 淫妇啪啪啪对白视频| 久久人人精品亚洲av| 日韩大码丰满熟妇| 色老头精品视频在线观看| 老汉色av国产亚洲站长工具| 色尼玛亚洲综合影院| 韩国精品一区二区三区| 青草久久国产| 亚洲一区中文字幕在线| 欧美日韩中文字幕国产精品一区二区三区 | 欧美绝顶高潮抽搐喷水| 亚洲熟妇中文字幕五十中出| 亚洲第一电影网av| 人人妻,人人澡人人爽秒播| 久久精品成人免费网站| 日韩精品中文字幕看吧| 国产精品精品国产色婷婷| 国产精品永久免费网站| 美女大奶头视频| 亚洲成av人片免费观看| 成年版毛片免费区| 91av网站免费观看| 淫妇啪啪啪对白视频| 国产在线观看jvid| 日本撒尿小便嘘嘘汇集6| 中文亚洲av片在线观看爽| 黄色视频,在线免费观看| 两个人看的免费小视频| avwww免费| 国产欧美日韩综合在线一区二区| 精品国产亚洲在线| а√天堂www在线а√下载| 国产99久久九九免费精品| 午夜福利欧美成人| x7x7x7水蜜桃| 亚洲国产高清在线一区二区三 | 日韩大码丰满熟妇| 亚洲色图av天堂| 精品久久久久久久毛片微露脸| 亚洲欧洲精品一区二区精品久久久| 亚洲伊人色综图| 禁无遮挡网站| 老熟妇乱子伦视频在线观看| 91av网站免费观看| 变态另类丝袜制服| 国产精品永久免费网站| 久久久国产欧美日韩av| 高清毛片免费观看视频网站| 日韩欧美在线二视频| aaaaa片日本免费| 在线观看免费日韩欧美大片| 9色porny在线观看| 不卡一级毛片| 国产欧美日韩精品亚洲av| tocl精华| 欧美久久黑人一区二区| 久久久久亚洲av毛片大全| av超薄肉色丝袜交足视频| 久久国产精品人妻蜜桃| 美国免费a级毛片| 变态另类成人亚洲欧美熟女 | 视频区欧美日本亚洲| 99精品久久久久人妻精品| 在线av久久热| 9色porny在线观看| 美国免费a级毛片| 黑人巨大精品欧美一区二区mp4| 久热爱精品视频在线9| 国产亚洲av高清不卡| 此物有八面人人有两片| 亚洲男人天堂网一区| 久久久水蜜桃国产精品网| 极品教师在线免费播放| 精品一区二区三区四区五区乱码| 欧美大码av| 看片在线看免费视频| 又紧又爽又黄一区二区| 999精品在线视频| 亚洲熟妇中文字幕五十中出| 国产av又大| 国产激情欧美一区二区| 国产欧美日韩一区二区三区在线| 国产一区二区三区在线臀色熟女| aaaaa片日本免费| 午夜福利影视在线免费观看| 亚洲国产精品sss在线观看| 欧美日韩亚洲综合一区二区三区_| 国产精品亚洲av一区麻豆| 中文字幕精品免费在线观看视频| 欧美黑人精品巨大| 成年女人毛片免费观看观看9| 999久久久国产精品视频| 青草久久国产| 精品国产亚洲在线| 夜夜爽天天搞| 亚洲午夜精品一区,二区,三区| 亚洲狠狠婷婷综合久久图片| 亚洲av片天天在线观看| 国产成人精品在线电影| 中文字幕久久专区| 长腿黑丝高跟| 淫妇啪啪啪对白视频| 国产精品久久久久久精品电影 | 好看av亚洲va欧美ⅴa在| 黄色视频不卡| 精品久久久久久久久久免费视频| 一进一出好大好爽视频| 国产人伦9x9x在线观看| 欧美午夜高清在线| 成人亚洲精品av一区二区| 满18在线观看网站| 亚洲专区字幕在线| 丝袜美腿诱惑在线| 无限看片的www在线观看| videosex国产| 精品久久久久久久人妻蜜臀av | 日本三级黄在线观看| 免费少妇av软件| 免费不卡黄色视频| 波多野结衣高清无吗| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人欧美| av福利片在线| 伊人久久大香线蕉亚洲五| 99在线视频只有这里精品首页| 久久香蕉精品热| 中亚洲国语对白在线视频| 国产精品国产高清国产av| 国产伦人伦偷精品视频| 欧美精品啪啪一区二区三区| 大型av网站在线播放| 欧美绝顶高潮抽搐喷水| 亚洲九九香蕉| 不卡av一区二区三区| 欧美激情高清一区二区三区| 美国免费a级毛片| 国产一区二区三区在线臀色熟女| 黑人操中国人逼视频| 久久久久久久午夜电影| 精品久久久久久久人妻蜜臀av | 日韩欧美国产在线观看| 亚洲视频免费观看视频| 久久影院123| 久久精品国产清高在天天线| 欧美久久黑人一区二区| 亚洲第一青青草原| 午夜免费鲁丝| 欧美激情极品国产一区二区三区| 麻豆一二三区av精品| 亚洲熟妇中文字幕五十中出| 可以在线观看毛片的网站| 久久精品影院6| 高潮久久久久久久久久久不卡| 美女高潮喷水抽搐中文字幕| 久久青草综合色| 亚洲久久久国产精品| 精品第一国产精品| 欧美黄色片欧美黄色片| 人妻丰满熟妇av一区二区三区| 亚洲,欧美精品.| 午夜久久久久精精品| 1024香蕉在线观看| 久久人人爽av亚洲精品天堂| 中文字幕人妻丝袜一区二区| 久久青草综合色| 欧美日本中文国产一区发布| 一级毛片精品| 精品久久久久久久毛片微露脸| 国产成人精品久久二区二区免费| 丝袜美腿诱惑在线| 岛国在线观看网站| 免费人成视频x8x8入口观看| 欧美精品亚洲一区二区| 高清毛片免费观看视频网站| 久久婷婷成人综合色麻豆| 久9热在线精品视频| 精品一区二区三区视频在线观看免费| 一区在线观看完整版| 少妇粗大呻吟视频| 国产成人av教育| 伊人久久大香线蕉亚洲五| 热99re8久久精品国产| 欧美+亚洲+日韩+国产| 久久这里只有精品19| 国产精品影院久久| 中文字幕色久视频| 日韩欧美一区视频在线观看| 国产日韩一区二区三区精品不卡| 天堂√8在线中文| 亚洲精品av麻豆狂野| 亚洲精品国产精品久久久不卡| 国产1区2区3区精品| 露出奶头的视频| 9191精品国产免费久久| а√天堂www在线а√下载| 亚洲一区中文字幕在线| 久久精品成人免费网站| 色哟哟哟哟哟哟| 久久香蕉激情| 人人妻,人人澡人人爽秒播| 夜夜夜夜夜久久久久| 亚洲成人精品中文字幕电影| 亚洲男人的天堂狠狠| 久久亚洲精品不卡| 亚洲情色 制服丝袜| 国产精品二区激情视频| 精品久久久久久久人妻蜜臀av | 亚洲 欧美一区二区三区| 国产一区二区在线av高清观看| 久久久久久久久中文| 国产精品av久久久久免费| 一本综合久久免费| av网站免费在线观看视频| 午夜激情av网站| 又大又爽又粗| 纯流量卡能插随身wifi吗| 国内精品久久久久久久电影| 欧美av亚洲av综合av国产av| svipshipincom国产片| 宅男免费午夜| 一级毛片精品| 欧美激情高清一区二区三区| 97人妻天天添夜夜摸| 亚洲精品久久国产高清桃花| 成人av一区二区三区在线看| av片东京热男人的天堂| avwww免费| 欧美大码av| АⅤ资源中文在线天堂| 在线国产一区二区在线| 两个人视频免费观看高清| 一级a爱视频在线免费观看| 十八禁人妻一区二区| 亚洲一区二区三区色噜噜| 麻豆国产av国片精品| 亚洲全国av大片| 国产精品秋霞免费鲁丝片| 亚洲国产高清在线一区二区三 | 久久久久九九精品影院| 亚洲专区国产一区二区| 一级片免费观看大全| 亚洲电影在线观看av| 操出白浆在线播放| 亚洲国产精品sss在线观看| √禁漫天堂资源中文www| 亚洲av美国av| 9191精品国产免费久久| 人成视频在线观看免费观看| 国产99白浆流出| 无人区码免费观看不卡| 中文字幕人成人乱码亚洲影| 两个人看的免费小视频| or卡值多少钱| 免费无遮挡裸体视频| 青草久久国产| 国产av在哪里看| 亚洲avbb在线观看| 亚洲av美国av| 咕卡用的链子| 麻豆av在线久日| 99久久久亚洲精品蜜臀av| 免费在线观看亚洲国产| 亚洲国产看品久久| 精品久久久久久久久久免费视频| 免费高清在线观看日韩| 9热在线视频观看99| 大香蕉久久成人网| 久久精品91蜜桃| 最近最新中文字幕大全免费视频| 国产亚洲av高清不卡| av福利片在线| 久久国产精品男人的天堂亚洲| 天堂影院成人在线观看| 国产精品日韩av在线免费观看 | 亚洲一码二码三码区别大吗| 美国免费a级毛片| 久久国产精品男人的天堂亚洲| 亚洲熟妇熟女久久| 免费在线观看日本一区| 日韩欧美三级三区| 国产精品永久免费网站| 亚洲片人在线观看| 精品日产1卡2卡| 最好的美女福利视频网| 国产欧美日韩综合在线一区二区| 欧美+亚洲+日韩+国产| 婷婷丁香在线五月| www.999成人在线观看| 搡老熟女国产l中国老女人| 亚洲五月天丁香| 国产成人欧美| 久久草成人影院| 久久精品国产亚洲av香蕉五月| 夜夜躁狠狠躁天天躁| 黄色视频,在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲一卡2卡3卡4卡5卡精品中文| 美女免费视频网站| 黄片播放在线免费| 国产精品自产拍在线观看55亚洲| 曰老女人黄片| 欧美人与性动交α欧美精品济南到| 亚洲成人精品中文字幕电影| 九色国产91popny在线| 亚洲人成77777在线视频| 人妻丰满熟妇av一区二区三区| videosex国产| 九色亚洲精品在线播放| 国产麻豆成人av免费视频| 日韩精品免费视频一区二区三区| 啦啦啦韩国在线观看视频| 亚洲av第一区精品v没综合| 激情视频va一区二区三区| 日本vs欧美在线观看视频| 日韩欧美一区二区三区在线观看| 黄色毛片三级朝国网站| 99国产精品一区二区蜜桃av| 亚洲国产欧美日韩在线播放| 中出人妻视频一区二区| 亚洲av成人av| 久久久久国产精品人妻aⅴ院| 真人一进一出gif抽搐免费| 自线自在国产av| 亚洲国产精品久久男人天堂| videosex国产| 在线观看66精品国产| 热99re8久久精品国产| 欧美黑人精品巨大| 精品高清国产在线一区| 夜夜爽天天搞| 此物有八面人人有两片| 亚洲欧洲精品一区二区精品久久久| 90打野战视频偷拍视频| 亚洲欧美日韩另类电影网站| 国产精品秋霞免费鲁丝片| 国产免费av片在线观看野外av| 中文字幕高清在线视频| 母亲3免费完整高清在线观看| 亚洲一区二区三区色噜噜| 成年女人毛片免费观看观看9| 中文字幕人妻丝袜一区二区| 少妇熟女aⅴ在线视频| 麻豆国产av国片精品| 99re在线观看精品视频| 久久 成人 亚洲| 久久久久久久久久久久大奶| 日韩精品青青久久久久久| 国产av一区在线观看免费| 90打野战视频偷拍视频| 国产精品亚洲一级av第二区| aaaaa片日本免费| 亚洲成人国产一区在线观看| 亚洲专区字幕在线| 一区在线观看完整版| 国产不卡一卡二| 免费一级毛片在线播放高清视频 | 亚洲国产精品合色在线| 国产精品美女特级片免费视频播放器 | 真人一进一出gif抽搐免费| 国产xxxxx性猛交| 日日摸夜夜添夜夜添小说| 国产精品久久久久久人妻精品电影| 女人被躁到高潮嗷嗷叫费观| √禁漫天堂资源中文www| or卡值多少钱| 日本黄色视频三级网站网址| 狂野欧美激情性xxxx| 日韩三级视频一区二区三区| 日韩精品中文字幕看吧| 国产成人av教育| 午夜免费观看网址| 国产伦一二天堂av在线观看| 如日韩欧美国产精品一区二区三区| 成人三级黄色视频| 成人欧美大片| 视频在线观看一区二区三区| 亚洲自拍偷在线| 精品国内亚洲2022精品成人| av网站免费在线观看视频| 亚洲精品在线观看二区| 国产在线精品亚洲第一网站| 男女下面插进去视频免费观看| 好男人电影高清在线观看| 久久久久精品国产欧美久久久| 久久中文字幕一级| 每晚都被弄得嗷嗷叫到高潮| 午夜久久久在线观看| 日韩视频一区二区在线观看| 国产99白浆流出| 亚洲国产中文字幕在线视频| 成在线人永久免费视频| 老汉色∧v一级毛片| 午夜成年电影在线免费观看| 国产精品自产拍在线观看55亚洲| 亚洲精品国产区一区二| 亚洲男人的天堂狠狠| 国产精品自产拍在线观看55亚洲| 不卡一级毛片| 国产精品 欧美亚洲| 久久人妻福利社区极品人妻图片| 99久久久亚洲精品蜜臀av| 亚洲第一电影网av| 黄频高清免费视频| 亚洲 国产 在线| 亚洲情色 制服丝袜| 免费人成视频x8x8入口观看| av中文乱码字幕在线| 欧美一级a爱片免费观看看 | 黄频高清免费视频| 大型黄色视频在线免费观看| 正在播放国产对白刺激| а√天堂www在线а√下载| 韩国精品一区二区三区| 50天的宝宝边吃奶边哭怎么回事|