High simulated surface ozone possibly caused by initial and boundary conditions

[Jaelyn Liu]

Hi,

 

I ran the WRF-Chem model using Chem_opt=202. And I used mozbc and CAM-Chem output data to generate initial and boundary conditions.

My WRF-Chem model simulated surface ozone much higher than observations. The averaged MDA8 ozone (both simulated and observations) over 3/20/2019-3/30/2019 are attached. The model spinup period is 3/1/2019-3/19/2019.

And I also found that the ozone concentration in the initial wrfinput file is high. I think the high simulated ozone may be affected by the initial (and/or boundary) conditions. My mozbc input is as follows:

Could you help check my mozbc process? Or could you help identify other factors that may cause high simulated ozone? My namelist input is also attached.

 

Thanks,

 

&control

do_bc           = .true.

do_ic           = .true.

domain          = 2

dir_wrf         = '//'

dir_moz         = '/mozbc/data/'

fn_moz          = 'ha0001.nc'

moz_var_suffix  = ''

def_missing_var = .true.

spc_map = 'o3 -> O3', 'n2o -> N2O', 'no -> NO',

          'no2 -> NO2', 'no3 -> NO3', 'nh3 -> NH3', 'hno3 -> HNO3', 'hno4 -> HO2NO2',

          'n2o5 -> N2O5','ho -> OH', 'ho2 -> HO2', 'h2o2 -> H2O2',

          'ch4 -> CH4', 'co -> CO', 'ch3ooh -> CH3OOH',

          'hcho -> CH2O', 'ch3oh -> CH3OH', 'c2h4 -> C2H4',

          'ald -> CH3CHO', 'ch3cooh -> CH3COOH', 'acet -> CH3COCH3', 'mgly -> CH3COCHO',

          'hcooh -> HCOOH', 'c2h2 -> C2H2',

          'pan -> PAN', 'mpan -> MPAN', 'macr -> MACR',

          'mvk -> MVK', 'c2h6 -> C2H6', 'c3h6 -> C3H6', 'c3h8 -> C3H8',

          'c2h5oh -> C2H5OH', 'apin -> 0.333*MTERP', 'bpin -> 0.333*MTERP', 'limon -> 0.333*MTERP',

          'onit -> ALKNIT+HONITR+NOA', 'onitr -> ONITR', 'isopr -> ISOP', 'bcary -> BCARY',

          'isopn -> ISOPNITA+ISOPNITB', 'acetol -> HYAC',

          'mek -> MEK',

          'bigene -> BIGENE', 'open -> BIGALD', 'bigalk -> BIGALK',

          'bigald1 -> BIGALD1', 'bigald2 -> BIGALD2', 'bigald3 -> BIGALD3', 'bigald4 -> BIGALD4',

          'gly -> GLYOXAL', 'ch4 -> CH4', 'pbznit -> PBZNIT',

          'tol -> TOLUENE', 'xyl -> XYLENES', 'benzene -> BENZENE', 'phen -> PHENOL',

          'cres -> CRESOL', 'dms -> DMS', 'so2 -> SO2',

          'oc_a01->0.1216*pom_a1+0.1216*pom_a4;1.e9',

          'oc_a02->0.7618*pom_a1+0.7618*pom_a4;1.e9',

          'oc_a03->0.1164*pom_a1+0.1164*pom_a4;1.e9',

          'oc_a04->0.0002*pom_a1+0.0002*pom_a4;1.e9',

          'bsoa1_a01->0.9886*soa1_a2+0.1216*soa1_a1;1.e9'

          'bsoa1_a02->0.0114*soa1_a2+0.7618*soa1_a1;1.e9',

          'bsoa1_a03->0.1164*soa1_a1;1.e9',

          'bsoa1_a04->0.0002*soa1_a1;1.e9',

          'bsoa2_a01->0.9886*soa2_a2+0.1216*soa2_a1;1.e9'

          'bsoa2_a02->0.0114*soa2_a2+0.7618*soa2_a1;1.e9',

          'bsoa2_a03->0.1164*soa2_a1;1.e9',

          'bsoa2_a04->0.0002*soa2_a1;1.e9',

          'bsoa3_a01->0.9886*soa3_a2+0.1216*soa3_a1;1.e9'

          'bsoa3_a02->0.0114*soa3_a2+0.7618*soa3_a1;1.e9',

          'bsoa3_a03->0.1164*soa3_a1;1.e9',

          'bsoa3_a04->0.0002*soa3_a1;1.e9',

          'bsoa4_a01->0.9886*soa4_a2+0.1216*soa4_a1;1.e9'

          'bsoa4_a02->0.0114*soa4_a2+0.7618*soa4_a1;1.e9',

          'bsoa4_a03->0.1164*soa4_a1;1.e9',

          'bsoa4_a04->0.0002*soa4_a1;1.e9',

          'bsoaX_a01->0.9886*soa5_a2+0.1216*soa5_a1;1.e9'

          'bsoaX_a02->0.0114*soa5_a2+0.7618*soa5_a1;1.e9',

          'bsoaX_a03->0.1164*soa5_a1;1.e9',

          'bsoaX_a04->0.0002*soa5_a1;1.e9',

          'bc_a01->0.1216*bc_a1+0.1216*bc_a4;1e9',

          'bc_a02->0.7618*bc_a1+0.7618*bc_a4;1e9',

          'bc_a03->0.1164*bc_a1+0.1164*bc_a4;1e9',

          'bc_a04->0.0002*bc_a1+0.0002*bc_a4;1e9',

          'so4_a01->0.9886*so4_a2+0.1216*so4_a1;1e9',

          'so4_a02->0.0114*so4_a2+0.7618*so4_a1+0.0002*so4_a3;1e9',

          'so4_a03->0.1164*so4_a1+0.0995*so4_a3;1e9',

          'so4_a04->0.0002*so4_a1+0.9003*so4_a3;1e9',

          'nh4_a01->0.1856*so4_a2+0.0050*so4_a1;1e9',

          'nh4_a02->0.0021*so4_a2+0.0930*so4_a1;1e9',

          'nh4_a03->0.0203*so4_a1+0.0186*so4_a3;1e9',

          'nh4_a04->0.1690*so4_a3;1e9',

          'na_a01->0.3889*ncl_a2+0.0479*ncl_a1;1.e9',

          'na_a02->0.0045*ncl_a2+0.2997*ncl_a1;1.e9',

          'na_a03->0.0458*ncl_a1+0.0391*ncl_a3;1.e9',

          'na_a04->0.3542*ncl_a3;1.e9',

          'cl_a01->0.5996*ncl_a2+0.0737*ncl_a1;1.e9',

          'cl_a02->0.0068*ncl_a2+0.4621*ncl_a1;1.e9',

          'cl_a03->0.0709*ncl_a1+0.0604*ncl_a3;1.e9',

          'cl_a04->0.0001*ncl_a1+0.5462*ncl_a3;1.e9',

          'oin_a01->0.9886*dst_a2+0.1216*dst_a1;1.e9',

          'oin_a02->0.0114*dst_a2+0.7618*dst_a1+0.0002*dst_a3;1.e9',

          'oin_a03->0.1164*dst_a1+0.0995*dst_a3;1.e9',

          'oin_a04->0.0002*dst_a1+0.9003*dst_a3;1.e9',

          'num_a01->0.9996*num_a2+0.7135*num_a1;1.0',

          'num_a02->0.0004*num_a2+0.2847*num_a1+0.0239*num_a3;1.0',

          'num_a03->0.0016*num_a1+0.6258*num_a3;1.0',

          'num_a04->0.3501*num_a3;1.0',

/