Modelling consortium for chemistry of indoor environments (MOCCIE): integrating chemical processes from molecular to room scales

Composition of indoor organic surface films in residences: simulating the influence of sources, partitioning, particle deposition, and air exchange

Indoor surfaces are coated with organic films that modulate thermodynamic interactions between the surfaces and room air. Recently published models can simulate film formation and growth via gas-surface partitioning, but none have statistically investigated film composition. The Indoor Model of Aerosols, Gases, Emissions, and Surfaces (IMAGES) was used here to simulate ten years of nonreactive film growth upon impervious indoor surfaces within a Monte Carlo procedure representing a sub-set of North American residential buildings. Film composition was resolved into categories reflecting indoor aerosol (gas + particle phases) factors from three sources: outdoor-originating, indoor-emitted, and indoor-generated secondary organic material. In addition to gas-to-film partitioning, particle deposition was modeled as a vector for organics to enter films, and it was responsible for a majority of the film mass after ∼1000 days of growth for the median simulation and is likely the main source of LVOCs within films. Therefore, the organic aerosol factor possessing the most SVOCs contributes most strongly to the composition of early films, but as the film ages, films become more dominated by the factor with the highest particle concentration. Indoor-emitted organics (e.g. from cooking) often constituted at least a plurality of the simulated mass in developed films, but indoor environments are diverse enough that any major organic material source could be the majority contributor to film mass, depending on building characteristics and indoor activities. A sensitivity analysis suggests that rapid film growth is most likely in both newer, more air-tight homes and older homes near primary pollution sources.

Source apportionment for indoor air pollution: Current challenges and future directions

Source apportionment (SA) for indoor air pollution is challenging due to the multiplicity and high variability of indoor sources, the complex physical and chemical processes that act as primary sources, sinks and sources of precursors that lead to secondary formation, and the interconnection with the outdoor environment. While the major indoor sources have been recognized, there is still a need for understanding the contribution of indoor versus outdoor-generated pollutants penetrating indoors, and how SA is influenced by the complex processes that occur in indoor environments. This paper reviews our current understanding of SA, through reviewing information on the SA techniques used, the targeted pollutants that have been studied to date, and their source apportionment, along with limitations or knowledge gaps in this research field. The majority (78 %) of SA studies to date focused on PM chemical composition/size distribution, with fewer studies covering organic compounds such as ketones, carbonyls and aldehydes. Regarding the SA method used, the majority of studies have used Positive Matrix Factorization (31 %), Principal Component Analysis (26 %) and Chemical Mass Balance (7 %) receptor models. The indoor PM sources identified to date include building materials and furniture emissions, indoor combustion-related sources, cooking-related sources, resuspension, cleaning and consumer products emissions, secondary-generated pollutants indoors and other products and activity-related emissions. The outdoor environment contribution to the measured pollutant indoors varies considerably (<10 %- 90 %) among the studies. Future challenges for this research area include the need for optimization of indoor air quality monitoring and data selection as well as the incorporation of physical and chemical processes in indoor air into source apportionment methodology.

Model Framework for Predicting Semivolatile Organic Material Emissions Indoors from Organic Aerosol Measurements: Applications to HOMEChem Stir-Frying

Cooking activities emit myriad low-volatility, semivolatile, and highly volatile organic compounds that together form particles that can accumulate to large indoor concentrations. Absorptive partitioning thermodynamics governs the particle-phase organic aerosol concentration mainly via temperature and sorbing mass impacts. Cooking activities can increase the organic sorbing mass by 1-2 orders of magnitude, increasing particle-phase concentrations and affecting emission rate calculations. Although recent studies have begun to probe the volatility characteristics of indoor cooking particles, parametrizations of cooking particle mass emissions have largely neglected these thermodynamic considerations. Here, we present an improved thermodynamics-based model framework for estimating condensable organic material emission rates from a time series of observed concentrations, given that adequate measurements or assumptions can be made about the volatility of the emitted species. We demonstrate the performance of this methodology by applying data from stir-frying experiments performed during the House Observations of Microbial and Environmental Chemistry (HOMEChem) campaign to a two-zone box model representing the UTest House. Preliminary estimates of organic mass emitted on a per-stir-fry basis for three types of organic aerosol factors are presented. Our analysis highlights that using traditional nonvolatile particle models and emission characterizations for some organic aerosol emitting activities can incorrectly attribute concentration changes to emissions rather than thermodynamic effects.

