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

    InteractionbetweenHumicAcidandp-ArsanilicAcidUsingFluorescenceQuenchingApproach

    2017-10-13 03:07:57GeSiyiZhouHongSongJiaYuJiaWeiWeiLiShiyin
    生態(tài)毒理學(xué)報(bào) 2017年3期
    關(guān)鍵詞:王萌毒理學(xué)腐殖酸

    Ge Siyi, Zhou Hong, Song Jia, Yu Jia, Wei Wei, Li Shiyin,*

    1. School of Environment, Nanjing Normal University, Nanjing 210023, China2. School of Geographical Science, Nanjing Normal University, Nanjing 210023, China

    InteractionbetweenHumicAcidandp-ArsanilicAcidUsingFluorescenceQuenchingApproach

    Ge Siyi1, Zhou Hong1, Song Jia2, Yu Jia1, Wei Wei1, Li Shiyin1,*

    1. School of Environment, Nanjing Normal University, Nanjing 210023, China2. School of Geographical Science, Nanjing Normal University, Nanjing 210023, China

    Received27 December 2016accepted13 March 2017

    Abstract: P-arsanilic acid (p-ASA), as a kind of organoarsenic feed additive, has been widely used in poultry and swine breeding. However, it has caused the arsenic pollution around the farm. Currently data shows humic acid (HA) is closely to the migration and transformation of p-ASA. Therefore, the interaction between p-ASA and HA was investigated by using the method of fluorescence quenching titration. The association constant changed from 2.74 to 4.88 L·mol-1at a pH varying from 5 to 9 and reached the maximum at pH 7. In addition, logK varied from 4.15 to 5.02 L·mol-1when the temperature increased from 15 ℃ to 35 ℃. The logK increased with an increase in the concentration of HA. The dominant mechanism between p-ASA and HA is static quenching. The primary interaction force was likely the hydrogen bond, and the binding behavior occurred on the As-O stretch of p-ASA and the carboxylic acid C=O stretch of HA. The results showed that dissolved organic matters could affect the fate and biogeochemical cycling of organoarsenic pharmaceuticals in aquatic ecosystems.

    Keywords: p-arsanilic acid; humic acid; fluorescence quenching; dissolved organic matter; FTIR

    As a speciation of organoarsenic feed additive, p-arsanilic acid [(4-aminophenyl) arsonic acid,p-ASA] has been widely used in poultry and swine breeding in developing countries, as well as in other organoarsenic drugs such as roxarsone (3-nitro-4-hydroxyphenylarsonic acid), to improve the poultry’s immunity and increase its weight gain[1-3]. Several places in China were detected that the concentration of arsenic in ground/underground water are out of limits. Most of them come from the abuse of organoarsenic drugs, which is scarcely assimilated in animal bodies, and unchangeably excrete with their litter into the aquatic or soil systems[4]. Some of them might be transformed into inorganic arsenic such as As(V) and As(III) which is more toxic to ecosystems and more carcinogenic to people[4]. Several studies have shown that the migration and transformation of arsenic (including inorganic arsenic and organic arsenic) is closed to dissolved organic matter (DOM)[5-7]. However, little is known about interaction between dissolved organic matter and p-ASA.

    Dissolved organic matter is a complex mixture widely existing in soils, aquatic environments, and sediments. It consists of heterogeneous organic macromolecules, and acts as a major reactant in nature because of its various functional groups, which can be combined with different contaminants within the environment[8-9]. As a collection of various substances with low molecular weight, humic substances (HS) play a significant role in biogeochemistry of the aquatic and soil ecosystems, and it is the most important component of DOM since it represents >70% DOM[8]. The HS can be divided into three categories: humic acid (HA), which can only be soluble in alkaline condition, fulvic acid (FA), which is soluble in all, conditions, and humin, which is insoluble in all conditions[10]. Many studies have found that HA exhibits most properties of DOM; for example, HA can combine with most organic contaminants (such as PAHs and metal ions) and affect the fate of the contaminants in the environment by adjusting the existing form, solubility, mobility, transformation, and toxicity[9,11]. Currently, little is known about the interaction mechanism and the possible role of HA in controlling the environmental and biogeochemical fate of p-ASA.

    The excitation-emission matrix (EEM) fluorescence quenching introduced by Ryan and Weber (1982) is a common method to analyze the interaction between DOM and contaminants such as PAHs, C60, and various metal ions like Cu2+, Fe3+, because of the productivity and sensitivity of the EEM[10,12-13]. In our work, to evaluate the interaction between the p-ASA and HA, we select the standard fluorescence quenching method, including the three-dimensional excitation-emission matrix fluorescence (3D-EEM) and synchronous-scan fluorescence (SSFS)[13]. Adding pollutants will change the fluorescent properties of HA. We also quantify the change of fluorescence intensity for HA due to the existence of p-ASA. This quantitative data is critical to evaluate the complexation between HA and p-ASA.

    The objectives of this study are to 1) explore the binding of p-ASA and HA using the fluorescence quenching titration, 2) measure its 3D-EEM and SSFS, 3) analyze the dominant quenching mechanism between p-ASA and HA, and 4) investigate the influence of environmental factors in this system, such as the HA concentration, pH values, and temperature.

    1 Materials and methods

    1.1 Materials

    Humic acid was purchased from Sigma-Aldrich. p-Arsanilic acid with a purity of 98% was obtained from Shanghai Jinchun Industrial Co. The buffer used in the experiment is 1.0 mmol·L-1of phosphate at pH 7. All chemicals used in this study were AR grade unless mentioned otherwise, and all solutions (e.g., 3 mol·L-1of NaOH, 3 mol·L-1HCl, the stock solution of p-ASA, and phosphate buffer) were dissolved by distilled Milli-Q water (Milli-Q, Molsheim, France). The concentration of HA was calculated by dissolved organic carbon (DOC) concentrations.

    1.2 Experiments

    The titration experiments between p-ASA and HA were measured by molecular fluorescence photometer (Lumina, Thermo Fisher, America). The p-ASA stock solutions from 0 to 20 mL were added into a series of 25 mL volumetric conical flasks. A certain volume of HA and 1 mL of 1 mmol phosphate buffer were mixed and diluted to 25 mL, and N2was inflated into the mixture solution for 15 minutes to avoid static quenching caused by oxygen. Independent kinetic experiments showed that the reaction could reach equilibrium within 24 h; therefore, all mixture solutions were shaken and stored in the dark for 24 h.

    Three independent experiments were carried out, and the results were reported as averages. First, for the HA concentration-dependent experiments, a series of concentrations of HA (6.32, 12.63, and 25.26 mg·L-1) were set to react with p-ASA at pH 7 and the temperature of 25 ℃. Second, for the pH-dependent experiments, 12.63 mg·L-1of HA was titrated by a series of volumes of p-ASA at the pH of 5, 6, 7, 8, 9, respectively, adjusted with 3 mol·L-1of NaOH and 3 mol·L-1of HCl, at the temperature of 25 ℃. Third, for the temperature-dependent experiments, the mixed system of 12.63 mg·L-1of HA and p-ASA were kept at several temperatures, 25 ℃, 30 ℃, 35 ℃, 45 ℃, 55 ℃, respectively, for 24 h using a water bath.

    The 3D-EEM spectrum was recorded with a speed of 600 nm·min-1and a 0.02 s delay of response time. The excitation wavelength of this scan was set from 260 to 550 nm, and the emission wavelength of this fluorescence spectra was from 360 to 600 nm. The slit width of both the excitation wavelength and the emission wavelength was 20 nm.

    To reveal the mechanism in the interaction between HA and p-ASA, the infrared radiation of HA, p-ASA, and HA-p-ASA was measured. All solutions before measured by Fourier transform infrared (FTIR, Nicolet iS5, Thermo Fisher, USA) were rotary evaporated and freeze-dried. The mg ratio of freeze-dried sample and KBr is 1:200, and it was mixed and grounded using agate mortar[14]. The ground sample was then put into pellets and measured by the FTIR spectrometer in the range of 4 000-400 cm-1with 16 scans and a set resolution of 8 cm-1.

    1.3 Data analyses

    Static quenching was used to describe the binding behavior between several organic contamination and humic substances (e.g., carbamazepine, PAHs, and pesticides) in diluted solutions. Nearly all of the solutions could be fitted by a 1:1 model[12]. Consequently, static quenching was assumed to describe the binding mechanism between p-ASA and HA in these experiments. The data measured by fluorescence spectroscopic was analyzed by the Ryan-Weber non-linear equation and the modified Stern-Volmer plot to fit the 1:1 model and static quenching[11-12].

    2 Results

    The binding behavior of HA and p-ASA resembled the interaction between HA and the other organic matters (e.g., PAHs and carbamazepine)[11]. Consequently, the p-ASA is assumed as the quencher, and the HA is quenched. Plots in Fig.1 shows that the fluorescence intensity of p-ASA is larger less than the fluorescence intensities of HA and HA-p-ASA. Therefore, HA is the fluorescence matter, and the p-ASA is the quencher, the assumption was confirmed. According to the fluorescence spectra of kinetic experiments, the reaction of HA and p-ASA can reach equilibrium within 24 h (data not shown here for simplicity). Therefore, after adding p-ASA into the HA solution, the sample should be stirred for 24 h to reach equilibrium.2.1 Fluorescence quenching of HA with the presence of p-ASA

    The fluorescence spectra of HA, HA with addition of p-ASA are measured by a three-dimensional excitation-emission matrix fluorescence (Fig.1). According to the contour plots of HA, four major peaks appeared at excitation/emission (Ex/Em) wavelengths of 300 nm/480 nm (peak A), Ex/Em wavelengths of 370 nm/480 nm (peak B), Ex/Em wavelengths of 420 nm/500 nm (peak C), and Ex/Em wavelengths of 460 nm/520 nm (peak D). The four peaks conform the fluorescence properties of dissolved organic matter. Peak A is the main peak, which is the visible fulvic-like fluorescence. Peaks B, C, and D are humic-like fluorescence[15]. Comparing the contour plot of HA without p-ASA (Fig.2a), HA with p-ASA (Fig.2b), we found that the fluorescence intensity of HA decreased with the existence of p-ASA, and the most obvious reduction occurs at peak A. In addition, the location of peak A shifted to the Ex/Em wavelengths of 320 nm/480 nm.

    Fig.3 (plot fitted by the Ryan-Weber model) shows that the fluorescence intensity of both peak A and peak D exhibited an obvious decrease with the augment of the concentration of p-ASA. The fluorescence of peak A quenched rapidly in the initial titration; however, after adding 460.72 μmol·L-1of p-ASA, the fluorescence quenched slowly. Instead, the quench rate of peak D remains low during the entire titration. When the concentration of p-ASA reached the maximum, the quenched ratio of peak A could reach more than 30%; meanwhile, the quenched ratio of peak D is less than 13%.

    The linear plot in Fig.4 (plot fitted by the modified Stern-Volmer model) indicates the quenching mechanism of HA and p-ASA is inconsistent with a single quenching mechanism, the linearity between F0/F, and the concentration of p-ASA is significant (R2>0.95,P<0.05). Toward the plot of peak D in Fig.4, the line is more parallel to the x-axis with an increase in p-ASA. The similar phenomenon is also detected by the previous study of HA and ROX[16].

    2.2 Transient fluorescence of HA and HA-p-ASA

    Table 1 is the data obtained by applying the tri-exponential compartment model to fit the fluorescence lifetime of HA and HA with p-ASA. It shows that the fluorescence lifetime of peak A, peak B, and peak D of HA changed slightly after adding p-ASA. Adding p-ASA increased fluorescence lifetime on peak A, but decreased fluorescence lifetime on peak B, peak C and peak D. What’s more, the ratio ofτ0/τof both four peaks is approximately equal to 1 indicates that the lifetime of HA did not change by adding p-ASA.

    Fig. 1 The fluorescence intensities of HA (A), HA-p-ASA (B), p-ASA (C)

    Fig. 2 3D-EEM of HA (a) and HA after addition of p-ASA (b) (conditions: HA 12.63 mg C·L-1, p-ASA 460.72 μmol·L-1, pH 7 and temperature 25 ℃)

    Fig. 3 Fluorescence data of HA peak A and peak D quenching by p-ASA, fitting by a modified Ryan-Weber model curve (conditions: HA 12.56 mg C·L-1, pH 7 and temperature 25 ℃. F0 is the fluorescence intensity of HA, F is the fluorescence intensity of HA with addition of p-ASA)

    Fig. 4 Fluorescence data of HA peak A and peak D quenching by p-ASA fitting by a modified Stern-Volmer model curve (conditions: HA 12.56 mg C·L-1, pH 7 and temperature 25 ℃)

    Fig.5 FTIR spectra of HA (line 1), p-ASA (line 2), HA-p-ASA (line 3)

    2.3 Interaction between p-ASA and HA at different HA concentration

    The logarithms of the association constant of HA and p-ASA at three different HA concentrations are shown in Table 2. All CLcalculated by the Ryan-Weber equation of this experiment is not given here because the value is very low and could contain large errors at a low HA concentration and conditional stability constants[10]. The log K decreased from 5.39 to 2.95 L·mol-1with an addition of the concentration of HA probably implies that the increase of the HA concentration leads to the decreased binding intensity of HA and p-ASA, and the combination degree of HA-p-ASA become larger at a low HA concentration. The value of Fendvaried from 71.85 to 77.85 (which remains greater than zero), probably implying that HA and p-ASA had not reacted completely.

    2.4 Influence of pH values and temperature on the reaction of p-ASA and HA

    The logarithms of association constant of HA-p-ASA, calculated by the Ryan-Weber model, are listed in Table 3. The logK fluctuated from 2.74 to 4.48 as pH increased from 5 to 9, and it reached the highest at pH 7. Previous studies in the dissolved organic matter and metal ions show that log K increased with pH augment[11], but there is no obvious tendency between pH and logK. The value of Fenddecreased slightly with pH varying from 5 to 6, and then increased when pH increases from 6 to 9. To investigate the influence of temperature to the binding strength of HA-p-ASA, five different temperatures were set, varying from 15 ℃ to 35 ℃.The logK of peak A varies from 4.15 to 5.02 (Table 3) at temperature decreased from 15 ℃ to 35 ℃. Meantime, the log K decreased with an increase in temperature (P<0.01), which means that the binding interaction is weaken with temperature increased.

    Table 1 Fluorescence lifetime of HA and HA-p-ASA

    Table 2 The logarithms of association constant (logK) of interaction between HA and p-ASA at varying HA concentration, pH 7 and 25 ℃

    Table 3 The logarithms of association constant (logK) of the interaction between HA (12.56 mg C·L-1) and p-ASA at varying pH and temperature

    The van’t Hoff equation is used to fit the relation of association constant and temperature:

    lnK=-ΔH0/RT+r

    Where ΔH0is the enthalpy change of the reaction between p-ASA and HA, and R is the gas constant calculated as 8.134 J·(mol·K)-1. The fitted line tendency shows that a positive correlation exists between the (natural logarithm of) association constant K and T-1(R2>0.95, P<0.05). The calculated ΔH0between p-ASA and HA is 75.99 kJ·mol-1, which is consistent with the enthalpy change of hydrogen bond[17].

    2.5 FTIR spectroscopic analysis of HA and p-ASA

    The FTIR spectrum was applied to analyze the mechanism of the interaction between HA and p-ASA. The IR spectra of HA, p-ASA, and HA-p-ASA were compared to each other (Fig.5). Results show that a peak (which assigns to As-O stretching vibration) appeared at 761 cm-1and 846 cm-1of p-ASA, respectively; however, this peak has not been discovered in the IR of HA and HA-p-ASA. This discrepancy probably indicates that the As-O of stretch p-ASA has participated the reaction. Two peaks at 1 600 cm-1and 1 390 cm-1of HA IR are attributed to the asymmetrical stretching vibration and symmetrical stretching vibration of carboxylic acid C=O stretch respectively[18]. However, in the IR of HA-p-ASA, the peak at 1 390 cm-1disappeared, which is likely caused by the reaction of carboxylic acid C=O stretch. To draw a conclusion, the interaction between HA and p-ASA probably occurs on carboxylic acid C=O stretch of HA and As-O stretch of p-ASA.

    3 Discussion

    The 3D-EEM of HA with or without the presence of p-ASA shows that the fluorescence intensity of HA can be quenched by p-ASA. Comparing the results calculated by fluorescence intensity, it shows that p-ASA fluoresced weakly, and the intensity of p-ASA is less than 3.34% of HA and less than 4.79% of the mixture of HA and p-ASA. This phenomenon confirms that p-ASA is the quencher, and HA is quenched. Meanwhile, the presence of p-ASA can lead to the red shift of peak A of HA. Previous study of the interaction between HA and Hg2+ions shows that the wavelength shift might indicate the change of conformation or rigidity[19].

    According to plots fitted by Ryan-Weber model, with the increase of p-ASA concentration, the trend of F/F0is similar to the composite action of metal ions and fulvic acid and confirms the trend of static quenching mechanism, which can be described by a simple 1:1 complexation model[13]. Also, previous study confirms that the dynamic quenching can decrease the fluorescence lifetime of fluorescent group, while the static quenching does not influence the lifetime. The fact that the ratio ofτ0/τis approximately equal to 1 indicates that the lifetime of HA do not change by adding p-ASA, and this phenomenon can further prove that static quenching is the dominant interaction between HA and p-ASA. As for the logarithm of the bimolecular dynamic quenching constant, logKD, calculated by the efficient quencher, O2, is less than 2; however, in HA-p-ASA system, the logK is more than 2. It shows that the logK of between HA and p-ASA is larger than the bimolecular dynamic quenching constant logKD, which can further confirm that the dominant quenching mechanism between HA and p-ASA was static quenching[20].

    Moreover, according to Fig.2 and Fig.3, we can found that p-ASA has a greater effect on peak A than peak D. This probably because p-ASA can interact with the fluorescent group represent by peak A. However, the fluorescence intensity of peak D is just slightly quenched means that p-ASA scarcely interacts with fluorescent group represent by peak D.

    The influence of different HA concentrations, pH and temperature show that there is still a little differences between the mechanism of metal ions with HA and HA with p-ASA. Experiments demonstrate that the concentration of HA can significant influence the binding degree of HA-p-ASA. With the increase of HA concentration, the binding degree of HA-p-ASA is decreased. This probably because the augment of HA concentration change the system equilibrium, so that decreases the association constant of HA-p-ASA[21]. Also, the Fendof a higher HA concentration is larger than that of a lower HA concentration, possibly caused by the inner filtration of HA and chemical equilibrium[22]. As for the influence of pH, experiments results show that the bonding strength of HA and p-ASA is larger when the system is under a neutral condition. This probably because pH can change the conformation of function group such as phenolic hydroxyl group and carboxyl group in HA, so that influence the binding degree of the system. For the effect of temperature, with an increasing temperature, the association constant of HA-p-ASA decreases. Instead of dynamic quenching, the increase of temperature cannot increase binding constant, on the contrary, it decreased the stability of HA-p-ASA, which could decrease the quench degree of HA by p-ASA, therefore, the logK was decreased[22].This phenomenon conforms to the tendency of statistic quenching[23].

    Meanwhile, the calculated △H0between p-ASA and HA is 75.99 kJ·mol-1, which is consistent with the enthalpy change of hydrogen bond[18]. The binding sites can be detected in plots measure by FTIR. Results show peaks which represent the C=O of HA and As-O of p-ASA has changed after interaction. This implies the carboxyle group of HA and As-O of p-ASA participant the reaction. In a word, the As-O of p-ASA is connected with C=O of HA through hydrogen bond during interaction.

    p-arsanilic acid as a kind of organoarsenic feed addictive,is exhausted to environment with feces of livestock, and cause arsenic pollution, which can do harm to human health. According to research, p-ASA have static reaction with HA. Hence, it is significant to degrade p-ASA with HA, and is important to eliminate pollution. The percentage of HA bound to p-ASA (α) can be written as:

    α = [QL]/[QL]+[L]=K[Q][L]/(K[Q][L]+[Q])

    (1)

    Where [QL] is the concentration of the compound of p-ASA and HA (mol·L-1), [Q] is the concentration of p-ASA at fluorescence-binding sites (mol·L-1), [L] is the concentration of HA at unbounded binding sites (mol·L-1). Defining CLas the total concentration of ligand L (HA), AssumingCL>>C0, CL≈ L, equation (1) can be simplified as:

    α= KCL/(KCL+1) (CL>>C0)

    (2)

    The percentage of bounded HA is calculated based on a simplified model by the association constant at pH 7, 12.56 mg·L-1HA concentration, and 25 ℃ temperature, which is most approaching the natural condition. The calculatedαaccounts for 69.7%, indicating a large number fraction of HA would interact with p-ASA. The association constant of combined complex HA-p-ASA approaches to pKa of weak acid such as formic acid (pKa 3.74), according previous research, this probably means bounded p-ASA will be released when HA in aquatic system decreased, however, this deduction still need to be proved in following study. What’s more, these data also can demonstrate that the migration and transformation of p-ASA connect tightly with aquatic dissolved organic matter by increase p-ASA solubility to motivate its mobility. In addition, the transformation of p-ASA from sediment to aquatic or soil can also be influenced by DOM by adsorption and complexation, which deserved further study.

    In summary, HA fluorescence intensity can be apparently quenched with the addition of p-ASA. The domain quench mechanism between HA and p-ASA is static quenching. The increased HA concentration promotes the logK of HA-p-ASA at peak A and the logK of HA-p-ASA is more than the bimolecular dynamic quenching constant. The complexation model between HA and p-ASA can probably be the 1:1 complexation model. The logK changes from 2.74 to 4.88 at a pH varying from 5 to 9, and it reaches the largest binding degree at a pH of 7. The logK decreases with the increase of temperature, and varies from 4.15 to 5.02 when temperature increases from 15 ℃ to 35 ℃. The enthalpy change of calculated ΔH0between HA and p-ASA is 75.99 kJ·mol-1, which means that the interaction force of HA and p-ASA was hydrogen bond. The binding sites of the interaction between HA and p-ASA are located at the As-O stretch of p-ASA and carboxylic acid C=O stretch of HA.

    Acknowledgements:This work was supported by the Major Project of the National Water Pollution Control (2012ZX07101-008), the National Natural Science Foundation of China (Grant No. 41373111).

    Corresponding author:Li Shiyin, Male, Assistant Dean and Professor in School of Environment, Nanjing Normal University, Nanjing 210046, China. Research direction: environmental toxicology.

    [1] 王素芬,寧建鳳,崔理華,等. 外源洛克沙胂對水生態(tài)系統(tǒng)砷蓄積及底泥微生物生長的影響 [J]. 生態(tài)毒理學(xué)報(bào), 2015, 10(3): 177-183

    Wang S, Ning J, Cui L, et al. Impacts of exogenous roxarsone on arsenic accumulation and growth of sediment microorganisms in aquatic microcosms [J]. Asian Journal of Ecotoxicology, 2015, 10(3): 173-183 (in Chinese)

    [2] Wang L L, Cheng H F. Birnessite (δ-MnO2) mediated degradation of organoarsenic feed additive p-arsanilic acid [J]. Environmental Science & Technology, 2015, 49(6):3473-3481

    [3] Jun J W, Tong M, Jung B K, et al. Effect of central metal ions of analogous metal-organic frameworks on adsorption of organoarsenic compounds from water: Plausible mechanism of adsorption and water purification [J]. Chemistry - A European Journal, 2015, 21(1): 347-354

    [4] Arts D, Abdus S M, Al-Abadleh H A. Surface interactions of aromatic organoarsenical compounds with hematite nanoparticles using ATR-FTIR: Kinetic studies [J]. Journal of Physical Chemistry A, 2013, 117(10): 2195-2204

    [5] Chen T C, Hseu Z Y, Jean J S, et al. Association between arsenic and different-sized dissolved organic matter in the groundwater of black-foot disease area, Taiwan [J]. Chemosphere, 2016, 159: 214-220

    [6] Mladenov N, Zheng Y, Simone B, et al. Dissolved organic matter quality in a shallow aquifer of Bangladesh: Implications for arsenic mobility [J]. Environmental Science & Technology, 2015, 49(18): 10815-10824

    [7] Zhu J P, Mei T, Peng Y, et al. Characterizing the interaction between roxarsone and humic acid by fluorescence quenching experiment [J]. China Journal of Environment Science, 2014, 35(7): 2620-2626

    [8] 王萌, 馬陶武, 龍奕, 等. 腐殖酸作用下沉積物中納米氧華銅對銅銹環(huán)棱螺生態(tài)毒性的影響[J]. 生態(tài)毒理學(xué)報(bào), 2014, 9(4): 803-808

    Wang M, Ma T, Long Y, et al. Impacts of humic acid on ecotoxicity of CuO nanoparticle in Bellamya aeruginosa under the sediment scenario [J]. Asian Journal of Ecotoxicology, 2014, 9(4): 803-808 (in Chinese)

    [10] Bai Y, Wu F, Liu C, et al. Interaction between carbamazepine and humic substances: A fluorescence spectroscopy study [J]. Environmental Toxicology and Chemistry, 2008, 27(1): 95-102

    [11] Esteves da Silva J C G, Machado A A S C, Oliveira C J S, et al. Fluorescence quenching of anthropogenic fulvic acids by Cu(II), Fe(III) and UO(2)(2+) [J]. Talanta, 1998, 45(6): 1155-1165

    [12] Ryan D K, Weber J H. Fluorescence quenching titration for determination of complexing capacities and stability constants of fulvic acid [J]. Cancer Biotherapy & Radiopharmaceuticals, 1982, 54(6): 12-22

    [13] Gauthier T D, Shane E C, Guerin W F, et al. Fluorescence quenching method for determining equilibrium constants for polycyclic aromatic hydrocarbons binding to dissolved humic materials [J]. Environmental Science & Technology, 1986, 20(11): 1162-1166

    [14] Rodríguez F J, Schlenger P, García-Valverde M. Monitoring changes in the structure and properties of humic substances following ozonation using UV-Vis, FTIR and 1 H NMR techniques [J]. Science of the Total Environment, 2016, 541: 623-637

    [15] Coble P G. Characterization of marine and terrestrial DOM in seawater using excitation-emission matrix spectroscopy [J]. Marine Chemistry, 1996, 51(4): 325-346

    [16] Fu Q L, He J Z, Blaney L, et al. Roxarsone binding to soil-derived dissolved organic matter: Insights from multi-spectroscopic techniques [J]. Chemosphere, 2016, 155:225-233

    [17] Bai Y C, Wu F C, Liu C Q, et al. Ultraviolet absorbance titration for determining stability constants of humic substances with Cu(II) and Hg(II) [J]. Analytica Chimica Acta, 2008, 616(1): 115-121

    [18] Bo P, Qiu M, Min W, et al. The opposite impacts of Cu and Mg cations on dissolved organic matter-ofloxacin interaction [J]. Environmental Pollution, 2012, 161(1): 76-82

    [19] Wu F, Cai Y, Evans D, et al. Complexation between Hg(II) and dissolved organic matter in stream waters: An application of fluorescence spectroscopy [J]. Biogeochemistry, 2004, 71(3): 339-351

    [20] Meissl K, Smidt E, Schwanninger M. Prediction of humic acid content and respiration activity of biogenic waste by means of Fourier transform infrared (FTIR) spectra and partial least squares regression (PLS-R) models [J]. Talanta, 2007, 72(2): 791-799

    [21] Pan B, Ghosh S, Xing B. Dissolved organic matter conformation and its interaction with pyrene as affected by water chemistry and concentration [J]. Environmental Science & Technology, 2008, 42(5): 1594-1599

    [22] Garg A, Manidhar D M, Gokara M, et al. Elucidation of the binding mechanism of coumarin derivatives with human serum albumin [J]. Plos One, 2013, 8(5): 1-17

    [23] Garbarino J R, Bednar A J, Rutherford D W, et al. Environmental fate of roxarsone in poultry litter. I. Degradation of roxarsone during composting [J]. Environmental Science & Technology, 2003, 37(8): 1509-1514

    10.7524/AJE.1673-5897.20161227002

    Foundation:National Natural Science Foundation of China (41373111)

    FirstauthorGe Siyi(1994-), female, postgraduate student of Nanjing Normal University, research direction is environmental toxicology, E-mail: 610530398@qq.com

    *Correspondingauthor,E-mail: lishiyin@njnu.edu.cn

    Ge S Y, Zhou H, Song J. Interaction between humic acid and p-arsanilic acid using fluorescence quenching approach [J]. Asian Journal of Ecotoxicology, 2017, 12(3): 667-674

    猜你喜歡
    王萌毒理學(xué)腐殖酸
    Scanning the optical characteristics of lead-free cesium titanium bromide double perovskite nanocrystals
    “學(xué)會道歉”寫作指導(dǎo)及例文展評
    Improved Five-Parameter Exponential-Type Potential Energy Model for Diatomic Molecules?
    豬糞中添加腐殖酸添加劑可降低糞便中的臭氣
    PM2.5毒理學(xué)實(shí)驗(yàn)染毒方法及毒理學(xué)效應(yīng)
    含腐殖酸固廢肥對菠菜產(chǎn)量及效益的影響
    災(zāi)害毒理學(xué)理論研究初探
    火災(zāi)毒理學(xué)若干問題的探討
    王萌作品選
    季銨型陽離子纖維素對水中腐殖酸的吸附
    久久影院123| 免费大片18禁| 精品久久蜜臀av无| 99国产综合亚洲精品| 欧美激情极品国产一区二区三区 | 黑人高潮一二区| 亚洲伊人久久精品综合| 久久影院123| 国产精品免费大片| 18禁在线播放成人免费| 美女内射精品一级片tv| 免费观看的影片在线观看| 婷婷色综合www| 国产精品久久久久久精品电影小说| 久久久久精品性色| 国产精品一二三区在线看| 97精品久久久久久久久久精品| 亚洲国产精品成人久久小说| 女的被弄到高潮叫床怎么办| 国产一区有黄有色的免费视频| 一区二区日韩欧美中文字幕 | 精品视频人人做人人爽| 国产探花极品一区二区| 欧美日韩综合久久久久久| 91精品一卡2卡3卡4卡| 日本黄色日本黄色录像| 日本午夜av视频| 3wmmmm亚洲av在线观看| 在线 av 中文字幕| 国产精品熟女久久久久浪| 黄色视频在线播放观看不卡| 亚洲人成网站在线播| 精品亚洲成a人片在线观看| 久热这里只有精品99| 久久这里有精品视频免费| 久久精品国产亚洲av涩爱| 久久精品人人爽人人爽视色| 只有这里有精品99| 日本猛色少妇xxxxx猛交久久| 久久ye,这里只有精品| 视频中文字幕在线观看| 一级毛片我不卡| 久久久久精品性色| 久久久a久久爽久久v久久| 国产精品国产av在线观看| 看十八女毛片水多多多| 超色免费av| 色婷婷av一区二区三区视频| 亚洲欧洲日产国产| av国产久精品久网站免费入址| 国产精品99久久99久久久不卡 | 国产精品国产av在线观看| 人妻人人澡人人爽人人| 999精品在线视频| 亚洲av成人精品一区久久| 麻豆成人av视频| 久久av网站| 大陆偷拍与自拍| 两个人的视频大全免费| 亚洲国产av新网站| 亚洲欧洲精品一区二区精品久久久 | 国产精品三级大全| 精品少妇内射三级| 一本一本综合久久| 精品99又大又爽又粗少妇毛片| 久久国产精品男人的天堂亚洲 | 黄色视频在线播放观看不卡| 久久久久久久国产电影| 成人午夜精彩视频在线观看| 成年美女黄网站色视频大全免费 | 我要看黄色一级片免费的| 简卡轻食公司| 欧美成人午夜免费资源| 国产午夜精品一二区理论片| 黄色一级大片看看| 亚洲av中文av极速乱| 亚洲av成人精品一二三区| 亚洲欧洲日产国产| 免费日韩欧美在线观看| 久久精品国产亚洲av天美| 久久久久久久国产电影| av视频免费观看在线观看| 亚洲成色77777| 一区二区三区精品91| 精品人妻偷拍中文字幕| 精品人妻熟女毛片av久久网站| 亚洲精品国产av蜜桃| 人妻系列 视频| 欧美 亚洲 国产 日韩一| 天堂8中文在线网| 日韩av免费高清视频| 街头女战士在线观看网站| 亚洲精华国产精华液的使用体验| 久久久精品区二区三区| 免费观看无遮挡的男女| 久久人人爽av亚洲精品天堂| 国产精品无大码| 欧美xxⅹ黑人| 国产在线免费精品| av视频免费观看在线观看| tube8黄色片| 欧美日韩视频高清一区二区三区二| 久久青草综合色| 亚洲国产精品成人久久小说| 一级黄片播放器| 91成人精品电影| 亚洲,欧美,日韩| 国产女主播在线喷水免费视频网站| 3wmmmm亚洲av在线观看| 亚洲欧洲日产国产| 一本久久精品| 纯流量卡能插随身wifi吗| 青春草视频在线免费观看| 不卡视频在线观看欧美| av在线观看视频网站免费| 国产精品秋霞免费鲁丝片| 欧美性感艳星| 91在线精品国自产拍蜜月| 一本大道久久a久久精品| 嘟嘟电影网在线观看| 免费少妇av软件| 久久久久久久久久久免费av| 亚洲第一区二区三区不卡| 街头女战士在线观看网站| 国产精品人妻久久久影院| 亚洲欧美日韩卡通动漫| 久久精品夜色国产| 日本欧美视频一区| 国产高清不卡午夜福利| 婷婷色综合大香蕉| 欧美少妇被猛烈插入视频| 午夜福利影视在线免费观看| 中文字幕亚洲精品专区| 不卡视频在线观看欧美| av播播在线观看一区| 自拍欧美九色日韩亚洲蝌蚪91| 91精品伊人久久大香线蕉| 欧美xxxx性猛交bbbb| 高清在线视频一区二区三区| 99re6热这里在线精品视频| 国产成人精品久久久久久| 亚洲熟女精品中文字幕| 免费av中文字幕在线| 亚洲国产最新在线播放| 成人无遮挡网站| 亚洲精品中文字幕在线视频| 亚洲性久久影院| 免费人妻精品一区二区三区视频| 精品一区在线观看国产| 日韩亚洲欧美综合| 久久精品国产a三级三级三级| 97超碰精品成人国产| 国产免费现黄频在线看| 69精品国产乱码久久久| 啦啦啦中文免费视频观看日本| 一级毛片我不卡| 精品少妇久久久久久888优播| 69精品国产乱码久久久| 亚洲美女视频黄频| 在线观看三级黄色| 爱豆传媒免费全集在线观看| 国产黄频视频在线观看| 日韩视频在线欧美| 秋霞在线观看毛片| 最近中文字幕高清免费大全6| 乱码一卡2卡4卡精品| 国产日韩欧美在线精品| 涩涩av久久男人的天堂| 熟妇人妻不卡中文字幕| 丝瓜视频免费看黄片| 99热这里只有是精品在线观看| 亚洲精品久久成人aⅴ小说 | 五月开心婷婷网| 日韩一本色道免费dvd| 亚洲综合色网址| videossex国产| 成年av动漫网址| 熟女av电影| 欧美精品一区二区大全| 国产精品三级大全| 亚洲国产色片| 精品一区二区三卡| 天天操日日干夜夜撸| 成人国产av品久久久| av.在线天堂| 最黄视频免费看| 成人18禁高潮啪啪吃奶动态图 | 免费人成在线观看视频色| 性高湖久久久久久久久免费观看| 女人精品久久久久毛片| 久久国产亚洲av麻豆专区| 啦啦啦啦在线视频资源| 亚洲少妇的诱惑av| 精品人妻熟女毛片av久久网站| 免费人妻精品一区二区三区视频| 久久99热这里只频精品6学生| 国产精品国产三级国产专区5o| 午夜免费男女啪啪视频观看| 成人手机av| 国产片内射在线| 国产片特级美女逼逼视频| 免费黄频网站在线观看国产| 色网站视频免费| 国产国语露脸激情在线看| 91久久精品国产一区二区三区| 国产片内射在线| 熟女人妻精品中文字幕| 久久久精品免费免费高清| 国产有黄有色有爽视频| 国产精品人妻久久久影院| 婷婷色综合www| 久久久久久久久久久丰满| 国产精品一区二区三区四区免费观看| 精品久久久久久电影网| 最后的刺客免费高清国语| 国产精品蜜桃在线观看| 在线精品无人区一区二区三| 这个男人来自地球电影免费观看 | 人人澡人人妻人| 国产 精品1| 久久这里有精品视频免费| 人成视频在线观看免费观看| 又粗又硬又长又爽又黄的视频| 2021少妇久久久久久久久久久| 免费黄网站久久成人精品| 亚洲中文av在线| 欧美日韩成人在线一区二区| 成人漫画全彩无遮挡| 亚洲精品中文字幕在线视频| 日本黄色日本黄色录像| 欧美人与性动交α欧美精品济南到 | 99热国产这里只有精品6| 在线观看三级黄色| 九色亚洲精品在线播放| 国产欧美日韩一区二区三区在线 | 国产乱人偷精品视频| 色94色欧美一区二区| 日韩精品有码人妻一区| 精品久久久精品久久久| 国产色婷婷99| 欧美日韩成人在线一区二区| 一级毛片电影观看| 国产精品久久久久久精品古装| 亚洲精华国产精华液的使用体验| 在线观看免费日韩欧美大片 | 国产av国产精品国产| 亚洲国产欧美日韩在线播放| 欧美日本中文国产一区发布| 99热国产这里只有精品6| 人妻人人澡人人爽人人| 久热这里只有精品99| 亚洲一区二区三区欧美精品| av天堂久久9| 麻豆乱淫一区二区| 国产成人精品无人区| 一区二区三区四区激情视频| 在线播放无遮挡| 最近最新中文字幕免费大全7| 丝袜在线中文字幕| 亚洲精品久久久久久婷婷小说| 最后的刺客免费高清国语| 欧美激情国产日韩精品一区| 午夜久久久在线观看| 国产黄色视频一区二区在线观看| av视频免费观看在线观看| 看免费成人av毛片| 视频区图区小说| 卡戴珊不雅视频在线播放| 亚洲精品一区蜜桃| 精品国产露脸久久av麻豆| 精品人妻偷拍中文字幕| 桃花免费在线播放| 美女cb高潮喷水在线观看| 欧美精品高潮呻吟av久久| 人成视频在线观看免费观看| 91久久精品国产一区二区成人| 国产精品久久久久久久久免| 国产在线免费精品| 爱豆传媒免费全集在线观看| 51国产日韩欧美| 中国美白少妇内射xxxbb| 国产免费一级a男人的天堂| 一级a做视频免费观看| 精品一区二区三区视频在线| 97超碰精品成人国产| 国产 一区精品| 亚洲精品av麻豆狂野| 内地一区二区视频在线| 国产精品99久久久久久久久| 五月伊人婷婷丁香| 国产精品熟女久久久久浪| 亚洲丝袜综合中文字幕| 成年美女黄网站色视频大全免费 | 日韩免费高清中文字幕av| 精品亚洲乱码少妇综合久久| 国产成人午夜福利电影在线观看| 国产午夜精品久久久久久一区二区三区| 伦理电影免费视频| 国产精品久久久久久久久免| 成年美女黄网站色视频大全免费 | 最近手机中文字幕大全| 亚洲三级黄色毛片| 国产永久视频网站| 观看美女的网站| av女优亚洲男人天堂| 国产又色又爽无遮挡免| 全区人妻精品视频| 男女边吃奶边做爰视频| 麻豆乱淫一区二区| 成人亚洲精品一区在线观看| 99热国产这里只有精品6| 亚洲婷婷狠狠爱综合网| 男女无遮挡免费网站观看| 精品人妻偷拍中文字幕| 免费观看性生交大片5| 男人爽女人下面视频在线观看| 亚洲精品国产av蜜桃| 久久精品国产鲁丝片午夜精品| 亚洲欧美中文字幕日韩二区| 国产乱来视频区| 十八禁高潮呻吟视频| 高清黄色对白视频在线免费看| 亚洲国产欧美在线一区| 建设人人有责人人尽责人人享有的| 久久av网站| 亚洲成色77777| 国产男人的电影天堂91| a级毛色黄片| 欧美激情极品国产一区二区三区 | 欧美人与善性xxx| 欧美xxⅹ黑人| 亚洲美女视频黄频| a 毛片基地| 国产精品人妻久久久久久| freevideosex欧美| av国产精品久久久久影院| 国产乱来视频区| 久久久国产精品麻豆| 国产欧美日韩综合在线一区二区| 老司机亚洲免费影院| av视频免费观看在线观看| 少妇被粗大猛烈的视频| 国产精品99久久99久久久不卡 | 免费看av在线观看网站| 国产成人一区二区在线| 人人妻人人澡人人看| 九色亚洲精品在线播放| 亚洲国产日韩一区二区| 大码成人一级视频| 欧美亚洲日本最大视频资源| 免费观看在线日韩| 久久综合国产亚洲精品| 伊人久久国产一区二区| 男的添女的下面高潮视频| 黑人猛操日本美女一级片| 纯流量卡能插随身wifi吗| 99re6热这里在线精品视频| 日韩精品有码人妻一区| 2021少妇久久久久久久久久久| 亚洲欧美成人精品一区二区| 免费观看a级毛片全部| 人妻人人澡人人爽人人| 秋霞伦理黄片| 如日韩欧美国产精品一区二区三区 | 韩国高清视频一区二区三区| 一区二区三区精品91| 国产精品一二三区在线看| 超色免费av| 麻豆精品久久久久久蜜桃| 高清黄色对白视频在线免费看| 我要看黄色一级片免费的| 人人妻人人澡人人爽人人夜夜| 日日摸夜夜添夜夜爱| 欧美人与善性xxx| 色视频在线一区二区三区| 亚洲国产精品999| 一级二级三级毛片免费看| 色5月婷婷丁香| 亚洲国产精品一区三区| 亚洲精品美女久久av网站| 国产av国产精品国产| 少妇高潮的动态图| 卡戴珊不雅视频在线播放| 国产成人精品福利久久| 亚洲国产最新在线播放| 寂寞人妻少妇视频99o| 国产亚洲精品久久久com| 亚洲欧美一区二区三区黑人 | 亚洲成人手机| av.在线天堂| 国产精品嫩草影院av在线观看| av.在线天堂| 免费日韩欧美在线观看| 男女边摸边吃奶| 久久精品夜色国产| 国产精品无大码| 如何舔出高潮| 色婷婷av一区二区三区视频| 一级爰片在线观看| 2018国产大陆天天弄谢| 日韩在线高清观看一区二区三区| 精品少妇久久久久久888优播| 久久毛片免费看一区二区三区| 国产成人91sexporn| 久久久久精品性色| 午夜免费男女啪啪视频观看| 五月伊人婷婷丁香| 国产在视频线精品| 少妇的逼水好多| 国产精品国产三级国产av玫瑰| 亚洲精品日本国产第一区| 精品人妻熟女毛片av久久网站| 成人国产av品久久久| 在线观看免费高清a一片| 亚洲人成网站在线观看播放| 草草在线视频免费看| 亚洲成人手机| 又黄又爽又刺激的免费视频.| 久久精品久久久久久久性| 黑丝袜美女国产一区| 天堂8中文在线网| 国产精品麻豆人妻色哟哟久久| 国产精品无大码| 久久久久精品性色| 精品人妻偷拍中文字幕| 国产精品久久久久久精品古装| 69精品国产乱码久久久| 精品亚洲成a人片在线观看| 久久久久久久久大av| 涩涩av久久男人的天堂| 精品久久国产蜜桃| 亚洲久久久国产精品| 亚洲色图综合在线观看| 在线天堂最新版资源| a级毛片黄视频| 2022亚洲国产成人精品| 亚洲国产欧美日韩在线播放| 亚洲av男天堂| 99久久精品一区二区三区| 亚洲精品国产色婷婷电影| 91精品三级在线观看| 免费观看性生交大片5| 夫妻性生交免费视频一级片| 国产片特级美女逼逼视频| 亚洲av国产av综合av卡| av.在线天堂| 一本久久精品| 日韩av不卡免费在线播放| 黄色视频在线播放观看不卡| 国产一区二区三区av在线| 在线亚洲精品国产二区图片欧美 | 日韩电影二区| 人人澡人人妻人| 久久久久久久久久久免费av| 国产精品嫩草影院av在线观看| 亚洲精品亚洲一区二区| 嘟嘟电影网在线观看| www.色视频.com| 精品一区二区三卡| 青青草视频在线视频观看| 国产免费现黄频在线看| 精品一区二区三区视频在线| 国产午夜精品久久久久久一区二区三区| 18禁观看日本| 国产成人精品婷婷| 高清视频免费观看一区二区| 午夜91福利影院| 大话2 男鬼变身卡| 性色av一级| 亚洲怡红院男人天堂| 成人综合一区亚洲| 国产亚洲午夜精品一区二区久久| 免费看光身美女| 好男人视频免费观看在线| 亚洲欧美精品自产自拍| 建设人人有责人人尽责人人享有的| 又大又黄又爽视频免费| 国产一区有黄有色的免费视频| 男女国产视频网站| 性高湖久久久久久久久免费观看| 国产老妇伦熟女老妇高清| 国产黄色免费在线视频| 欧美日韩亚洲高清精品| 日韩中文字幕视频在线看片| 国产伦理片在线播放av一区| 亚洲高清免费不卡视频| 永久免费av网站大全| 国产69精品久久久久777片| 菩萨蛮人人尽说江南好唐韦庄| 免费黄频网站在线观看国产| av天堂久久9| 少妇高潮的动态图| 亚洲国产色片| 观看av在线不卡| 乱码一卡2卡4卡精品| 亚洲av.av天堂| 国产日韩欧美亚洲二区| 欧美国产精品一级二级三级| 91精品三级在线观看| 亚洲综合精品二区| 久久99精品国语久久久| 中文字幕人妻熟人妻熟丝袜美| 亚洲av福利一区| 国产精品不卡视频一区二区| 搡老乐熟女国产| av有码第一页| av不卡在线播放| 国产精品国产三级专区第一集| 欧美成人午夜免费资源| 黑人猛操日本美女一级片| av有码第一页| 精品人妻熟女av久视频| av福利片在线| 久久国产精品男人的天堂亚洲 | 免费看av在线观看网站| 免费不卡的大黄色大毛片视频在线观看| 91精品一卡2卡3卡4卡| 美女内射精品一级片tv| 男的添女的下面高潮视频| 18禁动态无遮挡网站| 啦啦啦在线观看免费高清www| 99国产精品免费福利视频| 一本大道久久a久久精品| 一区二区三区精品91| 中文字幕人妻熟人妻熟丝袜美| 大片电影免费在线观看免费| www.av在线官网国产| 亚洲av免费高清在线观看| 国产黄频视频在线观看| 又黄又爽又刺激的免费视频.| 9色porny在线观看| 人体艺术视频欧美日本| 丰满饥渴人妻一区二区三| 日韩制服骚丝袜av| 久久这里有精品视频免费| 国产免费视频播放在线视频| 日韩人妻高清精品专区| 国产在线一区二区三区精| 伦理电影免费视频| 精品国产一区二区三区久久久樱花| 久久久久精品性色| 最新中文字幕久久久久| 午夜视频国产福利| 午夜福利视频精品| 日韩av在线免费看完整版不卡| 亚洲图色成人| 久久午夜综合久久蜜桃| 日本免费在线观看一区| 亚洲精品456在线播放app| 亚洲一级一片aⅴ在线观看| 亚洲性久久影院| av不卡在线播放| 观看av在线不卡| 日韩精品免费视频一区二区三区 | av女优亚洲男人天堂| 日韩伦理黄色片| 午夜视频国产福利| 亚洲国产最新在线播放| 中文欧美无线码| 丝袜喷水一区| 免费观看性生交大片5| 中文天堂在线官网| 亚洲图色成人| 99九九线精品视频在线观看视频| 亚洲第一区二区三区不卡| 嘟嘟电影网在线观看| 国产男女内射视频| av有码第一页| 99热全是精品| av又黄又爽大尺度在线免费看| 精品少妇黑人巨大在线播放| av女优亚洲男人天堂| 亚洲美女视频黄频| 尾随美女入室| 久久久久网色| 少妇丰满av| 国产成人午夜福利电影在线观看| 亚洲精品国产av成人精品| 哪个播放器可以免费观看大片| 国产色婷婷99| 亚洲成人手机| 国产精品成人在线| 国产日韩欧美在线精品| 91久久精品国产一区二区成人| 视频在线观看一区二区三区| 丝袜喷水一区| 黄色毛片三级朝国网站| 久久久午夜欧美精品| 91成人精品电影| av在线播放精品| 22中文网久久字幕| 国产亚洲精品第一综合不卡 | 晚上一个人看的免费电影| 亚洲av福利一区| 亚洲经典国产精华液单| 国产亚洲一区二区精品| 大片电影免费在线观看免费| .国产精品久久| 国产亚洲一区二区精品| 国产亚洲欧美精品永久| 久久久久久久久久久免费av| 99热这里只有是精品在线观看| 在线观看三级黄色| 亚洲成色77777| 国产无遮挡羞羞视频在线观看| 国产男人的电影天堂91| 波野结衣二区三区在线| a级毛色黄片| 国产国拍精品亚洲av在线观看| 久久国产精品大桥未久av| 欧美亚洲日本最大视频资源| 黑人猛操日本美女一级片| 欧美日韩一区二区视频在线观看视频在线| 久久久久久久大尺度免费视频| 国产高清国产精品国产三级|