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

    国产aⅴ精品一区二区三区波| 成年女人毛片免费观看观看9| 一个人免费在线观看的高清视频| 在线观看免费视频日本深夜| 女警被强在线播放| 久久亚洲精品不卡| 午夜福利18| 中文亚洲av片在线观看爽| 在线看三级毛片| 欧美在线一区亚洲| 亚洲18禁久久av| 国产精品久久久久久久电影 | 亚洲人成电影免费在线| 久久天躁狠狠躁夜夜2o2o| 在线观看美女被高潮喷水网站 | 日日干狠狠操夜夜爽| 99国产综合亚洲精品| 丁香欧美五月| 麻豆久久精品国产亚洲av| 一本大道久久a久久精品| 国产精品亚洲一级av第二区| 老熟妇仑乱视频hdxx| 国产视频一区二区在线看| 熟妇人妻久久中文字幕3abv| 国产单亲对白刺激| 欧美乱妇无乱码| av片东京热男人的天堂| a级毛片在线看网站| 久久香蕉精品热| 观看免费一级毛片| 国产人伦9x9x在线观看| 久久国产精品人妻蜜桃| 欧美日韩亚洲国产一区二区在线观看| 欧美日韩亚洲国产一区二区在线观看| 老司机午夜十八禁免费视频| 性欧美人与动物交配| 欧美一级毛片孕妇| 国产成人av教育| 少妇裸体淫交视频免费看高清 | 久久久国产欧美日韩av| 午夜福利在线在线| 亚洲人成网站在线播放欧美日韩| 国产单亲对白刺激| 日韩有码中文字幕| 成人av在线播放网站| 免费在线观看黄色视频的| 操出白浆在线播放| 成人国语在线视频| 久久精品91无色码中文字幕| 色老头精品视频在线观看| 97超级碰碰碰精品色视频在线观看| 免费观看人在逋| 男女下面进入的视频免费午夜| 一区二区三区激情视频| 亚洲熟女毛片儿| 亚洲美女黄片视频| 久久精品aⅴ一区二区三区四区| av欧美777| 亚洲精品美女久久av网站| 一个人免费在线观看电影 | 男女视频在线观看网站免费 | 不卡一级毛片| 国产一区二区三区视频了| 国产午夜福利久久久久久| 三级男女做爰猛烈吃奶摸视频| 级片在线观看| 午夜福利18| 国产精品久久电影中文字幕| xxx96com| 日韩大码丰满熟妇| 好男人在线观看高清免费视频| 深夜精品福利| 久久久久久人人人人人| 91在线观看av| 免费看日本二区| 这个男人来自地球电影免费观看| 黄片小视频在线播放| 丰满人妻一区二区三区视频av | 国产精品美女特级片免费视频播放器 | 999精品在线视频| 国产不卡一卡二| 日本免费a在线| 操出白浆在线播放| 免费观看精品视频网站| 欧美高清成人免费视频www| 中文字幕av在线有码专区| www.精华液| 成人国语在线视频| 亚洲成人久久爱视频| 午夜免费激情av| 国产成人精品无人区| 嫩草影院精品99| 久久精品影院6| 国产私拍福利视频在线观看| 国产精品久久久av美女十八| 亚洲 国产 在线| 日本黄大片高清| 变态另类成人亚洲欧美熟女| 色尼玛亚洲综合影院| 夜夜夜夜夜久久久久| 日韩av在线大香蕉| 国产麻豆成人av免费视频| 国产人伦9x9x在线观看| 在线十欧美十亚洲十日本专区| 午夜福利高清视频| 黄片小视频在线播放| 人成视频在线观看免费观看| 亚洲成a人片在线一区二区| 脱女人内裤的视频| 在线a可以看的网站| 久久久国产欧美日韩av| videosex国产| 日韩大码丰满熟妇| 真人一进一出gif抽搐免费| 老司机福利观看| 成人欧美大片| 亚洲色图av天堂| 国产高清视频在线播放一区| av在线播放免费不卡| 变态另类丝袜制服| 51午夜福利影视在线观看| av福利片在线| 蜜桃久久精品国产亚洲av| 男男h啪啪无遮挡| 久久香蕉精品热| 精华霜和精华液先用哪个| 男女床上黄色一级片免费看| 51午夜福利影视在线观看| 99热6这里只有精品| xxx96com| 亚洲自偷自拍图片 自拍| 变态另类成人亚洲欧美熟女| 亚洲国产欧美网| 亚洲一区二区三区色噜噜| 中文亚洲av片在线观看爽| xxx96com| 日本熟妇午夜| 中文亚洲av片在线观看爽| 欧美精品啪啪一区二区三区| 久热爱精品视频在线9| 欧美 亚洲 国产 日韩一| 久久精品91蜜桃| 成人三级黄色视频| 国产亚洲欧美在线一区二区| 日本 欧美在线| 夜夜夜夜夜久久久久| 淫妇啪啪啪对白视频| 欧美不卡视频在线免费观看 | 国产欧美日韩一区二区精品| 一卡2卡三卡四卡精品乱码亚洲| 色综合婷婷激情| 日韩免费av在线播放| 美女高潮喷水抽搐中文字幕| 亚洲专区中文字幕在线| 亚洲真实伦在线观看| 香蕉av资源在线| 亚洲精品美女久久久久99蜜臀| 色播亚洲综合网| 最新在线观看一区二区三区| xxx96com| 99国产精品一区二区三区| 亚洲 国产 在线| 99精品欧美一区二区三区四区| 国产精品 国内视频| 国产区一区二久久| 亚洲av中文字字幕乱码综合| 国产99白浆流出| 亚洲人成77777在线视频| 色av中文字幕| 中文亚洲av片在线观看爽| 亚洲在线自拍视频| 伊人久久大香线蕉亚洲五| 久久精品国产亚洲av高清一级| 亚洲天堂国产精品一区在线| 此物有八面人人有两片| 最近最新免费中文字幕在线| 亚洲成av人片在线播放无| 午夜a级毛片| 午夜精品在线福利| 日韩精品中文字幕看吧| videosex国产| xxx96com| 日韩 欧美 亚洲 中文字幕| 毛片女人毛片| 亚洲午夜精品一区,二区,三区| 99精品欧美一区二区三区四区| 亚洲欧美激情综合另类| 国产亚洲精品久久久久久毛片| 麻豆国产97在线/欧美 | 久久国产精品影院| www.999成人在线观看| 免费在线观看成人毛片| 全区人妻精品视频| 亚洲国产欧美人成| www.熟女人妻精品国产| 国产亚洲精品第一综合不卡| 俺也久久电影网| 色综合婷婷激情| 亚洲片人在线观看| 国产日本99.免费观看| 欧美日韩国产亚洲二区| 国产精品久久久人人做人人爽| 亚洲成人久久爱视频| 国产精品乱码一区二三区的特点| 国产午夜精品久久久久久| 国产69精品久久久久777片 | 国产91精品成人一区二区三区| 久久这里只有精品19| 国产精品乱码一区二三区的特点| 99在线人妻在线中文字幕| 特大巨黑吊av在线直播| 一进一出抽搐gif免费好疼| 性色av乱码一区二区三区2| 国产高清videossex| 97碰自拍视频| tocl精华| 国产精品98久久久久久宅男小说| 免费观看人在逋| 国产在线精品亚洲第一网站| 午夜福利18| 给我免费播放毛片高清在线观看| 欧美日韩瑟瑟在线播放| 美女午夜性视频免费| 免费搜索国产男女视频| 国产成人影院久久av| 狠狠狠狠99中文字幕| 精品欧美国产一区二区三| 亚洲精品av麻豆狂野| 亚洲真实伦在线观看| 久久亚洲真实| 亚洲欧洲精品一区二区精品久久久| 日本成人三级电影网站| 久久精品国产清高在天天线| 男人的好看免费观看在线视频 | 国产成人系列免费观看| 18禁裸乳无遮挡免费网站照片| 啦啦啦观看免费观看视频高清| 久久欧美精品欧美久久欧美| 亚洲狠狠婷婷综合久久图片| 99久久综合精品五月天人人| 中文亚洲av片在线观看爽| 国产一区二区在线av高清观看| 久久这里只有精品19| 亚洲一区二区三区不卡视频| 精品久久久久久,| 亚洲av日韩精品久久久久久密| 免费在线观看黄色视频的| 波多野结衣巨乳人妻| 人人妻人人澡欧美一区二区| 99国产综合亚洲精品| 欧美丝袜亚洲另类 | 很黄的视频免费| 黑人操中国人逼视频| 久久国产乱子伦精品免费另类| 日日干狠狠操夜夜爽| 日本三级黄在线观看| 麻豆成人av在线观看| 国产精品久久久久久精品电影| 男人的好看免费观看在线视频 | 国产精品自产拍在线观看55亚洲| 啦啦啦韩国在线观看视频| 日本 av在线| 99久久精品国产亚洲精品| 亚洲 欧美 日韩 在线 免费| 啪啪无遮挡十八禁网站| 欧美 亚洲 国产 日韩一| 久久香蕉国产精品| 九色国产91popny在线| 一级a爱片免费观看的视频| 亚洲人成电影免费在线| 久热爱精品视频在线9| 欧美日韩亚洲综合一区二区三区_| 两个人看的免费小视频| 欧美黑人精品巨大| 亚洲免费av在线视频| 亚洲成人久久性| 精品久久久久久久末码| 手机成人av网站| 国产精品爽爽va在线观看网站| 国产精品一区二区免费欧美| 国产私拍福利视频在线观看| svipshipincom国产片| 最好的美女福利视频网| 国产av又大| 免费无遮挡裸体视频| 99国产综合亚洲精品| 啦啦啦韩国在线观看视频| 国产精品久久久久久精品电影| 欧美中文日本在线观看视频| 日本五十路高清| 午夜精品在线福利| 亚洲成人免费电影在线观看| 一边摸一边抽搐一进一小说| 婷婷精品国产亚洲av| 香蕉丝袜av| 日日爽夜夜爽网站| 长腿黑丝高跟| 国产成人av激情在线播放| 欧美黑人精品巨大| 欧美另类亚洲清纯唯美| 午夜免费成人在线视频| 国产成人一区二区三区免费视频网站| 国产午夜福利久久久久久| av超薄肉色丝袜交足视频| 熟妇人妻久久中文字幕3abv| 91av网站免费观看| 香蕉久久夜色| 亚洲九九香蕉| 国产午夜精品论理片| 亚洲精品国产精品久久久不卡| 国产人伦9x9x在线观看| 精品国产美女av久久久久小说| 婷婷丁香在线五月| 亚洲乱码一区二区免费版| bbb黄色大片| 色哟哟哟哟哟哟| 一边摸一边抽搐一进一小说| 一级黄色大片毛片| 亚洲国产欧美一区二区综合| 成人永久免费在线观看视频| 亚洲九九香蕉| 国产伦一二天堂av在线观看| 久久中文字幕人妻熟女| 真人一进一出gif抽搐免费| 国产精品一区二区三区四区免费观看 | 老汉色∧v一级毛片| 免费看十八禁软件| 老司机在亚洲福利影院| 夜夜爽天天搞| 国产黄片美女视频| 亚洲欧美一区二区三区黑人| 超碰成人久久| 黄频高清免费视频| 国产免费男女视频| 男女之事视频高清在线观看| 亚洲自偷自拍图片 自拍| 精品久久久久久久毛片微露脸| 国产99久久九九免费精品| 深夜精品福利| 精品国产乱子伦一区二区三区| 久久中文看片网| 高清在线国产一区| 国产私拍福利视频在线观看| 91大片在线观看| 国产精品久久久人人做人人爽| 欧美另类亚洲清纯唯美| 久久这里只有精品中国| 蜜桃久久精品国产亚洲av| 美女午夜性视频免费| 国产单亲对白刺激| 啦啦啦免费观看视频1| 又爽又黄无遮挡网站| 日韩免费av在线播放| 1024手机看黄色片| 夜夜看夜夜爽夜夜摸| 舔av片在线| 久久久久久久久免费视频了| 亚洲成人国产一区在线观看| 亚洲av成人精品一区久久| 日韩欧美三级三区| 夜夜躁狠狠躁天天躁| 中文字幕人妻丝袜一区二区| 国产探花在线观看一区二区| 色综合站精品国产| 亚洲成av人片在线播放无| 亚洲av电影在线进入| 波多野结衣巨乳人妻| av有码第一页| 亚洲国产日韩欧美精品在线观看 | 久久精品影院6| 亚洲欧美日韩高清在线视频| 麻豆国产97在线/欧美 | 亚洲精品美女久久av网站| 久久热在线av| 免费高清视频大片| 国产精品电影一区二区三区| 久久精品综合一区二区三区| 黑人操中国人逼视频| 天天添夜夜摸| 亚洲va日本ⅴa欧美va伊人久久| 男女之事视频高清在线观看| 成人国产综合亚洲| 男女视频在线观看网站免费 | 日本一区二区免费在线视频| av天堂在线播放| 夜夜夜夜夜久久久久| 欧美日本亚洲视频在线播放| 91麻豆av在线| 草草在线视频免费看| 丰满人妻熟妇乱又伦精品不卡| 色综合亚洲欧美另类图片| 级片在线观看| 亚洲熟妇中文字幕五十中出| 波多野结衣巨乳人妻| 露出奶头的视频| 黑人操中国人逼视频| 亚洲九九香蕉| 国产精品一区二区三区四区免费观看 | 日本 av在线| 亚洲国产精品久久男人天堂| 日韩欧美三级三区| 可以在线观看毛片的网站| 久久精品综合一区二区三区| 香蕉久久夜色| 日韩av在线大香蕉| 成人特级黄色片久久久久久久| 后天国语完整版免费观看| 国产爱豆传媒在线观看 | 久久精品aⅴ一区二区三区四区| 国产精品永久免费网站| 男女床上黄色一级片免费看| 中文字幕av在线有码专区| 精品久久久久久久毛片微露脸| 亚洲一区二区三区不卡视频| 啪啪无遮挡十八禁网站| 精品国产超薄肉色丝袜足j| 真人一进一出gif抽搐免费| 日本免费一区二区三区高清不卡| www.www免费av| 天堂影院成人在线观看| 又大又爽又粗| 人人妻人人看人人澡| 亚洲精品美女久久av网站| 国产精品久久久久久久电影 | 成人国产综合亚洲| 亚洲精品粉嫩美女一区| 亚洲免费av在线视频| 神马国产精品三级电影在线观看 | 国产视频一区二区在线看| 日本免费a在线| 久久久久久久久免费视频了| 久久精品国产清高在天天线| 级片在线观看| 国产av一区二区精品久久| 亚洲成a人片在线一区二区| 亚洲成人久久性| 在线观看免费日韩欧美大片| 国产精品乱码一区二三区的特点| 久久久久久大精品| 一本精品99久久精品77| 国产av在哪里看| 久久午夜亚洲精品久久| 变态另类成人亚洲欧美熟女| 成人午夜高清在线视频| 一本大道久久a久久精品| 少妇裸体淫交视频免费看高清 | 777久久人妻少妇嫩草av网站| 国产精品精品国产色婷婷| 高潮久久久久久久久久久不卡| 久久精品国产清高在天天线| 操出白浆在线播放| 在线十欧美十亚洲十日本专区| 亚洲男人天堂网一区| 日韩欧美在线二视频| 美女免费视频网站| 日韩有码中文字幕| 久久久久国产精品人妻aⅴ院| 亚洲人成77777在线视频| 一级毛片女人18水好多| 国产精品一及| 精品久久久久久成人av| 亚洲成人精品中文字幕电影| 99热6这里只有精品| 亚洲成人久久性| 操出白浆在线播放| 免费在线观看影片大全网站| 又大又爽又粗| 国产乱人伦免费视频| 两个人视频免费观看高清| 色综合婷婷激情| 国产精品九九99| 亚洲人成网站高清观看| 国产精品一及| 国产三级在线视频| 国产v大片淫在线免费观看| 久久天堂一区二区三区四区| 波多野结衣巨乳人妻| 亚洲七黄色美女视频| 精品少妇一区二区三区视频日本电影| 99在线视频只有这里精品首页| 国产男靠女视频免费网站| 亚洲精品国产一区二区精华液| 88av欧美| 亚洲一区二区三区不卡视频| 免费看十八禁软件| 嫩草影视91久久| 成熟少妇高潮喷水视频| 波多野结衣高清作品| 欧美乱码精品一区二区三区| 欧美成人一区二区免费高清观看 | 国产av不卡久久| 好看av亚洲va欧美ⅴa在| 岛国在线免费视频观看| 精品国产超薄肉色丝袜足j| 在线看三级毛片| 国产一区在线观看成人免费| 淫妇啪啪啪对白视频| 听说在线观看完整版免费高清| 激情在线观看视频在线高清| 国产熟女午夜一区二区三区| 亚洲成人精品中文字幕电影| 久久国产精品人妻蜜桃| 亚洲九九香蕉| 亚洲中文字幕日韩| 老司机午夜十八禁免费视频| 国产av一区二区精品久久| 国产av麻豆久久久久久久| 性色av乱码一区二区三区2| 日韩中文字幕欧美一区二区| 欧美色欧美亚洲另类二区| 国产一区二区在线观看日韩 | 国产av又大| 亚洲成人精品中文字幕电影| 欧美大码av| 露出奶头的视频| av免费在线观看网站| 可以在线观看的亚洲视频| av国产免费在线观看| 午夜精品一区二区三区免费看| 少妇裸体淫交视频免费看高清 | 色播亚洲综合网| 成人av一区二区三区在线看| 老司机福利观看| 免费高清视频大片| 欧美日韩瑟瑟在线播放| 99在线视频只有这里精品首页| 亚洲国产精品999在线| 亚洲人成网站高清观看| 一个人免费在线观看电影 | 十八禁网站免费在线| 国产亚洲精品一区二区www| 日韩大码丰满熟妇| 国产野战对白在线观看| 欧美大码av| 免费看美女性在线毛片视频| 日韩中文字幕欧美一区二区| 国产精品亚洲av一区麻豆| 亚洲男人天堂网一区| 又黄又粗又硬又大视频| 国产成人影院久久av| 波多野结衣高清无吗| 欧美国产日韩亚洲一区| 免费高清视频大片| 国产精品国产高清国产av| 欧美日韩一级在线毛片| 亚洲片人在线观看| 人人妻人人看人人澡| 精品久久久久久久末码| 婷婷精品国产亚洲av在线| 天堂影院成人在线观看| 久久中文字幕一级| 国产一级毛片七仙女欲春2| 亚洲第一欧美日韩一区二区三区| 亚洲熟妇中文字幕五十中出| 亚洲一区中文字幕在线| 亚洲人与动物交配视频| 免费av毛片视频| 91在线观看av| 久久性视频一级片| 少妇人妻一区二区三区视频| 99热这里只有是精品50| 亚洲欧美一区二区三区黑人| 老熟妇仑乱视频hdxx| 法律面前人人平等表现在哪些方面| videosex国产| 国产视频一区二区在线看| 国产97色在线日韩免费| 男女床上黄色一级片免费看| 欧美激情久久久久久爽电影| 99久久国产精品久久久| 国产欧美日韩一区二区三| 女人爽到高潮嗷嗷叫在线视频| 两性夫妻黄色片| 午夜福利欧美成人| 一区二区三区国产精品乱码| 麻豆av在线久日| 亚洲美女黄片视频| 一级毛片女人18水好多| 亚洲国产欧美一区二区综合| 国产av在哪里看| 亚洲在线自拍视频| 人成视频在线观看免费观看| 波多野结衣高清无吗| xxx96com| 亚洲aⅴ乱码一区二区在线播放 | 精品欧美一区二区三区在线| 国产黄a三级三级三级人| 亚洲av五月六月丁香网| 午夜老司机福利片| 日韩欧美三级三区| 日韩国内少妇激情av| 在线观看免费午夜福利视频| 国产一区在线观看成人免费| 看片在线看免费视频| 18禁裸乳无遮挡免费网站照片| 国产在线观看jvid| 99国产精品99久久久久| a在线观看视频网站| 成人三级黄色视频| 真人做人爱边吃奶动态| 最近最新免费中文字幕在线| 亚洲自拍偷在线| 可以在线观看毛片的网站| 色精品久久人妻99蜜桃| 欧美日本亚洲视频在线播放| 久久久精品大字幕| 一个人观看的视频www高清免费观看 | 午夜激情av网站| 欧美色欧美亚洲另类二区| 国产精品一区二区三区四区久久| 老司机在亚洲福利影院| 亚洲一区高清亚洲精品| av有码第一页| 精品国产亚洲在线| 法律面前人人平等表现在哪些方面|