Spatial patterns and background ranges of naturally occurring radionuclides (NORs) (i.e. U-238, Th-232, K-40) and Cs-137 were studied in the urban soils of Yerevan, the capital city of Armenia. Multifractal Inverse Distance Weighting (MIDW) was used to generate and analyze distribution patterns of radionuclide activities. Based on Fourier transformation of radioactivity data, a spectral analysis was also applied to separate, where possible, background/baseline patterns from local anomalies: two ranges of background values were found to characterise the Yerevan territory. Specifically, in the south and south-east of Yerevan, the lower background ranges of U-238, Th-232 and K-40 comprised in the intervals 2.60–36.42 Bq/kg, 4.04–30.63 Bq/kg and 147.7–396.7 Bq/kg, respectively, were observed in association with the presence of sedimentary formations. In contrast, the higher ones were found, instead, in the central and northern parts of the city where andesite-basalt lavas and ignimbrite tuffs occur. Here, the background values rise to 142.4 Bq/kg, 138.76 Bq/kg and 1502 Bq/kg, respectively. As for the distribution of artificial Cs-137, its baseline levels in Yerevan seem to depend mostly on the global radioactive fallout and some local technogenic sources. Its distribution patterns partially differ from those of NORs. In the framework of this paper, Radium equivalent activity (RaEq), outdoor absorbed dose rate in air (ODRA) and annual effective dose equivalent (AEDEs) were also determined and mapped. They show a good coincidence of their spatial variations with those of NORs. The Monte Carlo simulation was used to assess excess lifetime cancer risk from a stochastic perspective. The related sensitivity analysis revealed that, among NORs, U-238 and Th-232 give the greatest contribution to the total variance (45.7% 42.8%, respectively). In comparison, K-40 has the lowest share (11.3%). Regarding Cs-137, a highly negligible contribution to the onset of health risks (accounting for 0.02%) was observed.

Breast milk samples were collected from 51 mothers in a seaside city Dalian, where the residents usually have higher dietary exposure to polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs) and polychlorinated biphenyls (PCBs) due to the larger consumption of seafood. The lipid-basis concentrations of ∑Cl₂–₈DD/Fs, ∑Cl₂–₁₀Bs, and total toxicity equivalent (WHO-TEQ) were measured to be in the ranges of 35.7–2727.8 pg/g, 4.91–52.64 ng/g, and 2.27–36.30 pg/g, respectively. The average proportion of ∑Cl₂–₃DD/Fs was higher than that of ∑Cl₄–₇DD/Fs in the collected human breast milk samples, suggesting that the health risk of Cl₂–₃DD/Fs should be especially concerned. The concentration data of PCDD/Fs and PCBs in human breast milk essentially followed a positive skew probability distribution. Women in high-level exposure scenarios exhibited a higher potential to accumulate homologues Cl₄DFs, Cl₇DFs, Cl₈DF, and Cl₆Bs in breast milk. Three PCDD/F congeners (1,2,3,6,7,8-Cl₆DF, 1,2,3,4,7,8-Cl₆DF, and 1,2,3,4,6,7,8-Cl₇DD) and three PCB congeners (PCB 126, PCB 138, and PCB 169) were identified as good indicators for the accumulation of PCDD/Fs and PCBs in human breast milk, respectively. The food-to-milk accumulation factors (FMAF) were calculated to evaluate the accumulation potentials of different PCDD/F and PCB congeners in human breast milk via dietary exposure. The calculated FMAF value presented a non-monotonic variation with the logarithm of n-octanol–water partition coefficient (log KOW) with a peak at a log KOW value of about 7.3 and a valley at a log KOW value of about 8. The mean value of the estimated daily intake (EDI) of total WHO-TEQ for breast-fed infants in Dalian, predicted by Monte Carlo simulation, was 10 folds higher than the upper range of the tolerable daily intake (TDI) value (4 pg WHO-TEQ/kg bw/d), suggesting continued and enhanced efforts should be made to reduce the exposure risk of infants to PCDD/Fs and PCBs.

