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

    Changes of Precipitation and Extremes and the Possible Effect of Urbanization in the Beijing Metropolitan Region during 1960–2012 Based on Homogenized Observations

    2015-06-09 21:30:01LIZhenYANZhongweiTUKaiandWUHongyi
    Advances in Atmospheric Sciences 2015年9期

    LI Zhen,YAN Zhongwei,TU Kai?,and WU Hongyi

    1Key Laboratory of Regional Climate-Environment in Temperate East Asia,Institute of Atmospheric Physics,Beijing 100029

    2Beijing Meteorological Bureau,Beijing 100089

    Changes of Precipitation and Extremes and the Possible Effect of Urbanization in the Beijing Metropolitan Region during 1960–2012 Based on Homogenized Observations

    LI Zhen1,YAN Zhongwei1,TU Kai?1,and WU Hongyi2

    1Key Laboratory of Regional Climate-Environment in Temperate East Asia,Institute of Atmospheric Physics,Beijing 100029

    2Beijing Meteorological Bureau,Beijing 100089

    Daily precipitation series at 15 stations in the Beijing metropolitan region(BMR)during 1960–2012 were homogenized using the multiple analysis of series for homogenization method,with additional adjustments based on analysis of empirical cumulative density function(ECDF)regarding climate extremes.The cumulative density functions of daily precipitation series,the trends of annual and seasonal precipitation,and summer extreme events during 1960–2012 in the original and final adjusted series atBeijing station were comparatively analyzed to show the necessity and efficiency ofthe new method.Results indicate that the ECDF adjustments can improve the homogeneity of high-order moments of daily series and the estimation of climate trends in extremes.The linear trends of the regional-mean annual and seasonal(spring,summer,autumn,and winter) precipitation series are?10.16,4.97,?20.04,5.02,and?0.11 mm(10 yr)?1,respectively.The trends over the BMR increase consistently for spring/autumn and decrease for the whole year/summer;however,the trends for winter decrease in southern parts and increase in northern parts.Urbanization affects local trends of precipitation amount,frequency,and intensity and their geographical patterns.For the urban-influenced sites,urbanization tends to slow down the magnitude of decrease in the precipitation and extreme amount series by approximately?10.4%and?6.0%,respectively;enhance the magnitude of decrease in precipitation frequency series by approximately 5.7%;reduce that of extremes by approximately?8.9%;and promote the decreasing trends in the summer intensity series of both precipitation and extremes by approximately 6.8%and 51.5%,respectively.

    homogenization,urbanization,precipitation extremes,MASH,ECDF

    1.Introduction

    Precipitation is an important climatic element,essential for studying water resources and the water cycle,and for quantifying climate change.However,long-term climate data series are usually interrupted by various non-natural factors, such as changes to observation locations,the environment, instrumentation,observation practices,orthe algorithms used for calculating any particular climate variable(Manton et al., 2001;Tuomenvirta,2001;Yan and Jones,2008).The inhomogeneity caused by these factors could severely distort the true climate signal and introduce considerable bias to estimates of climate trends,variability,and extremes.As a result,various methods have been introduced for detecting and adjusting inhomogeneity in different climatic variables (Potter,1981;Easterling and Peterson,1995;Lund and Reeves,2002;Li et al.,2014).

    To date,several helpful temperature series homogenization techniques have been developed for many regions(Vincent et al.,2002;Trewin,2013).However,the homogenization of precipitation series is rarely performed well, possibly because:(1)daily precipitation is a discrete variable with some non-normal distribution and stronger temporal and spatial variability;(2)the causes of inhomogeneity in precipitation series are more complex;and(3)the density of observation sites for precipitation is relatively sparse (Li et al.,2008b).Despite this,some meaningful studies have been conducted(Alexandersson,1986;Rhoades and Salinger,1993;Gonz′alez-Rouco et al.,2001;Beaulieu et al., 2008,2009;Wang et al.,2010).Tuomenvirta(2001)suggested that the precipitation adjustments could be±40%for 55 annual precipitation series in the North Atlantic region during 1890–1990,which were systematically biased by in-strumentalchanges.Begertetal.(2005)found thatthe adjustment factors ranged from 0.5 to 1.6 for 12 precipitation series in Switzerland from 1864 to 2000,and the introduction of automatic measuring equipment led to systematically lower measurements(approximately 5%less than the mean level). McRoberts and Nielsen-Gammon(2011)found that the remaining inhomogeneity from changes in gauge technology and station location may be responsible for an artificial trend of 1%to 3%per century,based on full network estimated precipitation data.Wang et al.(2010)developed a new algorithm(“transPMFred”)to detect shifts in nonzero daily precipitation series recorded at some stations across Canada and used a quantile-matching algorithm to adjustshifts in nonzero daily precipitation series.

    For China,although the importance of homogeneous climate data has been considered in earlier works(Tao et al., 1991;Yan et al.,2001),most previous climate change analyses associated with precipitation have been carried out based on unadjusted data.In recent years,some authors have pointed out the inhomogeneity in precipitation series(Liu and Sun,1995;Jiang et al.,2008;Li et al.,2008b).Li et al. (2012)applied the standard normal homogeneity test method to monthly precipitation series at 110 stations throughout China during 1900–2009 and ultimately completed two gridded datasets(5?×5?,2?×2?)of homogenized monthly precipitation.These homogenized datasets may be more reliable than the originals for climate change analyses.However,there is still potential for homogenization of daily precipitation series,which is important for evaluating changes of climate extremes.Daily meteorological observations have been used to assess changes in climate extremes or unusually anomalous weather fluctuations(Jones et al.,1999).As a result,a few homogenization methods,such as RHtests (Wang et al.,2010),have been developed for daily precipitation records.

    We have previously developed homogenized climate datasets based on multiple analysis of series for homogenization(MASH),including daily temperature(Li and Yan, 2010)and wind speed(Li et al.,2011a)series in the Beijing metropolitan region(BMR),and daily maximum/mean/ minimum temperature series over China(Li and Yan,2009). These results demonstrated MASH as a feasible technique for the homogenization of a large number of stations,especially in cases without complete metadata,such as in China.However,a problem arises from MASH where it is essential to adjust biases in the mean levels of the time series.To homogenize daily precipitation series,we developed an additional adjustment method to the relevant cumulative density function(CDF)dealing with climate extremes,which has been used in bias correction of model output.

    In the present study,we establish a homogenized daily precipitation dataset at 15 stations in the BMR for the period 1960–2012 based on a combination of the MASH and empirical cumulative density function(ECDF)method,quantifying any secular trends in precipitation series during the past 53 years based on homogenized observations.We subsequently identify the possible influence of urbanization on the amount,frequency,and intensity of summer precipitation and extremes,and their trends.The rest of this paper is organized as follows:The daily precipitation data and methods are described in section 2.The results,including the homogenization of daily precipitation series,changes in annual and seasonal precipitation,and summer precipitation extremes are illustrated,and the possible influences of urbanization in the BMR analyzed,in section 3.Further discussion and conclusions are presented in section 4.

    2.Data and methods

    2.1.Data

    Observed daily precipitation records from 20 meteorological stations in the BMR were collected and processed by the Information Center of the Beijing Meteorological Bureau. The observed values were quality controlled using some conventional procedures listed by Xu et al.(2013).Fifteen stations(Beijing–BJ,Shangdianzi–SDZ,Miyun–MY,Huairou–HR,Pinggu–PG,Changping–CP,Yanqing–YQ,Shunyi–SY, Chaoyang–CY,Tongxian–TX,Xiayunling–XYL,Daxing–DX,Fengtai–FT,Mentougou–MTG,Fangshan–FS)with few missing codes(a total of 184 days)were chosen,covering the study period of 1960–2012.Figure 1c shows the topography and geographical distribution of the 15 stations over the BMR.The metadata ofrelocation recordsateach station were also collected.All stations,apart from XYL,have relocated at least once and there were a total of 43 relocation records. There have been five relocations for BJ since 1960.The station was moved from Five Tower Temple to Daxing County on 1 January 1965,then to Zhanghua in Haidian County on 1 January 1969,to Daxing County on 1 July 1970,to Beiwa Road in Haidian District,a location closer to the city center on 1 January 1981,and back to the same suburban location in Daxing Country on 1 July 1997,with 38.8 km,44 km,24 km,22 km,and 22 km of horizontal moving distance,respectively.

    2.2.Methods

    2.2.1.Multiple Analysis of Series for Homogenization method

    The latest version(v3.03)of the MASH package(Szentimrey,1999)was used to detect and adjust possible break points due to relocations and other non-natural factors and homogenize the daily precipitation series in the BMR.Here we outline some of the basic ideas of MASH:(1)it is a relative homogeneity test procedure and does not assume that the reference series are homogeneous;(2)it includes a stepby-step iteration procedure where the role of series(candidate,reference)changes step by step in the course of the procedure;(3)additive(e.g.,for temperature)or multiplicative (e.g.,for precipitation)models can be used depending on the climatic elements;(4)it includes quality control and missing data completion procedures for monthly,seasonal,annual and daily data;(5)in the case of having monthly series for all the 12 months,the monthly,seasonal,and annual series canbe homogenized;(6)the daily inhomogeneity can be derived from the monthly inhomogeneity;and(7)metadata(possible dates of break points)can be used automatically.

    In the present study,the multiplicative model was applied to daily precipitation data underlying a quasi-lognormal distribution.A strict significance level for the test statistic (Monte Carlo method)was given at 0.01 to maintain the information in the original dataset as much as possible.The reference system(nine reference stations)of each candidate station was created based on their distances and correlation coefficients to the candidate station.The metadata of documented relocation records were used.

    2.2.2.Empirical Cumulative Density Function method

    The distribution-based homogenization method ECDF was introduced to the study.The idea of this method originates from bias correction techniques used in output adjustment of global or regional climate models for hydrological application or climate impact assessment(Piani et al.,2010; Yang et al.,2010;Dosio and Paruolo,2011).The basic concept of the method assumes that the series before and after the break point follows the same probability distribution. However,a non-parametric method,different from previous parametric methods cited above,was adopted to avoid uncertainties in the adjustment dealing with extremes,i.e.,only through the original ECDF,instead of parameter estimation by a certain probability distribution(e.g.,gamma).

    There are two circumstances in which MASH-based homogenized series can be additionally adjusted based on ECDF.If the break points detected by MASH agree with relocation records in each season,the MASH-based adjusted daily series will be additionally adjusted via ECDF,and thefinal adjusted daily series based on MASH+ECDF are obtained.In contrast,if no break point is detected,or the break points cannot be verified by metadata,the MASH-based series remain unchanged.

    The ECDF method includes the following steps:

    (1)Define the candidate sub-periods(CSPs)divided by break points,which are detected by MASH and consistent with metadata.CSPs are numbered starting from the most recent sub-period.

    (2)One or more reference station(s)is(are)chosen,if they have high correlation with the candidate station on the daily scale and no documented break points during the 10-year window centered on the relocation date of the candidate station,as suggested by Yan et al.(2010).Accordingly, the 10-year window is defined as the reference sub-periods (RSPs).RSPs are the same as CSPs,starting from the most recent sub-period.

    (3)Start the adjustment with the most recent inhomogeneity.The differences between the ECDFs of the reference and candidate series after the break point in RSP 1 are calculated after they are interpolated into the same intervals(1000 intervals between 0 and 1).If there is more than one reference station,their records will be combined together as one series.

    (4)The differences before the break point should be the same,if the candidate series is homogeneous.Based on this assumption,the differences are interpolated into the ECDF of the reference series in the RSP 2 before the break point and then interpolated back into the ECDF of the candidate series in the CSP 2,and the final adjustment values are obtained.

    (5)CSP 2 is homogeneous with respect to the latest subperiod(CSP 1)after subtracting the adjustment values.Repeat the process for all other CSPs sequentially.

    (6)Note that the adjustments are implemented separately for spring(March,April and May),summer(June,July and August;JJA),autumn(September,October and November; SON),and winter(December,January and February).

    To classify urban and rural sites regarding variability of summer precipitation,an empirical orthogonal function (EOF)analysis(Deng et al.,1989)was carried out on the standardized anomaly series of summer precipitation at the 15 stations.The linear trend and correlation were assessed to express long-term changes in each original and adjusted series during 1960–2012.A trend was considered statistically significant at the 95%confidence level and the slopes of trends were calculated by least-squares linear fitting.

    3.Results

    3.1.Homogenization of daily precipitation series

    First,the MASH method was applied to detect and adjust the inhomogeneity in daily precipitation series at the 15 stations in the BMR.The MASH method detected some documented changes associated with relocations.For annual series,the estimated inhomogeneity(ratio of original/adjusted) among the 15 stations varied between 0.92 and 1.15.Eight stations(53.3%of the total number)were inhomogeneous, with a total of 17 breakpoints and seven explained by relocation records.For the seasons,four,one,eight,and nine precipitation series were inhomogeneous and the estimated inhomogeneity ranged from 0.45 to 1.46,1 to 1.18,0.87 to 2.92,and 0.38 to 2.59 in spring,summer,autumn and winter, respectively.The number of break points in each season was different and 24.59%of them(of a total 51)were caused by relocations.For the BJ station,there were seven,zero,one, three,and four break points in the annual and seasonal series(Table 1).On the seasonal scale,only the break points in the summer of 1997 and the autumns of 1965 and 1969 coincided with relocation records.As a result,relocation was not a major reason for inhomogeneity in the precipitation series. The remaining break points may be due to factors including the new type of screens and instruments,introduction of automatic measurement equipment,and changes in observation times.

    Next,the ECDF method was used to additionally adjust the daily MASH-based precipitation series.To demonstrate how ECDF works in homogenization,the BJ station is selected as an example.Figure 2 shows the CSPs of summer and autumn precipitation series at BJ,and the corresponding ASPs of reference stations,according to the definition given in section 2.2.2.For the summer series,CSP 1(CSP 2)was July 1997 and JJA 1998–2012(JJA 1960–96 and June 1997);the reference stations were FT(0.85),CY(0.84),and DX(0.82);and ASP 1(ASP 2)was July and August 1997, JJA 1998–2001,and July 2002(July and August 1992,JJA 1993–96,and July 1997).For the autumn series,the CSPs were SON 1969–2012,1965–68,and 1960–64.For the break point in 1969/65,the reference stations were FT and CY/CP (0.68),and the corresponding RSPs were SON 1969–72 and 1965–68/SON 1965–69 and 1960–64,respectively.

    Figure 3 provides the CDFs of the original and final precipitation and precipitation>78 mm during 1960–2012 at BJ.It shows that the CDFs of daily precipitation≤0.8 mmwere almost identical in the original and adjusted daily precipitation series.The CDF of daily precipitation between 0.8 and 100.7 mm in the adjusted series is greater than that in the original series.However,the CDF of daily precipitation>100.7 mm in the adjusted series is less than that in the original series,indicating that the probability of extreme precipitation is adjusted towards the higher.For example,the CDF at 110 mm is 99.95%and 99.87%in the original and adjusted series,respectively.

    Table 1.Break points detected in the annual and seasonal series by multiple analysis of series for homogenization at Beijing station.

    To further understand the impact of inhomogeneity on trend estimation,we compared the annual/seasonal precipitation and summer precipitation extreme series(Fig.4)and their trends between the original and adjusted series(Table 2)at BJ.Here,the occurrence of an extreme precipitation (day)is recorded if daily precipitation exceeds the 90th percentile of the daily precipitation distribution during 1960–2012(Jones et al.,1999;Yan et al.,2002).The spring series was homogeneous with increasing trends of approximately 4.82 mm(10 yr)?1in the original and adjusted series(Fig. 4a).For the summer series,the precipitation was enhanced in some years during the early period after adjustment.The decreasing linear trend was slightly amplified from?23.22 to?24.01 mm(10 yr)?1due to some higher extreme values in the early period afteradjustment(Fig.4b).Forthe autumn series,the increasing trend was adjusted by 0.07 mm(10 yr)?1, which was less than the estimate of 5.84 mm(10 yr)?1in the original series(Fig.4c).In winter,both series show the same trends of approximately?0.23 mm(10 yr)?1,although there are some slight adjustments in the early period(Fig.4d).The adjusted annual series indicates a stronger decreasing trendof?13.66 mm(10 yr)?1than the original series[?12.79 mm(10 yr)?1](Fig.4e).The frequencies of summer precipitation extremes have slight differences of±1 day between the original and adjusted series(Fig.4f).The decreasing trend is?0.25 d(10 yr)?1after adjustment,more prominent than the estimate of?0.19 d(10 yr)?1in the original series.A moving two-sided t-test(n1=n2=10,degrees of freedom= 18)of the adjusted series indicates a sharp decrease(signif icant atα=0.05)in the mid-1990s in the annual and summer series,and the summer precipitation extremes series;and a rapid increase in the mid-1970s in the winter series.

    Table 2.Linear trends in annual and seasonal precipitation[units:mm(10 yr)?1]and summer precipitation extremes[units:d(10 yr)?1] at Beijing station during 1960–2012:Comparison between original and adjusted data.

    3.2.Geographical patterns of trends in the seasonal and annual precipitation

    Based on the final adjusted dataset,we calculated the trends in annual and seasonal precipitation and summer precipitation extremes averaged over the BMR.After adjustment,the spring and autumn precipitation series exhibited increasing trends of 4.97 and 5.02 mm(10 yr)?1,which were slightly less than the original series by?0.20 and?0.05 mm(10 yr)?1,respectively.However,the adjusted summer,winter,and annual series showed decreasing trends of?20.04,?0.11 and?10.16 mm(10 yr)?1,which were more prominent than the original values by?0.05,?0.03 and?0.33 mm(10 yr)?1,respectively.For the summer and annual precipitation,the inhomogeneity may be responsible for an artificial trend of 0.3%(10 yr)?1and 3.3%(10 yr)?1,respectively.The trend of the regional mean summer precipitation extremes was approximately?0.30 d(10 yr)?1in the original and adjusted series.The influence of homogenization on estimating the regional mean precipitation trend was notlarge,partly because some stationswere notinfluenced by relocation and partly because effects of local inhomogeneity compensated each other to some extent.The moving twosided t-tests indicated that there were more pronounced increasing trends in the annual and summer precipitation and summer precipitation extreme series during 1999–2012,fluctuating with linear variations of 173.1,127.8 and 3.2 d(10 yr)?1,respectively(Fig.5,left panels).

    The corresponding geographical patterns of annual and seasonal precipitation trends are also given(Fig.5,right panels).In spring,all 15 stations were found to have experienced increasing trends,especially in the northern part,indicating a southwest–northeast gradient.For TX,the increasing trend was reduced from 4.52 mm(10 yr)?1to 1.12 mm(10 yr)?1after adjustment,with the largest adjustment amplitude among all stations(Fig.5g,right panels).For the summer series,the pattern demonstrates the southeast–northwest gradient,with the lowest value of?26.23 mm(10 yr)?1at YQ and the highest value of?11.52 mm(10 yr)?1at CY,which is within the Central Business District(CBD)of Beijing, Chaoyang District,with the highest urbanization level among all districts and counties(Zhang,2006).Away from the CBD site,the decreasing trends become stronger,with larger values in the rural area,implying an underlying effect of urbanization in Beijing(Fig.5h,right panel).In autumn, the increasing trends vary from 1.65 mm(10 yr)?1to 9.19 mm(10 yr)?1;greater in the southeast(highest at DX)and lessin the west(lowestatMTG)(Fig.5i,rightpanel).In winter,the magnitudes of the trends are small,varying between?0.49 and 0.35 mm(10 yr)?1,showing a northern-increase–southern-decrease dipole pattern(Fig.5j,right panel).The annual geographical pattern is similar to the summer,because the precipitation in summer plays a dominant role in the overall annual precipitation amount.The adjusted annual data exhibit decreasing trends ranging from?16.22 to?0.73 mm(10 yr)?1,presenting a southeast–northwestgradient(Fig.5k, right panel).The magnitude of the decreasing trend is again lowest at CY and highest atYQ.Forthe summer precipitation extreme(Fig.5i,right panel),the spatial pattern indicates a northwest–southeast gradient.The trends are negative at 13 stations,especially at YQ(significant decreasing),ranging from?0.6 to?0.2 d(10 yr)?1;whereas,they are positive by approximately 0.1 and 0.2 d(10 yr)?1at CY and TX,respectively,with the highest urbanization levels and fastest urbanization rates over the past 10 years.This pattern suggests a possible influence of urbanization.

    3.3.Possible influence of urbanization on summer precipitation and extremes

    It is well known that the trend and pattern of regional precipitation is influenced by many factors,including global and regional water cycles,climate variation,topography,urbanization,and land use/cover change(Song et al.,2014). Hand and Shepherd(2009)found that the north-northeastern regions of the metropolitan Oklahoma City were wetter than other regions,and precipitation modification by the urban environment may be more significant than agricultural/topographic influences on weakly forced days.Over India,there is a significant increasing trend in the frequency of heavy rainfall climatology over urban regions during the monsoon season(Kishtawal et al.,2010).Over the Pearl River Delta metropolitan regions of China,the urbanization signatures in strong precipitation are significantly different from those in weak precipitation over the urban areas(Li et al.,2011b).For the BMR,the long-term decreasing trend in precipitation is consistent with that in northern China and mainly attributed to changes in regional atmospheric circulation(Xu et al.,2006;Wang et al.,2008;Cong et al.,2010). However,the trend and pattern of precipitation and extremes may be influenced by urbanization(Sun and Shu,2007;Li et al.,2008a;Zhang et al.,2009a,b;Miao et al.,2011;Yang et al.,2013,2014;You et al.,2014;Zhang et al.,2014).Here, we take summer precipitation as an example to quantify the possible influences of urbanization on the amount,frequency, and intensity of precipitation and extremes.A relevant study for low-temperature precipitation in the cold season can be found in Han et al.(2014).

    An EOF analysis was applied to the standardized anomaly time series ofsummerprecipitation atthe 15 stations during 1960–2012 to help identify whether there are different regional regimes of climate variability in summer precipitation in association with urban or rural divisions in the BMR. Figure 6 shows the first two EOF patterns and time coefficient series.The first EOF mode explains as much as 70.14% of the total variance,showing a coherent phase(negative coefficients)over the whole region,and indicating very similar inter-annual variation and consistent decreasing trends of summer precipitation among all stations in the BMR(Figs. 6a and c).However,the coefficients are larger in the southeast BMR,where is more urbanized and has a greater urbanization rate than the surrounding areas;and smaller in the northwest,where there is a mountainous area,suggesting a possible combined influence of urbanization and topography on the summer precipitation pattern.Based on the EOF 1 mode,four stations(BJ,DX,CY,and FT)are classified as the urban group,and the two stations(MTG,FS)away from BJ (in the upwind of the urban area)and with similar elevations to BJ are defined as the rural group.The second EOF mode explains 7.17%of the total variance and there is slight linear decreasing trend in the time coefficient series.The pattern clearly demonstrates a west–east gradient of summer precipitation across the BMR,suggesting the possible influence of large-scale factors,such as mean sea level pressure(Figs.6b and d).

    Figure 7 showsthe geographicaldistribution ofthe annual mean amount,frequency,and intensity of summer precipitation,as well as those of extremes.To further investigate the differences in these six indices between the urban and rural regimes,and quantify the possible influence of urbanization, we calculated the mean time series for the four urban sites,as classified by the first EOF mode,and the two rural stations. Figure 8 shows the amount,frequency and intensity series of summer precipitation and extremes in urban and rural areas during 1960–2012.Table 3 lists the trends derived from the series in Fig.8 and the corresponding possible contributions of urbanization,which are defined as

    where Turand Truare the trends of the urban and rural series during 1960–2012,respectively.

    Summer precipitation in the BMR usually occurs when a westerly trough approaches,and is attributed to water vapor transportation through the southwesterly jet in the mid and lower troposphere and/or to the easterly wind at the bottom ofthe troposphere.Figures7a and d show thatthe higherprecipitation occurs at the interface between the mountainous area and the plain areas(mainly on the windward side,especially in the southwest and northeast)and downwind of the urban area(especially in the northeast),while lower values occur in the urban core region.This confirms the combined influence of topography and urbanization on the pattern of precipitation amount,as well as extreme precipitation amount, in summer.This result is consistent with the work of Yin et al.(2011).The annual summer precipitation amount/extreme amount is 425.46/322.83 mm in the urban region,approximately 4.2%/3.0%higher than that in the rural area.The decreasing trendsof summer precipitation and extreme amounts are?17.94 and?16.50 mm(10 yr)?1,respectively,which are approximately 10.4%/6.0%less pronounced than those in the rural area(Figs.8a and d;Table 3).This suggests that urbanization tends to slightly enhance summer precipitation and extreme amounts,and reduce the magnitude of the decreasing trend in the urban area.

    The geographical pattern of precipitation frequencies indicates a southeast–northwestgradient.The lowestvalues occur mainly in the urban core area and the southeast part of the city(Fig.7b).This pattern may be shaped by topography and urbanization;while the frequencies of extremes are greater in the urban region,especially at CY,than those in the surrounding area.This shows that the extremes may be more likely to occur in the urban area(Fig.7e).The frequencies of precipitation in the rural area are greater than those in the urban area, but the reverse of the extremes in certain years(Figs.8b and e).The frequencies of precipitation and extremes exhibit declining trends of?0.57 and?0.16 d(10 yr)?1,respectively, in the urban area,which are slightly more and less prominent than the rural area by?0.03 and 0.02 d(10 yr)?1,respectively.Correspondingly,the influence of urbanization is estimated to be 5.7%and?8.9%,respectively,showing that urbanization may reduce the magnitude of the decreasing trend ofprecipitation frequency,butenhance thatof extremes in the urban area.

    Table 3.Trends in the amount[units:mm(10 yr)?1],frequency[units:d(10 yr)?1],and intensity[units:mm d?1(10 yr)?1]series of summer precipitation and extreme precipitation in the urban and rural areas during 1960–2012 and the possible contributions of urbanization to the trends of these six indices(*significance level:α=0.05).

    For intensity,the patterns of summer precipitation and extremes are similar to those of precipitation amount.The highest intensities are located in the northeast part(downwind of the urban area)(Figs.7c and f).The intensities of precipitation and extremes exhibit decreasing trends at each site,with a lower center in the core and eastern partof the city.For extreme intensity,the decreasing trend is significant at the 0.05 significance level at urban sites[?1.78 mm d?1(10 yr)?1],while insignificant in rural sites[?1.17 mm d?1(10 yr)?1].The differences in the trends[?0.02 and?0.61 mm d?1(10 yr)?1]between the rural and urban categories were used to estimate the influence of urbanization on the intensity of precipitation and extremes,accounting for around 6.8%and 51.5%,respectively,of the decreasing trend ofthe ruralcategory,which may representlarge-scale change. These results suggest that urbanization tends to promote a decrease in intensity,especially of extremes.

    4.Discussion and conclusions

    A dataset of daily precipitation observations in the BMR during 1960–2012 was homogenized based on MASH and additionally adjusted by the ECDF method regarding climate extremes.The new daily precipitation dataset provides a more reasonable basis for quantifying precipitation climate changes and the effect of urbanization on summer precipitation and extremes.The major conclusions from the present analysis are summarized with discussion as follows.

    Some break points associated with relocation records can be detected by MASH.However,the numbers are different in each season and significantly less than those in the temperature and wind speed series shown in our previous studies(Li and Yan,2010;Li et al.,2011a),suggesting that the precipitation series is less significantly influenced by the unnatural changes in the observation systems.

    The combination of the MASH and ECDF methods provides a new way to homogenize precipitation series on the daily scale from the perspective of both mean level and higher-order moment and shows the potential for improving and quantifying estimates of climate changes dealing with extreme precipitation and events,especially on the local scale.For BJ,the probability of extreme precipitation(>100.7 mm)increased after adjustment and the decreasing trend of summer precipitation extremes was amplified from?0.19 to?0.25 d(10 yr)?1.

    Inhomogeneity has little influence on the regional average long-term climate trends.When averaged over the BMR, the adjusted precipitation data exhibit an increasing trend of 4.82 and 5.77 mm(10 yr)?1for spring and autumn,respectively,and a decreasing trend of?24.01,?0.23 and?13.66 mm(10 yr)?1for summer,winter,and the whole year during 1960–2012,respectively.The summer extreme precipitation event series exhibits a declining trend of approximately?0.30 d(10 yr)?1.The increasing climatic trends since 1999 are notable,both in the summer/annual precipitation and summer precipitation extreme series.

    The geographic patterns of trends in precipitation for the whole year and the four seasons are different.In spring and autumn,precipitation of all the 15 stations shows a consistently increasing trend,with the highest values in the northern and southeast part of the BMR.For summer and the whole year,the patterns of decreasing trends demonstrate a southeast–northwest gradient,with the lowest value in the northwest and the highest value in the core urban area,implying an underlying effect of urbanization in Beijing.However, there are decreasing trends in the southern and increasing trends in the northern part for winter,revealing a northernincrease–southern-decrease pattern.For summer precipitation extremes,there are greater magnitudes of decreasing trends at the rural sites and smaller magnitudes of decreasing trends or even positive trends at the urban sites.

    Urbanization tends to slightly enhance the amount of summer precipitation and extremes,and reduce their magnitude of decrease in urban-influenced regions.Quantitatively, the influence of urbanization on the amount of summer precipitation and extremes during 1960–2012 was estimated to be approximately?2.09 and?1.04 mm(10 yr)?1,accounting for 10.4%and 6.0%,respectively,of the large-scale mean decreasing trends.For the frequencies of precipitation and extremes,the influence of urbanization was estimated to be 5.7%and?8.9%,respectively,suggesting that urbanization may reduce the magnitude of decrease of precipitation frequency,but enhance that of extremes in the urban areas.For precipitation intensity,urbanization tends to promote the decreasing trends in both summer precipitation and extremes, by approximately 6.8%and 51.5%,respectively.

    It is worth noting that the homogenization of daily precipitation series is a challenging topic,as daily precipitation is not a continuous variable and has larger variability.Hence, there is greater uncertainty in the homogenization of precipitation than that of other climatic variables,such as temperature.Nevertheless,the additional adjustments based on the ECDF method provide an experimental approach to improving the homogeneity of daily precipitation series in terms of probability distribution,showing potential for further application in homogenizing daily climate series.

    Acknowledgements.This study was supported by the State Programs of Science and Technology Development(Grant No. 2012BAC22B04),the National Natural Science Foundation of China(Grant No.41475078),and the National Department Public Benefit Research Foundation of China(Grant No. GYHY201206013).WU Hongyi was supported by the Science and Technology program of the Institute of Urban Meteorology(GrantNo.IUMKY201302PP0102).

    REFERENCES

    Alexandersson,H.,1986:A homogeneity test applied to precipitation data.J.Climatol.,6(6),661–675.

    Beaulieu,C.,O.Seidou,T.B.M.J.Ouarda,X.B.Zhang,G. Boulet,and A.Yagouti,2008:Intercomparison of homogenization techniques for precipitation data.Water Resour.Res.,44(2),W02425,doi:10.1029/2006WR005615.

    Beaulieu,C.,O.Seidou,T.B.M.J.Ouarda,and X.B.Zhang, 2009:Intercomparison of homogenization techniques for precipitation data continued:Comparison of two recent Bayesian change point models.Water Resour.Res.,45(8),W08410, doi:10.1029/2008WR007501.

    Begert,M.,T.Schlegel,and W.Kirchhofer,2005:Homogeneous temperature and precipitation series ofSwitzerland from 1864 to 2000.Int.J.Climatol.,25(1),65–80.

    Cong,Z.T.,J.J.Zhao,D.W.Yang,and G.H.Ni,2010:Understanding the hydrological trends of river basins in China.J. Hydrol.,388(3–4),350–356.

    Deng,A.J.,S.Y.Tao,and L.T.Chen,1989:The EOF analysis of rainfall in China during monsoon season.Chinese J.Atomos. Sci.,13(3),289–295.(in Chinese)

    Dosio,A.,and P.Paruolo,2011:Bias correction of the ENSEMBLES high-resolution climate change projections for use by impact models:Evaluation on the present climate.J.Geophys.Res.,116,D16106,doi:10.1029/2011JD015934.

    Easterling,D.R.,and T.C.Peterson,1995:A new method for detecting undocumented discontinuities in climatological time series.Int.J.Climatol.,15(4),369–377.

    Gonz′alez-Rouco,J.F.,J.L.Jim′enez,V.Quesada,and F.Valero, 2001:Quality control and homogeneity of precipitation data in the southwest of Europe.J.Climate,14(5),964–978.

    Han,Z.Q.,Z.W.Yan,Z.Li,W.D.Liu,and Y.C.Wang,2014:Impact of urbanization on low-temperature precipitation in Beijing during 1960–2008.Adv.Atmos.Sci.,31(1),48–56,doi: 10.1007/s00376-013-2211-3.

    Hand,L.M.,and J.M.Shepherd,2009:An investigation of warmseason spatial rainfall variability in Oklahoma City:Possible linkages to urbanization and prevailing wind.J.Appl.Meteor. Climatol.,48(2),251–269.

    Jiang,Z.H.,Q.Huang,and Q.X.Li,2008:Study of precipitation series homogeneous adjustment and their correction over China in the last 50 Years.Climatic and Environmental Research,13(1),67–74.(in Chinese)

    Jones,P.D.,E.B.Horton,C.K.Folland,M.Hulme,D.E. Parker,and T.A.Basnett,1999:The use of indices to identify changes in climatic extremes.Climatic Change,42,131–149.

    Kishtawal,C.M.,D.Niyogi,M.Tewari,R.A.Pielke,and J.M.Shepherd,2010:Urbanization signature in the observed heavy rainfall climatology over India.Int.J.Climatol.,30(13),1908–1916.

    Li.J.,R.C.Yu,and J.J.Wang,2008a:Diurnal variations of summer precipitation in Beijing.Chinese Science Bulletin,53(12),1933–1936.

    Li,Q.X.,Z.H.Jiang,Q.Huang,and Y.You,2008b:The experimental detecting and adjusting of the precipitation data homogeneity in the Yangtze Delta.Journal of Applied Meteorological Science,19(2),219–226.(in Chinese)

    Li,Q.X.,J.D.Peng,and Y.Shen,2012:Development of homogenized monthly precipitation dataset in China during 1900–2009.Acta Geographica Sinica,67(3),301–311.(in Chinese)

    Li,W.B.,S.Chen,G.X.Chen,W.M.Sha,C.Luo,Y.R.Feng, Z.P.Wen,and B.M.Wang,2011b:Urbanization signatures in strong versus weak precipitation over the Pearl River Delta metropolitan regions of China.Environmental Research Letters,6(3),034020.

    Li,Z.,and Z.W.Yan,2009:Homogenized daily mean/maximum/ minimum temperature series for China from 1960-2008.Atmos.Oceanic Sci.Lett.,2(4),237–243.

    Li,Z.,and Z.W.Yan,2010:Application of multiple analysis of series for homogenization to Beijing daily temperature series(1960-2006).Adv.Atmos.Sci.,27(4),777–787,doi: 10.1007/s00376-009-9052-0.

    Li,Z.,Z.W.Yan,K.Tu,W.D.Liu,and Y.C.Wang,2011a: Changes in wind speed and extremes in Beijing during 1960-2008 based on homogenized observations.Adv.Atmos.Sci.,28(2),408–420,doi:10.1007/s00376-010-0018-z.

    Li,Z.,Z.W.Yan,L.J.Cao,and P.D.Jones,2014:Adjusting inhomogeneous daily temperature variability using wavelet analysis.Int.J.Climatol.,34(4),1196–1207.

    Liu,X.N.,and J.A.Sun,1995:An inhomogeneity teststudy ofannual precipitation series.Meteorological Monthly,21(8),3–6. (in Chinese)

    Lund,R.,and J.Reeves,2002:Detection of undocumented changepoints:A revision of the two-phase regression model. J.Climate,15(17),2547–2554.

    Manton,M.J.,and Coauthors,2001:Trends in extreme daily rainfall and temperature in Southeast Asia and the South Pacific: 1961–1998.Int.J.Climatol.,21(3),269–284.

    McRoberts,D.B.,and J.W.Nielsen-Gammon,2011:A new homogenized climate division precipitation dataset for analysis ofclimate variability and climate change.J.Appl.Meteor.Climatol.,50(6),1187–1199.

    Miao,S.G.,F.Chen,Q.C.Li,and S.Y.Fan,2011:Impacts of urban processes and urbanization on summer precipitation:A case study of heavy rainfall in Beijing on 1 August 2006.J. Appl.Meteor.Climatol.,50(4),806–825.

    Piani,C.,J.O.Haerter,and E.Coppola,2010:Statistical bias correction for daily precipitation in regional climate models over Europe.Theor.Appl.Climatol.,99(1-2),187–192.

    Potter,K.W.,1981:Illustration of a new test for detecting a shift in mean in precipitation series.Mon.Wea.Rev.,109(9),2040–2045.

    Rhoades,D.A.,and M.J.Salinger,1993:Adjustment of temperature and rainfall records for site changes.Int.J.Climatol.,13(8),899–913.

    Song,X.M.,and Coauthors,2014:Rapid urbanization and changes in spatiotemporal characteristics of precipitation in Beijing metropolitan area.J.Geophys.Res.,119(19),11 250–11 271.

    Sun,J.S.,and W.J.Shu,2007:The effect of urban heat island on winter and summer precipitation in Beijing region.Chinese J. Atmos Sci.,31(2),311–320.(in Chinese)

    Szentimrey,T.,1999:Multiple Analysis of Series for Homogenization(MASH).Proc.Second Seminar for Homogenization of Surface Climatological Data,Budapest,Hungary,WMO, WCDMP-No.41,27–46.

    Tao,S.Y.,C.B.Fu,Z.M.Zeng,and Q.Zhang,1991:Two climate data sets of the long-term instrumental records in China.ORNL Technical Report CDIAC-47,Oak Ridge National Laboratory,Tennessee,198 pp.

    Trewin,B.C.,2013:A daily homogenized temperature data set for Australia.Int.J.Climatol.,33(6),1510–1529.

    Tuomenvirta,H.,2001:Homogeneity adjustments of temperature and precipitation series–Finnish and Nordic data.Int.J.Climatol.,21(4),495–506.

    Vincent,L.A.,X.Zhang,B.R.Bonsal,and W.D.Hogg,2002: Homogenization of daily temperatures over Canada.J.Climate,15(11),1322–1334.

    Wang,L.X.,H.F.Chen,Y.H.Wu,Y.Feng,and Q.Pu,2010:New techniques for the detection and adjustment of shifts in daily precipitation data series.J.Appl.Meteor.Climatol.,49(12), 2416–2436.

    Wang,X.R.,W.G.Wang,H.Z.Liu,and H.Wang,2008:Beijing region precipitation feature and some statistics of relationship between it and SHWP.Plateau Meteorology,27(4),822–829. (in Chinese)

    Xu,M.,C.P.Chang,C.B.Fu,Y.Qi,A.Robock,D.Robinson, and H.M.Zhang,2006:Steady decline of East Asian monsoon winds,1969–2000:Evidence from direct ground measurements of wind speed.J.Geophys.Res.,111,D24111,doi: 10.1029/2006JD007337.

    Xu,W.H.,Q.X.Li,X.L.Wang,S.Yang,L.J.Cao,and Y.Feng, 2013:Homogenization of Chinese daily surface air temperatures and analysis of trends in the extreme temperature indices.J.Geophys.Res.,118(17),9708–9720.

    Yan,Z.W,and P.D.Jones,2008:Detecting inhomogeneity in daily climate series using wavelet analysis.Adv.Atmos.Sci.,25(2),157–163,doi:10.1007/s00376-008-0157-7.

    Yan,Z.W,P.D.Jones,A.Moberg,H.Bergstrom,T.D.Davies, and C.Yang,2001:Recent trends in weather and seasonal cycles:an analysis of daily data from Europe and China.J. Geophys.Res.,106(D6),5123–5138.

    Yan,Z.W.,S.Bate,R.E.Chandler,V.Isham,and H.Wheater, 2002:An analysis of daily maximum wind speed in northwestern Europe using generalized linear models.J.Climate,15(15),2073–2088.

    Yan,Z.W.,Z.Li,Q.X.Li,and P.D.Jones,2010:Effects of site change and urbanisation in the Beijing temperature series 1977–2006.Int.J.Climatol.,30(8),1226–1234.

    Yang,L.,F.Q.Tian,J.A.Smith,and H.P.Hu,2014:Urban signatures in the spatial clustering of summer heavy rainfall events over the Beijing metropolitan region.J.Geophys.Res.,119(3),1203–1217.

    Yang,P.,G.Y.Ren,W.Hou,and W.D.Liu,2013:Spatial and diurnal characteristics of summer rainfall over Beijing Municipality based on a high-density AWS dataset.Int.J.Climatol.,33(13),2769–2780.

    Yang,W.,J.Andr′easson,L.P.Graham,J.Olsson,J.Rosberg,and F.Wetterhall,2010:Distribution-based scaling to improve usability of regional climate model projections for hydrological climate change impacts studies.Hydrology Research,41(3-4), 211–229.

    Yin,S.Q.,W.J.Li,D.L.Chen,J.H.Jeong,and W.L.Guo,2011: Diurnal variations of summer precipitation in the Beijing area and the possible effect of topography and urbanization.Adv. Atmos.Sci.,28(4),725–734,doi:10.1007/s00376-010-9240-y.

    You,H.L.,W.D.Liu,and G.Y.Ren,2014:Variation characteristics of precipitation extremes in Beijing during 1981-2010. Climatic and Environmental Research,19(1),69–77.(in Chinese)

    Zhang,C.L.,F.Chen,S.G.Miao,Q.C.Li,X.A.Xia,and C.Y.Xuan,2009a:Impacts of urban expansion and future green planting on summer precipitation in the Beijing metropolitan area.J.Geophys.Res.,114,D02116,doi: 10.1029/2008JD010328.

    Zhang,L.J.,T.J.Hu,F.Hu,L.Li,and C.Y.Xuan,2009b:Comparison of summer precipitation evolution over the last 30 years between urban and rural region of Beijing.Climatic and Environmental Research,14(1),63–68.(in Chinese)

    Zhang,S.Y.,2006:Beijing Modernization Report2006:Evaluation Research on Urbanization in Beijing.Beijing Science and Technology Press,1–8.[Available online at http://www.bjkw. gov.cn/n8785584/n8904761/n8904900/n8904975/n8934449/ 9028425.html.](in Chinese)

    Zhang,Y.,J.A.Smith,L.F.Luo,Z.F.Wang,and M.L.Baeck, 2014:Urbanization and rainfall variability in the Beijing metropolitan region.Journal of Hydrometeorology,15(6), 2219–2235.

    :Li,Z.,Z.W.Yan,K.Tu,and H.Y.Wu,2015:Changes of precipitation and extremes and the possible effect of urbanization in the Beijing metropolitan region during 1960–2012 based on homogenized observations.Adv.Atmos.Sci.,32(9),1173–1185,

    10.1007/s00376-015-4257-x.

    15 November 2014;revised 16 January 2015,accepted 22 January 2015)

    ?Corresponding author:TU Kai

    Email:tukai@tea.ac.cn

    九九热线精品视视频播放| 亚州av有码| 日日啪夜夜撸| 午夜福利欧美成人| 亚洲av成人av| 中文字幕熟女人妻在线| 毛片一级片免费看久久久久 | 国产午夜精品久久久久久一区二区三区 | 夜夜看夜夜爽夜夜摸| 国产老妇女一区| 亚洲成人免费电影在线观看| 大型黄色视频在线免费观看| 少妇丰满av| 国内精品一区二区在线观看| 午夜精品一区二区三区免费看| 禁无遮挡网站| 成年免费大片在线观看| 久久国产精品人妻蜜桃| netflix在线观看网站| 成人欧美大片| 成人永久免费在线观看视频| 亚洲国产高清在线一区二区三| 成人毛片a级毛片在线播放| 日韩精品青青久久久久久| 91av网一区二区| 国产高清激情床上av| 男人舔奶头视频| 精品人妻视频免费看| 免费av毛片视频| 欧美精品国产亚洲| 日本a在线网址| 亚洲va在线va天堂va国产| 精品人妻熟女av久视频| 亚洲国产精品久久男人天堂| 91狼人影院| av福利片在线观看| 日日啪夜夜撸| 精品乱码久久久久久99久播| 欧美性猛交╳xxx乱大交人| 欧美xxxx黑人xx丫x性爽| 国产淫片久久久久久久久| 91精品国产九色| 午夜激情福利司机影院| 久久国产精品人妻蜜桃| 国内精品一区二区在线观看| 国产三级在线视频| 有码 亚洲区| 亚洲aⅴ乱码一区二区在线播放| 最近在线观看免费完整版| 一级av片app| 99久久成人亚洲精品观看| 天堂av国产一区二区熟女人妻| 黄色一级大片看看| 精品国内亚洲2022精品成人| 久久九九热精品免费| 美女 人体艺术 gogo| 99热精品在线国产| 一区二区三区高清视频在线| 国产黄a三级三级三级人| 男人和女人高潮做爰伦理| 十八禁国产超污无遮挡网站| 老熟妇仑乱视频hdxx| 国产国拍精品亚洲av在线观看| 成年人黄色毛片网站| 国产色爽女视频免费观看| 国产亚洲精品久久久久久毛片| 精品一区二区免费观看| 午夜免费激情av| 午夜久久久久精精品| 国产精品亚洲一级av第二区| 国产国拍精品亚洲av在线观看| av女优亚洲男人天堂| 99热精品在线国产| 亚洲自拍偷在线| 少妇裸体淫交视频免费看高清| 最新在线观看一区二区三区| 老熟妇乱子伦视频在线观看| 日日撸夜夜添| 久99久视频精品免费| 国产精品电影一区二区三区| 亚洲真实伦在线观看| 极品教师在线免费播放| 日本三级黄在线观看| 亚洲专区国产一区二区| 亚洲熟妇中文字幕五十中出| 一级黄片播放器| 国产成人一区二区在线| 婷婷六月久久综合丁香| 午夜精品久久久久久毛片777| 99精品在免费线老司机午夜| 久久久午夜欧美精品| 亚洲性夜色夜夜综合| 两人在一起打扑克的视频| 国产精品自产拍在线观看55亚洲| 如何舔出高潮| 日本成人三级电影网站| 国产黄a三级三级三级人| 热99re8久久精品国产| 色综合婷婷激情| 亚洲精品一区av在线观看| 一级毛片久久久久久久久女| 乱码一卡2卡4卡精品| 亚洲第一区二区三区不卡| 精品午夜福利在线看| 麻豆成人午夜福利视频| 搡老熟女国产l中国老女人| 亚洲欧美日韩无卡精品| 午夜福利在线观看吧| 国产 一区 欧美 日韩| 亚洲经典国产精华液单| 亚洲自偷自拍三级| 午夜激情欧美在线| 国产午夜精品论理片| 91av网一区二区| 欧美成人a在线观看| 老师上课跳d突然被开到最大视频| 亚洲中文字幕日韩| 成年女人永久免费观看视频| 午夜爱爱视频在线播放| 国产午夜精品论理片| 亚洲自拍偷在线| 久久亚洲真实| 欧美潮喷喷水| 91av网一区二区| videossex国产| 又黄又爽又刺激的免费视频.| 我的老师免费观看完整版| 国产精品日韩av在线免费观看| 国产中年淑女户外野战色| 搡老熟女国产l中国老女人| 国产不卡一卡二| 午夜免费激情av| a级一级毛片免费在线观看| 一进一出抽搐gif免费好疼| 丰满乱子伦码专区| 女人被狂操c到高潮| 中国美白少妇内射xxxbb| 美女大奶头视频| 久久午夜福利片| 国产成人a区在线观看| 久久久久久九九精品二区国产| 噜噜噜噜噜久久久久久91| 亚洲精品456在线播放app | 久久婷婷人人爽人人干人人爱| 97热精品久久久久久| 久久久精品大字幕| 色在线成人网| 桃红色精品国产亚洲av| 韩国av一区二区三区四区| 久久人人精品亚洲av| 亚洲欧美日韩卡通动漫| 国产伦精品一区二区三区四那| av在线观看视频网站免费| 狠狠狠狠99中文字幕| 老师上课跳d突然被开到最大视频| 一区二区三区高清视频在线| 天天躁日日操中文字幕| 亚洲人成网站高清观看| 无遮挡黄片免费观看| 久久午夜福利片| 成人一区二区视频在线观看| 黄片wwwwww| a在线观看视频网站| 成人一区二区视频在线观看| 国内揄拍国产精品人妻在线| 日本a在线网址| 欧美极品一区二区三区四区| 身体一侧抽搐| 国产精品女同一区二区软件 | 国产精品久久久久久精品电影| 中国美女看黄片| 直男gayav资源| 国产一区二区三区视频了| 国产精品爽爽va在线观看网站| 老司机午夜福利在线观看视频| 欧美最新免费一区二区三区| 国产精品爽爽va在线观看网站| 亚洲人成网站在线播| 亚洲第一区二区三区不卡| 乱人视频在线观看| 真实男女啪啪啪动态图| 国产精品无大码| 中国美女看黄片| 国产在线男女| 日日撸夜夜添| 极品教师在线视频| 国产爱豆传媒在线观看| 久久久久国内视频| 欧美人与善性xxx| 黄色女人牲交| 色视频www国产| .国产精品久久| 变态另类丝袜制服| 日本在线视频免费播放| 精品人妻1区二区| 天天躁日日操中文字幕| 搡老熟女国产l中国老女人| 一进一出好大好爽视频| 国产69精品久久久久777片| 九九在线视频观看精品| 久久久久久久午夜电影| 中文在线观看免费www的网站| 中文字幕av成人在线电影| 校园人妻丝袜中文字幕| 中文字幕久久专区| 婷婷六月久久综合丁香| 夜夜夜夜夜久久久久| 在线观看66精品国产| 国产精品爽爽va在线观看网站| 女的被弄到高潮叫床怎么办 | 国产美女午夜福利| 亚洲经典国产精华液单| 色综合站精品国产| 亚洲国产精品久久男人天堂| 日韩一本色道免费dvd| 精品免费久久久久久久清纯| 国产探花在线观看一区二区| 久久久成人免费电影| 非洲黑人性xxxx精品又粗又长| 99久久九九国产精品国产免费| 国产成年人精品一区二区| 狂野欧美白嫩少妇大欣赏| 日韩精品中文字幕看吧| 国产精品三级大全| 精品久久久久久成人av| 此物有八面人人有两片| 一本久久中文字幕| 老司机深夜福利视频在线观看| 精品乱码久久久久久99久播| 极品教师在线视频| 亚洲精品亚洲一区二区| 国产成人影院久久av| 91在线观看av| 久久精品国产清高在天天线| 成年女人毛片免费观看观看9| av.在线天堂| 精品久久久久久久久av| 一卡2卡三卡四卡精品乱码亚洲| 免费观看在线日韩| 男女之事视频高清在线观看| 国产精品乱码一区二三区的特点| 床上黄色一级片| 精品99又大又爽又粗少妇毛片 | 午夜激情欧美在线| 永久网站在线| 国产高清有码在线观看视频| 久久国产精品人妻蜜桃| 日韩精品有码人妻一区| 亚洲欧美激情综合另类| av.在线天堂| 亚洲七黄色美女视频| 大又大粗又爽又黄少妇毛片口| 97人妻精品一区二区三区麻豆| 狠狠狠狠99中文字幕| 黄色欧美视频在线观看| 国产免费av片在线观看野外av| 日本黄色片子视频| 91午夜精品亚洲一区二区三区 | 中文字幕av成人在线电影| 最后的刺客免费高清国语| 老司机深夜福利视频在线观看| 国产成人福利小说| 国内精品久久久久精免费| 搡老岳熟女国产| 午夜爱爱视频在线播放| 国产熟女欧美一区二区| 欧美在线一区亚洲| 国产成人一区二区在线| 国产一区二区三区在线臀色熟女| 亚洲在线观看片| 男女下面进入的视频免费午夜| 欧美色欧美亚洲另类二区| 免费观看在线日韩| 春色校园在线视频观看| 亚洲精品一卡2卡三卡4卡5卡| 免费看美女性在线毛片视频| 色综合婷婷激情| 精品午夜福利在线看| 欧洲精品卡2卡3卡4卡5卡区| 日本色播在线视频| 狂野欧美白嫩少妇大欣赏| 免费观看在线日韩| 春色校园在线视频观看| 九九在线视频观看精品| 欧美成人一区二区免费高清观看| 婷婷精品国产亚洲av| 99九九线精品视频在线观看视频| 看片在线看免费视频| 国产91精品成人一区二区三区| 欧美成人免费av一区二区三区| 亚洲精品在线观看二区| 久99久视频精品免费| 欧美三级亚洲精品| 免费搜索国产男女视频| 亚洲最大成人中文| 深爱激情五月婷婷| 亚洲成av人片在线播放无| 少妇的逼好多水| 尤物成人国产欧美一区二区三区| 一区福利在线观看| 如何舔出高潮| 日韩中字成人| 老司机午夜福利在线观看视频| 欧美又色又爽又黄视频| 亚洲人成伊人成综合网2020| 波多野结衣高清无吗| 久久国内精品自在自线图片| 欧美三级亚洲精品| 给我免费播放毛片高清在线观看| 熟女电影av网| 久久久久国内视频| 国产精品亚洲一级av第二区| 性色avwww在线观看| 看十八女毛片水多多多| 久久6这里有精品| 中国美白少妇内射xxxbb| 欧美日韩中文字幕国产精品一区二区三区| 中文字幕久久专区| 亚洲自偷自拍三级| 免费观看精品视频网站| 国产黄a三级三级三级人| 亚洲av成人精品一区久久| 久久精品夜夜夜夜夜久久蜜豆| 蜜桃亚洲精品一区二区三区| 婷婷色综合大香蕉| 校园春色视频在线观看| 国产一区二区在线观看日韩| 国产单亲对白刺激| 亚洲美女黄片视频| 最新在线观看一区二区三区| 国产一区二区三区av在线 | 成人性生交大片免费视频hd| 久久久久久九九精品二区国产| 在现免费观看毛片| 亚洲一区二区三区色噜噜| 桃红色精品国产亚洲av| 女同久久另类99精品国产91| 亚洲欧美清纯卡通| 三级国产精品欧美在线观看| 日韩亚洲欧美综合| 人妻夜夜爽99麻豆av| 给我免费播放毛片高清在线观看| 国产精品一区二区性色av| 毛片女人毛片| 日韩欧美 国产精品| 禁无遮挡网站| 成人av一区二区三区在线看| 亚洲国产色片| 久久热精品热| 国产日本99.免费观看| 久久久午夜欧美精品| 精品无人区乱码1区二区| 男人舔奶头视频| 中文字幕精品亚洲无线码一区| 大型黄色视频在线免费观看| 一进一出抽搐动态| 12—13女人毛片做爰片一| 最近在线观看免费完整版| av在线蜜桃| 精品一区二区三区视频在线| 成人毛片a级毛片在线播放| 久久午夜亚洲精品久久| 午夜福利高清视频| 日本撒尿小便嘘嘘汇集6| av中文乱码字幕在线| 亚洲成a人片在线一区二区| 国产三级中文精品| 国产69精品久久久久777片| 亚洲美女黄片视频| 成人亚洲精品av一区二区| 在线观看一区二区三区| 亚洲不卡免费看| 成人国产麻豆网| 精品99又大又爽又粗少妇毛片 | 窝窝影院91人妻| 精品久久国产蜜桃| 欧美日本视频| 欧美日韩乱码在线| 亚洲色图av天堂| 亚洲五月天丁香| 国语自产精品视频在线第100页| 女人十人毛片免费观看3o分钟| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av一区综合| 日本熟妇午夜| 色播亚洲综合网| 国产黄片美女视频| 午夜老司机福利剧场| 在线a可以看的网站| 久久天躁狠狠躁夜夜2o2o| 欧美激情久久久久久爽电影| 色播亚洲综合网| 在线观看66精品国产| 久久久久九九精品影院| 三级国产精品欧美在线观看| 少妇人妻精品综合一区二区 | 91狼人影院| 一级毛片久久久久久久久女| 伦理电影大哥的女人| 91久久精品国产一区二区成人| 婷婷精品国产亚洲av| 午夜老司机福利剧场| 亚洲电影在线观看av| 久久久国产成人免费| 国产精品综合久久久久久久免费| 在线免费观看不下载黄p国产 | 天堂影院成人在线观看| 日本精品一区二区三区蜜桃| av专区在线播放| 国产精品日韩av在线免费观看| 18+在线观看网站| 在线观看美女被高潮喷水网站| 成人高潮视频无遮挡免费网站| 亚洲欧美日韩卡通动漫| 99视频精品全部免费 在线| 亚洲四区av| 国产女主播在线喷水免费视频网站 | 三级毛片av免费| 欧美又色又爽又黄视频| 女的被弄到高潮叫床怎么办 | a级一级毛片免费在线观看| 舔av片在线| 天天一区二区日本电影三级| 天天躁日日操中文字幕| 久久久久免费精品人妻一区二区| 一区二区三区四区激情视频 | 黄色配什么色好看| 嫩草影院入口| 免费在线观看影片大全网站| 国产美女午夜福利| 亚洲精品日韩av片在线观看| 又爽又黄a免费视频| 亚洲熟妇中文字幕五十中出| 亚洲国产日韩欧美精品在线观看| 男女之事视频高清在线观看| 欧美日韩综合久久久久久 | 日本熟妇午夜| 国产精品女同一区二区软件 | 噜噜噜噜噜久久久久久91| 成人欧美大片| 亚洲av.av天堂| 亚洲国产精品久久男人天堂| 亚洲精品成人久久久久久| 国产视频内射| 男人舔奶头视频| 少妇人妻精品综合一区二区 | 少妇猛男粗大的猛烈进出视频 | 露出奶头的视频| 人妻丰满熟妇av一区二区三区| 亚洲av美国av| 国内久久婷婷六月综合欲色啪| 精华霜和精华液先用哪个| 日本黄色片子视频| 无遮挡黄片免费观看| 两人在一起打扑克的视频| 亚洲精品亚洲一区二区| av女优亚洲男人天堂| 看十八女毛片水多多多| 国产黄色小视频在线观看| 看片在线看免费视频| 国产av在哪里看| 国产精品美女特级片免费视频播放器| 不卡一级毛片| 最后的刺客免费高清国语| 亚洲精品一卡2卡三卡4卡5卡| 久久6这里有精品| 欧美精品啪啪一区二区三区| 少妇裸体淫交视频免费看高清| 在线国产一区二区在线| 久久久久久久久大av| 超碰av人人做人人爽久久| 欧美最黄视频在线播放免费| 国产欧美日韩精品亚洲av| 99热这里只有是精品50| 久久这里只有精品中国| 老司机午夜福利在线观看视频| 悠悠久久av| 欧美色视频一区免费| 国产高潮美女av| 国产真实乱freesex| 男女啪啪激烈高潮av片| 三级毛片av免费| 村上凉子中文字幕在线| 欧美日韩亚洲国产一区二区在线观看| 波多野结衣高清无吗| 亚洲人成伊人成综合网2020| 欧美又色又爽又黄视频| 亚洲欧美激情综合另类| 97人妻精品一区二区三区麻豆| 欧美另类亚洲清纯唯美| 成人国产麻豆网| 欧美国产日韩亚洲一区| 国产伦精品一区二区三区视频9| 别揉我奶头 嗯啊视频| 九九久久精品国产亚洲av麻豆| 老司机午夜福利在线观看视频| 国产午夜精品久久久久久一区二区三区 | 99久久久亚洲精品蜜臀av| 毛片女人毛片| 成人av一区二区三区在线看| 好男人在线观看高清免费视频| 国产熟女欧美一区二区| 国产成人影院久久av| av在线亚洲专区| 1000部很黄的大片| 麻豆成人午夜福利视频| 国产亚洲精品久久久久久毛片| 亚洲欧美日韩卡通动漫| 看十八女毛片水多多多| 99久久无色码亚洲精品果冻| 白带黄色成豆腐渣| 中亚洲国语对白在线视频| 夜夜看夜夜爽夜夜摸| 午夜爱爱视频在线播放| av中文乱码字幕在线| 麻豆av噜噜一区二区三区| 午夜福利在线在线| 欧美色欧美亚洲另类二区| 亚洲av中文av极速乱 | 一个人看视频在线观看www免费| 日本爱情动作片www.在线观看 | 亚洲电影在线观看av| 成人性生交大片免费视频hd| 亚洲 国产 在线| 国产单亲对白刺激| 久久精品久久久久久噜噜老黄 | 神马国产精品三级电影在线观看| 在线免费观看不下载黄p国产 | 蜜桃亚洲精品一区二区三区| 18禁黄网站禁片免费观看直播| 精品日产1卡2卡| 国产一区二区在线观看日韩| 99热这里只有精品一区| 国产免费av片在线观看野外av| 亚洲在线自拍视频| 欧美一区二区国产精品久久精品| 看黄色毛片网站| 国产精品日韩av在线免费观看| 成人欧美大片| 国产伦人伦偷精品视频| 白带黄色成豆腐渣| 欧美成人一区二区免费高清观看| 日本黄色视频三级网站网址| 午夜爱爱视频在线播放| 国产爱豆传媒在线观看| 18禁黄网站禁片午夜丰满| 2021天堂中文幕一二区在线观| 亚洲avbb在线观看| 亚洲精品色激情综合| 久久久久久久久中文| 搡老妇女老女人老熟妇| av天堂中文字幕网| 久久国内精品自在自线图片| 97碰自拍视频| 成熟少妇高潮喷水视频| 亚洲一级一片aⅴ在线观看| 免费观看精品视频网站| 欧美性猛交╳xxx乱大交人| 少妇高潮的动态图| 中国美女看黄片| av在线蜜桃| 韩国av在线不卡| 欧美bdsm另类| 国产在视频线在精品| 亚洲av免费在线观看| 中文字幕久久专区| 偷拍熟女少妇极品色| 在线a可以看的网站| 高清毛片免费观看视频网站| 久久热精品热| 性欧美人与动物交配| 黄片wwwwww| 国产综合懂色| 综合色av麻豆| 99久久九九国产精品国产免费| 午夜免费激情av| 成人国产一区最新在线观看| 成人鲁丝片一二三区免费| 老司机福利观看| 在线观看一区二区三区| 国产aⅴ精品一区二区三区波| 丰满乱子伦码专区| 老熟妇乱子伦视频在线观看| 中亚洲国语对白在线视频| 1000部很黄的大片| 国产亚洲欧美98| 久久久成人免费电影| 国产成人a区在线观看| 欧美性猛交黑人性爽| 日韩一区二区视频免费看| 极品教师在线免费播放| 白带黄色成豆腐渣| 在现免费观看毛片| 嫩草影视91久久| 欧美中文日本在线观看视频| 久久精品久久久久久噜噜老黄 | 国产激情偷乱视频一区二区| 在线观看美女被高潮喷水网站| 国产视频一区二区在线看| 精品久久久久久久末码| 久久久久久九九精品二区国产| 最新在线观看一区二区三区| 九九久久精品国产亚洲av麻豆| 国产色婷婷99| 国内毛片毛片毛片毛片毛片| 久久久久国产精品人妻aⅴ院| 日韩欧美精品免费久久| 一区二区三区高清视频在线| 欧美3d第一页| 99久久精品热视频| 久久久精品大字幕| 91久久精品国产一区二区三区| 日日摸夜夜添夜夜添av毛片 | 成人欧美大片| 亚洲性夜色夜夜综合|