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

    Effect of Methane Emission Increases in East Asia on Atmospheric Circulation and Ozone

    2015-06-09 21:24:03SHANGLinLIUYiTIANWenshouandZHANGYuli
    Advances in Atmospheric Sciences 2015年12期

    SHANG Lin,LIU Yi,TIAN Wenshou,and ZHANG Yuli

    1Shandong Province Climate Centre,Jinan 250031

    2Key Laboratory of the Middle Atmosphere and Global Environmental Observation,Institute of Atmospheric Physics, Chinese Academy of Sciences,Beijing 100029

    3Key Laboratory for Semi-Arid Climate Change of the Ministry of Education,College of Atmospheric Sciences, Lanzhou University,Lanzhou 730000

    Effect of Methane Emission Increases in East Asia on Atmospheric Circulation and Ozone

    SHANG Lin1,LIU Yi?2,TIAN Wenshou3,and ZHANG Yuli2

    1Shandong Province Climate Centre,Jinan 250031

    2Key Laboratory of the Middle Atmosphere and Global Environmental Observation,Institute of Atmospheric Physics, Chinese Academy of Sciences,Beijing 100029

    3Key Laboratory for Semi-Arid Climate Change of the Ministry of Education,College of Atmospheric Sciences, Lanzhou University,Lanzhou 730000

    We used a fully coupled chemistry–climate model(version 3 of the Whole Atmosphere Community Climate Model, WACCM3)to investigate the effect of methane(CH4)emission increases,especially in East Asia and North America,on atmospheric temperature,circulation and ozone(O3).We show that CH4emission increases strengthen westerly winds in the Northern Hemisphere midlatitudes,accelerate the Brewer–Dobson(BD)circulation,and cause an increase in the mass flux across the tropopause.However,the BD circulation in the tropics between 10?S and 10?N at 100 hPa weakens as CH4emissions increase in East Asia and strengthens when CH4emissions increase in North America.When CH4emissions are increased by 50%in East Asia and 15%globally,the stratospheric temperature cools by up to 0.15 K,and the stratospheric O3increases by 45 ppbv and 60 ppbv,respectively.A 50%increase of CH4emissions in North America(with an amplitude of stratospheric O3increases by 60 ppbv)has a greater in fl uence on the stratospheric O3than the same CH4emissions increase in East Asia.CH4emission increases in East Asia and North America reduce the concentration of tropospheric hydroxyl radicals(4%and 2%,respectively)and increase the concentration of mid-tropospheric O3(5%and 4%,respectively)in the Northern Hemisphere midlatitudes.When CH4emissions increase in East Asia,the increase in the tropospheric O3concentration is largest in August.When CH4emissions increase in North America,the increase in the O3concentration is largest in July in the mid-troposphere,and in April in the upper troposphere.

    CH4,temperature,Brewer–Dobson circulation,O3

    1.Introduction

    Methane(CH4)is both an important greenhouse and chemically active gas.In the last glacial period(about 18 000 years ago),the concentration of CH4was 0.35 ppmv,but it reached 0.75 ppmv in the mid-19th century(e.g.,Rasmussen and Khalil,1984;Chappellaz et al.,1990).With the growth of human activities,the tropospheric CH4concentration reached 1.72 ppmv in 1995(e.g.,Lelieveld et al., 1998).Zhang et al.(2001)showed that from the end of the 20th century to the beginning of the 21st century,CH4emissions increased by 45%during three decades.Ramaswamy et al.(2001)pointed out that when the 1750 CH4mixing ratio of 0.72 ppmv increased to 1.77 ppmv in 2005,a radiative forcing of 0.48±0.05 W m?2was the result.CH4is the second most inf l uential radiative-forcing greenhouse gas after CO2, excluding water vapor.

    After the signing of the Montreal Protocol,atmospheric ozone(O3)depleting substances(ODSs)such as chlorides have been slowly reduced,while CH4may continue to increase significantly during the next 50 years.Satellite data show that in East Asia,with China’s rapid economic growth, the concentration of pollutants has continued to grow in recent years.Bergamaschi et al.(2013)found that from 2003 to 2008 CH4emissions caused by human activities grew fastest in China(18.6 Tg yr?1)—far faster than Brazil with the second fastest growth rate of 2.5 Tg yr?1.Zhang et al.(2011)analyzed Atmospheric Infrared Sounder satellite data and found that,compared with several major countries worldwide,the growth rate of CH4emissions has been significantly higher in China than other countries since 2007.

    In addition to changing the energy balance of the Earth’s atmosphere through the greenhouse effect,CH4emission increases also affect water vapor and the O3concentration in the stratosphere.Most CH4is oxidized in the troposphereby hydroxyl radicals(OH)(e.g.,Levy,1971),but about 7%–11%of CH4is transported into the stratosphere and oxidized (e.g.,Born etal.,1990).CH4oxidation is an importantsource of stratospheric water vapor(e.g.,Lelieveld and Crutzen, 1992;Fuglestvedt et al.,1996;Shindell et al.,2005).CH4emissions increase the O3concentration in the troposphere and stratosphere(e.g.,Owens et al.,1982,1985;Isaksen and Stordal,1986;Wuebbles and Hayhoe,2002).During the 20th century,the global background O3concentration rose by at least a factor of two,due mainly to increases in CH4and nitrogen oxides(NOx)emissions(e.g.,Marenco et al.,1994; Wang and Jacob,1998).

    Based on a series of model simulations,Butchart and Scaife(2001)found that mass flux transported from the tropical troposphere to stratosphere increases at a rate of 3%(10 yr)?1under conditions of increasing greenhouse gases,as supported by the IntergovernmentalPanelon Climate Change (IPCC)emissions scenarios.Increasing greenhouse gases enhances wave activities in the lower stratosphere and accelerates the Brewer–Dobson(BD)circulation(e.g.,Rind et al., 2001;Sigmond et al.,2004;Butchart et al.,2006).O3,water vapor and other trace gases in the troposphere are transported into the stratosphere faster as the BD circulation strengthens(e.g.,Austin et al.,2007).Austin and Li(2006)found that greenhouse gas and BD circulation increases jointly lead to increased CH4transport into the stratosphere,which is then oxidized into water vapor.The increase of stratospheric water vapor reduces the average age of the atmosphere in the stratosphere.The average age of the atmosphere indicates the transport time scales of substances leaving the troposphere and entering into the stratosphere(e.g.,Hall and Plumb,1994;Waugh and Hall,2002).

    The effects of a CH4concentration increase on atmospheric radiation and chemistry has been studied extensively (e.g.,Shindell et al.,2005;Chen et al.,2006;Shi,2006; Bi et al.,2007,2008;Guo et al.,2008;Bi,2009;Shi et al.,2009;Xie et al.,2013).Gruzdev and Brasseur(2005) used a two-dimensional model SOCRATES(Simulation of Chemistry,Radiation,and Transport of Environmentally Important Species,SOCRATES)and found that the stratosphere and middle atmosphere became cooler with an amplitude of less than 1 K when CH4was increased by 26%.Dyominov and Zadorozhny(2005)found CH4increases can cause an increase in stratospheric O3and lead to a warming in the stratosphere.However,the chemistry–radiation–dynamic interactions associated with CH4increases have not been well resolved in previous modeling studies,due to the limitation of model performance in the stratosphere.On the other hand, most previous studies have investigated the effects of globalscale CH4emission increases on climate and stratospheric O3,but the effect of a CH4emissions increase in East Asia on global temperature and O3is also still a popular and controversial issue.Wild and Palmer(2008)used a global chemistry transport model to demonstrate that the spatial extent of O3production and loss in the troposphere changes very little despite large projected increases in precursor emissions. Stevenson et al.(2006)used 26 atmospheric chemistry models to reveal that ensemble-mean changes in the tropospheric O3burden between 2000 and 2030 range significantly,but the effect of CH4increases on the tropospheric O3was not addressed.In this paper,we use a fully coupled chemistry–climate model to investigate the impact of a CH4emissions increase in East Asia on the stratospheric temperature,circulation and O3,and tropospheric O3.The results are compared with the impact of global CH4increases and a CH4emissions increase in North America,to assess the contribution of anthropogenic emissions in East Asia to climate change.

    2.Model description and numerical experiments

    We used version 3 of the Whole Atmosphere Community Climate Model(WACCM3),which is a global chemistry–climate model with 66 vertical levels extending from the surface to 4.5×10?6hPa(~160 km).WACCM3 is based on the software framework of the National Center for Atmospheric Research’s Community Atmosphere Model,version 3(CAM3),and includes all of the physical parameterizations of that model.The governing equations,physical parameterizations and numerical algorithms used in CAM3 are documented by Collins et al.(2004),and only the gravity wave drag and vertical diffusion parameterizations are modif i ed in WACCM3.The WACCM3 chemistry module is derived from the three-dimensional chemical transport model MOZART (Model for Ozone and Related Chemical Tracers)(e.g., Brasseur et al.,1998;Hauglustaine et al.,1998;Horowitz et al.,2003).This model resolves 51 neutral species,including all members of the Ox(oxides),NOx,HOx(hydrogen oxides),ClOx(chlorine oxides)and BrOx(bromine oxides) chemical families,along with tropospheric“source species”such as N2O(nitrous oxide),H2O,CH4,chlorof l uorocarbons (CFCs)and other halogenated compounds.WACCM3 performs well in simulating chemical and other atmospheric processes(e.g.,Garcia et al.,2007;Eyring et al.,2010;Liu and Liu,2009).Compared with observation,its simulations tend to be too warm in the Arctic winter,and produce cold temperatures and westerly winds that persist for too long in the Antarctic spring(e.g.,Garcia et al.,2007).These biases infl uence the temperature and O3concentration in the polar regions,so the simulation results in the polar regions are not discussed in detail in this study.

    We conducted six simulations at a horizontal resolution of 4?×5?with interactive chemistry,including one control experiment and fi ve sensitivity experiments,listed in Table 1.The surface emissions of CO2,CH4,Cl(chlorine),Br (bromine)and 14 other ODSs used in the control run(E0)are 12-month climatologies derived from IPCC’s A1B scenario (IPCC,2007)averaged over the period from 1979 to 2006. The A1 scenario meansa future world ofvery rapid economic growth,a global population that peaks mid-century and declines thereafter,and the rapid introduction of new and more ef fi cient technologies.Under the A1B scenario,technological changes occur in the energy system in a balanced manner across all sources.It has also been widely used in previouschemistry–climate model simulations in the literature to represent past and future ODS variations.The SST and sea ice fi elds used in the model are 12-month climatologies derived from Rayner et al.(2003).According to the IPCC 2007 A1B scenario,globally averaged CH4emissions increased by 15% in 2006,compared to 1979,and are expected to increase by 50%in 2050,meaning CH4emissions will increase observably in the future.In order to study the effect of different CH4emission increases in East Asia,in experiments E15, E30 and E50 we forced the surface CH4emissions in East Asia to increase by 15%,30%and 50%,respectively,compared to experiment E0(the increments were 7.5,15 and 25 Tg yr?1,approximately).East Asia was taken to range over the area(0?–90?N,60?–120?E).In order to study the effect of CH4emission increases in different areas,in experiment A50 the surface CH4emissions in North America were increased by 50%,the increment(25 Tg yr?1,approximately) and area of which were equal to E50 and the range of North America was taken as(0?–90?N,60?–120?W).In experiment G15 we forced global surface CH4emissions to increase by 15%(increment:90 Tg yr?1,approximately).The simulations were run for 50 years(from 2000 to 2050),with the first 10 years used as model spin-up,and the remaining 40 years (from 2010 to 2050)of model outputs used for analysis.

    3.Effect of CH4emission increases on stratospheric temperature and circulation

    Figure 1 shows the stratospheric temperature differences of E15,E30,E50,A50 and G15 relative to that in E0.When CH4emissions increase by 15%in East Asia,the temperature increases in the lower stratosphere in the Southern Hemisphere,but the temperature differences are not statistically signifi cant(Fig.1a).When CH4emissions increase by 30% in East Asia,there is a weaker temperature decrement in the lower stratosphere of the Southern Hemisphere compared with E15,but the temperature increases in the upper stratosphere(Fig.1b).The temperature difference is still not sta-tistically signifi cant in E30 either.When compared with the temperature differences induced by CH4emission increases of 15%and 30%,a signifi cant climate effect appears when CH4increases in East Asia by 50%,causing the middle stratosphere to cool significantly with a maximum cooling of 0.15 K at 30?N at 50 hPa(Fig.1c).Therefore,we base our conclusions on the results of E50 when studying the climate effect of CH4emission increases in East Asia in the following analysis.Figure 1 also indicates that a 50%CH4increase in North America causes a cooling in the middle stratosphere, with a maximum of 0.15 K at 30?N at 50 hPa.When global CH4emissions increase by 15%,the amplitude of the stratosphere maximum cooling is similar to that in E50 at 20?S at 60 hPa,but the range of cooling is larger than that in E50. That is to say,the greenhouse effect is much larger in experiment G15 than E50.

    Figure 2 shows the differences of the Eliassen–Palm (EP)flux and BD circulation of experiments E50,A50 and G15 compared with experiment E0.The BD circulation is a global-scale cell in the stratosphere in which air rises in the tropics and then moves polewards and downwards, mostly in the winter hemisphere.The BD circulation describes Lagrangian-mean transport and the transformed Eulerian mean residual velocities(e.g.,Andrews and McIntyre,1976,1978),which approximate the mean meridional mass transport for seasonally averaged conditions(e.g., Holton,1990).When CH4emissions increase,the stratospheric westerly weakens in the Southern Hemisphere midlatitudes and the mid-stratospheric westerly in the Northern Hemisphere strengthens.Compared with E0,the maximum westerly increment appears in the Northern Hemisphere in E50 and A50,and in the Southern Hemisphere in G15:the location is consistent with Fig.1,where the temperature gradient is evident.When the extratropical temperature gradient increases,the westerly and wave activity increase(e.g., Eichelberger and Hartmann,2005;Olsen et al.,2007;Garcia and Randel,2008).As shown in Fig.2,as the wave force becomes stronger in the Northern Hemisphere midlatitudes,the vertical velocity of BD circulation speeds up.So, the changes of the mean zonal wind,EP flux and BD circulation are consistent,as shown in Fig.2.Fomichev et al. (2007)found that the EP flux increases when the westerly is stronger at midlatitudes and the wave flux change in the range 20?–40?N is responsible for the tropical upwelling change. Garcia and Randel(2008)found that the maximum increase of EP flux appeared in the tropics during 1950–2003,and in the midlatitudes during 1980–2050,causing the BD circulation to change accordingly.It is worth noting that the vertical velocity of EP flux at 100 hPa in the Northern Hemisphere midlatitudes changes most significantly in all the sensitivity experiments,and the vertical velocity of BD circulation also increases significantly south of the equator in A50;but,in E50 and G15,the vertical velocity of BD circulation reduces south of the equator at 100 hPa.This is further analyzed laterin the paper.

    Using a series of models(including WACCM),Butchart et al.(2006)found that BD circulation upwelling regions extend roughly 30?either side of the equator,though there are variations in the width of about 10?between the models.In our simulations,the BD circulation rises between 22?S and 30?N under the SST and greenhouse gas concentrations used. So,the region between 22?S and 30?N is def i ned as the BD circulation rising zone.Lin and Fu(2013)found that the acceleration rate of the BD circulation is much larger in the lower stratosphere than upper stratosphere as greenhouse gas concentrations increase.As previous studies have used the upward mass flux across 100 hPa to infer the strength of the BD circulation(e.g.,Garcia and Randel,2008;Garny et al., 2011),Table 2 shows the changes of upward mass flux of E50,A50 and G15 at different altitudes.The increment of upward mass flux is larger at 100 hPa than at 70 hPa and 30 hPa,i.e.,when CH4emissions increase(in East Asia,North America and on the global scale)the increment of upward mass flux and the acceleration of the BD circulation is larger in the lower than in the upper stratosphere.It is worth noting that the change of upward mass flux is largest in experiment A50;so,compared with a 50%CH4emissions increase in East Asia,there is a larger effect of the CH4emissions increase in North America on the EP flux,BD circulation and mass flux.

    In order to study the BD circulation vertical velocity difference as CH4emissions increase by different concentrations and in differentareas,Fig.3 furthershows the difference of the BD circulation vertical velocity between sensitivityexperiments and the controlexperimentE0.When CH4emissions increase in East Asia and globally,the difference ofis rathersmallbetween 22?S and 10?S,while in the equatorial regionis stronger between 10?N and 30?N at 100 hPa.When CH4emissions increase in North America,the difference ofis similar to that in other experiments,but it is very different between 10?S–30?N,as in the equatorial region(10?S–10?N)becomes stronger and between 20?–30?Nis weaker. In the following analysis,the region between 22?S and 30?N is divided into three parts:10?–30?N,10?S–10?N and 10?–22?S.

    CH4emission increases accelerate the BD circulation, which can transport more air mass from the troposphere into the stratosphere.Table 3 shows the differences of upward mass flux at 100 hPa in experiments E15,E30,E50,A50 and G15,relative to that in experiment E0.Note that the upward mass flux between 10?S and 10?N decreases in all experiments except A50,i.e.,E15,E30,E50 and G15,compared to that in experiment E0.The changes of upward mass flux are consistent with BD circulation vertical velocity variations.

    Figure 4 shows the seasonal variations of the BD circulation vertical velocity.It is apparent that the vertical velocity variations caused by different levels of increase in CH4surface emissions are quite different,although the vertical velocity reduces in January and February in all the sensitivity experiments between 10?S and 10?N,and in this range the BD circulation ascending branch becomes stronger during March and June in E15.This change is not obvious in E30, and the vertical velocity reduces during May and August in E50.Table 3 shows that between 22?S and 10?S the upward mass flux decreases in E15 and E30,but increases in E50. When CH4emissions increase in East Asia the increment of upward mass flux between 10?N and 30?N is largest in E30. Above all,the change of upward mass flux at 100 hPa is nonlinear as the CH4emissions gradually increase in East Asia in E15,E30 and E50.In A50 the vertical velocity decreases slightly in January and February between 10?S and 10?N and increases obviously between 10?N and 20?N from January to May.The combined effect is that the upward mass flux increment at 100 hPa is largest in A50.There is no obvious change of mass flux between 10?S and 10?N in G15,but the mass flux increase is remarkable between 10?N and 30?N,asindicated by the noteworthy vertical velocity increase in January and February and the slight increase in other months. Therefore,there is also a large increase of the total mass flux in G15 between 22?S and 30?N.

    Table 3.The difference of mass flux between each sensitivity experiment and E0 at 100 hPa.Positive values indicate a mass flux increase and negative values indicate a mass flux decrease.

    4.Effect of CH4emission increases on stratospheric and tropospheric ozone

    Previous studies have shown that the tropospheric and stratospheric O3increases as CH4emissions increase(e.g., Owens et al.,1982,1985;Isaksen and Stordal,1986;Wuebbles and Hayhoe,2002).Figure 5 shows the difference of stratospheric O3in E50,A50 and G15.In the middle stratosphere,CH4reacts with active chlorine(CH4+Cl→CH3+ HCl)and decreases catalytic O3loss caused by Cl radicals; this is the main cause of O3increase in the middle stratosphere.

    When CH4emissions increase by 50%in East Asia(E50) and 15%on the global scale(G15),the maximum increase of stratospheric O3is 45 ppbv and 60 ppbv,respectively.Therefore,a CH4emissions increase is benef i cial for stratospheric O3recovery.It should be noted that although the increment of CH4emissions in East Asia and North America is the same,the stratospheric O3increase in A50 is larger,and the maximum increase reaches 60 ppbv.That is to say,when CH4emissions increase by the same amount in East Asia and North America,the effect on the stratospheric O3of the former is less than the latter.

    Around 90%of atmospheric CH4is removed through chemical reactions in the atmosphere,and 10%is removed by soil oxidation.Figure 6 shows the OH and O3difference in the troposphere between each sensitivity experiment and experiment E0.The reaction between CH4and OH radicals (CH4+OH→CH3+H2O)is the main removal mechanism of CH4in the troposphere.In the stratosphere,CH4reacts with the OH and O1(D)radical;however,the photolysis process of CH4can be ignored(e.g.,Qin and Zhao,2003).The photochemical reaction of the OH radical is important for oxidation and removal of numerous trace gases.In general,the concentration of the OH radical is highest in the tropics,and slightly higher in the Southern than the Northern Hemisphere [e.g.,Chameides and Davis(1982)].As shown in Figs.6a and b,OH in the Northern Hemisphere significantly reduces as CH4oxidation by OHmainly occurs in the Northern Hemisphere when CH4emissions increase in East Asia and North America.CH4is oxidized by OH on the global scale in G15,so the tropospheric OH reduces uniformly in both hemispheres.

    One quarterofthe hydroxylradicals reactwith CH4,starting several oxidation chains that strongly affect the atmospheric concentrations of OH and O3(e.g.,Crutzen,1973). The availability of NO plays an important role in determining the oxidation pathways:

    CH4oxidation by OH mainly occurs in the Northern Hemisphere troposphere when CH4emissions increase in East Asia and North America.This causes O3to increase significantly in the Northern Hemisphere troposphere.The concentration of CH4is higher at high latitudes than at low latitudes.With the same percentage increase of CH4in East Asia and North America,the increment of CH4at high latitudes is larger than at low latitudes.But,as shown in Figs.6dand e,the maximum percentage increase of the tropospheric O3appears in the low-and midlatitudes.This is because CH4oxidation by OH mostly occurs in the tropics and the concentration of NOxis largest in the midlatitudes(e.g.,Fiore et al.,2008).It should be noted that OH reduces by 4%,and O3in the Northern Hemisphere midlatitudes increases by 4% in E50.There is a small variation in A50,with OH reducing by 2%and O3in the Northern Hemisphere midlatitudes increasing by 2%.It is also interesting that the altitude of the maximum changes of OH and O3is higher in E50 than in A50.This is because the concentration of OH is higher in EastAsia than North America in the troposphere(notshown), and when CH4emissions increase in East Asia,CH4could be fully oxidized and generate more O3in the troposphere. When CH4emissions increase in North America,the oxidation reaction of CH4and OH is not suff i cient,so more CH4is transported into the stratosphere,where it is oxidized,and more O3is generated in the stratosphere.

    O3is generated through the CH4oxidation in the presence of NOxin the troposphere(e.g.,Crutzen,1973).During the 20th century,global background O3concentrations rose by at least a factor of two,due mainly to increases in CH4and NOxemissions(e.g.,Marenco et al.,1994;Wang and Jacob, 1998).Therefore,the distribution of CH4and NOxin the troposphere plays an important role in generating O3(e.g.,Fiore et al.,2008).

    Figure 7 shows O3differences at different heights.There is an obvious local impact of CH4emission increases,and the maximum increase of O3is 5%in E50 and 4%in A50 at 500 hPa.In the upper troposphere at 200 hPa,in E50 the O3concentration increases by 6%in East Asia,higher than surrounding regions.In A50 the O3concentration increment in North America is higher than surrounding regions,but the O3does not increase significantly in the other lower latitude regions,and the O3increment in North America is only 4%, much smaller than E50.The O3increase is symmetrically distributed along the equator in G15.

    Tropospheric O3is an important greenhouses gas,and as a chemically active gas it can react with numerous active substances.High concentrations of tropospheric O3are harmful to health,affecting the respiratory system.When CH4emissions increase by the same amount in East Asia and North America,the effect on tropospheric O3in the former is larger than the latter.

    Lang et al.(2012)found that the seasonal variation of O3change due to CH4emission increases in the troposphere is very weak,except near the tropical tropopause where O3has a maximum increase in July.As shown in Fig.7,the local change of tropospheric O3is signifi cant when CH4emissions increase in East Asia and North America;therefore,the O3changes in different months in East Asia and North America are interesting.The distribution of tropospheric O3concentrations in East Asia and North America are shown in Fig.8a, and are in good agreement with Wang et al.(2006).As shown in Fig.8b,there isan obviousseasonalvariation ofO3change in East Asia when CH4emissions increase in East Asia in the upper troposphere(200 hPa).The maximum increment appears in August.The seasonal variation of O3change inNorth America in A50 shows that the maximum increment appears in April.In G15,the variation is small in different monthsand the increment ofO3is large from May to September.

    O3increases are more consistent at 500 hPa than at 200 hPa,as the maximum increment of O3concentration appears in August in East Asia and in July in North America.However,the seasonal variation of O3change in G15 is not obvious(Fig.8c)in either East Asia or North America.This proves that seasonal variation of O3concentrations,shown in Fig.8,are affected by CH4emission increases locally in the lower troposphere.Comparing Figs.8a–c,the maximum O3increment appears when the original O3concentration is at its minimum in East Asia,so the percentage change of O3concentration is similar to Figs.8b and c.Although the O3concentration increase is smaller in the middle troposphere than in the upper troposphere,the seasonal variation is also signifi cant in the middle troposphere.The O3concentration is lower at 500 hPa than 200 hPa,so the percentage change of O3concentration is larger at 500 hPa than 200 hPa(not shown).

    Since CH4oxidation by OH depends on temperature, when CH4emissions increase the largest O3concentration increment should be in summer.But in A50,the maximum increment appears in April at 200 hPa.The seasonal variation of O3change in A50 is inconsistent with the seasonal variation of O3concentration simulated by the control experiment shown in Figs.8a and b.That is,the O3concentration change in A50 at 200 hPa is caused by transport processes. So,when CH4emissions increase in East Asia,its impact on the tropospheric O3can be noted throughout the troposphere. However,when CH4emissions increase in North America, the impact on the tropospheric O3can be noted only in the lower troposphere and the change of O3concentration in the upper troposphere is caused by transport processes.

    5.Summary and conclusions

    A fully coupled chemistry–climate model(WACCM3) was used to investigate the effect of CH4emission increases, especially in East Asia and North America,on atmospheric temperature,circulation and O3.

    The results show that CH4emission increases can strengthen the westerly winds in the Northern Hemisphere midlatitudes,and accelerate the BD circulation,increasing the mass flux across the tropopause.However,the tropical BD circulation between 10?S and 10?N at 100 hPa becomes weaker when CH4emissions increase in East Asia, and stronger when CH4emissions increase in North America. When CH4emissions increase by 50%in East Asia and 15% globally,the stratospheric temperature cools by up to 0.15 K, and the stratospheric O3increases by 45 ppbv and 60 ppbv, respectively.When CH4emissions increase by 50%in North America the stratospheric O3increases by 60 ppbv,so it has a greater inf l uence on stratospheric O3than the effect of a CH4emissions increase of the same amount in East Asia.CH4increases in East Asia and North America can reduce the tropospheric OH concentration(4%and 2%,respectively)and increase the tropospheric O3concentration(5%and 4%,respectively)in the Northern Hemisphere midlatitudes.When CH4increases in East Asia,the increase of the tropospheric O3concentration islargestin August.When CH4increasesin North America,the increase of the O3concentration reaches a maximum in July in the middle troposphere,and in April in the upper troposphere.When CH4emissions increase in East Asia,the impact on tropospheric O3can be noted throughout the troposphere.However,when CH4emissions increase in North America,the impact on tropospheric O3is notable only in the lower troposphere and the change of O3concentration in the upper troposphere is caused by transport processes.

    The emissions scenarios used in this paper are provided by the IPCC at a resolution of 4?×5?,so the CH4emissions source is coarse.CH4reacts with CO,NOxand other chemical substances,affecting tropospheric O3and stratospheric climate.Against the background of increasing pollutant emissions,the effect of a CH4emissions increase must be associated with CO,NOxand other substances.There-fore,much more work is needed to fully understand the effect of CH4emission increases on both tropospheric and stratospheric climate in East Asia.

    Acknowledgements.This work was supported by the China High Resolution Earth Observation Project(Grant No.1-113406) and the National Science Foundation of China(Grant Nos. 41225018 and 41175042).

    REFERENCES

    Andrews,D.G.,and M.E.McIntyre,1976:Planetary waves in horizontal vertical shear:The generalized Eliassen–Palm relation and the mean zonal acceleration.J.Atmos.Sci.,33, 2031–2048.

    Andrews,D.G.,and M.E.McIntyre,1978:Generalized Eliassen–Palm and Charney–Drazin theorems for waves on axismmetric mean fl ows in compressible atmospheres.J.Atmos.Sci., 35,175–185.

    Austin,J.,and F.Li,2006:On the relationship between the strength of the Brewer-Dobson circulation and the age of stratospheric air.Geophys.Res.Lett.,33,L17807,doi:10.1029/ 2006GL026867.

    Austin,J.,J.Wilson,F.Li,and H.Vo¨mel,2007:Evolution ofwater vapor concentrations and stratospheric age of air in coupled chemistry-climate model simulations.J.Atmos.Sci.,64(3), 905–921.

    Bergamaschi,P.,and Coauthors,2013:Atmospheric CH4in the first decade of the 21st century:Inverse modeling analysis using SCIAMACHY satellite retrievals and NOAA surface measurements.J.Geophys.Res.,118(13),7350–7369.

    Bi,Y.,2009:Study on the distributions,variations and climate impacts of the water vapor and methane in stratosphere.Ph.D. dissertation,Space Physics,University of Science and Technology of China,149 pp.(in Chinese)

    Bi,Y.,Y.J.Chen,L.Xu,S.M.Deng,and R.J.Zhou,2007:Analysis of H2O and CH4distribution characteristics in the middle atmosphere using HALOE data.Chinese Journal of Atmospheric Sciences,31(3),440–448.(in Chinese with English abstract)

    Bi,Y.,Y.J.Chen,R.J.Zhou,M.H.Fang,and L.Xu,2008:Study on H2O and CH4distributions and variations over Qinghai-Xizang Plateau using HALOE data.Plateau Meteorology, 27(2),249–258.(in Chinese with English abstract)

    Born,M.,H.Do¨rr,and I.Levin,1990:Methane consumption in aerated soils of the temperate zone.Tellus B,42(1),2–8.

    Brasseur,G.P.,D.A.Hauglustaine,S.Walters,P.J.Rasch,J.-F. Mu¨ller,C.Granier,and X.X.Tie,1998:MOZART:A global chemical transport model for ozone and related chemical tracers:1.Model description.J.Geophys.Res.,103,28 265–28 289.

    Butchart,N.,and A.A.Scaife,2001:Removal of chloro fl uorocarbons by increased mass exchange between the stratosphere and troposphere in a changing climate.Nature,410(6830), 799–802.

    Butchart,N.,and Coauthors,2006:Simulations of anthropogenic change in the strength of the Brewer-Dobson circulation.Climate Dyn.,27(7–8),727–741.

    Chameides,W.L.,and D.D.Davis,1982:The free radical chemistry of cloud droplets and its impact upon the composition of rain.J.Geophys.Res.,87(C7),4863–4877.

    Chappellaz,J.,J.M.Barnola,D.Raynaud,Y.S.Korotkevich,and C.Lorius,1990:Ice-core record of atmospheric methane over the past 160,000 years.Nature,345(6271),127–131.

    Chen,Y.J.,R.J.Zhou,C.H.Shi,and Y.Bi,2006:Study on the trace species in the stratosphere and their impact on climate.Adv.Atmos.Sci.,23(6),1020–1039,doi:10.1007/s00376-006-1020-3.

    Collins,W.D.,and Coauthors,2004:Description of the NCAR Community Atmosphere Model(CAM3).Tech.Note NCAR/ TN-464+STR,Natl.Center for Atmos.Res.,226 pp.

    Crutzen,P.,1973:A discussion of the chemistry of some minor constituents in the stratosphere and troposphere.Pure Appl. Geophys.,106(1),1385–1399.

    Dyominov,I.G.,and A.M.Zadorozhny,2005:Greenhouse gases and recovery of the Earth’s ozone layer.Advances in Space Research,35,1369–1374.

    Eichelberger,S.J.,and D.L.Hartmann,2005:Changes in the strength of the Brewer-Dobson circulation in a simple AGCM.Geophys.Res.Lett.,32,L15807,doi:10.1029/2005 GL022924.

    Eyring,V.,and Coauthors,2010:Multi-model assessment of stratospheric ozone return dates and ozone recovery in CCMVal-2 models.Atmospheric Chemistry and Physics, 10(19),9451–9472.

    Fiore,A.M.,J.J.West,L.W.Horowitz,V.Nail,and M.D. Schwarzkopf,2008:Characterizing the tropospheric ozone response to methane emission controls and the benefits to climate and air quality.J.Geophys.Res.,113,D08307,doi: 10.1029/2007JD009162.

    Fomichev,V.I.,A.I.Jonsson,J.de Grandpre,S.R.Beagley,C. McLandress,K.Semeniuk,and T.G.Shepherd,2007:Response of the middle atmosphere to CO2doubling:Results from the Canadian Middle Atmosphere Model.J.Climate, 20(7),1121–1144.

    Fuglestvedt,J.S.,I.S.A.Isaksen,and W.C.Wang,1996:Estimates of indirect global warming potentials for CH4,CO and NOx.Climatic Change,34(3–4),405–437.

    Garcia,R.R.,and W.J.Randel,2008:Acceleration of the Brewer-Dobson circulation due to increases in greenhouse gases.J. Atmos.Sci.,65,2731–2739.

    Garcia,R.R.,D.R.Marsh,D.E.Kinnison,B.A.Boville,and F.Sassi,2007:Simulation of secular trends in the middle atmosphere,1950–2003.J.Geophys.Res.,112,D09301,doi: 10.1029/2006JD007485.

    Garny,H.,M.Dameris,W.Randel,G.E.Bodeker,and R.Deckert,2011:Dynamically forced increase of tropical upwelling in the lower stratosphere.J.Atmos.Sci.,68,1214–1233.

    Gruzdev,A.N.,and G.P.Brasseur,2005:Long-term changes in the mesosphere calculated by a two-dimensional model.J. Geophys.Res.,110,304–321.

    Guo,S.C.,H.Zhou,D.R.Lv,Y.Q.Li,M.Dai,and Q.Li,2008: Temporal and spatial features of atmospheric methane and its relation to ozone variation in the stratosphere.Journal of Yunnan University(Natural Sciences Edition),30(4),381–387. (in Chinese with English abstract)

    Hall,T.M.,and R.A.Plumb,1994:Age as a diagnostic of stratospheric transport.J.Geophys.Res.,99(D1),1059–1070.

    Hauglustaine,D.A.,G.P.Brasseur,S.Walters,P.J.Rasch,J.-F. M¨uller,L.K.Emmons,and M.A.Carroll,1998:MOZART: A global chemical transport model for ozone and related chemical tracers:2.Model results and evaluation.J.Geophys. Res.,103,28 291–28 335.

    Holton,J.R.,1990:On the global exchange of mass between the stratosphere and troposphere.J.Atmos.Sci.,48,392–395.

    Horowitz,L.W.,and Coauthors,2003:A global simulation of tropospheric ozone and related tracers:Description and evaluation of MOZART,version 2.J.Geophys.Res.,108(D24), doi:10.1029/2002JD002853.

    IPCC,2007:Climate Change 2007:The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press,New York,996 pp.

    Isaksen,I.S.A.,and F.Stordal,1986:Antarctic ozone depletion: 2-D model studies.Geophys.Res.Lett.,13(12),1327–1330.

    Lang,C.,and Coauthors,2012:The impact of greenhouse gases on past changes in tropospheric ozone.J.Geophys.Res.,117, D23304,doi:10.1029/2012JD018293.

    Lelieveld,J.,and P.J.Crutzen,1992:Indirect chemical effects of methane on climate warming.Nature,355,339–342.

    Lelieveld,J.,P.J.Crutzen,and F.J.Dentener,1998:Changing concentration,lifetime and climate forcing of atmospheric methane.Tellus B,50(2),128–150.

    Levy,H.II.,1971:Normal atmosphere:Large radical and formaldehyde concentrations predicted.Science,173,141–143.

    Lin,P.,and Q.Fu,2013:Changes in various branches of the Brewer-Dobson circulation from an ensemble of chemistry climate models.J.Geophys.Res.,118(1),73–84.

    Liu,Y.,and C.X.Liu,2009:Simulation studies on seasonal variations of the stratospheric dynamics and trace gases using coupled chemistry-climate model WACCM-3.Chinese Journal of Space Science,29(6),580–590.(in Chinese with English abstract)

    Marenco,A.,H.Gouget,P.N′ed′elec,J.-P.Pag′es,and F.Karcher, 1994:Evidence of a long-term increase in tropospheric ozone from Pic du Midi data series:Consequences:Positive radiative forcing.J.Geophys.Res.,99(D8),16 617–16 632.

    Olsen,M.A.,M.R.Schoeberl,and J.E.Nielsen,2007:Response of stratospheric circulation and stratosphere-troposphere exchange to changing sea surface temperatures.J.Geophys. Res.,112,D16104,doi:10.1029/2006JD008012.

    Owens,A.J.,J.M.Steed,D.L.Filkin,C.Miller,and J.P.Jesson, 1982:The potential effects of increased methane on atmospheric ozone.Geophys.Res.Lett.,9(9),1105–1108.

    Owens,A.J.,C.H.Hales,D.L.Filkin,C.Miller,J.M.Steed, and J.P.Jesson,1985:A coupled one-dimensional radiativeconvective,chemistry-transport model of the atmosphere:1. Model structure and steady state perturbation calculations.J. Geophys.Res.,90,2283–2311.

    Qin,Y.,and C.S.Zhao,2003:Basic Atmospheric Chemistry.Meteorology Press,202 pp.(in Chinese)

    Ramaswamy,V.,and Coauthors,2001:Radiative forcing of climate change.Climate Change 2001:The Scientif i c Basis,J. T.Houghton et al.,Eds.,Cambridge University Press,New York,349–416.

    Rasmussen,R.A.,and M.A.K.Khalil,1984:Atmospheric methane in the recent and ancient atmospheres:Concentrations,trends,and interhemispheric gradient.J.Geophys.Res., 89(D7),11 599–11 605.

    Rayner,N.A.,and Coauthors,2003:Global analyses of sea surface temperature,sea ice,and night marine air temperature since the late nineteenth century.J.Geophys.Res.,108(D14), doi:10.1029/2002JD002670.

    Rind,D.,J.Lerner,and C.McLinden,2001:Changes of tracer distributions in the doubled CO2 climate.J.Geophys.Res., 106(D22),28 061–28 079.

    Shi,C.H.,2006:Study on the trends and chemical process of trace gases in stratosphere.PhD.dissertation,Space Physics, University of Science and Technology of China,191 pp.(in Chinese)

    Shi,C.H.,B.Zheng,Y.J.Chen,and Y.Bi,2009:The quasibiennial oscillation of water vapor in tropical stratosphere.Chinese Journal of Geophysics,52(10),2428–2435.(in Chinese with English abstract)

    Shindell,D.T.,G.Faluvegi,N.Bell,and G.A.Schmidt,2005: An emissions-based view of climate forcing by methane and tropospheric ozone.Geophys.Res.Lett.,32(4),L04803,doi: 10.1029/2004GL021900.

    Sigmond,M.,P.C.Siegmund,E.Manzini,and H.Kelder, 2004:A simulation of the separate climate effects of middleatmospheric and tropospheric CO2doubling.J.Climate, 17(12),2352–2367.

    Stevenson,D.,and Coauthors,2006:Multimodel ensemble simulations of present-day and near-future tropospheric ozone.J. Geophys.Res.,111(D8),doi:10.1029/2005JD006338.

    Wang,Y.H.,and D.J.Jacob,1998:Anthropogenic forcing on tropospheric ozone and OH since preindustrial times.J.Geophys.Res.,103(D23),31 123–31 135.

    Wang,W.G.,M.Yuan,J.Wu,W.X.Fan,H.Y.Wang,and X. L.Liu,2006:The variation of spatial temporal distribution of the global tropopause ozone.Journal of Yunnan University (Natural Sciences Edition),28(6),509–517.(in Chinese with English abstract)

    Waugh,D.,and T.Hall,2002:Age of stratospheric air:Theory, observations,and models.Rev.Geophys.,40,1-1–1-26.

    Wild,O.,and P.I.Palmer,2008:How sensitive is tropospheric oxidation to anthropogenic emissions?Geophys.Res.Lett., 35(22),doi:10.1029/2008GL035718.

    Wuebbles,D.J.,and K.Hayhoe,2002:Atmospheric methane and global change.Earth-Science Reviews,57(3–4),177–210.

    Xie,F.,W.S.Tian,J.P.Li,J.K.Zhang,and L.Shang,2013:The possible effects of future increase in methane emission on the stratospheric water vapor and global ozone.Acta Meteorologica Sinica,71(3),555–567.(in Chinese with English abstract)

    Zhang,R.J.,M.X.Wang,and Y.S.Wang,2001:Long-term trends of atmospheric methane and its future change.Climate and Environmental Research,6(1),53–57.(in Chinese with English abstract)

    Zhang,X.Y.,W.G.Bai,P.Zhang,and W.H.Wang,2011: Spatiotemporal variations in mid-upper tropospheric methane over China from satellite observations.Chinese Science Bulletin,56,3321–3327.

    :Shang,L.,Y.Liu,W.S.Tian,and Y.L.Zhang,2015:Effect of methane emission increases in East Asia on atmospheric circulation and ozone.Adv.Atmos.Sci.,32(12),1617–1627,

    10.1007/s00376-015-5028-4.

    22 January 2015;revised 10 May 2015;accepted 4 June 2015)?

    LIU Yi Email:liuyi@mail.iap.ac.cn

    男女视频在线观看网站免费| av欧美777| 啦啦啦韩国在线观看视频| 精品久久久久久久久久久久久| 久久精品国产99精品国产亚洲性色| 老司机午夜福利在线观看视频| 国产v大片淫在线免费观看| 国产91精品成人一区二区三区| 国产精品亚洲美女久久久| 国产精品日韩av在线免费观看| 国内久久婷婷六月综合欲色啪| 亚洲国产欧美人成| 欧美日韩精品网址| 好男人在线观看高清免费视频| 999久久久精品免费观看国产| 99国产精品一区二区蜜桃av| 男人舔女人的私密视频| 欧美成人免费av一区二区三区| 五月伊人婷婷丁香| 99久国产av精品| 宅男免费午夜| 美女被艹到高潮喷水动态| 国产精品久久久久久人妻精品电影| 久久久精品大字幕| 久久精品国产清高在天天线| 网址你懂的国产日韩在线| 国产亚洲精品久久久com| 国产男靠女视频免费网站| 欧美大码av| 日本a在线网址| 国产三级中文精品| 免费在线观看亚洲国产| 久久精品影院6| 国语自产精品视频在线第100页| 欧美中文日本在线观看视频| 麻豆国产97在线/欧美| 亚洲av五月六月丁香网| 嫁个100分男人电影在线观看| 两个人的视频大全免费| 在线观看午夜福利视频| 男人舔女人的私密视频| 欧美日本视频| 欧美中文日本在线观看视频| 午夜亚洲福利在线播放| 精华霜和精华液先用哪个| 12—13女人毛片做爰片一| 亚洲午夜精品一区,二区,三区| 国产av麻豆久久久久久久| 亚洲av成人一区二区三| 亚洲专区字幕在线| 国产在线精品亚洲第一网站| 日本黄色片子视频| 嫁个100分男人电影在线观看| 18禁美女被吸乳视频| 亚洲欧美日韩卡通动漫| 精品免费久久久久久久清纯| 国产精品九九99| 在线免费观看的www视频| 免费电影在线观看免费观看| 久久天堂一区二区三区四区| 人人妻,人人澡人人爽秒播| 国产精品久久久久久亚洲av鲁大| 国产伦在线观看视频一区| 日本在线视频免费播放| 国产高清三级在线| 国内精品久久久久精免费| 看黄色毛片网站| 91九色精品人成在线观看| 国产精品99久久久久久久久| 日韩国内少妇激情av| 特级一级黄色大片| 国产精品1区2区在线观看.| 国产成+人综合+亚洲专区| 国产精品99久久久久久久久| 亚洲最大成人中文| 国产乱人伦免费视频| 日本黄大片高清| 12—13女人毛片做爰片一| av天堂在线播放| 无遮挡黄片免费观看| 亚洲国产精品sss在线观看| 午夜福利高清视频| 亚洲精品在线观看二区| 国产精品综合久久久久久久免费| 久久精品国产综合久久久| 欧美日韩亚洲国产一区二区在线观看| 国产美女午夜福利| 日本a在线网址| 香蕉丝袜av| 久久久久久久午夜电影| 亚洲精品色激情综合| 麻豆国产97在线/欧美| 精品免费久久久久久久清纯| 两性夫妻黄色片| 欧美乱妇无乱码| 男女之事视频高清在线观看| 国产野战对白在线观看| 老汉色∧v一级毛片| 亚洲电影在线观看av| 国产精品1区2区在线观看.| 久久香蕉精品热| 国产一区二区三区视频了| 欧美在线黄色| 色吧在线观看| 色综合站精品国产| 90打野战视频偷拍视频| 国产精品久久久久久精品电影| 日本与韩国留学比较| 麻豆久久精品国产亚洲av| 午夜两性在线视频| 国产精品久久久久久久电影 | 久久国产乱子伦精品免费另类| 国内精品一区二区在线观看| 国产乱人伦免费视频| 国产黄片美女视频| 999久久久精品免费观看国产| www.自偷自拍.com| 久久国产乱子伦精品免费另类| 久久99热这里只有精品18| 香蕉av资源在线| 五月玫瑰六月丁香| 国产精品香港三级国产av潘金莲| 一二三四社区在线视频社区8| 嫩草影院精品99| 中国美女看黄片| 性色av乱码一区二区三区2| 51午夜福利影视在线观看| 亚洲精品国产精品久久久不卡| 他把我摸到了高潮在线观看| 全区人妻精品视频| 老熟妇乱子伦视频在线观看| 在线十欧美十亚洲十日本专区| 一级毛片精品| 亚洲欧美日韩高清专用| 久久中文字幕人妻熟女| 亚洲熟妇熟女久久| 国产主播在线观看一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 成人高潮视频无遮挡免费网站| 日日摸夜夜添夜夜添小说| 中出人妻视频一区二区| 午夜精品在线福利| 午夜免费激情av| 国产亚洲精品久久久com| 日本a在线网址| 婷婷六月久久综合丁香| 日韩欧美国产一区二区入口| 成年女人看的毛片在线观看| 成人永久免费在线观看视频| 一本精品99久久精品77| 成年女人毛片免费观看观看9| 国产又黄又爽又无遮挡在线| 精品欧美国产一区二区三| 男女视频在线观看网站免费| 亚洲,欧美精品.| 日本与韩国留学比较| 91字幕亚洲| 精品电影一区二区在线| 中文在线观看免费www的网站| 久久草成人影院| 日本 欧美在线| 中亚洲国语对白在线视频| 九九热线精品视视频播放| av天堂在线播放| ponron亚洲| 又紧又爽又黄一区二区| 国产精品野战在线观看| 欧美3d第一页| 色视频www国产| 亚洲五月婷婷丁香| a在线观看视频网站| 午夜福利免费观看在线| 久久久久久久久久黄片| 亚洲av免费在线观看| 日韩欧美在线乱码| 麻豆一二三区av精品| 法律面前人人平等表现在哪些方面| 国产精品久久视频播放| 国产黄色小视频在线观看| 丰满的人妻完整版| 中文字幕人妻丝袜一区二区| 日韩成人在线观看一区二区三区| 无遮挡黄片免费观看| 亚洲中文字幕日韩| 国产av一区在线观看免费| 欧美丝袜亚洲另类 | 在线播放国产精品三级| 国产单亲对白刺激| 亚洲精品456在线播放app | 啦啦啦免费观看视频1| 国产av一区在线观看免费| 国产成人精品无人区| а√天堂www在线а√下载| 91在线精品国自产拍蜜月 | 国产亚洲精品av在线| 制服人妻中文乱码| 老司机午夜十八禁免费视频| 精品熟女少妇八av免费久了| 国产淫片久久久久久久久 | 十八禁网站免费在线| 无遮挡黄片免费观看| 亚洲av片天天在线观看| 中文字幕最新亚洲高清| 99热这里只有精品一区 | 麻豆av在线久日| 美女黄网站色视频| 国内揄拍国产精品人妻在线| 日韩 欧美 亚洲 中文字幕| 亚洲国产看品久久| 亚洲真实伦在线观看| 日本五十路高清| www.自偷自拍.com| 欧美三级亚洲精品| 国产一区二区在线观看日韩 | 成人性生交大片免费视频hd| 性色avwww在线观看| 欧美成狂野欧美在线观看| e午夜精品久久久久久久| 狂野欧美白嫩少妇大欣赏| 1024手机看黄色片| 99精品久久久久人妻精品| 热99re8久久精品国产| 成年女人毛片免费观看观看9| 亚洲 欧美一区二区三区| h日本视频在线播放| 国产三级中文精品| 国产高清激情床上av| 综合色av麻豆| 亚洲av成人av| 少妇丰满av| 亚洲精品色激情综合| 欧美午夜高清在线| www.999成人在线观看| 亚洲国产精品999在线| 国产一区在线观看成人免费| 18禁国产床啪视频网站| 午夜福利免费观看在线| 免费大片18禁| 毛片女人毛片| 99在线视频只有这里精品首页| 禁无遮挡网站| 日本三级黄在线观看| 国产精华一区二区三区| 亚洲欧美日韩卡通动漫| 麻豆av在线久日| 在线观看免费午夜福利视频| 色视频www国产| 最新在线观看一区二区三区| 欧美色欧美亚洲另类二区| 亚洲一区高清亚洲精品| 蜜桃久久精品国产亚洲av| 亚洲欧洲精品一区二区精品久久久| 一区二区三区国产精品乱码| 欧美乱码精品一区二区三区| 亚洲欧美日韩东京热| 黑人巨大精品欧美一区二区mp4| 999精品在线视频| 女同久久另类99精品国产91| 亚洲成人久久性| 日韩成人在线观看一区二区三区| 久久中文字幕人妻熟女| 最近在线观看免费完整版| 黄色丝袜av网址大全| 日日夜夜操网爽| 一本综合久久免费| 人妻夜夜爽99麻豆av| 日本与韩国留学比较| av天堂在线播放| 免费电影在线观看免费观看| 精品人妻1区二区| 欧美av亚洲av综合av国产av| 婷婷六月久久综合丁香| 久久久久久久午夜电影| 叶爱在线成人免费视频播放| 久久中文字幕人妻熟女| 亚洲aⅴ乱码一区二区在线播放| 国产高清有码在线观看视频| 黄色片一级片一级黄色片| 国产aⅴ精品一区二区三区波| 99久久精品一区二区三区| 韩国av一区二区三区四区| 亚洲午夜精品一区,二区,三区| 日韩欧美 国产精品| 亚洲美女黄片视频| 天天躁狠狠躁夜夜躁狠狠躁| 人人妻,人人澡人人爽秒播| 在线观看一区二区三区| 色吧在线观看| 在线视频色国产色| 国产极品精品免费视频能看的| 精品国产超薄肉色丝袜足j| 国产av在哪里看| 99国产极品粉嫩在线观看| 男女床上黄色一级片免费看| 精华霜和精华液先用哪个| 熟妇人妻久久中文字幕3abv| 熟妇人妻久久中文字幕3abv| 欧美绝顶高潮抽搐喷水| 亚洲色图av天堂| www.自偷自拍.com| 99re在线观看精品视频| 少妇人妻一区二区三区视频| 午夜视频精品福利| 国产午夜福利久久久久久| 日韩欧美精品v在线| 午夜福利高清视频| 麻豆成人av在线观看| 一区二区三区激情视频| 久久久久国产精品人妻aⅴ院| 中文资源天堂在线| 国产亚洲av高清不卡| 亚洲中文字幕日韩| 老司机在亚洲福利影院| 日本免费一区二区三区高清不卡| 久久欧美精品欧美久久欧美| 看黄色毛片网站| 国产精品亚洲美女久久久| 亚洲欧美日韩无卡精品| 美女黄网站色视频| 免费观看的影片在线观看| 久久久色成人| 亚洲精品456在线播放app | 久久人妻av系列| 成年女人毛片免费观看观看9| 又爽又黄无遮挡网站| 成人18禁在线播放| 一个人看视频在线观看www免费 | 日本 欧美在线| 欧美日韩一级在线毛片| 又爽又黄无遮挡网站| 精品国产乱子伦一区二区三区| 欧美乱色亚洲激情| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲中文字幕一区二区三区有码在线看 | 亚洲熟妇中文字幕五十中出| 欧美日韩精品网址| 久久久久亚洲av毛片大全| 国产私拍福利视频在线观看| or卡值多少钱| 午夜激情福利司机影院| x7x7x7水蜜桃| aaaaa片日本免费| 国内精品久久久久久久电影| 成年免费大片在线观看| 久久草成人影院| 最近最新中文字幕大全电影3| 国内久久婷婷六月综合欲色啪| 高潮久久久久久久久久久不卡| 久久久精品大字幕| 日日夜夜操网爽| 亚洲色图av天堂| 亚洲成人精品中文字幕电影| 最近视频中文字幕2019在线8| 夜夜爽天天搞| 亚洲第一欧美日韩一区二区三区| 99久久综合精品五月天人人| 亚洲av中文字字幕乱码综合| 男人舔女人的私密视频| 深夜精品福利| 亚洲国产欧美一区二区综合| 成年女人毛片免费观看观看9| 国产综合懂色| 观看免费一级毛片| 免费观看精品视频网站| 黑人操中国人逼视频| 亚洲在线观看片| 久久久久免费精品人妻一区二区| 18禁黄网站禁片午夜丰满| 最新在线观看一区二区三区| 国产高清videossex| 国产高清videossex| 精品久久久久久久人妻蜜臀av| 好看av亚洲va欧美ⅴa在| 亚洲 欧美 日韩 在线 免费| 国产精品久久久久久亚洲av鲁大| 色在线成人网| 亚洲美女黄片视频| 亚洲国产欧洲综合997久久,| 国产精品久久久久久人妻精品电影| 亚洲国产欧美一区二区综合| 毛片女人毛片| 亚洲精品456在线播放app | 欧美一区二区精品小视频在线| av天堂在线播放| 亚洲激情在线av| 亚洲性夜色夜夜综合| 老鸭窝网址在线观看| 老鸭窝网址在线观看| 美女 人体艺术 gogo| 国产av一区在线观看免费| 国产成人精品久久二区二区免费| 老司机在亚洲福利影院| 午夜福利在线在线| 国产精品自产拍在线观看55亚洲| 国产美女午夜福利| 99精品在免费线老司机午夜| 国产精品99久久99久久久不卡| av欧美777| 国产一区二区在线av高清观看| 亚洲五月婷婷丁香| 夜夜爽天天搞| 亚洲第一欧美日韩一区二区三区| 高清在线国产一区| 国产又黄又爽又无遮挡在线| 香蕉久久夜色| 国产精品98久久久久久宅男小说| av福利片在线观看| 国产高清视频在线观看网站| 亚洲精品粉嫩美女一区| 嫩草影院入口| 天堂网av新在线| 国产精品98久久久久久宅男小说| 搞女人的毛片| 久久这里只有精品中国| 噜噜噜噜噜久久久久久91| 99riav亚洲国产免费| 在线观看舔阴道视频| 亚洲第一欧美日韩一区二区三区| 亚洲国产欧美网| 不卡一级毛片| 国产伦人伦偷精品视频| 女人被狂操c到高潮| 美女扒开内裤让男人捅视频| 精品久久久久久久人妻蜜臀av| 亚洲av片天天在线观看| 国产精品亚洲美女久久久| 两性午夜刺激爽爽歪歪视频在线观看| 999久久久精品免费观看国产| 久久国产乱子伦精品免费另类| 国产精品日韩av在线免费观看| 国产精品美女特级片免费视频播放器 | 午夜免费观看网址| 亚洲五月婷婷丁香| 精品国内亚洲2022精品成人| 黑人巨大精品欧美一区二区mp4| 岛国视频午夜一区免费看| 欧美成人一区二区免费高清观看 | 亚洲精品乱码久久久v下载方式 | 噜噜噜噜噜久久久久久91| 久9热在线精品视频| 免费一级毛片在线播放高清视频| 男人舔女人的私密视频| 亚洲av日韩精品久久久久久密| 在线观看日韩欧美| 两个人看的免费小视频| 午夜福利在线观看免费完整高清在 | 天堂av国产一区二区熟女人妻| 精品人妻1区二区| 国产淫片久久久久久久久 | 日韩精品中文字幕看吧| 亚洲国产精品999在线| 日本a在线网址| 亚洲欧美精品综合久久99| 国产精品98久久久久久宅男小说| 琪琪午夜伦伦电影理论片6080| 欧美日韩一级在线毛片| 成人精品一区二区免费| 国产精品1区2区在线观看.| 国产精品日韩av在线免费观看| av女优亚洲男人天堂 | 午夜福利18| 久久久国产精品麻豆| 中出人妻视频一区二区| 欧美三级亚洲精品| 九色成人免费人妻av| 熟女人妻精品中文字幕| 美女被艹到高潮喷水动态| 黄色成人免费大全| 久久久水蜜桃国产精品网| 免费人成视频x8x8入口观看| 日韩国内少妇激情av| 久久久久久久午夜电影| 国产麻豆成人av免费视频| 又黄又爽又免费观看的视频| 宅男免费午夜| 国产亚洲av嫩草精品影院| 国产亚洲精品一区二区www| 美女大奶头视频| 亚洲黑人精品在线| 日日夜夜操网爽| 精品一区二区三区视频在线观看免费| 亚洲av成人av| 欧美色欧美亚洲另类二区| 中文字幕最新亚洲高清| 精品国产超薄肉色丝袜足j| 视频区欧美日本亚洲| 高清在线国产一区| 亚洲中文字幕一区二区三区有码在线看 | 国产成人精品无人区| 成年免费大片在线观看| 国产综合懂色| 麻豆久久精品国产亚洲av| 女生性感内裤真人,穿戴方法视频| 婷婷丁香在线五月| 黄色成人免费大全| 免费av毛片视频| 久久天躁狠狠躁夜夜2o2o| 中文字幕久久专区| 亚洲国产欧美网| 亚洲一区二区三区色噜噜| 国产成年人精品一区二区| 国产欧美日韩一区二区精品| 无限看片的www在线观看| 88av欧美| 99久久精品国产亚洲精品| 亚洲av五月六月丁香网| 巨乳人妻的诱惑在线观看| 男插女下体视频免费在线播放| xxx96com| 色在线成人网| 欧美日韩黄片免| 一进一出抽搐gif免费好疼| 日韩三级视频一区二区三区| 精品久久久久久成人av| 婷婷亚洲欧美| 精品国产三级普通话版| 日韩有码中文字幕| 欧美日本视频| 男女午夜视频在线观看| 69av精品久久久久久| 两个人看的免费小视频| av中文乱码字幕在线| 久久久久久久精品吃奶| 91在线观看av| 男女午夜视频在线观看| 免费搜索国产男女视频| 国产成人精品久久二区二区免费| 日本与韩国留学比较| 久久亚洲真实| 伦理电影免费视频| 日本成人三级电影网站| 男女床上黄色一级片免费看| 可以在线观看的亚洲视频| 国产午夜精品论理片| 97人妻精品一区二区三区麻豆| 亚洲无线在线观看| 在线观看免费午夜福利视频| 午夜精品一区二区三区免费看| 九九在线视频观看精品| 欧美黄色片欧美黄色片| 国产成人一区二区三区免费视频网站| 在线观看舔阴道视频| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久大精品| x7x7x7水蜜桃| 国产精品电影一区二区三区| 成年女人看的毛片在线观看| 国产高清视频在线观看网站| 国产精品亚洲av一区麻豆| 99热6这里只有精品| 国产精品亚洲av一区麻豆| 国产乱人伦免费视频| 丰满人妻一区二区三区视频av | 美女cb高潮喷水在线观看 | 欧美中文日本在线观看视频| 午夜日韩欧美国产| 村上凉子中文字幕在线| 性色av乱码一区二区三区2| 美女扒开内裤让男人捅视频| 亚洲电影在线观看av| 久久久国产精品麻豆| 亚洲成人久久爱视频| 欧美高清成人免费视频www| 丁香欧美五月| 黄片大片在线免费观看| www日本在线高清视频| 亚洲,欧美精品.| 精品一区二区三区视频在线 | 久久久精品欧美日韩精品| 免费观看精品视频网站| 巨乳人妻的诱惑在线观看| АⅤ资源中文在线天堂| 午夜精品久久久久久毛片777| 成人永久免费在线观看视频| 国产三级中文精品| 嫁个100分男人电影在线观看| 亚洲一区二区三区不卡视频| 亚洲,欧美精品.| 中文字幕最新亚洲高清| 亚洲 欧美 日韩 在线 免费| 偷拍熟女少妇极品色| 精品日产1卡2卡| 国产99白浆流出| 男女视频在线观看网站免费| 中文字幕久久专区| 亚洲欧美日韩无卡精品| 久久久久九九精品影院| 看片在线看免费视频| 99国产极品粉嫩在线观看| 桃红色精品国产亚洲av| 老熟妇乱子伦视频在线观看| 亚洲,欧美精品.| 国产v大片淫在线免费观看| 日韩欧美 国产精品| 久久久久久久精品吃奶| 日本熟妇午夜| 色精品久久人妻99蜜桃| 亚洲电影在线观看av| 给我免费播放毛片高清在线观看| 免费av不卡在线播放| 麻豆av在线久日| 婷婷精品国产亚洲av| 亚洲七黄色美女视频| 日本三级黄在线观看| 97碰自拍视频| 久久性视频一级片| 女同久久另类99精品国产91| 亚洲中文av在线| 欧美3d第一页| 特级一级黄色大片| 日韩成人在线观看一区二区三区| 免费搜索国产男女视频| 免费在线观看日本一区| 18禁国产床啪视频网站| 99精品在免费线老司机午夜|