Next Article in Journal
Comparative Analysis of Neural Network Models for Predicting Ammonia Concentrations in a Mechanically Ventilated Sow Gestation Facility in Korea
Next Article in Special Issue
The Role of Industrial Structure Upgrading in Moderating the Impact of Environmental Regulation on Air Pollution: Evidence from China
Previous Article in Journal
Aerosol Properties and Their Influences on Marine Boundary Layer Cloud Condensation Nuclei over the Southern Ocean
Previous Article in Special Issue
In Vitro Bioaccessibility and Health Risk of Heavy Metals from PM2.5/PM10 in Arid Areas—Hotan City, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Extending Multi-Pathway Human Health Risk Assessment from Regional to Country-Wide—A Case Study on Kuwait

by
Mohammad Munshed
1,2,*,
Jesse Van Griensven Thé
1,2,
Roydon Fraser
1,
Bryan Matthews
2 and
Ashraf Ramadan
3
1
Department of Mechanical and Mechatronics Engineering, University of Waterloo, Waterloo, ON N2L 3G1, Canada
2
Lakes Environmental Software, Waterloo, ON N2L 3L3, Canada
3
Environmental Pollution and Climate Program, Environment & Life Sciences Research Center, Kuwait Institute for Scientific Research, P.O. Box 24885, Safat 13109, Kuwait
*
Author to whom correspondence should be addressed.
Atmosphere 2023, 14(8), 1247; https://doi.org/10.3390/atmos14081247
Submission received: 21 June 2023 / Revised: 31 July 2023 / Accepted: 3 August 2023 / Published: 5 August 2023
(This article belongs to the Special Issue Air Pollution Exposure and Health Impact Assessment)

Abstract

:
Air pollution has emerged as a pressing global issue in recent decades. While criteria pollutants and greenhouse gases contribute to the problem, this article explicitly addresses hazardous air pollutants (HAPs). This work estimates the country-wide cumulative human health impacts from exposure to HAPs. Kuwait is used as the case study due to data availability and non-fragmentation of data. At present, the evaluation of multi-pathway human health risks arising from exposure to HAPs is incomplete, as indirect pathways have not been considered. Furthermore, only a few HAPs, such as benzene, have established ambient air quality standards specifically intended to safeguard human health, leaving many HAPs unregulated. This study considers several pathways (both direct and indirect) and various environmental media (air, water, plants, soil, and animal tissue). The findings indicate that cumulative health risks in the coastal air quality zone are within acceptable limits but are notably higher when compared to the other air quality zones. For cancer risks, only the Ahmadi Hospital, with a cancer risk of 1.09 × 10−5 for the resident adult exposure scenario, slightly exceeds the acceptable risk level of 1 × 10−5. The proposed methodology integrates the results from a country-wide emissions inventory composed of different air quality zones, air dispersion and deposition modeling, multi-pathway transport-and-fate analysis, exposure quantification, and health risk and hazard characterization. It also extends and adapts EPA methodologies initially designed for hazardous waste combustion facilities to additional emission sources and provides a case study for a region seldom subjected to such human health risk assessments.

1. Introduction

Kuwait, a country in the Middle East, has been selected as the country of study, rather than, for example, the United States, due to the unfragmented nature of the data and the smaller size of the country. Kuwait is situated in the northeastern region of the Arabian Peninsula, bordered by Saudi Arabia to the south and Iraq to the north. It spans an area of approximately 17,818 square kilometers, slightly smaller than New Jersey in land area, and is located at the geographical coordinates of 29.3759° N, 47.9774° E [1].
Kuwait is a country that heavily relies on its abundant oil and gas resources, which have led to the development of various industries such as petrochemical production, refining, power generation, crude steel production using electric arc furnaces, and desalination using natural gas and heavy oil.
As a result of various oil, gas, and industrial activities, numerous pollutants are emitted from different emission source types. These include point sources, such as stacks associated with compressors, boilers, turbines, incinerators, flares, and exhaust vents related to glycol dehydration processes. In addition, area and volume sources represent other categories of emission sources. Area sources encompass fugitive emissions characterized by non-buoyant releases with negligible vertical extent, such as equipment leaks or landfill emissions. In contrast, volume sources involve non-point emissions with an initial vertical extent, including prescribed burns and marine vessel emissions.
Additionally, on-road and non-road mobile sources contribute to pollutant emissions. On-road mobile sources comprise motorcycles, passenger cars, light-duty trucks, buses, and heavy-duty trucks. Non-road mobile sources include agricultural equipment, aircraft jets, forklifts, and construction equipment, such as graders and backhoes.
Due to various industrial activities, oil and gas operations, traffic congestion, and frequent dust storms in Kuwait throughout the year [2,3], these factors lead to poor ambient air quality. Recent research indicates a decline in air quality across both urban and rural areas of Kuwait [4,5,6], which has subsequently been associated with increased mortality and morbidity rates [7,8,9].
Kuwait has three distinct air quality zones (AQZs): coastal, inland, and production, which are further divided into eleven subzones [10]. Each air quality zone and subzone presents a unique combination of emission source types, pollutants, and exposure scenarios that must be considered. The proposed approach employs the validated U.S. Environmental Protection Agency (EPA) Human Health Risk Assessment Protocol (HHRAP) risk assessment methodology [11] established in the US and Canada.

1.1. Literature Review, Objectives, and Additional Contributions

1.1.1. Literature Review

Exposure to hazardous air pollutants (HAPs), also known as air toxics, is associated with various health impacts [12], leading to short-term and long-term health problems. The specific health effects experienced by the population depend on the type, concentration, and mixture of air toxics, as well as the duration and frequency of exposure. For instance, respiratory issues may arise from exposure to methyl mercaptan, which can aggravate respiratory conditions [13], resulting in increased coughing, wheezing, and shortness of breath. Moreover, cardiovascular diseases, including atherosclerosis, have been linked to exposure to polycyclic aromatic hydrocarbons (PAHs) such as benzo(a)pyrene [14,15,16,17].
Neurological disorders may also be associated with exposure to air toxics, as pollutants like mercury and lead can cause cognitive impairments and memory loss [18], and in severe cases, conditions such as convulsions and coma [19,20,21,22]. Additionally, reproductive issues, including fertility problems [23], hormonal imbalances [24,25], and congenital disabilities [26], have been associated with exposure to certain air toxics, such as dioxins and some PAHs.
Lastly, long-term exposure to HAPs, such as benzene, formaldehyde, and 1,3-butadiene, can increase the risk of developing several types of cancer, including lung cancer [27,28,29,30,31,32] and lymphohematopoietic cancers [33]. It is important to note that vulnerable populations, including children, the elderly, and those with pre-existing health conditions, may be more susceptible to the adverse effects of air toxics. The HAPs summarized in Table 1 represent only a small fraction of air toxics released into the ambient environment from various industrial processes. These chemicals have been selected for their potency and the extensive availability of both cancer and non-cancer assessments in the literature.
To date, no country-wide studies have addressed the cumulative human health risks associated with exposure to hazardous air pollutants. Furthermore, no studies have yet considered multiple exposure (i.e., direct inhalation and scenario-relevant indirect pathways) in combination with numerous different sources (e.g., oil and gas operations; wastewater treatment plants; etc.). See Table 2 for a summary of the types of human health risk assessments that have and have not been conducted.

1.1.2. Objectives

The specific objectives of this research, using Kuwait as a case study, are as follows:
  • To implement a country-wide cumulative human health risk assessment incorporating multi-source and multi-pathway exposures.
  • To compare country-wide human health risk variability by region, for example, coastal versus inland.
  • To identify the chemical risk driver (defined by the authors as the dominant chemical of potential concern) for a particular sensitive receptor.
  • To identify if the direct or indirect pathway is the dominant pathway of risk for a particular sensitive receptor.
Though beyond the scope of this paper, these objectives support the long-term goals of targeted risk reduction and targeted mitigation resource allocations.

1.1.3. Additional Contributions

  • Extending and adapting EPA methodologies initially designed for hazardous waste combustion facilities to additional emission sources, including wastewater treatment plants, glycol dehydration units, but not including fugitives and mobile sources.
  • Providing a case study for a region of the world seldom subjected to such human health risk assessments.
Figure 1 presents the methodology framework used in this study.

2. Materials and Methods

The methodology consists of the following sequential steps:
  • Develop an emissions inventory (or augment an existing inventory) for each air quality zone.
  • Conduct air dispersion modeling using current US regulatory air dispersion models, such as AERMOD or similar, to estimate pollutant concentrations in the air and their deposition rates in various exposure media.
  • Estimate the concentrations of pollutants at the point-of-contact for receptor populations by conducting an environmental transport-and-fate analysis.
  • Identify realistic exposure scenarios to estimate the types and magnitudes of human exposure to pollutants.
  • Assess the levels, frequencies, and durations of contact between humans and pollutants.
  • Calculate multi-pathway and cumulative cancer risks and non-cancer hazards for each air quality zone (risk characterization).
  • Examine the contributing factors and underlying drivers of unacceptable risks (risk drivers).

2.1. Emissions Inventory (Step 1 in the Proposed Methodology)

An emissions inventory plays a vital role in managing air quality by enabling decisionmakers to:
  • Measure the contribution of each sector’s emissions.
  • Analyze emission trends both retrospectively and prospectively.
  • Provide guidance and assistance to policymakers and industry in regulating emissions and establishing achievable targets.
  • Supply the necessary data for accurate assessment of current or future emissions using air quality modeling tools, with an accurate emissions inventory being a prerequisite for air quality modeling.
  • Identify locations for monitoring hazardous air pollutant hotspots.
