A method for analysis of nuisance dust from integrated steel works: chemical and mineralogical characteristics of contributing sources

Integrated steelmaking is known to emit coarse airborne ‘nuisance’ dust (10–100 µm) to the production site and in the local environs. We present a method to quantitatively analyse the provenance, mineralogical and chemical attributes of the constituent particles in nuisance dust related to the integrated steelworks of Tata Steel, IJmuiden, the Netherlands. The dust is characterised per particle, using scanning electron microscopy with energy-dispersive spectrometry (SEM–EDS) microanalysis, and in bulk with quantitative X-ray diffraction (XRD) analysis. Based on mineralogical characteristics, particles in the dust are sorted into populations that can be related in detail to industrial processes and subsequent atmospheric weathering influence. The method is illustrated by application to a nuisance dust complaint sample from the neighbouring town Wijk aan Zee containing a large contribution of several dust sources from the integrated steelworks. Besides a background contribution from urban and natural dust, diverse sources from the integrated steelworks are identified in the nuisance dust sample, derived from coke-making, iron-ore agglomeration processes and blast furnace ironmaking, steelmaking slag processing, process fluxes, as well as steelmaking refractory materials. The most voluminous of these in the sample are directly verified by comparison with a set of reference source materials. The abundances, mineralogical and chemical attributes of the various dust particle populations in the sample are quantitatively examined including, specifically, the occurrence of the potentially toxic elements Mn and V. These elements occur with highest concentrations in dust derived from steelmaking converter slag: V is housed in dilute form (solid solution) in the phases di-calcium silicate and brownmillerite, and Mn chiefly in Mg–Fe-oxide (Mg-wustite ((Mg,Mn,Fe)O) and its oxidation product ((Mg,Mn,Fe)(Fe,Mn)2O4)). By treating nuisance dust as a particulate, multi-phase, multi-source material, the outlined method provides crucial information for toxicological evaluation and for mitigation of emissions, which is not obtainable by bulk chemical analyses alone. It also helps address the lack of adequate monitoring options for deposits of nuisance dust from integrated steel production, necessary to evaluate the relationship between deposition and monitored emissions that are regulated by the European Industrial Emissions Directive and by local permits based on this legislation.


Background
The dispersion of coarse airborne dust (10-100 µm) is experienced as a nuisance by those affected, as it soils surfaces where it is deposited, giving rise to the term 'nuisance dust' (used hereinafter to refer to this size fraction of dust). More importantly, depending on the nature of this nuisance dust, it can pose health-and environmental risks and damage property. Hence, the composition and extent of exposure to this dust need to be evaluated, certainly in the vicinity of industrial activity. Monitoring of nuisance dust is currently neither standardised nor widely practiced [1,2] unlike for the finer dust fractions PM10 (particulate matter < 10 µm in diameter) and PM2.5 (< 2.5 µm) as per European Ambient Air Directive [3,4] and standard EN 12341:2014 [5].
The health risks posed by nuisance dust (> 10 µm) are generally considered lower than those of inhalable PM10 and PM2.5 dust; however, dermal and oral exposure routes for nuisance dust ought to be considered. Nuisance dust can be inhaled but does not travel far beyond the uppermost respiratory tract, the oral/nasal cavity and larynx, into the tracheobronchial airway region [6]. By mucociliary clearance and coughing, the particles are brought back up to the oral cavity and they are subsequently swallowed [7,8]. If soluble along the digestive tract, the oral/nasal cavity, or in the lungs, dust constituents can pose health risks to the recipient [6]. Health risk evaluation should involve relevant particle-related metrics for dosimetry, such as solubility, particle count, size, surface area and composition [9].
Risks associated with dust contributions to the local soil should be assessed based on the dust constituents' bio-availability in the depositional environment [10][11][12]. To evaluate the risks associated both with digestion and soil contamination, it is insufficient to know only the inventory of elements present and their concentrations. Rather, the nature of the compounds needs to be known to evaluate the health and environmental impact of the dust.
This article presents a method to analyse nuisance dust regarding its constituent (mineral-)compounds and diverse contributing source materials. Our study relates specifically to integrated steel production where dust is associated with many processes, such as handling of raw materials (ores, fluxes and coals) and by-product streams (slags), high-temperature agglomeration (coke-making, sinter-and pellet production), as well as blast furnace ironmaking and steelmaking itself, see also, e.g. the study of Hleis et al. [13]. In the production environment, the risk of dispersal of nuisance dust is exacerbated by strong thermal updrafts associated with the processes at steelworks, which transports large particles to high elevations in open air with ample opportunity for dispersal by winds.
For risk assessment and effective mitigation, it is necessary to distinguish various dust types and their sources [14,15]. A methodology to do this has been developed at Tata Steel in IJmuiden, the Netherlands, using scanning electron microscopy and energy-dispersive spectrometry (SEM-EDS) and largely automated data processing. Based on mineralogical characteristics, the detected particles (grains) in a sample are classified under different populations, linked where possible with verified reference source materials. The abundances and chemical and mineralogical characteristics of these populations are then quantitatively analysed. This approach for analysing nuisance dust samples, and the type of information it yields, is discussed and illustrated here by application to a nuisance dust complaint sample from the town Wijk aan Zee (the Netherlands) containing a large contribution of several verified dust sources from the neighbouring Tata Steel IJmuiden integrated steelworks.

European legal framework for dust emissions
Dust emissions from industry in Europe are regulated by the Industrial Emission Directive [16]. According to this directive, industry applies for a permit providing the information on process and installation necessary for a country's local authority to set the permit conditions. The issued permit includes all the measures necessary to achieve a high level of protection for the environment as a whole (air, water and soil) set on the basis of Best Available Techniques (BAT) for emission prevention [16][17][18] with stipulation that the installation is properly operated. A permit includes emission limit values, or equivalent parameters or technical measures, and monitoring requirements. Limits for particulate matter in air are provided by the European Air Quality legislation [19]. Thus, with permits, the emissions from industrial production are regulated for a total mass allowance of emitted dust, and by legislation for concentration levels of PM10 and PM2.5, as well as for levels of arsenic, cadmium, mercury, nickel and polycyclic aromatic hydrocarbons.
In practice, for the integrated steel works of Tata Steel at IJmuiden, The Netherlands, the environmental permit limits the total emission of dusts from smokestacks as Keywords: Steelmaking, Coarse airborne dust, Nuisance dust, Characterisation, Analysis, X-ray diffraction, Microscopy, Mineralogy, Vanadium, Manganese well as from roofs/building tops with a maximum value currently set at 1750 ton/y total suspended particles (TSP) with separately defined limits for each process and operational dust scrubbing technique following BAT. The smokestack/roof-top emissions are quantified by continuous or regularly dust and PM10 monitoring. In 2019, the IJmuiden site emitted 688 ton TSP from these point sources, of which 481 ton PM10.
In addition, regarding nuisance dust, the European BREF Emissions from Storage [20] dictates the measures for raw materials, scrap and slag stockyards. In this regard, the Dutch Activities Decree applies to storage, handling and transportation. Currently, the permit is set for a maximum emission allowance of 1850 ton/y of total dust for the IJmuiden site. For these emissions, a methodology for quantification is used based on emission factors, size of storage areas, and activity intensity. Based on monitoring data, it is estimated that the site IJmuiden in 2019 emitted 1193 ton TSP from these open sources, of which 230 ton PM10. Emissions from the processing of steelmaking slags by a contracting company are regulated on a separate permit.
Total emissions (point, diffuse and open) from the IJmuiden site are publicly available information. PM10 is reported in the European Register EPRTR [21], with most recent entries from 2017. The TSP is reported in the Tata Steel Sustainability Report [22], with the most recent issue from 2018. Detailed data (Annual Emission Reports/Emissiejaarverslag) are reported to the Environmental Authorities in the Netherlands.

