Integrative Application of Life Cycle Assessment and Risk Assessment to Environmental Impacts of Anthropogenic Pollutants at a Watershed Scale.

Intense human activities have led to increasing deterioration of the watershed environment via pollutant discharge, which threatens human health and ecosystem function. To meet a need of comprehensive environmental impact/risk assessment for sustainable watershed development, a biogeochemical process-based life cycle assessment and risk assessment (RA) integration for pollutants aided by geographic information system is proposed in this study. The integration is to frame a conceptual protocol of "watershed life cycle assessment (WLCA) for pollutants". The proposed WLCA protocol consists of (1) geographic and environmental characterization mapping; (2) life cycle inventory analysis; (3) integration of life-cycle impact assessment (LCIA) with RA via characterization factor of pollutant of interest; and (4) result analysis and interpretation. The WLCA protocol can visualize results of LCIA and RA spatially for the pollutants of interest, which might be useful for decision or policy makers for mitigating impacts of watershed development.

An agent-based model for an air emissions cap and trade program: A case study in Taiwan.

To determine the actual status of individuals in a system and the trading interaction between polluters, this study uses an agent-based model to set up a virtual world that represents the Kaohsiung and Pingtung regions in Taiwan, which are under the country's air emissions cap and trade program. The model can simulate each controlled industry's dynamic behavioral condition with the bottom-up method and can investigate the impact of the program and determine the industry's emissions reduction and trading condition. This model can be used elastically to predict the impact of the trading market through adjusting different settings of the program rules or combining the settings with other measures. The simulation results show that the emissions trading market has an oversupply, but we find that the market trading amounts are low. Additionally, we find that increasing the air pollution fee and offset rate restrains the agents' trading decision, according to the simulation results of each scenario. In particular, NOx and SOx trading amounts are easily impacted by the pollution fee, reduction rate, and offset rate. Also, the more transparent the market, the more it can help polluters trade. Therefore, if authorities want to intervene in the emissions trading market, they must be careful in adjusting the air pollution fee and program rules; otherwise, the trading market system cannot work effectively. We also suggest setting up a trading platform to help the dealers negotiate successfully.

The influence of spatial resolution on human health risk co-benefit estimates for global climate policy assessments.

Assessment of the ability of climate policies to produce desired improvements in public health through co-benefits of air pollution reduction can consume resources in both time and research funds. These resources increase significantly as the spatial resolution of models increases. In addition, the level of spatial detail available in macroeconomic models at the heart of climate policy assessments is much lower than that available in traditional human health risk modeling. It is therefore important to determine whether increasing spatial resolution considerably affects risk-based decisions; which kinds of decisions might be affected; and under what conditions they will be affected. Human health risk co-benefits from carbon emissions reductions that bring about concurrent reductions in Particulate Matter (PM10) emissions is therefore examined here at four levels of spatial resolution (Uniform Nation, Uniform Region, Uniform County/city, Health Risk Assessment) in a case study of Taiwan as one of the geographic regions of a global macroeceonomic model, with results that are representative of small, industrialized nations within that global model. A metric of human health risk mortality (YOLL, years of life lost in life expectancy) is compared under assessments ranging from a "uniform simulation" in which there is no spatial resolution of changes in ambient air concentration under a policy to a "highly spatially resolved simulation" (called here Health Risk Assessment). PM10 is chosen in this study as the indicator of air pollution for which risks are assessed due to its significance as a co-benefit of carbon emissions reductions within climate mitigation policy. For the policy examined, the four estimates of mortality in the entirety of Taiwan are 747 YOLL, 834 YOLL, 984 YOLL and 916 YOLL, under Uniform Taiwan, Uniform Region, Uniform County and Health Risk Assessment respectively; or differences of 18%, 9%, 7% if the HRA methodology is taken as the baseline. While these differences are small compared to uncertainties in health risk assessment more generally, the ranks of different regions and of emissions categories as the focus of regulatory efforts estimated at these four levels of spatial resolution are quite different. The results suggest that issues of risk equity within a nation might be missed by the lower levels of spatial resolution, suggesting that low resolution models are suited to calculating national cost-benefit ratios but not as suited to assessing co-benefits of climate policies reflecting intersubject variability in risk, or in identifying sub-national regions and emissions sectors on which to focus attention (although even here, the errors introduced by low spatial resolution are generally less than 40%).

