TY - JOUR
T1 - A strategy for modelling heavy-tailed greenhouse gases (GHG) data using the generalised extreme value distribution
T2 - Are we overestimating GHG flux using the sample mean?
AU - Dhanoa, M. S.
AU - Louro, A.
AU - Cardenas, L. M.
AU - Shepherd, A.
AU - Sanderson, R.
AU - López, S.
AU - France, J.
PY - 2020/9/15
Y1 - 2020/9/15
N2 - In this study, we draw up a strategy for analysis of greenhouse gas (GHG) field data. The distribution of GHG flux data generally exhibits excessive skewness and kurtosis. This results in a heavy tailed distribution that is much longer than the tail of a log-normal distribution or outlier induced skewness. The generalised extreme value (GEV) distribution is well-suited to model such data. We evaluated GEV as a model for the analysis and a means of extraction of a robust average of carbon dioxide (CO2) and nitrous oxide (N2O) flux data measured in an agricultural field. The option of transforming CO2 flux data to the Box-Cox scale in order to make the distribution normal was also investigated. The results showed that average CO2 estimates from GEV are less affected by data in the long tail compared to the sample mean. The data for N2O flux were much more complex than CO2 flux data due to the presence of negative fluxes. The estimate of the average value from GEV was much more consistent with maximum data frequency position. The analysis of GEV, which considers the effects of hot-spot-like observations, suggests that sample means and log-means may overestimate GHG fluxes from agricultural fields. In this study, the arithmetic CO2 sample mean of 65.6 (mean log-scale 65.9) kg CO2–C ha−1 d−1 was reduced to GEV mean of 60.1 kg CO2–C ha−1 d−1. The arithmetic N2O sample mean of 1.038 (mean log-scale 1.038) kg N2O–N ha−1 d−1 was substantially reduced to GEV mean of 0.0157 kg N2O–N ha−1 d−1. Our analysis suggests that GHG data should be analysed assuming a GEV distribution of the data, including a Box-Cox transformation when negative data are observed, rather than only calculating basic log and log-normal summaries. Results of GHG studies may end up in national inventories. Thus, it is necessary and important to follow all procedures that contribute to minimise any bias in the data.
AB - In this study, we draw up a strategy for analysis of greenhouse gas (GHG) field data. The distribution of GHG flux data generally exhibits excessive skewness and kurtosis. This results in a heavy tailed distribution that is much longer than the tail of a log-normal distribution or outlier induced skewness. The generalised extreme value (GEV) distribution is well-suited to model such data. We evaluated GEV as a model for the analysis and a means of extraction of a robust average of carbon dioxide (CO2) and nitrous oxide (N2O) flux data measured in an agricultural field. The option of transforming CO2 flux data to the Box-Cox scale in order to make the distribution normal was also investigated. The results showed that average CO2 estimates from GEV are less affected by data in the long tail compared to the sample mean. The data for N2O flux were much more complex than CO2 flux data due to the presence of negative fluxes. The estimate of the average value from GEV was much more consistent with maximum data frequency position. The analysis of GEV, which considers the effects of hot-spot-like observations, suggests that sample means and log-means may overestimate GHG fluxes from agricultural fields. In this study, the arithmetic CO2 sample mean of 65.6 (mean log-scale 65.9) kg CO2–C ha−1 d−1 was reduced to GEV mean of 60.1 kg CO2–C ha−1 d−1. The arithmetic N2O sample mean of 1.038 (mean log-scale 1.038) kg N2O–N ha−1 d−1 was substantially reduced to GEV mean of 0.0157 kg N2O–N ha−1 d−1. Our analysis suggests that GHG data should be analysed assuming a GEV distribution of the data, including a Box-Cox transformation when negative data are observed, rather than only calculating basic log and log-normal summaries. Results of GHG studies may end up in national inventories. Thus, it is necessary and important to follow all procedures that contribute to minimise any bias in the data.
KW - Carbon dioxide
KW - Finney correction
KW - Generalised extreme value
KW - Heavy-tailed data
KW - Nitrous oxide
KW - Skewness correction
UR - http://www.scopus.com/inward/record.url?scp=85087308864&partnerID=8YFLogxK
U2 - 10.1016/j.atmosenv.2020.117500
DO - 10.1016/j.atmosenv.2020.117500
M3 - Article
AN - SCOPUS:85087308864
SN - 1352-2310
VL - 237
JO - Atmospheric Environment
JF - Atmospheric Environment
M1 - 117500
ER -