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

    Hydrodynamic and morphological processes in Yangtze Estuary: State-of-the-art research and its applications by Hohai University

    2012-08-11 15:03:27JinhaiZHENGYixinYANChaofengTONGZhiyiLEIChiZHANG
    Water Science and Engineering 2012年4期

    Jin-hai ZHENG*, Yi-xin YAN, Chao-feng TONG, Zhi-yi LEI, Chi ZHANG

    1. State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, P. R. China

    2. College of Harbor, Coastal and Offshore Engineering, Hohai University, Nanjing 210098, P. R. China

    Hydrodynamic and morphological processes in Yangtze Estuary: State-of-the-art research and its applications by Hohai University

    Jin-hai ZHENG*1,2, Yi-xin YAN2, Chao-feng TONG1,2, Zhi-yi LEI2, Chi ZHANG1,2

    1. State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, P. R. China

    2. College of Harbor, Coastal and Offshore Engineering, Hohai University, Nanjing 210098, P. R. China

    This paper presents a review of the state-of-the-art research and its applications developed at Hohai University relating to the hydrodynamic and morphological processes in the Yangtze Estuary. Longitudinal, lateral, and horizontal flow circulations have been revealed based on the measurements with acoustic Doppler current profilers (ADCP). The hydrodynamic mechanism at diversion points as well as the changing patterns of flow and sediment flux in the Yangtze Estuary has been investigated through long-term data analysis. A field survey has been carried out to detect the saltwater intrusion from the North Branch to South Branch. Different numerical models of flow motion, sediment transport, and saltwater intrusion have been developed to simulate the complicated processes and to evaluate the effects of engineering projects. The morphological processes of wetlands over a time scale of decades have been analyzed with an established database. Ideas for further research on the bio-geomorphological model system and long-term evolution mechanisms are put forward.

    Yangtze Estuary; hydrodynamic process; morphological process; saltwater intrusion; Hohai University

    1 Introduction

    The Yangtze Estuary has three diversion points and four outlets. With a 90-km entrance opening to the East China Sea, it accepts a huge tidal volume from the sea, and transmits tidal waves 642 km upstream to Datong, Anhui, where the tidal limit is located in the dry season, and 153 km upstream to Jiangyin, Jiangsu, where the tidal current limit is located in the flood season (Zhong 1985), as shown in Fig. 1.

    The Yangtze Estuary has a multi-annual average precipitation of 1 000 mm to 1 100 mm with quite a large amount of interannual variability. The precipitation in wet years is about1 200 mm, up to 1 400 mm at maximum, and in dry years it varies from 600 to 700 mm.

    According to statistical analysis of the measured data at the Datong Station from 1951 to 1988, the multi-annual average sediment concentration was 0.403 kg/m3, the highest annual average concentration was 3.14 kg/m3, and the lowest one was 0.28 kg/m3; the multi-annual average sediment load was 4.60 × 108t, with the largest annual load being 6.72 × 108t in 1964, and the smallest one being 3.11 × 108t in 1986.

    Fig. 1Tidal limit and tidal current limit of Yangtze Estuary

    According to statistical data from the Datong Hydrological Station from 1951 to 2005, the multi-annual average discharge was 28 226 m3/s; the multi-annual average runoff was 8.853 × 1011m3, with the maximum one being 1.358 9 × 1012m3in 1954, and the minimum one being 6.758 × 1011m3in 1978.

    In the Yangtze Estuary, tides are controlled by tidal waves propagating from the sea. The East China Sea’s progressive wave system includes the M2 partial tide as the major tide in this region, and is impacted by rotational tidal waves of the Yellow Sea, with the K1 and O1 tides being most notable. The tides in the estuary are non-formal neritic semidiurnal tides, with obviously different daily tides, especially at high tides. From the spring equinox to the autumnal equinox, night tides generally outdo diurnal tides, while from the autumnal equinox to the spring equinox of the following year, diurnal tides outdo night tides.

    Wind waves prevail in the Yangtze Estuary. For mixed waves, the cases in which wind waves prevail and swells prevail account for 77% and 23%, respectively. It is rare for only swells to occur. The dominant wave direction is NNE, with a frequency of 10.25%. The seasonal change of wave directions is significant: NW prevails in winter, SSE prevails in summer, and NE often prevails in spring and autumn.

    Since the Yangtze Estuary is characterized by these complicated dynamic processes, it has a significant value to scientific research and engineering studies. In recent decades, relevant studies on the Yangtze Estuary have been continuously conducted by various institutes in China. As one of the top universities in estuarine and coastal engineering, HohaiUniversity has played an important role in the development of fundamental theory and engineering applications in this region. This paper presents a general review of the state-of-the-art research and engineering studies conducted by Hohai University with respect to hydrodynamic and morphological processes, such as hydrodynamic characteristics, saltwater intrusion, and wetland evolution.

    2 Hydrodynamic characteristics

    2.1 Circulation feature analysis based on ADCP measurements

    The location of the Yangtze Estuary and the field measurement sections in the North Passage are shown in Fig. 2. The field data used herein were all primarily measured by acoustic Doppler current profilers (ADCP) in the Yangtze Estuary from September 22 to 30, 2002, and in the North Passage from October 21 to 23, 2002 (Yu 2005; Yan et al. 2007). The first step for circulation feature analysis was extracting the flow velocity in three directions: the west-east, north-south, and vertical directions. The next step was making some necessary transforms and comparisons among the relevant field data. Then, as the main task, the flow regime characteristics of the Yangtze Estuary were examined with a numerical model developed by Hohai University. Typical characteristics of circulation were found to exist in the flow regime of the Yangtze Estuary. The flow circulation also has significant effects on the transport and diffusion of dissolved pollutants in the estuary (Wu and Yan 2010). The possible generation mechanisms and some main factors were analyzed.

    Fig. 2Location of Yangtze Estuary and field measurement sections in North Passage

    It was found that there are three kinds of circulation in the Yangtze Estuary, shown in Figs. 3 through 5: during the low slack tide, longitudinal circulation with the surface layer seaward and the bottom layer landward (Fig. 3(b)), horizontal circulation with the north side landward and the south side seaward (Fig. 4(b)), and lateral circulation with the surface layersouthward and the bottom layer northward (Fig. 5(b)). Both the lateral circulation and horizontal circulation are mainly driven by Coriolis forces, and these three kinds of circulation always appear during the turning process of tides and move with tidal currents. If encountering a strong tidal current, the circulation will disappear. This indicates that, during the turning process of a tide, water flow moves in a complicated three-dimensional spiral formation under the influence of the interaction between river discharge and tidal currents. These three kinds of circulation have relationships with the density gradient induced by saltwater intrusion.

    Fig. 3Longitudinal velocity variations at section #1 of North Passage (The rightward arrowhead orients downstream)

    Fig. 4Velocity variations along horizontal lines on section #1 of North Passage (The upward arrowhead orients downstream, and the ADCP lines are arranged from the north side to south side)

    Fig. 5Cross-sectional velocity fields of section #1 of North Passage

    2.2 Hydrodynamic mechanism at diversion points

    Some methods such as theory analysis, laboratory experiments, and numerical simulation were used to study the hydrodynamic features and mechanisms at the diversion points of the Yangtze Estuary. It was shown that the diversion points rise and fall periodically, with the channels flourishing and withering alternatively, and the waterways silting periodically (Yan et al. 2000, 2001a, 2001b, 2001c, 2002, 2003; Tong 2005; Lei et al. 2009).

    The diversion point of the North Channel and South Channel is located at the Central Shoal. This is a changing diversion point, which swings between Liuhekou and Wusongkou. The Biandan Shoal and Liuhe Shoal are scoured and move downstream alternatively. The two shoals are cut frequently because they block the upper outlets of the North and South channels. The Biandan Shoal was cut in 1924 and 1978, with the sediment falling into the sea through the North Channel, and the Liuhe Shoal was cut in 1958 and 1963. Since the 1980s, the Yangtze Estuary has tended to be stable with the law of bed development still acting. At the outlets of the North Channel and South Channel, the Xinqiao Channel heading to the North Channel and the Nanshatou Shoal heading to the South Channel form an angle of 84°. The lower section of the Nanshatou Shoal is atrophied. Upstream, the Xinliuhe Shoal was cut and its flow was sped up by the flood in 1998, with the separated section moving downstream and uniting with the Central Shoal. Also, the diversion point of the North Channel and South Channel moved upstream.

    The river regime of the North Passage is healthy owing to the good boundary conditions: (1) the joint of the North Passage with the main channel of the South Channel is straight; (2)the diversion angle between the North Passage and South Passage is an acute angle; (3) the hydrodynamic features of the velocity field at the diversion point of the North Passage and South Passage are in favor of the bed sand discharge from the South Passage; and (4) a rotary current, in the Jigu Reef area outside the outlet of the North Passage, is advantageous to the discharge through the North Passage.

    In 1963, when the North Channel and South Channel were cut at their diversion point, a tremendous amount of sediment entered the South Channel, forming a central sand ridge in the South Channel during the 1970s. The central sand ridge was divided by the main flow from the South Channel during the 1980s, with its lower section forming the Jiangya’nan Shoal. The Jiangya’nan Shoal extended northward, moved downstream, and connected with the head of the Jiuduan Shoal during the 1990s. Because of this, the diversion point of the North Passage and South Passage moved upward to the head of the Jiangya’nan Shoal. During the 1990s, the Jiangya’nan Shoal was scoured downward by the ebb tide, and was restrained from moving downstream by a fish-mouth dike at the diversion point of the North Passage and South Passage in 1998 due to the deep waterway regulation work. The river regime has now entered a good period, with a diversion angle of 65°.

    On the east of the main channel of the South Channel, due to the shoreline deflecting southward, the watercourse becomes wide, and the main stream discharges to the 10-m isobath. The boundary layer may be separated, and an orchid vortex forms. The bed sand is brought downward by lateral circulation through the South Passage. The rotary current is intensive near the area of the Jigu Reef outside of the North Passage, with the circulation ratio of the M2 equinoctial tide being 0.7. Because of the lag of aggradations outside the North Passage, the vertical section of the North Passage maintains a larger gradient, and the flow discharge is expedited. This can be considered as a good boundary condition of the river regime.

    2.3 Patterns of changing flow and sediment flux

    Shao et al. (2011) carried out in situ measurements of the sediment settling velocity near the Baimao Shoal in the Yangtze Estuary. Using a three-dimensional baroclinic model, the fine silt particle pathline of sediment dredging in the Yangtze Estuary was studied (Xie et al. 2010).

    Based on linear regression analysis of the hydrological data at the Datong Station from 1950 to 2006, the discharge from the Yangtze River, including the annual discharge and monthly discharge, was basically stable (Shen 2006; Li 2007). According to the linear regression analysis, the sediment concentration in the downstream of the Yangtze River decreased gradually. The average sediment concentration in the flood season is obviously larger than that in the dry season. The variation coefficient of the maximum sediment concentration in the flood season relative to that in the dry season is 0.393, and that of average sediment concentration is 0.311. The annual sediment transport and sediment concentration have evidently been decreasing with the same trend since 1984 (Fig. 6).

    Fig. 6Linear regression analysis of annual average sediment concentration at Datong Station

    Like the characteristic flow in flood and dry seasons, the annual extreme flow as well as its cumulative frequency can be used to calculate the tidal wave and abstraction volume in lower reaches of the Yangtze River. Then locationships of tidal limit and tidal current limit can be determined, and their relations with the discharge of extreme flow can be obtained. The effect of abstraction volume on the locations of tidal limit and tidal current limit can be shown as well.

    The following was found in Li (2007): (1) With a flood volume of 92 600 m3/s in 1954, the tidal limit of the spring tide appeared near Wuhu City, 395.3 km upstream from Xuliujing, and the tidal current limit appeared near the Qigan River and the Chengtong Reach, 22 km upstream from Xuliujing. When an astronomical tide occurred, the tidal wave traced far upstream even if the discharge was very large, which could threat the warning level and flood embankment. (2) With the flood volume of 4 620 m3/s in 1979, the tidal limit of the spring tide was upward to Anqing City, about 600 km upstream from Xuliujing. The tidal current limit was near the Bamao Hill at Heishazhou, about 450 km upstream from Xuliujing. With the discharge decreasing in the dry reason, the location of the tidal current limit traced upstream. Thus, the duration and intensity of saltwater intrusion increased evidently. At the same time, the to-and-fro current was adverse to sewage discharge and self-purification. (3) Compared with the tidal limit, the tidal current limit has a higher sensitivity. The moving distance of the tidal current limit was larger than that of the tidal limit under the same variation of the discharge. (4) The locations of tidal limit and tidal current limit are to some extent affected by the abstraction volume in the lower reach of the Yangtze River when the Tongcheng Reach is mainly controlled by tides.

    2.4 Effects of engineering projects on hydrodynamic characteristics

    The effects of the Three Gorges Project (TGP) on the hydrodynamic condition, river regime, waterway stability, and estuary morphology were analyzed based on the sediment and shipping technique research for TGP (Yan 1991; Chen et al. 2008). The following conclusions were drawn: (1) The discharge amplitude diminishes because of the regulation of annualreservoir discharge. After the construction of TGP, the range of the stagnation zone is shortened during the period of reservoir water storage with the water level varying from 150 m to 175 m in October of dry years. (2) A decrease of sediment transport from the stream is advantageous to the channel stability. (3) In October of dry years, the decrease of discharge caused by reservoir regulation enhances the inverse flow from the North Branch into the South Branch, but from January to May the increase of discharge eases up this inverse flow.

    The following can be seen from the calculated results of the regulation work of the deep-water channel in the Yangtze Estuary: (1) The change of bed resistance and tidal wave distortion may cause the changes of the average tidal level, high tidal level, and low tidal level. At the diversion point of the North Passage, the average level has run up, the high level has fallen down, and the low level has run up. While for the South Passage, the average level has fallen down, the high level has run up, and the low level has fallen down. (2) Affected by jetties and groins, the discharges of the North Passage and South Passage have redistributed with some ebb currents transferring from the North Passage to South Passage and the diversion ratio decreasing. (3) The regulation work of the deep-water channel supports channel maintenance.

    Based on the research on saltwater intrusion in the North Branch and the construction of water source regions using the numerical model of Hangzhou Bay, five plans were investigated with respect to the narrowing of the North Branch, construction of a tidal barrier in the North Branch, and reclamation of the Qingcaosha Reservoir on the South Branch. It can be found that (1) the narrowing of the North Branch has little effect on the South Branch and South Channel; (2) the reclamation of the Qingcaosha Reservoir induces the ebb current to increase at the upstream outlet of the North Branch, and causes the flood current to decrease at the downstream outlet of the North Branch; and (3) the plans for the narrowing of the North Branch and reclamation of the Qingcaosha Reservoir are preferable to that for the construction of the tidal barrier in the North Branch.

    Delft3D was used to establish a two-dimensional tidal and salinity numerical model of both the Yangtze Estuary and Hangzhou Bay, with a model domain of 330 km from east to west and 340 km from south to north. The effects of the reclamation project and submerged dike and reservoir construction on hydrodynamics in the dry season were studied using this model, in which the discharge, flow velocity, flow direction, and water level were considered. It can be seen that the engineering works will have effects on the flood-ebb current in the Yangtze Estuary, with a decrease of 17.4% of the flood-ebb current in the dry season and a decrease of 18.8% of the flood-ebb current in the flood season. The diversion ratio will increase in the North Passage, whereas the flood-ebb current of other river sections will diminish. The tidal influx will be reduced under the effects of the narrowed North Passage, with an increase of the flood-ebb current and a decrease of the tidal current at the outlet. The estuary area will be more sensitive to the increasing water, whereas the area inside the estuarywill be less sensitive. The characteristics of water and sediment exchange between the Yangtze Estuary and Hangzhou Bay were also studied by Kong et al. (2007).

    3 Saltwater intrusion

    3.1 Motion laws of back-flowing saltwater from North Branch to South Branch

    Supported by the major program of the National Natural Science Foundation of China (Grant No. 50339010), for the purpose of understanding the motion patterns of saltwater back flowing to the South Branch and the Xuliujing Reach, a tracking test for the back-flowing saltwater was conducted from March 12, 2005 to March 14, 2005 (Yang 2006). Some motion characteristics of the back-flowing saltwater were also studied. The moving time and range of the salt wedge were comparatively small. Being mixed well, the back-flowing saltwater rushed downstream along with the flood-ebb current. Owing to a relatively little difference in salinity, the surface saltwater and bottom saltwater mixed completely at the Qiyakou cross-section, with a distance of 10 km to the diversion point of the North Branch and South Branch. The water source region of Shanghai City lies in the area of the south bank in the South Branch, and is under the disadvantageous effects of the back-flowing saltwater.

    3.2 Numerical simulation of saltwater intrusion in Yangtze Estuary

    In the early 1980s, research showed that the mixture of saltwater and freshwater in the Yangtze Estuary was basically mild. A one-dimensional salinity numerical model was established using the diffusion equation and conservation equation, and the salinity longitudinal distribution and saltwater intrusion range were calculated. In recent years, the unstable longitudinal analytic solution for the one-dimensional saltwater intrusion has been derived. The field data in the deep-water channel have also been analyzed.

    Since the 1990s, a two-dimensional vertical numerical model and a two-dimensional horizontal numerical model have been developed to study the saltwater intrusion in the Yangtze Estuary (Wang 1989; Wang and Zhu 1991). With the volume control method and the power function in the discrete format, the tidal level, flow velocity, and salinity have been validated. In recent years, supported by ELCIRC and Delft3D, a two-dimensional numerical model of saltwater intrusion has been established, and a numerical model combining a one-dimensional tidal numerical model covering the region from the Datong Hydrological Station to the Yangtze Estuary and a two-dimensional numerical model of the Hangzhou Bay has been established as well (Yang 2006).

    Since the late 1990s, for simulating freshwater and saltwater mixing in estuaries, a three-dimensional nonlinear baroclinic numerical model has been developed, in which the gradient of horizontal pressure contains the gradient of barotropic pressure arising from the gradient of tidal level and the gradient of baroclinic pressure caused by the gradient of salinity (Zhu et al. 2000; Zheng et al. 2002; Song et al. 2008).

    Based on the mode-splitting technique, through the σ-coordinate transformation, the three-dimensional motion can be divided into an external mode and an internal mode. An improved double-sweep-implicit (DSI) method was employed in the external mode. The Eulerian-Lagrangian method was employed to deduce both the momentum equations of tidal motion and the equation of saltwater diffusion so as to improve the computational stability and accuracy. Methods to provide the boundary conditions and the initial conditions were proposed. The calculated salinity fields were presented, as shown in Fig. 7, where the dotted and solid lines represent the depth and salinity contours, respectively. Computational results show that the salinity distribution has a partial mixing pattern, and that the model is suitable for the simulation of freshwater and saltwater mixing in the Yangtze Estuary.

    Fig. 7Salinity distribution in flood slack tide in South Passage and North Passage in September 1996

    3.3 Effects of engineering projects on saltwater intrusion

    Based on the National Key Scientific and Technological Project of the Seventh Five-YearPlan, a monographic study of sediment problems of TGP and their effects on navigation channels and saltwater intrusion induced by TGP were analyzed. In October of both average years and dry years, the average discharge of the Three Gorges Reservoir was 33 000 to 41 400 m3/s, and the water storage of the Three Gorges Reservoir had little effect on saltwater intrusion in the Yangtze Estuary. It had a distinct effect only in very dry years. However, the saltwater intrusion eased up when the freshwater discharge increased from January to May (Li et al. 2005).

    In the National Key Scientific and Technological Project of the Eighth Five-Year Plan, a study on the evolution characteristics of sand bars and the deep-water channel regulation schemes in the Yangtze Estuary, a two-dimensional numerical model for lateral-vertical salinity distribution was employed to research the effects of the deep-water channel regulation work on saltwater intrusion. It was found that 0.5% saltwater near the bed traced upward in the dry season, but the regulation work had no effect on the salinity at Wusong. A three-dimensional salinity numerical model was developed to study the influence of the regulation work on saltwater intrusion after stages I, II, and III from 1996 to 2000, and during the future phase. With the field data measured six times at seven survey points, the stability of the numerical model was validated.

    With the one-dimensional and two-dimensional salinity numerical models, the effects of some hydro-junctions on trunk streams, such as the South-to-North Water Diversion Eastern Route Project and TGP, were studied as well. From the research, it was found that the intensity of saltwater back flowing to the South Branch grew with the decrease of discharge. The discharge under the combined running of the above two projects was between those under the condition that the two projects ran separately, so was the intensity of saltwater intrusion. With a discharge of 25 000 m3/s at the Datong Station, the salinity of the South Branch did not exceed 0.2%, which was less than half the salinity there in the above three cases of project running.

    The effects on saltwater intrusion, induced by regulation work in the Yangtze Estuary, were studied with a two-dimensional salinity numerical model. The result showed that enclosing land can reduce the intensity of saltwater intrusion but has little effect on the saltwater back flowing to the South Branch. Other combined plans recommended can reduce the saltwater back flowing, including the regulation work of the South Branch combined with the work of North Branch narrowing, or the regulation work of the South Branch combined with the work of North Branch narrowing and the construction of a tidal barrier. With these projects, salinity around the Chenhang Reservoir will be reduced, while that near the Qingcaosha Reservoir region will remain basically unchanged.

    In order to forecast the salinity variations in the Qingcaosha Reservoir region, a hydrodynamic and saltwater intrusion model covering the whole Yangtze Estuary and Hangzhou Bay was developed with the TELEMAC model (Tong et al. 2010). According to thedifferent characteristic discharges corresponding to different cumulative frequencies, the salinity distribution in the Yangtze Estuary was calculated under different flow rates into the sea. The simulation results showed that the salinity values would decrease with the increase of the discharge. When the upstream discharge is below the average value of 16 700 m3/s in the dry season, the salinity response to the runoff is very sensitive.

    4 Wetland evolution

    4.1 Database and methodology

    Estuarine wetlands exist as transition zones among marine areas, land, and freshwater areas. The zones are characterized by coupling interactions among marine, atmosphere, biology, geology, and human activities, and represent the most dynamic region in the nearshore area. The complex processes occurring in the estuarine wetlands are of fundamental importance to climate change, water conservation, flood control, beach protection, land reclamation, biological diversity, and the ecological balance. The total area of wetlands in the Yangtze Estuary is about 3 052 km2, consisting of 2 506 km2of nearshore wetlands, 478 km2of permanent river wetlands, and 68 km2of permanent freshwater lake wetlands. The sustainable development of wetlands in the Yangtze Estuary is essential for engineering projects and economic growth in this region. Fig. 8 shows the spatial distribution of wetlands. Recent data show that the deposition of most wetlands is slowing down, and some areas are now suffering erosion.

    Fig. 8Spatial distribution of wetlands in Yangtze Estuary

    The remote sensing technique was applied in the Yangtze Estuary. The topography database was implemented using professional software, including MAPGIS, ArcGIS, and Oracle. Data sets are generally classified into two categories, referred to as the digital elevation data and the remote sensing images. These two types of data were combined to establish a time series of wetland topography in the Yangtze Estuary (Zheng et al. 2010). Using underwater data measured in different years, a digital elevation model (DEM) was developed. In addition, the Landsat ETM color composite images were used to determine alongshore landscape and channel regime. To obtain the best match with DEM, the WGS84 coordinate system and the enhanced ERDAS image format with a ground resolution of 30 m were used. The superposition of the remote sensing images and DEM was implemented on the basis of the raster data format, and the superposition of the buoy data and the channel line was based on the vector data format.

    4.2 Morphological processes of wetlands

    A numerical morphodynamic model, TIMOR3, which was coupled with a hydrodynamic model and a wave model was applied to simulation of the long-term morphological response to the water and sediment changes in a large area in the Yangtze Estuary. A detailed investigation was made of the South Branch, where the deep-water channel regulation project is under construction (Zhou et al. 2009).

    The wetland evolution over a time scale of decades was analyzed (Yang et al. 2011). Results show that on the east beach of the Chongming Island, the sections along the shoreline from the Baozhen Port to Xijia Port and from the Beisixiao Port to Beiliuxiao Port were stable, while the section from the Xijia Port to the Tuanjie Shoal suffered from erosion. Rapid and slow accretions were found in the migratory bird natural reserve and the section from there to the Beiliuxiao Port, respectively. The wetland in the north shoals of the North Channel kept migrating seaward along with the lower reach from the Baozhen Port to the Liuxiao sandy ridge. The outer Tuanjie Shoal kept growing, and the siltation promotion project on the eastern Hengsha Shoal shaped the artificial shoreline. The mouth bar reach of the North Channel developed towards a formed channel. The wetland area of the Jiuduan Shoal has been continuously increasing since the 1950s, and from a forecasting point of view, it will merge with the Jiangya’nan Shoal to form a new wetland region between the South Passage and North Passage of the Yangtze Estuary. After floods in 1980, 1982, 1983, and 1988, and ten continuous wet years from 1990 to 1999, the length, width, and elevation of the Baimao Shoal enhanced, and the number of individual sandy bodies reduced, while the Central Shoal formed. Then, the shoal head receded slowly due to the washout of conflux. During the period from May 1992 to September 2002, the –5-m depth contour of the Baimao Shoal head moved backward by 2 720 m, with an annual backward distance of 272 m. The total volume of scoured sediment during that period was about 1.3 × 106m3, which led to negative impacts onthe lower reach of the South Branch of the Yangtze Estuary. The Central Shoal in the lower reach of the South Branch was unstable. The shoal and inlets transferred frequently and regularly. During the 1980s, the sand ridge of the Momao Shoal joined to the south beach of the Yangtze Estuary beyond the depth contour of –2 m and maintained its embossed sandy body. The sand tail of the Momao Shoal extended downstream, and its direction has not changed till now.

    5 Conclusions

    A review has been presented of the state-of-the-art research and engineering studies conducted by Hohai University relating to hydrodynamic and morphological processes in the Yangtze Estuary. Longitudinal, lateral, and horizontal flow circulations have been revealed on the basis of ADCP measurements. During the low slack tide, there are longitudinal circulation with the surface layer seaward and the bottom layer landward, lateral circulation with the surface layer southward and the bottom layer northward, and horizontal circulation with the north side landward and the south side seaward. The hydrodynamic mechanism at diversion points and the laws of flow change and sediment flux in the Yangtze Estuary were demonstrated through long-term data analysis. It was shown that the diversion points rise and fall periodically, with the channels flourishing and withering alternatively, and the waterways silting periodically. A field survey was carried out to detect the saltwater back flowing from the North Branch to South Branch. Different numerical models of flow motion, sediment transport, and saltwater intrusion were developed to simulate complicated processes and to evaluate the effects of engineering projects (e.g., TGP and the deep-water channel regulation project). It was found that the large engineering projects have significant effects on the patterns of flow and sediment flux as well as saltwater intrusion. The morphological processes of wetlands over a time scale of decades were analyzed with an established database, which is considered important for the management and protection of wetlands and the continuous economic growth in the Yangtze Estuary. Future studies should include effects of varying river runoff and sediment load on estuarine dynamics, the micro-scale movement characteristics of fine sediment, the bio-geomorphological model system, the joint impacts of large engineering projects on the long-term morphological evolution, and an integrated framework of environmental impact assessment for construction projects in the Yangtze Estuary.

    Chen, X. Q., Yan, Y. X., Fu, R. S., Dou, X. P., and Zhang, E. F. 2008. Sediment transport from the Yangtze River, China, into the sea over the post-Three Gorge Dam period: A discussion. Quaternary International, 186(1), 55-64. [doi:10.1016/j.quaint.2007.10.003]

    Kong, J., Song, Z. Y., Xia, Y. F., and Zhang, W. S. 2007. Characteristics of water and sediment exchange between Yangtze Estuary and Hangzhou Bay. China Ocean Engineering, 21(2), 255-266.

    Lei, Z. Y., Zhang, J. S., and Kong, J. 2009. Numerical simulation of water level under interaction between runoff and estuarine dynamics in tidal reach of the Yangtze River. China Ocean Engineering, 23(3), 543-551.

    Li, J. Y. 2007. Study on Flow and Sediment Features of Datong-Xuliujing Channel and River-bed Evolution. Ph. D. Dissertation. Nanjing: Hohai University. (in Chinese)

    Li, T. L., Li, Y. C., Gao, X. Y., and Wang, Y. G. 2005. Effects of regulation project on salinity intrusion in Yangtze Estuary. The Ocean Engineering, 23(3), 31-38. (in Chinese)

    Shao, Y. Y., Yan, Y. X., and Maa, J. P. Y. 2011. In situ measurements of settling velocity near Baimao Shoal in Changjiang Estuary. Journal of Hydraulic Engineering, 137(3), 372-380. [doi:10.1061/(ASCE) HY.1943-7900.0000312]

    Shen, H. Y. 2006. Studies on the Impact of the Water Abstraction along the Long River Reach of Tong Cheng. M. E. Dissertation. Nanjing: Hohai University. (in Chinese)

    Song, Z. Y., Huang, X. J., Zhang, H. G., Chen, X. Q., and Kong, J. 2008. One-dimensional unsteady analytical solution of salinity intrusion in estuaries. China Ocean Engineering, 22(1), 113-122.

    Tong, C. F. 2005. Study on Flow and Sediment Movement of Bifurcation Area and Related Application of 3-D Numerical Flow Model. Ph. D. Dissertation. Nanjing: Hohai University. (in Chinese)

    Tong, C. F., Zheng, J. H., Zhang, C., and Claude, G. 2010. Salinity response to the runoff from Yangtze River Basin at Qingcaosha Reservoir area in Yangtze Estuary. Proceedings of the 29th International Conference on Ocean, Offshore and Arctic Engineering: Ocean Engineering, 1-9. Shanghai: ASME. [doi: 10.1115/OMAE2010-20422]

    Wang, Y. G. 1989. 2-D Vertical Numerical Calculation on Saltwater Intrusion in Estuary. M. E. Dissertation. Nanjing: Hohai University. (in Chinese)

    Wang, Y. G., and Zhu, L. Z. 1991. Vertical two-dimensional numerical model of salt intrusion in tidal estuaries. Journal of Hohai University (Natural Sciences), 19(4), 1-8. (in Chinese)

    Wu, D. A., and Yan, Y. X. 2010. Numerical simulation of the transport and diffusion of dissolved pollutants in the Changjiang (Yangtze) River estuary. Chinese Journal of Oceanology and Limnology, 28(3), 649-657. [doi:10.1007/s00343-010-9081-7]

    Xie, R., Wu, D. A., Yan, Y. X., and Zhou, H. 2010. Fine silt particle pathline of dredging sediment in the Yangtze River deepwater navigation channel based on EFDC model. Journal of Hydrodynamics, 22(6), 760-772. [doi:10.1016/S1001-6058(09)60114-1]

    Yan, Y. X. 1991. Numerical modeling of wave and current interaction in the tidal inlet area. China Ocean Engineering, 5(1), 65-74.

    Yan, Y. X., Gao, J., Mao, L. H., and Zheng, J. H. 2000. Calculation of diversion ratio of the North Channel in Yangtze Estuary. China Ocean Engineering, 14(4), 525-532.

    Yan, Y. X., Gao, J., Song, Z. Y., and Zhu, Y. L. 2001a. Calculation method for stream passing around the Jiuduan sandbank in the Yangtze River mouth. Journal of Hydraulic Engineering, 32(4), 79-84. (in Chinese)

    Yan, Y. X., Gao, J., Zhu, Y. L., Mao, L. H., and Zheng, J. H. 2001b. Mechanism of excellent regime of North Passage in Yangtze Estuary. Hydro-Science and Engineering, (3), 8-12. (in Chinese)

    Yan, Y. X., Gao, J., Zhu, Y. L., and Zheng, J. H. 2001c. Preliminary study on relationship between deepwater channel regulation and riverbed evolution. Journal of Hohai University (Natural Sciences), 29(5), 7-12. (in Chinese)

    Yan, Y. X., Gao, J., Zheng, J. H., and Tong, C. F. 2002. Hydrodynamic conditions for sediment movement in South Channel of Yangtze Estuary. Journal of Hohai University (Natural Sciences), 30(5), 1-6. (in Chinese)

    Yan, Y. X., Ge, L., and Gao, J. 2003. Application of the minimum energy dissipation rate theory in the braided river. Journal of Hydrodynamics, 18(6), 692-697. (in Chinese)

    Yan, Y. X., Tao, A. F., Yu, D. S., and Yang, J. Y. 2007. Analysis on the circulation of Yangtze Estuary based on ADCP measurements. China Ocean Engineering, 21(3), 485-494.

    Yang, J. Y. 2006. ELCIRC Model Applied in Yangtze Estuary. M. E. Dissertation. Nanjing: Hohai University. (in Chinese)

    Yang, Q., Zheng, J. H., Ju, Y., and Peng, C. 2011. Analysis on the recent wetland evolution in the north branch of Yangtze Estuary. Proceedings of the 8th Conference on Sediment Theory in China, 59-66. Nanjing: Hohai University Press. (in Chinese)

    Yu, D. S. 2005. The Flow and Sediment Movement Analysis Based on ADCP with 3D Numerical Flow Model. Ph. D. Dissertation. Nanjing: Hohai University. (in Chinese)

    Zheng, J. H., Yan, Y. X., and Zhu, Y. L. 2002. Three dimensional baroclinic numerical model for simulating fresh and salt water mixing in Yangtze Estuary. China Ocean Engineering, 16(2), 227-238.

    Zheng, J. H., Yan, Y. X., Yun, C. X., Tong, C. F., and Han, Z. 2010. Long-term morphological evolution of wetland in Yangtze River Estuary. Proceedings (Abstract) of the 11th International Symposium on River Sedimentation: Sedimentation and Sustainable Use of River Systems, 125. Stellenbosch: Sun Media.

    Zhong, X. C. 1985. Fluvial process of the fork and braided channel and its effects on each other in the South-North Channel of Yangtze Estuary. Acta Geographica Sinica, 40(1), 51-59. (in Chinese)

    Zhou, X. Y., Zanke, U., and Yan, Y. X. 2009. Morphodynamic simulation at South Branch of Yangtze Estuary. Proceedings of the 28th International Conference on Ocean, Offshore and Arctic Engineering, 299-305. Hawaii: ASME. [doi:10.1115/OMAE2009-79418]

    Zhu, Y. L., Zheng, J. H., Mao, L. H., and Yan, Y. X. 2000.Three-dimensional nonlinear numerical model with inclined pressure for saltwater intrusion at the Yangtze River Estuary. Journal of Hydrodynamics, Ser. B, 12(1), 57-66.

    (Edited by Ye SHI)

    This work was supported by the National Basic Research Program of China (973 Program, Grant No. 2010CB429002), the Fundamental Research Funds for the Central Universities (Grant No. 2012B06514), the Special Fund of the State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering at Hohai University (Grant No. 2009585812), and the Qing Lan Project of Jiangsu Province.

    *Corresponding author (e-mail: jhzheng@hhu.edu.cn)

    Received Mar. 6, 2012; accepted Aug. 17, 2012

    九九久久精品国产亚洲av麻豆| 亚洲国产色片| 大码成人一级视频| 最新的欧美精品一区二区| 亚洲经典国产精华液单| 美女视频免费永久观看网站| 久久6这里有精品| 亚洲一级一片aⅴ在线观看| 亚洲欧美精品专区久久| 国产高清不卡午夜福利| 在线观看一区二区三区激情| 精品人妻熟女av久视频| 亚洲精品国产av成人精品| 日韩强制内射视频| 久久久久精品久久久久真实原创| 高清黄色对白视频在线免费看 | 国产有黄有色有爽视频| 久久久a久久爽久久v久久| 妹子高潮喷水视频| 少妇人妻精品综合一区二区| 大香蕉97超碰在线| 久久韩国三级中文字幕| 激情五月婷婷亚洲| 亚洲精品乱久久久久久| 成人国产麻豆网| 边亲边吃奶的免费视频| 亚洲中文av在线| 两个人免费观看高清视频 | 成人毛片a级毛片在线播放| 国产av码专区亚洲av| 久久精品国产亚洲av涩爱| 国产91av在线免费观看| 国产精品久久久久久精品电影小说| 99九九线精品视频在线观看视频| 一级毛片 在线播放| 日韩精品有码人妻一区| 成人黄色视频免费在线看| 边亲边吃奶的免费视频| 免费大片黄手机在线观看| 欧美变态另类bdsm刘玥| 亚洲美女搞黄在线观看| 久久久久国产精品人妻一区二区| 国产精品熟女久久久久浪| 亚洲成人手机| 综合色丁香网| 亚洲,一卡二卡三卡| av又黄又爽大尺度在线免费看| 中文欧美无线码| 精品人妻偷拍中文字幕| 久久国产精品大桥未久av | 肉色欧美久久久久久久蜜桃| 日韩一区二区三区影片| 亚洲av免费高清在线观看| 天堂8中文在线网| 少妇的逼好多水| 国产午夜精品久久久久久一区二区三区| 九九在线视频观看精品| 久久av网站| 夜夜爽夜夜爽视频| 久热久热在线精品观看| 日韩亚洲欧美综合| 91在线精品国自产拍蜜月| 狂野欧美激情性xxxx在线观看| av又黄又爽大尺度在线免费看| 最新中文字幕久久久久| 国产精品成人在线| 亚洲精品视频女| 伊人亚洲综合成人网| 高清午夜精品一区二区三区| 成年人免费黄色播放视频 | 色网站视频免费| 麻豆成人av视频| 久久亚洲国产成人精品v| av视频免费观看在线观看| 一级爰片在线观看| 亚洲成人一二三区av| 欧美 亚洲 国产 日韩一| 一级a做视频免费观看| 美女主播在线视频| 久久久午夜欧美精品| 寂寞人妻少妇视频99o| 久久久久久久久久人人人人人人| 日本黄色日本黄色录像| 日本免费在线观看一区| 麻豆成人av视频| 99精国产麻豆久久婷婷| 91久久精品电影网| 国产午夜精品久久久久久一区二区三区| 亚洲av成人精品一二三区| 日韩伦理黄色片| 国产成人91sexporn| 亚洲成人手机| 9色porny在线观看| 日韩av免费高清视频| 一级毛片 在线播放| 亚洲精品日韩av片在线观看| a级毛片在线看网站| 欧美精品一区二区大全| 国产精品久久久久久久电影| 岛国毛片在线播放| 国产精品人妻久久久久久| 色婷婷久久久亚洲欧美| 美女福利国产在线| 日本av手机在线免费观看| 丰满乱子伦码专区| 亚洲精品色激情综合| 国内揄拍国产精品人妻在线| 观看av在线不卡| 一区二区三区精品91| 日韩人妻高清精品专区| 国产成人精品久久久久久| 亚洲成人av在线免费| www.av在线官网国产| 欧美人与善性xxx| 97超碰精品成人国产| 国产av一区二区精品久久| 亚洲国产毛片av蜜桃av| 国产乱人偷精品视频| 精品99又大又爽又粗少妇毛片| 亚洲精品成人av观看孕妇| 国产老妇伦熟女老妇高清| 婷婷色综合www| 亚洲婷婷狠狠爱综合网| 久久 成人 亚洲| 日本av手机在线免费观看| 一级二级三级毛片免费看| 日本猛色少妇xxxxx猛交久久| 我的女老师完整版在线观看| 又大又黄又爽视频免费| 蜜臀久久99精品久久宅男| 男的添女的下面高潮视频| 亚洲欧美精品自产自拍| 99精国产麻豆久久婷婷| 色94色欧美一区二区| 国产在线视频一区二区| 国产精品99久久久久久久久| 亚洲精华国产精华液的使用体验| 精品视频人人做人人爽| 夜夜爽夜夜爽视频| 国产成人精品一,二区| av黄色大香蕉| 涩涩av久久男人的天堂| 老熟女久久久| 亚洲精品乱码久久久v下载方式| 午夜精品国产一区二区电影| 国产69精品久久久久777片| 久久久久精品性色| 国产老妇伦熟女老妇高清| 久久热精品热| 日本免费在线观看一区| 久久鲁丝午夜福利片| 国产探花极品一区二区| 亚洲无线观看免费| 十八禁网站网址无遮挡 | 欧美日韩视频精品一区| 亚洲av中文av极速乱| 国产精品国产三级专区第一集| 日韩一本色道免费dvd| 天天躁夜夜躁狠狠久久av| 成人毛片a级毛片在线播放| 交换朋友夫妻互换小说| 亚洲精品亚洲一区二区| 亚洲精品日韩在线中文字幕| 亚洲av中文av极速乱| 国产黄频视频在线观看| 在线观看免费日韩欧美大片 | 久久久久久久久久久久大奶| 国产 一区精品| 十分钟在线观看高清视频www | 亚洲情色 制服丝袜| 国产av国产精品国产| 91精品伊人久久大香线蕉| 欧美日韩国产mv在线观看视频| 性色avwww在线观看| 免费黄色在线免费观看| 精品人妻一区二区三区麻豆| 少妇被粗大的猛进出69影院 | 婷婷色麻豆天堂久久| 一区在线观看完整版| 啦啦啦在线观看免费高清www| 国产日韩一区二区三区精品不卡 | www.色视频.com| 夫妻午夜视频| 青青草视频在线视频观看| 高清不卡的av网站| 久久97久久精品| 成人无遮挡网站| 少妇人妻精品综合一区二区| 91精品国产九色| 欧美丝袜亚洲另类| 国产91av在线免费观看| av在线老鸭窝| 蜜桃久久精品国产亚洲av| 亚洲不卡免费看| 少妇的逼水好多| 看免费成人av毛片| 欧美三级亚洲精品| 男女啪啪激烈高潮av片| 日韩不卡一区二区三区视频在线| 中文精品一卡2卡3卡4更新| 一级毛片久久久久久久久女| 人妻人人澡人人爽人人| 亚洲欧美一区二区三区国产| 国产爽快片一区二区三区| 亚洲丝袜综合中文字幕| 夫妻性生交免费视频一级片| 91午夜精品亚洲一区二区三区| freevideosex欧美| 亚洲av.av天堂| 亚洲内射少妇av| 天堂8中文在线网| 国产在视频线精品| 国产中年淑女户外野战色| 国产精品久久久久久久久免| 一级毛片aaaaaa免费看小| www.av在线官网国产| 免费高清在线观看视频在线观看| 国产精品国产三级国产av玫瑰| 久久久精品94久久精品| 国产免费一区二区三区四区乱码| 国产中年淑女户外野战色| 成年av动漫网址| 久久精品熟女亚洲av麻豆精品| 日韩一区二区视频免费看| 性高湖久久久久久久久免费观看| 高清黄色对白视频在线免费看 | 六月丁香七月| 一级黄片播放器| 99久久人妻综合| 午夜免费男女啪啪视频观看| 岛国毛片在线播放| 婷婷色麻豆天堂久久| 99热网站在线观看| 亚洲情色 制服丝袜| 一区二区三区乱码不卡18| 少妇熟女欧美另类| 99热这里只有是精品50| 蜜桃在线观看..| 国产成人精品无人区| a级毛片在线看网站| 99久久精品热视频| 中文字幕av电影在线播放| 精品人妻熟女av久视频| 一个人免费看片子| 欧美激情极品国产一区二区三区 | 最后的刺客免费高清国语| 人妻系列 视频| 中文字幕人妻熟人妻熟丝袜美| 亚洲性久久影院| 永久网站在线| 日本午夜av视频| 一级毛片久久久久久久久女| 熟女人妻精品中文字幕| 日本黄色片子视频| 黄色欧美视频在线观看| 大香蕉久久网| 欧美亚洲 丝袜 人妻 在线| 777米奇影视久久| 久久综合国产亚洲精品| av不卡在线播放| 女的被弄到高潮叫床怎么办| 国产在线男女| 99视频精品全部免费 在线| 国产午夜精品久久久久久一区二区三区| 免费久久久久久久精品成人欧美视频 | 日韩在线高清观看一区二区三区| av福利片在线观看| 亚洲av中文av极速乱| 各种免费的搞黄视频| videossex国产| av黄色大香蕉| 久久久久久久久久人人人人人人| 街头女战士在线观看网站| 国产亚洲午夜精品一区二区久久| 免费观看无遮挡的男女| 久久人人爽人人片av| 丝袜喷水一区| 人妻制服诱惑在线中文字幕| 日韩av不卡免费在线播放| 色吧在线观看| 国产av国产精品国产| 免费黄色在线免费观看| 丝袜喷水一区| 美女大奶头黄色视频| 伦理电影免费视频| 啦啦啦啦在线视频资源| 日本黄色日本黄色录像| 国产伦精品一区二区三区视频9| av卡一久久| 高清毛片免费看| 最近中文字幕2019免费版| 97超视频在线观看视频| 日韩免费高清中文字幕av| 成人综合一区亚洲| 国产成人精品无人区| av天堂久久9| 简卡轻食公司| 黄色欧美视频在线观看| 国产成人freesex在线| 婷婷色麻豆天堂久久| 久久久国产一区二区| 亚洲综合精品二区| 久久久久久久国产电影| 久久av网站| 国产永久视频网站| 最近最新中文字幕免费大全7| 91午夜精品亚洲一区二区三区| 亚洲欧美成人综合另类久久久| 国产日韩一区二区三区精品不卡 | 一区二区三区精品91| 免费不卡的大黄色大毛片视频在线观看| 国产高清有码在线观看视频| 最近的中文字幕免费完整| av专区在线播放| 夜夜看夜夜爽夜夜摸| 久久精品久久久久久噜噜老黄| 成人影院久久| 高清欧美精品videossex| 熟妇人妻不卡中文字幕| 人人妻人人爽人人添夜夜欢视频 | 国产伦精品一区二区三区视频9| 少妇裸体淫交视频免费看高清| 国产精品久久久久久精品电影小说| 又爽又黄a免费视频| 亚洲欧美日韩东京热| 亚洲精品日韩在线中文字幕| 多毛熟女@视频| 亚洲国产成人一精品久久久| 搡女人真爽免费视频火全软件| 搡老乐熟女国产| 国产无遮挡羞羞视频在线观看| 欧美日韩综合久久久久久| 日本猛色少妇xxxxx猛交久久| 下体分泌物呈黄色| 免费人成在线观看视频色| 精品熟女少妇av免费看| 97超视频在线观看视频| 亚洲精品日本国产第一区| 久久婷婷青草| 永久免费av网站大全| 国产亚洲欧美精品永久| 日韩亚洲欧美综合| 日本欧美视频一区| 国产熟女欧美一区二区| 色吧在线观看| 日韩免费高清中文字幕av| 精品亚洲乱码少妇综合久久| 国产成人freesex在线| 五月伊人婷婷丁香| 日韩成人av中文字幕在线观看| 少妇熟女欧美另类| 两个人免费观看高清视频 | 久久午夜综合久久蜜桃| 日韩 亚洲 欧美在线| 久久亚洲国产成人精品v| 国产真实伦视频高清在线观看| 天美传媒精品一区二区| av线在线观看网站| av.在线天堂| 激情五月婷婷亚洲| 免费看av在线观看网站| 欧美+日韩+精品| 欧美日韩视频精品一区| 亚洲精品国产av成人精品| 日本欧美视频一区| 亚洲激情五月婷婷啪啪| 成人18禁高潮啪啪吃奶动态图 | 91久久精品国产一区二区成人| 亚洲国产毛片av蜜桃av| 亚洲一区二区三区欧美精品| 国产av国产精品国产| 久久精品国产亚洲av天美| 色网站视频免费| 美女xxoo啪啪120秒动态图| 成人综合一区亚洲| 少妇人妻 视频| 热re99久久精品国产66热6| 三级国产精品欧美在线观看| 91精品伊人久久大香线蕉| 日韩在线高清观看一区二区三区| av又黄又爽大尺度在线免费看| 蜜桃久久精品国产亚洲av| 国产黄色免费在线视频| 春色校园在线视频观看| 精品人妻熟女毛片av久久网站| 一级毛片电影观看| 少妇精品久久久久久久| 99热国产这里只有精品6| 一级黄片播放器| 制服丝袜香蕉在线| 又大又黄又爽视频免费| 亚洲欧洲精品一区二区精品久久久 | 亚州av有码| 欧美xxxx性猛交bbbb| 青春草视频在线免费观看| 午夜福利网站1000一区二区三区| 免费av中文字幕在线| 国产乱来视频区| 亚洲久久久国产精品| 五月伊人婷婷丁香| 国产成人a∨麻豆精品| 日日啪夜夜爽| 三级国产精品欧美在线观看| 在线观看三级黄色| 高清欧美精品videossex| a级一级毛片免费在线观看| 免费人妻精品一区二区三区视频| 嫩草影院新地址| 精品99又大又爽又粗少妇毛片| 插逼视频在线观看| 日韩大片免费观看网站| 精品一区二区三区视频在线| 97超碰精品成人国产| 精华霜和精华液先用哪个| 热re99久久国产66热| 国产成人精品一,二区| 久久免费观看电影| 欧美日韩视频高清一区二区三区二| 晚上一个人看的免费电影| 成人亚洲欧美一区二区av| 女人精品久久久久毛片| 日本欧美国产在线视频| 80岁老熟妇乱子伦牲交| 街头女战士在线观看网站| 日日摸夜夜添夜夜爱| 中文在线观看免费www的网站| 国产熟女欧美一区二区| 久久久国产一区二区| 国产视频内射| 久久久久久久久久久免费av| 久久久久视频综合| 国产黄色视频一区二区在线观看| 日本免费在线观看一区| 五月玫瑰六月丁香| 三级国产精品片| 内射极品少妇av片p| 人人妻人人澡人人看| 春色校园在线视频观看| a级一级毛片免费在线观看| 亚洲精品日韩av片在线观看| 亚洲婷婷狠狠爱综合网| 精品久久久噜噜| 2022亚洲国产成人精品| 91久久精品国产一区二区成人| 大片免费播放器 马上看| 国产精品成人在线| 黄色毛片三级朝国网站 | 少妇熟女欧美另类| 蜜桃在线观看..| 日韩不卡一区二区三区视频在线| 性色avwww在线观看| 一级毛片aaaaaa免费看小| 国产精品无大码| 亚洲综合精品二区| 久久毛片免费看一区二区三区| 亚洲熟女精品中文字幕| 中文字幕精品免费在线观看视频 | 国产亚洲91精品色在线| 97精品久久久久久久久久精品| 岛国毛片在线播放| 欧美精品国产亚洲| 亚洲av综合色区一区| 国产熟女欧美一区二区| av天堂久久9| 老熟女久久久| 午夜福利网站1000一区二区三区| 嫩草影院新地址| 日本欧美视频一区| 一区在线观看完整版| 国产精品熟女久久久久浪| 伊人久久精品亚洲午夜| 九色成人免费人妻av| 国产欧美日韩一区二区三区在线 | 欧美三级亚洲精品| 91在线精品国自产拍蜜月| 99久久精品一区二区三区| 精品一区二区免费观看| 少妇 在线观看| 久久久久久久久久人人人人人人| 我要看日韩黄色一级片| 新久久久久国产一级毛片| 国产精品无大码| 亚洲国产精品成人久久小说| 性高湖久久久久久久久免费观看| 激情五月婷婷亚洲| 亚洲av成人精品一区久久| 99久久中文字幕三级久久日本| 又爽又黄a免费视频| 亚洲怡红院男人天堂| 99久久人妻综合| 人人妻人人澡人人看| 亚洲欧美日韩另类电影网站| 国产黄片美女视频| 99九九线精品视频在线观看视频| 欧美精品高潮呻吟av久久| 日韩中字成人| 99热全是精品| 成人毛片a级毛片在线播放| 王馨瑶露胸无遮挡在线观看| 日日摸夜夜添夜夜添av毛片| 爱豆传媒免费全集在线观看| 亚洲国产成人一精品久久久| 亚洲精品亚洲一区二区| 欧美 亚洲 国产 日韩一| 黄色一级大片看看| 噜噜噜噜噜久久久久久91| 乱码一卡2卡4卡精品| av免费观看日本| 免费在线观看成人毛片| 免费观看av网站的网址| 亚洲高清免费不卡视频| 亚洲无线观看免费| 有码 亚洲区| tube8黄色片| 日韩伦理黄色片| 在线观看www视频免费| 91久久精品电影网| 精品少妇内射三级| 91精品一卡2卡3卡4卡| 日韩av免费高清视频| av在线老鸭窝| 亚洲av.av天堂| 内地一区二区视频在线| 国产精品偷伦视频观看了| 欧美日韩av久久| 国产熟女欧美一区二区| 免费观看无遮挡的男女| 国产伦理片在线播放av一区| 大片电影免费在线观看免费| 精品亚洲成国产av| 自拍偷自拍亚洲精品老妇| 最近的中文字幕免费完整| 国产精品一区二区性色av| 精品久久久精品久久久| 免费大片黄手机在线观看| 亚洲国产精品一区三区| 国产伦精品一区二区三区四那| 99久久综合免费| 日日摸夜夜添夜夜添av毛片| 激情五月婷婷亚洲| 中文字幕久久专区| 日韩大片免费观看网站| 亚洲不卡免费看| 国产一区二区在线观看日韩| 亚洲人成网站在线播| 26uuu在线亚洲综合色| 丝袜在线中文字幕| 最新的欧美精品一区二区| a 毛片基地| 交换朋友夫妻互换小说| 春色校园在线视频观看| 大陆偷拍与自拍| 大片免费播放器 马上看| 少妇高潮的动态图| 亚洲国产精品一区三区| 纯流量卡能插随身wifi吗| 男女啪啪激烈高潮av片| 青青草视频在线视频观看| 国产色爽女视频免费观看| 最新中文字幕久久久久| 韩国高清视频一区二区三区| 边亲边吃奶的免费视频| 亚洲va在线va天堂va国产| 伦精品一区二区三区| 国产亚洲91精品色在线| 自线自在国产av| 亚洲精品日韩av片在线观看| 99热6这里只有精品| 国产免费一级a男人的天堂| 最近中文字幕高清免费大全6| h视频一区二区三区| 国产精品国产三级国产av玫瑰| 亚洲精品中文字幕在线视频 | 汤姆久久久久久久影院中文字幕| 精品一区二区三区视频在线| 插逼视频在线观看| 亚洲婷婷狠狠爱综合网| 岛国毛片在线播放| 久久久久久久精品精品| 国产精品人妻久久久久久| 国产亚洲欧美精品永久| 老女人水多毛片| 99热国产这里只有精品6| 国产精品久久久久久精品古装| 波野结衣二区三区在线| 色婷婷av一区二区三区视频| 少妇 在线观看| 永久免费av网站大全| 高清欧美精品videossex| 免费观看av网站的网址| 久久 成人 亚洲| 精品人妻熟女毛片av久久网站| 国产日韩一区二区三区精品不卡 | 国产精品偷伦视频观看了| 久久国内精品自在自线图片| 免费大片黄手机在线观看| 青青草视频在线视频观看| 少妇的逼水好多| 亚洲av福利一区| 亚洲不卡免费看| 亚洲美女视频黄频| 精品99又大又爽又粗少妇毛片| 最新中文字幕久久久久| 免费看av在线观看网站| 亚洲精品,欧美精品| 国产中年淑女户外野战色| 日韩av在线免费看完整版不卡| 乱系列少妇在线播放| 天天操日日干夜夜撸| 国产淫片久久久久久久久| 国产成人精品福利久久| 亚洲情色 制服丝袜| 亚洲精华国产精华液的使用体验| 亚洲欧洲日产国产| 日韩精品免费视频一区二区三区 |