Extraction of heavy metals from contaminated soil by Cinnamomum camphora.

83 acres of rice paddy fields in Taoyuan county, Taiwan, were polluted by cadmium (Cd), chromium (Cr) and copper (Cu) through a nearby irrigation channel, and rice plantation was ceased in 1987. Camphor trees (Cinnamomum camphora) have been planted in 2 acre of the above fields since 1991. Heavy metal accumulation of roots, leaves, branches and heartwood of camphor trees were analyzed during 20-year afforestation. Averaged Cd contents of the roots were found larger than the ones of the branches, leaves, sapwood and heartwood of camphor trees growing in three polluted plots. Averaged diameters at breast height (DBH) of the planted camphor trees were 13-15 cm. Cd pollution did not significantly impact the growth of camphor trees, as similar DBH's were found from both polluted and control sites. Annual growths of DBH were from 0.63 to 0.77 cm year(-1). Planting camphor trees sequestered 68.8 ton biomass per acre. During 20-year period, 0.69-1.98 ton C year(-1) ha(-1) were sequestered on three polluted plots. The above numbers exceeded IPCC LULUCF reference values 0.31-0.53 ton C year(-1) ha(-1) for activities at forest lands.

Development and application of dynamic hybrid multi-region inventory analysis for macro-level environmental policy analysis: a case study on climate policy in Taiwan.

We develop a novel inventory method called Dynamic Hybrid Multi-Region Inventory analysis (DHMRI), which integrates the EEMRIOA and Integrated Hybrid LCA and applies time-dependent environmental intervention information for inventory analysis. Consequently, DHMRI is able to quantify the change in the environmental footprint caused by a specific policy while taking structural changes and technological dynamics into consideration. DHMRI is applied to assess the change in the total CO2 emissions associated with the total final demand caused by the climate policy in Taiwan to demonstrate the practicality of this novel method. The evaluation reveals that the implementation of mitigation measures included in the existing climate policy, such as an enhancement in energy efficiency, promotion of renewable energy, and limitation of the growth of energy-intensive industries, will lead to a 28% increase in the total CO2 emissions and that the main driver is the export-oriented electronics industry. Moreover, a major increase in the total emissions is predicted to occur in Southeast Asia and China. The observations from the case study reveal that DHMRI is capable of overcoming the limitations of existing assessment tools at macro-level evaluation of environmental policies.

Integrating input output analysis with risk assessment to evaluate the population risk of arsenic.

Multimedia and site-specific risk assessments (RA) of major sources releasing arsenic (As) were converted into sector-based risk coefficients, which were integrated with the Input Output Table (IO) to analyze the association between sector activities and health risks. The developed IO-RA framework is a valuable tool for unfolding the risk chain linking the receptors, exposure pathways, emission sources, and production and consumption activities associated with various industrial sectors. The enlarged decision space along the chain can then be considered in planning risk management strategies. This case study estimates that air emissions of As result in 1.54 carcinogenic cases. Export is the primary driving force and accounts for approximately 48% of the final demand that leads to population risks of As. The ranking of the contribution of the five sectors in terms of total population risks is as follows: electricity supply (1.06E+00), steelmaking (2.2 × 10(-1)), cement kilns (1.50 × 10(-1)), semiconductor manufacturing (6.34 × 10(-2)) and incinerators (4.31 × 10(-2)). The electricity supply, steelmaking industry, and cement kilns are the major sectors, not only because their emissions directly cause risk but also because they have a stronger influence on the risk generated by other sectors.

