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

    Water characteristics of abyssal and hadal zones in the southern Yap Trench observed with the submersible Jiaolong*

    2020-06-08 05:33:56LIUXuehaiLIUYongzhiCAOWeiSUNChengjun
    Journal of Oceanology and Limnology 2020年3期

    LIU Xuehai , , , LIU Yongzhi , , , , CAO Wei , SUN Chengjun

    1 First Institute of Oceanography, Ministry of Natural Resources, Qingdao 266061, China

    2 Laboratory for Regional Oceanography and Numerical Modeling, Pilot National Laboratory for Marine Science and Technology (Qingdao), Qingdao 266237, China

    3 Key Laboratory of Marine Science and Numerical Modeling, Ministry of Natural Resources, Qingdao 266061, China

    Received Dec. 27, 2018; accepted in principle Feb. 20, 2019; accepted for publication Jun. 10, 2019 ? Chinese Society for Oceanology and Limnology, Science Press and Springer-Verlag GmbH Germany, part of Springer Nature 2020

    Abstract Using observations in an applied cruise of the submersible Jiaolong, water characteristics, geostrophic transport, and turbulent mixing in abyssal and hadal zones of the southern Yap Trench were studied. The spatial structures of deep water show that the abyssal water is cold, saline, and oxygen rich. The hadal water has very small changes in potential temperature and potential density, and a little decrease in salinity and obvious decrease in oxygen. The isotherm, isopycnal, and isohaline are depressed in abyss over the central trench. The turbulent mixing is enhanced in the near-bottom zone and the hadal water on the trench slope, especially at the steep slope, the dissipation rate and difference usivity is strong, which weakens the stratifi cation. The geostrophic fl ows move southward in the western region of the trench and northward in the eastern region, indicating cyclonic circulation. In the central region of the trench, the water transport is ~1.74 Sv southward. In the hadal zone, the northward and southward transports are balanced. Our analysis suggests that the abyssal water in the southern Yap Trench is from Lower Circumpolar Water (LCPW) and the hadal water seems to be of the isolated local water rather than LCPW.

    Keyword: Yap Trench; abyssal and hadal; the submersible Jiaolong; diapycnal mixing; geostrophic fl ows

    1 INTRODUCTION

    The Yap Trench is an important segment of the arcuate trench system in the Northwest Pacifi c, linking the waters in Mariana-Caroline basins (Kaneko et al., 1998). The length of the trench is ~700 km with a maximum depth ~8 527 m (Fujioka et al., 2000). Previous studies mainly focused on morphology and tectonics of the trench (Fujioka et al., 2000; Kobayashi, 2004), and so far, the hydrodynamic observation and studies are rare.

    Johnson and Toole (1993) and Kawabe and Taira (1998) proposed that there are three types of deep watermass in the Northwest Pacifi c. The bottom salty and oxygen-rich water is referred to as Lower Circumpolar Water (LCPW); the water lying above LCPW is referred to as North Pacifi c Deep Water (NPDW); and the water marked by a minimum salinity above the NPDW is Antarctic Intermediate Water. Usually the water exceeding 3 000 m is defi ned as abyssal water (Macdonald et al., 2009; Hautala, 2018). The abyssal water in Mariana-Caroline regions is from LCPW, with the infl ow of LCPW into the East Mariana Basin (EMB) and a partial outfl ow to the West Mariana Basin (Uehara and Taira, 1990; Johnson and Toole, 1993; Yanagimoto and Kawabe, 2007; Huang et al., 2018). Siedler et al. (2004) found that the channels at the Yap-Mariana Junction and Caroline Ridge were essential for the exchange of deep water such as water fl owing to East Caroline Basin (ECB) from EMB. Using observations by the submersible Jiaolong system, Liu et al. (2018) studied the watermass properties and deep currents in the northern Yap Trench, and suggested that the deep water originates from EMB and ECB as a part of the west propagating LCPW. However, the deep circulation and water properties in the southern Yap Trench remain unrevealed. Besides, the understanding in the hadal zone, which is below 6 000 m, is absent.

    Fig.1 Bathymetry and observation stations of the No. 38 Cruise

    The water linkage between the southern Yap Trench and the northern part is obstructed by a sill at a depth of ~5 000 m (as shown in Fig.1). Consequently, the properties of abyssal and hadal water in the two parts maybe have some difference erences. This means the local characteristics exists in the southern Yap Trench. Because the southern trench connects the Caroline Basin, the hydrodynamic investigation in the trench is helpful for identifi cation of the water transportation and deep circulation in Mariana-Caroline regions. Additionally, the investigation addressed here is also important for understanding the local biologicalchemical processes and material exchange with adjacent basins and trenches.

    Table 1 The observation time and instruments of No. 38 Cruise in the southern Yap Trench

    We conducted the investigation in the northern Yap Trench mainly by the submersible Jiaolong in No. 37 of Chinese Ocean Cruises (Liu et al., 2018). As the fi rst deep-sea manned submersible in China (Liu et al., 2010), Jiaolong has made hundreds of deep dives (maximum down to 7 062 m in Mariana Trench in 2012) in the South China Sea, Pacifi c, and Indian Ocean for multidisciplinary explorations (Li, 2017). Jiaolong and its mother ship equip tens of instruments including the conductivity-temperature-depth profi ler (CTD). The information of instruments used in this study is listed in Table 1.

    In this paper, using the observations in an applied cruise of the submersible Jiaolong (No. 38 of Chinese Ocean Cruises), we aim to study the water characteristics, geostrophic fl ow and turbulent mixing of abyssal and hadal zones in the southern Yap Trench. The paper is organized as follows. The data and methods used are described in Section 2. Section 3 displays the results and the analysis of water characteristics, geostrophic fl ow and turbulent mixing. A brief discussion and conclusion are presented in Section 4 and Section 5, respectively.

    2 DATA AND METHOD

    2.1 The study region

    The Yap Trench is separated into the northern and southern parts by a sill at 9°N at the depth of around 5 000 m (Fig.1). The southern part forms an important segment of an arc trench system in the Western Pacifi c. The maximum depth reaches 8 527 m in the central basin according to the topography from the General Bathymetric Chart of the Oceans (GEBCO_2014). This topography dataset has a resolution of 30 arc-second interval grid, access from the web (http://www.gebco.net/data_and_products/gridded_bathymetry_data/). Across the trench at around 8°N, the study area is located in the middle of the southern Yap Trench and in the southern margin of the Philippine Sea Plate.

    In this paper, only abyssal and hadal water is studied. The range for an abyss is 3 000-6 000 m, and below that, the hadal zone is defi ned.

    2.2 Observations and data processing

    The No. 38 Cruise was conducted in the southern Yap Trench from June 3 to June 14, 2017, which is an applied oceanic scientifi c expedition of the deep-sea manned submersible Jiaolong with the R/V Xiangyanghong 0 9 as its mother ship, lead by China Ocean Mineral Resources R&D Association and fi nancially supported by the National Key Basic Research Program of China.

    As shown in Fig.1b, nine observation sites were deployed across the abyssal trench along slope in a zonal direction overall. Among them, four casts (C01, C02, C03, and C04) were conducted using mainly ship-based CTD and dissolved oxygen sensors, and fi ve dives of Jiaolong (D148, D149, D150, D151, and D152) were observed using a set of CTD and dissolved oxygen sensor to obtain vertical profi les of temperature, salinity, and dissolved oxygen. Table 1 shows the information of the instruments, geographical positions, observation periods, and casting or diving depths. The maximum depth was 4 487-5 952 m for ship-based CTD observation and 4 187-6 685 m for Jiaolong dives, and all dives reached the sea bottom. C01, D148, D149, D150, and D151 were located in the west slope of the trench, C03 was in the central region, and C04, C04, and D152 were at the east slope. Situated in the fl ank of the trench, D148 and C02 had the shallowest bottom and observation depths.

    The observation data were reprocessed with SBE (Sea-Bird equipment) Data Processing Software V7.2 to obtain the values of temperature, salinity, density, and concentration of dissolved oxygen in a 1-m interval depth grid by the cubic spline interpolation in vertical after manually eliminating abnormal values. Only the data during lowering process are used in the study. The accuracy of SBE 911 Plus CTD is 0.001°C for temperature and 0.000 3 S/m for conductivity, and that of SBE 16 Plus CTD is 0.005°C for temperature and 0.000 5 S/m for conductivity. The equation of state is Thermodynamic Equation of Seawater 2010 (TEOS- 10), which is implemented with Gibbs Sea Water (GSW) Oceanographic Toolbox (GSW, http://www.teos-10.org/).

    2.3 Computation of geostrophic fl ows

    The meridional geostrophic fl ows at transects of Section-C and Section-D were computed from observed CTD data using the thermal wind relation:

    where v and v0are geostrophic velocities at the depth z and at the reference level z0, respectively; g is the gravitational acceleration; f is the Coriolis parameter; ρ and ρ0is the potential density and characteristic potential density, respectively.

    In this paper, the reference level is 3 000 m according to previous abyssal studies and adjacent trenches and basins (Johnson and Toole, 1993; Taira et al., 2005; Huang et al., 2018; Liu et al., 2018). At this weak circulation level, the meridional geostrophic velocity across transects Section-C and Section-D is set to be zero, hence we treated the calculated fl ows relative to the level as the geostrophic fl ows at abyssal and hadal depths. The geostrophic fl ow and volume transport were estimated from the observations by Jiaolong-carried CTD and ship-based CTD, respectively.

    2.4 Computation of diapycnal mixing

    Turbulent diapycnal mixing were computed using the Thorpe method (Thorpe, 1977) from CTD data. The turbulence dissipation rate ε and the vertical diapycnal difference usivity K are formulated by:

    where the ratio of the Ozmidov scale to the Thorpe overturn scale c=0.8 (Dillon, 1982), the mixing eき ciency m=0.2 (Osborn, 1980); N is the buoyancy frequency obtained from CTD data (Klymak et al., 2008); d is the Thorpe scale which is defi ned as the root-mean-square of all displacements between the real depth for a measured density and the depths in a reordered steady density profi le within a complete overturn. The overturn is detected by comparing the original and sorted density profi les. Here, the tests are implemented using the method of a previous study (Yang et al., 2016) to remove the spurious overturns.

    The average dissipation rate ε and diapycnal difference usivity K over a depth range H are as follows: where Li, εi, and Kiare the depth range, dissipation rate and diapycnal difference usivity of the i-th overturn within the range H, respectively. Here, H is set to be 300 m. The above turbulent mixing values are computed from 3 000 m to the maximum CTD casting or diving depths of nine stations.

    3 RESULT AND ANALYSIS

    3.1 Structure of abyssal water observed by shipbased CTD

    Figure 2 shows the vertical profi les of temperature and potential temperature θ3(referred to 3 000 dbar), salinity, potential density σ3(potential density referred to 3 000 dbar) and dissolved oxygen in the abyssal zone of the trench between 3 000 m and 6 000 m observed by ship-based instruments.

    Fig.2 Vertical profi les below 3 000 m at stations C01-C04 observed with ship-based instruments

    The temperature gradually decreases to the minimum values 1.453-1.463°C at around 4 500 m from 1.576-1.595°C, and then continually increases with depth due to adiabatic pressure. The maximum temperature is observed as 1.572°C at 5 716 m of C01, 1.607°C at 5 952 m of C03, and 1.558°C at 5 600 m of C04. Obviously, at the critical depth from the abyss to the hadal zone, the temperature surpasses that at 3 000 m. In a layer of 3 000-3 300 m, the order of the mean temperatures is C01>C04>C03>C02. From 3 500 to 4 300 m, the overall order turns to be C02>C04>C03>C01. This indicates that in the top 1 300 m abyssal water, the temperature on the west trench slope decreases more than that on the east slope. In the ~150 m layer below ~4 500 m where temperature steadily keeps the minimum temperature, C01, on a steeper slope, has the highest temperature, and C03, near the central of the trench, is the coolest. This suggests that the temperature over the central trench is depressed more than that at the fl ank. This is similar to the phenomena in the Challenger Deep of the Mariana Trench (Huang et al., 2018). At lower depths, temperature increases to the almost same in the three stations. On the other hand, the comparison for depths with the same temperature shows that C02 has the shallowest level for isotherm in the upper ~300 m layer but the deepest below the layer. From ~4 500 m to ~5 300 m, the isotherm level is shallower at C01 and deeper at C04.

    The potential temperature decreases continually with depth (Fig.2a). At the reference level, θ3is 1.572-1.599°C. The minimum is 1.243°C, appearing in the lowest layer of abyss at C03. The cooling is obvious in the top 1 500 m and very small in the next ~1 000 m. Below 5 500 m, temperature is steady. Similar to the real temperature, potential temperature varies difference erently at four stations. The profi les show that from 3 000 m to 4 500 m, potential temperature decreases less at a fl ank station C02 and more at a central station C03, which is same as the aforementioned for real temperature. In the top 300 m layer, C02 is colder than other stations. However, C01 is the coolest in the next 1 000 m where the slope is steep. Below 4 500 m, the potential temperature at C03 varies less than that at the fl ank stations. The pattern of the isotherm level for potential temperature is the same to that for the real temperature.

    The profi les of salinity present a similar pattern to those of potential temperature (Fig.2a & b). From the reference level to the maximum depth, salinity continually increases to 34.682 at more than ~5 500 m from 34.655-34.657. There is an obviously increase in the top ~1 600 m and a slight rise for deeper water. In addition, C02 is more saline than those of other stations at 3 000-3 400 m but evolves into less saline in the next ~1 000 m layer, in which C01 is more saline. This means that the isohaline at C02 is shallower in the fi rst 400 m and deeper in the layer from 400 m to 1 400 m than other stations, and in this layer C01 has shallower isohaline. Generally, the central station C03 is less saline at the isobar levels in abyss.

    Potential density has the similar patterns to those of salinity and potential temperature both in vertical and in horizontal. At 3 000 m, the density σ3is 41.441-41.447 kg/m3, and then increases to 41.510 kg/m3at ~4 700 m. Below the depth, the increase is very small (0.003 kg/m3totally). At the isobar levels in the upper 300 m, C02 is denser than other stations. Between 3 600 m and 4 300 m, the density is highest at C01 and least at C02. Below 4 700 m, C04 has higher densities and C03 has lower. From an isopycnal perspective, C02 is shallower for a density less than 41.459 kg/m3, but its isopycnal is deeper for the potential density exceeding 41.468 kg/m3. Overall, C3 has the relatively deep isopycnal. This means that the isopycnal of abyssal water is depressed in the central region of the trench.

    The content of dissolved oxygen increases with depth. In the upper 1 600 m, the continuous increase is apparent (from ~2.90 mL/L to ~3.41 mL/L). It is indicated from Fig.2a-d that the saline, cold, and dense water is rich in oxygen. Therefore, the pattern for profi les of oxygen is consistent with those of salinity, potential density, and potential temperature (oxygen-rich corresponding to cold). As shown in Fig.2d, C02 has more oxygen than other stations between 3 000-3 400 m but less content in the next ~1 000 m layer, in which C01 is relatively rich in oxygen. In addition, below ~4 800 m, the oxygen content decreases obviously with the depth increase at the stations.

    In summary, the conventional CTD observation displays the abyssal water characteristics for difference erent regions across the trench, and the spatial patterns are consistent for temperature, salinity, density, and oxygen content. Their spatial structures are directly related to geostrophic fl ow and turbulent mixing, which are investigated in following sections.

    3.2 Water characteristics in abyssal and hadal zones observed by Jiaolong

    Figure 3 presents the profi les of potential temperature θ3, salinity, potential density σ3and dissolved oxygen at abyssal and hadal depths in the trench observed by Jiaolong instruments.

    Above 5 500 m, four elements have the similar variation with depth. Potential temperatures decrease monotonously and the others increase continually. This means the deeper water is more saline, colder, denser, and oxygen-richer. The vertical changes are marked between 3 000-4 500 m, which is consistent with the observed by ship-based CTD. From ~1.60°C, θ3decreases to 1.25°C at ~5 300 m at D151 and ~5 000 m at D152. Below the depths, the value decreases only 0.007°C, indicating that potential temperature varies very little in hadal water and lower abyss. From 34.664 at 3 000 m, salinity increases to the maximum value 34.685 at ~5 300 m. Below the depth it keeps steady even though there is a very little drop at hadal depths. σ3increases to ~41.515 kg/m3at ~5 000 m from ~41.448 kg/m3, and keeps homogeneous in the hadal zone. From ~3.5 mL/L, dissolved oxygen increases to over 3.95 mL/L (with a maximum ~4.0 mL/L at the depth of ~4 800 m) below ~4 600 m, and then decrease obviously below ~5 400 m with a minimum of 3.5 mL/L.

    Fig.3 Vertical profi les below 3 000 m along Section-D observed by Jiaolong-carried instruments

    Additionally, although there is a difference erence in the oxygen content between the observed by casting and that by Jiaolong using difference erent sensors, their distribution patterns are identical. Two kinds of observations show that oxygen increases with depth, and the cold, saline, and dense water is relatively rich in oxygen. This has been observed in adjacent trenches and basins (Johnson and Toole, 1993; Kaneko et al., 1998; Kato and Kawabe, 2009; Huang et al., 2018; Liu et al., 2018). In particular, using Jiaolong-carried instruments, we observed that the oxygen is poorer in the hadal zone than that in the abyssal layer lying above the zone. This is probably related to the shortage of oxygen supplement from LCPW. Besides, the higher biomass in the near-bottom and hadal zones (Ichino et al., 2015) needs more consumption of oxygen. The observation demonstrates spatial variability. The contours of potential temperature, salinity, and potential density in Section-D are depressed over the trench between 3 000 m and ~4 200 m, in which the isotherm, isopycnal and isohaline at the western fl ank are shallower than those over the central trench. This was observed by shipbased CTD. Among the layers of ~4 500-5 500 m in the trench, the isolines are tilting and lift from west to east. Below 5 500 m, the difference erence is not obvious across the trench, so the potential temperature, salinity and potential density are uniformly distributed in the hadal zone.

    In summary, the observation of submersible Jiaolong presents similar spatial pattern with ship-based CTD, reveals that the isolines of potential temperature, salinity and potential density are depressed over the central trench. In hadal zone, these elements are homogeneous, but the oxygen content decreases.

    Fig.4 The θ- S scatter diagram with σ3 contours overlain

    3.3 Characteristics in θ- S diagram and abyssal watermass

    Figure 4a & b present the θ- S scatter diagrams below the reference level observed by ship-based and Jiaolong-carried CTDs. The diagrams show that the abyssal and hadal water becomes colder and more saline with the increasing of depth, with the exception of hadal and near-hadal depths. In Fig.4a for shipbased CTD, both of potential temperature and salinity have the very small rate of change when θ3drops to ~1.245°C below ~5 500 m (as shown in Fig.2a & b). This is also observed by Jiaolong as indicated in Figs.4b, 3a & 3b. Moreover, Fig.4b shows that in hadal zone, potential temperature keeps homogeneous and salinity decreases a little.

    The relation of potential temperature versus salinity with potential density contours overlain shows the spatial difference erence between stations zonally across the trench. For the fact that a lower potential temperature corresponds to a higher salinity for abyssal water, the difference erence of the relation at difference erent stations is not distinct in the diagrams. Even so, some small difference erence enlarges the spatial disparity in isopycnal. For the ship-based observation, at the same θ3between ~1.59°C and ~1.47°C, the salinity at C04 is highest, resulting in denser water and shallower isopycnal. In general, the central station C03 has a relatively low salinity at a same temperature, consequently lowering the density and deepening the isopycnal. For observations of Jiaolong (Fig.4b), at the same depth, the water at D148 is relatively cold and saline, and D150 has higher temperature and lower salinity. This indicates that above 4 200 m, the isopycnal is shallower at the fl ank than in the central trench.

    3.4 Geostrophic fl ows and water transport

    Figure 5a shows the meridional geostrophic fl ows at Section-C estimated from temperature and salinity observed by ship-based CTD. From C01 to C04, the fl ows basically travel southward except for the upper ~1 800 m layer (C01-C03) in which there is weak northward fl ow (<1 cm/s). Between ~4 000 m and ~5 000 m, the southward fl ows are relatively strong (in 2-3 cm/s). From C04 to C02, the east fl ank of the trench shows northward fl ows with a maximum velocity of 3.3 cm/s.

    Figure 5b presents the fl ows along Section-D estimated from observations of Jiaolong. At the west fl ank, the fl ows travel southward in the whole abyssal column. In the bottom layer at D149 on the steep slope, the velocity is the highest (with a maximum of 13.5 cm/s). From a west station D150 to D152 on the east slope, the southward fl ows gradually become to be in an opposite direction. Above ~3 700 m, the east of ~137.6°E shows a weak northward fl ow. Between ~3 700-4 700 m, the fl ows become northward at ~137.8°E. Below the layer, the fl ows are basically southward in the west of ~137.8°E and northward in the east with a maximum speed of 2.5 cm/s at depths close to the hadal zone.

    Fig.5 Meridional geostrophic fl ows relative to the reference level in abyss crossing transects

    Fig.6 Meridional geostrophic transports across transects estimated between 3 000 m and abyssal depths

    Both contours show that the meridional geostrophic fl ows are southward in the west and northward in the east. This opposing pattern suggests the cyclonic circulation in abyss. With consideration of the spatial difference erence between D148-D149 and C01, the strong fl ows occurred on the west slope according to Jiaolong’s observations, while this pattern was not presented in Section-C. Jiaolong can observe in deeper water than the conventional method to obtain hadal fl ows, but it has no observation in the central area, resulting in a certain discrepancy in estimated fl ows between two methods. Nevertheless, they show a similar pattern for the fl ows.

    With the results shown in Fig.5, we estimated the geostrophic transports from the reference level to abyssal depths at Section-C and Section-D (Fig.6a). Along section-C, the net volume in northward keeps increasing with depth and reaches the maximum of ~0.26 Sv at ~4 000 m. The volume reaches zero at ~4 500 m and increases to ~0.8 Sv (the positive means southward) at the maximum observation depth ~5 950 m. At Section-D, the net transport keeps close to zero in the fi rst 500 m and then continually increases to ~1.74 Sv at 6 000 m. In hadal zone, the transport reaches the maximum ~1.77 Sv at ~6 500 m and then decreases slightly to ~1.75 Sv at the maximum diving depth of 6 685 m.

    Fig.7 Profi les of average dissipation rate ε over a depth range H=300 m

    Fig.8 Profi les of average diapycnal difference usivity K over a depth range H=300 m

    For the difference erence in width and locations between two sections, their variations in the net transport are similar but somewhat difference erent. We computed the transport along two close transects: D149-D152 for Jiaolong and C01-D152 for ship-based CTD (Fig.6b). The panel shows two identical profi les of transport volumes. In the central part across the trench, the net transport is close to zero at the layer above 3 500 m, and below the level, it keeps increasing with depth. In the hadal zone, the southward transport is balanced by the northward.

    In summary, the meridional geostrophic fl ows in abyssal and hadal zones are southward in the western part of the southern Yap Trench and opposite in the eastern. In the main region of the trench, the net transport of abyssal water is southward. The net transport of the hadal water decreases to zero.

    3.5 Diapycnal mixing

    The estimated dissipation rate and eddy difference usivity of abyssal water at C01-C04 and D148-D152 are shown in Figs.7 & 8. The magnitude of dissipation rate and difference usivity is ε~O (10-8-10-6W/kg) and K~O (10-3-10-1m2/s), respectively.

    At the stations of ship-based observation, the average dissipation rate is almost constant in the upper ~1 400 m layer. Below 4 100 m, the average dissipation increases with depth. At the deepest water, the rate reaches ~7.8× 10-8W/kg. C03 has the highest dissipation rate above 4 700 m except for ~4 100-4 400 m. The diapycnal difference usivity presents the similar pattern to the dissipation rate. The average difference usivity increases a little in the fi rst ~700 m. After a slight decrease, the difference usivity strengthens continually from 3 700 m to 5 900 m and reaches the maximum of ~1.1× 10-2m2/s. Generally, the water on the steep slope and the deeper abyssal water in the central trench show relatively strong mixing.

    At stations of Jiaolong dives, the magnitude of dissipation rate and difference usivity is same as those of C01-C04. Generally, the mean values are higher in deeper water. Above 5 200 m, both the dissipation rate and difference usivity are relatively low at D151 where the topography is fl at, while those at locations on the steep slope are quite high and weaken the stratifi cation (shown in Fig.3). This corresponds to the depressions of isolines of potential temperature, salinity, and potential density over the central trench as shown in Fig.3. Compared with the results at stations C01-C04, we found that the mixing is active in the hadal water on the steep slope, especially at the near-bottom zone where the geostrophic fl ows are strong and the gradient of potential density is weak.

    In summary, the dissipation rate and difference usivity keep almost constant in the upper layer and then increase with depth. The turbulent mixing is strong on the steep slope and hadal zone where the strong geostrophic fl ows and small gradient of potential density exist.

    4 DISCUSSION

    4.1 Comparison of water properties in the southern Yap Trench with the northern

    The observations reveal that in the southern Yap Trench, the water obviously becomes colder, denser, more saline more oxygen with the depth increase between 3 000 m and ~4 500 m, and below this layer, the rate of change is very small but the oxygen content decreases obviously. The vertical patterns are consistent to the observed in the northern Yap Trench (Liu et al., 2018). This indicates that two parts of Yap Trench have the same or similar properties of watermass. In addition, the isolines of temperature, salinity, and density are depressed over the central trench, and in hadal zone. This phenomenon was observed in the Challenger Deep of the Mariana Trench (Huang et al., 2018). However, difference erent spatial patterns were shown in the north of the trench by Liu et al. (2018), where the isopycnal continually becomes deeper from the western side of the trench to the eastern side. These difference erences lead to the difference erent patterns for geostrophic transport. In the southern trench, the abyssal circulation is cyclonic with the southward fl ows in the western area and the northward in the eastern, and the net transport is southward. While in the northern trench, the route of LCPW is northward. Additionally, the circulation in the Challenger Deep is cyclonic (Taira et al., 2004; Huang et al., 2018). Johnson (1998) reported that the cyclonic sense of circulation is often evident in trenches of the Pacifi c.

    Comparison shows the similar pattern of diapycnal mixing in the south and north of the Yap Trench. The magnitude of dissipation rate and difference usivity in the southern trench is same to that in the north (Liu et al., 2018). In both areas, the dissipation changes little at the depth range of 3 000-4 500 m, in the deeper layer, the diapycnal mixing increases with depth, and the enhancement of mixing weakens water stratifi cation in the deeper abyssal and hadal zone. Besides, the mixing is strong along the slopes in the two areas of the trench. The strong mixing seems to be related to the geostrophic fl ows and internal tide dissipation analyzed in the north trench by Liu et al. (2018). Additionally, our fi ndings agree largely with the measurements by Huang et al. (2018) in the Challenger Deep.

    4.2 The origin of the abyssal and hadal water in the Southern Yap Trench

    Previous studies proved that the deep water in trenches and basins of Mariana-Caroline regions including the northern Yap Trench is from LCPW (Siedler et al., 2004; Huang et al., 2018; Liu et al., 2018). The same characteristics of the abyssal watermass is cold, saline, and rich in oxygen, indicating that the abyssal water in the southern Yap Trench is linked to adjacent water by deep fl ows, which is of LCPW. Although in the past, no observation was conducted in the abyss of the southern Yap Trench, Kaneko et al. (1998) found a broad isopycnal depression and suggested that a part of the westward fl ow at the Mariana Basin enters the Caroline Basin via the Yap Trench based on the observation of watermass at about 7°N and other sections in the west North Pacifi c. Studies by Kawabe et al. (2003) and Siedler et al. (2004) also proved the connection between the Yap Trench and two basins where the deep water below NPDW is of LCPW.

    On the other hand, the fl ow of LCPW is blocked by sills in trenches and basins, and the downwelling water is diき cult to reach the depths below 6 000 m. Furthermore, the salinity and oxygen content decrease in hadal water (seen in Figs.2d, 3c-d & 4b), so they are lower in the observed hadal zone than those of LCPW. This suggests that the hadal water in the Southern Yap Trench maybe is the isolated local water rather than LCPW. Nevertheless, based on this local observation and very few previous studies, we cannot make sure whether the hadal water in the trench is of LCPW or isolated local water.

    5 CONCLUSION

    Based on the observations of ship-based CTD and the manned submersible Jiaolong, we analyzed the water characteristics, geostrophic transport and turbulent mixing in abyssal and hadal zones of the southern Yap Trench. The conventional CTD observation presented the results for the depth less than 6 000 m at stations covering the trench slope and the central area, while the submersible conducted observation from the sea surface to the sea bed below 6 000 m along the slope of the hadal trench. Two kinds of observation obtained the consistent and mutually complementary results.

    The similar spatial patterns for temperature, salinity, density and oxygen content are shown. The abyssal water is cold, saline and oxygen rich. In hadal zone, the potential temperature and potential density keep almost homogeneous, with a little decrease in salinity and a distinct reduction in the content of oxygen. Results show that the isolines of potential temperature, salinity and potential density are depressed over the central trench in abyss. At the steep slope where the stratifi cation is weakened, the dissipation rate and difference usivity are relatively large. The diapycnal mixing is enhanced in the near-bottom zone and also active in the hadal water.

    In abyssal and hadal zones of the southern Yap Trench, the estimated meridional geostrophic fl ows are southward in the western part and northward in the eastern. In the main region of the trench, the water transport is ~1.74 Sv and southward in the range of 3 000-6 000 m. Below 6 000 m, the southward transport is balanced by the northward. This suggests existence of cyclonic circulation.

    The analysis of properties of watermass indicates that the abyssal water in the southern Yap Trench is of LCPW, and suggests that the hadal water seems to be of the isolated local water rather than LCPW. However, the existed studies are still not enough to verify the attribution of the watermass. Moreover, the connections and water exchanges of the southern Yap Trench to adjacent trenches and basins such as the northern Yap Trench, Caroline Basin, and Mariana Basin are to be studied with more observations, in order to reveal the origins and fl ow paths of the abyssal and hadal water in the trench.

    6 DATA AVAILABILITY STATEMENT

    The data that support the fi ndings of this study, described in Section 2.2, are available from the National Deep Sea Center, Ministry of Natural Resources of China but restrictions apply to the availability of the data, which were used under license for the current study, and so are not publicly available. However, data are available from the fi rst author of this paper upon reasonable request and with permission of the National Deep Sea Center, Ministry of Natural Resources of China.

    7 ACKNOWLEDGMENT

    We thank the crews and scientists aboard R/V Xiangyanghong 09 and the submersible Jiaolong for data acquisition and sample collections in the No. 38 of Chinese Ocean Cruises. We acknowledge the National Deep Sea Center, Ministry of Natural Resources of China for providing observation data of Jiaolong. We thank Dr. YANG Jichao, Dr. WANG Chunsheng, and Dr. LU Bo for their cooperation and helps in data sharing. We also thank the anonymous reviewers for providing valuable suggestions.

    References

    Dillon T M. 1982. Vertical overturns: a comparison of Thorpe and Ozmidov length scales. Journal of Geophysical Research, 87(C12): 9 601-9 613, https://doi.org/10.1029/JC087iC12p09601.

    Fujiwara T, Tamura C, Nishizawa A, Fujioka K, Kobayashi K, Iwabuchi Y. 2000. Morphology and tectonics of the Yap Trench. Marine Geophysical Research, 21(1-2): 69-86, https://doi.org/10.1023/A:1004781927661.

    Hautala S L. 2018. The abyssal and deep circulation of the Northeast Pacifi c Basin. Progress in Oceanography, 160: 68-82.

    Huang C J, Xie Q, Wang D X, Shu Y Q, Xu H Z, Xiao J G, Zu T T, Long T, Zhang T C. 2018. Seasonal variability of water characteristics in the Challenger Deep observed by four cruises. Scientifi c Reports, 8: 11 791, https://doi.org/10.1038/s41598-018-30176-4.

    Ichino M C, Clark M R, Drazen J C, Jamieson A, Jones D O B, Martin A P, Rowden A A, Shank T M, Yancey P H, Ruhl H A. 2015. The distribution of benthic biomass in hadal trenches: A modelling approach to investigate the efference ect of vertical and lateral organic matter transport to the seafl oor. Deep Sea Research Part I: Oceanographic Research Papers, 100: 21-33, https://doi.org/10.1016/j.dsr.2015.01.010.

    Johnson G C, Toole J M. 1993. Flow of deep and bottom waters in the Pacifi c at 10°N. D eep Sea Research Part I: Oceanographic Research Papers, 40(2): 371-394, https://doi.org/10.1016/0967-0637(93)90009-R.

    Kaneko I, Takatsuki Y, Kamiya H, Kawae S. 1998. Water property and current distributions along the WHP-P9 section (137°-142°E) in the western North Pacifi c. Journal of Geophysical Research: Oceans, 103(C6): 12 959-12 984.

    Kato F, Kawabe M. 2009. Volume transport and distribution of deep circulation at 165°W in the North Pacifi c. Deep Sea Research Part I: Oceanographic Research Papers, 56(12): 2 077-2 087, https://doi.org/10.1016/j.dsr. 2009. 08.004.

    Kawabe M, Fujio S, Yanagimoto D. 2003. Deep-water circulation at low latitudes in the western North Pacifi c. Deep Sea Research Part I: Oceanographic Research Papers, 50(5): 631-656, https://doi.org/10.1016/S0967-0637(03)00040-2.

    Kawabe M, Taira K. 1998. Water masses and properties of 165°E in the western Pacifi c. Journal of Geophysical Research: Oceans, 1031(C6): 12 941-12 958, https://doi.org/10.1029/97JC03197.

    Klymak J M, Pinkel R, Rainville L. 2008. Direct breaking of the internal tide near topography: Kaena Ridge, Hawaii. Journal of Physical Oceanography, 38(2): 380-399, https://doi.org/10.1175/2007JPO3728.1.

    Kobayashi K. 2004. Origin of the Palau and Yap trench-arc systems. Geophysical Journal International, 157(3): 1 303-1 315, https://doi.org/10.1111/j.1365-246X.2003. 02244.x.

    Li X Z. 2017. Taxonomic research on deep-sea macrofauna in the South China Sea using the Chinese deep-sea submersible Jiaolong. Integrative Zoology, 12(4): 270-282, https://doi.org/10.1111/1749-4877.12254.

    Liu F, Cui W C, Li X Y. 2010. China’s fi rst deep manned submersible, JIAOLONG. Science China Earth Science, 53(10): 1 407-1 410, https://doi.org/10.1007/s11430-010-4100-2.

    Liu Y Z, Liu X H, Lv X Q, Cao W, Sun C J, Lu J, Wang C S, Lu B, Yang J C. 2018. Watermass properties and deep currents in the northern Yap Trench observed by the Submersible Jiaolong system. Deep Sea Research Part I: Oceanographic Research Papers, 139: 27-42, https://doi.org/10.1016/j.dsr.2018.06.001.

    Macdonald A M, Mecking S, Robbins P E, Toole J M, Johnson G C, Talley L, Cook M, Wijfference els S E. 2009. The WOCEera 3-D Pacifi c Ocean circulation and heat budget. Progress in Oceanography, 82(4): 281-325.

    Osborn T R. 1980. Estimates of the local rate of vertical difference usion from dissipation measurements. Journal of Physical Oceanography, 10(1): 83-89, https://doi.org/10.1175/1520-0485(1980)010<0083:EOTLRO>2.0.CO;2.

    Siedler G, Holfort J, Zenk W, Müller T J, Csernok T. 2004. Deep-water fl ow in the Mariana and Caroline Basins. Journal of Physical Oceanography, 34(3): 566-581, https://doi.org/10.1175/2511.1.

    Taira K, Kitagawa S, Yamashiro T, Yanagimoto D. 2004. Deep and bottom currents in the challenger Deep, Mariana Trench, measured with super-deep current meters. Journal of Oceanography, 6 0(6): 919-926, https://doi.org/10.1007/s10872-005-0001-y.

    Taira K, Yanagimoto D, Kitagawa S. 2005. Deep CTD casts in the challenger deep, Mariana trench. Journal of Oceanography, 61(3): 447-454, https://doi.org/10.1007/s10872-005-0053-z.

    Thorpe S A. 1977. Turbulence and mixing in a Scottish Loch. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 286(1334): 125-181, https://doi.org/10.1098/rsta.1977.0112.

    Uehara K, Taira K. 1990. Deep hydrographic structure along 12°N and 13°N in the Philippine Sea. Journal of the Oceanographical Society of Japan, 46(4): 167-176.

    Yanagimoto D, Kawabe M. 2007. Deep-circulation fl ow at mid-latitude in the western North Pacifi c. Deep Sea Research Part I: Oceanographic Research Papers, 54(12): 2 067-2 081.

    Yang Q X, Zhao W, Liang X F, Tian J W. 2016. Threedimensional distribution of turbulent mixing in the South China Sea. Journal of Physical Oceanography, 46(3): 769-788, https://doi.org/10.1175/JPO-D-14-0220.1.

    免费人成视频x8x8入口观看| 51国产日韩欧美| 午夜视频国产福利| 欧美日韩综合久久久久久| 一级a爱片免费观看的视频| 国产高清视频在线观看网站| 国产高清有码在线观看视频| 麻豆精品久久久久久蜜桃| 最近的中文字幕免费完整| 俺也久久电影网| 日本黄大片高清| 国产精品亚洲美女久久久| 亚洲国产精品国产精品| 中文在线观看免费www的网站| 亚洲不卡免费看| 亚洲国产精品久久男人天堂| 精品日产1卡2卡| 老司机福利观看| 日本三级黄在线观看| 成年女人毛片免费观看观看9| 午夜福利成人在线免费观看| 内地一区二区视频在线| 国产aⅴ精品一区二区三区波| 看黄色毛片网站| 国产免费一级a男人的天堂| 久久久久久久久大av| 少妇猛男粗大的猛烈进出视频 | 国内精品宾馆在线| 美女被艹到高潮喷水动态| 成年女人看的毛片在线观看| 亚洲一级一片aⅴ在线观看| 99国产精品一区二区蜜桃av| 久久久久久久久久久丰满| 欧美不卡视频在线免费观看| 99久久精品热视频| 亚洲五月天丁香| 一卡2卡三卡四卡精品乱码亚洲| 欧美绝顶高潮抽搐喷水| 午夜福利18| 欧美一区二区亚洲| 欧美激情在线99| 久久久欧美国产精品| 欧美日韩乱码在线| 一卡2卡三卡四卡精品乱码亚洲| 日韩三级伦理在线观看| 一个人看视频在线观看www免费| 夜夜夜夜夜久久久久| 精品欧美国产一区二区三| 97超视频在线观看视频| 日韩一区二区视频免费看| 亚洲欧美日韩高清专用| 国产精品精品国产色婷婷| 日本与韩国留学比较| 97超碰精品成人国产| 精品人妻一区二区三区麻豆 | 99riav亚洲国产免费| 美女免费视频网站| 午夜精品国产一区二区电影 | or卡值多少钱| 18禁黄网站禁片免费观看直播| 一级a爱片免费观看的视频| 日日摸夜夜添夜夜添av毛片| 亚洲国产色片| a级毛片a级免费在线| 免费无遮挡裸体视频| 中文在线观看免费www的网站| 免费人成视频x8x8入口观看| 五月伊人婷婷丁香| 午夜精品在线福利| 久久久久国产网址| 久久99热6这里只有精品| av黄色大香蕉| 久久九九热精品免费| 高清日韩中文字幕在线| 高清午夜精品一区二区三区 | 高清毛片免费观看视频网站| av免费在线看不卡| av国产免费在线观看| 日本三级黄在线观看| 国产精品美女特级片免费视频播放器| 简卡轻食公司| 91在线精品国自产拍蜜月| 国产精品免费一区二区三区在线| www日本黄色视频网| 一区福利在线观看| av中文乱码字幕在线| 精品熟女少妇av免费看| 看免费成人av毛片| 中文字幕av在线有码专区| 真人做人爱边吃奶动态| 午夜日韩欧美国产| 精品国内亚洲2022精品成人| 久久精品国产鲁丝片午夜精品| 美女 人体艺术 gogo| 天天躁日日操中文字幕| 亚洲婷婷狠狠爱综合网| 亚洲精品粉嫩美女一区| 日韩中字成人| 99国产极品粉嫩在线观看| 大香蕉久久网| 大香蕉久久网| 在线观看免费视频日本深夜| 看免费成人av毛片| 国产在视频线在精品| 国产久久久一区二区三区| 国产视频内射| 别揉我奶头 嗯啊视频| 国产精品野战在线观看| 国产精品一区二区免费欧美| 男人狂女人下面高潮的视频| 免费搜索国产男女视频| 国产黄a三级三级三级人| 高清日韩中文字幕在线| 两个人的视频大全免费| 色综合色国产| 日日摸夜夜添夜夜添av毛片| 婷婷亚洲欧美| 麻豆久久精品国产亚洲av| 免费在线观看影片大全网站| 3wmmmm亚洲av在线观看| 大型黄色视频在线免费观看| 国产三级在线视频| 久久久成人免费电影| 黄色视频,在线免费观看| 亚洲熟妇熟女久久| 国产精品一区二区三区四区久久| 一区福利在线观看| 久久婷婷人人爽人人干人人爱| 久久热精品热| h日本视频在线播放| 97超碰精品成人国产| 亚洲欧美成人精品一区二区| 亚洲欧美精品综合久久99| 久99久视频精品免费| 中文字幕精品亚洲无线码一区| 国产一区二区三区在线臀色熟女| 成人三级黄色视频| 精品人妻偷拍中文字幕| 国产乱人视频| av福利片在线观看| 真人做人爱边吃奶动态| 亚洲人成网站高清观看| 精品人妻一区二区三区麻豆 | 久久99热6这里只有精品| 国产片特级美女逼逼视频| 色视频www国产| 熟女人妻精品中文字幕| 日日摸夜夜添夜夜爱| 午夜福利在线观看吧| 国产蜜桃级精品一区二区三区| 日韩一本色道免费dvd| 欧美三级亚洲精品| 久久精品人妻少妇| 成人永久免费在线观看视频| 国产在线精品亚洲第一网站| 国产伦精品一区二区三区四那| 国产精品嫩草影院av在线观看| 国产欧美日韩一区二区精品| 国产亚洲91精品色在线| 亚洲欧美日韩无卡精品| 亚洲精品日韩在线中文字幕 | 精品日产1卡2卡| 国产亚洲av嫩草精品影院| 亚洲va在线va天堂va国产| 午夜免费激情av| 久99久视频精品免费| 悠悠久久av| 日本免费一区二区三区高清不卡| 婷婷亚洲欧美| 亚洲专区国产一区二区| 国产女主播在线喷水免费视频网站 | 99国产极品粉嫩在线观看| 男女做爰动态图高潮gif福利片| 午夜福利在线在线| 少妇人妻一区二区三区视频| 97超视频在线观看视频| 校园人妻丝袜中文字幕| 亚洲精品国产av成人精品 | 一级a爱片免费观看的视频| 久久精品综合一区二区三区| 别揉我奶头~嗯~啊~动态视频| 国产在视频线在精品| 欧美不卡视频在线免费观看| 少妇人妻一区二区三区视频| 国产一区二区三区在线臀色熟女| 麻豆成人午夜福利视频| 丝袜喷水一区| 伦理电影大哥的女人| 欧美日本亚洲视频在线播放| 国产精品三级大全| 别揉我奶头~嗯~啊~动态视频| av在线老鸭窝| h日本视频在线播放| 国产精品福利在线免费观看| 欧美一区二区国产精品久久精品| 熟女人妻精品中文字幕| 人妻丰满熟妇av一区二区三区| 女同久久另类99精品国产91| 人人妻人人澡欧美一区二区| 国产一区二区亚洲精品在线观看| 久久午夜福利片| 夜夜夜夜夜久久久久| 午夜免费男女啪啪视频观看 | 99久久精品国产国产毛片| 亚洲不卡免费看| 熟妇人妻久久中文字幕3abv| 亚洲av中文av极速乱| 国产精品一及| av卡一久久| 亚洲成a人片在线一区二区| 国产av麻豆久久久久久久| 波野结衣二区三区在线| 亚洲自拍偷在线| 99在线人妻在线中文字幕| 国模一区二区三区四区视频| 欧美3d第一页| 亚洲成av人片在线播放无| 日本五十路高清| 亚洲熟妇熟女久久| 亚洲av免费高清在线观看| 欧美成人精品欧美一级黄| 欧美日本视频| 日本免费一区二区三区高清不卡| 一进一出好大好爽视频| 国产精品久久久久久av不卡| 蜜桃亚洲精品一区二区三区| 欧美日韩国产亚洲二区| 国产乱人视频| 噜噜噜噜噜久久久久久91| 久久人人爽人人片av| 免费一级毛片在线播放高清视频| 亚洲人成网站在线播放欧美日韩| 91麻豆精品激情在线观看国产| 精品欧美国产一区二区三| 亚洲欧美日韩高清在线视频| 老熟妇乱子伦视频在线观看| 国产精品福利在线免费观看| 女生性感内裤真人,穿戴方法视频| 俺也久久电影网| 老师上课跳d突然被开到最大视频| 白带黄色成豆腐渣| 丰满的人妻完整版| 精品乱码久久久久久99久播| 一边摸一边抽搐一进一小说| 欧美极品一区二区三区四区| 成人高潮视频无遮挡免费网站| 香蕉av资源在线| 欧美色视频一区免费| 日韩在线高清观看一区二区三区| 欧美日韩在线观看h| 嫩草影院入口| 亚洲国产精品久久男人天堂| 国产高清不卡午夜福利| 一进一出好大好爽视频| 99久久精品国产国产毛片| 亚洲精品国产av成人精品 | 在线天堂最新版资源| 日韩中字成人| 夜夜爽天天搞| 别揉我奶头 嗯啊视频| 高清午夜精品一区二区三区 | 99九九线精品视频在线观看视频| 女人被狂操c到高潮| 亚洲国产精品合色在线| 热99在线观看视频| 性插视频无遮挡在线免费观看| 久久久久久久久久久丰满| 九九在线视频观看精品| 精品国产三级普通话版| 欧美成人a在线观看| av在线老鸭窝| 97碰自拍视频| 观看美女的网站| 色尼玛亚洲综合影院| 国产熟女欧美一区二区| 亚洲国产欧美人成| 一级黄片播放器| 不卡视频在线观看欧美| 久久这里只有精品中国| 亚洲精品国产av成人精品 | 麻豆国产97在线/欧美| 99热全是精品| 国产成人freesex在线 | 久久久久精品国产欧美久久久| 国内精品美女久久久久久| 亚洲国产欧美人成| 国产v大片淫在线免费观看| 欧美日韩综合久久久久久| 国产成人福利小说| 久久久国产成人精品二区| 国内精品一区二区在线观看| 色吧在线观看| 99热全是精品| 亚洲色图av天堂| 1024手机看黄色片| 午夜久久久久精精品| 国产精品综合久久久久久久免费| 欧美性猛交黑人性爽| 蜜桃亚洲精品一区二区三区| 免费观看在线日韩| 91久久精品国产一区二区三区| 国产69精品久久久久777片| 久久久久国产网址| 国产精品一区二区三区四区久久| 少妇人妻精品综合一区二区 | 男人和女人高潮做爰伦理| 日日摸夜夜添夜夜添小说| 午夜福利18| 亚洲aⅴ乱码一区二区在线播放| 男女那种视频在线观看| 亚洲精品亚洲一区二区| 日本爱情动作片www.在线观看 | 成人国产麻豆网| 久久天躁狠狠躁夜夜2o2o| 最近最新中文字幕大全电影3| 国产乱人视频| av视频在线观看入口| 亚洲av成人精品一区久久| 亚洲av中文av极速乱| 成人午夜高清在线视频| 真实男女啪啪啪动态图| av中文乱码字幕在线| a级毛色黄片| 国产精品久久久久久亚洲av鲁大| 色吧在线观看| 男插女下体视频免费在线播放| 国产久久久一区二区三区| 国内精品一区二区在线观看| 一级a爱片免费观看的视频| 日本a在线网址| 亚洲成a人片在线一区二区| 看免费成人av毛片| 黄片wwwwww| 国产高清有码在线观看视频| 国产av一区在线观看免费| 熟女电影av网| 少妇丰满av| 亚洲精品影视一区二区三区av| 直男gayav资源| av福利片在线观看| 日韩制服骚丝袜av| 最后的刺客免费高清国语| 99久久成人亚洲精品观看| 免费在线观看影片大全网站| 久久久久久大精品| 桃色一区二区三区在线观看| 天堂av国产一区二区熟女人妻| 1000部很黄的大片| 亚洲精品色激情综合| 美女 人体艺术 gogo| 久久亚洲国产成人精品v| 亚洲真实伦在线观看| 亚洲国产欧美人成| 亚洲成人久久爱视频| 女人被狂操c到高潮| 国产成人a区在线观看| a级一级毛片免费在线观看| 亚洲电影在线观看av| 日韩中字成人| 久久天躁狠狠躁夜夜2o2o| 亚洲中文日韩欧美视频| av在线观看视频网站免费| 成人漫画全彩无遮挡| 精品久久久久久成人av| 亚洲专区国产一区二区| 久久久久久国产a免费观看| 成人av在线播放网站| 亚洲中文日韩欧美视频| 国产一区二区三区av在线 | 日韩精品有码人妻一区| 国产成人aa在线观看| 联通29元200g的流量卡| 人人妻人人澡人人爽人人夜夜 | 一个人看视频在线观看www免费| 人人妻人人澡欧美一区二区| 嫩草影院精品99| 亚洲av不卡在线观看| 成人特级黄色片久久久久久久| 精品久久国产蜜桃| 狂野欧美白嫩少妇大欣赏| 国产免费男女视频| 国产精品久久久久久精品电影| 国产精品嫩草影院av在线观看| 亚洲最大成人中文| 国产蜜桃级精品一区二区三区| 午夜爱爱视频在线播放| 国产精品亚洲美女久久久| 日本色播在线视频| 亚洲婷婷狠狠爱综合网| 国产av麻豆久久久久久久| 午夜影院日韩av| 狠狠狠狠99中文字幕| 久久精品国产亚洲av天美| 亚洲精品亚洲一区二区| 偷拍熟女少妇极品色| 成年av动漫网址| 日本黄大片高清| 久久久精品欧美日韩精品| 99热这里只有精品一区| 22中文网久久字幕| 亚洲av一区综合| 蜜桃久久精品国产亚洲av| 午夜福利高清视频| 中国国产av一级| 一卡2卡三卡四卡精品乱码亚洲| 日本-黄色视频高清免费观看| 国国产精品蜜臀av免费| 又爽又黄无遮挡网站| 国产老妇女一区| 中国国产av一级| 国产 一区 欧美 日韩| 国产男靠女视频免费网站| 久久6这里有精品| 午夜久久久久精精品| 日日啪夜夜撸| 精品免费久久久久久久清纯| 亚洲国产精品国产精品| 91久久精品电影网| 嫩草影院精品99| 中文字幕人妻熟人妻熟丝袜美| 色哟哟哟哟哟哟| 精品久久国产蜜桃| 免费看av在线观看网站| 老熟妇乱子伦视频在线观看| 久久久久久国产a免费观看| 中文字幕免费在线视频6| 一区二区三区高清视频在线| 国产精品女同一区二区软件| 成人美女网站在线观看视频| 久久欧美精品欧美久久欧美| 狠狠狠狠99中文字幕| 国产乱人偷精品视频| 国产精品综合久久久久久久免费| 卡戴珊不雅视频在线播放| 亚洲人成网站高清观看| 欧美又色又爽又黄视频| 亚洲五月天丁香| 亚洲欧美成人精品一区二区| 亚洲av成人av| 国产91av在线免费观看| 久久精品人妻少妇| 嫩草影院入口| 久久久久久久久久成人| 国产精品,欧美在线| 久久精品夜色国产| 你懂的网址亚洲精品在线观看 | 日韩人妻高清精品专区| 久久韩国三级中文字幕| 日韩大尺度精品在线看网址| 国产一区二区三区av在线 | 嫩草影院入口| 日韩,欧美,国产一区二区三区 | 国产在线精品亚洲第一网站| 国产男靠女视频免费网站| 久久久国产成人免费| 日韩国内少妇激情av| 女生性感内裤真人,穿戴方法视频| 久久久久久久久久久丰满| 亚洲成人久久性| 欧美日韩乱码在线| 午夜精品一区二区三区免费看| 午夜a级毛片| 99热6这里只有精品| 欧美一区二区精品小视频在线| 欧美xxxx黑人xx丫x性爽| 国产精品人妻久久久久久| 国产成人91sexporn| 成人鲁丝片一二三区免费| 伊人久久精品亚洲午夜| a级毛片免费高清观看在线播放| 国产探花在线观看一区二区| 日本精品一区二区三区蜜桃| 99热这里只有精品一区| eeuss影院久久| 欧美另类亚洲清纯唯美| 2021天堂中文幕一二区在线观| 亚洲国产色片| 免费观看在线日韩| 日日干狠狠操夜夜爽| 国产午夜精品论理片| 偷拍熟女少妇极品色| 中文字幕av在线有码专区| 久久精品影院6| 青春草视频在线免费观看| 久久这里只有精品中国| 成年免费大片在线观看| 亚洲一区二区三区色噜噜| 国产精华一区二区三区| 插逼视频在线观看| 国模一区二区三区四区视频| av天堂中文字幕网| 国产黄色小视频在线观看| 欧美潮喷喷水| 久久精品国产清高在天天线| 可以在线观看的亚洲视频| 又黄又爽又免费观看的视频| 久久久久久国产a免费观看| 亚洲无线观看免费| 18禁裸乳无遮挡免费网站照片| 青春草视频在线免费观看| 在线免费观看的www视频| 久久精品国产鲁丝片午夜精品| 欧美日本亚洲视频在线播放| 能在线免费观看的黄片| 日本成人三级电影网站| 免费大片18禁| 国产av麻豆久久久久久久| 久久精品影院6| 天堂网av新在线| 日本色播在线视频| 亚洲成人久久爱视频| 亚洲av五月六月丁香网| 噜噜噜噜噜久久久久久91| 日本一本二区三区精品| 国产蜜桃级精品一区二区三区| 高清毛片免费看| 亚洲欧美成人综合另类久久久 | 91在线观看av| 亚洲国产欧洲综合997久久,| 男女那种视频在线观看| 99热只有精品国产| 亚洲精品456在线播放app| 成人毛片a级毛片在线播放| 免费观看的影片在线观看| 国产伦在线观看视频一区| 十八禁网站免费在线| 黄色欧美视频在线观看| 亚洲内射少妇av| 精品久久久久久久人妻蜜臀av| АⅤ资源中文在线天堂| 欧美极品一区二区三区四区| 少妇人妻一区二区三区视频| 免费在线观看成人毛片| 久久久久久九九精品二区国产| 三级毛片av免费| 色在线成人网| 成人二区视频| 成人三级黄色视频| 99热全是精品| 久久久午夜欧美精品| 99在线视频只有这里精品首页| 色综合色国产| 91在线精品国自产拍蜜月| 97超碰精品成人国产| 亚洲最大成人av| 国产精品人妻久久久久久| 国产高清视频在线观看网站| 丝袜喷水一区| 日韩制服骚丝袜av| 校园春色视频在线观看| 久久人人精品亚洲av| 乱系列少妇在线播放| av.在线天堂| 亚洲国产色片| 国产精品久久久久久亚洲av鲁大| 精品无人区乱码1区二区| 91久久精品电影网| 1024手机看黄色片| 日本黄大片高清| 日韩亚洲欧美综合| 小说图片视频综合网站| 午夜福利在线观看吧| av天堂中文字幕网| 男人和女人高潮做爰伦理| 欧美+亚洲+日韩+国产| 亚洲成a人片在线一区二区| 少妇人妻一区二区三区视频| 人妻久久中文字幕网| 日本色播在线视频| 成人鲁丝片一二三区免费| 欧美+亚洲+日韩+国产| 国产淫片久久久久久久久| 成人毛片a级毛片在线播放| 在现免费观看毛片| 国产极品精品免费视频能看的| 麻豆成人午夜福利视频| av天堂在线播放| 日本黄色片子视频| 91午夜精品亚洲一区二区三区| 能在线免费观看的黄片| 亚洲中文字幕一区二区三区有码在线看| 一区二区三区免费毛片| 又黄又爽又免费观看的视频| 一区二区三区四区激情视频 | 色综合色国产| 欧美xxxx黑人xx丫x性爽| 国产精品乱码一区二三区的特点| 日本免费a在线| 婷婷精品国产亚洲av| 男人的好看免费观看在线视频| 亚洲中文字幕一区二区三区有码在线看| 日本一二三区视频观看| av在线蜜桃| 亚洲国产色片| 一进一出抽搐gif免费好疼| av女优亚洲男人天堂| 成人特级黄色片久久久久久久| 男人和女人高潮做爰伦理| 神马国产精品三级电影在线观看| 久久人人爽人人片av| 久久久久久九九精品二区国产| 熟女人妻精品中文字幕| 久久久久久九九精品二区国产| av女优亚洲男人天堂| 永久网站在线| 国产真实乱freesex| 此物有八面人人有两片| 国产视频内射| 亚洲国产精品sss在线观看| 在线播放国产精品三级| 一进一出好大好爽视频| 午夜久久久久精精品| 黑人高潮一二区| 亚洲av中文字字幕乱码综合| 露出奶头的视频|