hundreds of millions of people throughout the world are working under circumstances that foster ill health or unsafe . it is estimated that yearly over two million people worldwide die of occupational injuries and work - related diseases . in fact , more people die from diseases caused by work than are killed in industrial accidents ( 1 ) . according to american welding society ( aws ) and edison welding institution ( ewi ) ( 2 ) , welding will continue to be the preferred method of joining for world class product until 2020 . although there is a wide breadth of hazards that exist in welding operations , only 2% of occupational safety and health association ( osha ) general industry citations addressing on this matter ( 3 ) . previous researches had highlighted the challenges for developing countries in strategies for risk assessment and control in welding industries . transfer of technologies of welding from developed countries to developing countries which do not have similar infrastructures in terms of health and safety may be disastrous . uncritical adoption of new welding technologies by developing countries potentiates future health problems ( 4 , 5 ) . to this date , there is still very limited study that discussed the relationship between welding emissions with health risk of welder in automotive industry especially in malaysia ( 6 ) . currently the welding fumes exposure risk assessments were largely focused on single welding fumes constituents approach because the regulatory standard for compliance only caters for a single constituent . however , in reality , welders are simultaneously exposed to multiple welding fumes constituent at once . assessment of the hazards of multiple simultaneous exposures only had been done in limited study ( 7 ) . according to dominici et al . ( 8) , the shift from a single pollutant to multiple pollutant assessment was desirable by the scientific community and policy makers . karkoszka and sokovic ( 9 ) , developed the integrated risk estimation in welding process using qualitative method of assigning probability of occurrence , significance and risk involve in aspect of occupational and safety . yeo and neo ( 10 ) , on the other hand introduce the health hazard scoring system to quantify the environmental impact of the different welding process before choosing the most environmentally friendly welding processes . however , these models did not consider the quantitative data on welding fumes exposure and the developed tools had not been verified with actual data . on the other hand , leman et al . ( 11 ) had developed an environmental quality index ( eqi ) for industrial ventilation and occupational safety and health evaluation of welding processes in manufacturing plants . although the index has been developed based on actual data on welding exposure , there were no analysis had been done on the selection of the aggregation model used in this study . thus , there was still gap in developing a suitable risk assessment method relating to welding fume exposure to possible health risk of welder in quantitative manners . research needs had been highlighted to pursue a means of indexing exposure by job type or process by taking into account the intensity of the welding job and work practiced ( 12 ) . however , welders are not a homogeneous group , the potential adverse effect of welding fume exposures are oftentimes difficult to evaluate . differences exist in welder populations , such as industrial setting , types of ventilation , type of welding processes and materials used ( 13 ) . however , indexing exposure according to the location would be benefited as ranking tools between different locations on the same scale . welding risk assessment would be simpler if a single metric could embody all of the information in the measurement ( 14 ) . hence , this study aims to develop an index that can rank welding workplace that associate well with possible health risk of welders . the investigation was conducted in two automotive related industrial plants working on spot gun , spot weld and robotic metal inert gas ( mig ) weld . plant 1 consists of 53 male welders while plant 2 consists of 44 male welders . plant 1 had the average 12 hour working shifts while plant 2 has 14 hours average working shifts . fifty three non - welder male workers that did not have continuous exposure to welding fumes were selected from similar workplaces as control . another 30 male welders from plant 3 that work for average 8 hours at automotive assembly industries were investigated for index verification purposed . lung function test ( lft ) were performed on handheld spirometer ( micro medical dl , uk ) connected to spirometer software ( care fusion , san diego ) on a notebook computer . interviews were conducted before conducting ma neuvers to record demographic data , smoking habit and working experience . tests were conducted according to forced vital capacity procedure of the american thoracic society recommends ( 15 ) . measured parameters were forced vital capacity ( fvc ) , forced expiratory volume in 1 second ( fev1 ) and peak expiratory flow ( pef ) were all expressed as a percentage of the predicted value and fev1/fvc ratio . the predicted set used in this study interpretation and derivation of the value of normal , obstruction and restriction lung function result were done according to the american thoracic society ( ats ) ( 15 ) . personal samplings of welding fumes were conducted in plant 1 , plant 2 and plant 3 during march to june 2013 . personal sampling method was based on british standard guidelines bs en 689:1996 which stated at least one employee in ten of properly selected homogeneous group performing similar tasks must be sampled ( 17 ) . the filters media ( mixed cellulose ester 0.8 m pore sizes ) was used with sampling pump set to 2 l / min flow rate . personal sampling of welding fumes was done with the objective to get exposure on maximum risk workers . thus , in situation where more than one samples were obtained , the results with the highest concentration in most of the constituents were selected . in malaysia , under the occupational safety and health act 1994 , use and standards of exposure of chemical hazardous to health regulation ( usechh ) ( 18 ) , chemical classified hazardous to health with its permissible emission limits ( pel ) were listed and need to be comply by the employer . the analysis in certified laboratory was done based on american society for testing and materials ( astm ) d7439 - 08 method by using inductively coupled plasma mass spectrometry ( icp - ms ) ( agilent 7700 ) with microwave digestion ( nitric acid and hydrochloric acid ) . at least one field blanks were submitted together with batch of samples for each investigated plants to the accredited laboratory for analysis . duration of sampling was calculated and the concentrations of exposure were calculated in time weighted average 8 hours ( twa 8) . from a regulatory compliance perspective , threshold levels of controlled parameters are established in the context of possible adverse impacts to human health . it will be useful to relate the index to some acceptance parameters that are measurable . development of environmental index involves following four basic steps as shown in fig . 1 ( 1921 ) . basic steps in development of environmental index ( 19 ) each of these steps was explained in the next subsection . selection of relevant factors and parametersin this study , the analyses of welding fumes were conducted by icp - ms . currently there are only two standard method for determination of constituents in airborne particulate matters by using icp - ms , which is astm d7439 - 08 ( 22 ) and british standard ( bs ) international standardization organization ( iso ) bs iso 30011:2010 ( 23 ) . icp - ms has the advantage to analyze up to 25 multi constituents in a single sample . from these 25 constituents , only 15 constituents ( aluminum , an timony , arsenic , beryllium , cadmium , chromium , cobalt , copper , ferum , lead , manganese , molybdenum , nickel , silver and tin ) were shortlisted according to the constituents commonly associated with welding , cutting and brazing ( 24 , 25).transformation of selected parameters into sub - indexthe development of dose - effect information was often regarded as highly simplistic and not readily accepted by researches in epidemiology field . it is often not possible to identify the dose - effect information that applies to individual pollutant and properly covers all segments of the population . the dose - effect function must contend with the complexity of controlling extraneous factor that gave impact on the observed effect ( 26 ) . this scenario resulted in limited dose - effect information for pollutants available in the literatures . however , in risk assessment , the ideas on combining dose and health risk ( dose - risk ) were widely implemented ( 10 , 27 , 28 ) . the dose and risk effect in this study follow the dose - risk model for inhalable toxicity by ( 10 , 27 , 28 ) as shown in eq . 1 . sub  index  value  =  doses  rating      health  risk  rating  ( eq .  1)doses ratingthe 15 metal constituents ( aluminum , antimony , arsenic , beryllium , cadmium , chromium , cobalt , copper , ferum , lead , manganese , molybdenum , nickel , silver and tin ) were selected as parameters of sub - indexes . the doses rating values were derived from a segmented linear function or also known as staircase step function as shown in fig.2 ( 26 , 29 , 30 ) . 2 shows the relation between pollution concentration and doses rating values . by taking the pel as a reference line , the concentration of welding fumes below 5% of the reference line is considered as the moderately serious dose with a rating value of 2 . the concentration of welding fumes exceeding the reference line is considered as the lethal dose with a rating value of 3 . mild effect dose with rating value 1 is considered as 5% concentration of welding fumes below the reference line to limit of detection of value . for concentration below the limit of detection , 2relation between doses rating and pollution concentrationhealth risk ratingsfour health risks were considered in this study according to national institute for occupational safety and health ( niosh ) pocket guide to chemical hazard ( 24 ) ; sensitizer , respiratory toxins , target organ toxins and carcinogen . table 1 shows the how the health risk rating were categorized.table 1criteria for health risk ratinghealth risk ratingsensitizerrespiratory toxinstarget organ toxinscarcinogen0 ( no effect)no observed health riskno observed health riskno observed health riskno observed health risk1 ( mild)one sensitizer health risknose , nasal cavitiesone target organ health riskiarc 3,4 tlv a4,a52 ( moderately serious)two sensitizer health riskspharynx , larynx , tracheatwo target organs health risksiarc 2a,2b tlv a2,a33 ( lethal)more than two sensitizer health riskslower respiratory tracts : lung , bronchioles , alveolimore than two target organs health risksiarc 1 tlv a1the health risk ratings according to each investigated constituents were tabulated in table 2 . arsenic and cadmium had the highest rating score of 9 while aluminum , silver and tin had the lowest total rating score of 1 . the health risk information was extracted from niosh pocket guide to chemical hazard ( 24 ) for the investigated 15 welding fumes constituents.table 2health risk rating by constituentsmetal constituentssensitizerrespiratory toxinstarget organ toxinscarcinogentotal health risk ratingsaluminum ( al)01102antimony ( sb)02103arsenic ( as)03339beryllium ( be)03036cadmium ( cd)03339chromium ( cr)02013cobalt ( co)13026copper ( cu)02204iron ( fe)03003lead ( pb)01326manganese ( mn)03306molybdenum ( mo)02226nickel ( ni)23027silver ( ag)01001tin ( sn)01001derivation of weightthe weight of each metal constituent was derived according to their pel . metal constituents with lower pel has higher weight as shown in table 3 . as for manganese metal constituents , the pel is according to the ceiling value . the weight was selected so that their sum is unity.table 3weight according to constituentsno.constituentsusecch pel ( mg / m)weight1aluminum ( al)5.0 ( resp . ) 15.0 ( total)0.0252antimony ( sb)0.50.0503arsenic ( as)0.0100.1004beryllium ( be)0.002 c 0.0050.1905cadmium ( cd)0.0050.1906chromium ( cr)0.50.0507cobalt ( co)0.10.0508copper ( cu)1.00.0259iron ( fe)100.02010lead ( pb)0.050.10011manganese ( mn)c50.02512molybdenum ( mo)5.0 ( soluble ) 15 ( total insoluble)0.02513nickel ( ni)1.00.02514silver ( ag)0.010.10015tin ( sn)2.00.025total1.000aggregation modelthe aggregation process is the crucial part in calculation of environmental index . they affect the quality of results in many ways because aggregation process is where most of the simplifying process ( reduction of information ) takes place ( 26 , 31).aggregation model consist of ; additive form , multiplicative form and maximum or minimum operator form as shown in table 4 . in this study , 15 welding fumes constituents were considered as sub - indices . there was possibility that some of the constituents were below the detection limit leading to  0 value of sub index , thus a multiplicative form of aggregation model was not suitable for analysis in this study . maximum and minimum operators were also excluded from this study because these types of operators are biased towards extreme ( minimum or maximum ) sub index values . most of the air pollution indices reported in literatures use the additive form aggregation model and developed in the increasing scale form ( higher index portray the severe condition ) ( 3234 ) . following this , table 4list of aggregation model 
 penalty function of aggregation modelscombination of sub index by using aggregation models commonly arise issues such as ambiguity , eclipsing , compensation , and rigidity . these issues were explained in table 5 . in order to compare the aggregation models quantitatively , sadiq et . al ( 19 ) had proposed the usage of penalty functions in order to select the most appropriate aggregation models in a specific condition as shown in table 5 . it was highlighted by ( 19 ) that in the search of better aggregation model , a trade - off exists between properties ; a model may perform very well against one property , but perform poorly against another . therefore , in a selection of a model producing either ambiguous or eclipsed results , the index developer must consider to what degree the ambiguous or eclipsed result is acceptable . thus , the penalty function analysis must be done to compare and select the most suitable aggregation model for the developed index based on data collected in the case studies conducted . penalties are derived from the penalty function such that they are continuous over an interval , where maximum penalty. selection of relevant factors and parameters in this study , the analyses of welding fumes were conducted by icp - ms . currently there are only two standard method for determination of constituents in airborne particulate matters by using icp - ms , which is astm d7439 - 08 ( 22 ) and british standard ( bs ) international standardization organization ( iso ) bs iso 30011:2010 ( 23 ) . icp - ms has the advantage to analyze up to 25 multi constituents in a single sample . from these 25 constituents , only 15 constituents ( aluminum , an timony , arsenic , beryllium , cadmium , chromium , cobalt , copper , ferum , lead , manganese , molybdenum , nickel , silver and tin ) were shortlisted according to the constituents commonly associated with welding , cutting and brazing ( 24 , 25 ) . transformation of selected parameters into sub - index the development of dose - effect information was often regarded as highly simplistic and not readily accepted by researches in epidemiology field . it is often not possible to identify the dose - effect information that applies to individual pollutant and properly covers all segments of the population . the dose - effect function must contend with the complexity of controlling extraneous factor that gave impact on the observed effect ( 26 ) . this scenario resulted in limited dose - effect information for pollutants available in the literatures . however , in risk assessment , the ideas on combining dose and health risk ( dose - risk ) were widely implemented ( 10 , 27 , 28 ) . the dose and risk effect in this study follow the dose - risk model for inhalable toxicity by ( 10 , 27 , 28 ) as shown in eq . 1 . sub  index  value  =  doses  rating      health  risk  rating  ( eq .  1 ) the 15 metal constituents ( aluminum , antimony , arsenic , beryllium , cadmium , chromium , cobalt , copper , ferum , lead , manganese , molybdenum , nickel , silver and tin ) were selected as parameters of sub - indexes . the doses rating values were derived from a segmented linear function or also known as staircase step function as shown in fig.2 ( 26 , 29 , 30 ) . 2 shows the relation between pollution concentration and doses rating values . by taking the pel as a reference line , the concentration of welding fumes below 5% of the reference line is considered as the moderately serious dose with a rating value of 2 . the concentration of welding fumes exceeding the reference line is considered as the lethal dose with a rating value of 3 . mild effect dose with rating value 1 is considered as 5% concentration of welding fumes below the reference line to limit of detection of value . for concentration below the limit of detection , 2relation between doses rating and pollution concentration relation between doses rating and pollution concentration four health risks were considered in this study according to national institute for occupational safety and health ( niosh ) pocket guide to chemical hazard ( 24 ) ; sensitizer , respiratory toxins , target organ toxins and carcinogen . table 1criteria for health risk ratinghealth risk ratingsensitizerrespiratory toxinstarget organ toxinscarcinogen0 ( no effect)no observed health riskno observed health riskno observed health riskno observed health risk1 ( mild)one sensitizer health risknose , nasal cavitiesone target organ health riskiarc 3,4 tlv a4,a52 ( moderately serious)two sensitizer health riskspharynx , larynx , tracheatwo target organs health risksiarc 2a,2b tlv a2,a33 ( lethal)more than two sensitizer health riskslower respiratory tracts : lung , bronchioles , alveolimore than two target organs health risksiarc 1 tlv a1 criteria for health risk rating the health risk ratings according to each investigated constituents were tabulated in table 2 . arsenic and cadmium had the highest rating score of 9 while aluminum , silver and tin had the lowest total rating score of 1 . the health risk information was extracted from niosh pocket guide to chemical hazard ( 24 ) for the investigated 15 welding fumes constituents . table 2health risk rating by constituentsmetal constituentssensitizerrespiratory toxinstarget organ toxinscarcinogentotal health risk ratingsaluminum ( al)01102antimony ( sb)02103arsenic ( as)03339beryllium ( be)03036cadmium ( cd)03339chromium ( cr)02013cobalt ( co)13026copper ( cu)02204iron ( fe)03003lead ( pb)01326manganese ( mn)03306molybdenum ( mo)02226nickel ( ni)23027silver ( ag)01001tin ( sn)01001 health risk rating by constituents the weight of each metal constituent was derived according to their pel . metal constituents with lower pel has higher weight as shown in table 3 . as for manganese metal constituents , the pel is according to the ceiling value . table 3weight according to constituentsno.constituentsusecch pel ( mg / m)weight1aluminum ( al)5.0 ( resp . ) 15.0 ( total)0.0252antimony ( sb)0.50.0503arsenic ( as)0.0100.1004beryllium ( be)0.002 c 0.0050.1905cadmium ( cd)0.0050.1906chromium ( cr)0.50.0507cobalt ( co)0.10.0508copper ( cu)1.00.0259iron ( fe)100.02010lead ( pb)0.050.10011manganese ( mn)c50.02512molybdenum ( mo)5.0 ( soluble ) 15 ( total insoluble)0.02513nickel ( ni)1.00.02514silver ( ag)0.010.10015tin ( sn)2.00.025total1.000 weight according to constituents the aggregation process is the crucial part in calculation of environmental index . they affect the quality of results in many ways because aggregation process is where most of the simplifying process ( reduction of information ) takes place ( 26 , 31 ) . aggregation model consist of ; additive form , multiplicative form and maximum or minimum operator form as shown in table 4 . in this study , 15 welding fumes constituents were considered as sub - indices . there was possibility that some of the constituents were below the detection limit leading to  0 value of sub index , thus a multiplicative form of aggregation model was not suitable for analysis in this study . maximum and minimum operators were also excluded from this study because these types of operators are biased towards extreme ( minimum or maximum ) sub index values . most of the air pollution indices reported in literatures use the additive form aggregation model and developed in the increasing scale form ( higher index portray the severe condition ) ( 3234 ) . following this , list of aggregation model penalty function of aggregation models combination of sub index by using aggregation models commonly arise issues such as ambiguity , eclipsing , compensation , and rigidity . these issues were explained in table 5 . in order to compare the aggregation models quantitatively , sadiq et . al ( 19 ) had proposed the usage of penalty functions in order to select the most appropriate aggregation models in a specific condition as shown in table 5 . it was highlighted by ( 19 ) that in the search of better aggregation model , a trade - off exists between properties ; a model may perform very well against one property , but perform poorly against another . therefore , in a selection of a model producing either ambiguous or eclipsed results , the index developer must consider to what degree the ambiguous or eclipsed result is acceptable . thus , the penalty function analysis must be done to compare and select the most suitable aggregation model for the developed index based on data collected in the case studies conducted . penalties are derived from the penalty function such that they are continuous over an interval , where 1 refers to  maximum penalty. aggregation model s issues and penalty function adapted from sadiq et . al . ( 19 ) once the four penalty functions were calculated , a representative value of a cumulative penalty ( pc ) is derived to compare different model . the value of cumulative penalty that increase by an increase of  is suitable for the development of an index that considered compensation and rigidity as important characteristic . on the other hand , the value of a cumulative penalty that decrease by an increase of  is suitable for the development of an index that considered ambiguous and eclipsing as important characteristic . analysis of variance ( anova ) and multivariate analysis of variance ( manova ) was used to compare mean lung function parameters between welders and control groups . pearson correlation analysis was done to get association between working duration , smoking duration and type of welding with lung function . further analysis using multiple regression analysis was done to confirm the predictors of the lung functions increased / decreased value . the investigation was conducted in two automotive related industrial plants working on spot gun , spot weld and robotic metal inert gas ( mig ) weld . plant 1 consists of 53 male welders while plant 2 consists of 44 male welders . plant 1 had the average 12 hour working shifts while plant 2 has 14 hours average working shifts . fifty three non - welder male workers that did not have continuous exposure to welding fumes were selected from similar workplaces as control . another 30 male welders from plant 3 that work for average 8 hours at automotive assembly industries were investigated for index verification purposed . lung function test ( lft ) were performed on handheld spirometer ( micro medical dl , uk ) connected to spirometer software ( care fusion , san diego ) on a notebook computer . interviews were conducted before conducting ma neuvers to record demographic data , smoking habit and working experience . tests were conducted according to forced vital capacity procedure of the american thoracic society recommends ( 15 ) . measured parameters were forced vital capacity ( fvc ) , forced expiratory volume in 1 second ( fev1 ) and peak expiratory flow ( pef ) were all expressed as a percentage of the predicted value and fev1/fvc ratio . the predicted set used in this study was taken from pneumobile project , indonesia ( 16 ) . interpretation and derivation of the value of normal , obstruction and restriction lung function result were done according to the american thoracic society ( ats ) ( 15 ) . personal samplings of welding fumes were conducted in plant 1 , plant 2 and plant 3 during march to june 2013 . personal sampling method was based on british standard guidelines bs en 689:1996 which stated at least one employee in ten of properly selected homogeneous group performing similar tasks must be sampled ( 17 ) . the filters media ( mixed cellulose ester 0.8 m pore sizes ) was used with sampling pump set to 2 l / min flow rate . personal sampling of welding fumes was done with the objective to get exposure on maximum risk workers . thus , in situation where more than one samples were obtained , the results with the highest concentration in most of the constituents were selected . in malaysia , under the occupational safety and health act 1994 , use and standards of exposure of chemical hazardous to health regulation ( usechh ) ( 18 ) , chemical classified hazardous to health with its permissible emission limits ( pel ) were listed and need to be comply by the employer . the analysis in certified laboratory was done based on american society for testing and materials ( astm ) d7439 - 08 method by using inductively coupled plasma mass spectrometry ( icp - ms ) ( agilent 7700 ) with microwave digestion ( nitric acid and hydrochloric acid ) . at least one field blanks were submitted together with batch of samples for each investigated plants to the accredited laboratory for analysis . duration of sampling was calculated and the concentrations of exposure were calculated in time weighted average 8 hours ( twa 8) . from a regulatory compliance perspective , threshold levels of controlled parameters are established in the context of possible adverse impacts to human health . it will be useful to relate the index to some acceptance parameters that are measurable . development of environmental index involves following four basic steps as shown in fig . 1 ( 1921 ) . basic steps in development of environmental index ( 19 ) each of these steps was explained in the next subsection . selection of relevant factors and parametersin this study , the analyses of welding fumes were conducted by icp - ms . currently there are only two standard method for determination of constituents in airborne particulate matters by using icp - ms , which is astm d7439 - 08 ( 22 ) and british standard ( bs ) international standardization organization ( iso ) bs iso 30011:2010 ( 23 ) . icp - ms has the advantage to analyze up to 25 multi constituents in a single sample . from these 25 constituents , only 15 constituents ( aluminum , an timony , arsenic , beryllium , cadmium , chromium , cobalt , copper , ferum , lead , manganese , molybdenum , nickel , silver and tin ) were shortlisted according to the constituents commonly associated with welding , cutting and brazing ( 24 , 25).transformation of selected parameters into sub - indexthe development of dose - effect information was often regarded as highly simplistic and not readily accepted by researches in epidemiology field . it is often not possible to identify the dose - effect information that applies to individual pollutant and properly covers all segments of the population . the dose - effect function must contend with the complexity of controlling extraneous factor that gave impact on the observed effect ( 26 ) . this scenario resulted in limited dose - effect information for pollutants available in the literatures . however , in risk assessment , the ideas on combining dose and health risk ( dose - risk ) were widely implemented ( 10 , 27 , 28 ) . the dose and risk effect in this study follow the dose - risk model for inhalable toxicity by ( 10 , 27 , 28 ) as shown in eq sub  index  value  =  doses  rating      health  risk  rating  ( eq .  1)doses ratingthe 15 metal constituents ( aluminum , antimony , arsenic , beryllium , cadmium , chromium , cobalt , copper , ferum , lead , manganese , molybdenum , nickel , silver and tin ) were selected as parameters of sub - indexes . the doses rating values were derived from a segmented linear function or also known as staircase step function as shown in fig.2 ( 26 , 29 , 30 ) . 2 shows the relation between pollution concentration and doses rating values . by taking the pel as a reference line , the concentration of welding fumes below 5% of the reference line is considered as the moderately serious dose with a rating value of 2 . the concentration of welding fumes exceeding the reference line is considered as the lethal dose with a rating value of 3 . mild effect dose with rating value 1 is considered as 5% concentration of welding fumes below the reference line to limit of detection of value . for concentration below the limit of detection , 2relation between doses rating and pollution concentrationhealth risk ratingsfour health risks were considered in this study according to national institute for occupational safety and health ( niosh ) pocket guide to chemical hazard ( 24 ) ; sensitizer , respiratory toxins , target organ toxins and carcinogen . table 1 shows the how the health risk rating were categorized.table 1criteria for health risk ratinghealth risk ratingsensitizerrespiratory toxinstarget organ toxinscarcinogen0 ( no effect)no observed health riskno observed health riskno observed health riskno observed health risk1 ( mild)one sensitizer health risknose , nasal cavitiesone target organ health riskiarc 3,4 tlv a4,a52 ( moderately serious)two sensitizer health riskspharynx , larynx , tracheatwo target organs health risksiarc 2a,2b tlv a2,a33 ( lethal)more than two sensitizer health riskslower respiratory tracts : lung , bronchioles , alveolimore than two target organs health risksiarc 1 tlv a1the health risk ratings according to each investigated constituents were tabulated in table 2 . arsenic and cadmium had the highest rating score of 9 while aluminum , silver and tin had the lowest total rating score of 1 . the health risk information was extracted from niosh pocket guide to chemical hazard ( 24 ) for the investigated 15 welding fumes constituents.table 2health risk rating by constituentsmetal constituentssensitizerrespiratory toxinstarget organ toxinscarcinogentotal health risk ratingsaluminum ( al)01102antimony ( sb)02103arsenic ( as)03339beryllium ( be)03036cadmium ( cd)03339chromium ( cr)02013cobalt ( co)13026copper ( cu)02204iron ( fe)03003lead ( pb)01326manganese ( mn)03306molybdenum ( mo)02226nickel ( ni)23027silver ( ag)01001tin ( sn)01001derivation of weightthe weight of each metal constituent was derived according to their pel . metal constituents with lower pel has higher weight as shown in table 3 . as for manganese metal constituents , the weight was selected so that their sum is unity.table 3weight according to constituentsno.constituentsusecch pel ( mg / m)weight1aluminum ( al)5.0 ( resp . ) 15.0 ( total)0.0252antimony ( sb)0.50.0503arsenic ( as)0.0100.1004beryllium ( be)0.002 c 0.0050.1905cadmium ( cd)0.0050.1906chromium ( cr)0.50.0507cobalt ( co)0.10.0508copper ( cu)1.00.0259iron ( fe)100.02010lead ( pb)0.050.10011manganese ( mn)c50.02512molybdenum ( mo)5.0 ( soluble ) 15 ( total insoluble)0.02513nickel ( ni)1.00.02514silver ( ag)0.010.10015tin ( sn)2.00.025total1.000aggregation modelthe aggregation process is the crucial part in calculation of environmental index . they affect the quality of results in many ways because aggregation process is where most of the simplifying process ( reduction of information ) takes place ( 26 , 31).aggregation model consist of ; additive form , multiplicative form and maximum or minimum operator form as shown in table 4 . in this study , 15 welding fumes constituents were considered as sub - indices . there was possibility that some of the constituents were below the detection limit leading to  0 value of sub index , thus a multiplicative form of aggregation model was not suitable for analysis in this study . maximum and minimum operators were also excluded from this study because these types of operators are biased towards extreme ( minimum or maximum ) sub index values . most of the air pollution indices reported in literatures use the additive form aggregation model and developed in the increasing scale form ( higher index portray the severe condition ) ( 3234 ) . following this , only additive forms of aggregation model table 4list of aggregation model 
 penalty function of aggregation modelscombination of sub index by using aggregation models commonly arise issues such as ambiguity , eclipsing , compensation , and rigidity . these issues were explained in table 5 . in order to compare the aggregation models quantitatively , al ( 19 ) had proposed the usage of penalty functions in order to select the most appropriate aggregation models in a specific condition as shown in table 5 . it was highlighted by ( 19 ) that in the search of better aggregation model , a trade - off exists between properties ; a model may perform very well against one property , but perform poorly against another . therefore , in a selection of a model producing either ambiguous or eclipsed results , the index developer must consider to what degree the ambiguous or eclipsed result is acceptable . thus , the penalty function analysis must be done to compare and select the most suitable aggregation model for the developed index based on data collected in the case studies conducted . penalties are derived from the penalty function such that they are continuous over an interval , where  0 maximum penalty. selection of relevant factors and parameters in this study , the analyses of welding fumes were conducted by icp - ms . currently there are only two standard method for determination of constituents in airborne particulate matters by using icp - ms , which is astm d7439 - 08 ( 22 ) and british standard ( bs ) international standardization organization ( iso ) bs iso 30011:2010 ( 23 ) . icp - ms has the advantage to analyze up to 25 multi constituents in a single sample . from these 25 constituents , only 15 constituents ( aluminum , an timony , arsenic , beryllium , cadmium , chromium , cobalt , copper , ferum , lead , manganese , molybdenum , nickel , silver and tin ) were shortlisted according to the constituents commonly associated with welding , cutting and brazing ( 24 , 25 ) . transformation of selected parameters into sub - index the development of dose - effect information was often regarded as highly simplistic and not readily accepted by researches in epidemiology field . it is often not possible to identify the dose - effect information that applies to individual pollutant and properly covers all segments of the population . the dose - effect function must contend with the complexity of controlling extraneous factor that gave impact on the observed effect ( 26 ) . this scenario resulted in limited dose - effect information for pollutants available in the literatures . however , in risk assessment , the ideas on combining dose and health risk ( dose - risk ) were widely implemented ( 10 , 27 , 28 ) . the dose and risk effect in this study follow the dose - risk model for inhalable toxicity by ( 10 , 27 , 28 ) as shown in eq . 1 . sub  index  value  =  doses  rating      health  risk  rating  ( eq .  1 ) the 15 metal constituents ( aluminum , antimony , arsenic , beryllium , cadmium , chromium , cobalt , copper , ferum , lead , manganese , molybdenum , nickel , silver and tin ) were selected as parameters of sub - indexes . the doses rating values were derived from a segmented linear function or also known as staircase step function as shown in fig.2 ( 26 , 29 , 30 ) . 2 shows the relation between pollution concentration and doses rating values . by taking the pel as a reference line , the concentration of welding fumes below 5% of the reference line is considered as the moderately serious dose with a rating value of 2 . the concentration of welding fumes exceeding the reference line is considered as the lethal dose with a rating value of 3 . mild effect dose with rating value 1 is considered as 5% concentration of welding fumes below the reference line to limit of detection of value . for concentration below the limit of detection , 2relation between doses rating and pollution concentration relation between doses rating and pollution concentration four health risks were considered in this study according to national institute for occupational safety and health ( niosh ) pocket guide to chemical hazard ( 24 ) ; sensitizer , respiratory toxins , target organ toxins and carcinogen . table 1criteria for health risk ratinghealth risk ratingsensitizerrespiratory toxinstarget organ toxinscarcinogen0 ( no effect)no observed health riskno observed health riskno observed health riskno observed health risk1 ( mild)one sensitizer health risknose , nasal cavitiesone target organ health riskiarc 3,4 tlv a4,a52 ( moderately serious)two sensitizer health riskspharynx , larynx , tracheatwo target organs health risksiarc 2a,2b tlv a2,a33 ( lethal)more than two sensitizer health riskslower respiratory tracts : lung , bronchioles , alveolimore than two target organs health risksiarc 1 tlv a1 criteria for health risk rating the health risk ratings according to each investigated constituents were tabulated in table 2 . arsenic and cadmium had the highest rating score of 9 while aluminum , silver and tin had the lowest total rating score of 1 . the health risk information was extracted from niosh pocket guide to chemical hazard ( 24 ) for the investigated 15 welding fumes constituents . table 2health risk rating by constituentsmetal constituentssensitizerrespiratory toxinstarget organ toxinscarcinogentotal health risk ratingsaluminum ( al)01102antimony ( sb)02103arsenic ( as)03339beryllium ( be)03036cadmium ( cd)03339chromium ( cr)02013cobalt ( co)13026copper ( cu)02204iron ( fe)03003lead ( pb)01326manganese ( mn)03306molybdenum ( mo)02226nickel ( ni)23027silver ( ag)01001tin ( sn)01001 health risk rating by constituents the weight of each metal constituent was derived according to their pel . metal constituents with lower pel has higher weight as shown in table 3 . as for manganese metal constituents , table 3weight according to constituentsno.constituentsusecch pel ( mg / m)weight1aluminum ( al)5.0 ( resp . ) 15.0 ( total)0.0252antimony ( sb)0.50.0503arsenic ( as)0.0100.1004beryllium ( be)0.002 c 0.0050.1905cadmium ( cd)0.0050.1906chromium ( cr)0.50.0507cobalt ( co)0.10.0508copper ( cu)1.00.0259iron ( fe)100.02010lead ( pb)0.050.10011manganese ( mn)c50.02512molybdenum ( mo)5.0 ( soluble ) 15 ( total insoluble)0.02513nickel ( ni)1.00.02514silver ( ag)0.010.10015tin ( sn)2.00.025total1.000 weight according to constituents the aggregation process is the crucial part in calculation of environmental index . they affect the quality of results in many ways because aggregation process is where most of the simplifying process ( reduction of information ) takes place ( 26 , 31 ) . aggregation model consist of ; additive form , multiplicative form and maximum or minimum operator form as shown in table 4 . in this study , 15 welding fumes constituents were considered as sub - indices . there was possibility that some of the constituents were below the detection limit leading to  0 value of sub index , thus a multiplicative form of aggregation model was not suitable for analysis in this study . maximum and minimum operators were also excluded from this study because these types of operators are biased towards extreme ( minimum or maximum ) sub index values . most of the air pollution indices reported in literatures use the additive form aggregation model and developed in the increasing scale form ( higher index portray the severe condition ) ( 3234 ) . following this , list of aggregation model penalty function of aggregation models combination of sub index by using aggregation models commonly arise issues such as ambiguity , eclipsing , compensation , and rigidity . these issues were explained in table 5 . in order to compare the aggregation models quantitatively , sadiq et . al ( 19 ) had proposed the usage of penalty functions in order to select the most appropriate aggregation models in a specific condition as shown in table 5 . it was highlighted by ( 19 ) that in the search of better aggregation model , a trade - off exists between properties ; a model may perform very well against one property , but perform poorly against another . therefore , in a selection of a model producing either ambiguous or eclipsed results , the index developer must consider to what degree the ambiguous or eclipsed result is acceptable . thus , the penalty function analysis must be done to compare and select the most suitable aggregation model for the developed index based on data collected in the case studies conducted . penalties are derived from the penalty function such that they are continuous over an interval , where al . ( 19 ) once the four penalty functions were calculated , a representative value of a cumulative penalty ( pc ) is derived to compare different model . the value of cumulative penalty that increase by an increase of  is suitable for the development of an index that considered compensation and rigidity as important characteristic . on the other hand , the value of a cumulative penalty that decrease by an increase of  is suitable for the development of an index that considered ambiguous and eclipsing as important characteristic . analysis of variance ( anova ) and multivariate analysis of variance ( manova ) was used to compare mean lung function parameters between welders and control groups . pearson correlation analysis was done to get association between working duration , smoking duration and type of welding with lung function . further analysis using multiple regression analysis was done to confirm the predictors of the lung functions increased / decreased value . lung function data ( mean  standard deviation ( sd ) ) of welders on each plant and control groups was shown in table 6 . apparently lung function results showed that the mean of all lung function parameters of plant 1 and plant 2 were lower than the control group . the multivariate manova analysis reveals there was a significant mean difference of lung function values between groups f ( 10 , 282 ) = 6 . further separate univariate anovas on the outcome variable reveals welders on each plant had significant difference of the mean value of fev1/fvc , f ( 2 , 146 ) =3 . table 6mean values for control , plant 1 , plant 2 and plant 3criteriacontrol n=52 ( meansd)plant 1 n=53 ( meansd)plant 2 n=44 ( meansd)plant 3 n=30 ( meansd)age34.567.6530.625.9628.845.5529.739.04fvc ( % pred)88.33 12.1984.0915.7987.8613.2087.20 12.90fev1 ( % pred)94.5812.4088.5115.3091.1412.9990.8311.48fev1/fvc107.946.38105.919.89103.615.55104.877.82pef ( % pred)84.6711.9368.5816.0771.6814.5079.5315.64 mean values for control , plant 1 , plant 2 and plant 3 index verification was done by conducting a lung function and welding fumes investigation in an automotive assembly plant ( refer as plant 3 afterwards ) . lung function data ( mean  sd ) of welders on each plant and control groups was also shown in table 6 . apparently lung function results showed that the mean of all lung function parameters of plant 1 , plant 2 and plant 3 were lower than the control group . the multivariate manova analysis reveals there was a significant mean difference of lung function values between groups f ( 12 , 512 ) =3 . 84 , further separate univariate anovas on the outcome variable reveals welders on each plant had significant difference of the mean value of fev1/fvc , f ( 3 , 175 ) = 2 . table 7 shows the results of welding fumes personal sampling collected in plant 1 and plant 2 according to welding job type . iron was the highest constituent concentration in all plants with 0.602 ( mg / m ) detected in spot gun welding job process in plant 2 . welding fumes personal sampling several constituents such as antimony , beryllium , cadmium , cobalt , molybdenum , nickel and tin were detected below the limit of detection in all plants . the results of the field blanks did not show any contaminants or unusual concentration of metal constituents detected during sampling or handling of samples . iron constituents had the highest concentration with 0.633 mg / m3 . however , none of the constituents were exceeding the pel . welding fume concentration for plant 3 table 9 shows the correlation between lung functions of welders in each plant with working years , smoking years and type of welding . pearson correlation reveals there was a significant relationship in each plant as follows ; pearson correlation between lung functions of welders in each plant with working years , smoking years and type of welding plant 1 : significant relationship between fvc and the number of working years ( r=- . 40 , p ( two tailed ) < . 001 ) and significant relationship between fev1 and the number of working years ( r=- . 40 , p ( two tailed ) < . further analysis by multiple regressions ( backward stepwise method ) confirmed number of working years was the significant predictor to the decreased value of fvc . however , the number of smoking years and type of welding was not the significant predictors of decreased values of fvc . number of working years was also the significant predictor to the decreased value of fev1 . however , the number of smoking years and type of welding was not the significant predictors of decreased values fev1 . plant 2 : significant relationship between fev1/fvc and number of working years ( r=. 31 , p ( two tailed ) < . . further analysis by multiple regressions ( backward stepwise method ) confirmed number of working years was also the significant predictor to the increase values of fev1/fvc . however , the number of smoking years and type of welding was not the significant predictors of the increase values of fev1/fvc . plant 3 : significant relationship between fev1/fvc and number of smoking years ( r=. 58 , p ( two tailed ) < . . further analysis by multiple regressions ( backward stepwise method ) confirmed number of smoking years was the significant predictor to the increase values of fev1/fvc . however , working duration and type of welding were not the significant predictors of fev1/fvc . apparently high values fev1/fvc relates to a restrictive disorder which contradicts with smoking effects which was the low values of fev1/fvc ( obstruction disorder ) . the welding fume exposures often associate with restrictive disorders ( 4547 ) while smoking is associated with obstruction disorder ( 48 , 49 ) . there were also synergistic relations between the effects of smoking and welding exposure causing lung impairment reported by ( 5052 ) . to clarify these issues , analysis between smoker and nonsmoker welder thus , multiple regression analysis ( backward stepwise method ) was conducted again for fev1/fvc value as dependent variables , smoking status ( smoker and nonsmoker ) and working group ( less and more than 5 years working experience ) and as predictors . these results showed a synergistic relationship between the effects of number of smoking years and welding exposure for more and less than 5 years working experience causing restrictive disorder in plant 3 . table 10 shows the mean values of fvc , fev1 , fev1/fvc and pef in plant 1 , 2 and 3 . mean values for the control , plant 3 ( adjusting for smoker * ) , plant 1 and plant 2 the mean number of cigarette smoke by welder was 3 , 4 and 7 for plant 1 , 2 and 3 respectively . it is clear that welders in plant 3 smoke in average 2 times higher ( number of cigarette ) than plant 1 and 2 . this exclusion decreased the average number of cigarette smoke for plant 3 from 7 to 5 . for smoking years were not the significant predictors for plant 1 and 2 , no adjustment towards smoking years were made from these plants . data tabulated in table 8 were used for comparison between index value and percentage of predicted pulmonary function s value at the end of this section . table 11 shows the mean penalty functions for nine aggregation model ( explained in table 4 ) based on welding fume concentration in plant 1 and 2 . mean penalty functions values for aggregation models cumulative penalties pc for nine aggregation models were calculated and tabulated in table 12 . in the development of wfhi , ambiguity and eclipsing were considered more important because the index should not suffer from overestimation or underestimation of health risk related to welding fumes exposure . the values of the aggregation model of iam , iwam , iwrss and irmsa decreased by an increase of  thus suitable for index that considered ambiguity and eclipsing were more important ( 19 ) . from this four shortlisted aggregation models , iwam and iwrss were selected because the smallest penalty function values . cumulative penalty values for aggregation models table 13 shows the mean value of p1 , p2 , p3 and p4 for iwam and iwrss aggregation models . 3 , where the larger area would present poor performance , and vice versa ( 19 ) . these figures are plotted on a four quadrant point which represents the penalties of four characteristic properties ( ambiguity , eclipsing , compensation , rigidity ) as shown in fig . penalty function s value for iwam and iwrss aggregation model comparison of iwam and iwrss mean penalty values from the four quadrant analysis , iwrss had the smallest triangle area . table 14 shows the index value for plant 1 , 2 and 3 by implementing iwrss on each welding job type . index value for plant 1 and plant 2 with iwrss aggregation model figure 4 shows the relation between pulmonary functions with mean index value of each plant . the results of the study showed mean index value was directly proportional to percentage predicted of welder s lung functions in all the investigated plants . plant 1 has the highest index value which also has the lowest values of fvc , fev1 and pef . plant 3 has the lowest index value which also has the highest values of fvc , fev1 and pef . however , the values of , fev1/fvc were the lowest in plant 3 which in the agreement with restrictive disorder commonly exist in welders . the pulmonary functions of the investigated welders can not be compared directly with previous study from other researcher mainly due to the predicted set used were differs from one study to another . there were also differences existed in welding fumes analysis method , digestion and constituents being investigated which had limited the comparison to be made between this study and previous study by other researcher . relation between percentages predicted of lung function and index value for plant 1 , 2 and 3 table 11 shows the mean penalty functions for nine aggregation model ( explained in table 4 ) based on welding fume concentration in plant 1 and 2 . mean penalty functions values for aggregation models cumulative penalties pc for nine aggregation models were calculated and tabulated in table 12 . in the development of wfhi , ambiguity and eclipsing were considered more important because the index should not suffer from overestimation or underestimation of health risk related to welding fumes exposure . the values of the aggregation model of iam , iwam , iwrss and irmsa decreased by an increase of  thus suitable for index that considered ambiguity and eclipsing were more important ( 19 ) . from this four shortlisted aggregation models , iwam and iwrss were selected because the smallest penalty function values . cumulative penalty values for aggregation models table 13 shows the mean value of p1 , p2 , p3 and p4 for iwam and iwrss aggregation models . 3 , where the larger area would present poor performance , and vice versa ( 19 ) . these figures are plotted on a four quadrant point which represents the penalties of four characteristic properties ( ambiguity , eclipsing , compensation , rigidity ) as shown in fig . 3 . penalty function s value for iwam and iwrss aggregation model comparison of iwam and iwrss mean penalty values from the four quadrant analysis , iwrss had the smallest triangle area . table 14 shows the index value for plant 1 , 2 and 3 by implementing iwrss on each welding job type . figure 4 shows the relation between pulmonary functions with mean index value of each plant . the results of the study showed mean index value was directly proportional to percentage predicted of welder s lung functions in all the investigated plants . plant 1 has the highest index value which also has the lowest values of fvc , fev1 and pef . plant 3 has the lowest index value which also has the highest values of fvc , fev1 and pef . however , the values of , fev1/fvc were the lowest in plant 3 which in the agreement with restrictive disorder commonly exist in welders . the pulmonary functions of the investigated welders can not be compared directly with previous study from other researcher mainly due to the predicted set used were differs from one study to another . there were also differences existed in welding fumes analysis method , digestion and constituents being investigated which had limited the comparison to be made between this study and previous study by other researcher . relation between percentages predicted of lung function and index value for plant 1 , 2 and 3 the conclusions of this work are as follows 
 there was possibility that some of the metal constituents were below the detection limit leading to  0 value of sub index , thus the multiplicative form of aggregation model was not suitable for analysis in this study . on the other hand , maximum or minimum operator forms suffer from compensation issues and were not considered in this study.it is important that the developed index should not suffer from overestimation ( ambiguity ) or underestimation ( eclipsing ) of health risk related to welding fumes exposure . thus , four aggregation models ( iam , iwam , iwrss and irmsa ) were shortlisted according to the cumulative penalty value calculated.the penalty function analysis ( p1 , p2 , p3 and p4 ) calculated from metal fume concentration in plant 1 and 2 , suggest that iwrss had the minimal penalty function values and the best performance in describing wfhi.results from 3 case studies successfully relate mean index value with welder s lung functions in each investigated plant . there was possibility that some of the metal constituents were below the detection limit leading to  0 value of sub index , thus the multiplicative form of aggregation model was not suitable for analysis in this study . on the other hand , maximum or minimum operator forms suffer from compensation issues and were not considered in this study . it is important that the developed index should not suffer from overestimation ( ambiguity ) or underestimation ( eclipsing ) of health risk related to welding fumes exposure . thus , four aggregation models ( iam , iwam , iwrss and irmsa ) were shortlisted according to the cumulative penalty value calculated . the penalty function analysis ( p1 , p2 , p3 and p4 ) calculated from metal fume concentration in plant 1 and 2 , suggest that iwrss had the minimal penalty function values and the best performance in describing wfhi . results from 3 case studies successfully relate mean index value with welder s lung functions in each investigated plant . ethical issues ( including plagiarism , informed consent , misconduct , data fabrication and/or falsification , double publication and/or submission , redundancy , etc . ) have been completely observed by the authors .
abstractbackgroundthis study aimed to develop an index that can rank welding workplace that associate well with possible health risk of welders.methodswelding fumes health index ( wfhi ) were developed based on data from case studies conducted in plant 1 and plant 2 . 
 personal sampling of welding fumes to assess the concentration of metal constituents along with series of lung function tests was conducted . 
 fifteen metal constituents were investigated in each case study . 
 index values were derived from aggregation analysis of metal constituent concentration while significant lung functions were recognized through statistical analysis in each plant.resultsthe results showed none of the metal constituent concentration was exceeding the permissible exposure limit ( pel ) for all plants . however , statistical analysis showed significant mean differences of lung functions between welders and non - welders . 
 the index was then applied to one of the welding industry ( plant 3 ) for verification purpose . 
 the developed index showed its promising ability to rank welding workplace , according to the multiple constituent concentrations of welding fumes that associates well with lung functions of the investigated welders.conclusionthere was possibility that some of the metal constituents were below the detection limit leading to  0 value of sub index , thus the multiplicative form of aggregation model was not suitable for analysis . on the other hand , maximum or minimum operator forms suffer from compensation issues and were not considered in this study .