Receptor models such as the chemical mass balance (CMB) family and factor-analytic methods (PMF and relatives) apportion ambient PM2.5 to emission sources, but different models applied to the same chemical-composition dataset resolve different numbers of sources and different contribution magnitudes. Shi et al. (2018), comparing seven receptor models on Tianjin PM2.5 data, combined them by a simple unweighted average of the available model estimates. Weighted and Bayesian ensembles have been demonstrated at single sites (Atlanta, St. Louis), but the underlying methodology remains unresolved — how to weight member models, how to treat sources that some models fail to resolve, how to propagate uncertainty across correlated models, and how to validate the ensemble when no ground-truth source impacts exist. The open problem is to establish and validate a principled procedure for combining multiple receptor-model source apportionments of the same dataset, replacing the ad hoc unweighted average.
Background
PM2.5 designates airborne particulate matter with aerodynamic diameter below 2.5 micrometers, a dominant air pollutant in Chinese megacities and a leading health concern. Source apportionment (SA) estimates the contribution of each emission source category (e.g., crustal dust, coal combustion, vehicle exhaust, secondary sulfate and nitrate formed in the atmosphere from SO2 and NOx) to measured PM2.5 mass, and its results feed directly into air-quality control policy. Receptor models perform SA at a measurement ("receptor") site by solving a mass-conservation relation between the ambient chemical composition and source profiles, i.e., the relative chemical composition of PM emitted by each source category.
Two families dominate. Chemical mass balance (CMB) is a weighted least-squares fit of ambient species concentrations to known source profiles, applicable when the sources and their profiles are known; factor-analytic models such as positive matrix factorization (PMF) decompose the ambient data matrix into factors that the user interprets as sources without requiring profiles in advance. Variants used in the seven-model comparison include CMB-Iteration (iterative reweighting incorporating source-profile and receptor-data uncertainties); CMB-GC and CMBGC-Iteration (gas-constrained CMB, which adds ambient gas concentrations such as SO2, CO and NOx, through their ratios to PM2.5, as constraints on the solution, with a bound factor n, recommended as three, allowing for gas-phase chemical transformations during transport); WALSPMF (weighted alternating least squares PMF); and NCAPCA (non-negative-constrained absolute principal-component analysis).
A persistent structural difficulty is collinearity: chemically similar sources such as cement dust and crustal dust share marker species (e.g., Ca, Si, Al), so receptor models cannot reliably separate their contributions. Receptor-model SA is also not identifiable in a unique mathematical sense: different models, and even different runs of one model, yield different numbers of resolved sources and different contribution estimates for the same dataset, and no ground-truth time series of source impacts exists at a receptor to judge which model is correct.
Against this background, the idea of combining multiple SA results into an ensemble estimate has emerged. Shi, Liu, Wang, Tian, Wen, Shi, Feng, Ivey and Russell (2018) proposed CMBGC-Iteration, applied it and six other receptor models to 114 PM2.5 samples collected in Tianjin (April 2014 to January 2015), and computed ensemble-average source impacts as the unweighted arithmetic mean of whichever of the seven models provided an estimate for each source category, noting only that this average avoided negative and zero estimates that individual models produced. Earlier, related ensembles had been built with explicit weighting: Lee, Balachandran, Pachon, Shankaran, Lee, Mulholland and Russell (2009) combined four SA methods for health-effect studies using weights inversely proportional to method uncertainty; Balachandran, Pachon, Hu, Lee, Mulholland and Russell (2012) analyzed ensemble method uncertainty and sensitivity to which member methods are included; Balachandran, Chang, Pachon, Holmes, Mulholland and Russell (2013) formulated a Bayesian ensemble averaging of multiple SA methods; and Maier, Balachandran, Sarnat, Turner, Mulholland and Russell (2013) applied an ensemble-trained approach at a site impacted by multiple industrial point sources. None of these established a general, validated procedure for the combination step itself.
Problem Statement
Given several receptor-model source apportionments of the same ambient PM2.5 chemical-composition dataset — e.g., the seven models compared by Shi et al. (2018) on the Tianjin data, which resolve different numbers of source categories (four to six) and different contribution magnitudes — establish and validate a principled procedure for combining them into ensemble source-impact estimates. The procedure must specify: (i) a weighting scheme for member models and the criterion it optimizes, in a setting where member-model estimates are mutually correlated and no ground-truth source impacts exist at the receptor; (ii) how a source category that some member models fail to resolve (e.g., secondary organic carbon, i.e., particulate carbon formed in the atmosphere by oxidation of volatile organic compounds, or cement dust) is treated in the ensemble estimate and its uncertainty; and (iii) how uncertainty is propagated across member models to yield quantitative uncertainty bounds on ensemble source impacts. Validation must demonstrate, against defined performance criteria (e.g., mass closure, i.e., agreement of the sum of apportioned source contributions with the measured PM2.5 mass; agreement with independent constraints such as tracer species or emission-inventory changes; robustness to member inclusion and exclusion; comparison of ensemble uncertainty against individual-model uncertainties), whether and when the ensemble estimate is more accurate or better calibrated than both the individual member models and the unweighted average of available estimates that current practice uses, which the ensemble is meant to replace — with an executed comparison.
References
- Guoliang Shi, Jiayuan Liu, Haiting Wang, Yingze Tian, Jie Wen, Xurong Shi, Yinchang Feng, Cesunica E. Ivey, Armistead G. Russell (2018). Source apportionment for fine particulate matter in a Chinese city using an improved gas-constrained method and comparison with multiple receptor models. Environmental Pollution 233, 1058-1067. DOI: 10.1016/j.envpol.2017.10.007. https://doi.org/10.1016/j.envpol.2017.10.007
- Dongho Lee, Sivaraman Balachandran, Jorge Pachon, Roshini Shankaran, Sangil Lee, James A. Mulholland, Armistead G. Russell (2009). Ensemble-Trained PM2.5 Source Apportionment Approach for Health Studies. Environmental Science & Technology 43(18), 7023-7031. DOI: 10.1021/es9004703. https://doi.org/10.1021/es9004703
- Sivaraman Balachandran, Jorge E. Pachon, Yongtao Hu, Dongho Lee, James A. Mulholland, Armistead G. Russell (2012). Ensemble-trained source apportionment of fine particulate matter and method uncertainty analysis. Atmospheric Environment 61, 387-394. DOI: 10.1016/j.atmosenv.2012.07.031. https://doi.org/10.1016/j.atmosenv.2012.07.031
- Marissa L. Maier, Sivaraman Balachandran, Stefanie E. Sarnat, Jay R. Turner, James A. Mulholland, Armistead G. Russell (2013). Application of an Ensemble-Trained Source Apportionment Approach at a Site Impacted by Multiple Point Sources. Environmental Science & Technology 47(8), 3743-3751. DOI: 10.1021/es304255u. https://doi.org/10.1021/es304255u
- Sivaraman Balachandran, Howard H. Chang, Jorge E. Pachon, Heather A. Holmes, James A. Mulholland, Armistead G. Russell (2013). Bayesian-Based Ensemble Source Apportionment of PM2.5. Environmental Science & Technology 47(23), 13511-13518. DOI: 10.1021/es4020647. https://doi.org/10.1021/es4020647
- Zhenyu Wang, Haofei Yu, Weiqing Liang, Feng Wang, Gen Wang, Da Chen, Weichao Wang, Huan Zhao, Yinchang Feng, Zongbo Shi, Guoliang Shi (2022). Ensemble source apportionment of air pollutants and carbon dioxide based on online measurements. Journal of Cleaner Production 370, 133468. DOI: 10.1016/j.jclepro.2022.133468. https://doi.org/10.1016/j.jclepro.2022.133468
- Weiqing Liang, Zhenyu Wang, Dongheng Zhao, Jinxing Liu, Tong Ma, Yinchang Feng, Guoliang Shi (2025). Ensemble source apportionment of particulate matter and volatile organic compounds and quantifying ensemble source impacts on ozone. Journal of Environmental Sciences 157, 534-545. DOI: 10.1016/j.jes.2024.07.026. https://doi.org/10.1016/j.jes.2024.07.026
- Feng Wang, Haofei Yu, Zhenyu Wang, Weiqing Liang, Guoliang Shi, Jian Gao, Mei Li, Yinchang Feng (2021). Review of online source apportionment research based on observation for ambient particulate matter. Science of The Total Environment 762, 144095. DOI: 10.1016/j.scitotenv.2020.144095. https://doi.org/10.1016/j.scitotenv.2020.144095
- Qili Dai, Philip K. Hopke, Xiaohui Bi, Yinchang Feng (2020). Improving apportionment of PM2.5 using multisite PMF by constraining G-values with a priori information. Science of The Total Environment 736, 139657. DOI: 10.1016/j.scitotenv.2020.139657. https://doi.org/10.1016/j.scitotenv.2020.139657
- Elson Silva Galvão, Rita de Cassia Feroni, Marcos Tadeu D'Azeredo Orlando (2021). A review of the main strategies used in the interpretation of similar chemical profiles yielded by receptor models in the source apportionment of particulate matter. Chemosphere 269, 128746. DOI: 10.1016/j.chemosphere.2020.128746. https://doi.org/10.1016/j.chemosphere.2020.128746
- Yufei Zhang, Xin He, Chenchu Wang, Xuehan Wang, Lilai Song, Zhichao Lu, Xiaohui Bi, Yinchang Feng (2024). Methods and applications for quantitative assessment of uncertainty in atmospheric particulate matter source profiles. Atmospheric Environment 338, 120815. DOI: 10.1016/j.atmosenv.2024.120815. https://doi.org/10.1016/j.atmosenv.2024.120815
- Ran Huang, Zongrun Li, Cesunica E. Ivey, Xinxin Zhai, Guoliang Shi, James A. Mulholland, Robert Devlin, Armistead G. Russell (2022). Application of an improved gas-constrained source apportionment method using data fused fields: A case study in North Carolina, USA. Atmospheric Environment 276, 119031. DOI: 10.1016/j.atmosenv.2022.119031. https://doi.org/10.1016/j.atmosenv.2022.119031