In this study, emissions were calculated using the methods and algorithms summarized in Table 3 below.
Table 3 lists methods chosen for their applicability, accuracy in emissions estimations, and broad acceptance within the scientific, regulatory, and industrial community. These validated methods were developed by reputable organizations such as the U.S. Environmental Protection Agency (EPA) and American Petroleum Institute (API), based on research, empirical data, and engineering principles.
For each air quality zone, unique identification codes were assigned to each emission unit and to each process based on the fuel type consumed. These codes facilitated the tracking of results and thorough examination of emissions from individual units and processes across facilities and sectors. This coding procedure was instrumental in the study’s quality assurance process, enabling the review of assumptions, calculations, procedures, and the evaluation of the quality and representativeness of the inventory data [51].
It is important to note that emissions calculated using AP-42 emission factors were based on normal, steady-state conditions and did not take into account worst-case scenario conditions, such as during process start-up and shut-down, as well as process upsets. While these abnormal operating conditions have a relatively short time span, resulting emissions can easily exceed relevant emission limits and potentially cause ground-level concentrations that exceed relevant ambient air quality standards (e.g., National Ambient Air Quality Standards). According to Obaid et al. [52], various industry sectors, including power and/or heat generation, nuclear power generation, petrochemical production, energy-from-waste generation, and sulfuric acid production, exhibit different characteristics for start-up and shut-down emissions compared to steady-state conditions, and these characteristics differ for different source types and pollutants. The same authors note that attempts to characterize emissions during these conditions face numerous challenges, including the unavailability of manufacturers’ design data and the complex relationship between different dynamic process functions and timing of events. In addition to these challenges, other uncertainties associated with the calculation of criteria pollutants and HAPs such as benzene, toluene, ethylbenzene, and xylene (BTEX) emissions include data gaps, data quality, and inconsistencies across different power plants and sometimes even different emission units within the same power plant.
The general formula for estimating emissions is [46]:
E m i s s i o n s = A × E F × ( 1 E R 100 )
where A represents the activity rate, EF is the emission factor, and ER is the total percentage (%) of emission reduction efficiency.
It is important to note that the availability and completeness of an emissions inventory are critical factors in evaluating the risks to human health associated with exposure to hazardous air pollutants.

2.2. Atmospheric Dispersion and Deposition Modeling (Step 2 in the Proposed Methodology)

Air dispersion modeling is defined as the mathematical representation of pollutant transport in the atmosphere, along with the quantification of deposition rates in various exposure media [12]. Air dispersion modeling is performed to understand the location and magnitude and to predict the spatial distribution of pollutant concentrations. It is also used to compare modeled and monitored air concentrations to evaluate the air dispersion model performance and representativeness of the emissions inventory. Moreover, air dispersion and deposition modeling is a prerequisite for conducting a multi-pathway human health risk assessment, as such modeling is used to estimate pollutant concentrations in the air, which are crucial for calculating inhalation cancer risks, and to estimate deposition rates, which are critical for assessing indirect exposure pathways.
In this study, air dispersion and deposition modeling was conducted using the current EPA-preferred regulatory model, American Meteorological Society/United States Environmental Protection Agency Regulatory Model (AERMOD) [53], employing version 22112.
AERMOD, an extensively validated steady-state Gaussian plume model, integrates air dispersion principles based on planetary boundary layer turbulence structure and scaling concepts. It addresses surface and elevated emission sources and accommodates simple and complex terrain. Since 2005, the EPA has approved AERMOD as its regulatory model, and it is widely used in various applications, including permitting, air quality assessments, and human health risk evaluations [53].
Designed as a near-field model for distances up to 50 km from the source, AERMOD exhibits robust performance across diverse terrains and under a wide array of meteorological conditions. AERMOD calculates concentrations using mathematical principles, such as plume rise, atmospheric dispersion of buoyant (or neutrally buoyant) effluent, turbulence theory, and Gaussian distribution [54].
The AERMOD modeling system does not utilize stability classes, which were used in the previous Industrial Source Complex Short Term 3 (ISCST3) model but employs a more advanced turbulence scaling parameterization scheme. The model predominantly relies on surface weather observations to generate vertical profiles of variables, including temperature, temperature gradient, wind speed, wind direction, and turbulent velocities within the atmospheric boundary layer. The EPA’s AERMET model is the meteorological preprocessor responsible for producing AERMOD meteorological data files.
While validation of the AERMOD model is a desirable aspect of air dispersion modeling to ensure result credibility for a specific study area, such an evaluation was not conducted due to an absence of available measured air pollutant concentrations.

2.2.1. Unitized Emission Approach

The unitized emission approach in air dispersion modeling is a highly efficient method for evaluating the impact of multiple pollutants emitted from a single source [55]. This technique relies on the principle that air concentration levels and deposition fluxes of pollutants are linearly related to the source’s emission rate, meaning that a tenfold increase in emissions would result in a tenfold increase in receptor concentrations. By employing a standardized emission rate of 1 g per second (g/s) for the source, this approach streamlines the modeling process and significantly reduces computational time.
Using the unitized emission rate, a single air dispersion modeling run can provide concentration and deposition data for all emitted pollutants. This is achieved by multiplying the unit emission results by the actual emission rate of each individual pollutant. Consequently, the need for running the model separately for each pollutant is eliminated, enhancing the method’s efficiency.
The output of this approach includes unitized annual average and maximum one-hour average air concentrations and deposition rates for the vapor, particle, and particle-bound phases. These unitized values can be further used as inputs for risk modeling, wherein media-specific concentrations and deposition of each contaminant are calculated by scaling the unitized concentrations by the appropriate pollutant emission rate.
The unitized emission approach is particularly beneficial for modeling a vast number of pollutants, as it minimizes calculation imprecision due to low emission values and simplifies the process by avoiding multiple runs for each modeled pollutant.

2.2.2. Meteorological Data

The Weather Research and Forecasting (WRF) model [12,56] was used to generate two AERMOD meteorological data files (surface and profile files) for each of the three air quality zones in Kuwait (i.e., coastal, inland, and production). To ensure adequate coverage across the country, three pseudo-meteorological stations (pseudo-stations) were strategically positioned. Pseudo-stations define the locations for which WRF data are generated, with meteorological parameters representative of the selected location.

2.2.3. Air Dispersion Model General Options

The modeling options and settings used in AERMOD for the analysis are provided in Table 4.

2.2.4. Receptors (Calculation Points)

Receptors are specific locations with designated coordinates where pollutant concentrations are calculated. These receptors can represent human populations, ecosystems, or other areas sensitive to air pollution. AERMOD uses meteorological data, emission source parameterizations, and terrain datasets, such as the National Aeronautics and Space Administration (NASA) Shuttle Radar Topography Mission (SRTM) digital files, to predict the dispersion and concentration of pollutants at these receptors, thereby helping to evaluate air quality and assess the potential impact on health or the environment.
In each of the three air quality zones, two uniform Cartesian receptor grids were established centered around the facility of interest. These grids were configured with both fine (100 m spaced grid from the centroid of the emission sources out to a radius of 3 km) and coarse (500 m spaced grid extending from 3 km to 10 km) settings to cover all areas where emissions may have a significant impact. The selection of receptor locations typically depends on regulatory requirements and the specific objectives of the air quality assessment, as different projects may have varying goals and criteria for receptor placement. Additionally, sensitive receptors, defined as locations where people are more susceptible to adverse effects from exposure to hazardous air pollutants (such as schools, hospitals, residential areas, daycare facilities, care facilities, and places of worship), were included in the project. Further information about these sensitive receptors can be found in the case study section of this paper.

2.3. Transport, Fate, Exposure, and Risk Characterization (TFER)

2.3.1. Transport-and-Fate Modeling (Step 3 in the Proposed Methodology)

This step examines the behavior and movement of hazardous air pollutants (HAPs) following their release into the atmosphere, which is critical for determining their potential impact on human health and the environment. Transport refers to the physical movement of HAPs as they travel through the environment, while fate pertains to their ultimate destination within environmental media (e.g., soil, food, water). This analysis aims to calculate the concentrations of HAPs in air, soil, produce, milk, meats, eggs, fish, and drinking water. The focus is on ensuring accuracy in assessing the presence of HAPs in different environmental media, which is vital for understanding their potential impact on human health and the environment. The media concentrations were estimated using over forty mathematical equations [11], including a system of non-linear equations, such as the cumulative soil concentration for carcinogens and non-carcinogens, and the hazardous air pollutants’ loss constant due to runoff. The modeling approach accounted for numerous factors, such as meteorological conditions, physicochemical properties of the pollutants, and environmental characteristics of the study area, including annual evapotranspiration, annual irrigation, annual precipitation, and annual runoff. It is important to acknowledge that the modeling process may involve certain assumptions, limitations, or uncertainties that should be considered when interpreting the results. To facilitate transport-and-fate modeling, the computer application IRAP-h View [57] was used. An example of the output of the transport-and-fate modeling is listed below:
  • Aboveground exposed produce concentration due to direct deposition: [mg/kg]

2.3.2. Exposure Quantification (Steps 4 and 5 in the Proposed Methodology)

