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

    Thinner Sea Ice Contribution to the Remarkable Polynya Formation North of Greenland in August 2018

    2021-07-26 14:38:26XiaoyiSHENChangQingKEBinCHENGWentaoXIAMengmengLIXueningYUandHailiLI
    Advances in Atmospheric Sciences 2021年9期

    Xiaoyi SHEN ,Chang-Qing KE ,Bin CHENG ,Wentao XIA ,Mengmeng LI ,Xuening YU ,and Haili LI

    1School of Geography and Ocean Science,Nanjing University,Nanjing 210023,China

    2Jiangsu Provincial Key Laboratory of Geographic Information Science and Technology,Nanjing University,Nanjing 210023,China

    3Collaborative Innovation Center of South China Sea Studies,Nanjing University,Nanjing 210023,China

    4Collaborative Innovation Center of Novel Software Technology and Industrialization,Nanjing 210023,China

    5Finnish Meteorological Institute,Helsinki 00101,Finland

    ABSTRACT In August 2018,a remarkable polynya was observed off the north coast of Greenland,a perennial ice zone where thick sea ice cover persists.In order to investigate the formation process of this polynya,satellite observations,a coupled iceocean model,ocean profiling data,and atmosphere reanalysis data were applied.We found that the thinnest sea ice cover in August since 1978 (mean value of 1.1 m,compared to the average value of 2.8 m during 1978?2017) and the modest southerly wind caused by a positive North Atlantic Oscillation (mean value of 0.82,compared to the climatological value of?0.02) were responsible for the formation and maintenance of this polynya.The opening mechanism of this polynya differs from the one formed in February 2018 in the same area caused by persistent anomalously high wind.Sea ice drift patterns have become more responsive to the atmospheric forcing due to thinning of sea ice cover in this region.

    Key words:polynya,sea ice thickness,wind,sea ice drift,Greenland

    1.Introduction

    The term “polynya” refers to a large opening among the sea ice pack that remains ice-free or is only covered by thin ice for a time period (a minimum timespan of~15 days).Generally,there are two types of polynya,i.e.,latent heat polynya and sensible heat polynya.The former forms from persistent wind which drives the sea ice away from a barrier (e.g.,coast,fast ice and ice shelf);the latter is driven by thermodynamic factors,e.g.,the warm water upwelling which melts the sea ice cover.In the Arctic,polynyas usually appear around the Arctic shelf areas (e.g.,Siberian Shelf and Beaufort Shelf regions,

    Preu?er et al.,2016

    ),increase the heat extraction from the ocean to the atmosphere [e.g.,an increase of approximate 150 W min the heat extraction from the ocean is due to polynyas forming;

    Martin et al.(2004)

    ] and affect the sea ice production [e.g.70 kmper year;

    Martin et al.(2004)

    ].Thus,they can serve as a new proxy for weather and climate variations (

    April et al.,2019

    ).Sea ice north of Greenland is among the oldest multiyear ice (MYI) zones in the Arctic.In the past 30 years,a decline of more than 50% in MYI area has been observed in the Arctic (

    Kwok,2018

    ).In February 2018,extremely high wind (higher than 25 m s) led to the development of a polynya in the Wandel Sea.Such a pronounced polynya has never been observed before in this region at that time of year since 1979 (

    Moore et al.,2018

    ;

    Ludwig et al.,2019

    ).In August 2018,a polynya was observed again.It appeared north of Kap Morris Jesup,developed rapidly towards the Lincoln Sea (west of 40°W) and eventually extended along the north coast of Greenland (covering the area between 15°?45°W).The formation of polynyas in the Arctic is mostly determined by regional wind (

    Smith et al.,1990

    ),and it has been found that high wind played an essential role in the February polynya formation (

    Moore et al.,2018

    ).In summer,thermodynamic effects,e.g.,ocean-to-atmosphere heat fluxes(

    Minnett and Key,2007

    ),ocean interior salt fluxes (

    Martin et al.,2004

    ) and transport of upwelling oceanic heat (

    Campbell et al.,2019

    ),are also critical.However,a full explanation for the August polynya formation is still not established.Understanding the mechanism,including both dynamic and thermodynamic processes for sudden opening of summer polynyas is still a challenging task.In this study,we applied a set of comprehensive data to investigate 2018 August polynya formation,i.e.,to analyze the effects of atmosphere,sea ice and ocean on the polynya occurrence.Satellite and in situ observations were used to identify the evolution of this polynya.In situ observations of sea ice thickness (SIT) and ocean interior properties (e.g.,ocean temperature,salinity and density) are very limited to the north of Greenland (e.g.,ice mass balance buoy data,

    http://imb-crrel-dartmouth.org/archived-data

    ).Hence,SIT model data from the Pan-Arctic Ice-Ocean Modelling and Assimilation System (PIOMAS) were used to investigate the ice thickness.Ocean interior temperature and salinity during the polynya period were detected using weekly gridded ocean temperature and salinity profiling data (generated from European Union (E.U.) Copernicus Marine Service Information).European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 reanalysis data as well as a large-scale atmospheric circulation index (i.e.,North Atlantic Oscillation,NAO) were also applied here.

    2.Data and methods

    Sea ice concentration (SIC) data were obtained from the National Snow and Ice Data Center (NSIDC) Nimbus-7 Scanning Multichannel Microwave Radiometer (SMMR)and Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave Imager (SSM/I)?Special Sensor Microwave Imager/Sounder (SSMIS) passive microwave observations (

    Parkinson et al.,1996

    ;

    Cavalieri et al.,1997

    ),which were produced using a revised National Aeronautics and Space Administration (NASA) Team algorithm.Daily products at a grid cell size of 25 × 25 km for the period 1979?2018 were used.SIT data were taken from the PIOMAS model (

    Zhang and Rothrock,2003

    ),and a monthly product with a spatial resolution of 1° (latitudinal direction)× 2.67° (longitudinal direction) for the period 1978?2018 was used.NSIDC sea ice drift data (

    Tschudi et al.,2019

    )were used to calculate sea ice drift speed.Daily products with a spatial resolution of 25 × 25 km for the period 1979?2018 were used.Weekly gridded ocean temperature and salinity profiling data were obtained from the E.U.Copernicus Marine Service Information MULTIOBS_GLO_PHY_NRT_015_001 dataset.The weekly product with a spatial resolution of 0.25° × 0.25° was obtained for the period 2014?18.This dataset is a combination of satellite and in situ data [for details see

    Guinehut et al.(2012)

    and

    Mulet et al.(2012)

    ].

    Hourly 2-m air temperature (T2m),10-m level wind speed/direction and surface sensible heat flux data with a spatial resolution of 0.25° × 0.25° were taken from the ECMWF ERA5 reanalysis data.Daily mean values were calculated for analysis.The daily NAO index was obtained from the Climate Prediction Center of National Oceanic and Atmospheric Administration (NOAA) for the period 1972?2018.

    Meteorological data from the weather station at Kap Morris Jesup (83.65°N,33.37°W) were used to analyze the local weather conditions.We used observations of T2m and wind speed/direction.The time step of observed data was 6 hours,and each observation was averaged daily in the present study.

    Moderate Resolution Imaging Spectroradiometer(MODIS) true-color combination imagery from NASA Worldview application,part of the NASA Earth Observing System Data and Information System (EOSDIS) were used.These images are the combination of MODIS level-1B reflectance images with a spatial resolution of 250 m,using band 1 (approximately 0.62?0.67 μm wavelength),band 4 (approximately 0.54?0.57 μm wavelength) and band 3 (approximately 0.46?0.48 μm wavelength) as the red,green and blue channels,respectively.

    The region of interest (ROI) was defined as 55°W?5°E and 82°?85°N,in which the polynya occurred.Unless otherwise stated,average values calculated within the ROI for all gridded fields were used to analyze the ice,ocean and atmosphere conditions.We calculated the time-series on the original grid cells of the respective dataset and used a SIC threshold of 50% for a definition of the polynya region,i.e.,grid cells with SIC less than 50% were defined as polynya grid cells.The polynya extent was comprised of all grid cells below 50% SIC within the ROI.The polynya extent was defined as the surface area of all SIC grid cells below 50% (i.e.,not weighted by SIC).The average SIC and polynya extent were calculated for each day.

    3.Results and discussions

    3.1.Formation of polynya and local meteorological conditions

    Fig.1.Sea ice concentration evolution between late July and early September,north of Greenland in the Arctic Ocean.Black polygons indicate the region of interest where the polynya was observed.Black curves indicate the boundaries of the polynya.The location of the Kap Morris Jesup weather station is marked as *.The location of the Wandel Sea is marked as △.Sea ice concentration data are taken from NSIDC.

    On 29 July 2018,in the middle of the melting period,no polynyas or leads were observed (

    Fig.1a

    ).The average SIC in the ROI was almost 98% (

    Fig.2a

    ) and the T2m was about 0.28°C (

    Fig.2b

    ).On 5 August,a polynya developed in the Wandel Sea (

    Fig.1b

    ,the location is marked as triangle) with an extent of approximately 4.4 × 10km(

    Fig.2a

    ),accompanied by the continued above-freezing T2m since the beginning of August (

    Fig.2b

    ).On 13 August,this polynya was observed to cover areas west of 40°W (

    Fig.1c

    ) with an extent of 11.3 × 10kmand a T2m of 0.69°C (

    Fig.2b

    and

    2c

    ).The extent of this polynya subsequently expanded until it reached its maximum size on 18 August (about 17.5 × 10km,

    Fig.2a

    ).On 5 September,the SIC in the ROI had returned to climatological values(

    Fig.2a

    ) with a T2m of ?5.97°C (

    Fig.2b

    ).The daily-mean T2m and wind speed observed at Kap Morris Jesup weather station during the polynya period are shown in

    Fig.3

    .From 15 July to 31 August,the average T2m was 4.9°C with standard deviation of 4.4°C,and the maximum T2m was 12°C.From 1 September onward,the T2m dropped below freezing.The average wind speed was 3.6 m swhile the standard deviation was 1.9 m s,and the average wind direction observed at the Kap Morris Jesup weather station was 28° west of due south (offshore).The warm and moderate southerly wind are consistent with the formation of August polynya.

    3.2.Environment factors linked to polynya development

    To investigate the environmental factors linked to the polynya formation and development,we present a comprehensive analysis of the atmosphere,sea ice and ocean conditions during the polynya period.Although the polynya appeared suddenly on 3 August,to investigate the effect of various environment factors,we carried out the data analyses from the beginning of 2018.

    3.2.1.The role of the atmosphere

    Usually,the contribution of the atmosphere to the polynya formation can be attributed to two aspects.First,sea ice melts due to high air temperatures;second,sea ice drift occurs under regional wind (a combination of wind speed and direction drives the sea ice drift).Although the T2ms were high in summer,they barely reached the freezing point (

    Fig.2b

    ) and thus can only lessen the sea ice refreezing.Hence,air temperature is not enough to explain the polynya formation and development.It is expected that wind is dominant for the polynya initial formation and further development,while the high temperatures lessened refreezing and thus aided in maintaining the polynya development.

    3.2.1.1.Sea ice drift driven by wind

    Fig.2.(a) Daily sea ice concentration;(b) daily 2-m air temperature and (c) monthly sea ice thickness in the region of interest.Daily values through 2018 of sea ice concentration and 2-m air temperature,and monthly values of sea ice thickness are shown as solid black lines.The climatological values (1979?2018) are shown as the medians (grey lines),and the 25%?75% interquartile range is marked by grey shading.The polynya extent is shown as a dotted line in (a).It should be noted that the climatological values for SIC are also provided in (a),and that their values are quite close to 100%.The black horizontal line in (b) is the freezing point.The vertical black dashed lines in (a) and (b) indicate the polynya period.The black line in (c) shows the monthly sea ice thickness in 2018 in the region of interest,and the black dot denotes the monthly minimum sea ice thickness during 1979?2018.The sea ice concentration,air temperature and sea ice thickness are taken from NSIDC,ERA5 data and PIOMAS model,respectively.

    Winds can alter sea ice cover by causing sea ice drift,especially when the sea ice cover is quite thin.Here we examined the relationship between the polynya extent evolution and the cumulative wind-speed anomaly from the average value during the August polynya period (

    Fig.4a

    ).The cumulative wind-speed anomaly,which was calculated from the sum of the wind speed anomaly from the average value during the time period between 15 July and 13 September,can reveal the degree of the wind speed anomaly in a time period.A high correlation coefficient (0.68,p-value <0.01) indicates that the polynya expansion was associated with the local wind.On 3 August when the polynya started to open,a positive ocean?atmosphere sensible heat flux of 52 W m(

    Fig.4c

    ,vertical fluxes are positive downwards)associated with warm southerly wind (observed at the Kap Morris Jesup weather station) caused the heat exchange from the atmosphere to ocean.This heat exchange may have created more heat into the ocean,which would contribute to the maintenance of the polynya and further affect the water mass transformation and deep water formation(

    Moore and Pickart,2012

    ).Although local wind was found to be associated with the polynya opening,the wind speeds in the ROI during polynya period were moderate.According to the ERA5 data,the wind speeds during August 2018 (3.8 m s,

    Fig.5b

    )were mostly within the climatological ranges and at a relative minimum compared to other times of the year.

    Figure 5c

    shows the wind directions in August 2018 accompanied with the average wind direction during August 1979?2017.The average wind direction pattern is normally northward flow during August 1979?2017,and this is similar to that in August 2018.Both wind speeds and directions (

    Figs.5b

    and

    5c

    ) were comparable to these in previous years,hence regional wind is not the determining factor to cause the polynya opening.Instead,some other factors appear to have made sea ice drift stronger than expected,under the modest wind speeds.

    Fig.3.(a) 2-m air temperature and (b) wind speed from 15 July to 13 September 2018 observed at Kap Morris Jesup weather station.Black horizontal dashed line in (a) is the zerodegree temperature.Black vertical dashed lines indicate the polynya period.These data are sub-daily measurements and averaged daily here.

    Fig.4.Polynya extent,weather parameters and surface heat flux between 15 July and 13 September 2018.(a) Daily polynya extent (solid line) and cumulative wind-speed anomaly (dot line) based on (b);(b) wind speed anomaly from the average value at Kap Morris Jesup weather station and (c) The spatial average sensible heat flux (solid) and an average sensible heat flux during the polynya period (dashed) estimated from ERA5 data in the region of interest.The flux is positive downward.

    Fig.5.The atmospheric circulation index and wind speed/direction in the region of interest.(a) North Arctic Oscillation index;(b) 10 m wind speed;(c) wind direction in August 2018.For (a) and (b),daily values are shown by black lines,associated climatological values (1979?2018) are shown by the medians (grey lines)and the 25%?75% interquartile range are marked by grey shading.Black dashed lines in (a) and (b) indicate the polynya period.NAO index values are taken from NOAA,10 m wind speed and direction are taken from ERA5 data.The arrows in (c) only represents the direction,their length cannot be used to indicate the wind speed.

    The average August SIT in 2018 in the ROI was the minimum since 1979 (

    Fig.2c

    ),according to the PIOMAS model.The continued decrease in SIT during the past 40 years contributed to the minimum in 2018.The SIT anomaly may be responsible for the polynya opening as thinner sea ice tends to be easier to drift (

    Kwok et al.,2013

    ).In order to test this assumption off the north coast of Greenland,we examined the relationship between wind speeds and ice drift speeds in different periods with different SIT conditions (

    Table 1

    ).Correlations between wind speeds and ice drift speeds tend to be higher when ROI SIT decreases,which indicates that the sea ice drift has become more responsive to atmospheric forcing due to sea ice thinning in this area.In addition,it was also found that most of the MYI was replaced by first-year ice (FYI) due to the February polynya along the north coast of Greenland (

    Ludwig et al.,2019

    ),and this persisted throughout the winter.Hence,the unusual conditions during the winter of 2018 also played an important role in the August polynya formation,as thin FYI is easier to drift than thick MYI.

    Table 1.Sea ice thickness,wind speed and sea ice drift speed in the region of interest in different periods from August 1979?2018.Correlations between wind speeds and sea ice drift speeds are also shown in table,correlation coefficients that were significant at the 95% confidence level according to two-tailed Student’s t-tests are shown in bold.

    3.2.1.2.Spatial pattern of sea ice drift

    Sea ice drift can be caused by both wind and oceanic current.Due to the lack of the observations of oceanic current,it is impossible to estimate the quantitative contributions of oceanic current to the sea ice drift.

    Kwok et al.(2013)

    used a squared correlation coefficient to describe the fraction of the variance of sea ice drift which can be explained by the geostrophic wind.They found a correlation coefficient of 0.80 in summer from 1982?2009,i.e.,wind can explain all but about 36% of the variance of sea ice drift.As the sea icecover is trending towards being thinner and weaker,this coefficient near the ice edge is expected to be higher in the future (

    Kwok et al.,2013

    ).Hence,wind is an essential factor to affect the sea ice drift compared to the ocean currents.The spatial pattern of wind is usually inhomogeneous;thus,the sea ice drift is also not uniform.Regionally-averaged ice drift vectors may hide the details of the spatial pattern of sea ice drift.While the sea ice divergence (calculated from sea ice drift vectors) can reveal detailed sea ice drift patterns,the sea ice drift product is more uncertain in areas of low SIC (

    Sumata et al.,2014

    ).In the present study,time-series of MODIS images in August 2018 shown in

    Fig.6

    reveal the spatial pattern of sea ice drift,accompanied by the local wind vectors (red arrows,taken from ERA5 data).The sea ice drift vectors (blue arrows) were determined based on visual interpretation of MODIS image sequences.In general,quite consistent spatial patterns between the wind and sea ice drift fields can be seen in

    Fig.6

    .This also confirmed that wind played a more essential role in sea ice motion (

    Kwok et al.,2013

    ).On 11,14,19,21 and 28 August when northerly winds occurred,the polynya coverage shrank.On other days with a southerly wind,the polynya expanded.Thus,the changes in prevailing direction of the wind (and associated sea ice drift) were consistent with the polynya evolution.This indicates that wind is the direct factor for August polynya formation.However,as was mentioned before,neither the wind speed nor direction in August 2018 were anomalous,thus the sea ice drift became more responsive to the regional wind favoring polynya formation when sea ice was getting thinner.Hence,thinner sea ice cover played a relatively greater role in the polynya formation.

    3.2.1.3.The role of the large-scale atmosphere circulation

    Large-scale atmosphere records were also examined(hereafter,wind speed data were derived from the ERA5 reanalysis data,

    Fig.5

    ).During the polynya period,the climate records showed positive NAO phases (mean value of 0.82,compared to climatological value of ?0.02),which were accompanied by a modest wind speed (3.8 m s,comparing to climatological value of 4.2 m sand averaged value of 4.5±0.6 m sduring 1979?2017),with correlation coefficients of 0.51 (p-value <0.01) in the ROI.In addition,a positive phase of the NAO index in August is associated with a low pressure anomaly over Greenland and therefore southeasterly wind anomalies in the ROI (

    Fig.5c

    ).Although NAO showed positive phases during polynya periods,its associated wind speeds were not anomalous in comparison to the climatological values (

    Fig.5

    ).Thinner sea ice cover increased the responsiveness of sea ice drift to wind,which indicates that polynyas are more susceptible to forming in this region under normal atmospheric conditions (i.e.,the MYI area north of Greenland has become less resilient with respect to latent heat polynya formation).

    3.2.2.The role of the ocean

    The temperature and salinity mainly control ocean stratification under sea ice (e.g.,

    Campbell et al.,2019

    );hence,they can be used to reveal ocean interior property changes(e.g.,heat/salinity upwelling).The upwelling induced by wind stress and subsequent Ekman divergence/offshore transport in the upper layer may bring the sub-surface heat to the surface,which can contribute to maintaining the polynya.Ocean profiling data shows that mixed-layer temperatures during the August polynya period were about 0.1°C higher than climatological values (

    Fig.7a

    ),whereas the mixed-layer salinity and depth were within the climatological range (

    Figs.7b

    and

    7c

    ).Following the polynya opening,temperature dropped only to ?1.7°C from ?1.5°C and salinity remained at 33 psu at a depth of 50 m in the ocean interior based on ocean property profiling data (

    Fig.8

    ).Thus,there is no evidence for the occurrence of upward heat flux in the polynya period.Upwelling of intermediate Atlantic water would result in ocean warming and salinification.Oceanic heat played a weaker role in August polynya formation than wind forcing,and therefore this polynya is classified as a latent heat polynya.On the other hand,sea ice may be subject to windinduced wave impact and break off,especially under the influence of storm-generated waves (

    Kohout et al.,2014

    ).Although the significant wave height (SWH) values near the sea ice edge suggest there was potential to break up the sea ice cover north of Greenland due to the thinner thickness,no major storms were observed prior this polynya event,and the wind speed pattern was not unusual.We may conclude that the impact of ocean waves on polynya formation in this case was minor.The effect of waves on ice break-up largely depends on the SWH at the ice edge and the distance between the ice cover and ice edge (

    Squire and Moore,1980

    ).A better quantitative estimation of wave impacts on ice break up would require more accurate in situ long term observations.

    Fig.6.Time-series of MODIS true-color combination images and wind vectors in August 2018.The blue arrows indicate sea ice drift vectors and red arrows represent wind vectors.The wind vectors were derived from ERA5 data.The sea ice drift vectors can only be used to indicate direction.

    3.3.Comparison to the previous polynyas

    As was mentioned earlier,the main mechanism for previous polynya opening in the Arctic is regional wind,hence the timing,area and duration of a polynya are usually variable (

    Morales-Maqueda et al.,2004

    ).Arctic polynyas mostly appeared on confined locations near the coast,such as “fast ice edge” and ice shelf.Most previous studies focused on the winter polynyas to quantify the ice production rates and surface energy fluxes.In comparison to these winter polynyas,the August polynya extent is similar(

    Table 2

    ).Two studies focused on winter polynyas to the north of Greenland:1)

    Preu?er et al.(2019)

    found a polynya with an average extent of 2.5×10km,and 2) a polynya with a mean length of 5.7 km was found by

    Winsor and Bj?rk (2000)

    .Much larger polynyas are found in both February and August 2018,which is remarkable in comparison to previous reported polynyas along the north coast of Greenland (

    Table 2

    ).There are also some polynyas that have appeared in Arctic summer (e.g.,Wrangel Island polynya and North-East Water polynya);however,no studies reported the summer polynyas along the north coast of Greenland.Although some winter polynyas have been found to the north of Greenland,remarkable polynyas only occurred in February and August 2018.As these two polynyas appeared in different seasons,atmospheric and ice conditions are different,hence it is helpful to compare the different mechanisms,which will benefit further research and forecasting of polynyas north of Greenland.

    Fig.7.Mixed-layer properties in the region of interest in 2018.(a) Mixed-layer temperature;(b) Mixed-layer salinity and (c) Mixed-layer depth.Weekly values are shown by black lines.The 2014?18 means are shown by the medians (grey lines) and 25%?75% interquartile ranges are denoted by grey shading.Black dashed lines indicate the polynya period.All these data are taken from E.U.Copernicus Marine Service Information.

    The February polynya featured anomalously high T2m,and the wind speed reached an extremely high value(15.3 m s) during the polynya period.These results agree with those of previous studies (

    Moore et al.,2018

    ;

    Ludwig et al.,2019

    ).Upon comparing to the February event,some climate conditions were found to be different for the August polynya.On the one hand,NAO was positive before the February polynya event,it subsequently decreased during the polynya opening,and then increased during the refreezing.Although the NAO during the August polynya period featured a similar change pattern,all its values were positive.The NAO during the February polynya period was negative since it displayed the maximum extent on 26 February,and this situation was associated with a sudden stratospheric warming (

    Moore et al.,2018

    ).On the other hand,lower average wind speed was found during the August polynya(3.8 m scomparing to 6.1 m s,

    Fig.5b

    ).We have found that thinner sea ice increases its wind-driven drift (

    Table 1

    )and the wind did not have to be as strong as the storm conditions in February 2018 to induce the polynya.Hence,more polynyas are expected to occur in the future as the overall Arctic SIT is continuously decreasing (

    Kwok,2018

    ).In terms of the ocean,during the February polynya period,a stable ocean stratification within the climatological boundary was found (

    Fig.7

    ).This is similar to that during the August polynya period,as no evidence for the occurrence of upward heat flux was found.

    Table 2.Statistics of the Arctic polynyas in previous studies.Polynya time,study area and extent/length are included.

    In addition,the peak temperature during the event in February rose almost to the freezing point.The temperatures before and after the event,however,were significantly lower.The temperatures in August polynya were much warmer than those in February polynya (?0.6°C comparing to ?16.4°C,from ERA5 data,

    Fig.2b

    ).The August polynya occurred before the period of minimum sea ice extent when T2ms were partly above the freezing point,which means that temperatures prevented refreezing in August,while in February they only decreased refreezing.Usually,a narrow lead developed off the ice zone on the north coast of Greenland in summer,and seldom evolved into a polynya like the one that occurred in August 2018.This polynya appeared when SIT was low and thus its formation and closing may also be associated with the seasonal cycle,which is quite different from that in February 2018 when a polynya formed in winter with higher SIT.Thinner sea ice in August increased its wind-driven drift and hence modest wind (compared to the anomalously strong wind in February) caused the August polynya opening,hence it was responsible for the formation of the August polynya.

    Fig.8.Ocean property profiling in the region of interest in 2018.(a) potential temperature;(c) salinity.(b) and (d)provide vertically enlarged plots showing the ocean property profiling from the ocean surface to a depth of 100 m for(a) and (c),respectively.Mixed-layer depth values are shown by white lines.All these data are taken from E.U.Copernicus Marine Service Information.

    4.Conclusion

    In the present study,we conclude that two factors were mainly responsible for the August polynya development.First,the thinnest August sea ice cover since 1979 occurred in 2018,which was due to the continued decreases in SIT in the Arctic.Sea ice drift has become more responsive to the regional wind favoring polynya formation when sea ice cover is becoming thinner.It was found that most of the MYI was replaced by FYI due to the February polynya,a situation that persisted throughout the winter.Hence,the unusual conditions during the winter of 2018 may also haveplayed an important role as the FYI is easier to move.Second,moderate wind persisted through the polynya period accompanied by the positive NAO.The modest wind speed during August polynya was sufficient to cause sea ice divergence under the thinnest sea ice conditions since 1979,and the sea ice drift pattern was consistent with the regional wind field (southerly),thus wind had a direct effect on the August polynya formation.

    This study indicates that the August polynya was caused by the combination of these two factors.Furthermore,our analysis suggests that more polynyas will presumably occur in the Arctic as continuing thinning of sea ice is easier to be moved by wind.This will further influence the atmospheric-ocean heat exchange and the large-scale ice mass balance.

    Exact understanding and prediction of Arctic polynya formation requires improved knowledge of both the dynamic (e.g.,wind-induced sea ice drift) and thermodynamic (e.g.,upwelling heat) processes during the polynya evolution.Unfortunately,the detailed processes are not well known because of limitations imposed by the applied satellite or model data.Sea ice production was not assessed in the present study.Its effects on the fate (e.g.time of opening and extent) of the polynyas in the future are largely unknown.More in situ and satellite data are needed to investigate in detail August polynyas in future studies.

    .This work was supported by the National Key Research and Development Program of China (Grant No.2018YFC1407206),Academy of Finland (Grant No.317999),and European Union’s Horizon 2020 research and innovation programme (Grant No.727890-INTAROS).

    Data availability

    SIC and sea ice drift data were derived from NSIDC at

    https://doi.org/10.5067/8GQ8LZQVL0VL

    and at

    https://doi.org/10.5067/INAWUWO7QH7B

    ,respectively (last accessed date:20 July 2020);PIOMAS SIT data were obtained at

    http://psc.apl.washington.edu/zhang/IDAO/data_piomas.html

    (last accessed date:20 July 2020);ocean temperature and salinity profiling data were acquired using E.U.Copernicus Marine Service Information at

    http://marine.copernicus.eu/services-portfolio/access-to-products/?option=com_csw&view=details&product_id=MULTIOBS_GLO_PHY_NRT_015_001

    (last accessed date:20 July 2020);ERA5 reanalysis data were obtained at

    https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=overview

    (last accessed date:20 July 2020);NAO index data were derived from NOAA at

    http://www.cpc.ncep.noaa.gov

    (last accessed date:20 July 2020);meteorological data from the Kap Morris Jesup weather station can be downloaded at

    https://www.timeanddate.com/weather

    (last accessed date:20 July 2020).MODIS image was obtained from NASA worldview application at

    https://worldview.earthdata.nasa.gov

    (last accessed date:20 July 2020).This article is distributed under the terms of the Creative Commons Attribution License which permits any use,distribution,and reproduction in any medium,provided the original author(s) and the source are credited.This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (

    http://creativecommons.org/licenses/by/4.0/

    ),which permits unrestricted use,distribution,and reproduction in any medium,provided you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons license,and indicate if changes were made.

    黑人巨大精品欧美一区二区mp4| 国产成人av激情在线播放| 国产男女超爽视频在线观看| 三上悠亚av全集在线观看| 亚洲 欧美一区二区三区| 国产欧美日韩一区二区三| 久久久久精品国产欧美久久久| 热99久久久久精品小说推荐| 国产一区二区 视频在线| 亚洲精品美女久久久久99蜜臀| 丰满迷人的少妇在线观看| 人妻一区二区av| 怎么达到女性高潮| 高清视频免费观看一区二区| 国产高清videossex| 亚洲精品一二三| 丰满迷人的少妇在线观看| 国产欧美日韩精品亚洲av| 久久99热这里只频精品6学生| 深夜精品福利| 国产av国产精品国产| 国产精品一区二区在线观看99| 亚洲午夜理论影院| 俄罗斯特黄特色一大片| 欧美性长视频在线观看| 69av精品久久久久久 | 在线观看一区二区三区激情| 精品久久蜜臀av无| 精品视频人人做人人爽| av一本久久久久| 亚洲av欧美aⅴ国产| 我的亚洲天堂| 女性生殖器流出的白浆| 精品国产国语对白av| 天堂中文最新版在线下载| 亚洲色图av天堂| 黄片大片在线免费观看| 飞空精品影院首页| 国产精品电影一区二区三区 | 999久久久精品免费观看国产| 国产欧美日韩一区二区精品| 伊人久久大香线蕉亚洲五| 夜夜夜夜夜久久久久| 日韩精品免费视频一区二区三区| 国产精品久久电影中文字幕 | 黑人欧美特级aaaaaa片| 精品少妇黑人巨大在线播放| 99国产精品99久久久久| 中文字幕人妻丝袜制服| 国产精品一区二区精品视频观看| 中文欧美无线码| 欧美+亚洲+日韩+国产| 一本—道久久a久久精品蜜桃钙片| 日韩视频在线欧美| 国产高清国产精品国产三级| 日韩三级视频一区二区三区| 亚洲欧美一区二区三区久久| 男女下面插进去视频免费观看| 欧美成人午夜精品| 午夜精品久久久久久毛片777| 涩涩av久久男人的天堂| 免费久久久久久久精品成人欧美视频| 精品乱码久久久久久99久播| 精品少妇内射三级| 高清黄色对白视频在线免费看| 在线观看免费视频日本深夜| 18禁黄网站禁片午夜丰满| 亚洲第一av免费看| 青青草视频在线视频观看| 久久国产精品大桥未久av| 精品福利观看| 99精品久久久久人妻精品| 国产在线观看jvid| 精品第一国产精品| 国产有黄有色有爽视频| 国产精品亚洲一级av第二区| 国产精品久久久久久精品古装| 男男h啪啪无遮挡| 久久ye,这里只有精品| 桃花免费在线播放| 欧美乱妇无乱码| 久久狼人影院| 免费一级毛片在线播放高清视频 | 成年动漫av网址| 国产成人精品久久二区二区免费| 国产精品 国内视频| 成人国产一区最新在线观看| 丰满人妻熟妇乱又伦精品不卡| 成人亚洲精品一区在线观看| 91成年电影在线观看| 免费看十八禁软件| 日本vs欧美在线观看视频| 久久国产精品人妻蜜桃| av电影中文网址| 在线亚洲精品国产二区图片欧美| 色在线成人网| 亚洲精品一二三| 亚洲av电影在线进入| 国产亚洲精品久久久久5区| 国产av又大| 亚洲精品乱久久久久久| www.熟女人妻精品国产| 我的亚洲天堂| 国产精品电影一区二区三区 | 欧美日韩亚洲国产一区二区在线观看 | 国产亚洲欧美精品永久| 午夜激情久久久久久久| 国产男女内射视频| 日本av手机在线免费观看| 国产亚洲午夜精品一区二区久久| 91麻豆av在线| 精品欧美一区二区三区在线| 精品第一国产精品| 女性被躁到高潮视频| 搡老熟女国产l中国老女人| 中文亚洲av片在线观看爽 | 精品人妻在线不人妻| 精品国产一区二区三区久久久樱花| av福利片在线| 久久ye,这里只有精品| 成年版毛片免费区| 欧美在线黄色| 男女之事视频高清在线观看| 成人国产av品久久久| 国产av一区二区精品久久| 国产成人精品无人区| 视频区图区小说| 国产高清videossex| 91九色精品人成在线观看| 亚洲色图av天堂| 国产高清视频在线播放一区| 日本五十路高清| 国产一区有黄有色的免费视频| 国产伦理片在线播放av一区| 视频区欧美日本亚洲| 性色av乱码一区二区三区2| 国产精品成人在线| 欧美黄色淫秽网站| 高潮久久久久久久久久久不卡| a级毛片黄视频| 国产精品久久久久久精品古装| 久久精品熟女亚洲av麻豆精品| 视频区欧美日本亚洲| 中国美女看黄片| 国产成人免费观看mmmm| 最近最新中文字幕大全电影3 | 91精品国产国语对白视频| 欧美+亚洲+日韩+国产| 亚洲欧美一区二区三区黑人| 久久青草综合色| 亚洲精品在线观看二区| 一级,二级,三级黄色视频| 国产成人av激情在线播放| av线在线观看网站| 日韩大码丰满熟妇| 国产亚洲精品一区二区www | 国产欧美日韩一区二区精品| 国产老妇伦熟女老妇高清| 精品乱码久久久久久99久播| 狂野欧美激情性xxxx| 变态另类成人亚洲欧美熟女 | 一本大道久久a久久精品| 精品人妻在线不人妻| 一区在线观看完整版| 王馨瑶露胸无遮挡在线观看| 宅男免费午夜| 一本色道久久久久久精品综合| 九色亚洲精品在线播放| 亚洲性夜色夜夜综合| 午夜成年电影在线免费观看| 日韩欧美国产一区二区入口| 久久久久国内视频| av电影中文网址| 在线永久观看黄色视频| 女人精品久久久久毛片| 不卡一级毛片| 大型av网站在线播放| 黑人欧美特级aaaaaa片| 国产又色又爽无遮挡免费看| 亚洲国产av影院在线观看| 在线 av 中文字幕| 久久精品人人爽人人爽视色| 欧美精品av麻豆av| 91成人精品电影| 最近最新免费中文字幕在线| 搡老熟女国产l中国老女人| 国产精品久久久人人做人人爽| 欧美乱妇无乱码| 美国免费a级毛片| av视频免费观看在线观看| 在线观看免费午夜福利视频| 在线观看免费视频日本深夜| 亚洲精品中文字幕一二三四区 | 国产av精品麻豆| av一本久久久久| 欧美+亚洲+日韩+国产| av网站免费在线观看视频| 日本av手机在线免费观看| 波多野结衣av一区二区av| 高清毛片免费观看视频网站 | 菩萨蛮人人尽说江南好唐韦庄| 777久久人妻少妇嫩草av网站| 精品人妻熟女毛片av久久网站| 亚洲免费av在线视频| 亚洲五月色婷婷综合| 性高湖久久久久久久久免费观看| 最新美女视频免费是黄的| 97在线人人人人妻| 高清毛片免费观看视频网站 | 久久久久久亚洲精品国产蜜桃av| 国产成人欧美| 一边摸一边抽搐一进一小说 | 精品人妻1区二区| 一级,二级,三级黄色视频| 中文字幕高清在线视频| kizo精华| 亚洲精品国产一区二区精华液| 蜜桃国产av成人99| 亚洲国产欧美一区二区综合| 香蕉久久夜色| 亚洲熟妇熟女久久| 亚洲综合色网址| 亚洲精品在线美女| 色精品久久人妻99蜜桃| 高清欧美精品videossex| 男女床上黄色一级片免费看| 国产亚洲午夜精品一区二区久久| 波多野结衣av一区二区av| 中文欧美无线码| 不卡av一区二区三区| 亚洲熟妇熟女久久| 久久人妻av系列| 丝袜喷水一区| 一进一出抽搐动态| 热re99久久精品国产66热6| 中文字幕另类日韩欧美亚洲嫩草| 啦啦啦在线免费观看视频4| 久久精品熟女亚洲av麻豆精品| 国产欧美日韩综合在线一区二区| 亚洲三区欧美一区| 人成视频在线观看免费观看| 在线观看舔阴道视频| 宅男免费午夜| 午夜久久久在线观看| 天天躁夜夜躁狠狠躁躁| 久久性视频一级片| 国产老妇伦熟女老妇高清| 欧美日韩中文字幕国产精品一区二区三区 | 欧美精品人与动牲交sv欧美| 精品久久久精品久久久| 免费在线观看完整版高清| 桃花免费在线播放| 午夜福利视频精品| 好男人电影高清在线观看| 美女高潮喷水抽搐中文字幕| 制服诱惑二区| 成人18禁高潮啪啪吃奶动态图| 亚洲一区中文字幕在线| 国产一区有黄有色的免费视频| 日韩视频一区二区在线观看| 操美女的视频在线观看| 亚洲 欧美一区二区三区| 12—13女人毛片做爰片一| 亚洲精品美女久久av网站| 操出白浆在线播放| 热re99久久国产66热| 久久精品国产a三级三级三级| 黑丝袜美女国产一区| 精品一区二区三区四区五区乱码| 精品久久久精品久久久| 波多野结衣av一区二区av| 黑丝袜美女国产一区| 777米奇影视久久| 一边摸一边做爽爽视频免费| 国产欧美亚洲国产| 国产激情久久老熟女| 亚洲精品国产区一区二| 国产免费现黄频在线看| 欧美黄色淫秽网站| 建设人人有责人人尽责人人享有的| 国产一区二区激情短视频| 桃花免费在线播放| 麻豆成人av在线观看| 欧美激情极品国产一区二区三区| 中文字幕色久视频| 欧美成人午夜精品| 黄色视频不卡| 国产单亲对白刺激| 麻豆成人av在线观看| 欧美黄色片欧美黄色片| 日韩大码丰满熟妇| 啦啦啦在线免费观看视频4| 国产日韩欧美亚洲二区| 一夜夜www| 亚洲精品国产色婷婷电影| 日韩一区二区三区影片| 在线天堂中文资源库| 99热网站在线观看| 国产精品98久久久久久宅男小说| 欧美黑人欧美精品刺激| 中文字幕人妻熟女乱码| 99九九在线精品视频| 日韩精品免费视频一区二区三区| 又大又爽又粗| 亚洲精品av麻豆狂野| 丝袜喷水一区| 三级毛片av免费| 999精品在线视频| 18在线观看网站| 手机成人av网站| 国产有黄有色有爽视频| 一区二区三区乱码不卡18| 成年人免费黄色播放视频| 正在播放国产对白刺激| 一区福利在线观看| 亚洲精品国产一区二区精华液| 免费看a级黄色片| 18禁黄网站禁片午夜丰满| 老汉色∧v一级毛片| 视频在线观看一区二区三区| 国产精品1区2区在线观看. | 国产黄频视频在线观看| 午夜两性在线视频| 欧美日韩成人在线一区二区| 亚洲成人免费电影在线观看| 18在线观看网站| 欧美日韩精品网址| cao死你这个sao货| 精品福利观看| 一级,二级,三级黄色视频| 国产老妇伦熟女老妇高清| 欧美激情高清一区二区三区| 我要看黄色一级片免费的| 国产精品自产拍在线观看55亚洲 | 国产高清视频在线播放一区| 午夜精品久久久久久毛片777| 一区二区三区激情视频| 国产黄频视频在线观看| 999久久久精品免费观看国产| 一级,二级,三级黄色视频| 一边摸一边抽搐一进一出视频| 91av网站免费观看| 国产精品熟女久久久久浪| 丝袜美腿诱惑在线| 亚洲自偷自拍图片 自拍| 99精国产麻豆久久婷婷| 一区二区三区乱码不卡18| 一个人免费在线观看的高清视频| 99九九在线精品视频| 久久久国产精品麻豆| 久久青草综合色| 日本精品一区二区三区蜜桃| 久久久久久久久免费视频了| 亚洲久久久国产精品| 亚洲免费av在线视频| 亚洲五月色婷婷综合| a级毛片在线看网站| 精品少妇内射三级| 国产一卡二卡三卡精品| 丰满人妻熟妇乱又伦精品不卡| 久久亚洲精品不卡| 欧美精品av麻豆av| 一本一本久久a久久精品综合妖精| 国产aⅴ精品一区二区三区波| 两个人免费观看高清视频| 国产亚洲午夜精品一区二区久久| 他把我摸到了高潮在线观看 | 国产欧美亚洲国产| 美女福利国产在线| 国产av精品麻豆| 91九色精品人成在线观看| 亚洲av电影在线进入| 亚洲成人手机| 极品教师在线免费播放| 99国产精品一区二区三区| 久久人妻熟女aⅴ| 高清欧美精品videossex| 狠狠婷婷综合久久久久久88av| 天天躁日日躁夜夜躁夜夜| 免费高清在线观看日韩| 搡老熟女国产l中国老女人| 久久中文字幕人妻熟女| av不卡在线播放| 免费观看人在逋| 夜夜骑夜夜射夜夜干| 少妇裸体淫交视频免费看高清 | 少妇的丰满在线观看| 欧美av亚洲av综合av国产av| 欧美另类亚洲清纯唯美| 久久久欧美国产精品| 在线播放国产精品三级| 欧美黑人欧美精品刺激| 国产成人免费无遮挡视频| 欧美激情极品国产一区二区三区| 亚洲av成人一区二区三| 香蕉丝袜av| 免费在线观看视频国产中文字幕亚洲| 日韩欧美免费精品| 老司机在亚洲福利影院| 嫁个100分男人电影在线观看| 成人特级黄色片久久久久久久 | 国产精品99久久99久久久不卡| 亚洲黑人精品在线| 亚洲欧美日韩另类电影网站| 如日韩欧美国产精品一区二区三区| 999精品在线视频| 美女主播在线视频| 久久精品国产a三级三级三级| 亚洲成a人片在线一区二区| 一个人免费看片子| 久久久国产一区二区| 久热这里只有精品99| 正在播放国产对白刺激| 捣出白浆h1v1| 纯流量卡能插随身wifi吗| 欧美日韩福利视频一区二区| 亚洲专区中文字幕在线| 一本—道久久a久久精品蜜桃钙片| 丝袜喷水一区| 国产精品 欧美亚洲| 久久久欧美国产精品| 国产在线精品亚洲第一网站| 日韩三级视频一区二区三区| 在线观看舔阴道视频| 日本a在线网址| 国产精品久久久久久精品电影小说| 男女边摸边吃奶| 一区二区三区激情视频| 18禁黄网站禁片午夜丰满| 亚洲少妇的诱惑av| 国产男女内射视频| 人妻一区二区av| 窝窝影院91人妻| 久久国产精品大桥未久av| 老司机深夜福利视频在线观看| cao死你这个sao货| 激情视频va一区二区三区| 久9热在线精品视频| 亚洲 国产 在线| 国产成人免费观看mmmm| h视频一区二区三区| 国产人伦9x9x在线观看| 欧美激情极品国产一区二区三区| 日韩免费高清中文字幕av| 天天操日日干夜夜撸| 大片免费播放器 马上看| 一本一本久久a久久精品综合妖精| 久久国产精品人妻蜜桃| 久久精品人人爽人人爽视色| 热99国产精品久久久久久7| 色精品久久人妻99蜜桃| 悠悠久久av| 亚洲人成伊人成综合网2020| 国产淫语在线视频| 黄片播放在线免费| 国产精品一区二区在线观看99| av不卡在线播放| 1024视频免费在线观看| 老司机靠b影院| aaaaa片日本免费| 日韩视频在线欧美| 9色porny在线观看| videosex国产| 99国产精品一区二区三区| 中国美女看黄片| 12—13女人毛片做爰片一| 日本撒尿小便嘘嘘汇集6| 美女午夜性视频免费| 久久午夜综合久久蜜桃| 成人亚洲精品一区在线观看| 美女国产高潮福利片在线看| 制服诱惑二区| 在线观看人妻少妇| 夫妻午夜视频| 亚洲熟妇熟女久久| 超碰成人久久| 久久精品亚洲av国产电影网| 99国产极品粉嫩在线观看| 高清av免费在线| 亚洲精品一二三| 99国产精品99久久久久| 男女下面插进去视频免费观看| 建设人人有责人人尽责人人享有的| 女警被强在线播放| 久久狼人影院| 人人澡人人妻人| 色94色欧美一区二区| 久久精品亚洲熟妇少妇任你| 超色免费av| 91字幕亚洲| 色综合欧美亚洲国产小说| 国产伦理片在线播放av一区| 男男h啪啪无遮挡| 亚洲一码二码三码区别大吗| 国产精品久久久久久精品古装| 亚洲成a人片在线一区二区| 久热这里只有精品99| 久久99热这里只频精品6学生| 丝袜美腿诱惑在线| 欧美成人午夜精品| 精品久久久久久电影网| 国产免费av片在线观看野外av| 国精品久久久久久国模美| 中文字幕另类日韩欧美亚洲嫩草| 国产日韩欧美视频二区| 极品人妻少妇av视频| 人妻 亚洲 视频| 日本av手机在线免费观看| 国产有黄有色有爽视频| 日韩一区二区三区影片| 日本撒尿小便嘘嘘汇集6| 考比视频在线观看| 美女福利国产在线| 国产精品 国内视频| 三上悠亚av全集在线观看| 欧美成人免费av一区二区三区 | 19禁男女啪啪无遮挡网站| av不卡在线播放| 在线观看免费日韩欧美大片| 国产野战对白在线观看| 中文字幕最新亚洲高清| 日韩视频一区二区在线观看| 黑人巨大精品欧美一区二区mp4| 99久久国产精品久久久| 亚洲精品久久成人aⅴ小说| 精品国内亚洲2022精品成人 | 久久性视频一级片| 午夜福利视频在线观看免费| 国产一区二区 视频在线| 99国产综合亚洲精品| 变态另类成人亚洲欧美熟女 | 国产淫语在线视频| 亚洲色图综合在线观看| 国产成人精品久久二区二区91| 黄色 视频免费看| 亚洲国产欧美在线一区| 一二三四社区在线视频社区8| 欧美 亚洲 国产 日韩一| 两性夫妻黄色片| 国产一区有黄有色的免费视频| 国产在线精品亚洲第一网站| 国产伦人伦偷精品视频| 999久久久精品免费观看国产| 欧美黄色淫秽网站| 欧美黑人欧美精品刺激| 不卡av一区二区三区| 欧美日韩视频精品一区| 国产成人系列免费观看| 精品国产乱子伦一区二区三区| 精品免费久久久久久久清纯 | 成年人免费黄色播放视频| 在线 av 中文字幕| 操出白浆在线播放| 久久久精品免费免费高清| 免费一级毛片在线播放高清视频 | 建设人人有责人人尽责人人享有的| 桃花免费在线播放| 国产精品二区激情视频| 搡老熟女国产l中国老女人| 亚洲五月婷婷丁香| 国产在线观看jvid| 我的亚洲天堂| cao死你这个sao货| 韩国精品一区二区三区| 男女高潮啪啪啪动态图| 99国产综合亚洲精品| 精品国产超薄肉色丝袜足j| 可以免费在线观看a视频的电影网站| 天堂动漫精品| 亚洲午夜精品一区,二区,三区| 国产精品美女特级片免费视频播放器 | 俄罗斯特黄特色一大片| 精品一区二区三区视频在线观看免费 | 男女无遮挡免费网站观看| 脱女人内裤的视频| 国产精品 国内视频| 人人妻,人人澡人人爽秒播| 久久国产精品人妻蜜桃| 国产免费福利视频在线观看| 国产无遮挡羞羞视频在线观看| 婷婷丁香在线五月| 欧美日本中文国产一区发布| 美女国产高潮福利片在线看| 人人妻人人澡人人看| 国产高清激情床上av| netflix在线观看网站| 国产色视频综合| 久久狼人影院| 色综合婷婷激情| 在线观看免费日韩欧美大片| 国产一区有黄有色的免费视频| 99精品在免费线老司机午夜| 一级,二级,三级黄色视频| videos熟女内射| 久久人人爽av亚洲精品天堂| 大型av网站在线播放| av线在线观看网站| 国产aⅴ精品一区二区三区波| 男女之事视频高清在线观看| 老熟妇仑乱视频hdxx| √禁漫天堂资源中文www| 日韩视频一区二区在线观看| 国产精品久久久久久精品古装| 一级片'在线观看视频| 国产成人精品在线电影| 亚洲精品在线美女| 美女扒开内裤让男人捅视频| 90打野战视频偷拍视频| 久久久精品区二区三区| 久久久久视频综合| 啦啦啦 在线观看视频| 亚洲精品成人av观看孕妇| 日本vs欧美在线观看视频| 亚洲精品美女久久久久99蜜臀| 亚洲精品一卡2卡三卡4卡5卡|