Phase state of organic aerosols may limit temperature-driven thermodynamic repartitioning following outdoor-to-indoor transport

Ambient aerosols often experience temperature and humidity gradients following outdoor-to-indoor transport, causing organic aerosols (OA) to either gain or lose mass via gas-particle repartitioning. Recent models have sought to quantify these effects using equilibrium partitioning thermodynamics. However, evidence suggests some indoor OA may possess glassy or semisolid phase states with higher viscosities than liquid OA. Characteristic partitioning timescales of higher-viscosity particles are significantly longer than for liquid particles, which may either fully or partially inhibit repartitioning. For outdoor OA experiencing a temperature change during transport indoors, the ultimate repartitioning state depends on the relationship between the gas-particle partitioning rate coefficient (k_gp) of semivolatile organics and the indoor particle loss rate coefficient (l_p). That is, thermodynamic equilibrium partitioning may occur when semivolatile k_gp ≫ l_p, no repartitioning when semivolatile k_gp ≪ l_p, and partial repartitioning when their magnitudes are similar. Longer indoor particle lifetimes, higher particle number, and larger particle sizes all raise k_gp (driving repartitioning towards equilibrium). For simulated U.S. residences, equilibrium condensation was likely reached in humid climate zones during warm meteorological conditions. In colder regions, the degree of evaporative repartitioning depended on whether organics could repartition before the particle phase state adjusts to indoor conditions, which is uncertain. When an appreciable temperature gradient exists, this study not only confirmed that all outdoor-originating OA that is liquid indoors will reach thermodynamic equilibrium, but also concluded that a plurality (46% for this domain) of such OA that is semisolid may also achieve thermodynamic equilibrium during its indoor lifetime.

Human skin responses to environmental pollutants: A review of current scientific models

Whatever the exposure route, chemical, physical and biological pollutants modify the whole organism response, leading to nerve, cardiac, respiratory, reproductive, and skin system pathologies. Skin acts as a barrier for preventing pollutant modifications. This review aims to present the available scientific models, which help investigate the impact of pollution on the skin. The research question was "Which experimental models illustrate the impact of pollution on the skin in humans?" The review covered a period of 10 years following a PECO statement on in vitro, ex vivo, in vivo and in silico models. Of 582 retrieved articles, 118 articles were eligible. In oral and inhalation routes, dermal exposure had an important impact at both local and systemic levels. Healthy skin models included primary cells, cell lines, co-cultures, reconstructed human epidermis, and skin explants. In silico models estimated skin exposure and permeability. All pollutants affected the skin by altering elasticity, thickness, the structure of epidermal barrier strength, and dermal extracellular integrity. Some specific models concerned wound healing or the skin aging process. Underlying mechanisms were an exacerbated inflammatory skin reaction with the modulation of several cytokines and oxidative stress responses, ending with apoptosis. Pathological skin models revealed the consequences of environmental pollutants on psoriasis, atopic dermatitis, and tumour development. Finally, scientific models were used for evaluating the safety and efficacy of potential skin formulations in preventing the skin aging process or skin irritation after repeated contact. The review gives an overview of scientific skin models used to assess the effects of pollutants. Chemical and physical pollutants were mainly represented while biological contaminants were little studied. In future developments, cell hypoxia and microbiota models may be considered as more representative of clinical situations. Models considering humidity and temperature variations may reflect the impact of these changes.

Simulating indoor inorganic aerosols of outdoor origin with the inorganic aerosol thermodynamic equilibrium model ISORROPIA

