• <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/).

    国产精品一二三区在线看| a级毛片免费高清观看在线播放| 久久精品熟女亚洲av麻豆精品| av国产久精品久网站免费入址| 一区二区av电影网| 大香蕉久久网| 涩涩av久久男人的天堂| 亚洲四区av| 一二三四中文在线观看免费高清| 少妇裸体淫交视频免费看高清| 久久久久人妻精品一区果冻| 精品久久久久久久久av| 99久久中文字幕三级久久日本| 日本色播在线视频| 婷婷色综合大香蕉| 精品少妇黑人巨大在线播放| 性色av一级| 久久久久精品性色| 在线 av 中文字幕| 亚洲四区av| av国产久精品久网站免费入址| 亚洲自偷自拍三级| av播播在线观看一区| 人人妻人人澡人人爽人人夜夜| 亚洲三级黄色毛片| 亚洲精品国产av蜜桃| 51国产日韩欧美| 麻豆成人午夜福利视频| 日日撸夜夜添| 国产高清不卡午夜福利| 亚洲欧美一区二区三区国产| 中文欧美无线码| 国产一区二区三区综合在线观看 | 成人漫画全彩无遮挡| 在线观看一区二区三区| 亚洲精品,欧美精品| 街头女战士在线观看网站| 亚洲精品日韩在线中文字幕| 日韩av在线免费看完整版不卡| 男女边摸边吃奶| 精品少妇黑人巨大在线播放| 国产一级毛片在线| 熟女人妻精品中文字幕| 久久久久网色| 男人和女人高潮做爰伦理| 91精品国产国语对白视频| 又粗又硬又长又爽又黄的视频| 亚洲综合色惰| 久久精品人妻少妇| 黑丝袜美女国产一区| 国产有黄有色有爽视频| 精品熟女少妇av免费看| 人体艺术视频欧美日本| 日韩国内少妇激情av| 成人二区视频| 一区二区三区乱码不卡18| 免费av中文字幕在线| 3wmmmm亚洲av在线观看| 交换朋友夫妻互换小说| 国产爽快片一区二区三区| 日韩 亚洲 欧美在线| 99热6这里只有精品| 人妻系列 视频| 99热6这里只有精品| 青春草国产在线视频| 中文字幕免费在线视频6| 国产日韩欧美亚洲二区| 黄色欧美视频在线观看| 国产无遮挡羞羞视频在线观看| 亚洲av电影在线观看一区二区三区| 日韩三级伦理在线观看| xxx大片免费视频| 大片电影免费在线观看免费| 99久久中文字幕三级久久日本| 国产亚洲5aaaaa淫片| 只有这里有精品99| 中文天堂在线官网| 看十八女毛片水多多多| 久久久久久久久久成人| 伊人久久精品亚洲午夜| 亚洲中文av在线| 国产成人精品福利久久| 亚洲精品国产av蜜桃| a 毛片基地| 视频中文字幕在线观看| 高清av免费在线| 一区二区av电影网| 哪个播放器可以免费观看大片| 国产又色又爽无遮挡免| 日本黄色片子视频| 一区二区三区免费毛片| 青春草亚洲视频在线观看| 日韩三级伦理在线观看| 日本色播在线视频| 我的老师免费观看完整版| 久久人妻熟女aⅴ| 高清在线视频一区二区三区| 亚洲成人一二三区av| 欧美高清性xxxxhd video| 亚洲国产精品国产精品| 成人美女网站在线观看视频| 亚洲精品一二三| 一个人看视频在线观看www免费| 99国产精品免费福利视频| a级一级毛片免费在线观看| 成年人午夜在线观看视频| 婷婷色综合www| 简卡轻食公司| 亚洲成人一二三区av| 99久久人妻综合| 乱码一卡2卡4卡精品| 伊人久久国产一区二区| 成年av动漫网址| av黄色大香蕉| 久久久久久久国产电影| 乱系列少妇在线播放| 美女视频免费永久观看网站| 丝瓜视频免费看黄片| 观看免费一级毛片| 日韩电影二区| 青春草国产在线视频| 在线 av 中文字幕| 国模一区二区三区四区视频| 人人妻人人添人人爽欧美一区卜 | 我的老师免费观看完整版| 精品酒店卫生间| 久久av网站| 另类亚洲欧美激情| av在线老鸭窝| 日日啪夜夜爽| 内地一区二区视频在线| 蜜桃亚洲精品一区二区三区| 成年免费大片在线观看| 一区二区三区乱码不卡18| 在线观看免费高清a一片| 亚洲精品国产av成人精品| 99热全是精品| 又爽又黄a免费视频| 久久久欧美国产精品| 黄色日韩在线| 男女啪啪激烈高潮av片| 亚洲欧美精品专区久久| 一区在线观看完整版| 亚洲精华国产精华液的使用体验| 久久精品国产亚洲av天美| 成人综合一区亚洲| 欧美+日韩+精品| 成人亚洲欧美一区二区av| 国产永久视频网站| 亚洲色图av天堂| 成人亚洲欧美一区二区av| 亚洲精品久久久久久婷婷小说| 全区人妻精品视频| 身体一侧抽搐| 日产精品乱码卡一卡2卡三| 国产精品秋霞免费鲁丝片| 高清av免费在线| 青青草视频在线视频观看| 久久人人爽人人爽人人片va| 美女国产视频在线观看| 五月伊人婷婷丁香| 少妇高潮的动态图| 97超碰精品成人国产| 亚洲色图av天堂| 男的添女的下面高潮视频| 亚洲欧美成人综合另类久久久| 国产一区亚洲一区在线观看| 舔av片在线| 欧美成人一区二区免费高清观看| 女性被躁到高潮视频| av国产久精品久网站免费入址| 高清不卡的av网站| 好男人视频免费观看在线| 狂野欧美白嫩少妇大欣赏| 欧美极品一区二区三区四区| 精品视频人人做人人爽| 97超碰精品成人国产| 久久亚洲国产成人精品v| 免费大片黄手机在线观看| 人体艺术视频欧美日本| 欧美最新免费一区二区三区| 身体一侧抽搐| 久久精品久久久久久噜噜老黄| 2018国产大陆天天弄谢| 国产成人a∨麻豆精品| 在线免费十八禁| 精品少妇久久久久久888优播| 日本爱情动作片www.在线观看| 日韩强制内射视频| 中国国产av一级| 国产精品一区二区在线不卡| 国产高清国产精品国产三级 | 国产精品蜜桃在线观看| 国产精品伦人一区二区| 妹子高潮喷水视频| 免费看日本二区| 中文字幕免费在线视频6| 国产深夜福利视频在线观看| 国产欧美日韩一区二区三区在线 | 亚洲欧美精品专区久久| 五月伊人婷婷丁香| 搡女人真爽免费视频火全软件| 精品一区二区免费观看| 中国三级夫妇交换| 国产精品免费大片| 亚洲精品国产色婷婷电影| 偷拍熟女少妇极品色| 卡戴珊不雅视频在线播放| 婷婷色综合www| 久久久色成人| 五月天丁香电影| 美女视频免费永久观看网站| 日韩欧美一区视频在线观看 | 免费大片黄手机在线观看| 欧美精品人与动牲交sv欧美| 美女国产视频在线观看| 最近中文字幕高清免费大全6| 色婷婷久久久亚洲欧美| 午夜福利影视在线免费观看| 视频中文字幕在线观看| 色视频在线一区二区三区| 97精品久久久久久久久久精品| 直男gayav资源| 中文资源天堂在线| 日韩一本色道免费dvd| 亚洲欧美一区二区三区黑人 | 久久精品久久精品一区二区三区| 2018国产大陆天天弄谢| 成人影院久久| 亚洲av免费高清在线观看| 大片免费播放器 马上看| 亚洲欧洲日产国产| 爱豆传媒免费全集在线观看| 欧美bdsm另类| 国产成人免费观看mmmm| 啦啦啦在线观看免费高清www| 黑人猛操日本美女一级片| 少妇裸体淫交视频免费看高清| av卡一久久| 老司机影院成人| 秋霞在线观看毛片| 亚洲国产成人一精品久久久| 涩涩av久久男人的天堂| 国产欧美日韩精品一区二区| 麻豆成人av视频| 伊人久久精品亚洲午夜| 最近中文字幕2019免费版| 身体一侧抽搐| 在线观看免费日韩欧美大片 | 国产精品一区二区在线不卡| videossex国产| 亚洲精品成人av观看孕妇| 亚洲国产毛片av蜜桃av| 精品少妇黑人巨大在线播放| 22中文网久久字幕| 中国三级夫妇交换| 日韩免费高清中文字幕av| 中文字幕制服av| 国产精品久久久久久精品电影小说 | 久久国产乱子免费精品| 国产亚洲午夜精品一区二区久久| 热99国产精品久久久久久7| 亚洲国产精品成人久久小说| 日本黄色片子视频| 亚洲第一av免费看| 久久久久国产网址| .国产精品久久| 波野结衣二区三区在线| 欧美bdsm另类| 日本爱情动作片www.在线观看| 黄色视频在线播放观看不卡| 香蕉精品网在线| 在线亚洲精品国产二区图片欧美 | av在线观看视频网站免费| 九九久久精品国产亚洲av麻豆| 国产精品免费大片| 国产av国产精品国产| 久久久久网色| 亚洲av不卡在线观看| 又黄又爽又刺激的免费视频.| 国产成人freesex在线| 成人亚洲精品一区在线观看 | 最近中文字幕2019免费版| 亚洲精品aⅴ在线观看| 国产精品福利在线免费观看| 日韩不卡一区二区三区视频在线| 各种免费的搞黄视频| 美女xxoo啪啪120秒动态图| 国模一区二区三区四区视频| 99热6这里只有精品| 精品亚洲乱码少妇综合久久| 国产又色又爽无遮挡免| 色网站视频免费| 最新中文字幕久久久久| 日产精品乱码卡一卡2卡三| 肉色欧美久久久久久久蜜桃| 日本午夜av视频| 午夜福利视频精品| 久久久久久久久久久丰满| 女人久久www免费人成看片| 少妇被粗大猛烈的视频| 少妇猛男粗大的猛烈进出视频| 蜜桃久久精品国产亚洲av| 国产探花极品一区二区| 婷婷色综合大香蕉| 精品少妇久久久久久888优播| 国产乱人偷精品视频| 久久久久久久国产电影| 国产在线免费精品| 国产一区二区在线观看日韩| 亚洲,一卡二卡三卡| 极品教师在线视频| 中文资源天堂在线| 九草在线视频观看| 天堂中文最新版在线下载| 一级片'在线观看视频| 在线观看免费高清a一片| 免费在线观看成人毛片| av在线老鸭窝| 国产精品一区二区性色av| 久久久久久久久大av| 久久午夜福利片| 久久精品国产鲁丝片午夜精品| 亚洲中文av在线| 亚洲精品aⅴ在线观看| 亚洲精品,欧美精品| 日韩中文字幕视频在线看片 | 狠狠精品人妻久久久久久综合| 在线观看一区二区三区激情| 亚洲欧美日韩东京热| av国产免费在线观看| 美女福利国产在线 | 在线精品无人区一区二区三 | 最近的中文字幕免费完整| 欧美变态另类bdsm刘玥| 日韩av在线免费看完整版不卡| 丰满迷人的少妇在线观看| 插阴视频在线观看视频| 亚洲国产色片| 亚洲av不卡在线观看| 日韩,欧美,国产一区二区三区| 又大又黄又爽视频免费| 婷婷色综合www| av专区在线播放| 国产黄色免费在线视频| 黑人高潮一二区| 最近2019中文字幕mv第一页| 蜜桃在线观看..| 男人舔奶头视频| 亚洲av电影在线观看一区二区三区| 精品一区二区三区视频在线| 久久精品国产鲁丝片午夜精品| 色综合色国产| 精品久久久久久电影网| 亚洲综合精品二区| 精品一区二区三卡| 这个男人来自地球电影免费观看 | 国产午夜精品久久久久久一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 女人十人毛片免费观看3o分钟| 热re99久久精品国产66热6| 国产欧美日韩精品一区二区| 寂寞人妻少妇视频99o| 久久久久久久久大av| av卡一久久| 欧美高清性xxxxhd video| 久久鲁丝午夜福利片| 国产一区二区三区av在线| 亚洲国产精品一区三区| 少妇熟女欧美另类| h视频一区二区三区| 制服丝袜香蕉在线| 久久久成人免费电影| 大片免费播放器 马上看| 亚洲人成网站在线播| 一个人看视频在线观看www免费| 亚洲精品日韩在线中文字幕| 国产精品三级大全| 日韩一区二区三区影片| 亚洲不卡免费看| 大香蕉久久网| 亚洲国产精品专区欧美| 18禁裸乳无遮挡免费网站照片| 在线观看免费视频网站a站| 亚洲天堂av无毛| 亚洲va在线va天堂va国产| 国产伦精品一区二区三区四那| 免费人妻精品一区二区三区视频| 天天躁日日操中文字幕| 精品亚洲成国产av| 免费高清在线观看视频在线观看| av国产久精品久网站免费入址| 熟妇人妻不卡中文字幕| 久久久久国产网址| 精品久久久久久久久亚洲| 在线观看免费高清a一片| 青春草视频在线免费观看| 久久精品夜色国产| 精品人妻一区二区三区麻豆| 身体一侧抽搐| 中国国产av一级| 嘟嘟电影网在线观看| 欧美激情极品国产一区二区三区 | 永久网站在线| 免费av中文字幕在线| 晚上一个人看的免费电影| 国产av码专区亚洲av| 国产精品一区二区在线不卡| 在线亚洲精品国产二区图片欧美 | 六月丁香七月| 亚洲av在线观看美女高潮| 99视频精品全部免费 在线| 爱豆传媒免费全集在线观看| 中文字幕精品免费在线观看视频 | 欧美成人精品欧美一级黄| 十八禁网站网址无遮挡 | 婷婷色麻豆天堂久久| 久久人人爽av亚洲精品天堂 | 亚洲av国产av综合av卡| 国产精品女同一区二区软件| 免费看日本二区| 网址你懂的国产日韩在线| 日韩三级伦理在线观看| 老司机影院成人| 97超碰精品成人国产| 黄片wwwwww| 简卡轻食公司| 不卡视频在线观看欧美| 亚洲自偷自拍三级| 妹子高潮喷水视频| 欧美日韩一区二区视频在线观看视频在线| 亚洲av电影在线观看一区二区三区| 老熟女久久久| 天堂中文最新版在线下载| 嫩草影院新地址| 国产大屁股一区二区在线视频| 久久久久久久久久久免费av| xxx大片免费视频| 久热久热在线精品观看| 美女脱内裤让男人舔精品视频| 国产精品一区二区在线不卡| 看十八女毛片水多多多| 成人国产av品久久久| www.色视频.com| 在线免费十八禁| 午夜免费男女啪啪视频观看| 日韩精品有码人妻一区| 国产精品嫩草影院av在线观看| 伦精品一区二区三区| 色婷婷av一区二区三区视频| 免费黄色在线免费观看| 色综合色国产| 一边亲一边摸免费视频| 青春草国产在线视频| kizo精华| 国内揄拍国产精品人妻在线| 精品久久久精品久久久| 国语对白做爰xxxⅹ性视频网站| 国产黄片视频在线免费观看| videos熟女内射| 街头女战士在线观看网站| 中文字幕精品免费在线观看视频 | 美女中出高潮动态图| 婷婷色综合大香蕉| 精品人妻熟女av久视频| 国产探花极品一区二区| 亚州av有码| 亚洲av不卡在线观看| 午夜福利高清视频| av卡一久久| 国产伦理片在线播放av一区| 97超视频在线观看视频| 欧美极品一区二区三区四区| 男人舔奶头视频| 欧美成人a在线观看| 久久久精品免费免费高清| av免费观看日本| av不卡在线播放| 尾随美女入室| 国产成人免费观看mmmm| 日韩 亚洲 欧美在线| 国产在线免费精品| 中文精品一卡2卡3卡4更新| 99热全是精品| 色网站视频免费| av一本久久久久| 久久精品国产a三级三级三级| av网站免费在线观看视频| 日本爱情动作片www.在线观看| 国产高清三级在线| 在线观看免费日韩欧美大片 | 国产精品国产av在线观看| 观看av在线不卡| 男人爽女人下面视频在线观看| 国产成人freesex在线| xxx大片免费视频| 如何舔出高潮| 大话2 男鬼变身卡| 人人妻人人添人人爽欧美一区卜 | 国产av精品麻豆| 男男h啪啪无遮挡| 亚洲中文av在线| 色视频在线一区二区三区| 人人妻人人看人人澡| 777米奇影视久久| 久久久久久久大尺度免费视频| 两个人的视频大全免费| 能在线免费看毛片的网站| 看免费成人av毛片| 18禁在线播放成人免费| av专区在线播放| 国产欧美另类精品又又久久亚洲欧美| 中文字幕免费在线视频6| 欧美精品一区二区免费开放| 熟女人妻精品中文字幕| 天堂中文最新版在线下载| 美女内射精品一级片tv| 亚洲国产日韩一区二区| 美女xxoo啪啪120秒动态图| 日韩一区二区视频免费看| 五月玫瑰六月丁香| 男人狂女人下面高潮的视频| 亚洲av欧美aⅴ国产| 卡戴珊不雅视频在线播放| 中国美白少妇内射xxxbb| 有码 亚洲区| 成人美女网站在线观看视频| 久久久成人免费电影| 18禁动态无遮挡网站| av播播在线观看一区| 少妇 在线观看| 在线看a的网站| 精品久久久久久久末码| 国产成人freesex在线| 在线观看国产h片| 久久婷婷青草| 高清日韩中文字幕在线| 啦啦啦啦在线视频资源| 国产爱豆传媒在线观看| 女性被躁到高潮视频| 全区人妻精品视频| 亚洲高清免费不卡视频| 国产伦理片在线播放av一区| 伦理电影免费视频| 国产av一区二区精品久久 | 成人特级av手机在线观看| 最近最新中文字幕大全电影3| 观看av在线不卡| 国产精品国产三级国产专区5o| 大话2 男鬼变身卡| 好男人视频免费观看在线| 精品少妇黑人巨大在线播放| a级毛片免费高清观看在线播放| 丝袜喷水一区| 亚洲欧洲国产日韩| 妹子高潮喷水视频| 亚洲一区二区三区欧美精品| 免费观看的影片在线观看| 国产黄片美女视频| 日本爱情动作片www.在线观看| 18+在线观看网站| 少妇的逼好多水| 最近最新中文字幕大全电影3| 又大又黄又爽视频免费| 成人二区视频| 亚洲激情五月婷婷啪啪| 女性被躁到高潮视频| 精品酒店卫生间| 一级二级三级毛片免费看| 日韩强制内射视频| 国产视频首页在线观看| 一级爰片在线观看| 亚洲一区二区三区欧美精品| 97精品久久久久久久久久精品| 最近中文字幕2019免费版| 97热精品久久久久久| www.av在线官网国产| 成人18禁高潮啪啪吃奶动态图 | 干丝袜人妻中文字幕| 国内少妇人妻偷人精品xxx网站| 99久国产av精品国产电影| 肉色欧美久久久久久久蜜桃| 欧美国产精品一级二级三级 | 亚洲精品456在线播放app| 久久久亚洲精品成人影院| 最近的中文字幕免费完整| 久久久午夜欧美精品| 国产视频内射| 街头女战士在线观看网站| 成年免费大片在线观看| 汤姆久久久久久久影院中文字幕| 青春草国产在线视频| 国产午夜精品久久久久久一区二区三区| 性色av一级| 2022亚洲国产成人精品| 麻豆乱淫一区二区| 两个人的视频大全免费| 亚洲精品日本国产第一区| 国产免费福利视频在线观看| 国产精品久久久久成人av| 国产高清不卡午夜福利| 国产成人a∨麻豆精品| 亚洲色图av天堂| 91久久精品国产一区二区成人| 伊人久久精品亚洲午夜| 综合色丁香网| 麻豆成人av视频| 美女xxoo啪啪120秒动态图| 日韩成人伦理影院| 深爱激情五月婷婷| 亚洲国产精品999| 国产av国产精品国产| 黄色一级大片看看| 亚洲高清免费不卡视频| 黄色日韩在线|