To better evaluate the potential risks of aquatic polycyclic aromatic compounds (PACs), multiple approaches have been implemented in this study to assess the human health and ecological risks of parent, nitrated and oxygenated polycyclic aromatic hydrocarbons (PAHs, NPAHs and OPAHs) in the surface water of Jialing and Yangtze Rivers in downtown Chongqing in southwestern China. The concentrations of ∑PAHs (334 ± 125 ng L⁻¹) were much higher than those of ∑OPAHs (20.2 ± 7.49 ng L⁻¹) in the two rivers, while NPAHs were barely detected. Concentrations of detected PACs were higher in wet season than dry season, probably resulted from the elevated particle input due to heavy rainfall in wet season. Concentrations of PAHs were higher in the particulate phase than dissolved phase, while OPAHs levels showed a reverse pattern. The partition coefficients (Kₚ) of PACs in the water-SPM (suspended particulate matter) system were mainly affected by SPM concentrations and octanol/water partition coefficients of specific PACs. Human health risks calculated from non-probabilistic risk assessment model and probabilistic risk assessment model based on Monte Carlo simulation showed similar data pattern with slight difference in absolute values. Both models revealed potential or even severe human health risks contributed mainly by dermal exposure to aquatic PACs in this study. Furthermore, these models also manifested that infant stage was highly sensitive for PAC exposure. Sensitivity analysis indicated that health risk results was most sensitive to Benzo[a]pyrene equivalent toxic concentration (BaPₑq), followed by showering time and daily water intake volume. Levels of ecological risks and contributions of individual PACs differed from models based on different quality values. The adequacy of toxicity data was crucial for the reliability of ecological risk assessment.

To evaluate pesticide regulatory standards in agricultural crops, we introduced a regulatory modeling framework that can flexibly evaluate a population’s aggregate exposure risk via maximum residue levels (MRLs) under good agricultural practice (GAP). Based on the structure of the aggregate exposure model and the nature of variable distributions, we optimized the framework to achieve a simplified mathematical expression based on lognormal variables including the lognormal sum approximation and lognormal product theorem. The proposed model was validated using Monte Carlo simulation, which demonstrates a good match for both head and tail ends of the distribution (e.g., the maximum error = 2.01% at the 99th percentile). In comparison with the point estimate approach (i.e., theoretical maximum daily intake, TMDI), the proposed model produced higher simulated daily intake (SDI) values based on empirical and precautionary assumptions. For example, the values at the 75th percentile of the SDI distributions simulated from the European Union (EU) MRLs of 13 common pesticides in 12 common crops were equal to the estimated TMDI values, and the SDI values at the 99th percentile were over 1.6-times the corresponding TMDI values. Furthermore, the model was refined by incorporating the lognormal distributions of biometric variables (i.e., food intake rate, processing factor, and body weight) and varying the unit-to-unit variability factor (VF) of the pesticide residues in crops. This ensures that our proposed model is flexible across a broad spectrum of pesticide residues. Overall, our results show that the SDI is significantly reduced, which may better reflect reality. In addition, using a point estimate or lognormal PF distribution is effective as risk assessments typically focus on the upper end of the distribution.

This study aimed to assess the carcinogenic and non-carcinogenic risks of in-vehicle exposure in Tehran, Iran to formaldehyde and acetaldehyde for different models of taxis, and to explore the effects of city zone, taxi vehicle type, the taxi's age (<1, 1–5, 5–10), fuel type (gasoline, CNG, and LPG), and refueling activities on the estimated health risks based on previously measured concentrations. The overall and age-specific carcinogenic and non-carcinogenic risks of these compounds for taxi drivers and passengers were estimated separately using Monte Carlo simulations. Three scenarios of exposure frequency were defined for taxis commuting in different zones of city: Restricted Traffic Zone (RTZ) and Odd-Even Zone (OEZ) as two plans to reduce air pollution, and no-restriction zone (NRZ). The carcinogenic risks for drivers and passengers, the average risks of formaldehyde and acetaldehyde for most cases were above the 1 × 10⁻⁴. The health risks were greater in Restricted Traffic Zone (RTZ) and Odd-Even Zone (OEZ) in comparison to no-restriction zone (NRZ). The carcinogenic risk from formaldehyde exposures were higher than those for acetaldehyde in all cases. Taxis fueled with LPG showed lower cancer risks for both acetaldehyde and formaldehyde. Refueling increased the carcinogenic risk from both compounds. For non-carcinogenic risks from acetaldehyde, the average hazard ratios for both drivers and passengers were >1, indicating a non-negligible risk. Cancer and non-cancer risks for the taxi drivers were greater than the passengers given the higher time of occupancy. The present study showed that transportation in taxis can impose significant long-term health risks to both passengers and drivers. Development and investment in cleaner choices for public transportations are required.

