Skip to contents

Calculate Ozone Formation Potential (OFP) of VOC time series. Note: for Chinese VOC name, please also use English punctuation.

Usage

ofp(
  df,
  inunit = "ppbv",
  outunit = "ppbv",
  t = 25,
  p = 101.325,
  stcd = FALSE,
  sortd = TRUE,
  chn = FALSE,
  mtype = "usa",
  bvoc = TRUE
)

Arguments

df

dataframe contains time series.

inunit

input's unit for VOC concentration. A character vector from these options: "ugm" or "ppbv". "ugm" means ug/m3. "ppbv" means part per billion volumn. The default vaule is "ppbv".

outunit

output's unit for VOC concentration. A character from these options: "ugm" or "ppbv". "ugm" means ug/m3. "ppbv" means part per billion volumn. The default vaule is "ppbv".

t

Temperature, in Degrees Celsius, used to convert data in micrograms per cubic meter to standard conditions (25 Degrees Celsius, 101.325 kPa). By default, t equals to 25 Degrees Celsius.

p

Pressure, in kPa, used to convert data in micrograms per cubic meter to standard conditions (25 Degrees Celsius, 101.325 kPa). By default, p equals to 101.325 kPa.

stcd

logical. Does it output results in standard conditions? The default vaule is FALSE.

sortd

logical value. It determines whether the VOC species are sorted or not. By default, sortd has value "TRUE". If TRUE, VOC species in time series will be arranged according to VOC group, relative molecular weight, and MIR value.

chn

logical. Dose colnames present as Chinese? The default vaule is FALSE.

mtype

text. "usa" for MIR value from USA, "chn" for MIR value from CHINA.

bvoc

logical. Whether you want to list BVOC as a separate VOC group? The default vaule is TRUE.

Value

a list contains 5 tables: MIR_Result: the matched MIR value results; OFP_Result: the OFP time series of VOC by species; OFP_Result_stat: the statistics of OFP of VOC by species; OFP_Result_group: OFP the time series of VOC classified by groups; OFP_Result_group_stat: the statistics of OFP of VOC according to major groups.

Details

The CAS number is matched for each VOC speices (from column name), and the Maximum Incremental Reactivity (MIR) value is matched through the CAS number and used for time series calculation.
The MIR value comes from <https://ww2.arb.ca.gov/sites/default/files/classic/regact/2009/mir2009/mir10.pdf>, Zhang et al.(2021) <doi:10.5194/acp-21-11053-2021>.

Examples

voc_ofp=ofp(voc)
summary(voc_ofp)
#>                       Length Class      Mode
#> MIR_Result            8      data.frame list
#> OFP_Result            6      data.frame list
#> OFP_Result_stat       9      data.frame list
#> OFP_Result_group      4      data.frame list
#> OFP_Result_group_stat 9      data.frame list