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

    Relationship between South China Sea Precipitation Variability and Tropical Indo-Pacific SST Anomalies in IPCC CMIP5 Models during Spring-to-Summer Transition

    2015-06-09 21:30:01HUWentingandWURenguang
    Advances in Atmospheric Sciences 2015年9期

    HU Wentingand WU Renguang

    1State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029

    2Center for Monsoon System Research,Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029

    Relationship between South China Sea Precipitation Variability and Tropical Indo-Pacific SST Anomalies in IPCC CMIP5 Models during Spring-to-Summer Transition

    HU Wenting1and WU Renguang?2

    1State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029

    2Center for Monsoon System Research,Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029

    The present study evaluates the precipitation variability over the South China Sea(SCS)and its relationship to tropical Indo-Pacific SST anomalies during spring-to-summer transition(April–May–June,AMJ)simulated by 23 Intergovernmental Panel on Climate Change Coupled Model Intercomparison Project Phase 5 coupled models.Most of the models have the capacity to capture the AMJ precipitation variability in the SCS.The precipitation and SST anomaly(SSTA)distribution in the SCS,tropical Pacific Ocean(TPO),and tropical Indian Ocean(TIO)domains is evaluated based on the pattern correlation coefficients between model simulations and observations.The analysis leads to several points of note.First,the performance of the SCS precipitation anomaly pattern in AMJ is model dependent.Second,the SSTA pattern in the TPO and TIO is important for capturing the AMJ SCS precipitation variability.Third,a realistic simulation of the western equatorial Pacific (WEP)and local SST impacts is necessary for reproducing the AMJ SCS precipitation variability in some models.Fourth, the overly strong WEP SST impacts may disrupt the relationship between the SCS precipitation and the TPO–TIO SST. Further work remains to be conducted to unravel the specific reasons for the discrepancies between models and observations in various aspects.

    atmosphere–ocean interaction,South China Sea,CMIP5 models,spring-to-summer transition

    1.Introduction

    Utilizing climate models is a key approach to the understanding of historical climate changes and projections of climate in the future.The evaluation of model performances in various aspects is necessary for unraveling model biases, which then benefits the further development and application of models.Given the availability of the state-of-the-art models involved in the Coupled Model Intercomparison Project Phase 5(CMIP5)for the Intergovernmental Panel on Climate Change(IPCC)Fifth Assessment Report(AR5),we evaluate the precipitation variability over the South China Sea (SCS)and its relationship to tropical Indo-Pacific sea surface temperature anomalies during spring-to-summer transition(April–May–June,AMJ)simulated by 23 coupled models in this study.

    The SCS displays remarkable atmosphere–ocean interaction that is related to the climate variability in the surrounding regions(Wang et al.,1997;Wu and Wang,2001;Wu, 2002,2010;Lestari et al.,2011;Roxy and Tanimoto,2012; He and Wu,2013a,2013b),as it is a part of the tropical eastern Indian Ocean–western Pacific warm pool and the Asian monsoon region.The SCS is a region of frequent tropical cyclone/typhoon activity(Wang et al.,2007;Chen et al.,2012; Feng et al.,2013)and a pathway of moisture transport from the North Indian Ocean and the Southern Hemisphere to East Asia(Tao and Chen,1987;Ding,1994;Wu et al.,2006a), both of which have a great impact on the weather and climate of China,especially South China.The onset of the SCS summer monsoon mostly occurs in May(Wu and Wang,2000, 2001).So,the AMJperiod signifiesthe transition from spring to summer in the SCS and the surrounding regions.Therefore,it is important to evaluate model capacities in terms of the SCS climate variability during AMJ.

    Previous studies have revealed that the atmosphere–ocean relationship in the SCS shows obvious seasonality and regionality.He and Wu(2013b)pointed out that atmosphericforcing is dominant in the northern and central SCS during the local warm season,and that oceanic forcing occurs in the northern SCS during the local cold season.For AMJ in the SCS,which is the season of transition to the warm and rainy season,the atmospheric forcing in the central SCS is characterized by remarkable cloud-radiation and wind-evaporation effects across the central basin,and wind-driven oceanic upwelling effects along the west coasts.Furthermore,the SCS rainfall variability is subject to remote forcing from tropical Indo-Pacific SST anomalies(SSTAs).El Ni?no–Southern Oscillation(ENSO),one of the strongest signals in short-term climate variability,influences the climate variability in remote regions both directly and indirectly(Wu et al.,2009; Wu and Kirtman,2011;Wu et al.,2012;Yun et al.,2013). After El Ni?no dissipates during spring,the tropical Indian Ocean(TIO)SST warming persisting into summer(June–August)forces a Kelvin wave response into the northwest Pacific,which triggers suppressed convection and an anomalous anticyclone there(Xie et al.,2009).In addition,He and Wu (2014)further pointed out that North Indian Ocean SSTAs serve as a medium for an indirect impact of preceding equatorial central Pacific SSTAs on the SCS summer rainfall variability.

    The SCS rainfall variability during AMJ has been found to be influenced by SSTAsin the equatorialPacific(EP),TIO, and western North Pacific(WNP)(Hu et al.,2014).Abovenormal precipitation in the SCS corresponds to lower EP and TIO SST and higher WNP SST.Anomalous precipitation in the SCS occurs via anomalous circulations induced by different combinations of SSTAs in the above regions.Does a link between the SCS precipitation variability to SSTAsin the tropical Indo-Pacific regions during AMJ exist in the CMIP5 models?The main objective of the present study is to understand how climate models perform in this respect.

    The organization of the text is as follows:The data and methods used in this study are described in section 2.In section 3,the leading bias modes in the SCS and the tropical Indo-Pacific are examined.Then,we evaluate the model performance in simulating the AMJ precipitation variability in the SCS and the corresponding precipitation and SSTA pattern in the tropicalIndian and Pacific oceansin section 4.Following that,we classify the models into different types and analyze the relationship between SCS precipitation variability and tropical Indo-Pacific SSTAs during spring-to-summer transition in different types of models separately in section 5. A summary of the results is provided in section 6.

    2.Data and methods

    2.1.Data

    This study utilizes monthly mean precipitation,wind fields and skin temperature from historical simulations of 23 CMIP5 coupled models,which are available online at http:// pcmdi9.llnl.gov/.Due to the large number of experiments included in the CMIP5 framework,only one member for each of the 23 model outputs from the historical experiments is analyzed.Furthermore,the analysis in the present study is based on a common period from 1979 to 2005 to enable a fair comparison with the observations.The information for the selected models is given in Table 1.

    The National Oceanic and Atmospheric Administration (NOAA)extended reconstructed SST,version 3(ERSST3; Smith et al.,2008)is used in the present study.The ERSST3 is available on a 2?×2?grid from 1979 to 2005.We also use monthly mean rainfall from version 2 of the Global Precipitation Climatology Project(Adler et al.,2003).This dataset is available on a 2.5?×2.5?grid from 1979 to 2005.Monthly mean horizontal winds are adopted from the National Centers for Environmental Prediction–U.S.Department of Energy Reanalysis 2(Kanamitsu et al.,2002)on a regular 2.5?×2.5?grid for the period 1979–2005.The reanalysis product is provided by the NOAA/Cooperative Institute forResearch in Environmental Sciences Climate Diagnostics Center,Boulder, Colorado(available via anonymous ftp at ftp://ftp.cdc.noaa. gov/).

    2.2.Methods

    The rainfall and SSTAs in both observations and model outputs have been interpolated to 1?×1?grid for a better coverage of the SCS domain and for convenience of comparison between observations and model simulations.The significance level of correlation is determined based on the Student’s t-test.

    Li and Xie(2012,2014)introduced an intermodel Empirical Orthogonal Function(EOF)analysis to isolate the dominant patterns of model biases.The advantage of this method is maximizing the variance of variable differences among models instead of the temporal variance in conventional EOF analysis(Li and Xie,2014).Here,we carry out a similar analysis to examine the main model biases of the AMJ-mean precipitation.First,the AMJ-mean precipitation is normalized in two dimensions(x–y)for each CMIP5 model.Second,an EOF analysis is applied to each model’s precipitation deviations from the multimodel ensemble(MME)mean for a specific region.The obtained series of each intermodel principal component(IPC)represent the degree of deviation with respect to each intermodel EOF bias pattern in the individual models.

    We use the Taylor diagram(Taylor,2001)to provide a concise statistical summary of how well patterns match with each other in terms of their correlation and the ratio of their standard deviations,which are simply indicated by a single point on the diagram.Through comparison with observations in both pattern and magnitude,we can clearly distinguish the degree of accuracy in the simulation of each model.

    3.Leading bias modes

    First,we use the intermodel EOF analysis described above to obtain the dominant bias patterns of AMJ precipitation,SST and low-level winds over the SCS and the tropical oceans in the models.Figure 1 shows the regressed patterns of precipitation,SST and 850-hPa horizontal winds onto thefirst two IPCs as well as the IPCs with the intermodel EOF analysis applied to the SCS domain.

    The first intermodel EOF mode(Fig.1a),explaining 31.2%of total intermodel variability,features more precipitation and warmer SST in most regions of the SCS,accompanied by an anomalous cyclonic circulation.Models with higher or lower IPC1 values would have larger biases of precipitation in the central SCS and SST in the northern SCS than other normal ones(Fig.1c).The three models with the highest IPC1 values are CanCM4,CanESM2 and HadCM3, while the three models with the lowest IPC1 values are GISSE2-H,MIROC5 and NorESM1-M.The second intermodel EOF mode explains 17.0%of total intermodal variance(Fig. 1b).Thismode hasnegative precipitation anomaly centerslocated in the northeastern and southwestern SCS and a positive center in the southeastern SCS,accompanied by an anomalous cyclonic circulation and warmer SST in the central SCS. The three models with the largest IPC2 values are BCCCSM1.1,CSIRO-MK-3-6-0 and GFDL-ESM2G,and the three models featuring the opposite situation are ACCESS1-0,MIROC5 and MPI-ESM-LR.Both bias patterns show that the positive precipitation anomaly corresponds to a positive SSTA,indicative of the feedback of local SST change on the atmospheric change.However,He and Wu(2013b)pointed out that from April to June,an atmospheric forcing of the ocean with a significantly negative precipitation–SST tendency correlation is dominant in the northern SCS and the central SCS.The bias patterns suggest an inaccurate description of the precipitation–SST relationship over the SCS in some models.

    Hu et al.(2014)suggested a remote forcing of the tropical Indo-Pacific SST to the SCS AMJ rainfall variability. Thus,we further evaluate the main model biases in the tropical Indo-Pacific region during AMJ.Figure 2 shows the leading bias modes of precipitation,SST and low-level winds and the corresponding IPCs in the tropical Indo-Pacific domain (30?S–30?N,40?E–90?W).

    Similar to Fig.2a in the study of Li and Xie(2014),thefirst intermodel EOF mode(Fig.2a),explaining 23.2%of total variability,reflects the equatorial cold tongue SST bias.It features de ficient precipitation and colder SST in the EP,with excessive precipitation and warmer SST on the flanks of the cold tongue.A large area of negative SSTA covers the SCS, the northern Indian Ocean and the northwest of Australia.In contrast,a positive SSTA center is located in the southern Indian Ocean.The three models with the highest IPC1 values are GFDL-ESM2G,IPSL-CM5A-LR and IPSL-CM5AMR,and the three models with the lowest IPC1 values are ACCESS1-0,FGOALS-g2,and MIROC-ESM.

    Table 1.Information on the 23 IPCC AR5 coupled models used in this study.

    The second intermodel EOF mode(Fig.2b),explaining 14.7%of total variance,exhibits positive precipitation anomaly areas in the SCS and the Maritime Continent,in the tropics of the Southern Hemisphere,and in the eastern Paci fic.Furthermore,the most pronounced negative rainfall anomaly center is located in the central and western Paci fic. Positive SSTAs tend to be located in excessive rainfall areas, suggestive of an oceanic forcing of the atmosphere.The three models with the highest(lowest)IPC2 values are CanCM4, CSIRO-MK3-6-0 and HadCM3(FGOALS-s2,GISS-E2-H, and NorESM1-M).

    4.SCS rainfall variability

    Based on the intermodel EOF analysis,excessive precipitation in the central(Fig.1a)and southeastern(Fig.1b)SCS stands out as the most prominent error in some of the CMIP5 models.In this section,the SCS rainfall variability is investigated and evaluated through conventional EOF analysis and Taylor diagrams.The least squares linear trends are removed from all of the variables before the following analysis is conducted.Figure 3 shows the first EOF modes of the AMJ precipitation anomalies over the SCS simulated by 23 IPCC CMIP5 models compared with the observation.

    The first mode accounts for about 40.8%of the total variance in the observation(Fig.3a),which is separated from the other modes according to the rule of North et al.(1982).A large positive loading is located over the central SCS,while a negative loading is confined to the southern coast of China. The percentage variance explained by the leading EOF mode in the models varies from less than 20%(CNRM-CM5)to over 60%(GFDL-ESM2M).About two thirds of the models have the capacity to capture the main features of the observed leading mode(Figs.3b–x).However,some mod-els show remarkable deviations from the observed pattern. The CanESM2 and IPSL-CM5A-LR models display loading that extends to the northern SCS(Figs.3e and q).The distribution of loading in BCC-CSM1.1,IPSL-CM5A-MR, MIROC5,and MIROC-ESM is dissimilar to the observations (Figs.3c,r,t and u).In the GFDL-CM3 model,large loading shifts to both the northern and southern SCS are apparent(Fig.3k).The differences in the pattern of SCS rainfall variability among the models should be related to the AMJ-averaged model errors shown in the previous section.

    To quantify the similarity of the precipitation anomaly distribution,we use the Taylor diagram to compare models with observations based on both the pattern correlation coefficient(PCC)and the ratio of standard deviations(RSD)of the modeled and observed fields in the SCS domain[(0?–25?N,98?–122?E),with land grid points excluded)].The PCC and RSD values are calculated based on the precipitation anomalies obtained by regression with respect to the first conventional principal component(CPC1)of the observation and model,respectively.Since the pattern of observed AMJ precipitation variability may not be represented as EOF1,but may be captured by EOF2 in the models,we also calculate the PCC and RSD between the observation CPC1 and model CPC2.The related Taylor diagram is given in Fig.4.The radial distances from the origin to the points represent the RSD, while the azimuthal positions are determined by the PCC.

    Both the PCC and RSD for the model CPC1(solid dot) show an obvious diversity.For example,the PCC is only 0.136 for BCC-CSM1.1,but reaches 0.855 for NorESM1-M.The value of RSD varies from 0.551(MIROC-ESM)to 1.546(CanCM4).This indicates that the performance in the simulation of the SCS precipitation anomaly pattern is model dependent.The regressed precipitation of seven models(CanESM2,FGOALS-s2,GFDL-CM3,inmcm4,IPSLCM5A-LR,MIROC5,and MIROC-ESM)shows a higher PCC corresponding to the model CPC2(hollow hexagon) than to the model CPC1.For these seven models,the distribution of loading in EOF1 shows obvious differences from observations(Figs.3e,j,k,p,q,r,t and u).As such,the corresponding PCC is relatively low for the model CPC1. There are five models(CCSM4,CSIRO-Mk3-6-0,FGOALS-g2,GFDL-ESM2G,and HadCM3)that have both high a PCC for CPC1 of more than 0.6 and an RSD of around 1,suggestive of a better capacity to simulate the SCS rainfall variability.The PCC and RSD values of the MME for CPC1/CPC2 are 0.465/0.799 and 0.048/1.603,respectively.

    The SCS rainfall anomalies in AMJ are influenced by SST and convection anomalies in both the tropical Indian and Pacific oceans(Hu et al.,2014).Therefore,we evaluate the regressed SST and precipitation anomaly patterns in the tropical Pacific and Indian Ocean regions to the CPC1/2 in the models to understand which regional SSTAs contribute to the SCS AMJ rainfall variability.The evaluation is based on the PCC and RSD values calculated for precipitation and SST in both the tropical Pacific Ocean(TPO;30?S–30?N,120?–270?E)and the TIO(30?S–30?N,40?–100?E).The results for SST and precipitation are given in Fig.5.

    The PCCs of simulated SST and precipitation in the TPO and TIO vary widely among the models,indicative of diversity in model performance.Some models have a high positive PCC corresponding to the model CPC1 in both the TPO and TIO for both precipitation and SST,e.g.,CNRM-CM5, FGOALS-g2,GFDL-ESM2G,GFDL-ESM2M,GISS-E2-H, and IPSL-CM5B-LR.Some models show a high PCC for SST in the TPO(Fig.5d),but a low PCC in the TIO(Fig.5c), e.g.,ACCESS1-0,GFDL-CM3,HadCM3,and inmcm4.Two models(CSIRO-Mk3-6-0 and IPSL-CM5A-LR)show a high PCC for precipitation in the TIO(Fig.5a),but a low PCC in the TPO(Fig.5b).There are two models(CanCM4 and CanESM2)that have large negative PCCs for precipitation in both the TPO and TIO,with the PCC for SST being either negative or small in the TIO.Two models(FGOLAS-s2 and IPSL-CM5A-MR)have negative PCCs for SST in both the TPO and TIO.Two models(MIROC5 and MIROC-ESM)in which EOF2 is more similar than EOF1 to observation EOF1 display large positive PCCs corresponding to CPC2.As for the distribution of model RSD,almost all the models simulate a smaller(larger)magnitude for precipitation and SST in the TPO(TIO)than that in the observation.

    The above method is based on a comparison of the pattern of anomalies.It does not consider the value of SST or precipitation anomalies in each SST key region.An evaluation of SSTAs is also needed.Here,we briefly evaluate the SSTAs in some key regions.Hu et al.(2014)defined three key SST regions based on the distribution of the SST correlation with the SCS precipitation variability during AMJ(shown by the rectangular boxes in Fig.6b)in the observations;namely, the equatorial Indian Ocean(EIO),WNP,and EP.The SSTAs averaged over the EIO,WNP and EP in the observations are?0.130,0.168 and?0.356 K,respectively.

    The above analysesindicate the complexity ofmodelsimulations of the SCS precipitation anomaly pattern and associated SSTAs in the tropical Indo-Pacific regions.For illustrative purposes,we classify the models into different types based on both the PCCs of SST and precipitation in the TPO and TIO and the comparison of SSTAs averaged over the EIO,WNP and EP in models with the observations(Table 2).The classification is determined based on relative values of PCCs and SSTAs.The precise classification of each model may need a strict criterion,which is hard to reach due to thecomplexity involved.The five types of models and the corresponding values of PCCs and SSTAs are listed in Table 2. Note that the BCC-CSM1.1 model is not classified due to low PCCs for both SST and precipitation.

    Table 2.Classification of the models based on two factors:(1)PCCs of SST or precipitation(pr)regression anomalies with respect to CPCs between the observation and models for the SCS,TPO and TIO domains;(2)regressed SSTAs area-averaged over the EIO,WNP and EP.

    In the TPO–TIO type,the PCC for SST or precipitation in the TPO and TIO is large and positive,while the SSTAs in the EIO,WNP and EP are relatively close to the observations.It includes the following ten models:CCSM4,CNRM-CM5, GFDL-ESM2G,GFDL-ESM2M,IPSL-CM5A-LR,IPSLCM5B-LR,MIROC5(CPC2),MIROC-ESM(CPC2),MPIESM-LR,and NorESM1-M.Note that the CCSM4 and NorESM1-M models have relatively low PCCs for SST in the TIO,but they are classi fied as the TPO–TIO type because they have large PCCs for precipitation in the TIO,comparable to that in the TPO.We classify CNRM-CM5 as the TPO–TIO type because it has a large EIO SSTA and EP SSTA, comparable to those in the TPO type,though the SSTA in the WNP is small and negative.The IPSL-CM5A-LR model has a higher PCC corresponding to CPC2(0.319)than CPC1 (0.068)for precipitation in the SCS(Fig.4).Nevertheless,a visual inspection of the distribution of EOF1 and EOF2 loading tells us that the EOF1 resembles more than the EOF2 to the observed EOF1.Furthermore,the PCC corresponding to CPC2 is small for both precipitation and SST in both the TPO and TIO(Fig.5).Also note that the percentage variance explained by CPC2(15.8%)is less than half of that by CPC1(35.0%).Thus,the IPSL-CM5A-LR model is classified as the TPO–TIO type based on large positive PCCs for both precipitation and SST in both the TPO and TIO.The EOF2 loading of the MIROC5 and MIROC-ESM models is more similar than the EOF1 loading to the observed EOF1 loading(Fig.4).The PCC for SST in the TPO or TIO is larger corresponding to CPC2 than CPC1.Thus,MIROC5 and MIROC-ESM are classi fied to the TPO–TIO type based on CPC2.

    In the TPO type,the PCC for SST or precipitation is large in the TPO,but small in the TIO.It contains four models: ACCESS1-0,GFDL-CM3,HadCM3,and inmcm4.In the TIO type,the PCC is large in the TIO but small in the TPO. It includes CSIRO-Mk3-6-0 and GISS-E2-H.In the western EP(WEP)type,the PCC(SST)or PCC(precipitation)in the TPO and TIO are signi ficantly negative.It consists of two models:CanCM4 and CanESM2.As shown later,in the WEP type models,there is an overly strong impact of WEP SST that overcomes the eastern TPO and TIO SST influences on the AMJ SCS precipitation variability.In the local type(FGOALS-g2,FGOALS-s2,IPSL-CM5A-MR,and MRI-CGCM3),the PCCs for SST in both the TPO and TIO are relatively small or even negative,or the SSTAs in the EIO, WNP and EP are relatively small.

    5.Relationship between SCS precipitation variability and tropical Indo-Pacific SSTAs

    In this section,we analyze in detail the precipitation,850-hPa wind,and SSTA distribution in different types of models.For comparison,we first analyze the observed anomaly distribution.The regressed precipitation,850-hPa wind and SST with respect to the normalized CPC1 of AMJ precipitation anomalies for the observation are shown in Fig.6, which is similar to Fig.3 of Hu et al.(2014).Significantpositive precipitation anomalies cover most parts of the SCS, the Indo-China Peninsula,and the WNP,with a cyclonic circulation consisting of northeasterly winds from South China and westerly winds from the Indian Ocean(Fig.6a).Three negative precipitation anomaly centers are located over the western TIO,northwest of Australia,and the equatorial central Pacific.Strong easterlies prevail over the equatorial central Pacific.The SSTA field shows two negative regions over the TIO and the EP,and one positive region over the WNP (Fig.6b).The relevant physical processes of the relationship between AMJ SCS precipitation variability and tropical Indo-Pacific SSTAs are shown by Hu et al.(2014).Overall,both anomalous cross-equatorial flows from the southwestern TIO induced by negative SSTAs there,and an anomalous Walker circulation forced by negative EP SSTAs,contribute to enhanced convection over the SCS and the surrounding regions. Additionalcontribution comes from positive WNP SSTAs via a Rossby wave–type response(Hu et al.,2014).

    For the TPO–TIO type,the simulated patterns of rainfall,850-hPa wind,and SSTAs obtained by regression with respect to the normalized CPC of AMJ precipitation are given in Figs.7 and 8.Seven models(CCSM4,GFDLESM2G,GFDL-ESM2M,IPSL-CM5A-LR,IPSL-CM5BLR,MIROC5 and NorESM1-M)simulate the observed anomalies well(Figs.7a,c,d,e,f,g,j and 8a,c,d,e,f, g,j).For example,negative SSTAs in the TPO are accompanied by below-normal precipitation and anomalous easterlies overthe EP.Negative SSTAs in the southern TIOare followed by negative rainfall anomalies and cross-equatorial southerlies there.Furthermore,positive SSTAs in the WNP are accompanied by cyclonic wind anomalies over the SCS and the Philippine Sea due to the destabilization of the lower troposphere(Wu and Wang,2000,2001;Hu et al.,2014).These features are similar to observations,indicating the contribution of tropical Indian and Pacific Ocean SSTAs to the SCS precipitation variability(Hu et al.,2014).The other three models capture the precipitation and wind anomalies overand surrounding the SCS,but precipitation and wind anomalies over the broad TIO–TPO are relatively weak(Figs.7b,h and i),which is related to relatively small SSTAs(Figs.8b,h andi).The SCS precipitation variability in these three models appears to be mainly influenced by SSTAs in the TIO and WNP.

    In IPSL-CM5A-LR(Figs.7e and 8e),the precipitation, wind,and SSTA patterns are quite similar to observations. One notable discrepancy from observations is that the precipitation anomalies over the SCS shift to the northern part compared to observations.This may suggest a systematic shift in the response to tropical Indo-Paci fic SSTAs in this model.In MIROC5(Figs.7g and 8g),the SSTA pattern is similar to observations,as is the precipitation and wind anomaly distribution.For example,there are anomalous cross-equatorial flows over the TIO and anomalous easterlies over the EP.The precipitation,wind,and SSTA distribution in MIROC-ESM (Figs.7h and 8h)are also similar to observations,but with weaker magnitude.

    For the TPO type,precipitation,850-hPa wind,and SSTAs corresponding to the normalized CPC1 of AMJ precipitation are shown in Fig.9.The SSTA distribution similar to observations is visible in the TPO(right panels of Fig.9). However,the wind and precipitation anomalies over the EP are weak(left panels of Fig.9).Above-normal rainfall and the anomalous cyclone around the SCS may be contributed by positive SSTAs in the WNP via a Rossby wave–type response(Wu et al.,2014).Wind and precipitation anomalies over the TIO are similar to observations,but the SSTAs display an east–west pattern of contrast,leading to a low PCC (SST)in the TIO(Table 2).

    Two models,CSIRO-Mk3-6-0 and GISS-E2-H,are included in the TIO type.The precipitation and wind anomalies in the TIO and western tropical Pacific are similar to observations,as are the SSTAs in the TIO and WNP(Fig.10).However,negative SSTAs in the tropical Pacific are confined to north ofthe equatorin CSIRO-Mk3-6-0,and SSTAs in the EP are weak in both models(right panels of Fig.10).The SCS precipitation and wind anomalies appear to be contributed by both the TIO and WNP SSTAs.

    For the WEP type,prominent SSTAs are present in the TIO and TPO,but different from observations(right panels of Fig.11).The SSTAs in the equatorial central-eastern Pa-cific and TIO are opposite to observations,whereas the WEP SSTAs are similar to observations.Precipitation and wind anomalies over the SCS appear as a northwestward extension of the Rossby wave–type response to WEP SSTAs(left panels of Fig.11).

    For the local type(FGOALS-g2,FGOALS-s2,IPSLCM5A-MR,and MRI-CGCM3),the SSTAs in the TPO and TIO are small(right panels of Fig.12).In FGOALS-s2, prominent precipitation anomalies are present over the SCS and WNP(Fig.12a).These are accompanied by negative SSTAs over the northern SCS and subtropical WNP(Fig. 12b).Thus,it appears that the SSTA is a response to atmospheric change.A similar feature is seen in FGOALS-g2,IPSL-CM5A-MR and MRI-CGCM3,but with weaker anomalies.These results seem to suggest that the SCS precipitation variability is largely related to internal atmosphericdynamics.

    There are three points to note according to the above analysis of anomalies corresponding to SCS precipitation variability during AMJ in different model types.First,the signal of SCS precipitation variability in the TPO and TIO SST and precipitation varies largely among the models.Second,the SSTA pattern in the TPO and TIO is important for capturing the AMJ SCS precipitation variability.And third,a realistic simulation of the WEP and local SST impacts is necessary for reproducing the AMJ SCS precipitation variability.

    Previous studies have revealed that SST is influenced by the atmospheric changes through the cloud-radiation effect, wind-evaporation effect,and wind-induced oceanic processes (Klein et al.,1999;Wu and Kirtman,2007;Wu and Kin-ter III,2010).Meanwhile,SST change could modulate regional convection via lower-level moisture convergence,surface evaporation,and lower tropospheric stability(Lindzen and Nigam,1987;Wu and Wang,2000,2001).The performance of precipitation–SST correlation in the climate models could provide information about whether the physical processes of air–sea interaction in the models are realistic.Positive precipitation–SST correlation indicates an oceanic forcing of precipitation,with ocean surface warming inducing more precipitation(Wu et al.,2006b).By contrast,negative precipitation–SST correlation means an atmospheric forcing of SST,with decreased downward shortwave radiation reaching the surface induced by more precipitation and leading to surface cooling(Trenberth and Shea,2005;Wu et al.,2013). Thus,we compare the SST–precipitation correlation in models against observations during AMJ to understand the biases of the relationship between AMJ SCS precipitation variability and SST key regions.The point-wise and simultaneous precipitation–SST correlations for the observation and the 23 models are given in Fig.13.

    In the observations,positive correlation regions cover the eastern TPO,the WNP,and the southern TIO,while negative correlations appear in the eastern Bay of Bengal and most of the SCS(Fig.13a),which is consistent with He and Wu(2013b).This indicates the forcing of SST in the equatorial central-eastern Pacific,southern TIO,and WNP, supported by numerical experiments(Hu et al.,2014).In turn,anomalous convection and winds may affect the SCS precipitation and wind.In the SCS domain,multiple models(CanESM2,CCSM4,CSIRO-Mk3-6-0,GFDL-ESM2M, GISS-E2-H,HadCM3,inmcm4,IPSL-CM5A-LR,IPSLCM5A-MR,MIROC5,MIROC-ESM,MPI-ESM-LR,MRICGCM3,and NorESM1-M)show a sign of correlation opposite to the observed,indicating unrealistic oceanic forcing to the atmosphere in these models.Also of note is that overly strong WEP SST impacts are seen in CanCM4, CanESM2,CSIRO-Mk3-6-0,HadCM3,inmCM4,IPSLCM5A-LR,IPSL-CM5A-MR,and MIROC5.For example, the positive correlation overthe WEP exceeds 0.8 in CanCM4 (Fig.13d).The warm SST there induces large positive precipitation anomalies that extend northwestward to the SCS (Figs.11a and b).

    6.Summary

    Previous studies have indicated a remote forcing of the SCS precipitation variability during spring-to-summer transition in observations.In the present study,using 23 IPCC AR5 coupled models,we first examine the dominant bias patterns over the SCS and tropical oceans,and then evaluate the precipitation variability over the SCS and its relationship to tropical Indo-Pacific SSTAs in AMJ.The main results can be summarized as follows:

    Based on the intermodel EOF analysis,the leading model biases in the SCS feature excessive precipitation and warmer SST in the central SCS,with an anomalous cyclonic circulation.Moreover,deficient precipitation and colder SST in the EP is detected as the leading model bias mode,plus negative SSTA over the SCS,the northern Indian Ocean,and the northwest of Australia.

    Most of the models have the capacity to capture the AMJ precipitation variability in the SCS.This is related to the TPO–TIO SSTA pattern and its influence on the SCS precipitation and wind.When the models can simulate the SSTAs in the equatorial central-eastern Pacific Ocean,southern TIO, and WNP,either in all three regions or in a combination of two regions,they tend to produce a similar SCS precipitation variability as in observations.This indicates that the SSTA pattern in the TPO and TIO is important for capturing the AMJ SCS precipitation variability.

    There are,however,exceptions.One is the case when there is excessive SST forcing in the WEP that disrupts the relationship between the SCS precipitation variability and the tropical Indo-Pacific SST.In this case,the above-normal precipitation over the SCS is mainly due to the positive SSTA in the WEP.Another case is when there is no apparent SST influence.Here,the AMJ SCS precipitation variability is not affected by the SST conditions over the Indian Ocean and the Pacific,but is seemingly dominated by internal atmospheric dynamics.This indicates that a realistic simulation of regional SST impacts and internal atmospheric variability is necessary for reproducing the relationship between the AMJ SCS precipitation variability and tropical Indo-Pacific SST.

    The model evaluation presented here provides a general understanding about how these models perform in simulating the SCS precipitation variability and its relationship with the SSTAs in the tropical Indo-Pacific regions.Unraveling the reasons why the models perform well or badly calls for a detailed analysis of the physical connections in various aspects. Further studies should focus on the physical processes that connect the SCS to the tropical oceans,the relative impacts of SSTAs in different regions,and the relative contributions of internal atmospheric dynamics and SST forcing to the SCS precipitation variability.

    Acknowledgements.We thank the two anonymous reviewers for their valuable comments.The present study was supported by the National Key Basic Research Program of China(Grant No. 2014CB953902).RW acknowledges the support of the Hong Kong Research Grants Council(Grant No.CUHK403612),the National Natural Science Foundation of China(Grants Nos.41275081 and 41475081),and a Chinese University of Hong Kong direct grant (Grant No.4052057).WH acknowledges the support of a Chinese Academy of Sciences project(Grant No.XDA11010402),the National Natural Science Foundation of China(Grant Nos.41305065 and 41305068),and also the support of the State Key Laboratory of Tropical Oceanography,South China Sea Institute of Oceanology of Chinese Academy of Sciences(Grant No.LTO1203).

    REFERENCES

    Adler,R.F.,and Coauthors,2003:The version-2 Global Precipitation Climatology Project(GPCP)Monthly Precipitation Analysis(1979–present).Journal of Hydrometeorology,4,1147–1167.

    Chen,J.P.,R.G.Wu,and Z.P.Wen,2012:Contribution of South China Sea tropical cyclones to an increase in southern China summer rainfall around 1993.Adv.Atmos.Sci.,29(3),585–598,doi:10.1007/s00376-011-1181-6.

    Ding,Y.H.,1994:Monsoons over China.Kluwer Academic,Dordrecht,419 pp.

    Feng,X.,R.G.Wu,J.P.Chen,and Z.-P.Wen,2013:Factors forinterannual variations of September–October rainfall in Hainan, China.J.Climate,26,8962–8978.

    He,Z.Q.,and R.G.Wu,2013a:Coupled seasonal variability in the South China Sea.J.Oceanogr.,69,57–69.

    He,Z.Q.,and R.G.Wu,2013b:Seasonality of interannual atmosphere–ocean interaction in the South China Sea.J. Oceanogr.,69,699–712.

    He,Z.Q.,and R.G.Wu,2014:Indo-Pacific remote forcing in summer rainfall variability over the South China Sea.Climate Dyn.,42,2323–2337.

    Hu,W.T.,R.G.Wu,and Y.Liu,2014:Relation of the South China Sea precipitation variability to tropical indo-pacific SST anomalies during spring-to-summer transition.J.Climate,27,5451–5467.

    Kanamitsu,M.,W.Ebisuzaki,J.Woollen,S.K.Yang,J.J.Hnilo, M.Fiorino,and G.L.Potter,2002:NCEP-DOE AMIP-II Reanalysis(R-2).Bull.Amer.Meteor.Soc.,83,1631–1643,doi: 10.1175/BAMS-83-11-1631.

    Klein,S.A.,B.J.Soden,and N.C.Lau,1999:Remote sea surface temperature variations during ENSO:Evidence for a tropical atmospheric bridge.J.Climate,12,917–932.

    Lestari,R.K.,M.Watanabe,and M.Kimoto,2011:Role of air–sea coupling in the interannual variability of the South China Sea summer monsoon.J.Meteor.Soc.Japan,89A,283–290.

    Li,G.,and S.-P.Xie,2012:Origins of tropical-wide SST biases in CMIP multi-model ensembles.Geophys.Res.Lett.,39, L22703,doi:10.1029/2012GL053777.

    Li,G.,and S.-P.Xie,2014:Tropical biases in CMIP5 multimodel ensemble:The excessive equatorial Pacific cold tongue and double ITCZ problems.J.Climate,27,1765–1780.

    Lindzen,R.S.,and S.Nigam,1987:On the role of sea surface temperature gradients in forcing low-level winds and convergence in the tropics.J.Atmos.Sci.,44,2418–2436.

    North,G.R.,T.L.Bell,R.F.Cahalan,and F.J.Moeng,1982: Sampling errors in the estimation of empirical orthogonal functions.Mon.Wea.Rev.,110,699–706,doi:10.1175/1520-0493(1982)110<0699:SEITEO>2.0.CO;2.

    Roxy,M.,and Y.Tanimoto,2012:Influence of sea surface temperature on the intraseasonal variability of the South China Sea summer monsoon.Climate Dyn.,39(5),1209–1218.

    Smith,T.M.,R.W.Reynolds,T.C.Peterson,and J.Lawrimore,2008:Improvements to NOAA’s historical merged land-ocean surface temperature analysis(1880–2006).J.Climate,21,2283–2296,doi:10.1175/2007JCLI2100.1.

    Tao,S.Y.,and L.X.Chen,1987:A review of recent research on the East Asian summer monsoon in China.Monsoon Meteorology,C.P.Chang and T.N.Krishnamurti,Eds.,Oxford University Press,Oxford,60–92.

    Taylor,K.E.,2001:Summarizing multiple aspects of model performance in a single diagram.J.Geophys.Res.,106(D7), 7183–7192.

    Trenberth,K.E.,and D.J.Shea,2005:Relationships between precipitation and surface temperature.Geophys.Res.Lett.,32, L14703,doi:10.1029/2005GL022760.

    Wang,D.X.,Z.H.Qin,and F.X.Zhou,1997:Study on the air–sea interaction on the interannual time scale in the South China Sea.Acta Meteorologica Sinica,55(1),33–42.(in Chinese)

    Wang,G.H.,J.L.Su,Y.H.Ding,and D.K.Chen,2007:Tropical cyclone genesis over the South China Sea.J.Mar.Syst.,68, 318–326.

    Wu,B.,T.J.Zhou,and T.Li,2009:Contrast of rainfall-SST relationships in the western North Pacific between the ENSO-developing and ENSO-decaying summers.J.Climate,22, 4398–4405.

    Wu,G.X.,and Coauthors,2006a:The key region affecting the short-term climate variations in China:The joining area of Asia and Indian-Pacific Ocean.Advances in Earth Science,21(11),1109–1118.(in Chinese)

    Wu,R.,and B.Wang,2001:Multi-stage onset of the summer monsoon over the western North Pacific.Climate Dyn.,17(4), 277–289.

    Wu,R.G.,2002:Processes for the northeastward advance of the summer monsoon over the Western North Pacific.J.Meteor. Soc.Japan,80(1),67–83.

    Wu,R.G.,2010:Subseasonal variability during the South China Sea summer monsoon onset.Climate Dyn.,34(5),629–642.

    Wu,R.G.,and B.P.Kirtman,2007:Regimes of seasonal air-sea interaction and implications for performance of forced simulations.Climate Dyn.,29,393–410.

    Wu,R.G.,and B.Wang,2000:Interannual variability of summer monsoon onset over the western North Pacific and the underlying processes.J.Climate,13,2483–2501.

    Wu,R.G.,and J.L.Kinter III,2010:Atmosphere–ocean relationship in the midlatitude North Pacific:Seasonal dependence and east–west contrast.J.Geophys.Res.,115,D06101,doi: 10.1029/2009JD012579.

    Wu,R.G.,and B.P.Kirtman,2011:Caribbean Sea rainfall variability during the rainy season and relationship to the equatorial Pacific and tropical Atlantic SST.Climate Dyn.,37(7–8), 1533–1550.

    Wu,R.G.,B.P.Kirtman,and K.Pegion,2006b:Local air–sea relationship in observations and model simulations.J.Climate,19,4914–4932.

    Wu,R.G.,J.L.Chen,and W.Chen,2012:Different types of ENSO influences on the Indian summer monsoon variability. J.Climate,25(3),903–920.

    Wu,R.G.,J.P.Chen,and Z.-P.Wen,2013:Precipitation-surface temperature relationship in the IPCC CMIP5 models.Adv.Atmos.Sci.,30(3),766–778,doi:10.1007/s00376-012-2130-8.

    Wu,R.G.,G.Huang,Z.C.Du,and K.M.Hu,2014:Cross-season relation of the South China Sea precipitation variability between winter and summer.Climate Dyn.,43(1–2),193–207, doi:10.1007/s00382-013-1820-y.

    Xie,S.P.,K.M.Hu,J.Hafner,H.Tokinaga,Y.Du,G.Huang, and T.Sampe,2009:Indian Ocean capacitor effect on Indo-Western Pacific climate during the summer following El Ni?no.J.Climate,22,730–747.

    Yun,K.-S.,S.-W.Yeh,and K.-J.Ha,2013:Distinct impact of tropical SSTs on summer North Pacific high and western North Pacific subtropical high.J.Geophys.Res.,118,4107–4116, doi:10.1002/jgrd.50253.

    :Hu,W.T.,and R.G.Wu,2015:Relationship between South China Sea precipitation variability and tropical Indo-Pacific SST anomalies in IPCC CMIP5 models during spring-to-summer transition.Adv.Atmos.Sci.,32(9),1303–1318,

    10.1007/s00376-015-4250-4.

    14 November 2014;revised 26 February 2015;accepted 5 March 2015)

    ?Corresponding author:WU Renguang

    Email:renguang@mail.iap.ac.cn

    久久这里只有精品19| 久久精品久久精品一区二区三区| 又黄又粗又硬又大视频| 亚洲av成人不卡在线观看播放网 | 国产成人精品无人区| bbb黄色大片| 亚洲欧美日韩高清在线视频 | 好男人视频免费观看在线| 热re99久久精品国产66热6| 美女中出高潮动态图| 久久中文字幕一级| 菩萨蛮人人尽说江南好唐韦庄| 黄片播放在线免费| 韩国高清视频一区二区三区| 亚洲,欧美,日韩| 最新的欧美精品一区二区| 国产精品免费大片| 亚洲欧美日韩另类电影网站| 日韩精品免费视频一区二区三区| 亚洲欧美一区二区三区黑人| 婷婷丁香在线五月| 97人妻天天添夜夜摸| 搡老岳熟女国产| 国产精品av久久久久免费| 视频区图区小说| av国产精品久久久久影院| 啦啦啦啦在线视频资源| 成年人免费黄色播放视频| 久久精品亚洲av国产电影网| 啦啦啦中文免费视频观看日本| 黄网站色视频无遮挡免费观看| 午夜福利一区二区在线看| 久久99一区二区三区| 欧美黄色淫秽网站| 高清欧美精品videossex| 午夜免费成人在线视频| 中文字幕人妻丝袜制服| 七月丁香在线播放| av不卡在线播放| 精品人妻1区二区| 性高湖久久久久久久久免费观看| 少妇人妻久久综合中文| 欧美日韩国产mv在线观看视频| 国产欧美日韩一区二区三区在线| 成年人免费黄色播放视频| 亚洲天堂av无毛| 色视频在线一区二区三区| 国产精品一二三区在线看| 1024香蕉在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品久久午夜乱码| 晚上一个人看的免费电影| 国产欧美日韩精品亚洲av| 国产色视频综合| av天堂久久9| 国产日韩欧美视频二区| 日本黄色日本黄色录像| 免费不卡黄色视频| 一本久久精品| 欧美黑人精品巨大| 国产精品国产三级专区第一集| www.熟女人妻精品国产| 在线观看www视频免费| a级片在线免费高清观看视频| 亚洲精品自拍成人| 人人妻人人添人人爽欧美一区卜| 免费看av在线观看网站| 久久久久精品人妻al黑| 亚洲国产av影院在线观看| 热re99久久国产66热| 亚洲精品一区蜜桃| 男女下面插进去视频免费观看| 大片电影免费在线观看免费| 亚洲中文字幕日韩| 美女高潮到喷水免费观看| 日韩伦理黄色片| 久久国产亚洲av麻豆专区| 女人被躁到高潮嗷嗷叫费观| 成人亚洲欧美一区二区av| 99久久精品国产亚洲精品| 咕卡用的链子| 国产亚洲av片在线观看秒播厂| 色网站视频免费| 日韩大片免费观看网站| 亚洲五月色婷婷综合| 久久久欧美国产精品| 男女边摸边吃奶| 成人影院久久| 三上悠亚av全集在线观看| www.av在线官网国产| 视频区图区小说| 欧美乱码精品一区二区三区| 午夜免费观看性视频| 视频在线观看一区二区三区| 欧美精品av麻豆av| 国产爽快片一区二区三区| 成人亚洲欧美一区二区av| 可以免费在线观看a视频的电影网站| 久久精品亚洲熟妇少妇任你| 丝袜人妻中文字幕| 脱女人内裤的视频| 嫩草影视91久久| 国产高清不卡午夜福利| 国产精品 欧美亚洲| 亚洲成人免费电影在线观看 | 亚洲成人手机| 午夜福利视频在线观看免费| 最近最新中文字幕大全免费视频 | 啦啦啦 在线观看视频| 亚洲欧美成人综合另类久久久| 99久久99久久久精品蜜桃| 波多野结衣av一区二区av| 国产成人精品无人区| 国产片内射在线| 热re99久久精品国产66热6| 国产在线视频一区二区| 99热国产这里只有精品6| 好男人视频免费观看在线| 欧美av亚洲av综合av国产av| 汤姆久久久久久久影院中文字幕| 国产99久久九九免费精品| 丝瓜视频免费看黄片| 精品国产一区二区三区久久久樱花| 久久久久视频综合| 一边摸一边做爽爽视频免费| 人人妻人人澡人人爽人人夜夜| 女人精品久久久久毛片| 午夜两性在线视频| 国产一区二区激情短视频 | 一级毛片我不卡| 国产真人三级小视频在线观看| 中国美女看黄片| 黄片播放在线免费| 十八禁高潮呻吟视频| 国产免费福利视频在线观看| 久久ye,这里只有精品| 老司机靠b影院| 久久国产精品大桥未久av| 精品欧美一区二区三区在线| 黄网站色视频无遮挡免费观看| 黄色一级大片看看| 免费久久久久久久精品成人欧美视频| 欧美日韩亚洲高清精品| 亚洲国产精品一区二区三区在线| 精品第一国产精品| 男人舔女人的私密视频| 亚洲国产欧美一区二区综合| 真人做人爱边吃奶动态| 久久亚洲精品不卡| 亚洲av成人精品一二三区| 亚洲自偷自拍图片 自拍| 赤兔流量卡办理| 99香蕉大伊视频| 日韩免费高清中文字幕av| 亚洲精品久久成人aⅴ小说| 青草久久国产| 中文字幕精品免费在线观看视频| 国产又色又爽无遮挡免| 午夜老司机福利片| 一级黄片播放器| 欧美日本中文国产一区发布| 国产男人的电影天堂91| 午夜视频精品福利| 免费观看av网站的网址| 校园人妻丝袜中文字幕| 午夜免费观看性视频| 亚洲天堂av无毛| 亚洲精品国产av成人精品| 一区二区日韩欧美中文字幕| 日韩一卡2卡3卡4卡2021年| 黄色视频在线播放观看不卡| 在线精品无人区一区二区三| 国产成人啪精品午夜网站| 老汉色∧v一级毛片| 成人亚洲欧美一区二区av| 多毛熟女@视频| 国产一级毛片在线| 大型av网站在线播放| 久久综合国产亚洲精品| 天天影视国产精品| 在线观看免费午夜福利视频| 亚洲精品一卡2卡三卡4卡5卡 | 欧美xxⅹ黑人| 中文字幕av电影在线播放| 国产xxxxx性猛交| 久久99精品国语久久久| av不卡在线播放| 欧美黑人欧美精品刺激| av网站在线播放免费| 在线观看www视频免费| 亚洲九九香蕉| 校园人妻丝袜中文字幕| 精品人妻在线不人妻| 国产欧美日韩综合在线一区二区| 国产亚洲av高清不卡| 久久久国产精品麻豆| 免费久久久久久久精品成人欧美视频| 精品亚洲乱码少妇综合久久| 久久久亚洲精品成人影院| 国产av国产精品国产| 免费少妇av软件| 香蕉国产在线看| 欧美人与性动交α欧美软件| 国产成人一区二区在线| 色综合欧美亚洲国产小说| 天堂俺去俺来也www色官网| 国产高清不卡午夜福利| 高潮久久久久久久久久久不卡| 国产野战对白在线观看| 观看av在线不卡| 99国产精品一区二区三区| 亚洲久久久国产精品| 一本—道久久a久久精品蜜桃钙片| 一边摸一边做爽爽视频免费| 丝袜美腿诱惑在线| 亚洲,一卡二卡三卡| 一区福利在线观看| 国产一区二区在线观看av| 国产亚洲av高清不卡| 叶爱在线成人免费视频播放| 脱女人内裤的视频| 久久亚洲国产成人精品v| 欧美黄色淫秽网站| 婷婷色av中文字幕| 国产精品久久久久久人妻精品电影 | 国产熟女欧美一区二区| 亚洲,一卡二卡三卡| 久久亚洲精品不卡| 日韩中文字幕视频在线看片| 久久人人97超碰香蕉20202| 一级片免费观看大全| 啦啦啦视频在线资源免费观看| 欧美黑人精品巨大| 成年女人毛片免费观看观看9 | 这个男人来自地球电影免费观看| 99国产精品一区二区三区| 啦啦啦 在线观看视频| 国产成人啪精品午夜网站| 黑人巨大精品欧美一区二区蜜桃| 亚洲天堂av无毛| 国产精品久久久久久精品电影小说| 久久鲁丝午夜福利片| 午夜激情久久久久久久| 新久久久久国产一级毛片| 成人手机av| 老司机深夜福利视频在线观看 | 国产男人的电影天堂91| 亚洲精品中文字幕在线视频| 1024视频免费在线观看| 无限看片的www在线观看| 亚洲久久久国产精品| 熟女少妇亚洲综合色aaa.| 亚洲图色成人| 亚洲一卡2卡3卡4卡5卡精品中文| 一区二区三区激情视频| 男的添女的下面高潮视频| 亚洲欧美日韩高清在线视频 | 免费黄频网站在线观看国产| 视频区图区小说| 亚洲成人免费电影在线观看 | 午夜两性在线视频| av有码第一页| 91国产中文字幕| 欧美中文综合在线视频| 国产成人精品久久二区二区免费| 嫩草影视91久久| 美女高潮到喷水免费观看| 国产精品一国产av| 高清不卡的av网站| 国产成人一区二区在线| av天堂久久9| 精品人妻1区二区| 欧美精品人与动牲交sv欧美| 久久久精品94久久精品| 丰满少妇做爰视频| 美女午夜性视频免费| 久久精品国产亚洲av涩爱| 中文欧美无线码| 黑人欧美特级aaaaaa片| 好男人电影高清在线观看| 黄频高清免费视频| 色综合欧美亚洲国产小说| 九色亚洲精品在线播放| 午夜福利,免费看| 最新的欧美精品一区二区| 久久精品国产综合久久久| 精品人妻在线不人妻| 十八禁人妻一区二区| 亚洲av片天天在线观看| 国产亚洲精品第一综合不卡| 电影成人av| 亚洲 欧美一区二区三区| 精品亚洲乱码少妇综合久久| 交换朋友夫妻互换小说| 国产高清videossex| 国产精品99久久99久久久不卡| 亚洲人成网站在线观看播放| 国产伦人伦偷精品视频| av网站在线播放免费| 午夜激情av网站| 国产伦理片在线播放av一区| 亚洲人成电影免费在线| 亚洲av片天天在线观看| 免费少妇av软件| 天天躁夜夜躁狠狠久久av| 久久热在线av| 婷婷成人精品国产| 女人久久www免费人成看片| 成人国语在线视频| 国产精品久久久久久精品古装| 电影成人av| 欧美国产精品va在线观看不卡| 亚洲av日韩精品久久久久久密 | 国产在线观看jvid| 国产精品人妻久久久影院| 99久久人妻综合| 日本午夜av视频| 丰满饥渴人妻一区二区三| 一级片免费观看大全| 免费日韩欧美在线观看| 亚洲一码二码三码区别大吗| 欧美乱码精品一区二区三区| 亚洲九九香蕉| 又黄又粗又硬又大视频| 男女之事视频高清在线观看 | 日韩,欧美,国产一区二区三区| 国产视频首页在线观看| 久久天堂一区二区三区四区| 国产视频首页在线观看| 宅男免费午夜| 国产亚洲午夜精品一区二区久久| 国产精品香港三级国产av潘金莲 | 亚洲男人天堂网一区| 亚洲精品国产色婷婷电影| 久久狼人影院| 夫妻午夜视频| 久久天堂一区二区三区四区| 亚洲av电影在线观看一区二区三区| 久久久精品免费免费高清| 2018国产大陆天天弄谢| 欧美黑人精品巨大| 国产高清videossex| 久热这里只有精品99| 久久久久国产一级毛片高清牌| 91成人精品电影| 亚洲欧美成人综合另类久久久| 亚洲国产成人一精品久久久| 精品国产乱码久久久久久小说| 飞空精品影院首页| 波多野结衣一区麻豆| 婷婷色av中文字幕| 中文字幕人妻熟女乱码| 久久久久久久国产电影| 黄色片一级片一级黄色片| 每晚都被弄得嗷嗷叫到高潮| 又大又黄又爽视频免费| av线在线观看网站| 少妇人妻久久综合中文| av在线老鸭窝| 婷婷丁香在线五月| 精品人妻在线不人妻| 国产一区亚洲一区在线观看| 国产成人欧美| 啦啦啦 在线观看视频| 免费黄频网站在线观看国产| 免费女性裸体啪啪无遮挡网站| 免费久久久久久久精品成人欧美视频| 亚洲精品中文字幕在线视频| 国产福利在线免费观看视频| 免费av中文字幕在线| 乱人伦中国视频| 最新的欧美精品一区二区| 一级毛片我不卡| 看免费成人av毛片| 亚洲精品国产一区二区精华液| 免费看av在线观看网站| 香蕉丝袜av| 久久热在线av| 精品人妻1区二区| 久久久亚洲精品成人影院| 一级毛片女人18水好多 | 久久天躁狠狠躁夜夜2o2o | 岛国毛片在线播放| 国产男人的电影天堂91| 久久久精品区二区三区| 国产男女超爽视频在线观看| 男女国产视频网站| 新久久久久国产一级毛片| 中文字幕人妻丝袜制服| 亚洲色图 男人天堂 中文字幕| www日本在线高清视频| 国产成人系列免费观看| 最黄视频免费看| 黄网站色视频无遮挡免费观看| 国产成人精品久久久久久| 久久国产精品人妻蜜桃| 日本五十路高清| 99国产精品一区二区蜜桃av | 一区二区三区激情视频| 婷婷色综合大香蕉| 日本五十路高清| 亚洲精品乱久久久久久| 亚洲黑人精品在线| 男人爽女人下面视频在线观看| 国产女主播在线喷水免费视频网站| 久久久久网色| 国产黄色免费在线视频| 免费看av在线观看网站| 久久久精品国产亚洲av高清涩受| 少妇精品久久久久久久| 欧美另类一区| 十八禁网站网址无遮挡| 男男h啪啪无遮挡| 国产精品av久久久久免费| 日本欧美国产在线视频| 色精品久久人妻99蜜桃| 久久精品久久久久久久性| 国产成人影院久久av| 另类精品久久| 国产麻豆69| 精品人妻一区二区三区麻豆| 欧美日韩亚洲国产一区二区在线观看 | 中文乱码字字幕精品一区二区三区| 亚洲图色成人| 国产精品久久久久久人妻精品电影 | 在线观看人妻少妇| 成人黄色视频免费在线看| 亚洲欧美清纯卡通| 91精品三级在线观看| 国产精品 欧美亚洲| 国产日韩欧美亚洲二区| 国产一区二区三区av在线| 亚洲午夜精品一区,二区,三区| 欧美激情极品国产一区二区三区| 欧美精品av麻豆av| 国产免费又黄又爽又色| 亚洲伊人久久精品综合| 中文字幕亚洲精品专区| 免费高清在线观看视频在线观看| 亚洲精品中文字幕在线视频| 在线观看免费午夜福利视频| 老司机影院成人| 老鸭窝网址在线观看| 久久精品久久久久久久性| 亚洲人成电影免费在线| 欧美av亚洲av综合av国产av| 国产人伦9x9x在线观看| www.av在线官网国产| 亚洲国产最新在线播放| 欧美在线一区亚洲| 老汉色∧v一级毛片| 精品久久久久久电影网| 久久狼人影院| 色视频在线一区二区三区| 久久精品人人爽人人爽视色| 亚洲精品日本国产第一区| 视频在线观看一区二区三区| 亚洲精品国产av蜜桃| 国产精品亚洲av一区麻豆| 人妻 亚洲 视频| 50天的宝宝边吃奶边哭怎么回事| 午夜福利在线免费观看网站| 国产精品久久久久久精品电影小说| 欧美日韩成人在线一区二区| 亚洲av成人不卡在线观看播放网 | 一边摸一边抽搐一进一出视频| 亚洲激情五月婷婷啪啪| 999久久久国产精品视频| 亚洲精品乱久久久久久| 久久久欧美国产精品| 日日夜夜操网爽| 叶爱在线成人免费视频播放| 久久国产精品影院| 国产成人av教育| 免费少妇av软件| 国产高清国产精品国产三级| 女人久久www免费人成看片| av不卡在线播放| 欧美亚洲 丝袜 人妻 在线| 午夜精品国产一区二区电影| 色综合欧美亚洲国产小说| 黄片小视频在线播放| 99久久综合免费| 交换朋友夫妻互换小说| 丝袜在线中文字幕| 精品久久蜜臀av无| 欧美人与善性xxx| 欧美日韩国产mv在线观看视频| av网站免费在线观看视频| 超碰97精品在线观看| 如日韩欧美国产精品一区二区三区| 天堂俺去俺来也www色官网| 亚洲五月婷婷丁香| 成人影院久久| 韩国高清视频一区二区三区| 欧美精品一区二区免费开放| 日韩中文字幕欧美一区二区 | 国产99久久九九免费精品| 男人添女人高潮全过程视频| 十分钟在线观看高清视频www| 蜜桃国产av成人99| 久久免费观看电影| 美女国产高潮福利片在线看| 免费女性裸体啪啪无遮挡网站| 国产精品久久久久成人av| 国产伦人伦偷精品视频| 欧美性长视频在线观看| 黄频高清免费视频| 老司机亚洲免费影院| 成年美女黄网站色视频大全免费| 亚洲精品国产av成人精品| 免费观看av网站的网址| 不卡av一区二区三区| 欧美xxⅹ黑人| 亚洲精品日本国产第一区| 国产女主播在线喷水免费视频网站| 美女国产高潮福利片在线看| 久久99一区二区三区| 国产精品久久久久久精品电影小说| 纵有疾风起免费观看全集完整版| a级毛片黄视频| 伊人亚洲综合成人网| 秋霞在线观看毛片| 两人在一起打扑克的视频| 1024香蕉在线观看| 人妻人人澡人人爽人人| 国产黄频视频在线观看| 国产精品麻豆人妻色哟哟久久| 久久ye,这里只有精品| 天天添夜夜摸| 亚洲专区国产一区二区| 性少妇av在线| 热99久久久久精品小说推荐| 精品卡一卡二卡四卡免费| 中国美女看黄片| 精品亚洲乱码少妇综合久久| 在线观看国产h片| 精品久久久久久久毛片微露脸 | 亚洲国产成人一精品久久久| 国产男女内射视频| 成人亚洲精品一区在线观看| 一二三四在线观看免费中文在| 一区二区av电影网| 午夜福利一区二区在线看| 精品高清国产在线一区| 一级片免费观看大全| 日韩制服丝袜自拍偷拍| 一级毛片电影观看| 成人亚洲欧美一区二区av| 麻豆乱淫一区二区| 亚洲精品乱久久久久久| 亚洲精品美女久久av网站| 亚洲午夜精品一区,二区,三区| 免费在线观看黄色视频的| 人妻 亚洲 视频| 好男人视频免费观看在线| 嫁个100分男人电影在线观看 | e午夜精品久久久久久久| 国产精品久久久久久精品古装| 在线观看免费视频网站a站| 男女国产视频网站| 亚洲国产精品999| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产日韩一区二区| 欧美乱码精品一区二区三区| 亚洲一码二码三码区别大吗| 欧美精品亚洲一区二区| 18在线观看网站| 欧美日韩黄片免| 麻豆国产av国片精品| 亚洲精品国产区一区二| 国产欧美日韩综合在线一区二区| 国产在线一区二区三区精| 国产人伦9x9x在线观看| 狠狠婷婷综合久久久久久88av| 黄色 视频免费看| 亚洲精品国产一区二区精华液| 国产免费现黄频在线看| 亚洲三区欧美一区| 自拍欧美九色日韩亚洲蝌蚪91| 日本黄色日本黄色录像| 九色亚洲精品在线播放| 七月丁香在线播放| 久久久久精品人妻al黑| 日韩 亚洲 欧美在线| 两人在一起打扑克的视频| 在线精品无人区一区二区三| 欧美日韩精品网址| e午夜精品久久久久久久| www.999成人在线观看| a 毛片基地| 亚洲伊人久久精品综合| 亚洲,欧美,日韩| 下体分泌物呈黄色| av片东京热男人的天堂| 精品人妻熟女毛片av久久网站| 黄色视频在线播放观看不卡| 两人在一起打扑克的视频| 日本黄色日本黄色录像| 色视频在线一区二区三区| 永久免费av网站大全| 国产在线一区二区三区精| 黄色a级毛片大全视频| 伊人亚洲综合成人网| videos熟女内射| 国产av国产精品国产| 18在线观看网站| 日韩视频在线欧美| 王馨瑶露胸无遮挡在线观看| 亚洲色图 男人天堂 中文字幕| 一区二区日韩欧美中文字幕| 女人爽到高潮嗷嗷叫在线视频| 波多野结衣一区麻豆| 国产又色又爽无遮挡免| www日本在线高清视频|