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

    Tropical Cyclone Genesis Potential Index over the Western North Pacific Simulated by CMIP5 Models

    2015-09-04 02:49:05SONGYajuanWANGLeiLEIXiaoyanandWANGXidong
    Advances in Atmospheric Sciences 2015年11期

    SONG YajuanWANG LeiLEI Xiaoyanand WANG Xidong

    1College of Physical and Environmental Oceanography,Ocean University of China,Qingdao266100

    2First Institute of Oceanography,State Oceanic Administration,Qingdao266061

    3State Key Laboratory of Tropical Oceanography,South China Sea Institute of Oceanology, Chinese Academy of Sciences,Guangzhou510301

    4Key Laboratory of Marine Environment Information Technology,State Oceanic Administration,Tianjin300171

    Tropical Cyclone Genesis Potential Index over the Western North Pacific Simulated by CMIP5 Models

    SONG Yajuan1,2,WANG Lei3,LEI Xiaoyan?2,and WANG Xidong4

    1College of Physical and Environmental Oceanography,Ocean University of China,Qingdao266100

    2First Institute of Oceanography,State Oceanic Administration,Qingdao266061

    3State Key Laboratory of Tropical Oceanography,South China Sea Institute of Oceanology, Chinese Academy of Sciences,Guangzhou510301

    4Key Laboratory of Marine Environment Information Technology,State Oceanic Administration,Tianjin300171

    Tropical cyclone(TC)genesis over the western North Pacific(WNP)is analyzed using 23 CMIP5(Coupled Model Intercomparison Project Phase 5)models and reanalysis datasets.The models are evaluated according to TC genesis potential index(GPI).The spatial and temporal variations of the GPI are first calculated using three atmospheric reanalysis datasets (ERA-Interim,NCEP/NCAR Reanalysis-1,and NCEP/DOE Reanalysis-2).Spatial distributions of July–October-mean TC frequency based on the GPI from ERA-interim are more consistent with observed ones derived from IBTrACS global TC data.So,the ERA-interim reanalysis dataset is used to examine the CMIP5 models in terms of reproducing GPI during the period 1982–2005.Although most models possess deficiencies in reproducing the spatial distribution of the GPI,their multimodel ensemble(MME)mean shows a reasonable climatological GPI pattern characterized by a high GPI zone along 20°N in the WNP.There was an upward trend of TC genesis frequency during 1982 to 1998,followed by a downward trend.Both MME results and reanalysis data can represent a robust increasing trend during 1982–1998,but the models cannot simulate the downward trend after 2000.Analysis based on future projection experiments shows that the GPI exhibits no significant change in the first half of the 21st century,and then starts to decrease at the end of the 21st century under the representative concentration pathway(RCP)2.6 scenario.Under the RCP8.5 scenario,the GPI shows an increasing trend in the vicinity of 20°N,indicating more TCs could possibly be expected over the WNP under future global warming.

    tropical cyclone,genesis potential index,CMIP5,western North Pacific,global warming.

    1.Introduction

    Tropical cyclones(TCs),known as typhoons in the western North Pacific(WNP),are one of the most disastrous weather phenomena influencing agriculture,human life,and property(Chan,2005).TC landfall,accompanied by strong wind and heavy rainfall,can result in windstorms and flood damage.There has been growing interest in TCs because of their large socioeconomic impact.In particular,anticipated future changes in TC genesis frequency under global warming scenarios have raised concern.

    With the development of computing resources and numerical models,general circulation models(GCMs)are widely employed for simulating TCs.Recent GCMs with high resolution can simulate most aspects of typhoons(Zhao et al.,2009).There are several approaches for identifying TCs in GCMs,such as examining TC statistical features in the models and analyzing the large-scale environmental variables associated with TC activities.Clearly,not all models can simulate TC characteristics realistically.Therefore,some downscaling methods have been applied to represent the spatial and temporal characteristics of TCs(Villarini and Vecchi, 2013).

    A series of Coupled Model Intercomparison Projects (CMIPs),promoted by the World Climate Research Program, have provided a favorable platform to evaluate the performance of climate models in a systematic manner.Furthermore,these projects have provided a great opportunity to compare and understand the characteristics of natural phenomena under different climate conditions.Phase 5 of the CMIP(CMIP5)datasets have been employed to evaluate model abilities in simulating TCs,and to predict variation in TCs under global warming scenarios(Taylor et al.,2012). Although most of the models in CMIP5,especially for longterm integration,do not have sufficiently high resolution forreproducing specific TC activities,the large-scale environmental variables associated with TC genesis can be reproduced in the models.Therefore,we focus in this study on these large-scale environmental variables associated with TC generation,rather than on individual cyclones.

    Understanding the influence of large-scale environmental factors on tropical cyclone genesis is an important problem.In the absence of a comprehensive theory,empirical methods have been used to represent the relationship between large-scale environmental factors and tropical cyclone genesis.There are a few TC indices based on large-scale environmental factors,such as the yearly genesis parameter(Gray, 1979),the modified yearly convective genesis potential index (Royer et al.,1998),and the genesis potential index(GPI) (Emanuel and Nolan,2004).The GPI summarizes the environmental factors influencing TC genesis such as low-level vorticity,vertical wind shear,mid-level air humidity,ocean temperature,and a conditionally unstable atmosphere.Such an index represents climatological spatial and seasonal distribution of TC genesis,and the relative contribution of individual large-scale environmental factors on TC genesis may be different in different regions and on different timescales.To date,theGPIhasbeenwidelyusedtoinvestigatethevariation of TC genesis frequency related to climate changes(Menkes et al.,2012),as well as in the western Pacific(Zhang et al., 2010;Zhao et al.,2012;Deng et al.,2014).The GPI was developed to fit the spatial distribution and annual cycle of TC genesis at the hemispheric scale(Camargo et al.,2007),and has been employed to depict the temporal variation of TCs on the seasonal scale and in future climate(Yokoi et al.,2009).It has been shown that low-resolution models can simulate the GPI better than the genesis of TCs(Walsh et al.,2013).Considering the insufficient resolution of most CMIP5 models, evaluating the GPI is a more effective way to infer variations in TC activities.

    Some studies,using coarse-resolution models,have shown there to be a decreasing trend in tropical storm frequency and a slight increase in its intensity(e.g.,McDonald et al.,2005).However,another experiment with a coarseresolution atmospheric model(T42,approximately 300 km grid spacing)indicated that the change in TC frequency was insignificant(Tsutsui,2002).According to the Intergovernmental Panel on Climate Change Fourth Assessment Report, for models with coarse resolutions,there is no consistent evidence of large variations in either TC frequency or intensity in global warming simulation(Solomon et al.,2007). The variations in TC frequency and intensity under global warming scenarios rely on the model resolutions and physical parameterizations.Considering that CMIP5 models possess better parameterization schemes and higher horizontal resolutions than CMIP3 models,there is a need to assess CMIP5 model performances in reproducing TC activities.

    In this study,because different regions have different large-scale environmental factors and air–sea interaction processes determining TC frequency and intensity(Wang et al., 2012;Wang et al.,2014),we focus on the WNP,which is the most TC-active region globally.We examine how well GCMs simulate the GPI and evaluate the future changes in large-scale environmental factors linked to TCs.This study aims to evaluate model performances in simulating the climatological GPI related to TCs in the WNP,and to investigate the changes in large-scale environmental conditions for TC genesis frequency under global warming,based on the CMIP5 models.The geographic domain of the WNP used in this study is(0°–40°N,100°–160°E).The rest of the paper is organized as follows:Section 2 introduces the models, reanalysis datasets and the GPI used in this study.Section 3 describes the temporal and spatial variations of the GPI in CMIP5 models,in historical and future projection simulations.Finally,a discussion and conclusions are provided in section 4.

    2.Data and methods

    2.1.Observation and reanalysis datasets

    To evaluate model performances,we compare the simulated results with the International Best Track Archive for Climate Stewardship(IBTrACS)global tropical cyclone data for the 24-year period of 1982–2005(Knapp et al.,2010; www.ncdc.noaa.gov/oa/ibtracs/),which is endorsed by the World Meteorological Organization(WMO).The dataset contains the best-estimated TC positions,maximum sustained wind,and minimum central pressure.We also use the latest global reanalysis data produced by the European Centre for Medium-Range Weather Forecasts:ERA-Interim (Simmons et al.,2007);the National Center for Environmental Prediction–National Center for Atmospheric Research (NCEP–NCAR)Reanalysis-1(Kalnay et al.,1996)and NCEP/Department of Energy(DOE)Reanalysis-2 datasets (Kanamitsu et al.,2002);and the NOAA(National Oceanic and Atmospheric Administration)Optimum Interpolation (OI)Sea Surface Temperature(SST)V2(Reynolds et al., 2002).Atmospheric humidity,wind vectors,air temperature, and SST are derived from these datasets to construct the climatological summer mean(July–October)GPI.The monthly datasets used in this study include absolute vorticity at 850 hPa,relative humidity at 600 hPa,and vertical wind shear between 850 hPa and 200 hPa.In addition,the air temperature and humidity data at all levels in the troposphere are used to calculate the potential intensity.

    2.2.Model outputs

    The CMIP5 models are used in this study to evaluate model performances in simulating the climatological GPI related with TCs in the WNP.The 20th century simulations with 24-year periods from 23 state-of-the-art coupled models are compared with observational data.All models considered in this study are listed in Table 1,which provides information about their host centers and the resolutions of their atmospheric models.The historical simulations of the CMIP5 models are forced by observed greenhouse gas concentrations,solar radiation forcing,and aerosols during the last 150 years.The historical model outputs are availablefrom 1850 to 2005.To investigate changes of the GPI under future global warming scenarios,multi-model-based future projections from the representative concentration pathway(RCP)scenarios are analyzed.Two scenario datasets, denoted as RCP8.5 and RCP2.6,which indicate the highest and lowest radiative forcing values in the year 2100 relative to pre-industrial values,are analyzed to show the variation in TC genesis under the influence of atmospheric greenhouse gases,air pollutants,land use,and land cover(Meinshausen et al.,2011;Taylor et al.,2012).The GCMs have more than one member available;only the firstrealization ofeach model is used in this study.Finally,all the model outputs with different spatial resolutions are interpolated to the same grid (1°×1°)using the bilinear interpolation method.

    Table 1.Climate models whose outputs are analyzed in this study.

    2.3.Method

    The large-scale environmental conditions for TC generation can be evaluated by the GPI.The GPI has been widely employed to represent TC activities over the WNP,and it facilitates an easier comparison with the results from CMIP3 (Yokoi et al.,2009).The four terms of the GPI are shown below:

    whereηis the absolute vorticity(s?1)at 850 hPa,Vshear(m s?1)represents the magnitude of the vertical wind shear between 850 and 200 hPa,andH(%)is the relative humidity at 600 hPa.Vpot(m s?1),representing the maximum TC potential intensity(Emanuel,1986,1988,2000)is calculated as follows:

    whereTSis the ocean surface temperature,Tois the mean outflow temperature,CKis the exchange coef ficient for enthalpy, andCDis the drag coef ficient.is the saturation equivalent potential temperature at the ocean surface,andθeis the boundary layer equivalent potential temperature.The potential intensity is indicated as the upper bound of intensity that a TC may reach under appropriate thermodynamic and dynamic conditions.The four terms representing atmosphericvorticity,vertical wind shear,air humidity,and atmospheric potential intensity play crucial roles in TC formation.The absolute vorticity and vertical wind shear are considered as dynamic components,and the relative humidity and potential intensity are thermodynamic components.

    3.Results

    3.1.Observed TC frequency and the GPI in reanalysis datasets

    We select the TCs whose genesis locations are over the WNP based on the IBTrACS dataset and calculate the spatial distributions of observed TC frequency(Fig.1a).From the WNP to the South China Sea,the characteristics of the atmosphere and ocean provide favorable conditions for TC generation during July–October.From Fig.1a,the TC formation frequency shows a band of high values between 10°N and 20°N,with an elliptic pattern from southeast to northwest in the WNP.There are two centers of high values:one over the sea area east of the Philippines and the other over the South China Sea.Three reanalysis datasets are used to calculate the GPI.As shown in Fig.1b,the distribution of the GPI based on ERA-Interim has a high value patch over the WNP along 10°N.Correlation analysis shows that the GPI based on ERA-Interim corresponds to observed TC generation frequency with a correlation coefficient(CC)of 0.64. NCEP Reanalysis-1 and Reanalysis-2 seem to produce lower GPI values than ERA-Interim over the WNP,but show larger GPI offshore in Chinese seas.The GPI of NCEP Reanalysis-2 shows poor agreement with observed TC frequency,and its GPI has a lower CC with TC formation frequency.From these results,there are large inconsistencies in the distributions of GPI for different reanalysis datasets.The GPI of ERA-Interim has the highest correlation with observed TCs in terms of frequency distribution;therefore,we use ERAInterim for the model evaluations in this study.

    3.2.GPI from CMIP5 models

    The GPI values based on the 23 CMIP5 models in the period 1982–2005 are illustrated in Fig.2.We can use them to examine the ability of these models in reproducing the geographical pattern of large-scale variables.Most models can simulate the GPI pattern characterized by the highvalue zone along 20°N,except FGOALS-g2,HadGEM2-CC, HadGEM2-ES and MIROC5.The models present a wide range of GPI values.Seven models including ACCESS1.0, CanESM2,FGOALS-g2,GISS-E2-H,HadCM3,HadGEM2-CC,and HadGEM2-ES have smaller GPI values than the other models.As shown in Fig.2,the largest GPI values in most models are located east of the Philippines.The highvalue zones in FGOALS-g2,HadCM3,HadGEM2-CC and HadGEM2-ES are closer to the equator than those in the other models.Meanwhile,four models including CanESM2, HadGEM2-AO,MPI-ESM-LR and MPI-ESM-MRI exhibit relatively small GPI values near the equator.The last panelin Fig.2 is the multi-model ensemble(MME)mean of the 23 model simulations.The GPI of the MME mean exhibits a zonal distribution along 20°N,which is consistent with the observed TC distributions.The conspicuous maximum of the GPI indicates that the air–sea environment provides favorable conditions for active TC frequency over the WNP along 20°N.

    From Fig.2,it is noted that the spatial distributions of the GPI simulated by the 23 CMIP5 models have large differences.The CCs and root-mean-square errors(RMSEs) between the GPI from models and the GPI from ERAInterim–OISST reanalysis data are calculated for each of the 23 CMIP5 models(Fig.3).The models with CCs greater than 0.5 and RMSEs less than 0.37 are selected to produce the“best model ensemble”and to carry out analysis of future global warming projections.In total,11 climate models are selected:BCC-CSM1-1,CCSM4,FIO-ESM, GFDL-ESM2G,HadGEM2-AO,INMCM4,IPSL-CM5ALR,IPSL-CM2A-MR,MPI-ESM-LR,MPI-ESM-MR and MRI-CGCM3.The statistical analysis shows that BCCCSM1-1 can reproduce well the distribution of the GPI with the highest CC,but its RMSE is still larger than the MME mean.The RMSE in the MME is smaller and its CC is relatively higher than most individual models,except for BCCCSM1-1 and MPI-ESM-MR.This indicates that the MME performs reasonably well in simulating the GPI characteristics over the WNP.Although some differences in the GPI amplitude exist between the MME mean and the reanalysis data,the general patterns of the GPI in the MME mean are in agreement with those in the reanalysis data.

    The GPI differences between CMIP5 models and ERAInterim–OISST reanalysis data are shown in Fig.4.The GPI values in BCC-CSM1-1,CCSM4,FIO-ESM,HadGEM2-AO,IPSL-CM5A-LR and IPSL-CM5A-MR are much lower than those in the reanalysis data.Some models show negative differences in the northern part of the WNP and positive differences in the southern part,such as FIO-ESM,GFDLESM2G,INMCM4,IPSL-CM5A-LR,MPI-ESM-LR,MPIESM-MR and MRI-CGCM3.

    Both the GPI and TC number exhibit significant temporal variation(Fig.5).If we examine the long-term linear trends(as shown in Table 2),it is found that there was an upward trend of TC genesis frequency during 1982to the late 1990s,followed by a downward trend.Both ERA-Interim and NCEP Reanalysis-1showrobust increasing trends during 1982–1999,with linear trend values of 0.13 and 0.16,respectively.We also find a positive trend in the MME mean, with a trend value of 0.03.The 90%confidence interval of the MME in 1982–1999 was from?0.04 to 0.11.For the observed TC number,the 90%confidence interval is from?0.29 to 0.61during 1982–1999.However,the MME and ERA-Interim show upward trends in 2000–2005,which are not consistent with the trend of observed TC number.These results indicate that the variability of observed TC number and that of ERA-Interim GPI are not always consistent.The GPI is not always able to present the real TC number variations due to the limitation of the GPI and the complexity of TC activities.On the other hand,the GPI from the MME and individual models do not have to follow the observed real TC number variability,considering that the upward and downward trends of observed TC number may mainly be parts of long-term natural variability and the possible biases in individual model simulations.Due to the limitations of the GPI and the model biases,it should be noted that there are many uncertainties to assessing future change in TC number using the GPI from the CMIP5 models.For 1982–1999,7 out of 11 models show the same upward trend as observed:CCSM4, INMCM4,IPSL-CM5A-LR,IPSL-CM5A-MR,MPI-ESMLR,MPI-ESM-MR and MRI-CGCM3.The temporal variation of the GPI in the MME mean before 1999 is similar to observed,but it has a significant increasing trend after 1999 that is inconsistent with observations.

    So,what causes the GPI biases in these models?The GPIis composed of four terms;therefore,we examine those individual four terms in more detail in Fig.6.The first term, representing atmospheric absolute vorticity,is almost the same in the MME mean and the reanalysis data.The second term,associated with vertical wind shear,is smaller in the MME than in the two reanalysis datasets in the area around Taiwan Island.Since the second term of the GPI has an inverse relationship with the vertical wind shear,it means that the simulated vertical wind shear is larger along 20°N than those of ERA-Interim and NCEP Reanalysis-1,which hinders TC generation in this region.The third term,associated with mid-level humidity,shows a robustly high value in the tropical region from the equator to 20°N in the MME.Previous studies also suggest that models usually produce a higherrelative humidity than observed(Walsh et al.,2013).In this study,the humidity in the climate models is also higher than that in ERA-Interim.NCEP Reanalysis-1 shows a lower humiditycomparedwith theMMEmean and ERA-Interim.The fourth term is for potential intensity,with higher values in the low-latitude area in ERA-Interim and NCEP Reanalysis-1.The potential intensity in the MME mean is smaller than those in ERA-Interim and NCEP Reanalysis-1.In general, the second term,representing vertical wind shear,and the fourth term,associated with potential intensity,show lower values in the MME mean than the two reanalysis datasets; whereas,the atmospheric humidity in the MME mean is larger than those in the reanalysis data.By comparing the GPI from ERA-Interim and NCEP Reanalysis-1,lower GPI values in NCEP Reanalysis-1 are mainly due to lower midlevel humidity and weaker potential intensity.If we compare the GPI in the MME mean and ERA-Interim,smaller GPI values in the MME are mainly induced by weaker potential intensity.

    Table 2.Trends and 90%confidence intervals(in square brackets)of the GPI derived from different datasets and time periods.

    Potential changes in GPI due to global warming are evaluated using the RCP2.6 and RCP8.5 scenarios(Fig.7). We note that the GPI values in the two scenarios have strong interannual variations.In the RCP8.5 scenario,there is an obvious upward trend in the period 2006–2095;the trend value is 0.02 and the 90%confidence interval is from 0.01 to 0.033.This shows a high level of statistical significance.Eight out of eleven models have positive trends under RCP8.5:BCC-CSM1-1,FIO-ESM,GFDL-ESM2G, HadGEM2-AO,INMCM4,IPSL-CM5A-LR,IPSL-CM5AMR and MRI-CGCM3.It is suggested that global warming may possibly provide favorable large-scale environmental conditions for TC formation over the WNP.On the other hand,thetimeseriesofGPIundertheRCP2.6scenarioshows a negative trend of?0.002,but this is not statistically signif icant.Under the RCP2.6 scenario,the GPI exhibits no signif icant trend before 2060 in the 21st century,and only decreases at the end of the 21st century.

    A point of interest is the spatial characteristics of the GPI under global warming.We use 11 models including BCCCSM1-1,CCSM4,FIO-ESM,GFDL-ESM2G,HadGEM2-AO,INMCM4,IPSL-CM5A-LR,IPSL-CM2A-MR,MPIESM-LR,MPI-ESM-MR and MRI-CGCM3 under the highest scenario RCP8.5to demonstrate the change.As shown in Fig.8,the differences in GPI between RCP8.5(2076–2100) and historical simulations(1982–2005)in the WNP have two positive value patches along 15°N and 25°N.In the South China Sea,the GPI has a downward trend in the north(north of 17°N)and an upward trend in the south(between 10°N and 17°N).The changes in the GPI result in an out-of-phase relationship between TCs formed in the north and south.In the WNP,the GPI differences show a positive anomaly with two centers of large value in the region east of Taiwan from 20°N to 30°N and east of the Philippines from 10°N to 20°N, where typhoons frequently occur.The result suggests that the WNP,especially from 10°N to 30°N,may be expected to have more cyclone genesis under the global warming scenarios.The correlation coefficients between each model’s difference and the MME results are calculated.Eight of eleven models have CCs greater than 0.5,and the CCs of five models are greater than 0.8.Correlation analysis shows that most models project the same pattern as the MME under RCP8.5.The differences in the four terms between RCP8.5 and historical simulations are shown in Fig.9.For absolute vorticity,there is a positive anomaly north of 20°N and a negative anomaly to the south(Fig.9a).The distributional characteristics of the second term,related to vertical wind shear, show similar patterns,with a positive anomaly lying further south,extending from the South China Sea to the WNP.Both the absolute vorticity and vertical wind shear are out of phase in the zonal direction in the WNP.This means that TCs would increase in the northern part of the WNP and decrease in the southern part of the WNP.In addition,both relative humidity and potential intensity differences produce positive values over almost the whole of the WNP.The differences in the GPI between RCP8.5 and historical simulations have larger correlation coef ficients with the differences of the third term, indicating that the change in the GPI under RCP8.5 is mainlycaused by changes in relative humidity.

    4.Discussion and conclusions

    This study has assessed the GPI over the WNP in CMIP5 models.Three reanalysis datasets including ERA-Interim, NCEP Reanalysis 1,and NCEP/DOE Reanalysis 2 were first evaluated using the IBTrACS best track data.Based on the climatology of the 20th century,most of the 23 models still show deficiencies in terms of geographical patterns of the GPI in the WNP and produce a wide range of index values over time.Their MME mean reasonably represents the GPI distributions characterized by a large-value region from the South China Sea to the western Pacific,which are consistent with those in the observation and reanalysis products. Eleven out of 23 models with higher CCs and lower RMSEs were used to present the GPI variation over time.The GPI exhibits an increasing trend in both reanalysis data and the MME mean during the period 1982–1999.Analysis based on future projection experiments showed that the GPI under the RCP2.6 scenario has no significant change in the first half of the 21st century and then decreases at the end of the 21st century,while the GPI under RCP8.5 exhibits an increasing trend.The GPI differences between RCP8.5 and historical results showed that there are bands of positive values along 15°N and 25°N,suggesting that more TCs could possibly be expected over the WNP under future global warming,particularly over this two-band region.By examining the four terms in the GPI,it was found that the changes in the GPI under RCP8.5 are mainly caused by changes in relative humidity.

    Many models,owing to their low resolutions,produce much weaker and larger TCs(Camargo et al.,2005).Although the model resolutions in CMIP5 are higher compared with the GCMs in CMIP3,they are still insufficient for reproducing TC structures and spatiotemporal distributions.Considering the insufficient resolution of most CMIP5 models, evaluating the GPI is a more effective way to infer variations in TC activities.This study has provided an evaluation of CMIP5 model performances in simulating the GPI in the WNP and an investigation of the changes in large-scale environmental conditions associated with TC genesis under global warming.However,it should be noted that there are many uncertainties in future TC projection from the GPI.The GPI is not always able to present the actual TC variations due to the limitation of the GPI and the complexity of TC activities.In this study,information on the GPI under RCP2.6 and RCP8.5 has been presented,which provides only some clues and possible information on future TC projection.Although TC frequency cannot be determined by the GPI in anabsolutely accurate way,it could still provide some helpful information and clues to understanding the possible changes in TC activities.

    Acknowledgements.This study was supported by the National Basic Research Program of China(973 Program)(Grant No. 2013CB430304),the Scientific Research Foundation of the First Institute of Oceanography,the State Oceanic Administration(Grant No.GY0213G19),and the National Natural Science Foundation of China(GrantNos.41205026and41206026).WLwasalsosupported by the Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.XDA11010104)and the Knowledge Innovation Program of the Chinese Academy of Sciences(Grant No. SQ201208).The authors wish to thank the ESG-PCMDI for providing the CMIP5 model outputs.

    REFERENCES

    Bentsen,M.,and Coauthors,2013:The Norwegian Earth System Model,NorESM1-M-Part1:Descriptionand basicevaluation of the physical climate.Geoscientific Model Development,6, 687–720.

    Bi,D.H.,and Coauthors,2013:The ACCESS coupled model: Description,control climate and evaluation.Australian Meteorological and Oceanographic Journal,63,41–64.

    Camargo,S.J.,A.G.Barnston,and S.E.Zebiak,2005:A statistical assessment of tropical cyclone activity in atmospheric general circulation models.Tellus A,57,589–604.

    Camargo,S.J.,K.A.Emanuel,and A.H.Sobel,2007:Use of a genesis potential index to diagnose ENSO effects on tropical cyclone genesis.J.Climate,20,4819–4834.

    Chan,J.C.L.,2005:Interannual and interdecadal variations of tropical cyclone activity over the western North Pacific.Meteor.Atmos.Phys.,89,143–152.

    Chylek,P.,J.Li,M.K.Dubey,M.Wang,and G.Lesins,2011:Observed and model simulated 20th century Arctic temperature variability:Canadian Earth System Model CanESM2.Atmos. Chem.Phys.,11,22893–22907.

    Collins,W.J.,and Coauthors,2011:Development and evaluation of an Earth-System model-HadGEM2.Geoscientific Model Development,4,1051–1075.

    Deng,S.R.,L.G.Wu,R.F.Wang,and C.Wang,2014:Analysis of tropical cyclone genesis potential index with multiple reanalysis data.Chinese J.Meteor.Sci.,34,243–251.(in Chinese)

    Dunne,J.P.,and Coauthors,2012:GFDL’s ESM2 global coupled climate-carbon earth system models.Part I:Physical formulation and baseline simulation characteristics.J.Climate,25, 6646–6665.

    Emanuel,K.A.,1986:An air-sea interaction theory for tropical cyclones.Part I:Steady-state maintenance.J.Atmos.Sci.,43, 585–605.

    Emanuel,K.A.,1988:The maximum intensity of hurricanes.J. Atmos.Sci.,45,1143–1155.

    Emanuel,K,2000:A statistical analysis of tropical cyclone intensity.Mon.Wea.Rev.,128,1139–1152.

    Emanuel,K.A.,and D.S.Nolan,2004:Tropical cyclone activity and the global climate system.Preprints,26th Conf.on Hurricanes and Tropical Meteorology,Miami,FL,Amer.Meteor.Soc.,240–241.

    Gent,P.R.,and Coauthors,2011:The community climate system model version 4.J.Climate,24,4973–4991.

    Giorgetta,M.A.,and Coauthors,2013:Climate and carbon cycle changes from 1850 to 2100 in MPI-ESM simulations for the Coupled Model Intercomparison Project phase 5.Journal of Advances in Modeling Earth Systems,5,572–597.

    Gray,W.M.,1979:Hurricanes:Their formation,structure and likely role in the tropical circulation.Meteorology over the Tropical Oceans,D.B.Shaw,Ed.,Royal Meteor.Soc.,James Glaisher House,Grenville Place,Bracknell,Berkshire,155–218.

    Hourdin,F.,and Coauthors,2013:LMDZ5B:The atmospheric component of the IPSL climate model with revisited parameterizations for clouds and convection.Climate Dyn.,40, 2193–2222.

    Ji,D.,and Coauthors,2014:Description and basic evaluation of Beijing Normal University Earth System Model(BNU-ESM) version 1.Geoscientific Model Development,7,2039–2064.

    Kalnay,E.,and Coauthors,1996:The NCEP/NCAR 40-year reanalysis project.Bull.Amer.Meteor.Soc.,77,437–470.

    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.

    Knapp,K.R.,M.C.Kruk,D.H.Levinson,H.J.Diamond,and C. J.Neumann,2010:The International best Track archive for climate stewardship(IBTrACS).Bull.Amer.Meteor.Soc.,91, 363–376.

    Li,L.J.,and Coauthors,2013:The flexible global oceanatmosphere-landsystemmodel,Grid-pointVersion2: FGOALS-g2.Adv.Atmos.Sci.,30,543–560,doi:10.1007/ s00376-012-2140-6.

    McDonald,R.E.,D.G.Bleaken,D.R.Cresswell,V.D.Pope,and C.A.Senior,2005:Tropical storms:Representation and diagnosis in climate models and the impacts of climate change.Climate Dyn.,25,19–36.

    Meinshausen,M.,and Coauthors,2011:The RCP greenhouse gas concentrations and their extensions from 1765 to 2300.Climatic Change,109,213–241.

    Menkes,C.E.,M.Lengaigne,P.Marchesiello,N.C.Jourdain,E. M.Vincent,J.Lef`evre,F.Chauvin,and J.-F.Royer,2012: Comparison of tropical cyclogenesis indices on seasonal to interannual timescales.Climate Dyn.,38,301–321.

    Pope,V.D.,M.L.Gallani,P.R.Rowntree,and R.A.Stratton, 2000:The impact of new physical parametrizations in the Hadley Centre climate model:HadAM3.Climate Dyn.,16, 123–146.

    Qiao,F.,Z.Song,Y.Bao,Y.Song,Q.Shu,C.Huang,and W. Zhao,2013:Development and evaluation of an earth system model climate with surface gravity waves.J.Geophys.Res. Oceans,118,4514–4524.

    Reynolds,R.W.,N.A.Rayner,T.M.Smith,D.C.Stokes,and W.Q.Wang,2002:An improved in situ and satellite SST analysis for climate.J.Climate,15,1609–1625.

    Royer,J.-F.,F.Chauvin,B.Timbal,P.Araspin,and D.Grimal, 1998:A GCM study of the impact of greenhouse gas increase on the frequency of occurrence of tropical cyclones.Climatic Change,38,307–343.

    Shindell,D.T.,and Coauthors,2012:Interactive ozone and methane chemistry in GISS-E2 historical and future climate simulations.Atmos.Chem.Phys.,13,2653–2689.

    Simmons,A.J.,S.M.Uppala,D.P.Dee,and S.Kobayashi,2007: ERA-Interim:New ECMWF reanalysis products from 1989onwards.ECMWF Newsletter,110,25–35.

    Solomon,S.,D.Qin,M.Manning,Z.Chen,M.Marquis,K.B. Averyt,M.Tignor,and H.L.Miller,2007:Climate Change 2007:The Physical Science Basis.Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change,Cambridge University Press,Cambridge,UnitedKingdomandNewYork,NY,USA, 996 pp.

    Taylor,K.E.,R.J.Stouffer,and G.A.Meehl,2012:An overview of CMIP5 and the experiment design.Bull.Amer.Meteor. Soc.,93,485–498.

    Tjiputra,J.F.,C.Roelandt,M.Bentsen,D.M.Lawrence,T. Lorentzen,J.Schwinger,?.Seland,and C.Heinze,2013: Evaluation of the carbon cycle components in the Norwegian Earth System Model(NorESM).Geoscientific Model Development,6,301–325.

    Tsutsui,J.,2002:Implications of anthropogenic climate change for tropical cyclone activity:A case study with the NCAR CCM2.J.Meteor.Soc.Japan,80,45–65.

    Villarini,G.,and G.A.Vecchi,2013:Projected increases in North Atlantic tropical cyclone intensity from CMIP5 models.J. Climate,26,3231–3240.

    Volodin,E.M.,N.A.Diansky,and A.V.Gusev,2010:Simulating present-day climate with the INMCM4.0 coupled model of the atmospheric and oceanic general circulations.Izvestiya. Atmospheric and Oceanic Physics,46,414–431.

    Walsh,K.,S.Lavender,E.Scoccimarro,and H.Murakami, 2013:Resolutiondependenceoftropicalcycloneformationin CMIP3 and finer resolution models.Climate Dyn.,40,585–599.

    Wang,X.,W.Zhou,C.Y.Li,and D.X.Wang,2012:Effects of the East Asian summer monsoon on tropical cyclone genesis over the South China Sea on an interdecadal timescales.Adv.Atmos.Sci.,29,249–262,doi:10.1007/s00376-011-1080-x.

    Wang,X.,W.Zhou,C.Y.Li,and D.X.Wang,2014:Comparison of the impact of two types of El Nin?o on tropical cyclone genesis over the South China Sea.Int.J.Climatol.,34,2651–2660.

    Watanabe,M.,and Coauthors,2010:Improved climate simulation by MIROC5:Mean states,variability,and climate sensitivity.J.Climate,23,6312–6335.

    Xin,X.G.,T.W.Wu,J.L.Li,Z.Z.Wang,W.P.Li,and F.H.Wu, 2013:How well does BCC CSM1.1 reproduce the 20th century climate change over China?.Atmos.Oceanic Sci.Lett., 6,21–26.

    Yokoi,S.,Y.N.Takayabu,and J.C.L.Chan,2009:Tropical cyclone genesis frequency over the western North Pacific simulated in medium-resolution coupled general circulation models.Climate Dyn.,33,665–683.

    Yukimoto,S.,and Coauthors,2011:Meteorological Research Institute-Earth System Model Version 1(MRI-ESM1)-Model Description.Technical Reports of the Meteorological Research Institute No.64,Meteorological Research Institute, Tsukuba,Japan,83 pp.

    Zhang Y.,H.J.Wang,J.Q.Sun,and H.Drange,2010:Changes in the tropical cyclone genesis potential index over the western North Pacific in the SRES A2 scenario.Adv.Atmos.Sci.,27, 1246–1258,doi:10.1007/s00376-010-9096-1.

    Zhao,J.P.,L.G.Wu,and H.K.Zhao,2012:Improvement of tropical cyclone genesis potential index in the western North Pacific Basin.Journal of the Meteorological Sciences,32,591–599.(in Chinese)

    Zhao,M.,I.M.Held,S.-J.Lin,and G.A.Vecchi,2009:Simulations of global hurricane climatology,interannual variability,and response to global warming using a 50-km resolution GCM.J.Climate,22,6653–6678.

    :Song,Y.J.,L.Wang,X.Y.Lei,and X.D.Wang,2015:Tropical cyclone genesis potential index over the western North Pacific simulated by CMIP5 models.Adv.Atmos.Sci.,32(11),1539–1550,

    10.1007/s00376-015-4162-3.

    13 November 2014;revised 7 April 2015;accepted 5 May 2015)

    ?Corresponding author:LEI Xiaoyan

    Email:annylxy@fio.org.cn

    婷婷丁香在线五月| 观看免费一级毛片| 69人妻影院| 国产蜜桃级精品一区二区三区| 精品久久久久久久久亚洲 | 亚洲自偷自拍三级| 日韩欧美一区二区三区在线观看| 神马国产精品三级电影在线观看| 一区二区三区激情视频| 欧美精品啪啪一区二区三区| 欧美日韩国产亚洲二区| av黄色大香蕉| 极品教师在线视频| 夜夜夜夜夜久久久久| 国产精品美女特级片免费视频播放器| aaaaa片日本免费| 最新在线观看一区二区三区| 最新在线观看一区二区三区| 亚洲精品456在线播放app | 黄色日韩在线| 国产av麻豆久久久久久久| 精品99又大又爽又粗少妇毛片 | 国产精品98久久久久久宅男小说| 中文字幕高清在线视频| 国产激情偷乱视频一区二区| 午夜精品一区二区三区免费看| 国产大屁股一区二区在线视频| 国产aⅴ精品一区二区三区波| 精品午夜福利视频在线观看一区| 久99久视频精品免费| 亚洲av免费高清在线观看| 午夜福利在线观看吧| 国内少妇人妻偷人精品xxx网站| 99热这里只有是精品50| 在线免费观看不下载黄p国产 | 看黄色毛片网站| 国产精品一区二区免费欧美| 久久精品国产99精品国产亚洲性色| 国产亚洲av嫩草精品影院| av在线天堂中文字幕| 国产欧美日韩精品亚洲av| 3wmmmm亚洲av在线观看| 国产三级中文精品| 精品一区二区三区视频在线观看免费| 熟妇人妻久久中文字幕3abv| 精品久久久噜噜| 黄色欧美视频在线观看| av黄色大香蕉| 国产在视频线在精品| 99久久成人亚洲精品观看| 内地一区二区视频在线| 精品99又大又爽又粗少妇毛片 | 国产免费av片在线观看野外av| 欧美黑人欧美精品刺激| 国产成年人精品一区二区| 99国产极品粉嫩在线观看| 日韩精品中文字幕看吧| 欧美中文日本在线观看视频| 两人在一起打扑克的视频| 久久久国产成人免费| 春色校园在线视频观看| 高清毛片免费观看视频网站| 国产精品久久久久久av不卡| 丝袜美腿在线中文| 嫩草影院新地址| 听说在线观看完整版免费高清| 欧美日韩中文字幕国产精品一区二区三区| 亚洲熟妇中文字幕五十中出| 日韩中字成人| 老熟妇仑乱视频hdxx| 日本-黄色视频高清免费观看| 国产午夜福利久久久久久| .国产精品久久| 老熟妇乱子伦视频在线观看| 亚洲经典国产精华液单| 国产一区二区三区视频了| 麻豆久久精品国产亚洲av| 国产一级毛片七仙女欲春2| 色哟哟哟哟哟哟| 日韩人妻高清精品专区| 男人舔女人下体高潮全视频| 国产私拍福利视频在线观看| 国产精品,欧美在线| 深夜精品福利| 精品午夜福利在线看| 黄色视频,在线免费观看| 国产男人的电影天堂91| 久久久久久久午夜电影| 色综合站精品国产| 欧美性猛交黑人性爽| 淫秽高清视频在线观看| 搡老妇女老女人老熟妇| 免费电影在线观看免费观看| 少妇丰满av| 五月伊人婷婷丁香| 日日撸夜夜添| 国产私拍福利视频在线观看| 国产主播在线观看一区二区| 国产三级在线视频| 人人妻人人澡欧美一区二区| 亚洲乱码一区二区免费版| 91av网一区二区| 国产精品无大码| 美女大奶头视频| 国产一区二区三区av在线 | 欧美色视频一区免费| 国产真实乱freesex| 简卡轻食公司| 18禁黄网站禁片午夜丰满| 非洲黑人性xxxx精品又粗又长| 1024手机看黄色片| 亚洲va在线va天堂va国产| 给我免费播放毛片高清在线观看| 国产精华一区二区三区| 午夜爱爱视频在线播放| 欧美xxxx性猛交bbbb| 成人性生交大片免费视频hd| 成人av在线播放网站| 在线免费观看不下载黄p国产 | 99久久久亚洲精品蜜臀av| 国产三级中文精品| 久久精品久久久久久噜噜老黄 | 色在线成人网| 亚洲真实伦在线观看| 悠悠久久av| 22中文网久久字幕| 日韩欧美精品v在线| 午夜免费成人在线视频| 亚洲av中文av极速乱 | 在现免费观看毛片| 乱系列少妇在线播放| 国产成人aa在线观看| 国产精品美女特级片免费视频播放器| 久久久久九九精品影院| 亚洲精华国产精华精| 亚洲七黄色美女视频| 两个人视频免费观看高清| 搡老妇女老女人老熟妇| 制服丝袜大香蕉在线| 偷拍熟女少妇极品色| 中文字幕熟女人妻在线| 听说在线观看完整版免费高清| 噜噜噜噜噜久久久久久91| 欧美性猛交黑人性爽| 国产精品日韩av在线免费观看| 色综合站精品国产| 国产精品爽爽va在线观看网站| 嫩草影院精品99| 久久久久久久午夜电影| 欧美一区二区国产精品久久精品| 有码 亚洲区| 国产成年人精品一区二区| a级毛片免费高清观看在线播放| 成年女人永久免费观看视频| 国产一区二区亚洲精品在线观看| 久久久久久国产a免费观看| 偷拍熟女少妇极品色| 国产精品嫩草影院av在线观看 | 欧美绝顶高潮抽搐喷水| 99久久久亚洲精品蜜臀av| 亚洲 国产 在线| 男人舔奶头视频| 国产男靠女视频免费网站| 久久人人精品亚洲av| 少妇人妻一区二区三区视频| av国产免费在线观看| 婷婷色综合大香蕉| 日韩精品有码人妻一区| 国产精品久久久久久精品电影| 他把我摸到了高潮在线观看| 免费av毛片视频| 亚洲精品色激情综合| 国产国拍精品亚洲av在线观看| 亚洲av第一区精品v没综合| 一边摸一边抽搐一进一小说| 人妻少妇偷人精品九色| 亚洲第一电影网av| 又黄又爽又刺激的免费视频.| 成人欧美大片| 网址你懂的国产日韩在线| 波多野结衣高清无吗| 亚洲精品日韩av片在线观看| 久久中文看片网| 极品教师在线免费播放| 俺也久久电影网| 欧美成人一区二区免费高清观看| 成人无遮挡网站| 在线免费观看不下载黄p国产 | av在线老鸭窝| 国产 一区 欧美 日韩| 国产伦精品一区二区三区四那| 色5月婷婷丁香| 在线观看免费视频日本深夜| 麻豆久久精品国产亚洲av| av在线观看视频网站免费| 国产蜜桃级精品一区二区三区| 狂野欧美激情性xxxx在线观看| 91久久精品国产一区二区成人| 一区二区三区四区激情视频 | 91狼人影院| 日韩,欧美,国产一区二区三区 | 亚洲,欧美,日韩| 国产精品一区二区三区四区免费观看 | 亚洲精品久久国产高清桃花| 男人和女人高潮做爰伦理| 欧美丝袜亚洲另类 | 亚洲成人免费电影在线观看| 永久网站在线| 男人舔奶头视频| 国产不卡一卡二| 久久久成人免费电影| 在线观看av片永久免费下载| 国产高潮美女av| 国产视频一区二区在线看| 热99在线观看视频| 亚洲精品粉嫩美女一区| 欧美日韩中文字幕国产精品一区二区三区| 噜噜噜噜噜久久久久久91| 欧美xxxx黑人xx丫x性爽| 精品人妻熟女av久视频| 成年免费大片在线观看| 老司机福利观看| 久久久久精品国产欧美久久久| 午夜精品一区二区三区免费看| 亚洲性久久影院| 成年女人看的毛片在线观看| 人人妻,人人澡人人爽秒播| 综合色av麻豆| 精品无人区乱码1区二区| 中文字幕av在线有码专区| 深夜a级毛片| 亚洲性夜色夜夜综合| 在线免费十八禁| 99精品在免费线老司机午夜| 国语自产精品视频在线第100页| 日本与韩国留学比较| 中文资源天堂在线| 欧美人与善性xxx| 国产高清视频在线播放一区| 免费av毛片视频| 欧美在线一区亚洲| 中文字幕免费在线视频6| 亚洲精华国产精华液的使用体验 | 一级毛片久久久久久久久女| 在线看三级毛片| 午夜免费男女啪啪视频观看 | 日韩欧美在线乱码| 一个人看视频在线观看www免费| 国语自产精品视频在线第100页| 国产乱人视频| 丰满的人妻完整版| 一区二区三区免费毛片| 国产真实乱freesex| 午夜福利成人在线免费观看| 又黄又爽又刺激的免费视频.| 亚洲无线观看免费| 在线观看免费视频日本深夜| 国产av在哪里看| 我的女老师完整版在线观看| 国产精品自产拍在线观看55亚洲| videossex国产| 在线看三级毛片| 最近最新中文字幕大全电影3| 在线播放无遮挡| 久久香蕉精品热| 亚洲美女黄片视频| 亚洲国产欧洲综合997久久,| 亚洲七黄色美女视频| 麻豆成人午夜福利视频| 国产精品久久电影中文字幕| 天堂av国产一区二区熟女人妻| 欧美日韩瑟瑟在线播放| 两个人的视频大全免费| 男女下面进入的视频免费午夜| 国产一区二区亚洲精品在线观看| 性插视频无遮挡在线免费观看| 午夜激情欧美在线| 真实男女啪啪啪动态图| 日本熟妇午夜| 老熟妇仑乱视频hdxx| 久久精品国产自在天天线| 亚洲精品粉嫩美女一区| 夜夜爽天天搞| 两人在一起打扑克的视频| 九色国产91popny在线| 国产精品国产高清国产av| 国产精品久久久久久亚洲av鲁大| 久久久成人免费电影| 亚洲性夜色夜夜综合| 一本久久中文字幕| 黄片wwwwww| 777米奇影视久久| 毛片一级片免费看久久久久| a级毛色黄片| 日韩电影二区| 一个人看的www免费观看视频| 久久精品人妻少妇| 一本色道久久久久久精品综合| 亚洲精品久久午夜乱码| 下体分泌物呈黄色| 涩涩av久久男人的天堂| 亚洲最大成人中文| 日韩av在线免费看完整版不卡| 久久久久久伊人网av| 亚洲av电影在线观看一区二区三区| 成人无遮挡网站| 亚洲av不卡在线观看| a级一级毛片免费在线观看| 免费观看av网站的网址| 少妇人妻 视频| 黄色怎么调成土黄色| 亚洲欧美一区二区三区国产| 国产亚洲av片在线观看秒播厂| 久久久久国产网址| 汤姆久久久久久久影院中文字幕| 啦啦啦在线观看免费高清www| 亚洲高清免费不卡视频| 国产69精品久久久久777片| 亚洲人成网站高清观看| 国产黄片美女视频| 街头女战士在线观看网站| 一边亲一边摸免费视频| 国语对白做爰xxxⅹ性视频网站| 久久久久久久国产电影| 麻豆成人av视频| 国产精品国产三级专区第一集| 中文精品一卡2卡3卡4更新| 中文字幕精品免费在线观看视频 | 一级av片app| 久久毛片免费看一区二区三区| 夜夜看夜夜爽夜夜摸| 日韩制服骚丝袜av| 少妇猛男粗大的猛烈进出视频| 欧美3d第一页| 国产伦理片在线播放av一区| 久久精品人妻少妇| 丰满乱子伦码专区| 黄色视频在线播放观看不卡| 男女无遮挡免费网站观看| 欧美精品一区二区大全| 久久久久视频综合| 日本欧美国产在线视频| 又爽又黄a免费视频| 国产精品一区二区在线观看99| 国产一区二区在线观看日韩| 嫩草影院新地址| 热re99久久精品国产66热6| 日韩欧美一区视频在线观看 | 国产精品一区二区三区四区免费观看| 只有这里有精品99| 午夜免费鲁丝| 一级av片app| 伊人久久精品亚洲午夜| 日本av免费视频播放| 18禁在线播放成人免费| 亚洲av综合色区一区| 中文资源天堂在线| 免费人妻精品一区二区三区视频| 国产视频首页在线观看| 人妻少妇偷人精品九色| 亚洲欧美成人精品一区二区| 久久精品久久久久久噜噜老黄| 国产免费视频播放在线视频| 久久久久久久亚洲中文字幕| 我的老师免费观看完整版| 噜噜噜噜噜久久久久久91| 国产亚洲一区二区精品| 精品一区二区三卡| 五月开心婷婷网| 亚洲精品日韩av片在线观看| 各种免费的搞黄视频| 一本久久精品| 熟女av电影| 亚洲美女黄色视频免费看| 日本av手机在线免费观看| 中文字幕亚洲精品专区| 成人免费观看视频高清| xxx大片免费视频| 欧美xxxx性猛交bbbb| 国产伦理片在线播放av一区| 亚洲精品亚洲一区二区| 在线观看国产h片| 嘟嘟电影网在线观看| 高清视频免费观看一区二区| 亚洲真实伦在线观看| 久久97久久精品| 亚洲精品一区蜜桃| 人妻制服诱惑在线中文字幕| 日本av免费视频播放| 久久 成人 亚洲| 欧美高清性xxxxhd video| 欧美区成人在线视频| 五月玫瑰六月丁香| 亚洲精品第二区| 色5月婷婷丁香| av免费观看日本| 嫩草影院入口| 免费少妇av软件| 国产亚洲5aaaaa淫片| 亚洲欧美日韩卡通动漫| 亚洲国产精品国产精品| 大话2 男鬼变身卡| 国产高清有码在线观看视频| 涩涩av久久男人的天堂| 边亲边吃奶的免费视频| 久久久久久久精品精品| 尤物成人国产欧美一区二区三区| 国产高潮美女av| 嫩草影院新地址| 精品视频人人做人人爽| 久久精品国产自在天天线| 亚洲综合精品二区| 欧美高清成人免费视频www| 欧美日韩视频精品一区| 性色avwww在线观看| 成人一区二区视频在线观看| av免费观看日本| 毛片女人毛片| 国产黄片美女视频| 一级a做视频免费观看| 亚洲精品一区蜜桃| 亚洲精品aⅴ在线观看| 一区二区三区免费毛片| 成人一区二区视频在线观看| 国产午夜精品久久久久久一区二区三区| 99久久综合免费| 午夜免费观看性视频| 亚洲欧美日韩无卡精品| 春色校园在线视频观看| 欧美高清成人免费视频www| 久久精品久久久久久噜噜老黄| 亚洲精品国产av成人精品| 成人一区二区视频在线观看| 深夜a级毛片| 国产精品嫩草影院av在线观看| 国产伦在线观看视频一区| 99久久中文字幕三级久久日本| 日韩 亚洲 欧美在线| 精品酒店卫生间| 日本猛色少妇xxxxx猛交久久| 精品亚洲成国产av| 久久99热这里只有精品18| 三级国产精品欧美在线观看| 色网站视频免费| 亚洲国产欧美人成| 妹子高潮喷水视频| 亚洲成色77777| 一级毛片久久久久久久久女| 日韩视频在线欧美| 久久久久久久久大av| 中文字幕亚洲精品专区| av.在线天堂| 好男人视频免费观看在线| 国产成人午夜福利电影在线观看| 少妇人妻久久综合中文| 少妇熟女欧美另类| 一级毛片黄色毛片免费观看视频| 99九九线精品视频在线观看视频| 尤物成人国产欧美一区二区三区| 这个男人来自地球电影免费观看 | 视频中文字幕在线观看| 亚洲不卡免费看| 日韩av不卡免费在线播放| 久久国产亚洲av麻豆专区| 久久婷婷青草| 男女下面进入的视频免费午夜| 日韩亚洲欧美综合| 国产精品99久久久久久久久| 日韩大片免费观看网站| 狂野欧美激情性bbbbbb| 亚洲成人av在线免费| 国产精品一区二区在线不卡| 极品少妇高潮喷水抽搐| 3wmmmm亚洲av在线观看| 老司机影院成人| 伦理电影大哥的女人| 亚洲精品国产成人久久av| 极品少妇高潮喷水抽搐| 街头女战士在线观看网站| 99久久精品国产国产毛片| 少妇高潮的动态图| 男女国产视频网站| 中文在线观看免费www的网站| 国产精品国产三级专区第一集| 免费观看的影片在线观看| 国产在线免费精品| 少妇人妻 视频| 久久久久国产精品人妻一区二区| 成人亚洲欧美一区二区av| 国产精品精品国产色婷婷| 久久国产精品男人的天堂亚洲 | 日本色播在线视频| 51国产日韩欧美| 欧美日韩一区二区视频在线观看视频在线| 成人18禁高潮啪啪吃奶动态图 | 国产高潮美女av| 欧美zozozo另类| 亚洲欧美成人精品一区二区| 亚洲自偷自拍三级| 亚洲美女搞黄在线观看| 国产极品天堂在线| 一本一本综合久久| 女性生殖器流出的白浆| 成人影院久久| 国产成人精品一,二区| 精品国产三级普通话版| 国产视频内射| 中文天堂在线官网| 麻豆成人午夜福利视频| 精品熟女少妇av免费看| 亚洲,一卡二卡三卡| 在线播放无遮挡| 夫妻午夜视频| 久久国产精品男人的天堂亚洲 | www.色视频.com| 国产一级毛片在线| 亚洲美女黄色视频免费看| 精品人妻熟女av久视频| 制服丝袜香蕉在线| 好男人视频免费观看在线| 国产在线视频一区二区| 激情 狠狠 欧美| 成人午夜精彩视频在线观看| 一边亲一边摸免费视频| 亚洲丝袜综合中文字幕| 爱豆传媒免费全集在线观看| 又爽又黄a免费视频| 久久久久精品久久久久真实原创| 极品教师在线视频| 国产精品福利在线免费观看| 午夜激情福利司机影院| 亚洲av日韩在线播放| 国模一区二区三区四区视频| 国产精品一区二区三区四区免费观看| 欧美一区二区亚洲| 赤兔流量卡办理| 国产国拍精品亚洲av在线观看| 亚洲精品一区蜜桃| 深夜a级毛片| 亚洲国产精品专区欧美| 午夜福利影视在线免费观看| 女人久久www免费人成看片| 久久综合国产亚洲精品| 在线观看免费视频网站a站| 一本—道久久a久久精品蜜桃钙片| 纯流量卡能插随身wifi吗| 街头女战士在线观看网站| 欧美xxxx性猛交bbbb| 永久网站在线| 国产伦精品一区二区三区视频9| 亚洲精品第二区| 久久久久久九九精品二区国产| 日本一二三区视频观看| 国产精品一区二区在线观看99| 国产黄片视频在线免费观看| 日日摸夜夜添夜夜添av毛片| 男女边吃奶边做爰视频| 国产免费又黄又爽又色| 色吧在线观看| 尤物成人国产欧美一区二区三区| 男女下面进入的视频免费午夜| 免费看日本二区| 自拍偷自拍亚洲精品老妇| 亚洲精品第二区| 日韩av不卡免费在线播放| 一本一本综合久久| 在线天堂最新版资源| 精品久久久精品久久久| 国产免费福利视频在线观看| 人妻少妇偷人精品九色| 只有这里有精品99| 女的被弄到高潮叫床怎么办| 日本av免费视频播放| 久久精品久久久久久久性| 国产极品天堂在线| 99热全是精品| 欧美亚洲 丝袜 人妻 在线| 一级黄片播放器| 日韩av在线免费看完整版不卡| 精品一区二区三区视频在线| 寂寞人妻少妇视频99o| 大香蕉97超碰在线| 欧美xxⅹ黑人| 色婷婷久久久亚洲欧美| 身体一侧抽搐| 日韩亚洲欧美综合| 我的老师免费观看完整版| 国产淫片久久久久久久久| 久久人妻熟女aⅴ| 久久国内精品自在自线图片| 肉色欧美久久久久久久蜜桃| 夜夜骑夜夜射夜夜干| 亚洲第一区二区三区不卡| 99久久精品国产国产毛片| 日本-黄色视频高清免费观看| 亚洲四区av| 国产精品一区二区性色av| 麻豆成人午夜福利视频| 中文天堂在线官网| 欧美xxxx性猛交bbbb| 午夜福利在线在线| 久久99热6这里只有精品| 黄片无遮挡物在线观看| 国产精品伦人一区二区| 欧美日韩精品成人综合77777| 国产精品蜜桃在线观看| 日韩伦理黄色片| 久久精品国产亚洲av天美| 亚洲综合色惰| 亚洲国产毛片av蜜桃av| 多毛熟女@视频| 国产在线视频一区二区| 久久午夜福利片| 在线看a的网站| 亚洲精品第二区|