modeling)oxidant)chemistry)in)a)high3 isoprene)atmospheric...

12
Modeling oxidant chemistry in a high3 isoprene atmospheric environment: sensi7vity to model resolu7on and constraints from aircra: observa7ons Karen Yu, Daniel Jacob, Jenny Fisher, Sungshik Kim, Katherine Travis, Lei Zhu, Robert Yantosca, Melissa Sulprizio, Thomas Ryerson, Armin Wisthaler, Alan Fried, Paul Wennberg, and the SEAC 4 RS Science Team 7 May 2015 IGC7

Upload: others

Post on 20-Oct-2020

0 views

Category:

Documents


0 download

TRANSCRIPT

  • Modeling)oxidant)chemistry)in)a)high3isoprene)atmospheric)environment:)sensi7vity)to)model)resolu7on)and)

    constraints)from)aircra:)observa7ons!

    Karen!Yu,!Daniel!Jacob,!Jenny!Fisher,!Sungshik!Kim,!Katherine!Travis,!Lei!Zhu,!Robert!Yantosca,!Melissa!Sulprizio,!Thomas!Ryerson,!Armin!

    Wisthaler,!Alan!Fried,!Paul!Wennberg,!and!the!SEAC4RS!Science!Team!)

    7)May)2015)IGC7)

  • Computa7on)expense)of)tropospheric)chemistry)simula7ons)requires)evalua7on)of)benefits)of)increasing)grid)resolu7on)

    2!

    Resolu7on) Resource)requirements)(normalized)to)4x5))

    4ox5o!! 1!

    2ox2.5o! 8!

    0.25ox0.3125o! 1536!

  • Computa7on)expense)of)tropospheric)chemistry)simula7ons)requires)evalua7on)of)benefits)of)increasing)grid)resolu7on)

    2!

    •  Effect!of!grid!resoluTon!on!NOXXsaturated!vs.!NOXXlimited!regime!has!been!extensively!studied.!

    Resolu7on) Resource)requirements)(normalized)to)4x5))

    4ox5o!! 1!

    2ox2.5o! 8!

    0.25ox0.3125o! 1536!

    10!

    1!

    0.1!

    0.01!

    NOX!(ppb)!

    NOX!concentraTon!

    Ozone

    !prod

    ucTo

    n!

    NOX3saturated)HOX!loss:!NO2+OH!

    NOX3limited)HOX!loss:!HO2+HO2!

  • Computa7on)expense)of)tropospheric)chemistry)simula7ons)requires)evalua7on)of)benefits)of)increasing)grid)resolu7on)

    2!

    •  Effect!of!grid!resoluTon!on!NOXXsaturated!vs.!NOXXlimited!regime!has!been!extensively!studied.!

    •  As!NOX!emissions!decline,!it!becomes!increasingly!important!to!evaluate!model!ability!to!represent!interacTon!of!biogenic!and!anthropogenic!emissions.!

    Resolu7on) Resource)requirements)(normalized)to)4x5))

    4ox5o!! 1!

    2ox2.5o! 8!

    0.25ox0.3125o! 1536!

    10!

    1!

    0.1!

    0.01!

    NOX!(ppb)!

    High3NOX))RO2+NO!

    Low3NOX))RO2+HO2!

    How)does)grid)resolu7on)affect)model)representa7on)of)high3NOX)vs.)low3NOX)regimes?)

    VOC)+)OH)

    High!NOX:!nitrate!producTon,!ozone!producTon!

    Low!NOX:!chain!terminaTon,!sink!for!HOx!

    RO2)

  • Evaluate)GEOS3Chem)run)at)4ox5o,)2ox2.5o)and)0.25Ox0.3125o)resolu7ons)against)SEAC4RS)observa7ons)

    3!

    4x5!global,!starTng!January!2013!

    2x2.5!global,!start!January!2013!

    0.25x0.3125!N.!America!with!4x5!boundary!condiTons,!start!run!August!2013!

    GEOS!meteorology!at!0.25x0.3125!resoluTon!

  • Evaluate)GEOS3Chem)run)at)4ox5o,)2ox2.5o)and)0.25Ox0.3125o)resolu7ons)against)SEAC4RS)observa7ons)

    3!

    4x5!global,!starTng!January!2013!

    2x2.5!global,!start!January!2013!

    0.25x0.3125!N.!America!with!4x5!boundary!condiTons,!start!run!August!2013!

    GEOS!meteorology!at!0.25x0.3125!resoluTon!

    Instrument)(PI)) Measurement)

    CSD!CL!(T.!Ryerson)! Ozone!

    NO!

    NO2!PTRMS!(A.!Wisthaler)! Isoprene!

    CAMS!(A.!Fried)! Formaldehyde!

    Observa7ons)from)SEAC4RS)aircra:)campaign)

    Boundary!!layer!comparisons!

  • Emissions)of)NOX)and)isoprene)in)Southeast)US)show)high)granularity)and)spa7al)an73correla7on)

    YU ET AL.: IMPACT OF RESOLUTION ON TROP CHEM SIMULATIONS X - 15

    NOx emissions at 0.25x0.3125

    30oN

    35oN

    100oW 90oW 80oW

    0 30 60 90 120 molec/cm2/s * 1e9

    Isoprene (0.25x0.3125)

    30oN

    35oN

    100oW 90oW 80oW

    0.00 1.00 2.00 3.00 4.00 molec/cm2/s * 1e12

    Figure 1. Emissions of anthropogenic NOX, soil + fertilizer NOX, lightning NOX, and isoprene

    at each resolution for August 2013. Anthropogenic NOX emissions are from the NEI2008 inven-

    tory. Soil and fertilizer NOX are computed using the parameteriztation of Hudman et al. [2012].

    Lightning NOx computed based on satellite data as in Murray et al. [2012]. Isoprene emissions

    from MEGAN [Guenther et al., 2012].

    D R A F T April 24, 2015, 5:59pm D R A F T

    4!

  • Emissions)of)NOX)and)isoprene)in)Southeast)US)show)high)granularity)and)spa7al)an73correla7on)

    YU ET AL.: IMPACT OF RESOLUTION ON TROP CHEM SIMULATIONS X - 15

    NOx emissions at 0.25x0.3125

    30oN

    35oN

    100oW 90oW 80oW

    0 30 60 90 120 molec/cm2/s * 1e9

    Isoprene (0.25x0.3125)

    30oN

    35oN

    100oW 90oW 80oW

    0.00 1.00 2.00 3.00 4.00 molec/cm2/s * 1e12

    Figure 1. Emissions of anthropogenic NOX, soil + fertilizer NOX, lightning NOX, and isoprene

    at each resolution for August 2013. Anthropogenic NOX emissions are from the NEI2008 inven-

    tory. Soil and fertilizer NOX are computed using the parameteriztation of Hudman et al. [2012].

    Lightning NOx computed based on satellite data as in Murray et al. [2012]. Isoprene emissions

    from MEGAN [Guenther et al., 2012].

    D R A F T April 24, 2015, 5:59pm D R A F T

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●● ●

    ●●●

    ● ●

    ● ●

    ● ●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●●

    ● ●

    ●●●

    ● ●

    ●●

    ●●

    ●●●

    ●●

    ●●●

    ●●●

    ●●●

    ● ●●

    ●●

    ●● ●●●

    ●●

    ●●

    ●●

    ● ●●

    ●●

    ●●

    ●●

    ●● ●

    ●●

    ●●

    ●●

    ●●●●

    ●●●

    ●●●●

    ●●

    ●●●

    ●●

    ●●

    ●●●

    ●● ●●

    ●●

    ●● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●● ● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●● ●

    ● ●●

    ●●

    ● ●

    ●●

    ● ●●

    ●●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ●●

    ●●●●

    ●● ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●● ●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●●

    ●●

    ● ●

    ●●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ●● ●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ● ●●

    ●● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●●

    ●●

    ● ●

    ●●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ● ●●●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ● ●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●● ●

    ●●

    ● ●

    ● ●●●●

    ●●

    ●●

    ●●●

    ●●●

    ●●

    ●●

    ● ●●

    ● ●

    ●●●

    ●●

    ● ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ● ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ● ●

    r = −0.15

    1e+10

    1e+11

    1e+12

    1e+11 1e+12 1e+13Isoprene emiss. (molec cm−2s−1)

    NO e

    miss

    . (m

    olec

    cm−2

    s−1 )

    0.25x0.3125

    ●●●●●

    ●●

    ●●●●

    ● ●●

    ●●●●

    ●●

    ●●●●

    ●●

    ●● ●●●●●●

    ●●

    ●●●

    r = 0.10

    1e+10

    1e+11

    1e+12

    1e+11 1e+12 1e+13Isoprene emiss. (molec cm−2s−1)

    2x2.5

    ● ●●●

    ●●

    ●●●

    ●●●

    r = 0.76

    1e+10

    1e+11

    1e+12

    1e+11 1e+12 1e+13Isoprene emiss. (molec cm−2s−1)

    4x5

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●● ●

    ●●●

    ● ●

    ● ●

    ● ●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●●

    ● ●

    ●●●

    ● ●

    ●●

    ●●

    ●●●

    ●●

    ●●●

    ●●●

    ●●●

    ● ●●

    ●●

    ●● ●●●

    ●●

    ●●

    ●●

    ● ●●

    ●●

    ●●

    ●●

    ●● ●

    ●●

    ●●

    ●●

    ●●●●

    ●●●

    ●●●●

    ●●

    ●●●

    ●●

    ●●

    ●●●

    ●● ●●

    ●●

    ●● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●● ● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●● ●

    ● ●●

    ●●

    ● ●

    ●●

    ● ●●

    ●●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ●●

    ●●●●

    ●● ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●● ●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●●

    ●●

    ● ●

    ●●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ●● ●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ● ●●

    ●● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●●

    ●●

    ● ●

    ●●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ● ●●●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ● ●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●● ●

    ●●

    ● ●

    ● ●●●●

    ●●

    ●●

    ●●●

    ●●●

    ●●

    ●●

    ● ●●

    ● ●

    ●●●

    ●●

    ● ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ● ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●●

    ●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ● ●

    r = −0.15

    1e+10

    1e+11

    1e+12

    1e+11 1e+12 1e+13Isoprene emiss. (molec cm−2s−1)

    NO e

    miss

    . (m

    olec

    cm−2

    s−1 )

    0.25x0.3125

    ●●●●●

    ●●

    ●●●●

    ● ●●

    ●●●●

    ●●

    ●●●●

    ●●

    ●● ●●●●●●

    ●●

    ●●●

    r = 0.10

    1e+10

    1e+11

    1e+12

    1e+11 1e+12 1e+13Isoprene emiss. (molec cm−2s−1)

    2x2.5

    ● ●●●

    ●●

    ●●●

    ●●●

    r = 0.76

    1e+10

    1e+11

    1e+12

    1e+11 1e+12 1e+13Isoprene emiss. (molec cm−2s−1)

    4x5

  • Probability)distribu7ons)for)primary)and)secondary)species)

    5!

    •  Effect!of!grid!resoluTon!more!apparent!for!highly!skewed!distribuTons:!resoluTon!of!smallXscale!features!allows!for!more!extreme!tail!values!

    •  ResoluTon!differences!for!secondary!species!related!to!differences!in!chemical!pathways!

    ●●● ●

    ●●●

    ●●●●

    ● ●●

    ●●

    ●●

    ●●● ●

    ●●

    ●●

    ●●

    ●●

    ●● ●●

    ●●

    ●●

    ●●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●●●

    ●● ●

    ●●●

    ● ●

    ●●

    ●●●●●

    ●●

    ●●●

    ●●●●

    ●●●●●●

    ●●●●●

    ●●

    ●●●●

    ●●●●

    ●●

    ●●●● ●●

    ●●●●●●●

    ●●●

    ●●●●

    ●● ●●●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●●●

    ●●●

    ●●

    ● ●

    ●●

    ●●

    ●●

    ●●

    ●●●●●●●

    ●●●●●●

    ●●●

    ●●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●●●●●

    ●●●

    ●●●

    ●●

    ● ●● ●

    ●●●

    ●●

    ●●● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●●

    ●●●●● ●

    ●●●●

    ●●●●● ●●

    ●●

    ●●●

    ● ●●

    ●●

    ●●

    ●●

    ●●

    ●● ●●

    ●●●●

    ●●●●●●●

    ●●●

    ●●

    ●●●● ●

    ●●

    ●●●●

    ●●●●

    ● ●

    ●●●

    ●● ●

    ● ●●●●

    ●●●

    ●●●

    ●●

    ●●

    ●●●

    ●●

    ●●●

    ●● ●

    ●● ●●●

    ●●

    ●●●

    ●●●

    ●●●●

    ● ●

    ●●

    ●●●

    ●●●

    ●●

    ●●

    ●●

    ●●●● ●●

    ●●

    ●●●

    ●●

    ●●

    ● ●

    ● ●●

    ●●●●

    ●● ●

    ●●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●● ●●

    ●●●●

    ●●●

    ●●

    ●●

    ●● ●●●●

    ●●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●● ●

    ●●

    ●● ●

    ●●

    ●●

    ● ●●

    ●● ●

    ●●● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●●

    ●●●●

    ●●●●●●

    ●●

    ●●

    ●●

    ●●●●●

    ●●●

    ● ●●

    ●●

    ●●●

    ●●●

    ●●

    ●●●●●●●●

    ●●● ●●●●

    ●●●● ●

    ●●●

    ●●

    ●●●● ●●

    ●●

    ●●

    ●●●

    ●●● ●

    ●●●●●

    ●●

    ●●●●●●●

    ●●●

    ●●

    ●●

    ●●●●

    ●●● ●

    ●●●●●●●●●●

    ●●● ●

    ●●●

    ●●●

    ●●●●●●●●

    ●●●●●

    ●●

    ●●

    ●●●●●

    ●●●

    ●●

    ●●● ●●●

    ●●

    ●●●●●●

    ●● ●●

    ●●●

    ● ●●

    ●●●

    ●●

    ●●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●●

    ● ●●

    ●●

    ●●

    ●●

    ● ●● ●

    ●●●●

    ●●●●

    ●●●●●●

    ●●●●●●●●

    ●●

    ●●●●

    ●●●●

    ● ●●●

    ●●●

    ●●

    ●●

    ●●

    ●●●

    ●●●●

    ●●●●●

    ●●

    ●●

    ●●●●

    ●●

    ●●●●●●●●

    ●●

    ●●

    ●●●

    ●●●

    ●●●●

    ●●●●

    ●●●●●

    ●●

    ●●

    ●●

    ●●●●●

    ●●

    ●●

    ●●

    ●●●●●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●●●●●●●

    ●●●●

    ●●●●●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●●

    ●●●●

    ●●●

    ● ●

    ●●●●●

    ●●●

    ●● ●●

    ●●●●

    ●●

    ●●●

    ●●●●●●●

    ●●

    ●●

    ●● ●●

    ●●●

    ● ●

    ●● ●

    ●●

    ●●

    ● ●●

    ● ●●● ●●●

    ●●●●●●●●●

    ●●●

    ●●●●

    ●●●

    ●●●●

    ●●

    ●●

    ●●●●

    ●●

    ●●●

    ●●

    ●●●

    ●●● ●●

    ●●●

    ●●

    ●●●●

    ●●

    ●●

    ●●●

    ●●●●●●●

    ●●●

    ●●

    ●●

    ●●●

    ●●

    ●●●

    ●●●●●

    qqobs$x

    exp(

    qqob

    s$y)

    0.01

    0.1

    1

    10

    ●●●●●●●●●●●

    ●●●

    ●●●●

    ●●●●●●●●

    ●●●

    ●●●●●●

    ●●●●●●●

    ●●●●●●

    ●●●●●●●

    ●●●●●

    ●●●●●

    ●●

    ●●●●

    ●●●

    ●●

    ●●●

    ●●●●●●●●●●

    ●●●●

    ●●●●●●

    ● ●●●●●●●●●

    ●● ●●●●●●●

    ●●●●●●

    ●●●

    ●●●●●●●

    ●●●

    ● ●

    ●●●●●●●●●

    ●●●●●●●●●●

    ●●

    ●●●●●

    ●●●●●●

    ●●●

    ●●●●● ●●●●●

    ●●●

    ●●●●●●●●

    ●●●●●

    ●●

    ●● ●●●

    ●●●●●●●●●●●●

    ●●

    ●●

    ●●●●●●●●●

    ●●●●●

    ●●

    ● ●●

    ●●●●●●●●●

    ●●

    ●●

    ●●●● ●●●

    ●●●●

    ●●●●● ●

    ●●● ●●●●●●●●●●●●●●●●

    ●●●●●●●●●●●●●

    ●●●●●●●●●●●●●●●

    ●●●

    ●●●●●●●●●●●●●

    ●●●●●●●

    ●●●●●●● ●●●●

    ●●●●●● ●●●●

    ●●●●●●●●

    ●●●●●●●●

    ●●●●●●●●●

    ●●

    ●●

    ●●●●●●●●●●●●

    ●●●

    ●● ●●●●●●●●

    ●●●●●

    ●●

    ●●●

    ●●●

    ●●●●●●●

    ●●●●●

    ●●●●●●

    ●●●

    ●●●●

    ●●●

    ●●●●

    ●●●●●●●●●●●●●●

    ●●●●●●●●●●●

    ●●●

    ●●●●

    ●●

    ●●●●●●●

    ●●●

    ●●●●●●●●●

    ●●

    ●●●●

    ●●●●●●

    ●●

    ●●●

    ●●●

    ●●●●

    ●●●

    ●●●●●

    ●●

    ● ●●●

    ●●●●

    ●●●●

    ●●●

    ●●

    ●●●● ●●

    ●●●

    ●●

    ●●●●●●●

    ●●●●●●●●●●●

    ●●●●●●●●●●●●

    ●●●●●●●●●●●

    ●●●●●

    ●●●

    ●●●●●

    ●●

    ●●●●●● ●●●●●●

    ●●●●●●

    ●●●●●●●●●●●

    ●●●●●

    ● ●●

    ●●

    ●●●●●●●●●●

    ●●●●●●●●●●●●●

    ●●●●●

    ●●●●●

    ●●

    ●●●

    ●●●●

    ●●●●

    ●●●●

    ●●●●●●

    ●●●●

    ●●●●

    ●●●●

    ●●●●●●●●●●●

    ●●

    ●●

    ●●●

    ●●

    ●●●●

    ●●●

    ●●●

    ●●

    ●●●

    ●●●●●●●●●●

    ●●●●●●

    ●● ●

    ●●●

    ●●●●

    ●●

    ●●●●●

    ● ●

    ●●●●●●●●●●●●●●

    ●●●

    ●●●●●●●●●●●●●●

    ●●●● ●●

    ●●

    ●●●

    ●● ●

    ●●●●●

    ●●●●●●●

    ●●●●●

    ●●●●●●●●●●●●●●●●

    ●●●●●●●●●●● ●●●●●●●●

    ●●●●●●●●●

    ●●

    ●●

    ●●

    ●●●●

    ●●●●

    ●●●●●

    ● ●

    ●●●●

    ●●●● ●●●●

    ●●●●●●●●●●●●

    ●●●●

    ●●

    ●●●●●

    ●●●●

    ●● ●●●●●●●●●●

    ●●●●

    ●●

    ●● ●

    ● ● ●

    ●● ●

    ●●●● ●●●●●●●

    ●●●

    ●●

    ●●●●●●

    ●●●●●●●

    ●●●● ●●●●●●

    ●●●●●●●●●

    ●●

    ●●●●●●

    ●●●●●

    ●●●●●

    ●●●●

    ●●●●●●●

    ●●●●●●●

    ●●●●●●●

    ●●●●

    ●●●●

    ●● ●●●●●●●●●●

    ●●●

    ●●●●

    ●●●

    ●●●●●●

    ●●

    ●●●

    ●●

    ●●●

    ●●●●●●●● ●

    ●●●●●●

    ●●●●●● ●●●

    ●●●

    ●●●

    ●● ●●●

    ●●

    ●●●

    ●●●●●●●●●●●●●●● ●●●●●

    ●●●●●●●●●●●

    ●●●●

    ●●●●

    ●●●●●●

    ●●●●●●

    ●●●●●●●●●●●●●●●●●●●●●●●

    ●●

    ●●

    ●●●●●●●●●●●

    ●●●●●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●●●●●●●●●●

    ●●●●●●●●●●●●●●●●●●●

    ●●●●●●●●●●●●●●●

    ●●●●●●●●●●●●●

    ●● ●●●●●●●●●

    ●●

    ●●●●●●●

    ●●●●●●●●●●●●●●●●●●●●●

    ● ●●●●●●●

    ●●

    ●●●●●

    ●●

    ●● ●●●●●●●●

    ●●

    ●●

    ●●

    ●●●●●●● ●●●●●●

    ●●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●●●●●●●●●●●●●● ●●●●●●●●●●●

    ●●●●

    ●●●●●●●●●●●●●●●●●●●●●● ●●●●●●●●●●●●●●●●●●

    ●●●

    ●●●●●●●●●●●●

    ●●

    ●●

    ●●

    ●●

    ●●●●●●●●●

    ●●●●

    ●●

    ●● ●

    ●●●●●●●

    ●●●

    ●●

    ● ●●●●

    ●●●

    ●●●●●●●●●

    ●●●

    ●● ●●

    ●●●

    ●●●●●●●

    ●●●●●●●●●●●●●

    ●●● ● ●●●

    ●●●●●●●●●

    ●●

    ●●●●●●●●

    ●●

    ●●●●●●●●●●

    ●●

    ●●●●●●●●●●

    ●●●●●●●●●●

    ●●

    ●●●●●●●●●●●●●●●●●

    ●●●●●●●●●

    ●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●

    ● ●●●●●●●●●●●●●●●●● ●●

    ●●●●●●●●●●●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●●●●●●●●●●

    ●●●●●●●●●●●●●●

    ●●● ●●●●

    ●● ●●●●●●●

    ●●●●●●●●●●●●●●

    ● ●●●●●●●●●●●●

    ●●

    ●●

    ●●

    ●●●●●●●

    ●●●

    ●●

    ●●●●

    ●●●●●●●●●●●●●

    ●●●●●●●●●●●●

    ●● ●

    ●●●

    ●●●●●

    ●●

    ●●●●

    ●●●

    ●●●●●●●●●●●●●●●

    ●●●

    ●●●●●●●●●

    ●●●●

    ●●●●

    ●●

    ●●●●●●

    ●●●●

    ●●●●●●●●●●●●●●●●●●● ●●● ●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●●●●

    ●●●●●●●●●●●●●●●●●

    ●●●●●●● ●●●●●●●●●●●●●●●●●●●●●●●●●●●● ●●●●●●●

    ●●●●●●●●●●●●●● ● ●● ●●●●

    ●●●●●●●●●●●●● ●●●●●●●

    ●●● ● ●●●●●●●●●●● ●●●●●●●●●●●●●●●●● ●

    ●●

    ●●●●●●

    ●●●●●●●●

    ●●●●

    ●● ●●●●●●●●●●●●●●●● ●●●●●●●●

    ●●●●●●●●●●●●●●●●●●●●

    ●●●●●●●●●●●●●●●●

    ●●

    ●●●

    ●●●●●●●●●●● ●

    ●●●●●●●

    ● ●●●●● ●●●●

    ●●

    ●●

    ●●●

    ● ●● ●●●

    ● ●

    ●●●

    ●● ●●●●●●●●●●● ●

    ●●●●●●

    ●●●

    ●●●

    ●●●●●

    ●●

    ● ●●●●●●●●●●●●●

    ●● ●●●●

    ●●●●●●●●●●●●●●●●●●●●●●●

    ●●

    ●●

    ●●●●

    ●●●

    ●●●●●●●●●●

    ●●●●●●●●●● ●●●●●●●●●●●●●●●●●●●●●●●●

    ●●

    ●●●●●●●●●

    ●●●

    ●●●●●●●●●

    ●●●●●●●●●●●●●●●

    ●●

    ●●●●●● ●●●●●●●●●●●●●

    ● ●●●●●

    ●●●●

    ●●●

    ●●●●●●●●●●●

    ●●●●●●●●●●●●●●●●● ●●●●●●●●●

    ●●

    ●●

    ●●●●

    ●●

    ● ●●●●●●●●●●

    ●●

    ●●

    ●●●●●●●●●●●●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●●●●●●●●

    ●●●● ●●●●●●●●●●●●

    ●●●●●●●●●●●●●●●●

    ●●●●●●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●

    ●● ●

    ●●●

    ●●●●●●●●●●

    ●●●

    ●●

    ●●●●●●●●●

    ●●●●●●●●

    ● ●●●●●●●

    ●●

    ●●

    ●●●

    ●●●

    ●●

    ●●●●●●●●●

    ●●●●

    ●●●●●●●●●●●●●●

    ●●●●●●●

    ●●●●

    ●●●●●●●

    ●●●●●●●●●

    ●●

    ●●●●●●●●

    ●●

    ●●

    ●●●● ●●●●

    ●●●● ●●

    ●●

    ●●

    ●●●●●●●●

    ●●●

    ●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●

    ●●●●●

    ●●●●●●●●●●●●●●●●●●

    ●●●●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●●●●●●●●●

    ●●●●●●●●●

    ●●

    ●●●●●●●●

    ●●●

    ●●

    ●●●●●●●●●●●●●●●●●

    ●●●●●●●●●●●●●●●●●●●●●● ●●●●●●●●●●●●●●●●●●● ●●●●●●●

    ●●●●●●●●

    ●●●●●●

    ●●●●●●●●●●●●●●

    ●●

    ●●

    ●●●●●●●

    ●●●

    ●●

    ●●

    ●●

    ●●●●●●●●●●●●●●●●●●

    ●●●●●●

    ●●● ●

    ●●●

    ● ● ● ●●

    ●●

    ●●●●

    ●●

    ●●●●●●●●●●●●●●●

    ●●●

    ●●●●●●●● ●

    ●●●●●●●●

    ●●

    ●●●●

    ●●●●

    ●●

    ●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●● ●●●●●●●●

    ●●●●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●● ●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●●●

    ●●●

    ●●●●●●●●●●●●

    ●●●●●●●●●●●●●●●●●●●●●● ●●●

    ●●●●●●● ●●●●●●●●●●●●●●●●●●●●●●●●●

    ●● ●●

    ●●●●●●●●●●●●●

    ●●●●●●●●●●●●●

    ●●●●●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●●

    ●●●●

    ●●●●●●●●

    ●●●●●

    ●●●●●●

    ● ●●

    ●● ●● ●●

    ●●

    ●●●

    ●●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    1 2.5 16 50 84 97.5 99

    NOx●

    ●●

    ●●

    ●●●●●

    ●●

    ●●

    ●●

    ●●●

    ●●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ● ●

    ●●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●●●●

    ●●●●

    ●●●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ● ●●●

    ●●

    ●●

    ●●

    ●●●●●

    ●●

    ●●

    ●●●

    ●●

    ●●●●

    ●●

    ●●●●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●●●

    ●●

    ●●●●

    ● ●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●●●

    ●●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●●●●●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●●

    ●●

    ●●●

    ●●●●

    ●●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●●

    ●●●

    ●●

    ●●

    ●●●

    ●● ●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●●

    ●●

    ●●

    ●●

    ● ●

    ●●

    ●●●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●●●

    ●●

    ●●

    ●●

    ●●●●

    ●●●

    ●●

    ●●●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●●●

    ●●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●●

    ●●

    ●●

    ●●●●

    ●●●

    ●●●●

    ●●

    ●●

    ●●

    ●●●

    ●●●●●●

    ●●

    ●●●

    ●●

    ●●

    ●●●

    ●●●●●●

    ●●●●●

    ●●●

    ●●

    ●●

    ●●

    ●●●

    ●●

    ●●●

    ●●●

    ●●

    ●●

    ●●

    ●●●●

    ●●

    ●●

    ●●

    ●●●●●●●●●●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    ●●

    qqobs$x

    exp(

    qqob

    s$y)

    0.1

    1

    10

    ●●●●

    ●●●●

    ●●●●●●

    ●●●●

    ●●●●

    ●●●●●●●

    ●●●●●●

    ●●●●●●●

    ●●●●●●

    ●●

    ●●●●

    ● ●●●●

    ●●●●

    ●●●●●●

    ●●●

    ●●

    ● ●

    ●●

    ●●●

    ●●

    ●●●●●

    ●●●

    ●●●●

    ●●●●

    ●●

    ●●●●

    ●●●

    ●●●●●●

    ●●

    ●●●●●●

    ●●●●

    ●●

    ●●●

    ●●●●●●●●●●●●●

    ●●

    ●●●●●●●●●

    ●●●●

    ●●●●●●●●

    ●●●

    ●●●●●●

    ●●●●●● ●

    ●●●

    ●●●

    ●●

    ●●●●

    ●●●

    ●●●●

    ●●●●

    ●●●●

    ●●●●●●

    ●●●●●●●●

    ●●●●

    ●●

    ●●

    ●●●●

    ●●●●

    ●●

    ●●●●●●

    ●●●

    ●●●

    ●●●●

    ●●

    ●●●

    ●●●

    ●●●

    ●●●●

    ●●

    ●●

    ●●

    ●●

    ●●●

    ●●●●●

    ●●●●●●●●

    ●●●

    ●●●●

    ●●●●

    ●●●●●

    ●●●●

    ●●

    ●●●●

    ●●

    ●●

    ●●

    ●●●●

    ●●●●

    ●●●●

    ●●●

    ●●●●

    ●●●●●●●

    ●●●●

    ●●●●●●

    ●●●●

    ●●●●●●

    ●●

    ●●●●●

    ●●●●●●●

    ●●●●

    ●●●●

    ●●

    ●●●

    ●●●●

    ●●

    ●●

    ●●●●●

    ●●●●

    ●●●●

    ●●

    ●●●

    ●●●●

    ●●●●

    ●●

    ●●●●●●●●●●

    ●●●

    ●●●●●●●

    ●●●

    ●●●●●●●

    ●●●●

    ●●

    ●●●●●●●●●●●●

    ●●●●●●●

    ●●●●●●●●

    ●●●

    ●●●●●

    ●●●

    ●●

    ●●●●

    ●●●●●●

    ●●

    ●●●●●●●●●●

    ●●●

    ●●●●

    ● ●●

    ●●

    ●●

    ●●●●

    ●●●●

    ●●●●

    ●●●

    ●●

    ●●●●

    ●●●

    ●●

    ●●●●●●●

    ●●●●●●●●●●●●●●

    ●●●●●

    ●●●●

    ●●●●●●●

    ●●●●

    ●●●●●●●●●●●●

    ●●●●●●●●●●●●●●

    ●●●●

    ●●

    ●●●●

    ●●●

    ●●●●● ●

    ●●

    ●●●●●●

    ●●●

    ●●

    ●●●

    ●●

    ●●●●●●

    ●●●●●

    ●●●●

    ●●

    ●●●

    ●●●

    ●●●

    ●●●●

    ●●●

    ●●●

    ●●●●

    ●●●●

    ●●●●●●●●

    ●●●●

    ●●●●●●

    ●●●●●●●●●

    ●●●

    ●●●

    ●●●

    ●●

    ●●●●

    ●●●

    ●●●●●●●

    ●●●●●

    ●●

    ●●●●

    ●●

    ●●

    ●●●●●●

    ●●●●

    ●●●

    ●●

    ●●●

    ●●●●

    ●●●

    ●●●

    ●●

    ●●●

    ●●●●●●●●●●

    ●●●

    ●●●●

    ●●●

    ●●●

    ●●●●

    ●●

    ●●●●●●

    ●●●

    ●●●●●●

    ●●

    ●●●

    ●●●●●●●●●●●●●●●

    ●●

    ●●●●●●●●●

    ●●●●

    ●●●

    ●●●●●●

    ●●

    ●●●●●

    ●●●●●●●●

    ●●●●●●●●●●●●

    ●●●●

    ●●●●

    ●●●●

    ●●●●●●●●

    ●●●●

    ●●●●

    ●●

    ●●●

    ●●●●●

    ●●●●●●●

    ●●●

    ●●

    ●●●

    ●●●

    ●●●●●●

    ●●●●●

    ●●●●●●

    ●●●●●

    ●●

    ●●●●

    ●●●●

    ●●●●●

    ●●●

    ●●●●●●

    ●●●

    ●●●●●

    ●●●●●●

    ●●●●

    ●●●●

    ●●●●●

    ●●●

    ●●●●

    ●●●●●●●

    ●●●●

    ●●●●

    ●●

    ●●●●

    ●●

    ●●●●

    ●●●

    ●●●

    ●●●●●

    ●●

    ●●●

    ●●●

    ●●●●●●●●●●

    ●●●

    ●●●●

    ●●

    ●●

    ●●●●

    ●●●

    ●●

    ●●●

    ●●●

    ●●●

    ●●●●●

    ●●

    ●●

    ●●

    ●●

    ●●●●

    ●●●

    ●●●●●●●●●●●

    ●●●●●●

    ●●●●●●●●●●●

    ●●●●

    ●●●●

    ●●●●●

    ●●●●●●

    ●●●●●●●●●●●●●●●●●●●●●●●

    ●●●●

    ● ●

    ●●●●●

    ●●●●●●●

    ●●●●●●●

    ●●●●●●●

    ●●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●●●●●●●●●●●●

    ●●●●●●●●●●●

    ●●

    ●●●●●●●●●●

    ●●●●●●●●●●●●●●●●●●

    ●●●●

    ●●●●●●●●●

    ●●●●●●●●

    ●●●●●●

    ●●●●●●●●●●●●●●●●●●●●

    ●● ●● ● ● ●●

    ●●

    ●●●●●●●●

    ●●

    ●● ●●●●●●●●●●

    ●●

    ●●

    ●●●●●●●●●●●●●●

    ●●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●●●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●●

    ●●●

    ●●●●●●●●●

    ●●

    ●●●●

    ●●●●●●●

    ●●●●

    ●●

    ●●

    ●●●●●●

    ●●

    ●●●●●●●●●●

    ●●●●●●●●●

    ●●●

    ●●●

    ●●●●●●●

    ●●●●●●●●●●●●●●

    ●●●●●●●

    ●●●●●●●●●●

    ●●

    ●●●●●●●●

    ●●

    ●●●●●●

    ●●●●

    ●●

    ●●●●●●●●●●

    ●●●●●●●●●●

    ●●●

    ●●●●●●●●●●●●●●●●●

    ●●●●●●●●●

    ●●

    ●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●●

    ●●●●●●●●

    ●●●●●●●●●●●●

    ●●●●●●●●●