Outdoor aerosols can transform and have their composition altered upon transport indoors. Herein, IMAGES, a platform that simulates indoor organic aerosol with the 2-dimensional volatility basis set (2D-VBS), was extended to incorporate the inorganic aerosol thermodynamic equilibrium model, ISORROPIA. The model performance was evaluated by comparing aerosol component predictions to indoor measurements from an aerosol mass spectrometer taken during the summer and winter seasons. Since ammonia was not measured in the validation dataset, outdoor ammonia was estimated from aerosol measurements using a novel pH-based algorithm, while nitric acid was held constant. Modeled indoor ammonia sources included temperature-based occupant and surface emissions. Sensitivity to the nitric acid indoor surface deposition rate β g , HNO 3 , g was explored by varying it in model runs, which did not affect modeled sulfate due to its non-volatile nature, though the fitting of a filter efficiency was required for good correlations of modeled sulfate with measurements in both seasons. Modeled summertime nitrate well-matched measured observations when β g , HNO 3 , g = 2.75 h - 1 , but wintertime comparisons were poor, possibly due to missing thermodynamic processes within the heating, ventilating, and air-conditioning (HVAC) system. Ammonium was consistently overpredicted, potentially due to neglecting thirdhand smoke impacts observed in the field campaign, as well as HVAC impacts.

Understanding Ozone Transport and Deposition within Indoor Surface Boundary Layers

Ozone-initiated oxidation reactions on indoor surfaces meaningfully alter the chemical composition of indoor air and human exposure to air toxins. Ozone mass transport within the indoor surface boundary layer plays a key role in ozone-surface reaction kinetics. However, limited information is available on detailed ozone transport dynamics near realistic, irregular indoor surfaces. This paper presents a research framework to study the underlying mechanisms of ozone reactions with realistic indoor surfaces based on microscope scanning of surface material and detailed Computational Fluid Dynamics (CFD) simulation. The study results show that indoor surface topography can meaningfully affect ozone mass transport within a surface boundary layer, thereby modulating near-surface ozone concentration gradient and surface uptake. The results also reveal that the effective indoor surface area available for ozone reaction varies with indoor air speed and turbulent air mixing within the boundary layer. The detailed dynamic behaviors of ozone reactions with realistic indoor surfaces provide insights into the implications of pollutant-surface interactions on indoor chemistry and air quality.

Reactions and Products of Squalene and Ozone: A Review

This critical review describes the squalene-ozone (SqOz) reaction, or squalene ozonolysis. Ambient ozone penetrates indoors and drives indoor air chemistry. Squalene, a component of human skin oil, contains six carbon-carbon double bonds and is very reactive with ozone. Bioeffluents from people contribute to indoor air chemistry and affect the indoor air quality, resulting in exposures because people spend the majority of their time indoors. The SqOz reaction proceeds through various formation pathways and produces compounds that include aldehydes, ketones, carboxylic acids, and dicarbonyl species, which have a range of volatilities. In this critical review of SqOz chemistry, information on the mechanism of reaction, reaction probability, rate constants, and reaction kinetics are compiled. Characterizations of SqOz reaction products have been done in laboratory experiments and real-world settings. The effect of multiple environmental parameters (ozone concentration, air exchange rate (AER), temperature, and relative humidity (RH)) in indoor settings are summarized. This critical review concludes by identifying the paucity of available exposure, health, and toxicological data for known reaction products. Key knowledge gaps about SqOz reactions leading to indoor exposures and adverse health outcomes are provided as well as an outlook on where the field is headed.

Multiphase Ozonolysis of Oleic Acid-Based Lipids: Quantitation of Major Products and Kinetic Multilayer Modeling

