• <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

    大香蕉久久成人网| av天堂在线播放| 一级黄色大片毛片| 国产高清videossex| 国产成人av教育| 国产成人影院久久av| 黄色视频不卡| 80岁老熟妇乱子伦牲交| 国产视频一区二区在线看| 亚洲自偷自拍图片 自拍| 国产精品一区二区免费欧美 | 亚洲精品久久久久久婷婷小说| 国产亚洲欧美在线一区二区| 少妇被粗大的猛进出69影院| 精品亚洲成国产av| 男人舔女人的私密视频| 亚洲专区国产一区二区| 亚洲精品国产av成人精品| 色视频在线一区二区三区| 一级,二级,三级黄色视频| 日本黄色日本黄色录像| 久久天堂一区二区三区四区| av欧美777| 精品一区二区三卡| 国产亚洲精品一区二区www | 不卡av一区二区三区| 91大片在线观看| 亚洲国产av新网站| 久久久国产成人免费| 成人av一区二区三区在线看 | 建设人人有责人人尽责人人享有的| 美国免费a级毛片| 精品国内亚洲2022精品成人 | 亚洲精品自拍成人| 日韩,欧美,国产一区二区三区| 久久av网站| 久久精品熟女亚洲av麻豆精品| 少妇精品久久久久久久| 精品人妻一区二区三区麻豆| 久久久久久人人人人人| 欧美日韩黄片免| 亚洲国产精品成人久久小说| 真人做人爱边吃奶动态| 91麻豆av在线| 午夜精品久久久久久毛片777| 亚洲五月婷婷丁香| 免费观看人在逋| 亚洲男人天堂网一区| 国产在线观看jvid| 97在线人人人人妻| 九色亚洲精品在线播放| 欧美国产精品va在线观看不卡| 老司机深夜福利视频在线观看 | 男人爽女人下面视频在线观看| 免费一级毛片在线播放高清视频 | 建设人人有责人人尽责人人享有的| 久久精品亚洲熟妇少妇任你| 精品卡一卡二卡四卡免费| 人人澡人人妻人| 中国国产av一级| 黑人猛操日本美女一级片| 欧美另类亚洲清纯唯美| 日韩有码中文字幕| 在线观看免费日韩欧美大片| 中亚洲国语对白在线视频| a级毛片在线看网站| 亚洲av日韩在线播放| 午夜福利,免费看| av在线老鸭窝| 亚洲国产毛片av蜜桃av| 午夜免费成人在线视频| 在线观看人妻少妇| 国产精品久久久人人做人人爽| 久久国产精品男人的天堂亚洲| 巨乳人妻的诱惑在线观看| 欧美亚洲日本最大视频资源| 少妇 在线观看| 久久久久国产精品人妻一区二区| 亚洲精品国产色婷婷电影| 一级片'在线观看视频| av网站在线播放免费| 亚洲精品国产一区二区精华液| 国产欧美日韩综合在线一区二区| 欧美日韩视频精品一区| 中文欧美无线码| 免费不卡黄色视频| 精品国产一区二区三区久久久樱花| av福利片在线| 黄片小视频在线播放| 99久久人妻综合| 欧美成人午夜精品| 精品国产乱子伦一区二区三区 | 欧美一级毛片孕妇| av又黄又爽大尺度在线免费看| 51午夜福利影视在线观看| 午夜影院在线不卡| 久久久久久久大尺度免费视频| 亚洲成人手机| 免费av中文字幕在线| 亚洲天堂av无毛| 国产亚洲欧美在线一区二区| 啦啦啦免费观看视频1| 国产日韩欧美视频二区| 久久久久久久精品精品| 久久香蕉激情| 巨乳人妻的诱惑在线观看| 丰满少妇做爰视频| 欧美少妇被猛烈插入视频| 国产亚洲精品久久久久5区| 国产在线免费精品| 国产男女超爽视频在线观看| 一级片免费观看大全| 欧美亚洲 丝袜 人妻 在线| 国产有黄有色有爽视频| 宅男免费午夜| 久久av网站| 国产精品国产av在线观看| 国产伦理片在线播放av一区| 美女高潮到喷水免费观看| 天天添夜夜摸| 久久精品成人免费网站| 国产熟女午夜一区二区三区| 91麻豆精品激情在线观看国产 | 自线自在国产av| 国产视频一区二区在线看| 亚洲av男天堂| 国产成人av教育| 老司机午夜福利在线观看视频 | 狠狠婷婷综合久久久久久88av| 久久天堂一区二区三区四区| 久久狼人影院| 亚洲精品美女久久av网站| 国产一区二区三区综合在线观看| 香蕉丝袜av| 国产一区二区三区av在线| 亚洲国产av新网站| 免费不卡黄色视频| 欧美日韩亚洲国产一区二区在线观看 | 五月天丁香电影| 久久久久网色| 久久久久久免费高清国产稀缺| 我要看黄色一级片免费的| 久久综合国产亚洲精品| 欧美另类一区| 一区二区av电影网| bbb黄色大片| 麻豆国产av国片精品| 90打野战视频偷拍视频| 考比视频在线观看| 热re99久久精品国产66热6| 欧美日韩av久久| 少妇 在线观看| 一本综合久久免费| 免费一级毛片在线播放高清视频 | 亚洲精华国产精华精| 国产一卡二卡三卡精品| 飞空精品影院首页| 精品亚洲乱码少妇综合久久| 国产xxxxx性猛交| 国产深夜福利视频在线观看| 狠狠精品人妻久久久久久综合| 亚洲伊人色综图| 日本五十路高清| 久久狼人影院| 在线精品无人区一区二区三| 久久免费观看电影| av福利片在线| 国产区一区二久久| av福利片在线| 操美女的视频在线观看| 久久这里只有精品19| 国产精品成人在线| 午夜福利在线观看吧| 亚洲自偷自拍图片 自拍| 99国产精品99久久久久| 脱女人内裤的视频| 一区二区三区乱码不卡18| www日本在线高清视频| 丝袜脚勾引网站| 成年av动漫网址| 无遮挡黄片免费观看| 久久久久精品人妻al黑| 色精品久久人妻99蜜桃| 热99re8久久精品国产| 女性生殖器流出的白浆| av福利片在线| 窝窝影院91人妻| 12—13女人毛片做爰片一| 亚洲精品美女久久久久99蜜臀| a在线观看视频网站| 久久久精品94久久精品| 手机成人av网站| 久久久精品免费免费高清| 国产福利在线免费观看视频| 色婷婷av一区二区三区视频| 亚洲精品一二三| 亚洲av美国av| 久久精品国产亚洲av高清一级| 69精品国产乱码久久久| 午夜免费鲁丝| av电影中文网址| 中文字幕制服av| 久久精品aⅴ一区二区三区四区| 欧美xxⅹ黑人| 午夜老司机福利片| 99久久综合免费| 老司机靠b影院| 国产欧美日韩一区二区三 | 久久久国产精品麻豆| 超色免费av| 亚洲avbb在线观看| 国产成+人综合+亚洲专区| 韩国高清视频一区二区三区| 久久人人97超碰香蕉20202| 成年动漫av网址| 热99久久久久精品小说推荐| 性高湖久久久久久久久免费观看| 91av网站免费观看| 亚洲av国产av综合av卡| 最新的欧美精品一区二区| 一本色道久久久久久精品综合| 久久午夜综合久久蜜桃| 亚洲国产精品成人久久小说| 女人高潮潮喷娇喘18禁视频| 亚洲视频免费观看视频| 一级片免费观看大全| 欧美精品一区二区免费开放| 精品乱码久久久久久99久播| 午夜福利免费观看在线| 嫁个100分男人电影在线观看| 可以免费在线观看a视频的电影网站| 啦啦啦中文免费视频观看日本| 国产精品久久久久久精品电影小说| 午夜视频精品福利| 美女中出高潮动态图| 一本一本久久a久久精品综合妖精| 黄色视频,在线免费观看| 精品福利永久在线观看| 日韩有码中文字幕| 婷婷丁香在线五月| 高潮久久久久久久久久久不卡| 国产成人精品久久二区二区免费| 欧美日韩av久久| 免费在线观看影片大全网站| 999久久久精品免费观看国产| 在线亚洲精品国产二区图片欧美| 久久人人97超碰香蕉20202| h视频一区二区三区| 国产精品九九99| 男女免费视频国产| 精品国产乱码久久久久久小说| 熟女少妇亚洲综合色aaa.| 亚洲精品中文字幕一二三四区 | 无遮挡黄片免费观看| 成人亚洲精品一区在线观看| 中文欧美无线码| 欧美人与性动交α欧美精品济南到| 色视频在线一区二区三区| 国产精品1区2区在线观看. | 久久精品熟女亚洲av麻豆精品| 成人亚洲精品一区在线观看| 亚洲五月婷婷丁香| 夜夜夜夜夜久久久久| 高潮久久久久久久久久久不卡| 超碰成人久久| 男女午夜视频在线观看| 两个人免费观看高清视频| 黑丝袜美女国产一区| 亚洲一区二区三区欧美精品| 女警被强在线播放| 免费久久久久久久精品成人欧美视频| 欧美久久黑人一区二区| 丰满少妇做爰视频| 老司机靠b影院| 久久精品aⅴ一区二区三区四区| 叶爱在线成人免费视频播放| 国产精品成人在线| √禁漫天堂资源中文www| 欧美日韩国产mv在线观看视频| 九色亚洲精品在线播放| 欧美+亚洲+日韩+国产| 99国产精品99久久久久| 午夜日韩欧美国产| 黑人巨大精品欧美一区二区蜜桃| 日韩精品免费视频一区二区三区| 色94色欧美一区二区| 午夜福利在线观看吧| 精品久久久久久电影网| 国产精品 欧美亚洲| 丰满少妇做爰视频| 不卡一级毛片| 欧美变态另类bdsm刘玥| 欧美日韩亚洲综合一区二区三区_| 十八禁人妻一区二区| 日本猛色少妇xxxxx猛交久久| 亚洲国产欧美网| 欧美日韩黄片免| 欧美午夜高清在线| 精品国产乱码久久久久久男人| 日本av免费视频播放| e午夜精品久久久久久久| av天堂久久9| 夜夜骑夜夜射夜夜干| 成人国产一区最新在线观看| 亚洲精品一二三| 少妇裸体淫交视频免费看高清 | 精品久久蜜臀av无| 91老司机精品| 久久久久网色| 精品人妻一区二区三区麻豆| 深夜精品福利| 大型av网站在线播放| 美女主播在线视频| 久久ye,这里只有精品| 国产亚洲精品一区二区www | 欧美日韩亚洲综合一区二区三区_| 国产黄色免费在线视频| 十分钟在线观看高清视频www| 高潮久久久久久久久久久不卡| a级毛片黄视频| 日韩欧美国产一区二区入口| 亚洲成人国产一区在线观看| 人人妻人人爽人人添夜夜欢视频| 男女下面插进去视频免费观看| 国产成人一区二区三区免费视频网站| 亚洲av电影在线观看一区二区三区| 9191精品国产免费久久| 国产视频一区二区在线看| 久久久国产欧美日韩av| av有码第一页| 国产精品偷伦视频观看了| 国产男女内射视频| 9191精品国产免费久久| 飞空精品影院首页| 80岁老熟妇乱子伦牲交| 肉色欧美久久久久久久蜜桃| 美女主播在线视频| 国产日韩一区二区三区精品不卡| 欧美激情久久久久久爽电影 | 亚洲美女黄色视频免费看| 久久久久视频综合| 丝袜美足系列| 亚洲国产欧美日韩在线播放| 亚洲第一青青草原| 中亚洲国语对白在线视频| 777米奇影视久久| 亚洲中文日韩欧美视频| 亚洲三区欧美一区| 80岁老熟妇乱子伦牲交| 在线看a的网站| 成人国产av品久久久| 久久久久国内视频| 久久久久视频综合| 亚洲国产欧美网| 久久久久久免费高清国产稀缺| 国产主播在线观看一区二区| 狠狠狠狠99中文字幕| 亚洲国产欧美在线一区| 亚洲av成人一区二区三| 亚洲全国av大片| 亚洲黑人精品在线| 91大片在线观看| 一边摸一边抽搐一进一出视频| 精品亚洲成a人片在线观看| 91成年电影在线观看| 欧美精品一区二区免费开放| 成人影院久久| a级毛片黄视频| 国产精品久久久av美女十八| 成人免费观看视频高清| 久9热在线精品视频| 欧美大码av| 国产在视频线精品| 国产亚洲精品第一综合不卡| 亚洲人成77777在线视频| 18禁黄网站禁片午夜丰满| 国产av国产精品国产| 国产成人欧美| 久久av网站| 黄片大片在线免费观看| 国产精品免费大片| 国产精品免费视频内射| 一区二区日韩欧美中文字幕| 最近最新免费中文字幕在线| 91麻豆av在线| 欧美 亚洲 国产 日韩一| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲性夜色夜夜综合| 国产麻豆69| 国产1区2区3区精品| 不卡一级毛片| 国产黄色免费在线视频| 精品少妇内射三级| 另类亚洲欧美激情| 美女大奶头黄色视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲中文日韩欧美视频| 男人舔女人的私密视频| 成人免费观看视频高清| 少妇人妻久久综合中文| 久久国产亚洲av麻豆专区| 中文字幕制服av| 岛国毛片在线播放| 后天国语完整版免费观看| 高潮久久久久久久久久久不卡| 久久久精品区二区三区| 久久毛片免费看一区二区三区| 国产精品 国内视频| 熟女少妇亚洲综合色aaa.| a 毛片基地| 美女午夜性视频免费| 丝袜美腿诱惑在线| 国产一区二区在线观看av| 国产成人免费无遮挡视频| 一区二区三区激情视频| 人妻一区二区av| 中文字幕人妻熟女乱码| 老司机影院成人| 十八禁网站免费在线| 免费观看a级毛片全部| 欧美性长视频在线观看| 久久天堂一区二区三区四区| 99精品久久久久人妻精品| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品中文字幕在线视频| 性高湖久久久久久久久免费观看| 美女主播在线视频| 一本大道久久a久久精品| 美女午夜性视频免费| 中文字幕人妻熟女乱码| 日韩免费高清中文字幕av| 精品国产一区二区三区四区第35| 麻豆av在线久日| 国产av国产精品国产| 欧美精品高潮呻吟av久久| 久久久水蜜桃国产精品网| 美女主播在线视频| svipshipincom国产片| 美女视频免费永久观看网站| 久久久国产精品麻豆| 国产男人的电影天堂91| 亚洲精品一卡2卡三卡4卡5卡 | 天堂中文最新版在线下载| 法律面前人人平等表现在哪些方面 | 欧美亚洲日本最大视频资源| 亚洲av成人不卡在线观看播放网 | 一级毛片精品| 亚洲专区中文字幕在线| 黄频高清免费视频| 脱女人内裤的视频| 欧美av亚洲av综合av国产av| 久久久久国内视频| 国产亚洲精品一区二区www | 免费女性裸体啪啪无遮挡网站| 成年人黄色毛片网站| 一级毛片精品| 亚洲欧美色中文字幕在线| 亚洲av欧美aⅴ国产| av天堂久久9| 99久久综合免费| 国产精品二区激情视频| 欧美另类亚洲清纯唯美| 99久久99久久久精品蜜桃| 国产老妇伦熟女老妇高清| 考比视频在线观看| 日本撒尿小便嘘嘘汇集6| 激情视频va一区二区三区| 女人爽到高潮嗷嗷叫在线视频| 国产精品久久久av美女十八| 亚洲黑人精品在线| 中国国产av一级| 777久久人妻少妇嫩草av网站| www.av在线官网国产| 十分钟在线观看高清视频www| 岛国在线观看网站| a级片在线免费高清观看视频| 亚洲精品自拍成人| 脱女人内裤的视频| 久久精品aⅴ一区二区三区四区| 国产一区二区在线观看av| 国产精品一二三区在线看| 免费女性裸体啪啪无遮挡网站| 高清av免费在线| 嫩草影视91久久| 黑人猛操日本美女一级片| 水蜜桃什么品种好| 男女床上黄色一级片免费看| 一区二区av电影网| 999久久久精品免费观看国产| 国产成人精品久久二区二区91| 91老司机精品| 男人舔女人的私密视频| 99国产精品一区二区蜜桃av | av不卡在线播放| 国产xxxxx性猛交| 国产熟女午夜一区二区三区| 久久国产亚洲av麻豆专区| 日韩欧美国产一区二区入口| 国产成人精品无人区| 欧美大码av| 少妇精品久久久久久久| 久久精品久久久久久噜噜老黄| 啦啦啦在线免费观看视频4| av在线老鸭窝| 国产精品久久久久久人妻精品电影 | 国产一区有黄有色的免费视频| 麻豆av在线久日| 国产成人免费无遮挡视频| 国产精品麻豆人妻色哟哟久久| 日韩大码丰满熟妇| 搡老熟女国产l中国老女人| 两个人免费观看高清视频| 大香蕉久久成人网| 老司机在亚洲福利影院| 亚洲精品美女久久久久99蜜臀| 狂野欧美激情性xxxx| 韩国精品一区二区三区| 男女边摸边吃奶| 国精品久久久久久国模美| 国产淫语在线视频| 麻豆国产av国片精品| 久久青草综合色| 无遮挡黄片免费观看| 91国产中文字幕| 国产精品99久久99久久久不卡| 91成年电影在线观看| 91大片在线观看| 久久精品国产亚洲av香蕉五月 | av免费在线观看网站| 日韩,欧美,国产一区二区三区| 熟女少妇亚洲综合色aaa.| 涩涩av久久男人的天堂| 99久久精品国产亚洲精品| 国产精品1区2区在线观看. | 男女下面插进去视频免费观看| 啦啦啦视频在线资源免费观看| 国产一区有黄有色的免费视频| 亚洲人成电影免费在线| 美女视频免费永久观看网站| 美女主播在线视频| 9热在线视频观看99| 天天影视国产精品| 日本黄色日本黄色录像| 亚洲美女黄色视频免费看| 97人妻天天添夜夜摸| 亚洲九九香蕉| 亚洲成人手机| 一二三四社区在线视频社区8| 婷婷丁香在线五月| av线在线观看网站| 国产av又大| 97精品久久久久久久久久精品| 最近最新中文字幕大全免费视频| 午夜91福利影院| 国产精品九九99| 两个人看的免费小视频| 久久精品国产a三级三级三级| a在线观看视频网站| 国产精品1区2区在线观看. | 国产成人一区二区三区免费视频网站| 999久久久国产精品视频| 精品熟女少妇八av免费久了| 又紧又爽又黄一区二区| 欧美xxⅹ黑人| 久久精品亚洲熟妇少妇任你| 91成年电影在线观看| 91精品国产国语对白视频| 18禁国产床啪视频网站| 免费av中文字幕在线| 一区福利在线观看| 一本—道久久a久久精品蜜桃钙片| 欧美日韩亚洲高清精品| 大香蕉久久成人网| 秋霞在线观看毛片| 免费看十八禁软件| 一级a爱视频在线免费观看| 免费女性裸体啪啪无遮挡网站| 国产在视频线精品| 国产有黄有色有爽视频| 精品卡一卡二卡四卡免费| 美女高潮到喷水免费观看| 丁香六月欧美| 久久久水蜜桃国产精品网| 人妻人人澡人人爽人人| 国产福利在线免费观看视频| 久久精品亚洲熟妇少妇任你| 欧美精品一区二区大全| 日韩 亚洲 欧美在线| 亚洲激情五月婷婷啪啪| 免费在线观看影片大全网站| 99久久国产精品久久久| 欧美精品高潮呻吟av久久| 18禁裸乳无遮挡动漫免费视频| 国产真人三级小视频在线观看| 亚洲av电影在线观看一区二区三区| 国产主播在线观看一区二区| 午夜久久久在线观看| 国产又色又爽无遮挡免| 窝窝影院91人妻| 精品国产一区二区三区久久久樱花| 在线 av 中文字幕| 最近最新免费中文字幕在线| 下体分泌物呈黄色| 1024视频免费在线观看| 侵犯人妻中文字幕一二三四区| 伊人亚洲综合成人网| 五月天丁香电影| 亚洲精华国产精华精| 啦啦啦啦在线视频资源| 亚洲国产av新网站| 国产精品麻豆人妻色哟哟久久| 久久久久久亚洲精品国产蜜桃av| 久久精品国产a三级三级三级| 亚洲伊人久久精品综合| 国产在线一区二区三区精|