• <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)?
    跳高比賽中的意外
    成人欧美大片| 不卡一级毛片| 只有这里有精品99| 少妇的逼好多水| 一边摸一边抽搐一进一小说| 免费看日本二区| 看黄色毛片网站| av在线蜜桃| 又粗又爽又猛毛片免费看| 91av网一区二区| 久久精品夜夜夜夜夜久久蜜豆| 噜噜噜噜噜久久久久久91| 尤物成人国产欧美一区二区三区| 日本一二三区视频观看| 最近2019中文字幕mv第一页| 中文字幕久久专区| 高清午夜精品一区二区三区 | 国产精品爽爽va在线观看网站| 高清毛片免费观看视频网站| 精品久久久久久久久久久久久| 欧美日韩国产亚洲二区| 久久99精品国语久久久| 亚洲精品粉嫩美女一区| 日本一二三区视频观看| 亚洲人成网站高清观看| 小说图片视频综合网站| 欧美bdsm另类| 欧美在线一区亚洲| 中文字幕制服av| 97超碰精品成人国产| 亚洲国产高清在线一区二区三| 搡女人真爽免费视频火全软件| 97超视频在线观看视频| 日本av手机在线免费观看| 99热只有精品国产| 日韩高清综合在线| 久久午夜亚洲精品久久| 99久久无色码亚洲精品果冻| av在线播放精品| 日日摸夜夜添夜夜爱| 午夜福利视频1000在线观看| 国产精品美女特级片免费视频播放器| 成人鲁丝片一二三区免费| 别揉我奶头 嗯啊视频| 成人一区二区视频在线观看| 色综合色国产| 国产成人精品久久久久久| 最好的美女福利视频网| www日本黄色视频网| 国产蜜桃级精品一区二区三区| 男女做爰动态图高潮gif福利片| 国产成人福利小说| 亚洲七黄色美女视频| 99久久精品国产国产毛片| 国产中年淑女户外野战色| 一边摸一边抽搐一进一小说| 麻豆av噜噜一区二区三区| 蜜桃亚洲精品一区二区三区| 看黄色毛片网站| 此物有八面人人有两片| 淫秽高清视频在线观看| 国产69精品久久久久777片| 国产成人午夜福利电影在线观看| 国语自产精品视频在线第100页| 哪个播放器可以免费观看大片| 一级毛片久久久久久久久女| 人妻久久中文字幕网| 国产成人freesex在线| 亚洲欧美精品专区久久| 最近中文字幕高清免费大全6| 嘟嘟电影网在线观看| 久久久久久久久中文| 久久精品人妻少妇| 欧美一区二区国产精品久久精品| 一级黄片播放器| 日韩欧美在线乱码| 人人妻人人澡欧美一区二区| 欧美性感艳星| 少妇人妻精品综合一区二区 | 中文字幕av成人在线电影| 国产男人的电影天堂91| 国产淫片久久久久久久久| 精品国内亚洲2022精品成人| 你懂的网址亚洲精品在线观看 | 天美传媒精品一区二区| 性欧美人与动物交配| 精品少妇黑人巨大在线播放 | 国产视频首页在线观看| 国产69精品久久久久777片| av在线亚洲专区| 欧美日本亚洲视频在线播放| 最近中文字幕高清免费大全6| 中文字幕免费在线视频6| 看十八女毛片水多多多| 亚洲一级一片aⅴ在线观看| 成人无遮挡网站| 天堂影院成人在线观看| 不卡一级毛片| 综合色av麻豆| 久久精品国产99精品国产亚洲性色| 内射极品少妇av片p| 日韩欧美精品v在线| 一个人免费在线观看电影| 免费看日本二区| 国产精品.久久久| 一区福利在线观看| 99久久无色码亚洲精品果冻| 男女那种视频在线观看| 成人毛片60女人毛片免费| 深夜精品福利| 欧美xxxx性猛交bbbb| 波多野结衣巨乳人妻| 亚洲人成网站在线播放欧美日韩| 久久婷婷人人爽人人干人人爱| 99久久无色码亚洲精品果冻| 国产大屁股一区二区在线视频| 欧美三级亚洲精品| 国产私拍福利视频在线观看| 成人性生交大片免费视频hd| 久久久欧美国产精品| 免费看a级黄色片| 国产黄色小视频在线观看| 欧美3d第一页| 精品少妇黑人巨大在线播放 | 少妇猛男粗大的猛烈进出视频 | www日本黄色视频网| 亚洲五月天丁香| 亚洲天堂国产精品一区在线| 国产高潮美女av| 国产一区亚洲一区在线观看| 嫩草影院精品99| 毛片女人毛片| 午夜免费男女啪啪视频观看| 国产精品久久久久久亚洲av鲁大| 免费av毛片视频| 亚洲av熟女| 国产精品久久久久久久电影| 青春草国产在线视频 | 欧美3d第一页| 人妻久久中文字幕网| 成年av动漫网址| 亚洲第一电影网av| 久久久欧美国产精品| 伦精品一区二区三区| 久久久久久久久久久免费av| 插逼视频在线观看| 一个人看的www免费观看视频| 黄色一级大片看看| av.在线天堂| 亚洲av男天堂| 校园人妻丝袜中文字幕| 亚洲自偷自拍三级| 长腿黑丝高跟| 亚洲国产精品成人综合色| 51国产日韩欧美| 精品久久久噜噜| 老熟妇乱子伦视频在线观看| 亚洲欧美精品专区久久| 偷拍熟女少妇极品色| 免费在线观看成人毛片| 国产日本99.免费观看| 亚洲av不卡在线观看| 精品久久久噜噜| 成人特级黄色片久久久久久久| av在线观看视频网站免费| 亚洲国产欧美在线一区| 69人妻影院| 99热全是精品| 可以在线观看毛片的网站| .国产精品久久| 国产色婷婷99| 国产精品综合久久久久久久免费| 中文字幕熟女人妻在线| 99久国产av精品| 女同久久另类99精品国产91| 成人特级av手机在线观看| 国产一级毛片在线| 色哟哟哟哟哟哟| 中国国产av一级| 中出人妻视频一区二区| 小说图片视频综合网站| 中国国产av一级| 直男gayav资源| 欧美三级亚洲精品| 成人漫画全彩无遮挡| 午夜精品国产一区二区电影 | 少妇的逼好多水| 亚洲精品456在线播放app| 日本成人三级电影网站| 舔av片在线| 婷婷色av中文字幕| 一级黄色大片毛片| 久久久久久久午夜电影| 麻豆乱淫一区二区| 亚洲国产日韩欧美精品在线观看| 日韩欧美一区二区三区在线观看| 精品人妻熟女av久视频| 久久精品国产亚洲av涩爱 | 97超碰精品成人国产| 久久99热这里只有精品18| 在线免费十八禁| 搞女人的毛片| 国产一区亚洲一区在线观看| 老司机福利观看| .国产精品久久| 久久精品综合一区二区三区| 九九爱精品视频在线观看| 在线a可以看的网站| 欧美精品国产亚洲| 99在线视频只有这里精品首页| 国产av一区在线观看免费| 久久久精品大字幕| 成人特级黄色片久久久久久久| av在线蜜桃| 亚洲精品乱码久久久久久按摩| 看十八女毛片水多多多| 午夜福利在线观看免费完整高清在 | 99久久精品一区二区三区| 少妇熟女aⅴ在线视频| 日日啪夜夜撸| 国产一区二区三区av在线 | 午夜福利成人在线免费观看| 男女那种视频在线观看| 亚洲av二区三区四区| 变态另类成人亚洲欧美熟女| 99久久久亚洲精品蜜臀av| 国产精品一及| 免费观看的影片在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产大屁股一区二区在线视频| 国产黄片美女视频| 国产高清激情床上av| 国产精品久久久久久av不卡| 国产亚洲欧美98| 午夜久久久久精精品| 精品少妇黑人巨大在线播放 | 小说图片视频综合网站| 最新中文字幕久久久久| 久久韩国三级中文字幕| 99久国产av精品| 久久久久国产网址| 你懂的网址亚洲精品在线观看 | 中国美女看黄片| 欧美3d第一页| 亚洲欧美成人精品一区二区| 日产精品乱码卡一卡2卡三| 我要看日韩黄色一级片| 国产69精品久久久久777片| 伊人久久精品亚洲午夜| 五月伊人婷婷丁香| 免费搜索国产男女视频| 精品久久久久久久久久久久久| 免费观看在线日韩| 国产黄a三级三级三级人| 国产精品久久久久久精品电影小说 | 国产私拍福利视频在线观看| 麻豆成人午夜福利视频| 国产激情偷乱视频一区二区| 悠悠久久av| 亚洲婷婷狠狠爱综合网| 成人高潮视频无遮挡免费网站| 国产精品一区二区三区四区久久| 成人二区视频| 永久网站在线| 嫩草影院入口| 久久精品国产亚洲av涩爱 | 国产成人freesex在线| 少妇的逼好多水| 国产av一区在线观看免费| 国产精品无大码| 好男人在线观看高清免费视频| 国产探花极品一区二区| 男女边吃奶边做爰视频| 一本精品99久久精品77| 26uuu在线亚洲综合色| 免费在线观看成人毛片| 午夜福利高清视频| 国内精品宾馆在线| 国产亚洲精品久久久久久毛片| 国产精品女同一区二区软件| 成人午夜精彩视频在线观看| 亚洲av免费在线观看| 在线观看午夜福利视频| 亚洲自拍偷在线| 少妇的逼好多水| 亚洲国产精品成人久久小说 | 免费观看a级毛片全部| 99久久无色码亚洲精品果冻| 国产在视频线在精品| 国产熟女欧美一区二区| 亚洲内射少妇av| 非洲黑人性xxxx精品又粗又长| 国产成人freesex在线| 国产蜜桃级精品一区二区三区| 亚洲欧美日韩高清专用| 成人午夜精彩视频在线观看| 国产精品无大码| 日日摸夜夜添夜夜爱| 啦啦啦啦在线视频资源| 久久久久网色| 69av精品久久久久久| 免费看av在线观看网站| 亚洲在久久综合| 国产精品无大码| 日本爱情动作片www.在线观看| 九九爱精品视频在线观看| 三级国产精品欧美在线观看| 丰满人妻一区二区三区视频av| av女优亚洲男人天堂| 国产伦在线观看视频一区| 黄片wwwwww| 午夜免费激情av| 欧美潮喷喷水| 亚洲一区二区三区色噜噜| 亚洲av免费在线观看| 少妇人妻精品综合一区二区 | 亚洲色图av天堂| 日韩强制内射视频| 久久这里有精品视频免费| 亚洲精品色激情综合| 亚洲国产精品国产精品| 97人妻精品一区二区三区麻豆| 国产精品一区www在线观看| 亚洲国产精品合色在线| 免费大片18禁| 成人永久免费在线观看视频| 亚洲精品国产成人久久av| 国产精品电影一区二区三区| 99久久精品热视频| 高清毛片免费看| 91av网一区二区| 欧美丝袜亚洲另类| 精品人妻偷拍中文字幕| 亚洲人与动物交配视频| 国产女主播在线喷水免费视频网站 | 欧美日韩在线观看h| 精品一区二区三区视频在线| 不卡一级毛片| 亚洲最大成人手机在线| 久久精品91蜜桃| 18禁在线无遮挡免费观看视频| 国产亚洲精品av在线| 国产成人91sexporn| 边亲边吃奶的免费视频| 人人妻人人澡欧美一区二区| 免费大片18禁| 成熟少妇高潮喷水视频| 欧美+日韩+精品| 亚洲精品日韩av片在线观看| 国产精品女同一区二区软件| 青春草国产在线视频 | 国产亚洲av嫩草精品影院| 99热这里只有是精品在线观看| 欧美激情国产日韩精品一区| 国产一区二区在线观看日韩| 中文字幕熟女人妻在线| 成人av在线播放网站| 看片在线看免费视频| 亚洲精品乱码久久久久久按摩| 99久久中文字幕三级久久日本| 国产中年淑女户外野战色| 变态另类成人亚洲欧美熟女| 成人漫画全彩无遮挡| 天堂网av新在线| 亚洲欧洲国产日韩| 禁无遮挡网站| 看片在线看免费视频| 日本五十路高清| 亚洲av一区综合| 大又大粗又爽又黄少妇毛片口| 国产不卡一卡二| 欧美色视频一区免费| 国产精品国产高清国产av| АⅤ资源中文在线天堂| 中出人妻视频一区二区| 成人毛片a级毛片在线播放| 久久久久久伊人网av| 麻豆av噜噜一区二区三区| 日本成人三级电影网站| 综合色丁香网| 亚洲欧美中文字幕日韩二区| 久久久久久久亚洲中文字幕| 小说图片视频综合网站| 亚洲国产欧美在线一区| 极品教师在线视频| av.在线天堂| 99热只有精品国产| 精品人妻熟女av久视频| 免费一级毛片在线播放高清视频| av天堂在线播放| 欧美成人a在线观看| 日韩在线高清观看一区二区三区| 国产精品,欧美在线| 久久精品夜色国产| 夫妻性生交免费视频一级片| 中文字幕熟女人妻在线| 国产激情偷乱视频一区二区| 99久久中文字幕三级久久日本| 干丝袜人妻中文字幕| 精品一区二区三区人妻视频| 两个人视频免费观看高清| 天堂av国产一区二区熟女人妻| 色综合亚洲欧美另类图片| 日韩一区二区三区影片| 成年女人看的毛片在线观看| 我要搜黄色片| 久久这里有精品视频免费| www.色视频.com| 少妇熟女欧美另类| 亚洲成人久久性| 18禁在线无遮挡免费观看视频| 国产精品一区二区三区四区免费观看| 丝袜美腿在线中文| 全区人妻精品视频| 人妻少妇偷人精品九色| 亚洲精品日韩在线中文字幕 | 国产午夜精品久久久久久一区二区三区| 男人狂女人下面高潮的视频| 少妇被粗大猛烈的视频| 国产精品永久免费网站| 免费不卡的大黄色大毛片视频在线观看 | 狠狠狠狠99中文字幕| 久久鲁丝午夜福利片| 亚洲aⅴ乱码一区二区在线播放| 国产毛片a区久久久久| 美女内射精品一级片tv| 高清午夜精品一区二区三区 | 国产亚洲精品av在线| 亚洲欧洲国产日韩| 熟女人妻精品中文字幕| 伊人久久精品亚洲午夜| 久久久久久久午夜电影| 午夜精品国产一区二区电影 | 亚洲成a人片在线一区二区| 婷婷精品国产亚洲av| 欧美一区二区国产精品久久精品| a级毛色黄片| 亚洲国产高清在线一区二区三| 91久久精品国产一区二区三区| 日韩欧美国产在线观看| 26uuu在线亚洲综合色| 免费在线观看成人毛片| 人妻夜夜爽99麻豆av| 偷拍熟女少妇极品色| 国产日韩欧美在线精品| 国产乱人视频| 中国国产av一级| 国产精品蜜桃在线观看 | 色吧在线观看| 九九爱精品视频在线观看| 3wmmmm亚洲av在线观看| 亚洲乱码一区二区免费版| 国产精品美女特级片免费视频播放器| 日韩欧美精品v在线| 亚洲美女搞黄在线观看| 久久久久久伊人网av| 欧美最黄视频在线播放免费| 九九热线精品视视频播放| 久久精品国产亚洲av涩爱 | 色哟哟哟哟哟哟| 精品无人区乱码1区二区| 哪个播放器可以免费观看大片| 日本在线视频免费播放| 少妇的逼好多水| 欧美日本视频| 18禁在线无遮挡免费观看视频| 亚洲一区二区三区色噜噜| 一级黄片播放器| eeuss影院久久| 在线播放国产精品三级| 亚洲av成人av| 欧美人与善性xxx| 丝袜美腿在线中文| 精品一区二区三区人妻视频| 色视频www国产| 国产精品av视频在线免费观看| 99热这里只有是精品50| 国产精品精品国产色婷婷| 色综合色国产| 麻豆国产av国片精品| 亚洲乱码一区二区免费版| 日本三级黄在线观看| 久久人人爽人人爽人人片va| 中国美白少妇内射xxxbb| 夜夜爽天天搞| 国产精品一区二区性色av| 国产亚洲欧美98| 99久久精品热视频| 国产激情偷乱视频一区二区| 精品免费久久久久久久清纯| 中文亚洲av片在线观看爽| 国产片特级美女逼逼视频| 日韩一本色道免费dvd| 亚洲欧美日韩高清专用| 国产免费一级a男人的天堂| 非洲黑人性xxxx精品又粗又长| 国产精品无大码| 欧美3d第一页| 日本色播在线视频| 久久久久久久久久久丰满| 99在线人妻在线中文字幕| 内射极品少妇av片p| 一本久久中文字幕| 亚洲真实伦在线观看| 日本在线视频免费播放| 久久久久久久午夜电影| 亚洲最大成人中文| 能在线免费看毛片的网站| 又粗又爽又猛毛片免费看| 成人亚洲精品av一区二区| 国产男人的电影天堂91| 国产探花极品一区二区| 久久九九热精品免费| 久久99蜜桃精品久久| 亚洲av电影不卡..在线观看| 五月玫瑰六月丁香| 日韩成人伦理影院| 久久久久九九精品影院| 2021天堂中文幕一二区在线观| 高清毛片免费看| 男女边吃奶边做爰视频| 成年版毛片免费区| 熟女人妻精品中文字幕| avwww免费| 亚洲欧美日韩高清在线视频| 欧美性猛交╳xxx乱大交人| 亚洲经典国产精华液单| 少妇人妻精品综合一区二区 | 麻豆av噜噜一区二区三区| 久久久久久伊人网av| 日韩在线高清观看一区二区三区| 能在线免费看毛片的网站| 神马国产精品三级电影在线观看| 99热这里只有精品一区| 日韩制服骚丝袜av| 尤物成人国产欧美一区二区三区| 久久精品国产亚洲av香蕉五月| 亚洲av电影不卡..在线观看| 国产不卡一卡二| 九九爱精品视频在线观看| 亚洲欧美成人综合另类久久久 | 国产日本99.免费观看| 久久热精品热| 精品欧美国产一区二区三| 人人妻人人澡人人爽人人夜夜 | 可以在线观看的亚洲视频| 大型黄色视频在线免费观看| 热99re8久久精品国产| 波多野结衣高清无吗| 亚洲一区高清亚洲精品| 国产av在哪里看| 欧美成人精品欧美一级黄| 国产高清激情床上av| 国产又黄又爽又无遮挡在线| 亚洲自偷自拍三级| 成人一区二区视频在线观看| 丰满乱子伦码专区| 国产免费男女视频| 日韩亚洲欧美综合| 国产极品天堂在线| 级片在线观看| 少妇裸体淫交视频免费看高清| 久久精品91蜜桃| 99久国产av精品| 亚洲精品粉嫩美女一区| 亚洲经典国产精华液单| 人妻系列 视频| 爱豆传媒免费全集在线观看| 国产成人一区二区在线| 精品欧美国产一区二区三| 日韩中字成人| 午夜视频国产福利| 插阴视频在线观看视频| 少妇人妻一区二区三区视频| 99热这里只有是精品50| 亚洲av熟女| 又粗又爽又猛毛片免费看| 国产乱人偷精品视频| 亚洲精品国产成人久久av| 国产熟女欧美一区二区| 免费看日本二区| 麻豆国产97在线/欧美| 亚洲经典国产精华液单| 亚洲成a人片在线一区二区| av又黄又爽大尺度在线免费看 | 全区人妻精品视频| 午夜爱爱视频在线播放| 国产午夜精品久久久久久一区二区三区| 中文字幕免费在线视频6| 国产精品,欧美在线| 天堂av国产一区二区熟女人妻| 我要搜黄色片| 99久久九九国产精品国产免费| 一区二区三区免费毛片| 国产一级毛片七仙女欲春2| 精品欧美国产一区二区三| 亚洲精品亚洲一区二区| 少妇人妻精品综合一区二区 | 嫩草影院新地址| 国产精品一区二区在线观看99 | 九色成人免费人妻av| 色综合色国产| 国产日本99.免费观看| 亚洲国产欧洲综合997久久,| 网址你懂的国产日韩在线| 美女 人体艺术 gogo| 亚州av有码| 免费无遮挡裸体视频| 久久99热6这里只有精品| 爱豆传媒免费全集在线观看| 国产毛片a区久久久久| 国产日本99.免费观看| 又粗又硬又长又爽又黄的视频 | 日韩成人av中文字幕在线观看| ponron亚洲|