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

    久9热在线精品视频| 中文字幕精品免费在线观看视频| 真人一进一出gif抽搐免费| 极品教师在线免费播放| 中文字幕人妻丝袜一区二区| 一区二区日韩欧美中文字幕| 日韩欧美一区视频在线观看| 国产在线观看jvid| 国产一级毛片七仙女欲春2 | 99国产极品粉嫩在线观看| 身体一侧抽搐| 国产精品亚洲美女久久久| 真人做人爱边吃奶动态| 脱女人内裤的视频| 日日摸夜夜添夜夜添小说| 黑人欧美特级aaaaaa片| 国产久久久一区二区三区| 国产成人一区二区三区免费视频网站| 极品教师在线免费播放| 亚洲成av片中文字幕在线观看| 成人av一区二区三区在线看| 亚洲全国av大片| 麻豆久久精品国产亚洲av| 亚洲 欧美 日韩 在线 免费| 国产黄片美女视频| 曰老女人黄片| 亚洲国产欧美网| 成人18禁高潮啪啪吃奶动态图| 精品第一国产精品| 法律面前人人平等表现在哪些方面| 国产亚洲欧美98| 桃色一区二区三区在线观看| 首页视频小说图片口味搜索| 日本一区二区免费在线视频| 可以在线观看的亚洲视频| 热99re8久久精品国产| 国内精品久久久久久久电影| 久久九九热精品免费| 国产亚洲av嫩草精品影院| 成人手机av| 俺也久久电影网| 欧美乱妇无乱码| 欧美人与性动交α欧美精品济南到| 制服诱惑二区| 一级片免费观看大全| 母亲3免费完整高清在线观看| 午夜亚洲福利在线播放| 亚洲九九香蕉| 成人免费观看视频高清| 久久久久久大精品| 日韩欧美免费精品| 少妇裸体淫交视频免费看高清 | 黄色 视频免费看| 成人特级黄色片久久久久久久| 久久香蕉精品热| 精品一区二区三区视频在线观看免费| 国产精品综合久久久久久久免费| 麻豆久久精品国产亚洲av| 天天躁夜夜躁狠狠躁躁| 巨乳人妻的诱惑在线观看| 久久久国产欧美日韩av| 天天躁狠狠躁夜夜躁狠狠躁| 欧美av亚洲av综合av国产av| 亚洲国产欧美日韩在线播放| 国产精品久久久久久亚洲av鲁大| 免费看a级黄色片| 最好的美女福利视频网| 一进一出抽搐gif免费好疼| 91字幕亚洲| 自线自在国产av| 午夜免费成人在线视频| 久久精品国产亚洲av香蕉五月| 波多野结衣高清作品| 免费观看人在逋| 1024视频免费在线观看| 国产真实乱freesex| 老司机午夜十八禁免费视频| 久久伊人香网站| xxx96com| 一级毛片女人18水好多| 欧美最黄视频在线播放免费| 午夜福利一区二区在线看| 亚洲片人在线观看| 草草在线视频免费看| 久久久久久久久久黄片| 亚洲精品粉嫩美女一区| 国产一区二区在线av高清观看| 淫秽高清视频在线观看| 国产伦人伦偷精品视频| 人成视频在线观看免费观看| 一边摸一边抽搐一进一小说| 1024手机看黄色片| 1024手机看黄色片| 欧美黑人巨大hd| 18禁裸乳无遮挡免费网站照片 | 精品午夜福利视频在线观看一区| 狂野欧美激情性xxxx| 丰满人妻熟妇乱又伦精品不卡| 久久午夜综合久久蜜桃| 夜夜躁狠狠躁天天躁| 男男h啪啪无遮挡| 人人妻人人看人人澡| 日本熟妇午夜| 老司机在亚洲福利影院| 激情在线观看视频在线高清| 女同久久另类99精品国产91| 精品国产超薄肉色丝袜足j| 熟女电影av网| 日韩欧美 国产精品| 国产又爽黄色视频| 亚洲真实伦在线观看| 国产欧美日韩一区二区三| 久久久久久久午夜电影| 亚洲av第一区精品v没综合| 视频在线观看一区二区三区| 久久九九热精品免费| 久9热在线精品视频| 亚洲五月婷婷丁香| 亚洲色图 男人天堂 中文字幕| 国产亚洲av嫩草精品影院| 很黄的视频免费| 欧美中文综合在线视频| 国产aⅴ精品一区二区三区波| 久久国产亚洲av麻豆专区| 搞女人的毛片| 日本成人三级电影网站| 久久久久久人人人人人| 18禁黄网站禁片免费观看直播| 精品国产美女av久久久久小说| 国产av一区二区精品久久| 女生性感内裤真人,穿戴方法视频| 亚洲九九香蕉| 色精品久久人妻99蜜桃| 99精品久久久久人妻精品| 国产视频内射| 久久人妻av系列| 国产精品亚洲美女久久久| 黑人巨大精品欧美一区二区mp4| 亚洲三区欧美一区| 午夜视频精品福利| 国产又爽黄色视频| 他把我摸到了高潮在线观看| 久久久水蜜桃国产精品网| 男女之事视频高清在线观看| 桃红色精品国产亚洲av| 国产成人一区二区三区免费视频网站| 色av中文字幕| xxx96com| 熟妇人妻久久中文字幕3abv| 国产乱人伦免费视频| 久久久久久人人人人人| 97人妻精品一区二区三区麻豆 | 亚洲一码二码三码区别大吗| 日韩大码丰满熟妇| av福利片在线| 最新在线观看一区二区三区| 黄片播放在线免费| 亚洲国产欧美网| 国产亚洲精品久久久久久毛片| 免费看a级黄色片| 亚洲熟妇熟女久久| 久久亚洲真实| 很黄的视频免费| 精品国内亚洲2022精品成人| 熟妇人妻久久中文字幕3abv| 亚洲电影在线观看av| 在线观看一区二区三区| 成人亚洲精品一区在线观看| 91国产中文字幕| 国产成人欧美在线观看| 亚洲欧美日韩高清在线视频| 国产乱人伦免费视频| 久久久久亚洲av毛片大全| 久久久久久久久久黄片| 大型av网站在线播放| 美女 人体艺术 gogo| 国内精品久久久久精免费| 不卡av一区二区三区| 日韩精品青青久久久久久| 亚洲欧美日韩无卡精品| 亚洲国产欧洲综合997久久, | 中文字幕久久专区| 国产精品免费视频内射| 国产精品久久久久久亚洲av鲁大| 亚洲成人国产一区在线观看| 一二三四社区在线视频社区8| 在线观看午夜福利视频| 欧美色欧美亚洲另类二区| 一级黄色大片毛片| 曰老女人黄片| 国产精品久久久久久亚洲av鲁大| 亚洲精品色激情综合| 一夜夜www| 久99久视频精品免费| 午夜久久久久精精品| 国产亚洲精品av在线| 亚洲人成网站在线播放欧美日韩| 精品久久久久久,| 日韩三级视频一区二区三区| 精品久久久久久久末码| 久久午夜亚洲精品久久| 亚洲成av人片免费观看| 亚洲男人的天堂狠狠| 国产亚洲精品综合一区在线观看 | 午夜激情av网站| 国产亚洲精品第一综合不卡| 国产黄a三级三级三级人| 亚洲黑人精品在线| 狂野欧美激情性xxxx| 人人妻人人看人人澡| 12—13女人毛片做爰片一| 黄片小视频在线播放| 日本 欧美在线| 人妻久久中文字幕网| 91九色精品人成在线观看| 999精品在线视频| 1024香蕉在线观看| 色精品久久人妻99蜜桃| 婷婷丁香在线五月| 在线观看舔阴道视频| ponron亚洲| 久久99热这里只有精品18| 99国产精品一区二区蜜桃av| 国产亚洲欧美98| 国产成人精品久久二区二区免费| 国产一区二区三区视频了| 国产精品久久久久久精品电影 | 色综合婷婷激情| 国产视频一区二区在线看| 美女午夜性视频免费| 极品教师在线免费播放| 美女 人体艺术 gogo| 黄色毛片三级朝国网站| 亚洲精品在线美女| 久久精品国产清高在天天线| 亚洲黑人精品在线| 狂野欧美激情性xxxx| 看黄色毛片网站| 欧美一区二区精品小视频在线| 天堂动漫精品| 日韩高清综合在线| 在线播放国产精品三级| 女人爽到高潮嗷嗷叫在线视频| 伊人久久大香线蕉亚洲五| 变态另类丝袜制服| 亚洲成av人片免费观看| 观看免费一级毛片| 女同久久另类99精品国产91| 一进一出好大好爽视频| 两个人视频免费观看高清| e午夜精品久久久久久久| 亚洲午夜理论影院| 婷婷精品国产亚洲av在线| 黑丝袜美女国产一区| 久久这里只有精品19| 久久 成人 亚洲| 大型黄色视频在线免费观看| 天堂√8在线中文| 国产麻豆成人av免费视频| 国产精品电影一区二区三区| 大型黄色视频在线免费观看| 19禁男女啪啪无遮挡网站| 男人操女人黄网站| 精品一区二区三区av网在线观看| 国产视频一区二区在线看| 日本 av在线| 亚洲国产精品成人综合色| 色av中文字幕| 曰老女人黄片| 中文字幕人成人乱码亚洲影| cao死你这个sao货| 亚洲精品国产一区二区精华液| 老熟妇乱子伦视频在线观看| 日本 av在线| 亚洲精华国产精华精| 男女那种视频在线观看| 亚洲熟女毛片儿| 亚洲av美国av| 后天国语完整版免费观看| 淫妇啪啪啪对白视频| 亚洲人成网站在线播放欧美日韩| 亚洲第一青青草原| 精品欧美国产一区二区三| 怎么达到女性高潮| www.自偷自拍.com| 色综合站精品国产| 国产极品粉嫩免费观看在线| 99国产精品99久久久久| 欧美日本亚洲视频在线播放| 国产高清有码在线观看视频 | 中亚洲国语对白在线视频| 在线av久久热| 欧美绝顶高潮抽搐喷水| 精品久久久久久,| 国产成人精品久久二区二区91| 亚洲国产精品久久男人天堂| 精品少妇一区二区三区视频日本电影| 亚洲第一青青草原| 1024手机看黄色片| 国产精品电影一区二区三区| www.自偷自拍.com| 免费看美女性在线毛片视频| 欧美中文日本在线观看视频| 亚洲国产精品合色在线| 日韩欧美一区视频在线观看| 给我免费播放毛片高清在线观看| 91老司机精品| 一区二区三区高清视频在线| 亚洲av电影在线进入| 99热6这里只有精品| 好看av亚洲va欧美ⅴa在| 国产一区在线观看成人免费| 日本 欧美在线| 亚洲成av人片免费观看| 亚洲一区高清亚洲精品| 国产精品亚洲美女久久久| 国产野战对白在线观看| 午夜亚洲福利在线播放| 欧美丝袜亚洲另类 | 国产伦一二天堂av在线观看| 午夜免费成人在线视频| 午夜福利免费观看在线| 在线观看日韩欧美| netflix在线观看网站| 久久久久久人人人人人| 久久久水蜜桃国产精品网| 两性夫妻黄色片| 91大片在线观看| 国产午夜福利久久久久久| 国产欧美日韩精品亚洲av| 国产日本99.免费观看| 波多野结衣巨乳人妻| 精品人妻1区二区| 国产精品精品国产色婷婷| 国产日本99.免费观看| 无限看片的www在线观看| 欧美av亚洲av综合av国产av| 免费看美女性在线毛片视频| 国产成人欧美在线观看| 欧美日本视频| 久久久久亚洲av毛片大全| 日韩视频一区二区在线观看| 欧美色欧美亚洲另类二区| 久久精品亚洲精品国产色婷小说| 亚洲真实伦在线观看| 黄频高清免费视频| 欧美一级毛片孕妇| 超碰成人久久| 国产又爽黄色视频| 日本 欧美在线| 成人18禁在线播放| 亚洲 欧美 日韩 在线 免费| АⅤ资源中文在线天堂| 美女 人体艺术 gogo| 国产成+人综合+亚洲专区| 波多野结衣av一区二区av| 日本免费a在线| 久久久久国产一级毛片高清牌| 国产激情偷乱视频一区二区| 777久久人妻少妇嫩草av网站| 国产成人系列免费观看| 精品乱码久久久久久99久播| 俺也久久电影网| 成人精品一区二区免费| 男女午夜视频在线观看| 日韩欧美免费精品| 欧美一级毛片孕妇| 久热爱精品视频在线9| 国产精品久久久人人做人人爽| 视频在线观看一区二区三区| 天堂影院成人在线观看| 在线观看www视频免费| 99在线人妻在线中文字幕| 国产麻豆成人av免费视频| 老司机深夜福利视频在线观看| 巨乳人妻的诱惑在线观看| 69av精品久久久久久| xxx96com| 亚洲av成人不卡在线观看播放网| a级毛片a级免费在线| 久久精品国产99精品国产亚洲性色| 女同久久另类99精品国产91| 99在线人妻在线中文字幕| 亚洲一区高清亚洲精品| 变态另类成人亚洲欧美熟女| 国产精品二区激情视频| 亚洲成a人片在线一区二区| 国产精品99久久99久久久不卡| 国产一区二区在线av高清观看| 午夜福利成人在线免费观看| 精品国产亚洲在线| 一个人免费在线观看的高清视频| 女人被狂操c到高潮| 黄色成人免费大全| 亚洲中文av在线| 国产精品亚洲av一区麻豆| av免费在线观看网站| 国产精品日韩av在线免费观看| 十分钟在线观看高清视频www| 日韩中文字幕欧美一区二区| 日韩视频一区二区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 黑人巨大精品欧美一区二区mp4| 亚洲成av片中文字幕在线观看| 免费看日本二区| 亚洲成人免费电影在线观看| 黄色视频,在线免费观看| 国产av一区在线观看免费| 在线观看www视频免费| 日本免费a在线| 免费观看精品视频网站| 18禁黄网站禁片免费观看直播| tocl精华| 久久久久国产精品人妻aⅴ院| 99精品欧美一区二区三区四区| 成年免费大片在线观看| 不卡一级毛片| av视频在线观看入口| 免费一级毛片在线播放高清视频| 久热这里只有精品99| 精品无人区乱码1区二区| 在线观看66精品国产| 久久久久久国产a免费观看| 99国产精品99久久久久| 国产精品综合久久久久久久免费| 99国产精品一区二区蜜桃av| 国产精品野战在线观看| 久久亚洲精品不卡| 三级毛片av免费| 亚洲 欧美 日韩 在线 免费| 日本免费a在线| 日韩大尺度精品在线看网址| 亚洲片人在线观看| 无人区码免费观看不卡| 国产精品一区二区精品视频观看| 女性生殖器流出的白浆| svipshipincom国产片| 午夜福利在线在线| 久久久久九九精品影院| 亚洲一区中文字幕在线| av有码第一页| svipshipincom国产片| 国产三级在线视频| 一边摸一边抽搐一进一小说| 不卡av一区二区三区| 成人免费观看视频高清| 久久国产精品人妻蜜桃| 欧美av亚洲av综合av国产av| 国产真实乱freesex| 一二三四在线观看免费中文在| 免费高清在线观看日韩| 国产精品一区二区免费欧美| 搡老岳熟女国产| 亚洲精品一区av在线观看| 亚洲黑人精品在线| 精品卡一卡二卡四卡免费| 99精品欧美一区二区三区四区| www.熟女人妻精品国产| 亚洲成人精品中文字幕电影| 18禁黄网站禁片午夜丰满| 国产主播在线观看一区二区| 久久久水蜜桃国产精品网| 一本精品99久久精品77| 国产又爽黄色视频| 非洲黑人性xxxx精品又粗又长| 亚洲午夜理论影院| 亚洲精品粉嫩美女一区| 午夜久久久在线观看| 人人妻人人澡欧美一区二区| 精品欧美国产一区二区三| 侵犯人妻中文字幕一二三四区| 国产伦人伦偷精品视频| 中文字幕av电影在线播放| 久久中文看片网| 国产亚洲欧美精品永久| 18禁观看日本| 18禁裸乳无遮挡免费网站照片 | 欧美日韩瑟瑟在线播放| 琪琪午夜伦伦电影理论片6080| 我的亚洲天堂| 午夜免费观看网址| 校园春色视频在线观看| 国产成人系列免费观看| 国内精品久久久久久久电影| 999久久久精品免费观看国产| 精品久久久久久久末码| 在线十欧美十亚洲十日本专区| 亚洲av日韩精品久久久久久密| 91老司机精品| 亚洲精品久久成人aⅴ小说| 国产区一区二久久| 亚洲国产高清在线一区二区三 | 青草久久国产| 色婷婷久久久亚洲欧美| 他把我摸到了高潮在线观看| 免费在线观看亚洲国产| ponron亚洲| 色播亚洲综合网| 国产1区2区3区精品| 黄色视频不卡| 我的亚洲天堂| 国产主播在线观看一区二区| 国产v大片淫在线免费观看| 午夜福利在线在线| 亚洲人成电影免费在线| 国产一区二区三区视频了| 十八禁人妻一区二区| 一级a爱视频在线免费观看| 欧美成人性av电影在线观看| 成人永久免费在线观看视频| 久久久久久大精品| 成人一区二区视频在线观看| 变态另类丝袜制服| 日韩大尺度精品在线看网址| 亚洲第一青青草原| 高清在线国产一区| 日韩欧美一区视频在线观看| 国产成人欧美在线观看| av天堂在线播放| 很黄的视频免费| 欧美+亚洲+日韩+国产| 午夜老司机福利片| 日韩欧美 国产精品| 久99久视频精品免费| 欧美+亚洲+日韩+国产| 三级毛片av免费| 午夜福利高清视频| 成人手机av| 一区二区三区激情视频| 亚洲五月婷婷丁香| 亚洲中文字幕日韩| 久久九九热精品免费| 久久久久久免费高清国产稀缺| 久久久久九九精品影院| 欧美激情久久久久久爽电影| 香蕉丝袜av| av电影中文网址| 香蕉丝袜av| 搞女人的毛片| 国产成人精品无人区| 超碰成人久久| 最近最新中文字幕大全免费视频| 无限看片的www在线观看| 欧美日韩亚洲综合一区二区三区_| av电影中文网址| 狠狠狠狠99中文字幕| 又黄又爽又免费观看的视频| 大型av网站在线播放| 亚洲av五月六月丁香网| 亚洲国产欧美一区二区综合| 欧美精品亚洲一区二区| 99久久无色码亚洲精品果冻| 好看av亚洲va欧美ⅴa在| 国内毛片毛片毛片毛片毛片| 在线观看舔阴道视频| 制服丝袜大香蕉在线| 欧美绝顶高潮抽搐喷水| 国产黄色小视频在线观看| 天天一区二区日本电影三级| 51午夜福利影视在线观看| 欧美精品啪啪一区二区三区| 亚洲黑人精品在线| 一二三四在线观看免费中文在| 又大又爽又粗| 欧美三级亚洲精品| 高清毛片免费观看视频网站| 久久精品91蜜桃| 在线观看免费午夜福利视频| 91字幕亚洲| 精品无人区乱码1区二区| 自线自在国产av| 女人被狂操c到高潮| av欧美777| 久久人妻av系列| 亚洲av美国av| 久久国产亚洲av麻豆专区| 90打野战视频偷拍视频| 一本大道久久a久久精品| 日韩av在线大香蕉| av免费在线观看网站| 国产精华一区二区三区| 黄色 视频免费看| 国产成人一区二区三区免费视频网站| 国产亚洲精品第一综合不卡| 91麻豆av在线| 他把我摸到了高潮在线观看| 亚洲人成网站高清观看| 国产av不卡久久| 亚洲av片天天在线观看| 99国产综合亚洲精品| 天天一区二区日本电影三级| 日本 欧美在线| 妹子高潮喷水视频| 女警被强在线播放| 男女午夜视频在线观看| 亚洲av电影不卡..在线观看| 国产免费av片在线观看野外av| 久久久久久久久免费视频了| 视频区欧美日本亚洲| 欧美绝顶高潮抽搐喷水| 国产精品 国内视频| 不卡一级毛片| 国产国语露脸激情在线看| 欧美又色又爽又黄视频| 成熟少妇高潮喷水视频| 日韩欧美一区视频在线观看| 国产精品一区二区免费欧美| 变态另类丝袜制服| av免费在线观看网站| 无遮挡黄片免费观看| 欧美中文综合在线视频| 国产私拍福利视频在线观看| 脱女人内裤的视频| 免费高清视频大片| 午夜日韩欧美国产|