Next Article in Journal
Interstitial Ectopic Pregnancy: The Role of Mifepristone in the Medical Treatment
Previous Article in Journal
A Path Model for Burnout in Community Mental Health Professionals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Risk Assessment of Hazmat Road Transportation Considering Environmental Risk under Time-Varying Conditions

School of Business, East China University of Science and Technology, Shanghai 200237, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(18), 9780; https://doi.org/10.3390/ijerph18189780
Submission received: 30 June 2021 / Revised: 2 August 2021 / Accepted: 13 September 2021 / Published: 17 September 2021

Abstract

:
Hazardous materials shipments are integral to the development of industrial countries. Significant casualties and severe environmental pollution quickly ensue when accidents occur. Currently, relevant research on risk assessment of hazardous materials’ road transportation remains limited when both the population exposure risk and environmental risk are considered, especially in regard to analyzing the differences of accident impacts in different populations and environments. This paper adopts a Gaussian plume model to simulate dynamic areas at three levels of population exposure and assesses the pollution scope of air, groundwater, lakes, and rivers with a variety of diffusion models. Then, we utilize various costs to analyze the differences of accident impacts in population exposure and environmental pollution. Finally, a risk assessment model of hazardous materials road transportation under time-varying conditions is presented by considering the bearing capacity of the assessed area. Furthermore, this model is applied to a case study involving a risk assessment of hazardous materials transportation of a highly populated metropolitan area of Shanghai, China. The resulting analyses reveal that the safety of hazardous materials transportation could be effectively improved by controlling certain model parameters and avoiding road segments with a high risk of catastrophic accident consequences.

1. Introduction

Hazardous materials (hazmat) are flammable, explosive, toxic, corrosive and radioactive. Four billion tons of hazmat are transported worldwide, and approximately 400 million tons are shipped across China every year [1]. Hazmat accidents during road transport have caused catastrophic losses to humans and the environment worldwide. In July 2019, an oil truck in Nigeria lost control and turned over during driving, which caused heavy smog and led to 48 deaths and more than 90 injuries. In June 2020, a liquefied gas explosion accident occurred on an expressway in Zhejiang Province, China, causing 20 deaths and 172 hospitalizations. In September 2020, a truck carrying chlorine tanks exploded in Western Iran. The resulting explosion injured 217 people and released a large amount of chlorine into the air.
Therefore, risk assessment in hazmat transportation constitutes an important research area to enhance our understanding of hazmat transport management. In this paper, we analyze the differences in population exposure and environmental pollution due to accidents in different areas and propose a population exposure and environmental risk assessment model to better identify high-risk road segments, prevent the occurrence of catastrophic consequences, and provide a better guarantee for the safety of hazmat transportation.
The rest of this paper is organized as follows. Section 2 provides a literature review on hazmat transportation risk assessment and notes the differences between existing research and our contribution. Section 3 constructs the population exposure and environmental risk assessment model. The resulting analyses of the case study in Shanghai are presented in Section 4. Finally, conclusions are drawn in Section 5.

2. Literature Review

