Combined Model for IAQ Assessment: Part 1—Morphology of the Model and Selection of Substantial Air Quality Impact Sub-Models

: Indoor air quality (IAQ) is one of the most important elements a ﬀ ecting a building user’s comfort and satisfaction. Currently, many methods of assessing the quality of indoor air have been described in the literature. In the authors’ opinion, the methods presented have not been collected, systematized, and organized into one multi-component model. The application purpose of the assessment is extremely important when choosing IAQ model. This article provides the state-of-the-art overview on IAQ methodology and attempts to systematize approach. Sub-models of the processes that impact indoor air quality, which can be distinguished as components of the IAQ model, are selected and presented based on sensory satisfaction functions. Subcomponents of three potential IAQ models were classiﬁed according to their application potential: IAQ quality index, IAQ comfort index, and an overall health and comfort index. The authors provide a method for using the combined IAQ index to determine the indoor environmental quality index, IEQ. In addition, the article presents a method for adjusting the weights of particular subcomponents and a practical case study which provides IAQ and IEQ model implementation for a large o ﬃ ce building assessment (with a BREEAM rating of excellent). To determine the indoor air quality and indoor environmental quality (IEQ) indexes for the case study. (TC—thermal comfort, L—light quality, ACc—acoustic comfort). index ), indoor air quality (IAQ index ), acoustics ( ACc index ), and lighting quality ( L index ). certified). Standardized CEN and ISO analytical methods were used to determine the VOC concentrations and CO 2 and formaldehyde concentrations in the indoor air of the building. Selection of the sampling points was made with the BREEAM assessor in two representative office zones per tested floor and a minimum of two floors. The building was tested three days after formal final finishing works at the pre ‐ occupancy stage with no users inside. For this office building, the tests were conducted on the 55th and 47th floors. Air samples were collected using an active sampling procedure with an electronic mass flow controller, which controlled the air flow (10 dm 3 /h for VOC tests and up to 30 dm 3 for formaldehyde tests). Indoor samples were set up in selected representative office locations, approximately 1.5 m above the floor, away from windows, doors, potential emission sources, and direct sunlight. Air samples were tested in accordance with the ISO 16000 ‐ 6:2011 and ISO 16000 ‐ 3:2011 standards. The VOCs were assessed using tubes filled with Tenax