Commonly found in atmospheric aerosols, cooking oils, and human sebum, unsaturated lipids rapidly decay upon exposure to ozone, following the Criegee mechanism. Here, the gas-surface ozonolysis of three oleic acid-based compounds was studied in a reactor and indoors. Under dry conditions, quantitative product analyses by ¹H NMR indicate up to 79% molar yield of stable secondary ozonides (SOZs) in oxidized triolein and methyl oleate coatings. Elevated relative humidity (RH) significantly suppresses the SOZ yields, enhancing the formation of condensed-phase aldehydes and volatile C9 products. Along with kinetic parameters informed by molecular dynamics simulations, these results were used as constraints in a kinetic multilayer model (KM-GAP) simulating triolein ozonolysis. Covering a wide range of coating thicknesses and ozone levels, the model predicts a much faster decay near the gas-lipid interface compared to the bulk. Although the dependence of RH on SOZ yields is well predicted, the model overestimates the production of H₂O₂ and aldehydes. With negligible dependence on RH, the product composition for oxidized oleic acid is substantially affected by a competitive reaction between Criegee intermediates (CIs) and carboxylic acids. The resulting α-acyloxyalkyl hydroperoxides (α-AAHPs) have much higher molar yields (29-38%) than SOZs (12-16%). Overall, the ozone-lipid chemistry could affect the indoor environment through "crust" accumulation on surfaces and volatile organic compound (VOC) emission. In the atmosphere, the peroxide formation and changes in particle hygroscopicity may have effects on climate. The related health impacts are also discussed.

A holistic modeling framework for estimating the influence of climate change on indoor air quality

The IPCC 2021 report predicts rising global temperatures and more frequent extreme weather events in the future, which will have different effects on the regional climate and concentrations of ambient air pollutants. Consequently, changes in heat and mass transfer between the inside and outside of buildings will also have an increasing impact on indoor air quality. It is therefore surprising that indoor spaces and occupant well-being still play a subordinate role in the studies of climate change. To increase awareness for this topic, the Indoor Air Quality Climate Change (IAQCC) model system was developed, which allows short and long-term predictions of the indoor climate with respect to outdoor conditions. The IAQCC is a holistic model that combines different scenarios in the form of submodels: building physics, indoor emissions, chemical-physical reaction and transformation, mold growth, and indoor exposure. IAQCC allows simulation of indoor gas and particle concentrations with outdoor influences, indoor materials and activity emissions, particle deposition and coagulation, gas reactions, and SVOC partitioning. These key processes are fundamentally linked to temperature and relative humidity. With the aid of the building physics model, the indoor temperature and humidity, and pollutant transport in building zones can be simulated. The exposure model refers to the calculated concentrations and provides evaluations of indoor thermal comfort and exposure to gaseous, particulate, and microbial pollutants.

Using Real Time Measurements to Derive the Indoor and Outdoor Contributions of Submicron Particulate Species and Trace Gases

The indoor environment is usually more polluted than outdoors due to emissions of gas and particle-phase pollutants from multiple sources, leading to their accumulation on top of the infiltration of outdoor pollution. While it is widely recognized that negative health effects arise from the exposure to outdoor air pollution, exposure to indoor pollutants also needs to be well assessed since we spend most of our time (~90%) breathing indoors. Indoor concentrations of pollutants are driven by physicochemical processes and chemical transformations taking place indoors, acting as sources and/or sinks. While these basic concepts are understood, assessing the contribution of each process is still challenging. In this study, we deployed online instrumentation in an unoccupied room to test a methodology for the apportionment of indoor and outdoor pollutant sources. This method was successfully applied to the apportionment of PM₁ and VOCs, however, there are limitations for reactive gases such as O₃. The results showed that this unoccupied indoor environment acts as a source of VOCs and contributes 87% on OVOCs and 6% on C_xH_y, while it acts as a sink for particles, likely due to losses through volatilization up to 60%.

Kinetic multi-layer model of film formation, growth, and chemistry (KM-FILM): Boundary layer processes, multi-layer adsorption, bulk diffusion, and heterogeneous reactions

Large surface area-to-volume ratios indoors cause heterogeneous interactions to be especially important. Semi-volatile organic compounds can deposit on impermeable indoor surfaces forming thin organic films. We developed a new model to simulate the initial film formation by treating gas-phase diffusion and turbulence through a surface boundary layer and multi-layer reversible adsorption on rough surfaces, as well as subsequent film growth by resolving bulk diffusion and chemical reactions in a film. The model was applied with consistent parameters to reproduce twenty-one sets of film formation measurements due to multi-layer adsorption of multiple phthalates onto different indoor-relevant surfaces, showing that the films should initially be patchy with the formation of pyramid-like structures on the surface. Sensitivity tests showed that highly turbulent conditions can lead to the film growing by more than a factor of two compared to low turbulence conditions. If surface films adopt an ultra-viscous state with bulk diffusion coefficients of less than 10^-18  cm² s^-1 , a significant decrease in film growth is expected. The presence of chemical reactions in the film has the potential to increase the rate of film growth by nearly a factor of two.

