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

    A study of response of thermocline in the South China Sea to ENSO events*

    2018-08-02 02:50:44PENGHanbang彭漢幫PANAijun潘愛(ài)軍ZHENGQuanan鄭全安HUJianyu胡建宇
    Journal of Oceanology and Limnology 2018年4期
    關(guān)鍵詞:建宇

    PENG Hanbang (彭漢幫) PAN Aijun (潘愛(ài)軍) ZHENG Quan’an (鄭全安) HU Jianyu (胡建宇)

    1 State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen 361102, China

    2 Ocean Dynamics Laboratory, The Third Institute of Oceanography, State Oceanic Administration (SOA), Xiamen 361005, China

    3 Department of Atmospheric and Oceanic Science, University of Maryland, College Park 20742, USA

    Abstract This paper investigates the response of the thermocline depth (TD) in the South China Sea(SCS) to the El Ni?o-Southern Oscillation (ENSO) events using 51-year (from 1960 to 2010) monthly seawater temperature and surface wind stress data acquired from the Simple Ocean Data Assimilation(SODA), together with heat flux data from the National Centers for Environmental Prediction (NCEP),precipitation data from the National Oceanic and Atmospheric Administration (NOAA) and evaporation data from the Woods Hole Oceanographic Institution (WHOI). It is indicated that the response of the SCS TD to the El Ni?o or La Ni?a events is in opposite phase. On one hand, the spatial-averaged TDs in the SCS(deeper than 200 m) appear as negative and positive anomalies during the mature phase of the El Ni?o and La Ni?a events, respectively. On the other hand, from June of the El Ni?o year to the subsequent April, the spatial patterns of TD in the north and south of 12°N appear as negative and positive anomalies, respectively,but present positive and negative anomalies for the La Ni?a case. However, positive and negative TD anomalies occur almost in the entire SCS in May of the subsequent year of the El Ni?o and La Ni?a events,respectively. It is suggested that the response of the TD in the SCS to the ENSO events is mainly caused by the sea surface buoyancy flux and the wind stress curl.

    Keyword: South China Sea (SCS); thermocline depth; El Ni?o-Southern Oscillation (ENSO); buoyancy flux; wind stress curl

    1 INTRODUCTION

    El Ni?o-Southern Oscillation (ENSO) over the tropical Paci fic Ocean is the strongest climate fluctuation on interannual and planetary scale with worldwide impacts. Generally, the warming phase is known as El Ni?o while the cooling phase as La Ni?a(Trenberth, 1997). The South China Sea (SCS),which is located in the East Asian monsoon region and close to the equatorial warm pool (Fig.1), is greatly in fluenced by the ENSO through teleconnections (i.e. “atmospheric bridge”) with in fluencing factors such as sea surface temperature(SST), sea surface winds, air temperature and humidity, as well as cloud cover and surface heat flux(Klein et al., 1999; Wang, 2002).

    Fig.1 Bottom topography of the South China Sea

    The in fluence of the ENSO events on the SCS (e.g.its hydrometeorology, temperature, heat flux and winds) has been investigated in previous studies. The response of the SST in the SCS to the ENSO events has been discussed by several investigators (Klein et al., 1999; Wang, 2002; 2006; Liu et al., 2014; Qiu et al., 2016). For example, Qiu et al. (2016) showed that the SST anomaly is positive and negative during the El Ni?o events in the southwestern and northeastern SCS, respectively, but conversely during the La Ni?a events. Kuo et al. (2004) attributed the response of the Vietnam coastal upwelling to the El Ni?o and La Ni?a events to the manifestation of the sea surface winds.Fang et al. (2006) indicated that the sea surface winds in the SCS are closely associated with the ENSO events at the interannual timescale. Yan et al. (2010)suggested that upper ocean heat content anomaly is usually positive while wind stress curl and latent heat flux are usually positive in the SCS during El Ni?o,and conversely during La Ni?a. However, response of the thermocline depth (TD) in the SCS to the ENSO events remains unclear. Changes in ocean heat content and strati fication are the key to understanding the dynamics of the SCS and the TD is an important indicative measure for them. According to previous studies, the TD in the SCS shows remarkable monthly and seasonal variability (Liu et al., 2001; Xiong,2012; Peng et al., 2017). Based on historical observation data, Hao et al. (2012) investigated the climatological seasonal variations of the thermocline in the China seas and northwestern Paci fic Ocean and regarded the sea surface buoyancy flux as the dominant driving mechanism for the TD variations.Using a coarse-resolution model of the Paci fic Ocean,Miller et al. (1997) showed that the modeled largescale thermocline anomaly is signi ficantly driven by large-scale Ekman pumping by interannual wind stress curl variations. Wang et al. (1999) suggested that the deepening and rising of the thermocline occur in phase with the local anticyclonic and cyclonic wind stress forcing in the tropical Paci fic, respectively.Using the Princeton Ocean Model (POM), Chu et al.(1999) argued that the seasonal variability of the TD in the SCS is affected by upwelling and downwelling.Based on the Simple Ocean Data Assimilation(SODA) data, Peng et al. (2017) examined the mechanisms causing monthly variation of the SCS TD, and the results reveal that the TD variations are in a good agreement with the sea surface buoyancy flux and wind stress curl. Based on the SCS Physical Oceanographic Dataset 2014 (SCSPOD14), Zeng et al. (2016) showed the monthly distribution of mixed layer depth (MLD), which is closely related to the TD, and found the MLD shows remarkable seasonal variations. Soon after, Zeng and Wang (2017)analyzed the seasonal variations of the barrier layer(BL) de fined by the MLD and ILD (isothermal layer depth, which is also related to the TD) in the SCS and the mechanism of BL variation, and the results reveal that the sea surface buoyancy flux and wind stress curl are the main physical mechanisms for the BL change. According to these results, we can conclude that the sea surface buoyancy flux and wind stress curl are the main mechanisms that can in fluence the TD in the SCS. As previously mentioned, the sea surface heat flux (the main component of buoyancy flux (Duan et al., 2012)) and the sea surface wind in the SCS can be in fluenced by the El Ni?o and La Ni?a events. Accordingly, a scienti fic question of interest is whether the El Ni?o or La Ni?a events can affect the TD in the SCS through those mechanisms. Therefore,in this study, we investigate the response of the TD in the SCS to the ENSO events and the possible physical processes.

    The rest of the paper is organized as follows.Section 2 gives a brief description of the data and methods as well as the El Ni?o and La Ni?a records.Section 3 analyses the TD variability in the SCS associated with the ENSO events and the corresponding mechanisms. Section 4 shows the relationship between thermocline intensity in the SCS and Ni?o 3.4 index (hereafter N34). Section 5 gives the discussion and conclusions.

    2 DATA AND METHOD

    2.1 Data

    The monthly seawater temperature, sea surface salinity and wind stress products used in this paper are from SODA (Carton et al., 2000a, b, 2005). The SODA model is based on observations including almost all available hydrographic pro file data, as well as ocean station data, mooring temperature and salinity time series, surface temperature and salinity observations of various types. The basic temperature and salinity observation sets consist of approximately 7×106pro files over the global ocean, of which two thirds have been obtained from the World Ocean Database 2009 (Giese and Ray, 2011), and the others include observations from the Tropical Atmosphere-Ocean/Triangle Trans-Ocean Buoy Network (TAO/TRITON) mooring thermistor array and Argo drifters.Output variables are averaged every five days, and then mapped onto a uniform global 0.5°×0.5°horizontal grid (a total of 720×330×40-level grid points), spanning the latitude range from 79.25°S to 89.25°N and the longitude range from 179.75°W to 179.75°E, using the horizontal grid spherical coordinate remapping and interpolation package(Carton and Giese, 2008). In this study, we use 50×50 pro files bounded by 0.25°–24.75°N and 100.25°–124.75°E and use the upper 25 vertical layers of the seawater temperature (from 5 to 1 625 m). The version of SODA used in our study is SODA v2.2.4 (Giese and Ray, 2011).

    The sea surface buoyancy flux is calculated by the sea surface net heat flux and net fresh water flux (Gill,1982; Schmitt et al., 1989). The monthly net long wave radiation, net short wave radiation, sensible heat flux and latent heat flux data (1960–2010,1.875°×1.875°) used for the sea surface net heat flux calculation are downloaded from the National Centers for Environmental Prediction (NCEP) reanalysis(Kalnay et al., 1996) (ftp://ftp.cdc.noaa.gov/Datasets/ncep.reanalysis.derived/surface_gauss/). We use the monthly precipitation (1960–2010, 2.5°×2.5°) and evaporation data (1960–2010, 1°×1°) for the sea surface net fresh flux calculation. The monthly precipitation data are downloaded from the National Oceanic and Atmospheric Administration (NOAA)’s Precipitation Reconstruction (PREC) Dataset (Chen et al., 2002) (http://www.esrl.noaa.gov/psd/data/gridded/data.prec.html). The oceanic precipitation analysis (PREC/O) is produced by EOF (Empirical Orthogonal Function) reconstruction of historical gauge observations over islands and land areas. The monthly evaporation data are downloaded from Objectively Analyzed air-sea Fluxes (OAFlux)project (Yu et al., 2008) of Woods Hole Oceanographic Institution (ftp://ftp.whoi.edu/pub/science/oa flux/data_v3/monthly/evaporation/). The OAFlux project uses the objective analysis to obtain optimal estimates of flux-related surface meteorology and then computes the global fluxes by using the state-of-the-art bulk flux parameterizations.

    2.2 Method

    In this study, we adopt the gradient criterion method to determine the thermocline. The gradient criterion of thermocline is 0.05°C/m (General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China and Standardization Administration of the People’s Republic of China, 2008), i.e., the thermocline is de fined as the layer, from down to top, with the vertical gradient of temperature larger than 0.05°C/m for at least 5 consecutive meters. We de fine the depth of upper boundary of the thermocline as the TD.

    The sea surface buoyancy flux, a comprehensive index of thermodynamics, re flects the coherent role of air-sea heat and haline exchanges (Gill, 1982; Schmitt et al., 1989). The buoyancy flux B (kg/(m·s3)) is calculated by

    where Bqis the thermal buoyancy due to the net heat flux, Bpis the haline buoyancy due to the net fresh water flux, g is the gravitational acceleration, Q is the net heat flux (upward positive; W/m2), ρ is the reference water density and Cpis the speci fic heat of water, α is the thermal coefficient of expansion and β is the coefficient of haline contraction, S0, E, and P are surface salinity, evaporation and precipitation,respectively. If surface water loses heat (i.e., Q >0) or evaporation is greater than precipitation (i.e., E– P >0),the surface water will lose buoyancy that can enhance vertically convection and deepen the TD. On the other way, if surface water gets heat or precipitation is greater than evaporation, the surface water will get buoyancy that strengthens strati fication stability of the upper ocean and decreases the TD.

    The monthly N34 is downloaded from the Climate Prediction Center of NOAA (http://www.cpc.noaa.gov/data/indices/ersst4.nino.mth.81-10.ascii). N34 is the averaged SST anomaly in the region bounded by 5°S to 5°N, from 170°W to 120°W and is one of the several ENSO indicators. Based on the criterion developed by the Climate Prediction Center of NOAA and de finition from Trenberth (1997), a 5-month running average of N34 exceeding +0.4°C or -0.4°C for at least 6 consecutive months is de fined as an El Ni?o or a La Ni?a event. According to this de finition,there are 9 El Ni?o events and 12 La Ni?a events from 1960 to 2010 (Table 1). In order to investigate the interannual variability of TD, the seasonal cycle must be removed. In this study, the monthly TD anomaly(TDA) is obtained by subtracting the monthly climatological mean value (from January of 1960 to December of 2010) from the TD.

    Table 1 The El Ni?o and La Ni?a events from 1960 to 2010

    Fig.2 A comparison of the monthly TDA time series in the SCS (black line) with the El Ni?o events (shaded in orange) and La Ni?a events (shaded in blue) de fined by N34 from 1960 to 2010

    3 THE TDA IN THE SCS ASSOCIATED WITH THE ENSO EVENTS

    3.1 The spatial-averaged TDA in the SCS during the ENSO events

    To obtain the spatial average of the SCS TDA, the TDA is averaged over the region with bathymetry deeper than 200 m in the SCS (Fig.1). Figure 2 shows the TDA and N34 time series in all calendar months from 1960 to 2010, the correlation coefficient between them is -0.46. One can see that the TDA varies in opposition with the El Ni?o and La Ni?a events, i.e.,the TDAs are negative during most time of the El Ni?o events, but positive during the La Ni?a events(except the 1984/1985 La Ni?a event).

    Figure 3 shows comparisons between the TDA time series and N34 during the event year ([0]) and the event’s subsequent year ([1]) for 9 El Ni?o events and 12 La Ni?a events. As shown in Fig.3a, the SCS TDAs almost show the same evolution during 9 El Ni?o events, i.e., the TDAs are negative from October[0] to March [1] with the peak values in February [1],while are positive from April [1] to June [1] with the peak values in May [1]. On the contrary, during the La Ni?a events (except the 1984/1985 event) (Fig.3b) the TDAs are positive from September [0] to March [1]with the peak values in February [1], while are negative from May [1] to July [1] with the peak values in May [1]. Furthermore, from Fig.3c we can be more intuitive to see the opposite sign of TDA during El Ni?o and La Ni?a events.

    It is noted that, as Fig.3c shows, the peak of the mean TDA in February [1] is lagged by two months with respect to the peak of the mean N34 in December[0] during the El Ni?o or La Ni?a event. This may result from the “atmosphere bridge” that connects the SCS with the SST anomaly in the central equatorial Paci fic, with some delay.

    Fig.3 Comparisons between the spatial-averaged TDA in the SCS and N34 during 9 El Ni?o events (a) and 12 La Ni?a events (b), and the mean N34 and mean TDA averaged over all El Ni?o events and La Ni?a events (c) from 1960 to 2010 in the event year ([0]) and the event’s subsequent year ([1])

    Fig.4 Zonal mean of composites of (a) TDA (units: m), (b)anomaly of B (units: 10 -5 kg/(m·s 3)), (c) anomaly of WSC (units: 10 -7 N/m 3) in the SCS from the El Ni?o year ([0]) to the subsequent year of the El Ni?o year([1])

    3.2 The zonal average and the possible mechanisms of the SCS TDA

    In Section 3.1 we have shown that the TDA in the SCS (spatially averaged TDA in the SCS) is in opposition with El Ni?o and La Ni?a SSTs (N34). In this section, we will further investigate the zonally averaged SCS TDA and the possible physical processes responsible for the TDA using composites of El Ni?o and La Ni?a events. Hereafter, we consider El Ni?o and La Ni?a composites computed from events over 1960 to 2010. According to the composite N34 (the thick red curves in Fig.3), we de fine the period from June [0] to September [0], from October[0] to February [1] and from March [1] to May [1] as the onset, mature and decay phases of the El Ni?o event, and the period from May [0] to August [0],from September [0] to February [1] and from March[1] to June [1] as the onset, mature and decay phases of the La Ni?a event.

    As mentioned in the introduction, the main mechanisms to in fluence TD in the SCS are the buoyancy flux ( B) and wind stress curl (WSC) on the sea surface. Lombardo and Gregg (1989) and Lozovatsky et al. (2005) suggested that B and wind stress on the sea surface are the main dynamic factors to in fluence the structure of the upper ocean. The loss of buoyancy in the upper ocean, which triggers convection, would deepen the thermocline. On the contrary, if the upper ocean gains buoyancy, the seawater stability would be strengthened inhibiting thermocline from deepening due to the more difficult vertical convection. In this study, positive or negative B means the upper ocean loses or gains buoyancy. In addition, the WSC also in fluences the TD by Ekman pumping. Positive WSC, which triggers upwelling,can make the thermocline shallower. Conversely,negative WSC can deepen the thermocline through downwelling. In short, positive B tends to increase the TD while positive WSC tends to decrease the TD.

    Figure 4 shows composites of zonally-averaged thermocline depth anomalies (TDAs), buoyancy flux( B) and wind stress curl (WSC) anomalies during the El Ni?o events in the SCS, as a function of latitude and time. Unlike the spatial average from Fig.3, the zonally-averaged TDA is not consistently negative during the El Ni?o events or positive during the La Ni?a events over the entire SCS. Instead, the evolutions of TDA in the northern and southern SCS (taking the 12°N as the boundary) are different. As shown in Fig.4a, the TDA is negative in the northern SCS but positive in the southern SCS from June [0] to April [1]during the El Ni?o event, especially from December[0] to March [1], when the TDA reaches -30 to -20 m north of 16°N. In May [1], the TDA is positive almost in the entire SCS (except in the region north of 18°N).

    In order to explain the above phenomenon, we also show the composites of zonally-averaged anomaly of B (Fig.4b) and WSC (Fig.4c) during the El Ni?o event. The TDA tends to be positive when anomaly of B is positive or anomaly of WSC is negative, and to be negative when anomaly of B is negative or anomaly of WSC is positive. As Fig.4b and c show, at the onset of the El Ni?o event from June [0] to September [0],the anomaly of B is negative in the northern SCS and positive in the southern SCS, while the anomaly of WSC is conversely positive in the northern SCS and negative in the southern SCS. From October [0] to March [1], the mature phase of the El Ni?o event, the anomaly of B is negative in the northern SCS but positive in the southern SCS. In particular, the larger negative TDA north of 16°N from December [0] to March [1] may be due to the larger negative anomaly of B and positive anomaly of WSC. Later in May [1],the decay phase of the El Ni?o event, the positive anomaly of B and negative anomaly of WSC are the primary causes for the positive TDA nearly in the entire SCS (except in the region north of 18°N).

    Fig.5 The same as Fig.4, but for the La Ni?a event

    Figure 5 shows the composites of zonally-averaged TDA, the anomaly of B and WSC in the SCS during the La Ni?a events. One can see that the zonallyaveraged TDA, anomaly of B and WSC from June [0]to May [1] during the La Ni?a events are almost opposite of those during the El Ni?o events as shown in Fig.4, except in the region north of 18°N from June[0] to October [0] for the TDA, and the region north of 12°N from August [0] to November [0] and the region south of 12°N from February [1] to May [1]for the anomaly of B.

    3.3 Spatial distribution of TDA and its mechanism in the SCS during the ENSO events

    The TDA in most of the SCS, either spatially averaged or zonally averaged, appears almost in opposition with the El Ni?o and La Ni?a SSTs and has different values at different stages (the onset,mature and decay phases) of the El Ni?o or La Ni?a events. We select August [0], February [1] and May[1] as the representative month of the onset, mature and decay phases of the El Ni?o or La Ni?a event to investigate the spatial structures of TDA, together with the anomaly of B and WSC in the SCS.

    Figure 6 shows longitude-latitude maps of composite TDA, anomaly of B and WSC at different stages (onset, mature and decay phases) of the El Ni?o events. In August [0] of the El Ni?o event, the TDA is negative in the northern SCS, but positive in the southern SCS (Fig.6a1). At the same time, the pattern of anomaly of B (Fig.6b1), which is negative in the northern SCS but positive in the southern SCS,is nearly the same with that of the TDA. Furthermore,the anomaly of WSC (Fig.6c1) is positive in most of the northern SCS, but appears negative in the southern SCS. In February [1] of the El Ni?o event, the TDA becomes more intense than that in August [0].Meanwhile, the positive TDA in the southern SCS expands northward (Fig.6a2). The anomaly of B, like the TDA, is more intense than that in August [0] and the positive values in the southern SCS expand northward (Fig.6b2), and the WSC anomaly in the southeastern and northwestern SCS is negative and positive (Fig.6c2), respectively. In May [1] of the El Ni?o event, the TDA is positive almost in the entire SCS except in the northeastern SCS (Fig.6a3).Correspondingly, the anomaly of B is positive almost in the entire SCS (Fig.6b3). Besides, the anomaly of WSC is negative almost in the SCS except the positive anomaly on the continental shelf break (Fig.6c3).

    Figure 7 shows the spatial distributions of composite TDA, anomaly of B and WSC at different stages (onset, mature and decay phases) of the La Ni?a events. During the La Ni?a event, the TDA,anomaly of B and WSC nearly appear as opposite modes with those during the El Ni?o event, except in the region south of 12°N for the anomaly of B.

    4 THE RELATIONSHIP BETWEEN THERMOCLINE INTENSITY (TI) IN THE SCS AND N34

    As one of the indicative parameters of thermocline,the TD in the SCS presents pronounced variation during the ENSO events. Other indicative parameters of thermocline in the SCS, such as the thermocline thickness and intensity could also be in fluenced by the ENSO events. We analyze here the evolution of the thermocline intensity (TI) in the SCS and find that it shows a good correlation with N34.

    Fig.6 Composites of TDA (units: m) (a1–a3), anomaly of B (units: 10 -5 kg/(m·s 3)) (b1–b3), anomaly of WSC (units: 10 -7 N/m 3)(c1–c3) in the SCS in August [0], February [1] and May [1] of the El Ni?o event

    The TI is de fined as the average intensity over the thermocline, namely the ratio of the temperature difference and the depth difference between the upper and lower boundaries of the thermocline.

    Figure 8 shows the lag correlation coefficient between the TI anomaly (TIA) and N34, and the time series of spatially-averaged TIA in the SCS and N34 during the same period. The TIA averaged over the SCS shows a lag correlation with N34 and the highest correlation coefficient is 0.63, when the TIA lags N34 for 5 months (Fig.8a). Figure 8b shows the time series of the spatially averaged TIA and N34 at the same period and the correlation coefficient is 0.41 (see Fig.8a). From the spatial distribution of the correlation coefficients between TIA and N34 in the SCS (Fig.9),one can see that the correlation coefficients increase gradually and the larger ones become wider from -5 to 5 months that TIA lags N34. It is noted that the correlation coefficient in the southeastern SCS (0.30–0.50) is larger than that in the northwestern SCS (0–0.20). The reason for such a feature deserves further study in the near future.

    Fig.7 The same as Fig.6, but for the La Ni?a event

    5 CONCLUSION AND DISCUSSION

    In this study, the changes in thermocline depth(TD) in the South China Sea (SCS) associated with the El Ni?o-Southern Oscillation (ENSO) events during 1960 to 2010 and the mechanisms have been investigated. The main conclusions are as follows.

    Fig.8 Lag correlation between spatially averaged TIA in the SCS and N34 (a), a positive lag means N34 is leading TIA; time series of spatially averaged TIA in the SCS and N34 (shading) during the same period (b)

    Fig.9 Correlation between N34 and the SCS TIA at a lag of -5 months (a), 0 months (b) and 5 months (c)

    (1) There are 9 El Ni?o events and 12 La Ni?a events that happened from 1960 to 2010 according to the criterion proposed by the Climate Prediction Center of NOAA based on N34. The TD anomalies in the SCS are in opposite with El Ni?o and La Ni?a SSTs for all those events (except the 1984/1985 La Ni?a event).

    (2) By compositing the SCS TDA during 9 El Ni?o events, we find that the TDA averaged over the SCS is negative from the onset to mature phase of the El Ni?o event, but positive at the decay phase. The spatial patterns of the TDA are negative in the northern SCS and positive in the southern SCS during the onset phase, more intense during the mature phase, but become positive over the entire SCS in the decay phase.

    (3) By compositing the SCS TDA during 12 La Ni?a events (except the 1984/1985 La Ni?a event),we find that the TDA averaged over the SCS is positive from the onset to mature phase of the La Ni?a event, but negative in the decay phase. The TDA, for the spatial patterns, is positive and negative during the onset phase in the northern and southern SCS,respectively, more intense during the mature phase,but becomes negative over the entire SCS in the decay phase.

    (4) The main reason for the opposite sign of the TDA with respect to El Ni?o and La Ni?a SSTs is that the SST anomaly on the Tropical Paci fic Ocean,denoted by N34, can in fluence the buoyancy flux ( B)and the wind stress curl (WSC) on the sea surface of the SCS, which are the main mechanisms to affect the TD in the SCS, through the “atmospheric bridge”.

    The TDA, anomaly of B and WSC show the reverse mode between the northern and southern SCS during the mature phase of El Ni?o and La Ni?a events. That may be due to the following aspects. (1) Qu et al.(2004) argued that the Luzon Strait transport can play a role in conveying the import of ENSO to the SCS.This may in fluence the haline exchanges or the heat content, which can in fluence B, and then the TD in northern SCS. (2) Due to the in fluenced wind stress over the SCS during ENSO and the different topographic effect around SCS, the anomaly of WSC may also present different over the northern and southern SCS and hence the TD.

    Both B and WSC can in fluence TD in the SCS, but B is dominant (Hao et al., 2012; Peng et al., 2017).For the spatially average, TD matches B better than WSC during each individual ENSO event, except for the 1971/1972 La Ni?a event and the 2004/2005 and 1991/1992 El Ni?o events, of which both B and WSC are in a good agreement with TD. For the spatial distribution, the TDA is very similar to the composite mode during most of the ENSO events. The TDA north of SCS, however, is more intense and expands southward during the strong ENSO events (1997/1998 El Ni?o event and 1973/1974 La Ni?a event), anomaly of B and WSC are the same (not shown).

    猜你喜歡
    建宇
    Effects of confining pressure and pore pressure on multipole borehole acoustic field in fluid-saturated porous media
    Square grid pattern with direction-selective surface discharges in dielectric barrier discharge
    Formation of honeycomb-Kagome hexagonal superlattice pattern with dark discharges in dielectric barrier discharge
    Analysis of monthly variability of thermocline in the South China Sea*
    含 虛(一石二座)
    寶藏(2018年4期)2018-05-07 01:58:24
    酷蟲(chóng)學(xué)校
    糾紛的根源
    多悉善感的鹿角蟲(chóng)
    必勝嗎,狼蛛班長(zhǎng)?
    跳高比賽中的意外
    中文字幕高清在线视频| 成年人黄色毛片网站| 亚洲熟妇熟女久久| 国产av一区二区精品久久| 搡老岳熟女国产| 中文字幕制服av| 亚洲成人免费电影在线观看| 国产精品乱码一区二三区的特点 | 午夜影院日韩av| 天堂俺去俺来也www色官网| 在线观看免费视频日本深夜| 99国产精品一区二区三区| 久久青草综合色| 午夜成年电影在线免费观看| 桃红色精品国产亚洲av| 日韩欧美免费精品| 久久精品人人爽人人爽视色| 91精品三级在线观看| 另类亚洲欧美激情| 亚洲伊人色综图| 国产99白浆流出| 激情在线观看视频在线高清 | 他把我摸到了高潮在线观看| 午夜福利欧美成人| 日韩欧美国产一区二区入口| 建设人人有责人人尽责人人享有的| 下体分泌物呈黄色| 国产精品电影一区二区三区 | av免费在线观看网站| 国产免费男女视频| 成人永久免费在线观看视频| 精品高清国产在线一区| 免费观看精品视频网站| 日日爽夜夜爽网站| 看免费av毛片| 久久久久久久久免费视频了| 精品久久久久久久久久免费视频 | 午夜亚洲福利在线播放| 欧美日韩中文字幕国产精品一区二区三区 | 国产精品一区二区在线不卡| 99久久综合精品五月天人人| 成年女人毛片免费观看观看9 | 人妻丰满熟妇av一区二区三区 | 一本综合久久免费| 中文字幕色久视频| 国产成人av教育| 很黄的视频免费| 国产97色在线日韩免费| 国产精品秋霞免费鲁丝片| 国产精品久久久久久人妻精品电影| 国产精品国产av在线观看| 亚洲成人免费电影在线观看| 日韩中文字幕欧美一区二区| 老熟妇乱子伦视频在线观看| 欧美激情久久久久久爽电影 | 国产有黄有色有爽视频| 亚洲精品中文字幕一二三四区| 性色av乱码一区二区三区2| 夜夜爽天天搞| 国产亚洲精品久久久久5区| 久久人妻熟女aⅴ| 天天操日日干夜夜撸| 最近最新中文字幕大全免费视频| 18禁裸乳无遮挡免费网站照片 | 国产亚洲av高清不卡| 妹子高潮喷水视频| videos熟女内射| 十八禁高潮呻吟视频| 国产真人三级小视频在线观看| 天天影视国产精品| 欧美日韩亚洲高清精品| 午夜免费成人在线视频| 激情在线观看视频在线高清 | а√天堂www在线а√下载 | 少妇被粗大的猛进出69影院| 亚洲国产欧美日韩在线播放| 十八禁高潮呻吟视频| 91成人精品电影| 亚洲国产看品久久| 人人妻人人爽人人添夜夜欢视频| 亚洲,欧美精品.| 亚洲中文av在线| 亚洲视频免费观看视频| 黄片小视频在线播放| 少妇被粗大的猛进出69影院| 久久香蕉国产精品| 久久久久久久午夜电影 | 美女福利国产在线| 国产99白浆流出| 又黄又粗又硬又大视频| 国产亚洲精品久久久久久毛片 | 国产成人精品久久二区二区免费| 法律面前人人平等表现在哪些方面| 最近最新中文字幕大全电影3 | 99久久综合精品五月天人人| 久久久国产一区二区| 老司机福利观看| 伦理电影免费视频| 日日爽夜夜爽网站| 国产亚洲精品一区二区www | 国产亚洲欧美在线一区二区| 精品国产一区二区三区久久久樱花| 国内毛片毛片毛片毛片毛片| 男女之事视频高清在线观看| 亚洲aⅴ乱码一区二区在线播放 | 日本黄色视频三级网站网址 | 夜夜爽天天搞| av网站免费在线观看视频| 叶爱在线成人免费视频播放| 国产一区二区三区综合在线观看| 五月开心婷婷网| 一个人免费在线观看的高清视频| 久久久国产精品麻豆| 国产亚洲精品久久久久5区| 交换朋友夫妻互换小说| 亚洲一卡2卡3卡4卡5卡精品中文| 女人精品久久久久毛片| 十分钟在线观看高清视频www| 国内久久婷婷六月综合欲色啪| 亚洲在线自拍视频| 一级毛片高清免费大全| 又黄又爽又免费观看的视频| 手机成人av网站| 国产亚洲精品久久久久久毛片 | 少妇裸体淫交视频免费看高清 | 精品无人区乱码1区二区| 99香蕉大伊视频| 欧美黄色淫秽网站| 国产精品欧美亚洲77777| 日日爽夜夜爽网站| www日本在线高清视频| 国产精品1区2区在线观看. | 久久精品国产清高在天天线| 黄片小视频在线播放| 女人久久www免费人成看片| 天堂中文最新版在线下载| 母亲3免费完整高清在线观看| 久久国产精品大桥未久av| 老司机影院毛片| 夜夜爽天天搞| 男人操女人黄网站| 三级毛片av免费| 国产精品免费大片| 亚洲欧洲精品一区二区精品久久久| 久久人妻av系列| 一级毛片精品| 99精品在免费线老司机午夜| 久久久久久久久免费视频了| 精品久久久精品久久久| 中文亚洲av片在线观看爽 | 18禁裸乳无遮挡免费网站照片 | 国产精品九九99| 久久精品国产亚洲av高清一级| 欧美日韩成人在线一区二区| 久久这里只有精品19| 成熟少妇高潮喷水视频| 男女高潮啪啪啪动态图| 91九色精品人成在线观看| 午夜久久久在线观看| 搡老乐熟女国产| netflix在线观看网站| 国产精品1区2区在线观看. | 国产精品 国内视频| 狂野欧美激情性xxxx| 精品国产国语对白av| 久久久久久人人人人人| 在线免费观看的www视频| 亚洲精华国产精华精| 成人三级做爰电影| 欧美乱色亚洲激情| 丝瓜视频免费看黄片| 欧美黑人欧美精品刺激| 午夜福利在线观看吧| svipshipincom国产片| 97人妻天天添夜夜摸| 两个人免费观看高清视频| 1024视频免费在线观看| 丰满迷人的少妇在线观看| 国产精品影院久久| 电影成人av| 久久青草综合色| 亚洲国产精品sss在线观看 | 在线观看免费午夜福利视频| 欧美人与性动交α欧美软件| av欧美777| 国产熟女午夜一区二区三区| 国产在视频线精品| 一a级毛片在线观看| 欧美色视频一区免费| 黑人操中国人逼视频| 日韩精品免费视频一区二区三区| 69精品国产乱码久久久| 在线十欧美十亚洲十日本专区| 亚洲视频免费观看视频| 很黄的视频免费| 一二三四在线观看免费中文在| 国产精品二区激情视频| 操美女的视频在线观看| 手机成人av网站| 亚洲熟妇熟女久久| 少妇的丰满在线观看| 在线观看免费视频日本深夜| 黑丝袜美女国产一区| 夫妻午夜视频| 国产深夜福利视频在线观看| 午夜福利在线观看吧| 丝袜美足系列| 在线观看免费视频网站a站| 后天国语完整版免费观看| 老司机深夜福利视频在线观看| 国产91精品成人一区二区三区| 老司机午夜十八禁免费视频| 久久国产精品影院| 亚洲av日韩精品久久久久久密| 多毛熟女@视频| 亚洲av成人不卡在线观看播放网| 老鸭窝网址在线观看| 啦啦啦视频在线资源免费观看| 国产淫语在线视频| 在线观看日韩欧美| 手机成人av网站| av有码第一页| 伦理电影免费视频| 精品亚洲成国产av| 国产精品偷伦视频观看了| 一本大道久久a久久精品| 99re6热这里在线精品视频| 精品欧美一区二区三区在线| 久久久精品区二区三区| 狠狠狠狠99中文字幕| 人人妻人人添人人爽欧美一区卜| 波多野结衣av一区二区av| 久久午夜亚洲精品久久| 成年人午夜在线观看视频| 亚洲精品美女久久av网站| 欧美成人免费av一区二区三区 | 超色免费av| 午夜91福利影院| 日韩欧美三级三区| 777米奇影视久久| 悠悠久久av| 人成视频在线观看免费观看| 看片在线看免费视频| 精品人妻在线不人妻| 国产精品秋霞免费鲁丝片| 一边摸一边做爽爽视频免费| 国产有黄有色有爽视频| 人妻丰满熟妇av一区二区三区 | 18禁裸乳无遮挡动漫免费视频| 欧美日韩黄片免| 精品福利永久在线观看| 9191精品国产免费久久| 亚洲第一青青草原| 精品欧美一区二区三区在线| 亚洲全国av大片| 中亚洲国语对白在线视频| 搡老熟女国产l中国老女人| 日本精品一区二区三区蜜桃| 少妇的丰满在线观看| 纯流量卡能插随身wifi吗| 亚洲一区二区三区欧美精品| 免费在线观看视频国产中文字幕亚洲| 国产一区二区激情短视频| 久久ye,这里只有精品| 村上凉子中文字幕在线| 69精品国产乱码久久久| 国产亚洲av高清不卡| 久久人妻熟女aⅴ| 在线观看免费高清a一片| 国产精品国产av在线观看| 少妇的丰满在线观看| 999久久久国产精品视频| 国产精品久久视频播放| 不卡av一区二区三区| 黄色视频不卡| 免费av中文字幕在线| 精品亚洲成国产av| 女警被强在线播放| 亚洲三区欧美一区| 男女下面插进去视频免费观看| 亚洲av片天天在线观看| 捣出白浆h1v1| 搡老乐熟女国产| 国产淫语在线视频| 在线av久久热| 嫩草影视91久久| 亚洲五月婷婷丁香| 国产免费男女视频| 欧美黑人欧美精品刺激| 夫妻午夜视频| 巨乳人妻的诱惑在线观看| 日本一区二区免费在线视频| 国产一区二区三区视频了| 久久天躁狠狠躁夜夜2o2o| 久久精品亚洲精品国产色婷小说| 黄色视频不卡| 精品久久久久久久久久免费视频 | 美女高潮到喷水免费观看| 黄色丝袜av网址大全| 美女午夜性视频免费| 亚洲欧美激情综合另类| 丝袜在线中文字幕| 最新的欧美精品一区二区| 亚洲午夜精品一区,二区,三区| 国产精品av久久久久免费| 午夜免费鲁丝| 热99久久久久精品小说推荐| 中文字幕色久视频| 女人高潮潮喷娇喘18禁视频| 一级作爱视频免费观看| 精品国产亚洲在线| 久久久久国内视频| a级毛片在线看网站| 看黄色毛片网站| 女人高潮潮喷娇喘18禁视频| 一区在线观看完整版| 日韩中文字幕欧美一区二区| 亚洲一区二区三区欧美精品| 午夜激情av网站| 午夜亚洲福利在线播放| 女人爽到高潮嗷嗷叫在线视频| 日韩熟女老妇一区二区性免费视频| 人妻丰满熟妇av一区二区三区 | 国产精品免费大片| 女人高潮潮喷娇喘18禁视频| 无人区码免费观看不卡| 超碰97精品在线观看| 桃红色精品国产亚洲av| 精品少妇一区二区三区视频日本电影| 精品福利永久在线观看| 国产精品乱码一区二三区的特点 | 久久久国产精品麻豆| 国产激情欧美一区二区| 捣出白浆h1v1| 亚洲精品久久午夜乱码| 天堂动漫精品| 50天的宝宝边吃奶边哭怎么回事| 午夜福利,免费看| 91大片在线观看| 在线永久观看黄色视频| 国产精品久久电影中文字幕 | 夜夜夜夜夜久久久久| 乱人伦中国视频| 热99re8久久精品国产| 日韩免费av在线播放| 一级作爱视频免费观看| 欧美 日韩 精品 国产| 青草久久国产| 精品国产超薄肉色丝袜足j| 欧美亚洲 丝袜 人妻 在线| 变态另类成人亚洲欧美熟女 | 日本a在线网址| 高清欧美精品videossex| 免费高清在线观看日韩| 精品国产国语对白av| 免费高清在线观看日韩| 亚洲午夜理论影院| 天天添夜夜摸| 99久久99久久久精品蜜桃| 一个人免费在线观看的高清视频| 热re99久久国产66热| 成人亚洲精品一区在线观看| 一级毛片精品| 亚洲av日韩精品久久久久久密| 国产成人免费无遮挡视频| 老汉色av国产亚洲站长工具| 亚洲色图综合在线观看| 亚洲中文av在线| 一级毛片高清免费大全| 国产成人精品在线电影| 日日爽夜夜爽网站| videosex国产| 亚洲av成人一区二区三| 精品久久久久久久毛片微露脸| 亚洲av片天天在线观看| 99精国产麻豆久久婷婷| 9191精品国产免费久久| 国产精品av久久久久免费| 亚洲一区二区三区不卡视频| 俄罗斯特黄特色一大片| 中文字幕精品免费在线观看视频| а√天堂www在线а√下载 | 别揉我奶头~嗯~啊~动态视频| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩一级在线毛片| 黑人巨大精品欧美一区二区蜜桃| 99久久精品国产亚洲精品| 国产野战对白在线观看| 欧美丝袜亚洲另类 | 免费在线观看影片大全网站| 在线观看日韩欧美| 黑人巨大精品欧美一区二区蜜桃| 99久久精品国产亚洲精品| 亚洲美女黄片视频| 色综合欧美亚洲国产小说| 欧美在线黄色| 国产精品欧美亚洲77777| 亚洲视频免费观看视频| 精品一区二区三区av网在线观看| 啦啦啦 在线观看视频| 日日爽夜夜爽网站| 国产一区二区三区在线臀色熟女 | 一级a爱片免费观看的视频| av网站在线播放免费| 久热这里只有精品99| 国产精品一区二区在线不卡| www.999成人在线观看| 免费观看精品视频网站| 色94色欧美一区二区| 精品久久久精品久久久| 满18在线观看网站| 欧美日韩精品网址| 嫁个100分男人电影在线观看| 久久香蕉精品热| 老熟女久久久| 久久精品国产清高在天天线| 淫妇啪啪啪对白视频| 国产欧美日韩精品亚洲av| 少妇猛男粗大的猛烈进出视频| 天天添夜夜摸| 久久久久国产精品人妻aⅴ院 | 国产亚洲一区二区精品| 麻豆国产av国片精品| 欧美+亚洲+日韩+国产| 欧美丝袜亚洲另类 | 欧美av亚洲av综合av国产av| 国产精品久久电影中文字幕 | 国产精品免费大片| 久久中文看片网| 精品国内亚洲2022精品成人 | 脱女人内裤的视频| 黑丝袜美女国产一区| 超碰成人久久| 久久精品成人免费网站| 丰满饥渴人妻一区二区三| 制服诱惑二区| 999久久久国产精品视频| 两人在一起打扑克的视频| 亚洲国产精品sss在线观看 | 男男h啪啪无遮挡| 国产亚洲欧美在线一区二区| 一本大道久久a久久精品| 夜夜爽天天搞| 欧美日韩黄片免| 黄色片一级片一级黄色片| 女人高潮潮喷娇喘18禁视频| 熟女少妇亚洲综合色aaa.| 久久国产精品男人的天堂亚洲| 男女高潮啪啪啪动态图| 国产一区二区三区综合在线观看| 麻豆成人av在线观看| 大香蕉久久成人网| 韩国精品一区二区三区| 欧美一级毛片孕妇| 99热只有精品国产| 黑人操中国人逼视频| 99国产精品一区二区蜜桃av | 黄色片一级片一级黄色片| 日韩成人在线观看一区二区三区| 久久久久精品国产欧美久久久| 午夜福利,免费看| 一级毛片高清免费大全| 男女免费视频国产| 欧美精品亚洲一区二区| 亚洲色图 男人天堂 中文字幕| 91字幕亚洲| 亚洲成国产人片在线观看| 久久精品国产亚洲av香蕉五月 | 99国产综合亚洲精品| 老司机亚洲免费影院| 欧美亚洲 丝袜 人妻 在线| 天天躁夜夜躁狠狠躁躁| 精品久久蜜臀av无| 一进一出好大好爽视频| av有码第一页| 少妇猛男粗大的猛烈进出视频| 国产欧美日韩精品亚洲av| 久久久久久久久免费视频了| 精品久久久久久久毛片微露脸| 欧美激情久久久久久爽电影 | 亚洲五月天丁香| 精品一区二区三区视频在线观看免费 | 日本黄色日本黄色录像| 亚洲欧美一区二区三区久久| 色精品久久人妻99蜜桃| 免费在线观看影片大全网站| 老司机午夜福利在线观看视频| 国产不卡av网站在线观看| 窝窝影院91人妻| 女性生殖器流出的白浆| 手机成人av网站| 成人黄色视频免费在线看| 乱人伦中国视频| 欧洲精品卡2卡3卡4卡5卡区| 国产欧美日韩精品亚洲av| 日韩免费高清中文字幕av| 一边摸一边抽搐一进一小说 | 国产亚洲欧美精品永久| 精品乱码久久久久久99久播| 在线观看舔阴道视频| 啪啪无遮挡十八禁网站| 久久婷婷成人综合色麻豆| 成人国语在线视频| 日韩欧美三级三区| 午夜福利影视在线免费观看| 男人的好看免费观看在线视频 | 成在线人永久免费视频| 97人妻天天添夜夜摸| 99香蕉大伊视频| 国产又爽黄色视频| 精品少妇久久久久久888优播| 99国产综合亚洲精品| 亚洲五月色婷婷综合| 国产成人一区二区三区免费视频网站| 久久亚洲精品不卡| 真人做人爱边吃奶动态| 日韩人妻精品一区2区三区| 纯流量卡能插随身wifi吗| 亚洲综合色网址| 精品人妻熟女毛片av久久网站| 国产免费男女视频| 制服诱惑二区| 巨乳人妻的诱惑在线观看| 国产三级黄色录像| 在线观看舔阴道视频| 在线观看午夜福利视频| 婷婷精品国产亚洲av在线 | 久久久久国产精品人妻aⅴ院 | 色在线成人网| 黄片播放在线免费| 久久中文字幕一级| 精品国产亚洲在线| 老司机午夜福利在线观看视频| 亚洲男人天堂网一区| 咕卡用的链子| 青草久久国产| 免费少妇av软件| 久久青草综合色| 亚洲av熟女| 婷婷成人精品国产| 777米奇影视久久| 成年人免费黄色播放视频| 国产在视频线精品| 制服诱惑二区| 亚洲国产看品久久| 真人做人爱边吃奶动态| 亚洲一区二区三区欧美精品| 亚洲少妇的诱惑av| 国产精品一区二区免费欧美| 51午夜福利影视在线观看| 老司机靠b影院| 中出人妻视频一区二区| 久久久水蜜桃国产精品网| 曰老女人黄片| 日本wwww免费看| 精品一品国产午夜福利视频| 丁香六月欧美| 中国美女看黄片| 男女午夜视频在线观看| 淫妇啪啪啪对白视频| 在线看a的网站| 999久久久国产精品视频| 免费av中文字幕在线| 天天添夜夜摸| 动漫黄色视频在线观看| 亚洲综合色网址| 五月开心婷婷网| 侵犯人妻中文字幕一二三四区| 国产麻豆69| 午夜福利,免费看| 色婷婷久久久亚洲欧美| 99热国产这里只有精品6| 国产免费av片在线观看野外av| 成人精品一区二区免费| 亚洲人成电影观看| 色在线成人网| 精品国内亚洲2022精品成人 | 欧美 亚洲 国产 日韩一| 久久午夜亚洲精品久久| avwww免费| 91精品国产国语对白视频| 一边摸一边抽搐一进一小说 | 9色porny在线观看| 少妇裸体淫交视频免费看高清 | 免费少妇av软件| 少妇猛男粗大的猛烈进出视频| 国产黄色免费在线视频| 波多野结衣av一区二区av| 国产激情欧美一区二区| 麻豆av在线久日| 欧美av亚洲av综合av国产av| 日日爽夜夜爽网站| 日韩大码丰满熟妇| 巨乳人妻的诱惑在线观看| 精品高清国产在线一区| 大香蕉久久成人网| 精品视频人人做人人爽| 91在线观看av| 久9热在线精品视频| 国产精品久久久人人做人人爽| 久久精品亚洲精品国产色婷小说| 一级作爱视频免费观看| 欧美午夜高清在线| 一本大道久久a久久精品| 在线观看免费午夜福利视频| 80岁老熟妇乱子伦牲交| 成人永久免费在线观看视频| 欧美日韩福利视频一区二区| 国产一卡二卡三卡精品| 色婷婷久久久亚洲欧美| 很黄的视频免费| 一二三四社区在线视频社区8| 亚洲欧美精品综合一区二区三区| 国产精品亚洲av一区麻豆| 岛国毛片在线播放|