Applying probabilistic material flow analysis for quality control and management of waste recycling in steelmaking.

In modern steelmaking, multiple processes comprise a continuous manufacturing system, but not all phosphorus content data are connected or integrated into a holistic and systematic database. Disconnected data hinder the improvement of material management and resource efficiency in the industry. The objective of this study was to establish a method to evaluate material flows, reduce uncertainty, and perform quality control for waste recycling in the steelmaking industry. The results indicate that 10% of the phosphorus input is present in the final products, 30% accumulates in the slags, and more than 60% of the total mass remains in the processes. Comparing the material flow analysis results obtained using static and probabilistic approaches, the partition ratio of the phosphorus content in slags changes from 24.07% to 40.78%, but that in processes changes from 49.10% to 68.05%. This indicates that the variations in phosphorus content in slags and processes might affect the effectiveness of slag recycling and might increase the resource consumption required to maintain the quality of final products. The probability of forming substandard products in the baseline scenario is 0.43. Adopting a 50% removal rate, the probabilities of forming substandard products are reduced to 0.36 (waste removal scenario), 0.38 (slag reduction scenario), and 0.31 (raw material treatment scenario). The performance of raw material treatment and waste removal is more efficient for quality control. The method used in this study can be applied to evaluate the possible outcomes of waste recycling and reduce the probability of forming substandard products.

Health risk assessment of trace elements of ambient PM2.5 under monsoon patterns.

In order to identify the contribution to health risk derived from various emission sources, this study investigated monsoon variations in PM2.5 mass and concentrations of the associated trace elements in a region with complex pollution sources in central Taiwan. This study applied the Chemical Mass Balance model to analyze the source contribution of PM2.5. The source apportionment to obtain the risk contribution of different sources were conducted for different monsoon periods according to the monsoon patterns. In this way, the contributions of individual sources and chemicals to health risk under different monsoon types can be understood to support development of effective control strategies. Among the top contributors of PM2.5 during the north-east monsoon were Secondary Aerosol 28.93% >Coal Boiler 19.82% >Crustal Dust 15.99%; in south-west monsoon were Coal Boiler 37.29% >Traffic Emission 21.19% >Secondary Aerosol 17.84%. The total risk of cancer was above the acceptable risk (3.07 × 10-6), while the non-carcinogenic risk was within the acceptable range (0.262). The variation in the concentration and composition of PM2.5 was related to the change of monsoon type. During the north-east monsoon, the air mass had a long transmission distance and the PM2.5 concentration was relatively high. During the south-west monsoon, the air mass had a short transmission distance and the composition was mainly influenced by nearby emission sources, which resulted in higher risk due to chemical characteristics. To provide sound air quality management, attention should be paid to the composition of PM2.5 in addition to its concentration.

Applying multi-criteria decision-making to improve the waste reduction policy in Taiwan.

Over the past two decades, the waste reduction problem has been a major issue in environmental protection. Both recycling and waste reduction policies have become increasingly important. As the complexity of decision-making has increased, it has become evident that more factors must be considered in the development and implementation of policies aimed at resource recycling and waste reduction. There are many studies focused on waste management excluding waste reduction. This study paid more attention to waste reduction. Social, economic, and management aspects of waste treatment policies were considered in this study. Further, a life-cycle assessment model was applied as an evaluation system for the environmental aspect. Results of both quantitative and qualitative analyses on the social, economic, and management aspects were integrated via the Technique for Order Preference by Similarity to Ideal Solution (TOPSIS) method into the comprehensive decision-making support system of multi-criteria decision-making (MCDM). A case study evaluating the waste reduction policy in Taoyuan County is presented to demonstrate the feasibility of this model. In the case study, reinforcement of MSW sorting was shown to be the best practice. The model in this study can be applied to other cities faced with the waste reduction problems.

Uncertainty analysis of remediation cost and damaged land value for brownfield investment.

