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

    Variations of Summer Phytoplankton Community Related to Environmental Factors in a Macro-Tidal Estuarine Embayment, Hangzhou Bay, China

    2015-04-01 01:57:54ZHANGYuexiaYUJunJIANGZhibingWANGQinandWANGHui
    Journal of Ocean University of China 2015年6期

    ZHANG Yuexia, YU Jun, JIANG Zhibing, WANG Qin, and WANG Hui

    ?

    Variations of Summer Phytoplankton Community Related to Environmental Factors in a Macro-Tidal Estuarine Embayment, Hangzhou Bay, China

    ZHANG Yuexia1), 2), YU Jun2), JIANG Zhibing3), WANG Qin2), and WANG Hui4), *

    1),,266100,2),310007,3),,,310012,4),100081,

    To explore the distribution and composition of phytoplankton community and their responses to environmental changes, summer net-collected phytoplankton and physicochemical parameters in the Hangzhou Bay during 2004–2010 were investigated. A total of four phyla and 84 species were identified, including 67 diatom and 12 dinoflagellate species. The dominant species con- stantly consisted of the diatoms, although the dominance of dinoflagellate and cyanobacteria increased recently. Due to great spa- tio-temporal variations in environmental factors (salinity, suspended solids, and nutrient concentration), significant heterogeneities in community compositions among different years and subregions (inner and middle sections, and bay mouth) were found based on the analyses of multidimensional scaling and similarity. Canonical correspondence analysis showed that salinity and Si/N were the main variables associated with algal assemblage. Compared with the historical data since the 1980s, eutrophication (N, P, and N/P in- creased with decreasing Si/N) was exacerbated drastically. Moreover, climatic forcing and human activities resulted in a series of physical alterations, including sediment retention,temperature increase, and salinity decrease as well as reduction in water exchanges. All these changes induced obvious increases in cell density and Chl-while decreases in species diversity and diatom-dinoflagellate ratio as well as the shifting of dominant species. Therefore, the long-term phytoplankton variations were closely related to anthropo- genic and climatic perturbations in the Hangzhou Bay.

    Hangzhou Bay; phytoplankton; long-term changes; eutrophication; Qiantang River; Changjiang

    1 Introduction

    Hangzhou Bay (HZB) is a macro-tidal estuarine em- bayment located to the south of Changjiang Estuary (CE) in the East China Sea. Besides the Changjiang River, the bay receives other riverine loading, including Qiantang, Yong, and Cao’e River. The area (with cities, Shanghai, Hangzhou, Ningbo,) around the HZB is an important populous and economic center in China, and thus the bay has been suffering from intensive human activities, such as large-scale tidal flat reclamation (Ni and Lin, 2003; Zhang., 2010) and rapid industrial and agricultural development. Particularly since the 1980s, various all the nutrients and organic matters, and other pollutants in- flowing into the HZB have drastically increased (Zhang., 2001; Han., 2003; Gao., 2011). Moreover, recent precipitation increase with decreasing evaporation (Wang., 2011, 2013) as well as the successive fre-quent water conservancy projects in the Qiantang River and Changjiang basins (Xu., 2010; Yang., 2011) have significantly altered the freshwater and sediment-associated nutrient inputs (N, P, and Si) into the bay. All these anthropogenic activities and climatic forcing have significantly altered the geographical and physicochemi- cal environment (Xie., 2009, 2013). Consequently, HZB has confronted a series of ecological challenges caused by the anthropogenic and climate-driven physico- chemical changes.

    Phytoplankton is the important primary producer in estuarine/coastal ecosystem that plays a key role in the biogeochemical processes in riverine-seawater interface (Cloern and Dufford, 2005). It is also highly sensitive to environmental destabilizations and thus used in the as- sessment of ecological conditions and estuarine/coastal environmental changes (Zhou., 2008; Jiang., 2014). Therefore, its community variations are always paid great attention by researchers. Previous studies (Zhou., 2008; Wang and Wu, 2009; Li., 2010; Jiang., 2010; Jiang, 2013, 2014) have found signify- cant phytoplankton community alterations (, abundance increase, diatom-dinoflagellate ratio decrease, dominant species shifting, and frequent harmful algal blooms) in the CE and Xiangshan Bay adjacent to the HZB, owing to intensive human- and climate-driven changes. Phyto- plankton community in the HZB is also subjected to in- creased physicochemical stresses. However, distribution and composition of phytoplankton community associated with environmental factors in the HZB are not well documented (Zhang., 2001; Cai, 2006; Wang., 2008; Zhou., 2010). Moreover, the variations in phy- toplankton community response to environmental changes have not been examined.

    The present study focused on the variations of net-col- lected phytoplankton in the HZB during 2004–2010 com- bined with historical data since the 1980s, which allowed investigation into the long-term changes in the phyto- plankton community during the past 30 years. The influ- ences of environmental factors on phytoplankton were also examined. Our objectives are 1) to understand the distribution and composition of the phytoplankton com- munity in relation to environmental parameters in the HZB; 2) to explore the long-term changes in phytoplank- ton assemblage under increasing anthropogenic and cli- matic perturbations; and 3) to supply a scientific founda- tion for environmental protection and further study on this estuarine bay.

    2 Materials and Methods

    2.1 Study Area

    The HZB is a large funnel-shaped and shallow (with an average depth of about 10m) estuarine bay located in the East China Sea, immediately south of the CE and east of the Zhoushan archipelago (Fig.1). The semidiurnal tide is the main driving force behind the horizontal water flow in the HZB, with the M2constituent being the dominant tidal component (ECBCC, 1992). The average tide am- plitude is 3–4 m at the mouth and 4–6 m further upstream. The main river system discharging freshwater and sedi- ment directly into the HZB is Qiantang River, followed by Cao’e and Yong River. The Qiantang River has an average runoff of 163.7×108m3yr?1and an average sedi- ment load of 203.5×104tyr?1recorded at Lanxi station (data from the website: http://www.irtces.org/nishagb_2012. asp). The CE is situated immediately north of the HZB. The average water discharge from Changjiang is 8927×108m3yr?1, with a sediment load of 3.85×108tyr?1re- corded at Datong station (data from the website: http://www.irtces.org/nishagb_2012.asp). A part of water and sediment from the CE enters the HZB and imposes pro- found influences on both hydrodynamics and sedimenta- tion of the HZB (Su and Wang, 1989). The Changjiang Diluted Water (CDW) plays an important role in both the circulation and sediment transport inside the bay (Su and Wang, 1989; Xie., 2013). The plume front serves as a guide for sediment transport and it is closely related to the rapid accretion at the south shore of the bay (Xie., 2009, 2013).

    2.2 Sample Collection and Analysis

    A total of 12 stations (Fig.1) were set in HZB in the summer during 2004–2010. Phytoplankton samples were collected vertically from bottom to surface using a plank- ton net with a mesh size of 76μm (Jiang., 2013). All collection, preservation, and identification of the phyto- plankton samples were conducted according to the speci- fication GB 17378.7-2007. Surface salinity, dissolved in- organic nitrogen (DIN), phosphorus (DIP), silicate (DSi), chemical oxygen demand (COD), and suspended solids (SS) were measured according to the specification GB 17378.4-2007.

    2.3 Data Analysis

    The indices of Shannon diversity () and Pielou evenness () were calculated through the PRIMER 6.0. Dominance () was calculated according to the following equation:f×n/×100%,wherenis the number (cells) of individuals of species,is the total number of indi- viduals, andfis the frequency of speciesoccurring in all the samples. Only those species with≥2% were consid- ered to be the dominants. The eutrophication index () was calculated according to the following equation:=DIN×DIP×COD/4500, where the units of DIN, DIP, and COD are mgL?1. Two-way ANOVA was used to test sig- nificant differences in the phytoplankton community and environmental variables (, salinity, SS, COD, DIN, DIP, DSi, N/P, and Si/N) among different years and subregions (inner, middle section, and mouth of the bay). Prior to ANOVA, phytoplankton variables and physicochemical parameters were tested for normality (Kolmogorov-Smirnov test) and homogeneity (Levene test), and all were log transformed where necessary. Statistically significant differences (<0.05) among treatments were further as- sessed using the Student-Newman-Keuls test. The Kruskal-Wallis test was performed on data that did not satisfy the assumptions of normality and homogeneity.

    Non-metric multidimensional scaling (NMDS) was per- formed to determine spatial patterns in community struc- ture using PRIMER 6.0 software. Species abundance was log(+1) transformed and standardised before estimation using Bray-Curtis similarities between sample pairs. A two-way analysis of similarity (ANOSIM) was used to test significant differences among the phytoplankton com- munities in the different subareas and years. Phytoplank- ton species (excluded the rare species with the abundance contribution<0.01%andoccurrencefrequency<5%)abun- dance and environmental data were applied to canonical correspondence analysis (CCA) using CANOCO 4.5 soft- ware.

    3 Results

    3.1 Environmental Variables

    Both the salinity and SS decreased gradually from 2005 to 2010, while the COD, nutrients (DIN, DIP, and DSi), andincreased persistently in all the subregions of the HZB (Figs.2 and 3). Table 1 shows all the environmental parameters representing significant (<0.05) temporal differences. Except the COD, N/P, and Si/N, all the other variables showed significant (<0.01) spatial differences among the inner, middle section, and mouth of the HZB. The levels of SS, COD, and nutrients decreased with in-creasing salinity (Fig.2).

    3.2 Taxonomic Composition

    A total of four phyla (, Bacillariophyta, Dinophyta, Cyanophyta, and Chlorophyta), 43 genera, and 84 species of net-collected phytoplankton were indentified. Diatoms displayed the highest species number (67), followed by dinoflagellates (12). Diatoms dominating the phytoplank- ton community ranged in 75.0%–93.9% of the species number and 86.6%–99.0% of the total abundance during 2004–2010. Phytoplankton in HZB showed a complicated speciescomposition,includingfreshwater(,spp.), brackish (,), coastal (,), and offshore species (,) along the estuarine salinity gradient.

    Table 1 Results of the two-way ANOVA (F) or Kruskal- Wallis (H) test for phytoplankton community and environmental parameters

    Notes: SS: suspended solids; COD: chemical oxygen demand; DIN: dissolved inorganic nitrogen; DIP: dissolved inorganic phosphorus; DSi: dissolved silicate; Chl-: chlorophyll;: Shannon diversity index;: Pielou evenness index.nsNo significance,*<0.05,**<0.01,***<0.001.

    Fig.2 Summer environmental parameters in different subregions in the Hangzhou Bay during 2004–2010.

    3.3 Biomass and Species Diversity

    Both the cell density and chlorophyll(Chl-) increased gradually from 2004 to 2010, while the species diversity and evenness decreased constantly (Fig.4). Table 1 shows all these variables representing significant (<0.05) temporal differences, but there was no significant difference among different subregions in the HZB.

    Fig.4 Phytoplankton community parameters in different subregions of the Hangzhou Bay in summer during 2004–2010.

    3.4 Dominant Species

    A total of 13 dominant species were found in the HZB during 2004–2010 (Table 2).The dominant species in dif- ferent years always consisted of diatoms (,genera,,,,, and) except the year of 2009 (domi- nated by thespp.). The dominance ofincreased gradually with decreasing(Fig.5).

    Table 2 Dominance (%) of the dominant species in the Hangzhou Bay

    Note: ‘–’, not detected.

    3.5 Community Composition Analysis

    NMDS and ANOSIM showed a significant difference (Global=0.472;=0.001) in phytoplankton community in different years (Fig. 6). In addition, ANOSIM showed visible (Global=0.298;=0.001) assemblage differ- ences among different subareas in all years although there was no significant (=0.151;=0.065) difference be- tween the inner and middle sections.

    3.6 Canonical Correspondence Analysis (CCA)

    The CCA showed higher significant scores (=0.002) for Axis 1 and all canonical axes, thus the ordination re- sults are believable. These 8 environmental variables (sa- linity, SS, COD, DIN, DIP, DSi, N/P, and Si/N) in CCA explained the total variation of 25.8% in the phytoplankton community. The species-environment coefficients of Axis 1 and 2 were 0.796 and 0.829, indicating a significant relationship between the 8 environmental factors and the 75 samples. Fig.7 demonstrates that salinity was the primary factor associated with the phytoplankton assem- blage in 2004–2010.

    Fig.7 CCA ordination of phytoplankton samples with environmental variables. The symbols indicates sam- pling years which are the same as Fig.6a.

    4 Discussions

    4.1 Relationship Between Phytoplankton Distribution and Environmental Factors

    The HZB in summer is mainly influenced by the run- offs from Qiantang River and Changjiang, CDW, and Taiwan Warm Current (ECBCC, 1992). Therefore, the phytoplankton in HZB showed a complicated species composition, including the freshwater, brackish, coastal, and offshore species, based on the salinity gradient. How- ever, the use of a 76-μm net to collect phytoplankton samples in this study resulted in loss of the small-sized phytoplankton species (Jiang., 2013, 2014), and thus relatively low total species number.

    The present study found that the phytoplankton com- munity was absolutely dominated by diatoms (the domi- nant species mainly consisted of,, and), both in species number and cell density. These results are consistent with the compo- sition of net-collected microalgal community in the adja- cent Xiangshan Bay (Jiang., 2013). This finding supports the empirically based principle that diatoms re- spond rapidly to episodic high-light and high-nutrient conditions (Cloern and Dufford, 2005). However, the dominance of dinoflagellate and other taxonomic groups increased recently with decreasing diatom density. This might be attributed to the recent decreasing trend in Si/N.

    The SS concentrations in the inner and middle sections were much higher than those in the mouth of the HZB (Fig.2), due to the higher flow velocity, stronger tidal resuspension, and abundant sediment input from the Qiantang River (ECBCC, 1992; Xie., 2013). Other- wise, the freshwater discharges provided high levels of nutrients (DIN, DIP, and DSi) and COD in the HZB along the salinity gradient (Fig.2). However, the light limitation (due to the extremely high SS concentration) and change- able osmotic stresses in the upstream estuarine bay re- strained the phytoplankton growth and primary produc- tion (Fig.4), although the inner and middle sections pre- sented highlevel (Fig.3). Many studies found that high turbidity is an important physical process that limits light availability for estuarine primary productivity, such as the cases of the CE (Zhu., 2009) and Guadiana Estuary (Domingues., 2011). In contrast, high cell density and Chl-were found in the mouth in our study with rela- tively high nutrients and deep light penetration (Fig.2), due to exchanges with the external clear water impacted by the Taiwan Warm Current.

    CCA showed that salinity was the primary variable as- sociated with the phytoplankton community (Fig.7). This supports the empirically-based principle that salinity varying geographically is the most important factor which determines the phytoplankton composition and distribution in the estuarine waters (, Van Damme., 2009). The hydrological (circulation, SS, and salinity) and nu- trient conditions in the HZB were significantly different from the inner section to mouth, leading to the obvious spatial differences in assemblage compositions (Fig.6). According to the NMDS and ANOSIM, the phytoplank- ton also showed spatial heterogeneity along the salinity gradient, driven by the different environmental factors among different subareas (Cloern and Dufford, 2005; Jiang., 2013).

    4.2 Long-Term Changes in Phytoplankton Under Increasing Human and Climatic Perturbations

    The estuarine HZB, as a transition zone between river and ocean environments, is subject to riverine (, dis- charges of sediment and freshwater associated with huge nutrients) and marine (, tides, waves, and saline water intrusion) mutual impact. The changes in phytoplankton community in the HZB are profoundly impacted by the above-mentioned environmental factors, especially the runoff, sediment, nutrients, and water exchange (Van Damme., 2009). Yu. (2013) simulated that the phytoplankton dynamics in the estuarine Chesapeake Bay were controlled by the river flow associated with terri- genous materials, according to a 3-D physical-biogeo- chemical model based on ROMS.

    Large-scale dam construction in the Qiantang River and Changjiang basins (Xu., 2010; Yang., 2011) retained huge sediments, while the recent precipitation increase with decreasing evaporation (Wang., 2011, 2013) enhanced the freshwater discharges from Qiantang River (Fig.8). These anthropogenic and climatic changes resulted in an increase in nutrients and reduction in SS and salinity (Figs.2 and 9), which profoundly influenced the phytoplankton community variations. For example, the turbidity reduction (transparency enhancement) com- bined with eutrophication exacerbation (Figs.2 and 3) in this seriously light-limited HZB promoted the algal growth. Fig.10 shows that phytoplankton cell density and Chl-increased gradually during 1982–2012. In addition, the temperature increase caused by global warming might have amplified the microalgal biomass. Jiang. (2014) found that summer SST in the CE increased significantly at a rate of 0.49℃/decade during 1982 to 2009.

    Due to the extensive use of chemical fertilizers and the large discharge of industrial wastewater and domestic sewage, the riverine N and P loadings increased persis- tently, which induced continual elevation in DIN and DIP in the HZB, as well as drastic alteration in nutrient ratios (significant increase in N/P and decrease in Si/N) since 1980s (Fig.9). The results are similar to those in the pre- vious studies on the adjacent CE (Zhou., 2008; Jiang., 2010; Li., 2010; Gao., 2012; Jiang., 2014). Meanwhile, DSi concentration in the HZB did not change a lot, although the sediments from Changjiang and Qiantang River decreased sharply (Fig.8). Gao. (2012) concluded that the sharply reduced sediments do not nec- essarily mean the weathering rates of silicate rocks have also decreased. Moreover, large-scale tidal flat reclama-tion (Ni and Lin, 2003; Zhang., 2010) changed the hydrodynamic environment (tide, flow field, residual currents,), reduced the water exchange and residence time, prolonged the nutrient-lingering period, and thereby aggravated the eutrophication (Jiang., 2013). This adverse factor would exacerbate harmful algal blooms in the HZB mouth as well as the adjacent CE and Zhoushan archipelago in terms of frequency, range, persistence, and destructive capability (Zhou., 2008; Wang and Wu, 2009; Li., 2010).

    Fig.8 Long-term changes in annual average runoff and sediment of (a) Qiantang River and (b) Changjiang (data from the website: http://www.irtces.org/nishagb_2012.asp).

    Fig.9 Long-term changes of nutrient parameters in the Hangzhou Bay in summer (other data from personal data; ECBCC, 1992; Liu et al., 2001; Zhang et al., 2002; personal data).

    Fig.10 Long-term changes in summer phytoplankton community parameters in the Hangzhou Bay (other data from Liu and Ning, 1994; Liu et al., 2001; Cai, 2006; Wang et al., 2008; SOA, 2012, 2013; personal data).

    Fig.5 and Table 3 show that the dominance of non-group livingdecreased gradually, while that of small-size chain-forming species (,,, and) increased, especially the. This change might refer to eutrophication exacerbation and runoff increase (salinity decrease). As an extremely high siliceous group, the sharp decline of Si/N inevitably inhibited the growth and re- production of, while the relatively low siliceous diatoms (,, and) profited from this condition. These species, as-strategists, can rapidly multiply under increasing eutro- phication, particularly the(Liu., 2002). Besides, as estuarine species,spp. gradually takes the advantage of community competition in the tur- bid estuary with low salinity, such as the CE. Many field investigations (, Guo and Yang, 1992; Luan., 2006) found that this microalgae population generally dominates (>95% in total cell density) and blooms in the principal part of CDW with the salinity <23. The labora- tory experimental study by Li. (2005) supports this finding. They found that the optimal growth salinity ofis 19.2 with the suitable range of 14–23. Thus, on the conditions of decreasing salinity (<20), the domi- nance ofspp. in the HZB has increased ob- viously.

    Table 3 Long-term changes in summer dominant species in the Hangzhou Bay

    Under increasing Si-limitation (Si/N decrease) in the HZB, both the richness and cell density of the diatom species decreased gradually, while the dominance of dinoflagellate and other taxonomic groups (mainly the cyanobacteria) increased (Fig.11). Similarly, diatom-dinoflagellate ratio in species and cell number in the ad- jacent CE (Zhou., 2008; Li., 2010; Jiang., 2014) and Xiangshan Bay (Jiang., 2013) decreased significantly under the decreasing Si/N stresses. Besides, CDW shrinkage (caused by recent sharp reduction in runoffs from Changjiang in Fig.8; Jiang., 2014) combined with increasing intrusion of Taiwan Warm Current (Tang., 2009) brought about abundant off- shorein summer (Table 3), in despite of the salinity reduction in the HZB. However, due to the abso- lutely dominance by several chain-forming species, the phytoplankton community diversity and evenness de- creased evidently (Fig.10). This abnormal phenomenon was also found in the CE (Li., 2010).

    Fig.11 Long-term changes in summer phytoplankton community composition in the Hangzhou Bay.

    5 Conclusions

    The summer net-collected phytoplankton community in HZB was dominated by diatoms during 2004–2010, mainly consisting of genera,, and. Planktonic assemblages showed sig- nificant spatial heterogeneity along the salinity gradient. Increasing human activities and climatic forcing resulted in a series of physicochemical alterations, including eu- trophication exacerbation, light limitation alleviation, and salinity reduction. These changes induced an obvious increase in algal biomass with decreasing species diver- sity and diatom-dinoflagellate ratio as well as the shifting of dominant species. The present study confirmed that long-term changes in phytoplankton community in the HZB were closely related to the anthropogenic and cli- matic perturbations. However, the present use of a net with 76-μm mesh size resulted in loss of the small-sized phytoplankton species. We should obtain more informa- tion from water samples of microalgae with small cell sizes, and examine the data from net-collected samples to better understand the responses of phytoplankton to envi- ronmental changes in future.

    Acknowledgements

    The data that this paper used was provided by the Ma- rine Monitoring and Forecasting Center of Zhejiang. This work was funded by the National Basic Research Pro- gram of China (2010CB428903), the National Marine Public Welfare Research Project of China (201305009 and 201305043-3), the National Natural Science Founda- tion of China (41206103), and Basic Scientific Research of Second Institute of Oceanography, SOA (JG1222).

    Cai, Y. H., 2006. Study on the phytoplankton diversity in the Hangzhou Bay. Master thesis. Ocean University of China.

    Cloern, J. E., and Dufford, R., 2005. Phytoplankton community ecology: Principles applied in San Francisco Bay., 285: 11-28.

    Domingues, R. B., Anselmo, T. P., Barbosa, A. B., Sommer, U., and Galv?o, H. M., 2011. Light as a driver of phytoplankton growth and production in the freshwater tidal zone of a turbid estuary., 91 (4): 526-535.

    Editorial Committee of the Bay Chorography in China (ECBCC), 1992.. China Ocean Press, Beijing, 166-233.

    Gao, L., Li, D. J., and Zhang, Y. W., 2012. Nutrients and par- ticulate organic matter discharged by the Changjiang (Yang- tze River): Seasonal variations and temporal trends., 117 (G4), G04001.

    Gao, S. Q., Chen, J. F., Jin, H. Y., Wang, K., Lu, Y., Li, H. L., and Chen, F. J., 2011. Characteristics of nutrients and eutro- phication in the Hangzhou Bay and its adjacent waters., 29 (3): 36-47.

    Guo, Y. J., and Yang, Z. Y., 1992. Quantitative variation and ecological analysis of phytoplankton in the estuarine area of the Changjiang River., 33: 167-189.

    Han, Z. C., Dai, Z. H., and Li, G. B., 2003.. China Water Power Press, Beijing, 122-131.

    Jiang, T., Yu, Z. M., Song, X. X., Cao, X. H., and Yuan, Y. Q., 2010. Long-term ecological interactions between nutrient and phytoplankton community in the Changjiang Estuary., 28 (4): 887-898.

    Jiang, Z. B., Liu, J. J., Chen, J. F., Chen, Q. Z., Yan, X. J., Xuan, J. L., and Zeng, J. N., 2014. Responses of summer phyto- plankton community to drastic environmental changes in the Changjiang (Yangtze River) Estuary during the past 50 years., 54: 1-11.

    Jiang, Z. B., Zhu, X. Y., Gao, Y., Chen, Q. Z., Zeng, J. N., and Zhu, G. H., 2013. Spatio-temporal distribution of net-col- lected phytoplankton community and its response to marine exploitation in Xiangshan Bay., 31 (4): 762-773.

    Li, J. T., Zhao, W. H., Yang, D. F., and Wang, J. T., 2005. Effect of turbid water in Changjiang (Yangtze) Estuary on the growth of., 29 (1): 34-37.

    Li, Y., Li, D. J., Tang, J. L., Wang, Y. M., Liu, Z. G., and He, S. Q., 2010. Long-term changes in the Changjiang Estuary plankton community related to anthropogenic eutrophication., 13 (1): 66-72.

    Liu, D. Y., Sun, J., Chen, Z. T., and Wei, T. D., 2002. Effect of N/P ratio on the growth of a red tide diatom., 2002 (2): 39-44.

    Liu, Z. L., and Ning, X. R., 1994. Phytoplankton standing stock and primary production in the front of Hangzhou Bay., 12 (4): 58-66.

    Liu, Z. L., Ning, X. R., and Cai, Y. M., 2001. Primary produc- tivity and standing stock of the Hangzhou Bay to the Zhou- shan fishing ground during autumn., 23 (2): 93-99.

    Luan, Q. S., Sun, J., Shen, Z. L., Song, S. Q., and Wang, M., 2006. Phytoplankton assemblage of Yangtze River Estuary and the adjacent East China Sea in summer., 5 (2): 123-131.

    Ni, Y. Q., and Lin, J., 2003. The effects of regulation and recla- mation in Qiantang Estuary on Hangzhou Bay., 21 (3): 73-77.

    State Oceanic Administration (SOA), 2012..

    State Oceanic Administration (SOA), 2013..

    Su, J. L., and Wang, K. S., 1989. Changjiang River plume and suspended sediment transport in Hangzhou Bay., 9 (1): 93-111.

    Tang, X. H., Wang, F., Chen, Y. L., and Li, M. L., 2009. Warm- ing trend in northern East China Sea in recent four decades., 27 (2): 185-191.

    Van Damme, S., Struyf, E., Maris, T., Cox, T., and Meire, P., 2009. Characteristic aspects of the tidal freshwater zone that affect aquatic primary production In:. Barendregt, A.,, eds., Backhuys, 123-136.

    Wang, J. H., and Wu, J. Y., 2009. Occurrence and potential risks of harmful algal blooms in the East China Sea., 407 (13): 4012-4021.

    Wang, Y. J., Ding, Y. J., Ye, B. S., Liu, F. J., and Wang, J., 2013. Contributions of climate and human activities to changes in runoff of the Yellow and Yangtze rivers from 1950 to 2008., 56 (8): 1398-1412.

    Wang, Y. J., Liu, B., Su, B. D., Zhai, J. Q., and Gemmer, M., 2011. Trends of calculated and simulated actual evaporation in the Yangtze River Basin., 24: 4494-4507.

    Wang, Y. L., Yuan, Q., and Shen, X. Q., 2008. Distribution status and change tendency of phytoplankton during summer in Changjiang Estuary and adjacent waters., 27 (2):169-172.

    Xie, D. F., Pan, C. H., Cao, Y., and Zhang, B. H., 2013. Decadal variations in the erosion/deposition pattern of the Hangzhou Bay and their mechanism in recent 50a., 35 (4): 121-178.

    Xie, D. F., Wang, Z. B., Gao, S., and De Vriend, H. J., 2009. Modeling the tidal channel morphodynamics in a macro-tidal embayment, Hangzhou Bay, China., 29 (15): 1757-1767.

    Xu, K., Milliman, J. D., and Xu, H., 2010. Temporal trend of precipitation and runoff in major Chinese rivers since 1951., 73: 219-232.

    Yang, S. L., Milliman, J. D., Li, P., and Xu, K., 2011. 50000 dams later: Erosion of the Yangtze River and its delta., 75: 14-20.

    Yu, Q. Y., Wang, Y., Tang, X. X., and Li, M., 2013. River flow control on the phytoplankton dynamics of Chesapeake Bay., 12 (1): 103-114.

    Zhang, H. G., Sui, Y. Z., and Huang, W. G., 2010. Patterns of reclamation land use of Hangzhou Bay with remote sensing in the last two decades.. Michel, U., and Civco, D. L., eds., Proceedings of SPIE Vol. 7831, 78311E, 6pp.

    Zhang, J., Shi, Q. S., Wu, X. Y., Ye, X. R., and Zhu, G. H., 2002. Distribution characteristic analysis of main pollution factor in rainy season in the Hangzhou Bay., 20 (4): 35-41.

    Zhang, S. Y., Shao, J. B., and Dai, X. J., 2001. Studies on eu- trophication and phytoplankton diversity in Hangzhou Bay., 25 (6): 512-517.

    Zhou, M. J., Shen, Z. L., and Yu, R. C., 2008. Responses of a coastal phytoplankton community to increased nutrient input from the Changjiang (Yangtze) River., 28 (12): 1483-1489.

    Zhou, Y., Zhao, C. J., Gao, Y. S., Long, H., and Yu, J., 2010. Variation and distribution characteristics of phytoplankton in ecology-monitoring area of Hangzhouwan Bay from 2005 to 2008., 28 (2): 28-35.

    Zhu, Z. Y., Ng, W. M., Liu, S. M., Zhang, J., Chen, J. C., and Wu, Y., 2009. Estuarine phytoplankton dynamics and shift of limiting factors: A study in the Changjiang (Yangtze River) Estuary and adjacent area., 84 (3): 393-401.

    (Edited by Xie Jun)

    DOI 10.1007/s11802-015-2483-6

    ISSN 1672-5182, 2015 14 (6): 1025-1033

    ? Ocean University of China, Science Press and Springer-Verlag Berlin Heidelberg 2015

    (September 12, 2013; revised March 2, 2014; accepted October 1, 2015)

    * Corresponding author. E-mail: wangh@nmefc.gov.cn

    18禁美女被吸乳视频| 美女 人体艺术 gogo| 中出人妻视频一区二区| 欧美日韩国产mv在线观看视频| 欧美日韩黄片免| 亚洲黑人精品在线| 成人手机av| svipshipincom国产片| 两性午夜刺激爽爽歪歪视频在线观看 | 老司机午夜十八禁免费视频| 亚洲一卡2卡3卡4卡5卡精品中文| av网站在线播放免费| 我的亚洲天堂| 少妇 在线观看| 在线视频色国产色| 香蕉丝袜av| 久久香蕉激情| 久久 成人 亚洲| 动漫黄色视频在线观看| 免费在线观看完整版高清| 亚洲av熟女| 性少妇av在线| 亚洲一区二区三区不卡视频| 午夜免费成人在线视频| videosex国产| 精品无人区乱码1区二区| 亚洲欧美日韩高清在线视频| 999精品在线视频| 一进一出好大好爽视频| 这个男人来自地球电影免费观看| 久久久久精品国产欧美久久久| 亚洲美女黄片视频| 亚洲精品成人av观看孕妇| 老司机影院毛片| 国产三级黄色录像| 欧美激情高清一区二区三区| tocl精华| 欧美黄色淫秽网站| 久久久久久久久免费视频了| 搡老乐熟女国产| 国产成人免费无遮挡视频| 99国产精品一区二区三区| 欧美乱色亚洲激情| 久久久久久人人人人人| 国产不卡一卡二| 宅男免费午夜| 1024香蕉在线观看| 亚洲成av片中文字幕在线观看| 成人免费观看视频高清| 99国产精品一区二区蜜桃av | 看黄色毛片网站| 欧美日韩瑟瑟在线播放| 999精品在线视频| 在线观看一区二区三区激情| 91麻豆精品激情在线观看国产 | 成在线人永久免费视频| 1024视频免费在线观看| 日韩一卡2卡3卡4卡2021年| 午夜久久久在线观看| 国产成人啪精品午夜网站| 成人国产一区最新在线观看| 久热爱精品视频在线9| 成人特级黄色片久久久久久久| 国产亚洲一区二区精品| 亚洲少妇的诱惑av| 精品视频人人做人人爽| 亚洲av成人av| 99国产精品一区二区蜜桃av | 超碰成人久久| 女人被躁到高潮嗷嗷叫费观| 一区福利在线观看| 色94色欧美一区二区| 一个人免费在线观看的高清视频| 久99久视频精品免费| 丁香六月欧美| 老熟女久久久| 满18在线观看网站| 制服人妻中文乱码| 热re99久久精品国产66热6| 人人妻人人爽人人添夜夜欢视频| 18禁黄网站禁片午夜丰满| 99精品在免费线老司机午夜| 村上凉子中文字幕在线| 夜夜爽天天搞| 50天的宝宝边吃奶边哭怎么回事| 99国产极品粉嫩在线观看| 老司机影院毛片| 首页视频小说图片口味搜索| 一二三四在线观看免费中文在| 精品视频人人做人人爽| 欧美亚洲 丝袜 人妻 在线| 91大片在线观看| 亚洲精品在线美女| 亚洲av电影在线进入| 99热网站在线观看| 精品福利观看| 国产欧美日韩精品亚洲av| 99精品久久久久人妻精品| 亚洲 欧美一区二区三区| www日本在线高清视频| 在线天堂中文资源库| 一a级毛片在线观看| 少妇被粗大的猛进出69影院| a级片在线免费高清观看视频| 亚洲人成电影观看| 美女福利国产在线| 国产又爽黄色视频| 18在线观看网站| 两个人看的免费小视频| 国产男靠女视频免费网站| 精品国产乱码久久久久久男人| 国产单亲对白刺激| 国产视频一区二区在线看| 亚洲第一青青草原| 久久这里只有精品19| 咕卡用的链子| 两个人免费观看高清视频| 在线观看免费日韩欧美大片| 久久精品国产综合久久久| 国产亚洲精品第一综合不卡| a级毛片在线看网站| 欧美大码av| 婷婷成人精品国产| 久久中文字幕一级| 中出人妻视频一区二区| 国产无遮挡羞羞视频在线观看| 久久久久精品人妻al黑| 两人在一起打扑克的视频| 欧美性长视频在线观看| 丰满的人妻完整版| 免费在线观看黄色视频的| 亚洲国产欧美网| 国产一区有黄有色的免费视频| 亚洲av欧美aⅴ国产| 久久ye,这里只有精品| 久久久久久免费高清国产稀缺| 国产欧美日韩精品亚洲av| 精品久久久久久电影网| 一二三四在线观看免费中文在| 欧美黑人欧美精品刺激| 亚洲欧美激情综合另类| 麻豆成人av在线观看| 亚洲精品美女久久av网站| 夜夜夜夜夜久久久久| 曰老女人黄片| 日韩三级视频一区二区三区| 色综合欧美亚洲国产小说| 国产精品1区2区在线观看. | 亚洲欧美激情综合另类| 国产亚洲av高清不卡| 十八禁网站免费在线| 精品无人区乱码1区二区| 欧美国产精品va在线观看不卡| 久久天堂一区二区三区四区| 最新在线观看一区二区三区| 999久久久精品免费观看国产| 黄色视频不卡| 欧美日韩亚洲国产一区二区在线观看 | 啦啦啦免费观看视频1| 久久久久国内视频| 欧美乱色亚洲激情| 欧美日本中文国产一区发布| 大片电影免费在线观看免费| 久久久久久人人人人人| 亚洲五月天丁香| 亚洲在线自拍视频| 一a级毛片在线观看| 成人手机av| 亚洲欧美一区二区三区久久| 深夜精品福利| 精品久久蜜臀av无| 12—13女人毛片做爰片一| 国产精品av久久久久免费| 国产精品免费大片| 伊人久久大香线蕉亚洲五| 亚洲av第一区精品v没综合| 国产aⅴ精品一区二区三区波| 91字幕亚洲| 色综合欧美亚洲国产小说| 欧美精品亚洲一区二区| 色94色欧美一区二区| www.精华液| 12—13女人毛片做爰片一| 狠狠婷婷综合久久久久久88av| 久久久久精品国产欧美久久久| 日韩免费av在线播放| 日本a在线网址| 久久影院123| 日本五十路高清| 国产精品久久久久久人妻精品电影| 黄片播放在线免费| 日韩 欧美 亚洲 中文字幕| 国产精华一区二区三区| 不卡av一区二区三区| 久久中文看片网| 性少妇av在线| 精品久久久精品久久久| 黑人操中国人逼视频| 天堂中文最新版在线下载| 亚洲精品自拍成人| 免费在线观看完整版高清| 成人黄色视频免费在线看| 天天躁日日躁夜夜躁夜夜| 精品人妻在线不人妻| 午夜精品久久久久久毛片777| av不卡在线播放| 亚洲成人国产一区在线观看| 777米奇影视久久| 精品一区二区三区视频在线观看免费 | av中文乱码字幕在线| 老鸭窝网址在线观看| 免费一级毛片在线播放高清视频 | 亚洲美女黄片视频| 在线观看免费午夜福利视频| 老司机在亚洲福利影院| 亚洲性夜色夜夜综合| 多毛熟女@视频| 成人国产一区最新在线观看| 亚洲五月天丁香| 亚洲av片天天在线观看| 久久久久精品人妻al黑| av天堂久久9| 91成人精品电影| 十八禁人妻一区二区| 午夜福利影视在线免费观看| 亚洲熟妇熟女久久| 制服人妻中文乱码| 欧美黑人精品巨大| 18禁裸乳无遮挡动漫免费视频| 午夜两性在线视频| 免费人成视频x8x8入口观看| 欧美日韩国产mv在线观看视频| 在线免费观看的www视频| 亚洲成a人片在线一区二区| 国产亚洲精品久久久久久毛片 | 亚洲五月色婷婷综合| 久久中文字幕人妻熟女| 国产99久久九九免费精品| 99re在线观看精品视频| 女警被强在线播放| 中文字幕高清在线视频| 国产精品免费大片| 超色免费av| av网站在线播放免费| 免费在线观看日本一区| 精品亚洲成a人片在线观看| 最新的欧美精品一区二区| 久久精品亚洲精品国产色婷小说| 在线天堂中文资源库| 色尼玛亚洲综合影院| 免费看a级黄色片| 国产精华一区二区三区| 亚洲欧美精品综合一区二区三区| 欧美日韩乱码在线| 午夜91福利影院| 国产91精品成人一区二区三区| 两人在一起打扑克的视频| 欧美日韩黄片免| 亚洲av欧美aⅴ国产| 亚洲欧美一区二区三区久久| 五月开心婷婷网| 国产片内射在线| 午夜影院日韩av| 亚洲熟女精品中文字幕| 精品福利观看| 欧美在线一区亚洲| 欧美国产精品va在线观看不卡| 久久精品人人爽人人爽视色| 精品少妇久久久久久888优播| 国产av一区二区精品久久| 99热网站在线观看| 亚洲精品在线观看二区| 国产精品偷伦视频观看了| 免费日韩欧美在线观看| 精品国产一区二区三区久久久樱花| 久久久久久久久免费视频了| 欧美激情 高清一区二区三区| 成人三级做爰电影| 新久久久久国产一级毛片| 久久香蕉激情| 飞空精品影院首页| 两性夫妻黄色片| 亚洲精品一卡2卡三卡4卡5卡| 国产精品av久久久久免费| 欧美最黄视频在线播放免费 | 国产精品香港三级国产av潘金莲| 日本wwww免费看| 成人特级黄色片久久久久久久| 欧美在线一区亚洲| 国产精品国产高清国产av | 欧美日韩精品网址| 亚洲欧美一区二区三区黑人| 高清欧美精品videossex| 啦啦啦 在线观看视频| 中文字幕最新亚洲高清| 80岁老熟妇乱子伦牲交| 欧美一级毛片孕妇| 日本欧美视频一区| 女人被狂操c到高潮| 日韩制服丝袜自拍偷拍| 国产成人av教育| 动漫黄色视频在线观看| 久久性视频一级片| 亚洲熟女毛片儿| 大型av网站在线播放| netflix在线观看网站| 19禁男女啪啪无遮挡网站| 中文字幕最新亚洲高清| 91麻豆精品激情在线观看国产 | 麻豆成人av在线观看| 大型av网站在线播放| 成人永久免费在线观看视频| 香蕉丝袜av| 丝袜美足系列| 视频区欧美日本亚洲| 看免费av毛片| 久久久水蜜桃国产精品网| xxx96com| 天堂动漫精品| 久久青草综合色| 亚洲男人天堂网一区| 黄色丝袜av网址大全| 久久精品熟女亚洲av麻豆精品| 国产成人免费观看mmmm| 日韩成人在线观看一区二区三区| e午夜精品久久久久久久| 久久久久久人人人人人| 一进一出抽搐动态| 久久久久久久国产电影| 18禁黄网站禁片午夜丰满| 悠悠久久av| 午夜精品在线福利| 国产精品久久久久久精品古装| 高清毛片免费观看视频网站 | 午夜精品国产一区二区电影| 国产欧美日韩综合在线一区二区| 成人国产一区最新在线观看| 老司机午夜十八禁免费视频| 免费看a级黄色片| x7x7x7水蜜桃| 亚洲精品在线观看二区| 最新在线观看一区二区三区| 国产一区有黄有色的免费视频| 亚洲免费av在线视频| 日韩欧美国产一区二区入口| e午夜精品久久久久久久| 国产精品影院久久| 免费在线观看亚洲国产| 大型黄色视频在线免费观看| 极品人妻少妇av视频| 丝瓜视频免费看黄片| 女警被强在线播放| 午夜免费成人在线视频| 成人国语在线视频| 国产xxxxx性猛交| 91麻豆av在线| 人妻 亚洲 视频| 婷婷丁香在线五月| 国产亚洲一区二区精品| 国产欧美日韩一区二区三区在线| 国产成人欧美| 亚洲熟妇中文字幕五十中出 | 精品国产一区二区久久| 亚洲av日韩精品久久久久久密| 国产极品粉嫩免费观看在线| 国产日韩一区二区三区精品不卡| 人人妻,人人澡人人爽秒播| 亚洲五月婷婷丁香| 午夜福利影视在线免费观看| 美女高潮到喷水免费观看| 黄色丝袜av网址大全| videos熟女内射| 在线观看日韩欧美| 亚洲欧美日韩另类电影网站| 狠狠狠狠99中文字幕| 美女福利国产在线| 国产一卡二卡三卡精品| 亚洲欧美激情综合另类| 每晚都被弄得嗷嗷叫到高潮| 精品免费久久久久久久清纯 | 国产又爽黄色视频| 国产熟女午夜一区二区三区| 精品国产亚洲在线| 久久久久久久精品吃奶| 亚洲精品国产色婷婷电影| 女人高潮潮喷娇喘18禁视频| 国产熟女午夜一区二区三区| 日韩中文字幕欧美一区二区| 国产精品av久久久久免费| 国精品久久久久久国模美| 国产91精品成人一区二区三区| 亚洲国产中文字幕在线视频| 亚洲伊人色综图| 99香蕉大伊视频| 精品福利观看| 久久人妻福利社区极品人妻图片| 久久 成人 亚洲| 亚洲成av片中文字幕在线观看| 日韩欧美在线二视频 | 免费在线观看亚洲国产| 久久婷婷成人综合色麻豆| 老鸭窝网址在线观看| 久久精品亚洲精品国产色婷小说| 久久久水蜜桃国产精品网| 亚洲精品国产色婷婷电影| 新久久久久国产一级毛片| 久久久久久久久久久久大奶| 99久久99久久久精品蜜桃| 18禁国产床啪视频网站| 久久久久国产一级毛片高清牌| 日韩欧美国产一区二区入口| 18在线观看网站| 久久久久国产精品人妻aⅴ院 | 女人爽到高潮嗷嗷叫在线视频| 老汉色av国产亚洲站长工具| 欧美成人免费av一区二区三区 | 精品久久久久久久毛片微露脸| 午夜成年电影在线免费观看| videos熟女内射| 男女免费视频国产| 制服人妻中文乱码| 久久国产乱子伦精品免费另类| 男女下面插进去视频免费观看| 欧美另类亚洲清纯唯美| 乱人伦中国视频| 一区在线观看完整版| 久久国产精品影院| 性色av乱码一区二区三区2| 色播在线永久视频| 1024香蕉在线观看| 国产精品久久电影中文字幕 | 中亚洲国语对白在线视频| 欧美国产精品va在线观看不卡| 丝袜美腿诱惑在线| 久久 成人 亚洲| 久久天堂一区二区三区四区| 狂野欧美激情性xxxx| 国产高清视频在线播放一区| 电影成人av| 9191精品国产免费久久| 久久久久精品人妻al黑| 777米奇影视久久| 日本五十路高清| 成人黄色视频免费在线看| 亚洲五月色婷婷综合| cao死你这个sao货| 国产精品久久久久成人av| 国产一区二区三区在线臀色熟女 | 色播在线永久视频| 捣出白浆h1v1| 搡老乐熟女国产| 国产高清激情床上av| 欧美丝袜亚洲另类 | 亚洲视频免费观看视频| 日本撒尿小便嘘嘘汇集6| 天天躁狠狠躁夜夜躁狠狠躁| 久久 成人 亚洲| 久久香蕉国产精品| 中文字幕高清在线视频| 久久99一区二区三区| 欧美日韩瑟瑟在线播放| 成人手机av| 自线自在国产av| 国产单亲对白刺激| 1024香蕉在线观看| 亚洲成人手机| 亚洲av日韩在线播放| 久久香蕉精品热| 国产极品粉嫩免费观看在线| 国产真人三级小视频在线观看| 成人18禁在线播放| 国产一区二区激情短视频| 国产免费男女视频| 嫁个100分男人电影在线观看| 久久久久久久国产电影| 婷婷丁香在线五月| 高清视频免费观看一区二区| 韩国av一区二区三区四区| 啪啪无遮挡十八禁网站| 免费在线观看影片大全网站| 国产单亲对白刺激| 亚洲,欧美精品.| 黄色毛片三级朝国网站| 18禁美女被吸乳视频| 亚洲中文字幕日韩| 精品一区二区三区四区五区乱码| 一区二区日韩欧美中文字幕| 亚洲精品av麻豆狂野| 国内毛片毛片毛片毛片毛片| 欧美成人午夜精品| 国产麻豆69| av欧美777| 精品国产乱子伦一区二区三区| 大陆偷拍与自拍| 国产精品1区2区在线观看. | 精品国产超薄肉色丝袜足j| x7x7x7水蜜桃| 日韩欧美一区二区三区在线观看 | 人妻一区二区av| 久久精品国产综合久久久| 日韩人妻精品一区2区三区| 老汉色∧v一级毛片| 久久久国产成人精品二区 | 国产精品偷伦视频观看了| 婷婷精品国产亚洲av在线 | 国产亚洲欧美在线一区二区| 免费日韩欧美在线观看| 后天国语完整版免费观看| 无遮挡黄片免费观看| 最新美女视频免费是黄的| 一个人免费在线观看的高清视频| av免费在线观看网站| 国产欧美日韩一区二区精品| 啦啦啦免费观看视频1| 黄色a级毛片大全视频| 一二三四社区在线视频社区8| 夜夜躁狠狠躁天天躁| 亚洲全国av大片| 大香蕉久久网| 黄色怎么调成土黄色| 一本一本久久a久久精品综合妖精| 久久ye,这里只有精品| 一区二区日韩欧美中文字幕| 国产精品 欧美亚洲| 国产一区在线观看成人免费| 成人国语在线视频| 国产在线观看jvid| 老司机午夜福利在线观看视频| 国产一区有黄有色的免费视频| 久久久久国产一级毛片高清牌| av在线播放免费不卡| 国产精品 欧美亚洲| 亚洲成人国产一区在线观看| 黄片播放在线免费| 亚洲精品国产精品久久久不卡| 国产1区2区3区精品| 国产男女超爽视频在线观看| 黑人操中国人逼视频| 狠狠婷婷综合久久久久久88av| 女性被躁到高潮视频| 人成视频在线观看免费观看| 久久久久国产精品人妻aⅴ院 | 热re99久久国产66热| 99精国产麻豆久久婷婷| 热99国产精品久久久久久7| av网站免费在线观看视频| 黄色成人免费大全| 在线观看免费日韩欧美大片| 亚洲国产欧美日韩在线播放| 久久香蕉激情| 18禁观看日本| 国产极品粉嫩免费观看在线| 午夜日韩欧美国产| 国产高清激情床上av| av天堂久久9| 国产精品.久久久| 高清欧美精品videossex| 久久精品国产亚洲av高清一级| 麻豆乱淫一区二区| 制服诱惑二区| 亚洲精品中文字幕在线视频| 操美女的视频在线观看| 国产一区有黄有色的免费视频| av网站在线播放免费| 亚洲精品一二三| xxxhd国产人妻xxx| 视频区欧美日本亚洲| 午夜福利影视在线免费观看| 国产精品 国内视频| 免费在线观看影片大全网站| 日本vs欧美在线观看视频| 精品一区二区三卡| 9色porny在线观看| 成年人免费黄色播放视频| 成人精品一区二区免费| 女性生殖器流出的白浆| 在线观看免费高清a一片| 99精国产麻豆久久婷婷| 色老头精品视频在线观看| 亚洲精华国产精华精| 天堂俺去俺来也www色官网| 久99久视频精品免费| 欧美+亚洲+日韩+国产| 十分钟在线观看高清视频www| 国产淫语在线视频| 亚洲情色 制服丝袜| 91成年电影在线观看| www日本在线高清视频| 桃红色精品国产亚洲av| 久久久久久人人人人人| 最新在线观看一区二区三区| 国产高清激情床上av| 亚洲欧美日韩高清在线视频| 日韩视频一区二区在线观看| 9热在线视频观看99| tocl精华| 精品亚洲成a人片在线观看| 国产男女内射视频| 亚洲久久久国产精品| 搡老熟女国产l中国老女人| 最新的欧美精品一区二区| 日本欧美视频一区| 高潮久久久久久久久久久不卡| 亚洲国产欧美一区二区综合| 亚洲精品在线观看二区| 高潮久久久久久久久久久不卡| 日韩人妻精品一区2区三区| 午夜福利视频在线观看免费| 日日爽夜夜爽网站| 电影成人av| 久久国产乱子伦精品免费另类| 日日摸夜夜添夜夜添小说| 桃红色精品国产亚洲av| 在线视频色国产色| 成人精品一区二区免费|