Rationale of nuisance dust characterisation
Nuisance dust is a particulate material containing particles derived potentially from multiple diverse sources-industrial, urban, agricultural and natural. Its characteristics reflect the emission potential of the contributing source areas [23], the transfer efficiency of these, and the retention/remobilisation tendency of particles after their initial deposition. Characterisation of nuisance dust must provide relevant information about its constituent particles and (sub-)populations of these that correspond to identifiable sources, preferably verifiable with reference source materials. The most relevant information concerns 1) establishing the provenance of particle populations, which is critical for mitigating emissions, and 2) evaluating environmental and health risks associated with these particles.
Commonly, the chemical composition of a whole particle is measured in SEM-EDS studies on dust [24][25][26]. This reflects, but does not explicitly describe, the combination of (mineral-)phases present at its surface, which can be (heavily) altered by weathering reactions. Moreover, an individual particle is a sub-sample commonly too small to be representative of the bulk chemical composition of a multi-phase source material. For these reasons, analysing the (mineral-)phase constituents of particles is preferable for establishing provenance since the chance of encountering a diagnostic combination of phases in a small and potentially weathered particle is still high. It is also preferable for understanding the form of binding of elements in mineral phases, which is crucial for establishing elements' bio-availability [9,27].
To characterise both the (surface-) chemistry and mineralogy of dust particles, our approach combines SEM-EDS spectral imaging (SI) and an in-house-developed software package, PARC (PhAse Recognition and Characterisation) [28], to convert SI data into phase distribution maps. The phase distribution maps have a spatial resolution of between 1 and 2 µm (depending on settings for data acquisition) and reflect the surface mineralogy of the dust particles since the SEM electron beam generates characteristic element X-ray signal only to very limited depths of a few µm into the sample [29]. Quantitative X-ray diffraction (QXRD) analysis of the dust sample (as mounted for SEM-EDS analysis) is used to provide an independent, complementary characterisation of the crystalline phases present in the bulk material, albeit not spatially resolved. The definition of phases in QXRD is crystallographic; whereas in the SEM-EDS-SI-PARC (hereinafter simply PARC) approach phases are chemically defined. Certain phases which are chemically difficult to separate with the SEM-EDS-SI approach, such as the Fe-oxides hematite and magnetite, are easily distinguishable and quantifiable with QXRD. Conversely, the QXRD analyses cannot identify and quantify the abundances of different amorphous phases which are chemically identifiable with SEM-EDS, nor very lowabundance phases (with detection limits of 0.1 wt% to 5 wt% depending on the crystal structure and crystal statistics). Combining techniques ultimately provides a useful cross-verification of the abundance of certain phases or phase-groups.
An additional strength of the PARC approach is that it can readily detect even extremely low-abundance phases occurring in a dust sample as long as (1) the phase is coarse enough to be spatially resolved by the measurements and (2) the phase can be defined and detected based on a unique combination of its own major constituent elements (see the section "Phase distribution maps from spectral imaging data" for more detail). Such trace phases may contain high concentrations of elements that are present only at very low concentration in the bulk sample. The detection limit with EDS of an element in a phase varies with matrix material (e.g. oxide/metallic etc.), surface topography and analytical conditions, but as a rule of thumb ~ 0.1 wt% is widely accepted. The PARC approach is, therefore, able to detect the occurrence of low-concentration elements in the bulk where these are present in a concentration of > ~ 0.1 wt% in one or more distinctive phases, which provides valuable information on the speciation of such elements.
The approach described here is, in short, to detect individual dust particles and to sort these into interpretable populations based on characteristic (mineralogical) phase combinations by processing SEM-EDS-SI measurements with PARC. This yields detailed, quantitative mineralogical, (phase-)chemical and morphological information on individual particles and particle populations and the abundances of these populations in the bulk sample are quantified.

Sample collection and preparation
Nuisance dust samples are collected related to complaints on-and off-site. Regularly cleaned surfaces such as car panels, windowsills and garden furniture typically form the substrate on which nuisance dust is noticed. Dust is collected from these surfaces with a fine brush and stored in a glass vial for transport to the lab. Samples constitute normally at most a few milligrams of dust. In the lab, the dust is sprinkled onto a double-sided adhesive conductive carbon disc, mounted on an aluminium stub. All subsequent analyses are performed on the original particles, i.e. particles are not sectioned and polished. The primary reasons for this are (1) the relative speed and ease of preparation, given the extremely small quantities of available material; (2) the fact that the particles are fully exposed, enabling shape and size quantification; and (3) the ability to examine the surface of particles, including thin covering layers which are barely resolvable in cross-section. Up to 26 dust samples can be combined in a sample holder with fixed positions for viewing. This facilitates automation and cross-correlation of data acquired on the same mounted sample using multiple analytical techniques such as 3D light optical microscopy, SEM-EDS-SI-PARC and XRD analysis.
Much of the discussion of the dust characterisation approach in this article is illustrated with a nuisance dust sample collected 26/01/2020 as part of an ongoing timeseries sampling campaign at a residential location in Wijk aan Zee, near the IJmuiden Tata Steel site. This specific sample contains a very high contribution of several dust sources from the integrated steelworks which are verified by comparison with reference source materials and which are encountered consistently in other nuisance dust complaint samples.

Raw data acquisition 3D light optical microscopy
Light optical microscopy (LOM) is performed on samples directly after mounting on an aluminium stub using a ZEISS SmartZoom5 (PlanApo D 5.0x/0.3 objective, with optical resolution of 0.92 µm). Since the particles have considerable topography, Z-stack imagery is acquired to ensure sharpness at different focal levels, referred to hereafter as 3D-LOM. Such an image is shown in Fig. 1a for an exemplary nuisance dust sample.