Flexible land management that allows risk-based remediation approaches can increase local financial assistance in brownfield development investments. The evaluation of brownfields investment has always been mired in uncertainty because of the complexity of contamination and land policy. To clarify the uncertainty in investment costs resulting from a flexible risk-based land management, this research analyzes the uncertainty of remediation costs (Ir) and damaged land value (Ip). Sizable statistical data for 30,000 abandoned factories in Taiwan are used to derive the uncertainty of investment cost. The results show that the range and value of the uncertainty of Ip is more than Ir because the present land value and the change in flexible land management are the main sensitivity factors. The saved investment costs at the 90% confidence interval ranges five orders of magnitude. Moreover, investors could have saved approximately 5.74E+13 USD resulting from Ip and Ir. The study's two major findings with regard to investment risk in brownfield redevelopment are: (1) an explanation of the extent to which investment costs are to be reduced for investors by adjusting remediation goals depending on different land types; (2) characterization of the uncertainty of brownfields in the real investment market for brownfield redevelopment.

Twice the effort: Ineffectiveness of selecting air pollution control targets with emission quantity for risk reduction.

Protection of human health from air pollution has been typically pursued primarily via regulations of air quality standards and emission standards. Although reducing air pollution from the largest sources and placing more stringent emission limits on the industries of focus is a criterion used by decision makers to control air pollution, it is not clear whether this criterion is the most effective and efficient in improving health protection. Pollutants released from sources into the environment are spatially fluctuating rather than uniformly distributed, and hence, health risk is an issue of geographic variability. To address this issue, this study used a representative example of lead (Pb) in Taiwan. This study implemented an IO-RA methodology to redefine the effectiveness of air pollution management and rank the control priorities of target industries using different perspectives, i.e., environmental responsibility, economic benefit and repercussion potential. This study also considered the potential differences in policy effectiveness based on the air pollution control targets and ranked the industries according to their effectiveness in health risk improvement across the three perspectives and pure emission quantities. After determining the cause-effect chain of health risk through IO-RA, authorities can partner with specific industries according to the chosen effectiveness criteria and thus facilitate better policy performance.

Assessment of the value of reducing uncertainty by sampling in a groundwater remediation system.

The success of cleaning up a contaminated site depends on the degree of knowledge of the site's characteristics. When there is much uncertainty associated with the knowledge, the uncertainty regarding whether the remediation will work increases consequently. It is therefore essential to know how much reduction of uncertainty is needed for the purpose of designing a successful and reliable remediation system. The understanding of the site characteristics is basically increased by site investigation, and thus the uncertainty is decreased by sampling information. This study develops a method to evaluate the value of reducing uncertainty by sampling the hydrogeological parameters in a groundwater remediation system. Hydraulic conductivity being taken as an example of the site characteristics, random field generation and conditional simulation are coupled to obtain a range of hydraulic conductivity fields based on the sampling outputs. A multiple-realization management model that incorporates a chance constraint of health risk is used to find the lowest remediation cost under specific remediation criteria of risk through genetic algorithm. The remediation cost, which is expected to decrease with collection of more samples, serves as the measure of the value of uncertainty reduction by sampling. A case study shows that the variation of the hydraulic conductivity fields among the potential sites as well as the remediation cost is reduced as a result of increase of samples. It also shows that the risk after remediation decreases with the collection of more samples, which implies that the reduction of risk can also be used to assess the value of sampling.

A systemic health risk assessment for the chromium cycle in Taiwan.

Health risk assessment (HRA) has been recognized as a useful tool for identifying health risks of human activities. In particular, this method has been well applied to spatially defined units, such as a production plant, a treatment facility, and a contaminated site. However, the management strategies based on the risk information will be more efficient if the comprehensive picture of total risks from all kinds of sources is depicted. In principle, the total risks can be obtained when all risk sources are assessed individually. Apparently, this approach demands huge amount of efforts. This study develops a methodology that combines substance flow and risk estimation to facilitate examination of risk in a systemic way and provide comprehensive understanding of risk generation and distribution corresponding to flows of substances in the anthroposphere and the environment. Substance flow analysis (SFA) and HRA method is integrated to produce a systemic risk assessment method, from which substance management schemes can be derived. In this study, the chromium cycle in Taiwan is used as an example to demonstrate the method, by which the associated substance flow in the economy and the risk caused by the substance in the environmental system is determined. The concentrations of pollutants in the environmental media, the resultant risks and hazard quotients are calculated with the widely-used CalTOX multimedia model.

