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

    A Water Movement Study in Lianzhou Bay, Guangxi Province, China

    2014-05-02 05:41:56SUNTingAndreasMacranderandDavidKaiser
    Journal of Ocean University of China 2014年1期

    SUN Ting, Andreas Macrander, and David Kaiser

    1) College of Physical and Environmental Oceanography, Ocean University of China, Qingdao 266100, P. R. China

    2) Marine Research Institute, Skúlagata 4, 121 Reykjavík, Iceland

    3) Leibniz Centre for Tropical Marine Ecology, Fahrenheitstr, 6, D-28359 Bremen, Germany

    A Water Movement Study in Lianzhou Bay, Guangxi Province, China

    SUN Ting1),*, Andreas Macrander2), and David Kaiser3)

    1) College of Physical and Environmental Oceanography, Ocean University of China, Qingdao 266100, P. R. China

    2) Marine Research Institute, Skúlagata 4, 121 Reykjavík, Iceland

    3) Leibniz Centre for Tropical Marine Ecology, Fahrenheitstr, 6, D-28359 Bremen, Germany

    This study investigates the physical conditions (water depth, current speed, salinity, temperature) in Lianzhou Bay, a shallow coastal bay in southern China, during two expeditions in the dry and wet seasons of 2011. Based on these expedition data, basic hydrodynamic parameters like Brunt-V?is?l? Frequency, Richardson Number, Rossby radius, and Resonance Period are calculated. The results show that Lianzhou Bay is characterized by comparatively small quantity of freshwater input and weak stratification. Strong tides, which are spatially uniform within the bay, cause turbulent mixing. Residence time of the water is shorter in winter due to a stronger coastal current in that season. Consideration of the water movement may help to reduce the harmful ecological impact of aquaculture waste water discharge.

    water movement; water mixing process; Rossy radius; tidal resonance; residence time; Lianzhou Bay

    1 Introduction

    Water movement studies aim to characterize the physical oceanography of a given area regarding tidal regime, water mixing processes, water exchange processes,etc.Such studies have been carried out on different spatial scales,e.g., the Baltic Sea (Osińskiet al., 2010), the Gulf of Tonkin, which is also called the Beibu Gulf (Wuet al., 2008), the Tachibana Bay of Japan (Matsunoet al., 1999), Rio of Ferrol (deCastroet al., 2004), west coast of India (Unnikrishnan and Antony, 1990), and the estuary system of the Columbia River (Chawlaet al., 2008). In the management of aquaculture, sea water movement studies may be employed to minimize the ecological impact of wastewater. Moreover, it has been suggested that, an appropriate understanding of the water movement dynamics could be applied for any kind of anthropogenic activity in coastal areas (Páez-Osuna, 2001; Ulseset al., 2005; Matsunoet al., 1999).

    The site of this study–Lianzhou Bay–is one of the coastal regions of China which accommodate intense aquaculture activities. The study assesses the physical oceanography of Lianzhou Bay and its driving factors, usingin situdata (temperature, salinity, pH, dissolved oxygen, water depth, current speed,etc.) from several expeditions during August to October in 2010 and March to May in 2011.

    2 Methods

    2.1 Study Site

    Lianzhou Bay is a shallow coastal bay located on the northern coast of the Gulf of Tonkin, which is part of the South China Sea (Fig.1).

    The Nanliu River delta borders the Lianzhou Bay. In its vicinity are Beihai city in the east and Beibu Gulf in the south (Fig.1). Lianzhou Bay is among the five largest bays in Guangxi province, Southern China. The width of its mouth and the length of its coastal line are 25 and 79 km, respectively (Chen, 2003; Lvet al., 1995; Liang and Li, 2002). The bay is rather shallow (average depth 5 m; Chen, 1996), and has extensive tidal flats in its northern half. There is intensive oyster farming in the shallow subtidal area of the bay.

    Lianzhou Bay is typical for coastal bays with moderate freshwater input and large tides. Nanliu River is the largest fresh water source to the sea in Guangxi Province. It has three main branches: Nandong River in the west (which is the largest in terms of size and runoff), Nanxi River, and Nangan River in the east (which is the smallest). According to the data provided by the Changle Hydrological Station, the average fresh water discharge of Nanliu River from 1954–1982 reaches 168 m3 s?1. In general, the largest discharge occurs in summer during June to August (Gu and Wu, 2001; Jianget al., 2008). This runoff is small compared with estuaries of large riverslike,e.g., the Columbia River System which has a runoff between 2000 m3 s?1and 10000 m3 s?1(Chawlaet al., 2008). The sediments carried by the Nanliu River are the main material source for the Lianzhou Bay with the material transport direction from NE to SW (Gu and Wu, 2001).

    The bay is affected by strong tides. The mean tidal difference is 2.54 m, and the largest tidal difference reaches 5.36 m (Lai and Wei, 2003). Fig.2 shows a typical tidal time series spanning the time of the dry season expedition, demonstrating that in Lianzhou Bay diurnal tides are mostly dominant over semi-diurnal tides. Previous studies concluded that tidal currents dominate water movement in the bay, where the primary pattern are reversing tidal currents (Lai and Wei, 2003; Jianget al., 2008; Sun and Huang, 2001; Liang and Li, 2002). Sea water rushes toward N and NE in flood, and returns sea-ward to S and SW after reaching the head of the bay. The average tidal current velocity in Lianzhou Bay is larger on the ebb than on the flood tide, being around 104 cm s?1and 88 cm s?1, respectively (Chen, 1996; Chen and Shi, 1996).

    Fig.1 Map of Lianzhou Bay. The location is marked by the inset map. Local bathymetry data are obtained from a local bathymetry map in 1987 based on theoretical depth datum. Depths are referenced to mean sea level. Tidal flats are shaded.

    Fig.2 Predicted water level at Beihai city from 2011/4/1 to 2011/4/30 during the dry season expedition. Mean sea level is 255 cm.

    The ecology of Lianzhou Bay is affected by the economic activities of a population of about 320000 living within the catchment area of Nanliu River (Qiu and Lai, 2004). The land is mostly used for agriculture further upstream and for aquaculture closer to the coast. Local farmers discharge waste water from the aquaculture ponds into ditches, from where it flows into Nanliu River. The average annual runoff of nutrients (nitrogen, phosphorus, silicate) flushed from the river into Lianzhou Bay approaches 5900 t (Chen, 2001; Qiu and Lai, 2004). Consequently, the large load of nutrients and organic matter creates medium-polluted and eutrophicated conditions in Lianzhou Bay according to Sea water quality standard GB 3097-1997 and Environmental quality standard for surface water GB 3838-2002 (Chen, 2001). Algal (MicrocystisKuetz) blooms occurred in Lianzhou Bay in March 1995, and in February 2004 (Wei and He, 1998; Weiet al., 2004; Qiu and Lai, 2005).

    2.2 Data Collection

    Two expeditions have been carried out to collect in-situ data of physical parameters. To cover seasonal variability, expeditions have been done during wet season and dry season. The wet season extended from late August to early November in 2010. During this period, samples were taken once along each river tributary (not discussed in this paper), and twice on sections across the bay. During the dry season, an expedition was done during March and April in 2011 with samples taken once in each area. Sampling stations across the bay are shown in Fig.3. Typical parameters,i.e., temperature, salinity, pH, and dissolved oxygen were measured using a Hach HQ40d portable sensor and compatible probes (Hach Company, 2012). All measurements, except for DO, were made at both surface and bottom at each station of the bay and the three river branches. Water depths larger than 1 m were determined using a handheld echo sounder, while shallower depths were measured with a wooden stick and measuring tape. Current speed was measured with a simple and effective method recording the time for a handmade drift bottle to drift over a distance of 10 m along a moored boat or a fixed structure.

    In addition to sampling cruises, there were five fixed sampling sites to observe the 24-h tidal changes. Two were at the mouth of Nangan River during dry season (Fig.3C), and the rest were in a ‘Clam house’ built by local fishermen around 1 km away from the other 24-h sampling station, more out to the sea (Figs.3A, C). The serial number instead of the date in Table 1 will be used hereafter. Sampling was done every 2 h during the wet season (Station 1) and hourly during the dry season (Stations 2–5). Procedures and parameters are similar to sampling cruises but taken at the same location.

    Fig.3 Bay expeditions and 24 h sampling. In (A), (B), and (C), dark dots denote bay expedition cruises of 2010/08/ 30, 2010/10/05, and 2011/04/08, respectively. Numbers account for the direction of the cruise. Light squares denote the 24-h sampling stations for No.1 wet season, and Nos. 2 to 5 dry season.

    Table 1 Description of the 24-h samplings

    2.3 Data Processing

    Calculation of parameters such as Brunt-V?is?l? Frequency (BVF), Richardson Number (Ri), and water volume differences related to tides are based onin situdata of Lianzhou Bay. The potential density was calculated using the UNESCO equation of state (Fofonoff and Millard, 1983). Data were visualized using Ocean Data View (Schlitzer, 2012) and overlain with a modified map provided by colleagues from the Guangxi Mangrove Research Center (GMRC) in Beihai, Guangxi, China. Numerical modelling of the current dynamics in the bay is beyond the scope of this study, as no high-resolution dataset of the local bathymetry is available to set up a high- resolution ocean circulation model for Lianzhou bay. Nevertheless, it is possible to determine the general characteristics of the water movement in the bay by hydrodynamic key parameters based on the data of the expeditions.

    3 Results and Discussion

    3.1 Physico-Chemical Profiles of Lianzhou Bay

    Fig.4 shows the temperature and salinity distribution in Lianzhou Bay. Limited by the shallow water, the first wet season cruise on 2010/8/30 was conducted during a semidiurnal tidal flood and the first half of the ebb; the dry season cruise on 2011/4/8 was launched while a diurnal tidal flood started.

    Fig.4A gives a clear example of the seasonal variation of water temperature. The average temperate was 30℃during the wet season (A1) and 23℃ during the dry season (A2). No significant temperature gradient was observed between surface and bottom water. During wet season, water temperature is rather uniform across the bay with values up to 31.75℃. During the cooler dry season, temperatures decrease from the river mouth (approx. 25℃) to the deep bay area (approx. 21℃).

    Salinity within the bay is relatively high. During wet season, sea surface salinity ranges between 22 and 26, while at the bottom a salt wedge is found at the seaward edge with a salinity of 31 (Fig.4 B1). During dry season, fresh water at the surface does not spread as far, with a significant salinity gradient between the river mouths and the southern (outer) part of the bay, where the salinity increases to 31 (Fig.4 B2).

    3.2 Major Hydrodynamic Parameters

    3.2.1 Stratification/Brunt-V?is?l? frequency

    The Brunt-V?is?l? Frequency (BVF)Nis a measure of the stratification of a water body and is proportional to the stability of stratification (Houryet al., 1987; Osińskiet al., 2010).

    The Brunt-V?is?l? FrequencyNin a given depth is given by:

    In Lianzhou Bay,Nranges between 91 cph (cycles per hour) during wet season and 59 cph during dry season. These values are somewhat lower than in,e.g., the Weser Estuary (North Sea), which is under a stronger riverine influence (average runoff 323 m3s?1, Grabemann and Krause, 2001). There, a meanNof 127 cph was found (Macrander, 2009), which indicates a stronger stratification, especially during ebb tide when the fresh water flows out faster than the bottom salt water (Geyer, 1993). In another highly stratified water system, Baltic Sea,Nis found between 126 and 198 cph within thermo- and haloclines (Osińskiet al., 2010). Compared to these examples, stratification in Lianzhou Bay is weaker.

    Fig.5 showsNat the surface of the bay during wet season (2010/8/30) and dry season (2011/4/8). Towards the coast, salinity and density are lower, with the lowest salinities occurring close to the river mouths. Stratification in the rather shallow coastal areas is weak, too (smallNin,e.g., Station 3 of Fig.5A). In contrast, a higher saline bottom layer is found in the seaward Station 8 (Fig.5A), associated with a stronger stratification. Measurements at Stations 8–10 were taken during ebb tide; obviously, the lighter fresh water at the surface flows out faster, while close to the bottom, a salt wedge persists. In dry season, stratification is weaker also in the seaward part (Fig.5B), as the fresh river water does not spread that far. The data were taken during flood, hence, tidal current induces turbulent mixing which causes equally distributed density, consequently,Nat the bay mouth decreases to about 30 cph (Stations 1 and 2). The higherNin Stations 3–5 is caused mainly by higher surface water temperature rather than lower salinity. Measurements at these stations were taken during noontime, where the surface layer is additionally heated up by solar radiation, accounting for a larger vertical density gradient.

    Fig.4 Temperature (A) and salinity (B) distribution in Lianzhou Bay in 2010/8/30 (1) and 2011/4/8 (2).

    Fig.5 Distribution of N on the surface of the bay. (A) 2010/8/30, during wet season; (B) 2011/4/8, during dry season.

    3.2.2 Turbulent mixing/Richardson number

    A stable stratification with a fresh water surface layer and a more saline deeper layer can be removed by turbulent mixing, provided that the potential energy of the stratification is smaller than the kinetic energy which is associated with a sheared current. This ratio betweenEpotandEkinis defined as Richardson numberRi:

    whereNdenotes the BVF, and du/dzthe vertical gradient of the current speed. IfRi<0.25, turbulent mixing breaking up the stratification can be expected (Galperinet al., 2007). As current speed was measured only at the surface, we useusurface/H(with water depthH) as approximation for du/dz, assuming zero current speed right at the bottom.

    As current speed data under consideration are limited to the 24 h sampling,Riis accessed only for these timeseries. Turbulent mixing was prevalent during 66% of the observation time, whenRiwas smaller than 0.25, whereas during 34% of the time (mostly during slack tide),Riwas larger than 0.25. As an example, Fig.6 shows data from No. 5 24-h sampling. The potential density profile is largely determined by salinity (Figs.6A, 6B). The largest density difference between surface and bottom occurred around 14:00 on 2011/4/5 during the higher low water within the sampling period. The density difference implies maximumNat the same time (Fig.6C). The smallest measured current speed was 0.09 m s?1and normally occurred during slack tide at low water, while the highest speed generally occurred during late stage of flood tides (Fig.6D). Consequently,Riis largest during periods of largeNand low current speed (Fig.6E). While during the higher low tide, a stable stratification persisted for some time, current velocities were large enough to keepRibelow 0.25 at all other times. During the lower low tide the water depth was so small (less than 1 m), that the stratification was completely removed by turbulent mixing.

    Relying on the calculatedNandRi, we conclude that conditions in Lianzhou Bay are mostly dominated by turbulent mixing. This results in a vertically largely homogeneous distribution of salinity and density, while it weakens the baroclinic flow. Considering the shallowness, we conclude that currents in Lianzhou Bay are largely barotropic,i.e., in the same direction throughout the entire water column.

    Fig.6 24-h sampling No. 5: Temporal variation of salinity (A), potential density (B), N (C), current speed (D) and Ri (E). Tides are indicated by the measured water depth (grey bars).

    3.2.3 Earth rotation/Rossby radius

    On large scales, oceanic currents are deflected to theright (northern hemisphere) by the rotation of the earth. A characteristic length scale is given by the Rossby radius. As currents in Lianzhou bay are largely barotropic, the barotropic Rossby radiusLris assessed here:

    With the measured depth between 1 m and 5.6 m,Lrranges from 58 to 141 km. This is far larger than the size of the bay (24 km), hence, current dynamics in Lianzhou bay are not significantly affected by the rotation of the earth,i.e., currents will not be focused to the coast on the right-hand side.

    3.2.4 Tidal resonance

    Tides in bays can be amplified by interference between incoming and reflected tidal waves (e.g., in Bay of Fundy, Greenburg, 1979; in the Strait of Georgia, Forman, 2005), provided the resonance period of the bay is close to the tidal period.

    The shortest possible basin length L which produces resonances of tidal waves with period T (Leder and Orli?, 2004) is given by:

    with acceleration of gravity g and water depth H.

    With typical values for Lianzhou Bay, i.e., L=15 km and H=5 m a resonance period of T=2.38 h is found. Thisis much shorter than the main tidal periods, which are 12 h (S2), 12 h 25 min (M2), 24 h (S1), or 24 h 50 min (O2). To be in resonance with semi-diurnal tides (T=12 h 25 min), for instance, the length scale of the bay would need to be at least 78 km, which is by far larger than the actual size of Lianzhou Bay. Therefore, the tidal range will not be amplified within the bay. This is supported by the close agreement between predicted tides in Beihai city and the measured tides at the 24-hour samplings (Fig.7). Tidal waves are long gravity waves, which travel with a phase speed ofHence, the wave lengthλis given by

    Take againT=12 h 25 min as an example, thenλequals to 313 km. Referring to Sammariet al.’s work (2006), when the length scale is smaller than one quarter of the wave length (L<312/4=78.2 km), then the semi-diurnal tide will be distributed almost equally everywhere in the bay. This is a plausible parameter to support the prediction of tidal distribution dyamics in the next section.

    3.2.5 Tidal distribution dynamic

    Through comparing the predicted tides at Beihai City with the measured water depth at the 24-h samplings, it was confirmed that the tides arrive at the mouths of Nanliu River with a very short time lag compared to Beihai Port (Fig.7). The predicted tides closely match the observed ones; the only significant differences occurred during events of anomalous southerly wind stress.

    Fig.7 comparison between tidal sampling data and predicted water depth. (A) Station No.4 24-hour sampling. Note the higher water level around midnight when southerly wind was prevalent for several hours (B) Station No. 5 24-hour sampling.

    The nearly simultaneous arrival of tides within Lianzhou bay is further supported by tide models which show the phases of the dominant tides at the Guangxi coastal areas,i.e., O1, K1and M2(Sun and Huang, 2001). Generally, on the northern hemisphere, tidal waves rotate counterclockwise around amphidromic points. In Lianzhou Bay, all tides progress from SE to NW, and it needs only 8 min for O1and K1or 4 min for M2, respectively, to travel from Beihai city to the northwestern estuaries of Nanliu River. These are much shorter than the tidal periods, hence, high and low tides all occur almost simultaneously throughout Lianzhou bay.

    These findings confirm that due to the small scale of the bay, the current dynamics are not considerably affected by the earth rotation. Tide phase, amplitude and current are found to be rather uniform over the entire bay area.

    3.3 Coastal Currents

    The sea water input into Lianzhou Bay stems from the coastal currents which flow across the mouth of the bay. Those currents are part of the general wind driven water circulation in Beibu Gulf, of which Lianzhou Bay is a sub-embayment (Sunet al., 2001).

    Cyclonic current circulation prevails during winter. Though there is some debate about summer circulation (Yuan and Deng, 1999; Baoet al., 2005), newer studies indicate a (weaker) cyclonic circulation also during sum-mer due to wind and tidal forcing (Caiet al., 2003), and inflow through the Qiongzhou Strait between Hainan and China’s mainland (Wuet al., 2008). The circulation includes mainly two components,i.e., a tidal residual current and a wind driven current. The latter is evidently stronger (Sunet al., 2001), with estimated current speeds being 3.5 times larger than the residual current induced by tides. Due to the prevalent strong NE monsoon winds during winter, the wind-driven current in Beibu Gulf forms a cyclonic circulation. During summer, wind forcing and hence wind driven circulation is generally weaker and less distinguishable. Because of the shallow depth of the study site, and the largely barotropic character of the circulation in Beibu Gulf (Wuet al., 2008), it is feasible to assume that the surface current accounts for the water renewal in Lianzhou Bay. The velocity of the NW coastal currents at the mouth of Lianzhou Bay, summarized from Sunet al.(2001), is shown in Table 2.

    Observations and calculations reveal that the coastal current off Lianzhou bay is always directed to northwest, but is 4 times faster during winter, associated with the stronger wind stress during that season.

    Table 2, coastal current velocities of Lianzhou Bay?

    3.4 Residence Time

    Residence time is a useful value to estimate the contribution of fresh water input and coastal current to the water exchange processes in the bay. It is defined as the period of time that a water parcel needs to leave a particular area (Takeoka, 1984).

    If the water in Lianzhou bay were to be renewed only by the fresh water runoff from Nanliu river (average discharge 168 m3 s?1,i.e., 0.4 km3per month), it would take 4.5 months until the entire volume of the bay of 1.79 km3is replaced. However, water is also exchanged on the seaward side of the bay. Tidal currents in Lainzhou bay are mainly flowing back and forth, as tides arrive simultaneously everywhere across the bay, and the bay is too small for currents to be affected by the rotation of the earth. Hence, tidal currents themselves do not contribute to water renewal, carrying the same water into and out of the bay. The coastal current is not expected to turn right into the bay as the bay is smaller than the barotropic Rossby radius (as discussed above), and is shallower than the water in the offshore main pathway of the coastal current. Nevertheless, during ebb tide water from Lianzhou bay is displaced seaward, where it is carried northwestward with the coastal current (Fig.8). At the same time, ‘new’ water is advected by the coastal current to the southern part of the bay, where it enters the bay with the next flood. We predict that by this interaction between tidal currents and coastal current, water from Lianzhou bay is gradually displaced to the northwest, and replaced by ‘new’ water from the southeast.

    During winter, water of the coastal current flows with 0.2 m s?1, and thus needs about 1.4 d to pass across the 24 km width of the mouth of Lianzhou Bay. Hence, every 1.4 d, all the water flushing the bay due to the tides has been completely renewed by Gulf water from the coastal current. During summer, with the coastal current being weaker at 0.05 m s?1, renewal may take 6 d. Although some aspects are simplified in this concept, it suggests that water renewal by tides and coastal current occurs at the scale of days, and is faster in winter than in summer. In comparison, water exchange due to freshwater inflow is much smaller, with residence times in the order of months.

    The longer residence time within the bay during summer implies more accumulation of nutrients from aquaculture runoff.

    Fig.8 Conceptual sketch of coastal current along the mouth of Lianzhou Bay (thick arrows) and tidal current in Lianzhou Bay (thin arrows). Water from the southern part of the bay is gradually displaced to the northern part (open arrows), while ‘new’ water enters the southern part (black arrows).

    4 Conclusion

    This study assessed water movement processes in the shallow Lianzhou Bay, Guangxi Province, China. The findings are based onin situdata of temperature, salinity, current speed,etc., which are used to calculate several hydrodynamic parameters,i.e., Brunt-V?is?l? Frequency, Richardson Number, Rossby radius and Resonance Period. Owing to the small size of the bay, tides are uniform throughout the entire area. The hydrographic conditions in the bay are found to be barotropic and well mixed due to high tides, and with only moderate freshwater runoff. Water is primarily renewed by the combined action oftidal currents and coastal current. Residence time in summer is markedly longer, as the mainly wind-driven coastal current is weaker during that time. Particularly during summer, nutrient-rich waste water is affecting the ecology of the bay, as the coastal region in Guangxi is heavily utilized for aquaculture. Further research is necessary to find out strategies that minimize the harmful effects of aquaculture effluents.

    Acknowledgements

    We acknowledge the financial support from the German Federal Ministry for Education and Research (BMBF) to the Leibniz-Zentrum für Marine Tropen?kologie GmbH (ZMT), Germany, for the project ‘The Role of Mangroves for Biogeochemical Fluxes into the Coastal Ecosystems under the Influence of Anthropogenic Alterations’ as part of the Sino-German Research Project ‘BEIBU-Holocene environmental evolution and anthropogenic impact of Beibu Gulf, South China Sea’, WTZ China (Grant No. 03F0607B) is acknowledged. The enduring logistical support from our colleagues from the Guangxi Mangrove Research Center in Beihai, Guangxi, enabled us to carry out the field study and is greatly appreciated. The authors are sincerely grateful for the comments and suggestions of the anonymous reviewers.

    Bao, X. W., Hou, Y. J., Chen, C. S., Chen, F., and Shi, M. C., 2005. Analysis of characteristics and mechanism of current system on the west coast of Guangdong of China in summer. Acta Oceanologica Sinica, 24 (4): 1-9.

    Cai, S. Q., Huang, Q. Z., and Long, X. M., 2003. Three- dimensional numerical model study of the residual current in the South China Sea. Oceanologica Acta, 26 (5-6): 597-607.

    Chawla, A., Jay, D. A., Baptista, A. M., Wilkin, M., and Seaton, C., 2008. Seasonal variability and estuary-shelf interactions in circulation dynamics of a river-dominated estuary. Estuaries and Coasts, 31: 269-288.

    Chen, B., 1996. Preliminary analysis of the characteristics of the tidal and residual currents in Lianzhou Bay. Guangxi Sciences, 3 (1): 18-21, 37.

    Chen, B., and Shi, M. C., 1996. Numerical calculation about the tidal current and wind current in the Lianzhou Bay. Guangxi Sciences, 3 (3): 32-35.

    Chen, J. F., 2003. Analysis on temperature, salinity and pH value of seawater in Lianzhou Bay. Journal of Guangxi Academy of Sciences, 19 (1): 40-43 (in Chinese with English abstract).

    Chen, Q. Y., 2001. Assessment on situation of water quality in Lianzhou Bay of Guangxi. Marine Environmental Science, 20 (2): 56-58 (in Chinese with English abstract).

    deCastro, M., Gomez-Gesteira, M., Prego, R., and Alvarez, I., 2004. Ria-ocean exchange driven by tides in the Ria of Ferrol (NW Spain). Estuarine, Coastal and Shelf Science, 61: 15-24.

    Fofonoff, P., and Millard Jr, R. C., 1983. Algorithms for computation of fundamental properties of seawater. Unesco Technical Papers in Marine Sciences, 44: 53.

    Forman, M., 2005. Tidal resonance in Juan de Fuca Strait and the Strait of Georgia. Notes and Correspondence, 35: 1279-1286.

    Galperin, B., Snkoriansky, S., and Anderson, P. S., 2007. On the critical Richardson number in stably stratified turbulence. Atmospheric Science Letter, 8: 65-69.

    Geyer, W. R., 1993. The importance of suppression of turbulence by stratification on the estuarine turbidity maximum. Estuaries, 16 (1): 113-125.

    Grabemann, I., and Krause, G., 2001. On different time scales of suspended matter dynamics in the Weser Estuary. Estuaries, 24 (5): 688-698.

    Greenburg, D. A., 1979. A numerical model investigation of tidal phenomena in the Bay of Fundy and Gulf of Marine. Marine Geodesy, 2: 161-187.

    Gu, D. Q., and Wu, S. Y., 2001. Analysis of the sediment sources and the transport tendency in the Southern Lianzhou Bay. Journal of Oceanography of Huanghai and Bohai Seas, 19 (1): 25-31 (in Chinese with English abstract).

    Hach Company, 2012. P.O. Box 389, Loveland, Colorado 80539-0389.

    Houry, S., Dombrowsky, E., De, M. P., and Minster, J. F., 1987. Brunt-V?is?l? Frequency and Rossy Radii in the South Atlantic. American Meteorological Society, 17: 1619-1626.

    Jiang, L. M., Chen, B., and Qiu, S. F., 2008. The relationship of silt transportation and ocean dynamical condition in Delta of Lianhou Bay. Journal of Guangxi Academy of Sciences, 24 (1): 25-28 (in Chinese with English abstract).

    Lai, T. H., and Wei, M. X., 2003. The relation of change of five nutrients and environmental factors in Lianzhou Bay. Journal of Guangxi Academy of Sciences, 19 (1): 36-39 (in Chinese with English abstract).

    Leder, N., and Orli?, M., 2004. Fundamental Adriatic seiche recorded by current meters. Annales Geophysicae, 22: 1449-1464.

    Liang, W., and Li, G. Z., 2002. Sediment dynamic effects and sediment transport in Lianzhou Bay. Nanhai Yanjiu Yu Kaifa, 2: 12-16 (in Chinese with English abstract).

    Lv, R. H., Xia, B., and Mao, X. H., 1995. A study on the primary productivity of Lianzhou Bay waters of Guangxi Province. Journal of Oceanography of HuangHai and Bohai Seas, 13 (2): 53-60 (in Chinese with English abstract).

    Macrander, A., 2009. R/V Uth?rn ISATEC student’s cruise to outer Weser estuary on 07.04.2009. Unpublished cruise report, Universit?t Bremen.

    Matsuno, T., Shigeoka, M., Tamaki, A., Nagata, T., and Nishimura, K., 1999. Distributions of water masses and currents in Tachibana Bay, west of Arioake Sound, Kyushu, Japan. Journal of Oceanography, 55: 515-529.

    Osiński, R., Rak, D., Walczowski, W., and Piechura, J., 2010. Baroclinic Rossby radius of deformation in the southern Baltic Sea. Oceanologia, 52 (3): 417-429.

    Páez-Osuna, F., 2001. The environmental impact of shrimp aquaculture: A global perspective. Environmental Pollution, 112: 229-231.

    Qiu, S. F., and Lai, T. H., 2004. The composition characteristics of nutrient salts and their relations with the main environmental factors in the seawaters of Lianzhou Bay. Journal of Guangxi Academy of Sciences, 20 (3): 179-181 (in Chinese with English abstract).

    Qiu, S. F., and Lai, T. H., 2005. The changing characteristics of contents of physical and chemical factors in the sea water of Lianzhou Bay. Journal of Guangxi Academy Sciences, 21 (3): 151-154 (in Chinese with English abstract).

    Sammari, C., Koutitonsky, V. G., and Moussa, M., 2006. Sea level variability and tidal resonance in the Gulf of Gabes, Tu-nisia. Continental Shelf Research, 26: 338-350.

    Schlitzer, R., 2012. Ocean Data View. http://odv.awi.de.

    Sun, H. L., and Huang, W. M., 2001. Observation, analysis and numerical study of the tide and current in the Guangxi offshore area. Journal of Oceanography of Huanghai and Bohai Seas, 14 (4): 12-21 (in Chinese with English abstract).

    Sun, H. L., Huang, W. M., and Zhao, J. S., 2001. Three-dimensional numerical simulation of tide-induced, wind-driven and thermohaline residual currents in the Beibu Gulf. Oceanologia Et Limnologia Sinica, 21 (5): 561-569 (in Chinese with English abstract).

    Takeoka, H., 1984. Fundamental concepts of exchange and transport time scales I a coastal sea. Continental Shelf Research, 3 (3): 311-326.

    Ulses, C., Grenz, C., Marsaleix, P., Schaaff, E., Estournel, C., Meulé, S., and Pinazo, C., 2005. Circulation in a semi- enclosed bay under influence of strong freshwater input. Journal of Marine Systems, 56: 113-132.

    Unnikrishnan, A. S., and Antony, M. K., 1990. On vertical velocity fluctuations and internal tides in an upwelling region off the west coast of India. Estuarine, Coastal and Shelf Science, 31: 865-873.

    Yuan, S., and Deng, J., 1999. A numerical study on circulation in the Beibu Gulf. Nanhai Yanjiu Yu Kaifa, 12 (2): 41-46 (in Chinese with English abstract).

    Wei, M. X., and He, B. M., 1998. A study on eutrophication and red tide formation in Lianzhou Bay. Tropical Oceanology, 17 (4): 66-72.

    Wei, M. X., He, B. M., and Lai, T. H., 2004. The temporal and spatial distribution of pH value and DO and their relation with the environmental factors during the formation of the Algal Bloom in Lianzhou Bay. Guangxi Sciences, 11 (3): 221-224.

    Wu, D. X., Wang, Y., Lin, X. P., and Yang, J. Y., 2008. On the mechanism of the cyclonic circulation in the Gulf of Tonkin in the summer. Journal of Geophysical Research, 113: 1-10.

    (Edited by Xie Jun)

    * Corresponding author. E-mail: entergoing_9@hotmail.com

    (Received March 13, 2012; revised March 22, 2012; accepted July 26, 2013)

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

    免费看a级黄色片| 久久久久久久久久久丰满| 国产视频内射| 好男人在线观看高清免费视频| 国产日本99.免费观看| 国产黄色视频一区二区在线观看 | 免费电影在线观看免费观看| 长腿黑丝高跟| 亚洲av免费在线观看| 亚洲精品色激情综合| 有码 亚洲区| 色哟哟·www| 嫩草影院精品99| 中文字幕熟女人妻在线| 青青草视频在线视频观看| 毛片女人毛片| 国产亚洲av嫩草精品影院| 免费观看的影片在线观看| 午夜免费激情av| 亚洲精品国产av成人精品| 亚洲欧美中文字幕日韩二区| 欧美xxxx性猛交bbbb| 99热这里只有是精品在线观看| 如何舔出高潮| 精品99又大又爽又粗少妇毛片| 亚洲图色成人| 中文字幕久久专区| 中文资源天堂在线| 伊人久久精品亚洲午夜| 国产黄片美女视频| 日本黄大片高清| 久久99蜜桃精品久久| 成人国产麻豆网| 国产一区亚洲一区在线观看| 国产成年人精品一区二区| 亚洲欧洲日产国产| 爱豆传媒免费全集在线观看| 在线观看一区二区三区| 搡女人真爽免费视频火全软件| 99久久精品热视频| 亚洲电影在线观看av| 国产熟女欧美一区二区| 人体艺术视频欧美日本| 一级黄色大片毛片| 三级国产精品欧美在线观看| ponron亚洲| 国产一区二区三区av在线 | 99久久精品国产国产毛片| 亚洲国产精品久久男人天堂| 亚洲人成网站在线播放欧美日韩| 日韩三级伦理在线观看| 毛片女人毛片| 天堂av国产一区二区熟女人妻| 国产女主播在线喷水免费视频网站 | 九九热线精品视视频播放| 在现免费观看毛片| 欧美最黄视频在线播放免费| 欧美+亚洲+日韩+国产| 国产成人福利小说| 国产精品一二三区在线看| 欧美变态另类bdsm刘玥| 久久99精品国语久久久| 我要搜黄色片| 日本三级黄在线观看| 亚洲色图av天堂| 成人三级黄色视频| 91在线精品国自产拍蜜月| 久久精品国产99精品国产亚洲性色| 91精品一卡2卡3卡4卡| 91午夜精品亚洲一区二区三区| 亚洲精品粉嫩美女一区| 色视频www国产| 尾随美女入室| 亚洲自拍偷在线| 激情 狠狠 欧美| 亚洲国产精品成人久久小说 | 波多野结衣高清作品| 男人和女人高潮做爰伦理| 精品久久久久久久久久免费视频| 波多野结衣高清无吗| 色哟哟·www| 久久精品夜色国产| 亚洲av成人精品一区久久| 久久亚洲国产成人精品v| 国内少妇人妻偷人精品xxx网站| 日本五十路高清| 国产视频首页在线观看| 欧美极品一区二区三区四区| 波野结衣二区三区在线| 99久久精品热视频| 亚洲欧美中文字幕日韩二区| 国产免费男女视频| 日韩欧美一区二区三区在线观看| 亚洲欧美清纯卡通| 亚洲人成网站高清观看| 国产三级在线视频| 91在线精品国自产拍蜜月| 国产精品久久电影中文字幕| 成人亚洲欧美一区二区av| 99热网站在线观看| 禁无遮挡网站| 欧美变态另类bdsm刘玥| 一区福利在线观看| 在线观看午夜福利视频| 波多野结衣巨乳人妻| 日日啪夜夜撸| 午夜亚洲福利在线播放| 亚洲国产精品国产精品| ponron亚洲| 激情 狠狠 欧美| 免费看日本二区| 成年版毛片免费区| 亚洲国产欧美人成| 久久精品久久久久久久性| 欧美人与善性xxx| 91av网一区二区| 久久韩国三级中文字幕| 精品国内亚洲2022精品成人| 自拍偷自拍亚洲精品老妇| 国产黄色小视频在线观看| 精品日产1卡2卡| 国产黄片美女视频| 色综合站精品国产| 一卡2卡三卡四卡精品乱码亚洲| 网址你懂的国产日韩在线| 精品一区二区免费观看| 乱系列少妇在线播放| 日韩成人伦理影院| 国产蜜桃级精品一区二区三区| 春色校园在线视频观看| 亚洲欧美成人精品一区二区| 国产精品av视频在线免费观看| 日韩欧美精品免费久久| 18+在线观看网站| 91av网一区二区| 国模一区二区三区四区视频| 精品久久久久久久人妻蜜臀av| 国产av在哪里看| 亚洲av第一区精品v没综合| 国产91av在线免费观看| 久久精品人妻少妇| 日本av手机在线免费观看| 国产精品久久电影中文字幕| 日韩av不卡免费在线播放| 青青草视频在线视频观看| 亚洲精品乱码久久久v下载方式| 中文亚洲av片在线观看爽| 九色成人免费人妻av| 久久久久性生活片| h日本视频在线播放| 97超视频在线观看视频| 国产精品久久久久久久久免| 午夜福利在线观看免费完整高清在 | 九色成人免费人妻av| 国产精品女同一区二区软件| 一个人观看的视频www高清免费观看| 国产人妻一区二区三区在| 色视频www国产| 18+在线观看网站| 欧美又色又爽又黄视频| 亚洲欧洲日产国产| 免费观看的影片在线观看| 国产大屁股一区二区在线视频| 欧美高清成人免费视频www| 亚洲欧美日韩东京热| 我的老师免费观看完整版| 久久中文看片网| 嫩草影院入口| 精品日产1卡2卡| 久久久a久久爽久久v久久| 精品久久久久久成人av| 麻豆一二三区av精品| 国产高清不卡午夜福利| 男人舔奶头视频| 97超碰精品成人国产| 午夜福利在线观看免费完整高清在 | 久久精品国产亚洲av香蕉五月| 黄片wwwwww| 久久综合国产亚洲精品| 国产真实乱freesex| 尤物成人国产欧美一区二区三区| 日韩国内少妇激情av| 此物有八面人人有两片| 99热网站在线观看| 亚洲婷婷狠狠爱综合网| 18+在线观看网站| 国产片特级美女逼逼视频| 国产伦在线观看视频一区| 床上黄色一级片| or卡值多少钱| 亚洲国产欧美在线一区| 又粗又硬又长又爽又黄的视频 | 一个人看视频在线观看www免费| 日韩一区二区视频免费看| 麻豆成人av视频| av天堂中文字幕网| 22中文网久久字幕| 中文字幕av成人在线电影| 亚洲精品日韩av片在线观看| 男女做爰动态图高潮gif福利片| 亚洲内射少妇av| а√天堂www在线а√下载| 夜夜爽天天搞| 狂野欧美白嫩少妇大欣赏| 一卡2卡三卡四卡精品乱码亚洲| 免费一级毛片在线播放高清视频| 欧美+日韩+精品| 精品不卡国产一区二区三区| 久久精品91蜜桃| 国产极品天堂在线| 国产高清激情床上av| 精品久久久久久久久av| av卡一久久| 51国产日韩欧美| 色尼玛亚洲综合影院| 精品久久久噜噜| 看十八女毛片水多多多| 又粗又爽又猛毛片免费看| 久久精品夜夜夜夜夜久久蜜豆| 国产成人福利小说| 亚洲av熟女| 美女 人体艺术 gogo| 亚洲国产精品国产精品| 亚洲成av人片在线播放无| 国产午夜福利久久久久久| 亚洲精品日韩av片在线观看| 97在线视频观看| 日韩av在线大香蕉| 国产人妻一区二区三区在| 亚洲精品456在线播放app| 国产精品久久电影中文字幕| 国产伦精品一区二区三区四那| 国产精品.久久久| 在线观看av片永久免费下载| 久久亚洲精品不卡| 一级毛片电影观看 | 22中文网久久字幕| 国产一区二区三区av在线 | 老司机福利观看| 成人一区二区视频在线观看| 白带黄色成豆腐渣| 久久久午夜欧美精品| 美女国产视频在线观看| 国产精品女同一区二区软件| 一本久久精品| 又爽又黄无遮挡网站| 日韩大尺度精品在线看网址| 亚洲av二区三区四区| 看片在线看免费视频| 日日摸夜夜添夜夜添av毛片| 不卡一级毛片| www.色视频.com| 免费大片18禁| 简卡轻食公司| 丝袜喷水一区| 白带黄色成豆腐渣| 少妇人妻一区二区三区视频| 精品久久久久久久久av| 久久九九热精品免费| 亚洲欧美精品综合久久99| 国产伦精品一区二区三区四那| 一个人看的www免费观看视频| 亚洲乱码一区二区免费版| 波多野结衣高清作品| 好男人视频免费观看在线| 又爽又黄无遮挡网站| 亚洲人成网站在线观看播放| 成人特级黄色片久久久久久久| 婷婷色综合大香蕉| 国产一区二区在线av高清观看| 在线观看午夜福利视频| 精品久久久久久成人av| 国产伦一二天堂av在线观看| 人妻夜夜爽99麻豆av| av免费在线看不卡| 日本免费a在线| 午夜爱爱视频在线播放| 级片在线观看| av福利片在线观看| 精品久久久久久久久亚洲| 丰满乱子伦码专区| 亚洲欧美日韩高清专用| 亚洲图色成人| 国产精品乱码一区二三区的特点| 久久久午夜欧美精品| 一级毛片电影观看 | 国产精品av视频在线免费观看| 国产精品精品国产色婷婷| 日韩 亚洲 欧美在线| 亚洲人成网站高清观看| 国产白丝娇喘喷水9色精品| 亚洲乱码一区二区免费版| 中文字幕熟女人妻在线| 国产精品久久久久久精品电影小说 | 成年女人永久免费观看视频| 狠狠狠狠99中文字幕| 麻豆乱淫一区二区| 国产成人精品一,二区 | 少妇丰满av| 一个人免费在线观看电影| 热99re8久久精品国产| 亚洲人成网站高清观看| 欧美区成人在线视频| 深夜a级毛片| 欧美一区二区国产精品久久精品| 欧美在线一区亚洲| 婷婷色av中文字幕| 欧美高清成人免费视频www| 久久久久久九九精品二区国产| 中文精品一卡2卡3卡4更新| 18禁裸乳无遮挡免费网站照片| 免费黄网站久久成人精品| 日韩,欧美,国产一区二区三区 | 永久网站在线| 久久精品国产清高在天天线| 日韩av在线大香蕉| 女的被弄到高潮叫床怎么办| 国产成人一区二区在线| 亚洲精品久久久久久婷婷小说 | 在线观看免费视频日本深夜| 色噜噜av男人的天堂激情| 一级毛片久久久久久久久女| 中文字幕精品亚洲无线码一区| 日韩视频在线欧美| 国产熟女欧美一区二区| 久久久久久久久大av| 91久久精品国产一区二区成人| 国产日本99.免费观看| 亚洲无线观看免费| 波多野结衣高清无吗| 国产一区二区亚洲精品在线观看| 亚洲国产色片| 天堂av国产一区二区熟女人妻| 男的添女的下面高潮视频| 一边亲一边摸免费视频| 国产真实乱freesex| 啦啦啦观看免费观看视频高清| 国产成人精品久久久久久| 国产精品一及| 日韩亚洲欧美综合| 亚洲在线观看片| 国产精品综合久久久久久久免费| 国产av不卡久久| 在线免费十八禁| 日韩欧美在线乱码| 久久久午夜欧美精品| 中国国产av一级| 久久久国产成人精品二区| 国产成人freesex在线| www日本黄色视频网| 欧美xxxx黑人xx丫x性爽| 丝袜喷水一区| 欧美精品一区二区大全| 中文字幕熟女人妻在线| 桃色一区二区三区在线观看| 熟妇人妻久久中文字幕3abv| 26uuu在线亚洲综合色| 亚洲aⅴ乱码一区二区在线播放| 99热这里只有精品一区| 全区人妻精品视频| 国产国拍精品亚洲av在线观看| 可以在线观看的亚洲视频| 少妇的逼水好多| 97超碰精品成人国产| 天天躁夜夜躁狠狠久久av| 一级毛片电影观看 | 精品久久久久久久久av| 国产精品美女特级片免费视频播放器| 国产精品野战在线观看| 日产精品乱码卡一卡2卡三| 精品免费久久久久久久清纯| 中文字幕制服av| h日本视频在线播放| 美女黄网站色视频| 日韩亚洲欧美综合| 赤兔流量卡办理| 2021天堂中文幕一二区在线观| 在线观看美女被高潮喷水网站| 欧美成人精品欧美一级黄| 黄色日韩在线| 久久中文看片网| 好男人在线观看高清免费视频| av专区在线播放| 男的添女的下面高潮视频| 噜噜噜噜噜久久久久久91| 麻豆av噜噜一区二区三区| 国产成人freesex在线| 国产不卡一卡二| 免费人成视频x8x8入口观看| 观看免费一级毛片| 男的添女的下面高潮视频| 欧美性猛交╳xxx乱大交人| 美女高潮的动态| 日本熟妇午夜| 亚洲欧洲国产日韩| 久久午夜亚洲精品久久| 人人妻人人看人人澡| 黄片无遮挡物在线观看| 性插视频无遮挡在线免费观看| 黄色欧美视频在线观看| 日韩精品青青久久久久久| 国产精品女同一区二区软件| 欧美高清性xxxxhd video| 欧美最新免费一区二区三区| 丰满人妻一区二区三区视频av| 成人av在线播放网站| 成人高潮视频无遮挡免费网站| 免费电影在线观看免费观看| 亚洲欧美中文字幕日韩二区| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品av视频在线免费观看| 99在线视频只有这里精品首页| 亚洲精华国产精华液的使用体验 | 日韩高清综合在线| 日韩一区二区视频免费看| 国产三级在线视频| 久久午夜福利片| 国产激情偷乱视频一区二区| 日本五十路高清| 69av精品久久久久久| 青春草视频在线免费观看| 亚洲乱码一区二区免费版| 美女内射精品一级片tv| 国产探花在线观看一区二区| 亚洲乱码一区二区免费版| 久久午夜福利片| 国产激情偷乱视频一区二区| 久久国产乱子免费精品| 春色校园在线视频观看| 亚洲无线观看免费| 久久中文看片网| 韩国av在线不卡| 国产老妇女一区| 久久久精品欧美日韩精品| 成年av动漫网址| 91久久精品电影网| av视频在线观看入口| 老师上课跳d突然被开到最大视频| 国产免费一级a男人的天堂| 久久综合国产亚洲精品| 真实男女啪啪啪动态图| 日韩欧美精品v在线| 插阴视频在线观看视频| 18+在线观看网站| 久久精品久久久久久噜噜老黄 | 久久九九热精品免费| 国产亚洲精品av在线| 国产精品免费一区二区三区在线| 非洲黑人性xxxx精品又粗又长| 在线观看免费视频日本深夜| 99九九线精品视频在线观看视频| 国产在线精品亚洲第一网站| 国产精品福利在线免费观看| 毛片女人毛片| 国内精品一区二区在线观看| 久久精品综合一区二区三区| 欧美极品一区二区三区四区| 午夜福利高清视频| 国模一区二区三区四区视频| 午夜免费男女啪啪视频观看| 国产精品av视频在线免费观看| 非洲黑人性xxxx精品又粗又长| 国产老妇女一区| 美女xxoo啪啪120秒动态图| 99热这里只有是精品50| 99久久九九国产精品国产免费| 国产高清视频在线观看网站| 啦啦啦啦在线视频资源| 成年女人看的毛片在线观看| 成人漫画全彩无遮挡| 精品午夜福利在线看| 99riav亚洲国产免费| 美女被艹到高潮喷水动态| 少妇被粗大猛烈的视频| 国产精品一区www在线观看| h日本视频在线播放| 国产一区二区三区在线臀色熟女| 99久久成人亚洲精品观看| 久久99精品国语久久久| 亚洲婷婷狠狠爱综合网| 久久亚洲国产成人精品v| 免费不卡的大黄色大毛片视频在线观看 | 日本一二三区视频观看| 国产激情偷乱视频一区二区| 99久久人妻综合| 2022亚洲国产成人精品| 日韩高清综合在线| 亚洲无线观看免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩人妻高清精品专区| 性插视频无遮挡在线免费观看| 天堂av国产一区二区熟女人妻| 免费电影在线观看免费观看| 久久久国产成人精品二区| 天天躁夜夜躁狠狠久久av| 一区二区三区免费毛片| 久久久久国产网址| a级毛片免费高清观看在线播放| 色尼玛亚洲综合影院| 日本与韩国留学比较| 简卡轻食公司| 美女脱内裤让男人舔精品视频 | 青春草国产在线视频 | 久久婷婷人人爽人人干人人爱| 免费看av在线观看网站| 永久网站在线| 一本久久精品| 看免费成人av毛片| 丰满乱子伦码专区| 联通29元200g的流量卡| 99热网站在线观看| 午夜福利成人在线免费观看| 人人妻人人澡欧美一区二区| 男人和女人高潮做爰伦理| 变态另类成人亚洲欧美熟女| 亚洲av一区综合| 亚洲av中文字字幕乱码综合| www.色视频.com| 亚洲av电影不卡..在线观看| 欧美日本亚洲视频在线播放| 麻豆国产av国片精品| 免费观看在线日韩| 欧美成人a在线观看| 欧美日韩国产亚洲二区| a级毛片a级免费在线| 亚洲图色成人| 男女下面进入的视频免费午夜| 欧美丝袜亚洲另类| 国产白丝娇喘喷水9色精品| 青青草视频在线视频观看| 亚洲丝袜综合中文字幕| 亚洲欧美成人综合另类久久久 | 少妇裸体淫交视频免费看高清| 国产av麻豆久久久久久久| 99热只有精品国产| 亚洲美女搞黄在线观看| 日韩欧美精品v在线| 非洲黑人性xxxx精品又粗又长| 12—13女人毛片做爰片一| 久久国内精品自在自线图片| 色综合色国产| 99国产精品一区二区蜜桃av| 欧美激情国产日韩精品一区| 舔av片在线| 国产av麻豆久久久久久久| 97超视频在线观看视频| 欧美极品一区二区三区四区| 亚洲精品亚洲一区二区| 99在线视频只有这里精品首页| 国产精品久久久久久久电影| 日本av手机在线免费观看| 国产午夜精品论理片| 国产视频首页在线观看| 青春草视频在线免费观看| 国产成人精品婷婷| 亚洲在线自拍视频| 男人舔女人下体高潮全视频| 免费看av在线观看网站| 黄片无遮挡物在线观看| a级毛色黄片| 91午夜精品亚洲一区二区三区| av女优亚洲男人天堂| 亚洲中文字幕日韩| 久久久久九九精品影院| 成人亚洲欧美一区二区av| 毛片一级片免费看久久久久| 岛国毛片在线播放| 91在线精品国自产拍蜜月| 不卡视频在线观看欧美| 天美传媒精品一区二区| 成人永久免费在线观看视频| 色综合站精品国产| 精品久久久久久久人妻蜜臀av| 国产亚洲91精品色在线| 天堂网av新在线| 国产一区亚洲一区在线观看| 国产 一区 欧美 日韩| 在线免费观看的www视频| av在线老鸭窝| 国产精品免费一区二区三区在线| 国产在视频线在精品| 久久这里只有精品中国| 久久精品国产亚洲av天美| 少妇丰满av| 亚洲不卡免费看| 久久99热这里只有精品18| 亚洲国产精品国产精品| 久久99热这里只有精品18| 国产av一区在线观看免费| 久久精品夜夜夜夜夜久久蜜豆| 日本三级黄在线观看| 美女国产视频在线观看| 真实男女啪啪啪动态图| 国产成人精品婷婷| 国产伦精品一区二区三区四那| 亚洲最大成人av| 麻豆成人av视频| 日本av手机在线免费观看| 国产成人freesex在线| 日韩av在线大香蕉| 成人亚洲欧美一区二区av| 亚洲精品国产av成人精品| 国产精品一区二区性色av| 三级毛片av免费| 自拍偷自拍亚洲精品老妇| 又爽又黄无遮挡网站| 一区福利在线观看| 亚洲国产精品久久男人天堂| 亚洲av.av天堂| 欧美最黄视频在线播放免费| 成年女人永久免费观看视频| 午夜免费男女啪啪视频观看| 超碰av人人做人人爽久久| 看免费成人av毛片| 精品久久久久久久末码|