In previous studies, various risk assessment models have been presented to assess the potential risk of hazmat transportation.
In 1980, a traditional risk (TR) model was proposed [2]. Suppose a road network G = (U, V), consisting of a node set U to represent road interactions and a link set V to represent road segments; node i, j ϵ U, and segment (i, j) ϵ V. Then, the hazmat transportation risk R along a route is equal to the sum of the product of the accident probability Pij and accident consequence Cij of each road segment, as expressed in Equation (1).
R = ( i , j ) V P i j × C i j
Considering risk preference, Abkowitz et al. [3] proposed a perceived risk (PR) model by introducing risk preference parameter q to reflect the level of decision maker risk preference.
Then, Sivakumar et al. [4] proposed a conditional risk (CR) model, which evaluates the risk after the occurrence of the first hazmat accident. Erkut et al. [5] developed the minimum maximum (MM) model, mean-variance (MV) model, and disutility (DU) model.
Subsequently, Kang [6] and Kang et al. [7] applied the VaR risk measurement tool ubiquitous in financial risk management to hazmat shipments and proposed the value at risk (VaR) model to assess the risk of hazmat transportation. Toumazis et al. [8] improved the VaR model and developed the conditional value at risk (CVaR) model. The VaR model only focuses on the part of the risk not exceeding the VaR value, while the CVaR model focuses on the part of the risk exceeding the CVaR value [9].
Risk assessment under time-varying conditions is an emerging analysis topic to identify crash-prone traffic conditions and improve the safety of transportation. Wei [10] proposed a time-varying risk (TVR) model, which studies the change characteristics of accident probability and population density under time-varying conditions. Hossain et al. [11] applied a Bayesian belief net (BBN) to build a crash prediction model based on real-time traffic condition data retrieved from a high-resolution detector. Yu et al. [12] introduced a support vector machine (SVM) to evaluate real-time crash risk and analyze the relationship between risk and real-time traffic data detected by 30 RTMS radars. Sun et al. [13] proposed a dynamic Bayesian network (DBN) model of time sequence traffic data to investigate the relationship between crash occurrence and dynamic speed condition data. Wang et al. [14] presented a multilevel Bayesian logistic regression model of crash risk considering traffic, geometric and weather factors under time-varying conditions. Xu et al. [15] investigated the impacts of real-time traffic flow conditions on crash casualties under different collision types via geometric and traffic data. Yang et al. [16] introduced Bayesian dynamic logistic regression (Bayesian dynamic LR) to develop a real-time crash risk evaluation model with in-field streaming traffic data. Yu et al. [17] proposed a real-time crash risk analysis approach with multidimensional traffic flow input features to establish relationships between crash occurrence probability and precrash traffic operational conditions. Li et al. [18] built a real-time crash risk model of arterials with the long–short-term memory convolutional neural network (LSTM-CNN), which explicitly learns from various features, such as traffic flow characteristics, signal timing, and weather conditions.
In the risk assessment model of hazmat transportation, in addition to population exposure, environmental pollution should be considered [19].
On the one hand, most works in the literature have considered population exposure risk. Fabiano et al. [20] modified the accident probability of hazmat transportation taking into account the aspects of road characteristics, weather and traffic flow and analyzed the number of exposed people. Zhang et al. [21] applied the information diffusion theory to estimate the accident probability of hazardous chemical transportation and simulated the accident scope and number of people affected in liquefied petroleum gas leakage accidents at different levels. Fan et al. [22] adopted the product of the accident probability and number of affected populations within the determined scope along a given road as the risk. Zero et al. [23] associated the triangular fuzzy number (TFN) of population exposure in hazmat transportation accidents with the risk of the network. Dong et al. [24] and Sun [25] adopted a population exposure risk objective in their optimization model. In addition, Khanmohamadi et al. [26] and Orozco et al. [27] utilized ALOHA software to simulate the influence scope and identify the number of people exposed. Most studies have paid more attention to population. However, when assessing the accident consequences of population exposure, less consideration is given to the differences in the impact of accidents on populations in different areas.
On the other hand, there are few studies that consider both population exposure risk and environmental risk. Ren [28] and Xia et al. [29] analyzed the number of population exposures and the scope of environmental pollution. Saat et al. [30] estimated the environmental cleaning costs of soil and groundwater and the evacuation costs associated with population exposure caused by hazmat railway transportation accidents according to the accident scope specified in the Emergency Response Guidelines (ERG) of the United States. Cordeiro et al. [31] employed a multicriteria analysis method to assess accident consequences such as population vulnerability and environmental vulnerability and developed a risk assessment model. Huang et al. [32] examined the vulnerability of social and environmental functional areas along a road and established a risk assessment of the area along the highway (RAAH) model. Wang et al. [33] comprehensively calculated losses in terms of casualties and environmentally sensitive areas when assessing the time-varying risk of road segments. In addition, research on diffusion models and environmental bearing capacity has been relatively comprehensive in the field of environmental science. A diffusion model simulates the diffusion and migration of hazmat in the environment, so it can be adopted to evaluate the influence scope of hazmat transportation accidents. Zhao et al. [34] incorporated a Gaussian plume model into the box model to assess the air pollution area at facilities and routes. Yuan et al. [35] calculated the influence scope of the volatilization and diffusion of hazardous chemicals driven by wind with a Gaussian plume model. Zhao et al. [36] focused on three bodies of water, namely, rivers, lakes, and aquifers, and proposed a water diffusion model for each considering pipeline transportation accidents. The environmental bearing capacity is the natural purification capacity of the environment. Within the capacity limit, pollutants discharged into the environment may not cause harm to human health or the natural ecology through natural circulation [37]. Wang et al. [38] developed an evaluation index system for water environment carrying capacity (WECC) that considers the mutual interactions among six subsystems: industry, agriculture, population, water supply, water ecology, and water pollution. Peng et al. [39] built an indicator system of urban resources and environmental carrying capacity (URECC) based on ecological civilization, which included 18 indicators chosen from water, land, atmospheric environment, energy, and solid waste factors. Zhang et al. [40] included PM2.5 in their index system and established an evaluation system of regional resources and environmental carrying capacity (RECC). Thus, it is necessary to introduce the theory and method of environmental science into risk assessment of hazmat transportation. Moreover, most works in the literature consider the impact of hazmat accidents on a single environmental resource such as air or water, and seldom distinguish the impact of accidents on different environmental resources.
Furthermore, accident consequences may be affected by the bearing capacity of the assessed area. The bearing capacity of the assessed area includes the acceptable value of the population affected and the carrying capacity of the environment in the area. Within the capacity limit, the accident consequences may not cause serious casualties and harm to the natural ecology in the area. Although the accident consequences of a certain road segment are limited, they could be unacceptable if they exceed the bearing capacity range. In contrast, accident consequences could be significant, but could also be acceptable within the bearing capacity range. Hence, decision makers should appropriately process accident consequences according to the severity of the accident to better determine those road segments actually experiencing catastrophic consequences. More specifically, we set the population exposure risk preference parameter as higher than 1 if accident consequence exceeded the acceptable value of the population affected, or lower than 1 if it did not exceed it. We set the environmental pollution risk preference parameter as higher than 1 if accident consequence exceeded the carrying capacity of the environment, or lower than 1 if it did not exceed it.
Accounting for the above mentioned research gaps, this paper aims to develop a risk assessment model of hazmat road transportation under time-varying conditions. The main contributions of this paper can be summarized as follows:
(1)
A Gaussian plume model was adopted to simulate the dynamic areas at three levels of population exposure under time-varying conditions and utilized evacuation costs, inspection costs, medical costs and casualties to analyze the difference in the impact of accidents on populations in different areas.
(2)
We assessed the pollution scope of air, groundwater, lakes, and rivers with a variety of diffusion models under time-varying conditions and extended the pollution consequences to emergency disposal costs, monitoring costs and pollution remediation costs in order to scientifically evaluate the environmental pollution consequences of hazmat accidents.
(3)
Population exposure and environmental pollution risk preference parameters were set according to the severity of population exposure, considering the acceptable upper limit of population exposure and the seriousness of environmental pollution, accounting for the environmental bearing capacity. Then, population exposure and environmental risk assessment models were built under time-varying conditions.
Therefore, this paper improves the risk analysis system of hazmat transportation under time-varying conditions and provides a solid foundation for research on hazmat transportation location and path optimization.

3. Model

3.1. Leakage Accident Probability Model under Time-Varying Conditions

In this paper, Table 1 lists the weighted averages of truck accident probabilities in three major states determined by Harwood et al. [41]. We regard them as the initial truck accident probability P i j A . Table 2 shows the conditional leakage accident probability P i j B under the given accident types. Then, the product of P i j A and P i j B is applied to express the leakage accident probability P i j C [28]. In other words, P i j A and P i j B are multiplied to obtain P i j C .
Therefore, the leakage accident probability of hazmat road transportation is expressed in Equation (2). For the convenience of follow-up research, the unit is converted as indicated in Equation (3), where Pij(t) is the leakage accident probability in road segment (i, j) at time t; Mij(t) is the vehicle density of road segment (i, j) at time t, vehicle/km; and Lij is the driving length of road segment (i, j), km.
P i j C = P i j A × P i j B
P i j t = 10 6 × P i j A × P i j B × M i j t × L i j 2

3.2. Accident Consequence Assessment Model under Time-Varying Conditions

In this paper, the following assumptions apply when establishing the accident consequence assessment model: (1) hazmat immediately leaks after the vehicle overturns on the road, and the vehicle is regarded as a point risk source, where pollution spreads along all directions simultaneously into the external environment [2,34]; (2) the spilled hazmat does not degrade via the process of diffusion and spreads into air and water resources in equal amounts; (3) leakage occurs after an accident in a road segment, which does not affect the population and environment along other road segments.