SEM and EDS spectral imaging
Samples are coated with carbon prior to SEM analysis. Back-scattered electron (BSE) images and spectral images (SI) are normally collected for several adjacent fields (typically 3 × 3) of 1024 × 768 (BSE) and 512 × 384 (SI) pixels with a pixel step-size of 0.8 µm, for BSE, and 1.6 µm for SI, although higher resolutions can be obtained. Fig. 1b shows the BSE image of the same area as the LOM image in Fig. 1a.
The SI data acquisition for the analyses reported here was performed with a JEOL 7001 FEG-SEM instrument equipped with a Thermo Fisher Scientific microanalysis system, running Pathfinder 1.3 software (same manufacturer), with two SDD/EDS detectors (crystal area 30 mm 2 each) positioned roughly diametrically opposite one another. This detector configuration is essential for SI data collection of 3D objects such as dust particles because EDS detectors are typically mounted at a 35° take-off angle, and signal is obtained only from the parts of the object's surface facing towards the detector. With detectors viewing from two diametrically opposite positions and combining the signals, this mitigates to a large degree the signal-shadow effects from sample topography. An accelerating voltage of 15 kV and beam current of 10 nA were used, with a total acquisition time per field of 30 min.
Quantitative chemical analyses were derived from EDS spectra with the Pathfinder 1.3 software in standardless mode, using the inbuilt filter fit method for peakfitting and PROZA (Phi-Rho-Z) correction method.  Table 1 and d PARC grain populations in Table 3. Area proportions (area %) are for the whole analysed area as shown. The {empty spectra} count for image c includes the background pixels of the sample holder. Area percentages for d are normalised to 100% excluding the background shown in black. Note, a small number of loosely mounted particles visible in a had physically moved prior to acquisition of images b-d. Full resolution images and accompanying tabulated data given in Additional file 1 Element concentrations are quantified and reported using assumed oxide forms (oxygen concentration is not independently analysed). Carbon is excluded from the quantification (note, for the sake of completeness, both carbon and oxygen net count values are reported together with those of the other analysed elements in Additional file 1: F).

QXRD
Dust samples are analysed with QXRD for crystalline compounds on the same microscopy mount but a larger area of about 4 mm 2 is examined. Except for the converter-slag reference dust samples, a Bruker D8 Discover instrument was used with a 2D position-sensitive area detector (Våntec 500) and Co Kα radiation (30 kV/40 mA) with sample rotation during measurement. In total, three individual frames were recorded to cover a range from 10 to 100° 2θ with a total measurement time of 30 min. The individual frames are merged. The XRD patterns for the converter-slag reference dusts were recorded with Bruker D8 Endeavour diffractometer equipped with position-sensitive detector (Lynxeye XE-T). The radiation used was Cu Kα (40 kV/40 mA). The advantage of this measurement procedure is that more crystals are found in diffraction condition than with a conventional Bragg-Brentano diffractometer. The data are converted into a standard diffraction pattern of angle (2θ) vs. intensity as presented in Fig. 2. Phase proportions in weight percent (wt%) are quantified using the Rietveld approach and subsequently converted into volume percent (vol%) using phase densities derived from the ICSD database [30] (see Additional file 1: A). Volume percentage provides greater ease of comparison with microscopy images that also represent volumetric information.

Data analysis
The following sections describe the processing of SEM-BSE-and EDS-SI image data with PARC for phase identification, particle detection (segmentation) and subsequent assignment of particles to different populations. A schematic overview of the raw data, processing steps and data output from PARC is shown in Fig. 3.

Phase distribution maps from spectral imaging data
The PARC software is used to convert spectral imaging (SI) data into phase distribution maps following methods described in more detail in van Hoek et al. [28]. The raw SI dataset contains the EDS spectrum for each pixel in an image field. The software sorts all pixels into groups, hereafter referred to as PARC groups, based on the  Fig. 1), showing raw and modelled diffraction pattern, left, and quantitative phase analysis, right. Phase proportions are shown in terms of mass (wt%, with 2σ uncertainty) and volume (vol%). Coloured bars on right-hand side of table show typical association of phases with key source materials as outlined in Table 2. *Q-C-F-M = quartz-clayfeldspar-mica source/material category as in Table 2 detected combinations of spectral peaks. A peak is registered if its energy and intensity are above user-defined threshold values. The pixel sorting results in groups of unique peak, i.e. element, combinations. The distribution of these groups is shown in a colour-coded phase map (as in the example in Fig. 1c) alongside the SEM backscattered electron (BSE) images (Fig. 1b), for a single field or for multiple stitched fields. Area percentages of all groups are quantified (as displayed in the key to the phase map in Fig. 1c). The relationship between area and volumetric proportions is discussed in the section "Quantification of volumetric contributions". The naming convention used in this article for PARC groups is {p1, p2, …, pn}, referring to the spectral peaks by which the group is defined, and listed in order of increasing energy, i.e. a purely descriptive label rather than an interpretative one in terms of true mineralogical phases with which the group may correspond. The PARC groups occurring in materials discussed in this article are listed in Table 1. The chemical composition of a PARC group can be quantified using the sum spectrum of all pixels belonging to the group. Marginal pixels which are bordering with other groups in an image, where their spectrum is likely contaminated by signal from the adjacent phase, can be excluded in the quantification by applying a pixel-erosion filter. The thus-cleaned sum-spectra can be exported by PARC and converted to quantitative chemical analyses as described in "Raw Overview of raw (spectral-)image data acquired, processing and data generated by PARC. Processing steps are indicated with arrows and italic text in the flow-chart. Numbers in parentheses indicate relevant text sections providing more detail as follows: 4.2 = "Raw data acquisition"; 5 = "Data analysis"; 5.1 = "Phase distribution maps from spectral imaging data"; 5.2 = "Confirmation of phase identity with QXRD"; 5.3 = "Defining grain populations"; 6 = "Results and discussion"; 6.1 = "Source apportionment of grain populations". References to associated figures and tables are given. Key to outlines in flow-chart: single solid line = image/image stack; double solid line = segmented image (map); single dashed line = derived information. Data produced per PARC group are listed under the relevant grain/grain population/field mask used to collate pixels. The number of PARC groups, grain populations and source/material categories in the processing is indicated (including unclassified). Number of pixels and grains varies per dataset and the numbers shown are approximate indications only. 3D-LOM is included in the overview for completeness but is not further used in the data processing described here data acquisition" section and yield compositions with reasonable correspondence to known phase stoichiometry, sufficient for phase identification. Locally, if two or more phases are finely intergrown on a length-scale smaller than the spatial resolution (1 µm excitation volume), a mixed spectrum is generated that cannot be distinguished from a single phase. Such mixed pixels are nonetheless assigned to a PARC group based on their characteristic peak combination and, with appropriate account taken of this potential artefact, can also be used as a diagnostic feature for particle populations.
With grouping based on peak combinations only, separate phases containing the same major elements in different proportions will not be distinguished and are classified as one PARC group. The possibility exists within PARC to further sub-divide such groups based on continuous variations in element-peak intensity (ratios) using so-called "density plots" [28]. However, with 3D objects like dust particles, these intensities are generally too strongly affected by topography to robustly apply this approach. Hence, the groups here are based purely on element-peak combinations. As a result, there is a bracket of plausible/commonly occurring true phases corresponding with a PARC group within a dust sample or a given particle.
Peaks with energies below 0.9 keV are excluded here when defining PARC groups, due to (1) overlap of the K-line peaks of the light elements C, N, O and F with the L-line peaks of first series transition elements; and (2) strong absorption effects on low-energy signal due to the topography of the 3D particles. Consequently, no direct distinction is made between carbonates, hydrates, and various oxidation states of elements. The generally much higher intensity of the Fe Kβ peak in spectra measured in (verified) metallic Fe is, however, exploited to help distinguish dust particles rich in metallic Fe as opposed to FeOx.
The thus-defined pixel groups are saved in a PARC group model and our approach is to apply one standard model consistently to all datasets on nuisance dust samples (which may be periodically updated and re-applied to old datasets). All PARC groups with more than 1% area contribution to any dust particle (grain) population in the exemplary nuisance dust sample are shown in Table 1 (full listing of all PARC groups is given in Additional file 1: D). The spectral image converted to a PARC phase map is given in Fig. 1c.
A large number (around half the total) of PARC groups are defined by the presence of either chlorine-, sulphur-, or phosphorous peaks in addition to other elements. These {± , Cl, ±}, {± , S, ±} and {± , P, ±} groups are mostly not identifiable as genuine chemical compounds and they result from a thin layer of contaminant chlorides, sulphates or phosphates on top of a grain, giving mixed spectra. Table 1 shows the interpreted correspondence between PARC groups and phases detected by XRD in the example nuisance dust and in reference source materials (see the section "Reference source materials"). The PARC groups with the highest abundances in these samples all correspond clearly with one or more phases detected by XRD. PARC analyses in the samples shown here, and more generally in our analyses of nuisance dust samples, show higher abundance of certain phases than obtained from quantitative XRD analysis. This is consistently the case for the following PARC groups: {Si, Ca} and {Si, P, Ca} compared with XRD C2S phases; {Ca} compared with XRD calcite/lime/portlandite; {Na, Cl} compared  with XRD halite (NaCl); and {± , S, ±} PARC groups versus gypsum/anhydrite. In all these cases, we interpret the discrepancy as reflecting one or both of the following: (1) PARC analysis detecting material which is amorphous and, hence, not detected by XRD; (2) an over-representation of phases occurring dominantly in thin surface layers on grains in the PARC volumetric estimate, while these are barely detected with XRD due to the deeper penetration of X-rays c. 20-30 µm into the grains' interiors and the small contribution of the surface layers to the bulk grains. The mutual proportions of the Ca-ferrites, ironoxides and quartz are much more closely in agreement between PARC and XRD analyses, suggesting that these phases are fully crystalline and that they are not overrepresented in PARC analyses by occurrence as thin surface layers on grains.