Assessment of water quality management with a systematic qualitative uncertainty analysis.

Uncertainty is an inevitable source of noise in water quality management and will weaken the adequacy of decisions. Uncertainty is derived from imperfect information, natural variability, and knowledge-based inconsistency. To make better decisions, it is necessary to reduce uncertainty. Conventional uncertainty analyses have focused on quantifying the uncertainty of parameters and variables in a probabilistic framework. However, the foundational properties and basic constraints might influence the entire system more than the quantifiable elements and have to be considered in initial analysis steps. According to binary classification, uncertainty includes quantitative uncertainty and non-quantitative uncertainty, which is also called qualitative uncertainty. Qualitative uncertainty originates from human subjective and biased beliefs. This study provides an understanding of qualitative uncertainty in terms of its conceptual definitions and practical applications. A systematic process of qualitative uncertainty analysis is developed for assisting complete uncertainty analysis, in which a qualitative network could then be built with qualitative relationship and quantifiable functions. In the proposed framework, a knowledge elicitation procedure is required to identify influential factors and their interrelationship. To limit biased information, a checklist is helpful to construct the qualitative network. The checklist helps one to ponder arbitrary assumptions that have often been taken for granted and may yield an incomplete or inappropriate decision analysis. The total maximum daily loads (TMDL) program is used as a surrogate for water quality management in this study. 15 uncertainty causes of TMDL programs are elicited by reviewing an influence diagram, and a checklist is formed with tabular interrogations corresponding to each uncertainty cause. The checklist enables decision makers to gain insight on the uncertainty level of the system at early steps as a convenient tool to review the adequacy of a TMDL program. Following the instruction of the checklist, an appropriate algorithm in a form of probability, possibility, or belief may then be assigned for the network. Consequently, the risk or evidence of the success of outcomes will be obtained. The incorporation of the systematic consideration of qualitative uncertainty into water quality management is expected to refine the decision-making process.

Combining the cost of reducing uncertainty with the selection of risk assessment models for remediation decision of site contamination.

The decision as to whether a contaminated site poses a threat to human health and should be cleaned up relies increasingly upon the use of risk assessment models. However, the more sophisticated risk assessment models become, the greater the concern with the uncertainty in, and thus the credibility of, risk assessment. In particular, when there are several equally plausible models, decision makers are confused by model uncertainty and perplexed as to which model should be chosen for making decisions objectively. When the correctness of different models is not easily judged after objective analysis has been conducted, the cost incurred during the processes of risk assessment has to be considered in order to make an efficient decision. In order to support an efficient and objective remediation decision, this study develops a methodology to cost the least required reduction of uncertainty and to use the cost measure in the selection of candidate models. The focus is on identifying the efforts involved in reducing the input uncertainty to the point at which the uncertainty would not hinder the decision in each equally plausible model. First, this methodology combines a nested Monte Carlo simulation, rank correlation coefficients, and explicit decision criteria to identify key uncertain inputs that would influence the decision in order to reduce input uncertainty. This methodology then calculates the cost of required reduction of input uncertainty in each model by convergence ratio, which measures the needed convergence level of each key input's spread. Finally, the most appropriate model can be selected based on the convergence ratio and cost. A case of a contaminated site is used to demonstrate the methodology.

Site-specific health risk assessment of dioxins and furans in an industrial region with numerous emission sources.

