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

    日韩国内少妇激情av| 搡老岳熟女国产| 国产欧美日韩精品一区二区| 波多野结衣高清作品| 亚洲午夜精品一区,二区,三区| 伦理电影免费视频| 97超级碰碰碰精品色视频在线观看| 久久久久久国产a免费观看| 婷婷精品国产亚洲av| 99热这里只有是精品50| 欧美另类亚洲清纯唯美| 久久久国产成人免费| 国产亚洲av高清不卡| 国产麻豆成人av免费视频| 一边摸一边抽搐一进一小说| 亚洲av第一区精品v没综合| 久久久国产欧美日韩av| 黑人操中国人逼视频| 国产视频内射| 亚洲精品中文字幕一二三四区| 97人妻精品一区二区三区麻豆| 久久久久久人人人人人| 国产精品野战在线观看| 丰满人妻熟妇乱又伦精品不卡| 操出白浆在线播放| 嫩草影院入口| 国内毛片毛片毛片毛片毛片| 亚洲国产欧美人成| 巨乳人妻的诱惑在线观看| 首页视频小说图片口味搜索| 久久中文看片网| 黄色女人牲交| 亚洲国产日韩欧美精品在线观看 | 精品乱码久久久久久99久播| 神马国产精品三级电影在线观看| 老熟妇乱子伦视频在线观看| 99久久成人亚洲精品观看| 极品教师在线免费播放| 男女下面进入的视频免费午夜| 日本与韩国留学比较| 欧美三级亚洲精品| 亚洲午夜理论影院| 欧美中文综合在线视频| 精品一区二区三区视频在线 | 免费av毛片视频| 日韩欧美国产一区二区入口| 婷婷精品国产亚洲av| 一级毛片高清免费大全| 欧洲精品卡2卡3卡4卡5卡区| av女优亚洲男人天堂 | 亚洲av成人一区二区三| 国产久久久一区二区三区| 精品一区二区三区视频在线 | 国产一区二区三区在线臀色熟女| 亚洲片人在线观看| 色哟哟哟哟哟哟| 激情在线观看视频在线高清| 久久国产精品人妻蜜桃| 日韩欧美国产在线观看| 真人一进一出gif抽搐免费| 国产伦人伦偷精品视频| 宅男免费午夜| 国产高清激情床上av| 中国美女看黄片| 国产午夜精品久久久久久| 亚洲精品美女久久久久99蜜臀| 亚洲精品美女久久av网站| 国产主播在线观看一区二区| 亚洲人成网站高清观看| 久久精品91无色码中文字幕| 亚洲性夜色夜夜综合| 欧美日韩乱码在线| 欧美成人性av电影在线观看| 午夜福利欧美成人| 成人av一区二区三区在线看| 婷婷亚洲欧美| 99国产精品99久久久久| 日本三级黄在线观看| 午夜成年电影在线免费观看| 最新在线观看一区二区三区| 中国美女看黄片| 母亲3免费完整高清在线观看| 国产成人福利小说| 男人舔奶头视频| 嫩草影视91久久| 亚洲 欧美 日韩 在线 免费| 成年免费大片在线观看| 免费在线观看视频国产中文字幕亚洲| 别揉我奶头~嗯~啊~动态视频| 欧美在线一区亚洲| 偷拍熟女少妇极品色| 欧美日韩中文字幕国产精品一区二区三区| 黄片小视频在线播放| 黄色女人牲交| 草草在线视频免费看| 男女做爰动态图高潮gif福利片| 亚洲精品美女久久久久99蜜臀| 久久中文字幕一级| 九九热线精品视视频播放| 最近视频中文字幕2019在线8| 亚洲国产欧洲综合997久久,| 国产真人三级小视频在线观看| 999久久久国产精品视频| 女同久久另类99精品国产91| 欧美绝顶高潮抽搐喷水| 色综合婷婷激情| 国产精品久久电影中文字幕| 首页视频小说图片口味搜索| 亚洲欧美精品综合久久99| 18美女黄网站色大片免费观看| 久久久国产精品麻豆| 嫩草影院精品99| 性色av乱码一区二区三区2| cao死你这个sao货| 亚洲欧美日韩卡通动漫| 成年女人永久免费观看视频| 叶爱在线成人免费视频播放| 一个人观看的视频www高清免费观看 | 99国产综合亚洲精品| 黄色 视频免费看| 国产高潮美女av| 日日摸夜夜添夜夜添小说| 亚洲第一欧美日韩一区二区三区| 中出人妻视频一区二区| 亚洲天堂国产精品一区在线| 亚洲人成电影免费在线| 好男人在线观看高清免费视频| tocl精华| 美女 人体艺术 gogo| 午夜激情福利司机影院| 国产97色在线日韩免费| 亚洲男人的天堂狠狠| 国产精品 欧美亚洲| 夜夜夜夜夜久久久久| 女人被狂操c到高潮| 熟女少妇亚洲综合色aaa.| 欧美国产日韩亚洲一区| 夜夜躁狠狠躁天天躁| netflix在线观看网站| 一个人看视频在线观看www免费 | 午夜福利在线观看吧| 亚洲国产欧美一区二区综合| 99在线视频只有这里精品首页| 国产成人一区二区三区免费视频网站| 一本久久中文字幕| 亚洲av熟女| 日本a在线网址| cao死你这个sao货| 免费看十八禁软件| 国产精华一区二区三区| 男人和女人高潮做爰伦理| 后天国语完整版免费观看| www.自偷自拍.com| 欧美日韩黄片免| 99国产极品粉嫩在线观看| 99久久久亚洲精品蜜臀av| 成年女人看的毛片在线观看| 一区二区三区激情视频| 国产视频一区二区在线看| 变态另类成人亚洲欧美熟女| 在线a可以看的网站| 亚洲 欧美 日韩 在线 免费| 99国产精品一区二区三区| 国产成+人综合+亚洲专区| 国产精品亚洲一级av第二区| 两个人的视频大全免费| 久久久久国产精品人妻aⅴ院| 国产亚洲精品一区二区www| 国产伦精品一区二区三区四那| 日韩精品中文字幕看吧| 欧美日韩黄片免| 女人高潮潮喷娇喘18禁视频| 久久人人精品亚洲av| 午夜免费观看网址| 亚洲性夜色夜夜综合| 欧美成人一区二区免费高清观看 | 亚洲天堂国产精品一区在线| 久久午夜综合久久蜜桃| 免费在线观看日本一区| 国产免费男女视频| 亚洲成人免费电影在线观看| 又粗又爽又猛毛片免费看| 日韩精品青青久久久久久| 两个人的视频大全免费| 波多野结衣高清无吗| 国内揄拍国产精品人妻在线| 巨乳人妻的诱惑在线观看| 丰满人妻一区二区三区视频av | 12—13女人毛片做爰片一| 2021天堂中文幕一二区在线观| 欧美在线黄色| 岛国视频午夜一区免费看| 成在线人永久免费视频| 日本 欧美在线| 两性夫妻黄色片| 精品久久蜜臀av无| 淫妇啪啪啪对白视频| 国产人伦9x9x在线观看| 99久久精品一区二区三区| 少妇熟女aⅴ在线视频| 黄色 视频免费看| 麻豆成人av在线观看| 此物有八面人人有两片| 美女大奶头视频| 国产精品一及| 90打野战视频偷拍视频| 成人午夜高清在线视频| 成人无遮挡网站| 好男人在线观看高清免费视频| 欧美乱妇无乱码| 99久久精品国产亚洲精品| 国产精品永久免费网站| 精品久久久久久久人妻蜜臀av| 91老司机精品| 一级毛片高清免费大全| 亚洲自偷自拍图片 自拍| 免费无遮挡裸体视频| 日本黄大片高清| 午夜福利高清视频| 久久久国产欧美日韩av| 亚洲电影在线观看av| 一区二区三区高清视频在线| 亚洲成人久久性| 久久久久久九九精品二区国产| 国产av在哪里看| 草草在线视频免费看| 亚洲av成人av| 制服丝袜大香蕉在线| 亚洲精品久久国产高清桃花| 国产真实乱freesex| 少妇的丰满在线观看| 亚洲第一电影网av| 国产精品野战在线观看| 亚洲无线观看免费| 黄色成人免费大全| 国内少妇人妻偷人精品xxx网站 | 18禁美女被吸乳视频| 国产成人影院久久av| 老司机福利观看| 91在线观看av| 欧美日韩黄片免| 亚洲国产精品sss在线观看| 热99在线观看视频| 亚洲精品粉嫩美女一区| 免费观看人在逋| 无人区码免费观看不卡| 天堂网av新在线| 全区人妻精品视频| 日韩成人在线观看一区二区三区| 美女黄网站色视频| 黄频高清免费视频| 男人和女人高潮做爰伦理| 色尼玛亚洲综合影院| 天堂av国产一区二区熟女人妻| 99热这里只有是精品50| 久久久久久久久中文| 一个人看的www免费观看视频| 天天躁日日操中文字幕| 久久久久性生活片| 国产激情偷乱视频一区二区| 亚洲欧美精品综合一区二区三区| av在线蜜桃| 国内揄拍国产精品人妻在线| 噜噜噜噜噜久久久久久91| 久久久久九九精品影院| 成年免费大片在线观看| 窝窝影院91人妻| 动漫黄色视频在线观看| 国产私拍福利视频在线观看| 99精品欧美一区二区三区四区| 久久久久久大精品| 99riav亚洲国产免费| 熟女少妇亚洲综合色aaa.| 久久久久久久久中文| 一二三四在线观看免费中文在| 91av网一区二区| 噜噜噜噜噜久久久久久91| 香蕉国产在线看| 很黄的视频免费| 99热这里只有精品一区 | 久久精品国产综合久久久| 午夜福利高清视频| 啦啦啦免费观看视频1| 九九热线精品视视频播放| 这个男人来自地球电影免费观看| 欧美黄色淫秽网站| 国产激情偷乱视频一区二区| 小说图片视频综合网站| 成人三级黄色视频| 视频区欧美日本亚洲| 桃红色精品国产亚洲av| 亚洲精品乱码久久久v下载方式 | 丰满人妻熟妇乱又伦精品不卡| 国内精品一区二区在线观看| 国产精品影院久久| 久久精品aⅴ一区二区三区四区| 90打野战视频偷拍视频| 毛片女人毛片| 欧洲精品卡2卡3卡4卡5卡区| 9191精品国产免费久久| 三级国产精品欧美在线观看 | 国产成人一区二区三区免费视频网站| 亚洲在线观看片| 国产精品野战在线观看| 夜夜看夜夜爽夜夜摸| 91麻豆精品激情在线观看国产| 国产精品久久视频播放| 欧洲精品卡2卡3卡4卡5卡区| 99热精品在线国产| 757午夜福利合集在线观看| 三级国产精品欧美在线观看 | 婷婷丁香在线五月| 亚洲国产精品合色在线| 国产精品久久久av美女十八| 欧美黄色片欧美黄色片| 亚洲国产欧美一区二区综合| 国产成人啪精品午夜网站| 白带黄色成豆腐渣| 欧美在线黄色| 国产成人福利小说| 黄频高清免费视频| 黄色片一级片一级黄色片| 国产在线精品亚洲第一网站| 欧美日韩黄片免| 长腿黑丝高跟| 一本一本综合久久| 成人一区二区视频在线观看| 中文亚洲av片在线观看爽| 美女大奶头视频| 男女之事视频高清在线观看| 很黄的视频免费| 国产精品永久免费网站| 麻豆一二三区av精品| 18禁观看日本| 久久久久久久午夜电影| 成人av在线播放网站| 人妻丰满熟妇av一区二区三区| 精品不卡国产一区二区三区| 国产不卡一卡二| 91久久精品国产一区二区成人 | 亚洲精品美女久久av网站| www.www免费av| 一卡2卡三卡四卡精品乱码亚洲| 不卡av一区二区三区| avwww免费| 欧美成人一区二区免费高清观看 | 亚洲专区国产一区二区| 他把我摸到了高潮在线观看| 亚洲七黄色美女视频| 成人亚洲精品av一区二区| 麻豆成人av在线观看| 欧美在线黄色| 成人午夜高清在线视频| 婷婷丁香在线五月| 亚洲欧美日韩东京热| 搡老岳熟女国产| 91在线精品国自产拍蜜月 | aaaaa片日本免费| 亚洲成av人片在线播放无| av天堂在线播放| 亚洲狠狠婷婷综合久久图片| 97超级碰碰碰精品色视频在线观看| 中文字幕高清在线视频| 精品久久久久久,| 久久中文字幕人妻熟女| 欧美日韩福利视频一区二区| 国产探花在线观看一区二区| 三级国产精品欧美在线观看 | 亚洲国产中文字幕在线视频| 搞女人的毛片| 91av网一区二区| 日日干狠狠操夜夜爽| 亚洲欧美精品综合久久99| 国产精品乱码一区二三区的特点| 午夜福利在线在线| 欧美一区二区国产精品久久精品| 精品一区二区三区视频在线 | 午夜成年电影在线免费观看| 欧美日韩综合久久久久久 | 手机成人av网站| 亚洲乱码一区二区免费版| 亚洲国产欧美人成| 免费看日本二区| 免费无遮挡裸体视频| 国产精品免费一区二区三区在线| 一个人观看的视频www高清免费观看 | 亚洲五月天丁香| 精品久久久久久久久久久久久| 女同久久另类99精品国产91| 久久99热这里只有精品18| 精品免费久久久久久久清纯| 国产精品亚洲av一区麻豆| 欧美不卡视频在线免费观看| 丰满的人妻完整版| 99精品久久久久人妻精品| 国产午夜精品论理片| 久久久成人免费电影| 亚洲五月婷婷丁香| 欧美大码av| 此物有八面人人有两片| 香蕉丝袜av| 老司机福利观看| 麻豆成人av在线观看| 久久久久久久久久黄片| 国产真实乱freesex| 熟妇人妻久久中文字幕3abv| 国产精品久久久久久精品电影| 久9热在线精品视频| 在线播放国产精品三级| 国产极品精品免费视频能看的| 嫩草影院精品99| 一卡2卡三卡四卡精品乱码亚洲| 白带黄色成豆腐渣| 欧美黑人欧美精品刺激| e午夜精品久久久久久久| 变态另类成人亚洲欧美熟女| 久久久久久久午夜电影| 亚洲欧美精品综合久久99| 国产精品一区二区精品视频观看| 最好的美女福利视频网| 观看美女的网站| 国产人伦9x9x在线观看| 亚洲精品一卡2卡三卡4卡5卡| 一二三四社区在线视频社区8| 很黄的视频免费| 欧美激情久久久久久爽电影| 丰满的人妻完整版| 久99久视频精品免费| 国产熟女xx| 在线永久观看黄色视频| 熟女人妻精品中文字幕| 特级一级黄色大片| 国产91精品成人一区二区三区| 久久久色成人| 国产午夜精品论理片| 精品国产美女av久久久久小说| 亚洲欧美日韩高清专用| 欧美绝顶高潮抽搐喷水| 亚洲av美国av| 又爽又黄无遮挡网站| 国产精品亚洲一级av第二区| 国产成年人精品一区二区| 热99re8久久精品国产| 成人高潮视频无遮挡免费网站| www国产在线视频色| 手机成人av网站| 成年女人永久免费观看视频| 男女那种视频在线观看| 一a级毛片在线观看| 中文字幕精品亚洲无线码一区| 特级一级黄色大片| 男女床上黄色一级片免费看| 欧美3d第一页| 夜夜看夜夜爽夜夜摸| 成在线人永久免费视频| 国产aⅴ精品一区二区三区波| 日本与韩国留学比较| 亚洲成av人片免费观看| 国产熟女xx| 老汉色av国产亚洲站长工具| 两人在一起打扑克的视频| av福利片在线观看| 亚洲中文av在线| 成人三级做爰电影| 日韩免费av在线播放| 日本 欧美在线| 人人妻人人看人人澡| 国产私拍福利视频在线观看| x7x7x7水蜜桃| 变态另类成人亚洲欧美熟女| 免费看美女性在线毛片视频| 国产探花在线观看一区二区| 日本 av在线| 国内精品久久久久久久电影| 成人特级黄色片久久久久久久| 亚洲狠狠婷婷综合久久图片| 国产久久久一区二区三区| 国产一级毛片七仙女欲春2| 久久中文字幕人妻熟女| 亚洲国产中文字幕在线视频| 每晚都被弄得嗷嗷叫到高潮| 久久亚洲真实| 亚洲一区二区三区不卡视频| 亚洲精品美女久久久久99蜜臀| 国产成人啪精品午夜网站| 老熟妇仑乱视频hdxx| 免费搜索国产男女视频| 欧美高清成人免费视频www| 在线免费观看的www视频| 特大巨黑吊av在线直播| 久久久久久人人人人人| 99国产综合亚洲精品| 久久热在线av| 国产一区二区在线观看日韩 | 国产主播在线观看一区二区| 精品无人区乱码1区二区| 亚洲国产高清在线一区二区三| 亚洲欧美一区二区三区黑人| 国产成年人精品一区二区| 很黄的视频免费| 在线观看免费视频日本深夜| 亚洲九九香蕉| 国产精品野战在线观看| 国产精品免费一区二区三区在线| 99久国产av精品| 亚洲av成人不卡在线观看播放网| 99久久无色码亚洲精品果冻| 亚洲成av人片在线播放无| 99精品欧美一区二区三区四区| 99久久综合精品五月天人人| 大型黄色视频在线免费观看| 国产精华一区二区三区| 婷婷精品国产亚洲av| 日本a在线网址| 成人永久免费在线观看视频| 国产精品一区二区三区四区久久| 99久久综合精品五月天人人| 日韩精品青青久久久久久| 毛片女人毛片| 一个人看的www免费观看视频| 免费大片18禁| 久久久久久久久中文| 免费高清视频大片| 99久久成人亚洲精品观看| 91av网一区二区| 免费在线观看亚洲国产| 1000部很黄的大片| 免费搜索国产男女视频| 91在线观看av| 久久久久久大精品| 国产精品一区二区免费欧美| 午夜福利欧美成人| e午夜精品久久久久久久| 1000部很黄的大片| 亚洲精品久久国产高清桃花| 天天一区二区日本电影三级| 久久久久久大精品| 高潮久久久久久久久久久不卡| 午夜福利欧美成人| 夜夜夜夜夜久久久久| 999久久久国产精品视频| 午夜福利视频1000在线观看| 国产乱人伦免费视频| 国产黄色小视频在线观看| 一区二区三区国产精品乱码| 五月玫瑰六月丁香| 久久久久精品国产欧美久久久| 亚洲成人中文字幕在线播放| 亚洲欧美激情综合另类| 12—13女人毛片做爰片一| 午夜免费激情av| 免费无遮挡裸体视频| 亚洲精品中文字幕一二三四区| 91在线观看av| 琪琪午夜伦伦电影理论片6080| 淫妇啪啪啪对白视频| 久久天堂一区二区三区四区| 亚洲国产中文字幕在线视频| 久久国产乱子伦精品免费另类| 非洲黑人性xxxx精品又粗又长| 悠悠久久av| 一本一本综合久久| 亚洲熟女毛片儿| 亚洲精品久久国产高清桃花| 国产黄片美女视频| 狠狠狠狠99中文字幕| 特级一级黄色大片| 日韩欧美免费精品| 人人妻人人看人人澡| 久久久久久久午夜电影| 少妇的逼水好多| 天堂√8在线中文| 国内精品久久久久精免费| АⅤ资源中文在线天堂| 99热这里只有是精品50| 亚洲激情在线av| 国内久久婷婷六月综合欲色啪| 五月伊人婷婷丁香| 亚洲欧美精品综合久久99| 亚洲国产精品sss在线观看| 久久草成人影院| 一级a爱片免费观看的视频| 久久草成人影院| 18禁黄网站禁片免费观看直播| 最新在线观看一区二区三区| 亚洲精品在线观看二区| 成人性生交大片免费视频hd| 国产精品综合久久久久久久免费| 精品熟女少妇八av免费久了| 久久久久久人人人人人| 国产美女午夜福利| 日韩欧美国产一区二区入口| 午夜激情福利司机影院| 婷婷精品国产亚洲av在线| 非洲黑人性xxxx精品又粗又长| 成人国产一区最新在线观看| 99久久国产精品久久久| 91字幕亚洲| 女警被强在线播放| 18禁观看日本| aaaaa片日本免费| 亚洲欧美日韩无卡精品| 成人三级做爰电影| 久久精品国产亚洲av香蕉五月| 日韩欧美国产在线观看| 丰满人妻一区二区三区视频av | 在线观看美女被高潮喷水网站 | 国产又色又爽无遮挡免费看| 欧美中文综合在线视频| 国产欧美日韩一区二区三| 国产成人av激情在线播放| 国产成人福利小说| 国产成人精品久久二区二区免费| 中文字幕熟女人妻在线| 亚洲专区字幕在线|