Spatial and temporal scales of variability for indoor air constituents

Historically air constituents have been assumed to be well mixed in indoor environments, with single point measurements and box modeling representing a room or a house. Here we demonstrate that this fundamental assumption needs to be revisited through advanced model simulations and extensive measurements of bleach cleaning. We show that inorganic chlorinated products, such as hypochlorous acid and chloramines generated via multiphase reactions, exhibit spatial and vertical concentration gradients in a room, with short-lived ⋅OH radicals confined to sunlit zones, close to windows. Spatial and temporal scales of indoor constituents are modulated by rates of chemical reactions, surface interactions and building ventilation, providing critical insights for better assessments of human exposure to hazardous pollutants, as well as the transport of indoor chemicals outdoors.

Improving Predictions of Indoor Aerosol Concentrations of Outdoor Origin by Considering the Phase Change of Semivolatile Material Driven by Temperature and Mass-Loading Gradients

Outdoor aerosols experience environmental changes as they are transported indoors, including outdoor-to-indoor temperature and mass-loading gradients, which can reduce or enhance their indoor concentrations due to repartitioning driven by changes in thermodynamic equilibrium states. However, the complexity required to model repartitioning typically hinders its inclusion in studies predicting indoor exposure to ambient aerosols. To facilitate exposure predictions, this work used an explicit thermodynamic indoor aerosol model to simulate outdoor-to-indoor aerosol repartitioning typical for residential and office buildings across the 16 U.S. climate zones over an annual time horizon. Results demonstrate that neglecting repartitioning when predicting indoor concentrations can produce errors of up to 80-100% for hydrocarbon-like organic aerosol, 40-60% for total organic aerosol, 400% for ammonium nitrate, and 60% (typically 3 μg/m³) for the total PM_2.5 aerosol. Underpredictions were more likely for buildings in hotter than colder regions, and for residences than offices, since both cooler indoor air and more meaningful residential organic aerosol concentrations encourage condensation of semivolatile organics. Furthermore, a method for computing correction factors to more easily account for thermodynamic repartitioning is provided. Applying these correction factors to mechanical-only aerosol predictions significantly reduced errors to <0.5 μg/m³ for the total indoor PM_2.5 while bypassing explicit thermodynamic simulations.

Yields and Variability of Ozone Reaction Products from Human Skin

The skin of 20 human participants was exposed to ∼110 ppb O₃ and volatile products of the resulting chemistry were quantified in real time. Yields (ppb product emitted/ppb ozone consumed) for 40 products were quantified. Major products of the primary reaction of ozone-squalene included 6-methyl 5-hepten-2-one (6-MHO) and geranyl acetone (GA) with average yields of 0.22 and 0.16, respectively. Other major products included decanal, methacrolein (or methyl vinyl ketone), nonanal, and butanal. Yields varied widely among participants; summed yields ranged from 0.33 to 0.93. The dynamic increase in emission rates during ozone exposure also varied among participants, possibly indicative of differences in the thickness of the skin lipid layer. Factor analysis indicates that much of the variability among participants is due to factors associated with the relative abundance of (1) "fresh" skin lipid constituents (such as squalene and fatty acids), (2) oxidized skin lipids, and (3) exogenous compounds. This last factor appears to be associated with the presence of oleic and linoleic acids and could be accounted for by uptake of cooking oils or personal care products to skin lipids.

Assessing Human Exposure to SVOCs in Materials, Products, and Articles: A Modular Mechanistic Framework