While the sources of releasing dioxins have received much attention, the health risks of incinerators are studied considerably more than those of other sources, such as metallurgical industry. Furthermore, risk assessments focus primarily on the effects of single sources; few have addressed the aggregate risks from multiple sources in a region. When many sources of dioxins exist in an area, such as the Industrial Park located in the densely populated district-Siaogang District in southern Taiwan where 17 primary known PCDD/F stationary emission sources are clustered, is the aggregate exposure to these sources imposing high risks even though individual sources comply with emission standards? The study combined the multimedia and multiple pathway exposure modeling and site-specific exposure scenario to assess dioxin risks contributed by the 17 emission sources in the District, including municipal waste incinerators, medical waste incinerators, sinter plants, electric arc furnace, secondary aluminum smelters, cement kilns, etc. The average cancer risk of a resident living in the District was found to be 3.43E-04 under the site-specific exposure scenario. The top emission source is the sinter plant, followed by the electric arc furnace. The information has driven the local government to conduct more complete assessment and at the same time to consider enforcing a stricter local standard of dioxin emissions in the Siaogang District.

Oxidation of methyl methacrylate from semiconductor wastewater by O3 and O3/UV processes.

This study investigated the oxidation of methyl methacrylate (MMA) by sole ozonation and ozone/UV treatments. The semi-batch ozonation experiments were proceeded under different reaction conditions to study the effects of ozone dosage and UV radiation on the oxidation of MMA. The experimental results indicated that both the oxidation of MMA by the sole ozonation and O3/UV processes can completely decompose MMA to form the following intermediates within 30 min reaction time. To increase the applied ozone dosage can significantly raise the removal efficiency of MMA. However, the mineralization of MMA via the direct oxidation reaction of molecular ozone was slow, while introducing the UV radiation can promote the mineralization rate of MMA. In addition, the pH value of the oxidized solution in the O3/UV treatment decreased lower than that in the sole ozonation treatment of about 1 unit. The possible scheme of the decomposition pathway of MMA under the ozonation process is proposed in this study. Formic acid and acetic acid were found to be the main ozonated intermediates.

A novel sustainable decision making model for municipal solid waste management.

This paper reviews several models developed to support decision making in municipal solid waste management (MSWM). The concepts underlying sustainable MSWM models can be divided into two categories: one incorporates social factors into decision making methods, and the other includes public participation in the decision-making process. The public is only apprised or takes part in discussion, and has little effect on decision making in most research efforts. Few studies have considered public participation in the decision-making process, and the methods have sought to strike a compromise between concerned criteria, not between stakeholders. However, the source of the conflict arises from the stakeholders' complex web of value. Such conflict affects the feasibility of implementing any decision. The purpose of this study is to develop a sustainable decision making model for MSWM to overcome these shortcomings. The proposed model combines multicriteria decision making (MCDM) and a consensus analysis model (CAM). The CAM is built up to aid in decision-making when MCDM methods are utilized and, subsequently, a novel sustainable decision making model for MSWM is developed. The main feature of CAM is the assessment of the degree of consensus between stakeholders for particular alternatives. A case study for food waste management in Taiwan is presented to demonstrate the practicality of this model.

An investigation on the potential of metal recovery from the municipal waste incinerator in Taiwan.

This study aimed to identify distribution of metals and to estimate the amount of these metals that can be potentially recovered from incineration residues. First, the partitioning behavior of Cr, Cu, Fe, Cd, Al, Zn, and Pb in bottom ash and fly ash was investigated in one large municipal waste incinerator in Taiwan. In addition, the material flow analysis (MFA) method was used to estimate the material flux of metals within incinerator plant, and to calculate the amount of metal recovery. According to the findings of this study, six metals (Fe, Al, Cu, Zn, Cr, and Pb) concentrated in bottom ash mostly, while Cd existed primarily in fly ash. The weight percentages of Fe (4.49%), Al (5.24%), Cu (1.29%), Zn (2.21%), and Pb (0.58%) in incinerator ash are high, and even higher than the compositions of natural minerals. Finally, the amount of Cr, Cu, Fe, Cd, Al, Zn and Pb that can be potentially recovered from incineration residues will reach 2.69 x 10(2), 1.46 x 10(4), 4.91 x 10(4), 6.92 x 10(1), 5.10 x 10(4), 1.85 x 10(4) and 4.66 x 10(3) ton/yr, respectively.

