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

    美女xxoo啪啪120秒动态图| 热re99久久精品国产66热6| 国产69精品久久久久777片| 国产一级毛片在线| 十八禁高潮呻吟视频| 丝袜喷水一区| 国产日韩欧美视频二区| 自线自在国产av| 亚洲综合色惰| 女人精品久久久久毛片| 多毛熟女@视频| 黄色毛片三级朝国网站| 大香蕉97超碰在线| 亚洲欧美一区二区三区国产| 欧美成人午夜免费资源| 精品人妻在线不人妻| 大香蕉久久网| 人成视频在线观看免费观看| 69精品国产乱码久久久| 精品久久蜜臀av无| av天堂久久9| 黄色毛片三级朝国网站| 国产免费视频播放在线视频| 国产 精品1| 精品卡一卡二卡四卡免费| 精品一区二区三区视频在线| 秋霞在线观看毛片| 久久韩国三级中文字幕| 国产精品嫩草影院av在线观看| 一边摸一边做爽爽视频免费| 九色亚洲精品在线播放| 亚洲精品一区蜜桃| 国产av一区二区精品久久| 国产白丝娇喘喷水9色精品| 中文字幕人妻丝袜制服| 欧美日本中文国产一区发布| 满18在线观看网站| 日产精品乱码卡一卡2卡三| 亚洲av男天堂| 国产熟女午夜一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 免费观看无遮挡的男女| 久久久久国产精品人妻一区二区| 日韩 亚洲 欧美在线| 日韩制服丝袜自拍偷拍| 国产午夜精品一二区理论片| 成人免费观看视频高清| 18在线观看网站| 十分钟在线观看高清视频www| 黄色怎么调成土黄色| 亚洲国产色片| 男人舔女人的私密视频| 午夜免费观看性视频| 亚洲精品美女久久久久99蜜臀 | 满18在线观看网站| 99热国产这里只有精品6| 777米奇影视久久| 免费观看无遮挡的男女| 麻豆乱淫一区二区| 下体分泌物呈黄色| 亚洲少妇的诱惑av| 日韩成人伦理影院| 免费观看在线日韩| 观看av在线不卡| 伦理电影免费视频| 五月天丁香电影| 99热这里只有是精品在线观看| 一个人免费看片子| 久久精品国产a三级三级三级| kizo精华| 国产免费福利视频在线观看| 精品国产乱码久久久久久小说| 亚洲一级一片aⅴ在线观看| 视频中文字幕在线观看| 国产精品熟女久久久久浪| 18禁动态无遮挡网站| 久热这里只有精品99| 久久久久久久久久人人人人人人| 黄片无遮挡物在线观看| 国产亚洲欧美精品永久| av国产精品久久久久影院| 水蜜桃什么品种好| 久久久久久久亚洲中文字幕| 成年美女黄网站色视频大全免费| 国产男人的电影天堂91| 男人爽女人下面视频在线观看| 成人综合一区亚洲| 纵有疾风起免费观看全集完整版| 九色成人免费人妻av| 曰老女人黄片| 日韩欧美一区视频在线观看| 一级片免费观看大全| 国产亚洲午夜精品一区二区久久| 王馨瑶露胸无遮挡在线观看| 99热这里只有是精品在线观看| 日韩一区二区视频免费看| 精品国产露脸久久av麻豆| 看非洲黑人一级黄片| 久久久久人妻精品一区果冻| 精品一区二区三区视频在线| 一个人免费看片子| 免费高清在线观看日韩| 欧美日韩视频精品一区| 国产高清国产精品国产三级| 成年人午夜在线观看视频| 高清毛片免费看| 亚洲,欧美,日韩| 久久99一区二区三区| 一级片免费观看大全| www.av在线官网国产| 一级毛片我不卡| 两个人看的免费小视频| 亚洲精品国产色婷婷电影| 国产成人av激情在线播放| 精品卡一卡二卡四卡免费| 欧美 亚洲 国产 日韩一| 久久久久国产精品人妻一区二区| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美日韩卡通动漫| 黄色毛片三级朝国网站| 成人漫画全彩无遮挡| 亚洲av.av天堂| 久久99一区二区三区| 大话2 男鬼变身卡| 边亲边吃奶的免费视频| 精品人妻偷拍中文字幕| 男女免费视频国产| 成人毛片60女人毛片免费| 亚洲av国产av综合av卡| 视频区图区小说| 视频中文字幕在线观看| 在线观看www视频免费| 国产探花极品一区二区| 汤姆久久久久久久影院中文字幕| 成人国产麻豆网| av国产精品久久久久影院| 成人二区视频| 午夜福利,免费看| 免费人妻精品一区二区三区视频| 多毛熟女@视频| 曰老女人黄片| 精品一品国产午夜福利视频| 大话2 男鬼变身卡| 在线观看三级黄色| 一本大道久久a久久精品| 韩国精品一区二区三区 | 成人国语在线视频| 18禁在线无遮挡免费观看视频| 自线自在国产av| 色94色欧美一区二区| 99热国产这里只有精品6| 亚洲,一卡二卡三卡| 亚洲精品456在线播放app| 亚洲性久久影院| 国产精品国产三级专区第一集| 国产视频首页在线观看| 97精品久久久久久久久久精品| 亚洲三级黄色毛片| 香蕉丝袜av| 美女国产视频在线观看| 黄色毛片三级朝国网站| 久久久久久人人人人人| 日本黄色日本黄色录像| 亚洲色图 男人天堂 中文字幕 | 国产精品久久久久久精品古装| 国产av一区二区精品久久| 欧美精品一区二区大全| 我的女老师完整版在线观看| 国产亚洲欧美精品永久| a 毛片基地| 久久这里有精品视频免费| 日本黄大片高清| 国产精品人妻久久久影院| 美国免费a级毛片| 国产亚洲精品第一综合不卡 | xxxhd国产人妻xxx| 黄色 视频免费看| 青春草视频在线免费观看| 男女免费视频国产| 少妇的逼好多水| 国产片特级美女逼逼视频| 男女无遮挡免费网站观看| 成人毛片60女人毛片免费| 亚洲综合色网址| 国产精品久久久久久精品电影小说| 日韩电影二区| 日韩人妻精品一区2区三区| 亚洲综合色网址| 五月天丁香电影| 免费人成在线观看视频色| 伊人亚洲综合成人网| √禁漫天堂资源中文www| 狂野欧美激情性bbbbbb| 久久av网站| 久久99蜜桃精品久久| 最近中文字幕高清免费大全6| 女人被躁到高潮嗷嗷叫费观| 极品少妇高潮喷水抽搐| 日韩一区二区视频免费看| 热99久久久久精品小说推荐| 亚洲成av片中文字幕在线观看 | 亚洲,一卡二卡三卡| 丝袜美足系列| 欧美 亚洲 国产 日韩一| 女性被躁到高潮视频| 亚洲精品日本国产第一区| 如何舔出高潮| 在线观看www视频免费| 成年人免费黄色播放视频| 中文欧美无线码| 在线观看美女被高潮喷水网站| 久久久精品区二区三区| 亚洲国产日韩一区二区| 精品一区二区免费观看| 大香蕉97超碰在线| av在线app专区| 免费观看在线日韩| www.熟女人妻精品国产 | 日韩欧美一区视频在线观看| 蜜臀久久99精品久久宅男| 中国国产av一级| 99热全是精品| 国产高清三级在线| av福利片在线| 国产成人免费观看mmmm| 人人妻人人添人人爽欧美一区卜| 免费日韩欧美在线观看| 国产精品一二三区在线看| 成人午夜精彩视频在线观看| 寂寞人妻少妇视频99o| 桃花免费在线播放| 咕卡用的链子| 日本91视频免费播放| 夫妻性生交免费视频一级片| 久久午夜综合久久蜜桃| av在线老鸭窝| av福利片在线| 少妇 在线观看| 欧美日韩综合久久久久久| 久久ye,这里只有精品| 高清在线视频一区二区三区| 在线观看免费日韩欧美大片| 中文乱码字字幕精品一区二区三区| 国产在线免费精品| 欧美精品av麻豆av| 一个人免费看片子| 看非洲黑人一级黄片| 久久精品夜色国产| 精品熟女少妇av免费看| 久久精品久久精品一区二区三区| 考比视频在线观看| www.色视频.com| 秋霞在线观看毛片| 在线看a的网站| 国产免费视频播放在线视频| 18+在线观看网站| av片东京热男人的天堂| 香蕉精品网在线| 91午夜精品亚洲一区二区三区| 人妻人人澡人人爽人人| 黄色配什么色好看| 激情五月婷婷亚洲| 少妇精品久久久久久久| 日本午夜av视频| 最近的中文字幕免费完整| 在线天堂最新版资源| 黄色怎么调成土黄色| 久热这里只有精品99| 女人精品久久久久毛片| 三上悠亚av全集在线观看| 18在线观看网站| www日本在线高清视频| 国产成人精品福利久久| 狂野欧美激情性xxxx在线观看| 国产成人精品一,二区| 少妇被粗大猛烈的视频| 国产精品国产三级国产专区5o| 男女免费视频国产| 9191精品国产免费久久| 亚洲国产日韩一区二区| 制服人妻中文乱码| 夜夜爽夜夜爽视频| 桃花免费在线播放| 中文字幕免费在线视频6| 蜜臀久久99精品久久宅男| 成年人午夜在线观看视频| 伦理电影大哥的女人| 少妇人妻久久综合中文| a级毛片黄视频| 久久国产精品大桥未久av| 久久99热这里只频精品6学生| 精品视频人人做人人爽| 蜜臀久久99精品久久宅男| 91精品国产国语对白视频| 两个人免费观看高清视频| 这个男人来自地球电影免费观看 | 国产亚洲午夜精品一区二区久久| 精品亚洲成国产av| 国产一区二区在线观看日韩| 欧美丝袜亚洲另类| 日韩一本色道免费dvd| 一本久久精品| 人人妻人人澡人人爽人人夜夜| 欧美性感艳星| 午夜视频国产福利| 高清av免费在线| 乱码一卡2卡4卡精品| 少妇被粗大猛烈的视频| 免费看av在线观看网站| 亚洲激情五月婷婷啪啪| 在现免费观看毛片| 成人亚洲欧美一区二区av| 精品少妇久久久久久888优播| 夜夜爽夜夜爽视频| 18禁裸乳无遮挡动漫免费视频| 久久久久久久国产电影| 黄片无遮挡物在线观看| 在线免费观看不下载黄p国产| 国产免费一级a男人的天堂| 人妻人人澡人人爽人人| 欧美最新免费一区二区三区| 久久国产亚洲av麻豆专区| 精品一品国产午夜福利视频| videosex国产| 欧美成人精品欧美一级黄| 伦精品一区二区三区| 亚洲熟女精品中文字幕| 卡戴珊不雅视频在线播放| 免费黄网站久久成人精品| 人体艺术视频欧美日本| 日日撸夜夜添| 曰老女人黄片| 久久久久久久久久成人| 日本wwww免费看| 高清欧美精品videossex| 男女无遮挡免费网站观看| 国产极品天堂在线| 免费观看性生交大片5| 亚洲精品一区蜜桃| 男女啪啪激烈高潮av片| 国产成人免费无遮挡视频| 涩涩av久久男人的天堂| 亚洲激情五月婷婷啪啪| 视频中文字幕在线观看| 国产乱人偷精品视频| 久久久亚洲精品成人影院| 国产日韩一区二区三区精品不卡| 成年女人在线观看亚洲视频| av片东京热男人的天堂| av在线app专区| 国产一区二区在线观看av| 在线亚洲精品国产二区图片欧美| 欧美日韩视频精品一区| 婷婷成人精品国产| 亚洲美女黄色视频免费看| 久久 成人 亚洲| 久久鲁丝午夜福利片| 免费观看性生交大片5| 国产探花极品一区二区| xxxhd国产人妻xxx| 看免费成人av毛片| 91在线精品国自产拍蜜月| 精品人妻熟女毛片av久久网站| 国国产精品蜜臀av免费| 人妻 亚洲 视频| 国国产精品蜜臀av免费| 久久毛片免费看一区二区三区| 国产男女内射视频| 精品国产一区二区三区四区第35| 高清av免费在线| 国产淫语在线视频| √禁漫天堂资源中文www| 日韩一区二区视频免费看| 精品久久国产蜜桃| 久久这里有精品视频免费| 夫妻午夜视频| 日韩精品免费视频一区二区三区 | 97人妻天天添夜夜摸| 精品人妻偷拍中文字幕| 亚洲欧美日韩另类电影网站| 亚洲精品视频女| 国产老妇伦熟女老妇高清| 国产av精品麻豆| 综合色丁香网| 国产欧美亚洲国产| 丝袜美足系列| 熟女电影av网| 欧美人与性动交α欧美软件 | 国产片特级美女逼逼视频| 男女免费视频国产| 国产成人aa在线观看| 国产免费一级a男人的天堂| 亚洲欧美中文字幕日韩二区| 一级片免费观看大全| 少妇人妻精品综合一区二区| 亚洲精品自拍成人| 狂野欧美激情性bbbbbb| 一区二区三区精品91| 在线观看www视频免费| 亚洲色图 男人天堂 中文字幕 | 亚洲精品视频女| 国产日韩欧美视频二区| 国产精品一国产av| 久久人人97超碰香蕉20202| 菩萨蛮人人尽说江南好唐韦庄| 国产日韩欧美亚洲二区| 日韩中文字幕视频在线看片| 久久精品熟女亚洲av麻豆精品| 午夜久久久在线观看| 热99久久久久精品小说推荐| 亚洲情色 制服丝袜| 色5月婷婷丁香| 在线观看美女被高潮喷水网站| 国产一区二区三区综合在线观看 | 国产福利在线免费观看视频| 最近的中文字幕免费完整| av在线观看视频网站免费| 欧美丝袜亚洲另类| 国产欧美日韩综合在线一区二区| 精品人妻一区二区三区麻豆| 国产精品人妻久久久影院| 免费人成在线观看视频色| 国产片内射在线| 亚洲性久久影院| 高清不卡的av网站| 99精国产麻豆久久婷婷| a级毛片黄视频| 亚洲精品乱码久久久久久按摩| 亚洲精品第二区| 黄片无遮挡物在线观看| 九九爱精品视频在线观看| 女人精品久久久久毛片| 男人添女人高潮全过程视频| 久久久久人妻精品一区果冻| 国产欧美亚洲国产| 亚洲综合精品二区| 人成视频在线观看免费观看| 韩国高清视频一区二区三区| 精品国产乱码久久久久久小说| 十八禁高潮呻吟视频| 777米奇影视久久| 亚洲三级黄色毛片| 国产日韩一区二区三区精品不卡| 亚洲国产日韩一区二区| 综合色丁香网| 成人午夜精彩视频在线观看| 一本—道久久a久久精品蜜桃钙片| 啦啦啦中文免费视频观看日本| 国产一级毛片在线| av播播在线观看一区| 韩国高清视频一区二区三区| 蜜臀久久99精品久久宅男| 日韩 亚洲 欧美在线| 91午夜精品亚洲一区二区三区| 国产精品99久久99久久久不卡 | 男男h啪啪无遮挡| 欧美日韩亚洲高清精品| 大香蕉久久成人网| 免费日韩欧美在线观看| 精品人妻在线不人妻| 国产成人精品婷婷| 成年女人在线观看亚洲视频| 国产精品三级大全| 午夜免费观看性视频| 又粗又硬又长又爽又黄的视频| 伊人久久国产一区二区| 日本wwww免费看| 亚洲伊人色综图| 天天躁夜夜躁狠狠久久av| 国产亚洲精品第一综合不卡 | 成年av动漫网址| 日韩伦理黄色片| 超碰97精品在线观看| 日韩人妻精品一区2区三区| 老司机亚洲免费影院| 伦理电影大哥的女人| 日韩成人av中文字幕在线观看| 国产伦理片在线播放av一区| 久久热在线av| 赤兔流量卡办理| 亚洲欧洲精品一区二区精品久久久 | 亚洲欧美清纯卡通| 欧美+日韩+精品| 免费人妻精品一区二区三区视频| 国产极品粉嫩免费观看在线| 中文乱码字字幕精品一区二区三区| 黑丝袜美女国产一区| 波野结衣二区三区在线| 国产激情久久老熟女| 91成人精品电影| 美女主播在线视频| 熟女人妻精品中文字幕| 高清黄色对白视频在线免费看| 99久久精品国产国产毛片| 制服丝袜香蕉在线| 久久99热这里只频精品6学生| 人妻 亚洲 视频| 国产亚洲欧美精品永久| 国产亚洲精品第一综合不卡 | 国产黄色视频一区二区在线观看| 国产免费一级a男人的天堂| 亚洲国产最新在线播放| 五月开心婷婷网| 久久久久久久精品精品| 亚洲久久久国产精品| 亚洲国产精品一区三区| 不卡视频在线观看欧美| 亚洲色图综合在线观看| 少妇高潮的动态图| av在线播放精品| 亚洲精品美女久久久久99蜜臀 | 丁香六月天网| 中文精品一卡2卡3卡4更新| 在线观看免费视频网站a站| 亚洲精品国产色婷婷电影| 久久久欧美国产精品| 国国产精品蜜臀av免费| 熟女av电影| 人人澡人人妻人| 少妇的逼水好多| 成人影院久久| 黑丝袜美女国产一区| 精品少妇久久久久久888优播| 久久久久久久久久成人| 毛片一级片免费看久久久久| 制服丝袜香蕉在线| 国产成人aa在线观看| 成人毛片60女人毛片免费| 一二三四在线观看免费中文在 | 日本vs欧美在线观看视频| 国产免费视频播放在线视频| 午夜免费男女啪啪视频观看| 国产熟女午夜一区二区三区| 亚洲精品日本国产第一区| 女性生殖器流出的白浆| 丰满少妇做爰视频| 满18在线观看网站| av女优亚洲男人天堂| 国产成人精品无人区| 夜夜骑夜夜射夜夜干| 欧美xxⅹ黑人| 水蜜桃什么品种好| 久久99蜜桃精品久久| 女的被弄到高潮叫床怎么办| 精品国产乱码久久久久久小说| 新久久久久国产一级毛片| 一级a做视频免费观看| 深夜精品福利| 免费av中文字幕在线| 欧美+日韩+精品| 国产1区2区3区精品| 大陆偷拍与自拍| 亚洲国产精品成人久久小说| 视频在线观看一区二区三区| 久久久精品区二区三区| 日韩制服骚丝袜av| 日本欧美视频一区| 国产精品国产三级国产av玫瑰| 国产日韩欧美在线精品| 日本猛色少妇xxxxx猛交久久| 人人妻人人爽人人添夜夜欢视频| 亚洲欧美日韩卡通动漫| 成人亚洲欧美一区二区av| 一本大道久久a久久精品| 精品熟女少妇av免费看| 国产精品嫩草影院av在线观看| 最新的欧美精品一区二区| 另类精品久久| 少妇被粗大猛烈的视频| 国产亚洲精品第一综合不卡 | 久久国产精品男人的天堂亚洲 | 乱码一卡2卡4卡精品| 一边亲一边摸免费视频| 欧美激情 高清一区二区三区| 高清不卡的av网站| 日本-黄色视频高清免费观看| 美女国产高潮福利片在线看| 亚洲av日韩在线播放| 天堂俺去俺来也www色官网| 99香蕉大伊视频| 亚洲少妇的诱惑av| 国产国拍精品亚洲av在线观看| 秋霞伦理黄片| 国产高清国产精品国产三级| 欧美老熟妇乱子伦牲交| 视频中文字幕在线观看| 久久久久国产网址| 在线观看免费高清a一片| 亚洲成色77777| 9191精品国产免费久久| 国产有黄有色有爽视频| av不卡在线播放| 伊人亚洲综合成人网| 秋霞在线观看毛片| 制服诱惑二区| 18禁观看日本| 国产精品秋霞免费鲁丝片| 日本91视频免费播放| 又粗又硬又长又爽又黄的视频| 久热这里只有精品99| 丰满迷人的少妇在线观看| 九草在线视频观看| 免费日韩欧美在线观看| 一级,二级,三级黄色视频| 热re99久久精品国产66热6| 制服诱惑二区| 97在线人人人人妻| 亚洲内射少妇av| 欧美+日韩+精品| 高清不卡的av网站| 免费女性裸体啪啪无遮挡网站| 国产成人aa在线观看| 日本午夜av视频| 国产日韩欧美在线精品| 久久热在线av|