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

    Summery Intra-Tidal Variations of Suspended Sediment Transportation–Topographical Response and Dynamical Mechanism in the Aoshan Bay and Surrounding Area, Shandong Peninsula

    2021-12-22 11:20:28WANGYanDOUYanguangZOULiangGAOFeiSUDapengHURuiYUEBaojingandXUEBiying
    Journal of Ocean University of China 2021年6期

    WANG Yan, DOU Yanguang, 2), *, ZOU Liang, GAO Fei, SU Dapeng, HU Rui, YUE Baojing, and XUE Biying

    Summery Intra-Tidal Variations of Suspended Sediment Transportation–Topographical Response and Dynamical Mechanism in the Aoshan Bay and Surrounding Area, Shandong Peninsula

    WANG Yan1), DOU Yanguang1), 2), *, ZOU Liang1), GAO Fei1), SU Dapeng1), HU Rui1), YUE Baojing1), and XUE Biying1)

    1),,266071,2),,266071,

    As the second largest bay in Qingdao, the Aoshan Bay and its adjacent sea area play an important role in aquaculture deve- lopment and urban land and sea coordination for the eastern gulf type of city in the Qingdao Blue Silicon Valley Core Area (QBSVCA). Based onsedimentary dynamical observation and previous research results, the thermohaline structure, the transportation of suspended sediment and its mechanism, and the coastal geomorphic response were elaborated and analyzed in detail in this paper. The result indicated that the thermohaline and density distribution have obvious intra-tidal characteristics in the QBSVCA and the adjacent waters of the islands, during summer neap tide stage. The development of the bottom high suspended sediment concentration (SSC) layer was slightly enhanced in flood slack at each of the four stations. Suspended sediment transportation near the QBSVCA is related closely with the vertical mixing-stratification mechanism. Combined with previous research results, this study once again showed that submarine topography and the grain size of sea bed sediments would respond to hydrodynamic forces. The medians of the bottomand D50in the Aoshan Bay were the highest, followed by those in the Daguan Island and Xiaoguan Island, and the data in the Laoshan Bay were the lowest. This showed that the capacity of suspended sediment transportation in the bottom water layer of the Aoshan Bay was stronger than that of the adjacent sea area. The re-suspension and migration of fine sediments lead to the strong coarsening of sediments in this area.

    suspended sediment transportation; thermohaline structure; dynamical mechanism; topographical response; coastal environment

    1 Introduction

    The research on the bay and coastal zone is closely related to the urban economic development, so it is necessaryto closely follow the urban planning and the national strategy.The Qingdao Blue Silicon Valley Core Area (QBSVCA) is not only the new business card of Jimo District, but also a national key project in the whole Shandong Peninsula development. The Aoshan Bay is the second largest bay, onlynext to the Jiaozhou Bay, in Qingdao. The eastern gulf type of city relying on the Aoshan Bay Group is an important fulcrum in the ‘Three Bays’ of the gulf-based city pattern of Qingdao.

    The Aoshan Bay is a semi-closed bay located in the nor- theast of Aoshanwei sub-district at the eastern part of the QBSVCA. The mouth of the bay faces to the Southeast. There are 9 reefs in the bay, such as the Ganzui Island, the Zhanggong Island and the Beijiao Reef. Along the bank of the Aoshan Bay, there are several small seasonal rivers with small runoff, including the Daren River, the Wenquan River, the Xinsheng River, the Gaoyu River, the Daqiao River, the Wangcun River, and so on. Similarly, the Xiaodao Bay is also a semi-closed bay located in the southeast of Aoshanwei sub-district, which is a part of the Laoshan Bay. There’s no river flow into the Laoshan Bay, within which the Tuzi Island, the Small Island, the Lion Island, the Daguan Island and the Maer Island are scattered (Fig.1). The Aoshan Bay and the Laoshan Bay are dominated by wind and waves with small surges. Both strong and normal wind directions are southeast in these bays, with a frequency of 23% and a wave height of 2.2m. The tidal type in these bays is normal semidiurnal tide, and the mean sea level is low in winter while high in summer. The se- diment sources mainly include river discharging sediment, coastal erosion sediment and offshore transported sediment,., while the river discharging sediment and the erosion sediment near the Aoshan Bay are limited. Suspended sedi- ment in the offshore sea can spread to the Aoshan Bay through water exchange and then deposit on the seabed (Zhao and Chen, 2001; Chen and Zhao, 2004; Li., 2011). Marine environment quality in the inshore water of Jimo District is related to marine fishery, marine nature reserves and reserves of rare and endangered marine life. The ecosystem provincial marine special protection zones of the Daguan Island and the Xiaoguan Island group in Jimo District was established on July 10, 2015, with the approval of Shandong Provincial Government. The biocoenosis of the protection zone was stable and the habitat environment was good (Qingdao Municipal Ocean and Fisheries Administration, 2018). The adjacent waters of Xiaoguan Island and the Aoshan Bay are important sea areas for proliferating and releasing prawns in Qingdao. The Aoshan Bay is also an important culturing area for theraft culturing of scallops and the bottom seeding and conservation of variegated clam in Jimo District. The prolife- ration and release of the Japanese prawns will promote the development of marine fishing industry, aquatic pro- cessing industry, ship maintenance industry, catering industry and other industries, which produced high social benefits. Scallop breeding industry has also created huge economic benefits. The growing and the distribution of aqua-tic products, affected by marine thermohaline environmentand hydrodynamic conditions, also have seasonal distribution characteristics. The breeding industry of aquatic pro- ducts mainly implemented in summer. Thus, the investigation of thermohaline and hydrodynamic environment in summer could provide scientific basis for the development of aquaculture in the Aoshan Bay and the Laoshan Bay near the QBSVCA.

    This research focuses on the construction needs of the QBSVCA, in the offshore areas of the Aoshan Bay, the Xiaodao Bay and the Dingzi Bay, around the QBSVCA, carrying out the sedimentary dynamic investigation of coastal geological environment in coastal important engineering areas. Hydrodynamic data and hydrological environment data were collected systematically to reveal the characteristics of sedimentary dynamic environment in this area. Based on the investigation of hydrographic or thermohaline structures and hydrodynamic conditions in the QBSVCA and the adjacent waters of islands, the spatial and temporal distribution of sedimentary dynamic parameters and the thermohaline environment features in typical areas were determined during summer neap tides in Aoshan Bay and Laoshan Bay, which provides a scientific basis for the development of aquaculture related industries in this area.

    2 Materials and Methods

    2.1 Cruise and Sampling

    The observation of sedimentary dynamics was carried out in the survey area during neap tide from June 15 to 21, 2017. A total of 208 suspended sediment samples were collected after 25hcontinuous sedimentary dynamic observation at 4 stations (Fig.1). The ‘used to conduct the sedimentary dynamics survey. The Trimble SPS351 satellite/beacon differential GPS system and the Haida6.0 integrated navigation and positioning system were adopted to achieve navigation and positioning. Thermohaline environment observations wereaccomplished by the SBE19plus CTD of Seabird Company. Current observations adopted a direct reading sea current profiler of Acoustic Doppler Current Profiler (ADCP) of RDI Company (broadband=300kHz and blanking zone=1m), with measurement errors for the velocity and direction of <5mms?2and <1?, respectively. Suspended sediment samples were collected by the ECO-55 automatic water sampler of Seabird Company of the United States.

    Fig.1 Study area, observation stations, and previous beach profiles in the offshore areas of the Aoshan Bay, the Xiaodao Bay and the Dingzi Bay, around the QBSVCA. The re- mote sensing base map was obtained from Google Earth.

    CTD was dropped once per hour to obtain temperature, salinity and turbidity data, which can reflect the spatial and temporal variation of water thermohaline and turbid conditions within 25h at each of the 4 stations. The current data were stably and continuously measured within 25h, which can reflect the hydrodynamic characteristics of all the stations. The interval of measured layer is set as 1m, and the sampling frequency is set as 10s. Suspended sediment samples were collected once per hour with the dropping of the CTD instruments, obtained from the surface layer (0.5m–2m below the sea surface) and the bottom layer (0.5m–2m above the seabed).

    2.2 Calibration of Suspended Sediment Concentration (SSC) and Turbidity

    Suspended sediment concentration (SSC) was measured by filtering suspended sediment samples on double-layer filter membranes weighed previously (the weight of the upper membrane and the lower membrane in this link were represented by WandW, respectively), and then drying and weighing the filter membranes (the weight of the upper membrane and the lower membrane in this link were represented by WandW, respectively) (Eq. (1)). In this experiment, the filter process adopted cellulose acetate filter membranes with aperture of 0.45μm, while the weighing process was operated on the electronic balance with precision of 0.00001g. The filter membranes were dried for 24h before and after the filter process with drying temperature about 45℃. The calculation formula ofis as follow,

    whererespects the volume of filtered water.

    The turbidity values of seawater collected from the layer where suspended sediment samples were filtered were averaged and fitted to SSCs. As shown in Fig.2, four abnormal data accounting for 0.02% of total data were remov- ed and the coefficient of determination (2) was 0.902116, suggesting the validity of calibration to theand turbidity.

    Fig.2 Correlation between suspended sediment concentration (SSC) and turbidity.

    3 Results

    3.1 Temporal-Spatial Distribution of Thermohaline Parameters in the Tidal Cycle

    1) The central part of the Aoshan Bay

    The thermohaline parameters and the density of water column were vertically homogeneous in flood tide phase at Station CM01 (average depth of 7.2m), which is located at the central part of the Aoshan Bay (Figs.1 and 3a), with salinity<32.1, temperature<21℃,>22.2kgm?3.here represented the density of water column and its calculation formula is as follows,

    whererepresents the density of water column. Then, the temperature of the water column was gradually increased from the top layer in flood slack, while the salinity was simultaneously decreased from the top layer or the middle layer, with downward increasing vertical density gradient. In ebb tide, the temperature and the salinity of water co- lumn increased over time, with vertically homogeneous ther- mohaline structure and decreasing vertical density. Until ebb slack, the temperature and the salinity of the water column reached the highest and the density of water co- lumn dropped to the lowest in the tidal cycle. Overall, the salinity of the water column at Station CM01 was in the range of 31.95–32.45, while the temperature range is 20.4–23.5℃, with the characteristics of higher temperature and salinity and lower density in the upper-middle-layer water during ebb slack phase than other phases of the tidal cycle.

    2) North of the Laoshan Bay

    Compared with that at Station CM01, the vertical gradients of salinity and density in the water column were higher at Station CM02 (average depth of 13.2m), located at the northern part of the Laoshan Bay (Figs.1 and 3b), with the salinity range of 31.5–32.2 and the temperature range of 17.3–18.7℃. Another difference between these two stations is the intra-tidal vertical salinity variation characteristic of the water column, which showed increase trend in ebb tide phase and decrease trend in flood one at CM02 Station, while appeared to be a reverse trend at the other station. In ebb tide, the vertical temperature gradient of the water column is higher than that in flood tide, and the thermohaline environment of high temperature, low salt and low density was vertically gradually developed. In ebb slack, the temperature of the water column tends to be more homogeneous. In flood tide, the vertical temperature gradient of the water column increased, and the water column with high salinity, low temperature and high density developed until flood slack.

    3) East of the Maer Island

    Similar to that at Station CM02, the variation of thermohaline environment at Station CM03 (average depth of 17.3m) was characterized by high vertical gradient of salinity and density in the water column (Figs.1 and 3c), with the salinity range of 31.2–32.4 and the temperature range of 16.3–18.5℃. At the depth of less than 7m above the bottom of seabed (MAB), the vertical thermohaline and density structures of the water column were homogeneous. In flood tide, along with the homogeneously increase of salinity, the vertical temperature structure at the depth of more than 7MAB was initially homogeneous with lower temperature, and then rapidly occupied by enhancing stratification structure with higher temperature, accompanied by enhancing stratification of vertical sali- nity structure. Until the next flood tide, the vertical temperature structure tended to be homogeneous again.

    4) East channel of the Xiaodao Bay

    Similar to that at Station CM01, along with the homogeneous decrease of salinity, the temperature of the water column at Station CM04 (average depth of 8.5m) was uni-formly distributed with lower temperature (Figs.1 and 3d), with the salinity range of 31.7–32.2 and the temperature range of 17.8–21.3℃. The variation regularity was different from two other stations with the depth of more than 10m.

    Fig.3 Temporal-spatial distribution of thermohaline parameters in the tidal cycle. MAB, meters above bottom.

    In summary, the thermohaline and density distribution have obvious intra-tidal characteristics in the QBSVCA and the adjacent waters of islands, during summer neap tide stage. Water salinity here was basically around 32. Temperature decreased from the coast of Aoshanwei sub-district (about 22℃) to the Laoshan Bay (18℃ or so). Temperature in summer is higher in shallower waters, which is unfavorable to biological cultivation (Janet and Dudgeon, 1996; Naylor., 1998; Joos and Spahni, 2008). It can provide scientific basis for the development of aquaculture industry to make clear the characteristics of the thermohaline environment in the QBSVCA and the adjacent waters of islands.

    3.2 The Variation Characteristic of Current Velocity

    The horizontal current velocity at Station CM01 was 200–800mms?1, and the current velocities of surface water and bottom water were higher than that of middle water. The average current direction of water column during floodtide was NNW, while it was SSE in ebb tide, and the magnitude of horizontal velocity can be more than 410mms?1in surface water and bottom water during flood tide with the direction of SE (Fig.4a). The eastern component of current velocity () tended to the direction of west (the value ofwas negative) in flood slack and ebb slack, while to the east (the value ofwas positive) in other periods. The northern component of current velocity () tended to the north (the value ofwas positive) in flood tide, while to the south (the value of v wasnegative) during ebb tide (Fig.5a).

    In the northern area of the Laoshan Bay (Station CM02), the horizontal current velocity was 300–600mms?1, and the vertical distribution of current velocity was more homogeneous than that in the Aoshan Bay (Station CM01), with uniform structures ofand. The average current direction of water column during flood tide was SWW-NWW while it was NE in ebb tide, and the magnitude of horizontal velocity may be more than 500mms?1during ebb tide (Fig.4b). The velocitytended to the east in flood tide, while to the west in ebb tide. The velocitytended to the north in flood tide, while to the south during ebb tide (Fig.5b).

    Fig.4 The variation characteristic of current velocity.

    Fig.5 The variation characteristic of east, west, and vertical components of current velocity.

    In the eastern area of the Maer Island (Station CM03), the horizontal current velocity was 300–750mms?1, and the vertical distribution of current velocity was homogeneous with uniform structure ofand, similar to Station CM02. The average current direction of water column wasSSW-W during flood tide, while it was NE in ebb tide, and the horizontal velocity more than 500mms?1mainly appeared during ebb tide (Fig.4c). The componenttended to the east in flood tide, while to the west in ebb tide. The componenttended to the north in flood tide, while to the south during ebb tide (Fig.5c).

    On account of the restriction by topography, the average current direction of water column in the east channel of the Xiaodao Bay (Station CM04) was NNW-NNE in most of the time. And the range of horizontal current velocity was 300–600mms?1, with higher values in surface layer and bottom layer than those in middle layer (Fig.4d). The componenttended to the west in flood tide, while to the east in ebb tide. The componenttended to the north throughout tidal cycle due to the restriction of topography (Fig.5d).

    Overall, the regular semidiurnal tide was dominant in the coastal waters of the QBSVCA, the reversing current was dominant in the Aoshan Bay, and the rotary current was dominant outside the bay (Fig.4; Li, 2012). The peak values of the vertical average current velocity in the Ao- shan Bay (Station CM01) and the east channel of the Xiao- dao Bay (Station CM04) were in the range of 225–278mms?1, while those outside the bay (Station CM02 and Station CM03) were 174–621mms?1(Fig.6). In the Ao- shan Bay, the vertical average current velocity showed quickly increasing first then slowly decreasing trend in flood tide, while it showed slowly increasing first then quickly decreasing trend in ebb tide, with the peak values appearing in the middle stage of flood tide and ebb tide, respectively (Figs.4a and 6). Outside the Aoshan Bay, the gradients of the increasing stage of vertical average current velocity were similar to those of the decreasing stage, with the peak values of vertical average current velocity appearing in flood slack and ebb slack (Figs.4b, c, d and 6).

    3.3 The Intra-Tidal Distribution of Suspended Sediments

    The surface sediments at each station are mainly clay silt (Li., 2011; Li, 2012). Their resuspension and re- transportation affect the water transmittance and the distribution of benthic organisms in this sea area (Steyaert., 2003; Fu., 2012). Sediment sources in the coastal area of the QBSVCA mainly include river-carrying sediments, coastal erosion and offshore sediments,. The occurrence of near-shore coarse sediments is mainly caused by the settlement of fluvial sediments and coastal erosion, while fine-grain sediments in deep-water area was mainly derived from suspended sediment deposition (Li., 2011).

    Fig.6 The average current velocity of water column.

    Based on the fitting relationship betweenand turbidity (Fig.2), the spatial and temporal distribution of SSCs at each station was obtained, as shown in Fig.7. Basically, the variation of SSCs at the four stations was relevant to the variation of horizontal current velocities and thermohaline structures (Figs.3 and 7). The larger the velocity and density were, the greater thewould be. Therefore, the spatial and temporal distribution of SSCs also had distinct intra-tidal characteristics in and outside the Aoshan Bay. In the Aoshan Bay (Station CM01), the bottomrapidly increased in ebb slack and the middle stage of flood tide with the peak value more than 45mgL?1, while promptly decreased in the earlier stage of flood tide and the stage transiting from flood slack to next ebb slack, with peak value less than 25mgL?1(Fig.7a). In the east channel of the Xiaodao Bay (Station CM04), the trans- portation of bottom suspended sediments was inhibited in ebb tide and the middle stage of flood tide, while it was enhanced in ebb slack and flood slack (Fig.7b). At Station CM02 and Station CM03, located at the northeast and the southeast of the Laoshan Bay, the suspended sediment transportation in the bottom layer of the water column was enhanced in flood slack and ebb slack, while it was inhibited in other stages (Figs.7c and d).

    In general, as shown in Fig.7, SSCs in the water column were higher at Station CM01 and Station CM02 than two other stations. The bottom high SSC layers were developed or slightly enhanced in flood slack at each of the four stations. In addition, suspended sediment plumes occur- red at all of the four stations within two meters of surface water layer, which might be formed by river-carrying se- diment transportation or offshore sediment transportation, with the similar order of SSCs to the bottom water layer.

    4 Discussion

    4.1 The Spatial and Temporal Variation Characteristics of Suspended Sediment Flux (SSF)

    SSF were calculated by decomposing instantaneous material transport vector based on relative water depths (Dyer, 1974; Su and Wang, 1986; Wu., 2006). According to the calculation formula proposed by Wu. (2006), the calculation formula of instantaneous SSF () per unit meter width is as follows,

    Fig.7 The intra-tidal distribution of suspended sediments.

    The calculation results show that the direction ofwas mostly towards the north in flood tides at Station CM01 in the Aoshan Bay, while it was to the direction ranging from the south to the southeast in ebb tides (Fig.8a). The values ofcould reach up to 62.6g(ms)?1 in flood tides, while basically lower than 18.1g(ms)?1 in ebb tides. Thesediment transportation was restrained during flood slacks and ebb slacks, accompanied by the diversion of direction of suspended sediment transportation.

    At Station CM02 in the northern area of the Laoshan Bay, the transport direction of suspended sediment was counterclockwise (Fig.8b). The duration of suspended se- diment transporting northward was longer than that southward, and the values ofwere larger during the northward with the maximum of 145g(ms)?1than during the southward.

    In the eastern area of the Maer Island (Station CM03), the temporal variation of the direction and value ofwas similar to that in the northern area of the Laoshan Bay (Station CM02) (Figs.8b and 8c). At Station CM03, the values ofexceeded the maximum ofat Station CM02 (145g(ms)?1) for 3h, with the maximum of 252g(ms)?1.

    Combined with the temporal and spatial variation ofs and current velocities (Figs.4d and 7d), suspended sediment transport at Station CM04 in the east channel of the Xiaodao Bay was mainly towards NNW direction during flood tide and flood slack, with the maximum of 150g(ms)?1(Fig.8d).

    4.2 The Vertical Mixing and Stratification Mechanism of Suspended Sediment Transportation

    To reflect the dynamical mechanism of suspended sediment transportation in the water column, Richardson number () is usually calculated, which is closely related to the density effect of vertical stratification and turbulence in the water column (Zhang and Anthes, 1982; Kundu, andBeardsley, 1991; Wang., 2010; Wang., 2014). The formula ofis as follows:

    where ρis the bulk density of water, u is the horizontal current velocity, z is the vertical dimension, positive upward from seabed. With parameters of water salinity, temperature and pressure, ρ can be calculated from the UN- ESCO equation (Millero et al., 1980). The numeratorpart of the function represents the stratification stability of the water column, determined by density gradient. The denominator part indicates the turbulent intension driven by current shear stress. The high Ri value shows more stable status of the water column. On the contrary, the low value of Ri indicates turbulence intensification of the water column. Generally, Ri<0.25 is a necessary condition for turbulent mixing to overcome stratification structure. It also means when Ri>0.25, stratification structure of the water column is stable, and not easy to be broken by turbulent mixing (Turner, 1973).

    Fig.9 The vertical mixing and stratification mechanism of suspended sediment transportation.

    The hourly vertical-averaged values ofandat the four stations were shown in Fig.9. At Station CM01 and Station CM04, the increases ofvalues were strictly in accord with the situation of<0.25, while the decreases ofhappened to coincide with the situation of>0.25 (Figs.9a, d). In the inshore area of the QBSVCA (Station CM01 and Station CM04), dramatic increase inmostly occurred whenwas less than 0.25 and decreasing, while significant reduction inhappened whenwas more than 0.25 and increasing. The above phenomenon suggested that suspended sediment transportation near the QBSVCA were related closely with vertical mixing and stratification mechanism. When turbulent mixing of the water column overcame stratification structure, it could enhance the resuspension of seabed sediments. However, if stratification structure of the water column was stable, the resuspension of seabed sediments was inhibited. Mean- while, it also gave significant evidence that near-shore su- spended sediments of the QBSVCA were mainly derived fromresuspension of seabed sediments (Figs.9a, d). The poor correlation betweenandat Station CM02 and Station CM03 (Figs.9b, c) represented the suspended sediments beyond the Laoshan Bay originated from more complicated sources, which could be explained by future observational studies.

    4.3 The Topographical Response of Suspended Sediment Transportation

    For the project ‘Comprehensive survey and potential eva- luation of mineral resources in Shandong Peninsula coastal zone’, which this paper relies on, surface sediments were sampled in the Aoshan Bay, the Laoshan Bay, and the ocean area surrounding the Daguan Island and the Xiao- guan Island in May 2017. These surface sediment samples were analyzed by particle size compositions (Qingdao Institute of Marine Geology, 2018). In summer and winter of the year 2017, the project also monitored two beach topographic profiles along the west side of the Ao- shan Bay (beach profiles C1 and C2 in Fig.1).

    Combined with the previous observations of beach pro- files (QAM4, QAM5 and QAM6 in Fig.1, Bi., 2015), there are obvious seasonal differences in erosion and de- position in the western beach of the Aoshan Bay. The re- sults showed that the western beach of the Aoshan Bay was silted up in winter and eroded away in summer. In terms of particle size of sediments on the west coast of the Ao- shan Bay, it was generally finer in winter than that in sum- mer. The near-shore sediment transportation plays a significant role in coastal erosion and deposition. Thus, the hydrodynamic erosion on the western beach was stronger in summer than that in winter (Bi., 2015; Qingdao Institute of Marine Geology, 2018).

    Fig.10 The direction of sediment transportation obtained in this study and previous researches (Li, 2012).

    In aspect of topographic and geomorphic responses, the role of bottom suspended sediment transportation is more prominent. For the convenience of discussion, the study area was divided into three regions (Fig.10), namely, the Aoshan Bay, the Daguan Island and Xiaoguan Island, and the Laoshan Bay. The Station CM01 was located in the Aoshan Bay, the Station CM04 was located near the Daguan Island and Xiaoguan Island, the Station CM02 and Station CM03 were located in the Laoshan Bay. Compared with the median grain size (D50) of surface sediments in these three regions (Fig.11), the variation trend of the median of the bottomwas consistent with that of D50. Among the three regions, the medians of the bottomand D50in the Aoshan Bay were the highest, followed by those near the Daguan Island and Xiaoguan Island, and the data in the Laoshan Bay were the lowest. This result showed that the capacity of suspended sediment transportation in the bottom water layer of the Aoshan Bay was stronger than those of two other regions. The re-suspension and transportation of fine particles resulted in the coarsening of sediments in this region. The above conclusion is based on the beach profilesin previous observations and the correlation of the bottomand D50, which can be used to understand the relationship between the suspended sediment transportation and the particle size compositions of the bottom sediments. However, the grain size characteristic is controlled by many factors, and is the result of long-term hydrodynamic modulation. This result was obtained in a short duration. Thus, in order to verify the conclusion in this section, more observations should be carried out in the future.

    Fig.11 The bottom E and the median grain size (D50) of surface sediments in the three regions.

    Further, the direction of sediment transportation (Fig.10) obtained in this study and one previous research (Li, 2012) provided sedimentary dynamical evidences for the above results. The suspended sediments in the middle and the west of the Aoshan Bay were mostly derived fromresuspended sediments, and the directions of transportation were towards the bay both in neap tide and in spring tide. The suspended sediments in the east of the Aoshan Bay were delivered by coastal currents in spring tide, while they were poured out of the bay in neap tide. It can be seen that most sediments settled down in the Aoshan Bay were mainly delivered by coastal currents from the east of the bay in spring tide. The direction of sediment transportation in the Laoshan Bay in spring tide was almost opposite to that in neap tide.

    5 Conclusions

    Thesedimentary dynamical observations were conducted at 4 stations in the Aoshan Bay and the adjacent sea area. They showed obvious intra-tidal characte- ristics in thermohaline structure, current velocity and. The regular semidiurnal tide was dominant in the coastal waters of the QBSVCA, the reversing current was dominant in the Aoshan Bay, and the rotary current was dominant outside the bay. In the Aoshan Bay, the peak values of vertical average current velocity appeared in the middle stages of flood tide and ebb tide. SSCs in the water column were higher at two stations in northern part of the study area. The bottom highlayer was developed or enhanced in flood slack at the four stations. The variation ofshowed obvious intra-tidal characteristics, which was affected by the vertical mixing and stratification mechanism of the water column. When turbulent mixing overcame stratification structure in the water column, it could enhance the resuspension of seabed sediments. However, if the stratification structure of the water column was stable, the resuspension of seabed sediments was inhibited. Also, the near-shore suspended sediment transportation of the QBSVCA was derived from theresuspension of seabed sediments. The capacity of suspended sediment transportation in the bottom water layer of the Aoshan bay was stronger than those of the adjacent sea area. The resuspension and transportation of fine sediments resulted in the coarsening of sediments in this region. These results can provide scientific basis for the development of aquaculture industry and urban land and sea coordination in the QBSVCA.

    Acknowledgements

    This work was jointly supported by the National Natural Science Foundation of China (No. 41606082) and the China Geological Survey (Nos. DD20189230, DD20160 148).

    Bi, S., Zhang, Y., and Bie, J., 2015. Monitoring and analysis of beach topographic profiles in the southern Jiaodong Peninsula., 39 (3): 71-76 (in Chinese with English abstract).

    Chen, J., and Zhao, J., 2004. Distribution and variation of hydrographic factors in the Aoshan Bay., 25 (2): 66-72 (in Chinese with English abstract).

    Dyer, K., 1974. The salt balance in stratified estuaries., 2 (3): 273-281.

    Fu, S., Mu, F., Yang, S., and Wang, C., 2012. Study on the spatial-temporal distribution of Meiofauna in the intertidal zone of Cangkou, Qingdao.,42 (Sup.): 124-130 (in Chinese with English abstract).

    Janet, E., and Dudgeon, S., 1996. Temperature dependent change in the complexity of form offronds., 207 (1-2): 15-24.

    Joos, F., and Spahni, R., 2008. Rates of change in natural and anthropogenic radiative forcing over the past 20,000 years., 105 (5): 1425- 1430.

    Kundu, P., and Beardsley, R. C., 1991. Evidence of a critical Richardson number in moored measurements during the upwelling season off northern California., 96: 4855-4868.

    Li, S., 2012. Research on hydrodynamic environment and sediment movement characteristics in Aoshan Bay. PhD thesis. Ocean University of China, Qingdao, 49-50.

    Li, S., Song, Y., Bian, S., Liu, J., and Hu, Z., 2011. Research on hydrodynamic environment and sediment movement characteristics in Aoshan Bay., 28 (5): 58-62 (in Chinese with English abstract).

    Millero, F., Chen, C., Bradshaw, A., and Schleicher, K., 1980. A new high pressure equation of state for seawater.,27 (3/4A): 255-264.

    Naylor, R., Goldburg, R., Mooney, H., Beveridge, M., Clay, J., Folke, C.,., 1998. Nature’s subsidies to shrimp and salmon farming., 282 (5390): 883-884.

    Qingdao Institute of Marine Geology, 2018. Results report of the project ‘Comprehensive survey and potential evaluation of mineral resources in Shandong Peninsula coastal zone’. Qingdao Institute of Marine Geology, Qingdao, 68-102 (in Chinese).

    Qingdao Municipal Ocean and Fisheries Administration, 2018. Qingdao Marine Environment Bulletin (2017). Qingdao Municipal Ocean and Fisheries Administration, Qingdao, 33-34 (in Chinese).

    Steyaert, M., Vanaverbeke, J., Vanreusel, A., Barranguet, C., Lucas, C., and Vincx, M., 2003. The importance of fine-scale, vertical profiles in characterising nematode community structure., 58 (2): 353-366.

    Su, J., and Wang, K., 1986. The suspended sediment balance in Changjiang Estuary., 23 (1): 81-98.

    Turner, J., 1973. Buoyant convection from isolated sources. In:. Cambridge University Press, Cam- bridge, 165-206.

    Wang, H., Bi, N., Wang, Y., Saito, Y., and Yang, Z., 2010. Tide- modulated hyperpycnal flows off the Huanghe (Yellow River) Mouth, China., 35: 1315-1329.

    Wang, Y., Li, R., Wen, Z., Chen, X., and Chen, B., 2014. The summer spring-neap variation of the water thermohaline-tur- bidity structure and its dynamical mechanism in the southern Bohai Strait., 90: 52-59.

    Wu, D., Zhang, R., Yan, Y., and Li, R., 2006. Mechanism of suspended sediment transport in Dongdagang tidal channel of radial sand ridges., 34 (2): 216-222 (in Chinese with English abstract).

    Zhang, D., and Anthes, R., 1982. A high-resolution model of the planetary boundary layer-sensitivity tests and comparisons with SESAME-70 data., 21: 1594- 1609.

    Zhao, J., and Chen, J., 2001. Research on hydro-environment of Aoshan Bay., 22 (3): 59-63 (in Chinese with English abstract).

    July 2, 2020;

    December 4, 2020;

    June 1, 2021

    ? Ocean University of China, Science Press and Springer-Verlag GmbH Germany 2021

    . Tel: 0086-532-85731050

    E-mail: douyanguang@gmail.com

    (Edited by Chen Wenwen)

    国产精品久久久久久人妻精品电影| 日日夜夜操网爽| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 蜜桃久久精品国产亚洲av| 小蜜桃在线观看免费完整版高清| 九九在线视频观看精品| 99久久精品热视频| 欧美在线一区亚洲| 国产v大片淫在线免费观看| 久久久久久久久中文| 不卡一级毛片| 久久精品国产自在天天线| 国产精品精品国产色婷婷| 在线观看日韩欧美| 观看免费一级毛片| 亚洲人成网站在线播放欧美日韩| 男女做爰动态图高潮gif福利片| 18禁美女被吸乳视频| 成年免费大片在线观看| 91九色精品人成在线观看| 悠悠久久av| 18禁黄网站禁片免费观看直播| 乱人视频在线观看| 一级黄片播放器| 成人一区二区视频在线观看| 日韩欧美 国产精品| 18禁国产床啪视频网站| 欧美激情在线99| 日本黄色片子视频| 9191精品国产免费久久| 午夜免费成人在线视频| 亚洲性夜色夜夜综合| 18禁黄网站禁片午夜丰满| 亚洲欧美日韩高清在线视频| 999久久久精品免费观看国产| 久久久久久久久大av| 欧美激情在线99| 欧美一级毛片孕妇| 欧美乱色亚洲激情| 此物有八面人人有两片| 淫秽高清视频在线观看| 久久久久性生活片| 看黄色毛片网站| 中文字幕av成人在线电影| 九色成人免费人妻av| 国产精华一区二区三区| 精品久久久久久久久久免费视频| 动漫黄色视频在线观看| 88av欧美| 99在线视频只有这里精品首页| 淫妇啪啪啪对白视频| 日本黄色视频三级网站网址| 中文字幕人妻丝袜一区二区| 久9热在线精品视频| 人妻丰满熟妇av一区二区三区| 亚洲电影在线观看av| 99热只有精品国产| 久久久久国内视频| 男女那种视频在线观看| 久久精品91蜜桃| 国产私拍福利视频在线观看| 手机成人av网站| 12—13女人毛片做爰片一| 日本精品一区二区三区蜜桃| 麻豆国产av国片精品| 一级毛片女人18水好多| 精品电影一区二区在线| 淫秽高清视频在线观看| 桃色一区二区三区在线观看| 亚洲欧美日韩东京热| 在线视频色国产色| 啪啪无遮挡十八禁网站| 丰满人妻熟妇乱又伦精品不卡| 国语自产精品视频在线第100页| 一区福利在线观看| 中亚洲国语对白在线视频| 亚洲av电影不卡..在线观看| 成年版毛片免费区| 香蕉av资源在线| 日韩高清综合在线| 精品一区二区三区av网在线观看| 一进一出抽搐gif免费好疼| 亚洲内射少妇av| 操出白浆在线播放| 久久精品国产综合久久久| av中文乱码字幕在线| 制服丝袜大香蕉在线| 亚洲第一欧美日韩一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 97人妻精品一区二区三区麻豆| 91在线精品国自产拍蜜月 | 他把我摸到了高潮在线观看| 亚洲av免费高清在线观看| 高潮久久久久久久久久久不卡| 日本黄大片高清| 亚洲av成人av| 精品无人区乱码1区二区| 香蕉久久夜色| 久久亚洲精品不卡| 男插女下体视频免费在线播放| 天堂网av新在线| 久久久久久人人人人人| 日韩欧美免费精品| 免费看日本二区| www日本在线高清视频| 在线观看一区二区三区| a在线观看视频网站| 精品一区二区三区av网在线观看| 国产三级中文精品| 性色av乱码一区二区三区2| 婷婷丁香在线五月| 亚洲专区国产一区二区| 亚洲熟妇中文字幕五十中出| 国内精品久久久久久久电影| 少妇熟女aⅴ在线视频| 国产高清激情床上av| 在线观看午夜福利视频| 窝窝影院91人妻| 日本五十路高清| 搡老妇女老女人老熟妇| av黄色大香蕉| 一区二区三区高清视频在线| 中文字幕熟女人妻在线| 午夜福利成人在线免费观看| 三级男女做爰猛烈吃奶摸视频| 国产精品,欧美在线| 欧美成狂野欧美在线观看| 精品久久久久久久人妻蜜臀av| 午夜免费成人在线视频| 久久亚洲真实| 亚洲在线自拍视频| 人人妻人人澡欧美一区二区| 亚洲av熟女| www.熟女人妻精品国产| 成人av一区二区三区在线看| 午夜亚洲福利在线播放| 国产综合懂色| 日韩精品青青久久久久久| 欧美一区二区亚洲| 亚洲精品乱码久久久v下载方式 | 国产精品亚洲av一区麻豆| 国内精品一区二区在线观看| 欧美黄色片欧美黄色片| x7x7x7水蜜桃| 欧美日韩瑟瑟在线播放| 亚洲18禁久久av| 国产精品美女特级片免费视频播放器| 免费看光身美女| 黄色女人牲交| 三级毛片av免费| 搞女人的毛片| 丰满人妻一区二区三区视频av | 国产亚洲精品久久久久久毛片| 十八禁网站免费在线| 三级毛片av免费| 成人高潮视频无遮挡免费网站| 一进一出好大好爽视频| 少妇人妻精品综合一区二区 | 99视频精品全部免费 在线| 天堂√8在线中文| 国产中年淑女户外野战色| 亚洲va日本ⅴa欧美va伊人久久| 国产精品一及| 亚洲成a人片在线一区二区| 激情在线观看视频在线高清| 中文字幕人妻丝袜一区二区| 美女高潮的动态| 精品福利观看| 国产精品三级大全| 成人亚洲精品av一区二区| 精品一区二区三区视频在线观看免费| 香蕉丝袜av| 午夜福利成人在线免费观看| 一进一出抽搐gif免费好疼| 观看美女的网站| 尤物成人国产欧美一区二区三区| 欧美一区二区精品小视频在线| 亚洲av第一区精品v没综合| 国产又黄又爽又无遮挡在线| 在线观看av片永久免费下载| 精品免费久久久久久久清纯| 日韩大尺度精品在线看网址| 国产一区在线观看成人免费| 成人av在线播放网站| 无限看片的www在线观看| 久久久成人免费电影| 亚洲人成电影免费在线| 国产精品一区二区三区四区免费观看 | 桃红色精品国产亚洲av| 午夜日韩欧美国产| 欧美午夜高清在线| 一级黄片播放器| 亚洲精品一卡2卡三卡4卡5卡| 长腿黑丝高跟| 欧美成人免费av一区二区三区| 国产午夜福利久久久久久| 最近视频中文字幕2019在线8| 亚洲久久久久久中文字幕| 天天躁日日操中文字幕| 午夜福利成人在线免费观看| 岛国在线观看网站| 91av网一区二区| 午夜亚洲福利在线播放| 成年免费大片在线观看| 韩国av一区二区三区四区| 久久欧美精品欧美久久欧美| 成年女人永久免费观看视频| 国产精品美女特级片免费视频播放器| 国产欧美日韩精品一区二区| 悠悠久久av| 久久久久久国产a免费观看| 97超级碰碰碰精品色视频在线观看| 一级毛片高清免费大全| 人人妻,人人澡人人爽秒播| 最近视频中文字幕2019在线8| 色视频www国产| 日本熟妇午夜| 国产免费av片在线观看野外av| 最后的刺客免费高清国语| 黑人欧美特级aaaaaa片| 国产午夜福利久久久久久| 亚洲av电影在线进入| 色尼玛亚洲综合影院| 91字幕亚洲| 中亚洲国语对白在线视频| 小蜜桃在线观看免费完整版高清| 美女黄网站色视频| 成熟少妇高潮喷水视频| 俺也久久电影网| 91在线观看av| 国产高清视频在线观看网站| 熟女少妇亚洲综合色aaa.| 国产爱豆传媒在线观看| 国产成人影院久久av| 亚洲一区二区三区不卡视频| 99在线人妻在线中文字幕| 亚洲av电影不卡..在线观看| av天堂在线播放| 亚洲五月婷婷丁香| 亚洲最大成人中文| 露出奶头的视频| av黄色大香蕉| 3wmmmm亚洲av在线观看| 久久久久性生活片| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久久久久人妻精品电影| 波多野结衣巨乳人妻| 男女做爰动态图高潮gif福利片| 99热6这里只有精品| 国产老妇女一区| 少妇裸体淫交视频免费看高清| 夜夜躁狠狠躁天天躁| 啦啦啦免费观看视频1| 国产美女午夜福利| 亚洲电影在线观看av| 亚洲av熟女| 色综合婷婷激情| 亚洲激情在线av| 性欧美人与动物交配| 国产一区二区三区视频了| 久久精品影院6| 18美女黄网站色大片免费观看| 日本a在线网址| 日本黄色片子视频| 久久久久国产精品人妻aⅴ院| 精品国产美女av久久久久小说| 麻豆成人午夜福利视频| 窝窝影院91人妻| 国产乱人视频| 淫秽高清视频在线观看| 首页视频小说图片口味搜索| 国产精品国产高清国产av| 香蕉av资源在线| 欧美日韩福利视频一区二区| 18+在线观看网站| 一级a爱片免费观看的视频| 在线观看午夜福利视频| 啦啦啦观看免费观看视频高清| 男人舔奶头视频| 在线观看美女被高潮喷水网站 | 久久久久久久精品吃奶| 在线播放国产精品三级| 欧美三级亚洲精品| 国产一区二区亚洲精品在线观看| 久久久国产成人精品二区| av福利片在线观看| 日韩欧美三级三区| 国产v大片淫在线免费观看| 嫩草影视91久久| 亚洲av免费高清在线观看| 国产精品98久久久久久宅男小说| 欧美zozozo另类| 在线播放无遮挡| 国产国拍精品亚洲av在线观看 | 少妇人妻一区二区三区视频| xxxwww97欧美| 国产三级中文精品| 一进一出抽搐动态| 成年女人看的毛片在线观看| 亚洲五月天丁香| 国产探花极品一区二区| 亚洲av一区综合| 亚洲中文日韩欧美视频| 天堂动漫精品| 国产色爽女视频免费观看| 午夜福利欧美成人| 欧美中文日本在线观看视频| 免费看a级黄色片| 国产真实乱freesex| 夜夜看夜夜爽夜夜摸| 天天躁日日操中文字幕| 老司机在亚洲福利影院| 亚洲黑人精品在线| 色老头精品视频在线观看| 深爱激情五月婷婷| 久久久久免费精品人妻一区二区| 乱人视频在线观看| 亚洲午夜理论影院| 露出奶头的视频| 日日摸夜夜添夜夜添小说| 精品乱码久久久久久99久播| 黑人欧美特级aaaaaa片| 五月玫瑰六月丁香| 国产在线精品亚洲第一网站| 精品乱码久久久久久99久播| 国产在线精品亚洲第一网站| 日本精品一区二区三区蜜桃| 99热这里只有是精品50| 黄片小视频在线播放| www.熟女人妻精品国产| 国产欧美日韩一区二区三| 欧洲精品卡2卡3卡4卡5卡区| 欧美bdsm另类| 熟女少妇亚洲综合色aaa.| 亚洲精品日韩av片在线观看 | 国产高清视频在线播放一区| 午夜免费激情av| 国产伦一二天堂av在线观看| 丰满的人妻完整版| www.999成人在线观看| 天堂√8在线中文| 国产视频内射| 91九色精品人成在线观看| www.色视频.com| av女优亚洲男人天堂| 91在线观看av| av欧美777| 熟女人妻精品中文字幕| 亚洲最大成人中文| 天堂√8在线中文| 午夜福利成人在线免费观看| 12—13女人毛片做爰片一| 最近最新中文字幕大全免费视频| 在线观看美女被高潮喷水网站 | 国产欧美日韩一区二区三| 香蕉av资源在线| 在线天堂最新版资源| 女人十人毛片免费观看3o分钟| 18+在线观看网站| 97超视频在线观看视频| 老司机午夜十八禁免费视频| 午夜两性在线视频| 国产精品精品国产色婷婷| 免费观看精品视频网站| 久久亚洲精品不卡| 一进一出抽搐动态| 国产老妇女一区| 国产三级黄色录像| 亚洲欧美精品综合久久99| 国产亚洲精品久久久久久毛片| 99国产综合亚洲精品| 级片在线观看| 精品人妻1区二区| 90打野战视频偷拍视频| 欧美黄色片欧美黄色片| 国产视频内射| 亚洲 欧美 日韩 在线 免费| 亚洲aⅴ乱码一区二区在线播放| 窝窝影院91人妻| 日本 av在线| 三级毛片av免费| 白带黄色成豆腐渣| 亚洲av电影在线进入| 观看免费一级毛片| 中文字幕人成人乱码亚洲影| 成人性生交大片免费视频hd| 欧美最新免费一区二区三区 | 在线观看av片永久免费下载| 国模一区二区三区四区视频| 内射极品少妇av片p| 一二三四社区在线视频社区8| 国产精品久久电影中文字幕| 久久久久亚洲av毛片大全| 国产亚洲精品久久久com| 无限看片的www在线观看| 少妇的逼好多水| 国内久久婷婷六月综合欲色啪| 国产三级在线视频| 久久九九热精品免费| 成人国产综合亚洲| a级毛片a级免费在线| 最近最新中文字幕大全电影3| 中文字幕人成人乱码亚洲影| 亚洲欧美激情综合另类| 嫩草影视91久久| 亚洲国产色片| 99视频精品全部免费 在线| 亚洲国产精品成人综合色| 在线播放无遮挡| 国产成人福利小说| www.999成人在线观看| 美女 人体艺术 gogo| 又紧又爽又黄一区二区| 亚洲欧美一区二区三区黑人| 国产精品一区二区免费欧美| 国产真人三级小视频在线观看| 亚洲人成电影免费在线| 亚洲国产高清在线一区二区三| 叶爱在线成人免费视频播放| 久久久久精品国产欧美久久久| 久久欧美精品欧美久久欧美| 国产亚洲精品久久久com| 色老头精品视频在线观看| 极品教师在线免费播放| 亚洲人成网站高清观看| 久久亚洲精品不卡| 亚洲一区高清亚洲精品| 成人av一区二区三区在线看| 日韩欧美精品免费久久 | 成人精品一区二区免费| 亚洲激情在线av| 亚洲人成网站在线播| 日本与韩国留学比较| 亚洲国产色片| 全区人妻精品视频| 91久久精品国产一区二区成人 | 在线观看一区二区三区| 国产单亲对白刺激| 嫩草影视91久久| 欧美色视频一区免费| 亚洲av二区三区四区| 一个人免费在线观看的高清视频| 叶爱在线成人免费视频播放| 毛片女人毛片| 韩国av一区二区三区四区| 国产亚洲欧美在线一区二区| 男女那种视频在线观看| 国产精品久久久久久精品电影| 国产v大片淫在线免费观看| 欧美+亚洲+日韩+国产| 久久九九热精品免费| 亚洲专区国产一区二区| 天美传媒精品一区二区| 国产亚洲精品综合一区在线观看| АⅤ资源中文在线天堂| 最新美女视频免费是黄的| 老鸭窝网址在线观看| 一个人看视频在线观看www免费 | 黄色丝袜av网址大全| 欧美在线黄色| 国产亚洲欧美98| 日韩欧美国产一区二区入口| 51国产日韩欧美| 国产爱豆传媒在线观看| 在线观看免费午夜福利视频| 麻豆久久精品国产亚洲av| 最后的刺客免费高清国语| 亚洲精品一卡2卡三卡4卡5卡| 丝袜美腿在线中文| 亚洲国产精品合色在线| 免费看美女性在线毛片视频| 午夜日韩欧美国产| 国产av不卡久久| 欧美极品一区二区三区四区| 老熟妇仑乱视频hdxx| 在线观看午夜福利视频| 一个人看的www免费观看视频| 日本免费a在线| 好看av亚洲va欧美ⅴa在| 欧美性猛交黑人性爽| 啦啦啦韩国在线观看视频| 88av欧美| 欧美成人免费av一区二区三区| 人人妻人人澡欧美一区二区| bbb黄色大片| 国产高清视频在线播放一区| 免费人成视频x8x8入口观看| 欧美乱色亚洲激情| a级毛片a级免费在线| 丁香六月欧美| 天堂动漫精品| 国产aⅴ精品一区二区三区波| 久久婷婷人人爽人人干人人爱| 亚洲久久久久久中文字幕| 久久精品国产亚洲av香蕉五月| 法律面前人人平等表现在哪些方面| 国内揄拍国产精品人妻在线| 国产成人av激情在线播放| 极品教师在线免费播放| 国产精品1区2区在线观看.| 日韩欧美精品免费久久 | 精品日产1卡2卡| avwww免费| 国内久久婷婷六月综合欲色啪| 一a级毛片在线观看| 女人被狂操c到高潮| 成人国产一区最新在线观看| 欧美色视频一区免费| 国产99白浆流出| 亚洲最大成人手机在线| 国产老妇女一区| 欧美xxxx黑人xx丫x性爽| 99视频精品全部免费 在线| 51国产日韩欧美| 午夜a级毛片| 欧美最黄视频在线播放免费| 69人妻影院| 亚洲 欧美 日韩 在线 免费| 操出白浆在线播放| 一区福利在线观看| 看片在线看免费视频| 69av精品久久久久久| 12—13女人毛片做爰片一| 国产v大片淫在线免费观看| 美女高潮喷水抽搐中文字幕| 欧美日韩国产亚洲二区| 俺也久久电影网| 亚洲真实伦在线观看| 精品人妻一区二区三区麻豆 | 一卡2卡三卡四卡精品乱码亚洲| 男女床上黄色一级片免费看| 欧美乱色亚洲激情| 亚洲最大成人手机在线| 亚洲国产精品999在线| 99热这里只有是精品50| 搡老妇女老女人老熟妇| av在线蜜桃| 色综合欧美亚洲国产小说| 波多野结衣高清作品| 麻豆一二三区av精品| 国产精品精品国产色婷婷| 美女高潮的动态| 国语自产精品视频在线第100页| 亚洲欧美一区二区三区黑人| 亚洲欧美日韩东京热| 99国产精品一区二区蜜桃av| 亚洲av二区三区四区| 日本 欧美在线| 免费无遮挡裸体视频| 母亲3免费完整高清在线观看| 一级黄色大片毛片| 久久九九热精品免费| 最好的美女福利视频网| 亚洲天堂国产精品一区在线| 超碰av人人做人人爽久久 | 亚洲精品美女久久久久99蜜臀| 免费看光身美女| 国产淫片久久久久久久久 | 国产欧美日韩精品亚洲av| 一个人看视频在线观看www免费 | 亚洲精品日韩av片在线观看 | 宅男免费午夜| 成人永久免费在线观看视频| 日本一二三区视频观看| 久久国产乱子伦精品免费另类| 淫妇啪啪啪对白视频| 欧美成人性av电影在线观看| 日韩免费av在线播放| www.熟女人妻精品国产| 国产精品精品国产色婷婷| 一区二区三区激情视频| 亚洲精品456在线播放app | 久久久久久久久大av| 国产午夜精品论理片| 久久午夜亚洲精品久久| 男女下面进入的视频免费午夜| 午夜视频国产福利| a在线观看视频网站| 成人三级黄色视频| 国产av麻豆久久久久久久| 中亚洲国语对白在线视频| av中文乱码字幕在线| 天堂网av新在线| 色播亚洲综合网| 神马国产精品三级电影在线观看| 九色国产91popny在线| 欧美成人性av电影在线观看| 大型黄色视频在线免费观看| 午夜久久久久精精品| 18禁美女被吸乳视频| 日本免费一区二区三区高清不卡| 成人18禁在线播放| 亚洲狠狠婷婷综合久久图片| 亚洲中文字幕日韩| 欧美一级毛片孕妇| 色综合站精品国产| 中文亚洲av片在线观看爽| 91麻豆av在线| 久久久久久大精品| 精品一区二区三区视频在线观看免费| 成人亚洲精品av一区二区| a级毛片a级免费在线| 成年女人看的毛片在线观看| 欧美日韩亚洲国产一区二区在线观看| 露出奶头的视频| 国产极品精品免费视频能看的| 精品99又大又爽又粗少妇毛片 | 变态另类丝袜制服| 变态另类成人亚洲欧美熟女| www国产在线视频色| 欧美日韩中文字幕国产精品一区二区三区| 亚洲人成网站在线播放欧美日韩| 欧美一级毛片孕妇| 中文字幕熟女人妻在线| 三级国产精品欧美在线观看|