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

    Role of Extratropical Cyclones in the Recently Observed Increase in Poleward Moisture Transport into the Arctic Ocean

    2018-12-06 07:35:52GianVILLAMILOTEROJingZHANGJuanxiongHEandXiangdongZHANGDepartmentofPhysicsNorthCarolinaStateUniversityGreensboroNorthCarolina27401USA
    Advances in Atmospheric Sciences 2018年1期

    Gian A.VILLAMIL-OTERO,Jing ZHANG,Juanxiong HE,and Xiangdong ZHANGDepartment of Physics,North Carolina A&T State University,Greensboro,North Carolina 27401,USA

    2Department of Energy and Environmental Systems,North Carolina A&T State University,Greensboro,North Carolina 27401,USA

    3Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029,China

    4Department of Atmospheric Science,University of Alaska Fairbanks,Fairbanks,Alaska 99775,USA

    Role of Extratropical Cyclones in the Recently Observed Increase in Poleward Moisture Transport into the Arctic Ocean

    Gian A.VILLAMIL-OTERO?1,Jing ZHANG1,2,Juanxiong HE3,and Xiangdong ZHANG41Department of Physics,North Carolina A&T State University,Greensboro,North Carolina 27401,USA

    2Department of Energy and Environmental Systems,North Carolina A&T State University,Greensboro,North Carolina 27401,USA

    3Institute of Atmospheric Physics,Chinese Academy of Sciences,Beijing 100029,China

    4Department of Atmospheric Science,University of Alaska Fairbanks,Fairbanks,Alaska 99775,USA

    Poleward atmospheric moisture transport(AMT)into the Arctic Ocean can change atmospheric moisture or water vapor content and cause cloud formation and redistribution,which may change downward long wave radiation and,in turn,surface energy budgets,air temperatures,and sea-ice production and melt.In this study,we found a consistently enhanced poleward AMT across 60°N since 1959 based on the NCAR–NCEP reanalysis.Regional analysis demonstrates that the poleward AMT predominantly occurs over the North Atlantic and North Pacific regions,contributing about 57%and 32%,respectively,to the total transport.To improve our understanding of the driving force for this enhanced poleward AMT,we explored the role that extratropical cyclone activity may play.Climatologically,about 207 extratropical cyclones move across 60°N into the Arctic Ocean each year,among which about 66(32%of the total)and 47(23%)originate from the North Atlantic and North Pacific Ocean,respectively.When analyzing the linear trends of the time series constructed by using a 20-year running window,we found a positive correlation of 0.70 between poleward yearly AMT and the integrated cyclone activity index(measurement of cyclone intensity,number,and duration).This shows the consistent multidecadal changes between these two parameters and may suggest cyclone activity plays a driving role in the enhanced poleward AMT.Furthermore,a composite analysis indicates that intensification and poleward extension of the Icelandic low and accompanying strengthened cyclone activity play an important role in enhancing poleward AMT over the North Atlantic region.

    atmospheric moisture transport,cyclone activity,atmospheric circulation,Arctic,climate change

    1.Introduction

    Alteration of poleward atmospheric moisture transport(AMT)into the Arctic may have significant impacts on Arctic climate changes since moisture,or water vapor,represents one of the major greenhouse gases and critically determines cloud formation and distribution.Increased moisture transport into the Arctic Ocean may,therefore,enhance water vapor-or cloud-forced downward longwave radiation,contributing to an increase in surface air temperature and summer sea-ice melt and a decrease in winter sea-ice production.Wang and Key(2005)found that increased water vapor content appears to be influential in the enhancement of downward longwave radiation in the Arctic,and Zuidema et al.(2005)further emphasized the role of liquid-water-containing clouds.Francis and Hunter(2006)suggested that the location of the summer ice edge is strongly correlated to the variability in longwave energy emitted by the atmosphere.Kay et al.(2008)found that changes in clouds are an increasingly important control in sea-ice extent minima.In particular,when studying the extreme event of Arctic sea-ice loss in 2007,Kay and Gettelman(2009)found that,while summer cloud reductions enhanced ice-albedo feedback,an increase in cloudiness in fall over newly open water helps trap heat by increasing the downwelling longwave radiation but also reduces net absorbed shortwave radiation.In addition,climate model projections have shown that an increase in poleward AMT caused by a warming planet likely creates a positive feedback on the Arctic system by increasing latent heat release and the emission of longwave radiation to the surface(Skific et al.,2009).

    The impact of intense poleward AMT across 70°N during boreal winter on the Arctic-mean downward longwave radi-ation was recently quantified by Woods et al.(2013).They found that poleward moisture intrusions contribute to about 40%of the interannual variance in downward longwave radiation and about 30%of the interannual variance in skin temperature.Further evidence shows that the increase in downward longwave radiation associated with these moisture intrusions can impact the loss of sea ice(Woods and Caballero,2016)and result in rapid surface air temperature increases(Kim et al.,2017).Liu and Barnes(2015)investigated these extreme moisture intrusions and found that they are linked to Rossby wave breaking.These planetary-scale waves are capable of diverting moisture,sensible,and latent heat fluxes into the Arctic,enhancing the downward infrared radiation and surface warming(Baggett et al.,2016;Baggett and Lee,2017;Gong et al.,2017).These moisture intrusions also resemble atmospheric rivers(ARs;Newell et al.,1992;Ralph and Dettinger,2011;Newman et al.,2012;Liu and Barnes,2015;Baggett et al.,2016).ARs have been found to account for the majority of poleward moisture transport in the midlatitudes(Zhu and Newell,1998;Guan and Waliser,2015).Also,Dacre et al.(2015)attributed the moisture in ARs to local evaporation and convergence within the warm sectors of cyclones,instead of long-range transport as the source.

    Synoptic cyclones constitute the major part of atmospheric transient eddies,which may also play an important role in moisture transport,in addition to the stationary and planetary waves(e.g.,Peixoto and Oort,1992;Baggett et al.,2016;Vihma et al.,2016;Baggett and Lee,2017).Groves and Francis(2002)found that the movement of transient features into the Arctic account for 74%of the net precipitation in the Arctic basin.On the other hand,studies have suggested that cyclone activity plays a key role in the interaction between the surface and the atmosphere and has strong impacts on the Arctic sea ice(Murray and Simmonds,1995;Sorteberg and Kvingedal,2006).The strength,rather than the number,of cyclones in the Arctic basin,plays a more critical role in the decrease in September ice extent(Simmonds and Keay,2009).The impacts of cyclones on the surface and sea ice may be affected through changes in atmospheric moisture content and resulting longwave radiation.In particular,it has been found that cyclone tracks have shifted poleward and cyclone activity has intensified in the Arctic(Zhang et al.,2004;Simmonds et al.,2008;Sepp and Jaagus,2011;Bender et al.,2012),and the poleward shift may continue in the future as projected by climate models under global warming scenarios(Yin,2005).So,to improve our understanding of the changes in Arctic atmospheric moisture and the water cycle,we decided to quantitatively investigate the linkages between poleward AMT and extratropical cyclone activity.In this study,we used 60°N as the boundary between the midlatitudes and the Arctic to investigate the impact of extratropical cyclones that move into the Arctic on the AMT.The 60°N line stands as the division between the“polar cap”and the midlatitudes,as it generally corresponds to the zero contour between the polar and midlatitude centers of action of the Arctic Oscillation(AO)(Thompson and Wallace,1998).

    The goal of this study is to understand the role that extratropical cyclone activity plays in the poleward AMT.We seek to achieve this goal by answering the following specific questions:First,do AMT and cyclone activity exhibit any similar regional and interannual variations?Second,which components(e.g.,intensity,number,etc.)of cyclone activity correlate closely with the changes in poleward AMT?Andfinally,how does cyclone activity behave differently under the positive and negative phases of poleward AMT,and how does this difference impact the AMT?The paper is organized as follows:Section 2 describes the data and methods used to obtain the AMT and cyclone activity data.Sections 3 and 4 describe the results of the time series analysis of AMT and cyclone activity and the analysis of the relationship between them.Section 5 concludes the paper with a summary and discussion.

    2.Data and methodology

    Following Zhang et al.(2012),AMT is calculated with the original six-hourly spectral output of the NCAR–NCEP reanalysis(Kalnay et al.,1996)from January 1948 to December 2016,using the formula where psis the surface pressure,q the specific humidity,v the meridional component of the wind vector,g the gravity acceleration,ρ the density of water,and σ is the terrain-following coordinate of the data.is the longitudinal integral over the region,andis the time integral for the selected time range.

    In this study,only the AMT across the latitude of 60°N is used to measure the amount of moisture transported meridionally,in which positive numbers mean poleward transport of moisture,and negative numbers mean equatorial transport.Monthly and yearly integrations of AMT are calculated to facilitate understanding how the AMT behaves through the entire study period globally,as well as in four regions comprising North America(120°–70°W),the North Atlantic(70°W–20°E),Eurasia(20°–140°E),and the North Pacific(140°E–120°W).

    The cyclone activity parameters used for this study are calculated with sea level pressure(SLP)from the same sixhourly NCEP–NCAR reanalysis for the period from January 1948 to December 2016.Then,the cyclone activity parameters between 55°N and 65°N are used to analyze and compare the relationship between cyclone activity and the AMT crossing 60°N.Selecting this area spanning 10°meridionally allows us to track,and therefore analyze,only those cyclones that developed below 60°N and crossed into the Arctic.Cyclones are identified and tracked using empirical criteria developed by Zhang et al.(2004)after examining six-hourly SLP reanalysis snapshots.The criteria used for the identification and tracking of cyclones from Zhang et al.(2004)include:

    (1)If the SLP at one grid point is lower than the surround-ing eight grid points,a cyclone candidate is identified with its center at that grid point.

    (2)The minimum pressure gradient calculated between the candidate cyclone center and the adjacent eight grid points is required to be at least 0.15 hPa(100 km)-1.

    (3)The minimum SLP gradient between the cyclone candidate center and its four adjacent grid points must be negative inward for effectively removing any open lows from the data.

    (4)If two or more cyclone candidates are detected within a radius of 1200 km of each other at the same time,they are considered to be a single cyclone.

    (5)The cyclone candidate is required to have a minimum lifetime of 12 h.In addition,if the location of a cyclone is within 600 km of a cyclone center identified 6 h earlier,the center is considered to be a new location of the existing cyclone.Otherwise,a new cyclone is generated.

    To describe the variation characteristics for cyclone activity,we use the following three cyclone activity parameters,as defined by Zhang et al.(2004):

    (1)Cyclone intensity:defined by averaging the differences between the central SLP of cyclones and the climatological monthly mean SLP at corresponding grid points for all cyclones in a particular region during the month.

    (2)Cyclone trajectory count(number):defined as the number of cyclones over the entire duration occurring in a particular region during the month.

    (3)Cyclone Activity Index(CAI):defined as the sum of the differences between the cyclone central SLP and the climatological monthly mean SLP at corresponding grid points over all cyclone centers in a particular region during the month.Thus,CAI is an integrated parameter measuring the intensity,number and duration of cyclone activity.

    3.Climatology of AMT and cyclone activity across 60°N

    Fig.1.Annual variation of yearly AMT(units:km3)across 60°N(a)globally and(b–e)for the four regions of(b)the North Atlantic,(c)the North Pacific,(d)Eurasia,and(e)North America.Positive values indicate a poleward flux of moisture and negative values an equatorward flux.The thin solid straight lines are trend lines that indicate upward/downward trends for each of the three periods determined through a change-point analysis.

    Performing a change-point analysis on the values of AMT that have been integrated for each year along 60°N allows us to determine whether any changes in the behavior of the data are present(Taylor,2015).With this analysis,we identify three distinct periods at a 99%significance level by the Student’s t-test(Fig.1).The first period is denoted by a downward trend of-26.95 km3yr-1of poleward yearly AMT in the first 11 years.The second period in the yearly AMT trend is identified from 1959 to 1982 where an upward trend of 6.08 km3yr-1is present.The third and last period is defined to span from 1982 up to the end of the study period and continues the upward trend seen in the second period,albeit at a smaller magnitude of 4.46 km3yr-1.The year 1982 marks the year in which the average of yearly AMT after 1982 is higher than the average of yearly AMT across the entire period of study,whereas the average of yearly AMT before 1982 is lower than the average yearly AMT for the period of study.Poleward yearly AMT is the largest in the North Atlantic region,followed by the North Pacific and Eurasia regions.Yearly AMT in the North America region is mainly equatorward,which was also found by Cullather et al.(2000)and Rogers et al.(2001).Annually,the yearly AMT across 60°N is 5464 km3,in which 89%is from the oceanic regions(57%from North Atlantic and 32%from North Pacific)and 11%is from the continental regions.The Eurasia region accounts for 14%of the total poleward yearly AMT,balanced by a 3%equatorward transport in the North America region.

    Using the same three periods found in the global yearly AMT(Fig.1a)for each region(Figs.1b–e)helps to explain the region’s contributions to the trends of global yearly AMT observed in each period.During thefi rst period,all regions except the North Atlantic present downward trends,with the North Pacific contributing to 70%of the global downward trend.For the second period,the contributions to the global upward trend are mainly provided by the North Atlantic and Eurasia.Finally,all regions except Eurasia have upward trends in the last period.Eurasia exhibiting a downward trend during this last period explains why the global upward trend is smaller than that of the second period.

    Measuring the number of cyclones across 55°–65°N also demonstrates a clear regional distribution(Table 1).Annually,an average of 207 cyclones cross the latitudes of 55°–65°N,of which 66(32%)cyclones are from the North Atlantic region and 57(27%)from Eurasia.North America presents the smallest number with 37(18%)cyclones,and the remaining 47(23%)are from the North Pacific.The total cyclone number across 55°–65°N demonstrates an upward trend before 1971 and a downward trend after 1971,as indicated by the change-point analysis(Fig.2).Similar to Zhang et al.(2004),the cyclone number in the North Pacific region exhibits a downward trend throughout both periods,which helps promote Eurasia to become a bigger contributor than the North Pacific(Table 1).In addition,Eurasia also presents the largest upward trend in cyclone number amongst all the regions for the first period.

    Fig.2.Global yearly cyclone number(top),cyclone intensity(units:hPa)(middle),and CAI(units:hPa)(bottom)that crossed 60°N.The thin solid straight lines are trend lines that indicate upward/downward trends for each period determined through a change-point analysis.

    Globally,the intensity of extratropical cyclones crossing 55°–65°N is around 17 hPa stronger than the mean SLP in that area(Table 1).Again,similar to Zhang et al.(2004),cyclone intensity exhibits an upward trend during the entire study period,though the change-point analysis indicatesa slight change in the upward trend behavior in 1955(Fig.2).Of the various regions,we find that the North Atlantic presents the highest intensity,and the trends in each period are the closest to the global trends,which again illustrates the importance of this region to global cyclone activity(Table 1).In contrast to the trend in cyclone number,cyclone intensity in the North Pacific shows an upward trend for both periods.

    Table 1.Cyclone number,intensity and CAI values globally and for the four regions.The percentage of the global total is provided in parentheses.For each region,trend1 and trend2 represent the trends for each period as determined from the change-point analysis and are provided in square brackets.

    Variations in cyclone number and intensity can have opposite phases;therefore,an integrated index(i.e.,the CAI)is also used here to evaluate the overall variation in cyclone activity across 55°–65°N(Fig.2).The change-point analysis is also performed for the CAI and results in a downward trend before 1956 and an upward trend after 1956.Globally,the CAI has an average of 2573 hPa across 55°–65°N,from which the most significant contributions are provided by the North Atlantic and Eurasia regions with 922 hPa(36%)and 592 hPa(23%)respectively(Table 1).The region that contributes the least to the global CAI is North America,with 512 hPa(20%).The North Pacific provides the remaining amount of 547 mb(21%),in accordance with the lower values of cyclone number and intensity when compared to the other ocean basin.Because the biggest contribution of the CAI comes from the North Atlantic region,the changes in the North Atlantic CAI throughout the study period are reflected in the global CAI.The CAI behavior during the downward trend period can be attributed to the North Atlantic and North America regions,while the upward trend period arises from the North Atlantic,Eurasia,and North America regions(Table 1).The CAI trends in the North Pacific and North America mostly cancel each other out,especially during the global downward trend period.

    4.Relationship between cyclone activity and enhanced poleward AMT across 60°N

    With the observed increases in both AMT and cyclone activity across 60°N and the similarity of their regional distributions,we further explore the relationship between them in this section.To understand if the trends that have been observed in AMT,as well as cyclone activity,have any relationship between one another,we calculate the trends within the distinct time series of AMT,cyclone number,cyclone intensity,and CAI.Because the observed trends are inherently long-term,a 20-yr running window is used to calculate these trends.Figures 3a–c illustrate how these trend time series behave over time.The time series of the yearly AMT trend exhibits a negative trend during 1958–67,followed by a positive trend during 1967–99.During the positive trend period,a maximum is reached around 1982,the year that the change-point analysis denotes as the beginning of the third period in the yearly AMT.The trends in cyclone number,intensity,and CAI behave differently.For example,cyclone number has a positive trend throughout the majority of the time series before 1999,while cyclone intensity begins with a relatively flat trend,then becomes negative during 1961–65,and remains positive till 1995.When the trend correlations of yearly AMT to cyclone number,intensity,and CAI are calculated,we have-0.01,0.46 and 0.70,respectively,of which only the latter two are significant according to the Student’s t-test at the greater than 99%confidence level.The correlation results reveal that when the relations between cyclone activity and yearly AMT are explored,one cannot focus on simply the amount of storms or the intensity.Instead,an index such as the CAI that encompasses most major characteristics of cyclone activity including quantity,intensity,and duration is required.The trends in CAI display a similar pattern as the trends in yearly AMT,except at the beginning of the examined period when the CAI trend is relatively flat,and during 1995–2005 when the CAI trend is negative.Therefore,the overall cyclone activity,not intensity or number only,seems to be closely related to the change in yearly AMT.

    Fig.3.Comparison of trend time series for AMT and(a)cyclone number,(b)cyclone intensity,and(c)CAI,from 1948 to 2016,using a 20-yr running window.Correlations between each variable set are shown in the plot.

    Cyclone activity,particularly in the North Atlantic,is closely tied to the AO or North Atlantic Oscillation(NAO)index.For example,a positive AO/NAO index is related to an increase in storminess in the vicinity of Iceland and the Norwegian Sea,provided by a higher amount of storms in the North Atlantic(Zhang et al.,2004;Hurrell,2015).The variability of AMT is also well related to the AO index,mainly due to changes in the mean flow and storminess changing the moisture transport pattern(Dickson et al.,2000;Jakobson and Vihma,2010;Hurrell,2015).However,in recent decades,the linkage between AMT and AO,as well as the cyclone activity,have become complicated,which might be attributable to the climate pattern change from a conventional tripolar AO to an anomalous dipolar leading pattern,also known as the Arctic Rapid-change Pattern(ARP)discussed by Zhang et al.(2008).

    To better understand how the atmospheric circulation pattern and associated cyclone activity impact the poleward AMT,we perform a composite analysis of SLP associated with positive and negative AMT phases.To construct the composites of SLP,plus or minus one standard deviation of monthly AMT for the month analyzed is used to group the years into when strong or weak poleward AMT is present.Figure 4 depicts the interannual variations of monthly AMT in each month.Here,we only focus on the time when the increased poleward AMT occurs,so the period with upward trends,i.e.,the second and third periods(1959–2016)in the yearly AMT change-point analysis,are used to calculate the trend of monthly AMT in each month.The monthly AMT varies from month to month with greater monthlyAMT(~400–500 km3)from August to November,and weaker monthly AMT(up to~350 km3)from February to June.Of special interest is that from 1956 through 2016,all months show upward trends with the highest trends occurring in August(1.824 km3yr-1),September(1.297 km3yr-1),and November(1.015 km3yr-1).This suggests that any changes that affect poleward monthly AMT might be more visible in these months.Also of interest is that October,the month with the second highest mean monthly AMT(497.74 km3),is also one of the months with a relatively flat trend(0.495 km3yr-1).The presence of this flat trend indicates that a relatively stable poleward monthly AMT is present in October.

    Fig.4.Monthly AMT for(a–l)January–December,respectively,with one standard deviation above and below the mean(dotted)illustrated by dashed lines.The trend line from 1959 to 2016 is shown by a bold line.

    The composites of SLP and monthly AMT across 60°N associated with the positive and negative AMT phases for January,April,July,and October,representative of each season,are shown in Fig.5.During the winter month of January,we observe significant planetary-scale systems including the Siberian high,North America high,Aleutian low,and Icelandic low along the belt of 60°N(Figs.5a and b).The differences in the SLP composites between the positive and negative phases(significant at the 99%confidence level based on the Student’s t-test)indicate that the Icelandic low is much stronger and extends poleward into the Arctic during the positive phase.Cyclone activity crossing 60°N over the North Atlantic region,as represented by the CAI in Table 2,is also stronger during the positive phase.Enhanced poleward AMT in the North Atlantic region co-occurs with stronger cyclone activity there(Fig.5a and Table 2).

    During the spring month of April,similar differences in the SLP composite to that seen in winter are present.Again,the Icelandic low(though much weaker when compared to winter)extends poleward during the positive phase(Figs.5d and e).Similarly,there is stronger cyclone activity crossing 60°N through the North Atlantic region co-occurring with enhanced poleward AMT in the region(Table 2).

    Over summer(July),the differences in the SLP composite between the positive and negative phases are much smaller(Figs.5g–i).The low pressure over the Arctic is slightly stronger during the positive phase.Also,the low pressure over Eurasia is slightly stronger during the negative phase.Cyclone activity crossing 60°N in the North Atlantic region is much weaker in July but still relatively strong during the positive phase(Table 2).

    During fall(October),a pattern of differences in SLP composites reminiscent of that in spring and winter is again observed(Figs.5j and k).The Icelandic low does not extend as far south as it did during the positive phase.Cyclone activity crossing 60°N measured by the CAI(Table 2)is similar to that in April,even though the entire low system is much stronger in October.During the negative phase,the Icelandic low is centered around 60°N,suggesting significant cycloneactivity across there,and a relatively large CAI is observed(Table 2).

    Table 2.AMT(units:km3)and CAI(units:hPa)composite values for AMT larger than 1 standard deviation(positive phase)and smaller than 1 standard deviation(negative phase)for the North Atlantic region for January,April,July and October.

    Through the composite analysis performed here,it is possible to examine the relationship between the circulation pattern and associated cyclone activity and the poleward monthly AMT.Focus here is given to the North Atlantic region where the contribution to global monthly AMT is the greatest and the differences of SLP composites are significant at the 99%confidence level.A consistent relationship is found insofar as stronger cyclone activity across 60°N measured by the CAI there generally co-occurs with enhanced poleward monthly AMT in each representative seasonal month examined.This further supports why a strong correlation exists between the trends of the CAI and poleward yearly AMT.

    5.Summary and discussion

    The poleward AMT calculated from the NCAR–NCEP reanalysis data set has been analyzed along 60°N as a general boundary between the midlatitudes and the Arctic.An upward trend of poleward yearly AMT across 60°N since 1959 has been found,which has then consistently enhanced since 1982.The poleward yearly AMT demonstrates a distinct regional distribution,in which the maximum poleward yearly AMT occurs in the North Atlantic region.Poleward AMT also shows apparent seasonal variations,with greater values during August to November and smaller ones from February to June.An upward trend of monthly AMT is identified in all months since 1959,with the highest trend occurring in August,September,and November.It is also shown that greater poleward monthly AMT,as well as a greater upward trend,generally occurs during late summer and early fall,suggesting that the trend identified previously has continued and even amplified along with the rapidly changing Arctic system(Serreze et al.,1995;Comiso et al.,2008;Zhang et al.,2008;Skific et al.,2009).

    A poleward shift in cyclone tracks has been identified by several studies(Zhang et al.,2004;Simmonds et al.,2008;Sepp and Jaagus,2011;Bender et al.,2012),which provides an incentive to study the relationship between extratropical cyclone activity and enhanced poleward AMT.The activity of cyclones across 55°–65°N measured with cyclone number,intensity and the integrated CAI is analyzed.Similar to the upward trend in poleward yearly AMT,cyclone intensity and the CAI also demonstrate a clear upward trend since 1956,with the biggest contribution to the CAI trend coming from the North Atlantic and Eurasia regions.

    The similarity of the trends observed in both the poleward yearly AMT and the activity of cyclones across 60°N motivated us to further explore the relationships between them.From the time series of linear trends constructed by using a 20-yr running window,wefi nd a positive correlation of 0.70 between the poleward yearly AMT and CAI,an integrated measurement of cyclone activity including intensity,number,and duration.These consistent multidecadal variations be-tween yearly AMT and the CAI may suggest a driving role played by the overall cyclone activity in the poleward AMT,considering the fundamental nature of cyclones in holding and transporting moisture and energy.

    Fig.5.Composite SLP for AMT larger than 1 standard deviation(left column)and smaller than 1 standard deviation(middle column)and their difference[right column;only statistically significant(99%confidence level)grid points plotted]for(a–c)January,(d–f)April,(g–i)July,and(j–l)October.Poleward(equatorward)moisture transport at 60°N is represented by red(blue)dots,whose sizes represent the magnitude of the transport.

    Furthermore,a composite analysis of SLP and cyclone activity measured by the CAI and based on the positive and negative AMT phases reveals the impact of changed atmospheric circulation and cyclone activity on the poleward AMT.It is found that the poleward intrusions of the Icelandic low and accompanying enhancement of cyclone activity across 60°N would play an important role in driving the enhanced poleward AMT.The Icelandic low is generally defined on a monthly or seasonal time scale,as an integral representation of daily-scale storm activities.So,enhanced storm activity would consequently intensify the Icelandic low,perhaps through the wave breaking process.The poleward intrusion is consistent with the changes in the largescale atmospheric circulation,which are characterized by a spatial structure shift from the AO to ARP(Zhang et al.,2008,2012),and have decisively driven the rapid changes in the Arctic climate.

    Poleward AMT impacts the Arctic climate and environment in many different ways.Since water vapor content is a major greenhouse gas,an increase in poleward AMT can thus divert moisture in the Arctic,enhancing the downward infrared radiation and,in turn,causing surface warming(e.g.,Woods et al.,2013;Baggett et al.,2016;Baggett and Lee,2017;Gong et al.,2017).In addition,a greater amount of water vapor content may lead to greater amounts of condensation,which in turn warms the atmosphere even more due to latent heat release.This may generate a positive feedback between latent heat release and the emission of atmospheric longwave radiation to the surface,further enhancing the warming and sea-ice decrease,and leading to enlarged open water and increased local evaporation.Details and quantitative analysis of this feedback process associated with cyclone activity and poleward AMT would be an appropriate avenue for follow-up research.

    Acknowledgements.The authors thank the editor,Dr.Tido SEMMLER,and the anonymous reviewers,for providing insightful suggestions and comments that have greatly improved the manuscript.This work was supported by the NSF (Grant Nos.ARC-1023592,ARC-1107509,and PLR-1304684).The computing resources were provided by the Arctic Region Supercomputing Center at the University of Alaska Fairbanks.

    Baggett,C.,and S.Lee,2017:An identification of the mechanisms that lead to arctic warming during planetary-scale and synoptic-scale wave life cycles.J.Atmos.Sci.,74,1859–1877,http://dx.doi.org/10.1175/JAS-D-16-0156.1.

    Baggett,C.,S.Lee,and S.Feldstein,2016:An investigation of the presence of atmospheric rivers over the North Pacific during planetary-scale wave life cycles and their role in Arctic warming.J.Atmos.Sci.,73,4329–4347,http://dx.doi.org/10.1175/JAS-D-16-0033.1.

    Bender,F.A.-M.,V.Ramananthan,and G.Tselioudis,2012:Changes in extratropical storm track cloudiness 1983–2008:Observational support for a poleward shift.Climate Dyn.,38,2037–2053,https://doi.org/10.1007/s00382-011-1065-6.

    Comiso,J.C.,C.L.Parkinson,R.Gersten,and L.Stock,2008:Accelerated decline in the Arctic sea ice cover.Geophys.Res.Lett.,35,L01703,https://doi.org/10.1029/2007GL031972.

    Cullather,R.I.,D.H.Bromwich,and M.C.Serreze,2000:The atmospherichydrologic cycle over the Arctic basin from reanalyses.Part I:Comparison with observations and previous studies.J.Climate,13,923–937,http://dx.doi.org/10.1175/1520-0442(2000)013<0923:TAHCOT>2.0.CO;2.

    Dacre,H.F.,P.A.Clark,O.Martinez-Alvarado,M.A.Stringer,and D.A.Lavers,2015:How do atmospheric rivers form?Bull.Amer.Meteor.Soc.,96,1243–1255,http://dx.doi.org/10.1175/BAMS-D-14-00031.1.

    Dickson,R.R.,and Coauthors,2000:The Arctic Ocean response to the North Atlantic oscillation.J.Climate,13,2671–2696,http://dx.doi.org/10.1175/1520-0442(2000)013<2671:TAORTT>2.0.CO;2.

    Francis,J.A.,and E.Hunter,2006:New insight into the disappearing arctic sea ice.EOS,87,509–511,http://dx.doi.org/10.1029/2006EO460001.

    Gong,T.T.,S.Feldstein,and S.Lee,2017:The role of downward infrared radiation in the recent arctic winter warming trend.J.Climate,30,4937–4949,http://dx.doi.org/10.1175/JCLID-16-0180.1.

    Groves,D.G.,and J.A.Francis,2002:Moisture budget of the arctic atmosphere from TOVS satellite data.J.Geophys.Res.,107,ACL 11-1–ACL 11-21,http://dx.doi.org/10.1029/2001JD001191.

    Guan,B.,and D.E.Waliser,2015:Detection of atmospheric rivers:Evaluation and application of an algorithm for global studies.J.Geophys.Res.:Atmos.,120,12514–12535,http://dx.doi.org/10.1002/2015JD024257.

    Hurrell,J.W.,2015:Climate variability:North Atlantic and Arctic oscillation.Encyclopedia of Atmospheric Sciences.2nd ed.,G.R.North et al.,Eds.,Elsevier Ltd.,47–60.

    Jakobson,E.,and T.Vihma,2010:Atmospheric moisture budget in the Arctic based on the ERA-40 reanalysis.International Journal of Climatology,30,2175–2194,http://dx.doi.org/10.1002/joc.2039.

    Kalnay,E.,and Coauthors,1996:The NCEP/NCAR 40-year reanalysis project.Bull.Amer.Meteor.Soc.,77,437–470, http://dx.doi.org/10.1175/1520-0477(1996)077<0437:TNYRP>2.0.CO;2.

    Kay,J.E.,and A.Gettelman,2009:Cloud influence on and response to seasonal Arctic sea ice loss.J.Geophys.Res.,114,D18204,http://dx.doi.org/10.1029/2009JD011773.

    Kay,J.E.,T.L’Ecuyer,A.Gettelman,G.Stephens,and C.O’Dell,2008:The contribution of cloud and radiation anomalies to the 2007 Arctic sea ice extent minimum.Geophys.Res.Lett.,35,L08503,http://dx.doi.org/10.1029/2008GL033451.

    Kim,B.-M.,and Coauthors,2017:Major cause of unprecedented arctic warming in January 2016:Critical role of an Atlantic windstorm.Scientific Reports,7,40051,http://dx.doi.org/10.1038/srep40051.

    Liu,C.J.,and E.A.Barnes,2015:Extreme moisture transport into the Arctic linked to Rossby wave breaking.J.Geophys.Res.:Atmos.,120,3774–3788,http://dx.doi.org/10.1002/2014JD 022796.

    Murray,R.J.,and I.Simmonds,1995:Responses of climate and cyclones to reductions in Arctic winter sea ice.J.Geophys.Res.,100,4791–4806,http://dx.doi.org/10.1029/94JC02206.

    Newell,R.E.,N.E.Newell,Y.Zhu,and C.Scott,1992:Tropospheric rivers?—A pilot study.Geophys.Res.Lett.,19,2401–2404,http://dx.doi.org/10.1029/92GL02916.

    Newman,M.,G.N.Kiladis,K.M.Weickmann,F.M.Ralph,and P.D.Sardeshmukh,2012:Relative contributions of synoptic and low-frequency eddies to time-mean atmospheric moisture transport,including the role of atmospheric rivers.J.Climate,25,7341–7361,http://dx.doi.org/10.1175/JCLI-D-11-00665.1.

    Peixoto,J.P.,and A.H.Oort,1992:Physics of Climate.American Institute of Physics,520 pp.

    Ralph,F.M.,and M.D.Dettinger,2011: Storms,floods,and the science of atmospheric rivers.EOS,92,265–266,http://dx.doi.org/10.1029/2011EO320001.

    Rogers,A.N.,D.H.Bromwich,E.N.Sinclair,and R.I.Cullather,2001:The atmospheric hydrologic cycle over the Arctic basin from reanalyses.Part II:Interannual variability.J.Climate,14,2414–2429,http://dx.doi.org/10.1175/1520-0442(2001)014<2414:TAHCOT>2.0.CO;2.

    Sepp,M.,and J.Jaagus,2011:Changes in the activity and tracks of arctic cyclones.Climatic Change,105,577–595,http://dx.doi.org/10.1007/s10584-010-9893-7.

    Serreze,M.C.,R.G.Barry,and J.E.Walsh,1995:Atmospheric water vapor characteristics at 70°N.J.Climate,8,719–731, http://dx.doi.org/10.1175/1520-0442(1995)008<0719:AWVCA>2.0.CO;2.

    Simmonds,I.,and K.Keay,2009:Extraordinary September Arctic sea ice reductions and their relationships with storm behavior over 1979–2008.Geophys.Res.Let.,36,L19715,http://dx.doi.org/10.1029/2009GL039810.

    Simmonds,I.,C.Burke,and K.Keay,2008:Arctic climate change as manifest in cyclone behavior.J.Climate,21,5777–5796,http://dx.doi.org/10.1175/2008JCLI2366.1.

    Skific,N.,J.A.Francis,and J.J.Cassano,2009:Attribution of projected changes in atmospheric moisture transport in the Arctic:A self-organizing map perspective.J.Climate,22,4135–4153,http://dx.doi.org/10.1175/2009JCLI2645.1.

    Sorteberg,A.,and B.Kvingedal,2006:Atmospheric forcing on the Barents Sea winter ice extent.J.Climate,19,4772–4784,http://dx.doi.org/10.1175/JCLI3885.1.

    Taylor,W.A.,2015:Change-Point Analysis:A powerful new tool for detecting changes.Taylor Enterprises.,[Available online from http://www.variation.com/cpa/tech/changepoint.html]

    Thompson,D.J.W.,and J.M.Wallace,1998:The Arctic oscillation signature in the wintertime geopotential height and temperature fields.Geophys.Res.Lett.,25,1297–1300,http://dx.doi.org/10.1029/98GL00950.

    Vihma,T.,and Coauthors,2016:The atmospheric role in the Arctic water cycle:A review on processes,past and future changes,and their impacts.J.Geophys.Res.:Biogeosci.,121,586–620,http://dx.doi.org/10.1002/2015JG003132.

    Wang,X.J.,and J.R.Key,2005:Arctic surface,cloud,and radiation properties based on the AVHRR polar pathfinder dataset.Part II.Recent trends.J.Climate,18,2575–2593,https://doi.org/10.1175/JCLI3439.1.

    Woods,C.,and R.Caballero,2016:The role of moist intrusions in winter Arctic warming and sea ice decline.J.Climate,29,4473–4485,http://dx.doi.org/10.1175/JCLI-D-15-0773.1.

    Woods,C.,R.Caballero,and G.Svensson,2013:Large-scale circulation associated with moisture intrusions into the Arctic during winter.Geophys.Res.Lett.,40,4717–4721,http://dx.doi.org/10.1002/grl.50912.

    Yin,J.H.,2005:A consistent poleward shift of the storm tracks in simulations of 21st century climate.Geophys.Res.Lett.,32,L18701,http://dx.doi.org/10.1029/2005GL023684.

    Zhang,X.D.,J.E.Walsh,J.Zhang,U.S.Bhatt,and M.Ikeda,2004:Climatology and interannual variability of Arctic cyclone activity:1948–2002.J.Climate,17,2300–2317,http://dx.doi.org/10.1175/1520-0442(2004)017<2300:CAIVOA>2.0.CO;2.

    Zhang,X.D.,A.Sorteberg,J.Zhang,R.Gerdes,and J.C.Comiso,2008:Recent radical shifts of atmospheric circulations and rapid changes in Arctic climate system.Geophys.Res.Lett.,35,L22701,http://dx.doi.org/10.1029/2008GL035607.

    Zhang,X.D.,J.X.He,J.Zhang,I.Polyakov,R.Gerdes,J.Inoue,and P.L.Wu,2012:Enhanced poleward moisture transport and amplified northern high-latitude wetting trend.Nat.Clim.Change,3,47–51,http://dx.doi.org/10.1038/nclimate1631.

    Zhu,Y.,and R.E.Newell,1998:A proposed algorithm for moisture fluxes from atmospheric rivers.Mon.Wea.Rev.,126,725–735,http://dx.doi.org/10.1175/1520-0493(1998)126<0725:APAFMF>2.0.CO;2.

    Zuidema,P.,and Coauthors,2005:An Arctic springtime mixedphase cloudy boundary layer observed during SHEBA.J.Atmos.Sci.,62,160–176,http://dx.doi.org/10.1175/JAS-3368.1.

    5 May 2017;revised 14 October 2017;accepted 17 October 2017)

    :Villamil-Otero,G.A.,J.Zhang,J.X.He,and X.Zhang,2018:Role of extratropical cyclones in the recently observed increase in poleward moisture transport into the Arctic Ocean.Adv.Atmos.Sci.,35(1),85–94,https://doi.org/10.1007/s00376-017-7116-0.

    ?Corresponding author:Gian VILLAMIL-OTERO

    Email:gavillam@ncat.edu

    ?Institute of Atmospheric Physics/Chinese Academy of Sciences,and Science Press and Springer-Verlag GmbH Germany 2018

    女同久久另类99精品国产91| www日本黄色视频网| 成年版毛片免费区| www日本在线高清视频| 天堂av国产一区二区熟女人妻 | 观看免费一级毛片| 国产精品久久电影中文字幕| 亚洲第一欧美日韩一区二区三区| 狂野欧美激情性xxxx| 老司机福利观看| 精品欧美一区二区三区在线| aaaaa片日本免费| 欧美日韩一级在线毛片| 日韩有码中文字幕| 精品欧美国产一区二区三| 成人手机av| 一本一本综合久久| 禁无遮挡网站| a在线观看视频网站| 久久九九热精品免费| 久9热在线精品视频| 青草久久国产| 18禁裸乳无遮挡免费网站照片| 在线观看免费日韩欧美大片| 十八禁人妻一区二区| 久久午夜亚洲精品久久| 国产精华一区二区三区| 给我免费播放毛片高清在线观看| 国产精品久久视频播放| 色播亚洲综合网| a级毛片在线看网站| 精品久久久久久久末码| 一个人观看的视频www高清免费观看 | 大型av网站在线播放| 岛国在线免费视频观看| 日韩欧美免费精品| 亚洲av熟女| 这个男人来自地球电影免费观看| 欧美黑人巨大hd| 中亚洲国语对白在线视频| 欧美在线一区亚洲| 国模一区二区三区四区视频 | 亚洲黑人精品在线| 亚洲国产欧美一区二区综合| 亚洲黑人精品在线| 国产亚洲av嫩草精品影院| av福利片在线| 精品无人区乱码1区二区| 老司机午夜十八禁免费视频| 久久久久九九精品影院| 成年人黄色毛片网站| 18禁裸乳无遮挡免费网站照片| 午夜精品久久久久久毛片777| 三级毛片av免费| aaaaa片日本免费| 老司机午夜福利在线观看视频| svipshipincom国产片| 国产人伦9x9x在线观看| 欧美大码av| 欧美大码av| 窝窝影院91人妻| 国产精品亚洲美女久久久| 国产精品国产高清国产av| 中文亚洲av片在线观看爽| 亚洲电影在线观看av| 国产精品综合久久久久久久免费| 最好的美女福利视频网| 五月伊人婷婷丁香| 精品国产美女av久久久久小说| 亚洲精品国产精品久久久不卡| 久久中文字幕人妻熟女| 亚洲欧美一区二区三区黑人| 中国美女看黄片| 成人手机av| 精品一区二区三区四区五区乱码| 夜夜夜夜夜久久久久| 国产成人啪精品午夜网站| 亚洲在线自拍视频| 中文资源天堂在线| 欧美日韩精品网址| 精品国产亚洲在线| 欧美日韩中文字幕国产精品一区二区三区| 国产av不卡久久| 叶爱在线成人免费视频播放| 久久亚洲真实| 高清毛片免费观看视频网站| 久久久久九九精品影院| 午夜激情福利司机影院| 亚洲男人的天堂狠狠| 黑人欧美特级aaaaaa片| netflix在线观看网站| 少妇熟女aⅴ在线视频| 亚洲人成网站在线播放欧美日韩| 99久久久亚洲精品蜜臀av| 波多野结衣高清无吗| avwww免费| 亚洲熟妇熟女久久| 黄频高清免费视频| 久久久久国内视频| 非洲黑人性xxxx精品又粗又长| 首页视频小说图片口味搜索| 国产精品永久免费网站| xxx96com| 在线观看舔阴道视频| 国产精品98久久久久久宅男小说| 国产爱豆传媒在线观看 | 久99久视频精品免费| 久久精品国产99精品国产亚洲性色| 美女扒开内裤让男人捅视频| 天堂影院成人在线观看| 黄色毛片三级朝国网站| 日本熟妇午夜| 亚洲 欧美 日韩 在线 免费| 午夜免费成人在线视频| 国产又色又爽无遮挡免费看| 亚洲av美国av| 欧美一区二区国产精品久久精品 | 久久久国产成人精品二区| 亚洲男人的天堂狠狠| 免费在线观看日本一区| 在线十欧美十亚洲十日本专区| 国产亚洲精品久久久久久毛片| 老汉色∧v一级毛片| 天天躁狠狠躁夜夜躁狠狠躁| 日韩精品免费视频一区二区三区| 国产高清视频在线观看网站| 俄罗斯特黄特色一大片| 久久天躁狠狠躁夜夜2o2o| 后天国语完整版免费观看| 国产区一区二久久| 亚洲第一电影网av| 亚洲中文日韩欧美视频| 老司机在亚洲福利影院| 嫩草影视91久久| 国产精品香港三级国产av潘金莲| 国产蜜桃级精品一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 亚洲五月天丁香| bbb黄色大片| 成熟少妇高潮喷水视频| 日本一本二区三区精品| 91九色精品人成在线观看| 变态另类成人亚洲欧美熟女| 久久久久久九九精品二区国产 | 精品欧美国产一区二区三| 欧美中文日本在线观看视频| 级片在线观看| 一级毛片高清免费大全| 免费无遮挡裸体视频| 亚洲全国av大片| 国产一区二区在线观看日韩 | 欧美日韩亚洲综合一区二区三区_| 久久人妻福利社区极品人妻图片| 一本一本综合久久| 欧美黄色淫秽网站| 久久伊人香网站| 黑人欧美特级aaaaaa片| 亚洲成av人片在线播放无| 成人特级黄色片久久久久久久| 成人三级黄色视频| 国产精品九九99| 最近最新中文字幕大全免费视频| 在线观看免费午夜福利视频| 亚洲国产精品久久男人天堂| 亚洲成av人片在线播放无| 亚洲精品在线美女| 午夜激情av网站| 久久久久性生活片| 国产伦一二天堂av在线观看| 免费在线观看视频国产中文字幕亚洲| 色老头精品视频在线观看| 久久久久久久久免费视频了| 亚洲午夜理论影院| 色哟哟哟哟哟哟| 国产免费av片在线观看野外av| av在线天堂中文字幕| 国产av又大| 亚洲第一电影网av| 无限看片的www在线观看| 妹子高潮喷水视频| 国内揄拍国产精品人妻在线| 中文字幕熟女人妻在线| 国产欧美日韩一区二区精品| 国产精品永久免费网站| 久久久久性生活片| 欧美最黄视频在线播放免费| 一边摸一边做爽爽视频免费| av国产免费在线观看| 90打野战视频偷拍视频| 嫁个100分男人电影在线观看| 亚洲欧美精品综合久久99| 久久精品国产综合久久久| 三级国产精品欧美在线观看 | 亚洲精品久久成人aⅴ小说| 国产黄色小视频在线观看| 狂野欧美白嫩少妇大欣赏| 一级毛片女人18水好多| 手机成人av网站| 国产亚洲精品第一综合不卡| 午夜福利免费观看在线| 亚洲色图 男人天堂 中文字幕| 高清在线国产一区| 欧美激情久久久久久爽电影| 久久久久亚洲av毛片大全| www.999成人在线观看| a在线观看视频网站| 特大巨黑吊av在线直播| 韩国av一区二区三区四区| 午夜视频精品福利| 美女高潮喷水抽搐中文字幕| 日本 欧美在线| 精品午夜福利视频在线观看一区| 亚洲午夜理论影院| 999久久久国产精品视频| 99精品欧美一区二区三区四区| 在线观看午夜福利视频| 日本在线视频免费播放| 欧美 亚洲 国产 日韩一| 最好的美女福利视频网| 午夜福利免费观看在线| 日韩有码中文字幕| 一本精品99久久精品77| 亚洲色图av天堂| 无遮挡黄片免费观看| 亚洲av日韩精品久久久久久密| 真人一进一出gif抽搐免费| 久久久精品大字幕| 欧美成人午夜精品| 日韩精品免费视频一区二区三区| 岛国在线免费视频观看| 欧美中文日本在线观看视频| 亚洲av片天天在线观看| 久久久久久人人人人人| 99热只有精品国产| a在线观看视频网站| 欧美大码av| 黑人巨大精品欧美一区二区mp4| 国产主播在线观看一区二区| 51午夜福利影视在线观看| 国产真实乱freesex| 久久国产精品人妻蜜桃| 亚洲成人免费电影在线观看| 1024香蕉在线观看| 精品久久久久久久末码| 色av中文字幕| 欧美日韩瑟瑟在线播放| 精品久久久久久,| 我要搜黄色片| 亚洲成人久久爱视频| 人妻丰满熟妇av一区二区三区| 亚洲成人免费电影在线观看| 午夜免费观看网址| 色精品久久人妻99蜜桃| 国产久久久一区二区三区| 超碰成人久久| 久久久国产欧美日韩av| 成在线人永久免费视频| 18禁黄网站禁片午夜丰满| 精品久久久久久久久久免费视频| av免费在线观看网站| 男女做爰动态图高潮gif福利片| 精品久久久久久久末码| 国产av一区在线观看免费| 精品少妇一区二区三区视频日本电影| 两个人免费观看高清视频| 亚洲国产精品成人综合色| 欧美成人免费av一区二区三区| 久久久久国产精品人妻aⅴ院| 人妻夜夜爽99麻豆av| 日本 av在线| 国产黄色小视频在线观看| 国产成人一区二区三区免费视频网站| 日日夜夜操网爽| 一区二区三区高清视频在线| av有码第一页| 可以在线观看的亚洲视频| 亚洲av中文字字幕乱码综合| 久久精品人妻少妇| 视频区欧美日本亚洲| 99久久国产精品久久久| 日日夜夜操网爽| 日韩免费av在线播放| 老汉色av国产亚洲站长工具| 欧美乱妇无乱码| 欧洲精品卡2卡3卡4卡5卡区| 欧美成人午夜精品| 99精品欧美一区二区三区四区| 久久午夜亚洲精品久久| 一级毛片精品| 国产三级中文精品| 国产精品av视频在线免费观看| 久久天躁狠狠躁夜夜2o2o| 国产高清videossex| 麻豆成人av在线观看| 少妇被粗大的猛进出69影院| 1024视频免费在线观看| 久久热在线av| 精品国产超薄肉色丝袜足j| 亚洲18禁久久av| 制服诱惑二区| 日日夜夜操网爽| 99国产极品粉嫩在线观看| 男女之事视频高清在线观看| 国产精品 欧美亚洲| 国产精品久久久久久久电影 | 国产亚洲精品一区二区www| 亚洲成人中文字幕在线播放| 一级片免费观看大全| 中文字幕最新亚洲高清| 亚洲人成网站在线播放欧美日韩| 亚洲欧美日韩无卡精品| 最近视频中文字幕2019在线8| 免费看日本二区| 国产精品亚洲一级av第二区| 亚洲激情在线av| 日本一二三区视频观看| 男男h啪啪无遮挡| av国产免费在线观看| 日韩高清综合在线| 看免费av毛片| 成人国产一区最新在线观看| 欧美日韩亚洲综合一区二区三区_| 很黄的视频免费| 激情在线观看视频在线高清| 黄色女人牲交| 欧美成狂野欧美在线观看| 波多野结衣巨乳人妻| 18禁黄网站禁片免费观看直播| 午夜福利在线观看吧| 成人欧美大片| 久久人人精品亚洲av| 日韩av在线大香蕉| 中文字幕最新亚洲高清| 免费在线观看影片大全网站| videosex国产| 成人三级做爰电影| 日韩中文字幕欧美一区二区| 亚洲国产欧美人成| 老司机午夜福利在线观看视频| 一本综合久久免费| 他把我摸到了高潮在线观看| 久久久国产成人免费| netflix在线观看网站| 国产主播在线观看一区二区| 国产熟女午夜一区二区三区| 国产精品亚洲美女久久久| 久久99热这里只有精品18| 亚洲一卡2卡3卡4卡5卡精品中文| 中文字幕高清在线视频| 在线视频色国产色| 亚洲精品久久国产高清桃花| 最新美女视频免费是黄的| 久久久久亚洲av毛片大全| 男女视频在线观看网站免费 | 免费电影在线观看免费观看| 国产成人av教育| 亚洲va日本ⅴa欧美va伊人久久| 国产91精品成人一区二区三区| 国产真实乱freesex| 成人高潮视频无遮挡免费网站| 美女黄网站色视频| 夜夜夜夜夜久久久久| 国产亚洲精品一区二区www| 国产欧美日韩一区二区三| 操出白浆在线播放| 欧美绝顶高潮抽搐喷水| 欧美黑人欧美精品刺激| 成年免费大片在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲av日韩精品久久久久久密| 欧美另类亚洲清纯唯美| 老熟妇仑乱视频hdxx| 国产黄色小视频在线观看| 国产成人av激情在线播放| 成人三级黄色视频| 麻豆av在线久日| 高潮久久久久久久久久久不卡| 999久久久国产精品视频| 啪啪无遮挡十八禁网站| 国内精品久久久久精免费| 久久人人精品亚洲av| 中国美女看黄片| 国产亚洲精品av在线| 国产av一区二区精品久久| 麻豆国产av国片精品| 在线国产一区二区在线| 国产精品av视频在线免费观看| 18禁黄网站禁片午夜丰满| 成人国产综合亚洲| 19禁男女啪啪无遮挡网站| 成人亚洲精品av一区二区| 黄色 视频免费看| 69av精品久久久久久| 中文字幕久久专区| 少妇人妻一区二区三区视频| 给我免费播放毛片高清在线观看| 久久精品国产综合久久久| 黄片小视频在线播放| 久久久精品国产亚洲av高清涩受| 搞女人的毛片| 老汉色av国产亚洲站长工具| 欧美黄色淫秽网站| 亚洲欧美精品综合一区二区三区| 国产一区二区激情短视频| 国产激情久久老熟女| 午夜视频精品福利| 母亲3免费完整高清在线观看| 脱女人内裤的视频| 中文字幕人成人乱码亚洲影| 好男人在线观看高清免费视频| 不卡av一区二区三区| 男女床上黄色一级片免费看| 麻豆国产av国片精品| 又黄又粗又硬又大视频| 国产精品免费一区二区三区在线| 国产精品亚洲av一区麻豆| 亚洲人成网站在线播放欧美日韩| 亚洲人成77777在线视频| 九色国产91popny在线| 最好的美女福利视频网| 男插女下体视频免费在线播放| 亚洲一区二区三区不卡视频| 1024香蕉在线观看| 亚洲成a人片在线一区二区| 久久香蕉精品热| а√天堂www在线а√下载| 国产1区2区3区精品| 亚洲激情在线av| 成年女人毛片免费观看观看9| 欧美日韩黄片免| 亚洲精品中文字幕一二三四区| 国产精品电影一区二区三区| av欧美777| 97超级碰碰碰精品色视频在线观看| 午夜两性在线视频| 少妇裸体淫交视频免费看高清 | 男女午夜视频在线观看| 一级a爱片免费观看的视频| 成人手机av| 色老头精品视频在线观看| or卡值多少钱| 成人一区二区视频在线观看| 午夜影院日韩av| 精品欧美国产一区二区三| 老熟妇乱子伦视频在线观看| 亚洲av电影在线进入| 黄频高清免费视频| 日韩大尺度精品在线看网址| av在线播放免费不卡| 免费在线观看成人毛片| 国产亚洲精品久久久久久毛片| 日本免费a在线| 一级毛片高清免费大全| 免费av毛片视频| 色综合欧美亚洲国产小说| 又粗又爽又猛毛片免费看| 久久久久九九精品影院| 国产亚洲精品综合一区在线观看 | 久久人妻av系列| avwww免费| 欧美色欧美亚洲另类二区| 国产爱豆传媒在线观看 | 国产乱人伦免费视频| 俄罗斯特黄特色一大片| 啦啦啦韩国在线观看视频| 国产私拍福利视频在线观看| 久久久国产成人精品二区| 国产爱豆传媒在线观看 | av免费在线观看网站| 国产av不卡久久| 级片在线观看| 日韩 欧美 亚洲 中文字幕| 国产亚洲欧美98| 91字幕亚洲| 欧美色视频一区免费| 久久精品亚洲精品国产色婷小说| 亚洲熟妇熟女久久| 亚洲人成网站高清观看| 欧美黑人精品巨大| 国产日本99.免费观看| 午夜两性在线视频| 母亲3免费完整高清在线观看| 一本大道久久a久久精品| 91成年电影在线观看| 人妻丰满熟妇av一区二区三区| 香蕉av资源在线| 亚洲国产中文字幕在线视频| av中文乱码字幕在线| 精品无人区乱码1区二区| 成人av一区二区三区在线看| 这个男人来自地球电影免费观看| 日本熟妇午夜| 国产精品自产拍在线观看55亚洲| 一级毛片女人18水好多| 久久中文看片网| svipshipincom国产片| 日韩成人在线观看一区二区三区| x7x7x7水蜜桃| 日韩免费av在线播放| 麻豆成人午夜福利视频| 全区人妻精品视频| 99国产精品一区二区蜜桃av| 日韩国内少妇激情av| 亚洲九九香蕉| 亚洲中文日韩欧美视频| 天堂√8在线中文| 91老司机精品| 亚洲成人国产一区在线观看| 国产精品永久免费网站| 国产成人系列免费观看| 又爽又黄无遮挡网站| av天堂在线播放| 一进一出抽搐动态| 日韩免费av在线播放| 99精品久久久久人妻精品| 三级国产精品欧美在线观看 | 香蕉国产在线看| 妹子高潮喷水视频| 欧美乱妇无乱码| 国产99白浆流出| 十八禁人妻一区二区| 一级黄色大片毛片| 丝袜人妻中文字幕| 亚洲欧美精品综合一区二区三区| 搡老妇女老女人老熟妇| 欧美 亚洲 国产 日韩一| 日韩欧美免费精品| 91九色精品人成在线观看| 亚洲人成77777在线视频| 午夜日韩欧美国产| 色综合站精品国产| 看黄色毛片网站| 国产精品爽爽va在线观看网站| 国产1区2区3区精品| 国产精品爽爽va在线观看网站| 亚洲av片天天在线观看| 午夜福利免费观看在线| 麻豆一二三区av精品| 19禁男女啪啪无遮挡网站| 丁香六月欧美| 亚洲av电影在线进入| 亚洲午夜理论影院| 久久久国产精品麻豆| netflix在线观看网站| 俄罗斯特黄特色一大片| 久久久久免费精品人妻一区二区| 极品教师在线免费播放| 无限看片的www在线观看| 亚洲欧美一区二区三区黑人| 亚洲成av人片在线播放无| av在线天堂中文字幕| 成人永久免费在线观看视频| 性色av乱码一区二区三区2| 国产午夜精品论理片| 亚洲av第一区精品v没综合| 757午夜福利合集在线观看| 夜夜爽天天搞| 女同久久另类99精品国产91| 青草久久国产| 亚洲精品中文字幕在线视频| 午夜精品久久久久久毛片777| 大型黄色视频在线免费观看| 在线看三级毛片| 色在线成人网| 精品欧美一区二区三区在线| 日韩精品免费视频一区二区三区| 国产av一区在线观看免费| 久久久久久久久久黄片| 久久热在线av| 亚洲,欧美精品.| 脱女人内裤的视频| 亚洲av日韩精品久久久久久密| 亚洲午夜理论影院| 国产亚洲精品久久久久5区| 亚洲一区中文字幕在线| 麻豆av在线久日| 熟女电影av网| 成人av在线播放网站| 欧美又色又爽又黄视频| 久久精品国产99精品国产亚洲性色| 国产一区二区在线观看日韩 | 久久婷婷成人综合色麻豆| 黄色片一级片一级黄色片| 日本五十路高清| 精品日产1卡2卡| 久久精品国产亚洲av香蕉五月| 伦理电影免费视频| 中文字幕人妻丝袜一区二区| 国产激情偷乱视频一区二区| 极品教师在线免费播放| 亚洲国产日韩欧美精品在线观看 | 激情在线观看视频在线高清| 麻豆av在线久日| 日本黄大片高清| 可以在线观看的亚洲视频| 国产亚洲精品第一综合不卡| 老司机福利观看| 国产精品久久久人人做人人爽| 国内精品久久久久精免费| 90打野战视频偷拍视频| 老司机深夜福利视频在线观看| 一个人观看的视频www高清免费观看 | 国产一区二区激情短视频| 在线观看舔阴道视频| 亚洲国产精品成人综合色| 国产精品久久久久久久电影 | 国产精品香港三级国产av潘金莲| 午夜久久久久精精品| 国产欧美日韩精品亚洲av| 看黄色毛片网站| 两个人视频免费观看高清| 亚洲成人免费电影在线观看| 国产成人aa在线观看| 亚洲国产精品sss在线观看| www.熟女人妻精品国产|