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

    Impacts of channel morphology on residues and ecological risks of polychlorinated biphenyls in water and sediment in Chahe River

    2016-03-03 00:59:05ZhenhuaZhaoJieSunXiaokunFangLilingXiaJaviHussain
    Water Science and Engineering 2016年4期

    Zhen-hua Zhao*,Jie Sun,Xiao-kun Fang,Li-ling Xia,Javi Hussain

    aDepartment of Plant,Soil,and Microbial Sciences,Michigan State University,East Lansing,MI 48824,USA

    bKey Laboratory of Integrated Regulation and Resource Development on Shallow Lakes of Ministry of Education,Hohai University,Nanjing 210098,China

    cCollege of Environment,Hohai University,Nanjing 210098,China

    dNanjing Institute of Industry Technology,Nanjing 210016,China

    Impacts of channel morphology on residues and ecological risks of polychlorinated biphenyls in water and sediment in Chahe River

    Zhen-hua Zhaoa,b,c,*,Jie Sunb,c,Xiao-kun Fangb,c,Li-ling Xiad,Javid Hussainb,c

    aDepartment of Plant,Soil,and Microbial Sciences,Michigan State University,East Lansing,MI 48824,USA

    bKey Laboratory of Integrated Regulation and Resource Development on Shallow Lakes of Ministry of Education,Hohai University,Nanjing 210098,China

    cCollege of Environment,Hohai University,Nanjing 210098,China

    dNanjing Institute of Industry Technology,Nanjing 210016,China

    The impacts of channel morphology on the residues and ecological risks of 14 polychlorinated biphenyl(PCB)congeners in water and sediment were investigated in summer(July)and autumn(September)in the Chahe River,in Nanjing,China.The residual concentrations of trichlorobiphenyls(tri-CBs,PCB 18)and tetra-CBs(PCB 52)in water were significantly higher than those of penta-CBs to deca-CBs,and the average residual concentration of∑PCBs(sum of 14 PCB congeners)in summer was about six times higher than in autumn.However,the residues in sediment did not change significantly.Redundancy analysis(RDA)indicated that channel morphology and the corresponding environmental indices had significant impacts on PCB residues and their composition profiles in water and sediment.The overfl ow weir and lake-type watercourse may remarkably reduce the residual concentration and ecological risks of PCBs in water.The highest reduction percentages of the residualconcentration and ecological risks of∑PCBs induced by an overflow weir were 78%and 67%,respectively,and those induced by a lake-type watercourse were 36%and 70%,respectively.The watercourses with different channel morphologies were ranked by residual∑PCBs concentrations in the following descending order:the natural ecological watercourse,vertical concrete watercourse,and vegetation-type riprap watercourse.However,they were ranked by residual∑PCBs concentrations in sediment in the following descending order:the vertical concrete watercourse,vegetation-type riprap watercourse,and natural ecological watercourse.

    PCBs;Channel morphology;RDA;Ecological risk;Toxic equivalent concentration

    1.Introduction

    Polychlorinated biphenyls(PCBs)are one set of semivolatile organic compounds,which have been used widely as electric insulators in transformers,hydraulic fl uids,and paint additives since they began being produced commercially in 1929(Wang et al.,2011).A great deal of concern has been expressed overthe presence of PCBs in differentenvironmental matrices,due to their high toxicity,persistence,bioaccumulation,biomagnification,and long-range transportation (Zhang et al.,2002;Wan et al.,2005).Once released into the aquatic environment,PCBs are usually bound to suspendedparticles or absorbed by aquatic organisms.They can be biomagnifi ed to about200-70000 times along the food chain and cause potential risks to organisms and humans(Ashley et al., 2000;Fontenot et al.,2000;Pruell et al.,2000;UNEP,2004; Dodoo et al.,2012).It is now common knowledge that PCBs pose a major threat to humans and the environment even at very low concentrations(Baars et al.,2004;Dodoo et al., 2012).Exposure to PCBs can cause a series of health impacts,including reproductive disorders and immune suppression of various organisms,and even result in cancer, teratogenesis,and mutations(Van den Berg et al.,2006;Zhao et al.,2014).

    China produced a total of approximately 10000 tons of PCBs from 1965 to 1974,including 9000 tons of tri-CBs and 1000 tons of penta-CBs(Wu et al.,2011,2015).Although PCBs are notproduced any more,some of them remain in the environment.Many studies have already analyzed and reported the distribution,level,and fate of PCBs in air,water, soil(Zhang et al.,2007;Wang et al.,2008,2010),sediment, organisms,and the human body(Zhang et al.,2003;Nakata et al.,2005;Li et al.,2008;Weijs et al.,2009;Wu et al., 2009,2011;Hu et al.,2010;Bao et al.,2012;Consonni et al.,2012;Dodoo et al.,2012).The levels of PCBs in the Haihe and Huaihe rivers are 311-3110 ng/L(∑PCB12,sum of 12 PCB congeners)and 46.8-194 ng/L(∑PCB57,sum of 57 PCB congeners),respectively,which are higher than those in the Yangtze River(with a∑PCB18concentration ranging from 0.92 to 5.11 ng/L,sum of 18 PCB congeners)and the Pearl River(with a∑PCB20concentration ranging from 0.12 to 1.47 ng/L,sum of 20 PCB congeners)(Bao et al.,2012).

    In addition,some studies show that increased amounts of aquatic plants and phytoplankton productivity can reduce the concentrations of persistent organic pollutants(POPs)in aquatic ecosystems,and the species and biomass of aquatic plants and phytoplankton can affect the bioaccumulation and distribution of PCBs(Taylor et al.,1991;Cailleaud et al., 2007;Nizzetto et al.,2012;Frouin et al.,2013;Galbˊan-Malagˊon et al.,2013;Zhao et al.,2014).The physical and biogeochemical characteristics of the aquatic environment such as light,temperature,eutrophication,and nutrient stress can affect the bioaccumulation of PCBs in aquatic ecosystems (Magnusson and Tiselius,2010;Berrojalbiz et al.,2011;Zhao et al.,2014).However,few researchers have investigated the impact of channel morphology on the residual characteristics of PCBs.Recently,China has embarked on a pilot national monitoring program to assess ecological integrities of major watersheds since 2010.The components used in this monitoring program included hydrology,channel morphology, physico-chemical parameters,residues of pollutants(e.g., heavy metals,nutrients,and POPs),ecotoxicological aspects, types and numbers of biota and age,and growth of fi sh(Wang etal.,2014).The change of channelmorphology can affectthe fl ow rate,physico-chemical indices,such as pH and dissolved oxygen(DO),and the deposition of suspended particles, leading to a change in residualcharacteristics of POPs in water and sediment(Ko and Baker,2004;Wurl and Obbard,2006; Yan et al.,2008;Sandy,2010).

    The objectives of this study were to describe the impact of channel morphology on the residues and ecological risks of PCBs based on the investigation and redundancy analysis (RDA)of residual PCBs and environmentalindices atdifferent sampling sites.

    2.Materials and methods

    2.1.Study area

    The Chahe Watershed,a typical small agricultural watershed consisting mainly of small patches of uplands and paddy fields,is located in Baima Town(ata latitude of 31°34′N and a longitude of 119°10′E)in Nanjing,China(Fig.1).The study area is located in a subtropical monsoon climate zone with an average annualtemperature of 15.4°C and an annual rainfallof 1087.4 mm(Zhu etal.,2013).The totalarea is about4.087 km2and mostof the land is used for farming,of which paddy fi elds account for 50%and uplands account for 30%.The river mainly undertakes the functions of irrigation and drainage.In addition,most of the watercourses show a natural formation with small cross-sections and large curvature.For the purposes of irrigation and impoundment,a large number of sluices, dams,and other impounding buildings were built along the watercourses.Overfl ow weir 1 and Overfl ow weir 2 are located near Bridge 1 and Bridge 2,respectively.These structures affect the transportation of pollutants and sediment.There are mainly two kinds of pollution sources(mainly containing high concentrations of nitrogen,phosphorous,and trace POPs)along the river:(1)a large amount of field water,with agricultural non-point source pollutants,directly draining into the upper reaches of the Chahe River,especially at sampling site#1;and (2)a point source sewage outfall located between sampling sites#4 and#5.In addition,a paintfactory and a waste plastics recovery plant located near the Chahe River are also possible pollution sources(Bao et al.,2012).The Chahe River has been heavily polluted by the agricultural non-pointsource pollutants and sewage.In order to improve its water quality,some ecological engineering measures,such as overfl ow weirs, aquatic plants,and ecological riverbanks and watercourses, have been applied to the river since 2009.

    2.2.Sample collection

    The samples of water and sediment were collected according to the channel morphology along the Chahe River in the summer(on July 7,2013)and autumn(on September 30, 2013)(Fig.1).In July,the rainy season of this region,there is a lot of surface runoff fl owing into the river.Rice in July is at its tillering stage,and it needs plenty of pesticides and extensive irrigation.Therefore,fi eld water easily drained into the river and increased the pollution load of the river.Meanwhile,in September,rice is at its milk stage,and the climatic characteristics are entirely different from those in July.The sampling sites were set according to the channel morphology of the Chahe River from north to south:the#1,#2,#3,#4,#5, and#6 sites(Fig.1 and Table 1).Sites#1 and#2 were locatedin a watercourse with a natural ecological riverbank,and#1 was the representative site.Sites#3 and#4 were located in a watercourse with vegetation-type rock revetment,and#4 was the representative site.Sites#5 and#6 were located in an area with a vertical concrete riverbank,and#6 was the representative site.

    Fig.1.Chahe Watershed and locations of sampling sites.

    Three parallel water samples at a water depth of 0.5 m at each site were collected in brown glass bottles with a polytetrafl uoroethylene stopple.They were used for the determination of PCB concentrations(Zhao et al.,2014).These glass bottles were cleaned with K2Cr2O7-H2SO4solution and distilled water before being used.The samples were kept in a refrigerator at4°C afterthey were broughtto the laboratory.At each site,aboutten random plant samples in a 20 cm×20 cm area were collected.The fresh plant samples were washed, weighed,and used to assess the biomass of aquatic plants(Li et al.,2002).Sedimentsamples were collected with a stainless Peterson grab.After the sedimentsamples were dried in a cool and well-ventilated place,grinded,treated to remove impurities,and sieved through 100 meshes,they were mixed thoroughly for detection.Then,the samples were stored in brown glass bottles.

    2.3.Analysis of environmental parameters

    After the portable electrochemical meters were precalibrated,the temperature,DO concentration,pH,oxidationreduction potential(ORP),and conductivity were measured in situ using SX751-series portable electrochemical meters (San-Xin Instrumentation Inc.,Shanghai,China)(Zhao et al., 2014).The river flow rate was measured using a LSH10-1A portable ultrasonic wave Doppler flowmeter(Boyida Instrumentation Inc.,Xiamen,China).The turbidity was measured with a 2100Q portable turbidimeter(Hach Company,Loveland, Colorado,USA)(Zhao et al.,2013).

    2.4.Chemicals and materials

    In this study,the acetone,hexane,methanol,methylene chloride,anhydrous sodium sulfate,neutral alumina and copper powder,C18-bonded silica gel,and non-bonded silica gel were used.Their detailed information is the same as in Zhao et al.(2014).

    Fourteen PCBs were purchased from Dr.Ehrenstorfer GmbH(Augsburg,Germany),including tri-CBs(PCB 18, PCB 28,and PCB 31),tetra-CBs(PCB 44 and PCB 52),penta-CBs(PCB 101 and PCB 118),hexa-CBs(PCB 138,PCB 149, and PCB 153),hepta-CBs(PCB 170 and PCB 180),octa-CBs (PCB 194),and deca-CBs(PCB 209).These PCBs were diluted with n-hexane into standard solutions for experimental work(Zhao et al.,2014).

    2.5.Extraction of PCBs from water samples

    A 10-L water sample was enriched with a C18 solid-phase extraction column ateach site.The eluent was concentrated to approximately 1 mL with a vacuum rotary evaporator,and blown to near dry with N2.After the concentrated solution was adjusted to 1 mL with n-hexane,the PCBs were determined by a gas chromatograph/ion trap mass spectrometer(GC-IT/MS, Thermo Fisher Scientifi c Inc.,Billerica,MA,USA).The detailed extraction method of PCBs from water was described by Zhao et al.(2014).

    2.6.Extraction of PCBs from sediment samples

    5.0 g of sediment sample and 2.0 g of anhydrous sodium sulfate were weighed,put into a 50-mL glass centrifuge tube, and mixed.After the sediment was soaked in a 15-mL methylene chloride and hexane mixture(with a volumetric ratio of 2:1)for 1 h,the ultrasonic extraction lasted for 10 min.Then,the supernatantwas transferred into a pear-shaped bottle.After the processes above were repeated twice,the extraction liquid was concentrated to about 1 mL at40°C with a vacuum rotary evaporator.The solid-phase purification column(SPPC)was obtained by fi lling the upper layer with 0.5 g of copper powder and 1 g of anhydrous sodium sulfate,the middle layer with 1 g of silica gel,and the bottom layer with 1 g of alumina in a 6-mL empty column(Zhao et al.,2014).After the SPPC was rinsed and activated with 5 mL of hexane,concentrated extracting liquor was transferred into the SPPC.Then,the SPPC was eluted with a 25-mL mixture of hexane and methylene chloride(with a volumetric ratio of 3:1),and concentrated to about 1 mL.Finally,the volume was adjusted to 1 mL with n-hexane for determination.

    Table 1 Division of different channel morphologies for Chahe Watershed.

    2.7.Chromatographic conditions

    The qualitative detection of PCBs in samples was based on the chromatographic conditions described in Zacharis et al. (2012).A thermo GC-IT/MS system consisting of a TRACE GC Ultra gas chromatograph and a Polaris-Q ion-trap mass spectrometer(Thermo Fisher Scientifi c Inc.,Billerica,MA, USA)was used.The flow rate of the carrier gas(helium,with a purity of 99.999%)was 2.0 mL/min in constantflow,and the non-splitmode was selected.The column temperature and ion source temperature were 290°C and 220°C,respectively.The injection volume was 2μL.The mass spectrometer was operated in the electron impact(EI)mode at 70 eV.The detection was achieved in the selected ion monitoring(SIM) mode.Full-scan mass spectrometer(MS)data were acquired over the range of mass to charge ratios of 50-500 to obtain the fragmentation spectra of the analytes.Also,the U.S.NIST v.2.0g mass spectral database was used for the MS data interpretation.Qualitative detection of the relative retention time and peak were performed according to an internal standard(2,4,5,6-tetrachloro-m-xylene,TMX)(Zacharis et al., 2012).

    2.8.Quality control

    To verify the accuracy of the experimental data,a field blank and a laboratory blank for each batch of samples were set up.TMX was used as the internal standard to monitor the sample-preparation process.The recovery of the internal standard was 78.6%-104.4%.Each sample was analyzed six times,the relative standard deviation(RSD)was 2.1%-5.3%, and the standard error was between 0.14%and 13.52%(Zhao et al.,2014).Calibration was performed by extracting the spiked water-methanol samples at seven calibration levels (0.1,0.5,2.5,5,10,50,and 100 ng/L).The detection limitwas defi ned as the concentration,giving a signal-to-noise ratio of 3.Linearity was investigated within the range between the detection limit and 100 ng/L,and the values of the coeffi cient of determination(R2)were between 0.9950 and 0.9999.The detection limitof the target compounds(including the internal standard)was 0.25-1.0 ng/L(Popp et al.,2005).

    2.9.Data analysis

    The residual level and component characteristics of PCBs in water and sedimentwere graphed using GraphPad Prism 6.0 (Windows software package).Redundancy analysis(RDA)of residual quantities of PCBs,environmental indices,and samples were performed using the Canoco 5.0 software package. Values for residual PCBs and environmental indices were normalized by lg(x+1)transformation.Through comparisonof the RDA ordination results,we identifi ed the dominant environmental indices,which affected the residues of PCBs at different sampling sites(Zhao et al.,2014).

    3.Results and discussion

    3.1.Analysis of physicochemical indices

    Physicochemical indices of water samples showed the water environmental characteristics at each sampling site. Through comparison of physicochemical indices,we can identify the key factors affecting the residual characteristics of PCBs in the river(Table 2).

    The spatial varieties of each physicochemical index in the river were not significant within the same season.However, most of the environmental indices in summer were significantly higher than in autumn at each sampling site.The minimum DO concentration appeared at site#4 in both summer and autumn.This may be related to the poor reoxygenation ability of water due to the low velocity at this site,which was located in a typicallake-type watercourse.Along the fl ow direction of the river,the value of conductivity gradually decreased in summer.The change of turbidity is closely related to the fl ow rate of water,which implies that the fl ow rate affected the sedimentation of suspended particulate matter (SPM).The slightly acidic pH values at sites#2 and#3 in autumn may be related to the decrease of aquatic plant biomass,which leads to weakened photosynthesis and a high concentration of CO2in water.

    3.2.Residualleveland componentcharacteristics ofPCBs

    3.2.1.Residues of PCBs in water

    The residual concentrations of PCBs in water in summer and autumn at each sampling site are shown in Fig.2.The∑PCBs(sum of 14 PCB congeners)concentration in summer and autumn ranged from 5.81 to 276.82 ng/L(with an average of 108.73 ng/L)and from 6.65 to 30.60 ng/L(with an average of 18.13 ng/L),respectively.The average of∑PCBs concentration in summer was about six times higher than in autumn.The highest∑PCBs concentrations in water in both summer and autumn were found at site#5,which might be related to the presence of a point source sewage outfall upstream of this sampling site.The decrease of the∑PCBs concentration in autumn(the highest∑PCBs concentration was only 30.60 ng/L)may be due to the abundant absorption of PCBs by rice in paddy fi elds in the late growth stage of rice and aquatic plants in the river.The absorption of PCBs by rice may decrease the residualconcentration of PCBs in field water and indirectly reduce the runoff output of PCBs into the river. The absorption of PCBs by a large number of aquatic plants may also lead to a decrease of the PCB concentration in rivers (Richard et al.,1997;Huesemann et al.,2009).

    Fig.2.Distribution characteristics of PCBs in water in summer and autumn.

    As for the profile of PCB congeners,the concentrations of low-chlorinated PCBs,including tri-CBs and tetra-CBs in water(accounting for 40.47%-96.39%,with an average of 74.65%,mainly including PCB 18,PCB 44,and PCB 52), were significantly higher than those of highly chlorinated PCBs including penta-CBs to deca-CBs(accounting for 3.61%-59.53%,mainly including PCB 149,PCB 170,and PCB 194).In general,the dominant PCBs were tetra-CBs, followed by tri-CBs,penta-CBs,and hexa-CBs in different seasons at all sites except at site#3.The concentrations of other PCB congeners were much lower.

    Table 2 Physicochemical indices of samples in summer and autumn.

    3.2.2.Residues of PCBs in sediment

    The residualconcentrations of PCBs in surface sedimentin summer and autumn in the study region are shown in Fig.3. The∑PCBs concentration in sedimentin summer and autumn ranged from 12.91 to 65.11 ng/g(with an average of 38.69 ng/g)and from 4.88 to 79.41 ng/g(with an average of 32.10 ng/g),respectively.The∑PCBs concentration did not change significantly in sedimentin summeror autumn,and the principal components were similar to those in water.The highest∑PCBs concentrations in summer and autumn were detected at sites#2 and#6,respectively,which might be related to the relative deposition and accumulation of PCBs in sediments upstream of the overfl ow weir.

    The profi le of PCB congeners in sediment was similar to that in water.The sum of average percentages of tri-CBs and tetra-CBs accounted for about 87.70%of the concentration of∑PCBs(mainly including PCB 18,PCB 44,and PCB 52).In general,the dominant congeners were tri-CBs,followed by tetra-CBs,hexa-CBs,and penta-CBs in different seasons.

    The production and application of PCBs in China mainly consist of tri-CBs(90%)and penta-CBs(10%)(Que et al., 2006).However,whether in water or in sediment,the concentration of penta-CBs was always lower in this study.In sediment,the concentration of tri-CBs ranged from 1.04 to 58.81 ng/g,accounting for 2.07%-91.18%of the∑PCBs concentration,and the highest concentration of tri-CBs in summer and autumn occurred at site#6.The concentration of penta-CBs was only between 0 and 0.76 ng/g.This may be due to easy degradation dechlorination of penta-CBs in the environment(Sun et al.,2008).

    Fig.3.Distribution characteristics of PCBs in sediment in summer and autumn.

    3.3.Impact of channel morphology on PCB residues

    In the upper reach of the Chahe River(with sites#1 and#2), site#1 was considered to have the typicalchannel morphology of a long and narrow river with natural ecological riverbank. Although there were plenty of floating plants and submerged plants in the watercourse to reduce the PCB residues at site#1, the∑PCBs concentration atthis site was stillthe highestof all the sampling sites except site#5,due to the high flow rate, difficulty ofdeposition of SPM,and receiptof non-pointsource pollution from farmlands atsite#1.Hence,the pollution source and high fl ow rate played a dominant role in the increase of PCB concentrations in water.Ko and Baker(2004)reported that large amounts of suspended particles in the Susquehanna River resulted in elevated levels of hydrophobic organic compounds(HOC)and increased loads of these contaminants flowing to Chesapeake Bay during high flow periods.

    When the water flowed to site#2 in a lake-type watercourse, the presence of Overfl ow weir 1 near Bridge 1 and the broadening of river width caused the flow rate to greatly slow down and a large number of PCBs bound to SPM to be deposited easily to the bottom of the river.Therefore,the∑PCBs concentration in surface water decreased atsite#2(the mean value in summer and autumn was 46.00 ng/L).At the same time,the∑PCBs concentration in sediment increased greatly at this site(Fig.3).Hence,the overfl ow weir and low flow rate played a dominant role in the decrease of PCB concentrations in water.The reduction of turbidity also means the sedimentation of particles,and the overflow weir may play a key role in separating the water column into different layers. This phenomenon of vertical distribution profi les of PCBs in the water column has also been reported in the Black Sea and along Singapore's coast(Maldonado and Bayona,2002;Wurl and Obbard,2006).The types and sizes of SPMs in different water layers may affect the distribution and enrichment of contaminants due to the adsorption ability difference of different SPMs in relation to contaminants,and this phenomenon is usually affected by the resuspension processes of contaminated sediments and heavy rainfall events(Wurl and Obbard,2006).

    In the middle reach of the Chahe River(with sites#3 and #4),after Overfl ow weir 1,the flow rate increased sharply at site#3.The narrow watercourse and high flow rate caused the∑PCBs concentration to further decrease in water and sediment,which mainly comes from the overlying water upstream, due to the quick volatilization of low-chlorinated PCBs.Allthe possible reasons caused the∑PCBs concentration at this site to reach the lowest value(the mean value in summer and autumn was 10.24 ng/L).Although there is no evidence that high fl ow rate can accelerate the volatilization of lowchlorinated PCBs,several reports have shown that high wind speed and high temperature can accelerate the volatilization of PCBs(Jeremiason etal.,1994;Agrelletal.,2002;Dachs etal., 2002;Meijeretal.,2006,2009;Rowe and Perlinger,2012).For large water bodies,gas exchange with the atmosphere is a dominantprocess in the mass balance for relatively volatile and water-insoluble chemicals such as PCBs and mercury(Jeremiason et al.,1994;Rowe and Perlinger,2012).In the Hudson River,as well as many other water bodies,mass balances indicate that the major process for the removal of PCBs is volatilization,which accounted for about 60%of all PCB losses for congeners 2 through 6.Mono-CBs through tri-CBs congeners accounted for about half of the∑PCBs fl ux,with tetra-CBs through hexa-CBs congeners accounting forthe other half.Fluxes of congeners with more than six chlorines were not calculated.The average daily∑PCBs flux was 4μg/(m2·d), and the mass transfer velocity of PCBs was primarily under the control of water phase(Yan et al.,2008;Sandy,2010).

    Site#4 was considered to have the typical channel morphology of a lake-type watercourse with vegetation-type rock revetment.The absorption of vegetation near the watercourse and plenty of aquatic plants in the river could reduce the residual∑PCBs concentrations(Larsson et al.,2000; Berglund et al.,2001;Sobek et al.,2004;Nizzetto et al., 2012).In addition,the wider watercourse and lower fl ow rate were beneficial to the deposition of PCBs on SPM to the bottom of the river.The change of turbidity and the residual PCB concentrations in water and sediment also proved our hypothesis(Table 2 and Figs.2 and 3).Several studies have also shown that the interaction between eutrophication and uptake of POPs in aquatic biota can resultin reduced uptake of POPs in primary producers(growth dilution),and the increase of the downward fl ux of pollutants to the bottom.Therefore, lower amounts of pollutants are transported in the food chain(Larsson et al.,2000).Nizzetto et al.(2012)reported that the compounds with lg Kow>6.5(where Kowis the octanol-water partition coefficient)in the lake were mainly associated with the biomass of the primary producer,and the phytoplankton biologicalpump was a major driverof exportof POPs from the epilimnion,causing the decline of dissolvedphase concentrations.Sobek et al.(2004)suggested that partitioning of PCBs to particulate organic carbon in surface water was equilibrated and not kinetically limited by factors such as high growth rate of phytoplankton or large cell size. Berglund et al.(2001)showed that more PCBs were accumulated and buried in the sediment of eutrophic lakes than in oligotrophic lakes,where a greater fraction of∑PCBs load was dissolved in water(Berglund etal.,2001).Our research in the Qinhuai River indicated that the bioaccumulation of PCBs was significantly correlated to the composition of phytoplankton.Moreover,Chlorophyta,Bacillariophyta,and Euglenophyta had strong capacities to take up PCBs(Zhao et al.,2014).

    However,the presence of a sewage outfallbetween sites#4 and#5 and the effect of Overflow weir 2 near Bridge 2 resulted in higher∑PCBs concentration in water and sediment at site#4(the mean value in summer and autumn was 40.26 ng/L in water)than atsite#3.Moreover,this also led to the highest∑PCBs concentration at site#5 in water(the average was 153.71 ng/L).However,the∑PCBs concentration in sediment remained at a lower level,which may mean that the pollution sources contain fewer suspended particles.

    Site#6 was considered to have the typical channel morphology of a river with a vertical concrete riverbank with very few plants and poor permeability.Although the presence of Overfl ow weir 2 upstream of site#6(near Bridge 2)and the high fl ow rate can decrease the∑PCBs concentration in water due to the pre-sedimentation at site#5,the high flow rate and the aerobic water environmentnear the overfl ow weir were not conducive to the deposition and degradation of PCBs in sediment and water.Therefore,the residual∑PCBs concentration in water and sediment was higher at site#6(the mean value in water in summer and autumn was 59.02 ng/L)than at other sites.

    In general,the∑PCBs concentrations in water at different sites were ranked in the following descending order:site#5 (153.71 ng/L),site#1(71.36 ng/L),site#6(59.02 ng/L),site #2(46.00 ng/L),site#4(40.26 ng/L),and site#3(10.24 ng/L). The highest reduction percentages of the residual∑PCBs concentrations induced by the overflow weir and lake-type watercourse were 78%and 36%,respectively.However,the residual∑PCBs concentrations in sediment were ranked in the following descending order:verticalconcrete watercourse, vegetation-type riprap watercourse,and natural ecological watercourse.

    Through the analysis of the characteristics of PCB residuals,we find thatthe overflow weir may remarkably reduce the∑PCBs concentration in surface water,which implies that PCBs bound to SPM may be a dominant factor affecting the∑PCBs concentration in surface water.The sedimentation of SPM is advantageous to the reduction of the∑PCBs concentration in water due to the low flow rate upstream of the overflow weir.Downstream of the overfl ow weir,the higher volatilization rate of PCBs is also an important contribution to the decrease of the residual∑PCBs concentration in water because of the high fl ow rate.In addition,the absorption of plentiful aquatic plants and the degradation in a lake-type watercourse are also important reasons for the decrease of the residual∑PCBs concentration in water.

    3.4.Redundancy analysis between environmental indices and residual∑PCBs

    Redundancy analysis(RDA)was conducted to study the relationship between the residual∑PCBs and environmental indices.When the residual PCBs data and environmental indices data acted as the species matrix and environmental matrix,respectively,the eigenvalues of the fi rst two species axes were 0.194 and 0.165(Table 3),and the total eigenvalue was 0.875.The correlation coeffi cients between the fi rst two species axes and the first two environmental axes were 0.729 and 0.614,respectively(Zhao et al.,2014).Due to the factthat corresponding data are compositional data and have a gradient length of 1.3 standard deviation(SD)units,the linear method, RDA,is recommended.The ordination diagram of PCB residues,environmental indices,and sampling sites in water and sediment in summer and autumn are provided in Fig.4.In Fig.4,“sampling site name”and“sampling site name+*”refer to the water samples in summer and autumn,respectively (e.g.,#1 and#1*),and“sampling site name+S”and“sampling site name+*S”refer to the sediment samples in summer andautumn,respectively(e.g.,#1S and#1*S).RDA results show that the first two variable axes can explain 92.1%of the total fi tting cumulative variable,indicating that these results are credible.

    Table 3 Eigenvalues of axes and correlation coefficients of environmental axis and species axis.

    From the length of the arrows in Fig.4,we can see that the impacts of water environmentalindices including conductivity, ORP,pH,temperature,river width,flow rate,and DO concentration on the residual concentration and distribution of PCBs were greater than those of turbidity and biomass at the corresponding sampling sites in both summer and autumn.In particular,temperature and environmental media(water and sediment)caused the samples in the two seasons to be divided into Group 1,Group 2,and Group 3(Fig.4 and Table 4).Tetra-CBs and hexa-CBs in water showed a signifi cant positive correlation with the environmental indices of ORP,pH,temperature,DO,biomass,and river width in Group 1.This implies thatthe environmentalcharacteristics of high temperature, ORP,and DO concentration in summer were propitious to the generation of tetra-CBs through the degradation of highly chlorinated congeners.Therefore,Group 1 represents the effects of environmental indices on low-chlorinated PCBs in water samples.The sediment samples were divided into Group 2 and Group 3.In Group 3,the highly chlorinated congeners,including penta-CBs to deca-CBs,showed a negative correlation with environmental indices in autumn.This suggests that low temperature was not propitious to the degradation of highly chlorinated congeners,so the percentage of highly chlorinated PCBs notably increased(Laietal.,2015). In Group 2,the sedimentsamples are the summer samples,due to the degradation of highly chlorinated PCBs and the deposition of tri-CBs from the overlying water,and the tri-CBs become the dominant type of residual∑PCBs in sediment(Fig.3).

    The relationship between the samples and environmental factors indicated that the samples#1,#1*,#1S,#1*S,#3,#3*, #3S,#3*S,#6,#6*,#6S,and#6*S,taken from a long and narrow watercourse with a high fl ow rate,were symmetrically distributed near the three environmental indices of flow rate, conductivity,and turbidity.The samples#2,#2*,#2S,#2*S, #4,#4*,#4S,#4*S,#5,#5*,#5S,and#5*S,taken from a laketype watercourse with plentiful plants,were all located near the two index arrows of river width and biomass,and in the reverse direction of fl ow rate index.

    Therefore,channel morphology and environmental indices had a significant impact on PCB residues and their composition in the water.The high temperature,ORP,DO concentration,pH,fl ow rate,and narrow watercourse were propitious to the generation of low-chlorinated PCBs through the degradation of highly chlorinated congeners and their diffusion into the overlying water.The low temperature,laketype watercourse,and overflow weir accelerated the deposition of suspended particles and PCBs adsorbed by the particles into sediment(Table 4).

    3.5.Effect of channelmorphology on ecological risks of PCBs in water and sediment

    The toxic equivalency factor(TEF)method was employed to assess the effect of channel morphology on the potential ecological risks of PCBs in water of the study area.The 2,3,7,8-tetrachlorodibenzo-p-dioxin(TCDD)is the mostpotent congener within these groups of compounds,and several PCBs have been found to have toxic responses similarto those caused by TCDD(Giesy and Kannan,1998).As a result,the TEF method established by the World Health Organization(WHO) was used to evaluate the effects of these compounds on human and environmental health.The concentrations of these dioxinlike compounds can be converted into 2,3,7,8-TCDD toxic equivalent concentrations(TEQ)(Van den Berg et al.,1998; Wang et al.,2008).Then,the toxicity of the sample equals the summation of TEQ of each congener.The calculation was carried out as in Eq.(1):

    where Fiis the toxic equivalency factor of the i th PCB congener,and Ciis the content of the i th PCB congener in water(Zhao et al.,2010).

    According to relevant literature(Van den Berg et al.,1998, 2006;Yang et al.,2010),we selected seven indicative dioxinlike PCBs to calculate the TEQ values of water and sediment samples.The Fivalues were provided by the WHO:0.000002, 0.000005,0.000030,0.000030,0.000010,0.000100,and 0.000010 for PCB 28,PCB 52,PCB 101,PCB 118,PCB 153, PCB 170,and PCB 180,respectively.

    As shown in Fig.5,the average TEQ values of dioxin-like PCBs in water samples ranged from 0.04 to 0.39 pg/L,with an average of 0.22 pg/L.The largest reduction percentages of TEQ induced by the overfl ow weir and lake-type watercourse were 67%and 70%,respectively.The variation tendency of TEQ in water with different channel morphologies in summer and autumn was basically consistent.The average TEQ values firstdecreased upstream of site#4,reached the highestvalue at site#5,and then decreased again at site#6.

    Fig.4.Ordination diagram of PCB residues,environmental indices,and samples.

    Table 4 Relationship between residual PCBs and environmental indices in RDA analysis.

    The variation tendency of TEQ values in sediment was signifi cantly different from that in water,and sites#2 and#4 were the inflection points of a trend line.The overfl ow weir and lake-type watercourse play an important role in sediment precipitation and the reduction of TEQ.Based on the comprehensive analysis of the two kinds of change trends of TEQ values in water and sediment,we find that there are two kinds of sources of pollution in the river.The first pollution source,which is mainly field water,has a higher percentage of low-risk PCBs and large amounts of particulate matters,and suffers from more deposition of suspended particles(Figs.2 and 5).The second pollution source,located between sites #4 and#5,has a high content of PCB 52 and PCB 153,with larger Fivalues and less suspended particle content.Therefore, the TEQ value at site#5 reached the highest risk value (0.39 pg/L)in water and maintained a lower levelin sediment. At site#4,the TEQ value decreased mainly because of the degradation and absorption ofshoreside vegetation and aquatic plants in a lake-type watercourse on PCBs.Therefore,the lake-type watercourse is propitious to the decrease of ecological risks induced by PCBs.

    Fig.5.TEQ values in water and sediment under different channel morphologies.

    4.Conclusions

    The overfl ow weir,channelmorphology,and corresponding environmentalfactors had significanteffects on PCB residues. They affected the spatio-temporal distribution characteristics and the composition profile of PCBs in water and sediment. The high temperature,ORP,DO concentration,pH,flow rate, and narrow watercourse were propitious to the generation of low-chlorinated PCBs through the degradation of highly chlorinated congeners and their diffusion into the overlying water.The low temperature,lake-type watercourse,and overflow weir accelerated the deposition of suspended particles and PCBs adsorbed by the particles into sediment. Therefore,the residual∑PCBs concentrations in water were ranked in the following descending order:natural ecological watercourse,vertical concrete watercourse,and vegetationtype riprap watercourse.However,the residual∑PCBs concentrations in sediment were ranked in the following descending order:vertical concrete watercourse,vegetationtype riprap watercourse,and natural ecological watercourse. The highest∑PCBs concentrations in water and sediment samples appeared at sites#5 and#6,respectively.The concentration of low-chlorinated PCBs,including tri-CBs and tetra-CBs,was signifi cantly higher than the concentration of highly chlorinated PCBs including penta-CBs to deca-CBs. Moreover,the TEF method shows that the overflow weir and lake-type watercourse can remarkably reduce the residual concentration and ecological risks of PCBs.These study results can aid in the comprehensive management and pollution control of similar rivers.

    Agrell,C.,Larsson,P.,Okla,L.,Agrell,J.,2002.PCB congeners in precipitation,wash outratios and depositionalfl uxes within the Baltic Sea region, Europe.Atmos.Environ.36(2),371-383.http://dx.doi.org/10.1016/S1352-2310(01)00228-X.

    Ashley,J.T.F.,Secor,D.H.,Zlokovitz,E.,Wales,S.Q.,Baker,J.E.,2000. Linking habitat use of Hudson River striped bass to accumulation of polychlorinated biphenyl congeners.Environ.Sci.Technol.34(6), 1023-1029.http://dx.doi.org/10.1021/es990833t.

    Baars,A.J.,Bakker,M.I.,Baumann,R.A.,Boon,P.E.,Freijer,J.I., Hoogenboom,L.A.P.,Hoogerbrugge,R.,van Klaveren,J.D.,Liem,A.K.D., Traag,W.A.,et al.,2004.Dioxins,dioxin-like PCBs and non-dioxin-like PCBs in foodstuffs:Occurrence and dietary intake in The Netherlands. Toxicol.Lett.151(1),51-61.http://dx.doi.org/10.1016/j.toxlet.2004.01.028.

    Bao,L.J.,Maruya,K.A.,Snyder,S.A.,Zeng,E.Y.,2012.China's water pollution by persistent organic pollutants.Environ.Pollut.163,100-108. http://dx.doi.org/10.1016/j.envpol.2011.12.022.

    Berglund,O.,Larsson,P.,Ewald,G.,Okla,L.,2001.Influence of trophic status on PCB distribution in lake sediments and biota.Environ.Pollut. 113(2),199-210.http://dx.doi.org/10.1016/S0269-7491(00)00166-4.

    Berrojalbiz,N.,Dachs,J.,Del Vento,S.,Ojeda,M.J.,Valle,M.C.,Castro-Jimenez,J.,Mariani,G.,Wollgast,J.,Hanke,G.,2011.Persistent organic pollutants in Mediterranean seawater and processes affecting their accumulation in plankton.Environ.Sci.Technol.45(10),4315-4322.http:// dx.doi.org/10.1021/es103742w.

    Cailleaud,K.,Forget-Leray,J.,Souissi,S.,Hilde,D.,LeMenach,K., Budzinski,H.,2007.Seasonal variations of hydrophobic organic contaminant concentrations in the water-column of the Seine Estuary and their transfer to a planktonic species Eurytemora affinis(Calano?da, copepoda),Part1:PCBs and PAHs.Chemosphere 70(2),270-280.http:// dx.doi.org/10.1016/j.chemosphere.2007.05.095.

    Consonni,D.,Sindaco,R.,Bertazzi,P.A.,2012.Blood levels of dioxins,furans, dioxin-like PCBs,and TEQs in generalpopulations:A review,1989-2010. Environ.Int.44,151-162.http://dx.doi.org/10.1016/j.envint.2012.01.004.

    Dachs,J.,Lohmann,R.,Ockenden,W.A.,Mejanelle,L.,Eisenreich,S.J., Jones,K.C.,2002.Oceanic biogeochemicalcontrols on globaldynamics of persistent organic pollutants.Environ.Sci.Technol.36(20),4229-4237. http://dx.doi.org/10.1021/es025724k.

    Dodoo,D.K.,Essumang,D.K.,Jonathan,J.W.A.,Bentum,J.K.,2012.Polychlorinated biphenyls in coastaltropicalecosystems:Distribution,fate and risk assessment.Environ.Res.118,16-24.http://dx.doi.org/10.1016/ j.envres.2012.07.011.

    Fontenot,L.W.,Noblet,G.P.,Akins,J.M.,Stephens,M.D.,Cobb,G.P.,2000. Bioaccumulation of polychlorinated biphenyls in ranid frogs and northern water snakes from a hazardous waste site and a contaminated watershed. Chemosphere 40(8),803-809.http://dx.doi.org/10.1016/S0045-6535(99) 00329-X.

    Frouin,H.,Dangerfield,N.,Macdonald,R.W.,Galbraith,M.,Crewe,N., Shaw,P.,Mackas,D.,Ross,P.S.,2013.Partitioning and bioaccumulation of PCBs and PBDEs in marine plankton from the Straitof Georgia,British Columbia,Canada.Prog.Oceanogr.115,65-75.http://dx.doi.org/ 10.1016/j.pocean.2013.05.023.

    Galbˊan-Malagˊon,C.J.,Del Vento,S.,Berrojalbiz,N.,Ojeda,M.J.,Dachs,J., 2013.Polychlorinated biphenyls,hexachlorocyclohexanes and hexachlorobenzene in seawater and phytoplankton from the Southern Ocean (Weddell,South Scotia,and Bellingshausen Seas).Environ.Sci.Technol. 47(11),5578-5587.http://dx.doi.org/10.1021/es400030q.

    Giesy,J.P.,Kannan,K.,1998.Dioxin-like and non-dioxin-like toxic effects of polychlorinated biphenyls(PCBs):Implications for risk assessment. Crit.Rev.Toxicol.28(6),511-569.http://dx.doi.org/10.1080/1040844 9891344263.

    Hu,G.,Luo,X.,Li,F.,Dai,J.,Guo,J.,Chen,S.,Hong,C.,Mai,B.,Xu,M., 2010.Organochlorine compounds and polycyclic aromatic hydrocarbons in surface sedimentfrom Baiyangdian Lake,North China:Concentrations, sources profiles and potential risk.J.Environ.Sci.22(2),176-183.http:// dx.doi.org/10.1016/S1001-0742(09)60090-5.

    Huesemann,M.H.,Hausmann,T.S.,Fortman,T.J.,Thom,R.M.,Cullinan,V., 2009.In situ phytoremediation of PAH-and PCB-contaminated marine sediments with eelgrass(Zostera marina).Ecol.Eng.35(10),1395-1404. http://dx.doi.org/10.1016/j.ecoleng.2009.05.011.

    Jeremiason,J.D.,Hornbuckle,K.C.,Eisenreich,S.J.,1994.PCBs in Lake Superior,1978-1992:Decreases in water concentrations reflect loss by volatilization.Environ.Sci.Technol.28(5),903-914.http://dx.doi.org/ 10.1021/es00054a023.

    Ko,F.C.,Baker,J.E.,2004.Seasonaland annualloads of hydrophobic organic contaminants from the Susquehanna River basin to the Chesapeake Bay. Mar. Pollut. Bull. 48(9-10), 840-851. http://dx.doi.org/10.1016/ j.marpolbul.2003.10.014.

    Lai,Z.,Li,X.,Li,H.,Zhao,L.,Zeng,Y.,Wang,C.,Gao,Y.,Liu,Q.,2015.Residualdistribution and risk assessmentofpolychlorinated biphenylsin surface sediments of the Pearl River Delta,South China.Bull.Environ.Contam. Toxicol.95(1),37-44.http://dx.doi.org/10.1007/s00128-015-1563-z.

    Larsson,P.,Andersson,A.,Broman,D.,Nordback,J.,Lundberg,E.,2000. Persistent organic pollutants(POPs)in pelagic systems.Ambio 29(4-5), 202-209.http://dx.doi.org/10.1579/0044-7447-29.4.202.

    Li,W.,Huang,B.,Li,R.R.,2002.Assessing the effect of fisheries development on aquatic vegetation using GIS.Aquat.Bot.73(3),187-199.http:// dx.doi.org/10.1016/S0304-3770(02)00024-4.

    Li,X.M.,Gan,Y.P.,Yang,X.P.,Zhou,J.,Dai,J.Y.,Xu,M.Q.,2008.Human health risk of organochlorine pesticides(OCPs)and polychlorinated biphenyls(PCBs)in edible fish from Huairou Reservoir and Gaobeidian Lake in Beijing,China.Food Chem.109(2),348-354.http://dx.doi.org/ 10.1016/j.foodchem.2007.12.047.

    Magnusson,K.,Tiselius,P.,2010.The importance of uptake from food for the bioaccumulation of PCB and PBDE in the marine planktonic copepod Acartia clausi.Aquat.Toxicol.98(4),374-380.http://dx.doi.org/10.1016/ j.aquatox.2010.03.009.

    Maldonado,C.,Bayona,J.M.,2002.Organochlorine compounds in the northwestern Black Sea water:Distribution and water column process.Estuarine Coast.Shelf Sci.54(3), 527-540.http://dx.doi.org/10.1006/ecss. 2000.0672.

    Meijer,S.N.,Dachs,J.,Fernandez,P.,Camarero,L.,Catalan,J.,Del Vento,S., van Drooge,B.,Jurado,E.,Grimalt,J.O.,2006.Modelling the dynamic air-water-sediment coupled fluxes and occurrence of polychlorinated biphenyls in a high altitude lake.Environ.Pollut.140(3),546-560.http:// dx.doi.org/10.1016/j.envpol.2005.06.015.

    Meijer,S.N.,Grimalt,J.O.,Fernandez,P.,Dachs,J.,2009.Seasonalfl uxes and temperature-dependent accumulation of persistent organic pollutants in lakes:The role ofinternalbiogeochemicalcycling.Environ.Pollut.157(6), 1815-1822.http://dx.doi.org/10.1016/j.envpol.2009.01.024.

    Nakata,H.,Hirakawa,Y.,Kawazoe,M.,Nakabo,T.,Arizono,K.,Abe,S.I., Kitano,T.,Shimada,H.,Watanabe,I.,Li,W.,et al.,2005.Concentrations and compositions of organochlorine contaminants in sediments,soils, crustaceans,fishes and birds collected from Lake Tai,Hangzhou Bay and Shanghai city region,China.Environ.Pollut.133(3),415-429.http:// dx.doi.org/10.1016/j.envpol.2004.07.003.

    Nizzetto,L.,Gioia,R.,Li,J.,Borg?,K.,Pomati,F.,Bettinetti,R.,Dachs,J., Jones,K.C.,2012.Biological pump control of the fate and distribution of hydrophobic organic pollutants in water and plankton.Environ.Sci. Technol.46(6),3204-3211.http://dx.doi.org/10.1021/es204176q.

    Popp,P.,Keil,P.,Montero,L.,Ru¨ckert,M.,2005.Optimized method for the determination of 25 polychlorinated biphenyls in water samples using stir barsorptive extraction followed by thermodesorption-gas chromatography/ mass spectrometry.J.Chromatogr.A 1071(1-2),155-162.http:// dx.doi.org/10.1016/j.chroma.2005.01.066.

    Pruell,R.J.,Taplin,B.K.,McGovern,D.G.,McKinney,R.,Norton,S.B.,2000. Organic contaminant distributions in sediments,polychaetes(Nereis virens)and American lobster(Homarus americanus)from a laboratory food chain experiment.Mar.Environ.Res.49(1),19-36.http://dx.doi.org/ 10.1016/S0141-1136(99)00046-X.

    Que,M.,Wen,Q.,Liu,G.,Li,Y.,2006.Distribution and treatment process of polychlorinated biphenyls in natural waters.China Water Wastewater 22(24),10-14(in Chinese).

    Richard,M.R.,Fox,M.E.,Pick,F.R.,1997.PCB concentrations and congener composition in macrophytes and sediments in the St.Lawrence River near Cornwall,Ontario.J.Great Lakes Res.23(3),297-306.http://dx.doi.org/ 10.1016/S0380-1330(97)70913-4.

    Rowe,M.D.,Perlinger,J.A.,2012.Micrometeorological measurement of hexachlorobenzene and polychlorinated biphenyl compound air-water gas exchange in Lake Superiorand comparison to modelpredictions.Atmos.Chem. Phys.12(10),4607-4617.http://dx.doi.org/10.5194/acp-12-4607-2012.

    Sandy,A.L.,2010.The Application of a Micrometeorological Technique to Measure Air-Water Exchange of Polychlorinated Biphenyls.Ph.D. Dissertation.The State University of New Jersey,New Brunswick.

    Sobek,A.,Gustafsson,¨O.,Hajdu,S.,Larsson,U.,2004.Particle-water partitioning of PCBs in the photic zone:A 25-month study in the open Baltic Sea.Environ.Sci.Technol.38(5),1375-1382.http://dx.doi.org/10.1021/ es034447u.

    Sun,Z.,Qi,J.,Zeng,Z.,Zhang,Y.,Sun,J.,Hao,Y.,2008.Polychlorinated biphenyls in the water,soil and aquatic animals from the Jiuduansha wetland of Yangtze River Estuary.Res.Environ.Sci.21(3),92-97(in Chinese).

    Taylor,W.D.,Carey,J.H.,Lean,D.R.S.,McQueen,D.J.,1991.Organochlorine concentrations in the plankton of lakes in southern Ontario and their relationship to plankton biomass.Can.J.Fish.Aquatic Sci.48(10), 1960-1966.http://dx.doi.org/10.1139/f91-233.

    United Nations Environment Programme(UNEP),2004.PCB Management and Disposal under the Stockholm Convention on Persistent Organic Pollutants.UNEP,Geneva.

    Van den Berg,M.,Birnbaum,L.,Bosveld,A.T.C.,Bronstrom,B.,Cook,P., Feely,M.,Giesy,J.P.,Hanberg,A.,Hasegawa,R.,Kennedy,S.W.,et al., 1998.Toxic equivalency factors(TEFs)for PCBs,PCDDs,PCDFs for humans and wildlife.Environ.Health Perspect.106,775-792.http:// dx.doi.org/10.2307/3434121.

    Van den Berg,M.,Birnbaum,L.S.,Denison,M.,De Vito,M.,Farland,W., Feeley,M.,Fiedler,H.,Hakansson,H.,Hanberg,A.,Haws,L.,etal.,2006. The 2005 World Health Organization reevaluation ofhuman and Mammalian toxic equivalency factors for dioxins and dioxin-like compounds.Toxicol. Sci.93(2),223-241.http://dx.doi.org/10.1093/toxsci/kfl055.

    Wan,Y.,Hu,J.,Liu,J.,An,W.,Tao,S.,Jia,Z.,2005.Fate of DDT-related compounds in Bohai Bay and its adjacent Haihe Basin,North China. Mar.Pollut.Bull.50(4),439-445.http://dx.doi.org/10.1016/j.marpolbul. 2004.11.037.

    Wang,D.G.,Yang,M.,Jia,H.L.,Zhou,L.,Li,Y.F.,2008.Levels,distributions and profiles of polychlorinated biphenyls in surface soils of Dalian,China. Chemosphere 73(1),38-42.http://dx.doi.org/10.1016/j.chemosphere. 2008.05.055.

    Wang,H.,An,Q.,Dong,Y.H.,Li,D.C.,Velde,B.,2010.Contamination and congener profiles of polychlorinated biphenyls from different agricultural top soils in a county of the Tai Lake Region,China.J.Hazard.Mater. 176(1-3),1027-1031.http://dx.doi.org/10.1016/j.jhazmat.2009.11.143.

    Wang,H.,Du,J.,Leung,H.,Leung,A.O.W.,Liang,P.,Giesy,J.P., Wong,C.K.C.,Wong,M.,2011.Distribution and source apportionments of polychlorinated biphenyls(PCBs)in mariculture sediments from the Pearl River Delta,South China.Mar.Pollut.Bull.63(5-12),516-522.http:// dx.doi.org/10.1016/j.marpolbul.2011.02.009.

    Wang,Y.,Teng,E.,Liu,T.,Lv,Y.,Jin,X.,Giesy,J.P.,Hollert,H.,2014.A nationalpilotscheme formonitoring and assessmentof ecologicalintegrity of surface waters in China.Environ.Dev.10,104-107.http://dx.doi.org/ 10.1016/j.envdev.2014.02.003.

    Weijs,L.,Dirtu,A.C.,Das,K.,Gheorghe,A.,Reijnders,P.J.H.,Neels,H., Blust,R.,Covaci,A.,2009.Inter-species differences for polychlorinated biphenyls and polybrominated diphenyl ethers in marine top predators from the Southern North Sea:Part2.Biomagnification in harbourseals and harbour porpoises.Environ.Pollut.157(2),445-451.http://dx.doi.org/ 10.1016/j.envpol.2008.09.025.

    Wu,C.,Zhu,H.,Luo,Y.,Teng,Y.,Song,J.,Chen,M.,2015.Levels and potential health hazards of PCBs in shallow groundwater of an e-waste recycling area,China.Environ.Earth Sci.74(5),4431-4438.http:// dx.doi.org/10.1007/s12665-015-4427-2.

    Wu,J.,Luo,X.,Zhang,Y.,Yu,M.,Chen,S.,Mai,B.,Yang,Z.,2009.Biomagnification of polybrominated diphenyl ethers(PBDEs)and polychlorinated biphenyls in a highly contaminated freshwater food web from South China.Environ.Pollut.157(3),904-909.http://dx.doi.org/10.1016/ j.envpol.2008.11.001.

    Wu,S.,Xia,X.,Yang,L.,Liu,H.,2011.Distribution,source and risk assessmentof polychlorinated biphenyls(PCBs)in urban soils of Beijing, China. Chemosphere 82(5), 732-738. http://dx.doi.org/10.1016/j. chemosphere.2010.10.090.

    Wurl,O.,Obbard,J.P.,2006.Distribution of organochlorine compounds in the sea-surface microlayer,water column and sedimentof Singapore's coastal environment.Chemosphere 62(7),1105-1115.http://dx.doi.org/10.1016/ j.chemosphere.2005.06.002.

    Yan,S.,Rodenburg,L.A.,Dachs,J.,Eisenreich,S.J.,2008.Seasonal airwater exchange fluxes of polychlorinated biphenyls in the Hudson River Estuary.Environ.Pollut.152(2),443-451.http://dx.doi.org/ 10.1016/j.envpol.2007.06.074.

    Yang,J.M.,Salmon,A.G.,Marty,M.A.,2010.Developmentof TEFs for PCB congeners by using an alternative biomarker:Thyroid hormone levels. Regul.Toxicol.Pharmacol.56(2),225-236.http://dx.doi.org/10.1016/ j.yrtph.2009.12.008.

    Zacharis,C.K.,Rotsias,I.,Zachariadis,P.G.,Zotos,A.,2012.Dispersive liquid-liquid microextraction for the determination of organochlorine pesticides residues in honey by gas chromatography-electron capture and ion trap mass spectrometric detection.Food Chem.134(3),1665-1672. http://dx.doi.org/10.1016/j.foodchem.2012.03.073.

    Zhang,G.,Parker,A.,House,A.,Mai,B.,Li,X.,Kang,Y.,Wang,Z.,2002. Sedimentary records of DDT and HCH in the Pearl River delta,South China.Environ.Sci.Technol.36(17),3671-3677.http://dx.doi.org/ 10.1021/es0102888.

    Zhang,J.Y.,Qiu,L.M.,He,J.,Liao,Y.,Luo,Y.M.,2007.Occurrence and congeners specific of polychlorinated biphenyls in agricultural soils from Southern Jiangsu,China.J.Environ.Sci.19(3),338-342.http:// dx.doi.org/10.1016/S1001-0742(07)60055-2.

    Zhang,Z.L.,Hong,H.S.,Zhou,J.L.,Huang,J.,Yu,G.,2003.Fate and assessment of persistent organic pollutants in water and sediment from Minjiang River Estuary, Southeast China. Chemosphere 52(9), 1423-1430.http://dx.doi.org/10.1016/S0045-6535(03)00478-8.

    Zhao,Z.,Xia,L.,Liu,J.,Jiang,Y.,Liu,Y.,2010.PCBs residues characteristics and ecological risk assessment in a Chinese paddy fields.In:IEEE the International Workshop on Construction,Environment,Transportation (CET2010),pp.1-4.http://dx.doi.org/10.1109/ICEEE.2010.5660346.

    Zhao,Z.,Mi,T.,Xia,L.,Yan,W.,Jiang,Y.,Gao,Y.,2013.Understanding the patterns and mechanisms of urban water ecosystem degradation:Phytoplankton community structure and water quality in the Qinhuai River, Nanjing City,China.Environ.Sci.Pollut.Res.20(7),5003-5012.http:// dx.doi.org/10.1007/s11356-013-1464-9.

    Zhao,Z.,Jiang,Y.,Xia,L.,Mi,T.,Yan,W.,Gao,Y.,Jiang,X.,Fawundu,E., Hussain,J.,2014.Application of canonical correspondence analysis to determine the ecological contribution of phytoplankton to PCBs bioaccumulation in Qinhuai River,Nanjing,China.Environ.Sci.Pollut.Res. 21(4),3091-3103.http://dx.doi.org/10.1007/s11356-013-2265-x.

    Zhu,L.,Liu,D.,Zhu,C.,2013.Experimental study and numerical modeling of surface/subsurface flow at field scale.J.Softw.8(3),680-686.http:// dx.doi.org/10.4304/jsw.8.3.680-686.

    Received 1 August 2016;accepted 3 October 2016

    Available online 7 January 2017

    This work was supported by the National Natural Science Foundation of China(Grants No.41371307 and 51509129),the Open Foundation of State Key Laboratory of Pollution Control and Resource Reuse(Grant No. PCRRF12010),the State Key Laboratory of Soil and Sustainable Agriculture (Institute of Soil Science,Chinese Academy of Sciences)foundation(Grant No.0812201228),a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions,and the Top-notch Academic Programs Project(TAPP)of Jiangsu Higher Education Institutions.

    *Corresponding author.

    E-mail address:zzh4000@126.com(Zhen-hua Zhao).

    Peer review under responsibility of Hohai University.

    ?2016 Hohai University.Production and hosting by Elsevier B.V.This is an open access article under the CC BY-NC-ND license(http:// creativecommons.org/licenses/by-nc-nd/4.0/).

    午夜福利18| 啪啪无遮挡十八禁网站| 美女午夜性视频免费| 成年版毛片免费区| 国产片内射在线| 禁无遮挡网站| 色综合站精品国产| 国产一区在线观看成人免费| 国产精品1区2区在线观看.| 我的亚洲天堂| 可以在线观看的亚洲视频| 听说在线观看完整版免费高清| 国产真人三级小视频在线观看| 亚洲欧美精品综合一区二区三区| 精品高清国产在线一区| 狂野欧美激情性xxxx| 国产成人精品久久二区二区91| 欧美午夜高清在线| 亚洲 欧美 日韩 在线 免费| 精品久久久久久久人妻蜜臀av| 精品久久蜜臀av无| 久久香蕉激情| 精品免费久久久久久久清纯| 久久久久久久久久黄片| 亚洲av成人不卡在线观看播放网| 国产精品1区2区在线观看.| 免费一级毛片在线播放高清视频| 久久久久国产一级毛片高清牌| 在线观看午夜福利视频| 黄色视频,在线免费观看| 婷婷六月久久综合丁香| 一进一出抽搐gif免费好疼| 免费在线观看视频国产中文字幕亚洲| 正在播放国产对白刺激| 中文字幕精品亚洲无线码一区 | 欧美黑人巨大hd| 黄色视频不卡| av在线播放免费不卡| 黑人巨大精品欧美一区二区mp4| 精品一区二区三区av网在线观看| 国产精品 国内视频| 日韩精品青青久久久久久| 一本综合久久免费| 久久久久精品国产欧美久久久| 黄色片一级片一级黄色片| 一级片免费观看大全| 视频区欧美日本亚洲| 美女午夜性视频免费| 老熟妇仑乱视频hdxx| 国产三级在线视频| avwww免费| 免费在线观看影片大全网站| 欧美日韩精品网址| 亚洲中文日韩欧美视频| 午夜久久久久精精品| 99热6这里只有精品| 亚洲成av人片免费观看| 日本a在线网址| 女生性感内裤真人,穿戴方法视频| 中文字幕人成人乱码亚洲影| 亚洲精品一区av在线观看| 亚洲电影在线观看av| av电影中文网址| 欧美最黄视频在线播放免费| 亚洲avbb在线观看| 亚洲中文av在线| www日本黄色视频网| 亚洲精品一区av在线观看| 亚洲成人久久性| 老司机在亚洲福利影院| 热re99久久国产66热| 正在播放国产对白刺激| 亚洲精品国产一区二区精华液| x7x7x7水蜜桃| 欧美成狂野欧美在线观看| 日本a在线网址| 日韩欧美免费精品| АⅤ资源中文在线天堂| 欧美黑人精品巨大| 亚洲精品一区av在线观看| 亚洲 国产 在线| 黄网站色视频无遮挡免费观看| 欧美绝顶高潮抽搐喷水| 国内精品久久久久精免费| 成人国语在线视频| 国产精品久久视频播放| 欧美性猛交黑人性爽| 18禁国产床啪视频网站| 男女做爰动态图高潮gif福利片| 午夜影院日韩av| 精品电影一区二区在线| 夜夜看夜夜爽夜夜摸| 国产乱人伦免费视频| 一本综合久久免费| 一级毛片高清免费大全| 国产成人av激情在线播放| 国产成人啪精品午夜网站| 国产精品亚洲av一区麻豆| 久9热在线精品视频| 日本 欧美在线| 51午夜福利影视在线观看| 亚洲精品中文字幕一二三四区| 天天添夜夜摸| 天堂√8在线中文| 黄片小视频在线播放| 国产精品免费一区二区三区在线| 久久狼人影院| 国产主播在线观看一区二区| 天天躁狠狠躁夜夜躁狠狠躁| www.自偷自拍.com| 久久精品91蜜桃| 国语自产精品视频在线第100页| 超碰成人久久| 香蕉久久夜色| 国产亚洲精品一区二区www| 久久久久久久久中文| 黄色毛片三级朝国网站| 我的亚洲天堂| 免费在线观看日本一区| 一级毛片女人18水好多| 50天的宝宝边吃奶边哭怎么回事| 亚洲第一电影网av| 精品久久久久久成人av| 亚洲欧美精品综合久久99| 国产精品亚洲av一区麻豆| 亚洲一区高清亚洲精品| 很黄的视频免费| 老司机午夜十八禁免费视频| 老司机午夜十八禁免费视频| 无人区码免费观看不卡| 亚洲人成电影免费在线| av福利片在线| 亚洲色图 男人天堂 中文字幕| 嫁个100分男人电影在线观看| 色精品久久人妻99蜜桃| 欧美黑人巨大hd| 视频在线观看一区二区三区| 亚洲黑人精品在线| 在线观看午夜福利视频| 免费看十八禁软件| 午夜免费激情av| 一区二区三区精品91| 黄片大片在线免费观看| 国产精华一区二区三区| 久久久久九九精品影院| 女人被狂操c到高潮| 脱女人内裤的视频| 黄色a级毛片大全视频| 国产精品一区二区免费欧美| 男女那种视频在线观看| 人妻丰满熟妇av一区二区三区| 久久热在线av| 精品久久久久久久久久免费视频| 黑人欧美特级aaaaaa片| 精品久久久久久久久久免费视频| 欧美最黄视频在线播放免费| 啦啦啦韩国在线观看视频| 国产99久久九九免费精品| 黄色成人免费大全| 亚洲国产欧美日韩在线播放| 亚洲国产欧美日韩在线播放| 嫩草影视91久久| 亚洲全国av大片| 日韩欧美国产一区二区入口| 精品国产亚洲在线| 国产亚洲av高清不卡| 成人特级黄色片久久久久久久| 一级黄色大片毛片| 亚洲一区二区三区色噜噜| ponron亚洲| 国产精品综合久久久久久久免费| 女人高潮潮喷娇喘18禁视频| 欧美在线黄色| 婷婷六月久久综合丁香| 久久久水蜜桃国产精品网| 亚洲无线在线观看| 精品国产乱码久久久久久男人| 日本a在线网址| 高清毛片免费观看视频网站| 操出白浆在线播放| 十分钟在线观看高清视频www| 欧美日韩亚洲国产一区二区在线观看| 一本精品99久久精品77| 欧美大码av| 夜夜躁狠狠躁天天躁| 人妻丰满熟妇av一区二区三区| 国产亚洲精品综合一区在线观看 | 国产激情久久老熟女| 中文资源天堂在线| 午夜视频精品福利| 亚洲色图av天堂| 午夜福利高清视频| 黄色女人牲交| 天堂动漫精品| 欧美成人性av电影在线观看| 色精品久久人妻99蜜桃| 中文字幕精品亚洲无线码一区 | 中文字幕久久专区| 制服人妻中文乱码| 女人被狂操c到高潮| 亚洲色图 男人天堂 中文字幕| 国产精品久久久久久人妻精品电影| 午夜精品久久久久久毛片777| xxxwww97欧美| 精品久久久久久久久久免费视频| 亚洲av片天天在线观看| 麻豆成人av在线观看| 男人舔奶头视频| 国产精品,欧美在线| a级毛片a级免费在线| 亚洲精品久久国产高清桃花| 免费在线观看黄色视频的| 精品久久久久久久久久免费视频| 熟妇人妻久久中文字幕3abv| 精品国产亚洲在线| av在线播放免费不卡| 亚洲国产欧洲综合997久久, | 高清在线国产一区| 好看av亚洲va欧美ⅴa在| 夜夜看夜夜爽夜夜摸| 亚洲国产精品sss在线观看| 国产精品一区二区三区四区久久 | 中文字幕精品亚洲无线码一区 | 亚洲最大成人中文| 动漫黄色视频在线观看| 国产精品香港三级国产av潘金莲| 国产又爽黄色视频| 亚洲国产精品sss在线观看| 人人妻人人澡欧美一区二区| 草草在线视频免费看| 18禁黄网站禁片免费观看直播| 精品国产一区二区三区四区第35| 日韩 欧美 亚洲 中文字幕| 亚洲国产精品999在线| videosex国产| 亚洲aⅴ乱码一区二区在线播放 | 两性夫妻黄色片| 色av中文字幕| 色综合婷婷激情| 久久精品91蜜桃| 男人舔女人下体高潮全视频| 淫秽高清视频在线观看| 亚洲国产精品成人综合色| 国产色视频综合| 午夜福利视频1000在线观看| 美女扒开内裤让男人捅视频| 婷婷丁香在线五月| 黑人操中国人逼视频| 婷婷亚洲欧美| 老司机在亚洲福利影院| 国产精品影院久久| 91麻豆精品激情在线观看国产| 日韩视频一区二区在线观看| 男女视频在线观看网站免费 | а√天堂www在线а√下载| 两性夫妻黄色片| 久久久久国产一级毛片高清牌| 国产成人av激情在线播放| 又黄又粗又硬又大视频| 久久精品人妻少妇| 国产精品影院久久| 男女床上黄色一级片免费看| 亚洲七黄色美女视频| 非洲黑人性xxxx精品又粗又长| 久久国产精品人妻蜜桃| 欧美国产精品va在线观看不卡| 成人手机av| 精品午夜福利视频在线观看一区| 国产精品永久免费网站| 又黄又粗又硬又大视频| 91字幕亚洲| 国产精品二区激情视频| 十八禁人妻一区二区| 欧美中文综合在线视频| 怎么达到女性高潮| 国产黄片美女视频| 午夜福利在线在线| 99热只有精品国产| 亚洲国产欧美一区二区综合| 欧美成人免费av一区二区三区| 成人精品一区二区免费| 在线观看www视频免费| 亚洲片人在线观看| 国产高清有码在线观看视频 | 久久精品aⅴ一区二区三区四区| 久久中文字幕一级| 久99久视频精品免费| 国产aⅴ精品一区二区三区波| 国产黄片美女视频| 午夜免费观看网址| 日韩三级视频一区二区三区| 色精品久久人妻99蜜桃| 制服诱惑二区| 国产精品亚洲av一区麻豆| 成人三级做爰电影| xxx96com| 色老头精品视频在线观看| 久久香蕉精品热| 欧美av亚洲av综合av国产av| 黄色视频,在线免费观看| 久久香蕉激情| 亚洲欧美日韩高清在线视频| 午夜两性在线视频| 久久天堂一区二区三区四区| 99久久无色码亚洲精品果冻| 亚洲欧美精品综合久久99| 在线观看舔阴道视频| 此物有八面人人有两片| 美女高潮喷水抽搐中文字幕| 18美女黄网站色大片免费观看| 精品无人区乱码1区二区| 成人欧美大片| 午夜视频精品福利| 亚洲国产看品久久| 欧美激情高清一区二区三区| 欧美日韩一级在线毛片| 观看免费一级毛片| 手机成人av网站| 日韩成人在线观看一区二区三区| www.熟女人妻精品国产| 欧美日本视频| 丁香欧美五月| 性色av乱码一区二区三区2| 麻豆成人午夜福利视频| 国产精品一区二区免费欧美| 欧美精品亚洲一区二区| 国产不卡一卡二| 满18在线观看网站| 一卡2卡三卡四卡精品乱码亚洲| 两性午夜刺激爽爽歪歪视频在线观看 | 一个人免费在线观看的高清视频| 久99久视频精品免费| 久久青草综合色| 国内少妇人妻偷人精品xxx网站 | 亚洲欧洲精品一区二区精品久久久| 一本大道久久a久久精品| 国产av在哪里看| 欧美激情高清一区二区三区| 久久久久久亚洲精品国产蜜桃av| 韩国av一区二区三区四区| 欧美日本亚洲视频在线播放| 黑丝袜美女国产一区| 嫩草影视91久久| 亚洲av美国av| 搡老妇女老女人老熟妇| 亚洲免费av在线视频| 亚洲天堂国产精品一区在线| 欧美日韩精品网址| 欧美一区二区精品小视频在线| 丰满的人妻完整版| 极品教师在线免费播放| 1024手机看黄色片| 日韩大尺度精品在线看网址| 午夜福利高清视频| 18禁美女被吸乳视频| 午夜久久久在线观看| 久热爱精品视频在线9| 国产av一区二区精品久久| 最近最新免费中文字幕在线| 最新美女视频免费是黄的| 精品一区二区三区视频在线观看免费| 亚洲av美国av| 别揉我奶头~嗯~啊~动态视频| 免费一级毛片在线播放高清视频| 99国产极品粉嫩在线观看| 免费观看人在逋| 婷婷丁香在线五月| 婷婷亚洲欧美| 久久天堂一区二区三区四区| 69av精品久久久久久| www国产在线视频色| 亚洲最大成人中文| 亚洲第一av免费看| 母亲3免费完整高清在线观看| 午夜成年电影在线免费观看| 美女国产高潮福利片在线看| 精品国产一区二区三区四区第35| 国产国语露脸激情在线看| 精品久久久久久成人av| 婷婷精品国产亚洲av| 少妇被粗大的猛进出69影院| 十八禁人妻一区二区| 高潮久久久久久久久久久不卡| 天天躁夜夜躁狠狠躁躁| 久久久久久久午夜电影| aaaaa片日本免费| 国产黄片美女视频| 88av欧美| 精品一区二区三区四区五区乱码| 亚洲专区字幕在线| 午夜a级毛片| 免费看十八禁软件| 脱女人内裤的视频| 91字幕亚洲| 亚洲七黄色美女视频| 国产久久久一区二区三区| 在线观看www视频免费| 久久久久久大精品| 精品午夜福利视频在线观看一区| 日本黄色视频三级网站网址| 嫩草影院精品99| 亚洲欧美精品综合久久99| 亚洲一区中文字幕在线| 一进一出好大好爽视频| 欧美乱妇无乱码| 国产成人精品久久二区二区免费| 美女大奶头视频| 亚洲精品久久国产高清桃花| 欧美最黄视频在线播放免费| 午夜福利在线在线| 啪啪无遮挡十八禁网站| 我的亚洲天堂| 婷婷精品国产亚洲av在线| а√天堂www在线а√下载| 久久国产精品人妻蜜桃| 精品国产乱码久久久久久男人| 丰满人妻熟妇乱又伦精品不卡| 一级毛片精品| 久久人人精品亚洲av| videosex国产| 成人18禁高潮啪啪吃奶动态图| 日韩 欧美 亚洲 中文字幕| 三级毛片av免费| 男人舔女人的私密视频| 99在线视频只有这里精品首页| 免费无遮挡裸体视频| 欧美日韩亚洲综合一区二区三区_| 怎么达到女性高潮| av有码第一页| 9191精品国产免费久久| 精品免费久久久久久久清纯| 欧美国产精品va在线观看不卡| 亚洲aⅴ乱码一区二区在线播放 | 天堂动漫精品| 美女扒开内裤让男人捅视频| 91字幕亚洲| 男人舔奶头视频| 免费在线观看完整版高清| 国产97色在线日韩免费| av中文乱码字幕在线| 女警被强在线播放| 国产男靠女视频免费网站| 日韩一卡2卡3卡4卡2021年| 无遮挡黄片免费观看| 麻豆国产av国片精品| 日韩国内少妇激情av| 叶爱在线成人免费视频播放| 黄色a级毛片大全视频| 午夜视频精品福利| 91字幕亚洲| 视频在线观看一区二区三区| 人人妻人人看人人澡| 亚洲国产看品久久| 国产麻豆成人av免费视频| 国产男靠女视频免费网站| 婷婷亚洲欧美| 波多野结衣高清作品| 欧美色欧美亚洲另类二区| 国产99白浆流出| 国产成人精品久久二区二区91| 黑人操中国人逼视频| 久久久久久久久久黄片| 99久久综合精品五月天人人| 少妇粗大呻吟视频| 亚洲av第一区精品v没综合| 亚洲av中文字字幕乱码综合 | 美女 人体艺术 gogo| 又紧又爽又黄一区二区| 国产成人精品久久二区二区91| 亚洲专区国产一区二区| 1024手机看黄色片| 日本免费一区二区三区高清不卡| 夜夜夜夜夜久久久久| 美女大奶头视频| www.熟女人妻精品国产| 黄色视频不卡| 欧美 亚洲 国产 日韩一| 国产精品 国内视频| 亚洲精品在线美女| 非洲黑人性xxxx精品又粗又长| 国产在线观看jvid| 国产aⅴ精品一区二区三区波| 国产午夜精品久久久久久| 香蕉久久夜色| 视频在线观看一区二区三区| 91字幕亚洲| 亚洲av电影不卡..在线观看| 88av欧美| 香蕉久久夜色| 久久久久九九精品影院| 久久欧美精品欧美久久欧美| 日本撒尿小便嘘嘘汇集6| 国产又色又爽无遮挡免费看| 制服人妻中文乱码| 91成年电影在线观看| 757午夜福利合集在线观看| 免费看十八禁软件| 好男人电影高清在线观看| 中出人妻视频一区二区| 91麻豆精品激情在线观看国产| 99精品在免费线老司机午夜| 黄频高清免费视频| 黄网站色视频无遮挡免费观看| 香蕉国产在线看| 精品人妻1区二区| 国产成人啪精品午夜网站| 亚洲熟妇中文字幕五十中出| 国产亚洲精品综合一区在线观看 | 欧美一级毛片孕妇| 亚洲男人天堂网一区| 9191精品国产免费久久| 丝袜人妻中文字幕| 韩国精品一区二区三区| www.精华液| 老汉色av国产亚洲站长工具| 一本大道久久a久久精品| 亚洲精品久久国产高清桃花| 丝袜在线中文字幕| 老司机午夜福利在线观看视频| 久久中文看片网| 午夜日韩欧美国产| 一级片免费观看大全| 99久久久亚洲精品蜜臀av| 国产成人欧美在线观看| 91麻豆精品激情在线观看国产| 亚洲av成人一区二区三| 国产精品久久视频播放| 午夜免费成人在线视频| 精品久久久久久久久久久久久 | 亚洲人成网站高清观看| 叶爱在线成人免费视频播放| 精品国产超薄肉色丝袜足j| 亚洲成a人片在线一区二区| 少妇被粗大的猛进出69影院| 一区二区三区国产精品乱码| 一个人观看的视频www高清免费观看 | 淫秽高清视频在线观看| 一级黄色大片毛片| 桃色一区二区三区在线观看| 亚洲精品色激情综合| 国产伦在线观看视频一区| 国产色视频综合| 国产在线观看jvid| 久久青草综合色| 狠狠狠狠99中文字幕| 亚洲一区二区三区色噜噜| 桃色一区二区三区在线观看| 欧美黑人欧美精品刺激| 熟妇人妻久久中文字幕3abv| 香蕉国产在线看| 99re在线观看精品视频| 国产成人av教育| 男人舔女人下体高潮全视频| 窝窝影院91人妻| 亚洲aⅴ乱码一区二区在线播放 | 亚洲性夜色夜夜综合| 日本成人三级电影网站| 无人区码免费观看不卡| 久久99热这里只有精品18| 成人精品一区二区免费| 色精品久久人妻99蜜桃| 国产三级在线视频| 黄片播放在线免费| 国产黄色小视频在线观看| 妹子高潮喷水视频| 亚洲欧美精品综合久久99| 久久精品影院6| 色老头精品视频在线观看| 国产国语露脸激情在线看| 久久精品国产亚洲av香蕉五月| 一级毛片精品| 在线观看午夜福利视频| 2021天堂中文幕一二区在线观 | 俄罗斯特黄特色一大片| 一区二区三区激情视频| 久久久国产成人精品二区| 日本精品一区二区三区蜜桃| 色婷婷久久久亚洲欧美| 99精品久久久久人妻精品| 精品国产一区二区三区四区第35| 日日爽夜夜爽网站| 亚洲五月天丁香| av天堂在线播放| 久久热在线av| 精品免费久久久久久久清纯| 日韩免费av在线播放| svipshipincom国产片| 亚洲狠狠婷婷综合久久图片| 午夜两性在线视频| 可以免费在线观看a视频的电影网站| 亚洲精品粉嫩美女一区| 国产精品1区2区在线观看.| 亚洲欧美激情综合另类| 日韩欧美一区二区三区在线观看| 一区二区三区国产精品乱码| 一区二区三区精品91| 97人妻精品一区二区三区麻豆 | 精品乱码久久久久久99久播| 国产免费男女视频| 午夜精品久久久久久毛片777| 国产主播在线观看一区二区| av电影中文网址| 麻豆国产av国片精品| 日韩欧美一区二区三区在线观看| 免费在线观看视频国产中文字幕亚洲| 午夜精品久久久久久毛片777| 一进一出好大好爽视频| 又黄又粗又硬又大视频| 中文字幕人妻熟女乱码| 色综合亚洲欧美另类图片| 亚洲专区字幕在线| 在线观看日韩欧美| 国产蜜桃级精品一区二区三区| 搞女人的毛片| 免费在线观看黄色视频的|