State-of-the-Art Indoor Air Quality Measurement Systems
Approximately 30 years ago, people began to realize that buildings not only provide them with a sense of security, but can also significantly affect their health and well-being. This is particularly important due to the fact that people spend an increasing amount of time in closed indoor environment. Air quality and ventilation approaches were initially based on the users' dissatisfaction with the scent of the human body and, as such, the understanding of indoor air quality (IAQ) had serious limitations. Large quantities of pollutants and their sources clearly influence the indoor comfort of building inhabitants, as well as their health. In 1998, Fanger [1] presented an approach to the quantitative determination of perceived IAQ based on the level of dissatisfaction of residents caused by bad odors and irritants, smoke, and other sources of pollution. This approach provided two new measures of IAQ: the olf, which quantifies the pollution generated from a strong source of human bio-pollutant in the range of the impact of emitted odors on perceived air quality, and the decipol, measuring the perceived air quality in an indoor space with a source of pollution of one olf at a ventilation rate of 10 l/s. The number of emitted olfs per floor unit in different types of buildings and the amounts of pollutants from tobacco smoking (in olfs) can then be determined. Consideration of only the odour of the human body, without taking into account the influence of pollutants from various other sources (for example Smoking alone emits more than 7000 different compounds, many of which are harmful [2] for humans and animals and may transport biological pollutants that can act as allergens. People and household animals emit gases which are unpleasant, transfer pathogens, and cause diseases. These examples show that there are many paths for penetration of and exposure to the sources of pollution in indoor environments. In connection with the growing need to determine levels of indoor air pollution, new centers performing tests and new methods have been created considering the ability to analyze an increasing number of harmful substances. Measurement of pollutant concentrations in the air is generally a task performed by experts mainly in accredited laboratories and the results are published in scientific journals, technical reports, and, eventually, in guidelines, e.g., those of American Society of Heating Refrigerating and Air Conditioning Engineers ASHRAE [3]. The presence and concentrations of pollutants are often detected and measured without careful consideration of the significance of these measurements, and the pollutants measured may not be the most widespread or the most harmful. Some emissions are incorrectly grouped together; for example, more than one million volatile organic compounds (VOCs) are known and their toxicities are generally unknown, but they are often reported as a single value and referred to as the total VOCs (TVOCs) component. Frequently, carbon dioxide is used as an indicator of IAQ, although it does not have such a negative effect on the health of residents in the concentrations in which it is usually found in buildings. In our opinion, CO2 is rather a marker of human bioeffluents. Examples of different understandings of the set of typical pollutants in an indoor environment are shown in Smoking alone emits more than 7000 different compounds, many of which are harmful [2] for humans and animals and may transport biological pollutants that can act as allergens. People and household animals emit gases which are unpleasant, transfer pathogens, and cause diseases. These examples show that there are many paths for penetration of and exposure to the sources of pollution in indoor environments.
In connection with the growing need to determine levels of indoor air pollution, new centers performing tests and new methods have been created considering the ability to analyze an increasing number of harmful substances. Measurement of pollutant concentrations in the air is generally a task performed by experts mainly in accredited laboratories and the results are published in scientific journals, technical reports, and, eventually, in guidelines, e.g., those of American Society of Heating Refrigerating and Air Conditioning Engineers ASHRAE [3]. The presence and concentrations of pollutants are often detected and measured without careful consideration of the significance of these measurements, and the pollutants measured may not be the most widespread or the most harmful. Some emissions are incorrectly grouped together; for example, more than one million volatile organic compounds (VOCs) are known and their toxicities are generally unknown, but they are often reported as a single value and referred to as the total VOCs (TVOCs) component. Frequently, carbon dioxide is used as an indicator of IAQ, although it does not have such a negative effect on the health of residents in the concentrations in which it is usually found in buildings. In our opinion, CO 2 is rather a marker of human bioeffluents. Examples of different understandings of the set of typical pollutants in an indoor environment are shown in Table 1 [4]. This table provides recommended values from the results of the European project HealthVent [5], which aimed to develop health-based ventilation guidelines. Table 1 also includes recommendations provided by the Word Health Organization (WHO) on the acceptable levels of pollutant concentrations [6,7], as well as recommendations from other organizations, such as China's IAQ standard values [8]. Different approaches to the IAQ issue mean that the exposure limits assumed in the various source materials differ.  * PAH Para-Aminohippuric Acid-cyclic aromatic hydrocarbons; ** TVOC Total Volatile Organic Compounds; Reference [13] provides additional TVOC certification tests for new office buildings by determining (at the ppbv level) the content of carbon tetrachloride, chloroform, 1,2-and 1,4-dichlorobenzene, ethylbenzene, toluene, and o-, m-and, p-xylene.
Theoretical work on a combined IAQ model allowing aggregation of the results of the assessment of components affecting humans [14] is not yet well recognized in the literature. However, studies on IAQ indicators, which aim to provide a quantitative description of indoor air pollution, have been conducted since the nineties. In 2003, a significant study by Sekhar et al. [15] was published related to the standard indoor pollutant index (IPSI), the disease symptom index in the building symptom index (BSI), and to the often-cited works by Moschandres and Sofuoglu [16,17] on the indoor environmental index (IEI), indoor air pollution index (IAPI), and the indoor pollutant standard index (IPSI). The IAPI characterizes air pollution in an office with a single number: the index. The index value ranges between zero (lowest pollution level, i.e., best indoor air quality) and 10 (highest pollution level i.e., worst indoor air quality). The IAPI is a composite index; sub-indices ed are aggregated using the arithmetic mean in conjunction with a tree-structured calculation scheme. This scheme gives rise to some reservations, because at the top of the tree-structured calculation scheme is the IEI (calculated as the arithmetic mean of the IAPI and the IDI (indoor air discomfort index), and the combination of IAQ sensation and thermal conditions does not appear until later.
While considering the indicators for the quantitative description of pollution, the proposal of the IEA Working Group named "Defining the Metrics of IAQ" should also be mentioned. This group prepared, in 2017, the document entitled "In the Search of Indices to Evaluate the Indoor Air Quality of Low-Energy Residential Buildings" [18]. The group made the following assumption for the categorization of various indicators: there should be one index per individual pollutant and a dimensionless coefficient should be specified to evaluate the IAQ, provided that the current (observed) concentrations of a given pollutant c j are related to the ELVs (exposure limit values) concentration c j,ELV .
The index is calculated for each individual pollutant [18], which is specific only for this exact pollutant. The report showed that aggregation can be performed by addition, by taking the maximum value or by other methods, in an attempt to define metrics that can be used to evaluate IAQ.
The assumption was that the reference value usually refers to health risks (accounting for chronic or acute effects), but other metrics can also be used, (e.g., odor or irritation threshold). There are two important properties to be considered when aggregating sub-indices: ambiguity and eclipsing. As a result of the analysis, the authors concluded "that there are problems with model aggregation methods. In the aggregation model I agg = I 1 + I 2 , ambiguity creates a false alarm and in the aggregation model I agg = 1/2(I 1 + I 2 ), eclipsing underestimates the effect" [18]. Therefore, the discussion remains open [18]. The report also showed how there are large spreads of concentrations of individual pollutants (up to seven rows), even in the group of pollutants for which sub-indices were built. It determined the difficulties of building a weighted scheme based on the simplest percentage adjustment of the concentration shares and, thus, the share of the mass of pollutants to be removed by ventilation.
The current state of knowledge does not provide information authorizing the omission of certain pollutants. Hence, taking into account the lack of data on the characteristics of each chemical compound and consideration of the "removal efficiency" [19] requires us to abandon thinking about the adjustment of many individual pollutants, and to focus only on the creation of a model based on the representative and target components. In this state of knowledge, there are hopeful studies and proposals with a grey combined IAQ index model and the grey clustering model for IAQ indicators proposed by Zhu and Li in 2017 [20] is particularly interesting, especially when the relationships between system factors and the system's IAQ behavior and the interrelationships among the factors are uncertain. At first, all specific indoor air pollutants and related parameters should be measured. However, this is a very complex and time-consuming process. On the basis of the characteristics and correlations of the pollutants, the indoor air quality can be characterized by representative indicators. Studies [20] have pointed out that respirable particulates, CO 2 and TVOCs, were the three most representative and independent environmental parameters which can be used as an evaluation index of indoor air Appl. Sci. 2019, 9,3918 5 of 35 quality in office buildings. Since each indicator represents a class of pollutants with similar sources and dissemination characteristics, this index group avoids unreliability due to the fact that these indicators are "too small" because of critical concentration depression. A data pretreatment method must be used in the calculation procedure, reflecting the differences in concentration levels among different pollutants, but also expressing their influence on the comfort and health of the indoor occupants. Moreover, the measured pollutant concentrations can be used to predict the probable levels of other parameters, and good agreement was found between the predictions and measured values.

The Research Questions
The main research question contained in the paper concerned whether it was possible with the current state of knowledge to create and use in practice an IAQ model that was based on a unified and coherent approach for input indoor air parameters (such as pollutant concentrations, odor levels, and moisture content) and provided one output parameter (we proposed occupant satisfaction, IAQ index (in %)). The authors looked for physical equations for the IAQ index 's subcomponents and dependencies for their predicted occupant satisfaction functions with a pollutant concentration c j (PD = f (c j ) in %) which could be used as a model for subcomponents. This paper's intention was to provide an IAQ model with a step-by-step process which can be used to determine the value of the overall indoor environmental quality index (in %) including another three components: thermal comfort, acoustic comfort, and lighting quality. The innovative approach and added value of this article is in the use of the proposed IAQ model in practice and the relatively simple calculation of the overall IEQ value (with an uncertainty estimation) using the actual results of measurements in the Building Research Establishment Environmental Assessment Method BREEAM certified case study office building. The authors also provided occupant satisfaction functions for CO 2 , TVOCs, and formaldehyde HCHO in two variants: with experimental %PD values taken from the literature and for these pollutants' %PD values converted from an Air Quality Index system (see Section 2.2.).

Research Content and Strategy
The proposed IAQ model is presented in Sections 2.2-2.7. The model is later used to analyze the case study of an office building described in Section 2.9. Figure 2 presents the subsequent research steps from theory to practical application. Section 2.8 shows the method for determining indoor environmental quality IEQ index where IAQ index is a subcomponent/part of the IEQ index model. In order to determine the IAQ and IEQ, physical measurements of the indoor environment in the building were conducted using the experimental approach provided in Section 2.10. Based on these physical indoor measurements the IAQ and IEQ indexes (number of occupants satisfied with the indoor air and overall indoor quality, respectively) were assessed (see Section 3) and discussed (see Section 4).

The IAQ Model Proposal-Basic Assumptions
In the IAQ model construction process (our proposal), the commonly accepted approach is to transform individual concentrations of pollutants into subcomponents before they are aggregated into a single index (occupant satisfaction in %). However, summation of sub-indices can lead to situations in which all are under individual health thresholds, but the final indicator shows when the threshold has been exceeded. Conversely, the averaging of partial sub-indices can lead to an overall indicator showing an acceptable IAQ, even though one or more partial indicators are larger than their individual thresholds. One solution is to use the maximum value of all sub-indices to create the final form of the ∑IAQindex. Taking these issues into consideration, the authors created the ∑IAQ model with three complication levels adapted to the purposes of potential applications of the model, as presented in Figure 3; i.
Certification of a building, e.g., via the BREEAM system using (three sub-indices), called "quality"; ii.
Design, including perceptible contaminants affecting comfort and using the IAQ index when calculating the IEQ (five sub-indices), called "comfort"; iii.
Complex design, with the ∑IAQindex representing both comfort and health (seven sub-indices) called "comfort/health".
The simplest one, "quality", is an inner part of the ∑IAQcomfort model and can be used separately for simple applications with the main purpose of supporting a green building certification, e.g., via the BREEAM system using three components, i.e., CO2, HCHO, and TVOC. This model is used later on the case study of a BREEAM building.

The IAQ Model Proposal-Basic Assumptions
In the IAQ model construction process (our proposal), the commonly accepted approach is to transform individual concentrations of pollutants into subcomponents before they are aggregated into a single index (occupant satisfaction in %). However, summation of sub-indices can lead to situations in which all are under individual health thresholds, but the final indicator shows when the threshold has been exceeded. Conversely, the averaging of partial sub-indices can lead to an overall indicator showing an acceptable IAQ, even though one or more partial indicators are larger than their individual thresholds. One solution is to use the maximum value of all sub-indices to create the final form of the IAQ index . Taking these issues into consideration, the authors created the IAQ model with three complication levels adapted to the purposes of potential applications of the model, as presented in Figure 3; i.
Certification of a building, e.g., via the BREEAM system using (three sub-indices), called "quality"; ii.
Design, including perceptible contaminants affecting comfort and using the IAQ index when calculating the IEQ (five sub-indices), called "comfort"; iii.
Complex design, with the IAQ index representing both comfort and health (seven sub-indices) called "comfort/health".

The IAQ Model Proposal-Basic Assumptions
In the IAQ model construction process (our proposal), the commonly accepted approach is to transform individual concentrations of pollutants into subcomponents before they are aggregated into a single index (occupant satisfaction in %). However, summation of sub-indices can lead to situations in which all are under individual health thresholds, but the final indicator shows when the threshold has been exceeded. Conversely, the averaging of partial sub-indices can lead to an overall indicator showing an acceptable IAQ, even though one or more partial indicators are larger than their individual thresholds. One solution is to use the maximum value of all sub-indices to create the final form of the ∑IAQindex. Taking these issues into consideration, the authors created the ∑IAQ model with three complication levels adapted to the purposes of potential applications of the model, as presented in Figure 3; i.
Certification of a building, e.g., via the BREEAM system using (three sub-indices), called "quality"; ii.
Design, including perceptible contaminants affecting comfort and using the IAQ index when calculating the IEQ (five sub-indices), called "comfort"; iii.
Complex design, with the ∑IAQindex representing both comfort and health (seven sub-indices) called "comfort/health".
The simplest one, "quality", is an inner part of the ∑IAQcomfort model and can be used separately for simple applications with the main purpose of supporting a green building certification, e.g., via the BREEAM system using three components, i.e., CO2, HCHO, and TVOC. This model is used later on the case study of a BREEAM building.  The simplest one, "quality", is an inner part of the IAQ comfort model and can be used separately for simple applications with the main purpose of supporting a green building certification, e.g., via the Appl. Sci. 2019, 9, 3918 7 of 35 BREEAM system using three components, i.e., CO 2 , HCHO, and TVOC. This model is used later on the case study of a BREEAM building. Figure 1 shows the processes influencing the morphology of the IAQ model and Figure 3 provides a list of the pollutants for which three IAQ submodels were built, containing human-perceived contaminants (IAQ quality and IAQ comfort models), but also the IAQ comfort/health model for both perceptible and imperceptible pollutants, i.e., those that are not perceptible by humans but affect health and require additional energy for intensive ventilation for health reasons. There are potential sub-indices, such as IAQ(VOC non-odorous ) or IAQ(CO) [5,7]. Dust pollutants may have their sub-index both in the comfort model (if reliable curves of human sensory perception of PM concentrations are known) or in the IAQ comfort/health model if their health impact is considered to be the dominant feature. Considering the types of pollutants harmful to health assigned to the sub-indices of IAQ, we only consider the most important air pollutants (i.e., target emissions) that were given in the WHO guide in 2010 [5,7]. Submodels of processes that impact on air quality in indoor environments, which can be distinguished as components of the IAQ model, were based on sensory satisfaction functions (index of occupant dissatisfaction (PD) with the level of air pollution). Subcomponents of the three potential IAQ models were classified according to their future potential applications: in the assessment of environmental quality index IEQ (models IAQ quality and IAQ comfort ) or in the design of ventilation taking into account all possible harmful-to-health pollutants (model IAQ comfort/health ). In our opinion, such systematization creates order and has a practical dimension as presented later on in the case study. The following are the target pollutant groups: i.
In the air quality model (IAQ) quality , the IAQ index subcomponents were assigned to the selected three pollutants. The submodels for the IAQ were CO 2 , TVOC, and formaldehyde HCHO, as recommended by References [5,10,21,22]; ii.
In the (IAQ) comfort model, thepreviously provided simplified IAQ quality subcomponents for the three main pollutants were extended with a set of selected compounds VOC odorous , related to the collection of IAQ sub-indices (VOC odorous ) with an unknown cardinality, increased appropriately for the number of dominant pollutants. In addition, we provided a conditional deluge of two more components: (1) calculated using the enthalpy of hot and humid air (high enthalpy h > 55 kJ/kg [23]), the percentage of persons dissatisfied with respiratory cooling with humid air at relatively high temperatures and (2) the percentage of persons dissatisfied with indoor pollution with respect to dust pollution (PM 10 and PM 2.5 ), measured via panel tests. The introduction of a dust-pollution subcomponent to the IAQ model may be debatable, because some experimenters [24,25] underline the unique results of sensory tests of discomfort from dust, and the influence of "emissions" of respiratory dust particles on satisfaction is still under-researched. Considering the above, we expected two variants of the comfort model: with PD (PM 10 , PM 2.5 ) or without this factor; iii.
In the overall IAQ model, comfort (IAQ) comfort and health risk (IAQ) health indicators were used, and, hence, this model was called (IAQ) comfort/health . Models for subcomponents of IAQ not perceived by humans but influencing health, can be borrowed from the index set in the AQI (air quality index) system [26 -28], which was adapted to assess the quality of indoor air based on, and in accordance with, the concepts of the air quality assessment system used globally by the American EPA.
Values of AQI indices published on active EPA websites using the air quality index system were introduced for application in US federal regulations in 1999 [28]. Currently, the AQI system for outdoor air includes the following pollutants: ozone, particulate pollutants (PM 10 and PM 2.5 ), carbon monoxide CO, sulfur dioxide SO 2 , and nitrogen dioxide NO 2 . To convert a specific air pollutant concentration to an AQI, the EPA developed a tool called the AQI Calculator, which is an open resource [29]. This system (referring to the index from 2004 [16] and the indoor pollutant standard index (IPSI)) was further developed, and the proposed IAQI for indoor air presented by Wang et al. [30] in 2008 and a newer proposal [27] from 2017 for a similar but narrower set of indices, also for indoor air, were both modeled on it. The AQI and IAQI indicators showed an increase in the level of impact on human health with increasing concentrations of air pollution. There are some detected difficulties here, since "AQI is a piecewise linear function of the pollutant concentration" [27]. The calculated values of the AQI [31] or IAQI [30] indices, over the entire 0-500 scale calculated from the measured concentrations of selected contaminants or in the part of the scale corresponding to the IAQ rating, ranged from "good" to "unhealthy", and can be converted to PD% for use in the model equation, (IAQ) comfort/health . Concentrations will be significant when the uncertainties of scale conversions are estimated. Authors believe that their way of converting the AQI scale to PD% (which is similar to the method of conversion of the IEQ components' ordinal scales from the OFFICAIR EC project [32] to PD% scale; for example, the occupant percentage dissatisfied with noise [33,34] should be accepted in light of the expected results of a metrological analysis of the reliability of the combined IAQ model. Subcomponent models (physical functions, PD%) of the IAQ model for all individual air pollutants are presented later in this section.

IAQ Model Weighting Scheme Considering Air Pollution Ventilating
To obtain a comprehensive picture of IAQ in a building, it is necessary to measure the number of pollutants with different individual concentrations. There are methods that weigh sub-indices [21] but the problem is finding an effective weighting scheme and understanding how to adjust them in the overall model of all the pollutants in IAQ. For this reason, we proposed an adjustment method for the weights. In our opinion, provided in detail in References [33,34], and also according to Reference [22], the best weighting scheme, which would lead to a credibly aggregated model of IAQ composed of many extractable components (sub-indices), would be a system based on concentration values (the "excess masses" of pollutants to be regarded as loadings for the ventilation system). Therefore, the we aimed to determine the individual pollutants assigned to the IAQ model, their concentrations, c j , as the inputs of the IAQ submodels, and their "excess concentrations" originating from emissions or determined within indoor environments. Thus, it was possible to determine directly the energy requirements for ventilation purposes and the required minimum global ventilation rate. Determining the input concentration value, c j , for each IAQ sub-index enables the determination of the total mass of pollutants in the air, which is the basis for determining the air change rate N 1, . . . 7 (overall air change rate), assuming that the model includes all significant IAQ pollutants.
Currently, according to References [35][36][37], the most common assumption made is that pollution from VOC j compounds arises only from emissions due to the presence of construction or finishing materials (for j = 1, 2, . . . n) (it can be assumed that the source i of an emission is the entire indoor environment and then i = 1) from the zero state. The physical model for determining the ventilation rate in indoor environments polluted with VOC-type pollutants from building materials is given by EN 16798-1:2019 [10], assuming that design parameters for indoor air quality are derived using limit values for substance concentrations. In accordance with ECA Report Number 11 [36], the design ventilation rate required to dilute an individual substance emitted from building materials is calculated as: where Q h is the ventilation rate required for dilution in m 3 per second, G h is the emission rate of the substance in micrograms per second, C h,j is the guideline value of the substance in micrograms per m 3 , C h,o is the concentration of the substance in the supply air in micrograms per m 3 , and ε v is the ventilation effectiveness.
In fact, in a building with active "indoor chemistry" (see Figure 1), the use of this formula seems to be increasing. Taking into account the dynamic nature of the processes of generating various pollutants, the approach to IAQ and its components should be changed and subcomponents should be treated as pollution load processes, increasing in number not only due to the emission processes but also due to the generation of bio-pollution, water evaporation, and even dust infiltration from outside. It is also possible to set steady-state (initial) concentrations of pollutants and to determine the expected time courses of removal of these pollutants by means of ventilation (curves c j = f(τ)), at constant values of air change rate per hour ACH (h −1 ). Such ventilation rate calculations for CO 2 were developed in 1997 by Persily [38] and similar ones were provided in 2017 by Gyot [39] at the Berkeley National Laboratory. These calculations are not very accurate, as shown in the general demonstration graph in Figure 4.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 9 of 32` developed in 1997 by Persily [38] and similar ones were provided in 2017 by Gyot [39] at the Berkeley National Laboratory. These calculations are not very accurate, as shown in the general demonstration graph in Figure 4. Less accurate time-dependent curves of the total minimum ventilation rate ACH needed for "contaminant exhaustion" can be determined using programs [40] based on generic engineering equations for the sum of pollutants ∑cj. The generic equation for pollution concentration (the ratio of the amount of polluting product to the amount of fluid in the space (such as air in a room) can be calculated from the following equation: where c is the pollution concentration in the space (or in the room) with perfect mixing (m 3 /m 3 ) or (kg/kg), q is the amount of pollution added to the space (m 3 /h) (kg/h), N is the air change rate per hour (h −1 ), V is the volume or mass of the space (m 3 ) or (kg), e is the number 2.72, and t is time (h). If the initial concentration (at t = 0) in the space and the concentration in the supply fluid is zero, after some time the concentration in the room will stabilize. The ventilation rate graph for an amount of pollution q = 1 and a volume of space V =1, shows the values of ∑cj, similar to Figure 4. In order to obtain more exact values of the VOC concentrations remaining in the room for the air change rate function, it is possible to use the published dependencies or for the assumed volumes V of ventilated spaces with determinate concentrations, as they can be determined experimentally. A simplified method for determining the ventilation rate N from a simple formula for the time, t, course of a trial ventilation was provided by the Japanese researchers Noguchi et al. [41] in 2016. A description of this method is worth reading. Based on the temporal changes of the TVOC concentration measured using a PIDTVOC meter [42], the air change rate N or the ventilation rate F was estimated using the following method. Assuming perfect mixing of the air in the room and a constant TVOC emission rate E, the concentration change of TVOC in the room can be expressed by the following equation: Finally, Equation (5) can be expressed with one unknown parameter, the air exchange rate N as: Less accurate time-dependent curves of the total minimum ventilation rate ACH needed for "contaminant exhaustion" can be determined using programs [40] based on generic engineering equations for the sum of pollutants c j . The generic equation for pollution concentration (the ratio of the amount of polluting product to the amount of fluid in the space (such as air in a room) can be calculated from the following equation: where c is the pollution concentration in the space (or in the room) with perfect mixing (m 3 /m 3 ) or (kg/kg), q is the amount of pollution added to the space (m 3 /h) (kg/h), N is the air change rate per hour (h −1 ), V is the volume or mass of the space (m 3 ) or (kg), e is the number 2.72, and t is time (h). If the initial concentration (at t = 0) in the space and the concentration in the supply fluid is zero, after some time the concentration in the room will stabilize. The ventilation rate graph for an amount of pollution q = 1 and a volume of space V = 1, shows the values of c j , similar to Figure 4. In order to obtain more exact values of the VOC concentrations remaining in the room for the air change rate function, it is possible to use the published dependencies or for the assumed volumes V of ventilated spaces with determinate concentrations, as they can be determined experimentally. A simplified method for determining the ventilation rate N from a simple formula for the time, t, course of a trial ventilation was provided by the Japanese researchers Noguchi et al. [41] in 2016. A description of this method is worth reading. Based on the temporal changes of the TVOC concentration measured using a PID TVOC meter [42], the air change rate N or the ventilation rate F was estimated using the following method.
Assuming perfect mixing of the air in the room and a constant TVOC emission rate E, the concentration change of TVOC in the room can be expressed by the following equation: Finally, Equation (5) can be expressed with one unknown parameter, the air exchange rate N as: The initial concentration c j can be determined from the experimental results. After a long time, when the exponential term in Equation (5) can be assumed to be zero, the concentration C(t) becomes constant. The steady-state concentration C st can be determined from the temporal change in the experimental results where the concentration levels off.
The rule that the IAQ model should include a weighting scheme, referring to the variation in the share of pollutants in the IAQ, has been noted in References [35,37]. According to the first proposal, the weighting system is based on the differentiation of coefficients R j , which are the ratios of the real concentration values (or mass of pollutants) to the values of reference concentrations, representing the so-called relative masses of non-eliminated pollutants. This can also be represented by the desirable reduction in the level of pollution by means of ventilation and, thus, also by the energy requirements. For one emission source, the proposed system with a coefficient R j for a given pollutant C j has the formula: where I j is the ratio of the gas-phase concentration to the reference concentration value. For example, for the LCI (Lowest Concentration of Interest) value for the jth compound emitted from the building material, the factor R j is dimensionless, since y j is the gas-phase concentration for the jth compound in µg·m −3 , I j is the lowest concentration of interest (LCI) [41] for the jth compound in µg·m −3 and m is the number of all elected compounds. The weight coefficients W j are used as the weights of the equations in the IAQ index, according to: The authors of Reference [35] justified adjusting the coefficients where R j ≤ 1, but they did not explain the physical meaning of this condition. We believe that further discussion should include the issue of whether the "relative mass" of contamination expressed by Equation (6) has a proper place in the weighting scheme for the equations. The dimensionless quantity (7) does not have a sound physical meaning [37]. In our opinion, this type of calculation method is debatable.
One should strive to cover all the sub-indices of the combined IAQ ij model with a weighting scheme that would give VOCs a share in the total energy requirements for ventilation. The term "relative mass" should correspond to weights rationally proportional to the energy expenditure for ventilation of individual pollutants (IAQ sub-indices). Therefore, our future work will focus on introducing weights for all expressions in the overall IAQ model equation. However, since this is currently not possible without an adjustment method adapted to weight determination for very small concentrations, it was decided to present our model as an interim solution. We proposed the use of weights based on the "excess concentration" values within the pollutant categories only with similar and comparable orders of concentration values, for example the VOC odorous and VOC non-odorous categories. The rules of adjustment with boundary conditions will be provided and justified in the follow-up article to this report.

IAQ Model Scheme Morphology
According to the new proposal, for models with air quality sub-indices IAQ(P j ) (developed with the standard EN 16798-1:2019 as a reference for IEQ model creation [22,33] and with assumptions described in Reference [34]), in the case where indoor air has many pollutants, P 1 , . . . j, the combined ΣIAQ index equation is: (8) where the W P1 , . . . Pj weighting system for IAQ components is created on the basis of the arithmetic mean and the concept of "excess concentration" is introduced only for groups of pollutants with similar concentration values. There is a difference in concentration ∆c j between the observed concentration of pollutant c j and the reference concentration c ref (c ELV or c LCI ), which is below the current concentration in contaminated rooms. Thus, the excess concentration is: The weights W 1, . . . j for all three IAQ models are determined on the basis of arithmetic means or by adjusting all the values of ∆c j in a given model using Equation (10): where the sum of the adjusted weights W j of all ventilated pollutants described with sub-indices should be unity. The weight values for a given IAQ model, (e.g., IAQ comfort ) may be different, but the sum of the sub-index weights must be ≤1.0. The values of the reference concentrations are the concentration levels that are acceptable or recommended as limit values for various pollutants P j . In the case of the IAQ quality submodel as part of the IEQ index model, weights should be used for the VOC odorous (HCHO and TVOC) reference threshold concentrations of odors. The weights in the weighting system should be adjusted to unity according to the Equation (11): There are, however, non-typical cases in which the scales have different values. This is the case for formaldehyde, the concentration of which is many times lower in the building than the threshold level c th . According to the WHO [7], the admissible value c ref is also higher than the concentration in the building. In this case, the authors recommend taking the reference value as zero. Then, the weight W HCHO described by Equation (11) (for two pollutants), would not be negative (c HCHO − c ref ). The ASHRAE Guideline 10 (2011) [3] recommends that the IEQ model (and appropriate weights, W i ) should contain synergy effects of environmental parameters included in the subcomponents and their sensory perceptions. Figure 5 shows the extended IAQ index model with its sub-indices treated as components of the IEQ index , but also with sub-indices of the IAQ comfort/health type, i.e., pollutants that do not belong to the IEQ model but are important to health and the energy balance of a building with a mechanical ventilation system. The experimental dependencies of the percentage of persons dissatisfied, %PD, and the values of the concentrations of pollutants, c j , sensed in indoor air in the appropriate ranges are of fundamental significance in the sub-indices relevant to the IEQ model [43]. Appl. Sci. 2019, 9, x FOR PEER REVIEW 12 of 32` The scheme of the ∑IAQcomfort/health model consists of seven (or more) components or IAQ submodels and these are models for the various types of pollutants: IAQ(CO2), IAQ(TVOC), IAQ(HCHO), IAQ(VOCodorous), IAQ(h), IAQ (PM2.5, PM10), and the selected IAQ(VOCnon-odorous). The IAQ(VOCodorous) and IAQ(VOCnon-odorous) models should be multiplied, depending on the number of dominant VOC pollutants, and, hence, the ∑IAQcomfort/health model will, in practice, have more than seven components.
The inputs of each IAQ submodel are unit concentrations in air of a given pollutant,cj (in the case of IAQ(h). This is the moisture content x in well-known units "g of water vapor (gw) per kg of dry air (kga)", converted to a concentration of cj in μgwater/m 3 (or H-absolute humidity in gw/m 3 which is a measure of water-vapor density). In some cases, it is necessary to convert the pollution-derived parameter to VOC concentration (conversion of the odor intensity OI to VOC concentration is described later in this section).
From the concentration values, total air pollution can be calculated, and, subsequently, also the energy needed to ventilate the indoor air pollution. When the concentration levels of pollutants are variable and are increasing due to the presence of emissions, then the formulas given in Reference [36] and the amended standard EN 16798-1:2019 are used to calculate the required ACH ventilation rate. When the level of contamination is set (or quasi-fixed) and the volume and other parameters of the ventilated room are known, it is possible to calculate ventilation-time curves, i.e., maximum From the dependencies, expressed as the curves for PD(CO 2 ) or PD(VOC odorous ), the equations of the models are derived Equations (12)- (14): The scheme of the IAQ comfort/health model consists of seven (or more) components or IAQ submodels and these are models for the various types of pollutants: IAQ(CO 2 ), IAQ(TVOC), IAQ(HCHO), IAQ(VOC odorous ), IAQ(h), IAQ (PM 2.5 , PM 10 ), and the selected IAQ(VOC non-odorous ). The IAQ(VOC odorous ) and IAQ(VOC non-odorous ) models should be multiplied, depending on the number of dominant VOC pollutants, and, hence, the IAQ comfort/health model will, in practice, have more than seven components.
The inputs of each IAQ submodel are unit concentrations in air of a given pollutant,c j (in the case of IAQ(h). This is the moisture content x in well-known units "g of water vapor (g w ) per kg of dry air (kg a )", converted to a concentration of c j in µg water /m 3 (or H-absolute humidity in g w /m 3 which is a measure of water-vapor density). In some cases, it is necessary to convert the pollution-derived parameter to VOC concentration (conversion of the odor intensity OI to VOC concentration is described later in this section).
From the concentration values, total air pollution can be calculated, and, subsequently, also the energy needed to ventilate the indoor air pollution. When the concentration levels of pollutants are variable and are increasing due to the presence of emissions, then the formulas given in Reference [36] and the amended standard EN 16798-1:2019 are used to calculate the required ACH ventilation rate. When the level of contamination is set (or quasi-fixed) and the volume and other parameters of the ventilated room are known, it is possible to calculate ventilation-time curves, i.e., maximum ventilation curves for ACH ventilation rate to reach concentration levels ELV, LCI or the olfactory threshold level, according to Reference [40] or another adequate equation.
There are two outputs of each IAQ submodel as described below: 1.
The weights of the weighting system for the model IAQ quality and hypothetically W 1, . . . 5 for the model IAQ comfort or W 1, . . . , 7 for the model IAQ comfort/health (in a hypothetical model with a set adjustment method). These should reflect the energy load of the IAQ expressed by the theoretically assumed increase in the current concentration of pollutant c j relative to the reference concentration of c j , ref , which determines the level of this concentration intended to be obtained by ventilation.

2.
The PD% with the IAQ as a function of air pollution concentration. These values, determined in panel tests, reflect the impact of the interaction of air with a given pollutant at the actual level of concentration, estimated via panelists' sensations/perceptions (PD = f(c j ) in %).
Examples of measurable physical parameters for the purpose of IAQ and IEQ calculations (see case study) are given in the following section.
For the construction of the combined model IAQindex with a weighting scheme useful for aggregating sub-indices, we proposed the model presented in Figure 5. In this scheme, the combined IAQ model is shown as the basic assumption for the aggregation of all sub-indices. First, the model was cut by a cross-connected vertical connection regarding the inputs of submodels-the calculation of the sum of the masses of all pollutants c j in the ventilated space, using the values for the inputs of all submodels of IAQ concentration values of contaminants. The sum of the concentrations of all air pollutants expressed as mass units of pollution per m 3 of volume (which can be read after multiplication by V (m 3 ) as the mass to be displaced by ventilation), is the basis for calculating the "air change rate per hour" (ACH), the minimum air exchange rate needed to reduce the observed mass level of air pollution in a ventilated room (see Reference [40]). The second connection concerns the submodel outputs-the conversion of excess concentration to a dimensionless value, which allows to for the weighting scheme of the IAQ index combined model and the weights of individual IAQ submodels to be determined.
Additional assumptions were as follows: i. The sum ∆c 1, . . . 7 , which admittedly constitutes an excess mass increase of the sum of pollutants described by the submodels, was treated only as a "virtual energy load of the building" for ventilation conducted for the elimination of pollutants, and therefore, when adjusting the weights, the possibility of dividing the excess concentration by the sum of concentrations should be considered. ii.
The percentage of persons dissatisfied PD(IAQ component ) determined experimentally in sensory studies using panelists' sense of air quality during their exposure to an internal environment deteriorated by a given contamination component (PD = f (c j )), was derived from the literature or direct experiments.
Values of weights W 1, . . . j in sets of three, five or more components of the three IAQ index models (see Equations (12)- (14), were adjusted to a value of unity by dividing ∆c 1, . . . , j of each component by the sum of excess concentrations ∆c 1, . . . j in µ/m 3 . We proposed the use of c ref values, apart from the values of c LCI [18], c ELV , and the threshold values c th for odorous compounds, were as follows. i.
For IAQ(TVOC) and IAQ(VOC odorous ), the threshold concentrations, c th , for identified odorous compounds or mixtures are from Reference [44]; iii.
For the IAQ(h), the water-vapor concentration H (g w /m 3 ), recalculated from the moisture content x (g w /kg dry air ) using the gas constant for the water vapor and the actual temperature, the value of h up to the critical value for "high enthalpy of humid air" was evaluated using the formula: where h is the specific enthalpy of humid air (kJ/kg), which must be >55 kJ/kg. The EN 16798-1:2019 standard [10] recommends a limit for the dehumidification of air of 12g w /kg dry air (this value must be converted to a c value in g w /m 3 ).
i. For IAQ(VOC non-odorous ), the c ELV value in cases where no established LCI values were derived from the EN 16798-1:2019 standard. ii.
For IAQ(PM 2.5 , PM 10 ), the c ELV values were derived from the WHO [7] or other organizations (Tables 1 and 2).
The proposed reference values of pollutants forming the sub-indices of the IAQ model are given in Table 2. With reference to the concentration values of c LCI , it should be noted that according to References [18], this value is typically acquired by dividing occupational exposure limits by a safety factor (100 or 1000). Concentration c LCI is taken from Lowest Concentration of Interest (EU-LCI) from European Commission lists. However, the model values for exposure limit values (ELVs) of indoor air pollution, in accordance with the recommendations of the health-based ventilation guidelines [5], should be adopted in accordance with the current WHO guidelines given in the periodically issued WHO Air Quality Guidelines [7].

Selection of Submodels for Pollutant Components
Our overall selection of physical subcomponent equations and dependences for %PD = f(c j ) is presented in Table 3. The models presented were used to determine the IEQ for the sample building. The highlighted pollutants were taken into account in the case study building assessment.    Formulae for the conversion of odorant concentration c j to odor intensity OI for recalculation (4)  [53] [54] [55] [51] [34] The recommended criteria for dimensioning of humidification and dehumidification. It is recommended to limit the absolute humidity value to x = 12g w /kg a or H(g w /m 3 ) h = 1.006t + 0.622(2501 + 1.  1 TVOC, according to Reference [45], represents a narrow chromatographic picture that excludes, for example, the lower aldehydes, e.g., formaldehyde. 2 The measurement of the intensity of the odors in the building in which emissions from construction materials occur can be performed by a panel of participants, where the room is treated as a "test room for background odor" according to Section 6.8.1 of ISO 16000-28:2013, "Indoor Air-Part 28: Determination of odor emissions from building products using test chambers" (2013) [51]. The assessment of the 90% confidence level is possible through the use of a 15-pi odor intensity scale with a reading uncertainty of ±2 pi. 3 OI is perceived odor intensity on a six-level scale from 0 to 5 (no odor = 0, slight odor = 1, moderate odor = 2, strong odor = 3, very strong odor = 4, and overpowering odor = 5). 4 Based on the study Kim (5) volatile fatty acids. For example, for the group of reduced sulfur compounds for compound number 1, the conversion equation for H 2 S is Y = 0.950logX + 4.14, for the carbonyl compounds group for compound number 10, ammonium NH 3 , it is Y = 670logX + 2.38, and for the VOC group for styrene it is Y = 1.420logX + 3.10. 5 If one considers a cooling system that removes heat from a space but does not remove moisture unless condensation occurs, such as radiant cooling without dehumidification in a ventilation system, the importance of humidity is very clear. The sensible cooling of air in a room (no change in absolute humidity) from 25 • C and 60% RH to 20 • C (process a-b for x = 0.012 kg w /kg a ).This value, coupled with a high temperature, t a = 25 • C, is accepted as the critical limit value for dehumidification by EN16798-1:2019 [10] and can be converted into an absolute humidity H (g w /m 3 ). These processes are expected to significantly increase thermal comfort and air quality. Nevertheless, the same change in enthalpy of the air h can be achieved by simply reducing the humidity by 10% RH and keeping the temperature constant (processes a-c). Since IAQ is a function of enthalpy, these are expected to be the same. Here, a change in humidity of 10% RH at constant temperature is equivalent to a change in temperature of 5 or 6 • C at constant moisture content in air x (kg w /kg dry air ). 6 The IAQI system [30]  concentration value c using the interpolation method for each air pollutant. In the IAQI system, the index range 0-50 is "good" with a significance level of "little or no risk", 51-100 is "moderate" where "sensitive persons or those with respiratory symptoms are concerned", 101-150 is "unhealthy for sensitive groups", 151-200 is "unhealthy for all individuals", 201-300 is "very unhealthy-more serious health effects for everyone for short-term exposure", and 301-500 is "hazardous" with "health warning of emergency conditions for everyone".
Therefore, the comfort scale for IAQI is adequate for index values from zero to 200 and our proposal is to use the "hypothetical" scale of PD* of 0-100% in this range, converted from IAQI. 7 The method of conversion of the IAQI scale to the PD* scale is based on experiences gained during research [27]. The IAQI values in a health-risk scale can be given in % for persons giving a verbal answer of "no risk", "moderate", and "unhealthy" as their health risk evaluation. Therefore, when determining, for instance, The air quality indexes (i.e., AQI [31] and IAQI [30]) are piecewise linear functions of the pollutant concentrations. At the boundary between AQI categories, there is a discontinuous jump of one AQI unit. To convert from concentration c j to I j (in the converted scale index I j will be PD j , Equation (16) is used.
where I j is the air quality index I in the PD*% scale, c low is the pollutant concentration break point, which is ≤c j , c high is the pollutant concentration break point, which is ≥c j , I low is the index break point corresponding to c low , I high is the index break point corresponding to c high and c p is the truncated (to an integer) actual concentration for the pollutant. Little data exist on the AQI's metrological reliability for AQI and IAQI. Only in the EPA Air Program undertaken at Cornell University [26] is there a previous review of the quality assurance requirements for AQI.

The Representative VOCs for Indoor Environment
The time when IAQ studies focused on a class of contaminants referred to as volatile organic compounds (VOCs) is bygone. The analytical methodology available was the primary basis for this focus, but the recent broadening of analytical methods has led to growing realization that other compounds (i.e., SVOCs) beyond traditional VOCs are implicated in IAQ problems. The choice of VOCs remains a challenge in IAQ assessment. Moreover, VOCs is somewhat vague term, the definition of which is not universally agreed upon. It has been defined in terms of vapor pressures and boiling points, as well as molecular chain lengths detectable by chromatographic techniques. Due to the complexity of VOC emission profiles, it is tempting to simplify the analysis and reporting of emissions by grouping all detected compounds together. The first problem with this approach is that individual compounds have highly variable health and/or comfort effects, the result being that concentration alone is not predictive of IAQ effects. Levels of concern vary by orders of magnitude, so a collective concentration will not correlate with IAQ. Second, VOC detection and quantification are highly method dependent. A given sampling and analysis system cannot capture or respond to all the VOCs present in any indoor environment or in the test chamber for a given emitting material. Thus, the term "total" is misleading. The important aspect of IAQ submodel selection is the strategy defined by the US EPA as "VOCs-Total versus Target: Irritancy, Odor and Health Impact". The representative 90 target VOCs were presented by Canada's National Research Council's Institute for Research in Construction (NRG-IRC) in collaboration with several academic and governmental partners, including Health Canada. The compounds were selected based on health impact, occurrence in indoor air, known emission from building materials, as well as suitability for detection and quantification by gas chromatography-mass spectrometry (GC-MS) or high-performance liquid chromatography (HPLC). Our list of target VOCs was actually representative for indoor environment and recommended by the HealthVent project and is provided in Table 1.

Steps of ΣIAQ index Calculation
After selection of the IAQ model type ( IAQ quality or IAQ comfort ), the IAQ index evaluation was carried out using the complex model IAQ from Figure 5, which should contain the following stages.
(a) Calculation of the total concentration of pollutants in the ventilated space or the total mass of air pollutants per m 3 , the level of which is to be reduced by the ventilation process (taking into account the ventilated volume of the room and the emissions present). (b) Selection of the IAQ index model shape from the models defined by Equations (12)- (14), with the provision that due to the multiplication of the submodels for IAQ(VOC odorous ), the number of subcomponents of the IAQ comfort model will be more than five. (c) Processing the input data of the submodels to obtain the concentration value c j , e.g., converting the measured OI value into a concentration value for a given pollutant c j in µg/m 3 .
(d) Calculation of the excess concentration values for each identified contaminant (Table 2) ∆c j = c j −c ref .
(e) Calculation of the sum of excess concentrations (see Table 2), ∆c j . (f) Calculation of adjusted weights Wj for the selected model equations. IAQ quality and IAQ comfort , are determined on the basis of arithmetic means or by adjusting all the values of ∆c j in a given model using Equation (10), only for groups of pollutants with similar concentration values. (g) Calculation of the value of the ventilating air flow for the environment described in the IAQ index model in accordance with the requirements of the standard EN 16798-1:2019 (a method using the criteria for the ventilation required for the individual substance emitted) [10]. (h) Calculation of given IAQ environmental input parameters, including concentrations of pollutants c j assigned to submodels. The PD values from their sensory equations (Table 3) are presented as the dependence of the percentage of persons dissatisfied PD = f (c j , . . . ), from one of the formulas from Table 3, in order to determine this function. (i) Selection of the IAQ quality model equation (with weights W 1 , W 2 , and W 3 ) or the IAQ comfort model equation (with weights W 1 , . . . W 5 or more) and calculation from Equation (13) of the value with adjusted weights, followed by multiplication of IAQ submodels (Equation (8)) and insertion as a term of the IEQ index in Equation (18) [34].
When selecting the IAQ comfort/health model type, an IAQ index evaluation is carried out using the combined model IAQ from Figure 5, which should contain the following steps.
(a) Calculation of the total concentration of pollutants c j in the ventilated space or the total mass of air pollutants per m 3 , the level of which is to be reduced by the ventilation process (taking into account the ventilated volume of the room and the emissions present). (b) Choosing the IAQ index model (12) from among the models defined by Equations (12)- (14), with the provision that by multiplying the IAQ(VOC non-odorous ) submodels, the number of subcomponents of the IAQ comfort/health model will be more than seven.  (Table 3) depending on the percentage of persons dissatisfied, PD = f (c j , . . . ), or selected from the formulas for determining this function given in Table 3. (i) Development of IAQ PM2.5 , IAQ PM10 , and IAQ non-odorous submodels. When it is planned to use the indoor air quality index scale IAQI or a similar scale, it is necessary to convert these to PD* values in %, in two steps: (1) by reading from the standard curves of the IAQI = f(c j ) all IAQI values for the determined (measured) VOC non-odorous concentration values and using the converted scale calibration curve, PD = f(IAQI), by reading from the recalibration curve (Equation (16)), the PD = f(c j ) values on the dissatisfaction rating scale from zero to 100 in %, according to Footnote 7 in Table 3.

The IEQ Assessment Equation with IAQ As a Subcomponent
The proposed IAQ model can be a substantial component of the IEQ model; for example, in the case study shown later in the article. The indoor environmental quality index refers to the quality of a building's environment with respect to the occupants' satisfaction in %. The morphology of the IEQ index model used to assess buildings, to determine as an IEQ component the IAQ index and to determine other subcomponents TC index -thermal comfort, ACc index -acoustic comfort and L index -light quality based on measurements of physical properties in each of the submodels-in accordance with the scheme of the Piasecki-Kostyrko model, is presented in Figure 6 [22]. Lindex-light quality based on measurements of physical properties in each of the submodels-in accordance with the scheme of the Piasecki-Kostyrko model, is presented in Figure 6 [22]. The EN 16798-1:2019 is the reference for IEQ model creation [22,33]. The standard allows complex indoor information to be presented as one overall indicator of indoor environmental quality of the building-IEQindex. The model reliability, including the uncertainties of measurements and data for this model, was discussed clearly in Reference [34], where the authors also presented the internal incongruity in the IEQ model structure and the justification for using the crude weights method for each subcomponent. Originally, the IEQ model was expressed as a polynomial equation consisting of four terms by Wong [43]. The IEQindexis composed of the following subcomponents (SIi): thermal comfort (TCindex), indoor air quality (IAQindex), acoustics (ACcindex), and lighting quality (Lindex). Multiplying their weights, Wi, leads to Equation (17).
The authors adopted the crude weighting system, where all elements are weighted in the same way (0.25 for W1-W4), as shown in Equation (18).
As a consequence of the equation, the subcomponents SIi (the predicted percentage of those satisfied) can be calculated using Equation (19).
where PD is the predicted percentage dissatisfied (PPD) and PD(SIi) is the percentage of persons dissatisfied with the IEQ subcomponent (SIi) level. The authors' simulations for IEQindex sub-indices and preliminary metrological analysis of the overall IEQ model fitting were performed with Monte Figure 6. The research steps necessary to determine the IEQ index for buildings, including physical and design parameters of buildings and subcomponent models.
The EN 16798-1:2019 is the reference for IEQ model creation [22,33]. The standard allows complex indoor information to be presented as one overall indicator of indoor environmental quality of the building-IEQ index . The model reliability, including the uncertainties of measurements and data for this model, was discussed clearly in Reference [34], where the authors also presented the internal incongruity in the IEQ model structure and the justification for using the crude weights method for each subcomponent. Originally, the IEQ model was expressed as a polynomial equation consisting of four terms by Wong [43]. The IEQ index is composed of the following subcomponents (SI i ): thermal comfort (TC index ), indoor air quality (IAQ index ), acoustics (ACc index ), and lighting quality (L index ). Multiplying their weights, W i , leads to Equation (17).
The authors adopted the crude weighting system, where all elements are weighted in the same way (0.25 for W 1 -W 4 ), as shown in Equation (18). As a consequence of the equation, the subcomponents SI i (the predicted percentage of those satisfied) can be calculated using Equation (19).
where PD is the predicted percentage dissatisfied (PPD) and PD(SI i ) is the percentage of persons dissatisfied with the IEQ subcomponent (SI i ) level. The authors' simulations for IEQ index sub-indices and preliminary metrological analysis of the overall IEQ model fitting were performed with Monte Carlo tests. It is easy to show that the standard deviations of these values are equal:

A Case Study of a Building
The experimental part of this study was performed simultaneously with the BREEAM certification process, including determination of the three primary IAQ pollutants: formaldehyde concentration, CO 2 , and VOCs in the indoor air [22]. The building is a high tower, made of a convex concrete-steel structure with a glass facade. The basic information on the assessed building is presented in Table 4. At the time of the test, the building had a standard empty office without furniture (so-called pre-occupancy stage). The walls were plastered and painted, the suspended ceilings were in place, and the floors were finished with synthetic carpets. All building installations were active, including the mechanical ventilation controlled by the Building Management System BMS system with zonal CO 2 concentration sensors. The building was tested a few days after the formal end of finishing works.
The tests were made on the 55th and 47th floor. The experimental part of this study was performed simultaneously with the BREEAM certification process, including determination of the three primary IAQ pollutants: formaldehyde concentration, CO2, and VOCs in the indoor air [22]. The building is a high tower, made of a convex concrete-steel structure with a glass facade. The basic information on the assessed building is presented in Table 4. At the time of the test, the building had a standard empty office without furniture (so-called pre-occupancy stage). The walls were plastered and painted, the suspended ceilings were in place, and the floors were finished with synthetic carpets. All building installations were active, including the mechanical ventilation controlled by the Building Management System BMS system with zonal CO2 concentration sensors. The building was tested a few days after the formal end of finishing works. The tests were made on the 55th and 47th floor. Measurement points in the building were determined based on the analysis of frequencies of designed occupancies of the room and interior finish standards (open spaces). The sampling plan was prepared with the BREEAM assessors conducting the certification process of the facility. The main focus was on the IAQ index of open spaces in which the largest number of people may reside, and these represent the largest occupied usable floor space. According to the detailed design project documents, the building emphasizes the use of materials with known and low emission levels (BREEAM certified).

The Equipment, Measurements, and Experimental Approach
Standardized CEN and ISO analytical methods were used to determine the VOC concentrations and CO2 and formaldehyde concentrations in the indoor air of the building. Selection of the sampling points was made with the BREEAM assessor in two representative office zones per tested floor and a minimum of two floors. The building was tested three days after formal final finishing works at the pre-occupancy stage with no users inside. For this office building, the tests were conducted on the 55th and 47th floors. Air samples were collected using an active sampling procedure with an electronic mass flow controller, which controlled the air flow (10 dm 3 /h for VOC tests and up to 30 dm 3 /h for formaldehyde tests). Indoor samples were set up in selected representative office locations, approximately 1.5 m above the floor, away from windows, doors, potential emission sources, and direct sunlight. Air samples were tested in accordance with the ISO Appl. Sci. 2019, 9, x FOR PEER REVIEW 21 of 33` The experimental part of this study was performed simultaneously with the BREEAM certification process, including determination of the three primary IAQ pollutants: formaldehyde concentration, CO2, and VOCs in the indoor air [22]. The building is a high tower, made of a convex concrete-steel structure with a glass facade. The basic information on the assessed building is presented in Table 4. At the time of the test, the building had a standard empty office without furniture (so-called pre-occupancy stage). The walls were plastered and painted, the suspended ceilings were in place, and the floors were finished with synthetic carpets. All building installations were active, including the mechanical ventilation controlled by the Building Management System BMS system with zonal CO2 concentration sensors. The building was tested a few days after the formal end of finishing works. The tests were made on the 55th and 47th floor. Measurement points in the building were determined based on the analysis of frequencies of designed occupancies of the room and interior finish standards (open spaces). The sampling plan was prepared with the BREEAM assessors conducting the certification process of the facility. The main focus was on the IAQ index of open spaces in which the largest number of people may reside, and these represent the largest occupied usable floor space. According to the detailed design project documents, the building emphasizes the use of materials with known and low emission levels (BREEAM certified).

The Equipment, Measurements, and Experimental Approach
Standardized CEN and ISO analytical methods were used to determine the VOC concentrations and CO2 and formaldehyde concentrations in the indoor air of the building. Selection of the sampling points was made with the BREEAM assessor in two representative office zones per tested floor and a minimum of two floors. The building was tested three days after formal final finishing works at the pre-occupancy stage with no users inside. For this office building, the tests were conducted on the 55th and 47th floors. Air samples were collected using an active sampling procedure with an electronic mass flow controller, which controlled the air flow (10 dm 3 /h for VOC tests and up to 30 dm 3 /h for formaldehyde tests). Indoor samples were set up in selected representative office locations, approximately 1.5 m above the floor, away from windows, doors, potential emission sources, and direct sunlight. Air samples were tested in accordance with the ISO 16000-6:2011 and ISO 16000-3:2011 standards. The VOCs were assessed using tubes filled with Tenax pre-occupant 49 59,000 3000 2 Measurement points in the building were determined based on the analysis of frequencies of designed occupancies of the room and interior finish standards (open spaces). The sampling plan was prepared with the BREEAM assessors conducting the certification process of the facility. The main focus was on the IAQ index of open spaces in which the largest number of people may reside, and these represent the largest occupied usable floor space. According to the detailed design project documents, the building emphasizes the use of materials with known and low emission levels (BREEAM certified).

The Equipment, Measurements, and Experimental Approach
Standardized CEN and ISO analytical methods were used to determine the VOC concentrations and CO 2 and formaldehyde concentrations in the indoor air of the building. Selection of the sampling points was made with the BREEAM assessor in two representative office zones per tested floor and a minimum of two floors. The building was tested three days after formal final finishing works at the pre-occupancy stage with no users inside. For this office building, the tests were conducted on the 55th and 47th floors. Air samples were collected using an active sampling procedure with an electronic mass flow controller, which controlled the air flow (10 dm 3 /h for VOC tests and up to 30 dm 3 /h for formaldehyde tests). Indoor samples were set up in selected representative office locations, approximately 1.5 m above the floor, away from windows, doors, potential emission sources, and direct sunlight. Air samples were tested in accordance with the ISO 16000-6:2011 and ISO 16000-3:2011 standards. The VOCs were assessed using tubes filled with Tenax adsorbent. Then, they were thermally desorbed using a thermal desorption apparatus (TD-20, Shimadzu, Tokyo, Japan). The process of separation and analysis of volatile compounds was achieved using a gas chromatograph equipped with a mass spectrometer (GC/MS) (model: GCMS-QP2010, Shimadzu, Tokyo, Japan). The following GC oven temperature program was applied: initial temperature 40 • C for five min, 10 • C per min to 260 • C, and the final temperature of 260 • C for 1 min. The 1:10 split ratio injection mode was applied. The method used has a limit of quantification of 2 µg/m 3 . The volatile compounds were identified by comparing the retention times of chromatographic peaks with the retention times of reference compounds and by searching the NIST data base (National Institute of Standards and Technology, Gaithersburg, MD, USA) mass spectral database. Identified compounds were quantified using a relative identification factor obtained from standard solution calibration curves. TVOC was calculated by summing identified and unidentified compounds eluting between n-hexane and n-hexadecane. In order to determine volatile aldehydes, air samples were taken via cassettes using a solid absorbent silica gel coated with 2,4-dinitrilophenyl hydrazine (2,4-DNPH), and then subjected to a laboratory test using high-performance liquid chromatography (HPLC) with UV-Vis detection (Dionex 170S, Dionex, Sunnyvale, CA, USA) and an isocratic pump (Dionex P580A, Dionex, Sunnyvale, CA, USA). The described method has a limit of quantification at 2 µg/m 3 .
Other IEQ index components were tested as follows. The acoustic tests confirming the designed values were carried out by the measurement of the equivalent sound levels, LAeq, in the selected locations. The measurements were carried out during the daytime (starting at 11:00). The following equipment was used for the measurements: Brüel&Kjaer 4231 acoustic calibrator (Brüel&Kjaer, Naerum, Denmark), Nor-121 analyzer (Norsonic, Tranby, Norway), Brüel&Kjaer 4165 measuring microphones (Brüel&Kjaer, Naerum, Denmark), analyzer with microphone Norsonic-140 (Norsonic, Tranby, Norway). Before the tests were carried out, the calibration of the measuring path was conducted in accordance with the instructions to "check the acoustic measurement channel". The test results were evaluated in relation to the requirements considering permissible sound levels A in rooms intended for human dwellings. Thermal environmental measurements were provided using the microclimate multifunctional instrument HD32.1 and the tests were in accordance with ISO 7726 and ISO 7730. VOCs were tested simultaneously at all points. Visual comfort (Hea 01) was confirmed by using a MAVOLUX 5032C instrument (USB version) with a 3C15683 detector (Gossen, Nürnberg, Germany), in accordance with EN 12464 provisions.

Additional Explanations
The adaptation of the IAQ model to a practical casestudy was mainly for illustrative purposes in the context of the presented IAQ calculation/aggregation method. We did not focus deeply on discussing the technical or environmental issues of the presented building. Other IEQ subcomponents, such as thermal, acoustic, and visual satisfaction (in %), used to determine the IEQ index, were experimentally determined and partly presented in References [22,33]. Authors do not focus on these results in this article, as they have already been discussed in other papers [22].

Results for theIAQ index and IEQ index Prediction
A previous publication of ours [22] reported on IEQ and IAQ building assessments for a larger number of BREEAM buildings, where IEQ was assessed without calculating the combined IAQ index. The combined model of the IAQ index presented in this paper had not yet been previously developed, and we were limited in determining the IEQ index , thus we only took into account two of the most well-known pollutants (i.e., CO 2 and TVOC) separately. The assessment of the IEQ index was made by adaptation of the measured parameters (complying with the draft EN 16798-1:2019 standard for indoor environments) as the input values for the submodels of the IEQ index . The input values for the case study are presented in Table 5, which provides the input data for determining the IEQ index sub-indices of thermal comfort (TC index ), indoor air quality (IAQ index ), acoustics (ACc index ), and lighting quality (L index ) for an office building (47th floor) three days after completion of the finishing work before users were allowed in the building (i.e., pre-occupancy stage). Table 5. Physical parameters 1 and IEQ index results calculated using Equation (9) separately for an IAQ index with internal air pollution of CO 2 and an IAQ index with internal TVOC air pollution, assuming a realistic uncertainty of parameter measurement for the case study of a building (47th floor; open space) three days after the completion of finishing works.

IEQ TVOC
Second variant with c TVOC as an IAQ index parameter IEQ TVOC = 80.1% ± 10.7% 1 The IEQ and its measurement's uncertainty (with subcomponent standard deviation values) were calculated for IEQ physical parameter values, where t a is the air temperature ( • C), t r is the mean radiant temperature ( • C), v a is the relative air velocity (m/s), p a is the water-vapor partial pressure (Pa), M is the metabolic rate (met), and I cl, is the clothing insulation (clo). In addition, c CO2 is the concentration in ppm, c TVOC is the highest observed TVOC concentration in µg/m 3 , actual noise is in dB(A), and E min is the minimum daylight illuminance (lux).

Results for the ΣIAQ index and IEQ index Assessment Including Identified Pollutants (CO 2 , TVOC, and HCHO)
The example of a modified calculation of the collective submodel IAQ quality for three basic pollutants, as a component of the IEQ model for determining one project value for this indicator, is provided in two variants. The first variant uses sub-indices of IAQ for two pollutants, CO 2 and TVOC, which are described in Table 5, as well the sub-index of the third pollutant, HCHO (according to Reference [47]), where these differences in the approaches mean one must combine them into one submodel IAQ in order to be used in IEQ calculation. The second variant uses submodels of IAQ for TVOC and HCHO pollutants based on the IAQI system [30] and then converts them into percentages of persons dissatisfied PD* in %. According to the diagram of the model IAQ from Figure 5 and using Equation (12) of the IAQ quality model, the submodel weights are calculated as follows.
W CO 2 for the submodel IAQ(CO 2 ) = 0.5 is a component of the polynomial: W VOC for submodel IAQ(VOC) = 0.5 is a weight for combined submodel of the polynomial: IAQ(VOC) = W TVOC ·IAQ(TVOC) + W HCHO ·IAQ(HCHO) (22) with the terms W TVOC and W HCHO calculated from Equation (14) using the measured values c j (actual concentration of TVOC and HCHO) and the reference values c ref (Table 6). In our case study, the value of the submodel IAQ weight for the model (HCHO) ought to also be calculated from the measured value and the reference value c ref . However, formaldehyde is an unusual pollutant because, although it belongs to VOC odorous compounds, the concentrations found in buildings are many times lower than the HCHO threshold c th = 300 µg/m 3 according to the WHO [7] and lower than the threshold concentrations of HCHO from 60 µg/m 3 to 70 µg/m 3 issued in 2013 by the American Industrial Hygiene Association [44]. The permissible value of c ref = 100 µg/m 3 is also higher than the formaldehyde concentration found in buildings, according to Reference [5] and the standard EN16798-1:2019 [10]. Therefore, the authors propose that in such a case (to avoid a negative value of ∆c j ), the value modelling the reference should be taken as zero. Then, the form of the adjusted W HCHO weight in the model described by Equation (11) for air with three pollutants, would be as follows.
The results of the weights assessment for the two variants of the IAQ quality model are presented in Table 6.
According to the diagram of the model IAQ from Figure 5 and Equation (8), we proposed sensory equations for the percentage of persons dissatisfied %PD* in two variants.
The submodel IAQ's first variant includes the following: 1.
The IAQ submodels used so far in References [22,33] for CO 2 and TVOC pollutants, as shown in Table 5; 2.
The IAQ submodel for formaldehyde, using two types of equations depending on the range of HCHO concentrations measured in the building. Formaldehyde concentrations in the air with values above the threshold concentration, c th , for its odor, i.e., above 60 or even 300 µg/m 3 , can be used to create IAQ submodels for rooms with volatile and aromatic VOC compounds as well as for the HCHO equation [50].
PD HCHO = exp (2.14 · OI − 3.81) exp (2.14 · OI − 3.81) + 1 (24) However, in the case study building, the maximum concentration of HCHO was 18 µg/m 3 and, therefore, its concentration in the air was several times lower than the concentration of the odor threshold, c th [44]. The intensity of the formaldehyde odor was undetectable under these conditions, and the sensory equation PD = f(OI), which is appropriate for sensory detection of IEQ, is not applicable for odors below the threshold. Therefore, for small concentrations, we proposed the use of the equation taken from the work of Zhu and Li [47] based on the analysis of "health effects on the human body", derived from "indoor air quality comfort evaluation experiments and the literature".
PMV HCHO = 2log c HCHO 0.01 (25) This equation links the value of the new unit "the effect of formaldehyde on human comfort", called PMV HCHO , with its c HCHO concentration (µg/m 3 ) in the air. It covers the range from 10 µg/m 3 to 320 µg/m 3 and, as declared by the authors, this value has the same nature as PMV thermal comfort, which can be converted into a PD% unit according to the formula in Reference [48], experimentally confirmed for nearly zero energy buildings (NZEBs) by Reference [58].
The submodel IAQ's second variant includes the following.
i. The IAQ(CO 2 ) submodel used so far for CO 2 pollution, as shown in Table 5. ii.
The IAQ submodels for TVOC and HCHO types of pollution used as indoor air quality index ratio values borrowed from the IAQI system [30], which are then converted into percentages of persons dissatisfied (PD* in %) in the following way (a) The reference curves of IAQI = f (c j ) [30] for two dependencies of the IAQI index on TVOC and HCHO contamination values must be reconstructed. On the y-axis are the IAQI index values from zero to 200 in the range from "no risk" to "unhealthy" and on the x-axis, the c j values are presented. (b) In accordance with the measured values of c TVOC and c HCHO , the values IAQI TVOC and IAQI HCHO are determined from the functions IAQI TVO c = f(c) and IAQI HCHO = f(c). (c) Based on the IAQI system parameters [30] given in Footnote 6 of Table 3, which are presented as data for the functions for indexes, the IAQI TVOC and IAQI HCHO values appropriate for the break pointsi n perceived pollution concentration values are calculated in a range from zero (good) to 200 (unhealthy) using the ordinal scale IAQI = f(c) [27]. This function was converted to IAQI TVOC and IAQI HCHO scales using the concentration function scales PD*(TVOC) and PD*(HCHO) in the PD* range from 0 to 100% ( Figure 8). (d) The data used for calculation and conversion of IAQI and PD* scales are presented in Table 7. (e) Based on data determined for the new converted scales (Table 7)   (e) Based on data determined for the new converted scales (Table 7)  In the context of the results for the overall IEQ model index when treating both the main pollutants CO2 and TVOC separately, we present in Table 8 the transformed IEQ calculations with the sub-indices ∑IAQ. The results for the individual case study building IAQ subcomponents are taken from Table 5 (for cTVOC = 787 μg/m 3 and cHCHO = 18 μg/m 3 [22]). The IEQ index with ∑IAQquality values was calculated using two variants-the first conventional and the second borrowed from the IAQI scale [30] (Table 8).  which can be converted into a PD% unit according to the formula in Reference [48], experimentally confirmed for nearly zero energy buildings (NZEBs) by Reference [58].
The submodel ∑IAQ's second variant includes the following.
i. The IAQ(CO2) submodel used so far for CO2 pollution, as shown in Table 5. ii.
The IAQ submodels for TVOC and HCHO types of pollution used as indoor air quality index ratio values borrowed from the IAQI system [30], which are then converted into percentages of persons dissatisfied (PD* in %) in the following way (a) The reference curves of IAQI = f(cj) [30] for two dependencies of the IAQI index on TVOC and HCHO contamination values must be reconstructed. On the y-axis are the IAQI index values from zero to 200 in the range from "no risk" to "unhealthy" and on the x-axis, the cj values are presented.
(b) In accordance with the measured values of cTVOC and cHCHO, the values IAQITVOC and IAQIHCHO are determined from the functions IAQITVOc= f(c) and IAQIHCHO = f(c). (c) Based on the IAQI system parameters [30] given in Footnote 6 of Table 3, which are presented as data for the functions for indexes, the IAQITVOC and IAQIHCHO values appropriate for the break pointsi n perceived pollution concentration values are calculated in a range from zero (good) to 200 (unhealthy) using the ordinal scale IAQI = f(c) [27]. This function was converted to IAQITVOC and IAQIHCHO scales using the concentration function scales PD*(TVOC) and PD*(HCHO) in the PD* range from 0 to 100% (Figure 7). (d) The data used for calculation and conversion of IAQI and PD* scales are presented in Table 7.  In the context of the results for the overall IEQ model index when treating both the main pollutants CO 2 and TVOC separately, we present in Table 8a the transformed IEQ calculations with the sub-indices IAQ. The results for the individual case study building IAQ subcomponents are taken from Table 5 (for c TVOC = 787 µg/m 3 and c HCHO = 18 µg/m 3 [22]). The IEQ index with IAQ quality values was calculated using two variants-the first conventional and the second borrowed from the IAQI scale [30] ( Table 8a).
Concentrations c and SD(c) are in µg/m 3 ; PD* H , PD* L , and SD(PD*) are in %, and standard deviations of the HCHO concentration of 12% was adopted on the basis of reports from IAQ research conducted as part of BREEAM in 2016 [22]. The assumptions were that c H and PD* H are the coordinates of the upper break point (i.e., high break point) of the converted scale PD* = f(c), and c L and PD* L are the coordinates of the lower break point (i.e., low break point) of the converted scale PD* = f(c). Standard deviations were assumed for c meas , c TVOC , and c HCHO as well as for PD* = ±12%, as this is half of the transformed segment of the scale, which according to Table 7 covers a range of 25% of PD*, with one perceived category of air quality, e.g.,"no risk" or "moderate". Therefore, the maximum standard deviation was 12.5% and, according to the literature on AQI and IAQI values, it should be rounded to a total value. 2) SD vote (PD(SI i )) from the ±u overall (IEQ) equation was the standard deviation of a probability distribution of an each. (SI i ) vote and was calculated primary using the PD(SI i ) equation calibration curve [34].

Discussion of the IAQ index Theoretical Model
For years, the authors, as accredited laboratory personnel, have conducted IAQ pollution tests in indoor environments for various applications. Based on our experience, it was concluded that the general approach to assessing combined IAQ has not yet been systematized and that there is a global tendency to assess individual IAQ parameters separately or to group them without a justified aggregation method. This is not a good situation from the point of view of building users' needs. This, in our opinion, may lead to incorrect IAQ interpretations in specific building situations. In the context of analyzing the problem in this paper, authors presented a summary of the state-of-the-art methods and also provided a new approach for solving some of these problems. As presented, it is possible to create a IAQ index aggregating the results of indoor air analyses, taking into account various representative pollutants. Three levels of comprehensive air quality assessments (with three, five or seven subcomponents), depending on the application of the assessment, were proposed, together with step-by-step procedures. This may be practical, as shown in the evaluation of a case study on a building. We originally selected the main IAQ subcomponent equations and user satisfaction dependences, %PD = f(c j ), and provided them all in one place (Table 3). We then proposed and justified the weighting schemes for the IAQ total equation. In most of the studies in the literature, the weighting schemes used for IEQ or IAQ assessments are not physically justified or explained. There are known methods of weighting sub-indices, but the problem that was solved in this paper was an effective system for weight adjustments. For the construction of the combined model IAQ index , with a weighting scheme useful for aggregating sub-indices, we proposed the model scheme presented in Figure 5. According to the results, the advantage of the complex model IAQ index , in which the input quantities always constitute concentrations of given pollutants, is the ability to use these concentrations to calculate excess pollution concentrations from Equation (10) and generate weighting schemes W 1, . . . n for all three models by adjusting the weights based on the concentration values of excess air pollutants to a value ≤1.0 for each IAQ index model. The ∆c j values determine the masses of pollutants that must be removed by ventilation to eliminate the target pollutant effect. They can be determined as differences between the current concentrations of pollutants and the concentration of pollutants at the reference or standard level (e.g., c ELV or c LCI ), and in the case of VOC odorous , the odor threshold c th. The presented approach may allow planning of air quality for the building.
As discussed, it is important to identify those VOCs with comfort, health, and impacts and focus on the IAQ sub-model choice aspect, briefly defined as the strategy "VOCs-Total vs. Target: Comfort, Irritancy, Odor, and Health Impact". The model from Figure 5 has uniform inputs, i.e., concentration levels c j and two outputs: (1) weighed (adjusted) and (2) sensory equations, PD* = f (c TVOC ), constituting the IAQ submodel equations. These second outputs of submodels (PD* values) are coefficients of satisfaction from the comfort sensation or lack of "health risk". These are the terms of the equation describing "combined IAQ", which meets the requirements of the abovementioned strategy of selecting IAQ sub-models related to IAQ components that have the most impact on the resulting IEQ perception.
Models for subcomponents of IAQ not perceived by humans but influencing health are recommended to be used from the index set in the AQI system [26 -28], which was adapted by the authors to assess the quality of indoor air based on, and in accordance with, the concepts of the air quality assessment system used globally by the American EPA. In the context of the subcomponent of the TVOC concentration in Figure 7, the authors provided the relationships of PD* = f(c TVOC ) based on Jokl research [49] and resulting from the IAQI scale [30] as converted by the authors. The relationship between PD* and TVOC concentration in both approaches is strongly correlated, as shown in Figure 9.
of the TVOC concentration in Figure 8, the authors provided the relationships of PD* = f(cTVOC) based on Jokl research [49] and resulting from the IAQI scale [30] as converted by the authors. The relationship between PD* and TVOC concentration in both approaches is strongly correlated, as shown in Figure 9. The curves obtained from the conversion confirm Jokl's predictions provided in [49] and previously accepted idea. However, the final confirmation of these curves will be done experimentally as panel tests, as planned by the authors for the near future.

Discussion of Results for the Case Study on a Building
The experimental study was performed in a BREEAM certified building, and included the determination of formaldehyde concentration, CO2, and VOCs in the indoor air. The example calculation of the combined ∑IAQ model for three basic pollutants, as components of the IEQindex model, is presented in two variants but the calculated PD*TVOC values obtained with both calculation methodswere very similar. The first variant of ∑IAQ calculation used %PD = f(cj) curves in % sub-indices of IAQ for three pollutants, and the differences in approach in Table 5 and Table 8 meant combining them into one IEQ submodel: the ∑IAQ model intended for the IEQ calculation. The second variant used submodels of IAQ for TVOC and HCHO pollutants based on the IAQI system [30]which were then converted into percentages of persons dissatisfied (PD* in %).
The first conclusion is that CO2 concentration cannot be used separately for the IAQindex assessment, especially at the pre-occupancy stage ( Table 5). The building was polluted with VOC emissions and HCHO from the construction products directly after finishing works were completed.
The authors confirmed that all three pollutions should be a simultaneously integrated part of the IAQ model, because the importance of TVOC is much greater, representing the main source of pollution-the construction and finishing materials.
According to the results, we recognized two variants of the combined ∑IAQindex calculation. For the first variant [22], the combined ∑IAQ index of satisfied users was 69.1% and, for the second variant (new approach with converted AQI index), the ∑IAQ index was 70.0% satisfied. The results ofthe IEQindex(1) (for Variant 1) were within the interval of combined overall uncertainty,16.24%, Figure 9. PD TVOC based on conversion from IAQI to the PD* scale and study of Weber-Fechner theory.
The curves obtained from the conversion confirm Jokl's predictions provided in [49] and previously accepted idea. However, the final confirmation of these curves will be done experimentally as panel tests, as planned by the authors for the near future.

Discussion of Results for the Case Study on a Building
The experimental study was performed in a BREEAM certified building, and included the determination of formaldehyde concentration, CO 2 , and VOCs in the indoor air. The example calculation of the combined IAQ model for three basic pollutants, as components of the IEQ index model, is presented in two variants but the calculated PD* TVOC values obtained with both calculation methodswere very similar. The first variant of IAQ calculation used %PD = f(c j ) curves in % sub-indices of IAQ for three pollutants, and the differences in approach in Tables 5 and 8a meant combining them into one IEQ submodel: the IAQ model intended for the IEQ calculation. The second variant used submodels of IAQ for TVOC and HCHO pollutants based on the IAQI system [30] which were then converted into percentages of persons dissatisfied (PD* in %).
tests carried out immediately after finishing works gave results that significantly exceeded the BREEAM limits for TVOC at 300 µg/m 3 (twice as high). It should be expected that an acceptable level should be reached after a minimum of one month from the completion of the work.
For correctness of the obtained calculations, the authors are conducting a model credibility analysis that will be provided in the next article-Indoor Air Quality Model Part II: The Combined Model IAQ index Reliability Analysis. The model uncertainty estimate may be compromised because the model reproduces the discomfort level associated with the dominant component.