• <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)

    亚洲av在线观看美女高潮| 久久人人爽人人爽人人片va| 嫩草影院入口| 成人亚洲精品av一区二区| 美女被艹到高潮喷水动态| 夫妻性生交免费视频一级片| 亚洲av成人精品一二三区| 视频中文字幕在线观看| 免费在线观看成人毛片| 久热这里只有精品99| 三级经典国产精品| 精品国产三级普通话版| 卡戴珊不雅视频在线播放| 国产亚洲一区二区精品| 亚洲,欧美,日韩| 91午夜精品亚洲一区二区三区| 国产成年人精品一区二区| 好男人视频免费观看在线| 中文字幕制服av| 人人妻人人爽人人添夜夜欢视频 | 国产老妇伦熟女老妇高清| 亚洲国产精品国产精品| 成人亚洲欧美一区二区av| eeuss影院久久| 亚洲精品国产成人久久av| 日韩亚洲欧美综合| 国产69精品久久久久777片| 亚洲成人久久爱视频| 又黄又爽又刺激的免费视频.| h日本视频在线播放| 亚洲精品成人久久久久久| 久久影院123| 国产 精品1| 国产一区有黄有色的免费视频| 欧美日韩综合久久久久久| 日韩伦理黄色片| 久久久久久久精品精品| 久久久久久久国产电影| 看十八女毛片水多多多| 亚洲综合色惰| 日本欧美国产在线视频| 精品视频人人做人人爽| 午夜日本视频在线| 国产精品精品国产色婷婷| 国产黄色视频一区二区在线观看| 免费av观看视频| 99热国产这里只有精品6| 久久久久久国产a免费观看| 在现免费观看毛片| 久热这里只有精品99| 久久久久久九九精品二区国产| 中文精品一卡2卡3卡4更新| 亚洲欧美清纯卡通| 天堂网av新在线| 亚洲电影在线观看av| 日本wwww免费看| 欧美日韩视频高清一区二区三区二| 国产乱来视频区| 中文资源天堂在线| 精华霜和精华液先用哪个| 国产黄片视频在线免费观看| 日本三级黄在线观看| 三级经典国产精品| 国产欧美另类精品又又久久亚洲欧美| 久久久久精品性色| 97超碰精品成人国产| 亚洲av国产av综合av卡| 在线观看美女被高潮喷水网站| 国产精品一二三区在线看| 国产在线男女| 少妇人妻一区二区三区视频| 在线观看免费高清a一片| 成人亚洲精品av一区二区| 国产毛片a区久久久久| 3wmmmm亚洲av在线观看| 好男人在线观看高清免费视频| 高清在线视频一区二区三区| 成人毛片60女人毛片免费| 一级片'在线观看视频| 亚洲国产色片| 97超碰精品成人国产| 免费看a级黄色片| 中国国产av一级| 性色av一级| 国产伦在线观看视频一区| 中文欧美无线码| 欧美日韩视频高清一区二区三区二| 国产一级毛片在线| 久久久久久久久久人人人人人人| 自拍欧美九色日韩亚洲蝌蚪91 | 91在线精品国自产拍蜜月| 欧美成人精品欧美一级黄| 九草在线视频观看| 日韩三级伦理在线观看| 国产中年淑女户外野战色| 久久久欧美国产精品| 亚洲精品久久久久久婷婷小说| 青春草亚洲视频在线观看| 国产精品国产三级国产专区5o| av在线亚洲专区| 亚洲熟女精品中文字幕| 欧美日韩国产mv在线观看视频 | 亚洲国产高清在线一区二区三| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 特级一级黄色大片| 亚洲美女视频黄频| 色婷婷久久久亚洲欧美| 尾随美女入室| 欧美日韩视频高清一区二区三区二| 18禁动态无遮挡网站| 黄色一级大片看看| 国产成人免费观看mmmm| 97精品久久久久久久久久精品| 熟女av电影| 国产 一区精品| 卡戴珊不雅视频在线播放| 97在线人人人人妻| 男人舔奶头视频| 美女xxoo啪啪120秒动态图| 一级av片app| 国产黄片美女视频| 夫妻午夜视频| 在线亚洲精品国产二区图片欧美 | 国产精品久久久久久精品电影小说 | 久久精品国产a三级三级三级| 一级毛片黄色毛片免费观看视频| 国产高清国产精品国产三级 | 精品少妇黑人巨大在线播放| 欧美精品人与动牲交sv欧美| 精品久久久久久电影网| 男人爽女人下面视频在线观看| 欧美变态另类bdsm刘玥| 国产 一区精品| 王馨瑶露胸无遮挡在线观看| 赤兔流量卡办理| 亚洲自拍偷在线| 色播亚洲综合网| 老师上课跳d突然被开到最大视频| 亚洲欧美精品自产自拍| 少妇裸体淫交视频免费看高清| 热re99久久精品国产66热6| 少妇高潮的动态图| 国产成人福利小说| 女的被弄到高潮叫床怎么办| 亚洲电影在线观看av| 久久精品综合一区二区三区| 欧美高清性xxxxhd video| 亚洲av二区三区四区| 国产一区亚洲一区在线观看| 午夜爱爱视频在线播放| 伦理电影大哥的女人| 99久久中文字幕三级久久日本| 男人爽女人下面视频在线观看| 熟女av电影| 国产精品秋霞免费鲁丝片| 亚洲四区av| 免费av观看视频| 少妇丰满av| 婷婷色综合www| 免费高清在线观看视频在线观看| 国产成人免费无遮挡视频| 国产一区有黄有色的免费视频| 2022亚洲国产成人精品| 免费观看无遮挡的男女| 秋霞伦理黄片| 男女下面进入的视频免费午夜| 免费观看的影片在线观看| 久久这里有精品视频免费| 男女边摸边吃奶| 欧美激情在线99| 亚洲精品乱码久久久v下载方式| av.在线天堂| 噜噜噜噜噜久久久久久91| 熟妇人妻不卡中文字幕| 欧美精品一区二区大全| 色播亚洲综合网| 亚洲三级黄色毛片| 搡女人真爽免费视频火全软件| 好男人在线观看高清免费视频| 麻豆成人av视频| 久久久成人免费电影| 久久国内精品自在自线图片| 18+在线观看网站| 视频区图区小说| 中文乱码字字幕精品一区二区三区| 熟女av电影| 国产成人91sexporn| 国产精品99久久99久久久不卡 | 一区二区三区免费毛片| 中国美白少妇内射xxxbb| 草草在线视频免费看| 五月天丁香电影| 久久久久久久久久成人| 九色成人免费人妻av| 国产精品伦人一区二区| 91久久精品电影网| 97热精品久久久久久| 国产伦精品一区二区三区四那| 26uuu在线亚洲综合色| 男人和女人高潮做爰伦理| 91精品一卡2卡3卡4卡| 国内揄拍国产精品人妻在线| 久久久久性生活片| 制服丝袜香蕉在线| 精品久久久久久久末码| 搞女人的毛片| 欧美少妇被猛烈插入视频| 亚洲人成网站在线播| 永久免费av网站大全| 国产黄片视频在线免费观看| 小蜜桃在线观看免费完整版高清| 成人国产av品久久久| 一级av片app| 欧美三级亚洲精品| 女人十人毛片免费观看3o分钟| 日韩视频在线欧美| 久久99热这里只有精品18| 亚洲欧美一区二区三区国产| 日本免费在线观看一区| 天天躁日日操中文字幕| 国产毛片在线视频| 有码 亚洲区| 99热这里只有精品一区| 只有这里有精品99| 在线免费十八禁| 黄色配什么色好看| av在线天堂中文字幕| 老女人水多毛片| 国产真实伦视频高清在线观看| 久久久色成人| 国产高清不卡午夜福利| 国产高清有码在线观看视频| 国产成人午夜福利电影在线观看| 大码成人一级视频| 精品亚洲乱码少妇综合久久| 美女主播在线视频| 国产乱来视频区| 天堂中文最新版在线下载 | 免费看不卡的av| 亚洲av男天堂| 亚洲人成网站在线播| 久久精品人妻少妇| 亚洲国产最新在线播放| 日本熟妇午夜| 久久综合国产亚洲精品| 人人妻人人看人人澡| 久久久久久伊人网av| 成年免费大片在线观看| 中文字幕av成人在线电影| 大片免费播放器 马上看| 熟女av电影| xxx大片免费视频| 国产91av在线免费观看| 自拍偷自拍亚洲精品老妇| 亚洲最大成人手机在线| 黄色日韩在线| av在线播放精品| 久热这里只有精品99| 国产精品一及| 亚洲av不卡在线观看| 在线观看av片永久免费下载| 老司机影院成人| 亚洲av二区三区四区| 18禁在线播放成人免费| a级一级毛片免费在线观看| 久久久久久久久大av| 亚洲成人一二三区av| 我的老师免费观看完整版| 少妇猛男粗大的猛烈进出视频 | videossex国产| 尾随美女入室| 一级毛片电影观看| 日韩一区二区视频免费看| 亚洲国产精品专区欧美| 欧美极品一区二区三区四区| 国产精品无大码| 亚洲成人久久爱视频| 高清午夜精品一区二区三区| 好男人视频免费观看在线| 欧美人与善性xxx| 有码 亚洲区| 观看免费一级毛片| 七月丁香在线播放| 女的被弄到高潮叫床怎么办| 亚洲天堂av无毛| 欧美xxxx性猛交bbbb| 欧美日韩一区二区视频在线观看视频在线 | 97精品久久久久久久久久精品| 免费观看av网站的网址| 97热精品久久久久久| 久久鲁丝午夜福利片| 丰满乱子伦码专区| av专区在线播放| 日韩大片免费观看网站| 国产综合懂色| 免费av毛片视频| 香蕉精品网在线| 亚洲精品久久久久久婷婷小说| 国产老妇伦熟女老妇高清| 99视频精品全部免费 在线| 亚洲aⅴ乱码一区二区在线播放| 午夜免费观看性视频| 亚洲精品乱久久久久久| 真实男女啪啪啪动态图| av免费在线看不卡| 日日啪夜夜爽| 国产成人91sexporn| 在线亚洲精品国产二区图片欧美 | 中文在线观看免费www的网站| 3wmmmm亚洲av在线观看| 18禁动态无遮挡网站| 亚洲av国产av综合av卡| 少妇人妻精品综合一区二区| 亚洲精品视频女| 国产片特级美女逼逼视频| 国产高清不卡午夜福利| 日韩欧美精品免费久久| 99视频精品全部免费 在线| 少妇猛男粗大的猛烈进出视频 | 久久国内精品自在自线图片| 在线天堂最新版资源| 国产成人a∨麻豆精品| 免费看光身美女| 激情 狠狠 欧美| 久久这里有精品视频免费| 777米奇影视久久| 老女人水多毛片| 精品一区二区免费观看| 直男gayav资源| 日韩伦理黄色片| 国产精品伦人一区二区| 精华霜和精华液先用哪个| 国产69精品久久久久777片| 神马国产精品三级电影在线观看| 亚洲欧美日韩另类电影网站 | 在线观看一区二区三区| 欧美高清性xxxxhd video| 国产一区二区在线观看日韩| 在线亚洲精品国产二区图片欧美 | 99热这里只有是精品50| 国产精品人妻久久久影院| 涩涩av久久男人的天堂| 亚洲一级一片aⅴ在线观看| 久久精品国产a三级三级三级| 男女下面进入的视频免费午夜| 特大巨黑吊av在线直播| 成人国产麻豆网| 久久久久久久国产电影| 九九久久精品国产亚洲av麻豆| 午夜精品国产一区二区电影 | 国产中年淑女户外野战色| 精品人妻偷拍中文字幕| 亚洲欧美一区二区三区国产| 国产黄色免费在线视频| 青春草国产在线视频| 国产精品国产三级专区第一集| 亚洲人与动物交配视频| 欧美亚洲 丝袜 人妻 在线| 97精品久久久久久久久久精品| 国产乱人视频| 免费在线观看成人毛片| 国产午夜精品久久久久久一区二区三区| 三级国产精品片| 一区二区三区精品91| 国产 精品1| 国内少妇人妻偷人精品xxx网站| 下体分泌物呈黄色| 亚洲精品亚洲一区二区| 久久精品国产亚洲网站| 欧美日韩视频精品一区| 特大巨黑吊av在线直播| 草草在线视频免费看| av线在线观看网站| 97超视频在线观看视频| 精品一区二区三卡| 久久久久久久午夜电影| 少妇的逼水好多| 国产一区二区亚洲精品在线观看| 亚洲人成网站在线观看播放| 舔av片在线| 乱码一卡2卡4卡精品| 久久久久久久久大av| 中文在线观看免费www的网站| av免费在线看不卡| 夜夜爽夜夜爽视频| av一本久久久久| 久久久欧美国产精品| 嫩草影院精品99| 日日摸夜夜添夜夜添av毛片| 观看美女的网站| 免费看光身美女| 欧美激情在线99| 国产白丝娇喘喷水9色精品| 色播亚洲综合网| 人人妻人人澡人人爽人人夜夜| 大片免费播放器 马上看| 国产高清不卡午夜福利| 久久久色成人| 午夜福利高清视频| 韩国av在线不卡| 国产精品偷伦视频观看了| 欧美丝袜亚洲另类| 欧美日韩精品成人综合77777| 亚洲人成网站在线播| 日韩制服骚丝袜av| 久久久久久久亚洲中文字幕| 国产精品三级大全| 麻豆成人av视频| 免费观看性生交大片5| 成人免费观看视频高清| 高清欧美精品videossex| 国产成人精品久久久久久| 久久久久精品性色| 2021天堂中文幕一二区在线观| 超碰97精品在线观看| av免费观看日本| 内地一区二区视频在线| 99热网站在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲精品国产成人久久av| av在线亚洲专区| 国产高清有码在线观看视频| 国产老妇女一区| 毛片女人毛片| 亚洲精品日韩在线中文字幕| 欧美 日韩 精品 国产| 3wmmmm亚洲av在线观看| 老女人水多毛片| 啦啦啦在线观看免费高清www| 成年女人在线观看亚洲视频 | 久久久久久久久久久免费av| 毛片女人毛片| freevideosex欧美| 午夜精品国产一区二区电影 | 国产欧美日韩一区二区三区在线 | 久久久欧美国产精品| 在线观看美女被高潮喷水网站| 欧美激情久久久久久爽电影| 一区二区三区乱码不卡18| 1000部很黄的大片| 美女高潮的动态| 欧美日韩精品成人综合77777| 亚洲精品一区蜜桃| 亚洲内射少妇av| 日日摸夜夜添夜夜添av毛片| 午夜福利视频精品| 在线看a的网站| 久久久精品欧美日韩精品| 看黄色毛片网站| 精品少妇久久久久久888优播| 亚洲色图av天堂| 2021少妇久久久久久久久久久| 中文字幕av成人在线电影| 啦啦啦啦在线视频资源| 国产精品秋霞免费鲁丝片| 国产成人精品一,二区| 欧美国产精品一级二级三级 | 成人午夜精彩视频在线观看| 欧美日韩视频高清一区二区三区二| 白带黄色成豆腐渣| 一本一本综合久久| 国产极品天堂在线| 熟妇人妻不卡中文字幕| 大片电影免费在线观看免费| 久久久久久久久久成人| 精品久久久精品久久久| 22中文网久久字幕| 亚洲精品456在线播放app| 国产精品秋霞免费鲁丝片| 国产免费又黄又爽又色| 青春草亚洲视频在线观看| 白带黄色成豆腐渣| 国产淫片久久久久久久久| 2021天堂中文幕一二区在线观| 蜜桃久久精品国产亚洲av| 大码成人一级视频| 国产伦在线观看视频一区| 精品久久久久久久末码| 天天躁夜夜躁狠狠久久av| 亚洲精品456在线播放app| 99热网站在线观看| 国产精品国产av在线观看| 亚洲欧美日韩卡通动漫| 大话2 男鬼变身卡| 国产69精品久久久久777片| 99久久精品热视频| 久久这里有精品视频免费| 亚洲精品久久午夜乱码| 亚洲国产高清在线一区二区三| 午夜激情福利司机影院| 少妇丰满av| 欧美97在线视频| 在线观看美女被高潮喷水网站| 大陆偷拍与自拍| 国产精品一区www在线观看| 深夜a级毛片| 乱码一卡2卡4卡精品| 亚洲精品国产成人久久av| 白带黄色成豆腐渣| 午夜爱爱视频在线播放| 日韩伦理黄色片| 精品一区二区三卡| 91在线精品国自产拍蜜月| 成年女人看的毛片在线观看| 青春草国产在线视频| 精品一区二区免费观看| 午夜免费男女啪啪视频观看| 亚洲婷婷狠狠爱综合网| 亚洲精品乱码久久久久久按摩| 亚洲人成网站在线播| 亚洲内射少妇av| 国产一区二区三区综合在线观看 | 中文字幕久久专区| 久久久久久久精品精品| 寂寞人妻少妇视频99o| 韩国高清视频一区二区三区| 日韩大片免费观看网站| 国产毛片在线视频| 美女cb高潮喷水在线观看| 一区二区三区精品91| 国产日韩欧美亚洲二区| 激情 狠狠 欧美| 久久久久久久久久人人人人人人| 亚洲欧美日韩另类电影网站 | 国产亚洲午夜精品一区二区久久 | 国产精品熟女久久久久浪| 国产真实伦视频高清在线观看| tube8黄色片| 国产精品精品国产色婷婷| 欧美区成人在线视频| 日韩一区二区三区影片| 亚洲国产精品专区欧美| 欧美日本视频| 精品久久久精品久久久| 久久99蜜桃精品久久| 乱码一卡2卡4卡精品| 一级二级三级毛片免费看| 大陆偷拍与自拍| a级毛色黄片| 国产成人免费观看mmmm| 人妻一区二区av| 春色校园在线视频观看| 成人黄色视频免费在线看| 黄色欧美视频在线观看| 女人十人毛片免费观看3o分钟| 成人毛片a级毛片在线播放| 国产91av在线免费观看| 欧美另类一区| 最近中文字幕2019免费版| 麻豆乱淫一区二区| 七月丁香在线播放| 欧美激情国产日韩精品一区| 国产精品嫩草影院av在线观看| 精华霜和精华液先用哪个| 男人添女人高潮全过程视频| 女人十人毛片免费观看3o分钟| 男的添女的下面高潮视频| 亚洲人成网站在线播| 久久久精品94久久精品| 三级男女做爰猛烈吃奶摸视频| 麻豆乱淫一区二区| 爱豆传媒免费全集在线观看| 国产成人freesex在线| av女优亚洲男人天堂| 少妇被粗大猛烈的视频| 午夜视频国产福利| 日韩欧美精品免费久久| 日韩欧美精品v在线| 精品久久久噜噜| 91在线精品国自产拍蜜月| 男的添女的下面高潮视频| av线在线观看网站| av在线天堂中文字幕| 国产成人精品久久久久久| 在线亚洲精品国产二区图片欧美 | 久久久久久久久久久免费av| 久久人人爽人人片av| 18禁裸乳无遮挡动漫免费视频 | 婷婷色综合www| 免费观看a级毛片全部| 欧美潮喷喷水| 好男人视频免费观看在线| 国产黄频视频在线观看| 精品久久久噜噜| 91精品伊人久久大香线蕉| 2022亚洲国产成人精品| 国产一区有黄有色的免费视频| 一区二区三区四区激情视频| 精品人妻偷拍中文字幕| 中文在线观看免费www的网站| 精品少妇久久久久久888优播| 国产爱豆传媒在线观看| av在线app专区| 午夜亚洲福利在线播放| 国内少妇人妻偷人精品xxx网站| av黄色大香蕉| 天美传媒精品一区二区| 亚洲精品aⅴ在线观看| 精品视频人人做人人爽| 精品久久久久久久末码| 日韩欧美精品v在线| av专区在线播放| 一区二区三区乱码不卡18| 日本午夜av视频| 少妇的逼水好多| 2021少妇久久久久久久久久久| 日韩精品有码人妻一区| 国产伦理片在线播放av一区| 麻豆成人午夜福利视频| 国产伦精品一区二区三区视频9| 人人妻人人爽人人添夜夜欢视频 | 成人特级av手机在线观看| 欧美日韩一区二区视频在线观看视频在线 | 日韩国内少妇激情av| 最近2019中文字幕mv第一页|