Modeling the uncertainty of several VOC and its impact on simulated VOC and ozone in Houston, Texas

A WRF-SMOKE-CMAQ modeling system was used to study Volatile Organic Compound (VOC) emissions and their impact on surface VOC and ozone concentrations in southeast Texas during September 2013. The model was evaluated against the ground-level Automated Gas Chromatograph (Auto-GC) measurement data from the Texas Commission on Environmental Quality (TCEQ). The comparisons indicated that the model over-predicted benzene, ethylene, toluene and xylene, while under-predicting isoprene and ethane. The mean biases between simulated and observed values of each VOC species showed clear daytime, nighttime, weekday and weekend variations. Adjusting the VOC emissions using simulated/observed ratios improved model performance of each VOC species, especially mitigating the mean bias substantially. Simulated monthly mean ozone showed a minor change: a 0.4 ppb or 1.2% increase; while a change of more than 5 ppb was seen in hourly ozone data on high ozone days, this change moved model predictions closer to observations. The CMAQ model run with the adjusted emissions better reproduced the variability in the National Aeronautics and Space Administration (NASA)'s Ozone Monitoring Instrument (OMI) formaldehyde (HCHO) columns. The adjusted model scenario also slightly better reproduced the aircraft HCHO concentrations from NASA's DISCOVER-AQ campaign conducted during the simulation episode period; Correlation, Mean Bias and RMSE improved from 0.34, 1.38 ppb and 2.15 ppb to 0.38, 1.33 ppb and 2.08 ppb respectively. A process analysis conducted for both industrial/urban and rural areas suggested that chemistry was the main process contributing to ozone production in both areas, while the impact of chemistry was smaller in rural areas than in industrial and urban areas. For both areas, the positive chemistry contribution increased in the sensitivity simulation largely due to the increase in emissions. Nudging VOC emissions to match the observed concentrations shifted the ozone hotspots outside the industrial/urban region and enhanced the peaked ozone in the outflow region with consistent southerly/southeasterly winds during the afternoon time (1-5 pm). This study helps in the understanding of these processes which are critical to constrain high peak ozone values in the outflow regions. The results indicate that formation of ozone in the outflow could complicate attainment status in neighboring counties.