3.2.1. Accident Consequence Model of Population Exposure

In this paper, a Gaussian plume model is applied to simulate the accident influence area of population exposure. The Gaussian plume model simulates the diffusion process of hazmat into air, so the impact area of air pollution simulated with this model is also applicable to determine the population exposed to air at any time.
Based on relevant research on the Gaussian plume model in references [34,35], the initial model is as follows:
C i j ( x , y , z , H , t ) = Q H a z a r d π v i j W i n d ( t ) σ y σ z × exp y 2 2 σ y 2 × exp ( z H ) 2 2 σ z 2 + exp ( z + H ) 2 2 σ z 2
where Cij(x, y, z, H, t) is the concentration of hazmat at downwind point (x, y, z) in road segment (i, j) at time t, mg/m3; x, y and z are the transversal distance, longitudinal distance and vertical distance, respectively, between the measuring point and leakage source in the case of downwind diffusion, m; H is the height of the leakage source, m; QHazard is the release speed of hazmat from the leakage source, mg/s; v i j W i n d ( t ) is the wind speed along road segment (i, j) at time t, m/s; σy is the horizontal diffusion coefficient, m; and σz is the vertical diffusion coefficient, m.
In this model, σy and σz mainly depend on the atmospheric stability parameters a, b, c, d and x [42], where σy = axc and σy = bxd. Suppose that y, z, and H are not considered, the initial model is modified as follows:
C i j ( x , t ) = Q H a z a r d π a b v i j W i n d ( t ) x c + d
According to the three injury threshold levels of the AEGL standard, three regions are divided, which are denoted by u = 1, 2, and 3. C u P o p is set as the limit value of the hazmat concentration in region u, mg/m3. The radius r i j u P o p ( t ) of population exposure due to the point leakage source in road segment (i, j) at time t is expressed in Equation (6), and the accident influence area S i j u P o p ( t ) of population exposure along road segment (i, j) is shown in Figure 1 and Equation (7), where r i j 4 P o p ( t ) = 0. S i j u P o p ( t ) is a circle or ring based on the above assumption (1).
r i j u P o p ( t ) = Q H a z a r d π a b v i j W i n d ( t ) C u P o p 1 c + d
S i j u P o p ( t ) = π Q H a z a r d π a b v i j W i n d ( t ) C u P o p 2 c + d Q H a z a r d π a b v i j W i n d ( t ) C u + 1 P o p 2 c + d
AEGL standards [43] were developed by the National Advisory Committee for Acute Exposure Guideline Levels for Hazmat (NAC/AEGL Committee) and recommended to the Environmental Protection Agency (EPA). In this paper, the AEGL standard is adopted as the C u P o p value, and Table 3 summarizes the classification of AEGL injury threshold levels.
Considering the differences in the population affected by accidents, this paper divides the costs of population exposure into four parts: evacuation costs, inspection costs, medical costs and losses associated with casualties. The four parts are expressed as F v P o p , where v = 1, 2, 3, and 4, including 1—evacuation costs; 2—inspection costs; 3—medical costs; and 4—losses associated with casualties. Evacuation costs can be converted into costs associated with the mobilization of medical, security police, traffic police and fire control departments. Inspection costs include internal and external examination costs, routine blood analysis costs, routine urine analysis costs, chest X-ray costs and electrocardiogram costs. Medical costs comprise the costs of registration, examination, laboratory testing, diagnosis, medication, surgery, hospitalization and rescue. The losses associated with casualties encompass relevant compensations.
Tij(t) is the population density along road segment (i, j) in area u at time t. The probability of cost type v generated by population exposure is αuv(t). Therefore, the total accident consequence of population exposure along road segment (i, j) under time-varying conditions is defined in Equation (8).
C i j P o p ( t ) = u = 1 3 v = 1 4 S i j u P o p ( t ) × T i j ( t ) × α u v ( t ) × F v P o p

3.2.2. Accident Consequence Model of Environmental Pollution

