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

    Seasonal Variations in Phytoplankton Community Structure in the Sanggou, Ailian, and Lidao Bays

    2014-04-26 10:54:51YUANMingliZHANGCuixiaJIANGZengjieGUOShujinandSUNJun
    Journal of Ocean University of China 2014年6期
    關(guān)鍵詞:數(shù)據(jù)通訊交通管理數(shù)據(jù)結(jié)構(gòu)

    YUAN Mingli, ZHANG Cuixia, JIANG Zengjie, GUO Shujin, and SUN Jun,

    1) College of Marine Science and Engineering, Tianjin University of Science and Technology, Tianjin 300457, P. R. China

    2) Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071, P. R. China

    Seasonal Variations in Phytoplankton Community Structure in the Sanggou, Ailian, and Lidao Bays

    YUAN Mingli1), ZHANG Cuixia1), JIANG Zengjie2), GUO Shujin1), and SUN Jun1),*

    1) College of Marine Science and Engineering, Tianjin University of Science and Technology, Tianjin 300457, P. R. China

    2) Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Qingdao 266071, P. R. China

    The seasonal variations in phytoplankton community structure were investigated for the Sanggou Bay (SGB) and the adjacent Ailian Bay (ALB) and Lidao Bay (LDB) in Shandong Peninsula, eastern China. The species composition and cell abundance of phytoplankton in the bay waters in spring (April 2011), summer (August 2011), autumn (October 2011), and winter (January 2012) were examined using the Uterm?hl method. A total of 80 taxa of phytoplankton that belong to 39 genera of 3 phyla were identified. These included 64 species of 30 genera in the Phylum Bacillariophyta, 13 species of 8 genera in the Phylum Dinophyta, and 3 species of 1 genus in the Phylum Chrysophyta. During the four seasons, the number of phytoplankton species (43) was the highest in spring, followed by summer and autumn (40), and the lowest number of phytoplankton species (35) was found in winter. Diatoms, especially Paralia sulcata (Ehrenberg) Cleve and Coscinodiscus oculus-iridis Ehrenberg, were predominant in the phytoplankton community throughout the study period, whereas the dominance of dinoflagellate appeared in summer only. The maximum cell abundance of phytoplankton was detected in summer (average 8.08 × 103cells L-1) whereas their minimum abundance was found in autumn (average 2.60 × 103cells L-1). The phytoplankton abundance was generally higher in the outer bay than in the inner bay in spring and autumn. In summer, the phytoplankton cells were mainly concentrated in the south of inner SGB, with peak abundance observed along the western coast. In winter, the distribution of phytoplankton cells showed 3 patches, with peak abundance along the western coast as well. On seasonal average, the Shannon-Wiener diversity indices of phytoplankton community ranged from 1.17 to 1.78 (autumn > summer > spring > winter), and the Pielou’s evenness indices of phytoplankton ranged from 0.45 to 0.65 (autumn > spring > summer > winter). According to the results of canonical correspondence analysis, phosphate level was the major factor that limited the occurrence of P. sulcata and C. oculus-iridis, whereas optimal temperature and low salinity were responsible for Prorocentrum blooms in summer. The detailed description of seasonal variations in phytoplankton community structure in the three bays provide reference data for future studies on marine ecosystems and mariculture in adjacent areas.

    phytoplankton; seasonal variation; community structure; Sanggou Bay; Ailian Bay; Lidao Bay

    1 Introduction

    Phytoplankton are the primary producers in marine ecosystems, which absorb CO2and solar energy into organic compounds and produce energy through photosynthesis. Phytoplankton play an important role in global material cycle and energy flow (Sun, 2011), as well as carbon flux, cloud albedo, seawater heat flux, and luminous flux (Sun et al., 2011). In view of the fundamental link of marine food chain, phytoplankton are considered to be an important food resource of zooplankton and provide directly or indirectly bait for breeding aquatic animals, such as fish, shrimp and shellfish (Ren et al., 2010). The species composition and biomass conditions of phytoplankton directly affect their harvest production (Jin, 1965; Kang, 1986). Thus, investigation of phytoplankton community characteristics is of reference value for estimating the potential yield and managing reasonable development and utilization of the marine biotic resource (Cheng, 2004). Due to the small individual size and short life cycle, phytoplankton can rapidly respond to environmental changes. Thus, the phytoplankton community structure serves as an indicator of water quality (Shen and Shi, 2002).

    The Sanggou Bay (SGB, about 133.3 km2) is a semienclosed basin located in the eastern Shandong Peninsula and open to the Yellow Sea. The SGB is one of the most important mariculture bases in northern China (Song and Cui, 1996; Zhang et al., 2012) and is rich in economically important aquatic species such as seaweed, oysters, and scallops (Mao et al., 1988). In mariculture process, the excessive baits and the metabolites of farm-raised cultures deposited into sediment or suspended in water are the major components that cause water quality degradation in shallow waters (Wan, 2012). Along with the rapid development of local mariculture industry, the environ-ment of SGB and the adjacent water bodies are negatively impacted, as most obviously reflected by the changes in the phytoplankton community structure (Sun and Zhao, 1998; Zhu et al., 2000). Because of selective filter-feeding to dinoflagellates, the scallop aquaculture has changed the phytoplankton community composition in SGB (Zhang et al., 2005). Early surveys indicated that to some extent, scallop culture reduced the community diversity of phytoplankton in the SGB during 1983–2004 (Song et al., 2007). Additionally, the distribution of phytoplankton species was uneven while their biodiversity was at the medium level in the SGB during 2003–2005 (Mu et al., 2009). With the development of aquaculture, the phytoplankton community structure in the SGB has significantly changed in the past 30 years. Due to the significant role of phytoplankton as the primary marine producers, the changes in their community structure would inevitably impose a feedback effect on local farming. Deep understanding of phytoplankton community structure in the SGB would provide guiding information for reasonable management of aquaculture.

    Previous studies on phytoplankton in the SGB have primarily focused on tow-net samples (Song et al., 2007; Zhang et al., 2005), in which small-sized phytoplankton cells are easily omitted due to the limitation of the net sampling method with mesh aperture (Sun et al., 2002; Sun et al., 2003a). Despite the use of water samples in some surveys, such kind of work is limited to phytoplankton community characterization and rarely correlates phytoplankton parameters to local hydrological factors, thus blurring the relationship between the environment changes and the phytoplankton-related biological reactions (Mu et al., 2009). In the present study, we focused on the seasonal succession of phytoplankton communities and the changes in several environmental factors such as water temperature, salinity, and nutrient levels during the four seasons of a year in SGB and the adjacent Ailian Bay (ALB) and Lidao Bay (LDB). The results were analyzed to test the relationship between the phytoplankton community structure and the environmental factors, in order to provide guiding information for relevant studies on marine ecosystem and aquaculture.

    2 Material and Methods

    2.1 Study Area and Sampling Strategy

    A multidisciplinary research was carried out in the SGB and the adjacent ALB and LDB (37?00′–37?15′N, 122?30′–122?45′E) (Fig.1) during April 23–25, 2011 (spring), August 2–5, 2011 (summer), October 27–28, 2011 (autumn), 2011, and January 5–7, 2012 (winter). A total of 33 stations were planned and investigated during the 4 cruises.

    Water depth, temperature, salinity and dissolved oxygen (DO) were measured in situ with a calibrated automatic analyzer (YSI Professional Plus Multiparameters Water Quality Meter, US). Surface water samples (underwater 0.5 m) were collected using a Ruttner sampler (HYDRO-BIOS, Germany) at each station. Additional samples (depth 10 m) were collected at the stations with water depth >10 m. Phytoplankton samples (250 mL) were taken in polyethylene bottles and fixed with buffered formaldehyde (final conc. 2% by volume). 197 water samples were fixed in the investigation and stored at room temperature for laboratory analyses.

    Fig.1 Study area and sampling stations in the Sanggou, Ailian and Lidao bays.

    2.2 Laboratory Analyses

    Cell enumeration and species identification of phytoplankton were performed at ×200 or ×400 magnification under an inverted microscope after 24-h sedimentation of the fixed samples in a 25 mL Uterm?hl chamber (American Optical Ltd.) (Sun et al., 2002; Sun et al., 2003a). Phytoplankton species were identified according to Jin et al. (1965) and Isamu (1991).

    Water samples were filtrated with 0.45-μm-pore-sized acetate fiber membranes for chemical analyses. Ammonium (NH4-N), nitrate (NO3-N), and nitrite (NO2-N) levels were analyzed by using the sodium hypobromite/hypochlorite oxidation, zinc-alternate reduction, and naphthyl ethlenediamine spectrophotometric methods, respectively. Soluble reactive phosphate (SRP) and silicate concentrations were determined by phosphorus- and silicon-molybdenum blue reactions, respectively (Grassholf et al., 1983).

    2.3 Community Structure Analysis

    The phytoplankton community diversity was evaluated using the Shannon-Wiener (S-W) diversity index (H') (Shannon, 1951; Sun and Liu, 2004) as follows:

    where Piis the relative cell abundance of a species, and i and S are the numbers of the species. The evenness index (J) (Pielou, 1984) was calculated from H' (Sun and Liu, 2004) as follows:

    where H' is the S-W diversity index, and S is the number of the species in a sample. The phytoplankton dominance index (Y) was calculated using the following formula (Sun et al., 2003b):

    where niis the number of the individual species, N is the total number of all species, and fiis the occurrence frequency of the species in a sample. The Jaccard similarity index (P) was calculated using the following formula (Jaccard, 1908):

    where a and b are the numbers of species in 2 different seasons, and c is the number of the common species in the 2 seasons.

    Canonical correlation analysis (CCA) (ter Braak, 1986) was performed to analyze the relationship between phytoplankton cell abundance and related environmental factors. Biplots of different seasons were made with the MVSP3.1 program package (Kovach Computing Services, UK).

    3 Results

    3.1 Hydrology

    The hydrological characteristics of SGB, LDB, and ALB, including water depth, temperature, salinity, pH, and nutrient levels, were investigated. The water depths of the 3 bays range from 5 m to 40 m and average approximately 16.8 m. As located in a typical temperate coastal area, the 3 bays have an average water temperature of 4.4℃ in winter and 20.2℃ in summer. There are small differences in the average salinity among 4 seasons, ranging from 30.41 in summer to 31.54 in winter. The average transparency depth varies seasonally from 0.67 m in winter to 1.49 m in spring. The maximum pH value 8.68 occurs in spring and the minimum pH value 8.24 is in autumn. The nitrogen level is the highest in autumn (average 0.377 mg L-1), followed by summer (average 0.254 mg L-1); the lowest value appears in spring and winter (average 0.244 mg L-1). The phosphate levels of bay waters remain low from spring to winter, i.e., 0.009, 0.005, 0.003, and 0.002 mg L-1on average, respectively. The average silicate level of bay waters is approximately 10-fold higher in winter (up to 3.453 mg L-1) than in spring (0.322 mg L-1), summer (0.199 mg L-1), and autumn (0.278 mg L-1).

    3.2 Phytoplankton Species Composition

    A total of 80 taxa belonging to 39 genera of 3 Phyla were identified in 197 phytoplankton water samples. These included 64 species of 30 genera in the Phylum Bacillariophyta (76.92% of total species number); 13 species of 8 genera in the Phylum Dinophyta (20.51%); and 3 species of 1 genus were Chrysophyta (2.57%). The species composition of phytoplankton community followed the temperate coastal and widely distributed regimes. The total number of phytoplankton species was greater in spring than in the other 3 seasons. The number of diatom species accounted for >50% of the total species number in 4 seasons, showing similar seasonal changes to those with phytoplankton. By comparison, dinoflagellates were most diverse in summer but were still less diverse than diatoms (Fig.2). Chrysophyta only appeared in summer and autumn with very few species.

    Fig.2 Comparison of the numbers of phytoplankton species in the Sanggou, Ailian and Lidao bays among 4 seasons.

    There were 20–25 identical species among 4 seasons in the surveyed area, of which 15 species appeared in all of the four seasons, including 14 species of diatom and 1 species of dinoflagellate. The Jaccard similarity index values were 0.317–0.500 (Table 1) with the highest level detected between autumn and winter. The above results suggest that the phytoplankton community structure obviously changes with seasons.

    In the surveyed area, the diatoms Paralia sulcata (Ehrenberg) Cleve and Coscinodiscus oculus-iridis Ehrenberg were most common and predominant in four seasons (Table 2). P. sulcata was the primary dominant species that accounted for 36.60% of the phytoplankton abundance in summer and >50% in the other seasons (up to 72.01% in winter). Other dominant diatoms species such as Thalassiosira sp. and Nitzschia sp. frequently occurred in spring, autumn, and winter, whereas dinoflagellates rather than diatoms were predominant in summer. Two dominant dinoflagellates Prorocentrum gracile and Prorocentrum triestinum accounted for 17.41% and 20.92% of the phytoplankton abundance, respectively. Both species were mainly concentrated at the Station 19. In spring, a large number of Detonula pumila appeared at the StationL4 only, accounting for 96.23% of the phytoplankton abundance. In winter, the dominant species were diatoms. The occurrence frequency of P. sulcata in winter (90%) was the highest among seasons.

    Table 1 The Jaccard similarity index of phytoplankton community in the Sanggou, Ailian and Lidao bays between seasons

    Table 2 Dominant species of phytoplankton in the Sanggou, Ailian and Lidao bays

    3.3 Phytoplankton Cell Abundance

    The cell abundance of phytoplankton ranged (0.15–188.80) × 103cells L-1(average 5.84 × 103cells L-1) and peaked in summer (seasonal average 8.08 × 103cells L-1), followed by spring (6.62 × 103cells L-1) and winter (6.07 × 103cells L-1). The lowest cell abundance of phytoplankton was found in autumn (2.60 × 103cells L-1) (Fig.3A). The cell abundance of diatom ranged from 0 to 76.67 × 103cells L-1(average 4.77 × 103cells L-1) and peaked in spring (seasonal average 6.61 × 103cells L-1), followed by winter (6.07 × 103cells L-1) and summer (3.85 × 103cells L-1). The lowest cell abundance of diatom was detected in autumn (2.55 × 103cells L-1) (Fig.3B). The cell abundance of dinoflagellate was 0–185.87 × 103cells L-1(average 1.07 × 103cells L-1), with average value peaking in summer (seasonal average 4.22 × 103cells L-1) and being ≤0.05 × 103cells L-1in other seasons (lowest in winter, 0.003 × 103cells L-1) (Fig.3C). Over the four seasons, the relative abundance of diatom was 99.96% and 99.76% of total phytoplankton abundance in winter and spring, respectively, followed by 98.05% in autumn and only 47.67% in summer. The relative abundance of dinoflagellate was 52.27% of the total phytoplankton abundance in summer and <2% in other seasons. Chrysophyta merely appeared in summer and autumn in low cell abundance.

    Fig.3 Box-Whisker plots of phytoplankton cell abundance (× 103cells L-1) in the Sanggou, Ailian and Lidao bays in the four seasons. Whiskers mean minimum and maximum values; boxes display the lower and upper quartiles; and lines mean median-sets. A, phytoplankton; B, diatom; and C, dinoflagellate.

    3.4 Horizontal Distribution of Phytoplankton Cell Abundance

    3.4.1 Surface water

    The horizontal distribution of phytoplankton cell abundance in surface water of the bays over the four seasons is shown in Fig.4. In spring, high-abundance phytoplankton were mainly concentrated in the north of central LDB and sparsely distributed in coastal waters. The maximum cell abundance of phytoplankton (76.67 × 103cells L-1) was observed at Station L4, where Detonula pumila was predominant, followed by P. sulcata and C. oculus-iridis. In summer, high-abundance phytoplankton were mainly distributed in the south of inner SGB with the abundance of (0.31–188.80) × 103cells L-1. A peak of the phytoplankton abundance was observed at Station 19 where the dinoflagellates P. gracile and P. triestinum were dominant. In autumn, the phytoplankton were concentrated in southern SGB, where cell abundance decreased toward the north and the inner bay. Maximum cell abundance of phytoplankton was found at Station 5 (9.64 × 103cells L-1) where the dominant species P. sulcata accounted for 74.05% of total phytoplankton abundance. In winter, the phytoplankton abundance showed 3 patches with the peak level (19.41 × 103cells L-1) observed at Station 18 where Skeletonema sp. was the dominant species (77.67% of total phytoplankton abundance).

    Fig.4 Horizontal distribution of phytoplankton cell abundances (× 103cells L-1) in surface water of the Sanggou, Ailian, and Lidao bays.

    The distribution of diatom and dinoflagellate cell abundance in surface water of the bays during the four seasons is shown in Fig.5. The high-abundance diatoms were concentrated in middle LDB in spring and moved southward in summer, toward southern LDB, outer ALB, and northern SGB. In autumn, the diatoms continued moving southward and were concentrated in the south of outer SGB. In winter, the diatoms were distributed in 3 patches in inner SGB, outer SGB, and middle ALB, respectively. High-abundance dinoflagellates were concentrated in the LDB in spring and autumn, the south of inner SGB in summer, and northern SGB and the south of middle SGB in winter.

    In spring, the diatom cell abundance peaked at Station L4 (76.67 × 103cells L-1) where Detonula pumila was dominant. The corresponding dinoflagellate cell abundance averaged 0.02 ×103cells L-1and peaked at Station L10. In summer, the diatoms were concentrated from southern LDB to northern SGB (3.28 × 103cells L-1on average) with the peak cell abundance at Station L8 (10.27 × 103cells L-1). At this station, the dinoflagellate abundance was the highest among seasons (average 6.11 × 103cells L-1), being 52.27% of the phytoplankton cell abundance. At Station 19, the phytoplankton cell abundance reached the peak value and Prorocentrum was the major contributor. The distribution of dinoflagellate cell abundance was similar to that of total phytoplankton in surface water. In autumn, the diatom cell abundance wasthe lowest among four seasons, averaging 2.19 × 103cells L-1. The diatoms were concentrated in the south of outer SGB while the dinoflagellates were concentrated in the south of middle LDB. In winter, the distribution of diatom cell was relatively even in the surveyed area. The diatom cell abundance peaked at Station 18 (19.41 × 103cells L-1), whereas the corresponding dinoflagellate cell abundance was generally low and averaged 0.004 × 103cells L-1, that is, the lowest among the four seasons. In general, except in summer, diatoms had the absolute ad-vantage of growth and were characterized by the distribution of phytoplankton cell abundance in surface water.

    Fig.5 Horizontal distribution of diatom (Left) and dinoflagellate (Right) cell abundances (× 103cells L-1) in surface water of the Sanggou, Ailian, and Lidao bays.

    3.4.2 Bottom water (10 m)

    In spring, the phytoplankton in bottom water were concentrated in outer ALB and sparsely distributed in coastal waters. The phytoplankton abundance peaked at Station A3 where P. sulcata was predominant (83.92% of phytoplankton cell abundance). In summer, the abundance distribution of phytoplankton in bottom water was similar to that in surface water, that is, being within the rang of (0.58–188.80) × 103cells L-1(average 9.64 × 103cells L-1) and mainly concentrated in the south of inner SGB. In autumn, the phytoplankton were concentrated around the center of LDB and ALB, with the maximum abundance detected at Station L3 (6.58 × 103cells L-1) where P. sulcata and C. oculus-iridis were the dominant species. In winter, the phytoplankton cell abundance ranged (1.63–19.41) × 103cells L-1(average 6.00 × 103cells L-1), with the maximum level detected at Station 18 (Fig.6). At this station, Skeletonema sp. was the dominant species, consistent with the case of surface water.

    The horizontal distribution of diatom and dinoflagellate cell abundance in bottom water of the bays is shown in Fig.7. The distribution of diatom abundance in bottom water was different from that in surface water, especially in spring. The diatom abundance was higher in the outer bay than in the inner bay in spring. The diatoms were concentrated in the north of offshore SGB in summer, outer LDB and ALB in autumn, and inner SGB, outer SGB and outer ALB in winter. The distribution of dinoflagellate abundance in bottom water was similar to that in surface water in summer and winter. The dinoflagellates were concentrated mostly in the outer LDB in spring and autumn.

    In spring, the cell abundance of diatoms peaked at Station A3 where P. sulcata followed by C. oculus-iridis were the predominant species. The corresponding cell abundance of dinoflagellates peaked at Station L9 (0.01 × 103cells L-1), being lower than that in surface water. In summer, the diatom cell abundance was 3.53 × 103cells L-1on average, with the maximum level (13.67 × 103cells L-1) observed at Station 2 where P. sulcata and C. oculus-iridis were the predominant diatoms. The dinoflagellate cell abundance peaked at Station 19, consistent with the case of surface water. In autumn, the diatom cell abundance averaged 2.36 × 103cells L-1, being lower than that in surface water, whereas the dinoflagellate cell abundance averaged 0.04 × 103cells L-1, being higher than that in spring and winter. The diatom and dinoflagellate cell abundances peaked at stations L3 and L6, respectively. In winter, the diatom and dinoflagellate cell abundances showed similar distribution patterns in bottom and surface water layers. On the whole, the phytoplankton were dominated by dinoflagellates in summer and diatoms in other seasons.

    Fig.6 Horizontal distribution of phytoplankton cell abundance (× 103cells L-1) in bottom water of the Sanggou, Ailian, and Lidao bays.

    Fig.7 Horizontal distribution of diatom (Left) and dinoflagellate cell abundances (Right) (× 103cells L-1) in bottom water of the Sanggou, Ailian, and Lidao bays.

    3.5 S-W Diversity Index and Pielou’s Evenness Index

    The diversity indices of phytoplankton in four seasons ranged from 1.17 to 1.78 with an average of 1.41 (autumn > summer > spring > winter, seasonal average). The evenness indices of phytoplankton ranged from 0.45 to 0.65 with an average of 0.54 (autumn > spring > summer > winter, seasonal average). The S-W diversity index of phytoplankton community was evenly distributed in spring and was relatively lower in ALB and LDB than in SGB insummer and winter, whereas the evenness index of phytoplankton community was roughly lower in ALB and LDB than in SGB in the four seasons. Overall, the relatively higher S-W diversity index of phytoplankton appeared in the north of inner SGB, decreasing toward the outer bay. Pielou’s evenness index of phytoplankton was higher in coastal waters of the bay than in the offshore area (Fig.8). These trends suggested that phytoplankton were evenly distributed in high diversity in the surveyed area. Comparisons of the distribution between cell abundance and community diversity indices of phytoplankton in surface water showed that the high phytoplankton diversity index was associated with medium-level phytoplankton cell abundance.

    Fig.8 The distribution of Shannon-Wiener diversity index (Left) and Pielou’s evenness index (Right) of phytoplankton in surface water of the Sanggou, Ailian, and Lidao bays.

    3.6 CCA

    The relative positions of dominant phytoplankton species in the CCA bioplots (Fig.9) reflect the dependence of the species on environmental factors. Due to different ecological adaptabilities of phytoplankton species, the more similar species are on taxonomy, the closer their positions are in the diagram. In spring, the abundance of P. sulcata, C. oculus-iridis, T. rotula, and T. pacifica were positively correlated with phosphate, nitrite, nitrate, and ammonium levels. The abundance of Detonula pumila was most related to silicate level and depth. These results indicate that nutrient level was an important factor that impacted the phytoplankton community structure in spring. There were 2 categories of the correlation between the cell abundance of phytoplankton species and the levels of environmental variables in summer. The cell abundance of P. gracile and P. triestinum of dinoflagellate were positively correlated with nitrate level, pH, temperature, DO level, and silicate level and negatively correlated with salinity. The cell abundance of diatoms such as P. sulcata, Donkinia recta, and C. oculus-iridis was correlated with phosphate level, salinity, and water depth. The cell abundance of Nitzschia sp. correlated with nitrite and ammonium levels. In autumn, there were positive correlations between P. sulcata cell abundance and water depth, temperature, salinity, phosphate level, whereas the abundance of Donkinia recta, Nitzschia sp., Thalassiosira sp. and Pleurosigma angulatum were correlated with nitrite and nitrate levels, pH and DO level. In winter, P. sulcata and C. oculus-iridis cell abundance were positively correlated with temperature, salinity and phosphate level; Nitzschia sp. cell abundance was correlated with pH and DO level.

    Fig.9 Biplots of phytoplankton species and environmental factors in the Sanggou, Ailian, and Lidao bays in four seasons. T: temperature; S: salinity; DO: dissolved oxygen; dep: depth; NO2: nitrite; NO3: nitrate; NH3: ammonium; P: phosphate; and Si: silicate. Ps: Paralia sulcata; Co: Coscinodiscus oculus-iridis; Tp: Thalassiosira pacifica; Dp: Detonula pumila; Na: Navicula sp.; Dr: Donkinia recta; Tr: Thalassiosira rotula; Pg: Prorocentrum gracile; Pt: Prorocentrum triestinum; Th: Thalassiosira sp.; Ni: Nitzschia sp.; Pa: Pleurosigma angulatum; Pd: Pseudonitzschia delicatissima; Sp: Stephanopyxis palmeriana; and Sk: Skeletonema sp.

    4 Discussion

    4.1 Comparison with Historical Data

    The results from this study showed that diatoms, followed by dinoflagellates, were constantly dominant in the phytoplankton community in the SGB. Compared with previous data, the species composition of phytoplankton remained widespread and mainly in the temperate nearshore zone, while the number of phytoplankton species gradually declined in the past 30 years (Table 3). Human activities such as varieties of aquaculture model and their quantity expanding might have changed the marine ecological environment and further caused seasonal variations in the phytoplankton community structure.

    During early breeding in the SGB, the cell abundance of phytoplankton peaked in July (summer) and October (autumn) (Mao et al., 1988). There were 2 peaks in late spring and winter in the 1990s (Fang et al., 1996). In the present survey, the phytoplankton cell abundance peaked in summer, consistent with the findings of Mao et al.(1988). Previously, Chaetoceros, Coscinodiscu, Nitzschia of Bacillariophyta and Ceratium of Dinophyta were found predominant in the SGB (Mu et al., 2009), while in the present study, Paralia, Coscinodiscus, Thalassiosira of Bacillariophyta and Prorocentrum of Dinophyta were the predominant phytoplankton species.

    Table 3 Changes of the species composition of phytoplankton in the Sanggou, Ailian, and Lidao bays (1983–2012)

    4.2 Seasonal Changes in Phytoplankton Cell Abundance and Species Composition

    The phytoplankton community in the survey area was mainly composed of diatoms and dinoflagellates, consistent with previous findings (Mao et al., 1988; Mu et al., 2009; Song et al., 2007). At individual sites, there were a small number of Chrysophyta (Distephanus speculum). The diatoms such as P. sulcata, C. oculus-iridis and Thalassiosira sp. were dominant in terms of species number and cell abundance. Only in summer dinoflagellates had an explosive growth at a few stations, especially Station 19. P. gracile and P. triestinum were the major species with the maximum cell abundance 100.3 × 103cells L-1and 78.3 × 103cells L-1, respectively. Prorocentrum was the major genus of red tide in the world (Huang et al., 1997). Although it did not reach the threshold number of red tide, the increase in Prorocentrum cell abundance in summer indicated that the rapid growth and reproduction of dinoflagellates were likely to cause the outbreak of red tide.

    Over the four seasons, the highest cell abundance of phytoplankton was detected in summer, followed by spring and winter; the lowest cell abundance of phytoplankton was found in autumn (Table 4). The seasonal variations in the phytoplankton abundance could be explained by several reasons. In spring (April), the repression of phytoplankton growth and reproduction from the competition for nutrients was weak and there were nutriments accumulated in winter. Together the suitable environmental conditions allowed the cell abundance of phytoplankton to reach a second peak. In summer (August), the harvest of kelp reduced the absorption of nutrients and the decomposing of residual kelp added substantial nutrients to the water column (Liu et al., 2003; Shi et al., 2010a), plus that the light and temperature were suitable for phytoplankton growth (Shi et al., 2010b). These conditions might promote phytoplankton blooms. In autumn (October), shellfish grew rapidly and consumed a large amount of phytoplankton while the phytoplankton growth was limited by light. Thus, the phytoplankton cell abundance declined drastically in this period (Shi et al., 2010b). Coming to winter (January), the harvest of shellfish weakened the pressure of feeding on phytoplankton whose cell abundance slightly recovered.

    Table 4 Mean values and ranges of phytoplankton cell abundance (× 103cellsL-1) at the stations in the Sanggou, Ailian, and Lidao bays during the four seasons

    Analyses for the different stations with their both surface and bottom water showed that (Table 4): in winter, the changes in the phytoplankton cell abundance among different sites were the least with the highest cell abundance being 11.91-fold that of the lowest value. In spring and summer, the maximum-to-minimum ratio of cell abundance was up to 517.50–609.03, indicating uneven distribution of phytoplankton. Farming activities including seeding and harvesting of aquaculture species could have changed the physicochemical environment of the waters, and in different farming areas the aquaculture species limited the growth and reproduction of phytoplankton to different degrees. These could be the main causes for the uneven distribution of phytoplankton in the study area in spring and summer. Additionally, highdensity breeding might have impacted the flow of sea water and impacted the transport of biogenic elements, resulting in local deficiency of water (Shi et al., 2010b) and further affecting the distribution of phytoplankton.

    4.3 Seasonal Changes in Phytoplankton Community Diversity

    The variations in phytoplankton community structure to some extent reflect the changes in water conditions, and thus can be used to evaluate the status of the aquatic environment (Mu et al., 2009). The S-W diversity index indicates the complexity of the community structure; the evenness index reflects the uniformity degree of the distribution of species as well as the maturity and stability of the community (Sun and Liu, 2004). Although the cell abundance of phytoplankton peaked in summer, their diversity was not high at all because the predominant species Prorocentrum accounted for 52.27% of the total phytoplankton cell abundance. In autumn, the cell abundance was the lowest while the diversity index was thehighest, verifying that the low phytoplankton cell abundance was associated with high species diversity. The S-W diversity index <1 and the evenness index <0.3 were used as the standard of poor diversity (Ma et al., 2002). According to the above standard, the phytoplankton community structure in the surveyed area was in a relatively stable state, as reflected by the S-W diversity index >1 and the evenness index >0.3.

    Generally, anthropogenic activities in nearshore areas would decrease the phytoplankton diversity. However, in the present survey, the coastal area as more affected by human activities had higher phytoplankton diversity than nearshore area, possibly due to the effect of the specific aquaculture species. Shellfish’s non-selectively filtering of phytoplankton made each species in a dynamic equilibrium and no one species had absolute advantage. Therefore, the inner bay associated with shellfish breeding showed the higher diversity of phytoplankton (Shi et al., 2010a). In kelp farming area at the bay mouth, low nutrient levels related to kelp’s consumption limited the growth and reproduction of phytoplankton during most of the year. Such species at a competitive advantage occupied the entire phytoplankton community, leading to the decrease in the phytoplankton biodiversity.

    4.4 Relationship Between Phytoplankton and Environmental Factors

    The results of CCA showed that there were discrepancies in the correlation between the cell abundance of dominant phytoplankton species and the levels of specific environment factors over the four seasons. The dominant species of diatom, such as P. sulcata and C. oculus-iridis, were mainly restricted by phosphate level. P. Sulcata, as a benthic species, was well correlated with salinity. The dominant species of dinoflagellates that appeared only in summer were associated with nitrate level, silicate level and pH value. In summer, the dinoflagellate cell abundance greatly increased with the nutrient-competitive advantage at lower concentrations of phosphate and silicate. In winter, the low temperature was not suitable for the growth of phytoplankton. Thus, temperature was the main factor that limited phytoplankton growth in this season.

    During the investigation period, summer was favorable for explosive growth of Prorocentrum. This is understandable because water temperature and the salinity in the surveyed area were 18.1–24.2℃ and 26.57–31.43, respectively, both being in the ranges of optimal growth conditions (Zhang and Liu, 2007). The low salinity in summer was more conducive to reproduction of phytoplankton because the growth energy was relatively less under such condition (Chen et al., 2005; Chen et al., 2006; Zhang and Liu, 2007). The total nutrient level was the lowest in summer when Prorocentrum was more competitive in the low-nutrient environment. The sufficient irradiance, suitable temperature, and low salinity and nutrient levels were favorable for explosive growth of dinoflagellates, which increased the possibility of red tide outbreaks. Once the red tide breaks out, local fisheries and aquaculture would be threatened, and even serious environmental problems and economic losses would be caused. Environmental monitoring should be ready for preventing the occurrence of red tide during the summer period.

    5 Conclusions

    Our preliminary investigation of phytoplankton community structure showed that diatoms and dinoflagellates were the major groups of phytoplankton in the SGB and the adjacent ALB and LDB. The diatom species P. sulcata and C. oculus-iridis were common and dominant over the four seasons, whereas the dinoflagellate Prorocentrum dominated the phytoplankton community in summer only. Compared to previous data (1983–1984), the present study employed different sampling strategy and microscopic examination methods. Still, a similar decreasing trend was observed in the number of phytoplankton species in the SGB area. Results of the CCA suggested that the abundance of predominant species P. sulcata and C. oculus-iridis were mainly limited by phosphate level in water, while the warm temperature and low salinity were primarily responsible for the dinoflagellate bloom in summer. More recently, a brown-watered picoplankton has been widely spreading and thriving in the SGB, which could heavily affect the adjacent marine ecosystem. Further monitoring and research are needed for investigating the implication of this picoplankton that will undoubtedly be coupled with the mariculture and trigger new marine issues.

    制定智慧交通管理系統(tǒng)后臺系統(tǒng)的數(shù)據(jù)結(jié)構(gòu)、數(shù)據(jù)通訊協(xié)議、數(shù)據(jù)查詢接口等標(biāo)準(zhǔn)和規(guī)范,提高智慧交通管理系統(tǒng)的開放性和兼容性.

    Acknowledgements

    This research was supported by the National Program on Key Basic Research Project of China (Grant Nos. 2011CB409804 and 2015CB954002), Program for New Century Excellent Talents in University (NCET-12-1065) and the National Natural Science Foundation of China (Grant No.41176136) to J. Sun. The authors thank Dr. Y. Yu (Institute of Oceanology, Chinese Academy of Sciences) for sampling assistance in the cruises.

    Chen, B. Z., Wang, Z. L., Zhu, M. Y., and Li, R. X., 2005. Effects of temperature and salinity on growth of Prorocentrum dentatum and comparisons between growths of Prorocentrum dentatum and Skeletonema costatum. Advances in Marine Science, 23 (1): 60-64.

    Chen, H. L., Lv, S. H., Zhang, C. L., and Zhu, D. D., 2006. A Survey on the red tide of Prorocentrum doghaiense in East China Sea, 2004. Ecologic Science, 25 (3): 226-230.

    Cheng, J. S., 2004. The Ecological Environment and Biological Community in Inshore Waters of Huanghai and Bohai Seas. Press of Ocean University of China, Qingdao, 127-162.

    Fang, J. G., Kuang, S. H., Sun, H. L., Sun, Y., Zhou, S. L., Song, Y. L., Cui, Y., Zhao, J., Yang, Q. F., Li, F., Zhang, A. J., Wang, X. Z., and Tang, T. Y., 1996. Study on the carrying capacity ofSanggou Bay for the culture of scallop chlamys. Marine Fisheries Reseach, 17 (2): 18-31.

    Grassholf, K., Ehrhardt, M., and Kremling, K., 1983. Methods of Seawater Analysis. 2nd edition, Chemie Weinheim, 1-419.

    Huang, X. H., Shi, D. M., and Zhang, J. P., 1997. Study on nitrogen meiabolism of Prorocentrum micans with special referece to the occurrence of red tide. Oceanologia et Limnologia Sinica, 289 (1): 33-38.

    Isamu, Y., 1991. Illustrations of the Marine Plankton of Japan. 3rdedition. Hoikushapublishing, Tokyo, 1-158.

    Jaccard, P., 1908. Nouvelles recherches sur la distribution florale. Bulletin de la Société Vaudense des Sciences Naturelles, 44: 223-270.

    Jin, D. X., 1965. The introduction of nanoplankton. Journal of Fishery of Fujian, 1-8.

    Jin, D. X., Chen, J. H., and Huang, K. G., 1965. Chinese Marine Planktonic Diatoms. Shanghai Science Press, Shanghai, 1-230.

    Kang, Y. D., 1986. The ecological characteristics of phytoplankton and the relationship with the fishery in Yellow sea. Marine Fisheries Research, 7: 103-107.

    Liu, H., Fang, J. G., Dong, S. L., Liang, X. M., Jiang, W. W., Wang, L. C., and Lian, Y., 2003. Study on phytoplankton in raft-culture areas of Laizhou Bay and Sangou Bay II. Marine Fisheries Research, 24 (3): 20-28.

    Ma, J. X., Zheng, Z., Li, Y. P., Xing, H. Y., and Liu, Y. H., 2002. The distribution characteristics of phytoplankton in Laizhou Bay. Transaction of Oceanology and Limnology, 4: 64-67.

    Mao, X. H., Zhang, W. X., and Zhang, J. Z., 1988. Comprehensive Investigation of Environment for Propagation and Multiplication in Sanggou Bay. Qingdao Press, Qingdao, 8-46, 113-148.

    Mu, J. D., Dong, W., and Chen, B. J., 2009. Ecological characteristics of phytoplankton in Sanggou Bay. Progress in Fishery Sciences, 30 (3): 91-96.

    Pielou, E. C., 1984. The Interpretation of Ecological Data. John Wiley & Sons, New York, 1-230.

    Ren, Y., C., Wang, F., Dong, S., L., and Liu, F., 2010. Seasonal characteristics of primary production of sea cucumber (Apostichopus japonicus) culture ponds in Jinghai Bay, Rongcheng. Periodical of Ocean University of China, 40 (3): 24-28.

    Shannon, C. E., 1951. A mathematical theory of communication. Bell System Technical Journal, 27: 379-423.

    Shen, G. Y., and Shi, B. Z., 2002. Marine Ecology. 2nd edition. Science Press, Beijing, 191-192.

    Shi, J., Wei, H., Zhao, L., Fang, J. G., and Zhang, J. H., 2010a. Study on ecosystem model of multi-species culture in Sanggou Bay: ? Estblishment of the culture ecosystem model and sensitivity analyses of the parameters. Progress in Fishery Sciences, 31 (4): 26-35.

    Shi, J., Wei, H., Zhao, L., Fang, J. G., and Zhang, J. H., 2010b. Study on ecosystem model of multi-species culture in Sanggou Bay: II Simulation of ecosystem and the circulation of nutrient. Progress in Fishery Sciences, 31 (4): 36-42.

    Song, H. J., Li, R. X., and Wang, Z. L., 2007. Interannual variations in phytoplankton diversity in the Sanggou Bay. Advances in Marine Science, 25 (3): 332-339.

    Song, Y. L., and Cui, Y., 1996. The analysis of nutritional status and its affecting factors in cultured areas of Sanggou Bay. Marine Fisheries Research, 17 (2): 41-51.

    Sun, J., 2011. Marine phytoplankton and biological carbon sink. Acta Ecologica Sinica, 31 (18): 5372-5378.

    Sun, J., Li, X., and Tian, W., 2011. Summer phytoplankton community in the Beibu Gulf. Transactions of Oceanology and Limnology, 3: 1-7.

    Sun, J., Liu, D. Y., and Qian, S. B., 2002. A quantative research and analysis method for marine phytoplankton: An introduction to Uterm?hl method and its modification. Journal of Oceanography of Huanghai & Bohai Seas, 20 (2): 105-112.

    Sun, J., Liu, D. Y., and Feng, S. Z., 2003a. Preliminary study on marine phytoplankton sampling and analysis strategy for ecosystem dynamic research in coastal waters. Oceanologia et Limnologia Sinica, 34 (2): 140-148.

    Sun, J., Liu, D. Y., Ning, X. R., and Liu, C. G., 2003b. Phytoplankton in the Prydz Bay and the adjacent Indian sector of the Southern Ocean during the austral summer 2001/2002. Oceanologia et Limnologia Sinica, 34 (5): 519-532.

    Sun, J., and Liu, D. Y., 2004. The application of diversity indices in marine phytoplankton studies. Acta Oceanologica Sinica, 26 (1): 62-75.

    Sun, Y., and Zhao, J., 1998. Water environmental characteristics in cultured areas of Sanggou Bay. Journal of Fishery Sciences of China, 5 (3): 69-75.

    Ter Braak, C. J. F., 1986. Canonical correspondence analysis: A new eigenvector technique for multivariate direct gradient analysis. Ecology, 67 (5): 1167-1179.

    Wan, L., 2012. Effect of shellfish farming on nutrient salts of seawater in Sanggou Bay in spring. Environmental Science and Management, 37 (6): 62-66.

    Zhang, J. H., Wang, W., Han, T. T., Liu, D. H., Fang, J. G., Jiang, Z. J., Liu, X. J., Zhang, X. J., and Lian, Y., 2012. The distributions of dissolved nutrients in spring of Sangou Bay and potential reason of outbreak of red tide. Journal of Fisheries of China, 36 (1): 132-139.

    Zhang, L. H., Zhang, X. L., Li, R. X., Wang, Z. L., Li, Y., Wang, L. C., Lian, L., and Liu, Y., 2005. Impact of scallop culture on dinoflagellate abundance in Sanggou Bay. Advances in Marine Science, 23 (3): 342-346.

    Zhang, X. F., and Liu, Y. J., 2007. Advances in the biological study of dinoflagellate Prorocentrum donghaiense Lu. Ecology and Environment, 16 (3): 1053-1057.

    Zhu, M. Y., Zhang, X. L., Li, R. X., and Chen, S., 2000. Impacts of shellfish culture on the coastal ecosystem. Journal of Ocean University of Qingdao, 30 (2): 53-57.

    (Edited by Qiu Yantao)

    (Received March 1, 2013; revised April 24, 2013; accepted August 29, 2014)

    ? Ocean University of China, Science Press and Spring-Verlag Berlin Heidelberg 2014

    * Corresponding author. Tel: 0086-22-60601116

    E-mail: phytoplankton@163.com

    猜你喜歡
    數(shù)據(jù)通訊交通管理數(shù)據(jù)結(jié)構(gòu)
    主動交通管理在智慧高速中的應(yīng)用探討
    機(jī)器視覺系統(tǒng)與PLC的數(shù)據(jù)通訊
    “翻轉(zhuǎn)課堂”教學(xué)模式的探討——以《數(shù)據(jù)結(jié)構(gòu)》課程教學(xué)為例
    高職高專數(shù)據(jù)結(jié)構(gòu)教學(xué)改革探討
    中國市場(2016年45期)2016-05-17 05:15:48
    對設(shè)區(qū)城市交通管理立法的思考
    中國市場(2016年12期)2016-05-17 05:10:05
    基于C#的船舶電力推進(jìn)監(jiān)控系統(tǒng)的數(shù)據(jù)通訊設(shè)計
    限行不限購應(yīng)是小汽車交通管理的新常態(tài)
    智能交通是城市交通管理的有效途徑
    智能電能表數(shù)據(jù)通訊安全性分析
    TRIZ理論在“數(shù)據(jù)結(jié)構(gòu)”多媒體教學(xué)中的應(yīng)用
    又黄又爽又刺激的免费视频.| 亚洲美女视频黄频| 国产精品免费一区二区三区在线| 国产高清有码在线观看视频| 日韩av在线大香蕉| 麻豆成人午夜福利视频| 老师上课跳d突然被开到最大视频| 日韩,欧美,国产一区二区三区 | 高清毛片免费观看视频网站| 特大巨黑吊av在线直播| 不卡视频在线观看欧美| 国产三级在线视频| 偷拍熟女少妇极品色| 乱系列少妇在线播放| 精品久久久久久久久久久久久| 少妇熟女aⅴ在线视频| 午夜福利在线在线| 日韩强制内射视频| 欧美xxxx性猛交bbbb| 联通29元200g的流量卡| 成年女人看的毛片在线观看| 国产高清视频在线观看网站| 成人性生交大片免费视频hd| a级一级毛片免费在线观看| 国产午夜福利久久久久久| 九九久久精品国产亚洲av麻豆| 一夜夜www| 欧美色欧美亚洲另类二区| 免费人成在线观看视频色| 亚洲专区国产一区二区| 五月伊人婷婷丁香| 91久久精品国产一区二区成人| 国产精品久久久久久久久免| 成人永久免费在线观看视频| 村上凉子中文字幕在线| 亚洲av一区综合| 日本黄色片子视频| 亚洲国产精品合色在线| 男女做爰动态图高潮gif福利片| 欧美一区二区亚洲| 久久草成人影院| 91狼人影院| 日韩欧美国产一区二区入口| 在线观看美女被高潮喷水网站| 久久久久久久精品吃奶| 国产精品久久电影中文字幕| .国产精品久久| 午夜a级毛片| 一区二区三区四区激情视频 | 成人性生交大片免费视频hd| 国产伦精品一区二区三区四那| av在线蜜桃| 欧美精品国产亚洲| 国产精品嫩草影院av在线观看 | 精品不卡国产一区二区三区| 日本精品一区二区三区蜜桃| 熟女人妻精品中文字幕| 窝窝影院91人妻| 永久网站在线| 天堂动漫精品| 51国产日韩欧美| 国产av麻豆久久久久久久| 夜夜看夜夜爽夜夜摸| 日本色播在线视频| 一个人免费在线观看电影| 女人被狂操c到高潮| 色5月婷婷丁香| av中文乱码字幕在线| or卡值多少钱| 91麻豆av在线| 国产乱人伦免费视频| 国产一区二区在线观看日韩| 日韩高清综合在线| 国产探花在线观看一区二区| 亚洲精品在线观看二区| 天美传媒精品一区二区| 欧美日韩亚洲国产一区二区在线观看| 在线观看66精品国产| 99久久精品热视频| 精品欧美国产一区二区三| 人妻丰满熟妇av一区二区三区| 成年女人看的毛片在线观看| 久久久久久大精品| 久久精品国产亚洲网站| 两性午夜刺激爽爽歪歪视频在线观看| 精品人妻视频免费看| 精品免费久久久久久久清纯| 免费黄网站久久成人精品| 精品福利观看| 草草在线视频免费看| 69av精品久久久久久| 一卡2卡三卡四卡精品乱码亚洲| 又黄又爽又免费观看的视频| 国产精品国产三级国产av玫瑰| 桃红色精品国产亚洲av| 成人精品一区二区免费| 日本色播在线视频| 一本久久中文字幕| 亚洲中文日韩欧美视频| 91av网一区二区| 日本爱情动作片www.在线观看 | 成人国产综合亚洲| 国产精品,欧美在线| 久久国内精品自在自线图片| 麻豆一二三区av精品| 波多野结衣高清无吗| 免费看美女性在线毛片视频| 久久久久久久久久黄片| 一区二区三区四区激情视频 | 中出人妻视频一区二区| 老熟妇乱子伦视频在线观看| 男女做爰动态图高潮gif福利片| 日韩高清综合在线| 日韩欧美在线乱码| 亚洲成人精品中文字幕电影| 日本免费a在线| 亚洲精品456在线播放app | 在现免费观看毛片| av在线亚洲专区| 91av网一区二区| 免费人成视频x8x8入口观看| 精品人妻一区二区三区麻豆 | 午夜福利高清视频| 嫩草影院新地址| 亚洲人成网站在线播| 三级男女做爰猛烈吃奶摸视频| 他把我摸到了高潮在线观看| 欧美色视频一区免费| 国产激情偷乱视频一区二区| 亚洲精品一区av在线观看| 男女之事视频高清在线观看| 国产高清三级在线| 精品人妻1区二区| 免费看美女性在线毛片视频| 国产综合懂色| 又粗又爽又猛毛片免费看| 舔av片在线| 精品日产1卡2卡| 亚洲av二区三区四区| 国产精品,欧美在线| 国产欧美日韩一区二区精品| 国产午夜精品论理片| 精华霜和精华液先用哪个| 亚洲五月天丁香| 国产高清三级在线| 亚洲专区中文字幕在线| 色播亚洲综合网| 国产大屁股一区二区在线视频| 不卡一级毛片| 国产高清有码在线观看视频| 天堂网av新在线| 日日啪夜夜撸| videossex国产| 淫妇啪啪啪对白视频| 又粗又爽又猛毛片免费看| 国产伦人伦偷精品视频| 亚洲男人的天堂狠狠| 中文亚洲av片在线观看爽| 麻豆国产97在线/欧美| xxxwww97欧美| 在线a可以看的网站| 五月玫瑰六月丁香| 欧美中文日本在线观看视频| 欧美日韩黄片免| 综合色av麻豆| 直男gayav资源| 国产精品一区二区性色av| 综合色av麻豆| 中文在线观看免费www的网站| 久久久久久九九精品二区国产| av在线天堂中文字幕| 色av中文字幕| 日本免费一区二区三区高清不卡| 亚洲黑人精品在线| 黄色欧美视频在线观看| 色噜噜av男人的天堂激情| 禁无遮挡网站| netflix在线观看网站| avwww免费| 午夜福利视频1000在线观看| 一个人看视频在线观看www免费| 我的女老师完整版在线观看| 国产麻豆成人av免费视频| 日韩国内少妇激情av| 91麻豆精品激情在线观看国产| 麻豆成人午夜福利视频| 日韩,欧美,国产一区二区三区 | 国产成人一区二区在线| 国产91精品成人一区二区三区| bbb黄色大片| 国产精品乱码一区二三区的特点| 麻豆国产97在线/欧美| 国产v大片淫在线免费观看| 成人特级黄色片久久久久久久| 成年女人永久免费观看视频| 日本在线视频免费播放| 午夜福利18| 一区二区三区免费毛片| 国产av麻豆久久久久久久| 欧美bdsm另类| 美女 人体艺术 gogo| 色综合色国产| 国产精品女同一区二区软件 | 国产成人福利小说| 日本三级黄在线观看| 国产高潮美女av| 亚洲七黄色美女视频| 成人美女网站在线观看视频| 3wmmmm亚洲av在线观看| 三级毛片av免费| 亚洲国产精品合色在线| 亚洲性夜色夜夜综合| 尤物成人国产欧美一区二区三区| 亚洲精品一区av在线观看| 俺也久久电影网| 亚洲av成人精品一区久久| 亚洲美女黄片视频| 18禁裸乳无遮挡免费网站照片| 中文字幕高清在线视频| 国产精品,欧美在线| 国产亚洲欧美98| 亚洲人成网站在线播放欧美日韩| 99九九线精品视频在线观看视频| 日本三级黄在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲精品影视一区二区三区av| 男女那种视频在线观看| 在线天堂最新版资源| av天堂中文字幕网| 久久草成人影院| 九九爱精品视频在线观看| 免费av观看视频| 一区二区三区激情视频| 久久人人精品亚洲av| 99热只有精品国产| 有码 亚洲区| 制服丝袜大香蕉在线| 午夜激情欧美在线| 99热精品在线国产| 亚洲,欧美,日韩| 不卡一级毛片| 香蕉av资源在线| 欧美成人性av电影在线观看| 亚洲成人久久爱视频| 国产综合懂色| 国产精品一区二区三区四区免费观看 | 日韩欧美国产一区二区入口| 亚洲av免费在线观看| 欧美成人性av电影在线观看| 特大巨黑吊av在线直播| 国产熟女欧美一区二区| 亚洲av成人av| 亚洲成a人片在线一区二区| 91狼人影院| 欧美bdsm另类| 午夜免费激情av| 又黄又爽又刺激的免费视频.| 色哟哟哟哟哟哟| 床上黄色一级片| 亚洲欧美激情综合另类| 日本熟妇午夜| 国产高清激情床上av| 亚洲精品久久国产高清桃花| 欧美一级a爱片免费观看看| 久久久久久九九精品二区国产| or卡值多少钱| 成人二区视频| 亚洲无线观看免费| 在线观看av片永久免费下载| 亚洲18禁久久av| 亚洲五月天丁香| 久久草成人影院| 老司机深夜福利视频在线观看| 日韩欧美精品v在线| 日韩人妻高清精品专区| 亚洲av日韩精品久久久久久密| 麻豆成人午夜福利视频| 特级一级黄色大片| 五月伊人婷婷丁香| 男女视频在线观看网站免费| 亚洲第一电影网av| 999久久久精品免费观看国产| 日本精品一区二区三区蜜桃| 日本一二三区视频观看| 久久人妻av系列| 成人国产综合亚洲| 久久精品夜夜夜夜夜久久蜜豆| 天美传媒精品一区二区| 人妻久久中文字幕网| 特级一级黄色大片| 精品人妻视频免费看| 欧美成人性av电影在线观看| 国产视频内射| 国产精品不卡视频一区二区| 久久精品人妻少妇| 亚洲内射少妇av| 在现免费观看毛片| 精品久久久久久久久久免费视频| 国产精品一区二区性色av| 亚洲熟妇熟女久久| 国产精品伦人一区二区| 天天一区二区日本电影三级| 国产精品福利在线免费观看| 18禁在线播放成人免费| 波多野结衣巨乳人妻| 日韩欧美三级三区| 别揉我奶头 嗯啊视频| 美女免费视频网站| 亚洲成a人片在线一区二区| 免费在线观看影片大全网站| 国产亚洲精品久久久com| 国产精品1区2区在线观看.| 超碰av人人做人人爽久久| 国产激情偷乱视频一区二区| 变态另类丝袜制服| 国产探花在线观看一区二区| 久久久久久久久大av| 亚洲avbb在线观看| 亚洲精品一区av在线观看| 国产真实乱freesex| 日本成人三级电影网站| 一区二区三区免费毛片| 在线a可以看的网站| 成人午夜高清在线视频| 欧美精品啪啪一区二区三区| 久久6这里有精品| 搡老熟女国产l中国老女人| 精品久久国产蜜桃| 午夜视频国产福利| 欧美性猛交黑人性爽| 国产 一区 欧美 日韩| 国产一区二区在线观看日韩| 中文字幕免费在线视频6| 婷婷精品国产亚洲av在线| av女优亚洲男人天堂| 精品不卡国产一区二区三区| 午夜老司机福利剧场| 一个人看的www免费观看视频| 91在线精品国自产拍蜜月| 中国美白少妇内射xxxbb| 国产av不卡久久| 男女那种视频在线观看| 波多野结衣高清作品| 国产久久久一区二区三区| 一级a爱片免费观看的视频| 在线看三级毛片| 中文字幕av成人在线电影| 99精品久久久久人妻精品| 国产午夜福利久久久久久| 性色avwww在线观看| 亚洲经典国产精华液单| 男女视频在线观看网站免费| 精品不卡国产一区二区三区| 精品久久久久久久末码| 亚洲av二区三区四区| 国产精品久久久久久av不卡| 黄色日韩在线| 狠狠狠狠99中文字幕| 亚洲,欧美,日韩| 欧美性猛交╳xxx乱大交人| 看片在线看免费视频| 变态另类丝袜制服| 看片在线看免费视频| 国产精品三级大全| av专区在线播放| 久久精品影院6| 黄色欧美视频在线观看| 少妇的逼好多水| 国产高清三级在线| 变态另类丝袜制服| 窝窝影院91人妻| 国产国拍精品亚洲av在线观看| 成人美女网站在线观看视频| 搡老岳熟女国产| 久久久国产成人免费| 麻豆av噜噜一区二区三区| 狂野欧美白嫩少妇大欣赏| 国产精品不卡视频一区二区| 午夜a级毛片| 亚洲精品久久国产高清桃花| 淫妇啪啪啪对白视频| 日本五十路高清| 亚洲五月天丁香| 99久久成人亚洲精品观看| 九九在线视频观看精品| 久久久久久伊人网av| 51国产日韩欧美| 可以在线观看毛片的网站| 少妇高潮的动态图| 国内揄拍国产精品人妻在线| 日本 欧美在线| 亚洲在线观看片| 一区二区三区高清视频在线| 亚洲av免费在线观看| 人妻少妇偷人精品九色| 国产精品野战在线观看| 搡老岳熟女国产| 国产精品嫩草影院av在线观看 | 欧美中文日本在线观看视频| 麻豆久久精品国产亚洲av| 中国美女看黄片| 夜夜看夜夜爽夜夜摸| 国产亚洲av嫩草精品影院| 成年女人看的毛片在线观看| 可以在线观看的亚洲视频| 国产精品综合久久久久久久免费| 精品无人区乱码1区二区| 午夜影院日韩av| 欧美3d第一页| 少妇丰满av| 日本在线视频免费播放| 国产精品人妻久久久影院| 久久久精品大字幕| 精品久久久久久久久av| 色5月婷婷丁香| 一个人免费在线观看电影| 少妇熟女aⅴ在线视频| 99久久无色码亚洲精品果冻| 波野结衣二区三区在线| 毛片女人毛片| 欧美又色又爽又黄视频| 欧美3d第一页| 日本 av在线| 国产 一区 欧美 日韩| 欧美黑人欧美精品刺激| 成人一区二区视频在线观看| 狂野欧美激情性xxxx在线观看| 亚洲午夜理论影院| 中文字幕熟女人妻在线| 国产成人a区在线观看| 久久久久久久久大av| 免费看av在线观看网站| 变态另类成人亚洲欧美熟女| 1000部很黄的大片| 有码 亚洲区| 一级黄片播放器| 亚洲 国产 在线| 1024手机看黄色片| 久久香蕉精品热| www.www免费av| 国语自产精品视频在线第100页| 免费不卡的大黄色大毛片视频在线观看 | 亚洲乱码一区二区免费版| 亚洲成人久久爱视频| 老熟妇仑乱视频hdxx| 久久国内精品自在自线图片| 亚洲人成网站在线播放欧美日韩| 亚洲精品影视一区二区三区av| 日韩国内少妇激情av| 久久6这里有精品| 日韩精品中文字幕看吧| 一a级毛片在线观看| 亚洲精品456在线播放app | 日本欧美国产在线视频| 99视频精品全部免费 在线| 欧美成人性av电影在线观看| 91在线观看av| 国产成人福利小说| 欧美另类亚洲清纯唯美| 国产精品久久久久久久电影| 国产高清有码在线观看视频| 国产探花在线观看一区二区| 亚洲自拍偷在线| 午夜激情福利司机影院| 国产精品综合久久久久久久免费| 中亚洲国语对白在线视频| 中国美白少妇内射xxxbb| 五月玫瑰六月丁香| 一卡2卡三卡四卡精品乱码亚洲| 好男人在线观看高清免费视频| 嫁个100分男人电影在线观看| 又紧又爽又黄一区二区| 欧美日韩中文字幕国产精品一区二区三区| 97超级碰碰碰精品色视频在线观看| 搞女人的毛片| 制服丝袜大香蕉在线| 欧美绝顶高潮抽搐喷水| 在线观看美女被高潮喷水网站| 日日夜夜操网爽| 日本精品一区二区三区蜜桃| 日韩一区二区视频免费看| 成年女人永久免费观看视频| 人人妻人人看人人澡| 女生性感内裤真人,穿戴方法视频| 国产精品人妻久久久久久| 欧美精品国产亚洲| 国产高清三级在线| 久久精品国产亚洲av香蕉五月| 亚洲avbb在线观看| www日本黄色视频网| 亚洲avbb在线观看| 国产精品av视频在线免费观看| 免费黄网站久久成人精品| 深夜a级毛片| 久久久国产成人精品二区| 2021天堂中文幕一二区在线观| 日本撒尿小便嘘嘘汇集6| 亚洲av中文av极速乱 | 黄色一级大片看看| 免费大片18禁| 人人妻人人看人人澡| 国产成人福利小说| 亚洲人成网站高清观看| 美女黄网站色视频| 日本一本二区三区精品| 十八禁网站免费在线| 国产在线精品亚洲第一网站| 国产免费av片在线观看野外av| 午夜福利成人在线免费观看| 自拍偷自拍亚洲精品老妇| 中文字幕精品亚洲无线码一区| 国产高潮美女av| 午夜爱爱视频在线播放| 国产乱人视频| 少妇人妻精品综合一区二区 | x7x7x7水蜜桃| 免费在线观看影片大全网站| 麻豆成人午夜福利视频| 精品人妻偷拍中文字幕| 精品无人区乱码1区二区| 日韩精品青青久久久久久| 韩国av在线不卡| 午夜影院日韩av| 搡女人真爽免费视频火全软件 | 亚洲av第一区精品v没综合| 特级一级黄色大片| 岛国在线免费视频观看| 国产午夜精品论理片| 我要看日韩黄色一级片| 真人做人爱边吃奶动态| 露出奶头的视频| 日日摸夜夜添夜夜添av毛片 | 精品福利观看| 国产亚洲精品久久久com| 欧美日韩黄片免| 国产欧美日韩一区二区精品| 不卡一级毛片| 1024手机看黄色片| 欧美+亚洲+日韩+国产| 18禁黄网站禁片午夜丰满| 我的老师免费观看完整版| 国产亚洲av嫩草精品影院| 国产黄片美女视频| 美女xxoo啪啪120秒动态图| 午夜视频国产福利| 综合色av麻豆| 窝窝影院91人妻| 亚洲av美国av| 天堂影院成人在线观看| 国产精品免费一区二区三区在线| 男人的好看免费观看在线视频| 少妇熟女aⅴ在线视频| 欧美在线一区亚洲| 国产精品1区2区在线观看.| 日本一二三区视频观看| av.在线天堂| 亚洲人与动物交配视频| 国产精品久久久久久亚洲av鲁大| 中出人妻视频一区二区| 香蕉av资源在线| 国产av麻豆久久久久久久| 免费看美女性在线毛片视频| 成年版毛片免费区| 亚洲欧美日韩高清专用| 波野结衣二区三区在线| 日本欧美国产在线视频| 变态另类成人亚洲欧美熟女| 亚洲最大成人手机在线| ponron亚洲| 国产日本99.免费观看| 国产精品综合久久久久久久免费| 97超级碰碰碰精品色视频在线观看| 黄片wwwwww| 亚洲第一电影网av| 亚洲18禁久久av| 久久久久久久久中文| 九九久久精品国产亚洲av麻豆| 成人美女网站在线观看视频| 国国产精品蜜臀av免费| 18+在线观看网站| 综合色av麻豆| 免费无遮挡裸体视频| 国产成人影院久久av| 国产极品精品免费视频能看的| 18禁黄网站禁片午夜丰满| 日本与韩国留学比较| 国产成人av教育| 黄片wwwwww| 老女人水多毛片| 亚洲最大成人av| 久久久久久久久久成人| 毛片一级片免费看久久久久 | 亚洲av.av天堂| 啦啦啦观看免费观看视频高清| 我要看日韩黄色一级片| 成人国产综合亚洲| 最后的刺客免费高清国语| av在线蜜桃| 欧美色欧美亚洲另类二区| eeuss影院久久| 真实男女啪啪啪动态图| 丰满乱子伦码专区| 人妻制服诱惑在线中文字幕| 一个人观看的视频www高清免费观看| 亚洲国产高清在线一区二区三| 精品人妻熟女av久视频| 亚洲欧美日韩无卡精品| 亚洲美女搞黄在线观看 | 俄罗斯特黄特色一大片| 成人永久免费在线观看视频| 男女之事视频高清在线观看| 国产亚洲精品久久久久久毛片| 日韩亚洲欧美综合| 免费av不卡在线播放| 一进一出抽搐动态|