• <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视频的电影网站 | 视频在线观看一区二区三区| 亚洲一区中文字幕在线| 中文字幕高清在线视频| 国产亚洲av高清不卡| 波多野结衣av一区二区av| 黄色一级大片看看| 人人妻人人澡人人看| 精品一区在线观看国产| 老司机影院成人| 自线自在国产av| 亚洲,欧美,日韩| 久久天躁狠狠躁夜夜2o2o | 成人三级做爰电影| 日日爽夜夜爽网站| 久久久精品区二区三区| 亚洲成av片中文字幕在线观看| 成人亚洲精品一区在线观看| 久久久久久久大尺度免费视频| 日韩 亚洲 欧美在线| 日日啪夜夜爽| 捣出白浆h1v1| 日韩av在线免费看完整版不卡| 99热国产这里只有精品6| 亚洲国产欧美在线一区| 国产免费一区二区三区四区乱码| 久久久久久免费高清国产稀缺| 99久久人妻综合| 国产又色又爽无遮挡免| 午夜福利在线免费观看网站| 成人黄色视频免费在线看| 大码成人一级视频| 一区二区三区乱码不卡18| 国产精品免费视频内射| 日韩 欧美 亚洲 中文字幕| 蜜桃在线观看..| 黄频高清免费视频| 色婷婷av一区二区三区视频| 少妇猛男粗大的猛烈进出视频| 免费黄色在线免费观看| 国产免费又黄又爽又色| 日韩欧美一区视频在线观看| 国产精品av久久久久免费| 人人妻人人添人人爽欧美一区卜| 老司机亚洲免费影院| 黄色 视频免费看| 免费少妇av软件| 成人免费观看视频高清| 国产成人精品久久二区二区91 | 国产免费视频播放在线视频| 精品国产乱码久久久久久小说| 欧美亚洲日本最大视频资源| 久久97久久精品| 在线观看三级黄色| 精品久久蜜臀av无| 国产精品一国产av| 欧美亚洲 丝袜 人妻 在线| 欧美精品人与动牲交sv欧美| 男女免费视频国产| 最近中文字幕2019免费版| 亚洲精华国产精华液的使用体验| 午夜影院在线不卡| 欧美日韩亚洲综合一区二区三区_| 日韩av不卡免费在线播放| 美女扒开内裤让男人捅视频| 国产av一区二区精品久久| 自线自在国产av| 岛国毛片在线播放| av网站免费在线观看视频| 欧美日韩成人在线一区二区| 两性夫妻黄色片| 极品人妻少妇av视频| 久久久久久久精品精品| 妹子高潮喷水视频| 色网站视频免费| 91aial.com中文字幕在线观看| 久久热在线av| 免费黄网站久久成人精品| av有码第一页| 国产 一区精品| 观看美女的网站| 悠悠久久av| 成年人免费黄色播放视频| 哪个播放器可以免费观看大片| 日本爱情动作片www.在线观看| 亚洲国产成人一精品久久久| 欧美日韩亚洲综合一区二区三区_| 国产成人欧美在线观看 | 亚洲精品在线美女| 国产精品.久久久| 久久精品久久久久久噜噜老黄| 黄色毛片三级朝国网站| 国产成人系列免费观看| 亚洲精品自拍成人| 一区二区三区四区激情视频| 欧美精品高潮呻吟av久久| 国产女主播在线喷水免费视频网站| 在线观看三级黄色| 妹子高潮喷水视频| 香蕉丝袜av| 午夜激情久久久久久久| 韩国高清视频一区二区三区| 国产男女超爽视频在线观看| 国产精品免费大片| 999精品在线视频| 啦啦啦 在线观看视频| 中文字幕人妻熟女乱码| 亚洲成人手机| 成年动漫av网址| 色精品久久人妻99蜜桃| 亚洲,一卡二卡三卡| 九九爱精品视频在线观看| 欧美激情高清一区二区三区 | 午夜免费男女啪啪视频观看| 成人国产麻豆网| 一本久久精品| 亚洲人成网站在线观看播放| 国产麻豆69| 亚洲精品国产av成人精品| 天天影视国产精品| 91国产中文字幕| 亚洲精品久久午夜乱码| 无遮挡黄片免费观看| 国产精品久久久久久人妻精品电影 | 亚洲av在线观看美女高潮| 99国产综合亚洲精品| 国产精品人妻久久久影院| 国产精品久久久久久精品电影小说| 国产av一区二区精品久久| 青草久久国产| 男男h啪啪无遮挡| 国产亚洲精品第一综合不卡| 久久免费观看电影| 91老司机精品| 精品一区二区三卡| 成人黄色视频免费在线看| 大片免费播放器 马上看| 国产日韩欧美在线精品| 欧美精品一区二区免费开放| 午夜福利影视在线免费观看| 高清视频免费观看一区二区| 十分钟在线观看高清视频www| 在线 av 中文字幕| 久久久久精品人妻al黑| 2018国产大陆天天弄谢| 亚洲精品在线美女| 亚洲欧美精品自产自拍| 黄色怎么调成土黄色| 国产爽快片一区二区三区| 国产熟女欧美一区二区| 欧美久久黑人一区二区| 午夜福利视频在线观看免费| 国产精品成人在线| 亚洲自偷自拍图片 自拍| 少妇被粗大的猛进出69影院| 国产成人91sexporn| 曰老女人黄片| 午夜福利在线免费观看网站| 国产亚洲最大av| 国产精品久久久久久精品电影小说| 男女无遮挡免费网站观看| 超色免费av| av卡一久久| 日本欧美视频一区| 成人黄色视频免费在线看| 亚洲婷婷狠狠爱综合网| 久久av网站| 精品国产一区二区久久| 免费黄色在线免费观看| 色94色欧美一区二区| 日韩av免费高清视频| 亚洲熟女毛片儿| 久久久国产精品麻豆| 日韩中文字幕视频在线看片| 黄色 视频免费看| 中文字幕另类日韩欧美亚洲嫩草| 黑人欧美特级aaaaaa片| 国产免费又黄又爽又色| 视频区图区小说| 老司机亚洲免费影院| 一本一本久久a久久精品综合妖精| 日韩制服骚丝袜av| 亚洲熟女毛片儿| 夫妻性生交免费视频一级片| 丰满乱子伦码专区| videos熟女内射| 色精品久久人妻99蜜桃| 午夜福利影视在线免费观看| 国产欧美亚洲国产| 国产激情久久老熟女| 亚洲av综合色区一区| 亚洲成人av在线免费| 男女高潮啪啪啪动态图| 国产乱人偷精品视频| 天天躁夜夜躁狠狠久久av| 啦啦啦视频在线资源免费观看| 欧美精品高潮呻吟av久久| 汤姆久久久久久久影院中文字幕| 卡戴珊不雅视频在线播放| 搡老岳熟女国产| 欧美日韩视频精品一区| 九色亚洲精品在线播放| 亚洲av日韩精品久久久久久密 | 综合色丁香网| 秋霞伦理黄片| e午夜精品久久久久久久| 日韩 亚洲 欧美在线| 午夜福利免费观看在线| 成人漫画全彩无遮挡| 国产日韩欧美亚洲二区| 最新的欧美精品一区二区| 国精品久久久久久国模美| tube8黄色片| 免费人妻精品一区二区三区视频| 老司机亚洲免费影院| 久久鲁丝午夜福利片| 国产一区亚洲一区在线观看| 国产欧美日韩一区二区三区在线| 久久久久久久久久久久大奶| 又黄又粗又硬又大视频| 19禁男女啪啪无遮挡网站| 亚洲国产欧美在线一区| 97在线人人人人妻| 在线观看免费午夜福利视频| 久久人人爽人人片av| 大陆偷拍与自拍| 欧美久久黑人一区二区| 久久精品aⅴ一区二区三区四区| 少妇被粗大猛烈的视频| 一级黄片播放器| 黑丝袜美女国产一区| 精品一区在线观看国产| 欧美97在线视频| 观看av在线不卡| 久久久久网色| 黑丝袜美女国产一区| 国产精品女同一区二区软件| av视频免费观看在线观看| 久久久久精品久久久久真实原创| 卡戴珊不雅视频在线播放| 咕卡用的链子| 两性夫妻黄色片| 亚洲av电影在线观看一区二区三区| 久久人人爽人人片av| kizo精华| 哪个播放器可以免费观看大片| 91精品国产国语对白视频| 亚洲国产毛片av蜜桃av| av网站免费在线观看视频| 91老司机精品| 男女午夜视频在线观看| 纵有疾风起免费观看全集完整版| 两个人看的免费小视频| 久久久精品区二区三区| 国产精品国产av在线观看| 操出白浆在线播放| 我要看黄色一级片免费的| tube8黄色片| 一级片免费观看大全| 自拍欧美九色日韩亚洲蝌蚪91| 日日爽夜夜爽网站| 亚洲精品第二区| 成人影院久久| 美女大奶头黄色视频| 亚洲精品久久久久久婷婷小说| 久久影院123| 性色av一级| 国产成人系列免费观看| 永久免费av网站大全| 国产成人欧美在线观看 | 大话2 男鬼变身卡| 日日撸夜夜添| 亚洲欧美清纯卡通| 丝袜喷水一区| 国产成人精品久久二区二区91 | 最近的中文字幕免费完整| 亚洲精品美女久久久久99蜜臀 | 午夜久久久在线观看| 麻豆乱淫一区二区| 成人午夜精彩视频在线观看| 青草久久国产| 久久午夜综合久久蜜桃| 久久精品久久久久久久性| 热99久久久久精品小说推荐| 91国产中文字幕| 91精品伊人久久大香线蕉| 狠狠精品人妻久久久久久综合| 久久久国产精品麻豆| 亚洲美女视频黄频| 咕卡用的链子| 美女午夜性视频免费| 久久久亚洲精品成人影院| 久久人人爽av亚洲精品天堂| 亚洲国产欧美日韩在线播放| 亚洲中文av在线| 国产男人的电影天堂91| 久久久国产欧美日韩av| 国产一区二区在线观看av| 飞空精品影院首页| av免费观看日本| 亚洲欧美激情在线| 1024视频免费在线观看| 下体分泌物呈黄色| 精品人妻熟女毛片av久久网站| 国产精品欧美亚洲77777| 涩涩av久久男人的天堂| 久久久久精品久久久久真实原创| 女的被弄到高潮叫床怎么办| 美女主播在线视频| 99热全是精品| 1024视频免费在线观看| 侵犯人妻中文字幕一二三四区| 高清黄色对白视频在线免费看| 亚洲欧美一区二区三区久久| 亚洲人成77777在线视频| 久久久久久久精品精品| 国产在线视频一区二区| 我的亚洲天堂| 久久精品国产亚洲av涩爱| www.熟女人妻精品国产| 久久久欧美国产精品| 亚洲国产中文字幕在线视频| 欧美黑人精品巨大| 美女扒开内裤让男人捅视频| 精品人妻一区二区三区麻豆| 亚洲国产av新网站| 嫩草影视91久久| 波多野结衣一区麻豆| 建设人人有责人人尽责人人享有的| 老鸭窝网址在线观看| 女性被躁到高潮视频| 免费观看av网站的网址| 在现免费观看毛片| 永久免费av网站大全| 日本欧美国产在线视频| 男女下面插进去视频免费观看| 欧美人与性动交α欧美软件| 欧美日韩亚洲综合一区二区三区_| 亚洲成人手机| 最近中文字幕2019免费版| 大片免费播放器 马上看| 精品酒店卫生间| 人人妻人人添人人爽欧美一区卜| 国产一区二区三区综合在线观看| 最近中文字幕高清免费大全6| 中文字幕人妻丝袜一区二区 | 午夜福利免费观看在线| 国产免费一区二区三区四区乱码| 五月开心婷婷网| 精品久久蜜臀av无| 夜夜骑夜夜射夜夜干| 建设人人有责人人尽责人人享有的| 在线天堂中文资源库| 精品午夜福利在线看| 老汉色∧v一级毛片| 人成视频在线观看免费观看| 国产1区2区3区精品| 国产淫语在线视频| 日韩伦理黄色片| 九色亚洲精品在线播放| 大陆偷拍与自拍| 啦啦啦在线免费观看视频4| 人妻 亚洲 视频| 视频区图区小说| 麻豆乱淫一区二区| 黄频高清免费视频| 国产99久久九九免费精品| h视频一区二区三区| 午夜老司机福利片| 一区二区三区激情视频| 亚洲欧美精品自产自拍| 丝袜美足系列| 高清视频免费观看一区二区| 少妇的丰满在线观看| 飞空精品影院首页| 高清黄色对白视频在线免费看| 黄色毛片三级朝国网站| 久久精品aⅴ一区二区三区四区| 亚洲av成人精品一二三区| av女优亚洲男人天堂| 九草在线视频观看| 大片电影免费在线观看免费| 一区二区三区激情视频| 你懂的网址亚洲精品在线观看| 欧美精品av麻豆av| 亚洲精品一区蜜桃| 亚洲精品,欧美精品| 亚洲四区av| 19禁男女啪啪无遮挡网站| 熟女少妇亚洲综合色aaa.| 在线亚洲精品国产二区图片欧美| 人成视频在线观看免费观看| 高清不卡的av网站| 久久精品熟女亚洲av麻豆精品| 哪个播放器可以免费观看大片| 久热爱精品视频在线9| 精品免费久久久久久久清纯 | 欧美变态另类bdsm刘玥| 一区在线观看完整版| 精品少妇一区二区三区视频日本电影 | 中文字幕色久视频| 日韩欧美一区视频在线观看| 两个人看的免费小视频| av有码第一页| 无限看片的www在线观看| 久久人妻熟女aⅴ| a级毛片黄视频| 天天躁夜夜躁狠狠久久av| 亚洲五月色婷婷综合| 中文欧美无线码| 老汉色∧v一级毛片| 日本一区二区免费在线视频| av在线老鸭窝| 男女国产视频网站| 亚洲,欧美,日韩| 亚洲精品av麻豆狂野| 波野结衣二区三区在线| 国产精品嫩草影院av在线观看| 一级片免费观看大全| 国产乱来视频区| 一区在线观看完整版| 自线自在国产av| 十分钟在线观看高清视频www| 日本欧美视频一区| 国产精品一区二区在线不卡| www.av在线官网国产| av国产久精品久网站免费入址| 亚洲第一区二区三区不卡| 亚洲一区中文字幕在线| 久久99精品国语久久久| 在线亚洲精品国产二区图片欧美| 最近最新中文字幕免费大全7| 亚洲国产看品久久| 国产精品一区二区在线不卡| av有码第一页| 日本欧美国产在线视频| 免费高清在线观看视频在线观看| 亚洲精品国产av蜜桃| 国产97色在线日韩免费| 两个人免费观看高清视频| 久久人人爽av亚洲精品天堂| 亚洲av电影在线进入| 午夜福利影视在线免费观看| 亚洲激情五月婷婷啪啪| 女人被躁到高潮嗷嗷叫费观| 亚洲成国产人片在线观看| h视频一区二区三区| 9色porny在线观看| 午夜福利在线免费观看网站| 亚洲国产精品成人久久小说| 黄片无遮挡物在线观看| 伊人久久国产一区二区| 高清av免费在线| 久久av网站| 男男h啪啪无遮挡| 亚洲七黄色美女视频| 日韩 亚洲 欧美在线| 如日韩欧美国产精品一区二区三区| 纯流量卡能插随身wifi吗| 少妇人妻久久综合中文| 久久精品国产亚洲av涩爱| 日日爽夜夜爽网站| 中文字幕另类日韩欧美亚洲嫩草| 侵犯人妻中文字幕一二三四区| 国产熟女午夜一区二区三区| 一个人免费看片子| av天堂久久9| 国产精品.久久久| 亚洲四区av| 18在线观看网站| 欧美成人精品欧美一级黄| e午夜精品久久久久久久| 在线观看免费视频网站a站| 精品一区二区三区av网在线观看 | 亚洲,一卡二卡三卡| 在线天堂最新版资源| 一二三四在线观看免费中文在| 国产亚洲午夜精品一区二区久久| 亚洲第一青青草原| 毛片一级片免费看久久久久| 国产免费福利视频在线观看| 在线观看www视频免费| 菩萨蛮人人尽说江南好唐韦庄| 国产精品香港三级国产av潘金莲 | 日韩一区二区三区影片| 又大又黄又爽视频免费| 久久精品熟女亚洲av麻豆精品| 青青草视频在线视频观看| 色94色欧美一区二区| 亚洲一区中文字幕在线| 哪个播放器可以免费观看大片| 亚洲成人免费av在线播放| 国产精品三级大全| 黄色视频不卡| 亚洲一区中文字幕在线| 日本午夜av视频| 波多野结衣av一区二区av| 中文字幕最新亚洲高清| 亚洲人成77777在线视频| 黄色视频在线播放观看不卡| 免费观看a级毛片全部| 精品亚洲成a人片在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 人人妻人人爽人人添夜夜欢视频| av电影中文网址| 这个男人来自地球电影免费观看 | 自拍欧美九色日韩亚洲蝌蚪91| 久久av网站| 咕卡用的链子| 超色免费av| 精品亚洲乱码少妇综合久久| 精品一区在线观看国产| 人人妻人人澡人人爽人人夜夜| 老司机在亚洲福利影院| 69精品国产乱码久久久| 国产成人系列免费观看| 欧美日韩视频高清一区二区三区二| 香蕉丝袜av| 建设人人有责人人尽责人人享有的| 极品人妻少妇av视频| 国产精品国产av在线观看| 嫩草影院入口| 国产黄色免费在线视频| 日韩制服丝袜自拍偷拍| 激情五月婷婷亚洲| 免费高清在线观看日韩| 欧美成人午夜精品| 一级毛片 在线播放| 久久亚洲国产成人精品v| 午夜日韩欧美国产| a级毛片在线看网站| 伊人久久大香线蕉亚洲五| 男女边吃奶边做爰视频| 午夜福利免费观看在线| 亚洲欧美日韩另类电影网站| 99热网站在线观看| kizo精华| 夜夜骑夜夜射夜夜干| 中文欧美无线码| 欧美av亚洲av综合av国产av | 中文字幕另类日韩欧美亚洲嫩草| 在线观看www视频免费| 爱豆传媒免费全集在线观看| 国产成人一区二区在线| 国产精品熟女久久久久浪| av免费观看日本| 国产淫语在线视频| 国产成人a∨麻豆精品| 亚洲激情五月婷婷啪啪| 女人高潮潮喷娇喘18禁视频| 亚洲情色 制服丝袜| 男女之事视频高清在线观看 | 精品国产超薄肉色丝袜足j| 国产成人精品无人区| av福利片在线| 自线自在国产av| 精品久久久久久电影网| 国产一区有黄有色的免费视频| 一级毛片电影观看| 男女床上黄色一级片免费看| 毛片一级片免费看久久久久| 大话2 男鬼变身卡| 免费黄色在线免费观看| 操出白浆在线播放| 在线观看一区二区三区激情| av一本久久久久| 老汉色av国产亚洲站长工具| 咕卡用的链子| 伊人亚洲综合成人网| 亚洲精品国产区一区二| 免费不卡黄色视频| 欧美在线黄色| 精品一区在线观看国产| 宅男免费午夜| 国产精品女同一区二区软件| 99国产精品免费福利视频| 人人妻人人添人人爽欧美一区卜| 精品少妇黑人巨大在线播放| 最近中文字幕高清免费大全6| 黄色一级大片看看| 嫩草影院入口| 国产有黄有色有爽视频| 欧美国产精品va在线观看不卡| 亚洲第一区二区三区不卡| www日本在线高清视频| 伊人亚洲综合成人网| 日韩一区二区三区影片| 亚洲美女黄色视频免费看| 国产精品久久久久久精品古装| 亚洲天堂av无毛| 亚洲av在线观看美女高潮| 一边亲一边摸免费视频| a级片在线免费高清观看视频| 日本爱情动作片www.在线观看| 精品国产一区二区三区久久久樱花| 高清欧美精品videossex| 纵有疾风起免费观看全集完整版| 狂野欧美激情性bbbbbb| 欧美精品一区二区免费开放| 成年av动漫网址| 777米奇影视久久| 大香蕉久久网| 久热爱精品视频在线9| 亚洲精品乱久久久久久| 国产一区二区激情短视频 | 综合色丁香网| 电影成人av| 丝袜喷水一区| 少妇猛男粗大的猛烈进出视频| av在线观看视频网站免费| 亚洲国产看品久久| 99九九在线精品视频| 日本色播在线视频| 国产午夜精品一二区理论片| 热99国产精品久久久久久7|