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

    Wind induced turbulence caused colony disaggregation and morphological changes in the cyanobacterium Microcystis*

    2021-03-10 07:32:06ZhuWeiChenHuaiminXiaoManrlingMiquelLiMing
    湖泊科學(xué) 2021年2期
    關(guān)鍵詞:水柱生長率微囊

    Zhu Wei, Chen Huaimin, Xiao Man, Lürling Miquel & Li Ming

    (1: Ministry of Education Key Laboratory of Integrated Regulation and Resource Development on Shallow Lakes, College of Environment, Hohai University, Nanjing 210098, P.R.China)(2: Australian Rivers Institute, Griffith University, Nathan, QLD 4111, Australia)(3: Aquatic Ecology & Water Quality Management Group, Department of Environmental Sciences, Wageningen University, P.O. Box 47, 6700 AA, Wageningen, The Netherlands)(4: College of Natural Resources and Environment, Northwest A & F University,Yangling 712100, P.R.China)

    Abstract: Colony size and morphology influence the vertical migration capacity of Microcystis and therewith the occurrence of surface accumulations or blooms. To explore the influence of wind-induced turbulence on the colony size and morphology of Microcystis in field conditions, a high-frequency field investigation was conducted in an enclosure in Meiliang Bay, Lake Taihu, China, from 26 August to 7 September 2012. A Pearson’s correlation analysis indicated that the in situ growth rate of Microcystis was negatively related to surface cell density and positively related to wind speed. Strong wind speed stimulated Microcystis growth by enhancing light transmission due to the homodispersion of Microcystis in the water column. The Microcystis colony size was negatively correlated with wind speed, suggesting that wind-induced turbulence could break up colonies in shallow lakes. The results indicated that Microcystis colonies could be broken up by a turbulence intensity of 2.33×10-5 m2/s3, which corresponds to an average wind speed of 5.00 m/s in a reservoir with 30-m depth. Different Microcystis morphotypes were present and negative relationships were detected between the proportion of Microcystis ichthyoblabe and the proportion of Microcystis wesenbergii and between the proportion of Microcystis aeruginosa and the proportion of Microcystis wesenbergii out of all Microcystis throughout the water column, but more evidence is need to support the hypothesis that the morphology of Microcystis colonies changes over time. Altogether, the results suggest that declining wind speed, driven by climate change, will promote surface blooms of Microcystis due to the formation of larger colonies.

    Keywords: Microcystis; water blooms; turbulence; colony size; morphology; Lake Taihu

    The cyanobacteriumMicrocystismay form dense blooms throughout water systems around the world, which represents severe ecological and environmental issues[1-3]. The vertical migration capacity ofMicrocystisis an important trait underlying surface accumulations[4-5]. This migration ability depends on the density of the cyanobacteria-interplay between gas vesicles and carbohydrate ballast-, as well as the colony size and morphology, which are affected by turbulence[6-7]. Therefore, insight is needed in the influence of turbulence on colony size and morphology ofMicrocystisto understand the effect of turbulence on cyanobacterial blooms and surface accumulations.

    The effect of turbulence on the colony size ofMicrocystishas mainly been studied in the laboratory. Some scholars[8-10]found that small-scale turbulence could increase the colony size ofMicrocystisin laboratory experiments. However, this was mainly because small-scale turbulence promoted the mixing of carbon dioxide and nutrients in the culture equipment, as well as the colonies, which brought them continuously in higher light that, consequently, improvedMicrocystisgrowth. Li et al.[7]reported that artificial high turbulence could change the colony size and morphology of someMicrocystisspecies. Nevertheless, their work indicated that the natural wind-induced turbulence could not disaggregateMicrocystiscolonies, whereas O’Brien et al.[11]got the opposite result.

    Only one study has explored the effect of wind-waves on the colony size ofMicrocystisin field conditions[12]. These authors found that the mean size ofMicrocystiscolonies in the water column increased from 32.8 to 69.4 μm within 48 h during the passage of Typhoon Soulik, when the average wind speed was 6.63 m/s[12]. This result does not rule out the possibility that wind-waves changed the microenvironment so that the growth ofMicrocystiswas promoted, thereby increasing the colony size. Therefore, it is necessary to conduct high-frequency systematic monitoring in the field to study the relationship between wind-waves or turbulence and the colony size ofMicrocystis.

    The differentMicrocystiscolonies have been interpreted as being different species[13]. There are obvious seasonal successions ofM.ichthyoblabe,M.wesenbergii, andM.aeruginosain most eutrophic lakes[14-16]. Recently, Xiao et al.[13]proposed a conceptual model of the morphological transition of aMicrocystiscolony, pointing out phenotypic plasticity rather than species replacement to explain variousMicrocystismorphotypes. Subsequently, Li et al.[7]confirmed in an indoor experiment that turbulence could change the colony morphology ofMicrocystis. All of these results suggest thatMicrocystiswith different morphologies belong to the same species and that morphology is constantly changing. However, it is important to conduct field research to obtain evidence for the transition of colony morphology. According to the morphology change model proposed by Xiao et al.[13], turbulence could potentially driveM.ichthyoblabetoM.wesenbergii-like colonies and toM.aeruginosa-like colonies. Hence, reduced biomass inM.ichthyoblabewould be accompanied by increased biomass ofM.wesenbergii, and likewise, reduced biomass inM.wesenbergiiwould be accompanied by increased biomass ofM.aeruginosa. Thus, the proportion ofM.ichthyoblabeshould be negatively related to the proportion ofM.wesenbergiiand the proportion ofM.wesenbergiibiomass should be negatively related to the proportion ofM.aeruginosa. To verify the above hypothesis, it is necessary to conduct high-frequency measurements of the biomass of different morphologies ofMicrocystisin lakes.

    To this end, it is necessary to analyze colony size and morphological changes inMicrocystisunder different wind-induced turbulence conditions through high-frequency measurementsinsitu. Hence, the aims of this study were: (1) to clarify the influence of wind-induced turbulence on colony size ofMicrocystisin the field; (2) to intensively investigate morphology and morphological changes inMicrocystiscolonies; (3) to provide a systematic data set that can be used in subsequent numerical simulations ofMicrocystisbloom formation. Therefore, an enclosure was set up in Lake Taihu and the biomass, colony size and morphology ofMicrocystisat different depths in the enclosure were sampled at 2-h intervals to clarify the effect of turbulence on these parameters and the occurrence ofMicrocystisblooms.

    1 Materials and Methods

    1.1 Site description

    Our investigation was conducted in a 4.00 m × 4.00 m enclosure, which was established approximately 150 m away from the eastern shoreline of Meiliang Bay, Lake Taihu (31°25′N, 120°13′E; Fig.1). The enclosure was made of a flexible geotextile underwater with floating rubber tubing above the surface of the water column, and was open to the atmosphere and to the bottom sediment. The mean water depth in the enclosure during the investigation was approximately 2.00 m and there was no exchange of phytoplankton between the inside and outside of the enclosure.

    圖1 太湖梅梁灣處圍隔的位置及研究期間圍隔內(nèi)水華生消過程

    1.2 Field sampling

    Sampling was performed every 2 h from 14:00 on August 26 to 12:00 on September 7 in 2012. At each sampling, water samples (50 mL) were collected at varying depths through a slender plastic pipe (diameter 0.004 m) linked to a syringe. Water samples were taken from the surface (0 m) and at 0.10 m, 0.25 m, 0.40 m, 0.80 m, 1.20 m, 1.60 m, and 2.00 m depth. Formalin (2% (v/v)) was immediately added to the samples before measurements ofMicrocystiscell density, colony size, and morphology. Meanwhile, equal volumes of water samples were collected at the surface and middle layers of the water column, mixed, and then added into a 500-mL plastic bottle for subsequent nutrient concentration analysis.

    1.3 Measurements of environmental factors

    Air temperature, light intensity, and wind speed were measured in the field using an electronic thermometer (Mettler SG7, Toledo, OH, USA), a digital lux meter (ZDS-10, Shanghai Xuelian Instruments, China), and a wind recorder (TPJ-30, Zhejiang TOP Instrument Co., Ltd., China), respectively. The wind-driven current in the enclosure was measured using a ship-mounted Acoustic Doppler Current Profiler (ADCP; FlowQuest 2000, San Diego, USA) with its beams downward. The ADCP configuration was set to 0.13 m in bin cells (standard deviation of 0.001 m/s). The 3-dimensional current velocities and directions at different water layers were monitored synchronously every second. However, due to the unreceived reflection of 0.27 m near the ADCP[17-18], the 3-dimensional currents in the 0-0.27 m blind zone were unable to be measured in this study.

    The root mean square (RMS) velocity (m/s) was used to define the characteristic speed of the turbulence in each water layer, and was calculated as follows:

    (1)

    (2)

    whereμRMSxis the fluctuation of the current for Cartesian vectorx(which is similarly calculated for theyandzvectors) andnis the number of samples per measurement. TheRMSvelocities are expressed as averages in different water layers. The turbulent dissipation rate (ε in m2/s3) in each water layer, which describes the turbulence intensity, was calculated from theRMSvelocity (m/s)[19]as follows:

    (3)

    whereAis a dimensional constant of order 1[20]andhis the water depth (m) describing the size of the largest vortices.

    1.4 Analysis of cell density, colony size and morphospecies of Microcystis

    To estimate the cell density in eachMicrocystissample, a 10-mL centrifuge tube containing 5 mL sample was shaken in a water bath oscillator (100℃, 180 r/min) for approximately 5 min to disperse the colonies completely into single cells[23]. Next, cells were counted at least three times in a blood cell counting chamber under an optical microscope (Olympus CX31; Olympus Corp., Japan) at ×400 magnification. The average value of three results was used as the cell density until the difference between three calculated results was less than 10%.

    EachMicrocystissample was shaken well and then photos were taken using a digital camera (Olympus C-5050) coupled to an optical microscope (Olympus CX31). The photomicrographs were analyzed using UTHSCSA ImageTool v3.00 software[24]. TheMicrocystiscolonies were classified into five morphologies:M.ichthyoblabe,M.aeruginosa, sphericalM.wesenbergii, irregularM.wesenbergii, and others, according to the taxonomic methods of Yu et al.[25]and Li et al.[7]. Each individual colony was assumed to be spherical to calculate the biovolume, because it is difficult to accurately measure and calculate the diameter ofMicrocystiscolonies, especially those with irregular morphologies. The diameter (D) ofMicrocystiscolonies was calculated as follows[26]:

    (4)

    where,Lis the longest axis andWis the shortest axis (aligned perpendicular to the longest axis) ofMicrocystiscolonies. More than 200 colonies per sample were measured to determine the biovolume percentage of various morphospecies for each size group of each morphospecies. The median colony diameter (D50) value was used to estimate the average colony size in all measured samples, which indicated that 50% of the colonies were smaller than this size[27-28].

    1.5 Data analysis

    In the current study, the time interval was 24 hours (from noon-12:00 hours- on the first day until noon the next day) and was defined as a complete day. Thus, the sampling period was divided into August 26th, August 27th... September 6th.

    The following formula was used to calculate the average cell density ofMicrocystis(Cave) in the water column:

    (5)

    where,Ciis the cell density ofMicrocystisat the depthi, andhiis the height of the water column at depthi(h1=0.010 m;h2=0.165 m;h3=0.150 m;h4=0.275 m;h5=h6=h7=0.400 m;h8=0.200 m), andhis the water depth (2.00 m).

    There was a certain deviation in each measurement of cell density ofMicrocystis. If only twoCavewere used to calculate theinsitugrowth rate, there will be a large systematic error. Hence, a least square linear regression of the natural logarithm ofCavevs. time was performed for each complete day, and the slopes of the regression lines represented theinsitugrowth rates ofMicrocystis[29].

    The average colony size ofMicrocystis(D50ave) in the water column was calculated as follows:

    (6)

    where,D50iis the average colony size ofMicrocystisat depthi.

    The proportion of each morphospecies (M.ichthyoblabe,M.aeruginosa, sphericalM.wesenbergii, irregularM.wesenbergiiand others) out of all theMicrocystisthroughout water column (Pmorphospecies) was calculated as follows:

    (7)

    where,Piis the proportion of each morphospecies out of allMicrocystisat depthi.

    Correlation analyses between theinsitugrowth rate ofMicrocystisand environmental factors were performed using SPSS 19.0 software (IBM, Armonk, NY, USA). The relationships among the proportions of each morphospecies out of allMicrocystisthroughout the water column were analyzed via the interval maxima regression (IMR) method[30]. The proportions of two morphospecies were defined as the independent variable and the dependent variable, respectively. First, the independent variables along with the dependent variables were arranged in ascending order. Each independent variable was divided into equal increments, resulting in 18-32 intervals. Next, the maximum dependent variables of every interval were obtained and then fitted linearly with each independent variable. Then, to test the correlation between the fitted data and the original data, a bivariate Pearson’s correlation analysis was conducted using SPSS 19.0 software. The relationships between average colony size throughout the water column and wind speed and turbulence dissipation rate were determined by the IMR method. For all analyses, aPvalue of <0.05 was considered to indicate statistical significance.

    The minimum-scale of turbulence is determined by the Kolmogorov scale (LK: μm):

    (8)

    (9)

    where,υis kinematic viscosity (0.901×10-6m2/s at 25℃),γis shear rate (s-1),εis turbulent dissipation rate (m2/s3), which is calculated as follows[31]:

    (10)

    where,wis wind speed (m/s), andhis water depth (h=2.00 m in this study).

    2 Results

    2.1 Meteorological factors and nutrient concentrations

    The air temperature showed a distinct diurnal trend of increasing in the morning and decreasing in the afternoon (Fig.2a). The range of maximum daily air temperature was 24.9-39.1℃, and the average air temperature during the whole study period was 26.9℃. Light intensity showed a similar diurnal trend as air temperature, and the maximum daily light intensity ranged from 178 to 1023 μmol·photons/(m2·s).

    圖2 研究期間氣象因子的晝夜變化

    The average wind speed was 4.87 m/s from 12:00 on August 27thto 12:00 on August 28th, with a peak wind speed of 6.55 m/s(Fig.2b). Wind speed above 3.00 m/s could also be observed on September 1stand 4th. The mean turbulence dissipation rate (ε) ranged from 2×10-5to 19×10-5m2/s3and showed the same trend as the wind speed.

    圖3 研究期間營養(yǎng)鹽的變化規(guī)律

    2.2 Cell density, colony size and morphology of Microcystis

    Microcystiswas the dominant species, accounting for >99% of all phytoplankton during the study. The surface cell density ofMicrocystisfluctuated widely during the study period (Fig.4a).Microcystismainly accumulated in the water surface with a high cell density above 1×107cells/mL on August 26thand 30thand September 3rd. When wind speed was higher than 3.00 m/s on August 28thand September 1stand 4th, the surface cell density decreased to <2×105cells/mL. The average cell density in the whole water column varied from 3×104to 18×104cells/mL during the study (Fig.4a).

    The median colony diameter (D50) ofMicrocystisat the water surface decreased dramatically under strong wind on August 27th(Fig.4b). After August 28th, theD50ofMicrocystisat the surface first increased and then declined. The averageD50ofMicrocystisthrough the whole water column was slightly lower than that at the water surface with range of 66.2-768.0 μm. The trend in the variation inD50in the whole water column also decreased under strong wind on August 27th(Fig.4b). When wind got down, theD50in the whole water column increased and then declined.

    圖4 研究期間微囊藻細(xì)胞密度(a)及群體粒徑(b)的晝夜變化

    圖5 研究期間表層及整水柱不同形態(tài)微囊藻所占比例的晝夜變化

    The proportion of each morphospecies out of allMicrocystisat the water surface and in the water column varied over time (Fig.5). SphericalM.wesenbergiiand irregularM.wesenbergiiwere initially dominant in the surface layer, but the proportions ofM.aeruginosaandM.ichthyoblabeincreased over time. The trends in the variations of the proportions ofMicrocystismorphospecies throughout water column can be summarized as follows: the dominant species wasM.ichthyoblabein the early stage of the study, followed by an increase in the proportions of irregularM.wesenbergiiandM.aeruginosa, and finally, dominance ofM.ichthyoblabe. The proportion of sphericalM.wesenbergiiin the water column gradually decreased during the study.

    2.3 Relationships between Microcystis growth and environmental factors

    Theinsitugrowth rates ofMicrocystisobtained by linear fitting are shown in Fig.6. Except for the growth rate fitting results on August 31st, September 1st, and September 3rd, the other linear fitting results were significant. Theinsitugrowth rate ofMicrocystisranged from -1.074 to 1.134 day-1.

    圖6 研究期間微囊藻原位生長率的變化規(guī)律

    表1 研究期間微囊藻原位生長率與環(huán)境因子和表層藻密度的皮爾遜相關(guān)系數(shù)

    AT, air temperature; L, light intensity; W, wind speed;Csur, surface cell density ofMicrocystis;Cave, average cell density ofMicrocystisthrough water column; T, turbidity.*Significant at 0.05 level (two-tailed).**Significant at 0.01 level (two-tailed).

    2.4 Changes in colonial morphology of Microcystis

    Fig.7 shows the correlations among the proportions of each morphospecies throughout the water column. The proportion ofM.ichthyoblabeout of allMicrocystisthroughout the water column was negatively (P< 0.001) correlated with that ofM.wesenbergii. The proportion ofM.aeruginosawas negatively (P< 0.001) correlated with that of sphericalM.wesenbergii, irregularM.wesenbergii, and the sum of those two morphospecies. There were significant negative correlations between the averageD50and wind speed and meanεin the water column (Fig.8).

    圖7 不同形態(tài)微囊藻所占比例之間的關(guān)系(藍(lán)色點代表最大值,灰色點代表除最大值之外其余的值)

    圖8 整水柱微囊藻群體平均粒徑與風(fēng)速(a)和平均紊流耗散(b)之間的關(guān)系(深灰色點代表最大值,灰色點代表除最大值之外其余的值)

    3 Discussion

    The short-term dynamics ofMicrocystisbiomass in an enclosure ecosystem in a shallow lake was investigated in this study, with theinsitugrowth rate ofMicrocystisanalyzed by curve-fitting method. The horizontal migration ofMicrocystiswas avoided because of the enclosure ecosystem. TheMicrocystisbiomass was not affected by fish or zooplankton because fish were excluded from the enclosure and theMicrocystiscolonies in the enclosure ecosystem were too large (D50> 200 μm (Fig.4b)) to be predated on by zooplankton[32]. Nevertheless, the cell density ofMicrocystisin the 1-cm surface water layer was hundreds of times higher than the average cell density in the rest of the water column. Therefore, an accurate assessment of the thickness ofMicrocystisaggregation on the water surface is an important factor affecting the analytical results, but is almost impossible to achieve. In the current study, surface thicknesses of 0.5 cm, 0.4 cm, 0.3 cm and 0.2 cm were used in calculations. We found that there was no significant difference in the average cell density ofMicrocystisin the whole water column calculated with the surface thickness of 1 cm (ttest,P>0.05, Attached Tab.Ⅰ ). Therefore, theinsitugrowth rate calculated with the surface thickness of 1 cm was deemed reliable.

    The majority of the calculatedinsitugrowth rates in the current study were within the range of 0.069 and 0.642 day-1(Fig.6), similar to the results measured by the FDC method in other studies[33-34]. Theinsitugrowth rate ofMicrocystisdetermined by Cao et al.[35]using specially designed chambers in Lake Taihu was far lower than that determined in our study. This is because they conducted their experiment in spring, when the temperature was much lower than that in the current study (conducted in August). Theinsitugrowth rate ofMicrocystismeasured by the FDC method by Wu and Kong[36]was between 0.2 and 0.3 day-1in August. However, if the FDC was substituted into formula 1 proposed by Tsujimura[33], the calculatedinsitugrowth rate could reach 0.6 day-1, a value similar to those obtained in this study. Similarly, Wilson et al.[37]reported that the growth rate of individualMicrocystiscolonies ranged from 0.2 to 0.4 day-1. Li et al.[38]found that the maximuminsitugrowth rate ofMicrocystisin Lake Taihu measured by the RNA/TOC method was 0.6 day-1. Although Stolte and Garcés[39]noted that most cyanobacteria have aninsitugrowth rate of >0.5 day-1, and many show growth rates of >1.0 day-1in the laboratory, theinsitugrowth rate ofMicrocystisin a shallow lake would be <0.6 day-1most of the time. However, some of theinsitugrowth rates calculated in the present study were between -0.3 and -0.5 day-1(Fig.6), indicating that there was a considerable decline in theMicrocystiseven in midsummer. Further studies on the mechanism ofMicrocystisdecline are expected to provide a new ideas for the control ofMicrocystisblooms.

    In this current study, the dissolved phosphorus in the water column may have been absorbed by the rapidly growingMicrocystis, whose growth was promoted by increasing turbulence. Xie et al.[46]reported a similar result where the TDP concentration in a shallow lake (Lake Donghu, China) decreased due toMicrocystisgrowth. The large number ofMicrocystisaccumulated on the water surface will be under a state of nutrient stress[47-48]and photoinhibition[49-50], which are not conducive to its growth. Sufficient turbulence can causeMicrocystisto become evenly distributed in the water column, rather than aggregated at the surface[51-52]. The homodispersion ofMicrocystiscolonies in the water column is conducive to light transmission and nutrient uptake, both of which stimulateMicrocystisgrowth[53-54]. The existence of this pathway was verified by the negative relationship between theinsitugrowth rate and the surface cell density ofMicrocystisdetected in the current study (Tab.1). The nitrogen concentration was relatively high in our study, and was not significantly related to the growth rate ofMicrocystis.

    This study showed a significant negative relationship betweenMicrocystiscolony size and wind speed during the 12 days of continuous monitoring, indicating that wind-induced turbulence in shallow lakes could break upMicrocystiscolonies. A similar negative correlation has been demonstrated in the indoor oscillating grids test carried out withMicrocystiscolonies by O’Brien et al.[11]. Zhang et al.[55]also found the length ofDolichospermumflos-aquaefilaments decreased noticeably with the increasing turbulence energy dissipation rates in chemostats with self-designed mixing propellers. The turbulence dissipation rate in the current study ranged from 2×10-5to 19×10-5m2/s3, indicating that turbulence at the above intensities could break upMicrocystiscolonies under natural conditions. According to the relationship between turbulence intensity and wind speed as well as water depth given by MacKenzie and Leggett[31], a turbulence intensity of 2.33×10-5m2/s3is equivalent to that produced by an average wind speed of 5.00 m/s in a 30-m depth reservoir. Li et al.[7]showed that the minimum turbulence intensity required to break aMicrocystiscolony was 0.020 m2/s3in a laboratory experiment using submerged impellers. The minimum turbulence intensity in the field would be much lower than that in a laboratory because the laboratory experimental time (30 min) was too short to observe disaggregation ofMicrocystisat lower turbulence intensities. Therefore, the disaggregation ofMicrocystiscolonies by wind-induced turbulence will occur continuously under natural conditions.

    The process of colony formation by cell-adhesion is faster than that by cell-division[56]. Although the extracellular polysaccharides on the surface ofMicrocystiscolonies and cells are negatively charged[57-58], Chen and Lürling[59]showed that Ca2+, an important divalent electrolyte cation in the water, will counteract the electrostatic repulsion to promote the formation of largeMicrocystiscolonies by cell adhesion. Moreover, Qin et al[8]observed that the colony size ofMicrocystisincreased from 32.8 to 69.4 μm under the action of wind-waves, suggesting that strong turbulence could increase the chance of collisions betweenMicrocystiscolonies and cells. However, our field results did not show that strong turbulence promoted the adhesion ofMicrocystiscells or colonies to form larger colonies (Fig.8). This may be affected by the minimum-scale of turbulence. As shown in Fig.9, the maximum ofD50avewas limited by theLKunder different wind speeds. The opposite results between Qin et al.[12]and our study was because the colony size ofMicrocystiswas smaller thanLKduring their monitoring period. During August 28thto 29thand September 4thto 5thwhen the wind speed was less than 1.00 m/s,Microcystiscolonies aggregated at the water surface (Fig.4a), resulting in significantly increasedD50ave(Fig.4b). Other studies also found that individual colony size in the surface scums can be up to 2000 μm, but theD50is usually less than 800 μm[4,9,26]. Therefore, low turbulence would promoteMicrocystisadhesion to form large colonies, but the colony sizes does not exceed the minimum-scale of turbulence.

    圖9 整水柱微囊藻群體平均粒徑與柯爾莫哥洛夫尺度之間的關(guān)系

    No significant correlation was detected between the proportions of sphericalM.wesenbergiiand irregularM.wesenbergii(Fig.7c), although Li et al.[7]have shown that sphericalM.wesenbergiican change into irregularM.wesenbergiiunder the effect of turbulence. The lack of correlation may be caused byM.ichthyoblabecolonies continuously changing into sphericalM.wesenbergiicolonies, thus interfering with the relationship between sphericalM.wesenbergiiand irregularM.wesenbergii. However, both the proportions of sphericalM.wesenbergiiand irregularM.wesenbergiiwere negatively correlated with that ofM.aeruginosa(Fig.7d, e), indicating that whetherM.wesenbergiiwas deformed or not, their colonies were continuously dissolved and turned intoM.aeruginosa. The relationships found in the current study between the variousMicrocystiscolony morphotypes (Fig.7) support the morphology change model proposed by Xiao et al.[13]that the morphology ofMicrocystiscolonies changes over time. Other researchers[60-62], however, have shown that there are some differences in cell size, mucilage sheath, toxin production, and gene sequences in different morphotypes ofMicrocystiscolonies. In the future studies, it is necessary to measure the above characterization in the dynamic change ofMicrocystismorphotypes for verifying the hypothesis.

    Recent studies have forecasted that global climate change will reduce average wind speeds on land[63-65], resulting in lower turbulence intensity in lakes. The results of this study imply that the colony size ofMicrocystisin lakes will increase with lower wind speeds, and the probability of surface accumulation will increase. Another prediction from global climate change is an increase in extreme weather events, such as tropical typhoons with high rainfall and winds, which may not only transport nutrients into lakes, but will also stir up sediments[66-68]. In our study, no relation between water column nutrient concentrations and wind speed was found, but at higher winds, such as during typhoons, dissolved inorganic nutrients at the sediment-water interface increased with wind[66]. Consequently, cyanobacterial biomass increased in the next month after tropical cyclones passed[66]. Hence, the availability of nutrients and a more favorable volume-to-surface ratio of smaller colonies due to colony breakage under high winds will promote biomass build-up, while during subsequent stable weather conditions likely larger sizedMicrocystiscolonies will accumulate at the water surface.

    The aggregation ofMicrocystisat the water surface is not conducive to light penetration. Thus, other phytoplankton may receive even less light[13, 65]. Janatian et al.[69]indicated that a decrease in wind speed will shift phytoplankton groups fromr-selected coccal and colonial green algae and cyanobacteria toK-selected shade-adapted thin filamentous cyanobacteria in the large and shallow Lake V?rtsj?rv. Therefore, decreasing wind speed will likely reduce phytoplankton biodiversity in eutrophic lakes, and buoyant cyanobacteria such asMicrocystiswill become more dominant.

    4 Conclusions

    Strong turbulence in field conditions promotedMicrocystisgrowth due to homodispersion ofMicrocystiscolonies and release of nutrient from sediment.Microcystiscolonies can be disaggregated by strong wind-induced turbulence in the field. The negative relationships betweenM.ichthyoblabeandM.wesenbergiiand betweenM.wesenbergiiandM.aeruginosain the proportion of morphospecies out of allMicrocystiswere found in this study, but more evidence is need to support the hypothesis that the morphology ofMicrocystiscolonies changes over time. On the basis of our results, we conclude that a climatically modulated decline in wind speed will promote the surface accumulation ofMicrocystisdue to larger colony formation.

    5 Appendix

    The Attached Table Ⅰ is seen in the electronic version (DOI: 10.18307/2021.0205).

    猜你喜歡
    水柱生長率微囊
    探探鯨的水柱
    日本落葉松以林分平均木材積生長率代替林分蓄積生長率可行性驗證
    Run through the rain
    海南省馬占相思生長率模型研建
    水柱有“魔力”
    水柱測量中的水下滑翔機(jī)轉(zhuǎn)向性能
    基于海南省連清資料建立松樹和橡膠樹材積生長率模型
    微囊懸浮-懸浮劑和微囊懸浮劑不是同種劑型
    海洋浮游纖毛蟲生長率研究進(jìn)展
    微囊藻毒素-LR對秀麗線蟲精子形成的毒性作用
    狂野欧美白嫩少妇大欣赏| 亚洲 国产 在线| 少妇丰满av| 观看美女的网站| 很黄的视频免费| 欧美最黄视频在线播放免费| 国产一区二区三区视频了| 久久婷婷人人爽人人干人人爱| 两个人视频免费观看高清| 美女高潮的动态| 免费av观看视频| 午夜福利高清视频| 九色成人免费人妻av| 亚洲av免费高清在线观看| 午夜久久久久精精品| 老熟妇乱子伦视频在线观看| 久久亚洲精品不卡| 国产成人影院久久av| 天堂网av新在线| 一个人看视频在线观看www免费| 日日啪夜夜撸| 一个人免费在线观看电影| 午夜激情欧美在线| 欧美日韩国产亚洲二区| 国产不卡一卡二| 一个人免费在线观看电影| 国产亚洲欧美98| 欧美日本视频| 少妇人妻一区二区三区视频| 观看美女的网站| 啦啦啦啦在线视频资源| 国产精品国产高清国产av| 少妇人妻一区二区三区视频| 日本 欧美在线| 欧洲精品卡2卡3卡4卡5卡区| 成人特级黄色片久久久久久久| 久久久久久久久大av| 午夜福利在线观看吧| 欧美极品一区二区三区四区| 国语自产精品视频在线第100页| 黄片wwwwww| 亚洲欧美日韩卡通动漫| 国产又黄又爽又无遮挡在线| 国内少妇人妻偷人精品xxx网站| 深夜精品福利| 一级黄片播放器| 亚洲av中文字字幕乱码综合| 精品久久久久久,| 在线观看美女被高潮喷水网站| 熟女人妻精品中文字幕| 国产伦一二天堂av在线观看| 国产不卡一卡二| 成人午夜高清在线视频| 婷婷丁香在线五月| 亚洲中文字幕日韩| 琪琪午夜伦伦电影理论片6080| 乱系列少妇在线播放| 天堂av国产一区二区熟女人妻| 国产乱人视频| 一级毛片久久久久久久久女| 国产精品永久免费网站| 国产午夜福利久久久久久| 美女高潮喷水抽搐中文字幕| 国产真实乱freesex| 可以在线观看毛片的网站| 在线观看av片永久免费下载| 嫩草影院新地址| 国产v大片淫在线免费观看| 国内精品一区二区在线观看| 99久久中文字幕三级久久日本| 亚洲va在线va天堂va国产| 久久国产精品人妻蜜桃| 国内少妇人妻偷人精品xxx网站| 久久久久久国产a免费观看| 狂野欧美白嫩少妇大欣赏| 午夜日韩欧美国产| 一级a爱片免费观看的视频| 成人午夜高清在线视频| 国产在线精品亚洲第一网站| 神马国产精品三级电影在线观看| 久久久久久久精品吃奶| 亚洲不卡免费看| 欧美日韩综合久久久久久 | 99在线视频只有这里精品首页| 亚洲三级黄色毛片| 无遮挡黄片免费观看| 免费大片18禁| 久久精品国产亚洲av天美| 欧美区成人在线视频| 亚洲自偷自拍三级| 最新在线观看一区二区三区| 又紧又爽又黄一区二区| 久久精品影院6| 欧美xxxx性猛交bbbb| 午夜视频国产福利| 午夜影院日韩av| 一区二区三区四区激情视频 | 啦啦啦啦在线视频资源| 一级黄色大片毛片| 国产男人的电影天堂91| 桃色一区二区三区在线观看| 婷婷精品国产亚洲av| 色吧在线观看| 99热这里只有是精品在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国内揄拍国产精品人妻在线| 别揉我奶头~嗯~啊~动态视频| 可以在线观看毛片的网站| 国产亚洲精品久久久久久毛片| 97碰自拍视频| 国产一区二区三区av在线 | xxxwww97欧美| 九九爱精品视频在线观看| 欧美另类亚洲清纯唯美| 蜜桃久久精品国产亚洲av| 精品久久久噜噜| 国产综合懂色| 又爽又黄a免费视频| 黄色欧美视频在线观看| netflix在线观看网站| 人妻制服诱惑在线中文字幕| 国产视频一区二区在线看| 亚洲欧美日韩无卡精品| 欧美黑人巨大hd| 国产综合懂色| 午夜福利18| 天天躁日日操中文字幕| 真人做人爱边吃奶动态| 成人国产麻豆网| 午夜福利视频1000在线观看| 中文字幕高清在线视频| 国产午夜精品久久久久久一区二区三区 | 久久这里只有精品中国| 狂野欧美白嫩少妇大欣赏| 国产精品99久久久久久久久| 欧美黑人欧美精品刺激| 91精品国产九色| 欧美区成人在线视频| 国产一区二区在线av高清观看| 麻豆av噜噜一区二区三区| 亚洲色图av天堂| 天堂动漫精品| 夜夜夜夜夜久久久久| 欧美激情久久久久久爽电影| 国产av不卡久久| 两人在一起打扑克的视频| 国产精品一区二区三区四区免费观看 | 国产主播在线观看一区二区| 女生性感内裤真人,穿戴方法视频| 亚洲人成伊人成综合网2020| 熟女电影av网| 国内毛片毛片毛片毛片毛片| 在线播放国产精品三级| 日韩欧美国产一区二区入口| 色精品久久人妻99蜜桃| 国内精品久久久久精免费| 99国产极品粉嫩在线观看| 色播亚洲综合网| 亚洲欧美日韩高清在线视频| 国产男靠女视频免费网站| 免费搜索国产男女视频| 国内精品久久久久精免费| 欧美日本亚洲视频在线播放| 俺也久久电影网| 中文字幕av在线有码专区| 久久精品国产自在天天线| 毛片女人毛片| 亚洲av五月六月丁香网| 久久久精品欧美日韩精品| 中文字幕免费在线视频6| 夜夜看夜夜爽夜夜摸| 国产精品爽爽va在线观看网站| 国产v大片淫在线免费观看| 亚洲中文字幕日韩| 91在线精品国自产拍蜜月| 十八禁国产超污无遮挡网站| 国产一区二区在线观看日韩| 一级黄片播放器| 网址你懂的国产日韩在线| 哪里可以看免费的av片| 成人亚洲精品av一区二区| 欧美三级亚洲精品| 日韩欧美免费精品| 国产精品野战在线观看| 一进一出好大好爽视频| 女人被狂操c到高潮| 亚洲专区国产一区二区| 91久久精品国产一区二区三区| 亚洲avbb在线观看| 99国产极品粉嫩在线观看| 最近最新免费中文字幕在线| 亚洲成av人片在线播放无| 一区福利在线观看| 美女高潮的动态| 在线天堂最新版资源| 舔av片在线| 国产在视频线在精品| 日韩人妻高清精品专区| 禁无遮挡网站| 免费看光身美女| 亚洲国产欧洲综合997久久,| 特大巨黑吊av在线直播| 日韩一本色道免费dvd| 夜夜爽天天搞| 91狼人影院| 国产91精品成人一区二区三区| 在线看三级毛片| 国产亚洲精品久久久久久毛片| or卡值多少钱| 淫妇啪啪啪对白视频| 亚洲精品日韩av片在线观看| 中文字幕精品亚洲无线码一区| 波野结衣二区三区在线| 日本a在线网址| 精品国内亚洲2022精品成人| 国产免费av片在线观看野外av| 无人区码免费观看不卡| 亚洲一区二区三区色噜噜| 超碰av人人做人人爽久久| 一本精品99久久精品77| 日韩中字成人| 麻豆精品久久久久久蜜桃| 精品福利观看| 一级毛片久久久久久久久女| 小蜜桃在线观看免费完整版高清| 日日撸夜夜添| 特大巨黑吊av在线直播| 日日啪夜夜撸| 人妻少妇偷人精品九色| 全区人妻精品视频| 亚洲精品影视一区二区三区av| 18+在线观看网站| 国产精品免费一区二区三区在线| 欧美日韩亚洲国产一区二区在线观看| 自拍偷自拍亚洲精品老妇| 99久久精品国产国产毛片| 两性午夜刺激爽爽歪歪视频在线观看| 欧美性猛交黑人性爽| 中文字幕人妻熟人妻熟丝袜美| 久久人人精品亚洲av| 午夜福利视频1000在线观看| 色综合色国产| 日韩强制内射视频| 99热这里只有是精品50| 日韩精品有码人妻一区| 成年免费大片在线观看| 老熟妇仑乱视频hdxx| 午夜福利在线观看吧| 尾随美女入室| 综合色av麻豆| 国产av麻豆久久久久久久| 丰满的人妻完整版| 精品久久久久久成人av| 18禁黄网站禁片免费观看直播| 亚洲欧美日韩高清专用| 国产老妇女一区| 舔av片在线| 99久国产av精品| 日本黄色视频三级网站网址| 国产精品一区二区三区四区久久| 欧美zozozo另类| 麻豆av噜噜一区二区三区| 亚洲色图av天堂| 麻豆久久精品国产亚洲av| 日韩欧美三级三区| 亚洲国产色片| 18+在线观看网站| 婷婷色综合大香蕉| 91久久精品国产一区二区成人| 真实男女啪啪啪动态图| 日本-黄色视频高清免费观看| 国产男人的电影天堂91| 久99久视频精品免费| a级一级毛片免费在线观看| 亚洲真实伦在线观看| 精品免费久久久久久久清纯| 国产日本99.免费观看| 少妇熟女aⅴ在线视频| 午夜a级毛片| 国产亚洲精品久久久com| 国产在线精品亚洲第一网站| 两个人的视频大全免费| 亚洲一级一片aⅴ在线观看| 亚洲最大成人手机在线| 乱码一卡2卡4卡精品| 欧美激情国产日韩精品一区| 国产午夜精品久久久久久一区二区三区 | 久久99热6这里只有精品| 成人一区二区视频在线观看| 日韩一本色道免费dvd| 听说在线观看完整版免费高清| 中出人妻视频一区二区| 亚洲欧美日韩东京热| 国产精品福利在线免费观看| 偷拍熟女少妇极品色| 日本撒尿小便嘘嘘汇集6| 少妇的逼水好多| 日本黄色视频三级网站网址| 日韩精品有码人妻一区| 听说在线观看完整版免费高清| 九九爱精品视频在线观看| 午夜福利成人在线免费观看| 日日摸夜夜添夜夜添小说| 亚洲欧美日韩东京热| 国产午夜精品久久久久久一区二区三区 | 国产大屁股一区二区在线视频| 小说图片视频综合网站| 有码 亚洲区| 亚洲人成伊人成综合网2020| 联通29元200g的流量卡| av福利片在线观看| 日韩中字成人| 亚洲av一区综合| 国产v大片淫在线免费观看| 老熟妇乱子伦视频在线观看| 男人舔奶头视频| 国产精品久久久久久精品电影| 校园人妻丝袜中文字幕| 99热这里只有是精品在线观看| 亚洲av中文字字幕乱码综合| 成人av一区二区三区在线看| 欧美精品国产亚洲| 精品欧美国产一区二区三| 波多野结衣巨乳人妻| 国产精品无大码| 嫩草影院精品99| 高清毛片免费观看视频网站| 国产v大片淫在线免费观看| 麻豆一二三区av精品| 国产极品精品免费视频能看的| 91午夜精品亚洲一区二区三区 | 99久久精品国产国产毛片| 亚洲男人的天堂狠狠| 成人特级黄色片久久久久久久| 日本免费a在线| av在线老鸭窝| 成人毛片a级毛片在线播放| 日本五十路高清| 国产一区二区亚洲精品在线观看| 日韩高清综合在线| 日本在线视频免费播放| 老司机深夜福利视频在线观看| 国产乱人视频| 亚洲内射少妇av| 俺也久久电影网| 一卡2卡三卡四卡精品乱码亚洲| av福利片在线观看| 性插视频无遮挡在线免费观看| 欧美激情久久久久久爽电影| 少妇猛男粗大的猛烈进出视频 | 日日摸夜夜添夜夜添小说| 国产国拍精品亚洲av在线观看| 天堂网av新在线| 五月玫瑰六月丁香| 日日夜夜操网爽| 亚洲内射少妇av| 日日夜夜操网爽| 床上黄色一级片| 窝窝影院91人妻| 国产精品人妻久久久影院| 最新在线观看一区二区三区| 日日摸夜夜添夜夜添av毛片 | 97热精品久久久久久| 欧美中文日本在线观看视频| 日日摸夜夜添夜夜添小说| 国产熟女欧美一区二区| 国产精品av视频在线免费观看| 男女边吃奶边做爰视频| 久久久久国内视频| 久久久久久久久久久丰满 | 少妇猛男粗大的猛烈进出视频 | 搞女人的毛片| 精品免费久久久久久久清纯| 亚洲av五月六月丁香网| 夜夜爽天天搞| 亚洲中文字幕日韩| 亚洲国产色片| 午夜福利在线在线| 成年女人毛片免费观看观看9| 一级a爱片免费观看的视频| 亚洲欧美精品综合久久99| 欧美日韩精品成人综合77777| 非洲黑人性xxxx精品又粗又长| 伦精品一区二区三区| eeuss影院久久| 最近最新免费中文字幕在线| 久久久精品大字幕| 欧美最黄视频在线播放免费| 伊人久久精品亚洲午夜| 国产精品久久久久久亚洲av鲁大| 亚洲七黄色美女视频| 如何舔出高潮| 小蜜桃在线观看免费完整版高清| 久久久久久九九精品二区国产| 99在线人妻在线中文字幕| 亚洲aⅴ乱码一区二区在线播放| 99久国产av精品| 国产黄片美女视频| 久久午夜福利片| 亚洲人成网站在线播放欧美日韩| 男人和女人高潮做爰伦理| 国产极品精品免费视频能看的| 国产欧美日韩一区二区精品| 此物有八面人人有两片| 国产精品伦人一区二区| 九九久久精品国产亚洲av麻豆| 午夜久久久久精精品| 亚洲自偷自拍三级| 看黄色毛片网站| 国产精品女同一区二区软件 | 亚洲七黄色美女视频| 欧美不卡视频在线免费观看| a在线观看视频网站| 日本-黄色视频高清免费观看| 美女cb高潮喷水在线观看| 人妻夜夜爽99麻豆av| 18禁黄网站禁片午夜丰满| 免费看a级黄色片| 极品教师在线视频| 精品一区二区免费观看| 联通29元200g的流量卡| 俄罗斯特黄特色一大片| 在线国产一区二区在线| 亚洲av免费在线观看| 18禁黄网站禁片午夜丰满| 尾随美女入室| 美女大奶头视频| 男女之事视频高清在线观看| 久久草成人影院| 永久网站在线| 舔av片在线| 极品教师在线视频| 免费av不卡在线播放| 久久人人爽人人爽人人片va| 亚洲人成伊人成综合网2020| 白带黄色成豆腐渣| 少妇的逼好多水| 最新在线观看一区二区三区| 在线a可以看的网站| 精品久久久久久久久av| 午夜激情福利司机影院| 国产高清三级在线| 亚洲成人精品中文字幕电影| 午夜福利在线观看吧| 欧美+日韩+精品| 亚州av有码| 干丝袜人妻中文字幕| 午夜福利在线观看免费完整高清在 | 别揉我奶头 嗯啊视频| 久久人人精品亚洲av| 别揉我奶头~嗯~啊~动态视频| 变态另类成人亚洲欧美熟女| 五月玫瑰六月丁香| 国产视频内射| 亚洲色图av天堂| 亚洲国产欧洲综合997久久,| 日本黄色片子视频| 日本-黄色视频高清免费观看| 99九九线精品视频在线观看视频| 国产伦一二天堂av在线观看| 色av中文字幕| 欧美日韩黄片免| 嫩草影院入口| 日韩欧美在线乱码| 亚洲狠狠婷婷综合久久图片| 波多野结衣巨乳人妻| 精品午夜福利在线看| 国产aⅴ精品一区二区三区波| 亚洲人成网站在线播放欧美日韩| 色噜噜av男人的天堂激情| 亚洲人成伊人成综合网2020| 在线观看av片永久免费下载| 国产亚洲精品久久久com| 国产精品98久久久久久宅男小说| or卡值多少钱| 久久人妻av系列| 国产精品无大码| 热99re8久久精品国产| 日本黄色视频三级网站网址| 女人十人毛片免费观看3o分钟| 尾随美女入室| 久久久久久久午夜电影| 精品久久久久久久久久免费视频| 午夜a级毛片| 国产单亲对白刺激| 日日啪夜夜撸| 无人区码免费观看不卡| 成人毛片a级毛片在线播放| 18禁裸乳无遮挡免费网站照片| 国产午夜福利久久久久久| 免费在线观看成人毛片| 黄色女人牲交| 看黄色毛片网站| 日韩av在线大香蕉| 国产高清视频在线播放一区| 成人二区视频| 91麻豆精品激情在线观看国产| 成人精品一区二区免费| 久久亚洲真实| 在线观看av片永久免费下载| 亚洲人成网站在线播| 国产亚洲91精品色在线| 亚洲性久久影院| 性色avwww在线观看| 久久热精品热| 欧美日韩亚洲国产一区二区在线观看| АⅤ资源中文在线天堂| 内射极品少妇av片p| 美女cb高潮喷水在线观看| 国内精品美女久久久久久| 精品欧美国产一区二区三| 悠悠久久av| 别揉我奶头~嗯~啊~动态视频| av福利片在线观看| 亚洲av熟女| 国产探花极品一区二区| 成年女人永久免费观看视频| 日韩欧美国产在线观看| 精品久久久久久,| 在线看三级毛片| 91午夜精品亚洲一区二区三区 | 亚洲在线观看片| 精品人妻视频免费看| 蜜桃亚洲精品一区二区三区| 久久人人精品亚洲av| 欧美丝袜亚洲另类 | 美女 人体艺术 gogo| 日韩欧美国产一区二区入口| x7x7x7水蜜桃| 亚洲欧美日韩无卡精品| 女的被弄到高潮叫床怎么办 | 色尼玛亚洲综合影院| 99久久精品热视频| 日韩欧美在线乱码| 啪啪无遮挡十八禁网站| 91午夜精品亚洲一区二区三区 | 女生性感内裤真人,穿戴方法视频| 国产伦一二天堂av在线观看| 免费观看的影片在线观看| 亚洲成人精品中文字幕电影| 亚洲在线自拍视频| 日日干狠狠操夜夜爽| 99久久久亚洲精品蜜臀av| 欧美另类亚洲清纯唯美| 欧美日韩中文字幕国产精品一区二区三区| 婷婷精品国产亚洲av| 精品无人区乱码1区二区| 免费av毛片视频| 搡老熟女国产l中国老女人| 男插女下体视频免费在线播放| 人人妻,人人澡人人爽秒播| 日本黄色片子视频| 午夜福利欧美成人| 国产久久久一区二区三区| 免费av毛片视频| 又爽又黄a免费视频| 88av欧美| 啦啦啦韩国在线观看视频| 有码 亚洲区| 亚洲无线观看免费| 日本a在线网址| 波野结衣二区三区在线| www.www免费av| 人妻夜夜爽99麻豆av| 日韩 亚洲 欧美在线| 日本与韩国留学比较| 婷婷丁香在线五月| 国产免费av片在线观看野外av| av在线老鸭窝| 亚洲精品粉嫩美女一区| 国产亚洲精品综合一区在线观看| 国产精品自产拍在线观看55亚洲| 小蜜桃在线观看免费完整版高清| 黄色视频,在线免费观看| 免费av观看视频| 俄罗斯特黄特色一大片| 熟妇人妻久久中文字幕3abv| 男女边吃奶边做爰视频| 日日摸夜夜添夜夜添小说| 久久久久久国产a免费观看| 99精品在免费线老司机午夜| 亚洲精品色激情综合| 啦啦啦韩国在线观看视频| 天堂√8在线中文| 免费一级毛片在线播放高清视频| 国产精品三级大全| 真人做人爱边吃奶动态| 在线观看美女被高潮喷水网站| 国产精品国产高清国产av| 午夜福利在线观看吧| 国内精品美女久久久久久| 亚洲人成伊人成综合网2020| 五月伊人婷婷丁香| 亚洲天堂国产精品一区在线| 少妇的逼好多水| 欧美+日韩+精品| 欧美xxxx黑人xx丫x性爽| 国产高潮美女av| 99久久成人亚洲精品观看| 一个人免费在线观看电影| 国内少妇人妻偷人精品xxx网站| 精品人妻一区二区三区麻豆 | av在线观看视频网站免费| 最新在线观看一区二区三区| 美女免费视频网站| 欧美另类亚洲清纯唯美| 国产蜜桃级精品一区二区三区| 尤物成人国产欧美一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 内射极品少妇av片p| 嫩草影院新地址| 国产高清不卡午夜福利| 国产人妻一区二区三区在| 精品久久久久久,| av黄色大香蕉|