Metal contamination in aquatic environments is a severe global concern to human health and aquatic ecosystems. This study used several risk assessment indices, to evaluate metal (Cu, Zn, Pb, Cd and Cr) environmental risk of riverine surface sediments across the rural-urban interface of the rapidly developing Wen-Rui Tang River watershed in eastern China. Risk assessments were determined for 38 sites based on the potential ecological risk index (RI), consensus-based sediment quality guidelines (SQGs) and risk assessment code (RAC). Land-use cluster analysis showed that sediments were severely contaminated, especially for Cd, whose concentrations were ∼100 times higher than background levels and had a high proportion in the bioaccessible fraction. According to RI, ErCd was identified with extremely high risk potential, resulting in the highest ecological risk of Cluster 4 (industrial). Similarly, risk within Cluster 4 (industrial) was also ranked highest by SQGs assessment due to the high proportion of industrial land use. Zinc was determined with high risk due to its high concentration compared to its effect range medium (ERM) value. Discrepancies in predicting environmental risks from metals among the three indices were mainly attributed to the contrasting definitions of these metrics. Environmental risk uncertainty derived from spatial variation was further estimated by Monte Carlo simulation and ranked as: Zn > Cd > Cr > Pb > Cu. This comprehensive environmental risk assessment provides important information to guide remediation strategies for management of metal contamination at the watershed scale.

This study characterized spatio-temporal variations in the concentration of benzene, toluene, ethylbenzene, and xylene (BTEX) compounds in the vicinity of gas and compressed natural gas (CNG) stations in Tehran, Iran. Health risk assessment (HRA) was computed using Monte Carlo simulations (MCS) for evaluating inhalation lifetime cancer risk (LTCR), the hazard quotient (HQ), and sensitivity analysis (SA) for BTEX exposure in different age groups (birth to <81) and as a function of distance (0–250 m) from the center of the stations. For all monitoring stations, the average values of benzene, toluene, ethylbenzene, and xylene in winter were 466.09 ± 132.25, 873.13 ± 233.51, 493.05 ± 141.22, and 910.57 ± 145.40 μg m⁻³, respectively. The mean wintertime ratios of T/B for the 12 stations ranged from 1.69 to 2.04. Furthermore, there was no significant relationship between the concentration of BTEX with either the specific month or distance from the center of stations (p > 0.05). Factors promoting BTEX formation in the study region were fuel evaporation and gas/CNG station emissions. The LTCRs for the target compounds in the winter for different age groups and distances from the center of stations was limited to 2.11 × 10⁻⁴ to 1.82 × 10⁻³ and 2.30 × 10⁻⁴ to 2.01 × 10⁻³, respectively, which exceeded proposed values by U.S. EPA. Moreover, the HQs for BTEX for three age groups and distances were limited to between 2.89 × 10⁻⁵ and 9.33 × 10⁻², which were lower than the acceptable limit (HQs < 1). The results of this work are applicable to similar areas that are heavily populated with vehicular traffic. This study motivates a closer look at mitigation strategies to limit the health effects of carcinogenic emissions such as benzene and ethylbenzene from gas/CNG stations.

