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

    Record Low Sea-Ice Concentration in the Central Arctic during Summer 2010

    2018-12-06 07:35:55JinpingZHAODavidBARBERShugangZHANGQinghuaYANGXiaoyuWANGandHongjieXIE4KeyLaboratoryofPhysicalOceanographyOceanUniversityofChinaQingdao6600China
    Advances in Atmospheric Sciences 2018年1期

    Jinping ZHAO,David BARBER,Shugang ZHANG,Qinghua YANG,Xiaoyu WANG,and Hongjie XIE4Key Laboratory of Physical Oceanography,Ocean University of China,Qingdao,6600,China

    2Centre for Earth Observation Science,Faculty of Environment Earth and Resources,University of Manitoba,

    Winnipeg,Manitoba,R3T 2N2,Canada

    3National Marine Environmental Forecasting Center,China,Beijing,100081,China

    4Department of Geological Sciences,University of Texas at San Antonio,San Antonio,Texas 78284,U.S.A.

    Record Low Sea-Ice Concentration in the Central Arctic during Summer 2010

    Jinping ZHAO?1,David BARBER2,Shugang ZHANG1,Qinghua YANG3,Xiaoyu WANG1,and Hongjie XIE41Key Laboratory of Physical Oceanography,Ocean University of China,Qingdao,266100,China

    2Centre for Earth Observation Science,Faculty of Environment Earth and Resources,University of Manitoba,

    Winnipeg,Manitoba,R3T 2N2,Canada

    3National Marine Environmental Forecasting Center,China,Beijing,100081,China

    4Department of Geological Sciences,University of Texas at San Antonio,San Antonio,Texas 78284,U.S.A.

    The Arctic sea-ice extent has shown a declining trend over the past 30 years.Ice coverage reached historic minima in 2007 and again in 2012.This trend has recently been assessed to be unique over at least the last 1450 years.In the summer of 2010,a very low sea-ice concentration(SIC)appeared at high Arctic latitudes—even lower than that of surrounding pack ice at lower latitudes.This striking low ice concentration—referred to here as a record low ice concentration in the central Arctic(CARLIC)—is unique in our analysis period of 2003–15,and has not been previously reported in the literature.The CARLIC was not the result of ice melt,because sea ice was still quite thick based on in-situ ice thickness measurements.Instead,divergent ice drift appears to have been responsible for the CARLIC.A high correlation between SIC and wind stress curl suggests that the sea ice drift during the summer of 2010 responded strongly to the regional wind forcing.The drift trajectories of ice buoys exhibited a transpolar drift in the Atlantic sector and an eastward drift in the Pacific sector,which appeared to benefit the CARLIC in 2010.Under these conditions,more solar energy can penetrate into the open water,increasing melt through increased heat flux to the ocean.We speculate that this divergence of sea ice could occur more often in the coming decades,and impact on hemispheric SIC and feed back to the climate.

    sea ice concentration,central Arctic,Beaufort Gyre,Transpolar Drift,ice motion,divergence

    1.Introduction

    Arctic sea-ice extent has been on the decline since the late 1970s(Lindsay and Zhang,2005;Lemke et al.,2007),and has recently been assessed to be unique over at least a 1450-year period,including the relatively warm Medieval Period(Kinnard et al.,2011).The summer sea-ice extent underwent a decreasing trend in the past 30 years,at a rate of more than 70 000 km2yr-1(Stroeve et al.,2007),and reached historic minima in summer 2007(Comiso et al.,2008)and in 2012(Overland and Wang,2013).At this rate,a seasonally ice-free Arctic is expected sometime over the next few decades,according to a variety of models(e.g.,Overland and Wang,2007,2013;Tietsche et al.,2011;Kay et al.,2011).The current metamorphosis from a multiyear-ice-dominated to a first-year-ice-dominated Arctic marine system appears to be a key feature of a warming planet(Lemke et al.,2007).With the reduction in multiyear ice,first-year ice has increased significantly in spatial extent and seasonal periodicity(Maslanik et al.,2011;Tschudi et al.,2016).The thinner ice and low sea-ice concentration(SIC)allow more absorption of solar radiation,which can increase the heat content of the ocean surface mixed layer(Perovich et al.,2008),and feeds back positively to more ice melt and increased energy and mass exchange across the ocean–sea-ice–atmosphere interface(Holland et al.,2006;Kumar et al.,2010,Raddatz et al.,2013).

    The rapid decline in sea-ice extent and concentration has mostly been reported in terms of the significant change in the marginal ice zone(MIZ),which is a transitional zone between open ocean and pack ice(e.g.,Stroeve et al.,2007).When the sea ice retreats and ice thickness reduces rapidly,the main feature of the ice concentration is the spread of the MIZ.The MIZ was historically a narrow zone,but has recently become much larger in extent.The low SIC area has become so wide that it may be becoming hard to distinguish the division between the MIZ and pack ice,throwing this historical definition into question.The perennial ice-extent loss is mainly caused by:ice advection toward the Canadian Archipelago coast;ice loading into the Transpolar Drift;acceleration of the Transpolar Drift carrying ice out of Fram Strait;and ice export to Baffin Bay through the Nares Strait(Nghiem et al.,2007).

    Within the region north of 80°N—hereinafter referred to as the central Arctic—the SIC has remained high this century,with multiyear ice.However,the sea-ice age in the central Arctic has clearly responded to global warming.Most multiyear ice of age greater than ten years has been replaced by much younger ice(Belchansky et al.,2005).The loss of the oldest ice is even more extreme,with ice of greater than five years reaching a minimum in 2010 of just 6%of the 1983–2002 mean(Maslanik et al.,2011),illustrating that central Arctic ice has also been affected by this warming.

    Here,we show the appearance of a record low ice concentration in the central Arctic(CARLIC)in summer 2010 based on field observations and satellite remote sensing.Large areas of open water appeared in high latitude areas,north of 85°N,resulting in a SIC that was actually lower than that of the surrounding pack ice at lower latitudes.This strikingly low ice concentration in the central Arctic is unique in our analysis period of 2003–15,and has not been previously reported in the literature.The remarkable opening of the ice appears not to have been produced by local melt,but rather by sea-ice divergence.The objectives of this paper are to present new evidence for an observed reduction in SIC at high latitudes in the central high Arctic basin,and to ascertain what role regional-scale climate processes played in the observed processes.

    Fig.1.Survey routes of Xue Long during its 2010 Arctic cruise.The blue line in(a)is the northward route and the red line the return route.(b)Navigation speed of the ship from 26 July to 6 August 2010,starting from Station S26(72°42.04′N,153°33.12′W)and ending at Station BN11(86°04.85′N,176°05.88′W).The red lines are latitudinally averaged navigation speed

    2.Record low ice concentration zone

    Fig.2.Distribution of SIC in the Arctic Ocean during August 2010 from AMSR-E data.The white circle around the pole is defined here as the “North Pole blind zone”,due to satellite orbital geometry.

    Very low ice concentration was first identified by the navigation speed of the Chinese R/V ship Xue Long.The ship started its northward journey from station S26(72°42.04′N,153°33.12′W)on 27 July along the blue line in Fig.1a.Ten days later,on 6 August,the ship arrived at Station BN11(86°04.85′N,176°05.88′W),from where the ship started ice camp measurements drifting with sea ice for 12 days.Ice conditions strongly influenced the navigation speed as the ship is an ice-strengthened ship,not an icebreaker.The ship traveled across a low SIC zone between 72.5°N and 75.6°N,with an average speed of 8.6 knots.Between 75.6°N and 82.6°N,the average navigation speed reduced to 6.5 knots,traveling in flat first-year ice of high concentration.From 82.6°N northward,the SIC decreased quickly and large areas of open water frequently appeared.The ship then traveled with an average navigation speed of more than 8.2 knots(Fig.1b).

    The spatial distribution of SIC could be determined from 89-GHz Advanced Microwave Scanning Radiometer-EOS(AMSR-E)daily satellite microwave data at 6.25-km resolution(http://www.iup.uni-bremen.de:8084/amsr/amsre.html).AMSR-E stopped working on 4 December 2011.Its successor—Advanced Microwave Scanning Radiometer 2(AMSR2)—started to provide data from 18 May 2012.The SIC data are retrieved using the ARTIST Sea Ice algorithm(Spreen et al.,2008).The effects of melt ponds,wet snow and atmospheric water vapor can degrade SIC estimates,but SIC from AMSR-E reliably reflects the relative spatial difference of SIC(Meier,2005).For convenience,we express SIC as a fraction of unity 0–1.

    As shown in Fig.2,a region with very low SIC centered at(83°N,180°W)appeared on 31 July.It extended to a larger area up to the orbital“North Pole blind zone”of satellite coverage in the following 15 days.Although the low SIC area disappeared in the satellite record between 18 and 25 August,it subsequently reappeared,with the SIC reaching its minimum on 6 September.

    The variation in low SIC in the central Arctic can be seen in the daily averaged SIC(ASIC),

    where C(x,y,t)is the SIC at each grid point,and S is the area of the zone circled between a latitude(here,it is taken as 85°N)and the blind zone of satellite coverage(about 88.25°N).The variation in ASIC during 2003–15 from 1 August to 30 September is plotted in Fig.3.It is shown that the ASIC was very high(>0.9)in the central Arctic in most years.The lower ASICs appeared in 2007,2010,2012 and 2013.The minimum ASICs of the central Arctic dropped to about 0.87(2007),0.85(2012)and 0.86(2013).However,in 2010,the ASIC dropped abruptly to about 0.78,though the ice coverage of the whole Arctic was more than that in 2007 and 2012.The very low SIC in the central Arctic appeared in late July and existed continuously until late September.Figure 3 shows that the ASIC in early September was even lower than that encountered by the ship in early August 2010.As of the publication of this article,the record low of SIC in 2010 in the central Arctic has not been broken.

    Here,we define a useful integral to express the interannual difference in SIC in the central Arctic:

    Fig.3.Interannual features of low SIC in central Arctic.Daily average SIC north of 85°N from 1 August to 30 September from 2003 to 2015.

    where T is equal to 61 days within August and September,and AOW2Mis a two-month averaged area of open water.Higher values of AOW2Mindicate a joint effect of lower SIC and/or longer duration of the low ASIC.The AOW2Mfor the central Arctic surrounded by 85°N in Fig.4 clearly shows that the averaged area of open water reached its maximum in 2010 because of a longer-lasting low SIC.The AOW2Min 2007 was the second lowest ASIC this century(Fig.4).In 2012 and 2013,the minimum ASICs were close to that in 2007,but the AOW2Min these years was much lower than that of 2007.

    Fig.4.Two-month averaged area of open water in the central Arctic surrounded by 85°N and the blind zone from 2003 to 2015.

    The two plots of SICs using the regional minimum days of 5 September 2007 and 6 September 2010 are plotted in Figs.5a and b,respectively,to compare their differences.It is clear that the low overall concentration in 2007 was caused by an extreme retreat of the ice edge;whereas,within the pack ice,the concentration remained high(Fig.5a).The low concentration in 2010 was different,as it formed as an opening within the pack ice(Fig.5b).It can be seen from Fig.5 that in 2010 not only did the SIC north of 85°N decline,but the overall SIC in the Atlantic sector was also significantly reduced.

    3.Discussion on the driving factors of CARLIC

    Based on optical measurements,Zhao et al.(2009)studied the sea-ice melt rate in the central Arctic.They found that only 0.33 cm d-1of sea ice could be melted by absorption of solar radiation,even with the strong solar insolation in August.The solar radiation penetrating open water in the summer is the main heat source to the ocean surface mixed layer,as widely addressed by previous studies(e.g.,Kadko and Swart,2004).Work by Perovich et al.(2008)showed the strong influence of heat in the ocean surface mixed layer on the reduction of sea ice in the Beaufort Sea,while questioning whether the receipt of surface radiation would be able to melt ice at higher latitudes.Repeated in-situ ice thickness measurements by an electromagnetic induction(EM31)were conducted during the 12-day ice camp in 2010 started from(86°04.85′N,176°05.88′W).The sea ice surrounding the ice camp was first-year ice,based on ice-core analysis(Lei et al.,2012).The measurements of ice thickness along four repeated profiles suggested an average melt rate of 2 cm d-1,primarily bottom melt(Xieetal.,2013).The melting rate was 2.5 times greater than the average~0.8 cm d-1during a similar period of the SHEBA experiments in 1998,even though SHEBA was at much lower latitudes(70°–80°N)(Perovich et al.,2003).The high melt rate occurred due to the heat absorption through the large area of open water.However,even with this rate,it would still need 100 days to melt a 2-m-thick layer of ice.Therefore,the sea-ice melt could not be solely responsible for the record low SIC at the end of July.Instead,the rapid opening and closing shift of the low SIC region suggests that the divergence/convergence of ice drift might be more relevant to the CARLIC.

    Fig.5.Distribution of ice concentration in the Arctic Ocean with the lowest daily average SIC north of 85°N in(a)2007 and(b)2010.

    The transfer of momentum from the atmosphere to the ice is critical to sea-ice drift(Martin and Gerdes,2007).Ignoring the geopotential gradient and nonlinear interaction,the equation of ice motion as a continuum is(Lepp?ranta,2005)

    can be obtained by taking the curl and divergence for both sides of Eq.(3):

    Then,the equation for the divergence of ice drift is

    Considering that the time scale of ice drifting is larger than the inertial period,the two-order temporal derivative is small,and the variation of the divergence of the stresses is negligible,through analyzing the order of magnitude,D can be approximately expressed by

    The equation of ice concentration is as follows(Hibler,1979):

    where SAis related to ice growth,and ε represents diffusion terms.Because the new ice formation is negligible in this season,and the diffusion and advection are both small,the relationship between the divergence of sea-ice drift and concentration becomes

    Replacing D with a spatially averaged version in the region north of 85°N,and substituting Eq.(8)into Eq.(1),we obtain

    where S is the area north of 85°N,and the left-hand side of Eq.(9)is the relative rate of variation of ASIC.Because the internal ice stress and the drag stress of water are unknown,ASIC(t)cannot be obtained from Eq.(9).The response of the variation in SIC depends on the SIC itself.When sea ice is dense,the sea ice responds weakly to the wind stress curl,as the internal stress σ arising from the interaction of different parts of ice floe balance most of the wind forcing.Otherwise,when the sea ice is sparse,the ice drift becomes more responsive to wind(Kwok et al.,2013;Olason and Notz,2014).The water drag stress for sea ice usually responds to ice drift.The curl of wind stress in this equation is the only forcing factor,and the other terms are response factors and are expected to respond to the wind in different ways.

    We calculated the daily averaged wind stress curl(AWSC)north of 85°N in August and September 2010 using daily wind velocity data from NCEP Reanalysis 1(Kalnay et al.,1996).The relative rate of variation of ASIC was calculated using the daily ASIC.The total correlation coefficient between ASIC and AWSC was-0.54,at much higher than the 99.5%confidence level.It can be seen from Fig.6a that ASIC responds well to each event with high averaged wind stress curl.It verifies that the wind stress curl is one of the most important factors in producing CARLIC.

    Although the relative rate of variation of ASIC correlated well with wind stress curl, the response of ASIC itself to wind stress curl was related to the degree of sea-ice concentrating,as shown in Fig.6b.From early August,a positive AWSC lasted for a couple of weeks,driving the ASIC decline from 0.97 to its first minimum of 0.89.Then,a seven-day negative AWSC acted on the area to cause a convergence of the ice-drift field,and the ASIC recovered to 0.96.After this,a three-week period dominated by positive AWSC occurred.The ASIC declined again and reached its minimum of 0.78 on 6 September.A rapid increase in the ASIC occurred again,responding to the strong negative trend of AWSC since 13 September,and all of the open water in the central Arctic closed over this period.The positive AWSC occurred once again after 20 September,but the ASIC did not respond to it anymore because freezing of sea ice had begun at these high latitudes.

    The AWSCs of August in recent years are plotted in Fig.7.Since 2007,the wind stress curl in the central Arctic was negative,except in 2010.This may explain why the CARLIC only occurred in 2010.During the sea-ice minimum in 2007,the amount of multiyear ice was reduced in the Pacific sector of the Arctic Ocean,and replaced with an increased areal extent of first-year ice(Barber et al.,2012).With positive wind stress curl,the ASIC in 2007 reached the second low record.In 2012 summer,the sea-ice cover of the Arctic reached its recorded minimum,but the ice concentration in the central Arctic was still high.This means that the wind stress curl did not drive the occurrence of low SIC in the central Arctic after 2010.

    Fig.6.Influence of wind stress curl on SIC.The green bars present the time series of negative averaged wind stress curl(AWSC)north of 85°N.(a)Correlation of AWSC and relative rate of variation of ASIC(blue dots).The correlation coefficient is-0.54.(b)Comparison of AWSC and ASIC(red dots).

    Fig.7.Averaged wind stress curl in August each year.

    However,in 2003 and 2006,positive wind stress curls with magnitudes much larger than those in 2007 and 2010 were found,but no CARLIC events occurred.Therefore,it seems that the wind forcing is not the only factor in generating CARLIC.A reasonable explanation is that the SIC in 2003 and 2006 was higher(see Fig.3),responding poorly to wind forcing.The response of the ice to wind forcing has recently become more pronounced,since in summer the Arctic is now dominated by first-year ice types,smaller floe sizes and decreased concentration(e.g.,Asplin et al.,2009),allowing for increased ice speeds within the gyres(Galley et al.,2013).

    Besides the regional wind forcing by positive wind stress curl,the ice drift patterns on the basin scale might also contribute to the CARLIC.The trajectories of the buoys from the International Arctic Buoy Program(IABP;Rigor,2002)from 1 January to 30 September 2010(Fig.8a)were clustered into two regions,separated by the dashed purple line:one group went toward the Fram Strait within the Transpolar Drift,and the other went eastward into the Beaufort Gyre.No buoy went across the dashed line in the first 10 months of 2010,which means that the ice flowing out of the Arctic was partly compensated for by export from the Laptev Sea and western sector.The lack of full compensation favored a low ice concentration.

    The averaged SLPfi eld of April to August 2010(Fig.8b)also matches the drift patterns identified from the IABP buoys.This pattern drives the production of a divergence in the central Arctic,which is quite similar to the double-gyre ice drift pattern reported by Wang and Zhao(2012).Using Polar Pathfinder monthly 25-km EASE-Grid Sea Ice Motion Vectors(Fowler,2008),Wang and Zhao(2012)divided the ice drift pattern into four main types:Transpolar Drift plus Beaufort Gyre(TPDBG;38%of total occurrence),anticyclonic(15%),cyclonic(16%),and double-gyre(15%).The TPDBG type is the typical ice drift pattern in the Arctic,with the highest occurrence.The double-gyre drift type is quite similar to the TPDBG type,except that the sea ice in the central Arctic drifts to the Canadian Archipelago,not to the Fram Strait—quite similar to the typical characteristic of the drifting pattern in 2010.Although the double-gyre drift pattern benefits divergence in the central Arctic,its occurrence usually lasts a short time.However,in 2010,the double-gyre pattern at the basin scale lasted more than 10 months,which might have been the immediate cause of the CARLIC.

    4.Summary and conclusions

    A record low concentration of sea ice and large area of open water in the central Arctic during summer 2010 is reported in this paper.The lowest averaged SIC north of 85°N reached as low as 0.78,becoming the sparsest than at any time in the historical record.We conclude that,in this particular case,the low SIC was caused more by ice divergence than by in-situ melt,based on the temporally resolved measurement of ice melt at a nearby ice camp.A high correlation between SIC and wind stress curl is revealed to address the contribution of regional wind forcing on the divergence of ice drift.The high correlation coefficient suggests that regional wind forcing might have been a key driving factor of the sea-ice drift in summer2010.However, in2003and2006,the magnitude of the wind stress curls were much larger than those in 2007 and 2010,but no CARLIC events occurred,because the heavy ice seems to have prevented the occurrence of low SIC,as the region was still dominated then by multiyear ice forms.The drift trajectories of ice buoys(IABP)depicted a divergent transpolar drift in the Atlantic sector and an eastward drift in the Pacific sector.This feature illustrates a double-gyre drift pattern,cyclonic in the Transpolar Drift and anticyclonic in the Beaufort Gyre,which resulted in ice divergence in the central Arctic.In 2010,this drift pattern resulted in high concentrations of sea ice at lower latitudes,which decreased within the high latitudes toward the pole(Fig.1).This pattern was also observed in 2009(Barber et al.,2009),and again in 2012(Babb et al.,2013),while the occurrence in 2010 lasted more than 10 months and the persistent divergence drove CARLIC.

    Fig.8.Ice-drift pattern in 2010.(a)Drift trajectories of ice buoys in the Arctic Ocean from 1 January to 30 September 2010.Red dots are the start positions of each buoy.The dashed purple line is the division of the two drifting groups.(b)Averaged SLP of April–August 2010.

    An important question relates to the frequency of future occurrences of CARLIC-type events in the central Arctic,which is difficult to project given its rare occurrence and our limited knowledge about their formation and maintenance.Based on our results,however,a long-lasting positive wind stress curl favors the occurrence of CARLIC,which might occur again when such a wind condition reappears.This divergence of ice drift in the central Arctic might be a significant feature of sea-ice rapid decline at high latitudes in the future,due to the preconditioning that the open water would have on increasing melt from an enhanced ocean surface mixed layer temperature.The physical significance of CARLIC is that more solar energy can penetrate into open water,which can in turn enhance the ice melt and feed back to the atmosphere(Vihma,2014).The large area of open water in the ice pack potentially has substantial biological implications as well.Further investigation is necessary to reveal the climatic significance of this double-gyre pattern,and its coupling to sea-ice motion and melt.

    Acknowledgements.This study was funded by the Global Change Research Program of China(Grant No.2015CB953900),the Key Program of the National Natural Science Foundation of China(Grant Nos.41330960 and 41406208),the Canada Research Chairs Program,NSERC,and Canadian Federal IPY Program Office.This paper is a contribution to the Arctic Science Partnership(ASP)and ArcticNet research networks.

    Asplin,M.G.,J.V.Lukovich,and D.G.Barber,2009:Atmospheric forcing of the Beaufort Sea Ice Gyre:Surface pressure climatology and sea ice motion.J.Geophys.Res.,114,C00A06,https://doi.org/10.1029/2008JC005127.

    Babb,D.,R.J.Galley,M.G.Asplin,J.V.Lukovich,and D.G.Barber,2013:Multiyear sea ice export through the Bering Strait during winter 201112.J.Geophys.Res.,118,5489-5503,https://doi.org/10.1002/jgrc.20383.

    Barber,D.G.,and Coauthors,2009:Perennial pack ice in the southern Beaufort Sea was not as it appeared in the summer of 2009.Geophys.Res.Lett.,36,L24501,https://doi.org/10.1029/2009GL041434.

    Barber,D.G.,and Coauthors,2012:Change and variability in sea ice during the 2007-2008 Canadian International Polar Year Program.Climatic Change,115,115–133,https://doi.org/10.1007/s10584-012-0477-6.

    Belchansky,G.I.,D.C.Douglas,and N.G.Platonov,2005:Spatial and temporal variations in the age structure of Arctic sea ice.Geophys.Res.Lett.,32,L18504,https://doi.org/10.1029/2005GL023976.

    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.

    Fowler,C.,2008:Polar Pathfinder Daily 25 km EASE-Grid Sea Ice Motion Vectors.National Snow and Ice Data Center,Boulder,Colorado USA.Digital Media.

    Galley,R.J.,B.G.T.Else,S.J.Prinsenberg,and D.G.Barber,2013:Summer sea ice concentration,motion,and thickness near areas of proposed offshore oil and gas development in the Canadian Beaufort Sea-2009.Arctic,66(1),105–116.

    Hibler III,W.D.,1979:A dynamic thermodynamic sea ice model.J.Phys.Oceanogr.,9,815–846,https://doi.org/10.1175/1520-0485(1979)009<0815:ADTSIM>2.0.CO;2.

    Holland,M.M.,C.M.Bitz,E.C.Hunke,W.H.Lipscomb,and J.L.Schramm,2006:Influence of the sea ice thickness distribution on polar climate in CCSM3.J.Climate,19(11),2398–2414,https://doi.org/10.1175/JCLI3751.1.

    Kadko D.,and P.Swart,2004:The source of the high heat and freshwater content of the upper ocean at the SHEBA site in the Beaufort Sea in 1997.J.Geophys.Res.,109,C01022,https://doi.org/10.1029/2002JC001734.

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

    Kay,J.E.,M.M.Holland,and A.Jahn,2011:Inter-annual to multi-decadal Arctic sea ice extent trends in a warming world.Geophys.Res.Lett.,38,L15708,https://doi.org/10.1029/2011GL048008.

    Kinnard,C.,C.M.Zdanowicz,D.A.Fisher,E.Isaksson,A.de Vernal,and L.G.Thompson,2011:Reconstructed changes in Arctic sea ice over the past 1,450 years.Nature,479,509–512,Https://doi.org/10.1038/nature10581.

    Kumar,A.,and Coauthors,2010:Contribution of sea ice loss to Arctic amplification.Geophys.Res.Lett.,37,L21701,https://doi.org/10.1029/2010GL045022.

    Kwok,R.,G.Spreen,and S.Pang,2013:Arctic sea ice circulation and drift speed:Decadal trends and ocean currents.J.Geophys.Res.,118,2408–2425,https://doi.org/10.1002/jgrc.20191.

    Lei,R.B.,Z.H.Zhang,I.Matero,B.Cheng,Q.Li,and W.F.Huang,2012:Reflection and transmission of irradiance by snow and sea ice in the central Arctic Ocean in summer 2010.Polar Research,31,17325,https://doi.org/10.3402/polar.v31i0.17325.

    Lemke,P.,and Coauthors,2007:Observations:Changes in snow,ice and frozen ground.Climate Change 2007:The Physical Science Basis.Contribution of Working Groupito the Fourth Assessment Report of the Intergovernmental Panel on Climate Change,Soloon et al.,Eds.Cambridge University Press,Cambridge,United Kingdom and New York,USA,339–383.

    Lepp?ranta,M.,2005:The Drift of Sea Ice.Springer-Verlag,266 pp.

    Lindsay,R.W.,and J.Zhang,2005:The thinning of Arctic sea ice,1988-2003:have we passed a tipping point?J.Climate,18(22),4879–4894,https://doi.org/10.1175/JCLI3587.1.

    Martin,T.,and R.Gerdes,2007:Sea ice drift variability in Arctic Ocean Model Inter comparison Project models and observations.J.Geophys.Res.,112,C04S10,https://doi.org/10.1029/2006JC003617.

    Maslanik,J.A.,J.Stroeve,C.Fowler,and W.Emery,2011:Distribution and trends in Arctic sea ice age through spring 2011.Geophys.Res.Lett.,38,L13502,https://doi.org/10.1029/2011GL047735.

    Meier,W.N.,2005:Comparison of passive microwave ice concentration algorithm retrievals with AVHRR imagery in Arctic peripheral seas.IEEE Transactions on Geoscience and Remote Sensing,43(7),1324–1337,https://doi.org/10.1109/TGRS.2005.846151.

    Nghiem,S.V.,I.G.Rigor,D.K.Perovich,P.Clemente-Col′on,J.W.Weatherly,and G.Neumann,2007:Rapid reduction of Arctic perennial sea ice.Geophys.Res.Lett.,34,L19504,https://doi.org/10.1029/2007GL031138.

    Olason,E.,and D.Notz,2014:Drivers of variability in Arctic sea-ice drift speed.J.Geophys.Res.,119,5755–5775,https://doi.org/10.1002/2014JC009897.

    Overland,J.E.,and M.Y.Wang,2007:Future regional Arctic sea ice declines.Geophys.Res.Lett.,34,L17705,https://doi.org/10.1029/2007GL030808.

    Overland,J.E.,and M.Y.Wang,2013:When will the summer Arctic be nearly sea ice free?Geophys.Res.Lett.,40,2097–2101,https://doi.org/10.1002/grl.50316.

    Perovich,D.K.,T.C.Grenfell,J.A.Richter-Menge,B.Light,W.B.Tucker III,and H.Eicken,2003:Thin and thinner:Sea ice mass balance measurements during SHEBA.J.Geophys.Res.,108(C3),8050,https://doi.org/10.1029/2001JC001079.

    Perovich,D.K.,J.A.Richter-Menge,K.F.Jones,and B.Light,2008:Sunlight,water,and ice:Extreme Arctic sea ice melt during the summer of 2007.Geophys.Res.Lett.,35,L11501,https://doi.org/10.1029/2008GL034007.

    Raddatz,R.L.,R.J.Galley,L.M.Candlish,M.G.Asplin,and D.G.Barber,2013:Integral profile estimates of latent heat flux under clear skies at an unconsolidated sea-ice surface.Atmos.-Ocean.,51(3),239–248.

    Rigor,I.,2002:IABP drifting buoy pressure,temperature,position,and interpolated ice velocity.Compiled by the Polar Science Center,Applied Physics Laboratory,University of Washington,Seattle,in association with NSIDC.National Snow and Ice Data Center,Boulder,CO,https://dx.doi.org/10.7265/N53X84K7.

    Spreen,G.,L.Kaleschke,and G.Heygster,2008:Sea ice remote sensing using AMSR-E 89-GHz channels.J.Geophys.Res.,113,C02S03,https://doi.org/10.1029/2005JC003384.

    Stroeve,J.,M.M.Holland,W.Meier,T.Scambos,and M.Serreze,2007:Arctic sea ice decline:Faster than forecast.Geophys.Res.Lett.,34,L09501,https://doi.org/10.1029/2007 GL029703.

    Tietsche,S.,D.Notz,J.H.Jungclaus,and J.Marotzke,2011:Recovery mechanisms of Arctic summer sea ice.Geophys.Res.Lett.,38,L02707,https://doi.org/10.1029/2010GL045698.

    Tschudi,M.A.,J.C.Stroeve,and J.S.Stewart,2016:Relating the age of Arctic Sea ice to its thickness,as measured during NASA’s ICESat and IceBridge campaigns.Remote Sensing,8,457,https://doi.org/10.3390/rs8060457.

    Vihma,T.,2014:Effects of Arctic sea ice decline on weather and climate:A review.Surveys in Geophysics,35,1175–1214,https://doi.org/10.1007/s10712-014-9284-0.

    Wang,X.Y.,and J.P.Zhao,2012:Seasonal and inter-annual variations of the primary types of the Arctic sea-ice drifting patterns.Advances in Polar Science,23(2):72–81,https://doi.org/10.3724/SP.J.1085.2012.00072.

    Xie,H.,R.Lei,C.Ke,H.Wang,Z.Li,J.Zhao,and S.F.Ackley,2013:Summer sea ice characteristics and morphology in the Pacific Arctic sector as observed during the CHINARE 2010 cruise.The Cryosphere,7,1057–1072,https://doi.org/10.5194/tc-7-1057-2013.

    Zhao,J.P,T.Li,S.G.Zhang,and Y.T.Jiao,2009:The shortwave solar radiation energy absorbed by packed sea ice in the central Arctic.Advances in Earth Science,24(1),34–42,https://doi.org/10.3321/j.issn:1001-8166.2009.01.004.(in Chinese with English ABSTRACT)

    24 March 2017;revised 4 August 2017;accepted 5 September 2017)

    :Zhao,J.P.,D.Barber,S.G.Zhang,Q.H.Yang,X.Y.Wang,and H.J.Xie,2018:Record low sea-ice concentration in the central Arctic during summer 2010.Adv.Atmos.Sci.,35(1),106–115,https://doi.org/10.1007/s00376-017-7066-6.

    ?Corresponding author:Jinping ZHAO

    Email:jpzhao@ouc.edu.cn

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

    日韩av免费高清视频| 91精品三级在线观看| 久久久久国产精品人妻一区二区| 女性被躁到高潮视频| 免费在线观看视频国产中文字幕亚洲 | 国产 一区精品| 如日韩欧美国产精品一区二区三区| 日本午夜av视频| 最黄视频免费看| 婷婷色综合大香蕉| 亚洲欧美成人综合另类久久久| 亚洲人成电影观看| 99国产综合亚洲精品| 亚洲,欧美精品.| 国产色婷婷99| 尾随美女入室| 91精品三级在线观看| 成年人免费黄色播放视频| 成人亚洲精品一区在线观看| 制服丝袜香蕉在线| 中文天堂在线官网| 高清黄色对白视频在线免费看| 日韩制服丝袜自拍偷拍| 亚洲图色成人| 精品亚洲成a人片在线观看| 国产野战对白在线观看| 久久精品熟女亚洲av麻豆精品| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品久久久久久久性| 999久久久国产精品视频| 久久久久久久大尺度免费视频| 女人被躁到高潮嗷嗷叫费观| 自线自在国产av| 亚洲伊人色综图| 极品少妇高潮喷水抽搐| 丝袜喷水一区| 91成人精品电影| 国产免费又黄又爽又色| 99re6热这里在线精品视频| 亚洲欧美成人精品一区二区| 久久毛片免费看一区二区三区| 精品少妇黑人巨大在线播放| 亚洲中文av在线| 永久网站在线| 午夜日本视频在线| 日本午夜av视频| 国产成人免费观看mmmm| √禁漫天堂资源中文www| 国产成人av激情在线播放| 丝袜美足系列| 五月天丁香电影| 少妇熟女欧美另类| 少妇 在线观看| 久久久精品区二区三区| 老女人水多毛片| 日韩在线高清观看一区二区三区| 欧美日韩综合久久久久久| 最近最新中文字幕大全免费视频 | 亚洲av免费高清在线观看| 人妻人人澡人人爽人人| 夜夜骑夜夜射夜夜干| 下体分泌物呈黄色| 另类精品久久| 久久久精品免费免费高清| 欧美老熟妇乱子伦牲交| 18禁国产床啪视频网站| 亚洲成人av在线免费| 中文字幕人妻丝袜制服| 成人毛片a级毛片在线播放| 亚洲欧美一区二区三区黑人 | 一本久久精品| 欧美成人午夜精品| 国产日韩欧美亚洲二区| 男女午夜视频在线观看| 一区二区日韩欧美中文字幕| 国产淫语在线视频| 日韩在线高清观看一区二区三区| videosex国产| 夫妻性生交免费视频一级片| 国产精品av久久久久免费| 亚洲av日韩在线播放| 国产熟女午夜一区二区三区| 91aial.com中文字幕在线观看| 亚洲激情五月婷婷啪啪| 欧美精品一区二区免费开放| 国产男人的电影天堂91| 亚洲情色 制服丝袜| 好男人视频免费观看在线| 99热全是精品| 天堂俺去俺来也www色官网| 老女人水多毛片| av视频免费观看在线观看| 国产精品不卡视频一区二区| 在线天堂中文资源库| 亚洲国产精品一区二区三区在线| 男的添女的下面高潮视频| 午夜av观看不卡| 菩萨蛮人人尽说江南好唐韦庄| 久久久久视频综合| 看免费av毛片| 在线观看www视频免费| 日韩制服丝袜自拍偷拍| 久久狼人影院| 国产精品.久久久| 纵有疾风起免费观看全集完整版| 男女下面插进去视频免费观看| 国产精品国产三级专区第一集| av免费在线看不卡| 日本av手机在线免费观看| av电影中文网址| 国产乱来视频区| 999久久久国产精品视频| 亚洲精品久久成人aⅴ小说| 九草在线视频观看| 天天操日日干夜夜撸| 日本vs欧美在线观看视频| 多毛熟女@视频| 亚洲欧美清纯卡通| 最近的中文字幕免费完整| 老汉色∧v一级毛片| 青春草视频在线免费观看| 老汉色av国产亚洲站长工具| 欧美日韩一区二区视频在线观看视频在线| 久久免费观看电影| 国产97色在线日韩免费| 国产日韩一区二区三区精品不卡| 99久国产av精品国产电影| 国产精品一区二区在线不卡| 亚洲成国产人片在线观看| 国产精品免费大片| av在线播放精品| 国产精品99久久99久久久不卡 | 飞空精品影院首页| 久久久久久久亚洲中文字幕| 亚洲欧美成人精品一区二区| 国产黄色免费在线视频| 一本色道久久久久久精品综合| 国产高清国产精品国产三级| 蜜桃国产av成人99| av.在线天堂| 国产精品 欧美亚洲| 国产精品一二三区在线看| 久久鲁丝午夜福利片| 桃花免费在线播放| 国产精品免费视频内射| 黄色配什么色好看| 亚洲国产av影院在线观看| 女性被躁到高潮视频| 熟女av电影| 日本wwww免费看| 欧美 亚洲 国产 日韩一| 久久这里有精品视频免费| 免费人妻精品一区二区三区视频| 久久 成人 亚洲| 老女人水多毛片| 91午夜精品亚洲一区二区三区| 看免费av毛片| 天天影视国产精品| 波多野结衣一区麻豆| 色播在线永久视频| 人妻 亚洲 视频| 日韩中字成人| 久久久久精品性色| 国产欧美亚洲国产| 日本猛色少妇xxxxx猛交久久| 另类精品久久| 欧美xxⅹ黑人| 综合色丁香网| 日韩视频在线欧美| 看非洲黑人一级黄片| 午夜福利影视在线免费观看| 国产成人午夜福利电影在线观看| 少妇人妻 视频| 亚洲国产色片| 成人免费观看视频高清| 一级毛片电影观看| 韩国精品一区二区三区| 这个男人来自地球电影免费观看 | 欧美成人午夜免费资源| 亚洲精品视频女| 91精品国产国语对白视频| 久久久久久久亚洲中文字幕| 哪个播放器可以免费观看大片| 丰满少妇做爰视频| 深夜精品福利| 各种免费的搞黄视频| 亚洲五月色婷婷综合| 九草在线视频观看| 黑人猛操日本美女一级片| 国产成人欧美| 高清黄色对白视频在线免费看| 国产成人av激情在线播放| 91成人精品电影| 亚洲精品久久久久久婷婷小说| 黄色毛片三级朝国网站| 国产精品国产av在线观看| 又大又黄又爽视频免费| 九色亚洲精品在线播放| 国产免费现黄频在线看| 国产片内射在线| 国产深夜福利视频在线观看| 免费人妻精品一区二区三区视频| 热99国产精品久久久久久7| 18禁国产床啪视频网站| 大码成人一级视频| 久久久久久久大尺度免费视频| 亚洲av日韩在线播放| 日本爱情动作片www.在线观看| 国产日韩欧美亚洲二区| 国产精品一国产av| 久久久久网色| 性高湖久久久久久久久免费观看| 搡老乐熟女国产| 国产在线一区二区三区精| 美女大奶头黄色视频| 啦啦啦中文免费视频观看日本| 美女午夜性视频免费| 欧美人与性动交α欧美精品济南到 | 久久精品国产鲁丝片午夜精品| 男女边摸边吃奶| 亚洲国产欧美网| 一区二区日韩欧美中文字幕| 一级毛片 在线播放| 国产日韩欧美视频二区| 免费观看在线日韩| 毛片一级片免费看久久久久| 蜜桃在线观看..| 女人精品久久久久毛片| 女的被弄到高潮叫床怎么办| 满18在线观看网站| 毛片一级片免费看久久久久| 欧美人与善性xxx| 国产精品人妻久久久影院| 久久精品久久久久久久性| 欧美xxⅹ黑人| 王馨瑶露胸无遮挡在线观看| 久久精品国产亚洲av涩爱| 午夜激情久久久久久久| 99久久人妻综合| 国产黄色免费在线视频| 五月天丁香电影| 女性生殖器流出的白浆| 国产日韩一区二区三区精品不卡| 成年动漫av网址| 亚洲国产av影院在线观看| 日韩一区二区三区影片| 边亲边吃奶的免费视频| 嫩草影院入口| 亚洲久久久国产精品| 成年女人毛片免费观看观看9 | 中文字幕av电影在线播放| 久久婷婷青草| 女人精品久久久久毛片| 精品99又大又爽又粗少妇毛片| 欧美精品人与动牲交sv欧美| 美国免费a级毛片| 日本黄色日本黄色录像| 国产亚洲最大av| 女人高潮潮喷娇喘18禁视频| 一区在线观看完整版| 国产黄色视频一区二区在线观看| 你懂的网址亚洲精品在线观看| 综合色丁香网| 亚洲精品一二三| 亚洲成国产人片在线观看| 女的被弄到高潮叫床怎么办| 两性夫妻黄色片| 国产精品久久久久久精品电影小说| 久久女婷五月综合色啪小说| 日日摸夜夜添夜夜爱| 免费看不卡的av| 日韩大片免费观看网站| 中文字幕制服av| 天天躁夜夜躁狠狠久久av| 日韩免费高清中文字幕av| 91精品三级在线观看| 夫妻午夜视频| 巨乳人妻的诱惑在线观看| 在线看a的网站| 啦啦啦视频在线资源免费观看| 赤兔流量卡办理| 国产一区亚洲一区在线观看| 精品99又大又爽又粗少妇毛片| 日产精品乱码卡一卡2卡三| 欧美 亚洲 国产 日韩一| 女人久久www免费人成看片| 亚洲av电影在线观看一区二区三区| 日韩视频在线欧美| av片东京热男人的天堂| 18禁观看日本| 美女主播在线视频| 999精品在线视频| 亚洲av国产av综合av卡| 亚洲国产色片| 在线看a的网站| 美国免费a级毛片| 日韩熟女老妇一区二区性免费视频| 久久ye,这里只有精品| 青春草视频在线免费观看| 国产成人免费观看mmmm| 亚洲精品一区蜜桃| 最新中文字幕久久久久| 麻豆乱淫一区二区| 老鸭窝网址在线观看| 国语对白做爰xxxⅹ性视频网站| 国产精品二区激情视频| 丰满迷人的少妇在线观看| 精品一区二区三卡| 国产一区亚洲一区在线观看| 国产精品久久久av美女十八| 18禁裸乳无遮挡动漫免费视频| 中文字幕制服av| 亚洲美女黄色视频免费看| 欧美日韩精品网址| 成人毛片60女人毛片免费| 9191精品国产免费久久| 老司机影院毛片| 久久99热这里只频精品6学生| 欧美变态另类bdsm刘玥| 国产亚洲午夜精品一区二区久久| 亚洲国产看品久久| 在线观看三级黄色| 美女中出高潮动态图| 爱豆传媒免费全集在线观看| 欧美精品高潮呻吟av久久| 国产精品秋霞免费鲁丝片| 久久久久久久大尺度免费视频| 成人手机av| 国产乱来视频区| 日韩av在线免费看完整版不卡| 日本欧美国产在线视频| 欧美日韩综合久久久久久| av卡一久久| 国产在线视频一区二区| 啦啦啦在线免费观看视频4| 国产毛片在线视频| 黄色怎么调成土黄色| 人妻一区二区av| 91久久精品国产一区二区三区| 国产色婷婷99| 日本猛色少妇xxxxx猛交久久| 国产免费又黄又爽又色| 一边亲一边摸免费视频| 日韩精品免费视频一区二区三区| 男女边吃奶边做爰视频| 秋霞伦理黄片| 久久ye,这里只有精品| 两个人免费观看高清视频| 欧美日韩亚洲国产一区二区在线观看 | 欧美国产精品一级二级三级| 精品人妻在线不人妻| 在线观看三级黄色| h视频一区二区三区| 久久久国产一区二区| 美女午夜性视频免费| 精品午夜福利在线看| 精品亚洲成国产av| 桃花免费在线播放| 在线观看一区二区三区激情| av又黄又爽大尺度在线免费看| 国产一区有黄有色的免费视频| 国产精品亚洲av一区麻豆 | 国产av码专区亚洲av| 咕卡用的链子| 精品国产乱码久久久久久男人| 国产在视频线精品| 高清黄色对白视频在线免费看| 热99国产精品久久久久久7| 亚洲第一区二区三区不卡| 亚洲欧美清纯卡通| 纵有疾风起免费观看全集完整版| 国产精品免费视频内射| 在线观看免费日韩欧美大片| 妹子高潮喷水视频| 69精品国产乱码久久久| 中文字幕亚洲精品专区| 日韩精品有码人妻一区| 精品人妻一区二区三区麻豆| 亚洲一级一片aⅴ在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 国产男人的电影天堂91| 久久狼人影院| 亚洲精品国产色婷婷电影| 蜜桃在线观看..| 丰满迷人的少妇在线观看| 亚洲,欧美,日韩| 最近中文字幕高清免费大全6| 久久鲁丝午夜福利片| 高清不卡的av网站| 午夜福利在线免费观看网站| 亚洲欧美清纯卡通| 伊人亚洲综合成人网| 国产 精品1| 母亲3免费完整高清在线观看 | 欧美黄色片欧美黄色片| 精品午夜福利在线看| 99热国产这里只有精品6| 久久午夜综合久久蜜桃| 国产乱来视频区| 亚洲人成77777在线视频| 黄片播放在线免费| 97在线人人人人妻| 18在线观看网站| 午夜91福利影院| 中文字幕人妻丝袜一区二区 | 中文字幕人妻熟女乱码| 亚洲欧美清纯卡通| 捣出白浆h1v1| 菩萨蛮人人尽说江南好唐韦庄| 青春草国产在线视频| 熟女电影av网| 久久热在线av| 91成人精品电影| 成人黄色视频免费在线看| 老司机影院成人| 女性生殖器流出的白浆| 国产视频首页在线观看| tube8黄色片| 国产欧美日韩综合在线一区二区| 亚洲四区av| 久久女婷五月综合色啪小说| 人人澡人人妻人| 在线天堂中文资源库| 99久久综合免费| 中文字幕人妻丝袜制服| 国产精品免费视频内射| 哪个播放器可以免费观看大片| 自线自在国产av| 日韩大片免费观看网站| 午夜福利一区二区在线看| 久久久久久久久免费视频了| 色94色欧美一区二区| 亚洲国产最新在线播放| 成人毛片a级毛片在线播放| 韩国av在线不卡| 国产精品无大码| 欧美成人午夜免费资源| 国产男人的电影天堂91| 亚洲欧美成人精品一区二区| 久久国产亚洲av麻豆专区| 国产精品av久久久久免费| 亚洲一区中文字幕在线| 免费高清在线观看视频在线观看| 国产精品国产三级国产专区5o| 久久久亚洲精品成人影院| 热re99久久国产66热| 国产淫语在线视频| 成人国语在线视频| 日韩制服丝袜自拍偷拍| 欧美精品一区二区免费开放| 国产精品久久久久久精品电影小说| 另类亚洲欧美激情| 午夜免费男女啪啪视频观看| 亚洲欧洲国产日韩| 人妻少妇偷人精品九色| 欧美最新免费一区二区三区| 亚洲精品av麻豆狂野| 人人妻人人爽人人添夜夜欢视频| 麻豆精品久久久久久蜜桃| 久久鲁丝午夜福利片| 久久精品国产综合久久久| 高清黄色对白视频在线免费看| 欧美成人午夜精品| 人人妻人人澡人人看| 国产亚洲av片在线观看秒播厂| xxx大片免费视频| 亚洲视频免费观看视频| av卡一久久| 亚洲av日韩在线播放| 日韩一区二区三区影片| 国产麻豆69| 久久精品熟女亚洲av麻豆精品| 欧美精品av麻豆av| 国产成人精品福利久久| 免费在线观看视频国产中文字幕亚洲 | 熟女少妇亚洲综合色aaa.| 午夜福利,免费看| av卡一久久| 国产成人精品无人区| 国产亚洲一区二区精品| 日日撸夜夜添| 亚洲精华国产精华液的使用体验| 成人影院久久| 亚洲av电影在线进入| 老司机影院成人| 久久99蜜桃精品久久| 久久人人97超碰香蕉20202| 爱豆传媒免费全集在线观看| 国产伦理片在线播放av一区| 精品少妇内射三级| 精品人妻偷拍中文字幕| 国产爽快片一区二区三区| 777米奇影视久久| 亚洲精华国产精华液的使用体验| 一二三四中文在线观看免费高清| 成年美女黄网站色视频大全免费| 欧美成人午夜精品| 美女国产视频在线观看| 成人午夜精彩视频在线观看| 亚洲色图综合在线观看| 满18在线观看网站| 亚洲国产欧美日韩在线播放| 777米奇影视久久| 欧美激情 高清一区二区三区| 大片免费播放器 马上看| 交换朋友夫妻互换小说| 亚洲国产成人一精品久久久| 可以免费在线观看a视频的电影网站 | 在线天堂中文资源库| 国产精品久久久久成人av| 亚洲精品日韩在线中文字幕| 亚洲伊人色综图| 丰满乱子伦码专区| 午夜福利在线免费观看网站| 亚洲精品第二区| 久久久精品国产亚洲av高清涩受| 人妻 亚洲 视频| 在线观看免费高清a一片| 国产免费福利视频在线观看| 欧美亚洲日本最大视频资源| 精品午夜福利在线看| 亚洲成av片中文字幕在线观看 | 免费大片黄手机在线观看| 夜夜骑夜夜射夜夜干| 在线观看免费日韩欧美大片| 午夜福利视频在线观看免费| 青青草视频在线视频观看| 亚洲成人av在线免费| 秋霞在线观看毛片| 日日撸夜夜添| 午夜福利,免费看| 一本大道久久a久久精品| 久久99蜜桃精品久久| 免费黄频网站在线观看国产| 在线观看国产h片| 免费高清在线观看日韩| 日韩一区二区三区影片| 中国三级夫妇交换| 欧美黄色片欧美黄色片| 2021少妇久久久久久久久久久| 久久99精品国语久久久| 少妇精品久久久久久久| 青青草视频在线视频观看| 亚洲第一区二区三区不卡| 成人二区视频| av又黄又爽大尺度在线免费看| 99国产精品免费福利视频| 亚洲欧美清纯卡通| 国产 一区精品| 亚洲三级黄色毛片| 久久国产亚洲av麻豆专区| a级毛片在线看网站| 99久久人妻综合| 制服人妻中文乱码| 黄频高清免费视频| 欧美 日韩 精品 国产| 日本猛色少妇xxxxx猛交久久| 成人免费观看视频高清| 一级黄片播放器| 99久久中文字幕三级久久日本| 伦精品一区二区三区| 在现免费观看毛片| 久久精品久久久久久噜噜老黄| 考比视频在线观看| 精品卡一卡二卡四卡免费| 视频在线观看一区二区三区| 亚洲国产精品国产精品| 搡老乐熟女国产| 肉色欧美久久久久久久蜜桃| 免费看av在线观看网站| 一边亲一边摸免费视频| 国产精品一区二区在线观看99| 97精品久久久久久久久久精品| 亚洲av.av天堂| 在线精品无人区一区二区三| 在现免费观看毛片| 视频区图区小说| 两个人免费观看高清视频| 亚洲国产最新在线播放| 国产精品一国产av| 一区二区日韩欧美中文字幕| 日韩欧美精品免费久久| 波多野结衣一区麻豆| 午夜免费观看性视频| 咕卡用的链子| av免费观看日本| 人妻 亚洲 视频| 国语对白做爰xxxⅹ性视频网站| 欧美日韩视频精品一区| 午夜免费观看性视频| 一级毛片我不卡| 日本午夜av视频| av视频免费观看在线观看| 亚洲精品国产av成人精品| 亚洲色图 男人天堂 中文字幕| 亚洲精品久久成人aⅴ小说| 日韩欧美精品免费久久| av免费观看日本| 国产午夜精品一二区理论片| 又黄又粗又硬又大视频| 波多野结衣一区麻豆| 赤兔流量卡办理| av不卡在线播放| 搡女人真爽免费视频火全软件| 午夜福利,免费看| a 毛片基地| 99久久中文字幕三级久久日本| 亚洲国产日韩一区二区| 一级黄片播放器| 自拍欧美九色日韩亚洲蝌蚪91| 伦理电影大哥的女人| 久久久精品区二区三区| 欧美中文综合在线视频| 日本免费在线观看一区| 欧美日韩亚洲高清精品|