Confirmation of phase identity with QXRD
Besides discrepancies due to amorphous material and thin surface layers, the PARC and XRD approaches bring complementary information on the samples. One specific example of this concerns the diverse occurring Ca-ferrites, which in PARC are detected dominantly as the group {Ca, Fe}, while they are readily distinguished with XRD analysis (Table 1, Fig. 4). In iron-ore sinter, the group of Ca-ferrite minerals referred to collectively as silico ferrites of calcium and aluminium (SFCA) is abundant and diagnostic of this material (they do not occur in other industrial source materials). In contrast, the mineral brownmillerite (C2(A,F)) is the predominant Ca-ferrite phase in converter-slag reference material. These distinctive mineralogical signatures can be used to refine and independently support estimates of source material contributions to nuisance dust. Overall, the XRD analyses provide independent evaluation of the phases present in a sample and help confirm the detection of source material contributions to nuisance dust samples based on diagnostic phase occurrences.

Defining grain populations Grain segmentation
Since the dust characterisation approach is based on individual dust particles, these need first to be recognised and defined, i.e. segmented, as individual objects (grains) in the acquired (spectral-)imagery. Grain segmentation is currently performed on the SEM-BSE image by defining a background grey level to produce a binary mask, followed (where required) by a sequence of pixel-erosion and -dilation steps, and application of a watershed-criterium operation identical to routines in ImageJ [31]. The resultant grain segmentation mask is applied to all correlated images, i.e. SEM-BSE, SI and derived PARC group maps, as well as the quantitative (spectral) information contained therein. The grain segmentation for the exemplary nuisance dust sample is visible in Fig. 1d.

