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

    Beach Morphology and Coastline Evolution in the Southern Bohai Strait

    2015-03-15 01:43:44ZHANGWeiWUJianzhengLIWeiranZHULonghaiHURijunJIANGshenghuiSUNYonggenandWANGHuijuan
    Journal of Ocean University of China 2015年5期

    ZHANG Wei, WU Jianzheng, LI Weiran, ZHU Longhai, HU Rijun,, JIANG shenghui, SUN Yonggen, and WANG Huijuan

    1) College of Marine Geoscience, Ocean University of China, Qingdao 266100, P. R. China

    2) Qingdao Marine Geological Engineering Investigation Institute, Qingdao 266071, P. R. China

    3) Key Laboratory of Submarine Geosciences and Prospecting Techniques, Ministry of Education, Ocean University of China, Qingdao 266100, P. R. China

    4) First Institute of Oceanography, SOA, Qingdao 266061, P. R. China

    5) Qingdao Boyan Marine Environment Science & Technology Co. Ltd., Qingdao 266101, P. R. China

    Beach Morphology and Coastline Evolution in the Southern Bohai Strait

    ZHANG Wei1),2), WU Jianzheng1),3), LI Weiran1),3), ZHU Longhai1),3), HU Rijun1),3),*, JIANG shenghui1),3), SUN Yonggen4), and WANG Huijuan5)

    1) College of Marine Geoscience, Ocean University of China, Qingdao 266100, P. R. China

    2) Qingdao Marine Geological Engineering Investigation Institute, Qingdao 266071, P. R. China

    3) Key Laboratory of Submarine Geosciences and Prospecting Techniques, Ministry of Education, Ocean University of China, Qingdao 266100, P. R. China

    4) First Institute of Oceanography, SOA, Qingdao 266061, P. R. China

    5) Qingdao Boyan Marine Environment Science & Technology Co. Ltd., Qingdao 266101, P. R. China

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

    The beach studied in this paper spans a length of 51 km and is one of several long sandy beaches in the southern Bohai Strait. Due to the obstruction of islands in the northeast and the influence of the underwater topography, the wave environment in the offshore area is complex; beach types and sediment transport characteristics vary along different coasts. The coastlines extracted from six aerial photographs in different years were compared to demonstrate the evolving features. Seven typical beach profiles were selected to study the lateral beach variation characteristics. Continuous wind and wave observation data from Beihuangcheng ocean station during 2009 were employed for the hindcast of the local wave environment using a regional spectral wave model. Then the results of the wave hindcast were incorporated into the LITDRIFT model to compute the sediment transport rates and directions along the coasts and analyze the longshore sand movement. The results show that the coastline evolution of sand beaches in the southern Bohai Strait has spatial and temporal variations and the coast can be divided into four typical regions. Region (I), the north coast of Qimudao, is a slightly eroded and dissipative beach with a large sediment transport rate; Region (II), the southwest coast of Gangluan Port, is a slightly deposited and dissipative beach with moderate sediment transport rate; Region (III), in the central area, is a beach that is gradually transformed from a slightly eroded dissipative beach to a moderately or slightly strong eroded bar-trough beach from west to east with a relatively moderate sediment transport rate. Region (IV), on the east coast, is a strongly eroded and reflective beach with a weak sediment transport rate. The wave conditions exhibit an increasing trend from west to east in the offshore area. The distribution of the wave-induced current inside the wave breaking region and the littoral sediment transport in the nearshore region exhibit a gradual weakening tendency from west to east, which is opposite to the trend of the wave conditions outside the breaking region. The presence of submerged shoal (Dengzhou Shoal), deep trough (Dengzhou Channel), islands and irregular topography influnces the wave climate, beach types, wave-induced current features, littoral sediment transport trends and coastline evolution patterns in the southern Bohai Strait. Human activities, such as the sand exploitation of Dengzhou Shoal and other coastal engineering projects, also influence the beach morphology and coastline evolution.

    sand beach; coastline evolution; profile feature; wave hind-casting; wave-induced current; littoral sediment transport

    1 Introduction

    Sand beaches are charming and offer beautiful sceneries and hydrophilic entertainment; they are also an important coastal zone for human activities. Domestic and foreign scholars have conducted extensive studies on sand beaches around the world (Pilkey and Richter, 1964; Hallermeier, 1980–1981; Wright and Short, 1984; Dean, 1991; Benedet and List, 2008). At present, domestic research mainly focuses on beaches in the South China Sea(Dai et al., 2007; Cai et al., 2009; Qi et al., 2010) and East China Sea (Chen et al., 2011; Qu et al., 2013). In addition, numerous studies (Kuang et al., 2010; Zhang et al., 2012) were conducted on the beach of Beidaihe, which is located on the west side of the Bohai Sea, China. The coastline of Shandong Province is the second longest in China, but research on sand beaches is lacking. With the increase in coastal development activities, the research and protection of sand beaches are becoming more and more important. Previous studies of sand beach morphology focused on investigations of the profile changes of the cross shore (Masselink and Pattiaratchi, 2001; Jeanson et al., 2013) and the coastline evolution along the beach (Schwarzer et al., 2003; Robinson, 2004). Com-monly used research techniques include field surveys (Rao et al., 2009), remote sensing image interpretation (White and Asmar, 1999; Liu et al., 2013), numerical simulation (Zheng and Dean, 1997; Takaaki et al., 2010), and theoretical analysis (Chiara et al., 2011; van Rijn, 2011).

    In the present work, the study area in the southern Bohai Strait is characterized by a complex morphology (submerged shoal, a deep trough, and islands), coastal structures and irregular depth contours, as well as varied wave conditions and sand beach evolution. The relevant previous studies mainly focused on the sediment characteristics (Dong et al., 2011), hydrodynamics (Jiang and Gao, 2002), geomorphological features (Wu, 1992), and wave conditions (Chen and Hu, 1992) around Bohai Strait and Dengzhou Shoal. However, scholars only sparsely focused on the sandy beach in the southern Bohai Strait. Based on a field survey, Luan (2011) described the coastal sediment and geomorphological characteristics of the beach east of Luanjiakou Port and indicated that the beach is eroded and straightened. The evolution of Dengzhou Shoal and its association with the erosion on the west coast of Penglai has also attracted considerable attention (Wang, 1994; Wu et al., 1997; Li et al., 2004; Wang et al., 2005). These studies suggested that the sand exploitation of Dengzhou Shoal was the main reason underlying the coastal erosion along Penglai. However, these studies were mostly conducted before 2005 without considering the beach between Luanjiakou Port and Qimudao. Due to the complex geomorphological characteristics in the study area, the sediment dynamics and beach evolution process, which are useful for coastal protection and development, remain to be understood.

    In this study, we attempt to elucidate the evolution of a sand beach and the controlling factors. More specifically, the research aims to 1) use satellite images to analyze the coastline variation, 2) compare in-situ measurements with historical profiles to study the lateral variation of the beach, 3) make use of the SW wave model to study the wave conditions in offshore areas, 4) apply the coastal sediment transport model (the LITDRIFT model) to simulate the wave-induced current and littoral sediment movement in nearshore areas.

    2 Study Area

    The study area is located in the southern Bohai Strait from the west side of Qimudao to the east side of Penglai Port (Fig.1). The northern part of the study area includes the Miaodao Archipelago and Daheishan Island, and the central part covers Sangdao Island. Freshwater inflow from the Huangshui River the east of Gangluan Port has a small discharge. The coastline of the beach runs through the study area with an overall direction of ENE–WSW and extends approximately 51 km. Unlike a normal straight coast, the bathymetry in the area is highly complex due to the presence of submerged shoal (Dengzhou Shoal), a deep trough (Dengzhou Channel), islands and coastal structures (Fig.1). The minimum water depth over the shoal is less than 5 m but the depth increases eastwards to approximately 47 m in the trough. In the study area, Dengzhou Shoal extends approximately 13 km northwestward from the Penglai Port. Along the coast, there are many ports and scenic spots, constructed with groins, breakwaters, and revetments to prevent coastal erosion. The characteristic sand beach width is mostly less than 70 m, although it is 140 m at the groyne root. The tide is small in the nearshore area with an average rage of 1.06 m.

    Fig.1 The study area with bathymetry contours and the four main regions with seven typical beach profiles.

    3 Data and Methods

    3.1 Basic Data

    According to the distribution of the beaches and the data from a field investigation, the study area was divided into four main regions with seven typical profiles, which were measured in the shore-normal direction (Fig.1). By contrasting the depth data in the same positions of the profiles obtained from the field investigation with a large-scale nautical chart, the evolution of the selected beach profiles were characterized (Table 1). The field investigation used a Leica SR 500 Real Time Kinematic Global Positioning System (RTKGPS) with a position accuracy of ±1 cm and an elevation accuracy of ±2 cm as well as an SDE-28S Single Frequency Depth Sounder with a depth accuracy of ±1 cm ±0.1%D. A base reference station was established at Luanjiakou Port. The vertical coordinate was referenced to the mean sea level (MSL) for accurate comparisons. The nautical chart data are processed by the natural neighbor interpolation method (Yvonnet et al., 2006; Dinis et al., 2010).

    To obtain the grain size parameter for sediment transport modeling, the sediment samples collected from all of the beach profiles were subjected to sieve analysis (Fernlund, 1998). For wave hindcast, the required field data for the boundary conditions of the regional spectral wave model were obtained from in situ wave recorders. SBA-2 Optical Wave Recorders were deployed at the Beihuangcheng ocean station (Fig.2) and the wave parameters were recorded at three-hour intervals from 1 January to 31 December 2009 (Fig.3). The wave parameters during a high wave event from 9 to 11 November 2009 at the Longkou ocean station were used to calibrate the modeling result.

    3.2 Coastline Evolution and Comparison

    For the multi-temporal analysis, we used six aerial photographs (Table 2) from 1979 to 2013 which were downloaded for free from http://glovis.usgs.gov/. Geometric correction was used in these Landsat images with a WGS-84 coordinate system and a Gauss-Kruger projection (with a RMS less than 0.5 pixels). Path radiance correction was used in dark object subtraction (DOS) (Chavez Jr., 1988). The deviation of the tidal level was ignored due to the small tides and steep slopes of the beaches. The coastlines were obtained by the single band density slicing method based on the NIR band (Liu et al., 2013).

    Table 1 Location and data date of the beach profiles in the study area

    Fig.2 Nested type of flexible mesh used in the regional spectral wave model for the study area.

    Fig.3 Time series of offshore significant wave height (Hs), peak period (Tp) and peak wave direction (θ) at a 10 m depth from January 2009 to December 2009.

    Table 2 Acquisition time, resolutions and the local tidal levels of the aerial images

    3.3 Wave Hindcast

    Wave hindcast for the study area was conducted by the MIKE 21 Spectral Wave (MIKE 21 SW) model (Byrnes and Hammer, 2004; DHI, 2005). This model can simulate the growth, decay and transformation of wind-generated waves and swell in offshore and coastal areas by dissipation due to white-capping, bottom friction and depth-i n-duced wave breaking, and by refraction and shoaling due to depth variation, wave interaction and the effect of time-varying water depth (Sabique et al., 2011). All forcing sources except tidal variation are considered in the model. A fully spectral formulation and instationary time formulation are used (Komen et al., 1994). Diffraction is included using the phase-decoupled refraction-diffraction approximation proposed by Holthuijsen et al. (2003) with a smoothing factor of α=1. Wave breaking is based on the formulation of Battjes and Janssen (1978) where α=1 and γ=0.73 (SWAN team, 2012). The Nikuradse roughness coefficient is used for the dissipation due to bottom friction, where KN=0.04. The dissipation due to white capping based on the theory of Hasselmann (1974) is found to be 4.5.

    The wave simulation area spanned the region from western Qimudao to eastern Zhifudao Island, and the northern side extended to northernmost Bohai Strait island, i.e., Beihuangcheng Island (37?32′ to 38?25′N and 120?13′ to 121?25′E). The grid distribution is shown in Fig.2. The total simulated area is 8640.64 km2, which totals 26093 grids; the largest grid size is 0.90 km2. The area with concentrated grids (minimum size of 59.91 m2) was the region of our special concernment; the coastal zone extends approximately 63 km from Qimudao to the east of Penglai. Water depth data were obtained from the latest available nautical chart, and nearshore high-r esolution bathymetric data obtained from a field investigation.Wave heights, wave periods and wave directions were obtained from the numerical simulation results. The calibration results at Longkou station exhibited a strong correlation between the simulation results and the measured parameters (Fig.4). The annual wave distribution patternand the corresponding rose plots at the 10 m water depth of the in-situ measurement profiles are shown in Fig.9. The hindcast data on wave heights, wave periods and wave directions were utilized as input for the LITDRIFT modeling of sediment transport (DHI, 2005).

    Fig.4 Calibration of spectral wave model results for a) significant wave heights b) wave directions and c) wave periods at Longkou station (located on the 16.8 m depth contour).

    3.4 LITDRIFT Model

    The Annual Sediment Drift and Longshore Current model of LITDRIFT has been applied to calculate the longshore transport of non-cohesive sediment and waveinduced currents on a long uniform beach (DHI, 2009). LITDRIFT provides a detailed deterministic description of longshore sediment transport (LST) for an arbitrary bathymetry and calculates the net and gross littoral transport amounts for a section of coastline over a specific design period considering wave breaking. The input data used for the LITDRIFT simulation are the cross-shore profiles, the coast-normal angle, wave heights, wave directions and wave periods at a reference water depth of 10 m. For all profiles, the profiles’ length used to calculate the LST amount was 650 m with maximum boundary depths from 8.3 m to 10.0 m. The simulation results of a wave-induced current inside a breaking wave area of a 2.0 m depth are described in Section 4.4. The boundary wave parameters of the profiles are derived from wave hindcast in Section 3.3. The model considers Battjes and Janssen’s approach of wave propagation from deep water (Battjes and Janssen, 1978).

    4 Results and Discussion

    4.1 Coastline Evolution

    The changes of the coastlines since 1979 as extracted from the satellite images are shown in Fig.5 and Fig.6. The coastline of Region (I) is flat and extends from east to west. The coast was continuously eroded; the maximum erosion distance was 60.6 m (2.75 m yr-1) from 1979 to 2000. Then, the coast was in a dynamic equilibrium state with negligible variation between 2000 and 2009. With the construction of breakwaters and other influencing factors, the coast exhibited an erosional trend again in which the maximum erosion distance was 52.8 m (5.28 m yr-1). An accretion trend with a nearly equivalent distance can be observed on both sides of the eastern breakwater root.

    Fig.5 Each coastline was extracted from the aerial photo only by visual analysis. These are examples of extracted coastlines in 1979 (a), 1984 (b), 1989 (c), 2000 (d), 2009 (e), and 2013 (h).

    Fig.6 Coastline changes at various coasts since 1979.

    Overall, Region (II) presents a coastline extending from northeast to southwest. A continuous natural sand beach existed before 2000 that showed a gentle erosion trend at an average speed of approximately 0.45 m yr-1. Currently, many artificial buildings such as groynes, breakwaters, scenic spots, reclamation engineering and revetment have been constructed along the coast. The coastline distribution pattern seems to be significantly influenced by these buildings. For example, the presence of the east port led to the disappearance of the original beach coastline after 2009 (Fig.6). However, the establishment of a tourism landscape, built on the east side of the region, formed a typical artificial static equilibrium headland bay (Moon Bay) (Hsu et al., 1989; González and Medina, 2001) and sand beach. The beach at the root of the northeast groyne shows significant accretion, but erosion occurred along the coastline southwest of Gangluan Port, which implies a longshore sediment transport from northeast to southwest. The natural sand beach on the west estuary of the Yongwen River disappeared due to the construction of the revetment after 2009.

    The coastlines of Region (I) and Region (III) have similar distributions. Three elliptic reclamation engineering projects and Gangluan Port are located on the west coast. The sand beach that extends along the coast to the west of Luanjiakou Port showed continuous erosion, with a maximum erosion distance of 77 m (3.5 m yr-1) from 1979 to 2000. Then, the erosion speed slowed with a maximum erosion distance of 36 m (2.57 m yr-1) between 2000 and 2013. The reclamation engineering in the west remarkably modified the coastline, which caused serious erosion between the engineering projects and deposition in the shadow zones.

    The regional coastline distribution pattern of Region (IV) is similar to that of Region (II): it extends from northeast to southwest. However, the beach of Region (IV) is short and narrow, and its distribution is not continuous. Moreover, there are also many eroded cliffs and houses and roads on the verge are swallowed (Fig.7). Additionally, this coastal segment undergoes changes because of a large amount of coastal engineering, such as breakwaters, groynes, ponds and revetments (Fig.6). Previous studies have been conducted on Dengzhou Shoal in the offshore area and suggest that the sand exploitation at the end of the last century severely damaged the shoal, which led to serious coastline erosion (Li et al., 2004; Wang et al., 2005). Most of these coastal constructions were established during the last decade to prevent serious coastal erosion. The beaches scattered around the breakwaters and groynes are mainly present at a small scale with primary sand types of coarse sand to gravel. The comparison in Fig.6 shows a dynamic variation between erosion and accretion from 1979 to 1984. This trend has been confirmed to relate to the existence of Dengzhou Shoal, which can weaken the wave influence significantly. However, the shoal was generally eroded, with an average erosion distance of 30 m (5.0 m yr-1) from 1984 to 1989. This result correlates well with that of Wu et al. (1992, 1997). Since 2009, the coastline has changed insignificantly, indicating that these coastal engineering constructions have played a certain role in preventing coastal erosion.

    Overall, the coastline evolution of sand beaches in the southern Bohai Strait exhibit spatial and temporal variations. Prior to the year 2000, natural sand beaches experienced serious coastline erosion in Regions (I) and (III) at a maximum speed of 3.5 m yr-1. The coastline erosion speed in Region (III), which is located in the middle of the region, was greater than that in Region (I) located in the west. After 2000, the speeds decreased in the two re-gions. The coastal evolution in Region (II), which showed gentle erosion previously, changed due to the coastal engineering after 2000. The coast exhibited a dynamic variation with the existence of the Dengzhou Shoal which weakened the wave influence significantly before 1984. Then, the coast was seriously eroded, with a maximum speed of 5.0 m yr-1between 1984 and 1989, due to the crude sand exploitation of Dengzhou Shoal. Since 2009, the coastline has changed insignificantly due to the coastal protection projects.

    Fig.7 Field photographs showing severe erosion of coastline in Region (IV).

    4.2 Profile Features

    Based on the hydrodynamic processes and the relative contributions of different mechanisms to sediment transport and morphologic change, Wright and Short (1984) classified sand beach types, such as dissipative, reflective and several intermediate states (e.g., bar-trough, rhythmic bar-trough, transverse bar-rip and ridge and runnel). The temporal morphological variability of the specific crossshore transects are shown in Fig.8, where the upper figures indicate the profile obtained in different years, and the lower figures show the difference between temporally varying profiles.

    Comparisons revealed that profile 1 of Region (I) (Fig.8) was in a slightly eroded state and exhibited very low mobility between 2007 and 2012. Two or three gently sloping sandbars formed by the storm surge were located below the water in Profile 1 and Profile 2. The beach of Region (I) can be classified as a slightly eroded, dissipative beach (Ojeda et al., 2011; Thiebot et al., 2012).

    Profile 3 in Region (II) is characterized by a deposition trend (Fig.8) from 1975 to 2012, with a maximum deposition thickness of approximately 2.8 m (0.56 m yr-1). The occurrence of different trends in Region (I) and Region (II) may be due to the constructed groynes along the coast. These groynes may cause heterogeneous wave circumstances and thus block the sediment transport and enhance the shoreward sediment accumulation process.

    As observed in Region (III), the profile changes were similar for Profile 4 and Profile 5 from 1975 to 2004; both of the profiles were characterized by erosion with a maximum value of approximately 1 m (0.03 m yr-1) (Fig.8). However, the profiles have showed different changes after 2004. The beach having Profile 4 therein with a gentle slope and a sand bar formed below 3 m in a storm surge, is a micro-erosion dissipative type. For Profile 5, a sand bar also existed below 5 m, and a 150 m wide trough with a relatively steeper slope is present. Many smallsized beach cusps exist along the coastline. All of these features suggest that the beach under Profile 5 is a bar- trough beach, which is the intermediate state between a dissipative and reflective beach. The maximum erosion value above a 4.3 m water depth in Region (III) is approximately 2.1 m (0.3 m yr-1) after 2004 but less than 1 m (0.14 m yr-1) under the same water depth before 2004. In general, the sand beach has been transformed from a slightly eroded dissipative beach to a moderately or somewhat strong eroded bar-trough beach from west to east in Region (III), both of which are intermediate states between dissipative and reflective ones.

    The changes of Profiles 6 and 7 in Region (IV) indicated strongly erosion landward and accretion seaward at 6 m to 7 m depths, the maximum erosion value being 3.5 m (0.11 m yr-1). There is a strong wave energy environment in Region (IV) where the waves clearly influence the depth more than in other regions. During 2004 to2012, these profiles were eroded landward and accreted seaward at 8 m to 8.5 m depths with a maximum erosion value of 1.6 m (0.18 m yr-1) and an increasing wave influence depth. All of these features indicate that the beachin Region (IV) is a strongly eroded reflective beach.

    Fig.8 Variations of the selected profiles along coast of the study area.

    4.3 Wave Climate

    The wave simulation results indicate that there are distinct regional characteristics over different profiles in our study area (Fig.9 and Table 3). Most observed waves in Region (I) are in the directions of N to NNE, then in WNW and NNW, scarcely with the directions of W and E. Region (II) is affected by the northeastern Sangdao Island, and the frequencies of NW to W waves are significantly higher than those of the NNW to ENE waves; the highest frequency waves are in the W direction. In Region (III), the most frequent wave directions are NNW and N. The wave frequency in the W direction is relatively high as well, it being impacted by the southern waterway on the west side of Sangdao Island. In contrast, due to the blocking of the Miaodao Archipelago, the NE and ENE waves rarely occur. Also, due to the Miaodao Archipelago, N and NNW waves occur less often, and W waves are frequently observed in Region (IV).

    Wave height statistics suggests that the frequency of wave heights greater than 1.8 m is higher in Regions (III) and (IV) on the eastern side of the study area, than those in Regions (I) and (II). An eastward trend increasingly occurs, corresponding to a wave climate of a dissipative or reflective beach from west to east (Table 4).

    Fig.9 Rose plots for wave heights and directions for offshore locations at a 10 m depth off the coast.

    Table 3 Frequency statistics of waves for different stations and directions

    Table 4 Frequency statistics of waves of various heights and stations

    4.4 Wave-Induced Current

    The wave-induced currents generated in coastal regions are directly responsible for sediment transport and morphology evolutions (Xie, 2012). Based on the wave simulation results, the wave-induced current inside the breaking region of a 2.0 meter depth in the selected seven profiles can be calculated by the LITDRIFT model (Fig.10). Similar variation tendencies and strong currents occur in C1 and C2 of Region (I). The max-speed was higher than 4.0 m s-1. The current speed in C3 of Region (II) is negligible due to blocking by the breakwaters around it. The current curves of C4 and C5 in Region (III) are similar in that the maximum values are greater than 3.0 m s-1(with opposite coastal sloping directions). The current at the east station (C5) is stronger than that in the west. The similar variation tendencies in C6 and C7 of Region (IV) can also be seen in Fig.10. However, the max current speeds of 2.3 m s-1in C6 and 1.4 m s-1in C7 are obviously smaller than those in other regions. Meanwhile, the occurrence frequency of wave-induced current in C7 is higher than that at other points. In all, the distribution of the wave-induced currents inside the wave breaking region exhibits a gradual weakening tendency from west to east, which is opposite to the trend of wave conditions outside the breaking region. The non-correspondence is mainly caused by the complex morphology. In the west region, the morphology is flat, such that waves break closeto the coast with strong energy to generate a strong longshore current. In the east region, the morphology is complex and includes shoals, which can wake the waves break in the offshore area, weaken the energy and generate a weak longshore current.

    Fig.10 Time series plots of mode-simulated wave-induced current for nearshore locations (C1 to C7).

    4.5 Littoral Sediment Transport

    The littoral profile sediment transport amount is shown in Table 5. Strong transports take place at Profiles 2 and 5 with large net and gross LST amounts. Moderate transports can be observed at Profiles 1, 3 and 7; meanwhile, weaker transports occur at Profiles 2 and 6. Comparing with the wave-induced current conditions in Fig.10 shows an incomplete coincidence trend. The strong longshore current at C1 and moderate LST amount at Profile 1 can be ascribed to the deep water depth, which causes weak turbulence and transportation of sediment. Otherwise, the weak longshore current at C3 and moderate LST amount at Profile 3 are due to the low water depth, which enhances the turbulence and the transportation of the sediment. Thus, topography is also an important factor for the longshore sediment transport.

    The LST directions off different coasts also exhibit a regional character. The LST direction is from east to west at Profile 1 and from northeast to southwest at Profile 3; it from west to east in Profile 2 at which the wave condition is influenced by the Sangdao Island located in the northeast. Affected by the Miaodao Archipelago and Dengzhou Shoal locating in the northwest, the LST directions at Profiles 4 and 5 are from west to east. Meanwhile, the LST directions at Profiles 6 and 7 are from northeast to southwest due to the coast span direction and the presence of Dengzhou Channal.

    Comprehensively considering the wave conditions, wave-induced current and littoral sediment transport features, we can identify a very interesting pattern, in which the wave height is enhanced in the offshore area with a 10 m depth along the coast from west to east; however, opposite trends in the wave-induced current and littoral sediment transport are observed in the nearshore area. The enhancement of wave height in the east is mainly due to the deep water in the Dengzhou Channel. The weakening trends of the wave-induced current and littoral sediment transport in the nearshore area are mainly attributed to Dengzhou Shoal, which can weaken the wave energy in the off-shore area.

    Table 5 Grain size (D50), shore slope and simulated amount of annual sediment transport at the beach profiles

    Fig.11 Longshore and cross-shore transport directions in the study area obtained from the LITDRIFT model.

    5 Conclusions

    The coastline in the study area, which is affected by sand exploitation and complicated geomorphology (i.e., channels, shoals, islands, and ports), exhibits wide variations in its evolution, wave conditions, wave-induced current and longshore sediment transport patterns.

    The coastline evolution of sand beaches in the southern Bohai Strait shows spatial and temporal variations. In the natural state, the west coast shows a serious erosion trend, while the coast in the eastern area that is shielded by Dengzhou Shoal exhibits a dynamic variation. Due to the sand exploitation of Dengzhou Shoal after 1984 and the coastal engineering after 2000, the coastline variations exhibits an obvious change: the erosion speed slows down in the west but becomes severe in the east. Currently, the coastline in the east is changing insignificantly due to the existence of the coastal protection projects.

    The beaches eroded along the coastline, except in Region (II) where the sediment transport and wave field were affected by the groynes. The beaches can be divided into several types based on their profile morphologies and comparisons. The beach in Region (I), which is located in the west of the study area, is a slightly eroded, dissipative beach. A slightly deposited, dissipative beach exists in Region (II) southwest of Gangluan Port. Region (III) in the central area is a transition region in which the beach is transformed from a slightly eroded dissipative one into a moderately or somewhat strongly eroded bar-trough beach which is an intermediate state between dissipative and reflective beaches from west to east. The beach in Region (IV) located in the eastern area is a strongly eroded, reflective beach.

    The various evolving coastline features and types reflect different dynamic environments. The wave hindcast results suggest that the frequency of strong waves has been obviously higher in Regions (III) and (IV) in the east than that in Regions (I) and (II). This trend corresponds to the process for dissipative beaches with low wave energy to be transformed into reflective beaches with high wave energy from west to east in the study area.

    The distribution of wave-induced currents inside the wave breaking region and the littoral sediment transport in the nearshore region displayed a gradual weakening tendency from west to east, which is opposite to the trend of the wave condition outside the breaking region.

    With the aid of aerial photography comparisons, field investigation, wave handcast, wave-induced current and littoral sediment transport simulation, this study determines the evolving features of sand beaches in the southern Bohai Strait and therefore provides useful guidance for coastal exploitation, beach regulation and beach protection.

    Acknowledgements

    The authors wish to express sincere thanks to Drs. Yongchen Xu and Peng Wang, and Mrs Chao Dong, Dongxiao Yin, and Xianfeng Zhang at the Ocean University of China for their hard work in the field investigation. This work is supported by the National Natural Science Foundation for the Youth (No. 41106039).

    Battjes, J. A., and Janssen, J. P., 1978. Energy loss and setup due to breaking of random waves. In: Proceedings of the 16thInternational Conference on Coastal Engineering. ASCE, New York, 569-587.

    Benedet, L., and List, J. H., 2008. Evaluation of the physical process controlling beach changes adjacent to nearshore dredge pits. Coastal Engineering, 55: 1224-1236.

    Byrnes, M. R., and Hammer, R. M., 2004. Collection of environmental data within sand resource areas offshore North Carolina and the environmental implications of sand removal for coastal and beach restoration. US Department of the Interior, Minerals Management Science, Leasing Division, Sand and Gravel Unit, Herndon, VA: Final Report, MMS Study 2000-056, 75-106.

    Cai, F., Su, X., Liu, J., Li, B., and Lei, G., 2009. Coastal erosion in China under the condition of global climate change and measures for its prevention. Progress in Natural Science, 19 (4): 415-426.

    Chiara, F. S., Massimo, B., and Marco, F., 2011. Application of the parabolic bay shape equation to sand and gravel beaches on Mediterranean coasts. Coastal Engineering, 59: 57-63.

    Chavez Jr., P. S., 1988. An improved dark-object subtraction technique for atmospheric scattering correction of multispectral data. Remote Sensing of Environment, 24: 459-479.

    Chen, B., Chen, Z. Y., Stephenson, W., and Finlayson, B., 2011. Morphodynamics of a boulder beach, Putuo Island, SE China coast: The role of storms and typhoon. Marine Geology, 283 (1-4): 106-115.

    Chen, X. Y., and Hu, Z. J., 1992. On waves and coastal erosion in the sea area near West Village, Penglai county, Shandong province. Journal of Oceanography of Huanghai & Bohai Seas, 10 (1): 19-26 (in Chinese with English abstract).

    Dai, Z. J., Du, J. Z., Li, C. C., and Chen, Z. S., 2007. The configuration of equilibrium beach profile in South China. Geomorphology, 86: 441-454.

    Dean, R. G., 1991. Equilibrium beach profiles characteristics and applications. Journal of Coastal Research, 7 (1): 53-84.

    DHI (Danish Hydraulics Institute), 2005. DHI Software. DHI Water and Environment, Agern Alle 5, DK-2970, Harshoalm, Denmark.

    DHI (Danish Hydraulics Institute), 2009. LITDRIFT User’s Guide: Longshore Current and Littoral Drift. DHI Software, 108pp.

    Dinis, L., Jorge, R., and Belinha, J., 2010. A 3D shell-like approach using a natural neighbour meshless method: Isotropic and orthotropic thin structures. Composite Structures, 92 (5): 1132-1142.

    Dong, C., Wu, J. Z., Zhu, L. H., Hu, R. J., and Zhang, W., 2011. Distribution pattern of surface sediments in Dengzhou Shoal area. Marine Geology Frontiers, 27 (8): 8-15 (in Chinese with English abstract).

    Fernlund, J. M. R., 1998. The effect of particle form on sieve analysis: A test by image analysis. Engineering Geology, 50 (1-2): 111-124.

    González, M., and Medina, R., 2001. On the application of static equilibrium bay formulations to natural and man-madebeaches. Coastal Engineering, 43 (3-4): 209-225.

    Hallermeier, R. J., 1980-1981. A profile zonation for seasonal sand beaches from wave climate. Coastal Engineering, 4: 253-277.

    Hasselmann, K., 1974. On the spectral dissipation of ocean waves due to whitecapping. Bound-Layer Meteor, 6 (1-2): 107-27.

    Holthuijsen, L. H., Herman, A., and Booij, N., 2003. Phasedecoupled refraction-diffraction for spectral wave models. Coastal Engineering, 49: 291-305.

    Hsu, J. R. C., Silvester, R., and Xia, Y. M., 1989. Generalities on static equilibrium bays. Coastal Engineering, 12 (4): 353-369.

    Jeanson, M., Anthony, E. J., Dolique, F., and Aubry, A., 2013. Wave characteristics and morphological variations of pocket beaches in a coral reef-lagoon setting, Mayotte Island, Indian Ocean. Geomorphology, 182: 190-209.

    Jiang, D. H., and Gao, S., 2002. Relationship between the tidally-induced near-bed shear stress and the distribution of surficial sediments in the Bohai Strait. Acta Sedimentologica Sinica, 20 (4): 663-667 (in Chinese with English abstract).

    Komen, G. J., Cavaleri, L., Donelan, M., Hasselmann, K., Hasselmann, S., and Janssen, P. A. E. M., 1994. Dynamics and Modeling of Ocean Waves. Cambridge University Press, 556pp.

    Kuang, C. P., Pan, Y., Zhang, Y., Liu, S. G., Yang, Y. X., and Zhang, J. B., 2010. Shoreline change modeling on emergency beach nourishment project on west beach of Beidaihe, China. China Ocean Engineering, 24 (2): 277-289.

    Li, F. L., Xia, D. X., Wang, W. H., Wu, S. Y., and Wang, Y. J., 2004. Discussion on the evolution cause and its recovery for the Dengzhou Shoal. Acta Oceanologica Sinica, 26 (6): 65-73 (in Chinese with English abstract).

    Liu, Y. X., Huang, H. J., Qiu, Z. F., and Fan, J. Y., 2013. Detecting coastline change from satellite images based on beach slope estimation in a tidal flat. International Journal of Applied Earth Observation and Geoinformation, 23: 165-176.

    Luan, T., 2011. Current situation and Analysis of the Sandy Coast in the North of Shandong Peninsula. Master thesis. Ocean University of China (in Chinese with English abstract).

    Masselink, G., and Pattiaratchi, C., 2001. Seasonal changes in beach morphology along the sheltered coastline of Perth, Western Australia. Marine Geology, 172 (3-4): 243-263.

    Ojeda, E., Guillén, J., and Ribas, F., 2011. Dynamics of singlebarred embayed beaches. Marine Geology, 280 (1-4): 76-90.

    Pilkey, O. H., and Richter, D. M., 1964. Beach profiles of a Georgia barrier island. Southeastern Geology, 6: 11-19.

    Qi, H. S., Cai, F., Lei, G., Cao, H. M., and Shi, F. Y., 2010. The response of three main beach types to tropical storms in South China. Marine Geology, 275: 244-254.

    Qu, J. J., Han, Q. J., Dong, G. R., Zhang, K. C., and Zu, R. P., 2013. A study of the characteristics of aeolian sand activity and the effects of a comprehensive protective system in a coastal dune area in Southern China. Coastal Engineering, 77: 28-39. Rao, V. R., Murthy, M. V. R., Bhat, M., and Reddy, N. T., 2009. Littoral sediment transport and shoreline changes along Ennore on the southeast coast of India: Field observations and numerical modeling. Geomorphology, 112 (1-2): 158-166.

    Robinson, E., 2004. Coastal changes along the coast of Vere, Jamaica over the past two hundred years: Data from maps and air photographs. Quaternary International, 120 (1): 153-161.

    Sabique, L., Annapurnaiah, K., Balakrishnan Nair, T. M., and Srinivas, K., 2011. Contribution of Southern Indian Ocean swells on the wave heights in the Northern Indian Ocean–A modeling study. Ocean Engineering, 43: 113-120.

    Schwarzer, K., Diesing, M., and Larson, M., 2003. Coastline evolution at different time scales–examples from the Pomeranian Bight, southern Baltic Sea. Marine Geology, 194 (1-2): 79-101.

    SWAN team, 2012. SWAN User Manual. SWAN Cycle III version 40.85. Delft University of Technology, Technical Documentation, 94pp.

    Takaaki, U., Masumi, S., Takayuki, K., and Kazuya, S., 2010. A new model for predicting three-dimensional beach changes by expanding Hsu and Evans’ equation. Coastal Engineering, 57: 194-202.

    Thiebot, J., Idier, D., Garnier, R., Falque′s, A., and Ruessink, B. G., 2012. The influence of wave direction on the morphological response of a double sandbar system. Continental Shelf Research, 32: 71-85.

    Van Rijn, L. C., 2011. Coastal erosion and control. Ocean & Coastal Management, 54 (12): 867-887.

    Wang, W. H., 1994. An experimental study on beach drift using tracing sand on the beach near Xiazhuang Penglai county, Shandong province. Journal of Oceanography of Huanghai & Bohai Sea, 12 (1): 19-27 (in Chinese with English abstract).

    Wang, Z. D., Wu, S. Y., and Yu, J., 2005. Study of the recent geomorphological evolution in the Denzhou Shoal area. Coastal Engineering, 24 (1): 29-38 (in Chinese with English abstract).

    White, K., and El Asmar, H. M., 1999. Monitoring changing position of coastlines using Thematic Mapper imagery, an example from the Nile Delta. Geomorphology, 29 (1-2): 93-105.

    Wright, L. D., and Short, A. D., 1984. Morphodynamic variability of surf zones and beaches: A synthesis. Marine Geology, 56: 93-118.

    Wu, G. Q., Meng, X. D., Xia, D. X., and Cui, J. R., 1997. Research on the effect of Dengzhou Shoal. Oceanologia et Limnologia Sinica, 28 (4): 419-425 (in Chinese with English abstract).

    Wu, S. Y., 1992. The relation between geomorphic environment and coastal erosion in the west of Penglai county, Shandong province. Coastal Engineering, 11 (4): 46-52 (in Chinese with English abstract).

    Xie, M. X., 2012. Three-dimensional numerical modelling of the wave-induced rip currents under irregular bathymetry. Journal of Hydrodynamics, 24 (6): 864-872.

    Yvonnet, J., Coffignal, G., Ryckelynck, D., Lorongb, P., and Chinestab, F., 2006. A simple error indicator for meshfree methods based on natural neighbors. Computers and Structures, 84 (21): 1301-1312.

    Zhang, Y., He, L., Kuang, C., Pan, Y., and Gu, J., 2012. Numerical study of shoreline changes by emergency beach nourishment project at the middle beach of Beidaihe, China. Acta Oceanologica Sinica, 31 (1): 125-133.

    Zheng, J., and Dean, R. G., 1997. Numerical models and intercomparions of beach profile evolution. Coastal Engineering, 30: 169-201.

    (Edited by Xie Jun)

    (Received December 25, 2013; revised March 10, 2014; accepted December 18, 2014)

    J. Ocean Univ. China (Oceanic and Coastal Sea Research)

    DOI 10.1007/s11802-015-2569-1

    ISSN 1672-5182, 2015 14 (5): 803-815

    http://www.ouc.edu.cn/xbywb/

    E-mail:xbywb@ouc.edu.cn

    * Corresponding author. Tel: 0086-532-66781882

    E-mail: hrj@ouc.edu.cn

    肉色欧美久久久久久久蜜桃| 在线观看免费高清a一片| 国产男女超爽视频在线观看| 韩国精品一区二区三区| 久久久久网色| 大码成人一级视频| 国产极品天堂在线| 高清黄色对白视频在线免费看| 亚洲精品一二三| 亚洲精品久久成人aⅴ小说| 妹子高潮喷水视频| 国产免费视频播放在线视频| 国产精品人妻久久久影院| 国产黄色免费在线视频| 香蕉丝袜av| 国产在线免费精品| 在线亚洲精品国产二区图片欧美| 国产老妇伦熟女老妇高清| 亚洲国产欧美网| 精品卡一卡二卡四卡免费| 亚洲精品国产色婷婷电影| 9色porny在线观看| 一本色道久久久久久精品综合| 午夜老司机福利片| 欧美亚洲 丝袜 人妻 在线| 久久免费观看电影| 亚洲国产日韩一区二区| 搡老岳熟女国产| 狂野欧美激情性bbbbbb| 卡戴珊不雅视频在线播放| 亚洲综合色网址| 亚洲美女视频黄频| 欧美日韩视频精品一区| 精品国产露脸久久av麻豆| 久久鲁丝午夜福利片| 国产成人午夜福利电影在线观看| 亚洲av成人不卡在线观看播放网 | 久久人人爽人人片av| 99九九在线精品视频| 久久 成人 亚洲| 国产一区二区 视频在线| 制服人妻中文乱码| 少妇人妻久久综合中文| 男女边摸边吃奶| 午夜福利在线免费观看网站| 亚洲第一青青草原| 精品人妻一区二区三区麻豆| 男人舔女人的私密视频| 女人被躁到高潮嗷嗷叫费观| 久久久久精品性色| 国产成人免费无遮挡视频| av国产精品久久久久影院| 男女国产视频网站| 久久久久精品国产欧美久久久 | 国产亚洲最大av| 一本—道久久a久久精品蜜桃钙片| 大陆偷拍与自拍| 久久精品国产亚洲av涩爱| 午夜福利,免费看| 国产不卡av网站在线观看| avwww免费| 两个人免费观看高清视频| 午夜免费观看性视频| av线在线观看网站| 性色av一级| 国产老妇伦熟女老妇高清| 在线观看免费高清a一片| 两个人看的免费小视频| 成人18禁高潮啪啪吃奶动态图| 久久人人97超碰香蕉20202| 97人妻天天添夜夜摸| 美国免费a级毛片| 热re99久久国产66热| 成年人午夜在线观看视频| 在线观看www视频免费| 日韩 亚洲 欧美在线| 久久人人爽人人片av| 激情视频va一区二区三区| 一边摸一边做爽爽视频免费| 一二三四在线观看免费中文在| 性色av一级| 国产亚洲av高清不卡| 欧美国产精品va在线观看不卡| 国产在视频线精品| 欧美国产精品va在线观看不卡| 欧美 亚洲 国产 日韩一| 国产片特级美女逼逼视频| 精品国产国语对白av| 人人妻人人添人人爽欧美一区卜| 久久久久视频综合| av又黄又爽大尺度在线免费看| 亚洲精品国产一区二区精华液| 国产男人的电影天堂91| 精品一区二区三区av网在线观看 | 免费看不卡的av| 欧美变态另类bdsm刘玥| 免费少妇av软件| 人人妻,人人澡人人爽秒播 | svipshipincom国产片| 成人国语在线视频| 久久久精品国产亚洲av高清涩受| 91aial.com中文字幕在线观看| 男的添女的下面高潮视频| 天天添夜夜摸| 免费高清在线观看日韩| 捣出白浆h1v1| 欧美黄色片欧美黄色片| 亚洲久久久国产精品| 熟女av电影| 成年动漫av网址| 在线观看一区二区三区激情| 黑人巨大精品欧美一区二区蜜桃| 国产一级毛片在线| 国产精品国产三级国产专区5o| 亚洲欧美日韩另类电影网站| av女优亚洲男人天堂| 丝袜脚勾引网站| 天堂俺去俺来也www色官网| 99热全是精品| 99热全是精品| 高清在线视频一区二区三区| 新久久久久国产一级毛片| 日韩精品免费视频一区二区三区| 亚洲第一av免费看| 精品视频人人做人人爽| 久久久欧美国产精品| 免费看不卡的av| 一本一本久久a久久精品综合妖精| 欧美黄色片欧美黄色片| 国产毛片在线视频| 久久亚洲国产成人精品v| 亚洲精华国产精华液的使用体验| 国产精品久久久av美女十八| 亚洲专区中文字幕在线 | 国产精品av久久久久免费| 老司机影院毛片| 交换朋友夫妻互换小说| 在线看a的网站| 欧美日韩成人在线一区二区| 国产一区二区三区av在线| 哪个播放器可以免费观看大片| 久久久久精品久久久久真实原创| 啦啦啦在线观看免费高清www| svipshipincom国产片| 91老司机精品| 最近手机中文字幕大全| 亚洲综合精品二区| 我要看黄色一级片免费的| 免费在线观看视频国产中文字幕亚洲 | 中国三级夫妇交换| 一本色道久久久久久精品综合| 成年女人毛片免费观看观看9 | 欧美日韩av久久| 亚洲欧美中文字幕日韩二区| 亚洲精品日本国产第一区| 不卡av一区二区三区| 亚洲欧美一区二区三区久久| 国产精品熟女久久久久浪| 一级a爱视频在线免费观看| 亚洲成人手机| 一级毛片电影观看| 中文字幕最新亚洲高清| 久久精品久久精品一区二区三区| 最近中文字幕高清免费大全6| 9色porny在线观看| svipshipincom国产片| 伦理电影免费视频| 国产精品熟女久久久久浪| 高清欧美精品videossex| 亚洲精品国产区一区二| 97在线人人人人妻| 深夜精品福利| 亚洲第一区二区三区不卡| 精品福利永久在线观看| 成人国产麻豆网| 爱豆传媒免费全集在线观看| 日韩av免费高清视频| 久久久久精品国产欧美久久久 | 国语对白做爰xxxⅹ性视频网站| 2021少妇久久久久久久久久久| 国产精品一区二区精品视频观看| 久久久久久人妻| 久久婷婷青草| 亚洲第一青青草原| 久久久久久久久久久免费av| 满18在线观看网站| 日韩大码丰满熟妇| 这个男人来自地球电影免费观看 | 天天躁日日躁夜夜躁夜夜| 欧美日韩亚洲高清精品| 日韩一区二区视频免费看| 久久女婷五月综合色啪小说| 国产福利在线免费观看视频| 最黄视频免费看| 免费在线观看黄色视频的| 中文精品一卡2卡3卡4更新| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品久久久久久婷婷小说| 最黄视频免费看| 美女扒开内裤让男人捅视频| 亚洲成人手机| 国语对白做爰xxxⅹ性视频网站| 国产精品久久久久久精品电影小说| 熟妇人妻不卡中文字幕| av视频免费观看在线观看| 九草在线视频观看| 亚洲av日韩在线播放| 日韩av在线免费看完整版不卡| 哪个播放器可以免费观看大片| 国产精品 欧美亚洲| 亚洲成人手机| 韩国av在线不卡| 美女中出高潮动态图| 国产有黄有色有爽视频| 成年美女黄网站色视频大全免费| 18在线观看网站| 一区二区av电影网| 国产精品 欧美亚洲| 日日摸夜夜添夜夜爱| 久久久久网色| 国产精品久久久久久久久免| av在线app专区| 一本一本久久a久久精品综合妖精| 夫妻性生交免费视频一级片| 亚洲一码二码三码区别大吗| 欧美日韩视频精品一区| 中文字幕色久视频| 一本久久精品| 国产精品久久久久久精品古装| 精品人妻一区二区三区麻豆| 午夜福利,免费看| 丝袜美腿诱惑在线| 亚洲精品,欧美精品| 精品久久蜜臀av无| 久久国产亚洲av麻豆专区| 国产黄色视频一区二区在线观看| 欧美亚洲日本最大视频资源| 精品国产一区二区三区四区第35| 日本欧美视频一区| 亚洲男人天堂网一区| 亚洲人成77777在线视频| 亚洲国产精品一区二区三区在线| 久久国产精品大桥未久av| a级毛片在线看网站| 天天操日日干夜夜撸| 日本色播在线视频| 精品午夜福利在线看| 久久精品亚洲熟妇少妇任你| 亚洲欧美激情在线| 巨乳人妻的诱惑在线观看| 欧美黄色片欧美黄色片| 久久人妻熟女aⅴ| 看非洲黑人一级黄片| 美女扒开内裤让男人捅视频| 午夜福利视频精品| 少妇 在线观看| 涩涩av久久男人的天堂| 久久av网站| 国产成人精品久久久久久| 亚洲色图 男人天堂 中文字幕| 国产一区亚洲一区在线观看| av在线播放精品| 免费黄色在线免费观看| 国产免费又黄又爽又色| 国产一区二区三区综合在线观看| 亚洲av欧美aⅴ国产| 免费不卡黄色视频| 日韩制服骚丝袜av| 一级片'在线观看视频| 久久国产亚洲av麻豆专区| 亚洲欧美精品综合一区二区三区| 亚洲成人免费av在线播放| 在线天堂中文资源库| 这个男人来自地球电影免费观看 | 又黄又粗又硬又大视频| 国产精品久久久久久久久免| 精品一品国产午夜福利视频| 国产探花极品一区二区| 亚洲精品日韩在线中文字幕| 成人三级做爰电影| 无限看片的www在线观看| 999久久久国产精品视频| 国产成人一区二区在线| 亚洲精品日韩在线中文字幕| 亚洲国产欧美在线一区| 日本av手机在线免费观看| 亚洲欧美中文字幕日韩二区| 无遮挡黄片免费观看| 成年人免费黄色播放视频| 亚洲精品日韩在线中文字幕| 成人国语在线视频| av国产精品久久久久影院| 国产精品免费视频内射| 精品亚洲成国产av| 亚洲在久久综合| 最新的欧美精品一区二区| 777米奇影视久久| 一本久久精品| 在线亚洲精品国产二区图片欧美| 亚洲,欧美精品.| 亚洲一码二码三码区别大吗| 精品久久久久久电影网| 黑人欧美特级aaaaaa片| 日韩 亚洲 欧美在线| 亚洲国产毛片av蜜桃av| 美女脱内裤让男人舔精品视频| 色吧在线观看| 一级毛片我不卡| 亚洲精品中文字幕在线视频| 亚洲,欧美精品.| 国产黄色免费在线视频| 久久久久久久久免费视频了| 最近2019中文字幕mv第一页| avwww免费| 纵有疾风起免费观看全集完整版| 国产精品三级大全| 久久精品国产综合久久久| 久久久久久久大尺度免费视频| 亚洲男人天堂网一区| 1024视频免费在线观看| 校园人妻丝袜中文字幕| 黄网站色视频无遮挡免费观看| 亚洲国产精品999| 这个男人来自地球电影免费观看 | 欧美日韩精品网址| 国产亚洲av片在线观看秒播厂| 日韩熟女老妇一区二区性免费视频| 成人黄色视频免费在线看| 国产人伦9x9x在线观看| 欧美成人午夜精品| 深夜精品福利| 免费观看a级毛片全部| 精品第一国产精品| 国产日韩欧美视频二区| 尾随美女入室| 欧美日韩视频精品一区| 中文字幕最新亚洲高清| 亚洲五月色婷婷综合| 欧美老熟妇乱子伦牲交| 国产黄频视频在线观看| 国产99久久九九免费精品| 啦啦啦啦在线视频资源| 男人添女人高潮全过程视频| 精品免费久久久久久久清纯 | videos熟女内射| 久久精品国产亚洲av高清一级| e午夜精品久久久久久久| 热re99久久国产66热| 日韩制服丝袜自拍偷拍| 亚洲久久久国产精品| 黄片无遮挡物在线观看| 久久人人爽av亚洲精品天堂| 一级a爱视频在线免费观看| 久热这里只有精品99| 亚洲精品久久久久久婷婷小说| av一本久久久久| 水蜜桃什么品种好| 久久久精品免费免费高清| 亚洲国产欧美一区二区综合| 亚洲精品国产色婷婷电影| 97人妻天天添夜夜摸| 人妻 亚洲 视频| 久久精品亚洲av国产电影网| 久久综合国产亚洲精品| 激情视频va一区二区三区| 国产午夜精品一二区理论片| 日韩大码丰满熟妇| 天美传媒精品一区二区| 亚洲国产精品国产精品| 亚洲第一区二区三区不卡| 熟女少妇亚洲综合色aaa.| 建设人人有责人人尽责人人享有的| av有码第一页| 七月丁香在线播放| 色婷婷久久久亚洲欧美| 久久精品aⅴ一区二区三区四区| 亚洲精品aⅴ在线观看| 欧美国产精品一级二级三级| 亚洲欧美日韩另类电影网站| 久久av网站| 99精品久久久久人妻精品| 久久韩国三级中文字幕| 黑人巨大精品欧美一区二区蜜桃| 丝袜美腿诱惑在线| 少妇人妻久久综合中文| av不卡在线播放| 精品一品国产午夜福利视频| 国产熟女欧美一区二区| 操美女的视频在线观看| 亚洲精品日本国产第一区| 亚洲成人免费av在线播放| 制服丝袜香蕉在线| 欧美精品人与动牲交sv欧美| 国产精品亚洲av一区麻豆 | 欧美av亚洲av综合av国产av | 女性被躁到高潮视频| 久久天堂一区二区三区四区| 国产 一区精品| 欧美日韩国产mv在线观看视频| 日韩一区二区视频免费看| 狠狠婷婷综合久久久久久88av| 午夜激情久久久久久久| 久久精品久久久久久久性| 性高湖久久久久久久久免费观看| 成人三级做爰电影| 黄色视频在线播放观看不卡| 午夜福利,免费看| 亚洲av男天堂| 五月天丁香电影| 国产成人午夜福利电影在线观看| 免费人妻精品一区二区三区视频| 日本爱情动作片www.在线观看| 亚洲精品国产av蜜桃| 国产av精品麻豆| 一本色道久久久久久精品综合| 97精品久久久久久久久久精品| 亚洲精品美女久久av网站| 国产97色在线日韩免费| 女性被躁到高潮视频| 一区福利在线观看| 亚洲精品视频女| 嫩草影院入口| 九九爱精品视频在线观看| 欧美日韩综合久久久久久| 又大又黄又爽视频免费| 日韩欧美精品免费久久| 午夜福利一区二区在线看| 婷婷色综合大香蕉| 99久国产av精品国产电影| 国产色婷婷99| 亚洲av中文av极速乱| 最近中文字幕高清免费大全6| 国产欧美日韩一区二区三区在线| 精品国产露脸久久av麻豆| 男女边吃奶边做爰视频| 亚洲精品中文字幕在线视频| 大陆偷拍与自拍| 欧美少妇被猛烈插入视频| 999久久久国产精品视频| 亚洲欧美精品自产自拍| 欧美日韩亚洲国产一区二区在线观看 | 国产一级毛片在线| 色网站视频免费| 在线观看国产h片| av在线app专区| 精品卡一卡二卡四卡免费| 精品久久久精品久久久| 人妻 亚洲 视频| 国产熟女午夜一区二区三区| 一级毛片我不卡| 人妻一区二区av| 国产成人啪精品午夜网站| 免费黄网站久久成人精品| a级毛片在线看网站| 咕卡用的链子| 久久天堂一区二区三区四区| 一级毛片 在线播放| 韩国高清视频一区二区三区| 99国产综合亚洲精品| 亚洲精品av麻豆狂野| 在线观看一区二区三区激情| 日韩不卡一区二区三区视频在线| 曰老女人黄片| 免费高清在线观看日韩| 日韩成人av中文字幕在线观看| 日韩视频在线欧美| 菩萨蛮人人尽说江南好唐韦庄| 九色亚洲精品在线播放| 久久久国产精品麻豆| 国产男人的电影天堂91| 18禁观看日本| 街头女战士在线观看网站| 五月天丁香电影| 欧美日本中文国产一区发布| 亚洲综合色网址| 18禁观看日本| 交换朋友夫妻互换小说| 91精品伊人久久大香线蕉| 赤兔流量卡办理| 菩萨蛮人人尽说江南好唐韦庄| www日本在线高清视频| 精品少妇黑人巨大在线播放| 最近最新中文字幕免费大全7| 久久久久久免费高清国产稀缺| 51午夜福利影视在线观看| 黄片播放在线免费| 亚洲欧美色中文字幕在线| 波多野结衣一区麻豆| 久久久久网色| 成人亚洲精品一区在线观看| 免费观看性生交大片5| 秋霞伦理黄片| 色婷婷av一区二区三区视频| 妹子高潮喷水视频| 永久免费av网站大全| 人人妻人人添人人爽欧美一区卜| 在线观看一区二区三区激情| 操出白浆在线播放| 国产精品一区二区精品视频观看| 大香蕉久久网| 亚洲国产中文字幕在线视频| 国产在线一区二区三区精| 日韩大码丰满熟妇| 美女国产高潮福利片在线看| 美女脱内裤让男人舔精品视频| 久久99热这里只频精品6学生| 国产av一区二区精品久久| 国产成人欧美在线观看 | 日韩大片免费观看网站| 亚洲欧美中文字幕日韩二区| 欧美日韩亚洲综合一区二区三区_| 啦啦啦中文免费视频观看日本| 两个人免费观看高清视频| 最近最新中文字幕免费大全7| 亚洲精品国产区一区二| 中文天堂在线官网| 亚洲欧美清纯卡通| 波野结衣二区三区在线| 久久久久久久精品精品| 美女中出高潮动态图| 精品国产一区二区久久| 啦啦啦啦在线视频资源| 一级片免费观看大全| 天天躁日日躁夜夜躁夜夜| 日韩,欧美,国产一区二区三区| 免费高清在线观看视频在线观看| 啦啦啦在线免费观看视频4| 午夜免费男女啪啪视频观看| 国产成人91sexporn| 天堂俺去俺来也www色官网| 永久免费av网站大全| av国产久精品久网站免费入址| 欧美黄色片欧美黄色片| www.精华液| 国产野战对白在线观看| av福利片在线| 丰满饥渴人妻一区二区三| 国产成人精品久久二区二区91 | 老司机在亚洲福利影院| 香蕉丝袜av| 久久久久久久精品精品| a级片在线免费高清观看视频| 欧美日韩亚洲国产一区二区在线观看 | 最近的中文字幕免费完整| 久久精品人人爽人人爽视色| 成人国语在线视频| 在线免费观看不下载黄p国产| 免费女性裸体啪啪无遮挡网站| www.精华液| 精品国产超薄肉色丝袜足j| 电影成人av| 在线 av 中文字幕| 国产男女超爽视频在线观看| 亚洲色图综合在线观看| 人妻一区二区av| 国产人伦9x9x在线观看| 日韩视频在线欧美| 性少妇av在线| 伊人久久国产一区二区| 日本一区二区免费在线视频| 狠狠精品人妻久久久久久综合| 久久精品亚洲av国产电影网| 777久久人妻少妇嫩草av网站| 久久精品国产亚洲av高清一级| 中文字幕高清在线视频| 在线天堂中文资源库| 午夜免费鲁丝| 国产欧美亚洲国产| 卡戴珊不雅视频在线播放| 国产欧美亚洲国产| 久久久久久久久免费视频了| 午夜福利乱码中文字幕| 免费不卡黄色视频| 国产野战对白在线观看| xxx大片免费视频| 波多野结衣av一区二区av| 亚洲精品乱久久久久久| 男人添女人高潮全过程视频| av又黄又爽大尺度在线免费看| av在线老鸭窝| 久久久国产一区二区| 国产av一区二区精品久久| 成人影院久久| 人人妻人人爽人人添夜夜欢视频| 青春草亚洲视频在线观看| 欧美老熟妇乱子伦牲交| 国产高清国产精品国产三级| 熟女av电影| 日韩成人av中文字幕在线观看| 久久精品国产亚洲av涩爱| 高清不卡的av网站| 成人亚洲精品一区在线观看| 久久精品国产综合久久久| 丝袜美腿诱惑在线| 老司机在亚洲福利影院| 欧美精品人与动牲交sv欧美| 777久久人妻少妇嫩草av网站| 18在线观看网站| 亚洲欧美日韩另类电影网站| 欧美日韩成人在线一区二区| 亚洲精品国产av蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 丝瓜视频免费看黄片| 免费不卡黄色视频| 午夜福利一区二区在线看| 国产日韩一区二区三区精品不卡| 高清视频免费观看一区二区| 一区二区三区四区激情视频| 国产成人精品在线电影| 人妻人人澡人人爽人人| 成年美女黄网站色视频大全免费| www日本在线高清视频| 满18在线观看网站| 下体分泌物呈黄色| 视频在线观看一区二区三区|