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

    Monthly Variations of Atmospheric Circulations Associated with Haze Pollution in the Yangtze River Delta and North China

    2021-04-07 10:20:54XinyuZHANGZhicongYINHuijunWANGandMingkengDUAN
    Advances in Atmospheric Sciences 2021年4期

    Xinyu ZHANG,Zhicong YIN*,2,3,Huijun WANG,2,3,and Mingkeng DUAN

    1Key Laboratory of Meteorological Disaster, Ministry of Education/Joint International Research Laboratory of Climate and Environment Change/Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters,Nanjing University of Information Science and Technology, Nanjing 210044, China

    2Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Zhuhai 519080, China

    3Nansen-Zhu International Research Centre, Institute of Atmospheric Physics,Chinese Academy of Sciences, Beijing 100029, China

    ABSTRACT Haze pollution in early winter (December and January) in the Yangtze River Delta (YRD) and in North China (NC)are both severe;however,their monthly variations are significantly different.In this study,the dominant large-scale atmospheric circulations and local meteorological conditions were investigated and compared over the YRD and NC in each month.Results showed that the YRD (NC) is dominated by the so-called Scandinavia (East Atlantic/West Russia)pattern in December,and these circulations weaken in January.The East Asian December and January monsoons over the YRD and NC have negative correlations with the number of haze days.The local descending motion facilitates less removal of haze pollution over the YRD,while the local ascending motion facilitates less removal of haze pollution over NC in January,despite a weaker relationship in December.Additionally,the monthly variations of atmospheric circulations showed that adverse meteorological conditions restrict the vertical (horizontal) dispersion of haze pollution in December(January) over the YRD,while the associated local weather conditions are similar in these two months over NC.

    Key words:monthly haze pollution,atmospheric circulations,Yangtze River Delta,North China

    1.Introduction

    Heavy haze pollution often occurs in winter (December—February) in North China (NC) and the Yangtze River Delta (YRD) (Ding and Liu,2014).Haze pollution is mainly characterized by high concentrations of fine particulate matter and reduced visibility,causing serious socioeconomic problems (Wang and Chen,2016).Thus,haze pollution has aroused widespread concern in government and among the public.Early studies showed that haze pollution is not only affected by anthropogenic emissions (mainly the long-term trend of haze),but also by climate variability(mainly the interannual variations) (Wang et al.,2013;Yang et al.,2016;Zhang et al.,2016).

    Several researchers have investigated the pollution situation over the YRD from monthly,seasonal and annual mean perspectives (Yu et al.,2011;Zhao et al.,2017;Sun et al.,2019).Affected by urbanization,large cities in the YRD have higher aerosol optical depths than remote areas (Sun et al.,2019),as well as high anthropogenic heat release (Wang et al.,2015),as verified by high-resolution Weather Research and Forecasting model results.In addition to impacts of anthropogenic emissions on the long-term trend of haze pollution in the YRD,a few studies have investigated the roles of weather—climate conditions in the YRD.He et al.(2019b) pointed out that PMthroughout the Yangtze River Basin is negatively correlated with surface wind speed,precipitation and planetary boundary layer height (PBLH) but positively correlated with temperature.High wind speed,more precipitation and a well-developed boundary layer could facilitate the dispersion of aerosols and lead to a decrease in PM(He et al.,2019b).The positive relationship between temperature and PMis possibly due to the enhancement in oxidation of aerosol precursors at high temperatures (He et al.,2019b).Compared to the conditions in NC,He et al.(2019a) also indicated that wet removal of haze brought by precipitation over southern China (south of 30°N) is important.Furthermore,winter haze extremes in India at lower latitudes are even affected by El Ni?o sea surface temperature and signals from the Southern Hemisphere (Gao et al.,2019).

    NC is the most haze-polluted area in China,and the associated atmospheric circulations have already been examined in previous studies.The number of winter haze days (HDs)in NC (WHD) has been found to be directly dominated by local anticyclonic circulations (Chen and Wang,2015;Yin and Wang,2016),connected by two Rossby wave trains [i.e.,the East Atlantic/West Russia (EA/WR) and Eurasia patterns]over the continent and one Rossby wave train (i.e.,the West Pacific pattern) over the ocean (Wallace and Gutzler,1981;Barnston and Livezey,1987).These wave trains provide the favorable atmospheric conditions associated with haze pollution,such as a lower planetary boundary layer,lower wind speed,anomalous southerlies,and higher relative humidity (RH),which restrict the horizontal and vertical ventilation of haze particles (Yin et al.,2017).

    The number of HDs in December over the YRD(DHD),as well as those in January (JHD),slowly increased from 1981 to 2017,while those over NC (DHDand JHD) slightly decreased from 1981 to 2010 but increased from 2011 to 2017 (Fig.1).In this study,we pay more attention to the interannual variations of monthly haze pollution.The correlation coefficient (CC) between DHDand JHDis 0.71;however it decreases to 0.33 in the YRD,indicating significant differences in haze variations in December and January over the YRD.The time series of DHDand DHD(JHDand JHD) during 1981—2017 has a CC of 0.49 (0.2),which are not high for two adjacent economically developed regions (Figs.1a and b).Furthermore,we also compared the variations in emissions (i.e.,the difference of emissions in different months and regions) with the associated differences in HDs and found the CCs between them were all lower than 0.2 [non-significant;Fig.S1 in electronic supplementary material(ESM) ].The above statistical results illustrate the need for an analysis of the monthly variations of atmospheric circulations associated with haze pollution (i.e.,the impacts of meteorology) in the YRD and NC.

    2.Datasets and methods

    The monthly geopotential height,horizontal wind(UV),vertical wind (omega),temperature,and RH at different levels for the winters of 1981—2018,as well as surface RH,surface UV,surface air temperature,and sea surface level,are from the NCEP—NCAR reanalyses with a horizontal resolution of 2.5° × 2.5° (Kalnay et al.,1996).The precipitation data are from Climate Prediction Center Merged Analysis of Precipitation (CMAP),with a horizontal resolution of 2.5° × 2.5°,in winters from 1981 to 2018 (Xie and Arkin,1997).Other precipitation data,including the 1° × 1°Precipitation Reconstruction over Land of the National Oceanic and Atmospheric Administration,as well as a station dataset of the China Meteorological Administration(Chen et al.,2002),were used to verify the results.Monthly PBLH data are also used here (resolution:1° × 1°),but were derived from the ERA-Interim dataset with the same time period as the above datasets (Dee et al.,2011).In addition,the EA/WR and Scandinavia patterns are two of three prominent teleconnection patterns that affect Eurasia,and have been referred to as the Eurasia-2 and Eurasia-1 pattern,respectively (Barnston and Livezey,1987).The EA/WR pattern is composed of two positive height anomaly centers over Europe and northern China,and two negative height anomaly centers over the central North Atlantic and north of the Caspian Sea,when in positive phase.The Scandinavia pattern has four centers in Eurasia (i.e.,Spain and the adjacent Mediterranean and Atlantic area,Scandinavia,Northwest China or western Mongolia,and Japan),and sometimes some centers do not show up.The calculation procedures for these two indices are available at https://www.cpc.ncep.noaa.gov/data/teledoc/teleindcalc.shtml,and numerical values were downloaded from https://www.cpc.ncep.noaa.gov/data/teledoc/telecontents.shtml.The emissions inventory data are from http://inventory.pku.edu.cn/.The calculations for atmospheric circulations,including the anomalous anticyclone (AAC),surface air temperature (SAT),surface temperature difference between the north and south (SAT diff.),RH,PBLH,surface wind speed (UV speed),and precipitation,are as follows.The areas calculated as AAC indices are shown in Figs.2a and b for December and Figs.2c and d for January.The areas calculated as SAT diff.indices are[R1 (38°—70°N,85°—130°E) ? R2 (8°—23°N,95°—120°E)]and [R1 (36°—53°N,83°—115°E) ? R2 (15°—19°N,110°—115°E)]in December and January,respectively.Other variables’ indices are defined as area-averaged values of the studied areas (the YRD and NC).

    Fig.1.(a) Variation of December haze days (HDs) in the Yangtze River Delta (YRD) and in North China (NC) from 1981 to 2017.(b) Variation of January HDs in the YRD and in NC from 1982 to 2018.Solid (dotted) lines represent detrended (original) data.The correlation coefficients (CCs) of two solid lines are shown in the spaces.

    As for haze pollution,the number of HDs in each month in the winters of 1981—2017 were adopted owing to the lack of long-term PMdata.Haze was calculated by defining visibility less than a certain threshold,which was 10 km before 1 January 2014 for most stations and 7.5 km after because of the switch in visibility observation means(i.e.,from manual observation to automatic observation)according to a report of the China Meteorological Administration in 2014 (Yin et al.,2017) and RH less than 90% (Chen and Wang,2015;Yin and Wang,2016).One HD at a station was defined as a day with haze at any time,after excluding other weather phenomena that affected visibility.All the HDs of a month at all stations in the region were accumulated,and then the average value of all stations in the region was calculated as the HDs in the region.HDs over the YRD(27.05°—33.77°N,118.15°—121.95°E) and NC (34°—42°N,114°—120°E) were computed as the HD average value in these two areas,respectively.For more detail,readers are referred to Chen and Wang (2015),Yin and Wang (2016),and Yin et al.(2017).

    Fig.2.(a) Regressions of December geopotential height (contours;units:m) and the associated wave activity flux(WAF;arrows;units:m2 s?2) at 500 hPa on the December HDs in the YRD during 1981—2017.(b) As in (a) but in NC.(c) Regressions of the January geopotential height (contours;units:m) and the associated WAF (arrows;units:m2 s?2) at 500 hPa on the January HDs in the YRD during 1982—2018.(d) As in (c) but in NC.Solid lines and dotted lines represent positive and negative values of regression,respectively.Shaded regions indicate that the 500-hPa geopotential height anomalies are significant at the 95% confidence level,and yellow (green) color represents positive (negative) regression coefficients passing the Student’s t-test.The linear trend of the data was removed.The purple (gray) solid-line box indicates the location of the YRD (NC).The orange dotted-line boxes in (a,b)respectively indicate the regions of the AAC indices of the YRD (38°—60°N,90°—123°E) and NC (30°—49°N,100°—135°E) in December.The orange dotted-line boxes in (c,d) respectively indicate the regions of the AAC indices of the YRD (30°—49°N,85°—115°E) and NC (35°—40°N,115°—135°E) in January.

    3.Results

    3.1.Atmospheric circulations drive the monthly variations in haze

    The atmospheric circulations associated with December and January haze pollution in the YRD were investigated.As illustrated in Fig.2a,a Rossby wave—like train at middle and high latitudes of continental Eurasia strongly contributes to the variations in DHD.In comparison,in January,the Rossby wave train becomes non-significant,and the only significant centers are the local AAC located in the northwest of the YRD (Fig.2c).That is,the large-scale atmospheric circulation affecting HDs degenerate to the local anticyclone from December to January.Lower PBLH (Fig.3c)in December is apparent because of anomalous descending motions (Fig.3a) over the YRD.Weak surface wind speed,induced by anomalous southerly wind,is apparent in January over the YRD (Fig.3d).Thus,the vertical movement(Figs.3a and c) is restricted in December and the horizontal movement (Figs.3b and d) is restricted in January,over the YRD,leading to high contemporaneous HDs.Furthermore,the anomalous precipitation plays a critical role in the removal of haze particles over the YRD (Figs.4a and b).Precipitation and HDs over the YRD show a significant outof-phase relationship in both December (CC=?0.34) and January (CC=?0.46) (Figs.4a and b).

    Fig.3.(a) Cross-section (115°—125°E mean) CCs between December HDs in the YRD and December air temperature (shading),RH (contours),and wind (arrows).(b) As in (a) but in January.(c) The CCs between December HDs in the YRD and December PBLH (shading) and surface wind (arrows) from 1981 to 2017.(d) The CCs between January HDs in the YRD and January surface wind speed (shading) and surface wind (arrows) from 1982 to 2018.All the above CCs share the same colorbar at the bottom of the figure.The south wind component in the meridional direction (a—d),the east wind component in the latitudinal direction (c,d),and the ascending movement in the vertical direction (a,b;owing to omega multiplied by ?1) are positive,and vice versa.White slanted lines indicate that the CCs represented by the shading exceed the 95% confidence level.The linear trends of all data were removed.The purple box or purple line indicates the location of the YRD.

    Fig.4.(a) Variation of December HDs in the YRD and the December area-averaged precipitation from CMAP datasets,with their CC in the top-right corner.(b) Variation of January HDs in the YRD and January area-averaged precipitation from CMAP datasets,with their CC in the top-right corner.(c,d) As in (a,b) but in NC.The linear trends of all data were removed.

    Atmospheric circulations associated with haze pollution in NC in December and January were also examined.In the mid-troposphere,a Rossby wave—like train clearly affects the haze pollution over NC in December (Fig.2b).However,in January,the associated atmospheric circulations become relatively local (Fig.2d),which is different from the results revealed in many previous studies for the whole winter season (Chen and Wang,2015;Yin and Wang,2016).Anomalous ascending motion (Figs.5a and b),warming near the surface (Figs.5a and b),southerly wind (Figs.5c and d),weak surface wind speed (Figs.S2a in the ESM and 5d),low PBLH (Figs.5c and S3a),high potential of thermal inversion,i.e.,the air temperature at 925 hPa minus the air temperature at 1000 hPa (Figs.S2b and S3b),and high RH (Figs.5a and b,Figs.S2c and S3c),are favorable for both DHDand JHD.That is,the monthly variations of associated local weather conditions are similar.However,precipitation makes different contributions to DHDand JHD.December precipitation is almost unrelated to DHD(Fig.4c),while January precipitation has a significantly positive relationship with DHD(CC=0.33,exceeding the 95% confidence level;Fig.4d).

    From the data in Table 1,we can see that the AAC index is crucial for both DHDand DHD,while the SAT diff.index and precipitation index (SAT index,RH index,PBLH index and UV speed index) play important roles in the YRD (NC).To quantify the contributions of atmospheric circulations to haze pollution,the highly correlated variables (≥ 95% confidence level) were selected to regress on the number of HDs.In other words,the AAC,SAT diff.and precipitation indices (ACC,SAT,RH,PBLH and UV speed indices) in December were selected as regression factors of DHD(DHD).The predicted DHDand DHDdid not include the observed trend here,and in the last calculation step we added the trend as fitted values,shown in Fig.6.The fitted and observed DHDand DHDhad a CC of 0.77 and 0.67 for the YRD and NC,respectively (Figs.6a and b).That is,atmospheric circulations made variance contributions of 59% and 45%,respectively (Figs.6a and b),confirming that these circulation variables have large influences on DHDand DHD.DHDhad the same sign as the AAC (SAT difference between the north and the south of the YRD) in 33 (31) of 37 years,while it had opposite signs to precipitation in 25 of 37 years (Fig.7a).DHDshared the same mathematical sign as the AAC,local SAT and RH in 24,28 and 28 of 37 years,respectively,while it had opposite mathematical signs to the PBLH and surface wind in 26 and 25 of 37 years,respectively (Fig.7b).As can be seen from Fig.7a,three variables correspond to DHDin the minimum year of 2010,while it is two variables in the maximum year of 2013.In contrast,as shown in Fig.7b,all five variables correspond to DHDin the minimum (maximum) year of 2010 (2015),which highlights the importance of meteorological and climate conditions.

    Fig.5.(a) The cross-section (110°—120°E mean) CCs between December HDs in NC and December air temperature(shading),RH (contours),and wind (arrows).(b) As in (a) but in January.(c) The CCs between December HDs in NC and December PBLH (shading) and surface wind (arrows) from 1981 to 2017.(d) The CCs between January HDs in NC and January surface wind speed (shading) and surface wind (arrows) from 1982 to 2018.All the above CCs share the same colorbar at the bottom of the figure.The south wind component in the meridional direction(a—d),the east wind component in the latitudinal direction (c,d),and the ascending movement in the vertical direction (a,b;owing to omega multiplied by ?1) are positive,and vice versa.White slanted lines indicate that the CCs represented by the shading exceed the 95% confidence level.The linear trends of all the data were removed.The white box or gray line indicates the location of NC.

    Table 1.The correlation coefficients (CCs) between atmospheric circulation indices and haze days (HDs) of two regions (the YRD and NC) in early winter (December and January).The linear trends in all data were removed.Double asterisks (**) and triple asterisks (***)indicate the results pass the 95% and 99% confidence level,respectively.A forward-slash (/) indicates the significance level is below 95%.

    Fig.6.Variation in the standardized observed HDs (black bars) and fitted HDs (red solid line) in the YRD in (a)December from 1981 to 2017 and (b) January from 1982 to 2018.Fitted detrended HDs plus the observed HDs’trend are fitted HDs,represented here by the red solid line.Red dotted lines represent their same trend line.The number on the right is the CC of the observed HDs and fitted HDs plus the trend (i.e.,the CC between the black bars and the red solid line).(c,d) As in (a,b) but in NC.

    Similarly,CCs were calculated between circulation indices in January and JHD(JHD),as shown in Table 1,and the significant indices (exceeding the 95% confidence level) were selected to fit JHD(JHD).Three (five) circulation indices,including the AAC,SAT diff.and precipitation indices (the ACC,RH,PBLH,UV speed and precipitation indices),were demonstrated to be useful for the formation of JHD(JHD).The fitted circulation indices explain 56% (42%) of the variance of JHD(JHD)with the trend,indicating large contributions of these atmospheric circulations (Figs.6c and d).In Fig.7c,the AAC,SAT and precipitation all show corresponding results with JHDin the minimum (maximum) year of 1993 (2015).Meanwhile,in Fig.7d,five meteorological variables show corresponding results with JHDin the minimum year of 2011;however,only RH shows corresponding results with JHDin the maximum year of 2016.Regardless,meteorological and climate conditions play important roles in the interannual variabilities of HDs rather than other elements.

    3.2.Comparison of atmospheric circulations associated with DHDNCand DHDYRD

    The wave activity flux (WAF) demonstrates the wave energy propagation paths (Takaya and Nakamura,2001;Hsu and Lin,2007;Honda et al.,2009).A quasi-stationary planetary wave route (i.e.,the Scandinavia pattern) is apparent,with an anomalous negative center over Scandinavia and an anomalous positive center over Mongolia (Fig.2a).The CC between the Scandinavia pattern and DHDis?0.48 (significant at the 99% confidence level).In contrast,the WAF associated with DHD(Fig.2b) shows a positive EA/WR phase (CC=0.56).These results demonstrate that the Scandinavia and EA/WR patterns contribute to the variations of DHDand DHD,respectively.The AAC near Northeast Asia is located in the crossover regions of the Scandinavia and EA/WR patterns and might explain the close relationship between DHDand DHD.

    Local meteorological conditions associated with DHDand DHDwere also compared (Figs.3—5).The climate-mean vertical motion is weak in the downward direction over the YRD (not shown),and the AAC enhances the anomalous descending movement (Fig.3a).Besides,anomalous reversion of the north—south temperature is favorable for DHD(Fig.3a),which restricts the invasion of the East Asian winter monsoon (EAWM) into the YRD.Low PBLH also restricts the vertical dispersion (Fig.3c).Against the background of intensive descending motion over the YRD,dry air is carried from the top to the ground,causing the weak negative RH anomaly (Fig.3a) and less precipitation.As we know,some pollutants can be removed via wet removal,i.e.,precipitation,especially in southern areas with large precipitation.In Fig.4a,precipitation has a negative CC with DHD(CC=?0.34,≥ 95% confidence level).Local meteorological conditions related to DHDare correspondingly discussed (Figs.5a and c).Different from DHD,anomalous upward movement (Fig.5a),southerly wind near the surface (Fig.5a),high RH (Figs.5a and S2c),high surface temperature (Fig.5a),high potential of temperature inversion (Fig.S2b),and low surface wind speed(Fig.S2a),are favorable for DHD.In addition,precipitation has little influence on the elimination of haze pollution over NC (Fig.4c),as compared to that over the YRD (Fig.4a).As shown in Figs.3c and 5c,the PBLH shows negative anomalies over both NC and the YRD,indicating a negative contribution to the elimination of haze pollution.

    Fig.7.(a) Variation of December HDs (line) and associated significant meteorological variables (histogram) in the YRD from 1981 to 2017.(b) As in (a) but in NC.(c) Variation of January HDs (line) and associated significant meteorological variables (histogram) in the YRD from 1982 to 2018.(d) As in (c) but in NC.Dark green,pink,cyan,brown,blue,gray,and purple solid circles at the top indicate the same or opposite signs between HDs and meteorological variables (AAC,SAT,SAT diff.,RH,PBLH,UV speed,precipitation) according to the CCs in Table 1.The linear trends of all data were removed.

    It is common knowledge that the wintertime weather and climate of East Asia is affected by the EAWM,so the monthly EAWM [i.e.,the East Asian December (January)monsoon,abbreviated as EADM (EAJM)]was separately analyzed to study the effects on DHDand DHD.Here,the definitions of the EADM and EAJM refer to previous studies (Wang and Jiang,2004;He and Wang,2012) and are shown in Table S1 (in the ESM),including the sub-indices and synthetic index.The CCs of the EADM and DHD(DHD) show that specific members,such as the East Asian deep trough,westerly jet,as well as the synthetic index,have an influence on DHD,while all members of the EADM have a significant impact on DHD(Table 2).The question,then,is how does the EADM affect DHDand DHD? A strong EADM may weaken the temperature difference between the north and south of the YRD (Fig.8a)and bring cold air from the mid—high latitudes to China.Meanwhile,there is sufficient water vapor in the YRD,so the intrusion of cold air may enhance the precipitation over the YRD (Fig.8b).When precipitation increases,haze pollution in the YRD will be removed.In addition,a strong EADM may also increase the surface wind speed in theYRD (not shown),strengthening the conditions of horizontal dispersion.As for DHD,a strong EADM may lead to favorable local ventilation conditions and unfavorable growth conditions for haze particles,such as low SAT (Fig.8c),low RH (Fig.8d),high PBLH (Fig.8e),and high surface wind speed (Fig.8f).

    Table 2.The CCs between East Asian early winter (December and January) monsoon indices and HDs of two areas (the YRD and NC) in December and January,respectively.Double asterisks(**) and triple asterisks (***) indicate the results pass the 95% and 99% confidence level,respectively.

    Fig.8.Scatterplot of (a) December East Asian jet index and December SAT diff.index between the north and south of the YRD,and (b) the December SAT diff.index between the north and south of the YRD and the December precipitation index of the YRD.Scatterplot of the East Asian December monsoon synthetic index and December (c)SAT index,(d) RH index,(e) PBLH index,and (f) UV speed index of NC.The linear trends in all data were removed.The straight lines are the lines of best fit,and the CCs indicate the CCs between two variables.

    3.3.Comparison of atmospheric circulations associated with JHDNCand JHDYRD

    The CCs of JHDand JHDbecome non-significant (Fig.1b),indicating weak synergistic variations of JHDand JHD.The associated wave-like train is significantly weakened (Figs.2c and d) in January rather than that in December.Consequently,the influence of Eurasian wave trains on local HDs is weakened.Local anticyclones affect local HDs in both areas,with CCs of 0.40 and 0.37(Table 1),significant at the 95% confidence level.Additionally,the AAC in January is also relatively weaker than that in December (Fig.2,Table 1).

    It was shown in Zhong et al.(2019) that the directions of vertical motion vary according to the locations of the AAC.Because the YRD is located in the east of the bottom of the AAC (Fig.2c),descending motion was delivered to the region (Fig.3b).Although the descending motion only exists in the south of the YRD,other local meteorological conditions induced by the AAC,such as the SAT diff.(Fig.3b),weak UV speed (Fig.3d),anomalous southerly wind (Fig.3d),and less precipitation (Fig.4b),are crucial for high JHD.By contrast,for JHD,anomalous ascending motion(Fig.5b),high SAT (Fig.5b),high RH (Figs.5b and S3c in the ESM),anomalous southerly wind (Figs.5b and d),weak UV speed (Fig.5d),low PBLH (Fig.S3a),and high potential of temperature inversion (Fig.S3b),provide adverse conditions for the dispersion of haze particles.Besides,in January,precipitation makes different contributions to HDs over NC and the YRD,i.e.,more precipitation brings more JHD(Fig.4d).Light rain is dominant when precipitation happens over NC,which is conducive to the hygroscopic growth of haze particles.

    As for the EAJM,the significant impact factors are the Siberian high index and synthetic index (Siberian high index,Aleutian low index and meridional wind at 850 hPa index) for JHD(JHD).Additionally,the influential process of the EAJM is as follows.A strong EAJM can reduce the temperature difference between the north and south of the YRD and enhance the precipitation over the YRD(Figs.9a and b).Meanwhile,for JHD,the Siberian high can affect local meteorological conditions (Figs.9c—e),including RH (Fig.9c),PBLH (Fig.9d),and UV speed (Fig.9e).In other words,a strong EAJM,represented by the Siberian high,creates an unstable environment for high JHD.

    Fig.9.Scatterplot of (a) the East Asian January monsoon synthetic index and January SAT diff.index between the north and south of the YRD,and (b) the January SAT diff.index between the north and south of the YRD and the January precipitation index of the YRD.Scatterplot of the January Siberian high index and the January (c) RH index,(d) PBLH index,and (e) UV speed index of NC.The linear trends of all data were removed.The straight lines are the lines of best fit,and the CCs indicate the CCs between two variables.

    4.Conclusions and discussion

    Monthly haze pollution over the YRD and NC was investigated in early winter (i.e.,December and January),as well as the associated large-scale atmospheric circulations and local meteorological conditions (Fig.10).The DHD(DHD) is dominated by the Scandinavia (EA/WR) pattern in December,and these two wave trains degenerate into local AACs in January.In December,the DHD(DHD) is dominated by descending motion,the large temperature difference between the north and south,and less precipitation (ascending motion,high SAT,high RH,low PBLH,and weak UV speed).The precipitation is more important in January,and it makes a negative (positive) contribution to JHD(JHD).Most members of the EAWM make negative contributions to DHD and JHD over the two areas by affecting the local meteorological conditions.In fact,the locations of the study areas relative to the AAC might affect local circulations,according to Zhong et al.(2019).Because the YRD is located in the front of the AAC and NC at the back,it stimulates the descending (ascending)motion in the vertical direction for the YRD (NC) in both December and January,which results in the differences of local meteorological conditions,especially for the precipitation.

    Fig.10.Schematic diagram of the subseasonal variation of the atmospheric circulation associated with haze pollution in the YRD and NC.The schematic of the YRD and NC in December (January) is on the left (right).Ellipses represent anomalous cyclones and anticyclones at high levels.Black (purple) curves on the left indicate the EA/WR(Scandinavia) pattern.Vertical arrows indicate the directions of vertical movement.Warm (cold) colors at the surface indicate warm (cold) anomalous surface air temperature.Green (red) clouds indicate positive (negative)anomalies of precipitation.Yellow crescents indicate a low boundary layer height.The purple (gray) box indicates the location of the YRD (NC).

    Utilizing these significant circulation indices to establish regression models of HDs,results showed that these circulation indices make large contributions and can explained more than 40% of the variance of DHD,JHD,DHDand JHD.Owing to the switch in the observation methods of visibility in January 2014,the reliability of HDs after January 2014 has been tested.It was found that the conclusions are still robust if the data after January 2014 are removed.The collinearity of the indices can be ignored because the variance inflation factors of most selected factors are lower than five in the four regression equations(Fig.6).In Fig.11 in Zhang et al.(2019),the adverse weather conditions in autumn and winter might reduce the effect of emissions reduction in the Beijing—Tianjin—Hebei region in 2014 and 2015.Yin and Zhang (2020) investigated the pollution situation in the winters of 2017 and 2018 under almost the same emission conditions,finding that haze pollution rebounded in 2018 compared with that in 2017,indicating the important effect of climate and weather conditions.

    Acknowledgements

    .This research was supported by the National Key Research and Development Plan (Grant No.2016YFA0600703),the National Natural Science Foundation of China (Grant Nos.91744311,41991283 and 41705058),and the funding of the Jiangsu Innovation &Entrepreneurship Team.

    Electronic supplementary material:

    Supplementary material is available in the online version of this article at https://doi.org/10.1007/s00376-020-0227-z.

    男人舔女人下体高潮全视频| 免费无遮挡裸体视频| 韩国av一区二区三区四区| 国产亚洲精品av在线| 亚洲国产日韩欧美精品在线观看 | 日本黄色视频三级网站网址| 黄网站色视频无遮挡免费观看| 久久热在线av| 欧美久久黑人一区二区| 精品久久久久久久久久免费视频| 丁香六月欧美| 一区二区三区国产精品乱码| 美国免费a级毛片| 亚洲av中文字字幕乱码综合 | 88av欧美| 久久精品影院6| 亚洲第一青青草原| 巨乳人妻的诱惑在线观看| 国产精品久久久人人做人人爽| 亚洲,欧美精品.| 亚洲熟妇熟女久久| 99精品欧美一区二区三区四区| 美女高潮到喷水免费观看| 狠狠狠狠99中文字幕| 欧美成狂野欧美在线观看| 俺也久久电影网| 中文字幕精品亚洲无线码一区 | 精品国产国语对白av| 国内精品久久久久精免费| 精品卡一卡二卡四卡免费| 婷婷亚洲欧美| 999精品在线视频| 欧美日韩精品网址| 妹子高潮喷水视频| 国产爱豆传媒在线观看 | 欧美一区二区精品小视频在线| 婷婷精品国产亚洲av在线| 国产精品香港三级国产av潘金莲| 国产男靠女视频免费网站| 国产高清激情床上av| 哪里可以看免费的av片| 久久久国产精品麻豆| 日韩欧美三级三区| 欧美激情极品国产一区二区三区| 久久中文字幕人妻熟女| 大型av网站在线播放| 日韩av在线大香蕉| 大香蕉久久成人网| 非洲黑人性xxxx精品又粗又长| 91国产中文字幕| 十分钟在线观看高清视频www| 美女扒开内裤让男人捅视频| 久久久久久久久免费视频了| 欧美中文综合在线视频| 怎么达到女性高潮| 欧美 亚洲 国产 日韩一| 亚洲熟女毛片儿| 成人av一区二区三区在线看| av中文乱码字幕在线| 国产又色又爽无遮挡免费看| 亚洲第一青青草原| a级毛片a级免费在线| 大香蕉久久成人网| 国产日本99.免费观看| 97超级碰碰碰精品色视频在线观看| 精品无人区乱码1区二区| 亚洲男人的天堂狠狠| 99精品在免费线老司机午夜| 757午夜福利合集在线观看| 一本大道久久a久久精品| 久久久久久免费高清国产稀缺| 村上凉子中文字幕在线| 亚洲激情在线av| 久久久国产精品麻豆| 亚洲精华国产精华精| 三级毛片av免费| 正在播放国产对白刺激| 人人妻人人看人人澡| 久久国产精品影院| 国产亚洲精品一区二区www| 白带黄色成豆腐渣| 成人国产一区最新在线观看| 国产又色又爽无遮挡免费看| 国产一区二区激情短视频| 亚洲国产精品成人综合色| 麻豆成人午夜福利视频| 欧美一级毛片孕妇| 中文字幕久久专区| 亚洲专区字幕在线| 国产99白浆流出| 国产成人系列免费观看| 亚洲精品美女久久久久99蜜臀| 制服丝袜大香蕉在线| 桃红色精品国产亚洲av| 久久久国产成人精品二区| 黑人欧美特级aaaaaa片| 久久精品国产综合久久久| www日本在线高清视频| 12—13女人毛片做爰片一| 99re在线观看精品视频| 国产三级黄色录像| 国产精品久久视频播放| 无人区码免费观看不卡| 级片在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲色图 男人天堂 中文字幕| 国产真实乱freesex| 精品日产1卡2卡| 午夜日韩欧美国产| 国产一级毛片七仙女欲春2 | 亚洲色图 男人天堂 中文字幕| 亚洲avbb在线观看| 麻豆一二三区av精品| 国产精品免费视频内射| 一二三四在线观看免费中文在| 女人高潮潮喷娇喘18禁视频| 国产av不卡久久| 国产精品久久视频播放| 午夜精品久久久久久毛片777| 熟女少妇亚洲综合色aaa.| 精品免费久久久久久久清纯| 极品教师在线免费播放| 一边摸一边抽搐一进一小说| 国产真实乱freesex| 久久久精品国产亚洲av高清涩受| 高清在线国产一区| 欧美乱码精品一区二区三区| 女性被躁到高潮视频| 久久久久久人人人人人| 一进一出好大好爽视频| 免费观看人在逋| 免费观看人在逋| 欧美色视频一区免费| bbb黄色大片| 99在线视频只有这里精品首页| 国产精品美女特级片免费视频播放器 | 两人在一起打扑克的视频| 观看免费一级毛片| 国产免费男女视频| 女同久久另类99精品国产91| 久久久国产欧美日韩av| 日韩 欧美 亚洲 中文字幕| 日韩欧美免费精品| 久久久久久亚洲精品国产蜜桃av| 国产精品亚洲美女久久久| 午夜免费成人在线视频| 成人亚洲精品一区在线观看| 18禁黄网站禁片午夜丰满| 精品久久久久久久久久免费视频| 久久久久久免费高清国产稀缺| 久久精品影院6| 欧美日韩福利视频一区二区| 久久午夜综合久久蜜桃| 日本成人三级电影网站| 久久婷婷人人爽人人干人人爱| 叶爱在线成人免费视频播放| 18美女黄网站色大片免费观看| 性欧美人与动物交配| 久久国产精品人妻蜜桃| 亚洲午夜精品一区,二区,三区| 久久久国产欧美日韩av| 精品乱码久久久久久99久播| 国产高清有码在线观看视频 | 精品乱码久久久久久99久播| 亚洲男人天堂网一区| 欧美zozozo另类| 亚洲人成伊人成综合网2020| 99精品在免费线老司机午夜| 激情在线观看视频在线高清| 女警被强在线播放| 老司机深夜福利视频在线观看| 欧美一级毛片孕妇| 9191精品国产免费久久| 18禁观看日本| 亚洲中文av在线| 可以在线观看的亚洲视频| 成年人黄色毛片网站| 午夜免费鲁丝| 午夜福利在线在线| 国产亚洲欧美在线一区二区| av片东京热男人的天堂| 日本熟妇午夜| 亚洲片人在线观看| 免费搜索国产男女视频| 级片在线观看| ponron亚洲| 亚洲九九香蕉| 免费看日本二区| 亚洲午夜精品一区,二区,三区| 久久婷婷人人爽人人干人人爱| 久久久久九九精品影院| 高潮久久久久久久久久久不卡| 好男人电影高清在线观看| 久久人妻福利社区极品人妻图片| 国产精品 欧美亚洲| 日韩欧美 国产精品| 国产精品av久久久久免费| 欧美另类亚洲清纯唯美| 亚洲国产高清在线一区二区三 | 丁香六月欧美| 亚洲精品一区av在线观看| 大型av网站在线播放| 最近最新中文字幕大全电影3 | 最近最新免费中文字幕在线| 91老司机精品| 波多野结衣高清无吗| 国产亚洲欧美在线一区二区| 老司机福利观看| 中文亚洲av片在线观看爽| 国产三级黄色录像| 少妇的丰满在线观看| 精品不卡国产一区二区三区| 欧美激情极品国产一区二区三区| 欧美色欧美亚洲另类二区| 又黄又爽又免费观看的视频| 免费av毛片视频| bbb黄色大片| 99国产精品一区二区三区| 精品国产美女av久久久久小说| 99在线视频只有这里精品首页| 搡老熟女国产l中国老女人| 日本 av在线| 久久香蕉激情| 欧美乱色亚洲激情| 午夜福利一区二区在线看| 麻豆成人午夜福利视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产免费男女视频| 国产私拍福利视频在线观看| 老司机午夜十八禁免费视频| 成人国产综合亚洲| 亚洲熟妇中文字幕五十中出| 久久久久久亚洲精品国产蜜桃av| 国产精品99久久99久久久不卡| 亚洲精品粉嫩美女一区| 精品日产1卡2卡| 亚洲国产高清在线一区二区三 | 日日夜夜操网爽| 美女免费视频网站| a在线观看视频网站| 精品久久久久久久人妻蜜臀av| 国产不卡一卡二| 给我免费播放毛片高清在线观看| 精品国产国语对白av| 黄色女人牲交| 国产极品粉嫩免费观看在线| 久久久久国产精品人妻aⅴ院| 亚洲第一青青草原| 国内精品久久久久久久电影| 变态另类成人亚洲欧美熟女| 90打野战视频偷拍视频| 国产av一区在线观看免费| 亚洲专区中文字幕在线| 免费人成视频x8x8入口观看| 亚洲精品av麻豆狂野| 久热这里只有精品99| 无人区码免费观看不卡| 精品少妇一区二区三区视频日本电影| 一个人免费在线观看的高清视频| 制服人妻中文乱码| 国产野战对白在线观看| 青草久久国产| 日韩国内少妇激情av| 国产精品免费一区二区三区在线| x7x7x7水蜜桃| 黄网站色视频无遮挡免费观看| 黑丝袜美女国产一区| 91大片在线观看| 美女免费视频网站| 国产成人av教育| 国产亚洲av高清不卡| 日韩欧美免费精品| 国内精品久久久久精免费| 欧美在线黄色| 国产麻豆成人av免费视频| 日本成人三级电影网站| 俄罗斯特黄特色一大片| 精品一区二区三区av网在线观看| 搡老岳熟女国产| 国产国语露脸激情在线看| 国产亚洲精品久久久久久毛片| 精品国产乱码久久久久久男人| 亚洲天堂国产精品一区在线| 亚洲国产毛片av蜜桃av| 一级作爱视频免费观看| 亚洲专区国产一区二区| 国产伦人伦偷精品视频| 亚洲精品av麻豆狂野| 国产一卡二卡三卡精品| 性欧美人与动物交配| 黄色 视频免费看| 母亲3免费完整高清在线观看| 亚洲午夜理论影院| 美女高潮喷水抽搐中文字幕| 校园春色视频在线观看| 男人舔奶头视频| 青草久久国产| 黄色成人免费大全| 黄频高清免费视频| 午夜日韩欧美国产| 色播在线永久视频| 一二三四社区在线视频社区8| 久久伊人香网站| 男女午夜视频在线观看| 久9热在线精品视频| 欧美日韩精品网址| 丝袜美腿诱惑在线| 99精品在免费线老司机午夜| 久久天堂一区二区三区四区| 人妻久久中文字幕网| 欧美+亚洲+日韩+国产| 桃红色精品国产亚洲av| 曰老女人黄片| 国产黄色小视频在线观看| 久久久久久久久中文| 欧美黄色片欧美黄色片| www国产在线视频色| 成人国产一区最新在线观看| 禁无遮挡网站| 一级黄色大片毛片| 狠狠狠狠99中文字幕| 99久久99久久久精品蜜桃| 久久精品国产亚洲av香蕉五月| 国产97色在线日韩免费| 好男人电影高清在线观看| 制服诱惑二区| 亚洲专区国产一区二区| 大型av网站在线播放| 一本综合久久免费| or卡值多少钱| 亚洲 欧美一区二区三区| 国产精品,欧美在线| 国产精品久久电影中文字幕| 日韩成人在线观看一区二区三区| 热99re8久久精品国产| 99久久久亚洲精品蜜臀av| 大香蕉久久成人网| 日本成人三级电影网站| 日本 欧美在线| 亚洲va日本ⅴa欧美va伊人久久| 亚洲狠狠婷婷综合久久图片| 国产1区2区3区精品| www.999成人在线观看| 一本精品99久久精品77| 一区福利在线观看| 成人三级做爰电影| 精品久久久久久,| 99久久精品国产亚洲精品| 大型黄色视频在线免费观看| 免费在线观看黄色视频的| 国产一级毛片七仙女欲春2 | 日韩欧美国产一区二区入口| 国产极品粉嫩免费观看在线| 老熟妇仑乱视频hdxx| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲中文字幕一区二区三区有码在线看 | 国产免费av片在线观看野外av| 亚洲avbb在线观看| 亚洲第一av免费看| 亚洲免费av在线视频| 99国产极品粉嫩在线观看| 18禁黄网站禁片午夜丰满| 少妇粗大呻吟视频| 亚洲色图 男人天堂 中文字幕| 午夜两性在线视频| 一边摸一边抽搐一进一小说| 国产精品一区二区精品视频观看| av超薄肉色丝袜交足视频| 曰老女人黄片| 欧美黑人巨大hd| 婷婷亚洲欧美| 男女视频在线观看网站免费 | 日韩 欧美 亚洲 中文字幕| 国产成人欧美在线观看| 午夜影院日韩av| 国产视频内射| 亚洲成人精品中文字幕电影| 十分钟在线观看高清视频www| 最好的美女福利视频网| 夜夜夜夜夜久久久久| 亚洲精品在线观看二区| 亚洲人成伊人成综合网2020| 精品人妻1区二区| 国产日本99.免费观看| 日韩欧美国产一区二区入口| 日本撒尿小便嘘嘘汇集6| 日本五十路高清| e午夜精品久久久久久久| 午夜免费成人在线视频| e午夜精品久久久久久久| 久久久国产欧美日韩av| 成年女人毛片免费观看观看9| 精品人妻1区二区| 欧美丝袜亚洲另类 | 国产精品亚洲一级av第二区| 又大又爽又粗| 熟女电影av网| 亚洲五月婷婷丁香| 999久久久精品免费观看国产| 国产亚洲av嫩草精品影院| 国产午夜福利久久久久久| 成人手机av| 欧美日韩黄片免| 亚洲色图 男人天堂 中文字幕| 99久久久亚洲精品蜜臀av| 国产麻豆成人av免费视频| 午夜福利在线在线| 欧美中文日本在线观看视频| 丁香欧美五月| 欧洲精品卡2卡3卡4卡5卡区| 久久婷婷人人爽人人干人人爱| 中文在线观看免费www的网站 | 好男人在线观看高清免费视频 | 精品午夜福利视频在线观看一区| 男女下面进入的视频免费午夜 | 搡老岳熟女国产| 亚洲精品国产一区二区精华液| 国产av一区二区精品久久| 中文字幕最新亚洲高清| 十八禁网站免费在线| 亚洲 欧美一区二区三区| 午夜福利免费观看在线| av电影中文网址| 18禁美女被吸乳视频| 岛国在线观看网站| 最近最新免费中文字幕在线| 一区二区三区激情视频| 男人操女人黄网站| 好男人电影高清在线观看| av片东京热男人的天堂| 美女午夜性视频免费| 国产成人啪精品午夜网站| 色综合亚洲欧美另类图片| 国产又黄又爽又无遮挡在线| 国产精品一区二区三区四区久久 | 亚洲精品色激情综合| 可以在线观看毛片的网站| 久久久国产精品麻豆| 在线观看舔阴道视频| 亚洲精品中文字幕一二三四区| av免费在线观看网站| 午夜免费激情av| 老汉色∧v一级毛片| 午夜激情福利司机影院| 久久久国产欧美日韩av| √禁漫天堂资源中文www| 窝窝影院91人妻| 国产亚洲欧美在线一区二区| 99国产精品一区二区三区| 色精品久久人妻99蜜桃| 一区二区三区高清视频在线| 欧美色欧美亚洲另类二区| 精品国产一区二区三区四区第35| 国产人伦9x9x在线观看| 国产一区二区在线av高清观看| 亚洲国产欧美网| 午夜亚洲福利在线播放| 日本成人三级电影网站| 女警被强在线播放| 国产一区在线观看成人免费| 一本精品99久久精品77| 视频区欧美日本亚洲| 日韩国内少妇激情av| 午夜福利在线观看吧| 悠悠久久av| 97超级碰碰碰精品色视频在线观看| 麻豆国产av国片精品| 伊人久久大香线蕉亚洲五| 久久草成人影院| 国产精品久久久人人做人人爽| 一级毛片女人18水好多| av在线天堂中文字幕| 亚洲色图av天堂| 禁无遮挡网站| 2021天堂中文幕一二区在线观 | 国产精品二区激情视频| 精品国产国语对白av| 亚洲国产精品成人综合色| 黄网站色视频无遮挡免费观看| 在线观看日韩欧美| 无遮挡黄片免费观看| 色尼玛亚洲综合影院| 波多野结衣高清作品| netflix在线观看网站| 99国产精品一区二区蜜桃av| 亚洲精品国产区一区二| 首页视频小说图片口味搜索| 他把我摸到了高潮在线观看| 男男h啪啪无遮挡| 大型av网站在线播放| 国产一级毛片七仙女欲春2 | 在线观看一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 久久久精品欧美日韩精品| 欧美日韩精品网址| 欧美日韩福利视频一区二区| 国产在线观看jvid| 国产精品一区二区免费欧美| 久久久精品欧美日韩精品| 精品福利观看| 一进一出抽搐动态| 国产主播在线观看一区二区| 黄色毛片三级朝国网站| 日韩av在线大香蕉| 老汉色∧v一级毛片| 久久久国产成人精品二区| 91麻豆精品激情在线观看国产| 一区二区三区国产精品乱码| 人成视频在线观看免费观看| 午夜视频精品福利| 中出人妻视频一区二区| 亚洲男人天堂网一区| 香蕉久久夜色| 午夜福利视频1000在线观看| 亚洲真实伦在线观看| 日本在线视频免费播放| 亚洲精品国产精品久久久不卡| 熟妇人妻久久中文字幕3abv| 男人的好看免费观看在线视频 | 欧美日韩亚洲综合一区二区三区_| 少妇 在线观看| 国产极品粉嫩免费观看在线| 欧美黑人精品巨大| 免费一级毛片在线播放高清视频| 人成视频在线观看免费观看| 亚洲专区中文字幕在线| 亚洲国产精品久久男人天堂| 国产精品日韩av在线免费观看| 欧美精品啪啪一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜福利欧美成人| 免费在线观看日本一区| 69av精品久久久久久| 丝袜在线中文字幕| 精品久久蜜臀av无| 国产一卡二卡三卡精品| 美国免费a级毛片| 大香蕉久久成人网| 亚洲国产毛片av蜜桃av| 人妻丰满熟妇av一区二区三区| 日韩欧美免费精品| 999久久久精品免费观看国产| 国产成人欧美在线观看| 长腿黑丝高跟| 欧美日韩瑟瑟在线播放| 一本一本综合久久| 97人妻精品一区二区三区麻豆 | 久久精品aⅴ一区二区三区四区| 看片在线看免费视频| 欧美在线一区亚洲| 久久婷婷成人综合色麻豆| 精品日产1卡2卡| 久久香蕉国产精品| 亚洲专区国产一区二区| 97超级碰碰碰精品色视频在线观看| 亚洲成人国产一区在线观看| 国产免费男女视频| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品久久成人aⅴ小说| 国产精品野战在线观看| 久久精品91无色码中文字幕| av中文乱码字幕在线| 日本一本二区三区精品| 狠狠狠狠99中文字幕| 久久久久久九九精品二区国产 | 一级黄色大片毛片| 亚洲免费av在线视频| 国产精品一区二区三区四区久久 | 国产精品精品国产色婷婷| 久久香蕉精品热| 久久久国产成人精品二区| 国产欧美日韩精品亚洲av| 老熟妇乱子伦视频在线观看| 精品一区二区三区四区五区乱码| 午夜免费观看网址| 国产免费男女视频| 99精品久久久久人妻精品| 草草在线视频免费看| 色综合婷婷激情| 大香蕉久久成人网| 精品国产一区二区三区四区第35| 亚洲av中文字字幕乱码综合 | 日本免费一区二区三区高清不卡| www.熟女人妻精品国产| 听说在线观看完整版免费高清| 午夜两性在线视频| 国产v大片淫在线免费观看| 国产黄a三级三级三级人| 国产黄片美女视频| 精品一区二区三区四区五区乱码| 久久久久久久久久黄片| 757午夜福利合集在线观看| 女警被强在线播放| 免费av毛片视频| 久久精品亚洲精品国产色婷小说| 好男人在线观看高清免费视频 | 久久国产亚洲av麻豆专区| 日韩高清综合在线| 巨乳人妻的诱惑在线观看| 中文字幕高清在线视频| 99riav亚洲国产免费| 波多野结衣av一区二区av| 人妻丰满熟妇av一区二区三区| 母亲3免费完整高清在线观看| 国产三级在线视频| 国产精品电影一区二区三区| 国产精品久久视频播放| 2021天堂中文幕一二区在线观 | 亚洲精品中文字幕一二三四区| 亚洲 国产 在线| 精品第一国产精品| 亚洲av熟女| 人人妻,人人澡人人爽秒播| 国内精品久久久久精免费| 精品久久久久久久末码|