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

    A climatic environmental performance assessment method for ecological city construction:Application to Beijing Yanqi Lake

    2015-02-09 00:58:05DUWuPenXIAOXioJunDANGBin
    Advances in Climate Change Research 2015年1期

    ,*,DU Wu-Pen,XIAO Xio-Jun, DANG Bin

    aBeijing Municipal Climate Center,Beijing 100089,China

    bBeijing Institute of Architectural Design,Beijing 100045,China

    A climatic environmental performance assessment method for ecological city construction:Application to Beijing Yanqi Lake

    FANG Xiao-Yia,CHENG Chena,*,LIU Yong-Honaga,DU Wu-Penga,XIAO Xiao-Junb, DANG Bing

    aBeijing Municipal Climate Center,Beijing 100089,China

    bBeijing Institute of Architectural Design,Beijing 100045,China

    In contrast to the input perspective for evaluating planning metrics,this research takes the climatic environmental output effects as the starting point for assessing ecological city construction.Based on approaches such as observation data analysis,meteorological model simulation,and remote sensing,a set of climatic environmental performance assessment methods is developed and established.These methods mainly focus on surface ventilation assessment and urban thermal environment assessment.With the Yanqi Lake ecological development demonstration area located in Huairou district,Beijing as an example,the assessment of the local climatic environment before and after the construction are conducted,and relevant policy suggestions for urban planning and construction are presented.The results show that after development,the ventilation capacity will decrease overall and the ventilation potential index will decrease from 0.53 to 0.44.While this is not a large reduction,and is still at a favorable level,the ventilation potential in some local areas will markedly decrease.Furthermore,the thermal environment will become poorer to some extent;the urban heat island(UHI)area and intensity will increase compared with the current situation; continuous heat islands may occur in local areas;the UHI potential index of the core area will rise from 0.0878 to 0.1217(still a favorable level). Therefore,urban surfaces should be carefully developed and arranged during planning.It is suggested that the negative impacts of large areas of urban construction on the local climatic environment in the Yanqi Lake could be mitigated by 1)strengthening the airfow by introducing fresh, cold,northwesterly air via constructed ventilation corridors,2)increasing the number of ecological cold sources,particularly for water bodies and green belts to prevent the UHI in the southern region of Yanqi Lake from becoming linked with each other,and 3)considering a pre-program before sub-domain and building planning to obtain optimum building locations.Different construction standards should be developed for different ventilation potential and UHI intensity levels.For strong heat island areas,land areas should be reserved to serve as cold sources.

    Ecological city construction;Ventilation environment;Thermal environment;Performance assessment;Yanqi Lake ecological development demonstration area

    1.Introduction

    Urban environmental construction in China began with pollution treatment and ecological system optimization,and was followed by the comprehensive construction of ecological cities.Furthermore,knowledge of the environmental effects of urban construction has increased.Although the construction of ecological cities benefts the nation and its people,there are many activities that are implemented without consideration oftheir environmental effects.Therefore,scientifc quantifcation methods for evaluating ecological city construction effects have become relevant topics for academic research,policymaking,and plan implementation(Zhang and Zhang,2013). However,current evaluation methods conducted by urban planning departments usually derive from the planning indices of the relevant ministries and commissions,and mainly evaluate the construction and implementation from the perspective of the inputs.Meanwhile,evaluating the output effects of ecological city construction is still being tentatively explored. Whether the added input obtains the corresponding environmental effects directly refects the success of ecological city construction.Furthermore,when an ecological city enters the implementation phase,the focus is on the environmental effects of the output in terms of investment.Numerous problems related to the urban climate and atmospheric environment have become increasingly apparent,such as the heat island effect and air pollution.To minimize negative impacts on the local atmospheric environment,various forms of ecological city planning and construction focuson rational,scientifc methods.As a result,it is required that effective evaluation of the local climate environmental performance be conducted after the implementation of planning schemes.Based on these evaluations,we believe that the local climate environmental effects caused by ecological city construction will be evident. This information will aid the understanding of governments, scientifc research institutions,and the general public.The outcome of this research will be a technical tool to support planners for adjusting and optimizing ecological city planning schemes and enhancing input-output effciency.

    Urban ventilation refers to the manner in which rivers, lakes,valleys,and other natural ventilation channels and artifcially established urban ventilation channels are used to guide urban airfow and improve urban air quality(Lei et al., 2011).It is of great importance to enhance the self-purifcation capacity of the urban atmospheric environment,as this relieves the urban heat island(UHI)effect and improves the comfort of urban residents;it is also currently the main environmental element in urban architectural design and urban planning.At present,urban ventilation research is still at a nascent stage. The research methods mainly include feld tests,wind tunnel tests,and numerical simulation,and the research results are mainly expressed through wind feld values and visualizations (Peng,2008;Zhang et al.,2004;Birgit et al.,2007).There have been numerous overseas and domestic studies on monomer building ventilation,small-scale group building ventilation,and urban street valley ventilation(Morhmald and Shinsuke,2008;Zhao et al.,2009;Li et al.,2012),but urban ventilation evaluation is still at a preliminary stage.Therefore, it is necessary to establish a ventilation assessment method according to local conditions,and assess the ventilation environment before and after the construction of ecological cities to evaluate the effects of urban planning on the improvement of urban air quality and relief from UHI effects.

    The UHI is a climatic phenomenon in which the temperature withinanurbanareaishigherthanthatontheoutskirts.Itaffects the material and energy fows in an urban ecological system, changes the structure and functions of the urban ecological system,gives rise to a series of ecological effects,and can even affect the health of urban residents.Therefore,for the assessmentoftheeffectsofecologicalcityconstructionsonclimate,it is very important to assess the urban thermal environment, which is dominated by the UHI effect.There have been a number of studies on the UHI using meteorological data analysis,remote sensing monitoring,and numerical simulation methods.Using 30-year air temperature data observed at fve meteorological stations(three on the outskirts and two in the urban area)in Toronto,Tanzina and William(2012)researched the UHI using the time series analysis method.Using Thematic Mapper(TM)remote sensing data,Qin et al.(2001)conducted surface temperature inversion using a mono-window algorithm andanalyzedtheUHIeffectattheIsrael-Egyptborder.Usinga monolayer urban canopy model(UCM)of the Weather Research Forecast(WRF)model coupling,Li et al.(2013) researched the UHI in Singapore.

    In contrast to the input perspective for evaluating planning metrics,this paper takes the climatic environmental output effects as a starting point for assessing ecological city construction.A set of climatic environmental performance assessment methods has been established for pre-evaluating the ventilation and thermal environment after ecological city construction.Furthermore,as a pilot study,the methods have been applied in the ecological development demonstration area of the Yanqi Lake in the Huairou district.Based on the results,relevant construction policies are recommended. These,above all,will assist in the accumulation of practical experienceforamending and perfecting theevaluation methods.

    2.Research content

    2.1.Assessment content

    2.1.1.Ventilation assessment

    The ventilation environment of a city depends on the surface roughness,surrounding area openness,and underlying surface ventilation potential,and it is also affected by the background wind environment.Therefore,this study analyzed the surface ventilation potential and background wind environment,and based on the planned land functions,comprehensively assessed the ventilation environment after the construction of an ecological city.

    2.1.2.Urban thermal environment assessment

    Since the number of meteorological stations is currently limited,satellite remote sensing technology was used to obtain the current spatial distribution of heat islands.The heat island distribution after construction was forecast to determine heat island changes due to construction.

    2.2.Brief introduction to the research area

    The Yanqi Lake ecological development demonstration area,located in the southeast of the Huairou district,is a keyindustrial drive demonstration project in Beijing,and will be constructed as a high-end service,high-quality cultural and leisure area(Fig.1a).Land use types are indicated in Fig.1b, where yellow shading represents a change in land type from vegetation to buildings under future planning.The core area, which is called the ecological city area in Fig.1b,covers 21 km2,and the supporting town area covers 10 km2. Compared with the current land use,the public greenbelt and protective greenbelt area will increase around the Yanqi Lake, the decentralized residential and industrial land of the supporting towns in the south will be adjusted and arranged in a centralized manner,and the industrial and residential land area will be increased.

    3.Methods and data

    The study focused on the assessment of climate environmental changes before and after the ecological city construction.By integrating the existing technologies,a set of methods was established for assessing the effects of construction on the ventilation and thermal environment of the ecological city.

    3.1.Surface ventilation assessment

    Based on the analysis of the local background wind environment,a qualitative evaluation of the ventilation potential before and after the construction was conducted.Then,using the established ventilation potential index(VPI),a quantitative assessment was performed.

    3.1.1.Estimation of surface ventilation potential

    The surface ventilation potential is determined by ground coverage of buildings,surface roughness of natural vegetation, and access to the surrounding open area.The building coverage is the main negative factor that reduces air circulation,and the natural vegetation and access to the surrounding open area are the main positive factors that increase airfow (HKPD,2009).The surface roughness and open area can be quantitatively estimated by the aerodynamic roughness length and sky openness,respectively.

    3.1.1.1.Calculation of surface roughness.Surface roughness is quantitatively expressed by the dynamic roughness length, Z0.It is commonly calculated using the meteorological and morphological methods.Considering the data availability,Z0in this study was calculated with the morphological method, using remote sensing and geographic information data.As the ecological city area contains large areas of vegetation and town,Z0of the underlying surfaces of the vegetation and town areas should be estimated.

    In this study,Z0for the vegetation area was estimated based on the simple morphological method of Raupach(1994).

    Formula expression:

    Where Zdis zero plane displacement height(m);Z0is dynamic roughness length(m),Zhis roughness element height(m); Zd/Zhis normalized zero plane displacement height;Z0/Zhis normalized roughness length;λFis windward area ratio of roughness element on unit surface area;Λ is canopy area index of roughness element on unit surface area(for isotropy roughness element,Λ=2λF);ψhis roughness layer infuence function;Csis surface resistance coeffcient;CRis resistance coeffcient of roughness element;c is shading coeffcient of roughness element;b is canopy width of roughness element; Cdis an empirical coeffcient;Uhis wind speed;u*is friction velocity(or shearing velocity);and k is the Karman constant (=0.4).

    The seven model parameters,CR,Cs(or β),α,c,(u*/Uh)max(or γmax),ψh,and Λmax,were determined according to Jasinski et al.(2005).The vegetation canopy area index Λ was calculated based on Borak et al.(2005)and Yang et al.(2011) through the method defned by Zeng et al.(2001):

    Where Igis green leaf area index,namely,the common leaf area index;Isis withered leaf area index,including tree stem area index and withered and fallen leaf area index;Ingis the green leaf area index of month n;Is,minis the monthly minimum withered leaf area index of different vegetation coverage types;and 1-α is the monthly falling rate of withered leaves. The Is,minand α of different land use types in the research area were obtained according to the Is,minand α parameter tables for different MODIS IGBP vegetation cover types by Zeng et al.(2001)(Fig.1b).The green leaf area index Igwas directly obtained from the leaf area index via remote sensing inversion.In this study,the leaf area index was inverted from a Landsat-TM remote sensing image obtained during clear, summer conditions.

    The estimation of vegetation height,Zh,includes forest and plain farmland areas;the former changes slightly with the season,but for the latter,as the phenological phase changes, plant height changes relatively rapidly.In the study area, farmland in the plain area contains mainly wheat and corn,and the vegetation height is generally constant across the seasons; therefore,it can be obtained from the abundant measured data or literature.The forest vegetation height can be obtained via satellite radar inversion through the methods of Zhou et al. (2011).In the present study,we used forest vegetation height data with a 1-km spatial resolution inverted by the American satellite GLAS,and the corrections were performed based on the average height of broad-leaved forest,coniferous forest,bush wood,and other vegetation types determined by Cui et al.(2000)and Song(2010).

    For the urban area,Z0was estimated according to the urban morphological model established by Grimmond(1999):

    Where the different parameters are subjected to the vegetation morphological model above.The term λFis the windward area ratio of urban buildings per unit surface area.Also called the building sectional area index,λFis a complex variable involving building form;it is closely related to the urban building coverage,λP.In Grimmond(1999),the corresponding building coverage λPof the peak value of the normalized roughness length in the urban area simulated by multiple morphological methods is 0.3-0.4,and approximately 0.35 in most cases.According to Equations(9-11),and using 0.35 as the peak value of urban building coverage in the present study, the relationship between building sectional area index λFand urban building coverage λPwas back-deduced,and the obtained expression λF=0.8λPis relatively accurate.Thus,the change curve of the normalized roughness length,Z0/Zh,and the normalized zero plane displacement height,Zd/Zh,against building coverage,λP,can be simulated(refer to Fig.2).The normalized roughness length Z0/Zhreaches its peak value when λP=0.35.This is consistent with the dynamic roughness length in urban areas(Grimmond,1999).

    3.1.1.2.Calculation of sky openness.Sky openness is also called the sky view factor(SVF).It is a numerical value describing the 3D urban form;it refects the different geometrical forms of street channels in a city,affects the surface energy balance relationship,changeslocalaircirculation,acts as a comprehensive index for urban planning of street,community,and other local features,and serves as a reference for the scientifc planning of a city.More details can be found in Oke(1981).There are multiple methods for the acquisition of SVF,including the vector calculation and grid calculation models.G'al et al.(2008)compared the vector calculation and grid calculation models,and noted that the grid calculation model is more suitable for the rapid calculation of urban surface openness from large-scale and largevolume data.In this study,the grid calculation model based on the high-resolution digital elevation model(DEM)of Zakˇsek et al.(2011)was adopted to estimate the SVF.Fig.3 shows the calculation principle.The sky visual solid angle Ω and the normalized sky visual solid angle(i.e.,the sky openness,SVF),can be expressed by the following formula

    Where Ω is sky visual solid angle;γ is infuence of terrain height angle in azimuth angle i;n is number of calculated azimuth angles;and SVF is normalized sky visual solid angle, namely sky openness.In this article,SVF was estimated using a digital elevation grid image of Beijing with a 5-m spatial resolution.The azimuth angle number was n=36 and the radius of infuence was 20-grid(namely,20×5 m=100 m).

    3.1.2.Background wind feld analysis

    Wind direction and speed in the research area and its surroundings were obtained from meteorological stations.Statistical analysis was performed on the prevailing wind and mountain-valley breeze data.Land utilization data of the ecological city was updated in the numerical prediction model, and numerical simulation was performed for the background wind feld.The statistical analysis and numerical simulation results show that the background wind environment in the assessmentarea was,on the whole,fairly accurately represented.

    (1)Prevailing wind.The annual and seasonal prevailing winds were determined based on the observational data of wind direction and wind speed from the Huairou state meteorological station(30-year data from the year 1981 to 2010)and 16 automatic meteorological stations (2009-2013)in the Huairou district.

    (2)Mountain breeze and valley breeze.The predominant wind direction and start and end times of the mountain and valley breezes were determined using the data from the automatic meteorological station.The mountainous area presents a northeast-southwest direction.Wind with a northwest direction is considered a mountain breeze, and wind with a southeast direction is a valley breeze. The wind direction(eight-direction calculation method) with the maximum frequency during a certain period is considered the prevailing wind direction,i.e.,the period during which more than half the wind direction was southeast(northwest)was regarded as a valley(mountain)breeze period.

    (3)Wind speed.The annual wind speed was determined based on the monthly speed from 1981 to 2010,as recorded by the Huairou state-level meteorological station.The spatial distribution of wind speed was obtained through spatial interpolation of the observational data from the automatic meteorological stations.

    A 1 km×1 km spatial resolution wind feld covering a typical year and typical weather conditions in Beijing wassimulated using the WRF model via multi-nesting,and the wind feld and characteristics of the background wind feld near the ground were analyzed.

    (1)Model parameter setting.After the grid points were created,the urban land use data for Beijing in 2012 were input to the WRF model(the urban land is divided into high-density urban land,medium-density urban land, and low-density urban land based on the impervious surface percentage,and they correspond to Types 31,32, and 33,respectively,in the WRF model(refer to Fig.3). A 4-layer area nest of 27 km,9 km,3 km,and 1 km was adopted.The innermost nest area covers all of Beijing, and the physical process parameterization scheme, which has been investigated and widely used in Beijing, was selected.An urban canopy model was used to perform the wind feld numerical simulation with 1-km resolution.

    (2)Wind feld simulation in a typical weather case.One case is selected as mild wind pollution condition,that means low observed average wind speed and hazy weather;the other case is selected as strong wind diffusion condition,that means high observed average wind speed and no hazy weather;the third case is selected as ordinary weather condition,that means the observed wind speed is close to the monthly average wind speed and no hazy weather occurs.

    (3)Average wind feld simulation.In order to analyze the stable background wind environment in the Yanqi Lake region,10 m-high average wind felds for half a month in the winter of 2013 were simulated using the WRF model.Furthermore,the simulation results in the southeast of Huairou were compared with the average wind direction and speed during the simulation period recorded at the Huairou state-level meteorological station(Fig.4).

    3.1.3.Ventilation potential index

    The ventilation potential index(VPI)is a quantitative index used to compare the ventilation potential at different times and areas:

    Where m is ventilation potential grade,i is serial number of ventilation potential grade from none to high,n is grade of ventilation potential,wiis weight at grade i,and Piis the percentage at grade i.Generally,the value of VPI is between 0 and 1.Larger values indicate higher ventilation potential. Here m=5 and n=4.

    3.2.Urban thermal environment assessment

    We used UHI intensity to assess the thermal environment. Based on a satellite remote sensing method which compared spatial changes in the heat island intensity before and after constructionoftheecologicalcity,thethermalenvironmentwas qualitatively evaluated,the urban heat island potential index (UHPI)was established,and the thermal environment after construction of the ecological city was quantitatively assessed.

    3.2.1.Surface temperature inversion

    The surface temperature(120-m resolution)inversion was based on the methods of Jimenez-Munoz and Sobrino(2003) and Xu and Liu(2013)using Landsat-TM high-resolution remote sensing data.

    The 2011 land use is considered to be the current land utilization and the land use in the detailed plan in Fig.1b regarded as the planned land use.The surface temperature of both is determined using the above method.

    3.2.2.Estimation and grading of UHI intensity

    The average surface temperature in the farmland area is taken as the outskirts temperature,and its difference from the surface temperature in all other areas is defned as the UHI intensity.The specifc calculation is

    Where UHIIiis UHI intensity;Tiis temperature of pixel i; Tcropis temperature of a pixel in the farmland area;and n is total number of pixels in the farmland area.It is evident that the key point in the calculation of the heat island intensity is the selection of the outskirts farmland.Considering the general characteristics of the outskirts farmland,and the terrain and underlying surface features in the study area,the areas that met the following conditions were regarded as outskirts farmland in this study:

    (i)Plain(the altitude difference between city and plain is less than 50 m).

    (ii)Farmland type in the outskirts.

    (iii)Vegetation coverage≥80%.

    (iv)Imperviousness coverage≤20%.

    The imperviousness coverage and vegetation coverage were determined using the linear spectrum mixed decomposition method.The vegetation coverage was mainly estimated using the V-I-S-W model(Ridd,1995):

    The impervious stratum may be indicated by the sum of the proportion of low albedo end members and high albedo end members(Wu and Murry,2003):

    The key point in the above mentioned method is the selection of the end member(pure pixel).The minimum noise fraction(MNF)transformation was performed for the Landsat-TM image,and the pure pixel was selected based on the frst four MNF wave bands(89%of the information amount).Theimpervious surface coverage and vegetation coverage estimations were performed based on this method.

    After the outskirts farmland was determined,the UHI intensity was calculated and divided according to the categories in Table 1(Ye et al.,2011).

    3.2.3.UHPI

    The defnition of UHPI is Where UHPI is urban heat island potential index;m is heat island intensity grade;n is grade of temperature difference between the urban area and the outskirts;i is serial numberof n; viis weight at grade i;and qiis the percentage of grade i. Generally,the UHPI value is between 0 and 1.Larger values indicate more severe heat island effects.Here,m=7 and n=3.

    3.3.Essential data

    The essential data used in this study include urban planning,remote sensing,geographic information system,andmeteorological data.The specifc data,sources,and applications are shown in Table 2.

    Table 1 Grading basis of heat island intensity.

    4.Effect of the future construction on climate

    4.1.Ventilation assessment

    4.1.1.Assessment of surface ventilation potential

    Fig.5(left)shows the SVFs of the Yanqi Lake region at a spatial resolution of 100 m.It is obvious that the SVFs of the mountainous areas in the western and northern parts of the Yanqi Lake are lower than those of the plains,which are generally below 0.8.Most SVFs in the built-up regions of the plains are 0.5-0.8,and those of the remainder of the plains are generally above 0.9.The Landsat-TM(30-m)image for 2011 was used to estimate the leaf area index and extract the land use types.Then,the roughness length,Z0,of the regions surrounding the Yanqi Lake was determined at a resolution of 100 m(Fig.5,right).The roughness lengths of the surface in the northwestern mountainous areas are relatively long, generally above 0.5 m(above 0.8 m in some areas).In the built-up regions,the values are 0.6-0.7 m,and in the nonbuilt-up areas are generally less than 0.3 m.

    Based on the calculated roughness length and sky view factor,we calculated the surface ventilation potential in the Yanqi Lake region,and graded the calculation results based on the values in Table 3.Due to buildings in the mountainous areas,we did not consider,or grade,the ventilation potential in these areas.The ventilation potential after construction of the ecological city was determined based on roughness length and sky view factor.The distribution of the ventilation potential before and after construction shows that the overall ventilation potential of the ecological city will decrease due to the construction of groups of building over large areas.Therefore,it is necessary to construct ventilation corridors to guarantee ventilation of the environment.

    4.1.2.Background wind feld analysis

    Northwest and southeast winds prevail in the study area, and there are obvious mountain and valley winds.The durations of mountain and valley winds within a day in each season are different due to the different durations of sunshine.In winter,the duration of mountain winds is the longest;while in summer,the duration of valley winds is the longest.

    In winter,when the diffusion condition is relatively good, the area experiences prevailing northerly winds(although southerly winds occur during short periods of pollutedweather).In summer,when the diffusion condition is relatively good,the area mainly experiences southerly wind(although southerly wind with relatively low wind velocity also occurs during short periods of polluted weather).A mountain-valley breeze prevails under ordinary weather.The ecological city is located in the region where the plains and mountains meet, and the complicated terrain causes uneven wind velocity distribution around the ecological city,i.e.,the wind velocity in the plains is lower than that in the mountainous area in the northwest.

    4.1.3.Comprehensive assessment before and after construction

    Fig.6 shows the ventilation environment before and after the construction of the ecological city,as well as the difference between the two,based on the ventilation potential and background wind.Fig.6(right)shows that the ventilation potentials of three regions,location A in the northeast,the related town B in the south,and location C in front of the mountains around the Hongluo Lake,will decrease slightly. Location A is in the plains area,where northerly and northwesterly winds prevail.Therefore,the construction of manygroups of building will not only infuence the ventilation environment,but may prevent northern airfow from entering the core area of the Yanqi Lake,town B,and location C around the Hongluo Lake;thus deteriorating the overall ventilation in the ecological city.However,in town B,the regions with very high ventilation potential are signifcantly reduced,while the increase in the number of buildings in the core area of the Yanqi Lake is relatively small;after the construction of the ecological city,northwesterly airfow can also arrive at this small town.Therefore,ventilation corridors in the construction area can be created to effectively reduce the infuence on regions where the reserved ventilation potential is relatively high;thus,to some extent relieving the negative infuence of the construction on the ventilation environment of location B. In location C around the Hongluo Lake,there are large areas with very high ventilation potential before the construction of the ecological city.The construction will result in a large number of built-up regions in the southeastern part of the Hongluo Lake,which will infuence the ventilation potential. In particular,the buildings on the south shore are located along northwesterly wind channels and the lake-land breeze path, which may block their main airfow.

    Table 3 Grading of ventilation potential.

    The ventilation potentials before and after the construction of the ecological city area were graded according to the VPI (see Table 4 for the detailed grading).The results show that at present the VPI is 0.5345,and that the VPI after construction will be 0.4375.Although the reduction in ventilation potential is only slight after the construction,the ventilation potential is grade 3(general),which is relatively good.

    Table 4 Grading of the regional ventilation potential.

    4.2.Thermal environment assessment

    4.2.1.Heat island monitoring and prediction after construction

    According to the standards given in Table 1,the UHI intensity in different years is shown in Fig.7.The results indicate that the heat island effect in 2013 was slightly stronger than that in 2011.In particular,both heat island area and intensity increased slightly in the Yanqi industrial grouping region in the southeastern part of the core area.However,the UHIs in 2011 and 2013 were,on the whole,scattered. Compared with the current situation,both the heat island area and intensity will increase slightly after the construction. Furthermore,the two largest heat island regions,the northern part of the Huairou district and the Yanqi industrial grouping region,are likely to become linked together.Therefore,close attention needs to be paid to the changes in the heat islands in the Yanqi residential region,especially in the middle of both the regions mentioned above.The surface should be carefully developed and arranged,the number of water bodies and green belts should be increased,and ventilation corridors should be created to prevent the UHIs in the southern region of the Yanqi Lake from becoming linked together.However,after the construction of the ecological city,the heat island area and intensity change in the core area are relatively small,and heat islands are distributed,in a scattered manner,only in some regions of new building.

    4.2.2.Quantitative assessment of UHIs

    Based on the standards in Table 5,the results show that the UHPI of the core area in July 2011 and July 2013 were 0.0624 and 0.0878,respectively;placing them in the slight or none grade.Meanwhile,the UHPI predicted for heat islands after the construction is 0.1217;placing it in the general grade, which is still relatively good.However,we can imagine that in other areas,especially in the two largest heat island regions, conditions are likely to become worse.

    5.Conclusions and suggestions

    5.1.Effect on the climate

    5.1.1.Ventilation

    The VPI of the study area will decrease from the current 0.5345 to 0.4375 after the construction of the ecological city, and the ventilation capability will also weaken.However,the assessment grade of the ventilation potential will remain at the general grade,which is relatively good.It is necessary to point out that after construction the ventilation potential will markedly decrease for a relatively large area with buildings in the northeastern part of the Yanqi Lake,the related town in the southern part of the Yanqi Lake,and the area in front of the mountains around the Hongluo Lake.

    5.1.2.Thermal environment

    The UHPIs for July 2011,July 2013,and after the construction of the ecological city are 0.0624,0.0878,and 0.1217,respectively.The heat island grade in 2013 changed to general from the original slight or none grade;however, the general grade is still relatively good.Furthermore,the heat scope and intensity changes in the core area of the ecological city are relatively small.It is necessary to point out that after the construction of the ecological city,both the heat island area and intensity will increase slightly as compared with the present situation,and that the two largestheat island regions,i.e.,the northern part of the Huairou district and Yanqi industrial grouping region,may become linked together.

    5.2.Development suggestions

    5.2.1.Constructing urban ventilation corridors

    The land use plan for the regions around the Yanqi Lake ecological development demonstration area shows a relatively large open space in the northwestern region of the main construction land,and multiple main roads and river channels, oriented from south to north,are also present.These river channels,main roads,and green corridors can be utilized as ventilation corridors(Matzarakis and Mayer,1992)so that the fresh air generated in the open space in the northwestern region of the built-up area(cold source)and the northwesterly background wind can fow into the constructed region in the southeastern region of the ecological city to relieve the negative infuence of the construction on the local weather and atmospheric quality(as shown in Fig.8).

    Table 5 Grading of the regional heat islands.

    5.2.2.Relieving UHIs

    The thermal assessment shows that after construction of the ecological city both the area and intensity of the heat islands in the Yanqi industrial grouping region will increase slightly as compared with the current situation.However,the two largest heat island regions,the northern part of the Huairou district and the Yanqi industrial grouping region,may become linked together.This linking should be prevented to avoid the relatively strong heat islands developing into a single,strongerisland.Furthermore,land should be reserved to serve as a cold source.

    5.3.Future research

    This research started by evaluating the climatic environment output effects of ecological city construction inputs. Qualitative and quantitative assessments of the local climatic environment before and after the construction of the Yanqi Lake ecological development demonstration area were performed and relevant policy suggestions were formulated.In the future,the following items will require further study:

    5.3.1.Establish more comprehensive assessments

    This research mainly focused on assessing the ventilation and thermal environments before and after construction of the ecological city.In the future,it will be necessary to establish more comprehensive assessments of the climatic environment covering the urban atmosphere and urban meteorological disasters.

    5.3.2.Improve the assessment methods

    This study utilized remote sensing,geographical information systems,meteorological numerical simulation methods, and other technologies that are relatively well developed.As technology progresses,remote sensing with higher resolution and precision,and the latest assessment methods and technologies(e.g.,a 100-m resolution meteorological numerical model),will need to be employed.

    5.3.3.Expand the scope of application

    The Yanqi Lake ecological development demonstration area belongs to a new breed of newly-built ecological cities on the outskirts of an urban area.In the future,assessments should be made for the construction of ecological cities in central urban areas and the reconstruction of old ecological cities with improved research methods and technologies, enhancing the applicability of the evaluation methods.

    5.3.4.Form relevant assessment technical specifcations

    Relevant research data,methods,and techniques should be extracted and summarized into technical specifcations for the assessment of the effects of ecological urban construction on climate.Based on these specifcations,planners can visualize the basic climate environment of the city and optimize the design before construction.

    Acknowledgements

    This work was sponsored by Beijing Municipal Science and Technology Project(Z131100001113026),the Program of the ResearchandInnovationTeam on UrbanClimate AssessmentofBeijing MeteorologicalBureau,Climate Change Special Foundation of China Meteorology Administration(CCSF201506),Science and technology project of Beijing Meteorological Bureau(BMBKJ201402002)and National Natural Science Foundation of China(71473146).

    Birgit,K.,Malcolm,C.,Kevin,L.,2007.Environmental performance of a naturally ventilated city centre library.Energy Buildings 39,792-801.

    Borak,J.S.,Jasinski,M.F.,Crago,R.,2005.Time series vegetation aerodynamic roughness felds estimated from MODIS observations.Agric.For. Meteorol.135(1-4),252-268.

    Cui,G.-F.,Cheng,K.-W.,Shen,G.-Z.,et al.,2000.Status and classifcation of vegetation in Labagoumen forest region in Beijing.J.Beijing For.Univ.22 (4),46-51(in Chinese).

    G'al,T.,Lindberg,F.,Unger,J.,2008.Comporting continues sky view factor using 3D urban raster and vector databases:comparison and application to urban climate.Theor.Appl.Climatol.95(1-2),111-123.

    Grimmond,C.S.B.,1999.Aerodynamic properties of urban areas derived from analysis of surface form.J.Appl.Meteorol.38,1262-1291.

    HKPD(Hong Kong Planning Department),2009.Hong Kong planning standards and guidelines.CPLD49/05.Committee on Planning and Land Development(in Chinese).

    Jasinski,M.F.,Borak,J.,Crago,R.,2005.Bulk surface momentum parameters for satellite-derived vegetation felds.Agric.For.Meteorol.133,55-68.

    Jimenez-Munoz,J.C.,Sobrino,J.A.,2003.A generalized single channel method for retrieving land surface temperature from remote sensing data. J.Geophys.Res.108(D22),4688-4695.

    Lei,X.-P.,He,L.-B.,Wu,C.-L.,2011.Research methods and improvement measure research progress of urban ventilation.Build.Energy Environ.30 (1),11-17(in Chinese).

    Li,L.,Wu,D.,Zhang,L.-J.,et al.,2012.Ventilation assessment on urbanblock detailed planning based on numerical simulation.Acta Sci.Circumstantiae 32(4),946-953(in Chinese).

    Li,X.X.,Koh,T.Y.,Dara,E.,et al.,2013.A multi-resolution ensemble study of a tropical urban environment and its interactions with the background regional atmosphere.J.Geophys.Res.Atmos.118(17),9804-9818.

    Matzarakis,A.,Mayer,H.,1992.Mapping of urban air paths for planning in Munchen.WissBer Inst.Meteorol.Klimaforsch.Univ.Karlsruhe 16, 13-22.

    Morhmald,B.,Shinsuke,K.,2008.Assessment of wind ventilation performance inside a street located within a dense urban area using exceedance probability.Seisan-Kenkyu 60(1),5-13.

    Oke,T.R.,1981.Canyon geometry and the nocturnal heat island:comparison of scale model and feld observations.J.Climatol.1,237-254.

    Peng,R.-Q.,2008.Wind feld numerical simulation in building environment.J. Xi Hua Univ.Nat.Sci.Ed.27(1),50-54(in Chinese).

    Qin,Z.,Karnieli,A.,Berliner,P.,2001.A mono-window algorithm for retrieving land surface temperature from Landsat TM data and its application to the Israel-Egypt border region.Int.J.Remote Sens.22(18), 3719-3746.

    Raupach,M.R.,1994.Simplifed expressions for vegetation roughness length and zero-plane displacement as a function of canopy height and area index. Bound.-Layer Meteor.71,211-216.

    Ridd,M.K.,1995.Exploring a V-I-S(vegetation-impervious surface-soil) model for urban ecosystem analysis through remote sensing:comparative anatomy for cities.Int.J.Remote Sens.16(2),2165-2185.

    Song,X.-Y.,2010.Research on structural features of main arbors of forest for conservation of water supply at the upstream of Miyun reservoir:taking Chaoguan west ditch forest farm for example.Soil Water 4,54-57(in Chinese).

    Tanzina,M.,William,A.G.,2012.Characterization and estimation of urban heat island at Toronto:impact of the choice of rural sites.Theor.Appl. Climatol.108(1-2),105-117.

    Wu,U.C.S.,Murry,A.T.,2003.Estimating impervious surface distribution by spectral mixture analysis.Remote Sens.Environ.84(4),493-505.

    Xu,Y.-M.,Liu,Y.-H.,2013.Study on the thermal environment and its relationship with impervious surface in Beijing city using TM image.Ecol. Environ.Sci.22(4),639-643(in Chinese).

    Yang,A.-Q.,Sun,G.-Q.,Lu,L.-X.,et al.,2011.Deriving aerodynamic roughness length and zero-plane displacement height from MODIS product for Eastern China.J.Meteorol.Sci.31(4),516-524(in Chinese).

    Ye,C.-H.,Liu,Y.-.H.,Liu,W.-D.,et al.,2011.Research on urban surface heat environment monitoring indexes and its application.Meteorol.Sci.Technol.39(1),95-101(in Chinese).

    Zakˇsek,K.,Oˇstir,K.,Kokalj,ˇZ.,2011.Sky-view factor as a relief visualization technique.Remote Sens.3,398-415.

    Zeng,X.-B.,Shaikh,M.,Dai,Y.-J.,et al.,2001.Coupling of the common land model to the NCAR community climate model.J.Clim.15, 1832-1854.

    Zhang,B.-Y.,Sang,J.-G.,Wu,G.-C.,2004.Features and simulation of wind feld of building group environment:research on wind environment simulation.Part I.Mech.Eng.26(3),1-9(in Chinese).

    Zhang,L.,Zhang,Y.,2013.Discussion of the quantitative assessment method of ecological planning performance.Environ.Sci.Manage.38(3),27-30 (in Chinese).

    Zhao,L.-H.,Zheng,Z.-J.,Li,L.,2009.Discussion of the evaluation indexes of outdoor wind environment quality for the enclosed residential quarter in Guangzhou.In:Proceedings of the 5th International Workshop on Energy and Environment of Residential Buildings/3rd International Conference on Built Environment and Public Health,Vol.1&2,pp.1264-1270.

    Zhou,Z.-Q.,Yue,C.-R.,Xu,T.-S.,et al.,2011.Review on forest height estimation using remote sensing.Mod.Agric.Sci.Technol.2,198-199(in Chinese).

    Received 20 March 2015;revised 15 May 2015;accepted 6 August 2015

    Available online 21 August 2015

    *Corresponding author.

    E-mail address:ccyx0817@126.com(CHENG C.).

    Peer review under responsibility of National Climate Center(China Meteorological Administration).

    Production and Hosting by Elsevier on behalf of KeAi

    http://dx.doi.org/10.1016/j.accre.2015.08.001

    1674-9278/Copyright?2015,National Climate Center(China Meteorological Administration).Production and hosting by Elsevier B.V.on behalf of KeAi. This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    一进一出好大好爽视频| 午夜影院日韩av| 国产成人av激情在线播放| 制服诱惑二区| 久久人妻福利社区极品人妻图片| 国产av一区二区精品久久| 涩涩av久久男人的天堂| 亚洲 欧美 日韩 在线 免费| 在线国产一区二区在线| 午夜久久久在线观看| tocl精华| 黑丝袜美女国产一区| 亚洲第一av免费看| 欧美日韩黄片免| 9191精品国产免费久久| 久久草成人影院| 国产精品国产av在线观看| 动漫黄色视频在线观看| 亚洲免费av在线视频| 成人手机av| 啦啦啦 在线观看视频| 丁香欧美五月| 黄网站色视频无遮挡免费观看| 看片在线看免费视频| 91麻豆av在线| 午夜精品在线福利| 久久99一区二区三区| 亚洲男人天堂网一区| 天天躁夜夜躁狠狠躁躁| 欧美性长视频在线观看| 91老司机精品| 亚洲免费av在线视频| 色精品久久人妻99蜜桃| 国产精品一区二区三区四区久久 | 在线观看免费高清a一片| av网站在线播放免费| 国产精华一区二区三区| 极品教师在线免费播放| 日韩欧美一区二区三区在线观看| 日本黄色视频三级网站网址| 亚洲av片天天在线观看| 午夜久久久在线观看| 国产精品 欧美亚洲| 亚洲成a人片在线一区二区| 国产精品 国内视频| 欧美成人性av电影在线观看| 国产一区二区在线av高清观看| 美女午夜性视频免费| 欧美日韩亚洲高清精品| 男男h啪啪无遮挡| 午夜成年电影在线免费观看| 桃红色精品国产亚洲av| 人妻久久中文字幕网| 国产精品久久视频播放| 91成年电影在线观看| 亚洲成人国产一区在线观看| av天堂在线播放| 欧美日韩精品网址| 久久国产精品影院| 又紧又爽又黄一区二区| 高清欧美精品videossex| 精品卡一卡二卡四卡免费| 精品国产乱子伦一区二区三区| av天堂久久9| 黄片播放在线免费| 日韩欧美在线二视频| 很黄的视频免费| 大香蕉久久成人网| 日韩人妻精品一区2区三区| 午夜福利在线观看吧| 性色av乱码一区二区三区2| 在线免费观看的www视频| 久久国产精品男人的天堂亚洲| 少妇粗大呻吟视频| 色老头精品视频在线观看| 一级毛片精品| svipshipincom国产片| 97超级碰碰碰精品色视频在线观看| 亚洲精品中文字幕一二三四区| 亚洲精品粉嫩美女一区| 另类亚洲欧美激情| 亚洲美女黄片视频| 嫩草影视91久久| 精品一区二区三区视频在线观看免费 | 亚洲人成77777在线视频| 日韩三级视频一区二区三区| 黄色女人牲交| avwww免费| 国产在线观看jvid| 久久欧美精品欧美久久欧美| 免费在线观看影片大全网站| 亚洲精品中文字幕在线视频| 亚洲欧美日韩另类电影网站| 九色亚洲精品在线播放| 欧美性长视频在线观看| 久久精品亚洲精品国产色婷小说| 高清欧美精品videossex| 男女高潮啪啪啪动态图| 精品熟女少妇八av免费久了| 亚洲人成电影免费在线| 欧美乱妇无乱码| 午夜福利在线观看吧| 免费女性裸体啪啪无遮挡网站| 啦啦啦 在线观看视频| 1024香蕉在线观看| 亚洲熟妇中文字幕五十中出 | 大码成人一级视频| 精品久久久久久久久久免费视频 | 欧美性长视频在线观看| 欧美乱色亚洲激情| 老司机亚洲免费影院| 国产亚洲精品一区二区www| 窝窝影院91人妻| 黄色a级毛片大全视频| 黄色丝袜av网址大全| 国产亚洲欧美精品永久| 国产区一区二久久| 亚洲五月婷婷丁香| 最近最新免费中文字幕在线| 最近最新中文字幕大全免费视频| 1024视频免费在线观看| 叶爱在线成人免费视频播放| 久久人妻熟女aⅴ| 久久亚洲精品不卡| 99久久久亚洲精品蜜臀av| 超碰97精品在线观看| 久热爱精品视频在线9| 一二三四在线观看免费中文在| 国产精品1区2区在线观看.| 水蜜桃什么品种好| 日本黄色视频三级网站网址| 欧美日韩黄片免| 又紧又爽又黄一区二区| 国产精品1区2区在线观看.| a级毛片黄视频| 成人精品一区二区免费| 久久青草综合色| 天堂影院成人在线观看| 黑人猛操日本美女一级片| 亚洲avbb在线观看| 亚洲成国产人片在线观看| 亚洲人成电影免费在线| 久久精品亚洲av国产电影网| 性色av乱码一区二区三区2| 亚洲av片天天在线观看| 亚洲五月天丁香| 99国产精品免费福利视频| 国产免费av片在线观看野外av| 亚洲成人精品中文字幕电影 | av天堂在线播放| 日本免费一区二区三区高清不卡 | 老司机在亚洲福利影院| 久久人人爽av亚洲精品天堂| 亚洲欧洲精品一区二区精品久久久| 高清在线国产一区| 国产1区2区3区精品| 欧美 亚洲 国产 日韩一| 欧美乱码精品一区二区三区| 黑人操中国人逼视频| 在线观看舔阴道视频| 日韩精品中文字幕看吧| 久久国产精品影院| 久久国产精品影院| 亚洲三区欧美一区| 国产精品一区二区精品视频观看| 亚洲精品成人av观看孕妇| 一二三四在线观看免费中文在| 日韩国内少妇激情av| 老熟妇乱子伦视频在线观看| 欧美乱色亚洲激情| 欧美性长视频在线观看| 黄色毛片三级朝国网站| 窝窝影院91人妻| 国产精品野战在线观看 | 国产av在哪里看| 久久欧美精品欧美久久欧美| av天堂在线播放| 一a级毛片在线观看| 一区二区日韩欧美中文字幕| 久久天堂一区二区三区四区| 午夜免费观看网址| 深夜精品福利| 露出奶头的视频| 久久精品人人爽人人爽视色| 黄色视频,在线免费观看| 91成年电影在线观看| 亚洲av熟女| 欧美日韩视频精品一区| 黑人猛操日本美女一级片| 伦理电影免费视频| 高潮久久久久久久久久久不卡| xxx96com| 看黄色毛片网站| 欧美人与性动交α欧美软件| 一二三四在线观看免费中文在| 动漫黄色视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 一区二区三区激情视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲伊人色综图| av有码第一页| 满18在线观看网站| 制服诱惑二区| 免费不卡黄色视频| 美女扒开内裤让男人捅视频| 亚洲欧洲精品一区二区精品久久久| 黄色a级毛片大全视频| 亚洲国产欧美网| 午夜福利影视在线免费观看| 欧美日韩av久久| 亚洲 欧美 日韩 在线 免费| 国产精品综合久久久久久久免费 | 欧美黄色淫秽网站| 亚洲avbb在线观看| 久久久久久久久免费视频了| 色综合婷婷激情| 757午夜福利合集在线观看| 波多野结衣高清无吗| 午夜福利欧美成人| 亚洲成av片中文字幕在线观看| 亚洲中文av在线| 日本免费a在线| 亚洲五月天丁香| 别揉我奶头~嗯~啊~动态视频| 热99re8久久精品国产| 精品午夜福利视频在线观看一区| 日日摸夜夜添夜夜添小说| 变态另类成人亚洲欧美熟女 | 色播在线永久视频| 又紧又爽又黄一区二区| 每晚都被弄得嗷嗷叫到高潮| 国产野战对白在线观看| 国产免费av片在线观看野外av| 亚洲,欧美精品.| 免费高清在线观看日韩| 日韩av在线大香蕉| 成年人黄色毛片网站| 最新美女视频免费是黄的| 丝袜美足系列| 亚洲 欧美 日韩 在线 免费| 亚洲全国av大片| 久久香蕉精品热| 国产无遮挡羞羞视频在线观看| 精品一区二区三卡| 国产精品爽爽va在线观看网站 | 色婷婷av一区二区三区视频| 日日爽夜夜爽网站| 欧美国产精品va在线观看不卡| 欧美成人午夜精品| 嫩草影视91久久| 午夜福利一区二区在线看| 国产亚洲欧美在线一区二区| 男女下面插进去视频免费观看| 真人一进一出gif抽搐免费| 女人被狂操c到高潮| 91老司机精品| 美国免费a级毛片| 97人妻天天添夜夜摸| 97碰自拍视频| 午夜免费观看网址| 中文字幕另类日韩欧美亚洲嫩草| a在线观看视频网站| 少妇的丰满在线观看| 黄片播放在线免费| 色老头精品视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 午夜免费观看网址| 日日干狠狠操夜夜爽| 18禁黄网站禁片午夜丰满| 黑人欧美特级aaaaaa片| 欧美不卡视频在线免费观看 | 亚洲专区国产一区二区| 国产在线精品亚洲第一网站| 中文字幕av电影在线播放| 麻豆成人av在线观看| 日本欧美视频一区| 国产精品日韩av在线免费观看 | 丝袜美腿诱惑在线| 国产欧美日韩综合在线一区二区| 日本三级黄在线观看| 国产亚洲精品综合一区在线观看 | 在线天堂中文资源库| 色精品久久人妻99蜜桃| 老司机午夜福利在线观看视频| av国产精品久久久久影院| 九色亚洲精品在线播放| 最好的美女福利视频网| 欧美老熟妇乱子伦牲交| 性色av乱码一区二区三区2| 日本一区二区免费在线视频| 久久人人爽av亚洲精品天堂| 九色亚洲精品在线播放| 午夜影院日韩av| 国产欧美日韩精品亚洲av| 人成视频在线观看免费观看| 成人三级做爰电影| 中文字幕精品免费在线观看视频| 欧美av亚洲av综合av国产av| 亚洲aⅴ乱码一区二区在线播放 | 无人区码免费观看不卡| 黄片播放在线免费| 一二三四在线观看免费中文在| 国产精品久久久av美女十八| 亚洲人成77777在线视频| 久久九九热精品免费| 久久草成人影院| 一本大道久久a久久精品| 欧美日韩亚洲高清精品| 在线免费观看的www视频| 日日干狠狠操夜夜爽| 精品无人区乱码1区二区| 日日夜夜操网爽| 麻豆成人av在线观看| 精品一品国产午夜福利视频| 99国产综合亚洲精品| 高清毛片免费观看视频网站 | 国产高清videossex| 美国免费a级毛片| 好男人电影高清在线观看| 99精品欧美一区二区三区四区| 亚洲av熟女| 欧美中文综合在线视频| 一边摸一边抽搐一进一出视频| 人人妻,人人澡人人爽秒播| 日本精品一区二区三区蜜桃| 国产精品 国内视频| xxxhd国产人妻xxx| 12—13女人毛片做爰片一| 国产午夜精品久久久久久| 99久久国产精品久久久| 成人精品一区二区免费| 日韩免费av在线播放| 国产一区二区在线av高清观看| 国产国语露脸激情在线看| 欧美人与性动交α欧美软件| 亚洲成人精品中文字幕电影 | 满18在线观看网站| 欧美丝袜亚洲另类 | 日韩视频一区二区在线观看| 男人操女人黄网站| 丰满迷人的少妇在线观看| 国产精品99久久99久久久不卡| 天天躁狠狠躁夜夜躁狠狠躁| 视频区欧美日本亚洲| 久久久国产成人免费| 久久人人爽av亚洲精品天堂| 欧美日韩精品网址| 午夜a级毛片| 18禁美女被吸乳视频| 欧美黄色淫秽网站| 黑人操中国人逼视频| 99久久人妻综合| 国产精品二区激情视频| 大香蕉久久成人网| 国产日韩一区二区三区精品不卡| 国产精品免费视频内射| 久久精品aⅴ一区二区三区四区| 男女下面插进去视频免费观看| 少妇的丰满在线观看| 搡老熟女国产l中国老女人| 亚洲专区中文字幕在线| 天堂动漫精品| 久久精品亚洲av国产电影网| 亚洲欧洲精品一区二区精品久久久| 欧美人与性动交α欧美软件| avwww免费| 国产精品 国内视频| 91大片在线观看| 久久精品91蜜桃| 国产一区在线观看成人免费| 性色av乱码一区二区三区2| 美女 人体艺术 gogo| 欧美日韩国产mv在线观看视频| 日韩成人在线观看一区二区三区| 精品久久久久久久毛片微露脸| 国产欧美日韩综合在线一区二区| 国产99久久九九免费精品| 午夜视频精品福利| 嫩草影院精品99| 亚洲第一欧美日韩一区二区三区| 国产精品久久久久成人av| 天天躁狠狠躁夜夜躁狠狠躁| 免费观看人在逋| 男女高潮啪啪啪动态图| 亚洲欧美激情在线| 久久久水蜜桃国产精品网| 99在线视频只有这里精品首页| 91精品国产国语对白视频| 999久久久国产精品视频| 精品久久久久久久毛片微露脸| 国产麻豆69| 热re99久久国产66热| 99国产极品粉嫩在线观看| 国产精品av久久久久免费| 久久精品国产99精品国产亚洲性色 | 精品国产美女av久久久久小说| 一夜夜www| 亚洲国产中文字幕在线视频| 亚洲精品av麻豆狂野| 一进一出抽搐gif免费好疼 | 精品熟女少妇八av免费久了| 麻豆一二三区av精品| 亚洲成人精品中文字幕电影 | 欧美中文日本在线观看视频| av免费在线观看网站| 久久久久久久午夜电影 | 亚洲精品一区av在线观看| 最近最新中文字幕大全免费视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲三区欧美一区| 天堂√8在线中文| 精品一区二区三区四区五区乱码| 无限看片的www在线观看| www.999成人在线观看| 日韩欧美在线二视频| 精品国产亚洲在线| 亚洲七黄色美女视频| 日日爽夜夜爽网站| 成人三级黄色视频| 丝袜人妻中文字幕| 国产三级黄色录像| 在线av久久热| 国产精品98久久久久久宅男小说| 午夜福利,免费看| 99久久国产精品久久久| 一本综合久久免费| 日本 av在线| 国产熟女午夜一区二区三区| 色播在线永久视频| 99re在线观看精品视频| 国产黄色免费在线视频| 亚洲 欧美一区二区三区| 欧美+亚洲+日韩+国产| 精品免费久久久久久久清纯| 久久伊人香网站| 亚洲精品中文字幕在线视频| 国产精品自产拍在线观看55亚洲| 高清av免费在线| 亚洲熟女毛片儿| 国产免费现黄频在线看| 狂野欧美激情性xxxx| 国产一区二区三区视频了| 一夜夜www| 欧美最黄视频在线播放免费 | 亚洲精品久久成人aⅴ小说| 国产单亲对白刺激| 成人18禁在线播放| 天堂影院成人在线观看| 亚洲专区中文字幕在线| 日韩人妻精品一区2区三区| 99精品在免费线老司机午夜| 国产xxxxx性猛交| 午夜精品在线福利| 免费在线观看视频国产中文字幕亚洲| 亚洲av美国av| 午夜久久久在线观看| 男女午夜视频在线观看| 十八禁人妻一区二区| 欧美日韩国产mv在线观看视频| 日韩欧美一区二区三区在线观看| 精品熟女少妇八av免费久了| 欧美精品啪啪一区二区三区| 国产成年人精品一区二区 | 又黄又粗又硬又大视频| 欧美午夜高清在线| 国产熟女午夜一区二区三区| 欧美 亚洲 国产 日韩一| 亚洲,欧美精品.| 国产亚洲欧美98| 欧美成人午夜精品| 国产精品爽爽va在线观看网站 | av电影中文网址| 亚洲成人国产一区在线观看| 免费在线观看亚洲国产| 麻豆成人av在线观看| 免费高清在线观看日韩| 天堂影院成人在线观看| 亚洲精品一区av在线观看| 99riav亚洲国产免费| av电影中文网址| 亚洲av美国av| 国产成人欧美在线观看| 欧美日韩黄片免| 久久精品国产亚洲av高清一级| 淫妇啪啪啪对白视频| 在线观看一区二区三区激情| 亚洲五月婷婷丁香| 欧美 亚洲 国产 日韩一| 免费高清视频大片| 色综合欧美亚洲国产小说| 色哟哟哟哟哟哟| 一级片免费观看大全| 91九色精品人成在线观看| 国产午夜精品久久久久久| 黑人欧美特级aaaaaa片| av欧美777| 亚洲欧美激情综合另类| 亚洲精品中文字幕在线视频| 啪啪无遮挡十八禁网站| 欧美 亚洲 国产 日韩一| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲专区国产一区二区| 亚洲性夜色夜夜综合| 人妻丰满熟妇av一区二区三区| 亚洲国产看品久久| aaaaa片日本免费| 大型黄色视频在线免费观看| 日韩大码丰满熟妇| 国产亚洲精品一区二区www| 欧美黄色片欧美黄色片| 亚洲欧美精品综合久久99| 免费在线观看影片大全网站| 一个人观看的视频www高清免费观看 | 9色porny在线观看| 久久天堂一区二区三区四区| 久久久久精品国产欧美久久久| 久久欧美精品欧美久久欧美| 18禁裸乳无遮挡免费网站照片 | 亚洲 国产 在线| 男女之事视频高清在线观看| 黄色怎么调成土黄色| av视频免费观看在线观看| 十分钟在线观看高清视频www| 不卡一级毛片| 国产成人精品无人区| 热re99久久国产66热| 岛国视频午夜一区免费看| 一级片'在线观看视频| ponron亚洲| 99riav亚洲国产免费| 老熟妇乱子伦视频在线观看| 99国产综合亚洲精品| 午夜福利影视在线免费观看| 精品久久久久久久毛片微露脸| 淫妇啪啪啪对白视频| 99热只有精品国产| 久久久久久亚洲精品国产蜜桃av| 韩国精品一区二区三区| 国产熟女午夜一区二区三区| 国产高清视频在线播放一区| 桃红色精品国产亚洲av| av免费在线观看网站| 精品国产乱码久久久久久男人| 91九色精品人成在线观看| 欧美日韩亚洲综合一区二区三区_| 真人做人爱边吃奶动态| 亚洲精品美女久久久久99蜜臀| 成人影院久久| 国产av一区二区精品久久| 自拍欧美九色日韩亚洲蝌蚪91| 91成人精品电影| 久久这里只有精品19| 午夜精品国产一区二区电影| 亚洲色图 男人天堂 中文字幕| 国产一区二区三区视频了| 成年人黄色毛片网站| 色哟哟哟哟哟哟| 一区二区三区精品91| 一二三四在线观看免费中文在| 欧美成人性av电影在线观看| 无限看片的www在线观看| 嫩草影院精品99| 国产成人精品久久二区二区免费| 黄色视频不卡| 精品一区二区三区av网在线观看| 在线永久观看黄色视频| 校园春色视频在线观看| 日本欧美视频一区| 久久久国产成人精品二区 | 久久久久久久精品吃奶| 欧美日韩瑟瑟在线播放| 久久人妻熟女aⅴ| 亚洲专区字幕在线| 俄罗斯特黄特色一大片| 色综合婷婷激情| 亚洲av日韩精品久久久久久密| 亚洲 欧美一区二区三区| 一二三四社区在线视频社区8| 成人影院久久| 午夜免费观看网址| 中文字幕最新亚洲高清| 另类亚洲欧美激情| www.熟女人妻精品国产| 婷婷六月久久综合丁香| 校园春色视频在线观看| 男人舔女人下体高潮全视频| 亚洲成国产人片在线观看| 涩涩av久久男人的天堂| 极品人妻少妇av视频| av天堂在线播放| 桃色一区二区三区在线观看| 男女午夜视频在线观看| e午夜精品久久久久久久| 国产男靠女视频免费网站| 亚洲成a人片在线一区二区| 精品人妻1区二区| e午夜精品久久久久久久| 最好的美女福利视频网| 久久久久久亚洲精品国产蜜桃av| 国产野战对白在线观看| 最新在线观看一区二区三区| 亚洲成人精品中文字幕电影 | 免费在线观看日本一区| 欧美黄色片欧美黄色片| 黄色毛片三级朝国网站| 久久久久久大精品| 老熟妇仑乱视频hdxx| 欧美日韩福利视频一区二区| 亚洲 国产 在线| 国产麻豆69| 村上凉子中文字幕在线| 性欧美人与动物交配| 久久国产精品男人的天堂亚洲| 法律面前人人平等表现在哪些方面| av欧美777|