Exposure evaluation is the process of assessing the frequency, duration, and extent of human exposure to air pollutants, focusing on hazardous air pollutants in this study. Exposure can be estimated using one of three approaches [58]: direct measurement, estimation, or exposure reconstruction. The proposed methodology utilizes the estimation approach.
In addition to the exposure scenarios and pathways, it is essential to consider the demographic and lifestyle factors that may influence exposure. These factors include age, gender, occupation, socioeconomic status, dietary habits, and cultural practices. Understanding these factors can help identify susceptible populations and inform targeted risk reduction strategies.
This section describes the equations used for estimating exposure, as well as exposure scenarios and locations and scenario-relevant exposure pathways. It is important to note that separate equations are used for inhalation and ingestion exposure to account for differences in exposure pathways, mechanisms, and potential health effects.
Exposure was estimated using Equations (2) and (3) [58]:
E x p o s u r e   ( i n h a l a t i o n ) = ADD   ( inhalation ) = C a i r × I n h R × E T × E F × E D BW × A T
where ADD represents the average daily dose (mg/kg-day), C a i r is the concentration of hazardous air pollutant(s) in the air (mg/m3), I n h R is the inhalation rate (m3/h), E T is the exposure time (h/day), E F is the exposure frequency (days/year), E D is the exposure duration (years), BW is the body weight (kg), and A T is the averaging time (days).
E x p o s u r e   ( i n g e s t i o n ) = ADD   ( ingestion ) = C m e d i u m × I n g R × E F × E D BW × A T
where ADD represents the average daily potential dose (mg/kg-day), C m e d i u m is the concentration of hazardous air pollutant(s) in the medium (e.g., mg/L, mg/g), I n g R is the ingestion rate (e.g., L/day, g/day), E F is the exposure frequency (days/year), E D is the exposure duration (years), BW is the body weight (kg), and A T is the averaging time (days).
Exposure scenarios define the combination of exposure pathways and exposure parameters applied to complete the risk and hazard calculations described in the risk characterization section. For the purposes of this study, exposure scenarios include the resident adult and resident child scenarios [11], which differ in factors such as exposure duration, frequency, and body weight. The included exposure pathways are summarized in Table 5.
Exposure scenario locations refer to the physical places (i.e., receptors) within the study area chosen for evaluating residential adult and child exposure scenarios. To ensure thorough evaluation of all potential residential exposure locations throughout Kuwait, air dispersion modeling results were generated using two uniform Cartesian receptor grids and sensitive receptors. For more information on receptor grids and sensitive receptors, please refer to Section 2.2.4. This approach aligns with established best practices in the field [11]. Figure 2 illustrates the relationship between human exposure to HAPs through multiple pathways. HAPs exist in various environmental media, such as air, water, and soil, and can enter the human body via different exposure routes, including ingestion and inhalation. HAPs are deposited into environmental media, such as water and soil, through atmospheric deposition.

2.3.3. Quantitative Estimation of Cancer Risk and Non-cancer Hazard (Step 6 in the Proposed Methodology)

Risk characterization is a crucial component in human health risk assessments. It synthesizes information from preceding assessment steps, such as transport-and-fate modeling (Section 2.3.1) and exposure quantification (Section 2.3.2), to determine the presence or absence of risk and quantify expected outcomes. Furthermore, risk characterization identifies key factors contributing to risk, enabling regulators, engineers, analysts, and stakeholders to make informed decisions about appropriate actions for risk reduction. It is essential to recognize the uncertainties and limitations inherent in the risk characterization process, which may stem from factors like incomplete emissions data, variability in exposed populations, and assumptions made during the assessment to account for temporal variability in emissions. Acknowledging these uncertainties provides context for the results and informs future research or decision making.
Cancer risk is the incremental probability of an individual developing cancer over a lifetime due to exposure to a specific carcinogenic hazardous air pollutant or a mixture of carcinogenic hazardous air pollutants [11]. The predicted risk represents the additional risk of cancer from the exposure being analyzed, which is supplementary to other risks resulting from numerous factors (e.g., smoking). This risk is expressed as a probability ranging from zero (no risk) to one (certainty of developing cancer), without any units or dimensions. Cancer risk is represented as a numerical value, such as “1 in 100,000” or “1 × 10−5”, indicating that one individual in a population of 100,000 is expected to develop cancer due to exposure to the carcinogenic hazardous air pollutant over a lifetime. The calculation of cancer risk considers several factors, such as the concentration of the carcinogenic hazardous air pollutant in the environment, exposure duration, frequency, and exposure route (e.g., ingestion or inhalation), as well as the hazardous air pollutant’s cancer potency factor or slope factor.
For inhalation exposures, chronic cancer risks for individual pollutants are estimated by multiplying the long-term exposure concentration estimate by the corresponding pollutant-specific inhalation unit risk. This calculation determines the potential incremental cancer risk, as illustrated in Equation (4):
C a n c e r   R i s k   ( i n h a l a t i o n ) = I n h a l a t i o n   E x p o s u r e   C o n c e n t r a t i o n × I n h a l a t i o n   U n i t   R i s k
For each ingestion exposure pathway being evaluated, chronic cancer risks for individual pollutants are estimated by multiplying the lifetime average daily dose estimate for each ingestion exposure pathway by the pollutant-specific cancer slope factor. This calculation estimates potential incremental cancer risk, as illustrated in Equation (5):
C a n c e r   R i s k   ( i n g e s t i o n ) = L i f e t i m e   A v e r a g e   D a i l y   D o s e × C a n c e r   S l o p e   F a c t o r
Non-cancer hazard refers to potential adverse non-cancer health effects, such as respiratory issues, neurological complications, and neurotoxicity, resulting from exposure to a hazardous air pollutant or a mixture of hazardous air pollutants. Exposure duration plays a significant role in hazard classification, with two primary categories: acute hazards associated with short-term exposure and chronic hazards linked to long-term exposure. Evaluating non-cancer hazards involves comparing the estimated exposure level of a substance to its reference dose or reference concentration, which represent levels deemed safe without appreciable risk of adverse health effects throughout a lifetime of exposure [59]. The non-cancer hazard is calculated using Equation (6):
H a z a r d   Q u o t i e n t = E s t i m a t e d   E x p o s u r e   L e v e l R e f e r e n c e   C o n c e n t r a t i o n   ( o r   R e f e r e n c e   D o s e )
The hazard index, a conservative measure, is used in human health risk assessments involving multiple hazardous air pollutants or exposure pathways. This index represents the sum of the individual hazard quotients for each hazardous air pollutant or pathway.
In this research, we define cumulative risk (cancer or non-cancer) as the risk resulting from simultaneous exposure to multiple hazardous air pollutants from multiple sources through one or more pathways or exposure routes, such as inhalation or ingestion. The total cumulative cancer risk was computed using Equation (7):
R t o t a l =   ( R i , j ) +   ( R i , i n h a l a t i o n )
where R t o t a l represents the cumulative cancer risk estimate, considering both ingestion and inhalation pathways, R i , j represents the individual risk estimate for hazardous air pollutant i and ingestion pathway j, R i , i n h a l a t i o n represents the individual risk estimate for hazardous air pollutant i through the inhalation pathway, i denotes the different hazardous air pollutants (e.g., hazardous air pollutant 1, hazardous air pollutant 2, hazardous air pollutant 3, etc.), and j denotes the different ingestion pathways (e.g., vegetable ingestion, fish ingestion, egg ingestion, beef ingestion, etc.). On the other hand, the total cumulative non-cancer risk (hazard index) was estimated using Equation (8):
H a z a r d   I n d e x t o t a l =   ( H Q i , j ) +   ( H Q i , i n h a l a t i o n )
where H a z a r d   I n d e x t o t a l represents the total hazard index, considering both ingestion and inhalation pathways, H Q i , j represents the hazard quotient for hazardous air pollutant i and ingestion pathway j, H Q i , i n h a l a t i o n represents the hazard quotient for hazardous air pollutant i through the inhalation pathway, i denotes the different hazardous air pollutants (e.g., hazardous air pollutant 1, hazardous air pollutant 2, hazardous air pollutant 3, etc.), and j denotes the different ingestion pathways (e.g., vegetable ingestion, fish ingestion, egg ingestion, beef ingestion, etc.).
Lastly, the acute inhalation hazard quotient was calculated using Equation (9):
A c u t e   I n h a l a t i o n   H a z a r d   Q u o t i e n t = a c u t e   a i r   c o n c e n t r a t i o n a c u t e   i n h a l a t i o n   e x p o s u r e   c r i t e r i a
Regulatory agencies establish target risk levels to assess potential health risks from hazardous air pollutant exposure and ensure public safety. Carcinogenic risk levels typically range from 1 × 10−4 to 1 × 10−6, indicating a negligible increase in cancer risk due to exposure. For non-cancer health effects, a hazard quotient or hazard index of 1 or lower indicates no expected adverse effects, while values above 1 suggest potential risks.
In this study, based on the U.S. EPA Region 6 HHRAP Addendum [60], we adopted narrower target risk levels for hazardous air pollutant exposure to account for conservatism and background concentration exposure within the study area. Consequently, cancer risk levels should not exceed 1 × 10−5, and the hazard quotient/index should not exceed 0.25.
It is essential to note that target risk levels are not definitive indicators of adverse effects, as risk estimations inherently have uncertainties. These uncertainties can arise from exposure variability, limited toxicological data, temporal emissions’ variability, and assessment assumptions. If calculated risks are within target values, no further investigation is needed, and the evaluated conditions are not considered to present unacceptable risks. However, when target values are exceeded, the methodology entails further review by risk managers, risk analysts, or engineers. This review process involves examining the scientific basis and uncertainties associated with the calculation, including toxicity factors and uncertainty factors addressing study limitations and data quality issues, as well as investigating the contributing factors and underlying drivers of unacceptable risks.

2.3.4. Risk Driver Analysis: Forensic Analysis of Human Health Risks (Step 7 in the Proposed Methodology)

