Characteristics of sub-10 nm particle emissions from in-use commercial aircraft observed at Narita International Airport
- 1Department of Chemistry, Graduate School of Science, Tokyo Metropolitan University, Hachioji, Tokyo 192-0397, Japan
- 2National Institute of Advanced Industrial Science and Technology, Tsukuba, Ibaraki 305-8563, Japan
- 3National Institute for Environmental Studies, Tsukuba, Ibaraki 305-8506, Japan
- 4Environmental Science Analysis and Research Laboratory, Hachimantai, Iwate 028-7302, Japan
Correspondence: Nobuyuki Takegawa (email@example.com)
The characterization of ultrafine particle emissions from jet aircraft equipped with turbofan engines, which are commonly used in civil aviation, is an important issue in the assessment of the impacts of aviation on climate and human health. We conducted field observations of aerosols and carbon dioxide (CO2) near a runway at Narita International Airport, Japan, in February 2018. We used an ultrafine condensation particle counter (UCPC) and a condensation particle counter (CPC) with unheated and 350 ∘C heated operation modes to investigate the contributions of sub-10 nm size ranges to the total and the non-volatile particle number concentrations. The performance of the 350 ∘C heated mode was tested in the laboratory to verify the consistency with existing methods for non-volatile particle measurements. We also used a scanning mobility particle sizer with unheated and 350 ∘C heated modes and an engine exhaust particle sizer for the measurements of particle number size distributions. Spiked increases in the particle number concentrations and CO2 mixing ratios were observed to be associated with the directions of wind from the runway, which can be attributed to diluted aircraft exhaust plumes. We estimated the particle number emission indices (EIs) for discrete take-off plumes using the UCPC, CPC, and CO2 data. The median values of the total and the non-volatile particle number EIs for diameters larger than 2.5 nm as derived from the UCPC data were found to be 1.1×1017 and 5.7×1015 kg per fuel, respectively. More than half the particle number EIs were in the size range smaller than 10 nm for both the total and the non-volatile particles in most of the cases analyzed in this study. The significance of sub-10 nm size ranges for the total particles in the diluted plumes was qualitatively consistent with previous studies, but that for the non-volatile particles was unexpected. Possible factors affecting the similarities and differences compared with the previous findings are discussed.
Civil aviation has grown rapidly as a result of global economic development. Consequently, the impacts of aircraft emissions on climate and human health have been recognized as an important issue (ICAO, 2017; Masiol and Harrison, 2014; Stacey, 2019; and references therein). The characterization of ultrafine particles (UFPs; diameters of <100 nm) is key to understanding the environmental impacts of aircraft emissions because the particle number and mass emissions from jet engines are often dominated by the UFP size range (Masiol and Harrison, 2014; Stacey, 2019; and references therein).
The primary importance of aviation-produced aerosol particles in assessing the climate impacts is the formation of contrail cirrus clouds from soot or black carbon (BC) emitted at aircraft cruising altitudes (Kärcher and Voigt, 2017; Kärcher, 2018). Furthermore, aircraft emissions can significantly affect the number concentrations of Aitken mode particles in the upper troposphere (Wang et al., 2000; Lee et al., 2010; Righi et al., 2013, 2016). Righi et al. (2013, 2016) showed that the impacts of aviation on aerosols were sensitive to the parameterization of nucleation-mode particles ( nm) in their simulation model. In our view, the contribution of aircraft emissions to the number concentration of nucleation or Aitken mode particles relative to that of other sources (new particle formation and surface emission sources) is poorly understood.
The health impacts of UFPs, although they are not specific to aircraft emissions, have been extensively studied by many researchers (Oberdörster et al. (2005), Ohlwein et al. (2019), and references therein). UFPs can be efficiently deposited in the nasal, tracheobronchial, and alveolar regions in the human respiratory system, and the uptake and translocation (physical clearance) of solid UFPs such as soot into the blood and lymph circulation could be an important pathway (Oberdörster et al., 2005). Sub-10 nm particles could be efficiently deposited in the olfactory mucosa, and the subsequent translocation of solid particles along the axons of the olfactory nerve might be a concern (Oberdörster et al., 2005). A recent study suggested that the surface reactivity of aviation-induced soot particles may increase with decreasing particle size (Jonsdottir et al., 2019). However, the health impacts of UFPs emitted from aircraft have not been well established, as pointed out by Ohlwein et al. (2019).
A number of experiments have been performed at engine-test facilities and under real-world conditions to investigate gaseous and particulate emissions from aircraft equipped with turbofan engines, which are commonly used in civil aviation (e.g., Hagen et al., 1998; Petzold et al., 1999, 2005; Kärcher et al., 2000; Brock et al., 2000; Herndon et al., 2008; Westerdahl et al., 2008; Onasch et al., 2009; Kinsey, 2009; Timko et al., 2010, 2013; Lobo et al., 2015a, b; Moore et al., 2017a, b; Kinsey et al., 2019; Yu et al., 2017, 2019; Durdina et al., 2019). The key findings from the previous studies include the following two points. First, significant formation and evolution of volatile particles with diameters smaller than 10 nm can take place during plume expansion, depending on the sulfur content of the fuel, the age of the plume, and the ambient conditions (Petzold et al., 1999, 2005; Kärcher et al., 2000; Brock et al., 2000; Onasch et al., 2009; Timko et al., 2010, 2013). Second, non-volatile particles, which are assumed to be equivalent to soot or BC particles, primarily have sizes in the range larger than ∼10 nm (geometric mean diameter (GMD) of ∼20–60 nm in particle number size distributions) under various operating conditions, with the largest GMDs under ∼100 % engine-thrust conditions (Petzold et al., 1999; Lobo et al., 2015a, b; Moore et al., 2017a; Durdina et al., 2019). Zhang et al. (2019) recently proposed a global-scale aviation emission inventory for BC particles by integrating the existing datasets of non-volatile particle emission indices (EIs).
Alongside these scientific studies, the International Civil Aviation Organization (ICAO) has authorized a new regulatory standard for the mass and number emissions of particles emitted from aircraft engines (ICAO, 2017). In the method for measuring non-volatile particle number concentrations described in the Aerospace Recommended Practice (ARP) 6320, issued by the Society of Automotive Engineers (SAE-ARP6320) (SAE, 2018), the number concentrations of aerosol particles with diameters larger than 10 nm are measured downstream of a volatile particle remover (VPR) heated to 350 ∘C. The SAE-ARP6320 has been developed for the certification of jet engine emissions and may not be directly compared with ambient measurement data. Some of the abovementioned studies employed sampling methods equivalent to SAE-ARP6320 for field measurements of particulate emissions from in-use aircraft (Lobo et al., 2015b).
Direct measurements of UFPs behind jet engines (either in engine test cells or aircraft hangers) can provide systematic emission data as a function of the engine thrust and fuel types from selected jet engines under well-controlled conditions. A common issue in measuring UFPs behind jet engines in the previous scientific studies and in the regulatory standard is the significant loss of particles in long sampling tubes and/or the VPR when fresh aircraft exhaust plumes are sampled (Kinsey, 2009; Lobo et al., 2015a, b; Durdina et al., 2019). Although corrections for particle loss have been extensively evaluated and carefully considered for quantifying particle number concentrations, the absolute values of the correction factors and relative errors associated with the corrections tend to be larger for smaller particles (Kinsey, 2009; Lobo et al., 2015a, b; Durdina et al., 2019). Furthermore, the large uncertainty in measuring the particle number size distributions for diameters smaller than 20 nm by using mobility size spectrometers is also of great concern (Wiedensohler et al., 2012). These technical issues should be properly considered for better characterization of aircraft exhaust particles.
Field measurements of advected (diluted) aircraft exhaust plumes near runways are not optimal for obtaining systematic emission data, whereas potential artifacts associated with long sampling lines and/or high concentrations of condensable materials can be reduced by this approach. Furthermore, a variety of exhaust plumes from different types of in-use aircraft engines can be collectively characterized by field measurements near runways. Considering that the accessibility to platforms for sampling fresh engine exhausts (engine test cells, aircraft hangers, and runways) is generally restricted, these approaches should be complementarily selected for better characterizing UFP emissions from aircraft. Consistent integration of the data is also important for constructing reliable emission inventories from the aviation sectors for the global troposphere (cruising altitudes), where the accessibility to sampling platforms is extremely limited.
We conducted field measurements of UFPs near a runway at Narita International Airport (NRT), Japan. We used multiple instruments for the measurements of particle number concentrations and size distributions and carefully investigated the performance and consistency of these instruments. The purpose of the present study was to investigate the emission characteristics of sub-10 nm particles from commercial aircraft operating under real-world conditions.
2.1 Field observations
The field measurements were performed using two containers placed at an observation point ∼180 m from the centerline of runway A (∼140 m from the edge of the runway) at NRT between 5 and 26 February 2018 (Fushimi et al., 2019). Figure 1 shows an approximate layout of NRT with the location of the observation point. The aerosol instruments used for the field measurements consisted of an ultrafine condensation particle counter (UCPC; model 3776, TSI, Inc., Shoreview, MN, USA; d50=2.5 nm), a condensation particle counter (CPC; model 3771, TSI; d50=10 nm), a scanning mobility particle sizer (SMPS; model 3080, TSI) equipped with a differential mobility analyzer (long DMA; model 3081, TSI) and a CPC (model 3022A; TSI; d50=7 nm), an engine exhaust particle sizer (EEPS; model 3090, TSI), and two sets of cascade impactor samplers (NanoMOUDI-II; model 125B, MSP Corp., Shoreview, MN, USA). The UCPC and CPC used for the field measurements were the same as those used for airborne measurements in our previous studies (Takegawa and Sakurai, 2011; Takegawa et al., 2014, 2017, 2020a), except for the dilution–heater sections described below. The specifications of the UCPC and CPC were similar to those of the TSI CPCs used for aircraft emission measurements near active runways in earlier studies, which included CPC 3007 (d50=10 nm; Westerdahl et al., 2008), CPC 3022A (d50=7 nm; Westerdahl et al., 2008; Herndon et al., 2008; Moore et al., 2017a), and CPC 3775 (d50=4 nm; Moore et al., 2017a). It should be noted that the overall detection efficiencies for particles with diameters of ∼10 nm may significantly depend on the configurations of the sampling tubes rather than a small difference in the d50 values. The penetration efficiencies of particles through the sampling tubes in our system are described in Sect. 3.1 and the Supplement. The combination of CPCs of different detectable size ranges has been used for aircraft emission measurements at cruising altitudes (Brock et al., 2000; Moore et al., 2017b), whereas the combination of a CPC and nano-SMPS/EEPS has been used for measurements near runways (e.g., Westerdahl et al., 2008; Herndon et al., 2008; Moore et al., 2017a). We used the former method to characterize sub-10 nm particles considering its high temporal resolution and relatively large uncertainties in mobility size spectrometers for the sub-20 nm size range (Wiedensohler et al., 2012). Number concentrations of aerosol particles with diameters larger than 2.5 or 10 nm, as measured by the UCPC or CPC, are referred to as N2.5 or N10, respectively. Details of the performance of the UCPC and CPC are described in Sect. 3.1. The other instruments included a carbon dioxide (CO2) monitor (model LI-840, Li-Cor Biosciences, Lincoln, NE, USA), a nitrogen oxide (NOx) monitor (APNA-370, Horiba, Japan), meteorological sensors, and a video camera for monitoring aircraft passages. We used the data from the UCPC, CPC, SMPS, EEPS, and CO2 monitor recorded on 15–22 February 2018 for the present analysis. Although NOx data are useful for characterizing aircraft emissions, we did not use those data because the response time was not sufficient for capturing rapid changes in the concentration in aircraft plumes.
Figure 2a illustrates a schematic diagram of the sampling setup for the UCPC, CPC, SMPS, and CO2 monitor. Ambient air was drawn into the container through a stainless-steel tube (ID: 10.7 mm; length: ∼3 m) and was split into a bypass flow connecting to a piston pump and sample flows for aerosol and CO2 measurements. The total flow rate through the main tube (sum of the bypass and sample flows) was ∼20 L min−1. The aerosol sample flow was diluted by a factor of ∼5 by using particle-free air (2 L min−1) to extend the concentration range measured by the UCPC and CPC, which would otherwise have been limited by particle-coincidence effects. The diluted sample flow was then passed through a stainless-steel evaporation tube (ID: 7.5 mm; tube length: 320 mm; heated section: 200 mm) for heated sampling or through a bypass tube for unheated (room temperature) sampling. A thermocouple sensor was attached to the upstream part of the evaporation tube for temperature control. The three-way valve was switched between the unheated and heated modes every 8 h. We set the heater temperature to 250 ∘C on 7–9 February, 150 ∘C on 11–13 and 22–23 February, and 350 ∘C on 15–21 February. We used the data obtained during the unheated and 350 ∘C heated modes for the present analysis. The characterization of the 150 and 250 ∘C heated modes is ongoing and will be presented elsewhere. A copper tube (ID: 7.5 mm; length: ∼600 mm) between the evaporation tube and the valve was used to cool the heated sample air. The sampling method for N10 heated to 350 ∘C approximately corresponded to the method for measuring non-volatile particles described in SAE-ARP6320 (SAE, 2018), although we did not use a thermal denuder or a catalytic stripper so as to reduce the particle diffusion loss. The estimation of potential artifacts is described in Sect. 2.2. The tube downstream of the three-way valve was split into individual sample flows for the UCPC (1.4 L min−1), CPC (1 L min−1), and SMPS (0.3 L min−1). Electrically conductive tubes (ID: 4.8 mm; Part 3001788, TSI) were used for the connections between the splitter and the instruments. The SMPS was disconnected from the dilution–heater section after 18 February, and the flow rate settings were changed accordingly. Note that the flow rates for calculating the dilution factor were calibrated by using the standard flowmeter at the National Institute of Advanced Industrial Science and Technology (AIST). The dilution factor was estimated to be ∼3.8 before 18 February and ∼6.1 after that date. The residence time in the evaporation tube (320 mm, 350 ∘C) was ∼0.15 s before 18 February and ∼0.17 s after that. These values are shorter than the requirement for the residence time (>0.25 s) described in SAE-ARP6320 (SAE, 2018), although the actual temperature profile of the sample air inside the evaporation tube and downstream of it was not measured. The performance of the evaporation tube was tested using tetracontane (C40H82) particles, as described in Sect. 2.3. After the SMPS had been disconnected, the CPC 3022 was connected directly to the main sampling tube upstream of the dilution–heater section. This permitted an overall validation of the dilution method, because the maximum concentration range of the CPC 3022 was extended to 107 cm−3 by photometric detection. Although we used the dilution method to reduce the effects of particle coincidence, we frequently observed high particle number concentrations (>105 cm−3) downstream of the dilution section. The effects of particle coincidence were corrected using the methods described in our previous studies (Takegawa and Sakurai, 2011; Takegawa et al., 2017). The correction factors were as high as ∼40 % at nominal (uncorrected) particle number concentrations of 5×105 and 1.2×105 cm−3 for the UCPC and CPC, respectively. The uncertainty in the corrections would become larger at even higher concentrations.
As mentioned earlier, the measurements of particles below 20 nm using mobility size spectrometers might include large uncertainties (Wiedensohler et al., 2012). The major sources of uncertainties in the SMPS measurements originated from the corrections for the charging efficiency and Brownian diffusion, the latter generally being more significant for nanoparticles. The Aerosol Instrument Manager (AIM) software provides correction tools for these factors. We have found that the number size distributions at diameters below 20 nm show non-negligible differences for different versions of the software. We used AIM version 9.0 for the present analysis, with the diffusion loss correction enabled. The effects of the AIM diffusion correction on the derived number size distributions are described in Sect. S2.
The EEPS was operated independently from the UCPC–CPC–SMPS inlet system (Fig. 2b), and it measured the unheated particle number size distributions during the entire period. The sample flow rate of the EEPS was 10 L min−1. A copper tube (i.d.: 10 mm, length: ∼2 m), electrically conductive tubes (i.d.: 4.8 and 7.9 mm; total length: ∼1 m; part 3001788, TSI), and a glass manifold (inner diameter: 40 mm; total length: 600 mm) were used for the ambient sampling. A bypass pump (flow rate: 10 L min−1) was used to reduce the particle diffusion loss. We used the default instrument matrix for the EEPS, which may underestimate size and concentration of particles larger than ∼75 nm (Wang et al., 2016). Although the particle diameter range detectable by the EEPS extended from ∼6 to 520 nm, our laboratory experiments showed that the EEPS may significantly underestimate particle number concentrations below 10 nm. Reduced detection efficiencies of sub-10 nm particles can also be inferred from the EEPS data obtained by other researchers (Moore et al., 2017a). The evaluation of the EEPS is ongoing and will be presented elsewhere. In this study we used the EEPS data for particle diameters larger than 10 nm.
The CO2 instrument was calibrated by using two CO2 standards (397.2 and 1032 ppmv – parts per million by volume) twice a day during the measurement period. The injection of the CO2 standards was performed automatically by using solenoid valves. We found that the sensitivity of the instrument was generally stable during the measurement period.
The N2.5, N10, EEPS, and CO2 data were obtained every 1 s, and the SMPS data for particle diameters from 15 to 660 nm were obtained every 5 min (scanning time: 3 min). When we observed spiked increases in the particle number concentrations and CO2 mixing ratios in aircraft plumes, the timing of the detection of the concentration peaks did not exactly match between the individual instruments. Because this was likely caused by differences in the response time of the instruments and the delay time in the sampling tubes, the data were shifted accordingly (<10 s).
The overall penetration efficiencies of particles through the sampling tubes were estimated by using the theoretical formulae proposed by Gormley and Kennedy (1949). Details of the calculation procedures are given in Sect. S1 of the Supplement. The penetration efficiency for the UCPC and CPC sampling line (unheated mode; from the top of the inlet to the flow splitter) was estimated to be 70 %, 87 %, and 94 % for particle diameters of 5, 10, and 20 nm, respectively. The penetration efficiency through the evaporation tube and the detection efficiencies of the UCPC and CPC are important parameters for the data interpretation and are discussed in detail in Sect. 3.1. The penetration efficiency for the SMPS sampling line (unheated mode; from the top of the inlet to the SMPS) was estimated to be 79 % and 85 % for 15 and 20 nm, respectively, and that for the EEPS sampling line (from the top of the inlet to the EEPS) was estimated to be 94 % and 98 % for 10 and 20 nm, respectively.
The particle diffusion loss during sampling is an important issue for the quantification of UFPs, as mentioned in Sect. 1. The corrections for the penetration efficiencies through the sampling tubes and the detection efficiencies (see Sect. 3.1 for details) were not incorporated in the UCPC and CPC data presented in Sect. 3.2.1–3.2.4 because the actual size distributions in the sub-10 nm size range were uncertain. Furthermore, the corrections for the penetration efficiencies were not incorporated in the SMPS and EEPS data presented in Sect. 3.2.2 and 3.2.3 for consistency with the UCPC–CPC data: i.e., they include only the default (internal) corrections for the individual instruments. We considered the effects of particle diffusion loss as systematic uncertainties. We also considered the penetration efficiencies through the sampling tubes and the detection efficiencies for estimating the size distributions of particle number EIs in Sects. 3.2.5 and 4.
2.2 Potential artifacts
Potential artifacts due to the nucleation of gaseous compounds vaporized from particles in the evaporation tube (hereafter referred to as nucleation artifacts) were evaluated. Predicting the nucleation rates requires an estimate of the supersaturation of nucleating compounds, which is highly uncertain. Here we evaluate the growth rate of nucleated clusters under the given condition. The upper limit of this effect can be estimated by assuming that the number concentration of the vaporized compounds remains constant after a certain period of time (∼1 s) and that all the compounds are condensed into nucleated clusters (e.g., sulfuric acid). In real situations, only a small fraction of the vaporized compounds might contribute to the condensational growth of particles because of possible increases in their saturation vapor pressures by thermal decompositions and their deposition onto the inner surface of the sampling tube.
Let us assume a condensable material (number concentration of molecules: c; molecular weight: MW) in the gas-phase vaporized from particles. The original mass concentration in the particle phase before heating, m, is expressed as follows:
where NA is the Avogadro number. Assuming a kinetic regime, the growth of nucleated particles is governed by the following expression (Seinfeld and Pandis, 2006):
where ρ is the particle density, Dp is the diameter of nucleated particles, is the mean thermal velocity of condensable gas molecules, α is the accommodation coefficient, and ceq is the equilibrium number concentration. Assuming α=1 and ceq=0 (maximum molecular flux), we obtain the expression
If we consider sulfuric acid as a condensable material, dDp∕dt (m s−1) can be approximated as 0.1 m at 350 ∘C (maximum gas temperature), where m is expressed in units of kilogram per cubic meter. The residence time for particle growth after the evaporation tube is estimated to be <1 s. The residence time of 1 s and the dilution factor of 5 lead to the estimation that the maximum particle growth is 1 nm at an ambient mass concentration of 50 µg m−3. The minimum particle diameter at which the detection efficiency becomes zero is ∼2 and ∼7 nm for the UCPC and CPC, respectively (Takegawa and Sakurai, 2011; Takegawa et al., 2017). If we assume an initial cluster size of 1 nm (approximately the critical cluster size of sulfuric acid), the effects of the artifacts on N2.5 would be significant in aircraft exhaust plumes with relatively high concentrations (>50 µg m−3), and those on N10 (and SMPS) would be significant at very high concentrations (>300 µg m−3). This estimation does not change significantly if we assume high-molecular-weight organic compounds from jet engine lubrication oil (slower thermal velocities and smaller particle densities compared with sulfuric acid).
Other potential artifacts may originate from the condensational growth of non-volatile particles (or residual particles downstream of the evaporation tube) smaller than the detectable size range of the UCPC (diameter <2 nm). For example, residual particles with diameters of ∼2 nm may grow to ∼3 nm at an ambient mass concentration of 50 µg m−3. However, this effect is significant only in the presence of a large fraction of non-volatile particles with diameters below 2 nm.
2.3 Laboratory experiments
The accuracy of the measurements of particle number concentrations was the key issue in this study and thus was evaluated in the laboratory at AIST before and after the field measurements. We mainly used the data obtained after the field measurements because they were more comprehensive than those obtained before the measurements. The test items included the size-resolved detection efficiencies of the UCPC and CPC, the penetration efficiency of non-volatile particles through the dilution–heater section, and the removal efficiency of volatile compounds through the evaporation tube.
An electrospray aerosol generator (EAG; model 3480, TSI), a combustion aerosol standard (CAST; Matter Engineering, AG, Wohlen, Switzerland) with a tube furnace for thermal treatment at 350 ∘C, and a custom-made tube furnace for supplying condensable vapors were used to generate polydisperse aerosol particles for the calibrations. We used sucrose particles supplied from the EAG (particle diameter range: 6–15 nm) for the detection efficiency experiment, non-volatile propane soot particles supplied from the CAST (15–100 nm) for the penetration efficiency experiment, and tetracontane particles supplied from the tube furnace particle generator (15–30 nm) for the removal efficiency experiment. The tube furnace downstream of the CAST was used to remove organic compounds internally or externally mixed with soot. We used a volatility tandem DMA with a heater temperature of 350 ∘C to confirm that the thermal treatment efficiently removed organic compounds; i.e., there was no significant difference in the mobility diameter of soot particles between the two DMAs.
The experimental apparatus for measuring the size-resolved detection efficiencies of the UCPC and CPC was similar to that used in our previous studies (Takegawa and Sakurai, 2011; Takegawa et al., 2017), except that the sampling line was longer for the CPC. The length of the sampling line for the CPC was ∼90 cm, which was limited by the instrument layout in the observation rack. We set the particle diameters at 15, 30, 50, and 100 nm for testing the penetration efficiencies and at 15 and 30 nm for testing the removal efficiencies, in accordance with the SAE-ARP6320 protocol. The flow rate through the evaporation tube was set at 2.4 L min−1. We used an aerosol electrometer (AE; model 3068B, TSI) and a reference CPC (ref CPC; model 3775, TSI) as reference instruments. These were calibrated with the standard AE maintained by AIST for particle number concentrations at ambient pressures in the size range of 10–300 nm.
We also tested the removal efficiencies of tetracontane particles for diameters of 30 and 50 nm in the laboratory at Tokyo Metropolitan University (TMU). We used another CPC (model 3772; TSI), which is essentially the same as CPC model 3771, for the detection of residual particles downstream of the evaporation tube because the conditions for the CPC 3771 were not optimized during the experiments at TMU. The flow rate through the evaporation tube was set at either 2.4 or 2.7 L min−1 to test the dependency of the removal efficiencies on the flow rate. Considering that the polydisperse size distributions of the tetracontane particles generated at TMU were rather broad, we also measured the removal efficiencies of doubly charged particles (43 for 30 nm and 72 for 50 nm).
3.1 Laboratory experiments
Figure 3a shows the size-resolved detection efficiencies of the UCPC and CPC measured at AIST. The detection efficiencies for the UCPC and CPC were empirically estimated by using our previous calibration results for the CPC (Takegawa and Sakurai, 2011), the manufacturer specifications for the UCPC (Takegawa et al., 2017), and the penetration efficiencies in the instrument (for the UCPC; Wimmer et al., 2013) and the sampling lines. The penetration efficiencies were calculated by using the theoretical formulae proposed by Gormley and Kennedy (1949). Further details of the empirically estimated detection efficiencies are given in Sect. S1 in the Supplement. We found a good agreement between the experimental data and the estimated detection efficiencies for both the UCPC and CPC.
Figure 3b shows the penetration and detection efficiencies of non-volatile soot particles through the 350 ∘C heated mode sampling tubes as measured by the UCPC. The temperature was set at room temperature (∼19–21 ∘C) and at 350 ∘C for comparison. The penetration and detection efficiencies include the penetration efficiency through the dilution–heater section and the detection efficiency of the UCPC shown in Fig. 3a. Based on the specifications in the SAE-ARP6320 protocol, the penetration efficiency through a VPR should be equal to or higher than 30 %, 55 %, 65 %, or 70 % for particle diameters of 15, 30, 50, and 100 nm, respectively. Our data showed that the penetration efficiency was well above the requirements for all diameters tested. It was rather difficult to perform a simple theoretical prediction of the penetration efficiency at 350 ∘C because the actual temperature profile of the sample air was uncertain. Therefore, the penetration efficiency curve at 350 ∘C (not shown in Fig. 3) was determined by scaling the calculated curve at 20 ∘C with the experimental values at 350 ∘C and at room temperature for larger diameters (>30, 50, and 100 nm).
The SAE-ARP6320 protocol specifies that the removal efficiencies of tetracontane particles in a VPR should be higher than 99.9 % for particle diameters of 15 and 30 nm. For the removal efficiency test at AIST, we confirmed that the removal efficiencies of tetracontane particles for diameters of 15 and 30 nm were >99.9 % (remaining fraction of <0.1 %) for both the UCPC and CPC. For the removal efficiency test at TMU, we first confirmed that there was no significant difference exceeding the experimental uncertainties between the two flow rate settings (2.4 and 2.7 L min−1) for both the UCPC and CPC, and we mainly used the results for 2.7 L min−1 (shorter residence time). Table 1 presents the experimental results. The remaining fractions for 30 nm particles were found to be 0.3 () % and <0.1 % for the UCPC and CPC, respectively. The difference between the results from the AIST and TMU experiments for the UCPC was not well identified, and we took the results from the TMU experiments for conservative estimates. The remaining fractions for 50 nm particles were found to be 5 () % and <0.1 % for the UCPC and CPC, respectively. The remaining fraction for 73 nm was not quantified because of low particle number concentrations (thus the multiple-charge correction for 50 nm was minor). These results suggested that about 5 % of the 50 nm tetracontane particles might not have fully vaporized but shrunk to sizes between 2.5 and 10 nm downstream of the evaporation tube. Potential influences on the interpretation of the ambient data are discussed later.
3.2 Field observations
3.2.1 Particle number concentrations
Figure 4 shows time series of N2.5 and CO2 obtained on 15 February 2018. The inlet valve was switched from the unheated to the 350 ∘C heated mode at 16:00 LT (local time). The wind directions were east or east-southeast for Fig. 4a (average ± SD (1σ) of the wind speeds was 3.9±0.7 m s−1) and east or east-northeast for Fig. 4b (average ±1σ of the wind speeds was 4.3±0.6 m s−1), indicating that relatively stable winds from the direction of the runway were dominant. The ambient temperature was ∼10 ∘C for Fig. 4a and ∼8 ∘C for Fig. 4b. We did not observe a significant enhancement of aerosols or CO2 when the air parcels originated from the opposite direction from the runway during aircraft operating times (06:00–23:00 LT). The spikes in the concentrations of aerosol particles and in CO2 mixing ratios shown in Fig. 4 can be interpreted as resulting from diluted exhaust plumes from individual aircraft movements. The observed values of N2.5 (and N10) for the 350 ∘C heated mode (Fig. 4b) were significantly lower than those for the unheated mode (Fig. 4a) when comparing plumes with similar enhancement levels of CO2. The depletions of N2.5 (and N10) for the 350 ∘C heated mode were much larger than those expected from the difference in the penetration efficiencies of aerosol particles between the unheated and the 350 ∘C heated modes (see Fig. 3). Although the characteristics of the plumes shown in Fig. 4a and b cannot be directly compared, the similarity in the wind conditions suggests that the ages of the plumes were comparable (∼50 s). These results imply that most of the aerosol particles in the observed plumes were volatile. We have previously shown that jet engine lubrication oil was the major source of aerosol particles with diameters ranging from ∼10 to 30 nm in air parcels observed during the measurement period (Fushimi et al., 2019). A key point in Fig. 4 is that the number fraction of particles with diameters ranging from 2.5 to 10 nm (red area) was generally larger than that of particles with sizes above 10 nm (gray area), for both the unheated and 350 ∘C heated modes. The significance of sub-10 nm size ranges for the total particles in diluted plumes (10–102 m from the jet engines) was qualitatively consistent with the findings from previous studies conducted on the ground (e.g., Petzold et al., 2005; Lobo et al., 2012, 2015b) or in-flight at cruise altitudes (Petzold et al., 1999; Kärcher et al., 2000; Brock et al., 2000), but that for non-volatile particles was unexpected.
Figure 5a is a scatterplot of 1 s averaged values of N10 vs. those of N2.5 for the 350 ∘C heated mode for the time period from 17:00 to 23:00 LT on 15 February 2018, the same data as used in Fig. 4b. Data obtained on 21 February 2018 at slightly different flow rate settings (see Sect. 2.1) are also plotted to show the reproducibility of the results. The 1:1 correspondence line (N10=N2.5) represents the state in which all particles are included in the size range larger than 10 nm, and deviations below the 1:1 line indicate an increase in the fractions of sub-10 nm particles. These results show that many data points lie below the 1:2 line (sub-10 nm fraction of >50 %) when N2.5 exceeds 105 cm−3. This point is more quantitatively found in the median values for each concentration bin (Fig. 5b). Some data points in Fig. 5a were found on the 1:10 correspondence line (sub-10 nm fraction of 90 %). The deviations from the 1:1 line would have been even larger had we considered the size-dependent loss of particles in the sampling line shown in Fig. 3.
3.2.2 Size distributions of total particles
Figure 6 shows time series of total particle number concentrations measured by the CPC and EEPS and the particle number and volume size distributions (dN∕dlog D and dV∕dlog D) measured simultaneously by the EEPS and unheated SMPS for selected time periods on 15 February 2018. The total particle number concentrations measured by the CPC (d50=10 nm) and EEPS (10–520 nm) showed reasonable agreement, which confirms the overall consistency between these two independent measurements. The time periods were carefully selected for investigating the effects of rapid changes in the particle number concentration during each SMPS scan on the derived particle size distribution. The SMPS needed ∼30 s to scan the major population of particle number concentrations in the aircraft plumes (<30 nm). Considering possible delays in the detection timing due to the residence time in the SMPS sampling tubes (<5 s), the EEPS data were averaged over 40 s around the timing of the SMPS scan corresponding to the diameter range of 15–30 nm (i.e., averaged from 10 s before to 30 s after the onset of each SMPS scan). For the SMPS scan at 14:10 LT, the particle number concentration varied by a factor of ∼2 during the first 30 s scan time, but there was no systematic increase or decrease in the particle number concentration. We found a reasonable agreement between the EEPS and SMPS (within a factor of ∼2). For the SMPS scan at 14:20 LT, the particle number concentration decreased significantly during the 30 s scan time and then started to increase rapidly afterwards. The resultant number size distribution was likely affected by these concentration changes. These results demonstrate that the SMPS data can be used to investigate the characteristics of particle size distributions within specific size ranges by selecting the appropriate time windows. Note that the SMPS data tend to exhibit higher concentrations (by a factor of ∼2) at smaller diameters (<20 nm) compared with the EEPS data in both cases. The difference cannot be explained by the uncertainty due to the penetration efficiencies through the sampling tube (as indicated by error bars in Fig. 6b). The difference may be due to uncertainties in the default (internal) correction algorithms in the SMPS (see Fig. S2). Nevertheless, the overall size dependency (i.e., increasing particle number concentrations with decreasing particle diameters below 20 nm) is consistent between the SMPS and EEPS data.
In the enhancement events at 14:10 and 14:20 LT in Fig. 6, the EEPS data indicate that the particle number size distribution functions below 50 nm showed significant increases compared with those in the non-enhancement event (14:00 LT). The particle volume size distribution function below 50 nm exhibited moderate increases in the enhancement events, but the total integrated particle volume concentrations were largely affected by accumulation-mode particles (>50 nm) in the background air. A similar feature was also found for other time periods (Fushimi et al., 2019).
3.2.3 Size distributions of non-volatile particles
Figure 7 shows the particle number and volume size distributions measured simultaneously by the EEPS and the 350 ∘C heated SMPS for selected time periods with and without enhancements of aerosol particle number concentrations on 15 February 2018. The time periods were selected so that there would be no systematic increase or decrease in the particle number concentrations during the first 30 s scan time of the SMPS. Similarly to Fig. 6, the particle number size distribution functions below ∼50 nm in the enhancement events showed significant increases compared to those in the non-enhancement event for both the EEPS and the 350 ∘C heated SMPS. The particle number size distribution functions below ∼50 nm measured by the 350 ∘C heated SMPS were smaller by more than an order of magnitude compared to those measured by the EEPS, indicating that the aircraft exhaust particles were mostly volatile below 50 nm. Our laboratory experiments suggest that there might remain residues of >50 nm particles in the sub-10 nm size range after the evaporation tube (the remaining fraction was ∼5 % for 50 nm tetracontane particles). However, they were likely negligible compared with the observed sub-10 nm non-volatile particles because the particle number concentrations for a diameter range of >50 nm measured by the EEPS (unheated) were far below the observed values of the 350 ∘C heated N2.5–N10 (Fig. 7a). In addition, residues of 30–50 nm particles in the sub-10 nm size range after the evaporation tube were likely minor compared with the observed sub-10 nm non-volatile particles considering the sharp decrease in the dN∕dlog D values from 30 to 50 nm as measured by the EEPS (the remaining fractions were ∼0.3 % and 2 % for 30 and 43 nm tetracontane particles, respectively).
The effects of nucleation artifacts (Sect. 2.2) might be a major concern but were likely small under the observation conditions because the mass concentrations of aerosol particles in the aircraft plumes inferred from Fig. 7 were much lower than the threshold concentration (50 µg m−3) described in Sect. 2.2. In fact, we did not find any systematic increases in the N2.5∕N10 ratios for the 350 ∘C heated mode with increasing total particle volume concentrations derived from the EEPS (not shown). Furthermore, the particle number and volume size distributions derived from the 350 ∘C heated SMPS (Fig. 7b and c) exhibited no indication of the presence of an additional mode resulting from nucleation artifacts (Sect. 2.2). The 350 ∘C heated dN∕dlog D functions in Fig. 7b exhibited gradual increasing trends with decreasing particle diameters from 100 to 15 nm. This feature is unlikely to be explained by the artificial growth of particles downstream of the evaporation tube because the mass concentrations of aerosols in the plumes would not have been sufficient to yield particle growth exceeding 1 nm (Sect. 2.2).
3.2.4 Particle emission indices for take-off plumes
The temporal variations in aerosol particles and number size distributions of aerosol particles clearly indicate that the observed air parcels were significantly affected by aircraft emissions under appropriate wind conditions. However, aircraft emissions from various cycles of take-off, landing, and idling may have been mixed in the atmosphere. We calculated the enhancements of N2.5, N10, and CO2 above the background levels (referred to as ΔN2.5, ΔN10, and ΔCO2, respectively) and extracted discrete plumes during the take-off phases by setting some criteria for these parameters (for details, see Sect. S3 and Figs. S3–S5). By using the data obtained on 15, 16, 20, 21, and 22 February, we identified 132 discrete plumes for the unheated mode and 63 for the 350 ∘C heated mode. Particle number and volume size distributions (dN∕dlog D and dV∕dlog D) for the take-off plumes were derived from the EEPS data (unheated conditions), and the enhancements of dN∕dlog D and dV∕dlog D above the background levels were calculated using the same method as for ΔN2.5.
Note: The particle number EI values for unheated N2.5, unheated N10, 350 ∘C heated N2.5, and 350 ∘C heated N10 are referred to as total EI(N2.5), total EI(N10), non-volatile EI(N2.5), and non-volatile EI(N10), respectively.
The ΔN2.5∕ΔCO2 and ΔN10∕ΔCO2 ratios were converted into particle number EIs by assuming a CO2 EI value of 3160 g of CO2 per kilogram of fuel (Moore et al., 2017a), as shown in Table 2. The corrections for the penetration and detection efficiencies of the UCPC and CPC were not incorporated into the above estimates because the actual size distributions were uncertain; these values should therefore be regarded as lower limits. This point is further discussed in Sect. 4.1. The arrival time of the plumes was estimated to be ∼30–120 s by considering the wind directions and speeds, corresponding to the transport distance of ∼180–370 m from the engine exits to the observation point (Sect. S3). We did not find a systematic dependence of the ΔN10∕ΔCO2 and ΔN2.5∕ΔCO2 ratios on the estimated arrival time of the plumes (Fig. S5), suggesting that the variation in plume ages did not yield significant biases in our results. We did not perform a detailed classification of the plumes by jet engine types because information on the engines of the aircraft observed at NRT was not available. We performed a simple classification of the plumes by aircraft model (Table S2), but we did not find significant differences between the different types of aircraft.
The median total and non-volatile EI(N2.5) values, which likely cover the major size range of the aircraft emissions, were found to be 1.1×1017 and 5.7×1015 kg per fuel, respectively. The difference in these values (a factor of ∼20) is interpreted as the average contribution of volatile particles. We defined the sub-10 nm fraction as for the identified discrete plumes. The median and the central 50th percentile range of the sub-10 nm fraction for the unheated mode were found to be 0.63 and 0.53–0.70, respectively, and those for the 350 ∘C heated mode were 0.54 and 0.44–0.72, respectively (Table 2). The significance of sub-10 nm particles for the total and the non-volatile particles, which was shown in Sect. 3.2.1–3.2.3 (Figs. 4–7) as a case study, was also found in the statistical analysis of the take-off plumes. Considering that the penetration efficiencies of particles through the sampling tubes (Fig. S1) and the evaporation tube (Fig. 3) tended to decrease with decreasing particle diameters below 10 nm, these results suggest that more than half the total and the non-volatile particle number EIs in the aircraft take-off plumes were found in the size range smaller than 10 nm in most cases.
Previous studies have shown that the particle number EIs can vary significantly depending on the engine type, the engine thrust, the fuel sulfur content, the plume age, and the ambient conditions (Petzold et al., 1999, 2005; Kärcher et al., 2000; Brock et al., 2000; Onasch et al., 2009; Timko et al., 2010). The sampling setup (deployed instruments and sampling location relative to the runways) and the analysis procedures (discrete plume analysis) of this study are similar to those of Lobo et al. (2012, 2015b) and Moore et al. (2017a). Therefore, it is worth discussing the similarities and differences between those studies and our results. The fuel sulfur content is an important parameter for the comparison with other studies. We do not have information on the sulfur content of the fuel that was actually used at NRT. Instead, we analyzed fuel samples (Jet A-1) provided by a jet fuel company in Tokyo (Ishinokoyu, Co. Ltd.) (Saitoh et al., 2019b). We obtained a total of five samples between August 2017 and August 2018. The sulfur content of the fuel samples ranged from 30.4 to 440 ppmw (parts per million by weight). We assume that these values are representative of the sulfur content of jet fuels commercially available in Tokyo during the observation period.
Lobo et al. (2012) reported particle number and mass EIs measured 100–350 m downwind of the runways at Oakland International Airport in August 2005. They used a fast particulate spectrometer (DMS500, Cambustion) to measure the particle number size distributions for diameters ranging from 5 to 1000 nm. The total particle number EIs for various types of engines under take-off conditions ranged from 4×1015 to 2×1017 kg per fuel, which inclusively covered the 25–75th percentile range of the total EI(N2.5) from our measurements. The fuel sulfur content was estimated to be 240–395 ppmw.
Lobo et al. (2015b) reported the particle number and mass EIs measured near the jet engine exits and 100–350 m downwind of the runways at Hartsfield–Jackson Atlanta International Airport in September 2004. The former corresponds to the non-volatile particle EIs and the latter corresponds to the total particle EIs. They used a DMS500 to measure the particle number size distributions. The non-volatile particle number EIs for a JT8D-219 engine were kg per fuel under 85 % and 100 % thrust conditions, which were close to the 75th percentile value of the non-volatile EI(N2.5) from our measurements. The total particle number EIs for various types of engines under take-off conditions ranged from 7×1015 to 9×1017 kg per fuel, which again inclusively covered the 25–75th percentile range of the total EI(N2.5) from our measurements. The information on the fuel sulfur content was not provided.
Moore et al. (2017a) reported the particle number and volume EIs for take-off plumes measured 400 m downwind of the runway at Los Angeles International Airport (LAX) in May 2014. They used a CPC 3775 (d50=4 nm) for the measurements of total particles, a CPC 3022A (d50=7 nm) for non-volatile particles, and an EEPS for the measurements of size distributions. We calculated the median values (and the 25–75th percentile range) of the total and the non-volatile particle number EIs provided by Moore et al. (2017a) to be 4.6 (3.1–5.8)×1016 and 2.1 (1.1–3.6)×1015 kg per fuel, respectively. The median and 25–75th percentile range of the total and the non-volatile EI(N10) from our measurements showed good agreement with those from Moore et al. (2017a). Although Moore et al. (2017a) used different types of CPCs, the detection efficiency curves for the CPC 3771 and 3022A were likely similar according to the manufacturer specifications. Furthermore, the EEPS data from Moore et al. (2017a) showed that the contributions from sub-10 nm particles to the total particle number EIs were relatively small (see Sect. 3.2.5 for details). Therefore, the good agreement between their results and the EI(N10) values from our measurements would be reasonable. The sulfur content of the jet fuel samples collected at LAX ranged from 620 to 1780 ppmw (average value of 1180 ppm).
Timko (2010) showed that the total particle number EIs in moderately diluted plumes (measured by a CPC 3022A), which were dominated by volatile particles, exhibited a relatively weak dependence on the fuel sulfur content (< 1500 ppmw) under high-thrust conditions for various types of engines. Provided that the fuel sulfur content was likely below ∼1500 ppmw for the previous studies (Lobo et al., 2012, 2015b; Moore et al., 2017a) and this study, it might not be the major factor affecting the variability in the emissions of volatile particles among these studies (at least for diameters larger than ∼10 nm).
3.2.5 Size distributions of total particles for take-off plumes
The enhancements of dN∕dlog D and dV∕dlog D above the background levels derived from the EEPS data were converted to the corresponding EI values (dEI(N)∕dlog D and dEI(V)∕dlog D). As described in Sect. 2.1, the corrections for the penetration efficiencies through the sampling tubes (Fig. S1) were considered for the data conversion. Figure 8 shows the derived dEI(N)∕dlog D and dEI(V)∕dlog D values for the take-off plumes. The size distributions of the particle number and volume EIs for the take-off plumes observed at LAX (Moore et al., 2017a), which are characterized by bimodal lognormal distributions (nucleation and soot modes), are shown for comparison. Other previous studies also reported that aircraft emissions can be characterized by distinct bimodal size distributions under high-engine-thrust conditions (e.g., Kinsey, 2009; Lobo et al., 2012, 2015b; Yu et al., 2017, 2019).
The total and non-volatile particle number EIs derived from the UCPC and CPC fell in the same range as those from the previous studies for take-off plumes under real-world operating conditions (Lobo et al., 2012, 2015b; Moore et al., 2017a), as described in Sect. 3.2.4. However, the characteristics of the size distributions appeared to be significantly different. The mode diameters of the dEI(N)∕dlog D for the nucleation mode reported by the previous studies were 10–20 nm in most cases. We compared the dEI(N)∕dlog D values between each size bin (midpoint diameters of 10.8, 12.4, 14.3, … nm) for the individual take-off plumes shown in Fig. 8. We found that the dEI(N)∕dlog D at 10.8 nm exhibited the highest value for ∼98 % of the plumes. We also found that the dEI(N)∕dlog D at 10.8 nm was more than 2 times larger than that at 14.3 nm for 79 % of the plumes. These results suggest that the dEI(N)∕dlog D values tended to increase with decreasing particle diameters around 10–20 nm and that the mode diameters of the dEI(N)∕dlog D for the nucleation mode were smaller than 10 nm for the majority of the plumes observed in this study. Our results also suggest that the peak values of the dEI(N)∕dlog D for the nucleation mode were much larger than those reported by Moore et al. (2017a), in which the contributions from the sub-10 nm size range to the total particle number EI size distributions were relatively small.
The uncertainties in the dEI(N)∕dlog D and dEI(V)∕dlog D values at larger particle diameters (>100 nm) were considerably large in our data because they were significantly affected by the accumulation-mode particles in the background air and also because the default instrument matrix of the EEPS might underestimate size and concentrations of particles larger than ∼75 nm (Wang et al., 2016). Nevertheless, it is likely that the non-volatile particle number and volume EIs originating from soot-mode particles (>20 nm) were much smaller than those reported by the previous studies (Lobo et al., 2012, 2015b; Moore et al., 2017a). A correction for the EEPS size bin by a scaling factor of 1.1, as was done by Moore et al. (2017a), would increase the particle volume concentrations by a factor of 1.33, which does not fill the gap between the previous studies and our results.
4.1 Estimates of particle number EI size distributions for take-off plumes
We estimated the possible particle number size distributions for the total and the non-volatile particles constrained by the UCPC and CPC observations and investigated the consistency with the EEPS and SMPS measurements. We assumed lognormal number size distributions with various geometric mean diameters (GMDs) and geometric standard deviations (GSDs). We calculated the number fraction of sub-10 nm particles () by integrating the number size distributions weighted by the penetration and detection efficiency curves:
where ηCPC and ηUCPC are the overall penetration and detection efficiencies for the CPC and UCPC, respectively. We assumed monomodal size distributions for dN∕dlog D and dEI(N)∕dlog D considering the shapes of the observed size distributions for the total and the non-volatile particles (Figs. 6 and 7). The overall penetration efficiency, which included the sampling tubes from the rooftop and the dilution–heater sections (Figs. 2 and 3b), was calculated based on the theoretical formulae by Gormley and Kennedy (1949). The GMD and GSD values assumed in Eq. (4) were also used to calculate the dEI(N)∕dlog D in Eq. (5). The calculated EI(N2.5) was compared with the median EI(N2.5) derived from the observations to retrieve the absolute values of dEI(N)∕dlog D and EI(N).
Figure 9 shows the calculation results for the total particles. GMD values of <10 nm are needed to explain the median sub-10 nm fraction derived from the plume analysis (0.63) with a realistic range of GSDs (1.2–1.6). For example, the GMD value of 7.9 nm is obtained if we assume a GSD value of 1.6 (or the GMD value of 8.6 nm is obtained if we assume a GSD value of 1.4). The retrieved dEI(N)∕dlog D values were found to be consistent with those derived from the EEPS (Fig. 8), which supports the validity of our estimate. The retrieved EI(N) value was 1.7×1017 and 1.6×1017 kg per fuel for (GMD, GSD) = (7.9 nm, 1.6) and (8.6 nm, 1.4), respectively. These values were larger by a factor of ∼1.5 compared to the median total EI(N2.5) value, suggesting that the particle number EIs derived from the unheated UCPC might underestimate the true particle number EIs by a factor of
Figure 10 shows the calculation results for the non-volatile particles. Similarly to the total particles, GMD values of <10 nm would also be needed to explain the median sub-10 nm fraction derived from the plume analysis (0.54). The GMD value of 9.0 nm is obtained if we assume a GSD value of 1.6 (or the GMD value of 9.5 nm is obtained if we assume a GSD value of 1.4). The retrieved dEI(N)∕dlog D values were found to be consistent with those derived from the SMPS, although the size distribution data to support the validity of this estimate are limited. The retrieved EI(N) value was 1.0×1016 and 0.95×1016 kg per fuel for (GMD, GSD) = (9.0 nm, 1.6) and (9.5 nm, 1.4), respectively. These values were larger by a factor of ∼1.7 compared to the median non-volatile EI(N2.5) value, suggesting that the particle number EI derived from the 350 ∘C heated UCPC might underestimate the true non-volatile particle number EIs by a factor of
4.2 Interpretation and implications
The characteristics of the total particles observed in this study are qualitatively consistent with the findings from previous studies (e.g., Petzold et al., 1999, 2005; Lobo et al., 2012, 2015b; Kärcher et al., 2000; Brock et al., 2000): i.e., the total particle number EIs are dominated by volatile particles, and the sub-10 nm particles make significant contributions to the total particle number EIs in the diluted plumes (10–102 m from the jet engines). In contrast to the total particles, the significance of sub-10 nm size ranges for the non-volatile particle number EIs is unexpected. Currently, we do not have direct evidence that the observed sub-10 nm non-volatile particles are composed mainly of soot. Nevertheless, it is worth considering possible mechanisms for the formation of sub-10 nm non-volatile particles. Our estimate implies the presence of very small non-volatile particles with diameters down to a few nanometers. This is not consistent with the size of the primary soot particles from jet engines estimated by using transmission electron microscopy (TEM) or laser-induced incandescence (LII) methods (e.g., Liati et al., 2014; Boies et al., 2015; Saffaripour et al., 2017, 2020), although there still remain substantial uncertainties in estimating the size of primary particles by those methods, as pointed out by Boies et al. (2015).
An alternative possibility for the significance of sub-10 nm non-volatile particles includes the potential contributions of less volatile organic matter (compared with tetracontane) or metal compounds. Fushimi et al. (2019) showed that the mass contribution of the sum of trace elements (other than carbonaceous and sulfur compounds) was comparable to that of elemental carbon (soot) for UFP samples (∼10–30 nm) collected at NRT. Saitoh et al. (2019a) reported that metal elements including Ca, Fe, Si, Mg, K, Zn, Pb, and Ni were the major compositions of these trace elements. These metal compounds might have contributed to the non-volatile particles observed in this study.
The key point in our results is that the non-volatile particle number and volume EIs originating from soot-mode particles (>20 nm) were much smaller than those reported by the previous studies for take-off plumes under real-world operating conditions (Lobo et al., 2012, 2015b; Moore et al., 2017a). This feature might be related to the significance of sub-10 nm non-volatile particles. The data presented by Moore et al. (2017a), which were obtained in 2014, exhibited lower contributions of soot-mode particles compared to those by Lobo et al. (2012, 2015b), which were obtained in 2005 and 2004, respectively. A possible explanation for this tendency is that newer engines might emit smaller non-volatile particles compared with older engines.
Our results have an implication for jet engine exhaust measurements by the SAE-ARP6320 method. The non-volatile N10 data obtained from the 350 ∘C heated CPC approximately correspond to the definition of non-volatile particles specified by SAE-ARP6320, at least in terms of the detectable size range of the particle counter and the removal efficiency of tetracontane particles (Sect. 3.1). If a population of non-volatile particles having a similar size distribution as shown in Fig. 10 is measured by the SAE-ARP6320 method, subtle changes in the size-dependent penetration and detection efficiencies of the measurement system and/or shifts in the mode diameter of the particle population might lead to large variabilities in the results. Further investigations are needed to quantify the size distributions of non-volatile particles for various types of engines under different conditions.
Our results also have an implication for the emission inventories of the aviation sector. Although the potential contributions of sub-10 nm particles inferred from our results likely have negligible impacts on the mass concentrations of ambient aerosol particles, they may have non-negligible impacts on the number concentrations in and around airports and also at aircraft cruising altitudes. The lifetime of sub-10 nm particles would be short near the ground level due to evaporative loss (e.g., Fushimi et al., 2008) and/or coagulation scavenging onto pre-existing larger particles, but it could be much longer under conditions of low temperatures and reduced concentrations of pre-existing aerosol particles. We propose that emissions of sub-10 nm particles from aircraft under real-world conditions should be properly considered for understanding the impacts of aviation on human health and also for developing aviation emission inventories for regional and global models.
We conducted field measurements of aerosols at an observation point ∼180 m from the centerline of a runway at Narita International Airport. We investigated the characteristics of particle emissions from in-use commercial aircraft under real-world operating conditions, with specific focuses on the contributions of sub-10 nm size ranges to total and non-volatile (350 ∘C heated) particles. We used the UCPC, CPC, SMPS, and EEPS for the measurements of particle number concentrations and size distributions and carefully investigated the performance and consistency of these instruments. The major conclusions are summarized below.
The median values of the total and the non-volatile EI(N2.5), which likely cover the major size range of aircraft emissions, were found to be 1.1×1017 and 5.7×1015 kg per fuel, respectively. The difference in these values (a factor of ∼20) is interpreted as the average contribution of volatile particles. We did not find a systematic dependence of the total particle number EIs on the estimated plume age (∼30–120 s). The true particle number EIs for total and non-volatile particles might be larger by a factor of ∼1.5 and ∼1.7, respectively, compared to the abovementioned values considering the penetration efficiencies through the sampling tubes and the evaporation tube.
More than half the total and the non-volatile particle number EIs in the aircraft take-off plumes were found in the size range smaller than 10 nm for most of the cases analyzed in this study. The significance of sub-10 nm size ranges for the total particles was qualitatively consistent with previous studies, but that for the non-volatile particles was unexpected.
The unheated UCPC, CPC, and EEPS data consistently suggest that the mode diameters of the dEI(N)∕dlog D for nucleation-mode particles were smaller than 10 nm for the majority of the observed take-off plumes.
The EEPS data suggest that the non-volatile particle number and volume EIs originating from soot-mode particles (>20 nm) were much smaller than those reported by previous studies for take-off plumes under real-world operating conditions.
The characteristics of particle emissions may significantly depend on the type of jet engine, the maintenance conditions, and the fuel sulfur content, which are not available in this study. Particle emissions may also depend on other factors including ambient pressure and temperature. These factors should be carefully considered for a more systematic comparison of different studies.
The field measurement data used in this study are available at the Zenodo data repository (https://doi.org/10.5281/zenodo.4279160; Takegawa et al., 2020b).
The supplement related to this article is available online at: https://doi.org/10.5194/acp-21-1085-2021-supplement.
NT, AF, and HS designed the research; NT, AF, KM, YF, and KS performed field observations and collected data; NT, YM, and HS performed laboratory experiments; NT and YM performed data analysis; NT, AF, YF, and HS wrote the paper.
The authors declare that they have no conflict of interest.
We thank Narita International Airport Corporation and Narita International Airport Promotion Foundation for their help during the field observations at Narita International Airport. We also thank Takumi Saotome at the Research Institute for Environmental Strategies, Inc., Makiko Mine and Anna Nagasaki at Tokyo Metropolitan University for their help in the observations and data analysis, and Kenjiro Iida at AIST for useful advice on the evaluation of the UCPC and CPC.
This research has been supported by the Ministry of the Environment, Japan (grant no. Environment Research and Technology Development Fund (5-1709)).
This paper was edited by Andreas Petzold and reviewed by two anonymous referees.
Boies, A. M., Stettler, M. E. J., Swanson, J. J., Johnson, T. J., Olfert, J. S., Johnson, M., Eggersdorfer, M. L., Rindlisbacher, T., Wang, J., Thomson, K., Smallwood, G., Sevcenco, Y., Walters, D., Williams, P. I., Corbin, J., Mensah, A. A., Symonds, J., Dastanpour, R., and Rogak, S. N., Particle emission characteristics of a gas turbine with a double annular combustor, Aerosol Sci. Tech., 49, 842–855, https://doi.org/10.1080/02786826.2015.1078452, 2015.
Brock, C. A., Schröder, F., Kärcher, B., Petzold, A., Busen, R., and Fiebig, M.: Ultrafine particle size distributions measured in aircraft exhaust plumes, J. Geophys. Res., 105, 26555–26567, https://doi.org/10.1029/2000JD900360, 2000.
Durdina, L., Brem, B. T., Schönenberger, D., Siegerist, F., Anet, J. G., and Rindlisbacher, T.: Nonvolatile particulate matter emissions of a business jet measured at ground level and estimated for cruising altitudes, Environ. Sci. Technol., 53, 12865–12872, https://doi.org/10.1021/acs.est.9b02513, 2019.
Fushimi, A., Hasegawa, S., Takahashi, K., Fujitani, Y., Tanabe, K., and Kobayashi, S.: Atmospheric fate of nuclei-mode particles estimated from the number concentrations and chemical composition of particles measured at roadside and background sites, Atmos. Environ., 42, 949–959, https://doi.org/10.1016/j.atmosenv.2007.10.019, 2008.
Fushimi, A., Saitoh, K., Fujitani, Y., and Takegawa, N.: Identification of jet lubrication oil as a major component of aircraft exhaust nanoparticles, Atmos. Chem. Phys., 19, 6389–6399, https://doi.org/10.5194/acp-19-6389-2019, 2019.
Gormley, P. G. and Kennedy, M.: Diffusion from a stream flowing through a cylindrical tube, P. Roy Irish Acad. A, 52, 163–169, 1949.
Hagen, D. E., Whitefield, P. D., Paladino, J. D., and Trueblood, M. B.: Particulate sizing and emission indices for a jet engine exhaust sampled at cruise, Geophys. Res. Lett., 25, 1681–1684, https://doi.org/10.1029/97GL03504, 1998.
Herndon, S. C., Jayne, J. T., Lobo, P., Onasch, T. B., Fleming, G., Hagen, D. E., Whitefield, P. D., and Miake-Lye, R. C.: Commercial aircraft engine emissions characterization of in-use aircraft at Hartsfield-Jackson Atlanta International Airport, Environ. Sci. Technol., 42, 1877–1883, https://doi.org/10.1021/es072029+, 2008.
International Civil Aviation Organization: Doc 10075. Assembly Resolutions in Force (as of 6 October 2016), ICAO, available at: https://www.icao.int/Meetings/a39/Documents/Resolutions/10075_en.pdf (last access: 10 December 2019), 2017.
Jonsdottir, H. R., Delaval, M., Leni, Z., Keller, A., Brem, B. T., Siegerist, F., Schönenberger, D., Durdina, L., Elser, M., Burtscher, H., Liati, A., and Geiser, M.: Non-volatile particle emissions from aircraft turbine engines at ground-idle induce oxidative stress in bronchial cells, Commun. Biol., 2, 90, https://doi.org/10.1038/s42003-019-0332-7, 2019.
Kärcher, B.: Formation and radiative forcing of contrail cirrus, Nat. Commun., 9, 1824, https://doi.org/10.1038/s41467-018-04068-0, 2018.
Kärcher, B. and Voigt, C.: Susceptibility of contrail ice crystal numbers to aircraft soot particle emissions, Geophys. Res. Lett., 44, 8037–8046, https://doi.org/10.1002/2017GL074949, 2017.
Kärcher, B., Turco, R. P., Yu, F., Danilin, M. Y., Weisenstein, D. K., Miake-Lye, R. C., and Busen, R.: A unified model for ultrafine aircraft particle emissions, J. Geophys. Res., 105, 29379–29386, https://doi.org/10.1029/2000JD900531, 2000.
Kinsey, J. S.: Characterization of emissions from commercial aircraft engines during the Aircraft Particle Emissions eXperiment (APEX) 1 to 3, EPA-600, R-09,130, EPA, available at: https://cfpub.epa.gov/si/ (last access: 25 July 2019), 2009.
Kinsey, J. S., Squier, W., Timko, M., Dong, Y., and Logan, R.: Characterization of the fine particle emissions from the use of two Fischer–Tropsch fuels in a CFM56-2C1 commercial aircraft engine, Energy Fuels, 33, 8821–8834, https://doi.org/10.1021/acs.energyfuels.9b00780, 2019.
Lee, D. S., Pitari, G., Grewe, V., Gierens, K., Penner, J. E., Petzold, A., Prather, M. J., Schumann, U., Bais, A., Berntsen, T., Iachetti, D., Lim, L. L., and Sausen, R.: Transport impacts on atmosphere and climate: Aviation, Atmos. Environ., 44, 4678–4734, https://doi.org/10.1016/j.atmosenv.2009.06.005, 2010.
Liati, A., Brem, B. T., Durdina, L., Vögtli, M., Dasilva, Y. A. R., Eggenschwiler, P. D., and Wang, J.: Electron microscopic study of soot particulate matter emissions from aircraft turbine engines, Environ. Sci. Technol., 48, 10975–10983, https://doi.org/10.1021/es501809b, 2014.
Lobo, P., Hagen, D. E., and Whitefield, P. D.: Measurement and analysis of aircraft engine PM emissions downwind of an active runway at the Oakland International Airport, Atmos. Environ., 61, 114–123, https://doi.org/10.1016/j.atmosenv.2012.07.028, 2012.
Lobo, P., Durdina, L., Smallwood, G. J., Rindlisbacher, T., Siegerist, F., Black, E. A., Yu, Z., Mensah, A. A., Hagen, D. E., Miake-Lye, R. C., Thomson, K. A., Brem, B. T., Corbin, J. C., Abegglen, M., Sierau, B., Whitefield, P. D., and Wang, J.: Measurement of aircraft engine non-volatile PM emissions: Results of the Aviation-Particle Regulatory Instrumentation Demonstration Experiment (A-PRIDE) 4 campaign, Aerosol Sci. Tech., 49, 472–484, https://doi.org/10.1080/02786826.2015.1047012, 2015a.
Lobo, P., Hagen, D. E., Whitefield, P. D., and Raper, D.: PM emissions measurements of in-service commercial aircraft engines during the Delta-Atlanta Hartsfield Study, Atmos. Environ., 104, 237–245, https://doi.org/10.1016/j.atmosenv.2015.01.020, 2015b.
Masiol, M. and Harrison, R. M.: Aircraft engine exhaust emissions and other airport-related contributions to ambient air pollution: A review, Atmos. Environ., 95, 409–455, https://doi.org/10.1016/j.atmosenv.2014.05.070, 2014.
Moore, R. H., Shook, M. A., Ziemba, L. D., DiGangi, J. P., Winstead, E. L., Rauch, B., Jurkat, T., Thornhill, K. L., Crosbie, E. C., Robinson, C., Shingler, T. J., and Anderson, B. E.: Take-off engine particle emission indices for in-service aircraft at Los Angeles International Airport, Sci. Data, 4, 170198, https://doi.org/10.1038/sdata.2017.198, 2017a.
Moore, R. H., Thornhill, K. L., Weinzierl, B., Sauer, D., D'Ascoli, E., Kim, J., Lichtenstern, M., Scheibe, M., Beaton, B., Beyersdorf, A. J., Barrick, J., Bulzan, D., Corr, C. A., Crosbie, E., Jurkat, T., Martin, R., Riddick, D., Shook, M., Slover, G., Voigt, C., White, R., Winstead, E., Yasky, R., Ziemba, L. D., Brown, A., Schlager, H., and Anderson, B. E.: Biofuel blending reduces particle emissions from aircraft engines at cruise conditions, Nature, 543, 411–415, https://doi.org/10.1038/nature21420, 2017b.
Oberdörster, G., Oberdörster, E., and Oberdörster, J.: Nanotoxicology: An emerging discipline evolving from studies of ultrafine particles, Environ. Health Persp., 113, 823–839, https://doi.org/10.1289/ehp.7339, 2005.
Ohlwein, S., Kappeler, R., Joss, M. K., Künzli, N., and Hoffmann, B.: Health effects of ultrafine particles: a systematic literature review update of epidemiological evidence, Int. J. Publ. Health, 64, 547–559, https://doi.org/10.1007/s00038-019-01202-7, 2019.
Onasch, T. B., Jayne, J. T., Herndon, S., Worsnop, D. R., Miake-Lye, R. C., Mortimer, I. P., and Anderson, B. E.: Chemical properties of aircraft engine particulate exhaust emissions, J. Propul. Power, 25, 1121–1137, https://doi.org/10.2514/1.36371, 2009.
Petzold, A., Döpelheuer, A., Brock, C. A., and Schröder, F.: In situ observations and model calculations of black carbon emission by aircraft at cruise altitude, J. Geophys. Res., 104, 22171–22181, https://doi.org/10.1029/1999JD900460, 1999.
Petzold, A., Fiebig, M., Fritzsche, L., Stein, C., Schumann, U., Wilson, C. W., Hurley, C. D., Arnold, F., Katragkou, E., Baltensperger, U., Gysel, M., Nyeki, S., Hitzenberger, R., Giebl, H., Hughes, K. J., Kurtenbach, R., Wiesen, P., Madden, P., Puxbaum, H., Vrchoticky, S., and Wahl, C.: Particle emissions from aircraft engines – a survey of the European project PartEmis, Meteorol. Z., 14, 465–476, https://doi.org/10.1127/0941-2948/2005/0054, 2005.
Righi, M., Hendricks, J., and Sausen, R.: The global impact of the transport sectors on atmospheric aerosol: simulations for year 2000 emissions, Atmos. Chem. Phys., 13, 9939–9970, https://doi.org/10.5194/acp-13-9939-2013, 2013.
Righi, M., Hendricks, J., and Sausen, R.: The global impact of the transport sectors on atmospheric aerosol in 2030 – Part 2: Aviation, Atmos. Chem. Phys., 16, 4481–4495, https://doi.org/10.5194/acp-16-4481-2016, 2016.
Saffaripour, M., Tay, L.-L., Thomson, K. A., Smallwood, G. J., Brem, B. T., Durdina L., and Johnson, M.: Raman spectroscopy and TEM characterization of solid particulate matter emitted from soot generators and aircraft turbine engines, Aerosol Sci. Tech., 51, 518-531, https://doi.org/10.1080/02786826.2016.1274368, 2017.
Saffaripour, M., Thomson, K. A., Smallwood, G. J., and Lobo, P., A review on the morphological properties of non-volatile particulate matter emissions from aircraft turbine engines, J. Aerosol Sci., 139, 105467, https://doi.org/10.1016/j.jaerosci.2019.105467, 2020.
Saitoh, K., Fushimi, A., Fujitani, F., Morino, Y., Saotome, T., Takegawa, N., and Sera, K.: Elemental composition of jet lubricating oil and fuel, and particles (particle size 10 nm–10 µm) collected near the Narita International Airport runway, Nishina Memorial Cyclotron Center Annual Report, 24 (fiscal year 2017), Japan Radioisotope Association, Tokyo, 41–52, 2019a (in Japanese).
Saitoh, K., Fushimi, A., Takegawa, N., and Sera, K.: Quantification of major and trace elements contained in aircraft JET A-1fuel by PIXE analysis, 35th PIXE Symposium, November 2019, Tokyo, 2019b (in Japanese).
Seinfeld, J. H. and Pandis, S. N.: Atmospheric Chemistry and Physics: From Air Pollution to Climate Change, 2nd Edn., Wiley, Hoboken, NJ, 2006.
Society of Automotive Engineers: Procedure for the continuous sampling and measurement of non-volatile particulate matter emissions from aircraft turbine engines: – Aerospace recommended practice 6320, SAE International, available at: https://saemobilus.sae.org/content/ARP6320/#scope, last access: 29 September 2018.
Stacey, B.: Measurement of ultrafine particles at airports: A review, Atmos. Environ., 198, 463–477, https://doi.org/10.1016/j.atmosenv.2018.10.041, 2019.
Takegawa, N. and Sakurai, H.: Laboratory evaluation of a TSI condensation particle counter (Model 3771) under airborne measurement conditions, Aerosol Sci. Tech., 45, 272–283, https://doi.org/10.1080/02786826.2010.532839, 2011.
Takegawa, N., Moteki, N., Oshima, N., Koike, M., Kita, K., Shimizu, A., Sugimoto, N., and Kondo, Y.: Variability of aerosol particle number concentrations observed over the western Pacific in the spring of 2009, J. Geophys. Res.-Atmos., 119, 13474–13488, https://doi.org/10.1002/2014JD022014, 2014.
Takegawa, N., Iida, K., and Sakurai, H.: Modification and laboratory evaluation of a TSI ultrafine condensation particle counter (Model 3776) for airborne measurements, Aerosol Sci. Tech., 51, 235–245, https://doi.org/10.1080/02786826.2016.1261990, 2017.
Takegawa, N., Seto, T., Moteki, N., Koike, M., Oshima, N., Adachi, K., Kita, K., Takami, A., and Kondo, Y.: Enhanced new particle formation above the marine boundary layer over the Yellow Sea: Potential impacts on cloud condensation nuclei, J. Geophys. Res.-Atmos., 125, e2019JD031448. https://doi.org/10.1029/2019JD031448, 2020a.
Takegawa, N., Misawa, K., Fushimi, A., Fujitani, Y., Saitoh, K., Sakurai, H., and Murashima, Y.: Field measurements of aerosol particles near a runway at Narita International Airport, Japan (Version 1.0) (Data set), Zenodo, https://doi.org/10.5281/zenodo.4279160, 2020b.
Timko, M. T., Onasch, T. B., Northway, M. J., Jayne, J. T., Canagaratna, M. R., Herndon, S. C., Wood, E. C., Miake-Lye, R. C., and Knighton, W. B.: Gas turbine engine emissions: Part II: Chemical properties of particulate matter, J. Eng. Gas Turb. Power, 132, 061505, https://doi.org/10.1115/1.4000132, 2010.
Timko, M. T., Fortner, E., Franklin, J., Yu, Z., Wong, H.-W., Onasch, T. B., Miake-Lye, R. C., and Herndon, S. C.: Atmospheric measurements of the physical evolution of aircraft exhaust plumes, Environ. Sci. Technol., 47, 3513–3520, https://doi.org/10.1021/es304349c, 2013.
Wang, X. L., Grose, M. A., Avenido, A., Stolzenburg, M. R., Caldow, R., Osmondson, B. L., Chow, J. C., and Watson, J. G.: Improvement of engine exhaust particle sizer (EEPS) size distribution measurement - I. Algorithm and applications to compact aerosols, J. Aerosol Sci., 92, 95–108, https://doi.org/10.1016/j.jaerosci.2015.11.002, 2016.
Wang, Y., Liu, S .C., Anderson, B. E., Kondo, Y., Gregory, G. L., Sachse, G. W., Vay, S. A., Blake, D. R., Singh, H. B., and Thompson, A. M.: Evidence of convection as a major source of condensation nuclei in the northern midlatitude upper troposphere, Geophys. Res. Lett., 27, 369–372, https://doi.org/10.1029/1999GL010930, 2000.
Westerdahl, D., Fruin, S. A., Fine, P. L., and Sioutas, C.: The Los Angeles International Airport as a source of ultrafine particles and other pollutants to nearby communities, Atmos. Environ., 42, 3143–3155, https://doi.org/10.1016/j.atmosenv.2007.09.006, 2008.
Wiedensohler, A., Birmili, W., Nowak, A., Sonntag, A., Weinhold, K., Merkel, M., Wehner, B., Tuch, T., Pfeifer, S., Fiebig, M., Fjäraa, A. M., Asmi, E., Sellegri, K., Depuy, R., Venzac, H., Villani, P., Laj, P., Aalto, P., Ogren, J. A., Swietlicki, E., Williams, P., Roldin, P., Quincey, P., Hüglin, C., Fierz-Schmidhauser, R., Gysel, M., Weingartner, E., Riccobono, F., Santos, S., Grüning, C., Faloon, K., Beddows, D., Harrison, R., Monahan, C., Jennings, S. G., O'Dowd, C. D., Marinoni, A., Horn, H.-G., Keck, L., Jiang, J., Scheckman, J., McMurry, P. H., Deng, Z., Zhao, C. S., Moerman, M., Henzing, B., de Leeuw, G., Löschau, G., and Bastian, S.: Mobility particle size spectrometers: harmonization of technical standards and data structure to facilitate high quality long-term observations of atmospheric particle number size distributions, Atmos. Meas. Tech., 5, 657–685, https://doi.org/10.5194/amt-5-657-2012, 2012.
Wimmer, D., Lehtipalo, K., Franchin, A., Kangasluoma, J., Kreissl, F., Kürten, A., Kupc, A., Metzger, A., Mikkilä, J., Petäjä, T., Riccobono, F., Vanhanen, J., Kulmala, M., and Curtius, J.: Performance of diethylene glycol-based particle counters in the sub-3 nm size range, Atmos. Meas. Tech., 6, 1793–1804, https://doi.org/10.5194/amt-6-1793-2013, 2013.
Yu, Z., Liscinsky, D. S., Fortner, E. C., Yacovitch, T. I., Croteau, P., Herndon, S. C., and Miake-Lye, R. C.: Evaluation of PM emissions from two in-service gas turbine general aviation aircraft engines, Atmos. Environ., 160, 9–18, https://doi.org/10.1016/j.atmosenv.2017.04.007, 2017.
Yu, Z., Timko, M. T., Herndon, S. C., Miake-Lye, R. C., Beyersdorf, A. J., Ziemba, L. D., Winstead, E. L., and Anderson, B. E.: Mode-specific, semi-volatile chemical composition of particulate matter emissions from a commercial gas turbine aircraft engine, Atmos. Environ., 218, 116974, https://doi.org/10.1016/j.atmosenv.2019.116974, 2019.
Zhang, X., Chen, X., and Wang, J.: A number-based inventory of size-resolved black carbon particle emissions by global civil aviation, Nat. Commun., 10, 534, https://doi.org/10.1038/s41467-019-08491-9, 2019.