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

    The IOD-ENSO precursory teleconnection over the tropical Indo-Pacific Ocean: dynamics and long-term trends under global warming*

    2018-04-02 03:03:44YUANDongliang袁東亮HUXiaoyue胡曉悅XUPeng徐鵬ZHAOXia趙霞YukioMASUMOTOHANWeiqing韓衛(wèi)清
    Journal of Oceanology and Limnology 2018年1期
    關鍵詞:徐鵬

    YUAN Dongliang (袁東亮) , HU Xiaoyue (胡曉悅) XU Peng (徐鵬) ZHAO Xia (趙霞) Yukio MASUMOTO , HAN Weiqing (韓衛(wèi)清)

    1 Key Laboratory of Ocean Circulation and Waves (KLOCW), and Institute of Oceanology, Chinese Academy of Sciences, and Qingdao National Laboratory for Marine Science and Technology, Qingdao 266071, China

    2 University of Chinese Academy of Sciences, Beijing 100049, China

    3 Qingdao Collaborative Innovation Center of Marine Science and Technology, Qingdao 266071, China

    4 University of Tokyo, Tokyo, Japan

    5 Department of Atmospheric and Oceanic Sciences, University of Colorado, Boulder, Colorado, USA

    1 INTRODUCTION

    Due to its large longitudinal span, the tropical Pacific and Indian Oceans play an important role in the variability and predictability of the global climate. The El Ni?o-Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD) are two dominant modes of interannual variations over the tropical Pacific and Indian Oceans, respectively, the interactions of which have received intensive attention. Early studies suggest that the interactions are through the Walker Circulation(Wu and Meng, 1998; Klein et al., 1999; Alexander et al., 2002; Behera and Yamagata, 2003; Clarke and Gorder, 2003; Lau and Nath, 2003; Wu and Kirtman,2004; Annamalai et al., 2005; Lau et al., 2005; Behera et al., 2006; Kug et al., 2006; Izumo et al., 2010; Luo et al., 2010), which is called the atmospheric bridge.Recent studies suggest that the precursory relation between IOD and ENSO at the one-year lead time is dominated by the Indonesian Throughflow (ITF)variability, which is called the oceanic channel (Yuan et al., 2011, 2013; Xu et al., 2013; Zhou et al., 2015).

    The IOD-ENSO precursory relationship provides an opportunity for enhanced predictability of ENSO at the lead time of one year across the spring persistence barrier (Clarke and Gorder, 2003; Izumo et al., 2010; Luo et al., 2010; Yuan et al., 2011, 2013;Zhou et al., 2015). Statistically, about a half of the standard deviation of the cold tongue sea surface temperature anomalies (SSTA) can be predicted at the one-year lead if IOD is used as a precursor. Lag correlation analyses have suggested that sea level anomalies propagate from the eastern equatorial Indian Ocean to the western equatorial Pacific Ocean through the Indonesian seas, suggesting the oceanic channel dynamics leading to the enhanced predictability at the one-year lead time (Yuan et al.,2013, 2016; Xu et al., 2013; Xu and Yuan, 2015). In contrast, lag correlation analyses on the atmospheric bridge suggest no significant teleconnection at time lags longer than a season. So far, the oceanic channel dynamics have only been investigated in observations with scarce data in the Indonesian seas and in coarseresolution climate model simulations. No estimate of the penetration rate of the equatorial Kelvin waves from the equatorial Indian Ocean to the equatorial western Pacific Ocean has been provided. In this study, we use forced simulations of eddy-permitting global general ocean circulation models (OGCM) to diagnose the oceanic channel dynamics between IOD and ENSO. In addition, we use a continuously stratified ocean model to investigate the dynamics of the Kelvin wave propagation through the Indonesian seas and to estimate the penetration rate for the Kelvin wave propagation.

    Previous studies have suggested that the intraseasonal, semiannual, and annual Indian Ocean baroclinic Kelvin waves can propagate into the Indonesian Seas along the Sumatra-Java island chain(Molcard et al., 2001; Pujiana et al., 2009; Susanto et al., 2012). Drushka et al. (2010) suggest that the intraseasonal waves propagate downward along the ray path and will be blocked by the Lombok Strait sill depth if they travel across a long distance. It is not clear whether the interannual waves can reach and modulate the Pacific Ocean circulation. In this paper,we use the high-resolution OGCM simulations to investigate the IOD-ENSO teleconnection in the period of 1990 to 2008 to complement the observational analysis of Yuan et al. (2013).

    The IOD-ENSO teleconnection is found to have decadal and multi-decadal variations (Xu et al., 2013;Izumo et al., 2014; Xu and Yuan, 2015; Xu and Yuan,2016; Yuan et al., 2017). No investigation of the centennial trend of the teleconnection has been conducted so far. Analyses by Yuan et al. (2013, 2017)and Xu and Yuan (2015, 2016) have shown that the period since the 1990’s is in the positive phase (strong precursory relation) of the decadal variability of the teleconnection. We will use the observed SSTA and the Coupled Model Intercomparison Project phase 5(CMIP5) simulations to investigate the centennial trend of the IOD-ENSO teleconnection under global warming.

    The paper is organized as the following. Section 2 introduces the observational data and the model simulations used in the study. Section 3 presents the results of the lag correlation analyses of the Ocean General Circulation Model for the Earth Simulator(OFES) and Hybrid-Coordinate Ocean Model(HYCOM) simulations during 1990–2008, which are compared with the analyses based on observational data. Results from a simple linear continuously stratified ocean model experiments are presented in Section 4. The centennial trends of the IOD-ENSO teleconnections are investigated in Section 5.Conclusions are summarized in Section 6 of this paper.

    2 DATA AND MODEL

    2.1 Observation

    The Hadley Centre Global Sea Ice and Sea Surface Temperature (HadISST) data used in this study cover the global ocean with a 1° longitude by 1° latitude resolution from 1870 to near present (Rayner et al.,2003). The extended reconstruction of global SST(ERSST V3b) data covers the global ocean on a 2°×2°horizontal grid from 1854 to present (Smith et al.,2008). The Kaplan SST (Version 2.0) data cover the global ocean on a 5°×5° grid from 1856 to present(Kaplan et al., 1998). The merged altimeter sea level data from the Aviso project cover the global oceans on a grid of 1/3° longitude by 1/3° latitude from 1993 to present. The subsurface temperature data from the Joint Environmental Data Analysis Center of the Scripps Institution of Oceanography (White, 1995)are on a grid of 5° longitude by 2° latitude at 11 levels in the equatorial Pacific Ocean from 1990 to 2003.The National Center for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) reanalysis surface zonal wind data are on a 2.5°×2.5° global grid from 1948 to present (Kalnay et al., 1996).

    Fig.1 Amplitude penetration rates of the Kelvin wave in the idealized LCSM experiments (with the barrier extending to 8°S) depending on the baroclinic modes

    2.2 The OFES and HYCOM model simulations

    The monthly output from the OFES model simulations during 1990–2008 is used to calculate the lag correlations. The OFES model grid covers the global ocean from 75°S to 75°N with a horizontal resolution of 0.1° longitude by 0.1° latitude and stretched vertical coordinate of 54 levels from the sea surface to the depth of 6 065 m. The model was spun up with the climatological NCEP/NCAR reanalysis data for 50 years and was integrated with realistic wind and surface flux data since 1990 (Masumoto et al., 2004).

    The Hybrid Coordinate Ocean Model (HYCOM,Halliwell et al., 1998; Bleck, 2002) with a horizontal resolution of 0.33° longitude by 0.33° latitude over the tropical Indo-Pacific oceans and 33 vertical layers from sea surface to 5 500 m is configured to the Indo-Pacific ocean basins from 30°E to 70°W, 55°S to 54.89°N. The model was spun up for 30 years using monthly climatology data from the Comprehensive Ocean-Atmosphere Dataset (Slutz et al., 1985). The European Center for Medium Range Weather Forecasts (ECMWF) 40 years reanalysis (ERA-40)wind data were used as the forcing for 1958–2001.From 2002 through 2008, the wind stress based on the QuikSCAT satellite scatterometer measurements is used as the forcing. More detailed descriptions of the model configuration can be found in Trenary and Han(2012, 2013). The model output for the period of 1990–2008 is used for this analysis, to be consistent with the OFES time series.

    2.3 The linear continuously stratified model

    Besides the OGCMs, which contain complete dynamics of the ocean circulation, a simple linear continuously stratified model (LCSM) is used to study the Kelvin wave propagation through the Indonesian Seas. This model is based on an averaged density profile of the tropical Indo-Pacific Ocean between 5°S and 5°N and is similar to the one described by Yuan and Han (2006). The model basin covers the domain between 20°S and 20°N and from 60°E to 160°E, with a barrier wall of the grid-size thickness extending from the northern boundary to the equatorial area at 110°E that separates the model domain into western and eastern basins. The two basins are connected only through the ocean south of the barrier wall. Sponge layers of 3° wide with a relaxation coeffi cients of 10-4/s are specified near the western boundary of the western basin and the eastern boundary of the eastern basin so that there are no reflections at these boundaries. The model horizontal resolution is 0.25° latitude by 0.25° longitude. The coeffi cients of the horizontal Laplacian and vertical mixing are 400 m2/s, and 1.3×10-8m2/s, respectively,which are very small and can be neglected for wave propagation study. The bottom of the LCSM is assumed to be flat at a depth of 4 000 m. The nonormal-flow and nonslip conditions are used at all boundaries. The Arakawa-C grid is adopted. Total 10 vertical modes are calculated for the LCSM, the solution of which is compared with the first baroclinic mode solution in the following text. Experiments suggest that the Kelvin wave penetration rate is essentially unchanged when the number of baroclinic modes is larger than 10 (Fig.1).

    The model is also deployed over the tropical eastern Indian and western Pacific Oceans (60°E to 160°E, 20°S to 20°N) with a high horizontal resolution of 0.1° longitude by 0.1° latitude suffi cient to resolve the main passages within the Indonesian Seas. The sea-land mask of the Indonesian seas follows the 200 m isobath derived from the ETOPO1 bathymetry data. Sponge layers are specified at the western and eastern boundaries to suppress the reflections.

    2.4 The CMIP5 simulations

    The centennial model outputs used in this study are the historical simulations of CMIP5 (Meehl and Bony,2011), with both anthropogenic and natural forcing in the experiments. A set of HistoricalNat experiments have considered only the natural forcing. Both the Historical and HistoricalNat experiments have covered the period from the mid-19th century to the present. Depending on the available data, we chose to analyze 20 Historical and 7 HistoricalNat simulations in this study. Most of the models have a few ensemble experiments in the Historical and HistoricalNat simulations. For simplicity, only one of the ensemble experiments is analyzed for each model.

    2.5 Statistical analysis

    Indices in this study are as follows. The Ni?o3.4 index is the area-averaged SSTA in the equatorial central-eastern Pacific Ocean (5°S–5°N,170°W–120°W). The Dipole Mode Index (DMI) is the diff erence of area-averaged SSTA between the western (10°S–10°N, 50°–70°E) and the southeastern(10°S–Equator, 90°–110°E) tropical Indian Oceans(Saji et al., 1999).

    Throughout the paper, the lag correlation is calculated between the anomalies over the Indo-Pacific Oceans in the following winter-fall seasons and the area-averaged oceanic anomalies in the southeastern tropical Indian Ocean (STIO in abbreviation, 10°S–Equator, 90°–110°E) in fall or the surface zonal wind anomalies (SZWA) over the far western equatorial Pacific (5°S–5°N, 130°–150°E) in fall. Boreal seasons are defined as winter (December to the next February), spring (March to May), summer(June to August), and fall (September to November)throughout this paper. We will call the IOD year Year 0 and the year following it Year+1. The confidence levels of the correlations are computed based on the Student’s- t test.

    3 RESULT OF THE OGCM SIMULATIONS

    In this section, the lag correlations of the IODENSO are calculated based on the OFES and HYCOM model simulations, which are compared with those based on observational analyses. The lag correlations of the atmospheric bridge in the wind forcing are also calculated and compared with the oceanic channel dynamics of the lag correlations.

    3.1 Lag correlations of sea level anomalies

    The lag correlations between the STIO sea level anomalies (SLA) in fall and the Indo-Pacific SLA in the following winter through fall seasons in the OFES and HYCOM simulations are compared with those of the satellite altimeter SLA during 1993–2008 (Fig.2).The lag correlations of both model simulations show east-west dipole patterns in the equatorial Indian and Pacific Oceans in winter, in good agreement with those based on satellite altimeter data. The dipole patterns are believed to be associated with the ENSO and IOD in the Pacific and Indian Oceans, which peaked in winter and fall, respectively. The lag correlations in the following spring through summer suggest propagation of SLA from STIO through the Indonesian seas into the western equatorial Pacific both in the models and in the observations. The movement and the narrow meridional width of the significant positive correlations suggest changes of the thermocline depth and subsurface temperature anomalies in the western and central-eastern equatorial Pacific Ocean by the arrival of the Kelvin waves in the summer and fall of Year +1. The above comparisons suggest that the oceanic channel dynamics are reproduced by the OGCM experiments.

    The SLA lag correlations in the cold tongue in the fall of Year +1 are stronger and over a larger longitudinal domain in the OFES model, but weaker in the HYCOM model, than those in the observations.The diff erences may be due to errors of model parameterizations of the diff usion and sub-grid processes, which aff ect the amplitudes of the Kelvin wave propagation.

    The Kelvin wave propagation appears independent of the ENSO events, because the lag correlations in the equatorial band of the Pacific Ocean remain significant even if the ENSO signals associated with the Ni?o3.4 index are removed. Partial correlation analyses show essentially the same propagation of the Kelvin waves associated with the IOD events (not shown). These analyses suggest that the oceanic channel dynamics are independent of the oceanatmosphere coupling within the tropical Pacific basin.

    3.2 Subsurface lag correlations

    The subsurface temperature anomalies are dominated by the thermocline movement and can be used to examine the Kelvin wave propagation. The simulated lag correlations between the surface temperature anomalies in STIO in fall and the subsurface temperature anomalies in the equatorial Pacific vertical section in the following winter through fall show an area of significant positive correlations propagating from the subsurface equatorial western Pacific to the central-eastern equatorial Pacific,consistent with the observational analyses (Fig.3).The oceanic channel dynamics are suggested by the propagation of the lag correlations in the subsurface equatorial Pacific Ocean both in the OGCM simulations and in the observations. The correlations also suggest that the subsurface anomalies reach the surface in the eastern Pacific cold tongue in the summer through fall of Year +1 to induce further ocean-atmosphere coupled evolution. The potential predictability of ENSO beyond the spring barrier apparently lies in the subsurface connections across the equatorial Pacific basin.

    The propagation of the subsurface lag correlations remains unchanged even if the ENSO signals are removed from the anomalies (figure omitted),suggesting the independence of the oceanic channel dynamics from the tropical ocean-atmosphere coupling both in the model and in the observations.

    Fig.2 Comparisons of the observed and simulated lag partial correlations with the ENSO signals removed by OFES and HYCOM between the STIO SLA in fall and the tropical Indo-Pacific Ocean SLA in the following winter, spring,summer and fall

    3.3 Eff ects of the atmospheric bridge

    Lag correlations between the surface zonal wind anomalies (SZWA) over the western equatorial Pacific(130°–150°E, 5°S–5°N) in fall and the SLA over the Indian and Pacific Oceans in the following winter through fall are calculated to study the atmospheric bridge. The calculations in the observations and in the OFES and HYCOM model simulations are compared in Fig.4. The use of the SZWA instead of the DMI in the lag correlation analyses is to avoid incorporating the STIO upwelling anomalies associated with the IOD in the analysis. The results suggest that the NCEP/NCAR wind anomalies and the altimeter SLA in the cold tongue in the eastern equatorial Pacific are not significantly correlated at the time lag of one year,suggesting that the IOD precursory signal is not carried to the Pacific Ocean by the atmospheric bridge. The same insignificant lag correlations are reproduced in the OFES and HYCOM models, except that the simulated lag correlations are larger than in the observations. The diff erences suggest model defects.

    The above results of the lag correlation analyses are consistent with the extensive investigations of the oceanic channel dynamics in the observational data and in the numerical experiments in the published literature (Yuan et al., 2013, 2016; Xu et al., 2013,2015). Those analyses are not repeated in this paper.

    Fig.3 Comparisons of the lag partial correlation with the ENSO signals removed between STIO surface temperature anomalies in fall and the subsurface temperature anomalies in the Pacific equatorial vertical section in the following winter, spring, summer, and fall in the observations (left), OFES (middle), and HYCOM (right)

    4 LCSM EXPERIMENT

    4.1 Experiments with an idealized barrier

    A patch of wind stress is specified in the western basin of the LCSM over the domain of (5°S–5°N,65°–95°E) with cosine and Gaussian profiles in the zonal and meridional directions, respectively, which oscillates at various periods from 60 to 720 days with an amplitude of 0.15 N/m2. Each experiment is integrated for 4 years to reach the periodic steady state (owing to no reflections at the western and eastern boundaries). The penetration rate of the Kelvin wave is defined as the ratio of the sea level in the eastern basin (150°E, 0°) in comparison to the sea level in the western basin (95°E, 0°) in the experiment without the barrier, which represents the Kelvin wave sea level in the western basin without the barrier reflections.

    Experiments are conducted with the barrier extending from the northern boundary to the equator,4°S, 8°S, and 12°S, separately. The sea level amplitudes of the first baroclinic mode and of the total multimode simulations are compared in Fig.5. When the barrier extends to the equator, half of the wave energy propagates into the eastern basin due to symmetry with respect to the equator (Fig.5a, b). As the barrier extends further to the south, the sea level amplitudes in the eastern basin decrease (Fig.5c–h).This is because part of the Kelvin waves in the southern hemisphere is reflected into equatorial Rossby waves at the barrier, decreasing the sea level amplitudes reaching the eastern basin. The multimode reflection at the barrier is slightly higher than that of the first baroclinic mode waves.

    The dependence of the penetration rate on the barrier length and on the periods of the forcing is investigated (Fig.6). When the barrier southern tip is at the equator, about 73% (equivalent to 50% of the wave energy) of the first baroclinic mode sea level reaches the eastern basin. The penetration rate for the multimode sea level is slightly smaller due to the larger diff usion on the higher modes. If the barrier extends to 8°S, about 22% of the first baroclinic mode sea level and 12%–15% of the multimode sea level amplitudes depending on the periods propagate to the eastern basin. If the barrier extends beyond 12°S, less than 15% (20% for periods shorter than 60 days) of the first baroclinic mode sea level and less than 10%of the multimode sea level amplitudes will be able to penetrate into the eastern basin. The penetration rates of the first baroclinic mode are nearly independent of the periods of the waves, whereas those of the multimode are more dispersive with the periods,especially when the barrier extension is in the vicinity of the equator (Fig.6). This is due to the downward propagation of Kelvin waves that is dependent on the frequency of the waves. The higher the frequency, the steeper downward propagation of the waves is.

    The penetration rate of the interannual Kelvin wave sea level is found only sensitive to the latitude of the barrier southern tip, notwithstanding the shape of the barrier. In another set of the experiments, in which the barrier south of the equator is set to be at an angle of 45° from the equator, the penetration rate of the sea level amplitude is nearly the same as those in the experiments with a straight barrier, except for the intraseasonal waves, which are not sensitive to the barrier southern latitude beyond about 6°S due to the critical latitude at about 7°S. (Fig.7). The sole dependence of the penetration rate on the latitude ofthe barrier southern tip is further evidenced by the experiments with an L-shape barrier extending to diff erent southern latitudes (Fig.8). The penetration rate of the Kelvin sea level is nearly independent of the zonal barrier length for the interannual waves. In comparison, the multimode intraseasonal penetration rate is sensitive to the length of the barrier because of vertical propagation of these waves in the course of the propagation.

    Fig.4 Comparison of lag partial correlation with the ENSO signals removed between SZWA over the far western equatorial Pacific Ocean in fall and the tropical Indo-Pacific Ocean SLA in the following winter, spring, summer, and fall in the observations (left), OFES (middle), and HYCOM (right)

    When the L-shape barrier has a gap in the zonal part, part of the Kelvin wave energy will penetrate into the eastern basin through the gap. However, for the interannual and lower-frequency Kelvin waves,the penetration rate of the Kelvin wave sea level is not sensitive to the width of the gap (Fig.9), because the Kelvin wave speed is very fast so that the waves through the gap and those going around the end of the barrier interfere with each other constructively. For intraseasonal waves, these waves may interfere with each other destructively due to the phase diff erence of the diff erent routes of the propagation, hence the sensitivity to the partition of the wave energy through the gap. To examine the eff ects of the barrier length,yet another set of the experiments is conducted with diff erent lengths of the barrier east of the gap of 0.3°wide (Fig.10). The penetration rates of the intraseasonal Kelvin waves decrease significant with the eastern barrier length, suggesting destructive interference.

    Fig.5 The sea level amplitude in the experiments using the first baroclinic mode (a, c, e, g) and LCSM (b, d, f, h) with a barrier extending from the northern boundary to the equator (a, b), 4°S (c, d), 8°S (e, f), and 12°S (g, h)

    4.2 Experiments with realistic geometry of the Indonesian seas

    Fig.6 The penetration rates depending on the barrier extensions ( x-coordinate) and on the periods of the forcing (line type) for the first baroclinic mode and the multimode sea level amplitude

    Fig.7 The penetration rates depending on the barrier extensions with an angle of 45° ( x-coordinate) and on the periods of the forcing for the first baroclinic mode and the multimode sea level amplitude

    Fig.8 The penetration rates depending on the L-shape barrier extensions in degrees ( x-coordinate) and on the periods of the forcing for the first baroclinic mode and the multimode sea level amplitude

    Fig.9 The penetration rates depending on the width of the gap in the L-shape barrier ( x-coordinate) and on the periods of the forcing for the first baroclinic mode and the multimode sea level amplitude

    The LCSM has also been deployed over the tropical eastern Indian and western Pacific Oceans between 60°E and 160°E and between 20°S and 20°N with a high horizontal resolution of 0.1° longitude by 0.1°latitude to resolve the main passages of the Indonesian seas. The open and closed boundary conditions and the wind forcing of the experiments are the same as in the experiments with the idealized barrier. The results of the experiments have shown that the Indian Ocean Kelvin waves propagate to the western equatorial Pacific Ocean through both the Makassar Strait and the eastern Indonesian seas (Fig.11). The magnitudes of the Kelvin waves in the Makassar Strait are much larger than those in the eastern seas. These Kelvin waves must propagate across the Maluku Channel to reach the western Pacific Ocean. The numerical experiment in Fig.11a suggests clearly that the Makassar Strait Kelvin waves lose substantial amplitudes after they propagate across the Maluku Strait.

    The total penetration rate of sea level amplitudes in these experiments are estimated relative to the area averaged SLA over the STIO. The total penetration rate of the Kelvin wave sea level amplitudes at the interannual time scales is about 9%–15%, depending on the vertical modes (Fig.12). The penetration rate is nearly unchanged if the Lombok Strait or both the Flores Sea and the Lombok Strait are closed. The closing of the Flores Sea helps the Kelvin waves propagate into the Makassar Strait, hence the slight increase of the penetration rate. The closing of the Halmahera Sea appears to have little eff ects on the penetration rate. Some of the incoming Kelvin wave energy propagates back into the Indian Ocean along the West Australia coasts. When the Australia continent is removed, the penetration rate is increased by 3%–4%. These experiments using realistic coastal lines agree with those using the idealized barrier in the penetrating rate that around 9%–15% of the Kelvin wave sea level reaching the western equatorial Pacific Ocean. The Kelvin wave sea level amplitude is expected to be amplified by the shoaling thermocline in the eastern equatorial Pacific, which will induce ocean-atmosphere coupling to further amplify the oceanic channel signals in the year to follow.

    Figure 12 also shows the sea level penetration rates in the above experiments for diff erent forcing periods.For the first baroclinic mode waves, the Australia continent results in ~3% interannual sea level amplitudes propagating back to the Indian Ocean. For intraseasonal waves, this percentage is as high as 10%. The penetration rates for the intraseasonal waves are much higher due to the cancellation of the sea level in the STIO of the incoming and the reflected waves near the eastern boundary of the Indian Ocean.The closing of the Lombok Strait, the Flores Sea, and the Halmahera Straits has little eff ects on the penetration rate of the interannual waves (about 1%).For multimode sea level amplitudes, the penetration rate reduced by the Australia continent is about 3%for the interannual waves. The closing of the internal straits of the Indonesian seas results in little change of the penetration rate. The penetration rate of the multimode sea level through the realistic Indonesian sea channels is about 10% for the interannual Kelvin waves.

    It is worth mentioning that the depth of the LCSM is set at a constant of 4 000 m, which is not realistic.The interannual waves are expected to propagate through the major channels of the Indonesian seas,unblocked by the sill depth due to the flat ray paths of their vertical propagation, which suggests that the estimated penetration rates are robust. For intraseasonal and higher frequency waves, some of the waves may not be able to penetrate through all of the channels due to the depth limitation of the shallow straits. In addition, the Kelvin waves can be reflected at the surface and bottom (McCreary, 1984), but for the interannual waves, the energy will be trapped primarily above the pycnocline (Moore and McCreary,1990).

    Fig.10 The penetration rates depending on the length of the L-shape barrier east of a 0.3°gap ( x-coordinate) and on the periods of the forcing for the first baroclinic mode and the multimode sea level amplitude

    5 THE CENTENNIAL TRENDS OF THE IOD-ENSO TELECONNECTION

    Using the HadISST, ERSST, and Kaplan SST data,the lag correlations of the STIO SSTA in fall and the eastern equatorial Pacific Ni?o3.4 SSTA in the next fall are calculated in the running windows of 5, 11,and 21 years. The centennial trends of the teleconnection based on the 5-year running window are calculated to be 0.004/a, 0.004/a, and 0.005/a,respectively, for the three datasets (Fig.13). The centennial trends of the correlations after the filtering of the 11- and 21-year running windows are similar(Table 1). All calculations present positive trends,suggesting that the IOD-ENSO precursory teleconnection is getting stronger over the past 100 years or so. This kind of lag correlation calculations is repeated for the CMIP5 Historical simulated SSTA and SLA. Most of the CMIP5 simulations have generated positive trends of the lag correlations between the STIO SSTA in fall and the cold tongue SSTA in the next fall over the past 100 years or so,consistent with the analysis of the observed SSTA(Fig.14). So are the SLA lag correlations, suggesting the oceanic channel dynamics behind the trends. A few CMIP5 models have presented negative trends,which may be due to model deficiency. All except the CanESM2 model show the same trends for the SSTA and SLA lag correlations, suggesting that the IODENSO precursory teleconnection is controlled by the Kelvin wave propagation dynamics in these CMIP5 simulations.

    Fig.11 The sea level amplitude in the experiments using the LCSM with realistic coastal lines of the Indonesian seas

    Table 1 Trends of the lag correlations between the STIO SSTA and the Ni?o3.4 index after the filtering of the 5-, 11-, and 21-year running windows

    The dynamics behind the increased IOD-ENSO teleconnection are associated with the ITF changes.Negative trends of ITF transports, suggesting decreased ITF transports from the Pacific to the Indian Ocean over the past 100 years or so, have been produced in the majority of the CMIP5 simulations(Fig.15). The decreased ITF transports are conducive for Kelvin wave propagation from the Indian Ocean to the western equatorial Pacific Ocean, considering the Doppler eff ects and the nonlinear diff usion. In comparison, the SZWA in the far western equatorial Pacific have shown negative trends in most of the CMIP5 simulations in the past 100 year or so,suggesting weaker connection through the atmospheric bridge (Fig.16). The change in the atmospheric bridge is inconsistent with the SSTA lag correlation trends.

    Comparisons of the trends of the SZWA and the ITF transport anomalies in the Historical and HistoricalNat simulations of the CMIP5 models suggest that the anthropogenic forcing on global warming has significant impact on the oceanic channel and the atmospheric bridge, hence on the IOD-ENSO teleconnection (Figs.16 and 17). The detailed dynamics of this impact is beyond the scope of this study.

    Fig.12 The penetration rates depending on the periods of the forcing in the control experiments of the LCSM with some of the Indonesian straits closed( x-coordinate) for the first baroclinic mode and the multimode sea level amplitude

    6 CONCLUSION

    Lag correlations between the SLA in STIO in fall and the Indo-Pacific SLA in the following winter through fall seasons are calculated in the satellite altimeter sea level data of 1993–2008 and in the highresolution simulations of the OFES and HYCOM models for the period of 1990–2008. The propagation of significant positive lag correlations from STIO to the equatorial Pacific Ocean through the Indonesian seas in the observations is simulated successfully by the OGCMs. The dynamics are due to the propagation of the Kelvin waves, which induces ITF transport anomalies and subsurface temperature anomalies in the western equatorial Pacific Ocean. The propagation of subsurface temperature anomalies from the western equatorial Pacific Ocean to the eastern equatorial Pacific cold tongue is evidenced by the significant lag correlations between the SSTA in STIO in fall and the Indo-Pacific subsurface temperature in the equatorial vertical section of the Pacific Ocean in the following winter through fall, the observational analyses of which are also simulated successfully by the OGCMs.The elevation of the subsurface temperature anomalies to the cold tongue surface suggests that oceanatmosphere coupling amplifies the oceanic channel signals from the tropical Indian Ocean. All these dynamics are also simulated well by the OGCMS,providing evidence of the dynamical link between the tropical Indian Ocean and the tropical Pacific Ocean through the ITF. In contrast, the lag correlations between the SZWA over the western equatorial Pacific in fall and the Indo-Pacific SLA in the following spring through fall are all insignificant in the equatorial Pacific Ocean, suggesting that the atmospheric bridge is unlikely the link between the IOD and ENSO at the one year time lag. Evidently, the precursory relation between IOD and ENSO at the one-year time lag is due to the oceanic channel dynamics through the Indonesian seas, namely the ITF variability.

    Fig.13 Lag correlations between the southeastern tropical Indian Ocean SSTA in fall and the cold tongue SSTA at the one-year time lag in the past 140 years based on the HadISST (a), ERSST (b), and Kaplan (c) SST data

    Fig.14 Comparison of the long-term trends of the SSTA and SLA lag correlations between the STIO and the cold tongue at a one-year time lag in the CMIP5 Historical experiments

    Fig.15 Comparison of the long-term trends of the ITF transport anomalies and of the SSTA lag correlations between the STIO and the cold tongue at a one-year time lag in the CMIP5 Historical experiments

    Fig.16 Comparison of the western Pacific SZWA trends in the Historical and the HistoricalNat experiments

    Fig.17 Comparison of the ITF transport anomaly trends in the Historical and HistoricalNat experiments

    The penetration rates for the Kelvin waves from the equatorial Indian Ocean to the equatorial Pacific Ocean are estimated based on experiments using a LCSM. The penetration rate of the interannual Kelvin waves is found to be dependent only on the southern tip latitude of the barrier between the two equatorial basins, and independent of the lengths and gaps of the barrier. The experiments with realistic coastal lines of the Indonesian seas suggest that the Kelvin waves propagate into the equatorial western Pacific Ocean through both the Makassar Strait and the eastern Indonesian seas. The total penetration rates for the sea level amplitudes are about 10%–15%, depending on the baroclinic modes and not sensitive to the closing of the Lombok Strait, the Flores Sea, or the Halmahera Sea. The experiments of the LCSM and the diagnosis of the OGCM simulations are consistent in that the Kelvin waves from the equatorial Indian Ocean can propagate through the Indonesian seas to the equatorial western Pacific Ocean to influence the ENSO evolution in the year to follow, which is the essence of the oceanic dynamics of the IOD-ENSO precursory relation.

    The IOD-ENSO precursory teleconnection is found to have an increasing trend over the past 100 years or so. The diagnoses of the CMIP5 simulations have shown positive trends of the lag correlations between the STIO SSTA (SLA) and the cold tongue SSTA (SLA) at the one year time lag over the latest 100 years or so, suggesting positive trends of the IOD-ENSO teleconnection associated with the global warming. The centennial trends are consistent with the negative trends of the ITF transport, but are in contrast to the negative trends of the SZWA over the far western equatorial Pacific. The weaker ITF transports from the Pacific to the Indian Ocean facilitate the propagation of the Kelvin waves through the Indonesian seas. In comparison, the negative trends of the SZWA cannot explain the positive IODENSO teleconnection trends and suggest insignificant roles played by the atmospheric bridge. The IODENSO teleconnection is found to be sensitive to global warming forcing, as suggested by the comparisons of the ITF transport trends and the SZWA trends in the Historical and HistoricalNat experiments.

    Alexander M A, Bladé I, Newman M, Lanzante J R, Lau N C,Scot J D. 2002. The atmospheric bridge: the influence of ENSO teleconnections on air-sea interaction over the global oceans. J. Climate, 15(16): 2 205-2 231.

    Annamalai H, Xie S P, McCreary J P, Murtugudde R. 2005.Impact of Indian Ocean sea surface temperature on developing El Ni?o. J. Climate, 18(2): 302-319.

    Behera S K, Luo J J, Masson S, Rao S A, Sakuma H, Yamagata T. 2006. A CGCM study on the interaction between IOD and ENSO. J. Climate, 19(9): 1 688-1 705.

    Behera S K, Yamagata T. 2003. Influence of the Indian Ocean dipole on the southern oscillation. J. Meteor. Soc. Jpn.,81(1): 169-177.

    Bleck R. 2002. An oceanic general circulation model framed in hybrid isopycnic-Cartesian coordinates. Ocean Modell.,4(1): 55-88.

    Clarke A J, Gorder S V. 2003. Improving El Ni?o prediction using a space-time integration of Indo-Pacific winds and equatorial Pacific upper ocean heat content. Geophys. Res.Lett., 30(7): 1 399, https://doi.org/10.1029/2002GL016673.

    Drushka K, Sprintall J, Gille S T, Brodjonegoro I. 2010.Vertical structure of kelvin waves in the Indonesian throughflow exit passages. J. Phys. Oceanogr., 40(9):1 965-1 987.

    Halliwell G, Bleck R, Chassignet E. 1998. Atlantic Ocean simulations performed using a new hybrid-coordinate ocean model. In EOS Transactions, American Geophysical Union ( AGU), Fall 1998 Meeting, San Francisco, CA.

    Izumo T, Lengaigne M, Vialard J, Luo J J, Yamagata T, Madec G. 2014. Influence of Indian Ocean dipole and pacific recharge on following year’s El Ni?o: interdecadal robustness. Climate Dyn., 42(1-2): 291-310.

    Izumo T, Vialard J, Lengaigne M, de Boyer Montegut C,Behera S K, Luo J J, Cravatte S, Masson S, Yamagata T.2010. Influence of the state of the Indian Ocean Dipole on the following year’s El Ni?o. Nat. Geosci., 3(3): 168-172.

    Kalnay E, Kanamitsu M, Kistler R, Collins W, Deaven D,Gandin L, Iredell M, Saha S, White G, Woollen J, Zhu Y,Leetmaa A, Reynolds R, Chelliah M, Ebisuzaki W,Higgins W, Janowiak J, Mo K C, Ropelewski C, Wang J.1996. The NCEP/NCAR 40-year reanalysis project. Bull.Am. Meteor. Soc., 77(3): 437-471.

    Kaplan A, Cane M A, Kushnir Y, Clement A C, Blumenthal M B, Rajagopalan B. 1998. Analyses of global sea surface temperature 1856-1991. J. Geophy. Res., 103(C9):18 567-18 589.

    Klein S A, Soden B J, Lau N C. 1999. Remote sea surface temperature variations during ENSO: evidence for a tropical atmospheric bridge. J. Climate, 12(4): 917-932.

    Kug J S, Li T, An S I, Kang I S, Luo J J, Masson S, Yamagata T. 2006. Role of the ENSO-Indian Ocean coupling on ENSO variability in a coupled GCM. Geophys. Res. Lett.,33(9): L09710, https://doi.org/10.1029/2005GL024916.

    Lau N C, Leetmaa A, Nath M J, Wang H L. 2005. Influences of ENSO-induced Indo-western Pacific SST anomalies on extratropical atmospheric variability during the boreal summer. J. Climate, 18(15): 2 922-2 942.

    Lau N C, Nath M J. 2003. Atmosphere-ocean variations in the Indo-Pacific sector during ENSO episodes. J. Climate,16(1): 3-20.

    Luo J J, Zhang R C, Behera S K, Masumoto Y, Jin F F, Lukas R, Yamagata T. 2010. Interaction between El Ni?o and Extreme Indian Ocean Dipole. J. Climate, 23(3): 726-742.

    Masumoto Y, Sasaki H, Kagimoto T, Komori N, Ishida A,Sasai Y, Miyama T, Motoi T, Mitsudera H, Takahashi K,Sakuma H, Yamagata T. 2004. A fifty-year eddy-resolving simulation of the world ocean—preliminary outcomes of OFES (OGCM for the Earth Simulator). J. Earth Simulator, 1: 35-56.

    McCreary Jr J P. 1984. Equatorial beams. J. Mar. Res., 42(2):395-430.

    Meehl G A, Bony S. 2011. Introduction to CMIP5. Clivar Exchanges, 16(2): 4-5.

    Molcard R, Fieux M, Syamsudin F. 2001. The throughflow within Ombai Strait. Deep Sea Res. Part I Oceanogr. Res.Papers, 48(5): 1 237-1 253.

    Moore D W, McCreary J P. 1990. Excitation of intermediatefrequency equatorial waves at a western ocean boundary:with application to observations from the Indian Ocean. J.Geophys. Res., 95(C4): 5 219-5 231.

    Pujiana K, Gordon A L, Sprintall J, Susanto R D. 2009.Intraseasonal variability in the Makassar strait thermocline. J. Mar. Res., 67(6): 757-777.

    Rayner N A, Parker D E, Horton E B, Folland C K, Alexander L V, Rowell D P, Kent E C, Kaplan A. 2003. Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century. J.Geophys. Res., 108(D14): 4 407, https://doi.org/10.1029/2002JD002670.

    Saji N H, Goswami B N, Vinayachandran P N, Yamagata T.1999. A dipole mode in the tropical Indian Ocean. Nature,401(6751): 360-363.

    Slutz R J, Lubker S J, Hiscox J D, Woodruff S D, Jenne R L,Joseph D H, Steurer P M, Elms J D. 1985. Comprehensive Ocean-Atmosphere Data Set: Release 1. NOAA Environmental Research Laboratories, Climate Research Program, Boulder, Colorado. 268p.

    Smith T M, Reynolds R W, Peterson T C, Lawrimore J. 2008.Improvements to NOAA’s historical merged Land-Ocean surface temperature analysis (1880-2006). J. Climate,21(10): 2 283-2 296.

    Susanto R D, Ffield A, Gordon A L, Adi T R. 2012. Variability of Indonesian throughflow within Makassar Strait, 2004-2009. J. Geophys. Res., 117(C9): C09013, https://doi.org/10.1029/2012JC008096.

    Trenary L L, Han W Q. 2012. Intraseasonal-to-interannual variability of South Indian Ocean sea level and thermocline: remote versus local forcing. J. Phys.Oceanogr., 42(4): 602-627.

    Trenary L L, Han W Q. 2013. Local and remote forcing of decadal sea level and thermocline depth variability in the south Indian Ocean. J. Geophys. Res., 118(1): 381-398,https://doi.org/10.1029/2012JC008317.

    White W B. 1995. Design of a global observing system for gyre-scale upper ocean temperature variability. Prog.Oceanogr., 36(3): 169-217.

    Wu G X, Meng W. 1998. Gearing between the Indo-monsoon Circulation and the Pacific-Walker Circulation and the ENSO Part I: data analyses. Sci. Atmos. Sinica, 22(4):470-480. (in Chinese with English abstract)

    Wu R G, Kirtman B P. 2004. Understanding the impacts of the Indian Ocean on ENSO variability in a coupled GCM. J.Climate, 17(20): 4 019-4 031.

    Xu T F, Yuan D L, Yu Y Q, Zhao X. 2013. An assessment of Indo-Pacific oceanic channel dynamics in the FGOALS-g2 coupled climate system model. Adv. Atmos. Sci., 30(4):997-1 016, https://doi.org/10.1007/s00376-013-2131-2.

    Xu T F, Yuan D L. 2015. Why does the IOD-ENSO teleconnection disappear in some decades? Chin. J.Oceanol. Limnol., 33(2): 534-544, https://doi.org/10.1007/s00343-015-4044-7.

    Yuan D L, Han W Q. 2006. Roles of equatorial waves and western boundary reflection in the seasonal circulation of the equatorial Indian Ocean. J. Phys. Oceanogr., 36(5):930-944.

    Yuan D L, Wang J, Xu T F, Xu P, Hui Z, Zhao X, Luan Y H,Zheng W P, Yu Y Q. 2011. Forcing of the Indian Ocean dipole on the interannual variations of the tropical Pacific Ocean: roles of the Indonesian throughflow. J. Climate,24(14): 3 593-3 608.

    Yuan D L, Xu P, Xu T F. 2017. Climate variability and predictability associated with the Indo-Pacific Oceanic channel dynamics in the CCSM4 coupled system model.Chin. J. Oceanol. Limnol., 35(1): 23-28, https://doi.org/10.1007/s00343-016-5178-y.

    Yuan D L, Zhou H, Zhao X. 2013. Interannual climate variability over the tropical Pacific Ocean induced by the Indian Ocean dipole through the Indonesian throughflow.J. Climate, 26(9): 2 845-2 861.

    Zhou Q, Duan W S, Mu M, Feng R. 2015. Influence of positive and negative Indian Ocean Dipoles on ENSO via the Indonesian Throughflow: results from sensitivity experiment. Adv. Atmos. Sci., 32(6): 783-793.

    猜你喜歡
    徐鵬
    Enhanced measurement precision with continuous interrogation during dynamical decoupling
    我院老年患者多重用藥情況調查
    《喜悅》背后的烈火底色
    徐鵬飛漫畫
    Stabbing Pain with Words
    明文林郎徐鵬、傅氏夫婦墓志研究
    那年杏花微雨,我跟小三分住樓上樓下
    杏樹莊子
    伊犁河(2020年2期)2020-10-09 10:50:58
    用青春和汗水鑄就忠誠
    聽說第七街有人等你
    伴侶(2018年11期)2018-11-22 03:18:50
    色哟哟哟哟哟哟| 欧美区成人在线视频| 白带黄色成豆腐渣| 精品少妇黑人巨大在线播放 | 中出人妻视频一区二区| 啦啦啦啦在线视频资源| 久久久欧美国产精品| 变态另类成人亚洲欧美熟女| 婷婷色综合大香蕉| 在线观看美女被高潮喷水网站| 日本一二三区视频观看| 中文字幕人妻熟人妻熟丝袜美| 国产免费男女视频| 日韩精品青青久久久久久| 亚洲自偷自拍三级| 午夜福利在线观看吧| 久久久国产成人精品二区| 国内精品久久久久精免费| 夫妻性生交免费视频一级片| 成人三级黄色视频| 欧美性猛交╳xxx乱大交人| 一本久久中文字幕| 老女人水多毛片| 亚洲在久久综合| 国产精品综合久久久久久久免费| 亚洲内射少妇av| 亚洲精品国产av成人精品| 欧美极品一区二区三区四区| 国产伦理片在线播放av一区 | 国产成人a∨麻豆精品| 国产精品一二三区在线看| 国产探花在线观看一区二区| 91久久精品电影网| 色哟哟·www| 久久久久九九精品影院| 精品久久久久久久久久久久久| 一个人看视频在线观看www免费| 精品人妻偷拍中文字幕| 亚洲无线在线观看| 日本撒尿小便嘘嘘汇集6| 国产欧美日韩精品一区二区| 天天躁夜夜躁狠狠久久av| 亚洲一级一片aⅴ在线观看| 男女视频在线观看网站免费| 男人狂女人下面高潮的视频| 国产精品乱码一区二三区的特点| av国产免费在线观看| 最近视频中文字幕2019在线8| 国产精品一区www在线观看| 国产成人a区在线观看| 人体艺术视频欧美日本| 丰满人妻一区二区三区视频av| 一个人观看的视频www高清免费观看| 成人鲁丝片一二三区免费| 亚洲欧美成人精品一区二区| 亚洲自拍偷在线| 床上黄色一级片| 国内精品宾馆在线| 蜜桃久久精品国产亚洲av| 熟妇人妻久久中文字幕3abv| av天堂在线播放| 国产一区二区在线av高清观看| 禁无遮挡网站| 男女那种视频在线观看| 精品熟女少妇av免费看| 美女黄网站色视频| 亚洲va在线va天堂va国产| 亚洲不卡免费看| 免费无遮挡裸体视频| 欧美日韩精品成人综合77777| 国产黄色小视频在线观看| 国产精品一二三区在线看| 亚洲精华国产精华液的使用体验 | 成人二区视频| 不卡视频在线观看欧美| 久久久成人免费电影| 国产成人a∨麻豆精品| 99热这里只有精品一区| 国产高清不卡午夜福利| 黑人高潮一二区| 性色avwww在线观看| 99久久成人亚洲精品观看| h日本视频在线播放| 久久精品久久久久久噜噜老黄 | 精品99又大又爽又粗少妇毛片| 久久久久国产网址| 搡女人真爽免费视频火全软件| 亚洲av熟女| 嘟嘟电影网在线观看| 不卡视频在线观看欧美| 久久久久久久久大av| 一边摸一边抽搐一进一小说| 国产探花在线观看一区二区| 2021天堂中文幕一二区在线观| 国产精品爽爽va在线观看网站| 精品无人区乱码1区二区| 成人午夜精彩视频在线观看| 日本与韩国留学比较| 我要看日韩黄色一级片| 久久久久性生活片| 欧美激情久久久久久爽电影| 久久婷婷人人爽人人干人人爱| 中国美女看黄片| 日韩欧美国产在线观看| 99久久久亚洲精品蜜臀av| 国产成人精品婷婷| 麻豆乱淫一区二区| 国产白丝娇喘喷水9色精品| 亚洲欧洲国产日韩| 日本欧美国产在线视频| 国产午夜精品一二区理论片| eeuss影院久久| 天天躁日日操中文字幕| 国产日韩欧美在线精品| 日本免费一区二区三区高清不卡| 最近中文字幕高清免费大全6| 九九爱精品视频在线观看| 黄色配什么色好看| 99热网站在线观看| 99久久精品国产国产毛片| 国产亚洲av片在线观看秒播厂 | 国产在视频线在精品| 美女脱内裤让男人舔精品视频 | 女的被弄到高潮叫床怎么办| 在线a可以看的网站| 国产精品国产三级国产av玫瑰| 午夜福利视频1000在线观看| 中文字幕av在线有码专区| 日日啪夜夜撸| av卡一久久| 一本精品99久久精品77| 22中文网久久字幕| 国产精品精品国产色婷婷| 欧美3d第一页| 变态另类成人亚洲欧美熟女| 一级毛片电影观看 | 精品久久久久久成人av| 悠悠久久av| 国国产精品蜜臀av免费| 一卡2卡三卡四卡精品乱码亚洲| 波野结衣二区三区在线| 高清日韩中文字幕在线| 精品人妻熟女av久视频| 大香蕉久久网| 免费电影在线观看免费观看| 一级黄片播放器| 亚洲性久久影院| av免费观看日本| 99久久无色码亚洲精品果冻| 欧美一级a爱片免费观看看| 日本黄大片高清| 九九热线精品视视频播放| 久久久国产成人精品二区| 国产精品.久久久| 高清日韩中文字幕在线| 国产视频内射| 精品不卡国产一区二区三区| 狂野欧美激情性xxxx在线观看| 国产一级毛片在线| 夜夜看夜夜爽夜夜摸| 一边亲一边摸免费视频| 国产69精品久久久久777片| 老熟妇乱子伦视频在线观看| 国产精品,欧美在线| 麻豆成人午夜福利视频| 国产探花在线观看一区二区| 日韩,欧美,国产一区二区三区 | 成人亚洲精品av一区二区| 午夜福利在线观看吧| 国产成人影院久久av| 非洲黑人性xxxx精品又粗又长| 国产亚洲av片在线观看秒播厂 | 五月玫瑰六月丁香| 一级毛片电影观看 | a级毛片免费高清观看在线播放| 18禁在线播放成人免费| 在现免费观看毛片| 国产av一区在线观看免费| 久久久久久久久久成人| ponron亚洲| av免费在线看不卡| 欧美zozozo另类| 日本黄色片子视频| 久久久欧美国产精品| 美女脱内裤让男人舔精品视频 | 国产大屁股一区二区在线视频| 国产精品一区二区性色av| 成人特级av手机在线观看| 国产一级毛片在线| 久久久久久九九精品二区国产| 国产黄色小视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 久久热精品热| 日韩欧美精品免费久久| 国产伦一二天堂av在线观看| 亚州av有码| 国产亚洲欧美98| h日本视频在线播放| 国产av不卡久久| 干丝袜人妻中文字幕| 国产高清不卡午夜福利| 一级毛片久久久久久久久女| 99久久人妻综合| 欧美一区二区精品小视频在线| 99九九线精品视频在线观看视频| 色尼玛亚洲综合影院| 边亲边吃奶的免费视频| 日韩av在线大香蕉| 欧美最黄视频在线播放免费| 久久99热这里只有精品18| 久久久国产成人免费| 麻豆成人av视频| 国产精品免费一区二区三区在线| 国产v大片淫在线免费观看| 日本-黄色视频高清免费观看| 最近视频中文字幕2019在线8| 欧美高清性xxxxhd video| 亚洲一区二区三区色噜噜| 少妇裸体淫交视频免费看高清| 老师上课跳d突然被开到最大视频| 国产爱豆传媒在线观看| 成人无遮挡网站| 小蜜桃在线观看免费完整版高清| 色吧在线观看| 欧美另类亚洲清纯唯美| 国产精品久久久久久亚洲av鲁大| 亚洲最大成人手机在线| 久久精品国产亚洲av涩爱 | 人人妻人人看人人澡| 久久精品国产亚洲av香蕉五月| 日韩 亚洲 欧美在线| 看免费成人av毛片| 国产视频首页在线观看| 国产色爽女视频免费观看| 亚洲自拍偷在线| 春色校园在线视频观看| 国产成人freesex在线| 亚洲自偷自拍三级| 亚洲精品乱码久久久v下载方式| 天天躁夜夜躁狠狠久久av| 毛片一级片免费看久久久久| 久久亚洲国产成人精品v| 我的女老师完整版在线观看| 久久久国产成人精品二区| 亚洲欧美精品自产自拍| 日本免费a在线| 国产色爽女视频免费观看| 中文欧美无线码| 午夜老司机福利剧场| 免费观看精品视频网站| 色噜噜av男人的天堂激情| 能在线免费观看的黄片| 高清午夜精品一区二区三区 | 女同久久另类99精品国产91| 超碰av人人做人人爽久久| 国产高清三级在线| 国产 一区精品| 午夜精品在线福利| 亚洲经典国产精华液单| 淫秽高清视频在线观看| 欧美最黄视频在线播放免费| 成人国产麻豆网| 99热6这里只有精品| 国产亚洲精品久久久久久毛片| 中文字幕熟女人妻在线| 亚洲欧美日韩高清在线视频| 久久久久久久久大av| 国产69精品久久久久777片| 色哟哟哟哟哟哟| 91久久精品国产一区二区成人| 一级毛片aaaaaa免费看小| 26uuu在线亚洲综合色| 国产麻豆成人av免费视频| 精品不卡国产一区二区三区| 免费人成视频x8x8入口观看| 国产午夜精品论理片| 中文精品一卡2卡3卡4更新| 国产精品蜜桃在线观看 | 在线观看午夜福利视频| 白带黄色成豆腐渣| 国产精品精品国产色婷婷| 日韩视频在线欧美| 国产精品,欧美在线| 欧美成人精品欧美一级黄| 国产免费一级a男人的天堂| 女人十人毛片免费观看3o分钟| 中文字幕av在线有码专区| 国产精品国产三级国产av玫瑰| 精品日产1卡2卡| ponron亚洲| 高清毛片免费观看视频网站| 欧美人与善性xxx| 亚洲欧洲日产国产| 国产综合懂色| 免费av毛片视频| 亚洲精品456在线播放app| 天天躁夜夜躁狠狠久久av| 插阴视频在线观看视频| 在线观看66精品国产| 亚洲av一区综合| av专区在线播放| 亚洲综合色惰| 久久久久久久久久久丰满| 2021天堂中文幕一二区在线观| 亚洲国产色片| 两个人的视频大全免费| 国内精品宾馆在线| av.在线天堂| 欧美zozozo另类| av在线亚洲专区| 麻豆国产97在线/欧美| 国产亚洲av片在线观看秒播厂 | 白带黄色成豆腐渣| 日日摸夜夜添夜夜爱| 午夜免费男女啪啪视频观看| 精品久久国产蜜桃| 欧美潮喷喷水| 国产亚洲91精品色在线| 成人美女网站在线观看视频| 亚洲av免费在线观看| 久久99热6这里只有精品| 成人高潮视频无遮挡免费网站| 亚州av有码| 久久久久久久久久久丰满| www日本黄色视频网| 美女cb高潮喷水在线观看| 欧美最黄视频在线播放免费| 亚洲18禁久久av| 国产精品伦人一区二区| 色综合色国产| 69人妻影院| 亚洲欧美成人精品一区二区| 亚洲最大成人av| 一级毛片aaaaaa免费看小| 国产精品国产三级国产av玫瑰| 九九爱精品视频在线观看| av免费观看日本| 国产亚洲欧美98| 大型黄色视频在线免费观看| 国产片特级美女逼逼视频| 熟女电影av网| 97热精品久久久久久| 免费无遮挡裸体视频| 嫩草影院新地址| 中国美白少妇内射xxxbb| 国国产精品蜜臀av免费| 小蜜桃在线观看免费完整版高清| 一级毛片电影观看 | 国产午夜精品论理片| 欧美性猛交黑人性爽| 啦啦啦观看免费观看视频高清| 国产精品免费一区二区三区在线| 十八禁国产超污无遮挡网站| 美女黄网站色视频| 国产在线精品亚洲第一网站| 在线免费观看的www视频| 别揉我奶头 嗯啊视频| 久久精品久久久久久噜噜老黄 | 伊人久久精品亚洲午夜| а√天堂www在线а√下载| 中国美白少妇内射xxxbb| a级毛片a级免费在线| 成人一区二区视频在线观看| 久久精品影院6| 精品久久久噜噜| 免费不卡的大黄色大毛片视频在线观看 | 人人妻人人澡欧美一区二区| 非洲黑人性xxxx精品又粗又长| 午夜福利高清视频| 中文字幕精品亚洲无线码一区| 日本黄大片高清| 午夜精品在线福利| 国产高潮美女av| 久久久国产成人免费| 身体一侧抽搐| 亚洲成人久久爱视频| ponron亚洲| 简卡轻食公司| 大型黄色视频在线免费观看| 少妇熟女欧美另类| 蜜臀久久99精品久久宅男| 中文字幕精品亚洲无线码一区| 大型黄色视频在线免费观看| 麻豆乱淫一区二区| 精品人妻视频免费看| 久久99热6这里只有精品| 国产av在哪里看| a级毛色黄片| av在线观看视频网站免费| 国产日本99.免费观看| 亚洲激情五月婷婷啪啪| 亚洲成a人片在线一区二区| 国产亚洲5aaaaa淫片| 人妻系列 视频| 国产精品久久久久久亚洲av鲁大| 非洲黑人性xxxx精品又粗又长| 波多野结衣高清无吗| 免费在线观看成人毛片| 国产一区二区三区av在线 | 99久久精品热视频| 尤物成人国产欧美一区二区三区| 国产精品av视频在线免费观看| 色尼玛亚洲综合影院| 国产精品久久久久久精品电影| 国产久久久一区二区三区| 成人国产麻豆网| 直男gayav资源| 色综合亚洲欧美另类图片| 中文字幕人妻熟人妻熟丝袜美| 91麻豆精品激情在线观看国产| 老女人水多毛片| 国产成年人精品一区二区| 高清在线视频一区二区三区 | 一边亲一边摸免费视频| 中文字幕av在线有码专区| 我要看日韩黄色一级片| 久久精品国产亚洲av涩爱 | 午夜福利在线观看免费完整高清在 | 日韩国内少妇激情av| 成人午夜高清在线视频| 日本-黄色视频高清免费观看| 日本欧美国产在线视频| 久久热精品热| 色视频www国产| 91aial.com中文字幕在线观看| 日本黄色片子视频| 少妇丰满av| 高清毛片免费看| 日韩欧美一区二区三区在线观看| 国产成人a区在线观看| 天堂√8在线中文| 中文字幕久久专区| 亚洲最大成人av| 久久久久久久久久成人| 一个人观看的视频www高清免费观看| 黄色日韩在线| 男人和女人高潮做爰伦理| 精品一区二区三区人妻视频| 非洲黑人性xxxx精品又粗又长| 村上凉子中文字幕在线| 国产精品伦人一区二区| 成人三级黄色视频| 国产av一区在线观看免费| 免费不卡的大黄色大毛片视频在线观看 | 一边摸一边抽搐一进一小说| 国产精品一区二区性色av| 最近2019中文字幕mv第一页| 深夜a级毛片| 蜜桃久久精品国产亚洲av| 国产精品一及| 看片在线看免费视频| 欧美性猛交╳xxx乱大交人| 最近视频中文字幕2019在线8| 亚洲无线观看免费| 亚洲国产精品成人综合色| 美女脱内裤让男人舔精品视频 | 成年版毛片免费区| 99久久精品国产国产毛片| 成人毛片a级毛片在线播放| 一级黄片播放器| 欧美日本亚洲视频在线播放| 国产精品伦人一区二区| 久久精品夜夜夜夜夜久久蜜豆| av黄色大香蕉| 欧美精品国产亚洲| 国产午夜精品一二区理论片| 超碰av人人做人人爽久久| 亚洲在久久综合| 一个人观看的视频www高清免费观看| 国产精品1区2区在线观看.| 免费观看的影片在线观看| 精品熟女少妇av免费看| 一进一出抽搐动态| 欧美日韩在线观看h| 欧美成人免费av一区二区三区| 男人舔女人下体高潮全视频| 老师上课跳d突然被开到最大视频| 亚洲精品色激情综合| 欧美最黄视频在线播放免费| 色吧在线观看| 日本在线视频免费播放| 美女国产视频在线观看| 国产色婷婷99| 国国产精品蜜臀av免费| 久久久久久久久久成人| 18禁在线无遮挡免费观看视频| 精品不卡国产一区二区三区| 国产精品久久久久久精品电影小说 | 联通29元200g的流量卡| 乱人视频在线观看| 老熟妇乱子伦视频在线观看| 最近手机中文字幕大全| 久久精品91蜜桃| 日韩 亚洲 欧美在线| 久久人人精品亚洲av| 蜜桃久久精品国产亚洲av| 在线免费观看的www视频| 成人综合一区亚洲| 嫩草影院精品99| 看免费成人av毛片| 18禁裸乳无遮挡免费网站照片| 我的老师免费观看完整版| 少妇熟女欧美另类| avwww免费| 久久久久久久亚洲中文字幕| 成人特级av手机在线观看| 精品少妇黑人巨大在线播放 | 国产单亲对白刺激| 亚洲综合色惰| 国产av麻豆久久久久久久| 欧美性猛交╳xxx乱大交人| 婷婷亚洲欧美| 男人和女人高潮做爰伦理| 久久欧美精品欧美久久欧美| 成人亚洲精品av一区二区| 国产精品国产三级国产av玫瑰| 亚洲人成网站在线观看播放| 一区二区三区四区激情视频 | 久久久久网色| av天堂在线播放| 国产精品久久电影中文字幕| av又黄又爽大尺度在线免费看 | 99久久精品一区二区三区| 综合色丁香网| 日本成人三级电影网站| 国产美女午夜福利| 亚洲国产精品sss在线观看| 亚洲图色成人| 国产国拍精品亚洲av在线观看| 中文字幕av在线有码专区| 欧美激情在线99| 国产成人精品久久久久久| 一级黄片播放器| avwww免费| 日韩大尺度精品在线看网址| 18禁在线播放成人免费| 成人亚洲欧美一区二区av| 久99久视频精品免费| 18禁黄网站禁片免费观看直播| 中文字幕人妻熟人妻熟丝袜美| 久久久久久国产a免费观看| 少妇熟女aⅴ在线视频| 国产亚洲精品久久久com| 日产精品乱码卡一卡2卡三| 九九爱精品视频在线观看| 99在线人妻在线中文字幕| 久久精品国产亚洲av涩爱 | 一区二区三区高清视频在线| 免费看日本二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲国产精品合色在线| 亚洲av免费高清在线观看| 久久国产乱子免费精品| 啦啦啦观看免费观看视频高清| 亚洲av第一区精品v没综合| 国产极品精品免费视频能看的| 久久精品国产99精品国产亚洲性色| 亚洲av免费在线观看| 亚洲在久久综合| 日韩亚洲欧美综合| 国产一区二区在线观看日韩| 国产欧美日韩精品一区二区| 天天躁日日操中文字幕| 此物有八面人人有两片| 一本精品99久久精品77| 免费观看人在逋| 麻豆久久精品国产亚洲av| 中文字幕制服av| 国产av不卡久久| 白带黄色成豆腐渣| 非洲黑人性xxxx精品又粗又长| 久久久久久伊人网av| 亚洲无线观看免费| 久久午夜亚洲精品久久| 国产一区二区在线观看日韩| 免费一级毛片在线播放高清视频| 丝袜美腿在线中文| 22中文网久久字幕| 国产av一区在线观看免费| 美女xxoo啪啪120秒动态图| 亚洲欧美精品自产自拍| 中文字幕制服av| 丰满乱子伦码专区| 小说图片视频综合网站| a级毛色黄片| 亚洲18禁久久av| 婷婷亚洲欧美| 在线免费观看不下载黄p国产| 91麻豆精品激情在线观看国产| 伦精品一区二区三区| 久久久久网色| 婷婷六月久久综合丁香| 亚洲无线观看免费| 日韩亚洲欧美综合| 久久亚洲精品不卡| 日日摸夜夜添夜夜添av毛片| 久久午夜亚洲精品久久| 男女下面进入的视频免费午夜| 日韩一区二区视频免费看| 男女视频在线观看网站免费| av在线观看视频网站免费| 成人美女网站在线观看视频| 欧美zozozo另类| 亚洲精品日韩在线中文字幕 | 青春草亚洲视频在线观看| 亚洲色图av天堂| 99久国产av精品国产电影| 国产极品天堂在线| 成熟少妇高潮喷水视频| 热99在线观看视频| 国产成人a区在线观看| 中国国产av一级| 三级经典国产精品| 精品一区二区免费观看| 日韩视频在线欧美|