Risk driver analysis is a systematic, project-specific process that addresses the fundamental question, “Where does the risk come from?”. It is important to identify and understand the primary factors contributing to unacceptable cancer and non-cancer risks associated with hazardous air pollutant exposure in a specific project setting. This stage involves an examination of various exposure elements, including scenarios, such as farmers, residents, local fish consumers, construction workers, and outdoor workers, as well as pathways like air, water, soil, plants, and animals.
The analysis assesses the significance of specific exposure pathways, such as inhalation, ingestion, and dermal contact, to pinpoint the most critical factors leading to elevated risk levels. A vital aspect of this process is tracing back to the sources or chemicals contributing to the increased risks. Additionally, the analysis investigates the potential for bioaccumulation of hazardous air pollutants in tissues or agricultural produce to enhance understanding of the risks associated with human health in the project area. For instance, methylmercury is a lipophilic substance that tends to bioaccumulate in the food chain [61,62].
The project-specific risk assessment is re-evaluated upon completing the analysis to ensure accuracy and guide development of targeted risk reduction strategies. These strategies may include emission controls, exposure reduction measures, or regulatory actions aimed at addressing the identified risk drivers and protecting public health within the project’s scope. Figure 3 illustrates a simplified example of a risk driver analysis. As shown in the example, the risk driver analysis employs a systematic and methodological approach to investigate the various components of human health risk assessment, identifying the elements contributing to elevated risk and tracing them back to the emission sources or hazardous air pollutants driving the risk.
By implementing this methodology, we effectively evaluated human health risks associated with exposure to hazardous air pollutants in Kuwait. The methodology encompasses steps such as developing emissions inventories for each air quality zone, conducting air dispersion modeling using AERMOD, estimating pollutant concentrations at the point of contact, identifying exposure scenarios, and assessing levels and frequencies of contact between humans and pollutants.
Through the process of calculating multi-pathway and cumulative risks, our assessment considered both direct and indirect exposure pathways. This allowed us to examine the contributing factors and underlying drivers of unacceptable risks in each air quality zone, ultimately providing a consistent and robust risk characterization upon which risk mitigation strategies can be developed to protect human health in Kuwait.

3. Results

3.1. Case Study

This section presents the application of the human health risk assessment methodology outlined above to a focused case study in Kuwait, using data from the year 2017. The analysis concentrates on hazardous air pollutants, emphasizing the impact of direct emissions from industrial sources across Kuwait’s three air quality zones. Although mobile sources and fugitive emissions are significant sources of localized hazardous air pollutants [12,55,63], they were intentionally excluded from this case study. This approach enables us to concentrate on the effects of industrial emissions, providing a clear baseline for the analysis. In addition, industrial sources are regulated which provides a legal framework upon which emissions reduction actions can be developed to address specific issues identified through the health risk analysis process.
Emission sources, representative of Kuwait’s distinctive pollution landscape, were selected to ensure an accurate depiction. Our selection employed a multi-criteria approach, factoring in not only the type of industry, geographical location within air quality zones, and the intensity of emissions but also their unique contribution to Kuwait’s pollution landscape. The unique contribution refers to how certain industries in Kuwait are characterized by specialized processes and materials, resulting in distinct emissions profiles and pollution effects, which differ from other regions. In addition, we further enriched our criteria by considering the proximity to populated areas, meteorological conditions, and land-use characteristics (inland versus shoreline), acknowledging their potential impact on the dispersion and intensity of pollutants. This multi-faceted selection process ensures that our study is comprehensive, representative, and aligned with the actual pollution landscape in Kuwait.
The results of our assessments for each air quality zone, to be presented later in this paper, provide an understanding of the prevailing air quality conditions and associated health risks in Kuwait. This case study effectively underscores the practical applicability and potential impact of our proposed methodology in human health risk assessment. A map of the study region is provided in Figure 4 for a clearer geographical context.
In our efforts to accurately represent Kuwait’s unique pollution landscape, we selected several facilities of interest across the three air quality zones, each characterized by different types of emission sources. Table 6 presents a summary of the number of modeled facilities in each air quality zone, along with the number of their corresponding emission sources. Additionally, Table 7 lists the hazardous air pollutants that were considered in our study.
For this study, we selected a small, representative subset of hazardous air pollutants known for their adverse impacts on human health and their common presence in industrial emissions. The selected pollutants, drawn from a larger emissions inventory, have been chosen due to their potency, their prevalence in industrial emissions, and the extensive availability of both cancer and non-cancer assessments in the literature. Furthermore, a subset of pollutants is all that is required in order to demonstrate that by using a normalized pollutant that the human health impacts from any number of hazardous air pollutants can be easily found, meaning only one calculation of the time-consuming air dispersion modeling is required. By limiting the pollutants analyzed, this paper does not provide a comprehensive picture of air pollutant health risks in Kuwait; such a comprehensive impact assessment is outside the focus and scope of this paper. Table 1 (refer to Section 1.1) details the adverse health effects associated with these selected pollutants.

3.2. Cumulative Risk Results

This subsection presents the cumulative health risks from all pathways evaluated within our human health risk assessment study. This includes the analysis of cancer, non-cancer (hazard index), and acute non-cancer health effects within Kuwait’s three distinctive air quality zones: inland, production, and coastal. For a detailed definition of cancer risk and non-cancer risk, please refer to Section 2.3.3 of this paper. In the tables presented below, we provide numerical data to illustrate the calculated health risks for each receptor, under specific exposure scenarios, in the inland, production, and coastal air quality zones. Each value represents a statistical estimate of the risk posed to a “receptor”, which refers to an individual, either an adult or a child, residing in the area. These risk estimates are calculated based on potential exposure to various air pollutants, considering factors such as pollutant concentrations, exposure duration, frequency, and the inherent toxicity of the pollutants.

3.2.1. Inland

Table 8, Table 9 and Table 10 provide a summary of the cumulative cancer risk, cumulative non-cancer risk, and cumulative acute non-cancer risk for the inland air quality zone. The sensitive receptors for this zone were identified based on the grid nodes, which exhibited the maximum air concentration values. This approach ensures the analysis focuses on the areas with the highest potential exposure.
The inland zone demonstrates a low health risk profile, with all sensitive receptors displaying risks below the established thresholds, as outlined in Section 2.3.3, for cancer, non-cancer, and acute non-cancer effects. The identified pollutants do not appear to pose significant health hazards for either adults or children within this air quality zone, based on our receptor selection.

3.2.2. Production

Table 11, Table 12 and Table 13 provide a summary of the cumulative cancer risk, cumulative non-cancer risk, and cumulative acute non-cancer risk for the production air quality zone. Consistent with the inland air quality zone, the sensitive receptors for this zone were determined based on the grid nodes with the calculated maximum air concentration values. Though these receptors share the same designations as those in the inland zone, they represent distinct locations within the production air quality zone. This method of selection maintains the focus of the analysis on areas with the highest potential exposure.
Our study findings for the production zone align closely with those from the inland zone. No significant health risks were identified, as all cumulative risks across different population segments remained below the specified thresholds. This finding indicates that the current levels of selected pollutants may not contribute significantly to health impacts within the production zone.

3.2.3. Coastal

Table 14, Table 15 and Table 16 present a summary of the cumulative cancer risk, cumulative non-cancer risk, and cumulative acute non-cancer risk within the coastal air quality zone. The sensitive receptors identified in this zone have been selected due to their strategic placement within densely populated areas and locations where populations are likely to be more susceptible to adverse effects from exposure to hazardous air pollutants. Such locations include critical community hubs such as hospitals, schools, resorts, and places of worship. The chosen receptors are of particular interest due to the intersecting factors of high potential pollutant concentrations and significant population presence, which jointly contribute to an increased potential for exposure.
The cumulative health risks in the coastal air quality zone are within acceptable limits, however, they are notably higher when compared to our findings for the inland and production zones. For cancer risks, only the Ahmadi Hospital, highlighted in boldface in Table 14, with a cancer risk of 1.09 × 10−5 for the resident adult exposure scenario, slightly exceeds the specified threshold of 1 × 10−5.
To provide a practical understanding of the cancer risk value above, consider it in terms of population-scale risk. If 100,000 adults were exposed to the same level of carcinogenic hazardous air pollutants—such as benzene, formaldehyde, toluene, and benzo(a)pyrene—from multiple emission sources over a lifetime, approximately 1.09 individuals in this group could be expected to develop cancer due to this exposure. However, this is a statistical estimate and does not guarantee that any particular individual in this group will indeed develop cancer from this exposure.
To further contextualize, the cancer risk at the Ahmadi Hospital sensitive receptor can be roughly compared to the risk of dying from an encounter with a hornet, wasp, or bee sting, which is 1 in 54,516 [64]. This comparison serves to convey that, while non-zero, the additional risk posed by exposure to these pollutants is relatively manageable when viewed alongside everyday risks. However, it should be remembered that these risk assessments rely on models and certain assumptions, including constant, long-term exposure. Actual risk may vary depending on factors such as genetic predispositions, lifestyle choices, and actual exposure patterns.
For the non-cancer and acute non-cancer risks, all sensitive receptors displayed values below the defined thresholds of 0.25. This result suggests that, while the level of pollutants in the coastal air quality zone is higher, it does not exceed the threshold that would pose significant health risks to the population. However, the risk is more pronounced in this zone compared to the inland and production zones.

4. Discussion

In this section, we analyze the key risk drivers, including specific pollutants, their sources, and the main exposure pathways. Subsequently, we consider the strengths and limitations of our methodology. Finally, we conclude by examining the variability and uncertainty associated with our analysis.

4.1. Risk Driver Analysis

