Anthropogenic atmospheric emissions of typical toxic heavy
metals have caused worldwide concern due to their adverse effects on
human health and the ecosystem. By determining the best available
representation of time-varying emission factors with S-shape curves, we
establish the multiyear comprehensive atmospheric emission inventories of 12
typical toxic heavy metals (Hg, As, Se, Pb, Cd, Cr, Ni, Sb, Mn, Co, Cu, and
Zn) from primary anthropogenic activities in China for the period of
1949–2012 for the first time. Further, we allocate the annual emissions of
these heavy metals in 2010 at a high spatial resolution of 0.5
Heavy metals (HMs) is a general collective term which applies to the group
of metals (e.g., Hg, Pb, Cd, Cr, Ni, Sb, Mn, Co, Cu, Zn) and metalloids
(e.g., As, Se) with atomic density greater than 4.5 g cm
Since the 1980s, the United States, the United Kingdom, Australia and several other developed countries have begun to compile their national emission inventories of varied hazardous air pollutants (including HMs), such as the US National Emission Inventory (NEI), the UK National Atmospheric Emission Inventory (NAEI), and the Australian National Pollutant Inventory (NPI). Further, the quantitative assessments of global contamination of air by HMs from anthropogenic sources have been estimated in previous studies (Nriagu, 1979; Nriagu and Pacyna, 1988; Pacyna and Pacyna, 2001; Streets et al., 2011; Tian et al., 2014b). With the increasing environmental pollution associated with economic growth, some researchers have paid special attention to estimating China's HM emission inventory, especially for Hg, which is regarded as a global pollutant (Fang et al., 2002; Streets et al., 2005; Wu et al., 2006). Streets et al. (2005) and Wu et al. (2006) developed Hg emission inventories from anthropogenic activities of China for the year 1999 and 1995 to 2003, respectively. A research group led by Tian established the integrated emission inventories of eight HMs (Hg, As, Se, Pb, Cd, Cr, Ni and Sb) from coal combustion or primary anthropogenic sources on the provincial level during 1980 to 2009 (Cheng et al., 2015; Tian et al., 2010, 2012a–c, 2014a, b). However, comprehensive and detailed studies on anthropogenic atmospheric emissions of 12 typical toxic HMs with highly resolved temporal and spatial distribution information in China are still quite limited. Moreover, we have little knowledge on what the past and accelerated emission levels of HMs are like from anthropogenic sources during the historical period since the founding of the People's Republic of China in 1949, to the open-poor policy (1978).
In this study, for the first time, we have evaluated the historical trend and spatial distribution characteristics by source categories and provinces of atmospheric emissions of 12 typical HMs (Hg, As, Se, Pb, Cd, Cr, Ni, Sb, Mn, Co, Cu, and Zn) from primary anthropogenic activities during the period of 1949–2012. Particularly, we have attempted to determine the temporal variation profiles of emission factors for several significant sources categories (e.g., nonferrous metal smelting, ferrous metal smelting, cement production, and municipal solid waste (MSW) incineration) during the long period of 1949 to 2012, which includes the technological upgrade of the industrial process and the progress of the application rate for various air pollutant control devices (APCDs).
We estimate the atmospheric emissions of the 12 target HMs (Hg, As, Se, Pb, Cd, Cr, Ni, Sb, Mn, Co, Cu, and Zn) from primary anthropogenic sources by combining the specific annual activities and dynamic emission factors by source category in this study. Table S1 in the Supplement lists the target heavy metal species and the associated emission sources. Generally, we classify all sources into two major categories: coal combustion sources and non-coal combustion sources.
Currently, coal plays a dominant role in China's energy consumption, making up about 70 percent of its total primary energy consumption (Tian et al., 2007, 2012b). Consequently, tons of hazardous HM pollutants can be released into the atmospheric environment, although the concentration of heavy metals in Chinese coals is normally at parts per million (ppm) levels.
Atmospheric emissions of varied HMs from coal combustion are calculated by
combining the provincial average concentration of each heavy metal in feed
coals, the detailed coal consumption data, and the specific emission
factors, which are further classified into subcategories with respect to
different boiler configurations and the application rates and removal
efficiencies of various APCDs. The basic formulas can be expressed as
follows:
Previous studies have demonstrated that concentrations of HMs in Chinese coals vary substantially depending on the type of the feed coals and their origin, as well as the affinity of the particular element for pure coal and mineral materials (Tang et al., 2002; Ren et al., 2006).
In this study, we compile and summarize provincial-level test data of HM content in coal from published literature to date: Hg (879 samples), As (1018 samples), Se (472 samples), Pb (831 samples), Cd (616 samples), Cr (956 samples), Ni (863 samples), Sb (1612 samples), Mn (545 samples), Co (888 samples), Cu (765 samples), and Zn (828 samples), and then we calculate the average concentration of each heavy metal in coal as produced and coal as consumed on a provincial level by using bootstrap simulation and a coal transmission matrix (Tian et al., 2011a, 2013, 2014a). More details about the algorithms to determine HM concentrations in cleaned coals, briquettes and coke are given in our previous publications (Tian et al., 2010, 2012a). The brief introduction of bootstrap simulation as well as averaged concentration values of Hg, As, Se, Pb, Cd, Cr, Ni, Sb, Mn, Co, Cu, and Zn in feed coals on the provincial level can be seen in the Supplement, Sects. S1–S2 and Table S3–S9.
In this study, various coal combustion facilities are separated into five subcategories: pulverized-coal boilers, stoker-fired boilers, fluidized-bed furnaces, coke furnaces, and domestic coal-fired stoves. Therein, pulverized-coal boilers are predominant in coal-fired power plants in most of the provinces in China, representing over 85.0 % of the total installed capacities. The remaining share is divided between fluidized-bed furnaces and stoker-fired boilers, which are mainly used in relatively small unit-size coal-fired power plants. Different from the thermal power plants sector, stoker-fired boilers make up a large proportion of the coal-fired industrial sector and other commercial coal-fired sectors. The release rates of HMs in flue gas from various boiler categories vary substantially due to the different combustion patterns and operating conditions, as well as their genetic physical and chemical characteristics (Reddy et al., 2005). Therefore, it is necessary to develop a detailed specification of the methods by which the coals are fed and burned in China. In this study, we have compiled the release rates of Hg, As, Se, Pb, Cd, Cr, Ni, Sb, Mn, Co, Cu, and Zn for different combustion boilers from published literature (see Supplement Table S10). The arithmetic mean values of release rates of these 12 HMs from different combustion boilers reported in the literature are adopted to calculate the final emissions (see Table 1).
Averaged release rates and removal efficiencies of various HMs from coal-fired facilities and the installed APCDs.
Furthermore, the conventional APCDs used to reduce criteria air pollutants
(e.g., PM, SO
The residential sector is another important coal consumer in China. The
traditional cook stoves and improved cook stoves are major combustion
facilities for residential cooking and heating, both of which normally
do not have any PM and SO
HM emissions from non-coal categories are calculated as a product of annual
activity data (e.g., fuel consumption, industrial product yields) and
specific emission factors of varied HMs. The basic calculation can be
described by the following equation:
Notably, atmospheric emissions of Pb have significantly dropped in China, as
a result of unleaded gasoline introduction since the early 2000s. The
proportion of lead in leaded gasoline emitted to the air is estimated at
about 77 % (Biggins and Harrison, 1979) or 75 % (Hassel et al., 1987),
and thus this parameter is assumed to be at about 76 % for the period before
2000 in this study. Consequently, for leaded gasoline used by motor vehicles
in China, the total Pb emitted to the atmosphere is calculated according to
the following equation:
For brake and tyre wear, the atmospheric emissions of several HMs are
estimated by the following equation:
Because remarkable changes in products, devices, processes as well as practices (technology improvement) have imposed positive effects on emission reductions of pollutants with the growth of the economy and the increasing awareness of environment protection, the resulting pollutant emission level at any given time is a competition between technology improvement and production growth. Consequently, one of the major challenges in this study is to develop a reasonable representation of the time-varying dynamic emission factors of HMs associated with each primary industrial activity.
Taking into consideration the updated air pollutant control technologies,
and outdated enterprises that have shut down, the HM emission factors show a
gradually declining trend. Generally, the patterns of technologies' diffusion
across competitive markets are evident, and an S-shaped curve is a typical
result when plotting the proportion of a useful service or product supplied
by each major competing technology (Grübler et al., 1999). At the
earliest stage of industrialization, growth rates in the removal efficiency
of an air pollutant are slow as the advanced technology with high investment
and operation costs is only applied in specialized niche sectors.
Subsequently, along with the progress in technology and awareness of public
environmental protection, growth rates accelerate as early commercial
investments have resulted in standard-setting and compounding cost
reductions, which leads to the increased application of advanced technologies
for reduction in emissions of air pollutants in a wider array of settings.
Eventually, growth rate in the removal efficiency will gradually approach
near to zero as the potential market of optimal control technology of HM
emissions is saturated. By using an S-shaped curve, both historical and
future emissions of carbon aerosol and Hg to the atmosphere from human
activities have been evaluated by Bond et al. (2007) and Streets et al. (2004, 2011). Their results show that an S-shaped curve fits historical and
future trends better than polynomial or linear fits, even though it cannot
account for economic shocks because of the form of monotonous smooth
transitions. Therefore, S-shaped curves are applied to estimate the dynamic
HM emission factors from primary industrial process sources in this study.
The basic formulas can be expressed as follows:
Based on the above method, we build the dynamic representation of HM
emission factors to reflect the transition from uncontrolled processes
pre-1900 to the relatively high-efficiency abatement processes in 2012.
Parameters for some of these transitions are discussed throughout the paper,
and are summarized in the Supplement, Table S14. Actually, on the basis of Eq. (5), the specific values of the shape parameter of the curve (
By 2010, bath smelting (e.g., Ausmelt smelting, Isa smelting), flash
smelting, and imperial smelting process (ISP) represent the three most
commonly used techniques for copper smelting, for about 52, 34, and
10 % of Chinese copper production, respectively. For lead smelting,
sintering plus a blast furnace technique (traditional technique) and bath
smelting (e.g., oxygen side blowing, oxygen bottom blowing) plus a blast
furnace technique (advanced technique) are the two most commonly used
techniques in China, accounting for about 48 and 47 % of lead
production, respectively. With respect to zinc smelting, hydrometallurgy is
the predominant technique in China, for about 77 % of the zinc production
capacity. The remaining share is divided among vertical retort (VR)
pyrometallurgy (
Because of limited information and lack of field experimental tests on HM emissions in these source categories in China, some emission factors for this source category are cited from published literature, with only nationally averaged levels. Streets et al. (2011) indicate that China, eastern Europe and the former USSR can be regarded as a uniform region with similar levels of technology development, whose emission factor trajectories are identical. Therefore, we presume the emission factors of HMs with higher abatement implementation in eastern Europe, Caucasus and central Asian countries are equivalent to those in China for the same calendar year (see Fig. 1a). Based on the above assumptions and default abatement efficiencies of HMs in nonferrous metal smelting sectors (EEA, 2013), as well as other specific emission factors of HMs from published literature to date (Nriagu, 1979; Pacyna, 1984; Pacyna and Pacyna, 2001), the unabated emission factors are determined (see Supplement Tables S15–S16).
Time variation of arsenic emission factors for copper production and zinc emission factors for steel production in China (an example).
Presently, compared to primary smelting of Cu, Pb and Zn, there is much less information about emission factors of HMs for secondary metal smelting of Cu, Pb and Zn and other nonferrous metals (Al, Ni and Sb) smelting from the published literature. Hence, it is much more difficult to estimate the time-varying dynamic emission factors of HMs from the above sectors by the use of S-shaped curves due to a lack of necessary baseline information. We presume that the average emission factors for secondary metal (Cu, Pb and Zn) smelting, aluminum smelting, antimony smelting, and nickel smelting remain unchanged before the year 1996, at which the Emission Standard of Pollutants for Industrial Kiln and Furnace was first issued in China. We also presume that the average emission factors of HMs from secondary metal smelting and other nonferrous metal smelting for developing countries, referred to in Pacyna and Pacyna's report (2001) and eastern Europe, Caucasus and central Asian countries with limited abatement referred to in the EEA Guidebook (EEA, 2009; S. V. Kakareka, personal communication, 2008), as well as the United Kingdom applied in the NAEI emission factors database pre-1990 (UK, 1995), are reasonable for China before the year 1996. Subsequently, atmospheric emission factors of HMs from the industry of nonferrous metals in China decrease gradually with the implementation of tightened emission limits regulated by the gradually stricter Emission Standards of Pollutants from nonferrous metals industry (e.g., GB 9078–1996, GB 25465–2010, GB 25466–2010, GB 25467–2010).
With respect to gold smelting (large-scale) and mercury mining industries, the time-varying Hg emission factors from these two subcategories are determined by referring to studies carried out by Feng (2005), Streets et al. (2005, 2011), Pacyna and Pacyna (2006) and Pirrone et al. (2010). Specific emission factors of HMs from nonferrous metal smelting sectors can be seen in the Supplement, Table S16.
Currently, the blast furnace is the primary technique for pig iron production in China. For steel production, there are two main routes: (1) “ore BF-BOF steel-making route” based on blast furnace (BF) and basic oxygen furnace (BOF) techniques, and (2) “scrap-EAF steel-making route” based on an electric arc furnace (EAF) technique, using steel scrap or sponge iron as basic raw materials (Zhang and Wang, 2008). In spite of the environmental friendly nature of EAF and its flexibility to produce a variety of value-added grades of steel, the share of electric furnace steel in Chinese output of crude steel only accounts for about 8.9 % in 2012, mainly due to the shortness of steel scrap resources in China (CISA, 2013).
Comparing the national emission standard of air pollutants for the iron
smelting industry in China with those in certain European Union countries
(e.g., United Kingdom, Germany, Netherlands, Austria), we choose to use
the emission factors of HMs for the iron smelting industry obtained from the
emission factors database of NAEI in 2000 as the national average emission
factors for iron smelting of China in 2015. This is mainly because the PM
emission limit of existing facilities for iron smelting of China in 2015
(20 mg m
Cement, glass, and brick manufacturing are the major mineral commodity industries. During the manufacturing process, various HMs emitted from raw materials as vapor and feed fuels associated with fine particulate matter are emitted from the kiln system at high temperatures.
Currently, the new dry rotary kiln process is the dominant technology in cement manufacturing factories of China, representing over 92 % of the total national cement output. The emission ceilings of air pollutants for cement, glass or brick manufacturing specified in the present standards of China (e.g., GB 4915–2013, GB 26453–2011, GB 29620–2013) are less stringent compared with those of developed countries (see Supplement Table S18). By contrasting the emission limits of air pollutants from non-metallic minerals (cement, glass, and brick) manufacturing between China and developed countries, we presume the best emission factors of air pollutants achieved in China today are approximately identical to the average emission factors of developed countries at the end of the 1990s. With respect to cement production, the unabated emission factors of HMs can be obtained from the Web Factor Information Retrieval System (WebFIRE) (US EPA, 2012). Moreover, the average emission factors for glass and brick manufacturing are assumed to remain unchanged since pre-1996. Subsequently, atmospheric emission factors of HMs from these two sub-source categories decrease gradually with the implementation of gradually tightened emission limits from the above-mentioned Emission Standards of Pollutants from the non-metallic mineral manufacturing industry. Specific emission factors of various HMs from non-metallic mineral manufacturing can be seen in the Supplement, Table S16.
China is the biggest developing country in the world. The rural population still accounted for nearly 47.4 % of the total population in 2012 (NBS, 2013a), and it has had a long history of using agricultural residues and firewood to satisfy household energy demands for cooking and heating. Recently, crop residues have become more commonly burned in open fields during the harvest season. Abundant gaseous and particulate pollutants emitted by open biomass burning have caused severe regional air pollution and contributed to worsening of haze events in central and eastern China (Cheng et al., 2014; Li et al., 2014).
In this paper, a total mass of ten crop straws burned is calculated based on the method discussed in previous studies by Tian et al. (2011b) and Lu et al. (2011), including paddy, wheat, maize, other grains, legumes, tubers, cotton, oil plants, fiber crops, and sugar crops. Because of quite limited field test data about HM emission characteristics from crop straw for household use and firewood for open burning, we presume HM emission factors from biofuel for open burning are equal to those for household use. It is acknowledged that this simple assumption may introduce additional uncertainties, and thus relatively large uncertainty ranges for HM emission factors of biofuel combustion are applied in the analysis, which merits substantial investigation in the future. The average emission factors of HMs from these ten crop straws and firewood are summarized in the Supplement, Table S16.
Besides major conventional pollutants (PM, SO
Historically, leaded gasoline combustion by vehicles has been recognized as
the most significant contributor for the increase in lead level in human blood
(Robbins et al., 2010). Leaded gasoline has been forced out of the market
place in China since 1 July 2000 due to the adverse health effects on the
neurologic and/or hematologic systems (Xu et al., 2012). Compared to the Pb
content limits of 0.64 g L
For municipal solid waste (MSW) incineration, emission characteristics of HMs significantly depend on the concentration of metals in the feed wastes, the performance of installed APCDs, combustion temperatures, as well as composition of the gas stream (Chang et al., 2000).
Presently, stoke grate and fluidized-bed combustion are the major MSW
incineration technologies being used in China. Because of relatively high
costs and the heat content requirement for the feed MSW (
To estimate the hazardous air pollutant emission inventory from MSW incineration in China, Tian et al. (2012d) have compiled and summarized the comprehensive average emission factors of hazardous HMs (Hg, As, Pb, Cd, Cr, Ni, and Sb) for MSW incineration from published literature. Additionally, the emission ceiling of HMs for the existing incinerators in the newly issued standard (GB 18485–2014) which will be conducted in 2016 is approximately comparable to that in Directive 2000/76/EC (see Supplement, Table S18). Here, we presume the best emission factors of HMs in China for MSW in 2016 are almost equivalent to those in developed EU countries in 2000. Based on specific emission factors of HMs for MSW incineration from published literature (Nriagu, 1979; Pacyna, 1984; Nriagu and Pacyna, 1988) and certain emission factors of HMs with uncontrolled technology from AP42, Fifth Edition, Volume I, Chapter 2: Solid Waste Disposal (US EPA, 1996), the unabated emission factors of HMs from this source category are determined. Specific emission factors of various HMs from MSW incineration can be seen in the Supplement, Table S16.
Brake linings as well as tyre wear of vehicles are known as one of the
important emission sources of particulate matter to the surrounding
environment, particularly in urban areas (Hjortenkrans et al., 2007).
Notably, not all of the worn materials of brake lining and tyre will be
emitted into the atmosphere as airborne particulate matter (Hulskotte et al.,
2006). Here, we adopt the average emission factors of TSP from brake wear
and tyre wear for passenger cars (0.0075 and 0.0107 g km
In addition to steel as brake pad support material, the agents present in
brake linings usually consist of Sb, Cu, Zn, Ba, Sn, and Mo (Bukowiecki et
al., 2009). Further, antimony is present in brake linings as
Sb
Because of limited information and lack of field experimental tests on HM contents in brake linings and tyres in Chinese vehicles, and the substantial quantity of vehicles sold in China that are imported from foreign countries or manufactured by foreign-invested transnational vehicle companies, we presume the composition of worn materials from brake and tyre wear in terms of HMs are consistent with foreign countries (see Supplement Table S13).
Coal and liquid fuels consumption data by sectors at the provincial level (e.g., power plant, coal-fired industrial boiler, coal-fired residential sector, other coal-fired sectors) are collected from China Energy Statistical Yearbooks. Industrial production data by provinces (e.g., the output of ferrous/nonferrous metal products, production of cement/glass/brick, amount of municipal waste incineration, number of vehicles) are compiled from relevant statistical yearbooks, such as China Statistical Yearbooks, the Yearbook of Nonferrous Metals Industry of China, China Steel Yearbook etc. The detailed data sources for the main sectors are listed in the Supplement, Table S19. Furthermore, trends of activity levels by different sectors in China between 2000 and 2012 are summarized in the Supplement, Figs. S1–S5.
It is necessary to examine the potential uncertainty in emissions by sources and regions to quantify the reliability, and identify the direction and degree of improvement, of emission inventories in the future. A detailed uncertainty analysis is conducted by combining uncertainties of both activity levels and emission factors, through adopting a Monte Carlo simulation (Zhao et al., 2011; Tian et al., 2014a, b). Streets et al. (2003) indicate that there is no way to judge the accuracy of activity data estimates. Furthermore, uncertainties are still inevitable when representative values are selected for specific emission sources, countries, and regions in spite of emission factors adopted from detailed experiments.
Most of the input parameters of specific activity levels and emission factors, with corresponding statistical distributions, are specified on the basis of the data fitting, or referred to the related published references (Wu et al., 2010; Zhao et al., 2011; Tian et al., 2012a, b). Additionally, for parameters with limited observation data, the probability distributions such as normal distribution and triangular distribution are assumed by the authors for corresponding sources. Further details about the probability distribution for each source discussed in this study are listed in Table S20. Finally, all of the input parameters are placed in a Monte Carlo framework, and simulations are run 10 000 times to estimate the uncertainty ranges of varied HM emissions with a 95 % confidence interval.
HM emissions from primary anthropogenic sources in China,
1949–2012 (t year
Historical trends of atmospheric HMs (Hg, As, Se, Pb, Cd, Cr, Ni, Sb, Mn, Co, Cu, and Zn) emissions from primary anthropogenic sources in China, 1949–2012. CCPP, coal consumption by power plants; CCIB, coal consumption by industrial boilers; CCRS, coal consumption by residential sectors; CCOS, coal consumption by other sectors; LFC, liquid fuels combustion; NFMS, nonferrous metal smelting; FMS, ferrous metal smelting; NMMM, non-metallic minerals manufacturing; B&TW, brake and tyre wear; ONCS, other non-coal sources (including BB, biomass burning; and MSWI, municipal solid waste incineration).
The historical trend of atmospheric emissions of Hg, As, Se, Pb, Cd, Cr, Ni, Sb, Mn, Co, Cu, and Zn by different source categories from 1949 to 2012 are illustrated in Fig. 2. The total emissions of HMs from primary anthropogenic sources since 1949 have shown substantial shifts among varied source categories that reflect technological and economic trends and transition during this period of over 60 years. Within the year after the establishment of the People's Republic of China in 1949, the total emissions of Hg, As, Se, Pb, Cd, Cr, Ni, Sb, Mn, Co, Cu, and Zn from anthropogenic sources are estimated at about 11.5–312.6 t (see Table 2). The discharges of HMs on a national scale increased by 3–20 times from 1949 to 1960 due to the increasing demands for energy consumption and industrial production (especially for the period of the Great Leap Forward from 1958 to 1960, resulting in a remarkable increasing output of industrial products), then, a substantial decrease in 1961 and 1962 of 27.6–55.7 % compared to 1960 on account of the serious imbalance of economic structure and the Great Leap Forward famine caused by policy mistakes together with natural disasters (Kung and Lin, 2003). In spite of the negative growth of heavy metal emissions in individual years such as 1967, 1974, and 1976, the annually averaged growth rates of national emissions of HMs from primary anthropogenic sources were still as high as 0.2–8.4 % during the period of 1963 to 1977.
Subsequently, the policy of openness and reformation was issued by the Chinese central government. With the implementation of this policy from 1978 to 2012, China's GDP has been growing at an average annual growth rate of about 9.8 %, resulting in tremendous energy consumption and enormous output of industrial products. As can be seen from Fig. 2, historically there have been two periods during which the total emissions of HMs (except Pb) increased rapidly after 1978. The first one is the period of 1978 to 2000, except for one remarkable fluctuation from 1998 to 1999, which reflects a decrease in input of raw materials and output of industrial products mainly owing to the influence of the Asian financial crisis (Hao et al., 2002). The second one is the period of the 10th Five-Year Plan (10th FYP, from 2001 to 2005) in which a sharp increase in emissions of Hg, As, Se, Cd, Cr, Ni, Sb, Mn, Co, Cu, and Zn occurred, with the emissions increasing from about 268.0–11 308.6 t in 2001 to about 378.9–15 987.9 t in 2005, at an annual average growth rate of 4.8–12.0 % (see Table 2).
In terms of the lead content requirement in gasoline, the past 64 years
since the foundation of the PR China (1949 to 2012) can be divided into two
phases: the leaded gasoline period (1949 to 1990: gasoline with high lead
content (0.64 g L
Due to technological progress resulting in relatively low emission factors of HMs, and economic development bringing about high coal consumption and industrial products output, the trends of total atmospheric emissions for different HMs in China are diverse during the period of 2006 to 2012 (Cheng et al., 2015). Generally speaking, the national atmospheric emissions of Hg, Pb, Cd, Cr, Sb, Cu, and Zn increased at an annual average growth rate of 1.5–7.2 % from 2006 to 2012. In spite of the remarkable growth in coal consumption and gross industrial production, the national As, Se, Ni, Mn, and Co emissions are well restrained in this period. These are mainly due to the different volatility of these 12 elements during the high temperature process resulting in diverse release rates of furnaces and synergistic removal efficiencies of control measures (Xu et al., 2004).
Due to limited information about historical ground-level concentrations of 12 HMs in different cities in China, the temporal characteristics of atmospheric concentrations of four HMs (As, Pb, Cr, and Cu) in Beijing during 2000 to 2012 are used as valid index to verify whether or not the trend of historical HM emissions is reasonable (see Supplement Fig. S6). The data sources and specific values about atmospheric concentrations of As, Pb, Cr, and Cu in Beijing during 2000 to 2012 are listed in the Supplement, Table S21. It should be acknowledged that the verification method applied in this study has certain limitations on account of sampling discrepancies, including sampling time, sampling site, and detection method etc. Therefore, the historical variation trends of HM emissions may be inconsistent with those of ambient concentrations of HMs in some years.
As can be seen from Fig. S6, minimum values of the atmospheric concentrations of As, Pb, Cr, and Cu occur in 2008. This is mainly because most of the aerosol samples compiled from published papers are collected during August in that year, the time during which China hosted the Beijing Olympics under which a series of strict measures about energy-saving and pollution reduction were implemented, such as suspending production of high polluting industries in the Beijing and neighboring municipalities, restricting the driving of vehicles on alternate days under an even–odd license plate system, limiting pollutant emissions from coal combustion facilities in Beijing and the surrounding provinces etc. Consequently, the variation trends of atmospheric concentrations of As, Pb, Cr, and Cu have some discrepancies with those of historical emissions of the above four HMs in Beijing in 2008. However, the historical emission trends of As, Pb, Cr, and Cu are consistent with those of atmospheric concentrations of the above four HMs during 2000 to 2012 in general (see Fig. S6), which indicates that the historical trend of HM emissions estimated by this study is reasonable.
Estimate of annual Hg emissions from primary anthropogenic
sources among various studies (t yr
Until now, the comprehensive and special studies on various HM (except Hg) emissions in China are quite limited. Therefore, only a detailed comparison with Hg emission estimates from other studies is discussed in this study (see Fig. 3). Specifically, limited data on China's Hg emissions can be cited directly from the global Hg inventories estimated by Pacyna and Pacyna (2001), Pacyna et al. (2006, 2010) and Streets et al. (2011). Consequently, here, we mainly focus on comparing our estimations with the results from China's specialized Hg emission inventories estimated by Streets et al. (2005) and Wu et al. (2006).
Overall, the estimated Hg emissions from fuel combustion (except for the
subcategory of coal consumption by residential sectors) in this work are
substantially consistent with those reported by Streets et al. (2005) and Wu
et al. (2006), although the values calculated for the same year are somewhat
different. This may be mainly attributed to the difference in the averaged
provincial content of Hg in raw coal. In our study, according to a
comprehensive investigation of published literature, we determine the
national averaged Hg content in China to be 0.18 mg kg
The power plant sector represents the largest consumer of coal in China. The
thermal power generation increased from 3.6 TWh in 1949 to 3925.5 TWh in
2012 (NBS, 2013a). Meanwhile, coal burned by power plants has increased from
5.2 to 1785.3 Mt (NBS, 2013b), with an annual growth rate of 9.9 % and a
percentage share of the total coal consumption increasing from 22.7 to
50.6 %. For the period of 1949 to 2005, the emissions of HMs from coal
combustion by power plants increased in rough proportion to coal
consumption. However, this trend began to change after 2006 due to the
implementation of policies of energy-saving and pollution reduction,
especially the strengthening of SO
Presently, the combination of pulverized-coal boilers plus ESPs plus WFGD is the most common APCDs configuration in coal-fired power plants of China. By the end of 2012, the installed capacities of FGD in power plants increased by nearly 14 times compared with those in 2005, reaching about 706.4 GWe, accounting for approximately 86.2 % of the installed capacity of total thermal power plants (MEP, 2014a). Of all of the units with FGD installation, approximately 89.7 % adopt the limestone-gypsum WFGD process. The discharges of Hg, As, Se, Pb, Cd, Cr, Ni, Sb, Mn, Co, Cu, and Zn from coal combustion by power plants in 2012 are estimated at about 15.2–3038.9 t (see Fig. 2), which have decreased by 1.7–11.8 % annually since 2006. Moreover, the distinction of integrated co-benefit removal efficiencies of these elements for the typical APCD configurations is the primary reason for the obvious variations of the declining rates among varied HMs, as illustrated in Table 1 and Fig. 2.
In general, coal combusted by industrial boilers is used to provide hot water and heating for industrial production processes. With the development of China's economy (GDP increased from CNY (Chinese Yuan) 46.6 billion in 1949 to CNY 51 894.2 billion in 2012), coal consumption by industrial boilers increased at a relatively lower growth rate than the power sector, from 11.5 Mt in 1949 to 1205.6 Mt in 2012 (NBS, 2013b). According to the statistical data from China Machinery Industry Yearbook, the combination of stoker-fired boiler plus wet scrubber and cyclone is the most common configuration in coal-fired industrial sectors of China, especially for the small- and medium-scale boilers (CMIF, 2013).
As can be seen from Fig. 2, the emission trends of HMs from coal consumption
by industrial boilers are consistent with the total national emissions
trends between 1949 and 1997, and negative growth appears in 1998 and 1999
due to the decreased coal consumption resultant of the Asian financial
crisis (Hao et al., 2002; Tian et al., 2007, 2012b). Subsequently, the
emissions of different toxic HMs from coal consumption by industrial boilers
have distinct variation tendencies, mainly due to the different removal
efficiencies of HMs through typical APCDs. Generally, Hg and Pb emissions
from coal consumption by industrial boilers increased almost
monotonically from 85.1 and 3717.8 t in 2000 to 179.0 and 5770.0 t in 2012,
with an annual growth rate of about 6.4 and 3.7 %, respectively. However,
the discharges of Mn from coal consumption by industrial boilers
decreased about 1.2 times from 5866.0 to 4951.8 t during this period
(2000–2012). Moreover, the discharges of the other nine HMs (As, Se, Cd, Cr,
Ni, Sb, Co, Cu and Zn) from coal consumption by industrial boilers present a
trend of first an increase, and then a decrease overall, with the
implementation of policies of saving energy and pollution reduction in the
coal-fired industrial boilers sector, especially the growing application of
high-efficiency dust collectors and various types of combined dust and
SO
Historically, a sharp fluctuation of Hg discharges from the nonferrous metal smelting sector occurred in the period of the Great Leap Forward to the Great Leap Forward famine (an increase from 92.6 t in 1957 to 221.7 t in 1959, then a decrease rapidly to 104.0 t in 1963); this is mainly due to the rapid increase or decline of mercury mining outputs in this period (an increase from 1060 t in 1957 to 2684 t in 1959, then a decrease rapidly to 1345 t in 1963). Subsequently, a sharp increase of emissions of Hg occurred, with emissions from about 60.6 t in 1998 increasing to about 218.6 t in 2012, at an annually averaged growth rate of 9.6 %. Simultaneously, the primary contributor of Hg emissions from the nonferrous metal smelting sector has changed to the subsector of primary Zn smelting, which made up about 36.9–52.7 % of the sector during 1998 to 2012. Unlike Hg emissions, the emissions of As, Se, Pb, Cd, Ni, Sb, Cu, and Zn from the nonferrous metal smelting sector increased by approximately 7–15 times to 442.3, 1856.4, 251.8, 412.7, 140.6, 1240.9, and 4025.6 t in 2012, respectively. This is mainly because the reduced shares of HM emissions from the nonferrous metal smelting sector, caused by increasing advanced pollutants control devices installation, were partly counteracted by the rapid growth of nonferrous metals production.
A steady increase of HM emissions from the pig iron and steel industry accompanied by certain undulations occurred from 1949 to 1999 (see Fig. 2). Specifically, because of the emphasis on the backyard furnaces for steel production in the period of the Great Leap Forward Movement, a sharp fluctuation of emissions occurred during the period of 1958 to 1963, with the emissions of Hg, As, Se, Pb, Cd, Cr, Ni, Sb, Mn, Cu, and Zn almost doubling (NBS, 2013b). Although emission factors gradually reduced between 2000 and 2012, the output of pig iron and steel has rapidly increased from 131.0 and 128.5 Mt in 2000 to 663.5 and 723.9 Mt in 2012 and, as a result, the emissions of Hg, As, Se, Pb, Cd, Cr, Ni, Sb, Mn, Cu, and Zn from this sector have quadrupled or quintupled in the past 12 years. Especially, the share of Zn emissions from the ferrous metal smelting sector to the national emissions has increased from 13.1 to 32.2 %. Therein, the steel production industry represents the dominant contributor to Zn emissions, accounting for about 60.9–62.9 % during this period.
In order to facilitate understanding of historical HM emissions in China, details about temporal variation trends of HM emissions from liquid fuel combustion and brake and tyre wear are discussed in the Supplement Sect. S3.
Provincial HM emissions from anthropogenic sources and national composition by source categories in 2010.
The total emissions of Hg, As, Se, Pb, Cd, Cr, Ni, Sb, Mn, Co, Cu, and Zn from primary anthropogenic sources by provinces in China for the year 2010 are estimated at about 72 955.1 t. As can be seen in Fig. 4, coal combustion sources represent the major contributors of Hg, As, Se, Pb, Cr, Ni, Mn, Co, and Cu emissions and are responsible for about 50.6, 74.2, 64.6, 60.1, 90.4, 56.2, 80.9, 98.6, and 53.4 % of total emissions, respectively, while their contribution to the total Cd, Sb, and Zn emissions are relatively lower, at about 32.7, 39.3, and 39.8 %, respectively.
Among all the coal-consuming sub-sectors, coal consumption by industrial boilers ranks as the primary source of total national emissions of 12 HMs, with the average proportion about 57.7 % of the total emissions from coal combustion. This may be attributed to the significant coal consumption of industrial boilers (about 1117.3 Mt in 2010) and the relatively high share of boilers with inadequate APCDs (Cheng et al., 2015; NBS, 2013b).
As the largest coal consumer in China, coal consumption by power plants is
identified as the second largest contributor and accounts for about 14.0 %
in total national emissions of 12 HMs. In order to achieve the emission
reduction of PM, SO
China has been the world's largest producer of pig iron and steel by a
rapidly growing margin. By the end of 2012, the output of steel has amounted
to 723.9 Mt, accounting for about 46 % of worldwide steel production
(CISA, 2013). Despite the enormous achievement by China's iron and steel
industry, China is still featured as a steel producer with low energy
efficiency and high pollutant emissions level, compared with other major
steel-producing countries (Guo and Fu, 2010). Because of limited application
of FGD and de-NO
Regarding nonferrous metal smelting emissions, the primary smelting processes resulting in HM emissions discharge far more pollution than those emitted from the secondary smelting processes. Nonferrous metal smelting, as the fourth largest contributor, accounting for about 11.0 % of total emissions, represents the primary contributor to the discharges of Hg and Cd. Therein, primary Cu smelting contributes the largest part of most of these elements, including 89.5 % for As, 37.3 % for Pb, 74.8 % for Cd, 38.7 % for Ni and 76.6 % for Cu; primary Pb smelting is the major source of Sb and Pb; primary Zn smelting accounts for the largest proportion of Hg and Zn emissions among the nonferrous metal smelting category. In addition, with respect to Hg emissions from the nonferrous metal smelting sector, the mercury smelting industry is the other dominant subcategory source, with a share of about 33.0 % of nonferrous metal smelting emission in 2010.
It can be concluded that the emissions of HMs from brake wear are associated with the vehicle population, vehicle mileage as well as the content of HMs in brake linings and tyres. Currently, numerous studies have reported that airborne HMs (e.g., Sb, Cu, Zn) in urban areas are associated with road traffic and more definitely with emissions from brake wear (Gómez et al., 2005; Hjortenkrans et al., 2007). As can be seen from Fig. 4h and k, the brake and tyre wear sector is the largest source of national Sb and Cu emissions – 39.9 and 26.3 %, respectively. Brake wear is the dominant sub-contributor, accounting for over 99.9 and 99.6 % for Sb and Cu emissions from this sector in 2010, respectively. This is mainly due to the high content of Sb and Cu in the brake linings (see Table S13, Hjortenkrans et al., 2007) and the explosive expansion of the vehicle population in China (see Fig. S5). Nevertheless, the adverse effects of airborne PM originating from brake wear on human health and the ecosystem have still not received sufficient attention from the policymakers as well as the public.
Although the non-metallic mineral manufacturing sector is not the dominant source of most HMs, the discharge of Se from this sector makes it the largest contributor to the total. Within this category, the glass production sector discharges about 92.9 % of the total Se emissions due to the widespread application of selenium powder as a decolorizing agent in the glass production process and the huge output of glass production (Kavlak and Graedel, 2013).
As can be seen from Fig. 4a–l, the source contributions on the provincial
scale in 2010 vary substantially due to the difference of industrial
conformations and energy structures (Cheng et al., 2015; NBS, 2013a, b).
Among the provinces with high HM emissions, Shandong ranks as the largest
province with As, Se, Cd, Ni, Mn, and Cu emissions; accounting for about
8.1–10.6 % of the national emissions; Hebei contributes the largest part
of about 9.3 and 11.3 % to national Pb and Zn emissions respectively;
Guizhou represents the primary province with Hg and Sb emissions; the key
provinces with Cr and Co emissions are found in Yunnan and Shanxi,
respectively. These can be mainly attributed to the follow reasons (NBS,
2013a, b; Wu et al., 2008): (1) the enormous coal consumption of industrial
boilers, the considerable electric power generation, a substantial increase
of vehicle population, and the huge output of industrial products in
Shandong, (2) the flourishing pig iron and steel production in Hebei, (3) the dominant outputs of mercury and obviously high average concentration of
Sb in feed coals in Guizhou (about 6.0
Gridded HM emissions from anthropogenic sources for the year
2010 (0.5
The spatial distribution patterns of HM emissions from anthropogenic sources are illustrated in Fig. 5. In this study, 1796 power plants with capacity larger than 6000 kW, 566 copper/lead/zinc smelting plants, 33 large iron and steel plants, and 101 MSW incineration plants are identified as large point sources and their emissions are precisely allocated at their latitude/longitude coordinates (the geographical distribution of 2496 point sources in China is shown in the Supplement, Fig. S7). It should be noted here, the emissions from point sources of nonferrous metal smelting industry and ferrous metal smelting industry contain two parts: emissions originating from fuel combustion and emissions emitted from industrial production processes. Except for the emissions from point sources discussed above, the remaining anthropogenic sources in the provincial level are all treated as regional area sources. The specific method of geographical location for area sources has been discussed in our previous studies (Tian et al., 2012b, c).
The spatial variation is closely related with the unbalanced economic
development and population density in the Chinese mainland, so that these
12 typical HM emissions are distributed very unevenly from one area to
another, with the annual As emissions at the province level ranging from
0.009 kg km
Moreover, several provinces in the southwestern and central-southern regions also play a prominent role for these 12 HM emissions, especially for Guizhou, Sichuan, Yunnan, Hubei, and Hunan provinces. In general, Guizhou province starts out with high emissions of HMs from coal consumption by other sectors, mainly owing to both the high HM contents in the feed coals and the large magnitude of coal consumption by this sector. In addition, the nonferrous industries of Hunan and Yunnan provinces are flourishing, especially the copper and zinc smelting industries. Consequently, the nonferrous metal smelting sector is seen as one of the major sources of Cu and Zn emissions in these two provinces.
The situations of atmospheric HM concentrations in the aerosols of 44 major cities in China during the last 10 years have been reviewed comprehensively by Duan and Tan (2013). Their results indicate that the ambient concentrations of HMs (As, Pb, Cd, Cr, Ni, Mn, Cu, and Zn) are high in some cities, including Beijing, Tianjin, Shijiazhuang, Shenyang, Harbin, Jinan, Zhengzhou, Hangzhou, Nanjing, Hefei, Xian, Yinchuan, Urumqi, Wuhan, Changsha, Chongqing, Guangzhou, Shenzhen, Foshan, Shaoguan, and others. For HM emissions on the urban scale in 2010, these 20 cities with high HMs concentrations also represent the cities with the highest HM emissions in China (see Fig. 5). In general, the spatial distribution characteristics of gridded HM emissions from primary anthropogenic sources for the year 2010 in this study are reasonable and representative of the real situation of this HM pollution.
Uncertainties in the total emissions of HMs in China in 2010 (uncertainties in the emissions of HMs by source categories in China in 2010 can be seen in the Supplement, Table S22).
Emissions of varied HMs from primary anthropogenic sources with
uncertainties in 2010 are summarized in Fig. 6 and Supplement Table S22. As
can be seen, the overall uncertainties of the total emissions in our
inventories quantified by the Monte Carlo simulation are
The earlier statistical data for activity levels are considered to have high
uncertainty for developing countries (including China) with less developed
statistical systems. Unfortunately, we have to acknowledge that it is quite
difficult to accurately assess the specific uncertainty of activity data
from China's earlier official statistics. Akimoto et al. (2006) argue that
the data of energy consumption of China during 1996–2003 is not recommended
for use in the study of emission inventories due to the probable
underestimates. However, the discrepancies in coal consumption from the
power sector are considered to be less than
The combined uncertainty bounds for the national emissions of 12 HMs
during the historical period are shown in Fig. S8. In general, the range of
uncertainty has gradually diminished over time. For example, we calculate an
uncertainty level of
Generally speaking, emission inventories are never complete and perfect, and most emissions estimates possess a significant associated uncertainty, mainly owing to the lack of representativeness of specific emission factors and the reliability of the source-specific activity data. In this study, we have made great efforts to evaluate the historical trend of these HM emissions by collecting detailed activity levels for various source categories, adopting the best available dynamic emission factors for various anthropogenic sources in China today, and integrating publication literature and reports from developed countries and districts. Nevertheless, considerable uncertainties are still present, and this may lead to under- or overestimation of HM emissions from some source categories. Consequently, more detailed investigations and long-term field tests for all kinds of coal-fired facilities and industrial production processes are in great demand.
Presently, control of atmospheric HM emissions still has not received sufficient attention by the government and public in spite of the frequent occurrence of HM pollution in China (especially for provinces with high point sources of HM emissions). The implementation of more rigorous emission standards on primary anthropogenic sources (thermal power plant, coal-fired boiler, nonferrous metallurgy, pig iron and steel production etc.) and national ambient air quality standards (NAAQS) are regarded as important triggers to promote enterprises with diminished HM emissions. Therefore, the MEP should speed up the revision of the system of hazardous air pollutant (including HM) emission standards, and strengthen the amendment of NAAQS. Especially brake wear has been confirmed to be the main source of HM emissions from traffic, particularly in urban areas. However, there is no related emission standard of air pollutants for brake wear. In the near future, the promulgation of emission standards of brake wear should be expected, which will further strengthen the control of atmospheric HM emissions in China.
In addition, some specific actions are suggested as follows: (1) lower or stop mining and burning coal with high HM concentrations in certain provinces where the coals are mainly mined from small coal mines such as Zhejiang and Guangxi (or lower or stop using high-sulfur coal in corresponding provinces due to the high affinity between HMs and pyrite in coal) (Yuan et al., 2013; Zhu et al., 2015); (2) promote coal washing before combustion (the removal efficiencies of coal preparation to lower heavy metals can reach as high as approximately 30.0–60.0 %, see Supplement Table S9); (3) increase the application rate of advanced APCDs configurations in newly built or retrofitted coal-fired boilers; (4) initiate pilot tests or demonstration projects for specified mercury control (SMC) technologies in some sectors with high Hg emissions and develop comprehensive HM control technologies capable of simultaneously removing multiple heavy metals; (5) strengthen energy conservation and boost electricity and/or heat generation using cleaned energy and renewable energy, such as nuclear, wind and solar energy; (6) suspend small-scale coal-fired boilers and industrial production plants with backward emission control technologies (e.g., cement plants, ferrous smelting plants, nonferrous smelting plants); (7) eliminate outdated production technology, such as VR pyrometallurgy and ISP pyrometallurgy; and (8) improve cyclic utilization rate of nonferrous metals and ferrous metals during the period of 12th FYP etc.
We have calculated the historical emissions of 12 typical HMs from primary man-made activities during the period of 1949–2012, based on the detailed statistical data at a provincial level from various statistical yearbooks and adopting comprehensive time-varying dynamic emission factors from relevant research and literature. Undoubtedly, taking consideration of the economic transition and emission control technology improvement, the dynamic emission factors used in this study will enhance the accuracy and reliability of the estimation of HM emissions.
The total national atmospheric emissions of Hg, As, Se, Pb, Cd, Cr, Ni, Sb, Mn, Co, Cu, and Zn from anthropogenic sources increased by about 22–128 times during the period of 1949–2012, reaching about 526.9–22 319.6 t in 2012.
In spite of the increasing coal consumption and gross industrial production, the national emissions of certain HMs (e.g., As, Se, Ni, Mn, Co) have been well restrained with the implementation of energy-saving and pollution reduction policies during 2006 to 2012. Especially, the declining share of HM emissions from industrial process sources (e.g., nonferrous metal smelting, ferrous metal smelting, non-metallic mineral manufacturing) caused by increasing installation of advanced pollutants control devices, has been partially counteracted by the added industrial production yields. Additionally, both high contents of antimony and copper in brake lining and the rapid growth of civilian vehicle population are thought to be the primary reasons for continuous significant growth rate of Sb and Cu emissions from brake and tyre wear during 2000 to 2012.
The spatial distribution characteristics of HM emissions are closely related with the unbalanced regional economic development and population density in China. One notable characteristic is that HM emission intensities are much higher in central and eastern China than those in western China, and coastal regions are classed as the most polluted areas of HMs. Notably, because of the flourishing of the nonferrous metal smelting industry, the southwestern and central-southern provinces also play a prominent role in HM emissions.
The overall uncertainties in our bottom-up inventories are thought to be reasonable and acceptable with the adequate data availability. Nevertheless, to achieve more reliable estimations of HM emissions in China, much more detailed investigations and long-term field tests for all kinds of coal-fired facilities and industrial processes are still greatly needed in the future.
This work is funded by the National Natural Science Foundation of China (21177012, 21377012 and 40975061), the Special Program on Public Welfare of the Ministry of Environmental Protection (201409022), Open fund of State Environmental Protection Key Laboratory of Sources and Control of Air Pollution Complex (SCAPC201305), and the special fund of State Key Joint Laboratory of Environmental Simulation and Pollution Control (13L02ESPC). We thank Tami Bond and Freed Chuck for their help in improving the English writing in this paper.Edited by: T. Bond