A critical review of the current state of knowledge of chemical emissions from indoor sources, partitioning among indoor compartments, and the ensuing indoor exposure leads to a proposal for a modular mechanistic framework for predicting human exposure to semivolatile organic compounds (SVOCs). Mechanistically consistent source emission categories include solid, soft, frequent contact, applied, sprayed, and high temperature sources. Environmental compartments are the gas phase, airborne particles, settled dust, indoor surfaces, and clothing. Identified research needs are the development of dynamic emission models for several of the source emission categories and of estimation strategies for critical model parameters. The modular structure of the framework facilitates subsequent inclusion of new knowledge, other chemical classes of indoor pollutants, and additional mechanistic processes relevant to human exposure indoors. The framework may serve as the foundation for developing an open-source community model to better support collaborative research and improve access for application by stakeholders. Combining exposure estimates derived using this framework with toxicity data for different end points and toxicokinetic mechanisms will accelerate chemical risk prioritization, advance effective chemical management decisions, and protect public health.

Indoor Surface Chemistry: Developing a Molecular Picture of Reactions on Indoor Interfaces

Chemical reactions on indoor surfaces play an important role in air quality in indoor environments, where humans spend 90% of their time. We focus on the challenges of understanding the complex chemistry that takes place on indoor surfaces and identify crucial steps necessary to gain a molecular-level understanding of environmental indoor surface chemistry: (1) elucidate key surface reaction mechanisms and kinetics important to indoor air chemistry, (2) define a range of relevant and representative surfaces to probe, and (3) define the drivers of surface reactivity, particularly with respect to the surface composition, light, and temperature. Within the drivers of surface composition are the roles of adsorbed/absorbed water associated with indoor surfaces and the prevalence, inhomogeneity, and properties of secondary organic films that can impact surface reactivity. By combining laboratory studies, field measurements, and modeling we can gain insights into the molecular processes necessary to further our understanding of the indoor environment.

Spatial distributions of ozonolysis products from human surfaces in ventilated rooms

Ozone has adverse effects on human health. Skin oil on the human surface acts as an ozone sink indoors, producing oxidation products that can cause skin and respiratory irritations. Concentrations of ozone and oxidation products near human surfaces, including the breathing zone, can be modulated by indoor ventilation modes and human surface conditions. The objective of this study is to examine concentrations and spatial heterogeneity of ozone and ozonolysis products under representative ranges of indoor ventilation, clothing, and breathing conditions. Using computational fluid dynamics (CFD) simulation in conjunction with a chemical kinetic model, details of ozone reactions with the human surface and subsequent chemical reactions are examined. The results show that primary ozonolysis products are concentrated near the soiled clothing, while the secondary products are relatively well distributed throughout the room. Increasing indoor air mixing enhances the ozone deposition to the human surface, thereby resulting in higher emission rates of oxidation products in the room. Soiled clothing consumes more ozone than clean clothing and accordingly produces ~ 65% more primary products and ~15% more secondary products. The results also reveal that unsaturated hydrocarbons from the human breath, such as isoprene, contribute to only ~0.5% of ozone removal compared to ozone deposition to the human surface.

Understanding the Spatial Heterogeneity of Indoor OH and HO2 due to Photolysis of HONO Using Computational Fluid Dynamics Simulation

Indoor photolysis of nitrous acid (HONO) generates hydroxyl radicals (OH), and since OH is fast reacting, it may be confined within the HONO-photolyzing indoor volume of light. This study investigated the HONO-photolysis-induced formation of indoor OH, the transformation of OH to hydroperoxy radicals (HO₂), and resulting spatial distributions of those radicals and their oxidation products. To do so, a computational fluid dynamics (CFD) model framework was established to simulate HONO photolysis in a room and subsequent reactions associated with OH and HO₂ under a typical range of indoor lighting and ventilation conditions. The results showed that OH and HO₂ were essentially confined in the volume of HONO-photolyzing light, but oxidation products were relatively well distributed throughout the room. As the light volume increased, more total in-room OH was produced, thereby increasing oxidation product concentrations. Spatial distributions of OH and HO₂ varied by the type of artificial light (e.g., fluorescent versus incandescent), due to differences in photon flux as a function of light source and the distance from the source. The HO₂ generation rate and air change rate made notable impacts on product concentrations.