The risk driver analysis from our assessment identified a singular case of cumulative cancer risk exceeding the target risk levels in the coastal air quality zone. This elevated risk was specifically noted at Ahmadi Hospital for the resident adult exposure scenario. Our systematic examination of exposure scenarios, emission sources, and exposure pathways indicated that the primary risk driver was benzo(a)pyrene. Moreover, our analysis suggested that the risk is primarily associated with indirect exposure pathways.
This finding implies that the origin of this pollutant may be localized emissions within the coastal zone, potentially originating from industrial activities or combustion processes. Due to confidentiality restrictions, specific information identifying the sources cannot be disclosed. However, it is worth noting that the likely culprits include industrial flares and other combustion sources. Consequently, the insights provided by our method serve as a key tool for policymakers and regulators. Our findings can guide the consideration and implementation of targeted risk reduction strategies, such as emission control measures at identified sources or various mitigation strategies. The ultimate goal of these strategies, informed by our research, would be to uphold public health standards within the project’s scope while reducing the cancer risk to acceptable levels.

4.2. Strengths and Limitations

In this section, we explore the strengths and limitations of the research methodology employed in this case study. A key strength is the integration of previously validated methods to quantify the country-wide cumulative multi-pathway human health impacts associated with exposure to hazardous air pollutants in Kuwait. While these methods have been previously validated, our research combines them in a unique manner to address the specific challenges of assessing cumulative, multi-pathway risks in a region with complex industrial activities. Additionally, the proposed methodology is flexible in incorporating more sophisticated multi-layer, multi-species air dispersion models, such as non-steady-state puff models like CALPUFF or SCIPUFF. The advantage of non-steady-state puff models over steady-state plume models lies in their capacity to capture the temporal and spatial variability in meteorological conditions, leading to a more accurate representation of pollutant dispersion patterns, thereby improving the precision and reliability of the exposure and risk estimates.
However, we must also acknowledge certain limitations of our approach. First, building downwash effects were not evaluated due to lack of available parameterization. Second, our study did not account for long-range transport of pollutants from areas outside of Kuwait.

4.3. Variability and Uncertainty

Variability and uncertainty in human health risk assessments originate from biological differences, external environmental factors, and the presence of gaps in both data availability and scientific knowledge. Variability refers to the differences or range of responses to environmental exposure among individuals or populations and can be attributed to factors such as genetics, age, lifestyle, and pre-existing health conditions. Environmental factors, such as the duration and route of exposure, as well as spatial–temporal variations in hazardous air pollutant concentrations, further contribute to this variability.
On the other hand, uncertainty refers to the lack or incompleteness of data regarding the input parameters for the methods employed to evaluate human health impacts. This form of uncertainty is inherent at every stage of the assessment process, from limitations in data availability to the extrapolation of results from animal studies to humans. Enhancing the reliability of risk assessments requires effectively addressing both variability and uncertainty. This can be accomplished by employing tools such as stochastic modeling and sensitivity analysis and by emphasizing the necessity for continuous research, monitoring, and data collection to reduce uncertainty over time.

5. Conclusions

This paper proposes an innovative approach to quantify country-wide cumulative risk (cancer, non-cancer, and acute) resulting from simultaneous exposure to multiple hazardous air pollutants from multiple sources through one or more pathways or exposure routes, such as direct inhalation or indirect ingestion. The authors also introduce a new method for identifying chemical and pathway risk drivers which provides a deeper understanding of potential health risks and provides a starting point for targeted risk reduction. The conclusions drawn from this work are as follows:
  • The overall health risk profile across the inland, production, and coastal air quality zones of Kuwait is low to moderate, as most risk values lie beneath the established risk thresholds. This suggests that the current levels of pollutants quantified in this case study do not likely pose significant health threats to the adult and child residential population.
  • The coastal air quality zone has a higher risk profile compared to the inland and production zones, particularly for cancer risks. However, these values are mostly within acceptable limits. An exception is the Ahmadi Hospital for the resident adult exposure scenario, where the cancer risk slightly exceeds the target level.
  • The risk driver analysis identified benzo(a)pyrene as the primary risk driver contributing to the elevated cancer risk at the Ahmadi Hospital in the coastal zone, calculated to be 1.09 × 10−5, with the likely sources being local industrial emissions or combustion processes. However, it is important to note that our analysis assumes a resident adult exposure scenario that may not accurately represent the real-world exposure at the hospital. Specifically, it is unlikely that an individual would be exposed 24 h a day, 7 days a week over a lifetime in this location, suggesting that our estimates might overstate the actual risk.
  • The inherent variability and uncertainty in the risk estimates are recognized, emphasizing the need for careful interpretation and further research, such as the employment of stochastic modeling.
This study improves our understanding of potential health impacts from exposure to hazardous air pollutants in this region, providing valuable insights for future environmental management strategies and policy decisions. Moreover, the benefits of this research, such as forecasting health impacts, could be utilized to optimize regulatory-based permitting, public education, emergency room care, and availability of response specialists, with the end goal of improving overall healthcare management that directly considers the health impacts of air quality on the exposed population of Kuwait.

6. Future Work

The proposed methodology, as presented, sets the stage for future work. In this phase, a stochastic model will be developed, incorporating probabilistic distributions of variables such as body weight, inhalation rate, metabolic rates, and ingestion rates. This approach will address variability and uncertainty in estimating risks and will be made possible through the development and implementation of the proposed methodology. Additionally, Step 3 of the proposed methodology, which involves estimating the concentrations of hazardous air pollutants at the point-of-contact for receptor populations by conducting an environmental transport-and-fate analysis, will be used in part to estimate the ecological screening quotient for country-wide ecological risk assessments. This assessment will also consider the potential pathways through which ecological changes might affect human health, such as impacts on food chains, water quality, and biodiversity. By integrating these ecological aspects with the direct human health risk factors, the methodology aims to provide a more holistic approach to both human health and ecological integrity.

Author Contributions

Conceptualization, M.M., J.V.G.T., and R.F.; Data curation, M.M., J.V.G.T., and A.R.; Formal analysis, M.M.; Funding acquisition, J.V.G.T.; Investigation, M.M.; Methodology, M.M. and A.R.; Project administration, J.V.G.T. and R.F.; Resources, J.V.G.T. and R.F.; Software, M.M. and J.V.G.T.; Supervision, J.V.G.T. and R.F.; Validation, M.M.; Visualization, M.M.; Writing—original draft, M.M.; Writing—review and editing, M.M., J.V.G.T., R.F., and B.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research and the APC were funded by Lakes Environmental Software.

Informed Consent Statement

Not applicable.

Data Availability Statement

A subset of the data is available upon request to the corresponding author, due to confidentiality constraints on certain portions of the dataset.

Acknowledgments

The authors would like to thank Jeff Stoakes from the Air Toxics Risk Assessments department at the Office of Air Quality, Indiana Department of Environmental Management, for his helpful comments and review of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations and Nomenclature

Abbreviations/NomenclatureMeaning
µg/m3micrograms per cubic meter
ADDaverage daily dose (inhalation, ingestion)
APIAmerican Petroleum Institute
AQZair quality zone
ATaveraging time
BTEXbenzene, toluene, ethylbenzene, and xylene
BWbody weight
Cairconcentration of pollutant(s) in the air
Cmediumconcentration of pollutant(s) in the medium
EDexposure duration
EFexposure frequency
ETexposure time
ggram
HAPshazardous air pollutants
HHRAPHuman Health Risk Assessment Protocol
IngRingestion rate
InhRinhalation rate
ISCST3Industrial Source Complex Short Term 3
Lliter
mgmilligram
NASANational Aeronautics and Space Administration
PAHspolycyclic aromatic hydrocarbons
SRTMShuttle Radar Topography Mission
TFERtransport, fate, exposure, and risk characterization
U.S. EPAU.S. Environmental Protection Agency
VOCsvolatile organic compounds
WRFWeather Research and Forecasting