From a health prospective, it is critical to provide a comprehensive model which integrates all the parameters involved in virus transmission and its consequences on human body. In order to estimate the health risks, for workers and residents, associated with an exposure airborne viruses emitted from a wastewater treatment (WWTP), the concentration levels of viruses in emitted bioaerosols over a twelve-month period were measured by real-time polymerase chain reaction (RT-PCR). A combined Gaussian plum dispersion model and quantitative microbial risk assessment (QMRA) with Monte-Carlo simulation served as suitable explanatory tools to estimate the risk of acquiring gastrointestinal illness (GI) due to exposure to air containing Rotavirus (RoV) and Norovirus (NoV) bioaerosols. Additionally, DALY metric was applied to quantify the disability and mortality for workers and residents. RoV and NoV were detected above aeration tank with annual mean concentration 27 and 3099 (Viruses/m³.h), respectively. The medium calculated DALY indicator based on viral loads in contaminant source (RoV:5.76 × 10⁻² and NoV:1.23 × 10⁻¹) and estimated in different distances away (300–1000 m) (RoV:2.87 × 10⁻²- 2.75 × 10⁻² and NoV:1.14 × 10⁻¹-1.13 × 10⁻¹) were markedly higher than the threshold values recommended by US EPA (10⁻⁴ DALY pppy) and WHO (10⁻⁶ DALY pppy). The sensitivity analysis highlighted dose exposure and disease burden per case (DBPC) as two most influential factors for both workers and residents following exposure to two pathogens of concern. Due to high resistance and high concentration in the environment, the presence of RoV and NoV can intensify the consequences of diarrhea especially for children under five years of age; A comprehensible and transparent presentation of DALYs and QMRA can help decision makers and responsibilities to justify the priorities of exposure to wastewater in comparison with other risks of daily life.

In the present study, the daily dose in terms of particle surface area received by citizens living in five cities in Western countries, characterized by different lifestyle, culture, climate and built-up environment, was evaluated and compared. For this purpose, the exposure to sub-micron particle concentration levels of the population living in Barcelona (Spain), Cassino (Italy), Guilford (United Kingdom), Lund (Sweden), and Brisbane (Australia) was measured through a direct exposure assessment approach. In particular, measurements of the exposure at a personal scale were performed by volunteers (15 per each population) that used a personal particle counter for different days in order to obtain exposure data in microenvironments/activities they resided/performed. Non-smoking volunteers performing non-industrial jobs were considered in the study.Particle concentration data allowed obtaining the exposure of the population living in each city. Such data were combined in a Monte Carlo method with the time activity pattern data characteristics of each population and inhalation rate to obtain the most probable daily dose in term of particle surface area as a function of the population gender, age, and nationality.The highest daily dose was estimated for citizens living in Cassino and Guilford (>1000 mm²), whereas the lowest value was recognized for Lund citizens (around 100 mm²). Indoor air quality, and in particular cooking and eating activities, was recognized as the main influencing factor in terms of exposure (and thus dose) of the population: then confirming that lifestyle (e.g. time spent in cooking activities) strongly affect the daily dose of the population. On the contrary, a minor or negligible contribution of the outdoor microenvironments was documented.

Plant uptake factors (PUFs) are of great importance in human cadmium (Cd) exposure risk assessment while it has been often treated in a generic way. We collected 1077 pairs of vegetable-soil samples from production fields to characterize Cd PUFs and demonstrated their utility in assessing Cd exposure risks to consumers of locally grown vegetables. The Cd PUFs varied with plant species and pH and organic matter content of soils. Once normalized PUFs against soil parameters, the PUFs distributions were log-normal in nature. In this manner, the PUFs were represented by definable probability distributions instead of a deterministic figure. The Cd exposure risks were then assessed using the normalized PUF based on the Monte Carlo simulation algorithm. Factors affecting the extent of Cd exposures were isolated through sensitivity analyses. Normalized PUF would illustrate the outcomes for uncontaminated and slightly contaminated soils. Among the vegetables, lettuce was potentially hazardous for residents due to its high Cd accumulation but low Zn concentration. To protect 95% of the lettuce production from causing excessive Cd exposure risks, pH of soils needed to be 5.9 and above.