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

    国产精品影院久久| 国产成人啪精品午夜网站| 性色avwww在线观看| 在线观看美女被高潮喷水网站 | 美女大奶头视频| 午夜福利在线观看免费完整高清在 | 2021天堂中文幕一二区在线观| 欧美+日韩+精品| 亚洲色图av天堂| 窝窝影院91人妻| 中文在线观看免费www的网站| 精品日产1卡2卡| 精品人妻熟女av久视频| 毛片女人毛片| 免费人成视频x8x8入口观看| 男女之事视频高清在线观看| 久久久久久大精品| 日韩欧美三级三区| 国产高清激情床上av| 丁香欧美五月| 亚洲自拍偷在线| 免费在线观看影片大全网站| 久久香蕉精品热| 91狼人影院| 丁香六月欧美| 赤兔流量卡办理| 免费在线观看日本一区| 真实男女啪啪啪动态图| 国内久久婷婷六月综合欲色啪| 老司机午夜福利在线观看视频| 国产美女午夜福利| 18禁黄网站禁片午夜丰满| 男插女下体视频免费在线播放| 成人无遮挡网站| 午夜亚洲福利在线播放| 亚洲av电影在线进入| 亚洲无线在线观看| 成年女人永久免费观看视频| 五月玫瑰六月丁香| 精品午夜福利在线看| 精品乱码久久久久久99久播| 婷婷六月久久综合丁香| www日本黄色视频网| 亚洲精品一区av在线观看| 久久精品人妻少妇| 久久热精品热| 日韩人妻高清精品专区| 色哟哟哟哟哟哟| 悠悠久久av| 人妻夜夜爽99麻豆av| 日本黄色视频三级网站网址| 夜夜看夜夜爽夜夜摸| 麻豆成人av在线观看| 国语自产精品视频在线第100页| 757午夜福利合集在线观看| 亚洲一区二区三区不卡视频| 成人精品一区二区免费| 美女被艹到高潮喷水动态| 国产高清三级在线| 精品久久久久久久久久久久久| 亚洲精品456在线播放app | 午夜精品久久久久久毛片777| 精品欧美国产一区二区三| 男人的好看免费观看在线视频| 国产亚洲精品久久久com| 午夜福利视频1000在线观看| 欧美成狂野欧美在线观看| 性色avwww在线观看| 亚洲av一区综合| 亚洲无线在线观看| 嫩草影院精品99| 在线a可以看的网站| 天天一区二区日本电影三级| 久久久久免费精品人妻一区二区| netflix在线观看网站| 91在线观看av| netflix在线观看网站| 免费av观看视频| 欧美日本视频| 一本综合久久免费| 熟妇人妻久久中文字幕3abv| 免费无遮挡裸体视频| 亚洲色图av天堂| 麻豆成人av在线观看| 色哟哟哟哟哟哟| 中文字幕av成人在线电影| 国产三级在线视频| 18禁黄网站禁片免费观看直播| 日韩中字成人| 草草在线视频免费看| 久久中文看片网| 久久人人爽人人爽人人片va | 亚洲精华国产精华精| 男女视频在线观看网站免费| 亚洲国产精品成人综合色| www.www免费av| 中国美女看黄片| 久久午夜亚洲精品久久| 午夜精品久久久久久毛片777| 一级黄片播放器| 在线看三级毛片| 如何舔出高潮| 色尼玛亚洲综合影院| 亚洲国产精品久久男人天堂| 黄色丝袜av网址大全| 久久久久精品国产欧美久久久| 国产在视频线在精品| 国产色婷婷99| 真人做人爱边吃奶动态| 日韩欧美在线乱码| 一本综合久久免费| 伦理电影大哥的女人| 永久网站在线| 我的老师免费观看完整版| 91字幕亚洲| 成人国产一区最新在线观看| 国产三级在线视频| 亚洲片人在线观看| 五月玫瑰六月丁香| 欧美性猛交黑人性爽| 美女xxoo啪啪120秒动态图 | 国产亚洲精品av在线| 日本熟妇午夜| .国产精品久久| 亚洲中文字幕一区二区三区有码在线看| 国产一级毛片七仙女欲春2| 日韩欧美 国产精品| 午夜免费成人在线视频| 热99在线观看视频| 五月玫瑰六月丁香| 一个人免费在线观看的高清视频| 51午夜福利影视在线观看| 久久午夜福利片| 日韩有码中文字幕| 首页视频小说图片口味搜索| 高清在线国产一区| 日日夜夜操网爽| 91狼人影院| 亚洲精品一卡2卡三卡4卡5卡| 日本黄色视频三级网站网址| 国内少妇人妻偷人精品xxx网站| 欧美又色又爽又黄视频| 欧美在线黄色| 高清日韩中文字幕在线| 亚洲精品粉嫩美女一区| 51午夜福利影视在线观看| 美女cb高潮喷水在线观看| 在线免费观看不下载黄p国产 | 日本黄色片子视频| 一卡2卡三卡四卡精品乱码亚洲| 尤物成人国产欧美一区二区三区| 久久久精品大字幕| 亚洲,欧美,日韩| 三级男女做爰猛烈吃奶摸视频| 悠悠久久av| 美女大奶头视频| 我的女老师完整版在线观看| 亚洲无线在线观看| 免费人成视频x8x8入口观看| 麻豆国产av国片精品| 色视频www国产| www.熟女人妻精品国产| 淫妇啪啪啪对白视频| 熟妇人妻久久中文字幕3abv| 午夜两性在线视频| 国内少妇人妻偷人精品xxx网站| 日韩成人在线观看一区二区三区| 久久婷婷人人爽人人干人人爱| 欧美精品国产亚洲| 人人妻人人澡欧美一区二区| 久久精品夜夜夜夜夜久久蜜豆| 久久久久久久久大av| 国产一区二区激情短视频| 我要搜黄色片| 成年人黄色毛片网站| 美女高潮的动态| 国产aⅴ精品一区二区三区波| 久久久久久久久中文| 午夜两性在线视频| 自拍偷自拍亚洲精品老妇| 精品一区二区免费观看| 露出奶头的视频| 一区二区三区免费毛片| 久久久久性生活片| 99久久成人亚洲精品观看| 日日摸夜夜添夜夜添小说| 伊人久久精品亚洲午夜| 不卡一级毛片| 亚洲成人中文字幕在线播放| 久久中文看片网| 亚洲aⅴ乱码一区二区在线播放| 精品熟女少妇八av免费久了| 免费黄网站久久成人精品 | 国产私拍福利视频在线观看| 超碰av人人做人人爽久久| 在线观看午夜福利视频| 在线看三级毛片| 99热这里只有是精品50| 老司机深夜福利视频在线观看| 午夜两性在线视频| 搡老岳熟女国产| 哪里可以看免费的av片| 亚洲五月天丁香| 国产美女午夜福利| 国产大屁股一区二区在线视频| 在线观看舔阴道视频| 国产精品一区二区性色av| 午夜a级毛片| or卡值多少钱| 午夜福利视频1000在线观看| 91麻豆av在线| 免费一级毛片在线播放高清视频| 看十八女毛片水多多多| 午夜福利在线在线| 97超视频在线观看视频| 一区二区三区四区激情视频 | 午夜精品久久久久久毛片777| 特大巨黑吊av在线直播| 国产精品人妻久久久久久| 在线天堂最新版资源| 久久国产乱子免费精品| 男女床上黄色一级片免费看| 一个人观看的视频www高清免费观看| 精品熟女少妇八av免费久了| 国产精品久久电影中文字幕| 好看av亚洲va欧美ⅴa在| 1000部很黄的大片| 亚洲av.av天堂| 日本a在线网址| 亚洲人成电影免费在线| 亚洲天堂国产精品一区在线| 757午夜福利合集在线观看| 午夜福利在线在线| www.999成人在线观看| 免费搜索国产男女视频| 桃色一区二区三区在线观看| 男女那种视频在线观看| 成人国产综合亚洲| 亚洲精品在线观看二区| 免费高清视频大片| 亚洲精品久久国产高清桃花| 校园春色视频在线观看| 色综合婷婷激情| 免费在线观看影片大全网站| 久久99热这里只有精品18| 女人十人毛片免费观看3o分钟| 欧美日韩福利视频一区二区| 亚洲天堂国产精品一区在线| 国产精品乱码一区二三区的特点| 91午夜精品亚洲一区二区三区 | 色在线成人网| 丰满乱子伦码专区| 国产淫片久久久久久久久 | 免费看光身美女| 97碰自拍视频| 国产黄a三级三级三级人| 最近最新免费中文字幕在线| 成人亚洲精品av一区二区| 免费一级毛片在线播放高清视频| 身体一侧抽搐| 丰满的人妻完整版| 岛国在线免费视频观看| 国产69精品久久久久777片| 国产av麻豆久久久久久久| 欧美三级亚洲精品| 亚洲美女视频黄频| 亚洲av熟女| 麻豆av噜噜一区二区三区| 99久久精品一区二区三区| 三级国产精品欧美在线观看| 国产精品一及| 91久久精品国产一区二区成人| 亚洲第一欧美日韩一区二区三区| 韩国av一区二区三区四区| 一区二区三区高清视频在线| 亚洲国产欧洲综合997久久,| 久久九九热精品免费| 久久久久精品国产欧美久久久| 黄色一级大片看看| 精品久久久久久,| 九九热线精品视视频播放| 国产美女午夜福利| 成人精品一区二区免费| 国产乱人视频| 国产av不卡久久| 国产欧美日韩一区二区精品| 成人鲁丝片一二三区免费| 看免费av毛片| 国产精品久久久久久久电影| 91麻豆精品激情在线观看国产| 亚洲av第一区精品v没综合| 97超级碰碰碰精品色视频在线观看| 赤兔流量卡办理| 国产高清视频在线播放一区| 深夜精品福利| a在线观看视频网站| 最近最新中文字幕大全电影3| av福利片在线观看| 少妇丰满av| 啪啪无遮挡十八禁网站| 国产精品久久久久久久久免 | 嫩草影视91久久| 国产一区二区在线av高清观看| 99精品久久久久人妻精品| av在线老鸭窝| 国产精品久久久久久久电影| 久久精品国产清高在天天线| 中国美女看黄片| 又黄又爽又刺激的免费视频.| 全区人妻精品视频| 99久久精品热视频| 69av精品久久久久久| 精品99又大又爽又粗少妇毛片 | 看黄色毛片网站| 久久久久久九九精品二区国产| 少妇的逼好多水| 国产淫片久久久久久久久 | 女同久久另类99精品国产91| 欧美在线黄色| 丝袜美腿在线中文| 内地一区二区视频在线| 亚洲专区国产一区二区| 国产精品久久久久久精品电影| 啪啪无遮挡十八禁网站| 特大巨黑吊av在线直播| 看十八女毛片水多多多| 国产精品久久久久久精品电影| 性色avwww在线观看| 亚洲人成网站高清观看| 三级国产精品欧美在线观看| 亚洲熟妇中文字幕五十中出| 国产真实乱freesex| 成人精品一区二区免费| 最好的美女福利视频网| 亚洲七黄色美女视频| 亚洲avbb在线观看| 变态另类成人亚洲欧美熟女| 男女做爰动态图高潮gif福利片| 国产精品久久久久久人妻精品电影| 亚洲美女搞黄在线观看 | 亚洲精华国产精华精| 亚洲天堂国产精品一区在线| 亚洲中文字幕一区二区三区有码在线看| 天堂√8在线中文| 欧美xxxx性猛交bbbb| 国产精品99久久久久久久久| 国产一区二区在线观看日韩| 亚洲五月天丁香| 日本一本二区三区精品| 老熟妇仑乱视频hdxx| or卡值多少钱| 久久草成人影院| 亚洲av二区三区四区| 宅男免费午夜| 国产高清有码在线观看视频| 日韩欧美国产一区二区入口| 久久精品国产清高在天天线| 淫秽高清视频在线观看| 国产真实乱freesex| 久久草成人影院| 成人高潮视频无遮挡免费网站| 国内揄拍国产精品人妻在线| 在线观看午夜福利视频| av在线老鸭窝| 一进一出抽搐gif免费好疼| 亚洲av电影在线进入| 首页视频小说图片口味搜索| 国产毛片a区久久久久| 日韩欧美一区二区三区在线观看| 99riav亚洲国产免费| 国产探花在线观看一区二区| 国内少妇人妻偷人精品xxx网站| 99国产综合亚洲精品| 特大巨黑吊av在线直播| 在线观看66精品国产| 88av欧美| 熟女人妻精品中文字幕| 中文资源天堂在线| 中文字幕人成人乱码亚洲影| 中国美女看黄片| 淫妇啪啪啪对白视频| 日韩精品青青久久久久久| 我要看日韩黄色一级片| 99久国产av精品| 在线a可以看的网站| 一二三四社区在线视频社区8| 熟女人妻精品中文字幕| 国产精品亚洲美女久久久| 国产精品av视频在线免费观看| 午夜精品在线福利| 99久久99久久久精品蜜桃| 少妇被粗大猛烈的视频| 欧美乱妇无乱码| 在线国产一区二区在线| 成人三级黄色视频| 中亚洲国语对白在线视频| 99热这里只有是精品50| 永久网站在线| 淫秽高清视频在线观看| 国产精品伦人一区二区| 国产毛片a区久久久久| 深夜精品福利| 97超级碰碰碰精品色视频在线观看| 国产精品乱码一区二三区的特点| 午夜a级毛片| 给我免费播放毛片高清在线观看| 99精品在免费线老司机午夜| 国产男靠女视频免费网站| 99久久无色码亚洲精品果冻| 高清在线国产一区| 中文字幕av成人在线电影| 最近中文字幕高清免费大全6 | 制服丝袜大香蕉在线| 内地一区二区视频在线| 蜜桃亚洲精品一区二区三区| 欧美不卡视频在线免费观看| 国产欧美日韩一区二区三| 18禁黄网站禁片午夜丰满| 天天一区二区日本电影三级| 国产精品一及| 日韩有码中文字幕| 成人国产综合亚洲| 国产精品一区二区三区四区久久| 麻豆av噜噜一区二区三区| 国产精品98久久久久久宅男小说| 色哟哟·www| 精品午夜福利视频在线观看一区| 亚洲国产高清在线一区二区三| 欧美一区二区国产精品久久精品| 校园春色视频在线观看| 午夜两性在线视频| 欧美性感艳星| 直男gayav资源| 欧美+亚洲+日韩+国产| 色播亚洲综合网| 亚洲av中文字字幕乱码综合| 97人妻精品一区二区三区麻豆| 久久久色成人| 欧美激情国产日韩精品一区| 五月玫瑰六月丁香| 男人狂女人下面高潮的视频| 成人特级黄色片久久久久久久| 亚洲狠狠婷婷综合久久图片| av黄色大香蕉| 嫩草影视91久久| 午夜福利在线观看吧| 在线播放无遮挡| 色综合亚洲欧美另类图片| 能在线免费观看的黄片| 成人国产综合亚洲| 麻豆国产97在线/欧美| 成人亚洲精品av一区二区| 免费在线观看影片大全网站| 91在线精品国自产拍蜜月| 99热精品在线国产| 一区福利在线观看| 久99久视频精品免费| 少妇高潮的动态图| 丰满乱子伦码专区| 亚洲天堂国产精品一区在线| 欧美日本亚洲视频在线播放| 黄色日韩在线| 国产主播在线观看一区二区| 美女 人体艺术 gogo| 在线观看66精品国产| 亚洲精品粉嫩美女一区| 99热这里只有精品一区| 91午夜精品亚洲一区二区三区 | 免费看美女性在线毛片视频| 亚洲天堂国产精品一区在线| 欧美又色又爽又黄视频| 青草久久国产| 一个人看的www免费观看视频| 亚洲无线在线观看| av在线老鸭窝| 欧美在线黄色| 非洲黑人性xxxx精品又粗又长| 在线观看66精品国产| 天美传媒精品一区二区| 国产成人a区在线观看| 久久精品久久久久久噜噜老黄 | 国内精品美女久久久久久| 两个人视频免费观看高清| 亚洲av中文字字幕乱码综合| 国产伦精品一区二区三区视频9| 99热6这里只有精品| 精品一区二区三区视频在线观看免费| av在线天堂中文字幕| 成人欧美大片| www.色视频.com| 成人鲁丝片一二三区免费| 亚洲精品在线美女| 久久6这里有精品| 中文字幕av在线有码专区| 久久婷婷人人爽人人干人人爱| 国产亚洲欧美98| 12—13女人毛片做爰片一| 最近最新免费中文字幕在线| 国产精品亚洲美女久久久| 长腿黑丝高跟| 亚洲第一电影网av| 村上凉子中文字幕在线| 麻豆成人午夜福利视频| 99视频精品全部免费 在线| 狠狠狠狠99中文字幕| 国产精品人妻久久久久久| 日韩欧美三级三区| 亚洲avbb在线观看| 亚洲精品一区av在线观看| 欧美色视频一区免费| av天堂中文字幕网| 精品免费久久久久久久清纯| 日韩欧美在线乱码| 天堂√8在线中文| 高清在线国产一区| 中文字幕人成人乱码亚洲影| 久久99热6这里只有精品| 成人av一区二区三区在线看| 日韩欧美精品免费久久 | 国产一区二区三区在线臀色熟女| 久久人妻av系列| 久久精品夜夜夜夜夜久久蜜豆| 久久久国产成人精品二区| 1000部很黄的大片| 性色avwww在线观看| 黄色配什么色好看| 噜噜噜噜噜久久久久久91| 成人欧美大片| 欧美黄色片欧美黄色片| 亚洲久久久久久中文字幕| 亚洲精品成人久久久久久| 很黄的视频免费| 午夜福利高清视频| 十八禁人妻一区二区| 在线观看66精品国产| 91久久精品国产一区二区成人| 亚洲午夜理论影院| 亚洲av日韩精品久久久久久密| 欧美一级a爱片免费观看看| 色吧在线观看| 国产久久久一区二区三区| 18禁黄网站禁片免费观看直播| 亚洲成a人片在线一区二区| 91久久精品电影网| 久久久色成人| 一区二区三区激情视频| 夜夜爽天天搞| 又黄又爽又免费观看的视频| 超碰av人人做人人爽久久| 在线国产一区二区在线| 久久精品影院6| 欧美一级a爱片免费观看看| 免费电影在线观看免费观看| 99热这里只有精品一区| 3wmmmm亚洲av在线观看| 欧美中文日本在线观看视频| 村上凉子中文字幕在线| 最新在线观看一区二区三区| 神马国产精品三级电影在线观看| 精品人妻熟女av久视频| 免费在线观看影片大全网站| 欧美色视频一区免费| 两个人的视频大全免费| 男人的好看免费观看在线视频| 首页视频小说图片口味搜索| 91午夜精品亚洲一区二区三区 | 91久久精品国产一区二区成人| 少妇熟女aⅴ在线视频| 国产黄色小视频在线观看| 亚洲一区二区三区不卡视频| 久久久久久大精品| 国产精品人妻久久久久久| 亚洲av电影在线进入| 看黄色毛片网站| 亚洲最大成人av| 18+在线观看网站| 综合色av麻豆| 色综合婷婷激情| 亚洲最大成人av| 人人妻人人看人人澡| 国产亚洲欧美在线一区二区| 窝窝影院91人妻| 国产亚洲欧美在线一区二区| 午夜视频国产福利| 淫秽高清视频在线观看| 在线播放无遮挡| 极品教师在线免费播放| 国产伦一二天堂av在线观看| 偷拍熟女少妇极品色| 91狼人影院| 午夜亚洲福利在线播放| 色5月婷婷丁香| 观看美女的网站| 少妇熟女aⅴ在线视频| 岛国在线免费视频观看| 99久久九九国产精品国产免费| 国产伦人伦偷精品视频| 亚洲国产精品合色在线| 亚洲黑人精品在线| 国产视频一区二区在线看| 桃红色精品国产亚洲av| 国产一级毛片七仙女欲春2| 欧美色视频一区免费| ponron亚洲| 少妇人妻精品综合一区二区 | 久久国产精品人妻蜜桃| 欧美午夜高清在线| 男女那种视频在线观看| 国产精品一区二区三区四区久久| 99热精品在线国产| 亚洲无线在线观看| 国产亚洲精品av在线| 日本免费a在线| 欧美激情久久久久久爽电影| 一个人观看的视频www高清免费观看|