References

  1. Central Intelligence Agency. Kuwait—The World Factbook. Available online: https://www.cia.gov/the-world-factbook/countries/kuwait/ (accessed on 12 April 2023).
  2. Khalaf, F.I.; Al-Ajmi, D. Aeolian Processes and Sand Encroachment Problems in Kuwait. Geomorphology 1993, 6, 111–134. [Google Scholar] [CrossRef]
  3. Al-Hemoud, A.; Al-Enezi, A.; Al-Dashti, H.; Petrov, P.; Misak, R.; AlSaraf, M.; Malek, M. Hazard Assessment and Hazard Mapping for Kuwait. Int. J. Disaster Risk Sci. 2023, 14, 143–161. [Google Scholar] [CrossRef]
  4. Barkley, M.P.; González Abad, G.; Kurosu, T.P.; Spurr, R.; Torbatian, S.; Lerot, C. OMI Air-Quality Monitoring over the Middle East. Atmos. Chem. Phys. 2017, 17, 4687–4709. [Google Scholar] [CrossRef] [Green Version]
  5. Al-Hurban, A.; Khader, S.; Alsaber, A.; Pan, J. Air Quality Assessment in the State of Kuwait during 2012 to 2017. Atmosphere 2021, 12, 678. [Google Scholar] [CrossRef]
  6. Al-Awadhi, J.M. Measurement of Air Pollution in Kuwait City Using Passive Samplers. Atmos. Clim. Sci. 2014, 4, 253–271. [Google Scholar] [CrossRef] [Green Version]
  7. Al-Hemoud, A.; Gasana, J.; Al-Dabbous, A.; Alajeel, A.; Al-Shatti, A.; Behbehani, W.; Malak, M. Exposure Levels of Air Pollution (PM2.5) and Associated Health Risk in Kuwait. Environ. Res. 2019, 179 Pt A, 108730. [Google Scholar] [CrossRef]
  8. Achilleos, S.; Al-Ozairi, E.; Alahmad, B.; Garshick, E.; Neophytou, A.M.; Bouhamra, W.; Yassin, M.F.; Koutrakis, P. Acute Effects of Air Pollution on Mortality: A 17-Year Analysis in Kuwait. Environ. Int. 2019, 126, 476–483. [Google Scholar] [CrossRef]
  9. Cohen, A.J.; Brauer, M.; Burnett, R.; Anderson, H.R.; Frostad, J.; Estep, K.; Balakrishnan, K.; Brunekreef, B.; Dandona, L.; Dandona, R.; et al. Estimates and 25-Year Trends of the Global Burden of Disease Attributable to Ambient Air Pollution: An Analysis of Data from the Global Burden of Diseases Study 2015. Lancet 2017, 389, 1907–1918. [Google Scholar] [CrossRef] [Green Version]
  10. Freeman, B.; Gharabaghi, B.; Thé, J.; Munshed, M.; Faisal, S.; Abdullah, M.; Al Aseed, A. Mapping Air Quality Zones for Coastal Urban Centers. J. Air Waste Manag. Assoc. 2017, 67, 565–581. [Google Scholar] [CrossRef] [Green Version]
  11. U.S. Environmental Protection Agency. Human Health Risk Assessment Protocol (HHRAP) for Hazardous Waste Combustion Facilities, Final; U.S. Environmental Protection Agency: Washington, DC, USA, 2005.
  12. Munshed, M.; Van Griensven Thé, J.; Fraser, R. Methodology for Mobile Toxics Deterministic Human Health Risk Assessment and Case Study. Atmosphere 2023, 14, 506. [Google Scholar] [CrossRef]
  13. Agency for Toxic Substances and Disease Registry. Toxicological Profile for Methyl Mercaptan. Available online: https://wwwn.cdc.gov/TSP/ToxProfiles/ToxProfiles.aspx?id=224&tid=40 (accessed on 12 April 2023).
  14. Mallah, M.A.; Mallah, M.A.; Liu, Y.; Xi, H.; Wang, W.; Feng, F.; Zhang, Q. Relationship between Polycyclic Aromatic Hydrocarbons and Cardiovascular Diseases: A Systematic Review. Front. Public Health 2021, 9, 763706. [Google Scholar] [CrossRef] [PubMed]
  15. Burstyn, I.; Kromhout, H.; Partanen, T.; Svane, O.; Langård, S.; Ahrens, W.; Kauppinen, T.; Stücker, I.; Shaham, J.; Heederik, D.; et al. Polycyclic Aromatic Hydrocarbons and Fatal Ischemic Heart Disease. Epidemiology 2005, 16, 744–750. [Google Scholar] [CrossRef] [PubMed]
  16. Gustavsson, P.; Plato, N.; Hallqvist, J.; Hogstedt, C.; Lewné, M.; Reuterwall, C.; Schéele, P. A Population-Based Case-Referent Study of Myocardial Infarction and Occupational Exposure to Motor Exhaust, Other Combustion Products, Organic Solvents, Lead, and Dynamite. Epidemiology 2001, 12, 222–228. [Google Scholar] [CrossRef] [PubMed]
  17. Xu, X.; Cook, R.L.; Ilacqua, V.A.; Kan, H.; Talbott, E.O.; Kearney, G. Studying Associations between Urinary Metabolites of Polycyclic Aromatic Hydrocarbons (PAHs) and Cardiovascular Diseases in the United States. Sci. Total Environ. 2010, 408, 4943–4948. [Google Scholar] [CrossRef]
  18. Jackson, A.C. Chronic Neurological Disease due to Methylmercury Poisoning. Can. J. Neurol. Sci. J. Can. Des Sci. Neurol. 2018, 45, 620–623. [Google Scholar] [CrossRef] [Green Version]
  19. Lee, H.; Lee, M.W.; Warren, J.R.; Ferrie, J. Childhood Lead Exposure Is Associated with Lower Cognitive Functioning at Older Ages. Sci. Adv. 2022, 8, eabn5164. [Google Scholar] [CrossRef]
  20. Reuben, A.; Caspi, A.; Belsky, D.W.; Broadbent, J.; Harrington, H.; Sugden, K.; Houts, R.M.; Ramrakha, S.; Poulton, R.; Moffitt, T.E. Association of Childhood Blood Lead Levels with Cognitive Function and Socioeconomic Status at Age 38 Years and with IQ Change and Socioeconomic Mobility between Childhood and Adulthood. JAMA 2017, 317, 1244. [Google Scholar] [CrossRef] [Green Version]
  21. Mansouri, M.T.; Muñoz-Fambuena, I.; Cauli, O. Cognitive Impairment Associated with Chronic Lead Exposure in Adults. Neurol. Psychiatry Brain Res. 2018, 30, 5–8. [Google Scholar] [CrossRef]
  22. Sanders, T.; Liu, Y.; Buchner, V.; Tchounwou, P.B. Neurotoxic Effects and Biomarkers of Lead Exposure: A Review. Rev. Environ. Health 2009, 24, 15–45. [Google Scholar] [CrossRef] [Green Version]
  23. Eskenazi, B.; Ames, J.; Rauch, S.; Signorini, S.; Brambilla, P.; Mocarelli, P.; Siracusa, C.; Holland, N.; Warner, M. Dioxin Exposure Associated with Fecundability and Infertility in Mothers and Daughters of Seveso, Italy. Hum. Reprod. 2020, 36, 794–807. [Google Scholar] [CrossRef]
  24. Mocarelli, P.; Gerthoux, P.M.; Patterson, D.G.; Milani, S.; Limonta, G.; Bertona, M.; Signorini, S.; Tramacere, P.; Colombo, L.; Crespi, C.; et al. Dioxin Exposure, from Infancy through Puberty, Produces Endocrine Disruption and Affects Human Semen Quality. Environ. Health Perspect. 2008, 116, 70–77. [Google Scholar] [CrossRef] [Green Version]
  25. Chen, S.C.; Liao, T.L.; Wei, Y.H.; Tzeng, C.R.; Kao, S.H. Endocrine Disruptor, Dioxin (TCDD)-Induced Mitochondrial Dysfunction and Apoptosis in Human Trophoblast-like JAR Cells. Mol. Hum. Reprod. 2010, 16, 361–372. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Pan, X.; Liu, X.; Li, X.; Niu, N.; Yin, X.; Li, N.; Yu, Z. Association between Environmental Dioxin-Related Toxicants Exposure and Adverse Pregnancy Outcome: Systematic Review and Meta-Analysis. Int. J. Fertil. Steril. 2015, 8, 351–366. [Google Scholar] [CrossRef]
  27. Huff, J. Benzene-Induced Cancers: Abridged History and Occupational Health Impact. Int. J. Occup. Environ. Health 2007, 13, 213–221. [Google Scholar] [CrossRef]
  28. Infante, P.F. Benzene Exposure and Multiple Myeloma: A Detailed Meta-Analysis of Benzene Cohort Studies. Ann. N. Y. Acad. Sci. 2006, 1076, 90–109. [Google Scholar] [CrossRef] [PubMed]
  29. Yin, S.-N.; Hayes, R.B.; Linet, M.S.; Li, G.-L.; Dosemeci, M.; Travis, L.B.; Li, C.-Y.; Zhang, Z.-N.; Li, D.-G.; Chow, W.-H.; et al. A cohort study of cancer among benzene-exposed workers in China: Overall results. Am. J. Ind. Med. 1996, 29, 227–235. [Google Scholar] [CrossRef]
  30. Smith, M.T.; Jones, R.M.; Smith, A.H. Benzene Exposure and Risk of Non-Hodgkin Lymphoma. Cancer Epidemiol. Prev. Biomark. 2007, 16, 385–391. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Adegoke, O.J.; Blair, A.; Shu, X.O.; Sanderson, M.; Jin, F.; Dosemeci, M.; Addy, C.L.; Zheng, W. Occupational History and Exposure and the Risk of Adult Leukemia in Shanghai. Ann. Epidemiol. 2003, 13, 485–494. [Google Scholar] [CrossRef] [PubMed]
  32. Kang, D.S.; Kim, H.S.; Jung, J.-H.; Lee, C.M.; Ahn, Y.-S.; Seo, Y.R. Formaldehyde Exposure and Leukemia Risk: A Comprehensive Review and Network-Based Toxicogenomic Approach. Genes Environ. 2021, 43, 1–10. [Google Scholar] [CrossRef]
  33. Chen, W.-Q.; Zhang, X.-Y. 1,3-Butadiene: A Ubiquitous Environmental Mutagen and Its Associations with Diseases. Genes Environ. 2022, 44, 1–22. [Google Scholar] [CrossRef]
  34. Redmon, J.H.; Kondash, A.J.; Womack, D.; Lillys, T.; Feinstein, L.; Cabrales, L.; Weinthal, E.; Vengosh, A. Is Food Irrigated with Oilfield-Produced Water in the California Central Valley Safe to Eat? A Probabilistic Human Health Risk Assessment Evaluating Trace Metals Exposure. Risk Anal. 2020, 41, 1463–1477. [Google Scholar] [CrossRef] [PubMed]
  35. Shih, Y.-H.; Kasaon, S.J.; Tseng, C.-H.; Wang, H.-C.; Chen, L.-L.; Chang, Y.-M. Health Risks and Economic Costs of Exposure to PCDD/Fs from Open Burning: A Case Study in Nairobi, Kenya. Air Qual. Atmos. Health 2015, 9, 201–211. [Google Scholar] [CrossRef]
  36. Schreck, E.; Bonnard, R.; Laplanche, C.; Leveque, T.; Foucault, Y.; Dumat, C. DECA: A New Model for Assessing the Foliar Uptake of Atmospheric Lead by Vegetation, Using Lactuca Sativa as an Example. J. Environ. Manag. 2012, 112, 233–239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Bseibsu, A.; Madhuranthakam, C.M.R.; Yetilmezsoy, K.; Almansoori, A.; Elkamel, A. Numerical Simulation of Dispersion Patterns and Air Emissions for Optimal Location of New Industries Accounting for Environmental Risks. Pollutants 2022, 2, 444–461. [Google Scholar] [CrossRef]
  38. Iqbal, A.; Afroze, S.; Rahman, M.M. Probabilistic Health Risk Assessment of Vehicular Emissions as an Urban Health Indicator in Dhaka City. Sustainability 2019, 11, 6427. [Google Scholar] [CrossRef] [Green Version]
  39. Pratt, G.C.; Dymond, M. Multipathway Screening Factors for Assessing Risks from Ingestion Exposures to Air Pollutants. J. Air Waste Manag. Assoc. 2009, 59, 419–429. [Google Scholar] [CrossRef]
  40. Malakan, W.; Thepanondh, S.; Kondo, A. Modeling of Inhalation Health Risk of Volatile Organic Compounds in the Vicinity of Maptaphut Petroleum and Petrochemical Industrial Estate, Thailand. Sustainability 2022, 14, 12073. [Google Scholar] [CrossRef]
  41. Tsai, J.-H.; Gu, W.-T.; Chung, I.-I.; Chiang, H.-L. Airborne Air Toxics Characteristics and Inhalation Health Risk Assessment of a Metropolitan Industrial Complex. Aerosol Air Qual. Res. 2019, 19, 2477–2489. [Google Scholar] [CrossRef]
  42. Pratt, G.; Vadali, M.; Kvale, D.; Ellickson, K. Traffic, Air Pollution, Minority and Socio-Economic Status: Addressing Inequities in Exposure and Risk. Int. J. Environ. Res. Public Health 2015, 12, 5355–5372. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Ollson, C.A.; Knopper, L.D.; Whitfield Aslund, M.L.; Jayasinghe, R. Site Specific Risk Assessment of an Energy-From-Waste Thermal Treatment Facility in Durham Region, Ontario, Canada. Part A: Human Health Risk Assessment. Sci. Total Environ. 2014, 466–467, 345–356. [Google Scholar] [CrossRef] [Green Version]
  44. Morselli, L.; Passarini, F.; Piccari, L.; Vassura, I.; Bernardi, E. Risk Assessment Applied to Air Emissions from a Medium-Sized Italian MSW Incinerator. Waste Manag. Res. 2010, 29 (Suppl. S10), S48–S56. [Google Scholar] [CrossRef] [PubMed]
  45. Roberts, R.J.; Chen, M. Waste Incineration—How Big Is the Health Risk? A Quantitative Method to Allow Comparison with Other Health Risks. J. Public Health 2006, 28, 261–266. [Google Scholar] [CrossRef] [PubMed]
  46. U.S. Environmental Protection Agency. AP-42: Compilation of Air Emissions Factors; Air Emissions Factors and Quantification; U.S. Environmental Protection Agency: Research Triangle Park, NC, USA, 2023. Available online: https://www.epa.gov/air-emissions-factors-and-quantification/ap-42-compilation-air-emissions-factors#:~:text=AP%2D42%2C%20Compilation%20of%20Air (accessed on 24 April 2023).
  47. U.S. Environmental Protection Agency. WATER9: The Model for Water Treatment; Version 3.0; U.S. Environmental Protection Agency: Washington, DC, USA, 2006. Available online: https://www.epa.gov/chief/water9-version-30 (accessed on 26 April 2023).
  48. U.S. Environmental Protection Agency. Landfill Gas Emissions Model (LandGEM); Version 3.02; U.S. Environmental Protection Agency: Washington, DC, USA, 2020. Available online: https://www.epa.gov/sites/default/files/2020-06/landgem-v303.xlsm (accessed on 26 April 2023).
  49. GTI Energy. GRI-GLYCalc; Version 4.0; GTI Energy: Chicago, IL, USA, 2019; Available online: https://sales.gastechnology.org/000102.html (accessed on 28 April 2023).
  50. American Petroleum Institute. Amine Unit Emissions Model AMINECalc; Version 1.0; American Petroleum Institute: Washington, DC, USA, 1999; Available online: https://tinyurl.com/4tfmhhne (accessed on 28 April 2023).
  51. Radian Corporation. General QA/QC Methods—Final Report; U.S. Environmental Protection Agency: Washington, DC, USA, 1997. Available online: https://www.epa.gov/sites/default/files/2015-08/documents/vi03.pdf (accessed on 8 May 2023).
  52. Obaid, J.; Ramadan, A.; Elkamel, A.; Anderson, W. Comparing Non-Steady State Emissions under Start-up and Shut-down Operating Conditions with Steady State Emissions for Several Industrial Sectors: A Literature Review. Energies 2017, 10, 179. [Google Scholar] [CrossRef] [Green Version]
  53. U.S. Environmental Protection Agency. AERMOD Modeling System Development. Available online: https://www.epa.gov/scram/aermod-modeling-system-development (accessed on 3 March 2023).
  54. Cooper, C.D.; Alley, F.C. Air Pollution Control: A Design Approach; Waveland Press: Long Grove, IL, USA, 2011. [Google Scholar]
  55. Munshed, M. Mobile Toxics Human Health Risk Assessment Framework. Master’s Thesis, University of Waterloo, Waterloo, ON, Canada, 2018. Available online: https://uwspace.uwaterloo.ca/handle/10012/13293 (accessed on 8 May 2023).
  56. Skamarock, W.; Klemp, J.; Dudhia, J.; Gill, D.; Barker, D.; Duda, M.; Huang, X.; Wang, W.; Powers, J. A Description of the Advanced Research WRF Version 3; University Corporation for Atmospheric Research: Boulder, CO, USA, 2008. [Google Scholar]
  57. IRAP-h View, 5.1.1; Lakes Environmental Software: Waterloo. 2022. Available online: https://www.weblakes.com/software/risk-assessment/irap-h-view (accessed on 10 May 2023).
  58. U.S. Environmental Protection Agency. Guidelines for Exposure Assessment; U.S. Environmental Protection Agency: Washington, DC, USA, 1992. Available online: https://cfpub.epa.gov/ncea/risk/recordisplay.cfm?deid=15263 (accessed on 30 April 2023).
  59. Lee, C.C.; Lin, S.D. Handbook of Environmental Engineering Calculations, 2nd ed.; McGraw-Hill Education: New York, NY, USA, 2007; pp. 4.3–4.113. [Google Scholar]
  60. Region 6, U.S. EPA Multimedia Planning and Permitting Division. Region 6 Risk Management Addendum—Draft Human Health Risk Assessment Protocol for Hazardous Waste Combustion Facilities; U.S. Environmental Protection Agency: Dallas, TX, USA, 1998. [Google Scholar]
  61. Mason, R.P.; Reinfelder, J.R.; Morel, F.M.M. Bioaccumulation of Mercury and Methylmercury. Water Air Soil Pollut. 1995, 80, 915–921. [Google Scholar] [CrossRef]
  62. Hansen, J.C.; Danscher, G. Organic Mercury—An Environmental Threat to the Health of Exposed Societies? Rev. Environ. Health 1997, 12, 107–116. [Google Scholar] [CrossRef]
  63. Pratt, G.C.; Dymond, M.; Ellickson, K.; Thé, J. Validation of a Novel Air Toxic Risk Model with Air Monitoring. Risk Anal. 2011, 32, 96–112. [Google Scholar] [CrossRef]
  64. National Safety Council. Odds of Dying. Injury Facts. Available online: https://injuryfacts.nsc.org/all-injuries/preventable-death-overview/odds-of-dying/ (accessed on 6 May 2023).
