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

    Increases in Anthropogenic Heat Release from Energy Consumption Lead to More Frequent Extreme Heat Events in Urban Cities

    2021-04-20 04:01:44BinLIUZhenghuiXIEPeihuaQINShuangLIURuichaoLILonghuanWANGYanWANGBinghaoJIASiCHENJinboXIEandChunxiangSHI
    Advances in Atmospheric Sciences 2021年3期

    Bin LIU, Zhenghui XIE, Peihua QIN, Shuang LIU, Ruichao LI,3, Longhuan WANG,3,Yan WANG,3, Binghao JIA, Si CHEN,3, Jinbo XIE, and Chunxiang SHI

    1State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China

    2School of Software Engineering, Chengdu University of Information Technology, Chengdu 610103, China

    3University of Chinese Academy of Sciences, Beijing 100049, China

    4Key Laboratory of Mountain Hazards and Earth Surface Processes, Institute of Mountain Hazards and Environment,Chinese Academy of Sciences, Chengdu 610041, China

    5National Meteorological Information Center, China Meteorological Administration, Beijing 100081, China

    (Received 13 May 2020; revised 21 August 2020; accepted 9 October 2020)

    ABSTRACT With economic development and rapid urbanization, increases in Gross Domestic Product and population in fastgrowing cities since the turn of the 21st Century have led to increases in energy consumption. Anthropogenic heat flux released to the near-surface atmosphere has led to changes in urban thermal environments and severe extreme temperature events. To investigate the effects of energy consumption on urban extreme temperature events, including extreme heat and cold events, a dynamic representation scheme of anthropogenic heat release (AHR) was implemented in the Advanced Research version of the Weather Research and Forecasting (WRF) model, and AHR data were developed based on energy consumption and population density in a case study of Beijing, China. Two simulations during 1999-2017 were then conducted using the developed WRF model with 3-km resolution with and without the AHR scheme. It was shown that the mean temperature increased with the increase in AHR, and more frequent extreme heat events were produced, with an annual increase of 0.02-0.19 days, as well as less frequent extreme cold events, with an annual decrease of 0.26-0.56 days,based on seven extreme temperature indices in the city center. AHR increased the sensible heat flux and led to surface energy budget changes, strengthening the dynamic processes in the atmospheric boundary layer that reduce AHR heating efficiency more in summer than in winter. In addition, it was concluded that suitable energy management might help to mitigate the impact of extreme temperature events in different seasons.

    Key words: anthropogenic heat release, extreme temperature event, Weather Research and Forecasting model, Beijing

    1. Introduction

    The level of urbanization has increased from 39% to 55% in the last four decades globally (Chen et al., 2014b;United Nations, 2018). Vast numbers of people move from rural areas to urban areas across the world with about 2%energy consumption growth annually during the 21st century (Enerdata, 2018). Anthropogenic heat flux is mainly released from energy consumption supporting residents’accommodation and transportation, leading to changes in the urban thermal environment and extreme heat events(Feng et al., 2012; Chen et al., 2014a). Anthropogenic heat release (AHR) increases the occurrence of extreme temperature by heating the surface atmosphere directly (Salamanca et al., 2014) and functions as the main factor in the urban heat island effect, which has synergistic interactions with heat waves (Zhao et al., 2018; Rizvi et al., 2019). The number of extreme temperature events in the first decade of the 21st century grew faster than in the last 10 years of the 20th century (WMO, 2013). From the news report of Beijing of China on 3 August 2018, the number of hot days up to that date had already exceeded the total number of hot days in 2017, and the temperature reached a maximum of 42.1°C(Sina, 2018). Besides, increased temperatures can substantially increase the rate of temperature-related illnesses(Schr?ter et al., 2005; Jenerette et al., 2011; Tian et al.,2012; Zhang et al., 2014) and impose a heavy water supply burden on the urban ecosystem (Huang et al., 2010; Jenerette et al., 2011). Understanding and quantifying the influence of AHR on urban climate and extreme temperature events is an important component of urban energy management (Wang et al., 2013).

    Previous studies have focused mainly on the weather and short-term climate effects of AHR, including daily and seasonal impacts. AHR is significantly positively correlated with the urban heat island effect (Basara et al., 2010;Du et al., 2016) and raises temperatures by about 0.5°C-1°C in urbanized areas, an effect that is more pronounced in winter than in other seasons (Wilby, 2003; Feng et al.,2014; Wang et al., 2015). The influence of AHR is greater at night than in the daytime (Narumi et al., 2009; Mohan et al., 2012; Schrijvers et al., 2015), producing a greater increase in temperature, and the contribution of AHR among all warming factors increases to 75% at 0800 LST(Wang et al., 2013). For a specific intense heat wave event,AHR plays a key role in the synergies between the urban heat island effect and heat waves (Founda and Santamouris,2017; Zhao et al., 2018; Rizvi et al., 2019). According to studies of extreme temperature events, the changes in extreme temperature events are mostly related to greenhouse gas concentrations, anthropogenic aerosols, land-use change, and so on (Wen et al., 2013), but there has been a lack of consideration of AHR from energy consumption despite it being a key human activity in urban areas. Besides, extreme temperature events induced by AHR change the long-term frequency and trend of extreme temperature events, a phenomenon that deserves more attention but is rarely considered by meteorologists. Although AHR contributes only about 0.3% of the total heat flux at a global scale (Basara et al., 2010), its effect is increasing and very large in specific regions. Taking London and Tokyo as examples, the values of AHR have been estimated to be 10.9-1590 W m(Ichinose et al., 1999; Iamarino et al., 2012). In Beijing, rapid changes in AHR have altered the surface energy balance much more than the Chinese average (Chen et al., 2012; Yu et al., 2014). Therefore, it is appropriate to regard Beijing as a case for studying the influence of AHR on extreme temperature events, but there are some difficulties in using observational or model data based on case studies. Existing AHR data are either too low in resolution or too short-term (Chen et al., 2012; Dong et al., 2017). Determining the temporal and spatial variation of AHR in line with regional economic development and energy consumption is essential to reveal the effect on long-term climate and extreme temperature events. The Weather Research and Forecasting (WRF)model, which is based on fixed diurnal profiles of anthropogenic heat, is a good tool for research on climate effects but lacks a dynamic representation scheme for land-atmosphere interaction and does not consider abrupt temperature changes due to the effect of human activities (Dong et al.,2017).

    Recently, geoengineering strategies like increasing the albedo of urban roofs, implementation of green infrastructure, and using highly reflective paint for road pavements(Yaghoobian and Kleissl, 2012; Wang et al., 2013; Norton et al., 2015) have been regularly shown to be effective ways to cool cities, but may not help to warm cities when extreme cold events occur. Engineers focus on these operational methods to reduce heat or cold stress risks, but give less consideration to the influence of energy consumption on climate.Anthropogenic heat flux is closely related to energy consumption, which is the main factor in urban heat island intensity and heat wave events (Ryu and Baik, 2012). Energy consumption adjustments can provide another way to deal with extreme temperature events from the perspective of AHR.

    This work aims to reveal the long-term frequency and trend of extreme temperature events caused by AHR from energy consumption in a case study of Beijing during the 21st century, and reveal the influences of AHR on extreme temperature events over long periods, which might provide advice to urban planners trying to mitigate extreme temperature events. Then, we discuss whether the study methods can be applied to other cities. Here, temporal and spatial AHR data for Beijing from 1999 to 2017 were gathered for the case study, and the corresponding modifications to the AHR scheme in the WRF urban canopy model (UCM) were made to start the simulations.

    The rest of this paper is structured as follows: Section 2 describes the model development, including a new AHR scheme considering the influence of abrupt temperature changes on anthropogenic heat flux. Section 3 discusses the data and the experimental design. Section 4 presents the effects of AHR on extreme temperature events, including AHR variation, model validation, the frequency and trend of extreme heat events, and the mechanism relating AHR to extreme temperature events. Concluding remarks are given in section 5.

    2. Model development

    As part of the sensible heat flux, AHR upsets the raw energy balance. Equations (1) and (2) represent the energy balance equation without and with AHR, respectively. Once AHR increases are incorporated into the energy balance, the net radiation flux transferred from the urban land surface to the atmosphere also increases:

    and

    where Rn, LE

    ,

    H

    ,

    G

    and AHR represent the net radiation flux, latent heat flux, sensible and ground heat flux, and anthropogenic heat flux from the urban land surface to the near-surface atmosphere. The subscripts “no-AHR” and“AHR” represent the energy balance without and with AHR, respectively. In the default scheme of the WRF model, the AHR (denoted as AHR) is the product of the maximum AHR (A HR) and the empirical diurnal profile [

    f

    (

    h

    )] without spatial differences. AHR in this paper(A HR) is a spatiotemporal variable that changes according to latitude and longitude, as well as yearly, monthly and daily. The climatological and abrupt temperature changes are also considered, showing the increasing release of AHR when the weather is very cold or very hot:

    The WRF model is a limited-area, nonhydrostatic, mesoscale modeling system with a terrain-following eta coordinate, coupled with a UCM to provide a better representation of the physical processes involved in the exchange of heat,momentum, and water vapor in an urban environment. The UCM is a single-layer model with a simplified urban geometry. Some of its features include shadowing from buildings, reflection of short- and longwave radiation, the wind profile in the canopy layer, and multi-layer heat transfer equations for roofs, walls, and road surfaces (Tewari et al.,2007). Anthropogenic heat is optional in urban environmental simulations, which are based on fixed AHR values according to urban types. For example, the daily maximum AHR values are 90, 50 and 20 W mfor high-, middle- and low-density urban areas in the default AHR scheme. An empirical diurnal profile is used to represent the daily variation in anthropogenic heat flux in the single-layer model.The advanced AHR scheme mentioned above was coupled into the WRF model (see Fig. 1). AHR is regarded as part of the sensible heat flux and changes the energy interaction by upsetting the energy balance between the land surface and the atmosphere. In addition, climatological temperature and AHR gridded data were input into this scheme by a data interface of the WRF model when conducting initialization of real-time case, and the current temperature was from each step of the model output.

    Fig. 1. Schematic diagram of the anthropogenic heat scheme coupled with the WRF model.

    3. Data and experiments

    3.1. AHR

    AHR can be divided into four components (AHR from vehicles, the building sector, industry, and human metabolism) (Sailor and Lu, 2004), representing the major sources of waste heat in the urban environment, as shown in Eq. (7):

    where

    Q

    ,

    Q

    ,

    Q

    and

    Q

    represent the AHR from vehicles(V), the building sector (B), industry (I), and human metabolism (M). Human metabolic rates are not constant over the course of a day, and normally metabolic heat flux accounts for less than 1 W mof the total anthropogenic heat and can therefore be ignored (Grimmond, 1992; Sailor, 2011).Hence, in this paper, metabolic heat flux is not included in the total AHR. Because anthropogenic heat is closely tied to energy consumption (Sailor, 2011), the AHR from vehicles,buildings, and industry can be estimated by listing all kinds of energy consumption. Hence, AHR can be estimated by the following function:

    where

    E

    represents the

    k

    th kind of energy consumption(with the subscripts denoting coal, coke, crude oil, gasoline,kerosene, diesel, fuel oil, natural gas, and electric power)and

    a

    is the heat release coefficient of the

    k

    th kind of energy consumption. Hence, AHR can be represented by nine kinds of energy consumption. However, the energy consumption of Beijing is recorded by the National Bureau of Statistics and released yearly in a form that is incompatible with the grid resolution of the model simulation. Hence, the energy consumption must be downscaled to the grid size. Population density has been commonly used to map total energy consumption to different spatial scales (Washington,1972; Block et al., 2004), including the energy consumption of Beijing for nine kinds of energy from 1999 to 2017.Empirical monthly profiles of Beijing’s AHR were extracted from a global anthropogenic heat flux database with a high spatial resolution of 30 arc-seconds and a temporal resolution of 1 h (Dong et al., 2017), which was based on gridded monthly temperatures derived from various measurements and radiance-calibrated nighttime lights. The monthly profiles showed spatial variations and were convincing, but the time period of the raw data was limited to one year,which could not meet the research needs of this paper. The calculation was performed as in Eq. (9), and a schematic of the AHR methodology is shown in Fig. 2:

    Fig. 2. Flowchart of AHR data estimation.

    3.2. Data

    Air temperature and wind speed data for Beijing from observation and reanalysis datasets were used for model validation (Table 1). In this paper, observed temperature data from 20 national meteorological stations in Beijing and wind speeds from 205 non-watching automatic weather stations were used. The regional reanalysis data came from the China Meteorological Administration Land Data Assimilation System (CLDAS) as hourly outputs with a resolution of 0.0625° × 0.0625°. Regional and site observations of temperature comparisons were carried out to verify the effectiveness of the WRF physical schemes and determine whether these schemes were reasonable for long-term urban simulation with the new AHR scheme. This paper mainly focuses on comparing AHR simulation results with reanalysis data and in-situ data.

    Table 1. Dataset descriptions.

    3.3. Experimental design in the case study

    The Advanced Research version of the WRF model(Skamarock et al., 2008) dynamic core, version 3.9.1,coupled with a single-layer UCM was used in the experiments in this study. Two kinds of experiments were performed, with and without AHR (denoted “AHR” and “no-AHR”, respectively, in this paper). In the AHR experiment,the AHR data were regarded as a new input variable by adding them to the model initial input files with the same spatial resolution as the model and replacing the raw files once the model was running. Other variables, like many-year averaged air temperature and monthly profiles of the spatial distribution, were input from the initial files to calculate the hourly AHR for a specific year and day. AHR was treated as part of the sensible heat flux in this model, transferring to the near-surface atmosphere and upsetting the normal energy balance.

    In this study, two-layer nested domains with horizontal resolutions of 15 km (d01; mesh size 121 (lon) × 144 (lat),most of northern China), 3 km (d02; mesh size 206 (lon) ×251 (lat); almost all the Jing-Jin-Ji metropolitan area, and Beijing as the area of interest), were designed for these experiments (Fig. 3). NCEP FNL 6-h data (including soil water,moisture, and temperature) were used for the first-guess initial field and lateral boundary conditions. SSTs were needed to run long-term simulations, for which the OISST high-resolution dataset (daily, 0.25° × 0.25°) was chosen in this study.The terrestrial and MODIS land-use data provided in the WRF model described the real terrestrial and land-cover characteristics of the regions of interest, and the default static data were used in the experiments. The simulation period ran from 1 December 1999 to 31 December 2017. The first year, from December 1999 to December 2000, was considered as an initialization period. Two-meter air temperature data from 2001 to 2017 were used for further analysis of the model outputs for each hour. In terms of physical options, the schemes described in Table 2 were used.

    Fig. 3. Simulation area with two domains. The left-hand panel shows that D01 covers North China mostly,and the right-hand panel shows the position of D02 nested in D01, the color bar represents topography height.

    4. Effects of AHR on extreme temperature events

    4.1. Model validation

    Compared with the seven-year (from 2008 to 2014) average temperature in the CLDAS and AHR simulation results shown in Fig. 4, both temperatures in the city center were higher than in the suburbs. The similar spatial distributions show that the WRF physical scheme is reasonable. The correlation coefficients between the CLDAS temperatures and the AHR simulation results were mostly quite close to one, and the average root-mean-square error (RMSE) of the two datasets was 0.8°C. The gridded model results were interpolated to the site according to the coordinates of stations in both urban and suburban areas, and comparisons between site observations and AHR simulation results also showed that the simulation results for temperature were in good agreement with observations (Fig. 4b). Table 3 gives the regression coefficients, RMSE, and square of the correlation.There were nine stations on six ring roads and 11 stations in the suburbs. The results show that the AHR simulation results are reasonable, with simulated temperatures slightly higher than observations and RMSEs less than 2.5°C. Also,10-m wind speed comparisons were conducted (see Fig. 5).Compared with CLDAS analysis data, and 205 non-watching automatic weather stations, the correlation coefficients are mostly larger than 0.6, and the RMSE is less than 1.3 m s,showing the AH simulation had good agreement with both CLDAS data and observations.

    Table 2. Physical parameterization schemes.

    4.2. Frequency and trend of extreme temperature events in the case study

    In the approach proposed here, daily average AHR i s mainly based on energy consumption and population density, which both have spatial distributions. Most AHR is concentrated within the six ring roads of Beijing (see Fig. 6a),with maximum values greater than 80 W min the city center; in the suburbs, the AHR is less than 2 W m. In yearly terms, from Fig. 6b, due to Beijing’s energy management for environmental protection, high-energy industries have moved to other provinces or have moved away from using dirtier forms of energy, such as coal and coke. Instead, environmentally friendly energy such as natural gas and electric power, along with irreplaceable energy like gasoline, are increasing rapidly. Energy upgrades are mostly conducted in the city center, with impacts on urban central heating in winter, traffic on roads, and some factories. Although dirty or clean energy consumption may increase or decrease with urban regulation, the anthropogenic heat of Beijing increased from 2001 to 2010, but after 2010, AHR became stable and reached a maximum. From the monthly variations shown in Fig. 6c, which are based on the work of Dong et al. (2017), it is apparent that AHR is higher in winter than in other seasons, which makes sense because urban central heating is used in winter to counteract low temperatures, leading to more energy consumption. Spatial differences in AHR are much greater than monthly variations, as can be seen from the bar value of each month representing the logof spatial variations in Fig. 6c. In winter, the variations in spatial distribution are bigger than in other seasons. The estimated diurnal profiles are different due to changes in daily lifestyle pattern such as holidays and festivals. Diurnal profiles in this research were based on the default profiles, but took abrupt temperature changes into consideration according to Eq. (6), which was applied to the UCM.

    Fig. 4. Comparison between observed and simulated values: (a) correlation coefficients of yearly temperatures between CLDAS and AHR simulation; (b) RMSE and correlation coefficients of temperature between observation and AHR simulation. The observation stations can be seen in Table 3.

    Table 3. Regression coefficients (RCs), RMSE, and square of the correlation (R2) of observed and AHR simulated temperatures.

    Fig. 5. Wind speed comparisons between observed and simulated values: (a) spatial distribution of temporal correlation coefficients between CLDAS reanalysis data and AHR simulations; (b) correlation coefficients and RMSEs between station observations and AHR simulations, where the sizes of circles denote the values of correlation coefficients and the colors represent the RMSE.

    Fig. 6. (a) Spatial distributions of daily average AHR from 2001 to 2017. The closed curve from the inside to the outside represents the second to sixth ring roads of Beijing with high population density. (b) Yearly and (c) monthly profiles of AHR, with the bars representing the spatial standard deviation of lgAHR.

    Comparing the AHR and no-AHR simulation results,the temperature differences were mainly located in the city center, with a maximum increase of 1.14°C, and there was variation among the seasons, with average increases of 0.06°C, 0.01°C, 0.11°C and 0.18°C for spring, summer, fall and winter, respectively. The frequency distribution of the temperature increase (Figs. 7f-j) shows that in suburban areas the temperature showed no increase, but in the city center it increased more because more than 90% of AHR is located within the six ring roads of Beijing, where the area of high population density is situated. There, anthropogenic heat flux release to the near-surface atmosphere heats the air directly and makes the temperature rise. From Figs. 7a-e,both frequency distributions of temperature increase show that the city center experiences a stronger effect than the suburbs. The yearly and four seasonal average temperature increases were 0.525°C, 0.362°C, 0.191°C, 0.621°C and 0.936°C, respectively, with seasonal temperature increases similar to the monthly AHR profiles. However, the heating efficiency of AHR in summer is less significant than in other seasons [with the lowest Δ

    T

    /AHR in Table 4, where the temperature increase per unit of AHR (Δ

    T

    /AHR) is referred to as the heating efficiency]. The results of this study suggest that atmospheric dynamic processes in different seasons may transport more or less heat to the upper atmosphere. Because summer is more unstable than other seasons, a large amount of anthropogenic heat is transferred from the low atmosphere to upper layers, making temperature rises in summer not correspond exactly with a given amount of AHR. This will be confirmed in section 4.3.

    Fig. 7. Yearly and seasonal temperature rises due to the spatial distribution and frequency distribution of AHR: (a)yearly averaged temperature rises due to the spatial distribution of AHR; (b-e) as in (a) but for summer, autumn and winter; (f) yearly temperature rises due to the frequency distribution of AHR; (g-j) as in (f) but for spring, summer autumn and winter.

    Table 4. Average temperature rises and AHR in the urban area.

    Extreme climate and weather events are generally multifaceted phenomena. The present study pays particular attention to discussing climate extremes based on daily temperature without a spatial distribution (Sillmann et al., 2013),such as the hottest or coldest day of the year. Gridded datasets of indices (Alexander et al., 2006) have been developed to represent the spatial distribution of extreme heat events.Seven extreme temperature indices (ETIs) can be divided into two groups: extreme cold indices (ECIs) and extreme heat indices (EHIs). The ECIs include frost days (FD), cold nights (TN10p), and cold days (TX10p). EHIs include summer days (SD), warm nights (TN90p), warm days (TX90p),and diurnal temperature range (DTR), which is a measure of the difference between extreme cold events and extreme heat events. Here, these indices are introduced to illustrate the frequency and trend of extreme heat events, although these indices are commonly applied to climate periods or much longer periods with reoccurrence times of more than 20 years. Hourly temperature data from simulations were used to calculate the seven indices; the indices are described in Table 5.

    The ETI differences between the AHR and no-AHR experiments (noted as ΔTX10p, ΔTX90p, ΔTN10p,ΔTN90p, ΔFD, ΔSD, and ΔDTR), which eliminate the effect of climate change, can be used to show the spatial distribution of extreme temperature events with the effect of AHR. Figures 8a-g show that extreme cold events and heat events in the city center were stronger than in the suburbs due to greater AHR in the city center, similar to the effect of the temperature rise. Beijing’s six ring roads form a dividing line between significant and insignificant effects.

    In addition, a linear regression between time series and extreme event index differences can be used to quantify the temporal trend of extreme temperature events with the effect of AHR. The yearly ETI difference regression coefficients for the city center and the suburbs are given in Table 6,and the extreme temperature event indices of the time series and the temporal trend line can be seen in Figs. 8o-u. These results show that extreme cold events are decreasing, while extreme heat events are increasing, in the city center, with the regression coefficients of extreme cold event indices(ΔTX10p, ΔTN10p, and ΔFD) being negative and those of extreme heat event indices (ΔTX90p, ΔTN90p, and ΔSD)being positive from 2001 to 2017. However, the regression coefficients of extreme event indices in the suburbs were close to 0 and showed few tendencies. Figures 8o-u (blue line) show that extreme events were not significantly changing in the suburbs. Comparisons of regression coefficients between urban and suburban areas illustrate that the increase in the extreme heat event signal in urban areas is much stronger than in suburban areas. AHR changed the trend of extreme temperature events (Figs. 8h-n): cold events were annually decreasing by 0.26 days for cold days(ΔTX10p), 0.56 days for cold nights (ΔTN10p), and 0.39 days for frost days (ΔFD). The changes in extreme cold events and extreme heat events were not identical. For extreme heat events, the indices increased by 0.11 days for warm days (ΔTX90p), 0.02 days for warm nights(ΔTN90p), and 0.19 days for summer days (ΔSD) annually,or less strongly than for extreme cold events. Furthermore,the decrease in diurnal temperature range (ΔDTR) shows that the tendencies of extreme heat events and extreme cold events are not identical, with extreme cold events more severe and more frequent than extreme heat events.Through its heating effect, AHR makes hot days even hotter, but mitigates occurrences of cold weather.

    Table 5. Descriptions of extreme temperature event indices.

    Because AHR and ETIs are of different orders of magnitude, both AHR and ΔETIs were normalized, and the correlation coefficients of urban and suburban areas between the normalized AHR and the seven ETIs were calculated to find out the relationship between yearly AHR and extreme events. From Table 6, the absolute values of the correlation coefficients in urban areas are mostly greater than in suburban areas. ΔTX10p, ΔTN10p and ΔFD are negative, and ΔTX90p, ΔTN90p and ΔSD are positive, showing that the relationships between AHR and extreme cold events are different from those between AHR and extreme heat events,with AHR increasing from 2001 to 2017. This relationship ismuch stronger in urban than in suburban areas because the correlation coefficients are close to 0 in suburban areas.

    Fig. 8. Temporal and spatial distributions of ΔETIs: (a-g) spatial distribution of average ΔETIs from 2001 to 2017; (h-n)spatial distributions of the regression coefficients of temporal ΔETIs; (o-u) ΔETI time series and temporal trend lines from 2001 to 2017 in urban areas (red dotted line) and suburban areas (blue dotted line); solid red and blue lines are temporal trend lines.

    The results shown in Table 6 are corroborated by the spatial distribution of temporal correlation coefficients between AHR and the ETIs and the

    t

    -test statistical significance of the ETIs between the no-AHR and AHR experiments (Fig. 9).Both the

    t

    -test significance results and the spatial coefficient coefficients show that AHR has a great impact on extreme temperature events, although some suburban areas had larger correlation coefficients, but did not pass the statistical significance test and are considered spurious signals.Comparing the shaded areas, the region affected by extreme cold events is slightly larger than that affected by extreme heat events. It is also clear that those areas that passed the significance test also had a relatively high correlation in urban areas.

    Table 6. Regression coefficients of yearly ETIs and correlation coefficients between normalized AHR and ΔETIs.

    Fig. 9. Spatial distribution of temporal correlation coefficients between AHR and ΔETIs; the ETIs from (a) to (g) represent TX10p, TX90p, TN10p, TN90p, FD, SD and DTR. The shaded areas, which passed the t-test at a confidence level of 90%,indicate that the differences in extreme temperature events were significant when releasing anthropogenic heat flux.

    4.3. Mechanisms of extreme temperature events with AHR

    Fig. 10. The (a-c) PBLH, (d-f) specific humidity at 2 m (Q2), and (g-i) spatial distribution of the velocity of 10-m (w10)winds, from the (a, d, g) AHR and (b, e, h) no-AHR simulations, and (c, f, i) their differences. The shaded areas passed the ttest at a confidence level of 90%.

    Increases in AHR had an impact on dynamic processes in the atmospheric boundary layer. The planetary boundary layer height (PBLH) increased with rising AHR (Figs. 10a-c),with this change in PBLH amounting to more than a 100-m increase in the city center. The convergence of 10-m surface wind was also strengthened (Fig. 10g-i), with air converging in the city center and moving to the higher atmosphere,which also strengthens the vertical exchange of energy.Moreover, higher PBLH and stronger winds made the nearsurface atmosphere less stable. With stronger vertical transportation, heat transfers from the near-surface to the higher atmosphere are easier, and this will cool the near-surface temperature, which can be called the cooling effect. AHR is one of the biggest influencing factors on urban heat islands, and some researchers have already found a synergistic effect between urban heat islands and heat waves (Li and Bou-Zeid, 2013). Hence, AHR may have a relationship with extreme heat events, when the high temperature reaches some specified level. Some researchers have discovered that enhanced evaporation and lowering of near-surface humidity can lead to more heat convergence, making extreme heat events happen more easily (Freychet et al., 2017). However,in the results of this study, specific humidity had little change and did not pass the

    t

    -test at a confidence level of 90% (Figs. 10d-f), showing that AHR influences the extreme temperature events through heating directly.Extreme temperature events in this study were found to be controlled by the heating effect of AHR, but can be relieved slightly due to the cooling effect of the boundary layer dynamic process. Here, -ΔPBLH is used to represent the dynamic process, assuming that the dynamic process strength varies negatively with temperature and positively with PBLH. A multiple regression of ΔETI, AHR, and-ΔPBLH was performed; all the data were normalized firstly to remove the effect of the order of magnitude. As shown in Table 7, the boundary layer dynamic process depressed the number of occurrences of extreme temperature events and reduced the impact of AHR on these events.Extreme heat events are more frequent with greater AHR,but the boundary dynamic process damps down the tendency toward an increase. For cold extreme events, for example, the regression coefficients of anthropogenic heat for TX10p and TN10p were -0.081 and -0.319, and those of vertical transportation were 0.521 and 0.572. Extreme cold events are less frequent with greater AHR, but the boundary layer dynamic process damps down the tendency toward a decrease. Extreme heat events have the opposite tendency,being more frequent with greater AHR. Heating is the main effect of AHR, by increasing sensible heat flux and raising urban temperatures. It shows seasonal variation in heating efficiency with temperature, with winter temperatures increasing more than summer temperatures for each unit of AHR.Multiple linear regression of Δ

    T

    , AHR, and -ΔPBLH was performed in the urban area to reveal the direct heating effect of AHR and the cooling effect of atmospheric dynamic processes in different seasons (shown in Table 8). The regression coefficients of AHR from spring to winter are 0.21,-0.258, 0.128 and 0.381. The results show that winter has a larger regression coefficient of AHR than the other seasons and that the heating effect of AHR in winter is about 1.6 times stronger than over the whole year for the temperature rise. In addition, the negative regression coefficient in summer reveals that AHR may have an opposite effect to the trend towards a temperature rise. The regression coefficients of -ΔPBLH from spring to winter are -0.81, 0.90,-0.88 and -0.42. From comparisons of the absolute regression coefficients of -ΔPBLH over the four seasons, it was concluded that the impact of the boundary layer atmospheric dynamic process on the temperature rise, with the strongest effect in summer and the weakest in winter, is quite opposite to that of AHR. Note also that the absolute regression coefficients of -ΔPBLH were always larger than those of AHR, which means that -ΔPBLH is more related to trends in ΔT, but it does not necessarily make a larger contribution than AHR.

    Table 7. Multiple linear regression of ΔETI, AHR and -ΔPBLH.

    Table 8. Multiple linear regression of ΔT, AHR and ΔPBLH.

    5. Conclusions and discussion

    This paper has revealed the relationship between AHR from energy consumption and extreme temperature events in a case study of Beijing, China. Temporal and spatial variation data for AHR were processed, which were based on energy consumption, population density, and previous research. AHR increased to a steady state with declining use of coal and more widespread use of clean energy from 2001 to 2017. AHR is mainly concentrated within the six ring roads of Beijing, with a higher monthly profile in winter.Besides, its spatial variation was much larger than its variation over time. A new AHR scheme incorporating abrupt temperature changes was implemented in the WRF model.Two experiments, with and without AHR, were performed to simulate the effect of AHR on climate in the case study.Temperature rise due to the AHR heating effect was mainly concentrated in the urban center. In spring, summer, fall and winter, the temperature increased by 0.36°C, 0.19°C,0.62°C and 0.93°C, respectively, with similar seasonal variations in AHR. However, the heating efficiencies were different, with the smallest temperature increase per unit AHR in summer and the largest increase in winter. Both extreme heat events and extreme cold events could be represented using seven extreme temperature event indices. AHR reinforced the occurrence of extreme heat events, but restrained the development of extreme cold events, with an increase in extreme heat events of 0.11 days for warm days (ΔTX90p),0.02 days for warm nights (ΔTN90p), and 0.19 days for summer days (ΔSD), and an annual decrease in extreme cold events of 0.26 days for cold days (ΔTX10p), 0.56 days for cold nights (ΔTN10p), and 0.39 days for frost days (ΔFD)in the city center. In addition, the diurnal temperature range(ΔDTR) was narrower. It shows that the decrease in the frequency of annual extreme cold events is larger than the increase in the frequency of annual extreme heat events.The reason might be that coal and oil consumption release more anthropogenic heat flux in winter and weaken the effect of vertical transportation, which decreases the frequency of extreme cold events. Besides, the correlation coefficients between AHR and ETI differences (ΔETIs) and the regression coefficients both show the strong signal of this trend in the city center.

    Note also that AHR and the atmospheric dynamic process play opposite roles in urban temperature rises and the occurrence of extreme temperature events. The dynamic process restrains the increasing trend in extreme heat events and helps increase the number of extreme cold events. In addition, the effect of lower heating efficiency in summer than in other seasons could be making a contribution, with strong convection transferring heat to the higher atmosphere and reducing the warming ability.

    A possible limitation is that ignoring the typical routes taken by inhabitants, as well as the distributions of workplaces and thermal plants, will result in AHR estimation errors. However, larger regional analysis will weaken the influence of the uncertainty of AHR, and that is why we divided Beijing into three parts to analyze the relationships between AHR and extreme temperature events. Also, the results were obtained from a case study, but the study method can still be applied to other cities with different simulation areas. From our research, AHR heated the atmosphere which changed the frequencies and trend of extreme temperature events. According to the global high-resolution anthropogenic heat release data, cities in northern mid-latitudes(such as Dalian, Shanghai, Guangzhou, Taiwan and cities of Yangtze River Delta in East China, Seoul of South Korea,Tokyo of Japan, lots of cities in Europe like Amsterdam, London and Lodz, San Francisco, New York and Washington DC in United States) which release more anthropogenic heat shall have more significant influence than other places on extreme temperature events. Developing countries like China, India and Egypt have a stronger increasing tendency of extreme temperature events due to higher energy consumption growth. However, the results may be changed with different energy management practices, as well as dynamic processes in the atmospheric boundary layer.

    Finally, it is foreseeable that decreasing AHR by reducing energy consumption could be a way to mitigate extreme heat events, but might not be helpful in dealing with extreme cold events. For Beijing, it is reasonable to reduce energy consumption to deal with extreme heat events because of hot days in summer; meanwhile, energy consumption could be adjusted to keep the city warm and decrease extreme cold events in winter.

    . This work was supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDA23090102), the National Natural Science Foundation of China (Grant No. 41830967), the Key Research Program of Frontier Sciences, Chinese Academy of Sciences (Grant No. QYZDY-SSW-DQC012), and the National Key Research and Development Program of China (Grant Nos. 2018YFC1506602 and 2020YFA0608203). We also thank the National Meteorological Information Center, China Meteorological Administration,for data support.

    亚洲av综合色区一区| 寂寞人妻少妇视频99o| 久久精品国产自在天天线| av天堂久久9| 日韩免费高清中文字幕av| 婷婷成人精品国产| 在线观看免费视频网站a站| 熟女电影av网| 国产精品无大码| 国产亚洲欧美精品永久| 最后的刺客免费高清国语| 纯流量卡能插随身wifi吗| 91精品一卡2卡3卡4卡| 人妻夜夜爽99麻豆av| 国产伦精品一区二区三区视频9| 一级二级三级毛片免费看| 丝袜脚勾引网站| 国产在线一区二区三区精| 国产成人freesex在线| 国产精品女同一区二区软件| 老司机亚洲免费影院| 中文字幕人妻熟人妻熟丝袜美| 欧美激情极品国产一区二区三区 | 又粗又硬又长又爽又黄的视频| 久久午夜福利片| 亚洲成人一二三区av| 国产成人午夜福利电影在线观看| 国产精品久久久久久久电影| 少妇人妻 视频| 国产黄片视频在线免费观看| 青春草国产在线视频| 在线观看免费视频网站a站| 亚洲av日韩在线播放| av免费在线看不卡| 美女福利国产在线| 丝袜在线中文字幕| 欧美激情 高清一区二区三区| 22中文网久久字幕| 中文精品一卡2卡3卡4更新| 成人无遮挡网站| 九九爱精品视频在线观看| 国产一区二区三区综合在线观看 | 国产精品女同一区二区软件| 久久精品国产亚洲网站| 午夜视频国产福利| 纵有疾风起免费观看全集完整版| 中国三级夫妇交换| av国产久精品久网站免费入址| 在线亚洲精品国产二区图片欧美 | 婷婷色综合www| 亚洲av在线观看美女高潮| 国产在线免费精品| 69精品国产乱码久久久| 中国国产av一级| 欧美xxxx性猛交bbbb| 成年女人在线观看亚洲视频| 亚洲第一区二区三区不卡| 久久久久久久久大av| 毛片一级片免费看久久久久| 纯流量卡能插随身wifi吗| 大香蕉久久网| av免费观看日本| 人人妻人人添人人爽欧美一区卜| 久久99热6这里只有精品| 午夜免费男女啪啪视频观看| 天天躁夜夜躁狠狠久久av| 有码 亚洲区| 男人爽女人下面视频在线观看| 天天影视国产精品| 狠狠婷婷综合久久久久久88av| 美女国产视频在线观看| 王馨瑶露胸无遮挡在线观看| 黑人巨大精品欧美一区二区蜜桃 | 亚洲国产av新网站| 国产精品久久久久久av不卡| 亚洲美女搞黄在线观看| 亚洲中文av在线| 少妇丰满av| 啦啦啦视频在线资源免费观看| 蜜桃国产av成人99| 伊人久久精品亚洲午夜| 日本色播在线视频| 久久精品久久精品一区二区三区| 嘟嘟电影网在线观看| 91精品三级在线观看| 在线观看免费日韩欧美大片 | 黄色一级大片看看| 制服诱惑二区| 男人操女人黄网站| 黄片播放在线免费| 91精品伊人久久大香线蕉| 亚洲成人一二三区av| 久久鲁丝午夜福利片| 少妇被粗大猛烈的视频| 啦啦啦视频在线资源免费观看| 亚洲图色成人| 欧美日韩综合久久久久久| 亚洲国产精品专区欧美| 欧美日韩国产mv在线观看视频| 精品人妻一区二区三区麻豆| 在线天堂最新版资源| 97超视频在线观看视频| 亚洲av免费高清在线观看| 久久国产精品男人的天堂亚洲 | 精品久久国产蜜桃| 久久久久国产精品人妻一区二区| kizo精华| 免费观看a级毛片全部| 啦啦啦视频在线资源免费观看| 久久精品国产a三级三级三级| 99精国产麻豆久久婷婷| 十分钟在线观看高清视频www| 国产精品嫩草影院av在线观看| av专区在线播放| 亚洲精品乱码久久久久久按摩| 国产精品女同一区二区软件| 久久久久久久大尺度免费视频| 国产精品.久久久| 少妇人妻 视频| 日韩免费高清中文字幕av| 美女中出高潮动态图| 22中文网久久字幕| 9色porny在线观看| 亚洲五月色婷婷综合| 亚洲性久久影院| 日韩人妻高清精品专区| 伦理电影免费视频| 另类精品久久| 国产在线免费精品| 亚洲一级一片aⅴ在线观看| 曰老女人黄片| 熟女av电影| 欧美人与性动交α欧美精品济南到 | 亚洲精品一二三| 大香蕉久久成人网| 啦啦啦中文免费视频观看日本| 九九爱精品视频在线观看| 免费观看在线日韩| 亚洲欧洲日产国产| 国产无遮挡羞羞视频在线观看| 岛国毛片在线播放| 国产成人免费观看mmmm| 2021少妇久久久久久久久久久| 久久婷婷青草| 欧美精品人与动牲交sv欧美| 国产成人av激情在线播放 | 婷婷色综合大香蕉| 国产成人a∨麻豆精品| 新久久久久国产一级毛片| 三级国产精品欧美在线观看| 精品一区在线观看国产| 婷婷色av中文字幕| 国产成人精品在线电影| 精品午夜福利在线看| 狂野欧美白嫩少妇大欣赏| 国产色婷婷99| 亚洲av中文av极速乱| 久久精品久久精品一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 欧美精品人与动牲交sv欧美| 日韩成人av中文字幕在线观看| 精品亚洲乱码少妇综合久久| 寂寞人妻少妇视频99o| 中国国产av一级| 内地一区二区视频在线| 亚洲国产精品一区二区三区在线| av一本久久久久| 国产乱来视频区| 黄色一级大片看看| 亚洲精品色激情综合| 久久精品夜色国产| 国产在线视频一区二区| 亚洲一级一片aⅴ在线观看| 少妇 在线观看| 亚洲精品美女久久av网站| 日本vs欧美在线观看视频| 午夜激情av网站| 欧美激情国产日韩精品一区| 我的女老师完整版在线观看| 亚洲精品第二区| 欧美激情极品国产一区二区三区 | 一本一本综合久久| 日韩在线高清观看一区二区三区| 亚洲欧洲国产日韩| 热99国产精品久久久久久7| 国产不卡av网站在线观看| 亚洲欧洲国产日韩| 人妻少妇偷人精品九色| 不卡视频在线观看欧美| 桃花免费在线播放| 三级国产精品片| 久久久国产精品麻豆| 欧美精品国产亚洲| 热re99久久精品国产66热6| 免费高清在线观看日韩| 亚洲精品国产色婷婷电影| 在线亚洲精品国产二区图片欧美 | 亚洲av综合色区一区| 国产精品久久久久久久电影| 欧美性感艳星| 高清在线视频一区二区三区| 成人毛片a级毛片在线播放| 人人妻人人添人人爽欧美一区卜| 国产色爽女视频免费观看| 免费高清在线观看视频在线观看| 精品卡一卡二卡四卡免费| 80岁老熟妇乱子伦牲交| 搡女人真爽免费视频火全软件| 欧美日韩av久久| 国产在线视频一区二区| 久久人人爽人人爽人人片va| 欧美丝袜亚洲另类| 国产有黄有色有爽视频| 久久ye,这里只有精品| 蜜桃久久精品国产亚洲av| 久久久久久久久久久久大奶| 久久97久久精品| 成年女人在线观看亚洲视频| 99热国产这里只有精品6| 下体分泌物呈黄色| 久久鲁丝午夜福利片| 日韩av免费高清视频| 晚上一个人看的免费电影| 亚洲欧美精品自产自拍| 午夜免费鲁丝| 国产成人精品一,二区| 91国产中文字幕| 国产在视频线精品| 少妇被粗大的猛进出69影院 | 日本午夜av视频| 观看美女的网站| 在线观看一区二区三区激情| 美女脱内裤让男人舔精品视频| 中文精品一卡2卡3卡4更新| 亚洲国产成人一精品久久久| 国产精品 国内视频| 91精品国产国语对白视频| 精品一品国产午夜福利视频| 最新中文字幕久久久久| 水蜜桃什么品种好| 老司机影院毛片| 黑人欧美特级aaaaaa片| 又黄又爽又刺激的免费视频.| 国产 一区精品| 色婷婷久久久亚洲欧美| 成人免费观看视频高清| 亚洲欧洲精品一区二区精品久久久 | 国产一区二区三区综合在线观看 | 97精品久久久久久久久久精品| 亚洲成色77777| 青春草亚洲视频在线观看| 涩涩av久久男人的天堂| 制服诱惑二区| 国产精品一区www在线观看| .国产精品久久| 日韩强制内射视频| 国产亚洲av片在线观看秒播厂| 街头女战士在线观看网站| 国产精品久久久久成人av| 精品少妇黑人巨大在线播放| 高清av免费在线| 蜜桃久久精品国产亚洲av| 日产精品乱码卡一卡2卡三| 日韩免费高清中文字幕av| 日韩,欧美,国产一区二区三区| 国产国拍精品亚洲av在线观看| 三上悠亚av全集在线观看| 亚洲av国产av综合av卡| 人妻人人澡人人爽人人| 亚洲在久久综合| 日韩av免费高清视频| 伊人久久国产一区二区| 一区在线观看完整版| 最近2019中文字幕mv第一页| 国产av精品麻豆| 欧美日韩视频高清一区二区三区二| 青春草视频在线免费观看| 人人妻人人澡人人看| 麻豆乱淫一区二区| 成年美女黄网站色视频大全免费 | 国产免费福利视频在线观看| 嫩草影院入口| 久久精品国产亚洲av涩爱| 九色成人免费人妻av| 乱人伦中国视频| av播播在线观看一区| 国产片特级美女逼逼视频| 国产高清不卡午夜福利| 色视频在线一区二区三区| 80岁老熟妇乱子伦牲交| 黑人高潮一二区| 99热这里只有精品一区| 在线观看免费日韩欧美大片 | 五月伊人婷婷丁香| 国产精品国产三级国产专区5o| 我的老师免费观看完整版| 亚洲成人手机| 91精品国产九色| 大话2 男鬼变身卡| 能在线免费看毛片的网站| 免费大片黄手机在线观看| 国产视频首页在线观看| 国产免费又黄又爽又色| 女人精品久久久久毛片| 亚洲不卡免费看| 亚洲精品日韩av片在线观看| 欧美日韩av久久| 久久人妻熟女aⅴ| 美女福利国产在线| 日本黄色日本黄色录像| 国产视频首页在线观看| 精品一品国产午夜福利视频| 五月玫瑰六月丁香| 老女人水多毛片| 久久精品久久久久久噜噜老黄| a级毛色黄片| 久久热精品热| 日韩亚洲欧美综合| 中文字幕人妻丝袜制服| 精品少妇内射三级| 少妇高潮的动态图| 欧美日韩视频高清一区二区三区二| 国产av一区二区精品久久| 各种免费的搞黄视频| 高清毛片免费看| 亚洲精品日韩av片在线观看| 亚洲国产精品成人久久小说| 熟女电影av网| videossex国产| 女性生殖器流出的白浆| 国产精品女同一区二区软件| 在线播放无遮挡| 欧美另类一区| 亚洲精品久久成人aⅴ小说 | 亚洲不卡免费看| 国产一区二区在线观看日韩| 欧美丝袜亚洲另类| 超色免费av| 亚洲五月色婷婷综合| 22中文网久久字幕| 国产精品久久久久久精品电影小说| 国产av码专区亚洲av| 18禁在线无遮挡免费观看视频| 亚洲国产日韩一区二区| 久久久久久久精品精品| 中国美白少妇内射xxxbb| 男女国产视频网站| av视频免费观看在线观看| 欧美精品国产亚洲| 国产日韩欧美在线精品| 日韩三级伦理在线观看| 涩涩av久久男人的天堂| 七月丁香在线播放| 亚洲人成网站在线播| 男人操女人黄网站| 久热久热在线精品观看| 国产精品无大码| 久久精品久久久久久久性| 只有这里有精品99| 亚洲欧美日韩另类电影网站| 国产黄色免费在线视频| 中文乱码字字幕精品一区二区三区| 亚洲欧美精品自产自拍| 免费日韩欧美在线观看| 亚洲第一区二区三区不卡| 中文字幕精品免费在线观看视频 | 日日撸夜夜添| videos熟女内射| 亚洲伊人久久精品综合| 精品人妻一区二区三区麻豆| 我的女老师完整版在线观看| tube8黄色片| 桃花免费在线播放| 丝袜喷水一区| 亚洲av福利一区| 精品一区二区三卡| 午夜激情久久久久久久| 日韩视频在线欧美| 亚洲国产欧美在线一区| 午夜影院在线不卡| 一级片'在线观看视频| 日韩电影二区| 一边亲一边摸免费视频| 一区二区日韩欧美中文字幕 | 久久综合国产亚洲精品| 国产伦理片在线播放av一区| 久久ye,这里只有精品| 国产精品不卡视频一区二区| 久久久精品94久久精品| 午夜免费男女啪啪视频观看| 国产日韩一区二区三区精品不卡 | 国产成人精品婷婷| av专区在线播放| 男女高潮啪啪啪动态图| 国产精品国产av在线观看| 久久久久人妻精品一区果冻| 国产精品不卡视频一区二区| 天堂8中文在线网| 国产在线视频一区二区| 卡戴珊不雅视频在线播放| 少妇高潮的动态图| 黄色毛片三级朝国网站| 国产高清不卡午夜福利| 久久国内精品自在自线图片| 永久免费av网站大全| 熟女人妻精品中文字幕| 蜜臀久久99精品久久宅男| 少妇熟女欧美另类| 久久99蜜桃精品久久| 国产精品久久久久久久久免| 少妇人妻 视频| av在线app专区| 我的女老师完整版在线观看| 欧美97在线视频| 精品久久久久久久久av| 赤兔流量卡办理| 9色porny在线观看| 国产在线免费精品| 日韩欧美精品免费久久| 夫妻性生交免费视频一级片| 丰满少妇做爰视频| 国产熟女欧美一区二区| 亚洲国产色片| 亚洲精品视频女| 亚洲av综合色区一区| 精品国产乱码久久久久久小说| 久久精品国产a三级三级三级| 中国三级夫妇交换| 久久婷婷青草| 99九九线精品视频在线观看视频| 人妻制服诱惑在线中文字幕| av国产久精品久网站免费入址| 日本-黄色视频高清免费观看| 日韩亚洲欧美综合| 视频在线观看一区二区三区| 国产精品久久久久久精品古装| 少妇丰满av| 涩涩av久久男人的天堂| 少妇的逼好多水| 在线播放无遮挡| 亚洲色图综合在线观看| 国产精品一区二区在线观看99| 春色校园在线视频观看| 久久久亚洲精品成人影院| 韩国av在线不卡| 五月伊人婷婷丁香| 精品亚洲成a人片在线观看| 爱豆传媒免费全集在线观看| 午夜视频国产福利| 天美传媒精品一区二区| 老司机影院毛片| 五月天丁香电影| 国产欧美另类精品又又久久亚洲欧美| 午夜影院在线不卡| 精品一区在线观看国产| www.av在线官网国产| 国产免费一区二区三区四区乱码| 国产成人精品无人区| 一级黄片播放器| 五月伊人婷婷丁香| 国产高清三级在线| 久久午夜综合久久蜜桃| 色94色欧美一区二区| 99久国产av精品国产电影| 亚洲精品自拍成人| 亚洲国产成人一精品久久久| 久久久久视频综合| 久久亚洲国产成人精品v| 91久久精品国产一区二区成人| 欧美日韩综合久久久久久| 久久久久精品性色| 欧美精品一区二区免费开放| 老熟女久久久| 国产精品三级大全| 看免费成人av毛片| 亚洲精品国产av成人精品| 欧美精品一区二区大全| 热re99久久精品国产66热6| 亚洲怡红院男人天堂| 国产成人精品一,二区| 亚洲精品国产色婷婷电影| 久久久久国产网址| 一本大道久久a久久精品| 亚洲精品视频女| 三上悠亚av全集在线观看| 久久久久视频综合| 国产毛片在线视频| 插阴视频在线观看视频| 18在线观看网站| 制服丝袜香蕉在线| 国产视频内射| 国产伦精品一区二区三区视频9| 大又大粗又爽又黄少妇毛片口| 9色porny在线观看| 最近中文字幕2019免费版| 亚洲国产精品成人久久小说| 日本午夜av视频| 一个人免费看片子| 亚洲情色 制服丝袜| 亚洲成人av在线免费| 亚洲人与动物交配视频| 一区二区三区精品91| 亚洲精品一区蜜桃| 夫妻午夜视频| 免费大片黄手机在线观看| 99久国产av精品国产电影| 亚洲精品,欧美精品| 国产乱人偷精品视频| 黑人猛操日本美女一级片| 一个人看视频在线观看www免费| 欧美日韩国产mv在线观看视频| 天天躁夜夜躁狠狠久久av| 特大巨黑吊av在线直播| 久久久久视频综合| 久久久久久久久久成人| 久久国内精品自在自线图片| 国产午夜精品久久久久久一区二区三区| 高清在线视频一区二区三区| 国产亚洲最大av| 国产精品一区二区三区四区免费观看| 亚洲av免费高清在线观看| 亚洲美女黄色视频免费看| 亚洲美女搞黄在线观看| 97在线视频观看| 女性生殖器流出的白浆| 国产精品一区www在线观看| 日韩中字成人| 夜夜爽夜夜爽视频| 黄色怎么调成土黄色| 午夜福利在线观看免费完整高清在| 一边摸一边做爽爽视频免费| 国产精品秋霞免费鲁丝片| 两个人免费观看高清视频| 美女福利国产在线| 亚洲精品视频女| av线在线观看网站| 成人影院久久| 亚洲成色77777| 午夜激情福利司机影院| 亚洲av成人精品一二三区| 久久久精品免费免费高清| 毛片一级片免费看久久久久| a级毛片在线看网站| 国产欧美日韩一区二区三区在线 | 亚洲精品美女久久av网站| 老司机影院成人| 午夜av观看不卡| 国产免费又黄又爽又色| 亚洲精品av麻豆狂野| 人妻夜夜爽99麻豆av| 在线亚洲精品国产二区图片欧美 | 日韩,欧美,国产一区二区三区| 日韩亚洲欧美综合| 看非洲黑人一级黄片| 大又大粗又爽又黄少妇毛片口| 国产av国产精品国产| av网站免费在线观看视频| 成年av动漫网址| a级毛片免费高清观看在线播放| 看十八女毛片水多多多| 亚洲五月色婷婷综合| 一级毛片 在线播放| 伦精品一区二区三区| 夜夜骑夜夜射夜夜干| 男人操女人黄网站| 一级毛片aaaaaa免费看小| 久久韩国三级中文字幕| 日韩人妻高清精品专区| 在线亚洲精品国产二区图片欧美 | 免费大片黄手机在线观看| 人妻少妇偷人精品九色| 菩萨蛮人人尽说江南好唐韦庄| 久久亚洲国产成人精品v| 高清欧美精品videossex| 街头女战士在线观看网站| 成人综合一区亚洲| 国产成人a∨麻豆精品| 汤姆久久久久久久影院中文字幕| videos熟女内射| 黄色怎么调成土黄色| 成人毛片a级毛片在线播放| 久久精品国产亚洲av天美| 亚洲精品国产色婷婷电影| 夜夜看夜夜爽夜夜摸| 妹子高潮喷水视频| 亚洲av福利一区| 国产白丝娇喘喷水9色精品| 国产日韩欧美视频二区| 中国三级夫妇交换| 人妻制服诱惑在线中文字幕| 97精品久久久久久久久久精品| 十分钟在线观看高清视频www| 精品少妇黑人巨大在线播放| 午夜激情av网站| 一级毛片 在线播放| 国产探花极品一区二区| 曰老女人黄片| 亚洲国产精品专区欧美| 亚洲在久久综合| 水蜜桃什么品种好| 成人国产av品久久久| 国产精品免费大片| 黑丝袜美女国产一区| 婷婷成人精品国产| √禁漫天堂资源中文www| 插逼视频在线观看| 亚洲av不卡在线观看| 成人影院久久| 满18在线观看网站| 女的被弄到高潮叫床怎么办| 精品一区在线观看国产| 国产免费福利视频在线观看| 午夜精品国产一区二区电影| 亚洲精品国产色婷婷电影| 制服丝袜香蕉在线|