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

    亚洲精品成人av观看孕妇| 色综合色国产| 欧美3d第一页| 亚洲综合精品二区| 亚洲成人手机| 伦理电影免费视频| 日韩成人av中文字幕在线观看| 97在线人人人人妻| 欧美极品一区二区三区四区| 国产视频内射| 男人爽女人下面视频在线观看| 欧美成人一区二区免费高清观看| 嫩草影院新地址| 欧美日韩视频高清一区二区三区二| 丰满乱子伦码专区| 亚洲av男天堂| 亚洲一级一片aⅴ在线观看| 一本色道久久久久久精品综合| 日本vs欧美在线观看视频 | 国产成人精品一,二区| 亚洲精品一区蜜桃| 免费黄频网站在线观看国产| 成人美女网站在线观看视频| 大又大粗又爽又黄少妇毛片口| 免费黄网站久久成人精品| 精品亚洲乱码少妇综合久久| 中文字幕制服av| 中文在线观看免费www的网站| 18禁裸乳无遮挡免费网站照片| 亚洲内射少妇av| 国产精品欧美亚洲77777| 久久久久视频综合| 久久久久视频综合| 亚洲av电影在线观看一区二区三区| 久久99热这里只有精品18| 日韩制服骚丝袜av| 少妇人妻 视频| 欧美日韩视频精品一区| 免费人妻精品一区二区三区视频| 高清黄色对白视频在线免费看 | h日本视频在线播放| 日韩国内少妇激情av| 日本-黄色视频高清免费观看| 你懂的网址亚洲精品在线观看| 肉色欧美久久久久久久蜜桃| 成人毛片a级毛片在线播放| 永久网站在线| 永久免费av网站大全| 蜜臀久久99精品久久宅男| 亚洲真实伦在线观看| 1000部很黄的大片| av在线蜜桃| 亚洲精品久久久久久婷婷小说| 蜜臀久久99精品久久宅男| 国产一区二区三区av在线| 亚洲av成人精品一二三区| 日韩三级伦理在线观看| 亚洲欧洲日产国产| 一二三四中文在线观看免费高清| 国产成人免费无遮挡视频| 国精品久久久久久国模美| 久久精品久久精品一区二区三区| 午夜福利影视在线免费观看| 亚洲国产最新在线播放| 伦理电影免费视频| 亚洲精品国产成人久久av| 卡戴珊不雅视频在线播放| 春色校园在线视频观看| 秋霞在线观看毛片| 观看免费一级毛片| av在线app专区| 日韩在线高清观看一区二区三区| 欧美最新免费一区二区三区| 中文字幕av成人在线电影| 在线观看一区二区三区| 欧美日韩视频精品一区| 少妇人妻久久综合中文| 国产男女内射视频| av卡一久久| 久久精品久久久久久久性| 秋霞在线观看毛片| 又爽又黄a免费视频| 另类亚洲欧美激情| 人体艺术视频欧美日本| 久久久午夜欧美精品| xxx大片免费视频| 久久 成人 亚洲| 毛片女人毛片| 国产精品一区www在线观看| 国产精品久久久久久精品古装| 内地一区二区视频在线| 老女人水多毛片| 亚洲第一av免费看| 日韩在线高清观看一区二区三区| 国产有黄有色有爽视频| 纵有疾风起免费观看全集完整版| 久久久精品免费免费高清| 国产成人aa在线观看| 亚洲av在线观看美女高潮| 久久婷婷青草| 中文资源天堂在线| 亚洲色图av天堂| 国产精品人妻久久久久久| 大又大粗又爽又黄少妇毛片口| 18+在线观看网站| 久久鲁丝午夜福利片| 欧美日韩一区二区视频在线观看视频在线| 亚洲精品aⅴ在线观看| av网站免费在线观看视频| 人体艺术视频欧美日本| 看十八女毛片水多多多| 一级av片app| 最近最新中文字幕大全电影3| 亚洲美女搞黄在线观看| 国产av一区二区精品久久 | 亚洲精品视频女| 久久久亚洲精品成人影院| 亚洲精品国产av成人精品| 精品久久久精品久久久| 久久国产亚洲av麻豆专区| 人人妻人人添人人爽欧美一区卜 | 国产高潮美女av| 久久精品国产鲁丝片午夜精品| 久久久a久久爽久久v久久| 免费看光身美女| 91久久精品电影网| 国产精品福利在线免费观看| 少妇被粗大猛烈的视频| 亚洲精品乱久久久久久| 黄色怎么调成土黄色| 插阴视频在线观看视频| 少妇熟女欧美另类| 色视频www国产| 少妇精品久久久久久久| 久久久欧美国产精品| 色综合色国产| 男的添女的下面高潮视频| 久久精品国产自在天天线| 亚洲欧洲日产国产| 日日摸夜夜添夜夜爱| 97精品久久久久久久久久精品| 国产亚洲一区二区精品| 日本欧美视频一区| 伦理电影大哥的女人| av在线蜜桃| 菩萨蛮人人尽说江南好唐韦庄| 成人特级av手机在线观看| 久久韩国三级中文字幕| 免费在线观看成人毛片| 新久久久久国产一级毛片| 熟女av电影| 国产高清不卡午夜福利| 人妻夜夜爽99麻豆av| 国产精品福利在线免费观看| 少妇的逼好多水| 成人漫画全彩无遮挡| 亚洲欧美一区二区三区国产| 亚洲精品乱码久久久v下载方式| 亚洲人成网站在线观看播放| 国产av码专区亚洲av| 精品国产乱码久久久久久小说| 欧美精品国产亚洲| 亚洲国产精品一区三区| 久久久久网色| 色视频www国产| 国产欧美日韩精品一区二区| 人人妻人人澡人人爽人人夜夜| 交换朋友夫妻互换小说| 国产成人免费观看mmmm| 中国国产av一级| 婷婷色综合www| 国产爽快片一区二区三区| 午夜免费男女啪啪视频观看| 久久久久久久久久久免费av| 欧美精品人与动牲交sv欧美| 色吧在线观看| 国产精品三级大全| 美女cb高潮喷水在线观看| 51国产日韩欧美| 国产乱人偷精品视频| 99re6热这里在线精品视频| 国产爱豆传媒在线观看| 国产精品福利在线免费观看| 搡老乐熟女国产| 一级黄片播放器| 日日啪夜夜爽| 国产白丝娇喘喷水9色精品| 这个男人来自地球电影免费观看 | 亚洲国产精品一区三区| 久久精品夜色国产| 欧美一级a爱片免费观看看| 国产av码专区亚洲av| 久久久久视频综合| 一级毛片我不卡| 免费在线观看成人毛片| 久热这里只有精品99| 熟女av电影| 国产女主播在线喷水免费视频网站| 亚洲精品,欧美精品| 精品久久国产蜜桃| 少妇猛男粗大的猛烈进出视频| 国产欧美日韩一区二区三区在线 | 卡戴珊不雅视频在线播放| 国产精品伦人一区二区| 国产精品人妻久久久久久| 亚洲第一区二区三区不卡| 九九久久精品国产亚洲av麻豆| 亚洲色图综合在线观看| 性色av一级| 欧美日本视频| 在线看a的网站| 亚洲av男天堂| av国产精品久久久久影院| 麻豆成人午夜福利视频| 精品人妻熟女av久视频| 亚洲不卡免费看| 蜜桃亚洲精品一区二区三区| 欧美日本视频| 欧美亚洲 丝袜 人妻 在线| 色哟哟·www| 国产精品熟女久久久久浪| 嘟嘟电影网在线观看| 久久精品熟女亚洲av麻豆精品| 国产一区二区三区综合在线观看 | 成人国产麻豆网| 国产黄片视频在线免费观看| 色视频www国产| 18+在线观看网站| 只有这里有精品99| 搡老乐熟女国产| 亚洲成人一二三区av| 国产一区二区三区综合在线观看 | 久久99热这里只有精品18| 久久国产亚洲av麻豆专区| 成年av动漫网址| 免费看光身美女| 婷婷色综合大香蕉| 深爱激情五月婷婷| 美女视频免费永久观看网站| 岛国毛片在线播放| 亚洲精品日韩av片在线观看| 在线观看美女被高潮喷水网站| 国产精品国产三级国产专区5o| 亚洲精品一二三| 男女啪啪激烈高潮av片| 各种免费的搞黄视频| 插阴视频在线观看视频| 亚洲国产毛片av蜜桃av| 国产女主播在线喷水免费视频网站| 午夜激情久久久久久久| 韩国高清视频一区二区三区| 久久久久国产网址| 99久久中文字幕三级久久日本| 亚洲av免费高清在线观看| 精品久久久噜噜| 日韩人妻高清精品专区| 亚洲精品国产成人久久av| 少妇熟女欧美另类| 日日摸夜夜添夜夜添av毛片| 只有这里有精品99| 免费大片18禁| 美女视频免费永久观看网站| 内射极品少妇av片p| av免费观看日本| 国产免费视频播放在线视频| 丝瓜视频免费看黄片| 大陆偷拍与自拍| av又黄又爽大尺度在线免费看| 亚洲欧美日韩另类电影网站 | 肉色欧美久久久久久久蜜桃| 少妇人妻久久综合中文| 亚洲va在线va天堂va国产| 亚洲精品乱码久久久v下载方式| 插逼视频在线观看| 搡女人真爽免费视频火全软件| 国产免费一区二区三区四区乱码| 国产精品久久久久久精品古装| 中文资源天堂在线| 3wmmmm亚洲av在线观看| 99久久综合免费| 新久久久久国产一级毛片| 制服丝袜香蕉在线| 亚洲欧洲日产国产| 亚洲精品成人av观看孕妇| 蜜桃在线观看..| 亚洲激情五月婷婷啪啪| 在线观看三级黄色| 亚洲精品国产av成人精品| 日韩电影二区| 亚洲精品乱码久久久久久按摩| 成人漫画全彩无遮挡| 亚洲欧美日韩另类电影网站 | 免费观看a级毛片全部| av在线播放精品| 国产色爽女视频免费观看| 久久人妻熟女aⅴ| 国产精品欧美亚洲77777| 欧美 日韩 精品 国产| 国产黄片美女视频| 777米奇影视久久| 免费黄色在线免费观看| 亚洲国产日韩一区二区| 美女国产视频在线观看| 国产精品久久久久久av不卡| 精品人妻一区二区三区麻豆| 久久久久久久久久人人人人人人| 日韩不卡一区二区三区视频在线| 国产精品久久久久成人av| 国产成人a区在线观看| 人体艺术视频欧美日本| 中国三级夫妇交换| 欧美日韩国产mv在线观看视频 | 久久久久久久久大av| 日韩 亚洲 欧美在线| 亚洲在久久综合| 日韩大片免费观看网站| 老司机影院成人| 国产精品不卡视频一区二区| 日韩电影二区| 免费人妻精品一区二区三区视频| 久久99热这里只频精品6学生| 中文字幕亚洲精品专区| 亚洲成人av在线免费| 99久久精品热视频| 国产精品熟女久久久久浪| 久久婷婷青草| 久久久午夜欧美精品| 观看美女的网站| av线在线观看网站| 国产免费福利视频在线观看| 噜噜噜噜噜久久久久久91| 国产亚洲一区二区精品| 麻豆成人av视频| 久久精品久久久久久噜噜老黄| 99久久精品国产国产毛片| h视频一区二区三区| 国产极品天堂在线| 国产成人91sexporn| 毛片女人毛片| 91精品一卡2卡3卡4卡| 又大又黄又爽视频免费| 亚洲欧美清纯卡通| 国产精品国产av在线观看| 久久精品国产自在天天线| 一本久久精品| kizo精华| 国产伦精品一区二区三区视频9| 在现免费观看毛片| 久久毛片免费看一区二区三区| 亚洲欧美精品专区久久| 一区二区三区乱码不卡18| 久久精品国产亚洲av天美| av天堂中文字幕网| 国产成人免费观看mmmm| 水蜜桃什么品种好| 女性生殖器流出的白浆| 国产精品人妻久久久影院| 久久久久性生活片| 三级经典国产精品| 日韩制服骚丝袜av| 在现免费观看毛片| 国产一区二区三区综合在线观看 | 狂野欧美激情性xxxx在线观看| 亚洲图色成人| av网站免费在线观看视频| 久久人人爽av亚洲精品天堂 | 日本wwww免费看| 下体分泌物呈黄色| 欧美人与善性xxx| 亚洲av.av天堂| 韩国高清视频一区二区三区| 黄色日韩在线| 国产精品国产三级国产av玫瑰| 女性被躁到高潮视频| 内射极品少妇av片p| 最近最新中文字幕免费大全7| 日韩一区二区三区影片| 亚洲国产高清在线一区二区三| 免费大片18禁| 国产成人精品一,二区| 最黄视频免费看| 简卡轻食公司| 国产黄片美女视频| 色视频在线一区二区三区| 免费人成在线观看视频色| 中文欧美无线码| 成年人午夜在线观看视频| 免费大片黄手机在线观看| 中文精品一卡2卡3卡4更新| 国产精品蜜桃在线观看| 高清黄色对白视频在线免费看 | 久久 成人 亚洲| 精品视频人人做人人爽| 日本午夜av视频| 乱系列少妇在线播放| 看十八女毛片水多多多| 国产精品三级大全| 丝瓜视频免费看黄片| 精品人妻偷拍中文字幕| 高清不卡的av网站| 乱码一卡2卡4卡精品| 久久国产乱子免费精品| 国产精品一区www在线观看| 少妇高潮的动态图| 国产国拍精品亚洲av在线观看| 国产69精品久久久久777片| 你懂的网址亚洲精品在线观看| 校园人妻丝袜中文字幕| 免费人成在线观看视频色| 少妇丰满av| 午夜激情福利司机影院| 黄色日韩在线| 国产精品精品国产色婷婷| 国产高清三级在线| 91aial.com中文字幕在线观看| 亚洲精品日本国产第一区| 久久鲁丝午夜福利片| 久久人人爽人人爽人人片va| 亚洲精品456在线播放app| 1000部很黄的大片| 91精品伊人久久大香线蕉| 蜜桃在线观看..| av又黄又爽大尺度在线免费看| 亚洲av不卡在线观看| 在线天堂最新版资源| 欧美性感艳星| 久久久久久久久久久丰满| 国产精品一区二区性色av| 91精品国产国语对白视频| 熟女av电影| 国产国拍精品亚洲av在线观看| 亚洲不卡免费看| av免费观看日本| 麻豆乱淫一区二区| 国产永久视频网站| 欧美xxxx黑人xx丫x性爽| 日韩国内少妇激情av| 日韩大片免费观看网站| 国产精品蜜桃在线观看| 日日啪夜夜撸| 深夜a级毛片| 亚洲精品成人av观看孕妇| 日本黄大片高清| 能在线免费看毛片的网站| 国产成人精品福利久久| 国产视频首页在线观看| 一个人免费看片子| 91aial.com中文字幕在线观看| 水蜜桃什么品种好| 在线天堂最新版资源| 久久韩国三级中文字幕| 日韩一本色道免费dvd| 日本黄色片子视频| 九草在线视频观看| 在线观看国产h片| 干丝袜人妻中文字幕| 男女啪啪激烈高潮av片| 亚洲四区av| 国产av精品麻豆| 在线观看美女被高潮喷水网站| av专区在线播放| 国产男人的电影天堂91| 精品人妻一区二区三区麻豆| 久久精品国产自在天天线| 久久精品国产亚洲网站| 高清午夜精品一区二区三区| 久久亚洲国产成人精品v| 高清黄色对白视频在线免费看 | 国产国拍精品亚洲av在线观看| 直男gayav资源| 亚洲欧美日韩东京热| 十分钟在线观看高清视频www | 国语对白做爰xxxⅹ性视频网站| 精品熟女少妇av免费看| 亚洲美女视频黄频| 韩国av在线不卡| 日韩av在线免费看完整版不卡| 在现免费观看毛片| 你懂的网址亚洲精品在线观看| 欧美精品亚洲一区二区| 夫妻性生交免费视频一级片| 啦啦啦视频在线资源免费观看| 熟女av电影| 国产精品熟女久久久久浪| 亚洲精华国产精华液的使用体验| 国产亚洲91精品色在线| av在线app专区| 免费不卡的大黄色大毛片视频在线观看| 欧美极品一区二区三区四区| 亚洲精品视频女| 亚洲av在线观看美女高潮| 国产精品av视频在线免费观看| 国产精品蜜桃在线观看| 国产成人精品婷婷| 女的被弄到高潮叫床怎么办| 欧美亚洲 丝袜 人妻 在线| 狂野欧美激情性bbbbbb| 日韩不卡一区二区三区视频在线| 国产美女午夜福利| 日本黄大片高清| 高清欧美精品videossex| 大陆偷拍与自拍| 国产精品久久久久成人av| 国产成人精品福利久久| 亚洲色图综合在线观看| 美女国产视频在线观看| 婷婷色av中文字幕| av播播在线观看一区| 少妇熟女欧美另类| 亚洲av男天堂| 亚洲国产成人一精品久久久| 黄片无遮挡物在线观看| 啦啦啦啦在线视频资源| 欧美精品国产亚洲| 亚洲不卡免费看| 搡女人真爽免费视频火全软件| 国产午夜精品一二区理论片| 女人久久www免费人成看片| 99久久综合免费| 一个人看的www免费观看视频| 国产成人a∨麻豆精品| 国产精品一区二区在线观看99| 人妻少妇偷人精品九色| 国产精品一区二区在线观看99| 麻豆乱淫一区二区| 中文天堂在线官网| 日本欧美视频一区| www.色视频.com| 国产淫语在线视频| 亚洲成人中文字幕在线播放| 女的被弄到高潮叫床怎么办| 黄片wwwwww| 亚洲国产成人一精品久久久| 久久久国产一区二区| 乱码一卡2卡4卡精品| 高清不卡的av网站| 亚洲人成网站高清观看| 中国国产av一级| 久久这里有精品视频免费| 狂野欧美激情性xxxx在线观看| 日韩欧美 国产精品| 五月天丁香电影| 国产精品欧美亚洲77777| videos熟女内射| 午夜激情久久久久久久| 天堂俺去俺来也www色官网| 国产在线男女| 亚洲精品日韩av片在线观看| 深夜a级毛片| 国产免费福利视频在线观看| 哪个播放器可以免费观看大片| 一本一本综合久久| 亚洲欧美一区二区三区黑人 | 又粗又硬又长又爽又黄的视频| 夜夜爽夜夜爽视频| 国产高清有码在线观看视频| 精品酒店卫生间| 久久午夜福利片| 国产一区二区在线观看日韩| 麻豆精品久久久久久蜜桃| 日韩av在线免费看完整版不卡| 秋霞在线观看毛片| 一级a做视频免费观看| 精品视频人人做人人爽| 成人特级av手机在线观看| 亚洲三级黄色毛片| 国产精品女同一区二区软件| 午夜免费观看性视频| a级毛片免费高清观看在线播放| 日本免费在线观看一区| 啦啦啦视频在线资源免费观看| 国产黄片美女视频| 国产久久久一区二区三区| 国产免费福利视频在线观看| 女性被躁到高潮视频| 又大又黄又爽视频免费| 成年免费大片在线观看| 国产极品天堂在线| 一级av片app| 草草在线视频免费看| 99久久人妻综合| 草草在线视频免费看| 老师上课跳d突然被开到最大视频| 精品一区二区免费观看| 亚洲av不卡在线观看| 国产成人freesex在线| 亚洲av电影在线观看一区二区三区| 尾随美女入室| 亚洲色图综合在线观看| 又大又黄又爽视频免费| 亚洲av不卡在线观看| 最新中文字幕久久久久| 丝瓜视频免费看黄片| 激情 狠狠 欧美| 欧美激情国产日韩精品一区| av在线播放精品| 国产一区二区在线观看日韩| 国产精品久久久久久精品古装| 一区在线观看完整版| 国产伦理片在线播放av一区| 久久久精品94久久精品| 高清视频免费观看一区二区| 六月丁香七月| 大话2 男鬼变身卡| 久久久久久人妻| xxx大片免费视频| 另类亚洲欧美激情| 青青草视频在线视频观看| 伦理电影大哥的女人| 国语对白做爰xxxⅹ性视频网站| 午夜福利在线观看免费完整高清在| 久久久久久久大尺度免费视频| 亚洲av成人精品一二三区| 国产免费又黄又爽又色| 国产高清国产精品国产三级 |