General approach to grain classification
Following grain segmentation (as listed in Fig. 3), data are generated per grain on its size and shape, sub-grains (spatially discrete domains comprising one PARC group), overall chemical composition, constituent PARC group area proportions and their compositions. A (self-developed) decision tree comprising sorting criteria (filters) is then used within PARC to classify individual grains under defined populations (Fig. 5). The filters used in this article are based on the area proportions of certain PARC groups in each grain (including summation, mutual ratios, etc.), which reflect genuine mineralogical phases (and phase intergrowths). The grain size/shape and sub-grain information are not discussed further in this article but can potentially be used to further sub-divide grain populations based on more than solely PARC group abundances (mineralogy). Each filter contains one or more thresholding rules where the result of a given filter is 'true' if and only if all the rules within a given filter are satisfied. Individual grains follow a path through the decision tree filters until they are classified under a population. The decision tree and filters together comprise the grain population model, which is applied consistently to all nuisance dust samples. The model is constructed and refined manually based on a combination of the following approaches: 1) Using characteristic/diagnostic properties of known potential source materials for nuisance dust emission 2) Identifying coherent populations, and their characteristic properties, within nuisance dust samples, which may or may not yet be confidently attributable to specific source materials Table 2 lists a number of interpretive source materials and distinctive material types identified, referred to henceforth as 'source/material' categories, to which grain populations are assigned (Table 3). Some categories comprise groupings of several grain populations, while others correspond directly with a single grain population. Their key defining mineralogical characteristics are outlined together with the translation of these to PARC group (pixel) abundances as measured in this study. Moreover, independent mineralogical characterisation of a number of potential source materials from XRD analysis (Fig. 4) is used to support the definitions in Table 2. Several key populations can readily be cross-identified by simply inspecting the PARC group maps of both nuisance dust samples and reference samples of potential source materials. Examples include the populations assigned to source/material categories 1-5, 7-8 and 10-13 (Table 2). These are the starting point in defining an initial set of populations, filters and decision tree structure for the grain population model. The grain population model is subsequently expanded and refined in an iterative manner by applying it both to reference samples of potential source materials and nuisance dust samples, assessing the results of the grain classification, modifying the model manually and repeating the cycle. The key aims are to (1) maximise the correspondence between the defined grain populations and specific reference samples where applicable (the evaluation of this is discussed further in the section "Application of grain population model to reference source materials and example nuisance dust sample") and (2) minimise the amount of unclassified grains in nuisance dust samples by defining coherent populations even where no (direct) match is found with any reference source material. The model is regularly updated in response to the results of applying it on new nuisance dust and potential source material samples. Hence, the grain population model shown in this article  Table 3. Arrows indicate flow direction of connectors between nodes. Key shows assignment of populations to interpretive source/material categories listed in Table 2. Note, alphabetic labelling of populations reflects order in which these appear in decision tree and hence deviates from order in key (See figure on next page.) Table 2 Interpretive source/material categories to which grain populations have been assigned, together with diagnostic or common true phase assemblage and translation to PARC groups (as listed in Table 1    is necessarily a snapshot of its state of development at time of writing. The main advantages of a manual decision tree approach are related to the transparent, 'white box' nature of the model, which allows full insight into how particles are classified. Disadvantages of this approach include the possibility of overlooking genuine but less-strikingly obvious populations while also potentially generating artefact populations by defining them too broadly or too tightly. The PARC data output lends itself readily to deeper statistical analysis and further work is in progress on interrogating the data, including on strengthening grain population classification and correspondence with source materials.

Reference source materials
Ultimately, the grain population classification approach requires validation on well-defined reference materials where the provenance of the material is independently known. Several reference samples are discussed in this article: • Dusts from conveyor belts transporting iron-ore sinter (sinter dust 1) and iron-ore pellet (pellet dust 1) material from their respective production sites to the charging point of the blast furnace; • Dusts sampled at the sintering plant (cooling section; sinter dust 2) and pelletising plant (abrasion-testing facility; pellet dust 2); • Metal-rich slivers from a scrap processing facility (uilenballen); • Dust from a scrap park location (scrap park); • Three dust samples (sieved to < 125 µm) taken from the pits used to discharge and cool converter slag prior to further processing (converter-slag pit dusts 1-3).
The setup and refinement of the grain population model were based to a large extent on using one of each type of the reference materials listed above, while in the case of the sinter, pellet and converter-slag reference materials, additional similar samples of each were treated purely as unknown materials for validation purposes (samples 2 and 3 of each material type in the list, where applicable). Figure 4 shows the volume percentages of phases from XRD analysis for the above reference materials, alongside the example nuisance dust sample (full tabulated data and indexed diffraction patterns are given in Additional file 1: A). The reference samples pertaining to sinter-, pellet-and converter-slag pit dust are clearly distinguished by their bulk mineralogical compositions. These mineralogical characteristics are taken directly into account in Table 2 for defining the diagnostic/ common true phase assemblage for source/material categories (1) converter slag, (3) iron-ore sinter (sintered material) and (4) iron-ore pellet (fired material). The results of applying the grain population model to the above listed reference materials are shown in the section "Application of grain population model to reference source materials and example nuisance dust sample". Figure 5 shows the full decision tree structure for the developed grain population model applied here. Table 3 lists the 49 grain populations shown in the flow-chart and provides brief text descriptions of the defining characteristics of each population, the preferred interpretation of each population and assignment to the 26 source/ material categories described in Table 2 (including that of unassigned grains). Complete filter definitions for Fig. 5 are tabulated in Additional file 1: B.

Quantification of volumetric contributions
Ultimately, the volumetric and mass contributions of nuisance dust particles are of interest, while additionally the sizing and (specific) surface area of particles are highly relevant for understanding environmental behaviour and bio-availability of toxic species. The nuisance dust particles (grains) have very variable shapes and tend to rest on their largest cross-section in the mounted sample. Translating the directly measured silhouette area fraction (reported in the figures and tables here) to volume fraction without any further stereographic correction implicitly assumes perfectly plate-shaped particles of equal thickness (rising above the sample mount). This is clearly not a realistic assumption for more spherical particles, while even genuinely platy-shaped particles have a wide range of aspect ratios and are not guaranteed to be oriented resting on their largest flat face. This is currently the main source of error in estimating volumetric contributions of grains and grain populations to the bulk samples. Our expectation is that (ongoing) further quantitative exploitation of 3D microscopic imagery will substantially improve estimates of volumetric and surface area contribution from different particle populations. The combination of QXRD and PARC analyses, meanwhile, provides an important internal consistency check regarding the volumetric proportions of mineralogical phases and can help flag up the occurrence of some phases disproportionately concentrated in thin superficial layers on particles.

Source apportionment of grain populations
Broadly speaking, the grain populations can be split between IJmuiden steelworks-derived, urban and natural/agricultural contributions in the dust. The source/ material categories listed in Table 2 most unequivocally attributable to the steelworks are 1-6 and 9-14, since for these, there are no plausible alternative sources to steelmaking raw materials (ores, fluxes), intermediate products (iron-ore agglomerates), by-products (slags) and highly specific installation materials (refractory materials). Source/material categories 6-8, 15 and 16 are found in samples collected on-site but are not yet uniquely attributed to a specific process or raw material. Ca-carbonate (most commonly calcite) is a natural rock-forming mineral which can occur as seashells in the local natural environment, while as limestone it is a raw material (flux) for steelmaking processes. It can also form as a weathering product of converter (and other calcic) slag. As such category 17 (Ca-carbonate/ (hydr)oxide) is a potentially ambiguous category of dust particles, except in cases where steelmaking-diagnostic phases are observed in conjunction with these particles. Currently, the carbon-rich grain populations are not optimally sub-divided and source/material category 24 may include coal, cokes, graphite and peat and pollen (natural). Urban dust sources include various building materials, glass, and paints: cement (category 18) and paint (category 21, potentially category 20) are identified as potential source materials for grain populations in the presented model based on widely known characteristics of these materials [32,33]. Gypsum (category 19) can be both urban building material and directly steelmaking process/site related. Natural dust sources include sand, rich in quartz with lesser quantities of feldspar, and weathering products such as claywhich are currently grouped together under the source/ material category 22. The surface contaminant NaCl is traceable to sea spray and has been observed partially covering the surface of diverse particles.

Application of grain population model to reference source materials and example nuisance dust sample
Following classification of grains to populations, data are generated per grain population as specified in Fig. 3. Here, we focus on the relative abundances of grain populations, their mineralogical makeup and (micro-)chemical attributes, which are of relevance for further understanding their provenance as well as the distribution of specific (potentially toxic) chemical elements on the population and phase level. Table 4 reports the area proportions of all defined grain populations for the example nuisance dust sample of Fig. 1 and for the reference source material samples referred to earlier in this section. The area proportions of the interpretive source/material categories, to which the populations are assigned, are shown in Fig. 6 (tabulated in Additional file 1: C; original SEM-BSE and PARCgenerated images are given in Additional file 1: D, along with tabulated area percentages of all PARC groups per sample). This set of results gives an important indication of the performance of the grain population model in distinguishing three key potential source materials-converter slag, iron-ore sinter and pellets-for which their respective particle (grain) populations are expected to have overlapping apparent mineralogy in terms of PARC group abundances. It can be seen that there is a very good correspondence between the defined grain populations and their respective interpretive source materials, with minimal overlap of 1% or less between the converter-slag and sinter dust samples and < 5% overlap between sinter and pellet dust samples, where overlap refers to misallocation regarding source such as 'pellet' grains in a sinter reference dust or vice versa. Similarly, the discrimination between particles rich in FeOx, prevalent in ironore sinter and pellets, and Fe-metal-rich particles such as encountered in the uilenballen sample, is also strong. Only around 2% or less (by area) of the pellet dust is classified as potentially metal rich (source/material 7), while in the uilenballen sample > 90% is classified this way. The remainder of the uilenballen grains are dominantly classified as source/material 5, with around 7% area contribution, though they are upon detailed inspection not considered to be ore/sinter/pellet related. Hence, there is in this grain population clearly some potential overlap in the classification of grains that are rich in FeOx from different sources. The scrap-park sample contains a significant population of grains (source/material 8) rich in Zn-oxide/Zn-Fe-oxide, while a small amount (c. 1% area) of clearly similar grains is also found in converter-slag pit dust sample 2. The scrap-park sample is a case in point of the challenges in finding appropriate reference samples for source materials: while the sampling location is clearly defined, the material collected there is highly heterogeneous, as can be seen in the lack of any single dominant grain population or source/material category. Detailed examination of this sample confirms that the grains classified as ore/sinter/pellet (5), pellet (4), sinter (3), converter (1) strongly resemble their counterpart populations in their respective reference source materials, while a significant portion of the material appears to be natural silicates (22: quartz-clay-feldspar-mica). The  Key to population labels and explanation in Table 3   Table 4  relatively large area fraction of unassigned grains contains numerous diverse silicate phases also considered likely to be of natural origin in most cases. The only grain populations considered likely to be associated with scrap material are those falling under source/material category 8 (Zn rich). The interpretation is that this sampling location captured a heterogeneous population of dust particles of which only small portion is derived from the material hoped to be sampled there.
Another key source material is not verified directly with regard to a reference sample here-(10) Ca-aluminate slag-since this source material is mineralogically distinctive and indicative of a well-defined set of steelmaking process slags which are mixed and processed together at their suspected emission point.

Mineralogy of example nuisance dust sample
Applying this grain population model, about 90% of the analysed area of the nuisance dust sample is found to be composed of source/material categories 1 (converter), 2 (slopping slag), 3 (sinter), 4 (pellet), 5 (ore/sinter/pellet) and 10 (Ca-aluminate slag); while converter-and sinterrelated populations are clearly the two dominant contributors (Fig. 6). Figure 7 provides a graphical summary of the mineralogical makeup of the bulk sample and the diversity of its constituent grain populations with > 0.5% area contribution (full data are tabulated in Additional file 1: E).
Clearly, the PARC group proportions chiefly reflect the filters used to classify the grains under their respective populations. Nevertheless, the distinctive characteristics of a grain population regarding its constituent PARC groups help to understand its provenance and its potential health implications in terms of distribution of hazardous elements in the constituent mineral phases. A few points merit some discussion here.
The converter-slag populations (source/material category 1) vary widely in their PARC group abundances in Fig. 6 Area percentages of interpretive source/material categories for the example nuisance dust sample together with a number of reference source materials from the Tata Steel works in IJmuiden. Key shows source/material number and shortened name-see Table 2 for complete names and descriptions. (Tabulated data given in Additional file 1: C) line with the series of filters used in each case. Essentially, the populations are split into three types of relatively 'pristine' converter-slag particles based on the predominance of each of the three major converter-slag phases C2S (or C3S; populations y and z), C2(A,F) (populations aa and ab) and Mg-wustite (populations ac and ad); plus 'weathered' converter-slag material (w and x) characterised by the presence of distinctive P-and Al-enriched (hydrated/carbonated) silicate phases formed by alteration of the original slag phases. These populations are further sub-divided based on the presence/absence of significant abundance of sulphate phases on-at leastthe surface of the grains, which is clearly seen in the PARC group proportions in the figure. Compared with the reference dust materials sampled in the converterslag cooling pits, the converter-attributed grains in the nuisance dust sample are markedly different in their detailed breakdown into populations and their mineralogy. Closely inspecting Table 4 reveals that the grains in the reference converter-slag pit dust samples are overwhelmingly classified as population y, i.e. pristine C2S/ C3S-dominated grains with no significant covering of sulphate, followed by minor (< 5%) population aa (C2(A,F)rich grains) and a small but clearly verified abundance (upon detailed inspection) of the 'weathered' slag population w. The example nuisance dust sample contains much more weathered slag population (w and x) as a fraction of the total converter-attributed grains and, moreover, a large portion of all the different converter-attributed grains have a significant covering of sulphate group pixels. While these weathered-and sulphate-covered grain populations are scarce in the reference converter-slag pit dust samples, they are still clearly recognisable there as being equivalent to the grains encountered in the nuisance dust. The preferred interpretation is that the converter-attributed material in the nuisance dust sample is, on the whole, significantly more weathered as well as having more sulphate covering it than the relatively pristine material encountered in the three reference dust samples from the cooling pits. Such detailed, nuanced information is being used to glean a better understanding of the formation mechanisms of the dust from converter slag and its further propagation into the surroundings.
Annotations 1 and 2 in Fig. 7 highlight the abundances of specific PARC group pixels chiefly associated with mixing of two or more true phases where these are finely intergrown with one another. The converter-slag populations contain abundant {Mg, Ca, Fe}  Tables 2 and 3, respectively. Source/ material short names shown here. PARC groups labelled as Table 1. PARC groups are clustered (see brackets) based on interpreted mineralogical association and ordered according to the abundance of these clusters in the bulk sample. Mineralogical identity/interpretation of PARC groups labelled in figure. Annotations 1 and 2 discussed further in text (see section "Mineralogy of example nuisance dust sample"). (Tabulated data given in Additional file 1: E) pixels (annotation 1), which is a characteristic feature of converter-slag reference materials and has been used to strengthen the filters to discriminate between converter-and sinter-related grains (hence, the results here strongly reflect that filtering). These are interpreted as mixed pixels reflecting finely intergrown Mg-wustite {Mg, Fe} and C2(A,F) {Ca, Fe}. Both converter-slag and iron-ore sinter/pellet populations contain an abundance of {± Mg, ± Ca, Fe, Si} pixels (annotation 2), which reflect different fine phase intergrowths in each case. For the converter-slag grain populations, these pixels represent intergrowth of C2S/C3S {Si, Ca} with Mg-wustite {Mg, Fe}. For the sinter/pellet grain populations, they represent mixed pixels of variously C2S, silicate glass, Mg-ferrite, Mg(-Ca)-olivine and diverse Ca-Mg-Fe-silicates, which occur abundantly in sintered or pelletised material but not in unreacted ore material. Population al (labelled 5-al in Fig. 7) is notably lacking in these pixels, consistent with it being largely unreacted ore material.

Quantified chemistry of grain populations
Quantitative chemical analyses per grain population have been derived for the example nuisance dust sample by processing their respective EDS sum-spectra (as described in the section "SEM and EDS spectral imaging"). The results are shown in Table 5 and Fig. 8 for grain populations contributing 0.5% or more (shown by the superimposed line in the figure). These chemical analyses are broadly in line with expectations for the respective source/material categories. The converter-related populations (source/material 1) vary rather widely, which in the case of populations y and z (C2S/C3S rich), aa-ab (C2(A,F) rich) and ac (Mg-wustite rich) chiefly reflects the relative abundance of the three main converter-slag phases. The 'weathered' converter-slag populations w and x resemble the C2S/C3S-rich populations y and z, while deviating in the Al and P concentrations-reflecting the abundance of the interpreted weathering product phases represented by the PARC groups {Al, Si, Ca} and {Si, P, Ca}, respectively. The sulphur concentration of the converter-related populations varies strongly according Table 5 Quantitative chemical analysis from EDS sum-spectra per bulk grain population for populations with > 0.5 area % contribution, reported as oxide (or element) wt% and normalised to 100% of the elements above detection limit (values > 3-sigma error in quantification result) Elements are reported with assumed oxide forms (except Cl). See Tables 3 and 2 for key to, respectively, populations and source/material to whether the population was defined as being 'sulphate rich' or not. Clearly, the bulk chemistry of the grains comprising the various converter-related populations varies very widely, corresponding with large variations in the mineralogy of individual grains. This illustrates the difficulty of classifying grains based only on their bulk chemical composition, versus the ability to draw together chemically disparate grains based on their related mineralogy. A key additional reason for basing grain classification on mineralogical phase proportions (via PARC groups) is the commonly heavy contamination and alteration of the surface by/to other phases. This may lead to significant deviations in chemistry between the analysed surface and the true bulk chemistry of the grains. By classifying grains based on characteristic proportions of specific PARC groups, it is possible to screen out likely surface contaminant and alteration products. These are classified as distinct PARC groups and can be excluded from the expressions in the filters. Based only on bulk chemistry per grain, such surface artefacts are much more difficult to screen out.

Mode of occurrence of identified potentially toxic elements-V, Mn
In a recent study [34] by the Dutch National Institute for Public Health and the Environment (RIVM) on dust deposits in the town Wijk aan Zee, three elements V, Mn and Pb were detected with potential exposure health risk. Their occurrence was not coupled to any specific source of dust other than that implied by the vicinity to the steel production facility. To estimate the health risk of dust particles, the dosimetry of its constituent elements needs to be determined [9]. For the chemical analysis as performed by the RIVM study, aqua regia digestion will likely dissolve most material, with the exception of quartz and other relatively insoluble minerals. The elemental composition analysed in this way can serve to determine a maximum health risk in the scenario that all elements were to become available in solution [34]. In the case of human ingestion of dust, dissolution of coarse particles may occur to a lesser extent during passage through the digestive tract than during the aqua regia preparation used for analysis. For those phases hosting V, Mn and Pb, knowing their solubility and dissolution kinetics in  various environments is important to come to a realistic health risk assessment.
In the present study, dust particle populations are identified with mineralogy and phase chemistry specific to them, providing direct information on which dust sources contribute the potentially toxic elements mentioned above and on their mode of occurrence. No Pb was detected in the EDS sum-spectra of any PARC group. V and Mn, however, were clearly detected in specific dust particle populations. Figure 9 compares the analysed V and Mn concentrations of dust grain populations with > 0.5 area% contribution to the example nuisance dust sample (Table 5). To shed light on the specific mineralogical phases containing high concentrations of these elements, Table 6 lists all PARC groups with concentrations > 0.5 wt% V reported as V 2 O 5 and > 10 wt% Mn reported as MnO quantified for their sum-spectra (collated per group per population, with full data tabulated in Additional file 1: F; note that listed populations in the table need not have a large bulk concentration of V or Mn if they contain only minor amounts of those PARC groups with the high concentrations).
Clearly, the highest bulk concentrations of V are measured in the populations associated with converter slag as well as slopping slag (also converter process derived). Significant concentrations are also measured in the grain populations associated with dolomite (11-n 1 ) and Ca-(hydr)oxide/-carbonate (17-t, 17-u). Iron-ore sinter and pellet grain populations (source/material categories 3, 4 and 5) have much lower concentrations around detection limits, while no detectable concentration is measured in the Ca-aluminate slag population (10-af ). Within the converter-slag grain populations with the highest overall V concentrations, the PARC groups with the highest V-concentrations in the {Si, Ca} (0.9-3.9 wt% V 2 O 5 ) and {Al, Ca, Fe} (0.5-1.2 wt% V 2 O 5 ) groups, interpreted, respectively, as C2S and C2(A,F) and weathering    Table 5. Left-hand axis = wt% component; right-hand axis = area % population (shown as superimposed line). Upper plot, V reported as V 2 O 5 wt% and, lower plot, Mn as MnO wt% products of these. This is consistent with earlier studies showing a high affinity of these phases for V [35]. Significant concentrations of V (0.9-1.5 wt%) are also measured in the {Si, Ca} group-interpreted as C2S-in the sinter-related grain populations, while the bulk concentration of V in these populations is much lower than in the converter-related populations. This corresponds with a much lower abundance of C2S in the sinter material compared to converter-slag material. Other grain populations also contain individual PARC groups-hence specific phases or phase intergrowths-with V concentration of V 2 O 5 > 0.5 wt% (listed in Table 6), including potentially other phases besides C2S and C2(A,F), though the V concentration in the bulk grain populations is relatively low. The grain populations with the highest measured concentrations of Mn are those attributed to converter slag and rich in C2(A,F) (1-aa and 1-ab) and especially the Mg-wustite-rich population (1-ac). Notably, the most Ca-ferrite (SFCA-)rich iron-ore sinter grain population  (Fig. 4) show that these contain both wustite and magnetite (spinel solid solution) but no direct estimate can be made with QXRD of the mutual proportions of these in the converter-slag dust particle populations within the nuisance dust sample since the latter contains a mixture of populations and abundant magnetite from iron-ore, sinter and pellets.
In the case of the sinter population (3-ao), an explanation is required for the relatively high bulk Mn concentration compared both with the other iron-ore-related grain populations and the chemistry of bulk sinter material (which typically contains c. 0.1-0.5 wt% MnO depending on the ore blend). This grain population contains a low but clearly detectable presence of PARC groups (listed in Table 5) with Mn concentrations greatly exceeding those of any groups encountered in the converter-slag grain populations. Detailed inspection of these Mn-rich pixel groups' distribution in the PARC phase map (Fig. 10a) reveals two individual grains (annotations 1 and 2), classified as sinter (3-ao), containing a high abundance of {± S, ± Ca, Mn, ± Fe} pixels together with the groups typically abundant in this grain population. The observed microstructural relations of the PARC groups (shown in Fig. 10b  To summarise, applying the PARC methodology indicates that the most V-rich material in the example nuisance dust sample is derived from steelmaking converter slag (and process) and that it is housed in dilute form (solid solution) in the phases C2S and brownmillerite, as well as weathering/alteration products of these. Similarly, the same material derived from converter slag has  Fig. 10 a Region of interest in SEM-BSE image of nuisance dust example (Fig. 1) showing distribution of PARC groups {± S, ± Ca, Mn, ± Fe} highlighted yellow. Note localised concentrations of these pixels in grains 1-3. Grains 1 and 2 are classified as population ao (iron-ore sinter rich in Ca-ferrites), grain 3 as population aw (unassigned, low empty spectra; Table 3 Both Mn and V are known to be heavily concentrated in converter slag during the production of steel from original raw materials (ores, fluxes, coal). An unexpectedly high Mn concentration detected in one of the ironore sinter grain populations flagged up the presence of scarce particles (< < 0.1% area contribution) containing abundant Mn-Fe-oxide (or hydroxide/carbonate), where this phase contains at least c. 80 wt% Mn reported as MnO. Further investigation is required to fully pin down the identity of this phase but it is clearly associated with material from the iron-ore sintering process and present in very low abundance in sinter reference materials. This example illustrates how with a micro-analysis approach-based on treating individual particles and particle populations in terms of their constituent mineral chemistry-it is possible to investigate simultaneously the mode of occurrence, distribution and origin (source) of potentially toxic elements occurring in bulk nuisance dust material. While bulk chemical analyses of nuisance dust may flag up any elevated concentrations of potentially toxic elements in the bulk material, an adequate characterisation must go beyond this and provide at minimum the type of information shown here on constituent particles' mineralogy and phase chemistry to support toxicological evaluation, establish provenance and mitigate emissions at source.

Conclusions
A robust method has been developed to quantitatively analyse the provenance, mineralogical and chemical attributes of the constituent particles in (10-100 µm) coarse airborne nuisance dust deriving from the integrated steelworks of Tata Steel IJmuiden, the Netherlands. The method involves detecting individual dust particles (grains) in a sample and sorting these into populations based on the mineralogical characteristics of their surfaces as analysed with SEM-EDS spectral imaging and processing with PARC software. The abundances, and mineralogical and chemical characteristics, of each sample's grain populations are then quantitatively analysed. QXRD provides independent complementary analysis of bulk sample mineralogy. A total of 48 grain populations have been defined based on comparison of nuisance dust grains with the reference characteristics or samples of 23 source materials of industrial (steelworks) (14), urban (3), natural (2) and potentially ambiguous origin (4), while a further two distinct materials are identified whose origin is not known. The method is demonstrably successful at distinguishing several key source materials derived from iron-ore agglomeration and the steelmaking converter process. Applying the characterisation method to a nuisance dust complaint sample from the town of Wijk aan Zee with a very high contribution from several verified source materials attributable to the steelworks, we demonstrate that detailed information can be gleaned pertaining to provenance and potential environmental/health risks associated with the encountered dust particles. The occurrence of potentially toxic elements can be traced to specific dust particle populations of known origin and the mode of occurrence in (mineral) phases can be identified to support robust toxicological evaluation. Specifically, for the nuisance dust sample shown, the most V-rich material is derived from steelmaking converter slag (and process) and is housed in dilute form (solid solution) in the phases C2S and brownmillerite, and in weathering/alteration products of these. Similarly, the same material derived from converter slag has the highest general concentrations of Mn, housed chiefly in the Mg-Fe-oxide solid solution phases Mg-wustite ((Mg,Mn,Fe)O) and magnetite ((Mg,Mn,Fe) (Fe,Mn) 2 O 4 ). The developed method is currently being applied to characterise larger sample sets of nuisance dust in the surroundings of the integrated steelworks of IJmuiden to aid efforts to mitigate nuisance dust emissions at source and provide more detailed and appropriate information for toxicological evaluation than that available from solely bulk chemical analyses. Moreover, the characterisation approach also helps address the lack of adequate monitoring options for deposits of nuisance dust from integrated steel production, necessary to evaluate the relationship between deposition and monitored emissions that are regulated by the European Industrial Emissions Directive and by local permits based on this legislation.
Additional file 1. Supplementary information is provided in the following additional files (prefixed with 'additional file'): (A) fully indexed diffraction patterns and tabulated quantitative results for QXRD analyses; (B) complete details of filters used in grain population classification model; (C) tabulated data plotted in Fig. 6; (D)-(i) full tabulated data for PARC group and grain population area percentages per sample, and (D)-(ii) original SEM-BSE and PARC-generated imagery for all samples (files labelled by sample name and, where relevant, by image type); (E) tabulated data as plotted in Fig. 7; (F) tabulated quantitative SEM-EDS analyses as applicable to sections "Quantified chemistry of grain populations" and "Mode of occurrence of identified potentially toxic elements-V, Mn" (PARC groupand grain population sum-spectra for example nuisance dust dataset).

Abbreviations
BAT: Best available techniques; BREF: Best available techniques reference document; BSE: Back-scattered electron; EDS: Energy-dispersive spectrometry; LOM: Light optical microscopy; PARC : PhAse Recognition and Characterisation (software); PM10/PM2.5: Particulate material with diameter < 10 or 2.5 µm; QXRD: Quantitative X-ray diffraction; SEM: Scanning electron microscopy; SI: Spectral imaging; SFCA: Silico ferrite(s) of calcium and aluminium; TSP: Total suspended particles; Note on cement chemist notation: Several phases are referred to in this article by their nominal chemical formulae expressed using this notation form, where C=CaO, M=MgO, A=Al 2 O 3 , S=SiO 2 , W=FeO and F=Fe 2 O 3 . E.g. C2S=2CaO·SiO 2 =Ca 2 SiO 4 =di-calcium silicate. To avoid confusion with standard chemical formulae, subscripts are not used in the cement chemist notation here.