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

    黑人欧美特级aaaaaa片| 亚洲国产精品sss在线观看| 欧美中文综合在线视频| 久久精品成人免费网站| 午夜a级毛片| 久久国产精品影院| 国产精品野战在线观看| 美国免费a级毛片| 男男h啪啪无遮挡| 在线av久久热| 亚洲九九香蕉| 两个人免费观看高清视频| 亚洲熟女毛片儿| 国产精品国产高清国产av| 国产亚洲av高清不卡| 悠悠久久av| 久久人人精品亚洲av| 亚洲精品久久成人aⅴ小说| 黄网站色视频无遮挡免费观看| 国产又色又爽无遮挡免费看| 国产精品一区二区三区四区久久 | 久久久久国内视频| 国产亚洲欧美98| 色av中文字幕| 一级毛片女人18水好多| 一进一出抽搐gif免费好疼| 国产精品久久久久久亚洲av鲁大| 免费av毛片视频| 国产成人精品久久二区二区91| 成人手机av| 在线看三级毛片| 国产av一区在线观看免费| 欧美精品啪啪一区二区三区| 国产精品影院久久| 一进一出好大好爽视频| 免费人成视频x8x8入口观看| 99热6这里只有精品| 美女扒开内裤让男人捅视频| 国产亚洲精品综合一区在线观看 | 12—13女人毛片做爰片一| 亚洲精品在线观看二区| 美女高潮喷水抽搐中文字幕| 听说在线观看完整版免费高清| 一本一本综合久久| 少妇的丰满在线观看| 国产成人av激情在线播放| 香蕉久久夜色| 欧美在线黄色| 久久香蕉国产精品| a在线观看视频网站| 欧洲精品卡2卡3卡4卡5卡区| 看黄色毛片网站| 少妇 在线观看| 久久精品人妻少妇| 国产一卡二卡三卡精品| 满18在线观看网站| 18禁黄网站禁片午夜丰满| 极品教师在线免费播放| 亚洲精品国产精品久久久不卡| 国产一区二区三区在线臀色熟女| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩有码中文字幕| 久久精品国产综合久久久| 亚洲三区欧美一区| 一区福利在线观看| 俺也久久电影网| e午夜精品久久久久久久| 日韩欧美在线二视频| 亚洲午夜精品一区,二区,三区| 亚洲成人久久性| 欧美在线黄色| 18美女黄网站色大片免费观看| 国产真实乱freesex| 久久中文字幕一级| 成人18禁在线播放| 久9热在线精品视频| 久久亚洲真实| 欧美性猛交黑人性爽| 天天躁狠狠躁夜夜躁狠狠躁| 国产视频一区二区在线看| 老司机深夜福利视频在线观看| 国产亚洲精品久久久久久毛片| 啦啦啦免费观看视频1| 十八禁网站免费在线| 国产色视频综合| 法律面前人人平等表现在哪些方面| 日日夜夜操网爽| 亚洲最大成人中文| 精品乱码久久久久久99久播| av视频在线观看入口| 精品熟女少妇八av免费久了| 国产精品 欧美亚洲| 国产亚洲av高清不卡| а√天堂www在线а√下载| 香蕉av资源在线| 91成人精品电影| 精品电影一区二区在线| 一区二区三区高清视频在线| 夜夜看夜夜爽夜夜摸| 无限看片的www在线观看| 美国免费a级毛片| 不卡av一区二区三区| 欧美又色又爽又黄视频| 久久国产精品男人的天堂亚洲| 午夜福利18| 中文字幕人妻丝袜一区二区| a级毛片在线看网站| 国产日本99.免费观看| 在线观看日韩欧美| 久久这里只有精品19| 精品一区二区三区四区五区乱码| 国产1区2区3区精品| 无限看片的www在线观看| 成人国产一区最新在线观看| 丰满的人妻完整版| 美女 人体艺术 gogo| 国产人伦9x9x在线观看| 精品福利观看| 国产精品免费一区二区三区在线| 国产高清有码在线观看视频 | 国产av在哪里看| 男女床上黄色一级片免费看| 国产黄色小视频在线观看| 亚洲欧美精品综合久久99| 成人精品一区二区免费| 18禁观看日本| 黄频高清免费视频| 国产v大片淫在线免费观看| 人人澡人人妻人| 国产午夜福利久久久久久| 久久精品国产亚洲av香蕉五月| or卡值多少钱| 给我免费播放毛片高清在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲最大成人中文| 日本黄色视频三级网站网址| 婷婷六月久久综合丁香| 欧美在线一区亚洲| 亚洲一卡2卡3卡4卡5卡精品中文| 91九色精品人成在线观看| 悠悠久久av| 精品国产美女av久久久久小说| 午夜福利一区二区在线看| 国产成人精品久久二区二区91| xxxwww97欧美| 男男h啪啪无遮挡| 亚洲一区二区三区色噜噜| 亚洲国产精品sss在线观看| 久久婷婷成人综合色麻豆| 1024香蕉在线观看| 大香蕉久久成人网| 50天的宝宝边吃奶边哭怎么回事| 天天躁夜夜躁狠狠躁躁| 亚洲av日韩精品久久久久久密| 欧美色欧美亚洲另类二区| 麻豆一二三区av精品| 亚洲精品在线美女| 国产精品自产拍在线观看55亚洲| 成年女人毛片免费观看观看9| 一本综合久久免费| 亚洲国产欧洲综合997久久, | 欧美另类亚洲清纯唯美| 亚洲三区欧美一区| 亚洲国产欧洲综合997久久, | 成人特级黄色片久久久久久久| 欧美av亚洲av综合av国产av| 美女高潮到喷水免费观看| 亚洲成国产人片在线观看| 亚洲国产日韩欧美精品在线观看 | 香蕉丝袜av| 中文字幕人妻熟女乱码| 国产av一区在线观看免费| 亚洲欧美精品综合一区二区三区| 久久香蕉国产精品| 国产蜜桃级精品一区二区三区| 色综合站精品国产| 中文字幕另类日韩欧美亚洲嫩草| 日韩欧美免费精品| 免费看美女性在线毛片视频| 俄罗斯特黄特色一大片| 亚洲人成77777在线视频| 人成视频在线观看免费观看| 午夜福利成人在线免费观看| 中文字幕精品免费在线观看视频| 午夜福利一区二区在线看| 欧美zozozo另类| 欧美日韩一级在线毛片| 白带黄色成豆腐渣| 国产免费av片在线观看野外av| 大型av网站在线播放| 国产成人av教育| 日本熟妇午夜| 香蕉av资源在线| 国产精品美女特级片免费视频播放器 | 可以免费在线观看a视频的电影网站| 成人精品一区二区免费| 在线观看66精品国产| 亚洲,欧美精品.| 1024视频免费在线观看| 91麻豆av在线| 亚洲一区中文字幕在线| 欧美激情 高清一区二区三区| 长腿黑丝高跟| 两人在一起打扑克的视频| 亚洲国产精品999在线| 国产亚洲欧美在线一区二区| 午夜老司机福利片| 国产成人av激情在线播放| 女人高潮潮喷娇喘18禁视频| 天天添夜夜摸| 精品欧美一区二区三区在线| 久久久久久人人人人人| 亚洲熟妇熟女久久| 久久久久久亚洲精品国产蜜桃av| 日本成人三级电影网站| 18禁国产床啪视频网站| 大香蕉久久成人网| 亚洲av片天天在线观看| 国产在线观看jvid| 亚洲精品久久国产高清桃花| 国产精品电影一区二区三区| 99国产极品粉嫩在线观看| 国产伦人伦偷精品视频| 国产激情久久老熟女| 又黄又粗又硬又大视频| 老熟妇仑乱视频hdxx| 日日干狠狠操夜夜爽| 亚洲欧美精品综合久久99| 欧美性长视频在线观看| 亚洲中文av在线| 日本一区二区免费在线视频| 亚洲av中文字字幕乱码综合 | 久久久久九九精品影院| 国产人伦9x9x在线观看| 国产三级在线视频| 人成视频在线观看免费观看| 中文在线观看免费www的网站 | av在线播放免费不卡| 51午夜福利影视在线观看| 可以免费在线观看a视频的电影网站| 侵犯人妻中文字幕一二三四区| 国产精品免费一区二区三区在线| 久久午夜亚洲精品久久| 欧美精品亚洲一区二区| 亚洲国产欧洲综合997久久, | 日本免费一区二区三区高清不卡| 久久久久久国产a免费观看| 熟女少妇亚洲综合色aaa.| 少妇 在线观看| 久久精品国产亚洲av香蕉五月| 精品一区二区三区四区五区乱码| 老熟妇乱子伦视频在线观看| av在线播放免费不卡| 午夜老司机福利片| 手机成人av网站| 久久久久久久午夜电影| 老司机福利观看| 国产精品一区二区精品视频观看| 成人亚洲精品一区在线观看| 久久久国产成人精品二区| 美国免费a级毛片| 亚洲精品美女久久久久99蜜臀| 成人手机av| 中文字幕久久专区| 欧美黑人精品巨大| 后天国语完整版免费观看| 亚洲人成网站高清观看| www日本黄色视频网| 免费搜索国产男女视频| 欧美一级a爱片免费观看看 | 天堂影院成人在线观看| 欧美日本亚洲视频在线播放| 老鸭窝网址在线观看| 一本精品99久久精品77| 香蕉久久夜色| 午夜免费观看网址| 无人区码免费观看不卡| 黑人欧美特级aaaaaa片| 人人妻,人人澡人人爽秒播| 丝袜美腿诱惑在线| 美女扒开内裤让男人捅视频| 国产精品影院久久| 1024视频免费在线观看| 日韩国内少妇激情av| 18禁黄网站禁片免费观看直播| 精品久久久久久久末码| 国产亚洲精品综合一区在线观看 | 少妇粗大呻吟视频| 99国产综合亚洲精品| 天堂影院成人在线观看| av中文乱码字幕在线| 日韩欧美免费精品| e午夜精品久久久久久久| 在线十欧美十亚洲十日本专区| 一区二区日韩欧美中文字幕| 午夜免费激情av| 曰老女人黄片| 国产精品日韩av在线免费观看| 国产成人av教育| 不卡av一区二区三区| 他把我摸到了高潮在线观看| 久久香蕉激情| 免费观看精品视频网站| 国产成人系列免费观看| 18禁黄网站禁片午夜丰满| 美国免费a级毛片| 国产av一区二区精品久久| 久久精品国产清高在天天线| 无人区码免费观看不卡| 99国产精品一区二区三区| 亚洲片人在线观看| 成人国语在线视频| 悠悠久久av| 18美女黄网站色大片免费观看| 中文字幕av电影在线播放| 欧美又色又爽又黄视频| 亚洲色图 男人天堂 中文字幕| 国产亚洲精品av在线| 亚洲国产毛片av蜜桃av| 久久这里只有精品19| 亚洲色图av天堂| 十八禁网站免费在线| 女警被强在线播放| 不卡一级毛片| 亚洲aⅴ乱码一区二区在线播放 | 久久久久久亚洲精品国产蜜桃av| 大型黄色视频在线免费观看| 母亲3免费完整高清在线观看| e午夜精品久久久久久久| 精品一区二区三区视频在线观看免费| 亚洲九九香蕉| 午夜久久久久精精品| 国产又色又爽无遮挡免费看| 别揉我奶头~嗯~啊~动态视频| 淫妇啪啪啪对白视频| 男女午夜视频在线观看| 高清毛片免费观看视频网站| 一本大道久久a久久精品| 国产精品爽爽va在线观看网站 | 国产在线精品亚洲第一网站| 一a级毛片在线观看| 99riav亚洲国产免费| www.www免费av| 岛国在线观看网站| 51午夜福利影视在线观看| 人妻丰满熟妇av一区二区三区| 日韩中文字幕欧美一区二区| 露出奶头的视频| 欧美日韩亚洲国产一区二区在线观看| 国产精品一区二区精品视频观看| 免费观看人在逋| 亚洲人成网站高清观看| 巨乳人妻的诱惑在线观看| 国产成人一区二区三区免费视频网站| 亚洲一区中文字幕在线| 成在线人永久免费视频| 黄色丝袜av网址大全| 国产一区二区三区视频了| 免费女性裸体啪啪无遮挡网站| 亚洲avbb在线观看| 精品国产一区二区三区四区第35| 日韩一卡2卡3卡4卡2021年| 亚洲欧美日韩高清在线视频| 中文字幕人妻丝袜一区二区| 精品午夜福利视频在线观看一区| 一边摸一边做爽爽视频免费| 黄色女人牲交| x7x7x7水蜜桃| 麻豆成人av在线观看| 人成视频在线观看免费观看| netflix在线观看网站| 日韩欧美 国产精品| 欧美乱妇无乱码| 欧美绝顶高潮抽搐喷水| 深夜精品福利| 日日夜夜操网爽| 午夜福利欧美成人| 久久久久国产一级毛片高清牌| 动漫黄色视频在线观看| 88av欧美| 精品久久久久久久久久免费视频| 一区福利在线观看| 日韩有码中文字幕| 天天躁夜夜躁狠狠躁躁| 成年人黄色毛片网站| 亚洲精品中文字幕在线视频| 99国产精品一区二区蜜桃av| 亚洲 欧美一区二区三区| 国产精品一区二区三区四区久久 | 三级毛片av免费| 一区二区三区高清视频在线| 亚洲黑人精品在线| 成人三级做爰电影| 午夜a级毛片| 国产精品久久视频播放| 国产亚洲av高清不卡| 嫁个100分男人电影在线观看| 一级a爱视频在线免费观看| av超薄肉色丝袜交足视频| 女性被躁到高潮视频| 久久精品国产亚洲av高清一级| 欧美中文综合在线视频| 亚洲精华国产精华精| 精品福利观看| 国产成人影院久久av| 女同久久另类99精品国产91| 国产av又大| 老熟妇仑乱视频hdxx| 人人妻人人澡人人看| 成人国产一区最新在线观看| 女人被狂操c到高潮| 亚洲国产精品999在线| 精品日产1卡2卡| 国产不卡一卡二| 亚洲男人的天堂狠狠| 婷婷精品国产亚洲av在线| 国产97色在线日韩免费| 啦啦啦免费观看视频1| 国产国语露脸激情在线看| 嫩草影视91久久| 又黄又粗又硬又大视频| 亚洲专区字幕在线| 村上凉子中文字幕在线| 一区二区三区精品91| 国产又黄又爽又无遮挡在线| 久久香蕉精品热| 亚洲av美国av| 久久午夜亚洲精品久久| 日日摸夜夜添夜夜添小说| 国产爱豆传媒在线观看 | 国产精品免费视频内射| 亚洲激情在线av| 国产欧美日韩一区二区精品| 可以在线观看毛片的网站| 男人舔女人的私密视频| 99久久无色码亚洲精品果冻| 夜夜躁狠狠躁天天躁| 久久国产精品男人的天堂亚洲| 国产精华一区二区三区| 久热爱精品视频在线9| 亚洲第一青青草原| 中文字幕久久专区| 亚洲成av片中文字幕在线观看| 侵犯人妻中文字幕一二三四区| 人人妻人人澡人人看| 国产精品野战在线观看| 男人舔女人的私密视频| 欧美黄色片欧美黄色片| 大香蕉久久成人网| 伦理电影免费视频| 精品少妇一区二区三区视频日本电影| 久久人人精品亚洲av| 淫秽高清视频在线观看| 别揉我奶头~嗯~啊~动态视频| 国产激情偷乱视频一区二区| 欧美+亚洲+日韩+国产| 听说在线观看完整版免费高清| 亚洲成国产人片在线观看| 熟妇人妻久久中文字幕3abv| 91成人精品电影| 亚洲精品国产一区二区精华液| 精品人妻1区二区| 一本一本综合久久| 久久精品人妻少妇| 亚洲一卡2卡3卡4卡5卡精品中文| www.自偷自拍.com| 99热6这里只有精品| 国产99白浆流出| 99热6这里只有精品| 一进一出抽搐gif免费好疼| 免费在线观看成人毛片| www.自偷自拍.com| 欧美成狂野欧美在线观看| 亚洲全国av大片| 两性夫妻黄色片| 好男人电影高清在线观看| 在线观看免费午夜福利视频| 成人特级黄色片久久久久久久| 日本一区二区免费在线视频| 亚洲国产欧美日韩在线播放| 免费看十八禁软件| 亚洲国产欧美日韩在线播放| 国产真实乱freesex| 久久久国产欧美日韩av| 午夜激情福利司机影院| 欧美日韩瑟瑟在线播放| 午夜免费激情av| 亚洲五月色婷婷综合| 欧美日韩一级在线毛片| 12—13女人毛片做爰片一| www.www免费av| 18禁观看日本| 日本一区二区免费在线视频| 久久国产精品男人的天堂亚洲| 欧美黄色淫秽网站| 国产免费av片在线观看野外av| 日本三级黄在线观看| 亚洲精品一区av在线观看| 亚洲一区高清亚洲精品| 日韩精品青青久久久久久| 一级作爱视频免费观看| 亚洲专区中文字幕在线| 精品午夜福利视频在线观看一区| 深夜精品福利| 亚洲精品av麻豆狂野| 久久精品影院6| 真人做人爱边吃奶动态| а√天堂www在线а√下载| 日本黄色视频三级网站网址| 国产亚洲欧美98| 一进一出好大好爽视频| 中文字幕高清在线视频| 夜夜爽天天搞| 91字幕亚洲| 久久久久国产一级毛片高清牌| 香蕉av资源在线| 色播在线永久视频| 母亲3免费完整高清在线观看| 99久久久亚洲精品蜜臀av| 成人国语在线视频| 亚洲男人天堂网一区| 午夜日韩欧美国产| 久久精品成人免费网站| 身体一侧抽搐| 免费在线观看黄色视频的| 亚洲成人免费电影在线观看| 日韩av在线大香蕉| www日本在线高清视频| 亚洲av电影在线进入| 白带黄色成豆腐渣| 中亚洲国语对白在线视频| 欧美亚洲日本最大视频资源| 午夜精品久久久久久毛片777| 午夜激情av网站| 别揉我奶头~嗯~啊~动态视频| 男人舔女人的私密视频| 老司机福利观看| 日本 av在线| 桃色一区二区三区在线观看| 免费在线观看完整版高清| 欧美黄色淫秽网站| 宅男免费午夜| 91大片在线观看| 欧美午夜高清在线| 在线国产一区二区在线| 成人免费观看视频高清| 757午夜福利合集在线观看| 18禁美女被吸乳视频| 麻豆一二三区av精品| 欧美性猛交╳xxx乱大交人| 亚洲午夜精品一区,二区,三区| 精品久久久久久久末码| 岛国在线观看网站| 亚洲成人久久爱视频| 欧美又色又爽又黄视频| 色尼玛亚洲综合影院| 91麻豆精品激情在线观看国产| 国产精品 欧美亚洲| 一个人免费在线观看的高清视频| 热re99久久国产66热| 亚洲国产精品成人综合色| 久热爱精品视频在线9| xxxwww97欧美| 国产精品久久久av美女十八| 黄片播放在线免费| 啦啦啦 在线观看视频| 99久久99久久久精品蜜桃| 18禁美女被吸乳视频| 亚洲最大成人中文| 免费在线观看黄色视频的| 欧美又色又爽又黄视频| 最近最新中文字幕大全电影3 | 夜夜爽天天搞| 熟妇人妻久久中文字幕3abv| 免费搜索国产男女视频| 在线十欧美十亚洲十日本专区| 大型av网站在线播放| 在线观看免费视频日本深夜| 丝袜美腿诱惑在线| 岛国在线观看网站| 国产精品精品国产色婷婷| 好看av亚洲va欧美ⅴa在| 69av精品久久久久久| 中文亚洲av片在线观看爽| 国内毛片毛片毛片毛片毛片| 成人精品一区二区免费| 国产黄色小视频在线观看| 国产高清视频在线播放一区| 精品一区二区三区四区五区乱码| 黄片小视频在线播放| 国产精品一区二区免费欧美| 最近最新中文字幕大全免费视频| 手机成人av网站| 亚洲欧美激情综合另类| 无遮挡黄片免费观看| 国产激情欧美一区二区| 国产精品九九99| 精品不卡国产一区二区三区| a在线观看视频网站| 久久国产乱子伦精品免费另类| 国产v大片淫在线免费观看| 国产黄a三级三级三级人| 国产伦在线观看视频一区| 日韩三级视频一区二区三区| 少妇 在线观看| 大香蕉久久成人网| 国产不卡一卡二| 丰满人妻熟妇乱又伦精品不卡| 嫩草影院精品99| 少妇 在线观看| 99精品久久久久人妻精品| 精品国产美女av久久久久小说| 午夜免费鲁丝| 久久天躁狠狠躁夜夜2o2o| 999精品在线视频|