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

    Regional Features and Seasonality of Land–Atmosphere Coup ling over Eastern China

    2018-04-08 10:59:11ChujieGAOHaishanCHENShanleiSUNBeiXUVictorONGOMASiguangZHUHediMAandXingLI
    Advances in Atmospheric Sciences 2018年6期

    Chujie GAO,Haishan CHEN?,Shanlei SUN,Bei XU,Victor ONGOMA,Siguang ZHU,Hedi MA,and Xing LI

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

    2College of Atmospheric Science,Nanjing University of Information Science and Technology,Nanjing 210044,China

    3Department of Meteorology,South Eastern Kenya University,P.O.Box 170,Kitui 90200,Kenya

    1. Introduction

    In the early 21st century,the Global Land–Atmosphere Coupling Experiment presented a basic concept for studying the effects of land surface factors on atmospheric predictability based on model intercomparison(Koster et al.,2004,2006).The strength of land–atmosphere coupling is used to denote the potential contribution of land surface states to the variabilities of atmospheric circulation and climate variables such as precipitation and temperature.Soil moisture(SM)acts as a reserve of water and energy because soil retains wet and dry conditions longer than the atmosphere.Therefore,the soil has the potential to affect the atmospheric conditions through different coupling mechanisms(Koster and Suarez,2001;Wu and Dickinson,2004;Seneviratne et al.,2006b;Wei et al.,2008;Spennemann and Saulo,2015).Land–atmosphere coupling,in terms of the degree of SM that affects precipitation and surface temperature(ST),has been extensively studied(Seneviratne et al.,2006a;Zhang et al.,2008a,2008b,2009,2011;Wei and Dirmeyer,2012;Ruscica et al.,2014;Guillod et al.,2015;Tuttle and Salvucci,2016).In fact,we need to consider the impact of SM on evapotranspiration(ET)and further ST,rather than the atmospheric control of ET and ST,to better understand the land–atmosphere coupling.

    The impact of SM on ET is considered to be the first segment of the connection from land surface states to the atmosphere(Dirmeyer,2011).First of all,SM is an important source of the atmospheric water.Through ET,a major component of the continental water cycle,SM returns approximately 60%of the total precipitation that falls on land back to the atmosphere(Oki and Kanae,2006).Secondly,SM also plays an important role in the land surface energy balance,since ET consumes more than half of the total solar energy absorbed by the land surface(Trenberth et al.,2009).Thus,SM can evidently influenc the thermal condition of the land surface through ET,consequently affecting the atmosphere(Zhang and Zuo,2011;Gao et al.,2014).So,the relationship between SM and ET is a key process in the land–atmosphere interface and is also the focus of our study.

    In addition,the key regions or“hot spots”of land–atmosphere coupling are unevenly distributed.By analyzing the mean states of multi-model simulations in summer(Koster et al.,2004),those“hot spots”mainly appear in transitional climates between dry and wet zones,such as the central Great Plains of North America,the Sahel,and India.For each of the individual models,the land–atmosphere coupling intensity and spatial pattern are distinctly different,indicating that there are also numerous uncertainties in the SM feedback.More studies have further con fi rmed that the SM effects on climate change are mainly detected in the dry–wet transition zones(Wei et al.,2008;Dirmeyer,2011;Hua et al.,2013).The reason is that SM and ET are not simply linearly related.For example,SM does not constrain ET if soil water content fully meets the land surface evaporation and demand of plants(i.e.,an“extreme wet condition”).On the other hand,if the soil is too dry(i.e.,an“extreme dry condition”),then ET stops completely(Seneviratne et al.,2010;Bellucci et al.,2015).SM does not control ET under either extreme wet or dry conditions.For this reason,most of the inferred effects of SM on climate systems(strong land–atmosphere coupling)are induced by its constraints on ET under the transitional conditions between dry and wet zones.

    Compared toSM–precipitation coupling,SM–temperature coupling is more distinct(Zhang et al.,2011;Zhang and Zuo,2011;Wu and Zhang,2013).For instance,lower SM usually results in less ET and weaker evaporative cooling,which leads to higher ST(Mei and Wang,2012;Zittis et al.,2013).Meanwhile,another interpretation of this negative correlation exists:high ST acts as a primary forcing factor,thereby increasing ET and drying the soil(Nicholls,2004;Cai and Cowan,2008).This hypothesis suggests that high(low)ST causes dry(wet)soil.In summary,SM–ST relationships can be used to classify the conditions for ET processes as water-lim ited(SM constrains ET)or energy-lim ited(ET is controlled by atmospheric backgrounds,and no land–atmosphere coupling exists under such a condition)(Koster et al.,2009;Seneviratne et al.,2010;Yin et al.,2014a),which is also closely related to dry–wet climate conditions.

    Previous studies have emphasized the spatial pattern of land–atmosphere coupling and identi fi ed“hot spots”with strong SM–climate feedback.However,less attention has been paid to seasonal changes of land–atmosphere coupling(Wei and Dirmeyer,2012;Wei et al.,2012),especially in China,where complex spatial and seasonal climate variations exist.In this study,we attempt to identify hot spots of land–atmosphere coupling and explore the seasonal changes of such coupling over China.

    The remaining parts of this work are organized as follows:Section 2 introduces the datasets and methods;section 3 presents the relationships among SM,ST and ET over the study area;sections 4 and 5 discuss and summarize the fi ndings,respectively.

    2. Data and methods

    2.1. Study area

    China has a vast territory and is a typical monsoon region,exhibiting varying climatic conditions(Guo et al.,2003;He et al.,2007).Figure 1 shows eastern China(east of 95°E)and its major geographical regions.In general,the elevation decreases from west to east,and a typical East Asian monsoon climate system prevails across this area(Jiang et al.,2015).Additionally,previous studies(e.g.,Ge and Feng,2009;Hu et al.,2015)have pointed out that more than 96%of China’s population lives in this area.Therefore,the local land–atmosphere coupling influences human life and needs to be explored extensively.

    2.2. Reanalysis dataset

    Because of a lack of observations,studies related to SM on larger spatial and longer time scales are still very limited.However,the recent development and availability of various satellite-derived and reanalysis SM datasets make it possible to understand relevant scientific questions more comprehensively.The uncertainties in satellite-derived SM data and the weaknesses of current quality controls(Yin et al.,2014b)are likely to introduce some uncertainties into our results.According to Zuo and Zhang(2007,2009),the spring SM of the European Centre for Medium-Range Weather Forecasts(ECMWF)ERA-40 dataset(1948–2002)can reproduce the temporal and spatial features of observed SM in eastern China reasonably well.Moreover,Zhang et al.(2008c)compared the ERA product with other multi-source datasets(e.g.,observations,NCEP/NCAR,GSWP2 and CLM)and found that ERA-40 can better capture the interannual and interdecadal variabilities of observed SM.

    Recently,the ECMWF released a new product:ERAInterim(Balsamo et al.,2015).In this dataset,the Tiled ECMWF Scheme used for the surface exchanges within the Land Surface Model(van den Hurk et al.,2000)is the same as that employed by ERA-40.The assimilation method,however,was upgraded from the 3D variational method of ERA-40 to a 4D variational approach in ERA-Interim(Dee et al.,2011).To further reduce the biases between the ERA-simulated variables and observations,the observed relative humidity and temperature are continuously used to adjust the SM(Douville et al.,2000;Mahfouf et al.,2000).By comparing several sets of reanalysis data(e.g.,ERA-Interim,MERRA,JRA,CFSR,and NCEP),Liu et al.(2014)noted that ERA-Interim performs best in reproducing the spatial and temporal features of the SM over eastern China.This dataset also captures the major characteristics of precipitation and evaporation,which are the two main factors affecting SM.Overall,ERA-Interim is preferable to other datasets in representing the spatial and temporal character-istics of the actual SM.Thus,this SM dataset is selected for this study,which is available at http://www.ecmw f.int/en/research/climate-reanalysis/era-interim/land.The soil consists of four layers of thickness,i.e.,7,21,72 and 189 cm from top to bottom,with a horizontal resolution of 1°×1°.Considering the greater accuracy of the near-surface ERA SM and its stronger interaction with ET(Zhang et al.,2008c;Zuo and Zhang,2009),we select the first layer of ERA SM for further analyses.Besides,the surface latent heat flux and top layer soil temperature are chosen to represent ET and ST in our study.The study period is from March 1979 to February 2014.

    Fig.1. Topography(units:m)and major geographical regions in eastern China.The geographical regions are based on Zhao(1995).The elevation data are from the Global 30 Arc-Second Elevation(GTOPO30-10min)dataset:https://lta.cr.usgs.gov/GTOPO30.

    2.3. Model simulation

    To enhance confidence in our results,a group of numerical experiments simulated by version 4.0 of NCAR’s Community Land Model(CLM 4.0,http://www.cesm.ucar.edu/models/cesm1.1/clm/)are used.This model is a welldeveloped land surface model and has been widely used for studies related to land surface processes(Chen et al.,2010;Xiong et al.,2011;Zhu et al.,2013).Recent studies(e.g.,Li et al.,2011;Lai et al.,2014)have stated that CLM 4.0 is capable of capturing the spatiotemporal characteristics of SM from in situ observations over China.The land cover type,soil color,texture,sand ratio and other information come from the land characteristic parameter data of the model,and the soil column is vertically represented by 15 layers(0.007,0.028,0.062,0.119,0.212,0.366,0.620,1.038,1.728,2.865,4.739,7.830,12.925,21.327,and 35.178 m).The global land surface conditions(e.g.,SM,ST,ET and heat fluxes)in the last five decades are produced by offline simulations,which are based on the global 1°×1°and 3-h atmospheric forcing dataset including near surface meteorological elements(i.e.,2-m temperature,wind speed,specific hum idity,precipitation,surface pressure,downward longwave and shortwave radiation)from 1948 to 2010 developed at Princeton University(Sheffield et al.,2006).The forcing dataset is available at http://hydrology.princeton.edu/data.pgf.php.The third layer SM(6.2 cm—close to the selected layer of ERA SM), first layer soil temperature and ET from March 1979 to February 2010 are selected for our study.

    2.4. Methods

    Our analyses are mainly based on Pearson correlation analysis.The monthly correlation coefficient between SM and ET(SM–ET relationship—hereafter referred to asRSM?ET)can be used to describe the extent to which SM affects ET(Dirmeyer et al.,2009;Dirmeyer,2011).A positiveRSM?ETsuggests that ET changes are primarily controlled by SM,e.g.,higher(lower)SM causes more(less)ET.This usually happens in arid regions under water-lim ited conditions.On the contrary,a negativeRSM?ETimplies that ET is mainly affected by the atmospheric environmental variables(e.g.,temperature,humidity and wind speed)and controls the changes of SM.In other words,increased(decreased)ET tends to decrease(increase)SM.Such a negativeRSM?ETtypically occurs in hum id regions under energy-limited conditions.For this reason,RSM?ETis often employed to determine the direction of land–atmosphere interactions,i.e.,how the land affects the atmosphere,and vice versa.Another widely used approach to diagnose land–atmosphere coupling is the correlation coefficient between ST and ET(RST?ET)(Zittis et al.,2013).This is also a measure of the interaction between soil and temperature through ET(Ruscica et al.,2014).RST?ETis generally positive and negative underenergy-limited and water-limited conditions,respectively.For example,whenRST?ETis positive,higher(lower)temperature(rather than SM)leads to more(less)ET.In contrast,a negativeRST?ETindicates that increased(decreased)ET results in lower(higher)temperature,which is very likely induced by an SM anomaly.To sum up,when both of the positiveRSM?ETand negativeRST?ETare significant,it is likely that SM regulates the ET process,thereby affecting ST.This situation suggests that there is a strong land–atmosphere coupling.

    In this study,all data are linearly detrended before perform ing statistical analyses.

    3. Results

    3.1.Seasonal mean SM and ST

    Figure 2 shows the seasonal mean states of SM in the two datasets.In general,both the ERA and CLM datasets are in agreement that the soil dries from southeast to northwest,and some differences(e.g.,the magnitude,location and tim ing of the maximum SM)exist among the four seasons,i.e.,spring(March–May,MAM),summer(June–August,JJA),autumn(September–November,SON)and winter(December–February,DJF).For ERA,the SM over North China is generally lower than 22%,with a dry center located west of Hetao(bend of the Yellow River),and partially reaches 25%during summer in the east.Northeast China maintains a relatively wet condition(over 25%)throughout the year.Meanwhile,the seasonal variability is larger in the southern part of the country.For the upstream region of the Yangtze River and central China,the SM in spring–autumn and winter are above and below 31%,respectively.The SM over East China and South China ranges from 28%and 25%during winter to 31%during summer,respectively.In particular,Southwest China is a dry center(<22%)during spring but a wet center(>31%)during summer and autumn.CLM basically reproduces the spatial patterns of seasonal SM.However,there are some systematic biases:wetter conditions in wet regions and drier conditions in dry regions,which is consistent with the fi ndings of Lai et al.(2014).Besides,the seasonal variability of SM for CLM is relatively low.Most obviously,there is a dry center near Hetao(<10%)and two wet centers in East China and Northeast China(>37%)that maintain their magnitudes and positions during the four seasons.Moreover,other regions demonstratecorresponding seasonalvariations,butmuch less so than in the ERA dataset.A ll in all,the dry–wet patterns from the two datasets match the dry and wet climate divisions in China over recent decades,insofar as humid/semi-hum id and arid/sem i-arid regions are located in southern and northern parts of China,and the Northeast China is an exceptional sem i-hum id region(Ma and Fu,2005;Wu et al.,2005;Zhang et al.,2016).

    Fig.2. Seasonal mean SM for the period 1979–2013 in the ERA(7 cm,top row)and CLM(6.2 cm,bottom row)datasets(units:volume%).

    Fig.3. As in Fig.2 but for ST(units:°C).

    Figure 3 presents the climatological ST of each season.The spatial distributions and seasonal changes of the ST in the two datasets are in close agreement.Generally,ST decreases from south to north.It is warmer in summer(ranging from 10°C to30°C)butcolderin winter(ranging from?15°C to 20°C).Comparing Figs.2 and 3,the ability of CLM to reproduce the ST climatology is apparently better than that of SM.

    The above results show that complex spatial and seasonal changes exist in both SM and ST.In the follow ing sections,the connections between SM and ST(through ET)are further investigated.These investigations are expected to reveal the possible pathway of land–atmosphere interactions and the“hot spots”of land–atmosphere coupling over eastern China.

    3.2. Spatial distributions and seasonal changes of RSM?ETand R ST?ET

    As shown in Fig.4,the extent to which SM affects ET,i.e.,theRSM?ET,is mainly opposite to the SM distribution in each season(Fig.2).Generally,positive and negativeRSM?ETare roughly distributed in the dry north and wet south,respectively.significant positive values are found over North China during the four seasons for both datasets,indicating that SM evidently affects ET because of the relatively dry soil condition.For other regions,differences exist not only in different seasons,but also in the different datasets.Firstly,for ERA,Northeast China exhibits lowerRSM?ETvalues,especially during summer and winter.In contrast,the seasonal changes in negative values over hum id/sem i-hum id regions are much more complex.For example,the Yangtze River basin shows significant negativeRSM?ETduring spring and summer,whereas its downstream region(East China)shows positive values during autumn and winter.Such a reversal in the changes also appears in South and Southwest China:theRSM?ETover South China is negative in summer but positive in winter,and theRSM?ETover Southwest China changes from positive to negative between winter/spring and summer/autumn.For the CLM dataset,Northeast China shows significant negativeRSM?ETthroughout the year,except for the winter season.Besides,the seasonal variations ofRSM?ETfor hum id/sem i-hum id regions in CLM are much weaker:East China and South China maintain negative values during autumn/winter.In particular,the significant positive values over Southwest China are only shown in spring.The negativeRSM?ETindicates ET is not constrained by SM under wet soil conditions,and these differences between the two datasets may be due to the systematic biases and weaker seasonal changes in the SM simulation,as discussed in section 3.1.Overall,the positiveRSM?ET,a necessary(but not suffi-cient)condition for land–atmosphere coupling(Dirmeyer et al.,2009),is found over North China throughout the year and over Southwest China during spring in both datasets.

    Figure 5 demonstrates the spatial distribution ofRST?ET,which is basically opposite to that ofRSM?ET(Fig.4).For ERA,positive values are mostly found in the humid/sem ihum id regions over the south and northea st,implying that the energy supply is the primary factor causing the ET anomaly.If the ST is higher than normal,the ET increases and dries the soil.Most of this positiveRST?ETcan be sustained for a whole year.On the contrary,theRST?ETover North China exhibits relatively larger seasonality.The values are mainly negative during the four seasons,but the most significant and robustly negative ones only appear in summer.Moreover,a reverse change inRST?ETis also shown over Southwest China:it is negative in spring and changes to positive during summer and autumn.CLM basically reproduces theRST?ETpattern of ERA.

    3.3. Seasonal evolution of land–atmosphere coupling

    Fig.4. Spatial distributions of correlation coefficients between SM and ET(RSM?ET)in the ERA(top row)and CLM(bottom row)datasets.The dotted grid points are significant at the 95%confidence level.All data are detrended.

    Fig.5. As in Fig.4 but for the correlation coefficients between ST and ET(RST?ET).

    To further investigate the relationships betweenRSM?ET/RST?ETand SM,we produce scatterplots of the two correlation coefficients at each grid point over eastern China for the ERA and CLM datasets(Fig.6).To reduce the influenc of cold conditions(e.g.,ice,snow and frozen soil)in high latitude areas(Fig.2),only the summer pattern is chosen for analysis.It is found that bothRSM?ETandRST?ETare linearly related to SM.Under dry soil conditions,the grid points with positiveRSM?ETare primarily accompanied by negativeRST?ET,which implies that SM affects ET and ST,and a strong land–atmosphere coupling exists.With increased SM,the positiveRSM?ETis reduced and even becomes negative.Meanwhile,the negativeRST?ETis turning positive.These fi ndings indicate that the influenc of SM on ET and ST is reduced,and the land–atmosphere coupling intensity is weakened,with the transition from arid to hum id regions.It is worth mentioning that when SM is relatively lower,a close linear relationship exists betweenRSM?ETandRST?ET,implying that the coupling changes with SM can be re fl ected by each of them.However,with higher SM,this relationship becomes weak,and suggests that ST may have different infl uences onRSM?ETandRST?ET,which can be found in their seasonal cycles and is further discussed in the follow ing text.

    Fig.6. Relationship between RSM?ET and RST?ET in summer at each grid point over eastern China for the(a)ERA and(b)CLM dataset.The different colors denote different SM ranges(units:volume%).A ll data are detrended before the correlation coefficient is calculated.

    Fig.7. SM intra-annual variability(multi-year average of maximum SM minus minimum SM within a year)in the(a)ERA and(b)CLM dataset(units:volume%).The four selected sub-areas are shown in(a):area I(21°–28°N,97°–104°E);area II(21°–33°N,106°–122°E);area III(35°–50°N,96°–121°E);and area IV(40°–53°N,123°–135°E).

    Additionally,SM differs not only by region but also by season(Fig.2),as does the pathway of land–atmosphere interactions(i.e.,SM affects ST and vice versa;Figs.4 and 5).To further exam ine the intra-annual variability of SM,the mean difference between maximum and m inimum SM within a year is shown in Fig.7.Large variations(>4%)are demonstrated in the southern half,southeastern North China and Northeast China for the ERA dataset.This is possibly due to the influenc of the transformation of the monsoon system(Jiang et al.,2015).Notably,Southwest China is a large SM-variation center.For the CLM dataset,as discussed in section 3.1,the seasonal variability is relatively lower,especially in the Yangtze River basin and Northeast China.Still,the largest intra-annual change is exhibited in Southwest China.Therefore,considering the evident seasonal variations inRSM?ET,RST?ETand SM,we first take area I(Southwest China:21°–28°N,97°–104°E)as an example to explore the monthly change at the regional scale.We also select the humid/semi-humid area II(including South China,East China and Central China:21°–33°N,106°–122°E),arid/sem i-arid area III(North China:35°–50°N,96°–121°E)and sem i-hum id area IV(Northeast China:40°–53°N,123°–135°E)for the comparative analyses.These four sub-areas are shown in the boxed areas in Fig.7a,and the dry–wet definitions are based on Ma and Fu(2005).

    In Southwest China,precipitation during the wet season(May–October)accounts for about 70%to 80%of the annual total(Zhang et al.,2014).Accordingly,as shown in Figs.8a and b,the SM in area I rapidly increases beginning in May and remains at a high level(>30%)from June to November.During the months with higher SM,negativeRSM?ETand positiveRST?ETare significant(p<0.05),indicating that energy primarily controls ET.During the dry season(November–April),the soil does not receive enough water;thus,SM decreases to slightly below 25%in March.Meanwhile,theRSM?ETandRST?ETbecome positive and negative,respectively.The alternating relationship signs imply that the feedback from SM to the atmosphere is enhanced(shifts from an energy-lim ited condition during the wet season to a water-lim ited condition during the dry season),particularly in early spring(March and April),characterized by themaximumRSM?ET.Sim ilarly,in area II,RSM?ETincreases andRST?ETdecreases as the soil turns from wet to dry in the ERA dataset(Fig.8c).For the CLM dataset,however,the weaker intra-annual SM variability in area II leads to a smaller corresponding change inRSM?ET(Fig.8d).

    Fig.8. Monthly variations of SM(blue lines;units:volume%),RSM?ET(red lines)and RST?ET(green lines)in areas(a,b)I,(c,d)II,(e,f)III and(g,h)IV from March to February for the ERA(left)and CLM(right)datasets.The dashed lines are significant at the 95%confidence level.A ll data are detrended before the correlation coefficient is calculated.

    In area III,the significant(p<0.05)positiveRSM?ETand lowerSM throughouttheyear indicate thatthisareaisundera water-lim ited condition in all seasons for both datasets(Figs.8eand f).However,thesigni fi cant(p<0.05)negativeRST?ET(strong land–atmosphere coupling)occurs only from May to August and becomes positive from September to April.The conditions in area IV are more complicated(Figs.8g and h),as compared with the other three areas.RSM?ETchanges with the ERA SM,butRST?ETis basically positive throughout the year.The CLM SM,by contrast,is high,and weaker variability exists during the 12 months,with theRSM?ETandRST?ETchanging differently.

    To further understand the monthly changes of the correlations in areas III and IV(Fig.8),we additionally exam ine the regional differences in ST,which also exhibit seasonal variations among the four sub-areas(Fig.3).The annual ST cycles based on the ERA and CLM datasets are shown in Fig.9.Areas I and II show low intra-annual ST variabilities,whereas areas III and IV show larger ones(bottom row of Fig.9).In areas III and IV,the ST is close to and often below 0°C from November to March,indicating that an insufficient energy supply exists for ET.In contrast,areas I and II are relatively warm throughout the year.These fi ndings provide a potential explanation for the abnormal positiveRST?ETobserved during the cold seasons in areas III and IV(Fig.8).

    To sum up,the above results show that differences exist in the land–atmosphere coupling in both space and seasonality.The differences are re fl ected in the spatial patterns ofRSM?ETandRST?ETand their seasonal variations.In consideration of the strong coupling only existing under the condition with positiveRSM?ETand negativeRST?ET,we additionally specify an indicator,which is expressed as the absolute value of the positiveRSM?ETmultiplied by the negativeRST?ET,to identify the“hot spots”.The results for both the ERA and CLM datasets are shown in Fig.10.For spring,the hot spots are in North China and Southwest China;during summer and autumn,however,these hot spots are primarily in North China.In addition,land–atmosphere coupling barely exists in winter.

    4. Discussion

    In this study,we explore the land–atmosphere coupling by using the correlations between SM and ET and those between ST and ET.North China(area III)is characterized by arid/sem i-arid regions;here,SM is maintained at a low level,and the SM–ET relationship is positive throughout the year(Figs.2,8e and f),which indicates that water-limited conditions predominate throughout the year.However,the coupling is strong in summer,relatively weak in spring and autumn,and eventually disappears in winter(Fig.10).These features may be associated with the ST seasonality,which is high in summer and low in winter(area III,Fig.9).During the cold seasons(e.g.,winter),the land surface is primarily coveredw th ice,snow and frozen soil,due to the low temperatures(<0°C,Figs.3 and 9).This situation i sbound to hinder the direct physical process between SM and the atmosphere.As a result,it is not difficult to understand why stronger coupling is always accompanied by higher temperatures(i.e.,sufficient energy supply),especially during summer.This suggests that a water-limited condition is necessary but not sufficient for land–atmosphere coupling in the relatively dry regions of China(i.e.,those at high latitude);and seasonal coupling variations are closely related to energy supply(as reflected by ST).In Southwest China(area I),the water supply is not sustainable because of the special climatic conditions(Zhang et al.,2014).Thus,the soil is generally drier during the spring than in other seasons(Figs.2,8a and b),which leads to a strong land–atmosphere coupling(Fig.10).However,in other humid/semi-hum id regions,over the southern part of eastern China(e.g.,area II)and Northeast China,the soil is moist throughout the year and the temperature is extremely low during local dry periods(October–February;Figs.8e,g and 9),respectively.Therefore,these areas exhibit no coupling(Fig.10).The results from the CLM 4.0 simulation are generally consistent with the regional features and seasonal changes from the ERA dataset.However,differences exist over some areas.Speci fi cally,CLM is based on a different model,which means the forcing dataset and descriptions of physical processes are different from ERA.This situation makes the differences in the results between the two datasets inevitable.

    Fig.9. Annual cycle(top row)of ST in areas I,II,III and IV,from March to February,in the(a)ERA and(b)CLM dataset.The solid dots and short lines(bottom row)are the annual mean,maximum and m inimum ST,respectively.Units:°C.

    Fig.10. Absolute value of positive RSM?ET multiplied by negative RST?ET.The dotted grid points indicate that both of the two correlations are significant at the 95%confidence level.All data are detrended.

    Previous studies have demonstrated that abnormal SM can affect the atmospheric general circulation mainly through altering the surface energy balance and land surface thermal forcing of the atmosphere,and further lead to precipitation anomalies in China(Zuo and Zhang,2007,2016;Zhang and Zuo,2011;Meng et al.,2014).Zhang and Zuo(2011)found that anomalous spring SM from the lower and middle reaches of the Yangtze River Valley to North China is highly correlated with the summer precipitation in China.Abnormally wet spring soil tends to increase the surface evaporation and lower the surface air temperature.As a result,the reduced temperature in late spring reduces the land–sea thermal contrast,resulting in a weakened East Asian monsoon.SM is de fi nitely an important factor re fl ecting the land–atmosphere interaction,and ET is a key component of the land–atmosphere feedback cycle.Our results further map the“hot spots”of land–atmosphere coupling over eastern China.In addition,we emphasize the seasonal variabilities of these hot spots.In the hum id/sem i-hum id regions of China,strong land–atmosphere coupling mainly happens in the relatively drier seasons,which is consistent with the conclusions of Wei et al.(2012)and Wei and Dirmeyer(2012),who showed that SM–precipitation coupling is stronger during dry periods for climatologically wet regions.However,in arid/sem i-arid regions of China,studies should focus on the relatively warmer seasons when significant negative SM feedbacks exist with relatively higher temperature(Zhang and Dong,2010;Zhang et al.,2011).

    Our study has certain lim itations.For example,land cover exerts an important effect on the exchanges of energy,water,carbon dioxide and other greenhouse gases between the land surface and atmosphere(Yin et al.,2015);thus,land-use and land-cover change(LUCC)m ight affect processes between the land and atmosphere.Moreover,remarkable LUCC has occurred alongside rapid econom ic developmentover thepast 30 years.Our study does not consider LUCC,which might lead to some uncertainties in the results.Although the analytical process may be subjective,it still provides a reference to further study land–atmosphere interactions.It is difficult to directly measure the strength of land–atmosphere coupling in the real world,and the specific land feedback on atmospheric circulation(e.g.,monsoon intensity)still needs to be explored in depth.With an increasing number of tools avail-able to determ inehow the land affects theclimate,researchers will have more opportunities to improve their understanding of climate predictions.

    5. Conclusions

    This study attempts to identify the“hot spots”of land–atmosphere coupling over eastern China using reanalysis data and simulation output.Besides the spatial distribution,we also examine the seasonal evolution of the coupling.Based on correlation analysis,we identify two hot spots:North China and Southwest China.The strongest couplings are found in summer over North China and in spring over Southwest China.

    Further analysis suggests that the coupling over North China is closely associated with temperature(partly re fl ecting the energy supply).In North China,the dry soil condition(water-limited)is maintained throughout the year;and summer is relatively warmer comparing to other seasons,which means a sufficient energy supply for ET,consequently leading to the strongest land–atmosphere coupling.On the other hand,the SM(water supply)plays a decisive role in the coupling over humid/semi-hum id regions.In Southwest China,the seasonal climate conditions are warm due to the lower latitude,meaning that the energy supply is sufficient throughout the year.However,the precipitation varies across the four seasons.By the end of the dry season(i.e.,spring),the SM is evidently lower,constraining ET.

    The uncertainties of this study should be noted.We use a high-quality reanalysis dataset(ERA-Interim)and model simulation(CLM 4.0)to maxim ize the robustness of the results,and the results from both datasets generally agree with each other,albeit with some differences observed in some areas.The methods applied in this study also have some limitations.Therefore,the fi ndings of this study should be regarded as prelim inary,and additional studies are needed.The specific feedback mechanisms from the land to the atmosphere over the aforementioned“hot spots”require more thorough investigation.

    Acknowledgements.This work was jointly supported by the National Natural Science Foundation of China(Grant Nos.41625019 and 41605042),the Natural Science Foundation of Jiangsu Province,China(Grant No.BK20151525),and the Priority Academ ic Program Development of Jiangsu Higher Education Institutions.

    Balsamo,G.,and Coauthors,2015:ERA-Interim/Land:A global land surface reanalysis data set.Hydrology and Earth System Sciences,19,389–407,https://doi.org/10.5194/hess-19-389-2015.

    Bellucci,A.,and Coauthors,2015:Advancements in decadal climate predictability:The role of nonoceanic drivers.Rev.Geophys.,53(2),165–202,https://doi.org/10.1002/2014RG 000473.

    Cai,W.J.,and T.Cowan,2008:Evidence of impacts from rising temperatureon inflows to theMurray-Darling Basin.Geophys.Res.Lett.,35,L07701,https://doi.org/10.1029/2008GL 033390.

    Chen,H.S.,M.M.Xiong,and W.Y.Sha,2010:Simulation of land surface processes over China and its validation Part I:Soil temperature.Scientia Meteorologica Sinica,30(5),621–630,https://doi.org/10.3969/j.issn.1009-0827.2010.05.008.(in Chinese)

    Dee,D.P.,and Coauthors,2011:The ERA-Interim reanalysis:Con fi guration and performance of the data assimilation system.Quart.J.Roy.Meteor.Soc.,137,553-597,https://doi.org/10.1002/qj.828.

    Dirmeyer,P.A.,2011:The terrestrial lmoistureclimatecoupling.Geophys.Res.Lett.,38,L16702,https://doi.org/10.1029/2011GL048268.

    Dirmeyer,P.A.,C.A.Schlosser,and K.L.Brubaker,2009:Precipitation,recycling,and land memory:An integrated analysis.Journal of Hydrometeorology,10,278–288,https://doi.org/10.1175/2008JHM 1016.1.

    Douville,H.,P.Viterbo,J.-F.Mahfouf,and A.C.M.Beljaars,2000:Evaluation of the optimum interpolation and nudging techniques for soil moisture analysis using FIFE data.Mon.Wea.Rev.,128,1733–1756,https://doi.org/10.1175/1520-0493(2000)128<1733:EOTOIA>2.0.CO;2.

    Gao,C.J.,H.S.Chen,B.Xu,and G.Zeng,2014:Possible relationships among South China Sea SSTA,soil moisture anomalies in southwest China and summer precipitation in eastern China.Journal of Tropical Meteorology,20(3),228–235,https://doi.org/10.16555/j.1006-8775.2014.03.005.

    Ge,M.L.,and Z.M.Feng,2009:Population distribution of China based on GIS:Classi fi cation of population densities and curve of population gravity centers.Acta Geographica Sinica,64(2),202–210,https://doi.org/10.3321/j.issn:0375-5444.2009.02.007.(in Chinese)

    Guillod,B.P.,B.Orlowsky,D.G.Miralles,A.J.Teuling,and S.I.Seneviratne,2015:Reconciling spatial and temporal soil moisture effects on afternoon rainfall.Nature Communications,6,6443,https://doi.org/10.1038/ncomms7443.

    Guo,Q.Y.,J.N.Cai,X.M.Shao,and W.Y.Sha,2003:Interdecadal variability of East-Asian summer monsoon and its impact on the climate of China.Acta Geographica Sinica,58(4),569–576,https://doi.org/10.3321/j.issn:0375-5444.2003.04.011.(in Chinese)

    He,J.H.,J.H.Ju,Z.P.Wen,J.M.Lu,and Q.H.Jin,2007:A review of recent advances in research on Asian monsoon in China.Adv.Atmos.Sci.,24(6),972–992,https://doi.org/10.1007/s00376-007-0972-2.

    Hu,L.L.,Y.L.Liu,Y.H.Ren,L.J.Yu,and C.Qu,2015:Spatial change of population density boundary in mainland China in recent 80 years.Journal of Remote Sensing,19(6),928–934,https://doi.org/10.11834/jrs.20155016.(in Chinese)

    Hua,W.J.,H.S.Chen,S.G.Zhu,S.L.Sun,M.Yu,and L.M.Zhou,2013:Hotspots of the sensitivity of the land surface hydrological cycle to climate change.Chinese Science Bulletin,58(30),3682–3688,https://doi.org/10.1007/s11434-013-5846-7.

    Jiang,J.,D.B.Jiang,and Y.H.Lin,2015:Monsoon area and precipitation over China for 1961–2009.Chinese Journal of Atmospheric Sciences,39(4),722–730,https://doi.org/10.3878/j.issn.1006-9895.1410.14195.(in Chinese)

    Koster,R.D.,and M.J.Suarez,2001:Soil moisture mem-ory in climate models.Journal of Hydrometeorology,2,558–570,https://doi.org/10.1175/1525-7541(2001)002<0558:SMM ICM>2.0.CO;2.

    Koster,R.D.,S.D.Schubert,and M.J.Suarez,2009:Analyzing the concurrence of meteorological droughts and warm periods,with implications for the determination of evaporative regime.J.Climate,22,3331–3341,https://doi.org/10.1175/2008JCLI2718.1.

    Koster,R.D.,and Coauthors,2004:Regions of strong coupling between soil moisture and precipitation.Science,305,1138–1140,https://doi.org/10.1126/science.1100217.

    Koster,R.D.,and Coauthors,2006:GLACE:The global land atmosphere coupling experiment.Part I:Overview.Journal of Hydrometeorology,7,590–610,https://doi.org/10.1175/JHM 510.1.

    Lai,X.,J.Wen,S.X.Cen,H.Q.Song,H.Tian,X.K.Shi,Y.He,and X.Huang,2014:Numerical simulation and evaluation study of soil moisture over China by using CLM 4.0 model.Chinese J.Atmos.Sci.,38(3),499–512,https://doi.org/10.3878/j.issn.1006-9895.1401.13194.(in Chinese)

    Li,M.X.,Z.G.Ma,and G.Y.Niu,2011:Modeling spatial and temporalvariationsin soilmoisture in China.ChineseScience Bulletin,56(17),1809–1820,https://doi.org/10.1007/s11434-011-4493-0.

    Liu,L.,R.H.Zhang,and Z.Y.Zuo,2014:Intercomparison of spring soil moisture among multiple reanalysis data sets over eastern China.J.Geophys.Res.Atmos.,119,54–64,https://doi.org/10.1002/2013JD020940.

    Ma,Z.G.,and C.B.Fu,2005:Decadal variations of arid and semi-arid boundary in China.Chinese Journal of Geophysics,48(3),519–525,https://doi.org/10.3321/j.issn:0001-5733.2005.03.008.(in Chinese)

    Mahfouf,J.-F.,P.Viterbo,H.Douville,A.C.M.Beljaars,and S.Saarinen,2000:A revised land-surface analysis scheme in the integrated forecasting system.ECMWF Newsletter,88,8–13.

    Mei,R.,and G.L.Wang,2012:Summer land-atmosphere coupling strength in the United States:Comparison among observations,reanalysis data,and numerical models.Journal of Hydrometeorology,13,1010–1022,https://doi.org/10.1175/JHM-D-11-075.1.

    Meng,L.,D.Long,S.M.Quiring,and Y.J.Shen,2014:Statistical analysis of the relationship between spring soil moisture and summer precipitation in East China.International Journal of Climatology,34(5),1511–1523,https://doi.org/10.1002/joc.3780.

    Nicholls,N.,2004:The changing nature of Australian droughts.Climatic Change,63(3),323–336,https://doi.org/10.1023/B:CLIM.0000018515.46344.6d.

    Oki,T.,and S.Kanae,2006:Global hydrological cycles and world water resources.Science,313,1068–1072,https://doi.org/10.1126/science.1128845.

    Ruscica,R.C.,A.A.Srensson,and C.G.Men′endez,2014:Hydrological links in Southeastern South America:Soil moisture memory and coupling within a hot spot.InternationalJournalofClimatology,34,3641–3653,https://doi.org/10.1002/joc.3930.

    Seneviratne,S.I.,D.Luthi,M.Litschi,and C.Schar,2006a:Landatmosphere coupling and climate change in Europe.Nature,443,205–209,https://doi.org/10.1038/nature05095.

    Seneviratne,S.I.,T.Corti,E.L.Davin,M.Hirschi,E.B.Jaeger,I.Lehner,B.Orlowsky,and A.J.Teuling,2010:Investigating soil moisture-climate interactions in a changing climate:A review.Earth-Science Reviews,99,125–161,https://doi.org/10.1016/j.earscirev.2010.02.004.

    Seneviratne,S.I.,and Coauthors,2006b:Soil moisture memory in AGCM simulations:Analysis of Global Land-Atmosphere Coupling Experiment(GLACE)Data.Journal of Hydrometeorology,7,1090–1112,https://doi.org/10.1175/JHM 533.1.

    Sheffield,J.,G.Goteti,and E.F.Wood,2006:Development of a 50-year-high-resolution global dataset of meteorological forcing for land surface modeling.J.Climate,19(13),3088–3111,https://doi.org/10.1175/JCLI3790.1.

    Spennemann,P.C.,and A.C.Saulo,2015:An estimation of the land-atmosphere coupling strength in South America using the Global Land Data assimilation System.InternationalJournalofClimatology,35,4151–4166,https://doi.org/10.1002/joc.4274.

    Trenberth,K.E.,J.T.Fasullo,and J.Kiehl,2009:Earth’s global energy budget.Bull.Amer.Meteor.Soc.,90(3),311–323,https://doi.org/10.1175/2008BAMS2634.1.

    Tuttle,S.,and G.Salvucci,2016:Empirical evidence of contrasting soil moisture-precipitation feedbacks across the United States.Science,352(6287),825–828,https://doi.org/10.1126/science.aaa7185.

    van den Hurk,B.J.J.M.,P.Viterbo,A.C.M.Beljaars,and A.K.Betts,2000:Offline validation of the ERA40 surface scheme.ECMWF Technical Memorandum No.295,1–42.

    Wei,J.F.,and P.A.Dirmeyer,2012:Dissecting soil moistureprecipitation coupling.Geophys.Res.Lett.,39,L19711,https://doi.org/10.1029/2012GL053038.

    Wei,J.F.,P.A.Dirmeyer,and Z.C.Guo,2008:Sensitivities of soil wetness simulation to uncertainties in precipitation and radiation.Geophys.Res.Lett.,35,L15703,https://doi.org/10.1029/2008GL034494.

    Wei,J.F.,P.A.Dirmeyer,M.G.Bosilovich,and R.G.Wu,2012:Water vapor sources for Yangtze River Valley rainfall:Climatology,variability,and implications for rainfall forecasting.J.Geophys.Res.,117,D05126,https://doi.org/10.1029/2011JD016902.

    Wu,L.Y.,and J.Y.Zhang,2013:Asymmetric effects of soil moisture on mean daily maximum and m inimum temperaturesovereastern China.Meteor.Atmos.Phys.,122,199–213,https://doi.org/10.1007/s00703-013-0284-2.

    Wu,S.H.,Y.H.Yin,D.Zheng,and Q.Y.Yang,2005:Aridity/hum idity status of land surface in China during the last three decades.Science in China Series D:Earth Sciences,48(9),1510–1518,https://doi.org/10.1360/04yd0009.

    Wu,W.R.,and R.E.Dickinson,2004:Time scales of layered soil moisture memory in the context of land-atmosphere interaction.J.Climate,17,2752–2764,https://doi.org/10.1175/1520-0442(2004)017<2752:TSOLSM>2.0.CO;2.

    Xiong,M.M.,H.S.Chen,and M.Yu,2011:Simulation of land surface processes over China and its validation.Part II:Soil moisture.Scientia Meteorologica Sinica,31(1),1–10,https://doi.org/10.3969/j.issn.1009-0827.2011.01.001.(in Chinese)

    Yin,D.Q.,M.L.Roderick,G.Leech,F.B.Sun,and Y.F.Huang,2014a:The contribution of reduction in evaporative cooling to higher surface air temperatures during drought.Geophys.Res.Lett.,41,7891–7897,https://doi.org/10.1002/2014 GL062039.

    Yin,J.F.,X.W.Zhan,Y.F.Zheng,J.C.Liu,C.R.Hain,and L.Fang,2014b:Impact of quality control of satellite soil moisture data on their assimilation into land surface model.Geo-phys.Res.Lett.,41(20),7159–7166,https://doi.org/10.1002/2014gl060659.

    Yin,J.F.,and Coauthors,2015:An assessment of impacts of land-cover changes on root-zone soil moisture.Int.J.Remote Sens.,36(24),6116–6134,https://doi.org/10.1080/01431161.2015.1111539.

    Zhang,C.J.,Y.M.Liao,J.Q.Duan,Y.L.Song,D.P.Huang,and S.Wang,2016:The progresses of dry-wet climate divisional research in China.ClimateChangeResearch,12(4),261–267,https://doi.org/10.12006/j.issn.1673-1719.2015.191.(in Chinese)

    Zhang,J.Y.,and W.J.Dong,2010:Soil moisture influenc on summertime surface air temperature over East Asia.Theor.Appl.Climatol.,100,221–226,https://doi.org/10.1007/s00704-009-0236-4.

    Zhang,J.Y.,W.C.Wang,and J.F.Wei,2008a:Assessing landatmosphere coupling using soil moisture from the Global Land Data Assimilation System and observational precipitation.J.Geophys.Res.,113,D17119,https://doi.org/10.1029/2008JD009807.

    Zhang,J.Y.,W.C.Wang,and L.R.Leung,2008b:Contribution of land-atmosphere coupling to summer climate variability over the contiguous United States.J.Geophys.Res.,113,D22109,https://doi.org/10.1029/2008JD010136.

    Zhang,J.Y.,W.C.Wang,and L.Y.Wu,2009:Land-atmosphere coupling and diurnal temperature range over the contiguous United States.Geophys.Res.Lett.,36,L06706,https://doi.org/10.1029/2009GL037505.

    Zhang,J.Y.,L.Y.Wu,and W.J.Dong,2011:Land-atmosphere coupling and summer climate variability over East Asia.J.Geophys.Res.,116,D05117,https://doi.org/10.1029/2010 JD014714.

    Zhang,R.H.,and Z.Y.Zuo,2011:Impact of spring soil moisture on surface energy balance and summer monsoon circulation over East Asia and precipitation in East China.J.Climate,24,3309–3322,https://doi.org/10.1175/2011JCLI4084.1.

    Zhang,W.J.,T.J.Zhou,and R.C.Yu,2008c:Spatial distribution and temporal variation of soil moisture over China Part I:Multi-data inter-comparison.Chinese Journal of Atmospheric Sciences,32,581–597,https://doi.org/10.3878/j.issn.1006-9895.2008.03.15.(in Chinese)

    Zhang,W.L.,J.Y.Zhang,and G.Z.Fan,2014:Dom inant modes of dry-and wet-season precipitation in southwestern China.Chinese Journal of Atmospheric Sciences,38(3),590–602,https://doi.org/10.3878/j.issn.1006-9895.2013.13156.(in Chinese)

    Zhu,S.G.,H.S.Chen,and J.Zhou,2013:Simulations of global land surface conditions in recent 50 years with three versions of NCAR Community Land Models and their comparative analysis.Transactions of Atmospheric Science,36(4),434-446,https://doi.org/10.13878/j.cnki.dqkxxb.2013.04.007.(in Chinese)

    Zittis,G.,P.Hadjinicolaou,and J.Lelieveld,2013:Land-Atmosphere coupling:The feedback of soil moisture into surface temperature in Eastern Mediterranean and middle East.Advances in Meteorology,ClimatologyandAtmospheric Physics,C.G.Helm is,and P.T.Nastos,Eds.,Springer,833–839,https://doi.org/10.1007/978-3-642-29172-2 117.

    Zuo,Z.Y.,and R.H.Zhang,2007:The spring soil moisture and the summer rainfall in eastern China.Chinese Science Bulletin,52,3310–3312,https://doi.org/10.1007/s11434-007-0442-3.

    Zuo,Z.Y.,and R.H.Zhang,2009:Temporal and spatial features of the soil moisture in boreal spring in eastern China.Science in China Series D:Earth Sciences,52(2),269–278,https://doi.org/10.1007/s11430-009-0011-5.

    Zuo,Z.Y.,and R.H.Zhang,2016:influenc of soil moisture in Eastern China on the East Asian summer monsoon.Adv.Atmos.Sci.,33(2),151–163,https://doi.org/10.1007/s00376-015-5024-8.

    色吧在线观看| 日本熟妇午夜| 天堂av国产一区二区熟女人妻| 男女那种视频在线观看| 在线播放国产精品三级| 久热久热在线精品观看| 成人鲁丝片一二三区免费| 最近手机中文字幕大全| 国内精品宾馆在线| 春色校园在线视频观看| 国内精品宾馆在线| 午夜激情福利司机影院| 干丝袜人妻中文字幕| 午夜激情福利司机影院| 久久精品夜色国产| 免费搜索国产男女视频| 国产精品日韩av在线免费观看| 久久久久久国产a免费观看| 99在线人妻在线中文字幕| 床上黄色一级片| 国产亚洲午夜精品一区二区久久 | 又黄又爽又刺激的免费视频.| 亚洲精华国产精华液的使用体验| 最近2019中文字幕mv第一页| 在现免费观看毛片| 国产黄a三级三级三级人| 成人国产麻豆网| 午夜日本视频在线| av免费观看日本| 国产一区二区在线观看日韩| 亚洲精品乱码久久久v下载方式| 国产精品一二三区在线看| 在现免费观看毛片| 久久久久九九精品影院| 超碰av人人做人人爽久久| 天堂√8在线中文| www日本黄色视频网| 午夜日本视频在线| 亚洲乱码一区二区免费版| 日本免费一区二区三区高清不卡| 成人毛片a级毛片在线播放| 麻豆成人av视频| 97超碰精品成人国产| 尤物成人国产欧美一区二区三区| 高清在线视频一区二区三区 | 久99久视频精品免费| 能在线免费观看的黄片| 国产v大片淫在线免费观看| 亚洲精品456在线播放app| 国产精品一二三区在线看| 国产人妻一区二区三区在| av在线播放精品| 久久精品久久精品一区二区三区| 国产男人的电影天堂91| 九九久久精品国产亚洲av麻豆| 国产三级中文精品| 亚洲国产高清在线一区二区三| 嫩草影院入口| 国产成人精品婷婷| 亚洲av成人av| 国产精品一区二区三区四区久久| 精品久久久久久电影网 | 91久久精品国产一区二区三区| 国产精品三级大全| 男女下面进入的视频免费午夜| 国产精品福利在线免费观看| ponron亚洲| 国产老妇女一区| 久久精品熟女亚洲av麻豆精品 | 有码 亚洲区| 爱豆传媒免费全集在线观看| 精华霜和精华液先用哪个| 亚洲精品乱码久久久久久按摩| 三级男女做爰猛烈吃奶摸视频| 日韩成人伦理影院| 午夜免费男女啪啪视频观看| 啦啦啦韩国在线观看视频| 日韩欧美国产在线观看| 久久精品国产亚洲av天美| 特大巨黑吊av在线直播| 欧美丝袜亚洲另类| 一级黄片播放器| 精品一区二区三区视频在线| 国内精品宾馆在线| 精品久久久久久久末码| 毛片一级片免费看久久久久| av视频在线观看入口| 观看美女的网站| 国产精品日韩av在线免费观看| 麻豆国产97在线/欧美| 免费av毛片视频| 午夜a级毛片| 在线观看一区二区三区| 久久久午夜欧美精品| 色综合色国产| 伊人久久精品亚洲午夜| 国产伦在线观看视频一区| 中文字幕av在线有码专区| 午夜久久久久精精品| 国产精品蜜桃在线观看| 建设人人有责人人尽责人人享有的 | 亚洲人成网站在线观看播放| 国产成人a区在线观看| 免费黄网站久久成人精品| 国产乱人偷精品视频| 直男gayav资源| 久久精品久久久久久久性| 国产人妻一区二区三区在| 久久久久网色| 成人亚洲精品av一区二区| 秋霞在线观看毛片| av女优亚洲男人天堂| 乱系列少妇在线播放| 在现免费观看毛片| 色视频www国产| 2021少妇久久久久久久久久久| 丝袜美腿在线中文| 精品免费久久久久久久清纯| 搡老妇女老女人老熟妇| 91在线精品国自产拍蜜月| 综合色丁香网| 亚洲av电影在线观看一区二区三区 | 国产亚洲av嫩草精品影院| 大香蕉97超碰在线| 亚洲四区av| 一本一本综合久久| 中文乱码字字幕精品一区二区三区 | 午夜福利在线在线| 三级国产精品欧美在线观看| 国产免费又黄又爽又色| 天堂影院成人在线观看| 亚洲av成人精品一区久久| 国产精品麻豆人妻色哟哟久久 | 非洲黑人性xxxx精品又粗又长| 18禁动态无遮挡网站| 久久久国产成人免费| 男人和女人高潮做爰伦理| 热99re8久久精品国产| 久久久国产成人精品二区| 久久精品国产亚洲av天美| 精品久久久久久久久久久久久| 亚洲人成网站在线播| 色尼玛亚洲综合影院| 蜜桃久久精品国产亚洲av| 亚洲美女搞黄在线观看| av黄色大香蕉| 精品人妻一区二区三区麻豆| 亚洲精品一区蜜桃| 国产久久久一区二区三区| 中文字幕免费在线视频6| 亚洲国产成人一精品久久久| 亚洲精品乱码久久久v下载方式| 久久久久久久久久久丰满| 十八禁国产超污无遮挡网站| 国产大屁股一区二区在线视频| 亚洲,欧美,日韩| 国产黄片美女视频| 精品少妇黑人巨大在线播放 | 97超视频在线观看视频| 国产精品av视频在线免费观看| 国产精品av视频在线免费观看| 国产私拍福利视频在线观看| 亚洲自偷自拍三级| 亚洲自偷自拍三级| 国产精品三级大全| 欧美一区二区亚洲| 18禁在线无遮挡免费观看视频| 丰满少妇做爰视频| 99热6这里只有精品| 又爽又黄a免费视频| 18+在线观看网站| 搞女人的毛片| 国语自产精品视频在线第100页| 国产精品一区二区三区四区免费观看| 麻豆一二三区av精品| 黄色一级大片看看| 全区人妻精品视频| 22中文网久久字幕| 日韩中字成人| 国产亚洲5aaaaa淫片| 欧美日本亚洲视频在线播放| 亚洲国产精品sss在线观看| 九九久久精品国产亚洲av麻豆| 久久久精品94久久精品| 免费大片18禁| 亚洲精品影视一区二区三区av| 亚洲精品影视一区二区三区av| 简卡轻食公司| 麻豆乱淫一区二区| 国产v大片淫在线免费观看| 嘟嘟电影网在线观看| 成人亚洲精品av一区二区| 欧美变态另类bdsm刘玥| 国产成人一区二区在线| 亚洲在久久综合| 国产成人精品久久久久久| 亚洲乱码一区二区免费版| 最新中文字幕久久久久| 精品不卡国产一区二区三区| 国产女主播在线喷水免费视频网站 | 久久6这里有精品| 91精品伊人久久大香线蕉| 男插女下体视频免费在线播放| 99热这里只有是精品在线观看| 国产精品久久久久久久久免| 级片在线观看| 美女被艹到高潮喷水动态| 日日摸夜夜添夜夜添av毛片| 国产精品电影一区二区三区| 黄色配什么色好看| 国产av码专区亚洲av| 欧美又色又爽又黄视频| 97超视频在线观看视频| 成人特级av手机在线观看| 天天躁日日操中文字幕| 99久国产av精品国产电影| 国产乱来视频区| 99久久人妻综合| 啦啦啦韩国在线观看视频| 亚洲在久久综合| 久久鲁丝午夜福利片| 欧美高清成人免费视频www| 久久久久久久久久成人| 国产精品国产高清国产av| 国产欧美日韩精品一区二区| 三级经典国产精品| 精品国内亚洲2022精品成人| 国产激情偷乱视频一区二区| 人妻制服诱惑在线中文字幕| 精品国产一区二区三区久久久樱花 | 日日撸夜夜添| 日本猛色少妇xxxxx猛交久久| 久久6这里有精品| 好男人视频免费观看在线| 狂野欧美白嫩少妇大欣赏| 内地一区二区视频在线| 欧美成人免费av一区二区三区| 国产欧美日韩精品一区二区| 亚洲精品影视一区二区三区av| 亚洲av熟女| 日本免费a在线| 国产精品嫩草影院av在线观看| 国产成人aa在线观看| 成人亚洲精品av一区二区| 国产成人a区在线观看| 日日摸夜夜添夜夜添av毛片| 成人亚洲欧美一区二区av| 男女啪啪激烈高潮av片| 菩萨蛮人人尽说江南好唐韦庄 | av又黄又爽大尺度在线免费看 | 亚洲国产精品成人久久小说| 国产视频内射| 国产综合懂色| 国产又色又爽无遮挡免| 三级经典国产精品| 亚洲精品自拍成人| 少妇被粗大猛烈的视频| 97在线视频观看| 国产精华一区二区三区| 麻豆av噜噜一区二区三区| 亚洲精品一区蜜桃| 欧美人与善性xxx| 天堂中文最新版在线下载 | 色哟哟·www| 欧美成人一区二区免费高清观看| 国产视频内射| 国产免费又黄又爽又色| 国产精品.久久久| 男的添女的下面高潮视频| 美女黄网站色视频| 免费看美女性在线毛片视频| av国产久精品久网站免费入址| 成年免费大片在线观看| 可以在线观看毛片的网站| 国产伦理片在线播放av一区| 99视频精品全部免费 在线| 久久国产乱子免费精品| 亚洲精品国产成人久久av| 18禁在线无遮挡免费观看视频| 国产亚洲一区二区精品| 天堂网av新在线| 国产毛片a区久久久久| 亚洲激情五月婷婷啪啪| 亚洲伊人久久精品综合 | 一边摸一边抽搐一进一小说| 韩国av在线不卡| 女的被弄到高潮叫床怎么办| 人妻系列 视频| 国产成人精品一,二区| 亚洲国产精品合色在线| 久久精品人妻少妇| 亚洲av免费在线观看| 国产黄片视频在线免费观看| 欧美色视频一区免费| 又黄又爽又刺激的免费视频.| 特级一级黄色大片| 欧美精品一区二区大全| 99久久精品热视频| 69人妻影院| 超碰97精品在线观看| 国产淫语在线视频| 99视频精品全部免费 在线| eeuss影院久久| 欧美一区二区国产精品久久精品| 一区二区三区免费毛片| 亚洲一区高清亚洲精品| 波多野结衣巨乳人妻| 亚洲第一区二区三区不卡| 日韩制服骚丝袜av| 色播亚洲综合网| 欧美成人精品欧美一级黄| 看片在线看免费视频| 我的老师免费观看完整版| 亚洲欧美成人综合另类久久久 | 日韩一区二区视频免费看| 男的添女的下面高潮视频| 国产精品一区二区三区四区久久| videossex国产| 少妇人妻精品综合一区二区| 成人综合一区亚洲| 小说图片视频综合网站| 欧美bdsm另类| 两个人视频免费观看高清| 久久精品国产亚洲网站| 看片在线看免费视频| 插逼视频在线观看| 亚洲av成人精品一二三区| 深爱激情五月婷婷| 精品一区二区三区视频在线| 日韩欧美国产在线观看| 成年版毛片免费区| 国产午夜福利久久久久久| 色噜噜av男人的天堂激情| 久久6这里有精品| 久久久久久久久久久丰满| 国产淫片久久久久久久久| 国产精品一二三区在线看| 日韩欧美在线乱码| 国产精品一区二区性色av| eeuss影院久久| 成人特级av手机在线观看| 一区二区三区免费毛片| av黄色大香蕉| 国产探花极品一区二区| 国产av在哪里看| 99在线人妻在线中文字幕| 人妻系列 视频| 伦精品一区二区三区| 精品国产三级普通话版| 亚洲最大成人手机在线| 欧美变态另类bdsm刘玥| 高清午夜精品一区二区三区| 国产精品一区二区在线观看99 | 国内精品一区二区在线观看| 日韩制服骚丝袜av| 国产精品一区www在线观看| 国产精品.久久久| av国产免费在线观看| 亚洲精品,欧美精品| 国产黄色小视频在线观看| av黄色大香蕉| 亚洲精品久久久久久婷婷小说 | 亚洲av电影在线观看一区二区三区 | 国产视频内射| 赤兔流量卡办理| 色综合亚洲欧美另类图片| 简卡轻食公司| 亚洲在久久综合| 亚洲美女视频黄频| 亚洲精品影视一区二区三区av| 天堂影院成人在线观看| 大话2 男鬼变身卡| 国产亚洲91精品色在线| 中文乱码字字幕精品一区二区三区 | 国产黄片视频在线免费观看| 中文字幕制服av| 国产老妇女一区| 久久亚洲精品不卡| 国产真实乱freesex| 赤兔流量卡办理| 少妇人妻一区二区三区视频| 欧美zozozo另类| 高清毛片免费看| 亚洲欧美精品自产自拍| 超碰av人人做人人爽久久| 女的被弄到高潮叫床怎么办| h日本视频在线播放| 国产一区亚洲一区在线观看| 一级av片app| 国产av码专区亚洲av| 成人特级av手机在线观看| 女人被狂操c到高潮| h日本视频在线播放| 亚洲熟妇中文字幕五十中出| 免费在线观看成人毛片| 亚洲欧美日韩东京热| 国产伦精品一区二区三区视频9| 在线免费十八禁| 国产亚洲精品av在线| 精品一区二区免费观看| 少妇人妻精品综合一区二区| 日韩国内少妇激情av| 我的女老师完整版在线观看| 麻豆成人午夜福利视频| 久久亚洲国产成人精品v| 日韩精品有码人妻一区| 99久国产av精品国产电影| 久久久久久久久中文| 国产三级在线视频| 女人十人毛片免费观看3o分钟| 精品人妻一区二区三区麻豆| 三级国产精品欧美在线观看| 欧美性猛交黑人性爽| 亚洲成人久久爱视频| 免费av观看视频| 亚洲成人精品中文字幕电影| 天天躁日日操中文字幕| 又黄又爽又刺激的免费视频.| 午夜福利高清视频| 少妇裸体淫交视频免费看高清| 如何舔出高潮| 精品国产露脸久久av麻豆 | 免费看光身美女| 国产成人免费观看mmmm| 秋霞在线观看毛片| 国产成人freesex在线| 亚洲av免费在线观看| 日韩 亚洲 欧美在线| 熟妇人妻久久中文字幕3abv| 免费播放大片免费观看视频在线观看 | 成人漫画全彩无遮挡| 欧美另类亚洲清纯唯美| 欧美区成人在线视频| 久久久久久国产a免费观看| 久久精品国产亚洲网站| 看十八女毛片水多多多| 亚洲自偷自拍三级| 久久99热6这里只有精品| 一个人看的www免费观看视频| 黑人高潮一二区| 建设人人有责人人尽责人人享有的 | 少妇的逼好多水| 国产视频内射| 夫妻性生交免费视频一级片| or卡值多少钱| 欧美另类亚洲清纯唯美| 身体一侧抽搐| 成人国产麻豆网| 久久久欧美国产精品| 香蕉国产在线看| videossex国产| 久久狼人影院| 日韩不卡一区二区三区视频在线| 欧美国产精品va在线观看不卡| 国产精品不卡视频一区二区| 制服诱惑二区| 大香蕉97超碰在线| 人人妻人人爽人人添夜夜欢视频| 亚洲伊人久久精品综合| 国产精品99久久99久久久不卡 | 最近手机中文字幕大全| 国产一区二区激情短视频 | 久久国产精品大桥未久av| 欧美性感艳星| 国产乱人偷精品视频| 亚洲精品日韩在线中文字幕| videossex国产| 蜜桃在线观看..| 日本vs欧美在线观看视频| 性色av一级| 国产av码专区亚洲av| 制服人妻中文乱码| 老熟女久久久| 亚洲精品乱码久久久久久按摩| 国产精品久久久久久av不卡| 久久婷婷青草| 久久久久精品久久久久真实原创| 亚洲欧美一区二区三区国产| 国产一区二区激情短视频 | 亚洲中文av在线| 免费久久久久久久精品成人欧美视频 | 在线免费观看不下载黄p国产| 成人亚洲欧美一区二区av| 国产综合精华液| 男人添女人高潮全过程视频| 亚洲av成人精品一二三区| 99热全是精品| 插逼视频在线观看| 91在线精品国自产拍蜜月| 免费黄频网站在线观看国产| 日韩中字成人| 少妇熟女欧美另类| 肉色欧美久久久久久久蜜桃| 免费大片18禁| 人妻 亚洲 视频| 亚洲伊人久久精品综合| 丰满饥渴人妻一区二区三| 久久免费观看电影| 午夜福利影视在线免费观看| 伦精品一区二区三区| 日本猛色少妇xxxxx猛交久久| 亚洲国产精品国产精品| 日本与韩国留学比较| 夜夜骑夜夜射夜夜干| 亚洲人成77777在线视频| 亚洲av电影在线观看一区二区三区| 精品国产一区二区三区四区第35| 日本wwww免费看| 午夜精品国产一区二区电影| 亚洲精品久久久久久婷婷小说| 国产精品久久久久久精品电影小说| 国产一区二区激情短视频 | 久热这里只有精品99| www.色视频.com| 在线亚洲精品国产二区图片欧美| 免费高清在线观看视频在线观看| 欧美人与性动交α欧美软件 | www.色视频.com| 亚洲四区av| 免费大片黄手机在线观看| 国产国拍精品亚洲av在线观看| 午夜影院在线不卡| 蜜臀久久99精品久久宅男| 啦啦啦中文免费视频观看日本| 午夜免费男女啪啪视频观看| 成人国产麻豆网| 麻豆精品久久久久久蜜桃| 国产黄频视频在线观看| 久热久热在线精品观看| 美女主播在线视频| 少妇人妻精品综合一区二区| 永久网站在线| 少妇猛男粗大的猛烈进出视频| 日韩,欧美,国产一区二区三区| 欧美日韩视频高清一区二区三区二| 国产精品不卡视频一区二区| 色哟哟·www| 国产日韩欧美在线精品| 中国三级夫妇交换| 亚洲av在线观看美女高潮| 99热国产这里只有精品6| 亚洲精品成人av观看孕妇| 天天操日日干夜夜撸| 国产亚洲精品久久久com| www.av在线官网国产| 久久人人爽av亚洲精品天堂| 久久婷婷青草| 国产精品三级大全| 亚洲第一av免费看| 日韩中文字幕视频在线看片| 日本91视频免费播放| 国产成人免费观看mmmm| 在线观看www视频免费| 母亲3免费完整高清在线观看 | 大香蕉久久成人网| 亚洲欧美清纯卡通| 国产探花极品一区二区| 9色porny在线观看| 日韩在线高清观看一区二区三区| 青青草视频在线视频观看| 尾随美女入室| 狂野欧美激情性xxxx在线观看| 国产一区有黄有色的免费视频| 菩萨蛮人人尽说江南好唐韦庄| 天天躁夜夜躁狠狠躁躁| 日本与韩国留学比较| 少妇熟女欧美另类| 午夜日本视频在线| 国产精品一二三区在线看| 大香蕉97超碰在线| 夫妻性生交免费视频一级片| 国语对白做爰xxxⅹ性视频网站| 国产精品国产av在线观看| 日韩中文字幕视频在线看片| 丁香六月天网| 99热国产这里只有精品6| 香蕉精品网在线| 91国产中文字幕| 激情视频va一区二区三区| 亚洲精品久久久久久婷婷小说| 亚洲五月色婷婷综合| 丝袜人妻中文字幕| 51国产日韩欧美| 日韩av不卡免费在线播放| 成年人午夜在线观看视频| 大话2 男鬼变身卡| 侵犯人妻中文字幕一二三四区| 男女边摸边吃奶| 国产日韩欧美视频二区| 天堂中文最新版在线下载| 五月天丁香电影| 国产有黄有色有爽视频| 99久国产av精品国产电影| 亚洲国产精品一区二区三区在线| 婷婷色综合www| 久久人妻熟女aⅴ| 丝袜脚勾引网站| 亚洲一码二码三码区别大吗| 一本久久精品| 国产精品久久久久久久电影| av电影中文网址| 国产女主播在线喷水免费视频网站| 日本欧美国产在线视频| 国产成人精品久久久久久| 婷婷色综合www| 中文字幕人妻丝袜制服| 日本欧美国产在线视频| 久久免费观看电影| 99热这里只有是精品在线观看| 亚洲精品,欧美精品| 国产精品蜜桃在线观看| 久久这里只有精品19| 久久久久久久亚洲中文字幕| 一边亲一边摸免费视频| 18+在线观看网站|