Identification of pollution source of cadmium in soil: application of material flow analysis and a case study in Taiwan.

BACKGROUND: Since the 1970s, at least 200 hectares (ha) of farm-land has been polluted by the heavy metal cadmium (Cd). Consequently, the Cd pollution has led to contaminate the rice production and caused acute social panic. According to the recent investigation results performed by the Taiwan Environmental Protection Administration (TEPA), it is indicated that most of the Cd pollution incidents in Taiwan resulted from the waste-water discharge of stearate Cd factories. To prevent the Cd pollution incidents from spreading, the TEPA has either forced these factories to close down or assisted them in improving their production processes since the 1980s. Unfortunately, accidental incidents of Cd pollution still emerge in an endless stream, despite the strict governmental controls placed on these questionable factories. Whether this pollution has resulted from undetected or hidden pollution sources stemming from two decades ago or comes from some new source, will be an outstanding issue. Therefore, this study attempts to identify the pollution sources of Cd in soil in Taiwan as well as to find the solution to the above-mentioned, outstanding issue by way of a methodology termed Material Flow Analysis (MFA). METHODOLOGY: The MFA has proved to be a useful tool on providing quantitative information of the flow of substances through an economic to an environmental system. Based upon the supply-and-demand theory of MFA, researchers have successfully conducted an overview of the use of materials in many industries, the construction industry being one of these. Therefore, this study tries to establish a set of analytical processes by way of MFA for identifying the pollution source of Cd in soil in Taiwan. In addition, the spirit of Life Cycle Assessment (LCA) technique was also employed to identify the materials, and products should be ignored as a crucial pollution source in this study. RESULTS AND DISCUSSION: According to the MFA methodology applied in this study and on the basis of related studies performed by Taiwanese governmental authorities, we arrive at the following analysis results: (1) the total amount of Cd from the economic perspective of material and product flow was approximately 441.2 tons; (2) the wastewater directly discharged into irrigation water can be concluded to be the major pollution route of Cd in farmland soil in Taiwan; (3) material plastic stabilizer (cadmium oxide, CdO), Zn-Pd compounds and Cu compounds should be the crucial pollution sources to contaminate environment through the route of wastewater in Cd flow analysis; (4) the crucial pollution sources to contaminate environment through the route of wastewater in Cd flow analysis were five factories, Coin, Jili, Taiwan Dye, Guangzheng and Mingguan, and they were all categorized as stearate Cd industries; (5) the typical source of the Cd pollution in soil in Changhua County through the pollution route of wastewater should be the metal surfacing process industries. CONCLUSIONS: This study proved that MFA can be a good tool for identifying Cd flow as well as for recognizing the crux of the problem concerning incidents of Cd pollution. This study led to the conclusion that the causal relationship between farmland pollution caused by Cd and stearate Cd factories in Taiwan seemed quite close by way of MFA methodology. In addition, this study also found that the wastewater discharged from a single metal surfacing process factory will not cause remarkable farmland pollution. However, the wastewater simultaneously discharged from a group of pollution factories can result in a significant pollution incident. RECOMMENDATIONS AND OUTLOOK: This case study is only a small contribution to the understanding of the toxic material flow related to Cd in the environment. This study recommends that Taiwanese governmental authorities should not deal with problems on an ad hoc basis, but should instead deal with Cd pollution problems overall employing control measures. Finally, the more accurate information or data we can collect, the more reliable results we can identify. Therefore, the quality and quantity of related data used in this MFA model should be closely scrutinized in order to ensure the most correct and comprehensive investigation on the toxic material flow.