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

    BCC-ESM1 Model Datasets for the CMIP6 Aerosol Chemistry Model Intercomparison Project (AerChemMIP)

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

    Jie ZHANG, Tongwen WU, Fang ZHANG, Kalli FURTADO, Xiaoge XIN, Xueli SHI, Jianglong LI,Min CHU, Li ZHANG, Qianxia LIU, Jinghui Yan, Min WEI, and Qiang MA

    1Beijing Climate Center, China Meteorological Administration, Beijing 100081, China

    2Met Office Hadley Centre, Exeter EX1 3PB, UK

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

    ABSTRACT

    Key words:BCC-ESM1,CMIP6,AerChemMIP,climate sensitivity,precipitation,surface air temperature,sulfate

    1.Introduction

    Climate impacts of atmospheric chemistry, especially the contribution of tropospheric ozone precursors to radiative forcing, have been estimated in all the Intergovernmental Panel on Climate Change (IPCC) reports. However,even the models used in the IPCC’s fourth report were mostly offline chemistry transport models, and none of the climate models included tropospheric ozone chemistry(Collins et al., 2017). The Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) was launched in 2011 and aimed to document changes in atmospheric composition and the associated radiative forcings in CMIP5 models (Lamarque et al., 2013). However, the model setups for ACCMIP were generally different from CMIP5 (with lower resolution but more complex chemistry)and only four of the models could be used to estimate future climate changes in surface ozone (Schnell et al., 2016). To build on the experiences of previous programs, the CMIP6-endorsed Aerosol Chemistry Model Intercomparison Project(AerChemMIP) is a series experiments with the goal of quantifying the effects of near-term climate forcers (NTCFs: methane, tropospheric aerosols and ozone, and their precursors)and reactive well-mixed greenhouse gases (WMGHGs) on climate and air quality (Collins et al., 2017). AerChemMIP is designed to answer four scientific questions: (1) how have anthropogenic emissions contributed to global radiative forcing and affected regional climate over the historical period? (2) How might future policies (on climate, air quality, and land use) affect the abundances of NTCFs and their climate impacts? (3) How do uncertainties in historical NTCF emissions affect radiative forcing estimates? (4) How important are climate feedbacks to natural NTCF emissions,atmospheric composition, and radiative effects?

    BCC-ESM1 is a new version of the Beijing Climate Center’s ESM that includes interactive atmospheric chemistry and aerosols. Its details are described in Wu et al. (2020a).BCC-ESM1, along with 14 other ESMs, have participated in AerChemMIP, including CESM2-WACCM, CNRMESM2-1, EC-Earth3-AerChem, EMAC-2-54-AerChem,GFDL-ESM4, GISS-E2.1, IPSL-CM5A2-INCA, IPSLCM6A-AER-LR, MIROC6, MIROC-ES2L, MPI-ESM-1-2-HAM, MRI-ESM2-0, NorESM2, and UKESM1-0-LL. All of these ESMs are run in AerChemMIP and the other programs under CMIP6 [DECK (Diagnostic, Evaluation and Characterization of Klima) and the other CMIP6-endorsed MIPs] with their maximum available complexity. Therefore, their results can be shared across MIPs and used for more comprehensive assessments. The AerChemMIP experiments are classified into three tiers by their priority. The Aer-ChemMIP is the first international MIP that BCC-ESM1 has been involved in, and thus far it has conducted all its Tier 1 experiments and seven experiments in lower priorities. BCC-ESM1 has also run the DECK and CMIP historical simulations, which serve as the “entry cards” for models participating in CMIP6 (Eyring et al., 2016).

    AerChemMIP is the first chemistry program fully involved in the CMIP structure. Most of its experiments are newly designed and of special scientific interest. This paper introduces the experiment settings as well as the BCCESM1 outputs, to facilitate their use. Following this introduction, a brief overview of BCC-ESM1 is provided in section 2. The experiments run by BCC-ESM1 are also summarized in that section. Section 3 presents a technical validation of BCC-ESM1, including the model climate sensitivity and evaluations of surface air temperature (SAT), precipitation, and surface sulfate concentrations. Section 4 provides a summary, and section 5 shows some usage notes.

    2.Model and experiments

    2.1.Model configuration

    BCC-ESM1 is based on the climate system model BCC-CSM2 (Wu et al., 2019a), but includes interactive representations of tropospheric chemistry and aerosols. As the successor of BCC-CSM1.1 (Wu et al., 2013), which participated in CMIP5, BCC-CSM2 included updates to the deep convection parameterization, considered the indirect effects of aerosols and gravity wave drag generated by deep convection, and introduced a new scheme for the cloud fraction.BCC-CSM2-MR, the medium resolution version of BCCCSM2, shows an overall improvement in many aspects,such as the tropospheric air temperature at the global and regional scale in East Asia, the stratospheric quasi-biennial oscillation, and the Madden–Julian Oscillation.

    The atmospheric component of BCC-ESM1 is BCCAGCM3-Chem (Wu et al., 2020a), which includes an interactive tropospheric gas-phase and aerosol chemistry scheme.The model includes 66 gas-phase species and 13 prognostic aerosol species. The main atmospheric chemistry built in BCC-AGCM3-chem is based on MOZART2 (Model for OZone And Related chemical Tracers, version 2), a troposphere-only chemistry transport model developed by the National Center for Atmospheric Research (Horowitz et al.,2003). SO, dimethyl sulfide (DMS), and NHare the gas tracers newly added to BCC-ESM1. A bulk aerosol scheme is also newly introduced in BCC-ESM1, including the relative reactions, transport, emissions, and depositions. The aerosol species involved are sulfate (SO), hydrophobic and hydrophilic organic carbon (OC1 and OC2), hydrophobic and hydrophilic black carbon (BC1 and BC2), four categories of soil dust (DST01, DST02, DST03, DST04), and four categories of sea salt (SSLT01, SSLT02, SSLT03, SSLT04).Dust and sea salt aerosols are classified into four size bins by diameter. The mass-mixing ratios of each aerosol type are prognostic variables. Prognostic aerosol masses are used to estimate the liquid cloud droplet number concentration and the indirect effects of aerosols are considered.

    To save computing resources and tune the model more efficiently, BCC-AGCM3-Chem has a relatively low horizontal resolution of T42 (about 280 km) and has 26 vertical levels with the top level at 2.914 hPa. The land component of BCC-ESM1 is the BCC Atmosphere and Vegetation Interaction Model, BCC-AVIM2.0, with terrestrial carbon cycle(Li et al., 2019). The oceanic component is the Modular Ocean Model, version 4, with 40 levels included (MOM4-L40; Griffies et al., 2005). The sea-ice component is the Sea Ice Simulator (SIS), with three vertical layers, one snow cover layer, and two ice layers (Winton, 2000). MOM4-L40 and SIS use a tripolar grid with a horizontal resolution of 1°longitude by 1/3° latitude between 10°S and 10°N, increasing to 1° at 30°S and 30°N.

    2.2.Experiments

    The experiments conducted by BCC-ESM1 are summarized in Table 1. The simulation data have been archived on the Earth System Grid Federation (ESGF) node (https://esgf-node.llnl.gov/projects/cmip6/). DECK includes four basic experiments: (a) a historical Atmospheric Model Intercomparison Project (amip) simulation; (b) a control simulation under pre-industrial conditions (piControl); (c) a simulation that branches from piControl in 1850 and is forced by an abruptly quadrupled global-mean concentration of CO(abrupt-4×CO2); and (d) a simulation that branches from piControl and is forced by a gradually increased COconcentration at a rate of 1% per year (1pctCO2). Amip is run by an atmospheric general circulation model (AGCM) with prescribed sea surface temperature (SST) and sea-ice concentration (SIC) and include three ensemble members. The other three experiments are all conducted by a coupled atmosphere–ocean general circulation model (AOGCM) and include only one ensemble member. The CMIP historical experiment (historical) is also conducted by the AOGCM. It branches from piControl and is forced by evolving external forcing. There are three ensemble members for historical with slightly varying initial conditions. Amip, piControl,and historical are baselines for many CMIP6-Endorsed MIPs. Abrupt-4×CO2 and 1pctCO2 are used to estimate model climate sensitivities and feedbacks.

    Table 1. Experiments conducted by BCC-ESM1.

    The AerChemMIP experiments are run with the same complexity as in DECK and historical. The experiments are classified into three tiers by their priority. The Tier 1 experiments mainly focus on the impacts of NTCFs and reactive well-mixed gases on historical and future climate, as well as future NTCFs under different air quality policies. Tier 2 and 3 will examine the effects of NTCF emissions and climate feedback more specifically, to provide detailed uncertainty estimates.

    BCC-ESM1 has run all the Tier 1 experiments and seven Tier 2 and 3 experiments:

    ● hist-piNTCF in Tier 1 parallels historical but with the anthropogenic emissions of aerosols, and ozone precursors fixed to their 1850 values. hist-piAer in Tier 3 is similar to hist-piNTCF but with the anthropogenic emissions of aerosol fixed to the 1850 values.By comparing with historical, these experiments can be used to estimate the climate impacts of anthropogenic emissions of NTCFs and aerosols. Three ensemble members are conducted for both histpiNTCF and hist-piAer with the same initial conditions as historical.

    ● piClim-control in Tier 1 is run under pre-industrial conditions for 31 years with prescribed SST and SIC from 30 years of piControl climatology. It is designed following the protocols of the Radiative Forcing Model Intercomparison Project (Pincus et al.,2016). Eight more piClim-X experiments are also conducted by BCC-ESM1 that branch from piClim-control: (a) piClim-NTCF and piClim-CH4 in Tier 1,which are the same as piClim-control but set NTCF emissions and the CHconcentration to present-day values, respectively; (b) piClim-aer, piClim-BC and piClim-O3 in Tier 2, which set the emissions of aerosol precursors, black carbon and ozone precursors to present-day values, respectively; (c) piClim-NOX,piClim-VOC and piClim-SO2 in Tier 3, which set the NOx, CO and VOCs, and SOemissions to presentday values, respectively. These experiments can help to quantify the effective radiative forcings (ERFs) at pre-industrial levels due to different drivers, as in Thornhill et al. (2020). The piClim-X experiments include only one ensemble per experiment in BCCESM1.

    ● histSST in Tier 1 also branches from piClim-control but with prescribed time-evolving SST and SIC from historical. All other forcings are the same as in historical. histSST-piNTCF and histSST-piCH4 in Tier 1 are the same as histSST and branch from the same point as histSST but with the NTCF precursors and methane fixed to their 1850 values, respectively.These experiments can help to diagnose the transient ERFs for NTCF and methane. The histSST-X experiments include only one ensemble per experiment.

    ● ssp370 in Tier 1 is a future simulation with the highest levels of short-lived climate pollutants without climate policy. It achieves forcing levels of 7.0 W min 2100. ssp370-lowNTCF in Tier 1 is a future simulation with strong levels of air quality control measures. Both experiments branch from the historical experiments in 2015 and are run for 50 years from 2015 to 2055 instead of 85 years from 2015 to 2100 as suggested by AerChemMIP, to reduce the computational expense. These experiments can help to estimate the impacts on future climate and air quality, as in Allen et al. (2020) and Turnock et al.(2020). Three ensemble members are conducted for both ssp370 and ssp370-lowNTCF.

    ● ssp370SST in Tier 1 is a transient future prescribed-SST simulation following the ssp370 scenario. It branches from the first member of amip in year 2015 with prescribed SST and SIC from the first member of the ssp370 experiment. ssp370SST-lowNTCF in Tier 1 branches at the same point as ssp370SST with prescribed SST and SIC from the first member of ssp370-lowNTCF following the ssp370-lowNTCF scenario. Different ERF changes between these future prescribed SST simulations and the relevant coupled model projections (ssp370 and ssp370-lowNTCF) help to quantify the efficacy of NTCFs to affect climate. Both experiments include only one ensemble.

    3.Technical validation

    3.1.Model climate sensitivities to idealized CO2 forcing

    Climate sensitivity is an important metric for model evaluation. It is one of the major factors responsible for different model climate responses and uncertainty in future climate projections (Grose et al., 2018). Here, we estimate the global mean SAT response to idealized COforcing by using the 1pctCO2 and abrupt-4×CO2 experiments. Both experiments are initialized from the piControl experiment,which is used as the reference experiment.

    The transient climate response (TCR) is calculated from the 1pctCO2 simulation in which COis gradually increased at a rate of 1% per year. The TCR is defined as the annual mean SAT change relative to the piControl equilibrium climate state, averaged over a 20-year period centered on the time at which the COconcentration was doubled(year 70 of the simulation). The 1pctCO2 forcing resembles the COforcing in the 21st century but is more idealized.As shown in Fig.1a, the evolution of SAT anomaly is linearly correlated with the COconcentration in the 1pctCO2 experiment. A high linear correlation is also evident in observations at both global and regional scales (Wu et al.,2019b). The TCR in BCC-ESM1 is about 1.78 K, which is 1.40 K for BCC-CSM2-MR (Shi et al., 2020). The TCRs in both BCC models lie within the “l(fā)ikely” range of 1.0–2.5 K suggested by IPCC AR5 (Bindoff et al., 2013). In the 1pctCO2 experiment and the real world, the distribution of heat flux at the atmosphere–ocean interface has not yet reached equilibrium.

    The equilibrium climate sensitivity (ECS) is measured as the SAT changes in a double COexperiment when the climate system reaches its equilibrium climate state. The ECS in BCC-ESM1 is estimated by the abrupt-4xCO2 experiment. SAT changes in the abrupt-4xCO2 experiment are shown in Fig.1b. The model response to quadrupled COcan be separated into a fast-response stage in the first two decades and a slow-response stage afterward. Following the approach of Gregory et al. (2004), we regress the global annual mean net radiative-flux changes at the top of the atmosphere (TOA) on SAT changes and then extrapolate to zero net flux (Fig.2). The ECS is then defined as half of the xintercept concerning a COdoubling. This method yields an ECS estimate of 3.34 K for BCC-ESM1, which is 3.03 K for BCC-CSM2-MR (Wu et al., 2019a). ECSs in both BCCESM1 and BCC-CSM2-MR are within the “l(fā)ikely” range of 1.5–4.5 K suggested by IPCC AR5.

    3.2.SAT in piControl and historical experiments

    Fig.1. (a) Temporal evolution of globally averaged SAT changes relative to the reference experiment (piControl) for the 151 model years in the 1%CO2 experiment. (b) As in (a)but for SAT changes in the abrupt-4×CO2 experiment. Units: K.

    Fig.2. Annual changes of the net TOA radiative flux and the global-mean SAT for the first 151 years of the abrupt CO2 quadrupling experiment (abrupt-4×CO2) relative to the preindustrial control run (piControl). The regression line is also calculated to estimate the ECS of BCC-ESM1.

    Figure 3a shows the evolution of the global mean annual averaged SAT in the piControl run. In piControl,SAT reaches a quasi-equilibrium state with a slight increasing trend of about 0.016°C (100 yr), which is much smaller than the standard deviation of SAT (0.10°C yr). Time series of annual mean SAT anomalies from the three members of the CMIP6 historical simulations are shown in Fig.3b, as well as the observation from HadCRUT4(https://www.metoffice.gov.uk/hadobs/hadcrut4/). The anomalies are relative to the 1850–1950 mean. The SAT simulations show good agreement with HadCRUT in the first 100 years. However, there is a significant cooling bias in the second half of the 20th century, especially over the period from the 1960s to 1990s. We colloquially refer to this significant cooling bias as the “pot-hole” bias (because it represents an undesirable dip in the “road” of an accurately simulated SAT trend). It is a common problem in most of the ESMs participating in CMIP6. One possible reason for the pot-hole bias is that the prognostic aerosol schemes produce overly strong aerosol cooling (Flynn and Mauritsen, 2020; Meehl et al., 2020).

    The spatial structure of SAT in the historical experiment is examined over the period from 1986 to 2005 (Fig.4),as well as its biases compared with ERA5 reanalysis (Copernicus Climate Change Service, 2017). The spatial structure of SAT is well captured with a spatial correlation of 0.99. Simulated SAT is about 0.35 K higher than that in ERA5, attributable to the warmer tropics and subtropics. Generally, the biases in BCC-ESM1 show great similarity with those in BCC-CSM2-MR (Wu et al., 2019a), including the slight systematic warm biases over oceans, as well as the relatively large warming/cooling biases over the polar regions. The biases over the oceans and polar regions may partly come from the ocean and sea-ice modules, which are the same in BCC-ESM1 and BCC-CSM2-MR. Both models also show systematic cool biases over most land regions. The cooling bias may be partly attributable to the land surface modeling component, since it is also evident in AMIP simulations, or to other processes like errors in cloud properties. The cooling biases in BCC-ESM1 may also be related to the overestimated sulfate concentration in the atmosphere shown in subsection 3.4.

    Fig.3. (a) Time series of the global mean SAT in the piControl simulation. (b) Time series of the global mean SAT anomalies in the historical simulations (three members, black lines) and the observation (red line) from HadCRUT4. The anomalies are relative to the 1850–1950 mean.

    3.3.Precipitation and associated near-surface horizontal wind

    Annual mean precipitation in the historical experiment and its comparison with GPCP observation (Adler et al.,2003) are examined and shown in Fig.5. The biases are generally within 1 mm dand the global mean precipitation in BCC-ESM1 is very close to the observation: 2.27 mm dversus 2.24 mm d. BCC-ESM1 still suffers from the double ITCZ problem: precipitation south of the equator in the eastern Pacific is overestimated and the SPCZ is too zonally elongated. Zhang et al. (2015) found that the double ITCZ related biases include sufficient precipitation in the southeastern Pacific, warmer SST, weaker easterly flow,and stronger meridional wind divergences away from the equator. However, no progress can be identified from the multimodel ensemble (MME) mean in CMIP3 to that in CMIP5.The double ITCZ measures are also examined for BCC-ESM1, as in Zhang et al. (2015). The SST bias in BCC-ESM1 closely resembles the SAT bias (not shown).As indicated by Fig.4c, the warmer SST south of the equator amplifies the precipitation over the Southeast Pacific and leads to the eastward-extended SPCZ. Over the eastern tropical Pacific, the near-surface zonal wind weakens and the meridional wind converges south of the equator(Fig.6). Both zonal and meridional wind biases contribute to the SPCZ biases in BCC-ESM1. Generally, biases of the double ITCZ measures in BCC-ESM1 are quite similar to those in the CMIP3 and CMIP5 MMEs. With great improvements in atmospheric physical parameterizations, the double ITCZ bias is much weaker in the newly developed high-resolution climate system model, BCC-CSM3-HR (Wu, 2020b).

    Fig.4. Annual mean SAT in (a) BCC-ESM1 and (b) ERA5 reanalysis from 1986 to 2005. (c) The SAT biases. Units: °C.

    Other model deficiencies in precipitation reproductions are evident, including the overestimations over eastern Indonesia, western Indian Ocean, southern Africa, and East Asia. Meanwhile, underestimations are pronounced north of the equator in the Pacific and Atlantic. Such biases seem to be common in CMIP3 and CMIP5 models (Oueslati and Bellon, 2015). The sufficient precipitation over eastern Indonesia and the western Indian Ocean accompanies stronger trade winds in the east. The anomalous easterlies converge toward the biased precipitation, west of the Maritime Continent and the eastern African coast. The underestimations north of the equator in the Pacific and Atlantic are asymmetric with the overestimations in the south, which may be attributable to the asymmetric warming amplitudes by the equator as well as the local meridional wind divergences(Zhang et al., 2015). The Southern Africa and East Asia are greatly affected by the global monsoon system, and the precipitation biases can be related to many physical and dynamic processes, but investigating these is beyond the scope of the present study.

    Fig.5. Annual mean precipitation in (a) BCC-ESM1 and (b)GPCP from 1986 to 2005. (c) The precipitation biases. Units:mm d?1.

    3.4.Surface sulfate and its variation

    Sulfate is a major species of aerosol. Sulfate aerosols can cool the atmosphere by reflecting solar irradiation directly or by making clouds more reflective via aerosol indirect effects. The anthropogenic emissions of SOfrom fuel combustion and industrial activities constitute the main precursor of atmospheric sulfate. A large amount of sulfate in the atmosphere is believed to have significant cooling effects at regional and global scales, although large uncertainties remain in the magnitude of these effects. BCC-ESM1 can capture the major features of aerosols reasonably well in comparison with observed aerosol concentrations and their optical properties (Wu et al., 2020a). In this paper, we use ground-based monitoring networks to extend the evaluation of sulfate. Simulation from the historical experiment is examined here.

    The sulfate observations are from three monitoring networks: the European Monitoring and Evaluation Program(EMEP; T?rseth et al., 2012); the United States Interagency Monitoring of Protected Visual Environment (IMPROVE;Malm et al., 2004) extensive ground-based networks; and the Acid Deposition Monitoring Network in East Asia(EANET, 2019). For each observation site, model data are selected only when the monthly observation is available. To maximize the amount of data available for comparison but also ensure the climatological representativeness of the observation, all 12 monthly means are required to calculate an annual mean, as in Mulcahy et al. (2020).

    Fig.6. Annual mean near-surface zonal wind (UAS) in (a) BCC-ESM1 and (b) ERA5 from 1986 to 2005. (c) The UAS biases. (d–f) As in (a–c) but for the near-surface meridional wind (VAS). Units: m s?1.

    The annual mean surface sulfate concentrations from 274 monitoring sites are calculated and shown in Fig.7a. Surface sulfate concentrations at EMAP, EANET and IMPROVE sites are about 3.33 μg m, 3.20 μg mand 1.35 μg m, respectively. Generally, sulfate is higher in southern Europe than in northern Europe and higher in the east of the U.S. than in the west. BCC-ESM1 tends to overestimate the observations by about 0.95 μg mbut shows a high spatial correlation of 0.77, significant at the 1% level using a Student’s t-test (Fig.7b). The surface sulfate concentrations are comparable with the observation at the EMEP sites (3.18 μg m) but overestimated by approximately 2 μg mat the EANET and IMPROVE sites. BCC-ESM1 can capture the observed sulfate centers well over the eastern U.S., southern Europe, and Japan.

    Regional statistical evaluations are shown in Fig.8. Considering the relatively low model resolution and the high heterogeneous distribution of sulfate, simulations falling with the range of 50% to 200% of the observations (gray dashed lines in Fig.8) are considered to be reasonable. Over the EMEP sites (Europe), the model shows a significant correlation of 0.66, with a slight underestimation over the southern EMEP sites, where the measured sulfate concentration is relatively large. Over EANET sites (East Asia), the comparison shows a high degree of scattering, with a low correlation with the observation of 0.39. Over the western U.S., the sulfate concentration is about twice that of the observation.As one of the major emission centers, the eastern U.S. is more polluted than the west. Among the four regions, the simulations perform best over the eastern U.S., with the highest correlation coefficient with the observation (0.87) and the smallest root-mean-square error (RMSE; 1.04 μg m). Generally, the performance of BCC-ESM1 in reproducing the surface sulfate concentration is comparable to instrumental measurements, and its performance is comparable with that of UKESM1 (Mulcahy et al., 2020). The major discrepancy is that BCC-ESM1 tends to overestimate sulfate levels over the western U.S. and East Asia, which may partly contribute to the regional cooling biases as shown in Fig.4c.

    Fig.7. (a) Annual mean surface sulfate concentration in 2000–09 at 274 monitoring sites in the three ground-based networks: EMEP in Europe, IMPROVE in the U.S., and EANET in East Asia. (b) Corresponding simulation in BCCESM1. All 12 monthly means are required to calculate annual mean results. Units: μg m?3.

    Fig.8. Comparison of simulated annual mean SO42? from BCC-ESM1 against ground-based measurements from (a)EMEP (Europe), (b) EANET (East Asia), (c) West IMPROVE(western U.S.), and (d) East IMPROVE (eastern U.S.)networks. Observations and model data cover the years 1978 to 2009 for EMEP, 2000 to 2014 for EANET, and 1990 to 2014 for IMPROVE. NMB stands for normalized mean bias,RMSE for root-mean-square error, and COR for correlation coefficient. Linear regression statistics are also included.Units: μg m?3.

    Figure 9 shows the evolution of the simulated and observed annual mean sulfate concentration, averaged across all the available ground-based observations, for all available years in the four typical regions. The uncertainty is calculated as ±1 standard deviation of the observed annual means across all the measurement sites. Simulations are suggested to be comparable with observations if the evolutions of the mean series of all sites are similar and the bands of uncertainty ranges are well overlapped. Generally, despite the systematic overestimation in East Asia and the western U.S., BCC-ESM1 can reproduce the changes in sulfate concentration at regional scales reasonably well.

    Over Europe (EMEP), the sulfate concentration decreases by about 4 μg mfrom 1978 to 2009. BCC-ESM1 captures the reduction and the changes are well within the observed uncertainty range, albeit the negative trend is slightly smaller. This underestimation of the annual mean surface sulfate concentration over Europe is also evident in UKESM1 simulations (Mulcahy et al., 2020) and is suggested to be related to a low bias in wintertime in UKESM1(Turnock et al., 2015).

    Over EANET sites, the surface sulfate concentration increases slightly with large interannual variability and uncertainties (about 3–4 μg m). It increases in the early 2000s,decreases by about 1.5 μg mfrom 2006 to 2012, and then bounces back afterward. The model generally tracks the observation quite well, despite a systematic overestimation by about 2 μg m. The simulation also shows large variations across all the measurement locations. Note that the EANET sites here are mostly in Japan and the southern countries.China is also a main sulfate emissions center with large spatial heterogeneity. Local evaluation over China will be interesting when high coverage monitoring data are available.

    Fig.9. Time series of annual mean observed (red lines) and simulated (blue lines) sulfate concentrations averaged across all measurement locations in each network for a particular year. Shaded areas show ±1 standard deviation of the observed and modeled annual mean values across all the measurement locations. Units: μg m?3.

    Over the western U.S., the sulfate concentration is relatively small with little to no trend. The simulated sulfate concentration is about twice that of the observation with a small decreasing trend. This region is generally remote from source regions and differences may be attributable to the biases in transportation from the source region, oxidants, or removal rates (Mulcahy et al., 2020).

    Over the eastern U.S., one of the major anthropogenic source regions, a significant negative trend is observed due to reductions in emissions. Although the model generally overestimates concentrations by about 0.8 μg m, the decreasing rate and the uncertainty range in BCC-ESM1 show reasonable agreement with the observations.

    4.Summary and relevant studies

    This paper summarizes the experiments and datasets of BCC-ESM1, specifically for the CMIP6-endorsed AerChem-MIP. The basic performance of BCC-ESM1 is evaluated,revealing that: the climate sensitivities in BCC-ESM1 are well within the likely ranges suggested by IPCC AR5; SAT and precipitation can be reasonably represented (the biases of which are discussed); and the present-day surface sulfate concentrations also compare well with observations, including their spatial structure and regional trends.

    BCC-ESM1 is one of the few ESMs participating in CMIP6 with interactive atmospheric chemistry, and its Aer-ChemMIP outputs have been widely used in recent atmospheric chemistry studies (i.e., Allen et al., 2020; Griffiths et al., 2020; Thornhill et al., 2020; Turnock et al., 2020).For example, in the special issue for AerChemMIP, BCCESM1 data are used to estimate the climate and air quality impact of mitigation policies. Turnock et al. (2020) used the outputs from BCC-ESM1 and nine other ESMs to make a first assessment of the impact on historical and future changes to air pollutants (Oand PM). CMIP6 models show a consistent overestimation of surface Oconcentration but an underestimation of surface PMwith a large diversity in the historical period. For scenario ssp370,which encompasses the highest levels of short-lived climate pollutants without climate policy, both surface Oand PMincrease across most regions. In BCC-ESM1, surface Ochanges are about 3 ppbv higher than the multimodel mean value for the historical period and the surface Oprediction is large in scenario ssp370. The large Oresponse in BCC-ESM1 may be related to the absence of NOtitration,which will reduce surface Oconcentrations. Allen et al.(2020) quantified the 2015–2055 climate and air quality effects of non-methane NTCFs with “weak” and “strong”levels of air quality controls. Their findings suggest air quality improvement with non-methane NTCF reductions, but additional surface warming and wetting may occur, particularly in Asia and the Arctic. Generally, results in BCCESM1 are consistent with the seven other ESMs participating in CMIP6, both in terms of global mean climate anomalies and air pollution anomalies due to NTCF mitigation.

    5.Usage notes

    The variables of Priority 1 in AerChemMIP have been generated, including meteorological variables at various time scales (monthly, daily, six-hourly, and hourly) and chemical variable outputs at the monthly scale. All the variables are listed at http://clipc-services.ceda.ac.uk/dreq/tabs03/expt_AerChemMIP_AerChemMIP_1_1.html. For brevity,we summarize the chemical outputs in Table 2, including the volume mixing ratios of gas-phase chemical species in the atmosphere, the mass mixing ratios of aerosols, atmospheric optical thickness, etc.

    The horizontal resolution of all the variables is T42(approximately 280 km). Most of the 3D atmospheric outputs are on model atmosphere levels (26 levels in a hybrid sigma-pressure vertical coordinate system). To facilitate multimodel comparisons, 11 atmospheric variables are interpolated to 19 fixed pressure levels for monthly outputs, and 6 of them are interpolated to 8 fixed pressure levels for daily outputs (Table 3).

    Acknowledgements.

    This work was jointly supported by the National Key Research and Development Program of China (Grant No. 2016YFA0602103), the National Key Research and Development Program of China CERC-WET Project (Grant No.2018YFE0196000) and the National Natural Science Foundation of China (Grant No. 41805063).

    Table 2. Chemical outputs of BCC-ESM1 in CMIP6. Variables marked with an asterisk (*) are only provided for the first member of historical experiments.

    Table 2. (Continued.)

    Table 3. Monthly outputs of 11 variables that have been interpolated to 19 fixed pressure levels at 1000, 925, 850, 700, 600, 500, 400,300, 250, 200, 150, 100, 70, 50, 30, 20, 10, 5, and 1 hPa. Daily outputs of 6 variables have been interpolated to 8 fixed pressure levels at 1000, 850, 700, 500, 250, 100, 50, and 10 hPa.

    Data availability statement

    The data for the BCC-ESM1 simulations (Zhang et al.,2019) are available on the ESGF nodes (https://esgfnode.llnl.gov/projects/cmip6/). The dataset is in Network Common Data Form (NetCDF), version 4.

    Disclosure statement

    No potential conflicts of interest are reported by the authors.

    Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source,provide a link to the Creative Commons license, and indicate if changes were made.

    两性午夜刺激爽爽歪歪视频在线观看 | 天堂动漫精品| 丝袜美腿诱惑在线| 1024手机看黄色片| 日本三级黄在线观看| 国产精品野战在线观看| 成人18禁在线播放| 在线永久观看黄色视频| avwww免费| 99精品欧美一区二区三区四区| АⅤ资源中文在线天堂| 91成年电影在线观看| 男女之事视频高清在线观看| 久久性视频一级片| 午夜激情福利司机影院| 日韩有码中文字幕| 男插女下体视频免费在线播放| 国产av又大| 欧美zozozo另类| 久久久久国内视频| 在线观看日韩欧美| 欧美一级毛片孕妇| 亚洲欧洲精品一区二区精品久久久| 久久人妻av系列| 中文字幕av在线有码专区| 久久精品亚洲精品国产色婷小说| 国产精品免费视频内射| 麻豆久久精品国产亚洲av| 国产三级中文精品| 久久精品夜夜夜夜夜久久蜜豆 | 久久久水蜜桃国产精品网| 在线观看美女被高潮喷水网站 | 亚洲国产高清在线一区二区三| 亚洲美女视频黄频| 18美女黄网站色大片免费观看| 精品久久久久久成人av| 777久久人妻少妇嫩草av网站| 亚洲成人精品中文字幕电影| 51午夜福利影视在线观看| 成年免费大片在线观看| 久久久久免费精品人妻一区二区| 国产av在哪里看| 久久久久久国产a免费观看| 好男人在线观看高清免费视频| 亚洲第一电影网av| 国产一区二区三区视频了| 国产精品免费一区二区三区在线| 亚洲五月天丁香| 国产成人啪精品午夜网站| 亚洲最大成人中文| 国产一区二区在线观看日韩 | 无遮挡黄片免费观看| 欧美性猛交黑人性爽| 国产熟女午夜一区二区三区| 两个人看的免费小视频| 日本黄大片高清| 国产91精品成人一区二区三区| 色播亚洲综合网| 国产单亲对白刺激| 成年版毛片免费区| 搡老妇女老女人老熟妇| 欧美丝袜亚洲另类 | 欧美一区二区精品小视频在线| 亚洲成人国产一区在线观看| 俺也久久电影网| 正在播放国产对白刺激| 99国产精品99久久久久| 久久久久久九九精品二区国产 | 夜夜躁狠狠躁天天躁| 婷婷精品国产亚洲av| 国产高清激情床上av| 制服丝袜大香蕉在线| 欧美色欧美亚洲另类二区| 可以在线观看毛片的网站| 亚洲成人免费电影在线观看| 国产蜜桃级精品一区二区三区| 亚洲七黄色美女视频| 夜夜夜夜夜久久久久| 国产精品野战在线观看| 日韩欧美国产一区二区入口| 91成年电影在线观看| 午夜福利视频1000在线观看| 久久久久久久精品吃奶| 亚洲av成人精品一区久久| 制服诱惑二区| 久久精品成人免费网站| 国产av又大| 少妇的丰满在线观看| 久久午夜综合久久蜜桃| 99精品在免费线老司机午夜| 岛国在线观看网站| 国产激情欧美一区二区| av超薄肉色丝袜交足视频| 一个人免费在线观看的高清视频| 久久久国产精品麻豆| 国产精品爽爽va在线观看网站| 两个人视频免费观看高清| a级毛片在线看网站| 黄色视频不卡| 变态另类丝袜制服| 变态另类成人亚洲欧美熟女| 亚洲中文字幕一区二区三区有码在线看 | 在线观看免费视频日本深夜| 三级国产精品欧美在线观看 | 日韩欧美一区二区三区在线观看| www.999成人在线观看| 亚洲九九香蕉| 欧美成人一区二区免费高清观看 | 国产69精品久久久久777片 | 欧美最黄视频在线播放免费| 又黄又爽又免费观看的视频| 午夜视频精品福利| 久久天堂一区二区三区四区| 成人国产一区最新在线观看| 欧美性猛交╳xxx乱大交人| 在线看三级毛片| 日日爽夜夜爽网站| 午夜免费成人在线视频| 欧美黑人巨大hd| 亚洲国产精品久久男人天堂| 免费在线观看影片大全网站| 宅男免费午夜| 国产日本99.免费观看| 高潮久久久久久久久久久不卡| 999精品在线视频| 五月玫瑰六月丁香| 欧美乱妇无乱码| 免费在线观看日本一区| 午夜免费成人在线视频| 亚洲精品久久国产高清桃花| 一本精品99久久精品77| 每晚都被弄得嗷嗷叫到高潮| 麻豆国产97在线/欧美 | 日韩欧美 国产精品| 国产精品九九99| 狂野欧美白嫩少妇大欣赏| 波多野结衣巨乳人妻| 神马国产精品三级电影在线观看 | 一级黄色大片毛片| 亚洲熟妇熟女久久| 国内精品久久久久久久电影| 国内精品一区二区在线观看| 欧美黑人欧美精品刺激| 给我免费播放毛片高清在线观看| 啦啦啦韩国在线观看视频| 亚洲av五月六月丁香网| 老熟妇仑乱视频hdxx| 国产精品影院久久| 久久中文看片网| 中文字幕熟女人妻在线| 麻豆成人午夜福利视频| 久久久久精品国产欧美久久久| 免费一级毛片在线播放高清视频| 久久国产精品人妻蜜桃| 成人永久免费在线观看视频| 夜夜看夜夜爽夜夜摸| 五月伊人婷婷丁香| 叶爱在线成人免费视频播放| 成人手机av| 午夜久久久久精精品| 夜夜躁狠狠躁天天躁| 亚洲国产高清在线一区二区三| 国产在线观看jvid| 舔av片在线| 国产成人av教育| 一进一出抽搐动态| 欧洲精品卡2卡3卡4卡5卡区| 91老司机精品| 亚洲国产欧美网| 午夜成年电影在线免费观看| 免费观看精品视频网站| 麻豆成人av在线观看| 日韩欧美精品v在线| 国产探花在线观看一区二区| 国产精品免费视频内射| 国产一区二区在线观看日韩 | 免费一级毛片在线播放高清视频| 免费在线观看视频国产中文字幕亚洲| 国产亚洲欧美98| 免费无遮挡裸体视频| 欧美性长视频在线观看| 男女午夜视频在线观看| 波多野结衣高清作品| a级毛片在线看网站| 成熟少妇高潮喷水视频| 欧美精品啪啪一区二区三区| 久久久国产欧美日韩av| 国产精品日韩av在线免费观看| 国产成人影院久久av| 亚洲avbb在线观看| 亚洲全国av大片| 欧美成人性av电影在线观看| 国产激情欧美一区二区| 久久婷婷人人爽人人干人人爱| 亚洲一卡2卡3卡4卡5卡精品中文| 成年版毛片免费区| 天天添夜夜摸| 国产成人av激情在线播放| 亚洲专区中文字幕在线| 亚洲av熟女| 亚洲天堂国产精品一区在线| 老熟妇仑乱视频hdxx| 欧美zozozo另类| 欧美一级a爱片免费观看看 | 久久久久久亚洲精品国产蜜桃av| 麻豆久久精品国产亚洲av| 黄色视频,在线免费观看| 熟女少妇亚洲综合色aaa.| 在线观看www视频免费| 天天添夜夜摸| 久久久久久久久免费视频了| 免费看a级黄色片| 一级片免费观看大全| 精品国内亚洲2022精品成人| 村上凉子中文字幕在线| 久久久久九九精品影院| 亚洲av五月六月丁香网| АⅤ资源中文在线天堂| 国产久久久一区二区三区| 午夜视频精品福利| 麻豆av在线久日| 亚洲av片天天在线观看| 美女高潮喷水抽搐中文字幕| 人妻久久中文字幕网| 不卡一级毛片| 性色av乱码一区二区三区2| 国产又黄又爽又无遮挡在线| 国产亚洲精品第一综合不卡| 制服诱惑二区| 欧美性长视频在线观看| 日韩三级视频一区二区三区| 一本精品99久久精品77| 亚洲国产看品久久| 中文字幕精品亚洲无线码一区| 亚洲色图av天堂| 变态另类丝袜制服| 精品第一国产精品| 国产又色又爽无遮挡免费看| 久久香蕉激情| 在线观看免费视频日本深夜| 一个人免费在线观看电影 | 国内精品一区二区在线观看| 日本免费一区二区三区高清不卡| 一级毛片女人18水好多| 国产精品久久久久久久电影 | 久久精品成人免费网站| 国产亚洲欧美98| 99国产精品一区二区三区| 亚洲国产精品sss在线观看| 嫩草影院精品99| 免费在线观看成人毛片| 可以在线观看毛片的网站| 伊人久久大香线蕉亚洲五| 国产在线观看jvid| 亚洲av电影在线进入| a在线观看视频网站| 久久久久性生活片| 亚洲一区二区三区不卡视频| 久久人人精品亚洲av| 欧美日韩亚洲综合一区二区三区_| 正在播放国产对白刺激| 国产不卡一卡二| 欧美日韩一级在线毛片| 精品少妇一区二区三区视频日本电影| 亚洲一码二码三码区别大吗| 搡老熟女国产l中国老女人| 久久久久久久久久黄片| 欧美一级毛片孕妇| 欧美日本亚洲视频在线播放| 黄色片一级片一级黄色片| 国产精品98久久久久久宅男小说| 午夜亚洲福利在线播放| 成人一区二区视频在线观看| 日韩欧美三级三区| 亚洲第一电影网av| 亚洲国产精品sss在线观看| 一级作爱视频免费观看| 欧美精品亚洲一区二区| 丝袜美腿诱惑在线| 久久久久免费精品人妻一区二区| 亚洲精品美女久久久久99蜜臀| 婷婷丁香在线五月| 欧美日本亚洲视频在线播放| 非洲黑人性xxxx精品又粗又长| 变态另类成人亚洲欧美熟女| a级毛片在线看网站| 国产精品亚洲美女久久久| 五月玫瑰六月丁香| e午夜精品久久久久久久| 丝袜人妻中文字幕| 久热爱精品视频在线9| 丰满的人妻完整版| 日本一二三区视频观看| 午夜福利成人在线免费观看| 无限看片的www在线观看| 亚洲全国av大片| 少妇裸体淫交视频免费看高清 | 国产精品电影一区二区三区| 大型黄色视频在线免费观看| 最近最新中文字幕大全免费视频| 99热这里只有是精品50| 成人国产一区最新在线观看| 51午夜福利影视在线观看| www.熟女人妻精品国产| 亚洲 欧美 日韩 在线 免费| 变态另类丝袜制服| 最新在线观看一区二区三区| 老熟妇乱子伦视频在线观看| 国产精品久久久人人做人人爽| 精品欧美国产一区二区三| 窝窝影院91人妻| 最好的美女福利视频网| 久久精品国产亚洲av高清一级| 国产精品永久免费网站| 亚洲国产欧美网| 99热这里只有是精品50| 国产精品综合久久久久久久免费| 一a级毛片在线观看| 亚洲国产精品成人综合色| 亚洲男人的天堂狠狠| 久久香蕉精品热| 国产亚洲欧美在线一区二区| 97超级碰碰碰精品色视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | av欧美777| 真人一进一出gif抽搐免费| 午夜日韩欧美国产| 少妇人妻一区二区三区视频| 啪啪无遮挡十八禁网站| 黑人欧美特级aaaaaa片| 亚洲av成人不卡在线观看播放网| cao死你这个sao货| 亚洲一区二区三区不卡视频| 久久中文字幕人妻熟女| netflix在线观看网站| 亚洲熟妇熟女久久| 国产99久久九九免费精品| 久久精品aⅴ一区二区三区四区| 狂野欧美激情性xxxx| 亚洲精品在线观看二区| 亚洲人成电影免费在线| 老司机靠b影院| 这个男人来自地球电影免费观看| 在线a可以看的网站| 色播亚洲综合网| 国产精品乱码一区二三区的特点| 欧美日韩亚洲综合一区二区三区_| 日日夜夜操网爽| 久久99热这里只有精品18| 欧美黄色淫秽网站| 亚洲国产精品999在线| 成人手机av| av福利片在线观看| 日韩av在线大香蕉| 国产黄片美女视频| 亚洲aⅴ乱码一区二区在线播放 | 丝袜人妻中文字幕| 欧美性长视频在线观看| 极品教师在线免费播放| 久久久久久久久中文| 老汉色av国产亚洲站长工具| 99精品在免费线老司机午夜| 国产av不卡久久| 50天的宝宝边吃奶边哭怎么回事| 男女床上黄色一级片免费看| 最新美女视频免费是黄的| 99久久精品热视频| www.熟女人妻精品国产| 少妇被粗大的猛进出69影院| 蜜桃久久精品国产亚洲av| 亚洲国产看品久久| 精品欧美一区二区三区在线| 日韩欧美精品v在线| 老司机午夜十八禁免费视频| 婷婷精品国产亚洲av在线| 一进一出抽搐动态| www.自偷自拍.com| 无人区码免费观看不卡| 97碰自拍视频| 很黄的视频免费| 18禁观看日本| 一进一出抽搐动态| 亚洲全国av大片| 一进一出抽搐动态| 女人高潮潮喷娇喘18禁视频| 国产成人欧美在线观看| 国产一区二区在线观看日韩 | 亚洲一码二码三码区别大吗| 色老头精品视频在线观看| 少妇裸体淫交视频免费看高清 | 19禁男女啪啪无遮挡网站| 91老司机精品| 天堂√8在线中文| 免费在线观看日本一区| 亚洲精品粉嫩美女一区| 亚洲精品av麻豆狂野| 哪里可以看免费的av片| 国产精品野战在线观看| 久久久国产精品麻豆| 天堂av国产一区二区熟女人妻 | 久久伊人香网站| 久久99热这里只有精品18| 日韩欧美在线二视频| 在线观看一区二区三区| 亚洲国产精品久久男人天堂| 一进一出好大好爽视频| 欧美大码av| 美女高潮喷水抽搐中文字幕| 国产激情欧美一区二区| 久久久精品欧美日韩精品| 91九色精品人成在线观看| 91av网站免费观看| 国内精品一区二区在线观看| 午夜福利在线观看吧| 国产精品亚洲一级av第二区| 禁无遮挡网站| 两个人的视频大全免费| 免费在线观看黄色视频的| 嫩草影视91久久| 中文字幕精品亚洲无线码一区| 国产三级中文精品| 欧美3d第一页| 美女午夜性视频免费| 国产单亲对白刺激| 在线观看一区二区三区| 中文字幕最新亚洲高清| 欧美黑人欧美精品刺激| 波多野结衣巨乳人妻| 最新在线观看一区二区三区| 免费看a级黄色片| 亚洲免费av在线视频| 国产三级在线视频| 久久久久国产一级毛片高清牌| 五月伊人婷婷丁香| 国产av一区在线观看免费| 国产av一区二区精品久久| 色尼玛亚洲综合影院| 一进一出好大好爽视频| 50天的宝宝边吃奶边哭怎么回事| 国产三级黄色录像| 国产亚洲精品av在线| 亚洲精品美女久久av网站| 成熟少妇高潮喷水视频| 中文字幕久久专区| 嫩草影视91久久| 国产野战对白在线观看| 香蕉丝袜av| 午夜精品一区二区三区免费看| 亚洲国产精品sss在线观看| 亚洲欧美日韩高清专用| 欧美日韩精品网址| a级毛片a级免费在线| 中文在线观看免费www的网站 | 久久天躁狠狠躁夜夜2o2o| 91九色精品人成在线观看| 香蕉av资源在线| 少妇的丰满在线观看| 欧美一区二区精品小视频在线| a级毛片在线看网站| 国产免费男女视频| 长腿黑丝高跟| 日韩欧美国产在线观看| 97人妻精品一区二区三区麻豆| bbb黄色大片| 亚洲国产高清在线一区二区三| 亚洲欧美一区二区三区黑人| 可以在线观看毛片的网站| 久久精品aⅴ一区二区三区四区| 色综合亚洲欧美另类图片| 欧美色视频一区免费| 精品一区二区三区视频在线观看免费| 又大又爽又粗| 国内精品一区二区在线观看| 免费观看人在逋| 精品久久久久久久人妻蜜臀av| 久久精品国产清高在天天线| 国产单亲对白刺激| 久久人人精品亚洲av| 美女 人体艺术 gogo| 亚洲 欧美 日韩 在线 免费| 久久精品成人免费网站| 亚洲国产日韩欧美精品在线观看 | 亚洲成av人片在线播放无| 国产91精品成人一区二区三区| 国产久久久一区二区三区| 美女大奶头视频| 成人av在线播放网站| 麻豆成人av在线观看| 久久亚洲真实| 亚洲专区字幕在线| 九九热线精品视视频播放| 精品久久久久久久久久免费视频| 国产蜜桃级精品一区二区三区| 女人被狂操c到高潮| 国产一区二区三区视频了| 黄色视频,在线免费观看| 美女 人体艺术 gogo| 欧美成人午夜精品| 天天添夜夜摸| 丁香六月欧美| 日本精品一区二区三区蜜桃| √禁漫天堂资源中文www| 国产男靠女视频免费网站| 麻豆成人午夜福利视频| 亚洲成人免费电影在线观看| 久久亚洲真实| 少妇熟女aⅴ在线视频| 精品乱码久久久久久99久播| 在线播放国产精品三级| 天堂动漫精品| 精品国产乱子伦一区二区三区| 99热只有精品国产| 特大巨黑吊av在线直播| 国产蜜桃级精品一区二区三区| 淫妇啪啪啪对白视频| 国产精品爽爽va在线观看网站| 国内精品一区二区在线观看| 日本一二三区视频观看| 久久欧美精品欧美久久欧美| 精品欧美一区二区三区在线| 国产精品免费一区二区三区在线| 精品高清国产在线一区| 国产成人精品久久二区二区91| 精品一区二区三区视频在线观看免费| x7x7x7水蜜桃| 男女下面进入的视频免费午夜| 99热这里只有是精品50| 久久久国产精品麻豆| 757午夜福利合集在线观看| 俺也久久电影网| 久久久久亚洲av毛片大全| 欧洲精品卡2卡3卡4卡5卡区| 国产av不卡久久| 久久精品亚洲精品国产色婷小说| 老鸭窝网址在线观看| 亚洲精品在线美女| aaaaa片日本免费| 久久亚洲精品不卡| 淫妇啪啪啪对白视频| 亚洲精品中文字幕一二三四区| 亚洲国产高清在线一区二区三| 欧美色欧美亚洲另类二区| 男人舔女人下体高潮全视频| 欧美一级a爱片免费观看看 | 欧美日韩精品网址| a在线观看视频网站| 婷婷精品国产亚洲av在线| 欧美国产日韩亚洲一区| 亚洲熟妇熟女久久| 欧美精品啪啪一区二区三区| 久久午夜亚洲精品久久| 国产亚洲精品第一综合不卡| 级片在线观看| 亚洲中文av在线| 曰老女人黄片| bbb黄色大片| 亚洲精品久久成人aⅴ小说| 国产av麻豆久久久久久久| 人成视频在线观看免费观看| 97碰自拍视频| 亚洲欧洲精品一区二区精品久久久| 国产精品爽爽va在线观看网站| 变态另类成人亚洲欧美熟女| 老司机午夜十八禁免费视频| 亚洲成人国产一区在线观看| 免费在线观看视频国产中文字幕亚洲| 精品欧美国产一区二区三| 搡老岳熟女国产| 淫妇啪啪啪对白视频| 51午夜福利影视在线观看| 国产精品综合久久久久久久免费| 欧美国产日韩亚洲一区| 亚洲九九香蕉| 日本黄色视频三级网站网址| 又爽又黄无遮挡网站| 国产高清视频在线观看网站| 亚洲欧美日韩东京热| 国产av又大| 久久亚洲精品不卡| 欧美乱码精品一区二区三区| 欧美成狂野欧美在线观看| 巨乳人妻的诱惑在线观看| 成人欧美大片| 久久中文字幕人妻熟女| 黄色 视频免费看| 欧美另类亚洲清纯唯美| 亚洲狠狠婷婷综合久久图片| 国产午夜福利久久久久久| 狂野欧美激情性xxxx| av有码第一页| 欧美成人午夜精品| 99久久综合精品五月天人人| 亚洲av日韩精品久久久久久密| 精品少妇一区二区三区视频日本电影| 国产成人欧美在线观看| 女生性感内裤真人,穿戴方法视频| 黄色成人免费大全| 亚洲人成电影免费在线| 亚洲激情在线av| 免费在线观看视频国产中文字幕亚洲| 久99久视频精品免费| 国产免费av片在线观看野外av| 亚洲成av人片在线播放无| 制服诱惑二区| 国产蜜桃级精品一区二区三区| 无限看片的www在线观看| 可以在线观看的亚洲视频| 国产精品,欧美在线| 在线十欧美十亚洲十日本专区| 美女 人体艺术 gogo| 五月玫瑰六月丁香| 欧美日韩精品网址| 很黄的视频免费| 一二三四社区在线视频社区8| 一二三四在线观看免费中文在|