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

    Future Changes in Extreme High Temperature over China at 1.5°C–5°C Global Warming Based on CMIP6 Simulations

    2021-02-26 08:22:26GuweiZHANGGangZENGXiaoyeYANGandZhihongJIANG
    Advances in Atmospheric Sciences 2021年2期

    Guwei ZHANG, Gang ZENG, Xiaoye YANG, and Zhihong JIANG

    Key Laboratory of Meteorological Disaster of Ministry of Education (KLME), Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters (CIC-FEMD), Nanjing University of Information Science and Technology, Nanjing 210044, China

    ABSTRACT

    Key words:extreme high temperature,China,CMIP6,1.5°C–5°C global warming

    1.Introduction

    The Intergovernmental Panel on Climate Change(IPCC, 2013) in its 5th Assessment Report reported that from 1880 to 2010, the global mean temperature (GMT) has risen about 0.85°C (0.65°C–1.06°C). In order to reduce the adverse impacts caused by a rapid global warming scenario,the Conference of the Parties of the United Nations Framework Convention on Climate Change (UNFCCC), in December 2015, established a goal to limit global warming to“well below 2.0°C” and to persist with efforts to limit the warming under 1.5°C above pre-industrial levels(UNFCCC, 2015). Limiting global warming under 1.5°C is deemed to be important by many researchers (King et al.,2017; Nangombe et al., 2018; Zhang et al., 2020a; Zhao et al., 2020). The World Meteorological Organization (WMO)stated that 2010–2019 was the hottest decade on record(WMO, 2020). Occurrences of extreme high temperature(EHT) events, which has enormous socio-economic and human health impacts, were becoming more frequent, consistent with the observed global warming over the last century(Robine et al., 2008; Ding et al., 2010; Seneviratne et al.,2016; Mora et al., 2017; Zhang et al., 2020b). The heatwave in Europe in 2003 established a new record for the hottest summer and killed more than 4000 people (Barriopedro et al., 2011). Long-term EHT events will damage infrastructure, overwhelm power and water facilities, and have a significant socio-economic impact (Wilbanks et al., 2012). In developing countries, under a scenario of a global temperature rise of 2°C, the increase in thermal stress will lead to the loss of productivity for one month per year (Yu et al.,2019).

    China is the most populous country and a huge energy consumer. Against the background of global warming, the rate of temperature rise across China has been 0.9°C–1.5°C since 1909, slightly higher than the global average (The Third National Assessment Report on Climate Change,2015). Most previous studies are based on CMIP5 simulations and the Community Earth System Model (CESM)low-warming projections, and they mainly focused on the EHT changes, specific to the 1.5°C and 2.0°C global warming scenarios (Sanderson et al., 2017; Lin et al., 2018; Shi et al., 2018a; Yang et al., 2018; Yu et al., 2018). For example,Li et al. (2018) and Zhang et al. (2020b) both found that the annual temperature in China was projected to increase at a rate that is approximately 10% higher than the global average level in the 1.5°C and 2.0°C warming scenarios. According to the conclusion from Yu et al. (2018) and Shi et al.(2018a), EHT events in China will increase significantly at the 1.5°C and 2.0°C thresholds, especially for the indices affected by maximum temperature. In addition, some studies also focused on even more extreme scenarios, such as 3°C, 4°C, and 5°C (Wang et al., 2018; Weber et al., 2018).Xu et al. (2017) and concluded that the mean surface temperature over Asia would increase about 4.6°C and 6.0°C at 3°C and 4°C global warming, and more substantial warming would occur in high latitudes than in low latitudes. Overall, most studies conclude that China will be warmer than global averages (Hu et al., 2017; Li et al., 2018; Shi et al.,2018b). However, due to the earlier modeling framework not accounting for population changes, previous studies have rarely considered the number of people in China that will suffer due to China’s warming exceeding that of the global average level.

    The Coupled Model Intercomparison Project Phase 6(CMIP6) provides the latest outputs of many climate models developed by institutions around the world for Scenario Model Intercomparison Project designed for climate projection in different emission scenarios. The new phase of the CMIP is designed to improve climate simulations and provide for additional modeling groups that are expected to be more reliable (Eyring et al., 2016). Some researchers have already examined the newly released simulations(Chen et al., 2020; Jiang et al., 2020; Yang et al., 2020;Zhou et al., 2020; Zhu et al., 2020). Zelinka et al. (2020) pointed out that for the implied social scenarios in CMIP6, the models responded in such a way that suggested higher equilibrium climate sensitivity. Chen et al. (2020) have compared the simulating climate extremes in CMIP6 and CMIP5 models, and they found that the models in CMIP6 had finer resolution that ultimately resulted in improved dynamical processes. Zhu et al. (2020) concluded that, compared with CMIP5, the CMIP6 multi-model ensemble mean not only simulates the spatial pattern of air temperature well, but out performs the CMIP5 in modeling China's climate index. Since the CMIP6 simulations were recently released and designed to incorporate population changes, they are well-suited to studying future changes in China's EHT and determining how many people in China will experience warming above the global average level.

    This study aims to analyze future changes in EHT based on the recently released CMIP6 outputs, and addresses following questions at the 1.5°C, 2°C, 3°C, 4°C,and 5°C global warming scenarios; (1) How many people in China will experience warming at higher than global average levels? (2) What is the future trend regarding the frequency of EHT in China? (3) How many impacts might be avoided by limiting warming under 1.5°C compared to the higher scenarios?

    2.Data

    Sixteen climate models of the CMIP6 archives which contain daily surface air temperature (SAT) were used in this study (Eyring et al., 2016). Their basic information is listed in Table 1, and further details can be found at https://esgfnode.llnl.gov/projects/cmip6/. The historical simulations(1850–2014) and 21st-century projections (2015–2100)under three future scenarios (SSP1-2.6, SSP2-4.5, and SSP5-8.5) are used, and only the first member (e.g.,r1i1p1f1) is selected for each model. These scenarios are the combination of Shared Socio-economic Pathways(SSPs; O’Neill et al., 2017) and forcing levels of the Representative Concentration Pathways (RCP). The SSP1-2.6,SSP2-4.5, and SSP5-8.5 approximately follow RCP2.6,RCP4.5, and RCP8.5 global forcing pathways with SSP1,SSP2, and SSP5 socio-economic conditions. In this study,the outputs of each CMIP6 model are bi-linearly interpolated onto the resolution of 1.0° × 1.0° grid. The EHT indicators were calculated using the original grids of each model and then bi-linearly interpolated to the 1.0° × 1.0° spatial resolution to maintain the internal consistency.

    The population datasets under three SSPs (SSP1, SSP2,and SSP5), named Global One-Eighth Degree Population Base Year and Projection Grids Based on the SSPs (V1.01),were used in this study (Jones and O’Neill, 2016, 2020).They can be considered as a sustainable development scenario, a business-as-usual scenario, and a fossil-fueled development scenario. The model uses global total population data at 10-year intervals from 2010 to 2100, with a horizontal resolution of 0.125° × 0.125°. The datasets were developed by the Center for International Earth Science Information Network (CIESIN), Columbia University, and were obtained from the NASA Socio-economic Data and Applications Center (SEDAC). Spatially explicit population projections of the SSP project are up-scaled to a 1.0° × 1.0° spatial resolution to match the resolution of the climate model.

    Table 1. Basic information (name and group) and atmospheric resolution (latitude × longitude) of sixteen CMIP6 global climate models.

    3.Methods

    3.1.Periods of preindustrial, present-day, and specific warming thresholds

    The period spanning 1995–2014 in the historical simulation represents the present-day, which is the same as the definition in Jiang et al. (2020). Consistent with previous researchers that used CMIP5, the pre-industrial period is defined to be 1861–1900 (Xu et al., 2017; Fu et al., 2018). The specific global warming periods in this study (i.e., 1.5°C, 2°C,3°C, 4°C, and 5°C) use the pre-industrial level as a baseline for comparison. For each climate model, we calculated the 20-year moving average of the GMT anomaly (by subtracting it from the average of the pre-industrial period) and then found the particular year in this time series that reached a specific warming threshold, which is then defined as the final year. We designate the final year as the temporal mid-point and take 10 years forward and nine years backward in order to establish the specific 20-year warming period. Consistent with previous studies (Collins et al., 2013; Zhou et al.,2018; Jiang et al., 2020), the changes in this study are deemed to be robust if at least 2/3 of total models agree on the sign of change.

    3.2.Extreme high temperature indices

    This study used three EHT indices, warmest day (TXx),warm days (TX90p), and warm spell duration indicator(WSDI) to analyze the impacts of global warming on EHT in China. Their definitions are shown in Table 2 (Meehl and Tebaldi, 2004; Smith et al., 2013). The present-day period(1995?2014) is referred to as the baseline. The EHT indices are calculated based on all days of the present-day period and the future period (2015?2100). The calculations of TX90p and WSDI use the 90th percentile of the baseline maximum temperature as their threshold. These respective values are then tallied individually for each calendar day. For example, the threshold for January 1st should exceed 90th percentile of the days on the same date during the present-day period. Therefore, TX90p is defined as the number of days with the daily maximum SAT higher than the threshold, and WSDI is defined as the number of occurrences with at least six consecutive days that had a maximum SAT higher than the threshold (Table 2). Furthermore, we calculated the indices separately for each model and averaged them to display as the multi-model ensemble mean (MME).

    3.3.Response to global warming

    To investigate the response of EHT to global warming,the projections of EHT indices and GMT are averaged over 10-year periods to eliminate interannual variability (Collins et al., 2013; Jiang et al., 2020). From 2016 to 2100, a fiveyear moving average is calculated to obtain a running-mean time series of sixteen values (i.e., the averages of 2016?2025, 2021?2030, and up to 2091?2100). The linear regression between these running-mean EHT indices and GMT is then computed, and the resulting regression coefficient is referred to as the response rate of EHT to global warming, which is recognized as the long-term forced signal. It is used to figure out the amount by which EHT indices will change in response to a global mean temperature rise of 1°C.

    3.4.Avoided Impacts

    The avoided impacts (AI) of EHT caused by additional warming were studied by using formula 1 (Li et al., 2018).

    Here AI represents avoided impacts, Cand Crepresent the changes in EHT indices at the x°C (x can be 2, 3, 4 and 5) and 1.5°C global warming relative to the present-day level. The AI is used to infer the degree of impact that canbe avoided by maintaining global warming at 1.5°C, as opposed to other higher warming thresholds (2°C, 3°C, 4°C and 5°C).

    Table 2. Definitions of extreme high temperature indices.

    In addition, we also calculated the avoided impacts of population-weighted changes (AIpop), by applying formula 2.

    Here, CPand CPrepresent the changes in the population-weighted EHT indices at the x°C (x can be 2, 3, 4 and 5) and 1.5°C scenarios compared with the present-day level.For example, CPis the result of EHT indices multiplied by the population at 1.5°C global warming minus the EHT indices multiplied by the population at the present-day level. The calculations above are performed at each specific grid point..

    4.Results

    4.1.Future changes in surface air temperature

    As shown in Figures 1 and 2a, the results of CMIP6 indicate that global warming following the high-emission pathway (e.g., SSP5-8.5) will increase more dramatically than that following the low-emission pathway (e.g. SSP1-2.6),which is consistent with the previous conclusions derived from CMIP5. For example, the specific warming threshold under SSP5-8.5 will be reached sooner compared to the threshold times indicated for the SSP1-2.6 and SSP2-4.5 scenarios (Fig.1). The multi-model ensemble mean (MME) results show that the annual GMT anomalies relative to preindustrial levels under SSP5-8.5 will exceed 5°C before 2100, while under SSP1-2.6 and SSP2-4.5 it will only reach 2°C and 3°C. Although there are some differences between the time sensitivity of each model, it can be concluded that global warming will reach 1.5°C (2°C) before 2030 (2050)under the aforementioned three future scenarios. As displayed in Fig.2a, after 2050, global warming under the SSP1-2.6 scenario projects the warming to become stable,while it will continue to increase under SSP2-4.5 and SSP5-8.5. Under SSP2-4.5, global warming will reach 3°C at about 2070. Under SSP5-8.5, global warming will reach 3°C, 4°C, and 5°C before 2060, 2080, and 2100, respectively.

    Fig.1. The timing of reaching the specific warming threshold under SSP1-2.6(green), SSP2-4.5 (blue) and SSP5-8.5 (red). The MME indicates the multimodel ensemble mean results.

    Fig.2. (a) Annual global mean surface air temperature (SAT) anomalies relative to pre-industrial levels (1861?1900) from 1850 to 2100. The black line is from 1850 to 2014 (the historical period), and the green, blue and red lines are from 2015 to 2100 (future period) under SSP1-2.6, SSP2-4.5, and SSP5-8.5, respectively. Shading indicates the uncertainties of the multimodel ensemble. (b) China’s regional mean SAT changes relative to the pre-industrial level at 1.5°C–5°C global warming.Green, blue and red bars represent results of SSP1-2.6, SSP2-4.5, and SSP5-8.5. The perpendicular black lines at the top of the bars represent the variation between the 10th and 90th percentiles of the multi-model ensemble.

    Many previous studies have pointed out that at 1.5°C and 2°C global warming, China’s temperature rise (relative to pre-industrial levels) will exceed the global average level.The CMIP6 results also show that not only at 1.5°C and 2°C, but also at the 3°C, 4°C, and 5°C global warming scenarios, the increases in China’s regional mean SAT will be higher than the global average (Fig.2b). Under the 1.5°C and 2°C scenarios, the annual SAT over China is projected to increase by 1.55°C and 2.12°C, under the SSP1-2.6, by 1.59°C and 2.18°C under the SSP2-4.5, and by 1.60°C and 2.20°C under SSP5-8.5. Intuitively, the temperature increases are projected to be greatest under SSP5-8.5 and least pronounced under SSP1-2.6. Furthermore, when the warming threshold increases, the average SAT in China will increase more than the global average level. For example,the SSP5-8.5 scenario, which suggests when global warming reaches 3°C, the average increase of regional SAT in China is approximately 3.40°C, which is about 13% higher than the global level, while the increase at 2°C global warming (2.20°C) is about 10% higher than the global level.

    The spatial patterns of changing SAT relative to the pre-industrial level, displayed in Fig.3, are used to display which regions will experience warming above the global average. In order to reduce the complicated description, we named the areas that will experience warming above the global average level as overheating areas. For example, at 1.5°C global warming, the areas where the increases in SAT exceed 1.5°C are the overheating regions. As shown in Fig.3,the overheating areas at 1.5°C–5°C are mostly located in Tibet and northern China (i.e., Northwest China, North China and Northeast China). This is similar to the conclusions drawn by previous studies that northern China and Tibet will experience more warming in the future (Zhang et al., 2020b). Analogous results, under different emission scenarios, found that the high-emission scenario (SSP5-8.5) will have more overheating area. Under SSP1-2.6/SSP2-4.5/SSP5-8.5 (Fig.4a), about 56%/57%/58% and 60%/61%/64% areas of China will experience warming above the global level at 1.5°C and 2°C global warming. At 3°C global warming, the land fraction of overheating areas in China will be 52% under SSP2-4.5, while the value under SSP5-8.5 will reach 69%. Meanwhile, when global warming increase to 4°C and 5°C under SSP5-8.5, the land fraction of overheating areas in China will reach nearly 70%.

    The population in these overheating areas of China(Fig.4b) indicate that there will be more people living in the overheating regions under SSP1-2.6 and SSP5-8.5 than under SSP2-4.5. For example, at 1.5°C (2°C) global warming, the population in the overheating areas under SSP1-2.6 and SSP5-8.5 will be 0.29 (0.40) billion and 0.30 (0.40) billion, accounting for 21% (30%) and 22% (30%) of China’s total population, while the number under SSP2-4.5 will be 0.27 (0.31) billion, accounting for 19% (30%) of China’s total population. The number at 3°C global warming under SSP2-4.5 will decrease to 0.19 billion (16% of China’s total population). While under SSP5-8.5, when global warming rises to 3°C/4°C/5°C, there will be about 0.44/0.41/0.36 billion people in the overheating areas, accounting for 37%/41%/34% of China’s total population. Overall, approximately 0.19 to 0.44 billion people in China will experience warming that is higher than the global level in the future.

    4.2.Future changes in extreme high temperature

    In an effort to quantify the response of EHT events to global warming, we calculated the response rate of the EHT indices to GMT (Fig.5). The results indicate that EHT events in China are expected to continue to increase as global warming continues. In most parts of China, the response rates of TXx (Figs. 5a–c), TX90p (Figs. 5d–f),and WSDI (Figs. 5g–i) to global warming are higher than 0.4°C °C, 10 d °C, and 5 d °C, respectively. The response rates of TX90p and WSDI demonstrate remarkable spatial consistency as opposed to the TXx. For TX90p and WSDI, there will be a significant response in southern China, with increased rates exceeding 30 d °Cand 15 d °C. This is particularly true in southern Tibet and South China, TX90p and WSDI will increase by more than 40 and 20 days annually when the GMT goes up by 1°C. The spatial distribution of the TXx shows almost an inverse distribution, that is, the high response areas are mostly located in northern China. For illustration, under SSP5-8.5 (Fig.5c),the TXx in most parts of northern China will increase by at least 1.2°C for every 1°C rise in GMT, while the increases in the south will be lower than this.

    Fig.3. Spatial patterns of SAT changes relative to pre-industrial levels. (a), (b) and (c) are for 1.5°C global warming under SSP1-2.6, SSP2-4.5 and SSP5-8.5, respectively. (d), (e) and (f) are for 2°C global warming under SSP1-2.6, SSP2-4.5, and SSP5-8.5,respectively. (g) and (h) are for 3°C global warming under SSP2-4.5 and SSP5-8.5. (i) and (j) are for 4°C and 5°C global warming under SSP5-8.5. Dotted areas denote where at least 2/3 models agree on the sign of the change.

    We used the present-day period (1995–2014) as the baseline to discuss future changes in EHT at each specific global warming scenario (1.5°C–5°C). As shown in Fig.6,the spatial distribution of the TXx changes indicates that the increase in northern China will be greater than in southern China. For example, under the 1.5°C global warming scenario, the increased TXx in southern China will be less than 1°C, while the increases in some parts of Northeast China will be higher than 1.5°C (Figs. 6a–c). The increase under the SSP5-8.5 scenario is slightly higher than those under SSP1-2.6 and SSP2-4.5. For example, at 3°C global warming, the areas where increased TXx exceed 3°C are larger under SSP5-8.5 (Fig.6g) than those under SSP2-4.5 (Fig.6h). For the regional mean results of China (Fig.9a), the increased TXx at 1.5°C and 2°C global warming is similar in each scenario, ranging from 0.2°C–2.1°C and 1.1°C–2.9°C respectively. When global warming reaches 3°C, 4°C, and 5°C, China's regional mean TXx will increase by 2.1°C–4.2°C, 3.6°C–5.4°C, and 4.1°C–6.8°C relative to present-day levels, respectively (Fig.9a).

    Fig.4. (a) The land fraction of China where warming is higher than the global average (%), and (b) the population of these regions (billion). Green, blue and red bars represent results of SSP1-2.6, SSP2-4.5, and SSP5-8.5. The perpendicular black lines at the top of the bars represent the variation between the 10th and 90th percentiles of the multi-model ensemble.

    Fig.5. The responses of EHT changes to global warming under SSP1-2.6, SSP2-4.5 and SSP5-8.5. (a)–(c) are for TXx (°C °C?1);(d)–(f) are for TX90p (d °C?1); (g)–(h) are for WSDI (d °C?1). Dotted areas denote where at least 2/3 models agree on the sign of the change.

    Fig.6. Spatial patterns of TXx changes relative to the present-day. (a), (b) and (c) are for 1.5°C under SSP1-2.6, SSP2-4.5 and SSP5-8.5, respectively. (d), (e) and (f) are for 2°C under SSP1-2.6, SSP2-4.5, and SSP5-8.5, respectively. (g) and (h) are for 3°C under SSP2-4.5 and SSP5-8.5. (i) and (j) are for 4°C and 5°C under SSP5-8.5. Dotted areas denote where at least 2/3 models agree on the sign of the change.

    Interestingly, the spatial distribution of changing TX90p (Fig.7) and WSDI (Fig.8), is similar to the results of the response rate to global warming, yet they are still almost the reverse of TXx. Large areas of increased frequency for TX90p and WSDI will occur in southern China,especially in southern Tibet and South China. For example,when global warming reaches 3°C, the increase in TX90p(WSDI) in southern China will exceed 75 (40) days per year, with parts of southern Tibet and southern China exceeding 100 (60) days per year. Furthermore, the increases under SSP5-8.5 will be higher than those under SSP1-2.6 and SSP2-4.5. For instance, compared to SSP1-2.6 and SSP2-4.5, the WSDI will increase by 10 days or more over more areas than under SSP5-8.5 (Figs. 8a?c). Overall, compared to the present-day level, China’s TX90p (WSDI) will increase by 7–45 (3–23) days, 31–64 (15–33) days, 31–64(30–52) days, 58–105 (30–52) days, 99–146 (48–68) days,and 121–183 (60–87) days at 1.5°C, 2°C, 3°C, 4°C, and 5°C global warming, respectively (Figures 9b and c).

    4.3.Avoided Impact at 1.5°C global warming

    Fig.7. The same as Fig.6, but for the changes in TX90p (units: d).

    Formula (1) quantifies the impacts avoided at 1.5°C global warming compared with other higher warming thresholds (2°C–5°C). As shown in Figures 10a–c, in comparison to 2°C–5°C global warming, the reduced warming in 1.5°C global warming scenario will serve to avoid approximately 36%–87%,47%–89%, and 46%–86% of the increases in TXx, TX90p, and WSDI, respectively. The avoided impacts for a 2°C warming scenario project a 20% reduction compared to the results for the 3°C–5°C scenarios, all of which exceed 65%. For example, when compared to the 3°C scenario, limiting the global warming to under 1.5°C in SSP2-4.5 will help avoid 68%, 80%, and 76% of the increases in TXx, TX90p, and WSDI, respectively. Formula (2) was used to investigate the population-weighted avoided impacts. When population changes are considered,the avoided impacts of TXx will increase significantly,while the results of TX90p and WSDI will change marginally (Figures 10e–f). For the TXx, the population-weighted avoided impacts are about 80%–100% for 2°C–5°C, which is 30%–40% higher than the results without considering population changes. While for the TX90p and WSDI (Figures 10e and 10f), the results are about 10% lower than the results displayed in Figures 10b and 10c. For instance, the populationweighted avoided impacts of WSDI for 2°C–5°C are 35%–80%.

    Overall, constraining global warming to 1.5°C rather than 2°C, 3°C, 4°C, and 5°C will help avoid approximately 36%, 73%, 80%, and 87% of the EHT indices increases in China. When considering the population changes, the population-weighted avoided impact will increase to 30%–40% in the TXx.

    5.Discussions

    Fig.8. The same as Fig 6, but for the changes in WSDI (units: d).

    Our analysis yields some new insights. Given that most CMIP6 climate models projected that global warming would reach 1.5°C and 2°C before 2030 and 2050 under three future scenarios (SSP1-2.6, SSP2-4.5 and SSP5-8.5).Under SSP5-8.5, global warming will likely reach 5°C by 2100. In China, the overall temperature increases will be greater than the global average, which will result in about 0.19 to 0.44 billion people suffering above global levels of warming. For the changes in the EHT, the TXx will increase much in northern China, while the TX90p and WSDI will increase much in southern China.

    Fig.9. Compared to present-day, China’s regional mean increases in TXx (a), TX90p(b), and WSDI (c) at 1.5°C–5°C global warming. Green, blue and red box whisker plots represent results of SSP1-2.6, SSP2-4.5, and SSP5-8.5. The box-whisker plots show the 10th, 25th, 50th, 75th, and 90th percentiles of the multi-model ensemble.

    However, the merits of this study are not without its limitations. For example, in an attempt to match the same resolution, both the climate model data and population data were interpolated. The difference is that the model data was mostly processed by the bi-linearly downscaled method and the population data was processed by the bi-linearly upscaled method. These differential methods may add some uncertainties, however it is worth noting that dynamic downscaling is not only a viable methodology to improve the horizontal resolution of the models (Liang et al., 2019), but this process has also been shown to improve the model performance in China in some cases (Liang et al., 2019). Furthermore, there are still some issues worthy of further study.For instance, the time needed to reach the specific global warming thresholds (time points) seems slightly earlier than those concluded by previous studies (Xu et al., 2017; Shi et al., 2018a). This is likely attributed to the higher average climate sensitivity in CMIP6, especially regarding the response to aerosols (Flynn and Mauritsen, 2020; Zelinka et al., 2020). Limiting global warming to 1.5°C or 2.0°C will require significant reductions in anthropogenic greenhouse gas emissions. Consequently, anthropogenic aerosol emissions are expected to decline as a result of reduced greenhouse gas emissions and air quality improvement measures.Samset et al. (2018) showed that aerosol removal would cause a global surface warming of 0.5°C–1.1°C and EHT indices to increase. EHT events in major aerosol emission regions (mainly China and the U.S.) are more sensitive to aerosol emission reductions. Therefore, an in-depth study concerning the impact of aerosol removal on China is necessary. Additionally, our results indicate that the TXx would increase most notably in northern China while the increases in TX90p and WSDI would be most pronounced in southern China. Since northern China has both the largest forest and crop areas in China (Tao and Zhang, 2013; Huang et al., 2017) in addition to large expanses of arid areas,enhanced high temperatures in northern China may increase the frequency and intensity of forest fires and also lead to severe droughts (Allen et al., 2004; Chai et al., 2018). Increasing drought will affect agricultural production, such as wheat and cotton in Northwest China (Wang et al., 2008).Since southern China is an important industrial region (e.g.Shanghai, Shenzhen, Hong Kong), high temperatures are likely to result in the consumption of more energy, therefore it is necessary to investigate the impacts of high temperatures on the economy (Jiang et al., 2017; Yu et al., 2019;Zhu et al., 2020), and to quantify the future heat-related risks in China, such as energy consumption (Yu et al.,2019), forest fires (Li et al., 2019), and droughts (Huang et al., 2017; Su et al., 2018).

    In summary, this study projected future EHT changes in China under five specific warming thresholds (1.5°C–5°C), which can provide the basis for an in-depth future projection study.

    6.Conclusions

    Fig.10. The avoided impacts (%) of TXx (a), TX90p (b), and WSDI (c) at 2°C, 3°C, and 4°C global warming relative to 1.5°C global warming. (d)–(f) are the same as (a)–(c) but for the population-weighted results (%). Green, blue and red bars represent the results of SSP1-2.6, SSP2-4.5, and SSP5-8.5, respectively. The perpendicular black lines at the top of the bars represent the variation between the 10th and 90th percentiles of the multi-model ensemble.

    Understanding how climate extremes respond to global warming is critical for climate change adaptation and international climate negotiations. In this study, we use the newly released CMIP6 simulations to project the future changes of extreme high temperature (EHT) in China at 1.5°C, 2°C,3°C, 4°C, and 5°C global warming. The conclusions are summarized as follows:

    (1) CMIP6 simulations project that, under three future scenarios (SSP1-2.6, SSP2-4.5 and SSP5-8.5), the global mean temperature will increase by 1.5°C and 2°C relative to pre-industrial levels (1861–1900) before 2030 and 2050,respectively. Furthermore, the timing of the 1.5°C/2°C global warming threshold under SSP5-8.5 is reached slightly earlier than that under SSP1-2.6 and SSP2-4.5.Under SSP 5-8.5, global warming will eventually exceed 5°C by 2100. While under SSP1-2.6, global warming will stabilize at around 2°C after 2050.

    (2) Under 1.5°C, 2°C, 3°C, 4°C, and 5°C global warming, China’s regional mean temperature increase relative to pre-industrial levels will be higher than the global average.For example, when global warming reaches 1.5°C and 2°C,China’s regional mean temperature is projected to increase by 1.55°C/1.59°C/1.60°C and 2.12°C/2.18°C/2.20°C relative to the preindustrial level under SSP1-2.6/SSP2-4.5/SSP5-8.5, respectively. In China, most of the areas where warming exceeds global average levels will occur in Tibet and northern China (Northwest China, North China and Northeast China), more than half the whole country. Especially under SSP5-8.5, this land fraction will reach nearly 70% of China. Furthermore, an estimated 0.19 to 0.44 billion people, accounting for 16% to 41% of the national population, will experience warming that is higher than the global average.

    (3) Regarding changes in EHT compared to the present day, China will experience an increase in the TXx, TX90p,and WSDI by 0.2°C–6.8°C, 7–183 days, and 3–87 days at 1.5°C–5°C global warming, respectively. The highest increases will occur at 5°C global warming under SSP 5-8.5. Furthermore, the TXx will increase most dramatically in northern China, while the TX90p and WSDI increases will be greatest in southern China. For example, relative to the present-day (1995–2014), TXx in northern China will increase by at least 1°C–5°C, and TX90p (WSDI) in southern China will increase by more than 25–150 (10–80) days at 1.5°C–5°C global warming. Compared to 2°C, 3°C, 4°C,and 5°C, limiting global warming to 1.5°C will likely prevent about 36%, 73%, 80%, and 87% of the EHT increases in China. When considering population changes, the population-weighted avoided impacts of TX90p and WSDI will change not much, but the avoided impacts of TXx will increase by 30%–40%.

    Acknowledgements.

    This research is supported by the National Key Research and Development Program of China(2017YFA0603804), the National Natural Science Foundation of China (41831174 and 41430528), and the Postgraduate Research& Practice Innovation Program of Jiangsu Province (KYCX19_1026). Guwei ZHANG was supported by the China Scholarship Council (NO. 201908320503). We acknowledge the High Performance Computing Center of Nanjing University of Information Science & Technology for their support of this work. We sincerely thank the editors and reviewers for their constructive critique and positive review.

    小蜜桃在线观看免费完整版高清| 晚上一个人看的免费电影| 性插视频无遮挡在线免费观看| 淫秽高清视频在线观看| 在线免费观看的www视频| 男人舔奶头视频| 日本成人三级电影网站| 久久久久久久久久成人| 日韩视频在线欧美| 女人被狂操c到高潮| 少妇人妻精品综合一区二区 | 亚洲av.av天堂| 特大巨黑吊av在线直播| 1024手机看黄色片| videossex国产| av视频在线观看入口| 国产高清有码在线观看视频| 91麻豆精品激情在线观看国产| 欧美bdsm另类| 亚洲真实伦在线观看| 成人毛片60女人毛片免费| 国产伦在线观看视频一区| 青春草视频在线免费观看| 少妇熟女欧美另类| 午夜福利视频1000在线观看| 亚洲av第一区精品v没综合| 国产视频首页在线观看| 五月伊人婷婷丁香| 国产成人精品久久久久久| 国产精品一区二区在线观看99 | 真实男女啪啪啪动态图| 成人一区二区视频在线观看| 亚洲人成网站在线播| 国产中年淑女户外野战色| 免费观看的影片在线观看| 六月丁香七月| 综合色av麻豆| 亚洲内射少妇av| 精品午夜福利在线看| 婷婷六月久久综合丁香| 免费观看a级毛片全部| 男女啪啪激烈高潮av片| 国内少妇人妻偷人精品xxx网站| kizo精华| 波多野结衣巨乳人妻| 99精品在免费线老司机午夜| 最近最新中文字幕大全电影3| 91精品国产九色| 99热这里只有是精品在线观看| 成人二区视频| 色吧在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 校园春色视频在线观看| 麻豆精品久久久久久蜜桃| 亚洲成人久久爱视频| 在线观看午夜福利视频| 亚洲av一区综合| 国产av不卡久久| 久久99热这里只有精品18| 少妇被粗大猛烈的视频| 亚洲精品成人久久久久久| 五月玫瑰六月丁香| 观看美女的网站| 最近手机中文字幕大全| 久久精品国产亚洲网站| 高清毛片免费看| 成人午夜精彩视频在线观看| 久久鲁丝午夜福利片| 午夜福利在线观看免费完整高清在 | 一个人看的www免费观看视频| 日本免费一区二区三区高清不卡| 少妇被粗大猛烈的视频| 国产一区二区在线av高清观看| 色哟哟·www| 久久久成人免费电影| 99热网站在线观看| 欧美日韩在线观看h| 日韩欧美精品v在线| 特级一级黄色大片| 联通29元200g的流量卡| 日本一本二区三区精品| 天堂网av新在线| 国产在线精品亚洲第一网站| 九九爱精品视频在线观看| 女人被狂操c到高潮| 久久精品国产99精品国产亚洲性色| 女人被狂操c到高潮| 国产麻豆成人av免费视频| 国产美女午夜福利| 国产麻豆成人av免费视频| 亚洲国产精品成人久久小说 | 久久亚洲国产成人精品v| 午夜视频国产福利| 国产av麻豆久久久久久久| 国内精品久久久久精免费| 国产色婷婷99| 成人一区二区视频在线观看| 国产在线精品亚洲第一网站| 青春草视频在线免费观看| av专区在线播放| 免费一级毛片在线播放高清视频| 可以在线观看毛片的网站| 免费一级毛片在线播放高清视频| 中文在线观看免费www的网站| 亚洲图色成人| 一本精品99久久精品77| 国产大屁股一区二区在线视频| 国产大屁股一区二区在线视频| 十八禁国产超污无遮挡网站| 精品久久久久久久久久免费视频| 国产 一区 欧美 日韩| 又黄又爽又刺激的免费视频.| 真实男女啪啪啪动态图| 尤物成人国产欧美一区二区三区| 日韩欧美国产在线观看| 亚洲第一电影网av| 久久人妻av系列| 啦啦啦观看免费观看视频高清| 可以在线观看毛片的网站| 国产成人午夜福利电影在线观看| 国产精品不卡视频一区二区| 久久亚洲国产成人精品v| 两性午夜刺激爽爽歪歪视频在线观看| 不卡视频在线观看欧美| 看十八女毛片水多多多| 国产精品一区二区在线观看99 | 在线播放国产精品三级| 亚洲人成网站高清观看| 可以在线观看毛片的网站| 看片在线看免费视频| 小蜜桃在线观看免费完整版高清| 亚洲无线观看免费| www.av在线官网国产| 亚洲av二区三区四区| 在线播放无遮挡| 老熟妇乱子伦视频在线观看| 人人妻人人澡人人爽人人夜夜 | 99国产精品一区二区蜜桃av| 一级黄片播放器| 国产成人精品一,二区 | 美女内射精品一级片tv| 国产69精品久久久久777片| 啦啦啦韩国在线观看视频| 99久久成人亚洲精品观看| 久久久欧美国产精品| 不卡一级毛片| 国产成人91sexporn| 精品免费久久久久久久清纯| 在线免费观看不下载黄p国产| 老女人水多毛片| 少妇熟女aⅴ在线视频| 亚洲激情五月婷婷啪啪| 亚洲国产欧美在线一区| 久久午夜福利片| 亚洲久久久久久中文字幕| 欧美人与善性xxx| 国产三级中文精品| 天堂√8在线中文| 亚洲va在线va天堂va国产| 精品久久久久久成人av| 国产免费男女视频| 成人毛片a级毛片在线播放| 丰满的人妻完整版| 午夜亚洲福利在线播放| 免费在线观看成人毛片| 精品久久久噜噜| 欧美最新免费一区二区三区| 美女大奶头视频| 日韩强制内射视频| 在线免费十八禁| 久久人妻av系列| 国产精品一区二区在线观看99 | 日本五十路高清| 中文在线观看免费www的网站| 欧美日韩国产亚洲二区| 日韩人妻高清精品专区| 亚洲国产欧美在线一区| 日韩高清综合在线| 成年av动漫网址| 日韩国内少妇激情av| 久久人人精品亚洲av| 乱人视频在线观看| 秋霞在线观看毛片| 国产伦精品一区二区三区四那| 午夜福利在线观看免费完整高清在 | 久久久精品94久久精品| 亚洲国产精品成人久久小说 | 天堂√8在线中文| 午夜福利视频1000在线观看| 国产乱人偷精品视频| 欧美不卡视频在线免费观看| 精华霜和精华液先用哪个| 欧美性感艳星| 免费在线观看成人毛片| 麻豆一二三区av精品| 1000部很黄的大片| 亚洲精品色激情综合| 99riav亚洲国产免费| 久久精品国产亚洲av天美| 少妇熟女aⅴ在线视频| 深夜a级毛片| 乱人视频在线观看| 99热全是精品| 国产午夜精品一二区理论片| 12—13女人毛片做爰片一| 久久久久久久久久黄片| 欧美成人一区二区免费高清观看| 亚洲av电影不卡..在线观看| 少妇熟女aⅴ在线视频| 免费观看a级毛片全部| 天美传媒精品一区二区| 免费在线观看成人毛片| 国产一区亚洲一区在线观看| 综合色丁香网| 一级av片app| 欧美激情国产日韩精品一区| 国产精品福利在线免费观看| 国产精品日韩av在线免费观看| 麻豆成人午夜福利视频| 日韩大尺度精品在线看网址| 久久精品91蜜桃| 在线免费观看的www视频| 午夜免费激情av| 国产高清三级在线| 亚洲国产高清在线一区二区三| 久久中文看片网| 国产成人影院久久av| 精品久久久久久久久久久久久| 一个人看的www免费观看视频| 午夜a级毛片| 日韩av在线大香蕉| 免费看日本二区| 熟妇人妻久久中文字幕3abv| 久久国产乱子免费精品| 校园人妻丝袜中文字幕| 一级毛片我不卡| 中文字幕制服av| 岛国毛片在线播放| 精品人妻偷拍中文字幕| 国产高清激情床上av| 国产精品人妻久久久影院| 深爱激情五月婷婷| av卡一久久| 热99re8久久精品国产| 蜜桃亚洲精品一区二区三区| 亚洲图色成人| 国产亚洲精品av在线| 夜夜看夜夜爽夜夜摸| 中文欧美无线码| 日韩欧美三级三区| 搡老妇女老女人老熟妇| 免费大片18禁| 卡戴珊不雅视频在线播放| 老司机影院成人| 亚洲欧美中文字幕日韩二区| 性欧美人与动物交配| 国内精品久久久久精免费| 黄色配什么色好看| 午夜福利在线观看免费完整高清在 | 尤物成人国产欧美一区二区三区| 日韩大尺度精品在线看网址| 一个人看视频在线观看www免费| 高清午夜精品一区二区三区 | 99久久精品热视频| 噜噜噜噜噜久久久久久91| 久久国内精品自在自线图片| 干丝袜人妻中文字幕| 国产色婷婷99| 欧美激情久久久久久爽电影| 久久欧美精品欧美久久欧美| 国产精品国产三级国产av玫瑰| 一级毛片久久久久久久久女| 国产淫片久久久久久久久| 欧美另类亚洲清纯唯美| 精品久久久久久久人妻蜜臀av| 2022亚洲国产成人精品| 97超视频在线观看视频| 欧美在线一区亚洲| 级片在线观看| 成人毛片60女人毛片免费| a级毛片a级免费在线| 亚洲成a人片在线一区二区| 不卡视频在线观看欧美| 国产伦精品一区二区三区视频9| 中国美女看黄片| 日韩大尺度精品在线看网址| 亚洲精品自拍成人| 日日干狠狠操夜夜爽| 午夜精品国产一区二区电影 | 欧美高清性xxxxhd video| 网址你懂的国产日韩在线| 丰满的人妻完整版| 精品欧美国产一区二区三| 成人一区二区视频在线观看| 国产极品天堂在线| 亚洲av熟女| 久久久久久久久久久丰满| 久久精品91蜜桃| 久久久久久久午夜电影| 色视频www国产| 网址你懂的国产日韩在线| 亚洲av免费在线观看| 能在线免费观看的黄片| 国产伦在线观看视频一区| 日韩精品有码人妻一区| 18禁裸乳无遮挡免费网站照片| 午夜精品一区二区三区免费看| 此物有八面人人有两片| 91在线精品国自产拍蜜月| 天堂影院成人在线观看| 国产精品综合久久久久久久免费| 男女边吃奶边做爰视频| 久久久国产成人精品二区| 亚洲经典国产精华液单| 亚洲精品色激情综合| 狂野欧美白嫩少妇大欣赏| 两性午夜刺激爽爽歪歪视频在线观看| 九九久久精品国产亚洲av麻豆| 好男人在线观看高清免费视频| 午夜免费男女啪啪视频观看| 男人舔女人下体高潮全视频| 日韩高清综合在线| 99久久精品一区二区三区| 精品人妻视频免费看| 少妇猛男粗大的猛烈进出视频 | 熟妇人妻久久中文字幕3abv| 欧美一区二区精品小视频在线| 成人漫画全彩无遮挡| 国产91av在线免费观看| 老熟妇乱子伦视频在线观看| 国产亚洲精品久久久久久毛片| 我的女老师完整版在线观看| 久久久国产成人免费| 亚洲精品日韩在线中文字幕 | 国产色婷婷99| 欧美又色又爽又黄视频| 岛国在线免费视频观看| 精品久久久久久成人av| 国产精品久久久久久久电影| 中文字幕熟女人妻在线| 国产精品久久电影中文字幕| 欧美区成人在线视频| 麻豆国产97在线/欧美| 特大巨黑吊av在线直播| 黄色日韩在线| 亚洲一级一片aⅴ在线观看| 在线免费十八禁| 国产免费男女视频| 搡女人真爽免费视频火全软件| 99久久精品国产国产毛片| 日韩成人伦理影院| 久久精品久久久久久噜噜老黄 | 一边亲一边摸免费视频| 男人狂女人下面高潮的视频| 亚洲欧美成人精品一区二区| av在线观看视频网站免费| 国国产精品蜜臀av免费| 亚洲成a人片在线一区二区| 成人高潮视频无遮挡免费网站| 丰满人妻一区二区三区视频av| 国产中年淑女户外野战色| h日本视频在线播放| 哪里可以看免费的av片| 欧美成人免费av一区二区三区| 插逼视频在线观看| 成人漫画全彩无遮挡| 一本精品99久久精品77| 亚洲av中文av极速乱| 亚洲婷婷狠狠爱综合网| 精品久久久久久久人妻蜜臀av| 午夜a级毛片| 久久精品夜色国产| av福利片在线观看| 国产一区二区三区在线臀色熟女| 亚洲国产精品久久男人天堂| 大又大粗又爽又黄少妇毛片口| 午夜福利在线观看免费完整高清在 | 国产日本99.免费观看| 一级黄色大片毛片| 免费大片18禁| 大又大粗又爽又黄少妇毛片口| 欧美bdsm另类| 少妇人妻精品综合一区二区 | 国产成人福利小说| 国产精品久久电影中文字幕| 亚洲av中文av极速乱| 亚洲欧美清纯卡通| 免费看光身美女| 长腿黑丝高跟| 一卡2卡三卡四卡精品乱码亚洲| 两个人的视频大全免费| 亚洲av熟女| 国产成人午夜福利电影在线观看| 99久久中文字幕三级久久日本| 成人二区视频| 日本免费一区二区三区高清不卡| 欧美xxxx黑人xx丫x性爽| 国产伦理片在线播放av一区 | 中文字幕精品亚洲无线码一区| 婷婷精品国产亚洲av| 黄色视频,在线免费观看| 午夜精品国产一区二区电影 | 久久精品夜夜夜夜夜久久蜜豆| 在线免费十八禁| 女同久久另类99精品国产91| 一夜夜www| 小蜜桃在线观看免费完整版高清| 高清日韩中文字幕在线| 能在线免费观看的黄片| 小说图片视频综合网站| 欧美一区二区亚洲| 免费人成在线观看视频色| 国产成人91sexporn| 干丝袜人妻中文字幕| 国产精品综合久久久久久久免费| 在线观看午夜福利视频| 看免费成人av毛片| 国产精品爽爽va在线观看网站| 国产在线精品亚洲第一网站| 精华霜和精华液先用哪个| 国产成人aa在线观看| 亚洲人成网站在线播放欧美日韩| 一夜夜www| 国产精品人妻久久久久久| 国产一区亚洲一区在线观看| 蜜桃久久精品国产亚洲av| 欧美日本亚洲视频在线播放| 精品一区二区三区视频在线| 夫妻性生交免费视频一级片| 深爱激情五月婷婷| 国产欧美日韩精品一区二区| 免费av毛片视频| 免费看日本二区| 99国产极品粉嫩在线观看| 成人av在线播放网站| 日本成人三级电影网站| 国产高清有码在线观看视频| 97超碰精品成人国产| 国产男人的电影天堂91| 亚洲不卡免费看| 久久国产乱子免费精品| 少妇人妻精品综合一区二区 | 我的老师免费观看完整版| 嫩草影院新地址| 国产成人午夜福利电影在线观看| 亚洲熟妇中文字幕五十中出| 国产综合懂色| 深夜a级毛片| 三级毛片av免费| av在线蜜桃| 午夜免费激情av| 亚洲av二区三区四区| 亚洲av第一区精品v没综合| 国产伦精品一区二区三区四那| 亚洲一级一片aⅴ在线观看| 久久精品国产亚洲av天美| 五月伊人婷婷丁香| 久久精品久久久久久久性| 日韩人妻高清精品专区| 久久九九热精品免费| 一个人看的www免费观看视频| 中国国产av一级| 日韩一区二区视频免费看| 亚洲av成人av| 国产精品永久免费网站| 亚洲自拍偷在线| 中文字幕精品亚洲无线码一区| 亚洲中文字幕一区二区三区有码在线看| 国产成人a∨麻豆精品| 国产精品乱码一区二三区的特点| 偷拍熟女少妇极品色| 可以在线观看毛片的网站| 国产精品,欧美在线| 日本在线视频免费播放| 五月伊人婷婷丁香| 久久久a久久爽久久v久久| 亚洲美女视频黄频| 亚洲av不卡在线观看| 久久国产乱子免费精品| 青春草国产在线视频 | 一级二级三级毛片免费看| 中文字幕av在线有码专区| 国产精品日韩av在线免费观看| 亚洲三级黄色毛片| 亚州av有码| 十八禁国产超污无遮挡网站| 色吧在线观看| 十八禁国产超污无遮挡网站| 色综合亚洲欧美另类图片| 九九爱精品视频在线观看| 日韩精品有码人妻一区| 亚洲av不卡在线观看| 亚洲精品日韩av片在线观看| 99热这里只有是精品在线观看| 欧美又色又爽又黄视频| 在线免费观看不下载黄p国产| 成人毛片a级毛片在线播放| 简卡轻食公司| 久久热精品热| 狂野欧美白嫩少妇大欣赏| 男女啪啪激烈高潮av片| 最近2019中文字幕mv第一页| 一进一出抽搐gif免费好疼| 精品欧美国产一区二区三| 成年女人看的毛片在线观看| 一边摸一边抽搐一进一小说| 国产高清不卡午夜福利| 一区二区三区免费毛片| www.av在线官网国产| av在线播放精品| 国产精品国产高清国产av| 亚洲aⅴ乱码一区二区在线播放| www.av在线官网国产| 日韩av在线大香蕉| 网址你懂的国产日韩在线| 少妇猛男粗大的猛烈进出视频 | 国产精品无大码| 国产精品伦人一区二区| 亚洲精品日韩av片在线观看| 成人鲁丝片一二三区免费| 国产一区二区三区av在线 | 自拍偷自拍亚洲精品老妇| 国产成人aa在线观看| 久久久午夜欧美精品| 久久99热6这里只有精品| 波多野结衣高清无吗| 久久人人爽人人片av| 99国产极品粉嫩在线观看| 久久久国产成人精品二区| 哪个播放器可以免费观看大片| 国产精品久久久久久久久免| www日本黄色视频网| 午夜精品国产一区二区电影 | 免费av毛片视频| 观看免费一级毛片| 九色成人免费人妻av| 精品久久久久久久人妻蜜臀av| 欧美性猛交黑人性爽| 日韩成人伦理影院| 日本免费a在线| 国产午夜精品论理片| 女人十人毛片免费观看3o分钟| 色综合色国产| 嫩草影院入口| 成年女人看的毛片在线观看| 精品一区二区三区视频在线| 日韩人妻高清精品专区| 麻豆乱淫一区二区| 国产 一区精品| 成人亚洲精品av一区二区| 久久这里只有精品中国| 国产大屁股一区二区在线视频| 国产精品人妻久久久影院| 精品国内亚洲2022精品成人| 成人特级av手机在线观看| 中国美白少妇内射xxxbb| 搡老妇女老女人老熟妇| 老熟妇乱子伦视频在线观看| videossex国产| 寂寞人妻少妇视频99o| 最近2019中文字幕mv第一页| 久久99精品国语久久久| 亚洲欧美精品综合久久99| 国产高清激情床上av| 日本一本二区三区精品| 精品久久久久久久久av| 亚洲精华国产精华液的使用体验 | 国产伦一二天堂av在线观看| 最近中文字幕高清免费大全6| 网址你懂的国产日韩在线| 男插女下体视频免费在线播放| 一进一出抽搐gif免费好疼| 成年免费大片在线观看| 欧美bdsm另类| av在线老鸭窝| 91久久精品国产一区二区三区| 亚洲在久久综合| videossex国产| 99久久久亚洲精品蜜臀av| 青春草国产在线视频 | 午夜福利视频1000在线观看| 两个人的视频大全免费| 韩国av在线不卡| 26uuu在线亚洲综合色| 国产色爽女视频免费观看| 一级毛片电影观看 | 国产精品电影一区二区三区| 黄色视频,在线免费观看| 国产精品国产高清国产av| 性色avwww在线观看| 日韩人妻高清精品专区| 国产精品伦人一区二区| 久久精品国产清高在天天线| 可以在线观看毛片的网站| 欧美性感艳星| 人人妻人人澡人人爽人人夜夜 | 两个人的视频大全免费| 国产日本99.免费观看| 婷婷色综合大香蕉| av在线老鸭窝| 中文字幕熟女人妻在线| 免费无遮挡裸体视频| 99久久精品国产国产毛片| 2022亚洲国产成人精品| 精品99又大又爽又粗少妇毛片| 免费不卡的大黄色大毛片视频在线观看 | 丰满人妻一区二区三区视频av| a级一级毛片免费在线观看| 综合色丁香网| 国产高清不卡午夜福利| 日韩 亚洲 欧美在线| 国产一区二区在线观看日韩| av在线老鸭窝| 成人漫画全彩无遮挡| 国产高清不卡午夜福利|