In this paper, the concentrations, influence lengths and areas of the various kinds of environmental pollution caused by accidents in road segment (i, j) at time t are expressed as C r E n v , r i j r E n v ( t ) and V i j r E n v ( t ) , respectively, where r = 1, 2, 3, and 4, for 1—air; 2—groundwater; 3—lakes; and 4—rivers.
We assessed the air pollution area. C 1 E n v is the limit value of the hazmat concentration in air (refers to the AEGL-2 standard). Based on Equation (6), the radius of the air pollution influence area of a given point leakage source in road segment (i, j) at time t is expressed in Equation (9). The accident influence area of air pollution along road segment (i, j) is shown in Figure 2 and Equation (10), which is a hemisphere [34].
r i j 1 E n v ( t ) = Q H a z a r d π a b v i j W i n d ( t ) C 1 E n v 1 c + d
V i j 1 E n v ( t ) = 2 3 π Q H a z a r d π a b v i j W i n d ( t ) C 1 E n v 3 c + d
Then, we simulated the water resource pollution area. In this paper, water resources are divided into groundwater and surface water, and surface water can be divided into static surface water (lakes) and dynamic surface water (rivers).
First, C 2 E n v is the limit value of the hazmat concentration in groundwater (refer to [45], the same below). Based on the groundwater diffusion model presented in [36], the influence radius of groundwater pollution due to the point leakage source in road segment (i, j) at time t is expressed in Equation (11), and the accident influence area (hemisphere) of groundwater along road segment (i, j) is shown in Figure 3 and Equation (12).
r i j 2 E n v ( t ) = v i j G w a t e r ( t ) t G w a t e r + H i j G w a t e r u t G w a t e r × ln Q C 2 E n v S i j G w a t e r π H i j G w a t e r u t G w a t e r
where v i j G w a t e r ( t ) is the average flow velocity of groundwater along road segment (i, j) at time t, m/s; tGwater is the time after hazmat release into groundwater, s; H i j G w a t e r is the average depth of groundwater along road segment (i, j), m; u* is the friction velocity of water flow [46], m/s, with the assumed value occurring in the interval of [0.02, 0.553]; Q is the release amount of hazmat into water, mg; and S i j G w a t e r is the wet cross-sectional area of groundwater along road segment (i, j), m2.
V i j 2 E n v ( t ) = 2 3 π v i j G w a t e r ( t ) t G w a t e r + H i j G w a t e r u t G w a t e r × ln Q C 2 E n v S i j G w a t e r π H i j G w a t e r u t G w a t e r 3
Second, C 3 E n v is the limit value of the hazmat concentration in lake water. According to the lake diffusion model reported in [36], the lake pollution influence radius of the point leakage source in road segment (i, j) at time t is defined in Equation (13).
r i j 3 E n v ( t ) = 0.58 H i j L a k e u t L a k e × ln Q 0.58 π H i j L a k e 2 u t L a k e C 3 E n v
where tLake is the diffusion time of hazmat into lake water, s; and H i j L a k e is the average depth of the lake along road segment (i, j), m.
Third, C 4 E n v is the limit value of the hazmat concentration in river water. Based on the river diffusion model proposed in [36], the influence radius of lake pollution of the point leakage source in road segment (i, j) at time t is expressed in Equation (14).
r i j 4 E n v ( t ) = v i j R i v e r ( t ) t R i v e r + 0.58 H i j R i v e r u t R i v e r × ln Q 2 + exp W i j R i v e r 2 0.145 H i j R i v e r u t R i v e r 0.58 π H i j R i v e r 2 u t R i v e r C 4 E n v
where v i j R i v e r ( t ) is the longitudinal velocity of river flow along road segment (i, j) at time t, m/s; t R i v e r is the time after hazmat release into the river, s; H i j R i v e r is the average depth of the river along road segment (i, j), m; and W i j R i v e r is the average width of the river along road segment (i, j), m.
Then, the accident influence area of lake and river pollution can be divided into two situations. One situation indicates that road segment (i, j) occurs along the lake or river, and the influence area (one-fourth of a sphere) is shown in Figure 4(left). The other situation states that road segment (i, j) passes through the lake or river, and the influence area (hemisphere) is shown in Figure 4(right) [36]. Then, the influence area (hemisphere) of lake and river pollution along road segment (i, j) is expressed in Equations (15) and (16), respectively.
V i j 3 E n v ( t ) = 1 3 o r 2 3 π 0.58 H i j L a k e u t L a k e × ln Q 0.58 π H i j L a k e 2 u t L a k e C 3 E n v 3 2
V i j 4 E n v ( t ) = 1 3 o r 2 3 π v i j R i v e r ( t ) t R i v e r + 0.58 H i j R i v e r u t R i v e r × ln Q 2 + exp W i j R i v e r 2 0.145 H i j R i v e r u t R i v e r 0.58 π H i j R i v e r 2 u t R i v e r C 4 E n v 3
In this paper, the pollution of air and water resources (groundwater, lakes and rivers) caused by accidents is comprehensively considered. According to the Recommended Method of Environmental Damage Appraisal and Assessment (version II) issued by the Ministry of Ecology and Environment of China (http://www.mee.gov.cn/gkml/hbb/bgt/201411/W020141105395741560668.pdf, accessed on 5 October 2020), the environmental pollution treatment costs associated with accidents in each road segment are divided into emergency disposal costs, monitoring costs and pollution remediation costs. These costs are expressed as F r w E n v , where r = 1, 2, 3, and 4, for 1—air pollution, 2—groundwater pollution, 3—lake pollution, and 4—river pollution; w = 1, 2, and 3, for 1—emergency disposal costs, 2—monitoring costs, and 3—pollution remediation costs. Emergency disposal costs include pollution control costs, on-site rescue costs, site cleaning costs, personnel transfer and resettlement costs, and emergency monitoring costs. Monitoring costs consist of the expenses for environmental monitoring, information disclosure, on-site investigation and supervision implementation. Pollution remediation costs comprise the expenses of scheme preparation and remediation and subsequent monitoring and supervision.
Therefore, the total consequence of environmental pollution in road segment (i, j) under time-varying conditions is as follows:
C i j E n v ( t ) = r = 1 4 w = 1 3 V i j r E n v ( t ) × F r w E n v

3.3. Transportation Risk Assessment Model under Time-Varying Conditions

3.3.1. Population Exposure Risk Assessment Model

In this paper, the risk preference parameter q i j P o p in the population exposure risk model is set according to the ratio p i j P o p of the actual population affected by the accident to the acceptable upper limit of the population affected along segment (i, j). This parameter is related to the severity of the accident impact on the population, which is set as follows:
p i j P o p = u = 1 3 S i j u P o p ( t ) × T i j ( t ) 10 3 × L i j × d P o p × T i j ( t )
where dPop is the width of the bearing area (referring to 600 m, as specified in the Emergency Response Guide [47]), m.
q i j P o p = 0.9 , 0 < p i j P o p 0.5 0.95 , 0.5 < p i j P o p < 1 1.0 p i j P o p = 1 1.05 1 < p i j P o p 1.5 1.1 p i j P o p > 1.5
Therefore, based on the leakage accident probability and accident consequence models of population exposure, the population exposure risk assessment model of hazmat road transportation under time-varying conditions is as follows:
E R i j P o p ( t ) = P i j ( t ) × C i j P o p ( t ) q i j P o p
where P i j t = 10 6 × P i j A × P i j B × M i j t × L i j 2
C i j P o p ( t ) = u = 1 3 v = 1 4 S i j u P o p ( t ) × T i j ( t ) × α u v ( t ) × F v P o p = v = 1 4 π × Q H a z a r d π a b v i j W i n d ( t ) C 1 P o p 2 c + d Q H a z a r d π a b v i j W i n d ( t ) C 2 P o p 2 c + d × α 1 v ( t ) × T i j ( t ) × F v P o p + π × Q H a z a r d π a b v i j W i n d ( t ) C 2 P o p 2 c + d Q H a z a r d π a b v i j W i n d ( t ) C 3 P o p 2 c + d × α 2 v ( t ) × T i j ( t ) × F v P o p + π × Q H a z a r d π a b v i j W i n d ( t ) C 3 P o p 2 c + d × α 3 v ( t ) × T i j ( t ) × F v P o p

3.3.2. Environmental Risk Assessment Model

In this paper, the risk preference parameter q i j E n v in the environmental risk model is introduced according to the ratio p i j E n v of the hazmat leakage amount in road segment (i, j) to the environmental bearing capacity of the area. This parameter is related to the severity of environmental pollution caused by the accident, as follows:
p i j E n v = Q r = 1 5 V i j r E n v ( t ) × C r E n v
q i j E n v = 0.9 , 0 < p i j E n v 0.5 0.95 , 0.5 < p i j E n v < 1 1.0 p i j E n v = 1 1.05 1 < p i j E n v 1.5 1.1 p i j E n v > 1.5
Therefore, based on the leakage accident probability model and accident consequence model of environmental pollution, the environmental risk assessment model of hazmat road transportation under time-varying conditions is expressed in Equation (23).
E R i j E n v ( t ) = P i j ( t ) × C i j E n v ( t ) q i j E n v
where P i j t = 10 6 × P i j A × P i j B × M i j t × L i j 2
C i j E n v ( t ) = r = 1 5 w = 1 3 V i j r E n v ( t ) × F r w E n v = w = 1 3 2 3 π Q H a z a r d π a b v i j W i n d ( t ) C 1 E n v 3 c + d × F 1 w E n v + 2 3 π v i j G w a t e r ( t ) t G w a t e r + H i j G w a t e r u t G w a t e r × ln Q C 2 E n v S i j G w a t e r π H i j G w a t e r u t G w a t e r 3 × F 2 w E n v + 1 3 o r 2 3 π 0.58 H i j L a k e u t L a k e × ln Q 0.58 π H i j L a k e 2 u t L a k e C 3 E n v 3 2 × F 3 w E n v + 1 3 o r 2 3 π v i j R i v e r ( t ) t R i v e r + 0.58 H i j R i v e r u t R i v e r × ln Q 2 + exp W i j R i v e r 2 0.145 H i j R i v e r u t R i v e r 0.58 π H i j R i v e r 2 u t R i v e r C 4 E n v 3 × F 4 w E n v

4. Case Study

Shanghai is one of the most highly populated metropolitan cities in the world, and it experiences heavy traffic, especially during the morning and afternoon rush hours. We apply the proposed models to the newly constructed real-world transportation network of the district of Songjiang, Shanghai. The network is simplified and contains 38 nodes and 52 arcs, as shown in Figure 5. We obtained the lengths of the arcs according to Baidu maps and retrieved the basic traffic flows along the arcs from the Shanghai Public Data Open Platform (https://data.sh.gov.cn/, accessed on 6 March 2021). The traffic flow divided by the length is the traffic density.
Population density information was acquired from the 2020 Shanghai Statistical Yearbook issued by the Municipal Bureau of Statistics in Shanghai. We collected wind speed data in different periods of Songjiang from the Shanghai Meteorological Service (http://sh.cma.gov.cn/, accessed on 7 March 2021). The costs of population exposure are based on the price standards of medical services in Shanghai (https://max.book118.com/html/2018/1009/5313144021001321.shtm, accessed on 7 March 2021) and the Interpretation of Several Issues Concerning the Application of Law in the Trial of Personal Injury Compensation Cases (http://www.court.gov.cn/fabu-xiangqing-282621.html, accessed on 7 March 2021). The calculation of the environmental pollution treatment costs refers to the instructions in the Recommended Method of Environmental Damage Appraisal and Assessment (http://www.mee.gov.cn/gkml/hbb/bgt/201411/W020141105395741560668.pdf, accessed on 8 March 2021).
A truck transporting 6 tons of ammonia (gaseous) is considered in the network. The results of the case study are analyzed.

4.1. Analysis of the Accident Probability

Figure 6 shows the accident probabilities during 24 periods in road segments (6, 8), (17, 18), (25, 26) and (33, 34). It is clear that the accident probabilities of hazmat transportation are low, and the accident probabilities increase during the morning and afternoon rush hours but decrease during other periods (especially late at night).
For the same road segment, only the vehicle density Mij(t) in Equation (3) is affected by time-varying conditions. As vehicle densities on roads tend to increase during the morning and afternoon rush hours (Figure 7), the accident probabilities also exhibit a similar trend. Hence, we can effectively reduce accident probability by controlling vehicle density in practice.

4.2. Analysis of the Accident Consequences

We determine population exposure accident consequences and wind speeds during 24 periods in road segments (15, 16) and (30, 31), as shown in Figure 8 and Figure 9, respectively. It is observed that there are significant population exposure accident consequences during hazmat transportation. In addition, for the same road segment, if αuv(t) remains fixed, only the wind speed v i j W i n d ( t ) and population density Tij(t) in Equation (7) are affected by time-varying conditions.
Although we simulate changes in population density in residential areas (population density is lower during the daytime and higher at night), commercial office areas (population density is higher during the daytime and lower at night), and boarding schools (population density during the daytime and at night is similar), the impact on population exposure accident consequences is limited and irregular.
However, Figure 8 and Figure 9 show that wind speed is inversely proportional to accident consequences. In other words, the higher the wind speed, the less significant the accident consequence is. Suppose that the diffusion time is not considered, and the amount of leaked hazmat remains the same; the higher the wind speed, the faster the diffusion process of hazmat occurs, which leads to the phenomenon of hazmat dilution. Because we set the hazmat concentration (referring to the AEGL standard) in advance, the influence area reaching the specified concentration is smaller, and the accident consequence is less remarkable. In practice, we should consider wind speed to avoid any catastrophic accident consequences for population exposure.
Then, we determine the accident influence areas of air, groundwater, lake and river pollution during 24 periods along road segment (17, 18), as shown in Figure 10. The variation in the accident influence area of environmental pollution follows the order of air > river > groundwater > lake. Under time-varying conditions, only the wind speed v i j W i n d ( t ) in Equation (10) is affected. Equation (15) indicates that lakes remain motionless, so the accident influence area of lake pollution is independent of time. Then, under time-varying conditions, only the wind speed v i j W i n d ( t ) in Equation (10) is affected. Additionally, in Equations (12) and (16), v i j G w a t e r ( t ) and v i j R i v e r ( t ) are affected by time-varying conditions, but they are much lower than the wind speed, and the former is lower than the latter. The lower the speed, the smaller the variation difference is. Thus, the variation area of the accident influence area of groundwater pollution is smaller than that of river pollution and is obviously smaller than that of air pollution. Therefore, we should pay attention to the impact of time-varying conditions on the influence area of air pollution in practice.
In addition, choosing road segment (32, 33) as an example, the accident consequences during 24 periods are shown in Figure 11 and Figure 12. The accident consequences of environmental pollution are similar to those of population exposure, and they are also quite notable. Hence, hazmat transportation accidents exhibit the characteristics of low probability and significant consequences.
In terms of population exposure and environmental pollution accident consequences, based on Equations (18), (19), (21) and (22), we set the population exposure risk preference parameter q i j P o p as higher than 1 if accident consequence exceeded the acceptable value of the population affected, or lower than 1 if it did not exceed it; then, we set the environmental pollution risk preference parameter q i j E n v as higher than 1 if accident consequence exceeded the carrying capacity of the environment, or lower than 1 if it did not exceed it. This extends the consequences beyond the bearing capacity or reduces the consequences within the bearing capacity, respectively, as shown in Figure 11 and Figure 12, which clearly reveal the period when the accident actually causes catastrophic consequences. In addition, if we analyze the accident consequences of each road segment, the insights are the same.

4.3. Analysis of the Transportation Risk

Comparing Figure 12 and Figure 13, the accident consequences considering the risk preference parameter and risk reveal the same trend during each period, and the corresponding transportation risks in the road segments with major accident consequences are higher. Because the risk preference parameter enhances the accident consequences to encompass the periods and road segments with serious accidents, this avoids the problem that certain periods and road segments with catastrophic accident consequences are overlooked due to a low probability when calculating risk.
When planning transportation routes or schemes, we should avoid segments with a higher transportation risk, and at the same time, we should reasonably avoid road segments with more notable accident consequences according to the characteristics of low probability and significant consequences in hazmat transportation accidents. In other words, we should consider transportation risk and avoid catastrophic consequences to identify road segments with a high risk of catastrophic consequences and improve the security of hazmat transportation.

5. Conclusions

Road accidents involving hazmat transportation have caused notable casualties and excessive environmental pollution. Therefore, it is necessary to avoid high-risk road segments and prevent catastrophic accidents through in-depth research and assessment of hazmat transportation risk. Based on the probability of leakage accidents, this paper explores transportation risk in two aspects of population exposure and environmental pollution and establishes a transportation risk assessment model under time-varying conditions. This compensates for certain shortcomings of previous models and further improves the risk assessment system of hazmat transportation under time-varying conditions.
Therefore, this paper not only promotes the development of the risk analysis theory, technology and methods of hazmat transportation to a certain extent, but also provides support for the scientific decision-making process in regard to hazmat road transportation locations and path optimization. More specifically, the results retrieved from the proposed model could be used as basics for hazmat transportation location and path optimization to help decision makers scientifically formulate decisions regarding hazmat transportation location and path optimization.
In addition, through the case study, it was clearly revealed that in practice, controlling certain parameters in the transportation risk assessment model could reduce the probabilities and consequences of hazmat road transportation accidents. We should also avoid road segments with a high risk of catastrophic accident consequences to improve the safety of hazmat transportation.
Several directions of research should be explored in the future. This paper should be applied to a larger case to test and improve the theoretical methods, and risk assessment of hazmat multimodal transportation under time-varying conditions should be further studied to more closely reflect reality.

Author Contributions

Conceptualization, L.L. and Q.W.; methodology, L.L. and Q.W.; validation, S.L., Y.L. and T.F.; data curation, Q.W.; writing—original draft preparation, Q.W.; writing—review and editing, L.L.; visualization, Q.W.; supervision, T.F.; project administration, S.L. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is funded by the National Natural Science Foundation of China (72032001, 72074076, 71302043); Ministry of Education, Humanities and Social Sciences Research Planning Foundation (21YJA630057); and Key Reserch and Development Projects of Shanghai Science and Technology Commission(21DZ1204104).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Xin, C.L.; Meng, J.; Zhang, J.W.; Feng, Q.R. Review on the research of hazardous materials transportation route optimization. Chin. J. Saf. Sci. 2018, 28, 102–107. [Google Scholar]
  2. Zhao, J.H. Study on Location Routing Problem of Hazardous Waste Recycling System Considering Environmental Risk; Chinese Southwest Jiaotong University: Chengdu, China, 2015. [Google Scholar]
  3. Abkowitz, M.; Lepofsky, M.; Cheng, P. Selecting criteria for designating hazardous materials road routes. Transp. Res. Rec. J. Transp. Res. Board 1992, 1333, 30–35. [Google Scholar]
  4. Sivakumar, R.A.; Batta, R.; Karwan, M.H. A multiple route conditional risk model for transporting hazardous materials. INFOR Inf. Syst. Oper. Res. 1995, 33, 20–33. [Google Scholar] [CrossRef]
  5. Erkut, E.; Ingolfsson, A. Catastrophe avoidance models for hazardous materials route planning. Transp. Sci. 2000, 34, 165–179. [Google Scholar] [CrossRef]
  6. Kang, Y.Y. Value-at-Risk Models for Hazardous Material Transportation; State University of New York at Buffalo: New York, NY, USA, 2011. [Google Scholar]
  7. Kang, Y.Y.; Batta, R.; Kwon, C. Value-at-Risk model for hazardous material transportation. Ann. Oper. Res. 2014, 222, 361–387. [Google Scholar] [CrossRef]
  8. Toumazi, I.; Kwon, C. Routing hazardous materials on time-dependent networks using conditional value-at-risk. Transp. Res. Part C Emerg. Technol. 2013, 37, 73–92. [Google Scholar] [CrossRef] [Green Version]
  9. Cao, H.; Fan, J.; Liu, L.P. Research on route selection of road rail combined transport of hazardous materials based on CVaR. Oper. Res. Manag. 2017, 26, 41–48. [Google Scholar]
  10. Wei, H.; Li, J. Population exposure risk analysis of hazardous materials transportation under time varying conditions. Chin. J. Saf. Sci. 2004, 10, 98–101+2. [Google Scholar]
  11. Hossain, M.; Muromachi, Y. A Bayesian network based framework for real-time crash prediction on the basic freeway segments of urban expressways. Accid. Anal. Prev. 2012, 45, 373–381. [Google Scholar] [CrossRef]
  12. Yu, R.; Abdel-Aty, M. Utilizing support vector machine in real-time crash risk evaluation. Accid. Anal. Prev. 2013, 51, 252–259. [Google Scholar] [CrossRef]
  13. Sun, J.; Sun, J. A dynamic Bayesian network model for real-time crash prediction using traffic speed conditions data. Transp. Res. Part C Emerg. Technol. 2015, 54, 176–186. [Google Scholar] [CrossRef]
  14. Wang, L.; Abdel-Aty, M.; Shi, Q.; Park, J. Real-time crash prediction for expressway weaving segments. Transp. Res. Part C Emerg. Technol. 2015, 61, 1–10. [Google Scholar] [CrossRef]
  15. Xu, C.C.; Wang, Y.; Liu, P.; Wang, W.; Bao, J. Quantitative risk assessment of freeway crash casualty using high-resolution traffic data. Reliab. Eng. Syst. Saf. 2018, 169, 299–311. [Google Scholar] [CrossRef]
  16. Yang, K.; Wang, X.; Yu, R. A Bayesian dynamic updating approach for urban expressway real-time crash risk evaluation. Transp. Res. Part C Emerg. Technol. 2018, 96, 192–207. [Google Scholar] [CrossRef]
  17. Yu, J.; Wang, Y.Y.; Zou, Z.H.; Wang, L.Q. Convolutional neural networks with refined loss functions for the real-time crash risk analysis. Transp. Res. Part C Emerg. Technol. 2020, 119, 102740. [Google Scholar] [CrossRef]
  18. Li, P.; Abdel-Aty, M.; Yuan, J. Real-time crash risk prediction on arterials based on LSTM-CNN. Accid. Anal. Prev. 2020, 135, 105371. [Google Scholar] [CrossRef] [PubMed]
  19. Holeczek, N. Hazardous materials truck transportation problems: A classification and state of the art literature review. Transp. Res. Part D Transp. Environ. 2019, 69, 305–328. [Google Scholar] [CrossRef]
  20. Fabiano, B.; Currò, F.; Palazzi, E.; Pastorino, R. A framework for risk assessment and decision-making strategies in dangerous good transportation. J. Hazard. Mater. 2002, 93, 1–15. [Google Scholar] [CrossRef]
  21. Zhang, J.H.; Zhao, L.J. Risk analysis of hazardous chemicals transportation. Theory Pract. Syst. Eng. 2007, 27, 117–122. [Google Scholar] [CrossRef]
  22. Fan, T.J.; Chiang, W.C.; Russell, R. Modeling urban hazmat transportation with road closure consideration. Transp. Res. Part D Transp. Environ. 2015, 35, 104–115. [Google Scholar] [CrossRef]
  23. Zero, L.; Bersani, C.; Paolucci, M.; Sacile, R. Two new approaches for the bi-objective shortest path with a fuzzy objective applied to HAZMAT transportation. J. Hazard. Mater. 2019, 375, 96–106. [Google Scholar] [CrossRef] [PubMed]
  24. Dong, S.; Zhou, J.B.; Ma, C.X. Design of a network optimization platform for the multivehicle transportation of hazardous materials. Int. J. Environ. Res. Public Health 2020, 17, 1104. [Google Scholar] [CrossRef] [Green Version]
  25. Sun, Y. A Fuzzy Multi-objective routing model for managing hazardous materials door-to-door transportation in the road-rail multimodal network with uncertain demand and improved service level. IEEE Access 2020, 8, 172808–172828. [Google Scholar] [CrossRef]
  26. Khanmohamadi, M.; Bagheri, M.; Khademi, N.; Ghannadpour, S.F. A security vulnerability analysis model for dangerous goods transportation by rail-Case study: Chlorine transportation in Texas-Illinois. Saf. Sci. 2018, 110, 230–241. [Google Scholar] [CrossRef]
  27. Orozco, J.L.; Van Caneghem, J.; Hens, L.; Gonzalez, L.; Lugo, R.; Diaz, S.; Pedroso, I. Assessment of an ammonia incident in the industrial area of Matanzas. J. Clean. Prod. 2019, 222, 934–941. [Google Scholar] [CrossRef]
  28. Ren, C.X. Research on Route Optimization Method of Dangerous Goods Road Transportation Based on Risk Analysis; Nankai University: Tianjin, China, 2007. [Google Scholar]
  29. Xia, Q.; Qian, Y.; Liu, M.F. Optimal route selection of dangerous goods road transportation based on environmental risk assessment—A case study of Zhangjiagang City. Chin. Environ. Sci. 2014, 34, 266–272. [Google Scholar]
  30. Saat, M.R.; Werth, C.J.; Schaeffer, D.; Yoon, H.; Barkan, C.P.L. Environmental risk analysis of hazardous material rail transportation. J. Hazard. Mater. 2014, 264, 560–569. [Google Scholar] [CrossRef] [PubMed]
  31. Cordeiro, F.G.; Bezerra, B.S.; Peixoto, A.S.P.; Ramos, R.A.R. Methodological aspects for modeling the environmental risk of transporting hazardous materials by road. Transp. Res. Part D Transp. Environ. 2016, 44, 105–121. [Google Scholar] [CrossRef]
  32. Huang, X.F.; Wang, X.H.; Pei, J.J.; Xu, M.; Luo, X.W. Risk assessment of the areas along the road due to hazardous material transportation accidents. Nat. Hazards 2018, 93, 1181–1202. [Google Scholar] [CrossRef]
  33. Wang, H.X.; Liang, Q.Q. Risk analysis and route optimization of dangerous goods transportation based on the empirical path set. J. Adv. Transp. 2020, 2020, 1–13. [Google Scholar] [CrossRef]
  34. Zhao, J.H.; Verter, V. A bi-objective model for the used oil location-routing problem. Comput. Oper. Res. 2014, 62, 157–168. [Google Scholar] [CrossRef]
  35. Yuan, W.Y.; Peng, Y.; Yang, F.M. Two level emergency center location model based on hazardous chemical accidents. Theory Pract. Syst. Eng. 2015, 35, 728–735. [Google Scholar]
  36. Zhao, J.H.; Verma, M.; Verter, V. Pipeline transportation of crude oil in Canada: Environmental risk assessment using modified diffusion models. Hum. Ecol. Risk Assess. Int. J. 2020, 27, 1206–1226. [Google Scholar] [CrossRef]
  37. Wang, H.D.; Xia, Q. Research progress of environmental capacity. Environ. Sci. Technol. 1983, 01, 32–36. [Google Scholar]
  38. Wang, Y.M.; Zhou, X.D.; Engel, B. Water environment carrying capacity in Bosten Lake basin. J. Clean. Prod. 2018, 199, 574–583. [Google Scholar] [CrossRef]
  39. Peng, T.; Deng, H. Evaluating urban resource and environment carrying capacity by using an innovative indicator system based on eco-civilization—A case study of Guiyang. Environ. Sci. Pollut. Res. 2020, 28, 6941–6955. [Google Scholar] [CrossRef]
  40. Zhang, G.; Luo, S.; Jing, Z.; We, S.; Ma, Y. Evaluation and forewarning management of regional resources and environment carrying capacity: A case study of Hefei city, Anhui province, China. Sustainability 2020, 12, 1637. [Google Scholar] [CrossRef] [Green Version]
  41. Harwood, D.W.; Viner, J.G.; Russell, E.R. Procedure for developing truck accident and release rates for HAZMAT routing. J. Transp. Eng. 1993, 119, 189–199. [Google Scholar] [CrossRef]
  42. Arya, S.P. Air Pollution Meteorology and Dispersion; Oxford University Press: Cambridge, UK, 1999. [Google Scholar]
  43. Krewski, D.; Bakshi, K.; Garrett, R.; Falke, E.; Rusch, G.; Gaylor, D. Development of acute exposure guideline levels for airborne exposures to hazardous substances. Regul. Toxicol. Pharmacol. 2004, 39, 184–201. [Google Scholar] [CrossRef]
  44. National Research Council. Impact of Revised Airborne Exposure Limits on Non-Stockpile Chemical Materiel Program Activities; National Academies Press: Washington, WA, USA, 2005. [Google Scholar]
  45. Hu, P.; Shuai, B.; Zhao, J.H. Environmental risk analysis and measurement of volatile hazardous waste road transportation. Ind. Saf. Environ. Prot. 2014, 40, 8–12. [Google Scholar]
  46. Seo, I.W.; Cheong, T.S. Predicting longitudinal dispersion coefficient in natural streams. J. Hydraul. Eng. 1998, 124, 25–32. [Google Scholar] [CrossRef]
  47. Alrukaibi, F.; Alrukaibi, D.; Alkheder, S.; Alojaiman, S.; Sayed, T. Optimal route risk-based algorithm for hazardous material transport in Kuwait. J. Loss Prev. Process. Ind. 2018, 52, 40–53. [Google Scholar] [CrossRef]
Figure 1. Accident influence area of population exposure.
Figure 1. Accident influence area of population exposure.
Ijerph 18 09780 g001
Figure 2. Accident influence area of air pollution.
Figure 2. Accident influence area of air pollution.
Ijerph 18 09780 g002
Figure 3. Accident influence area of groundwater pollution.
Figure 3. Accident influence area of groundwater pollution.
Ijerph 18 09780 g003
Figure 4. Accident influence area of lake and river pollution.
Figure 4. Accident influence area of lake and river pollution.
Ijerph 18 09780 g004
Figure 5. Road network in Songjiang District (simplified).
Figure 5. Road network in Songjiang District (simplified).
Ijerph 18 09780 g005
Figure 6. Accident probabilities.
Figure 6. Accident probabilities.
Ijerph 18 09780 g006
Figure 7. Vehicle densities on the roads.
Figure 7. Vehicle densities on the roads.
Ijerph 18 09780 g007
Figure 8. Population exposure accident consequences and wind speed in road segment (15, 16).
Figure 8. Population exposure accident consequences and wind speed in road segment (15, 16).
Ijerph 18 09780 g008
Figure 9. Population exposure accident consequences and wind speed in road segment (30, 31).
Figure 9. Population exposure accident consequences and wind speed in road segment (30, 31).
Ijerph 18 09780 g009
Figure 10. Environmental pollution accident influence areas.
Figure 10. Environmental pollution accident influence areas.
Ijerph 18 09780 g010
Figure 11. Population exposure accident consequences.
Figure 11. Population exposure accident consequences.
Ijerph 18 09780 g011
Figure 12. Environmental pollution accident consequences.
Figure 12. Environmental pollution accident consequences.
Ijerph 18 09780 g012
Figure 13. Transportation risk.
Figure 13. Transportation risk.
Ijerph 18 09780 g013
Table 1. Weighted average truck accident probability (accidents per million vehicle/km) in three states of the United States.
Table 1. Weighted average truck accident probability (accidents per million vehicle/km) in three states of the United States.
RegionRoad Type P i j A
CityTwo-lane road5.38
Undivided multilane road8.65
Divided multilane road7.75
Cycle lane6.03
Expressway1.35
Table 2. Leakage accident probability under certain accidents.
Table 2. Leakage accident probability under certain accidents.
Accident Type P i j B
Single-vehicle noncollision accidentLeaves the road0.331
Overturned on the road0.375
Other noncollision accident0.169
Single-vehicle collision accidentCollision with stopped vehicle0.031
Collision with a train0.455
Collision with a nonmotorized vehicle0.015
Collision with fixed objects0.129
Other collision accident0.059
Table 3. Classification of the AEGL injury threshold levels.
Table 3. Classification of the AEGL injury threshold levels.
GradeDefinition
AEGL-1The airborne concentration of a substance above which it is predicted that the general population, including susceptible individuals, could experience notable discomfort, irritation, or certain asymptomatic nonsensory effects. However, the effects are not disabling and are transient and reversible upon cessation of exposure.
AEGL-2The airborne concentration of a substance above which it is predicted that the general population, including susceptible individuals, could experience irreversible or other serious, long-lasting adverse health effects or an impaired ability to escape.
AEGL-3The airborne concentration of a substance above which it is predicted that the general population, including susceptible individuals, could experience life-threatening health effects or death.
Note: the information in Table 3 refers to [44].
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liu, L.; Wu, Q.; Li, S.; Li, Y.; Fan, T. Risk Assessment of Hazmat Road Transportation Considering Environmental Risk under Time-Varying Conditions. Int. J. Environ. Res. Public Health 2021, 18, 9780. https://doi.org/10.3390/ijerph18189780

AMA Style

Liu L, Wu Q, Li S, Li Y, Fan T. Risk Assessment of Hazmat Road Transportation Considering Environmental Risk under Time-Varying Conditions. International Journal of Environmental Research and Public Health. 2021; 18(18):9780. https://doi.org/10.3390/ijerph18189780

Chicago/Turabian Style

Liu, Liping, Qing Wu, Shuxia Li, Ying Li, and Tijun Fan. 2021. "Risk Assessment of Hazmat Road Transportation Considering Environmental Risk under Time-Varying Conditions" International Journal of Environmental Research and Public Health 18, no. 18: 9780. https://doi.org/10.3390/ijerph18189780

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