Figure 1. Flow diagram illustrating the application of a multi-pathway health risk assessment methodology in Kuwait’s air quality zones.
Figure 1. Flow diagram illustrating the application of a multi-pathway health risk assessment methodology in Kuwait’s air quality zones.
Atmosphere 14 01247 g001
Figure 2. A visual representation of the Multi-pathway Human Exposure to Hazardous Air Pollutants.
Figure 2. A visual representation of the Multi-pathway Human Exposure to Hazardous Air Pollutants.
Atmosphere 14 01247 g002
Figure 3. Simplified Example of a Risk Driver Analysis.
Figure 3. Simplified Example of a Risk Driver Analysis.
Atmosphere 14 01247 g003
Figure 4. Kuwait Air Quality Zones. Map Data: Lakes Satellite.
Figure 4. Kuwait Air Quality Zones. Map Data: Lakes Satellite.
Atmosphere 14 01247 g004
Table 1. Human Health Effects of Selected Hazardous Air Pollutants (HAPs).
Table 1. Human Health Effects of Selected Hazardous Air Pollutants (HAPs).
Hazardous Air
Pollutants
Emission Source
Type
Major Exposure RouteOrgan/System
Affected
Citation
Methyl mercaptanPoint, area, volumeInhalationRespiratory[13]
Polycyclic aromatic hydrocarbons (PAHs)Point, area, volume, on-road, and non-road mobile sourcesInhalation; oralCardiovascular[14,15,16,17]
MercuryPoint, area, volumeInhalation; oralCentral nervous and peripheral nervous systems[18]
LeadPoint, area, volumeInhalation; oralCentral nervous and peripheral nervous systems[19,20,21]
DioxinsPoint, area, volume, on-road, and non-road mobile sourcesOralReproductive[23,24,25,26]
BenzenePoint, area, volume, on-road, and non-road mobile sourcesInhalation; oralHematological and respiratory[27,28,29,30,31]
FormaldehydePoint, area, volume, on-road, and non-road mobile sourcesInhalation; oralHematological[32]
1,3-ButadienePoint, area, volume, on-road, and non-road mobile sourcesInhalationHematological and immune[33]
Table 2. Summary of Human Health Risk Assessments in Literature.
Table 2. Summary of Human Health Risk Assessments in Literature.
Multi-SourceMulti-PathwayCumulativeCitation
YesNoNo[34] 2020; [35] 2015; [36] 2012
NoYesNo[37] 2022; [38] 2019; [39] 2009
YesNoYes[40] 2022; [41] 2019; [42] 2015
NoYesYes[12] 2023; [43] 2014; [44] 2010; [45] 2006
YesYesYesThis Work
Table 3. Overview of EPA and Industry Methods/Algorithms for Estimating Emissions from Various Sources.
Table 3. Overview of EPA and Industry Methods/Algorithms for Estimating Emissions from Various Sources.
Method/AlgorithmObjectiveCitation
EPA 1 AP-42 emission factorsEstimate emissions from various source categories, including point sources (e.g., industrial stacks)[46]
AP-42, Fifth Edition, Volume I,
Chapter 7: Liquid Storage Tanks
Calculate volatile organic compound (VOC) and hazardous air pollutant (HAP) emissions from floating- and fixed-roof storage tanks[46]
EPA WATER9Estimate emissions from wastewater treatment plants[47]
EPA LandGEMCalculate emissions from landfills, utilizing the methane generation rate and potential methane generation capacity parameters[48]
GRI-GLYCalcEstimate emissions from glycol dehydration units[49]
API 2 AMINECalcEstimate emissions from amine gas treatment plants[50]
1 U.S. Environmental Protection Agency. 2 American Petroleum Institute.
Table 4. AERMOD Model Options and Settings Used in This Study.
Table 4. AERMOD Model Options and Settings Used in This Study.
Model Option/SettingSetting
AERMOD ExecutableVersion 22112
Dispersion OptionsNon-default regulatory option selected
Fast all sources (FASTALL)
Flat (FLAT)
Calculation TypeUnitized (unit emission rate concept representing the µg/m3 impact per 1 g/s of emissions) (refer to Section 2.2.1)
OutputConcentration, total deposition, dry deposition, and wet deposition
Dispersion CoefficientRural
PollutantsBenzene, formaldehyde, toluene, and benzo(a)pyrene
Averaging Periods1 hour and annual
Source TypesPoint
ReceptorsUniform Cartesian grid and discrete Cartesian receptors (sensitive receptors)
TerrainTerrain in Kuwait can be approximated as flat; however, terrain files were used for completeness in the model
Meteorological Data Files2017 hourly meteorological data, contained in surface and profile files, for Kuwait’s three air quality zones were generated utilizing the WRF model
Table 5. Exposure Pathways Used in Multi-pathway Exposure Assessments [11].
Table 5. Exposure Pathways Used in Multi-pathway Exposure Assessments [11].
Exposure Pathways
Direct inhalation of vapors and particles
Incidental ingestion of soil
Ingestion of drinking water from treated surface water sources
Ingestion of homegrown produce 1,2
1 Applies to farms in Wafra, Al-Abdali, and Al-Sulaibiya areas. 2 https://www.trade.gov/market-intelligence/kuwait-agriculture (accessed on 12 April 2023).
Table 6. Summary of the Number of Modeled Facilities and Corresponding Number of Emission Sources in the Air Quality Zones of Kuwait.
Table 6. Summary of the Number of Modeled Facilities and Corresponding Number of Emission Sources in the Air Quality Zones of Kuwait.
Air Quality ZoneNumber of Modeled
Facilities
Number of Emission Sources within
Modeled Facilities
Inland421
Production721
Coastal1047
Table 7. Hazardous Air Pollutants Considered in the Study.
Table 7. Hazardous Air Pollutants Considered in the Study.
Hazardous Air Pollutants
Benzene
Formaldehyde
Toluene
Benzo(a)pyrene
Table 8. Cumulative Cancer Risk.
Table 8. Cumulative Cancer Risk.
Sensitive Receptor/Exposure ScenariosResident AdultResident Child
Receptor 14.79 × 10−119.63 × 10−12
Receptor 24.74 × 10−99.55 × 10−10
Receptor 31.02 × 10−82.05 × 10−9
Receptor 43.77 × 10−111.76 × 10−11
Table 9. Cumulative Non-cancer Risk.
Table 9. Cumulative Non-cancer Risk.
Sensitive Receptor/Exposure ScenariosResident AdultResident Child
Receptor 18.83 × 10−78.94 × 10−7
Receptor 29.30 × 10−59.86 × 10−5
Receptor 32.02 × 10−42.16 × 10−4
Receptor 43.45 × 10−88.93 × 10−8
Table 10. Cumulative Acute Non-cancer Risk.
Table 10. Cumulative Acute Non-cancer Risk.
Sensitive Receptor/Exposure ScenariosResident Adult/Child
Receptor 14.47 × 10−5
Receptor 22.21 × 10−4
Receptor 33.00 × 10−4
Receptor 48.61 × 10−10
Table 11. Cumulative Cancer Risk.
Table 11. Cumulative Cancer Risk.
Sensitive Receptor/Exposure ScenariosResident AdultResident Child
Receptor 13.11 × 10−116.26 × 10−12
Receptor 22.34 × 10−94.73 × 10−10
Receptor 31.20 × 10−82.41 × 10−9
Receptor 43.83 × 10−111.76 × 10−11
Table 12. Cumulative Non-cancer Risk.
Table 12. Cumulative Non-cancer Risk.
Sensitive Receptor/Exposure ScenariosResident AdultResident Child
Receptor 15.76 × 10−75.86 × 10−7
Receptor 24.45 × 10−54.64 × 10−5
Receptor 32.38 × 10−42.55 × 10−4
Receptor 43.39 × 10−88.85 × 10−8
Table 13. Cumulative Acute Non-cancer Risk.
Table 13. Cumulative Acute Non-cancer Risk.
Sensitive Receptor/Exposure ScenariosResident Adult/Child
Receptor 12.49 × 10−5
Receptor 22.71 × 10−4
Receptor 32.71 × 10−4
Receptor 48.78 × 10−10
Table 14. Cumulative Cancer Risk.
Table 14. Cumulative Cancer Risk.
Sensitive Receptor/Exposure ScenariosResident AdultResident Child
General Ahmadi Hospital6.12 × 10−62.31 × 10−6
Fatima Bint Asad High School for Girls7.56 × 10−62.87 × 10−6
Ahmadi Hospital1.09 × 10−54.33 × 10−6
Ahmadi Zoo7.92 × 10−63.04 × 10−6
Adan Hospital4.05 × 10−61.52 × 10−6
Hilton Kuwait Resort3.63 × 10−61.31 × 10−6
Mosque (North Ahmadi)5.79 × 10−62.23 × 10−6
Table 15. Cumulative Non-cancer Risk.
Table 15. Cumulative Non-cancer Risk.
Sensitive Receptor/Exposure ScenariosResident AdultResident Child
General Ahmadi Hospital9.67 × 10−31.56 × 10−2
Fatima Bint Asad High School for Girls1.19 × 10−21.93 × 10−2
Ahmadi Hospital1.26 × 10−22.43 × 10−2
Ahmadi Zoo1.13 × 10−21.92 × 10−2
Adan Hospital6.02 × 10−39.92 × 10−3
Hilton Kuwait Resort6.84 × 10−31.00 × 10−2
Mosque (North Ahmadi)7.76 × 10−31.36 × 10−2
Table 16. Cumulative Acute Non-cancer Risk.
Table 16. Cumulative Acute Non-cancer Risk.
Sensitive Receptor/Exposure ScenariosResident Adult/Child
General Ahmadi Hospital1.19 × 10−1
Fatima Bint Asad High School for Girls1.16 × 10−1
Ahmadi Hospital9.69 × 10−2
Ahmadi Zoo1.05 × 10−1
Adan Hospital6.81 × 10−2
Hilton Kuwait Resort9.65 × 10−2
Mosque (North Ahmadi)8.29 × 10−2
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Munshed, M.; Van Griensven Thé, J.; Fraser, R.; Matthews, B.; Ramadan, A. Extending Multi-Pathway Human Health Risk Assessment from Regional to Country-Wide—A Case Study on Kuwait. Atmosphere 2023, 14, 1247. https://doi.org/10.3390/atmos14081247

AMA Style

Munshed M, Van Griensven Thé J, Fraser R, Matthews B, Ramadan A. Extending Multi-Pathway Human Health Risk Assessment from Regional to Country-Wide—A Case Study on Kuwait. Atmosphere. 2023; 14(8):1247. https://doi.org/10.3390/atmos14081247

Chicago/Turabian Style

Munshed, Mohammad, Jesse Van Griensven Thé, Roydon Fraser, Bryan Matthews, and Ashraf Ramadan. 2023. "Extending Multi-Pathway Human Health Risk Assessment from Regional to Country-Wide—A Case Study on Kuwait" Atmosphere 14, no. 8: 1247. https://doi.org/10.3390/atmos14081247

APA Style

Munshed, M., Van Griensven Thé, J., Fraser, R., Matthews, B., & Ramadan, A. (2023). Extending Multi-Pathway Human Health Risk Assessment from Regional to Country-Wide—A Case Study on Kuwait. Atmosphere, 14(8), 1247. https://doi.org/10.3390/atmos14081247

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop