Development of a new correction algorithm applicable to any filter-based absorption photometer

Among the various measurement approaches to quantify the light absorption coefficient (Babs), filter-based absorption photometers are dominant in monitoring networks around the globe. Numerous correction algorithms have been introduced to minimize the artifacts due to the presence of the filter in these instruments. However, from our recent studies conducted during the Fire Influence on Regional and Global Environments Experiment (FIREX) laboratory campaign, corrected filter-based Babs remains biased high by roughly a factor of 2.5 when compared to a reference value using a photoacoustic instrument for biomass burning emissions. Similar overestimations of Babs from filterbased instruments exist when implementing the algorithms on 6 months of ambient data from the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) user facility from 2013 (factor of roughly 3). In both datasets, we observed an apparent dependency on single-scattering albedo (SSA) and the absorption Ångström exponent (AAE) in the agreement between Babs based on existing correction factors and the reference Babs. Consequently, we developed a new correction approach that is applicable to any filter-based absorption photometer that includes light transmission from the filter-based instrument as well as the derived AAE and SSA. For the FIREX and SGP datasets, our algorithm results in good agreement between all corrected filter-based Babs values from different filter-based instruments and the reference (slopes ≈ 1 and R2 ≈ 0.98 for biomass burning aerosols and slopes ≈ 1.05 and R2 ≈ 0.65 for ambient aerosols). Moreover, for both the corrected Babs and the derived optical properties (SSA and AAE), our new algorithms work better or at least as well as the two common correction algorithms applied to a particle soot absorption photometer (PSAP). The uncertainty of the new correction algorithm is estimated to be ∼ 10 %, considering the measurement uncertainties of the operated instruments. Therefore, our correction algorithm is applicable to any filter-based absorption photometer and has the potential to “standardize” reported results across any filter-based instrument.

Abstract. Among the various measurement approaches to quantify the light absorption coefficient (B abs ), filter-based absorption photometers are dominant in monitoring networks around the globe. Numerous correction algorithms have been introduced to minimize the artifacts due to the presence of the filter in these instruments. However, from our recent studies conducted during the Fire Influence on Regional and Global Environments Experiment (FIREX) laboratory campaign, corrected filter-based B abs remains biased high by roughly a factor of 2.5 when compared to a reference value using a photoacoustic instrument for biomass burning emissions. Similar overestimations of B abs from filterbased instruments exist when implementing the algorithms on 6 months of ambient data from the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) user facility from 2013 (factor of roughly 3). In both datasets, we observed an apparent dependency on single-scattering albedo (SSA) and the absorption Ångström exponent (AAE) in the agreement between B abs based on existing correction factors and the reference B abs . Consequently, we developed a new correction approach that is applicable to any filter-based absorption photometer that includes light transmission from the filter-based instrument as well as the derived AAE and SSA. For the FIREX and SGP datasets, our algorithm results in good agreement between all corrected filter-based B abs values from different filter-based instruments and the reference (slopes ≈ 1 and R 2 ≈ 0.98 for biomass burning aerosols and slopes ≈ 1.05 and R 2 ≈ 0.65 for ambient aerosols). Moreover, for both the corrected B abs and the derived optical properties (SSA and AAE), our new algorithms work better or at least as well as the two common correction algorithms applied to a particle soot absorption photometer (PSAP). The uncertainty of the new correction algorithm is estimated to be ∼ 10 %, considering the measurement uncertainties of the operated instruments. Therefore, our correction algorithm is applicable to any filter-based absorption photometer and has the potential to "standardize" reported results across any filter-based instrument.

Introduction
Light-absorbing atmospheric aerosols directly affect the Earth's energy budget by absorbing solar radiation, leading to a warming effect when they are suspended in the atmosphere and to the melting of snow and ice following deposition (Bond and Bergstrom, 2006;Boucher, 2015;Horvath, 1993). For decades, scientists have conducted field experiments around the globe to investigate how absorbing aerosols influence the atmospheric radiative balance and interact with clouds (e.g., Andrews et al., 2011;Cappa et al., 2016;Lack et al., 2008b;Rajesh and Ramachandran, 2018;Schwarz et al., 2008). These experiments may be performed at fixed stations (e.g., observation sites maintained by the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) program or the National Oceanic and Atmospheric Administration (NOAA) Global Monitoring Division -GMD) or on mobile platforms (e.g., car trailer, aircraft, ship), typically involving measurements of aerosol chemical, physical, and optical properties. Crucial to the quantification of the radiative forcing of absorbing aerosols are measurements of the absorption coefficient (B abs ). For example, long-term monitoring of B abs provides essential data to evaluate chemistry-climate model simulations (e.g., Chen et al., 2019;Vignati et al., 2010), while intensive measurements of B abs during short-term field campaigns allow for the investigation of optical properties that govern features of aerosol forcing (e.g., McMeeking et al., 2014;Olson et al., 2015).
A variety of instruments have been used to measure B abs , which are generally classified into two large categories: filterbased techniques and photoacoustic techniques (Lack et al., 2014;Moosmüller et al., 2009). The major difference between the two categories of techniques is that B abs is measured after the aerosols are deposited on the filter media in filter-based instruments, while the aerosols are characterized within an airstream in photoacoustic instruments. Compared to filter-based instruments, photoacoustic instruments have the advantage of avoiding potential artifacts due to the contact of aerosols with filters; therefore, they are often used as the reference instruments in intercomparison studies of aerosol absorption (e.g., Arnott et al., 2005;Davies et al., 2019;Jiang et al., 2018;Li et al., 2019;Schmid et al., 2006;Sheridan et al., 2005).
Filter-based absorption photometers have been widely used at observational sites around the world due to their ease of operation and relatively low cost. Numerous instruments can be classified as filter-based absorption photometers including the Radiance Research Particle Soot Absorption Photometer (PSAP), the NOAA Continuous Light Absorption Photometer (CLAP), the Brechtel Manufacturing Tricolor Absorption Photometer (TAP), the Magee Scientific Aethalometer (AETH), and the Thermo Scientific Multi-Angle Absorption Photometer (MAAP). Operationally, all of these instruments are similar in that aerosols are deposited onto a filter and the reduction in the transmission (Tr) of light by the particles (sometimes called attenuation -ATN) is used to infer B abs . Where the instruments may differ is that some are multiwavelength (multi-λ) instruments (e.g., 3λ-PSAP, CLAP, TAP, 7λ-AETH models), while others are not (e.g., 1λ-PSAP, other AETH models, MAAP).
One challenge with filter-based absorption photometers is that biases can arise due to the presence of the filter. For example, light scattering by particles loaded onto the filter or by the filter itself may affect the transmission of light (e.g., Arnott et al., 2005;Bond et al., 1999); non-absorbing material may result in absorption enhancement (e.g., Cappa et al., 2008), or organic vapors adsorbed to the filter may absorb light (e.g., Subramanian et al., 2007). Consequently, various correction algorithms exist to minimize these biases, but they are often specific only to certain instruments. For example, some are applicable to the PSAP, CLAP, and TAP (e.g., Bond et al., 1999;Müller et al., 2014;Ogren, 2010;Virkkula, 2010;Virkkula et al., 2005), while others are applicable to the AETH (e.g., Arnott et al., 2005;Collaud Coen et al., 2010;Drinovec et al., 2017;Kirchstetter and Novakov, 2007;Schmid et al., 2006;Virkkula et al., 2007Virkkula et al., , 2015Weingartner et al., 2003).
Although the equations associated with these existing correction algorithms are different, they share some commonalities. For example, filter-based absorption photometers are assessed using laboratory (e.g., ammonium sulfate, fullerene soot) or ambient aerosols during experiments, which include reference measurements of B abs . These reference measurements often include either direct photoacoustic B abs or inferred B abs as the difference between the extinction coefficient (B ext ) and the scattering coefficient (B scat ). Correction equations are developed by comparing data between the filter-based instrument and the reference instrument, and the equations often contain one term that accounts for filter loading effects and another that accounts for multiple-scattering effects. Consequently, the correction equations frequently incorporate both Tr and either B scat or the single-scattering albedo (SSA) to account for these effects. However, even when the correction algorithms are applied, potential issues can remain, such as the following.

Comparisons between the reference instrument and
4. Derived products (such as absorption Ångström exponents -AAEs) may differ based on the implemented correction algorithm (e.g., Backman et al., 2014;Davies et al., 2019).
5. The agreement between measurements of B abs and estimates of B abs by chemistry-climate models may vary based on the implemented correction algorithm (e.g., Alvarado et al., 2016).
The first three issues in this list may arise due to differences in aerosol optical properties between those used in deriving the correction equation and those associated with a given aerosol sample, and these issues can propagate through to the fourth issue. The final issue is arguably most important because the evaluation of chemistry-climate models may be severely affected by the differences between different correction algorithms, which may inhibit the modeling community from providing accurate projections of future temperature and precipitation response.
In this work, we seek to address some of these issues. First, we evaluate the CLAP, TAP, and PSAP using two common PSAP-based correction algorithms, namely Bond et al. (1999) as updated by Ogren (2010) and  as updated by Virkkula (2010). For brevity, we refer to these corrections as B1999 and V2005 for Bond et al. (1999) and Virkkula et al. (2005), respectively, incorporating their respective updates. In addition, we propose new correction algorithms that are applicable to any filter-based absorption photometer (e.g., CLAP, TAP, PSAP, and AETH) across multiple wavelengths by combining observed filterbased B abs with B scat (e.g., from a colocated nephelometer -NEPH) and reference B abs (e.g., from a colocated photoacoustic instrument). However, in reality (e.g., at longterm observatories), reference values of B abs are rare, and in some cases, complementary B scat measurements may not exist; consequently, we also provide methods to correct filterbased B abs data in these scenarios. To our knowledge, this is the first study to simultaneously evaluate B1999 and V2005 corrections on PSAP "successors" (i.e., CLAP and TAP) and to present a correction algorithm that is broadly applicable to any filter-based absorption photometer. Regarding the latter, even if our correction algorithm has its own limitations, its use can nevertheless standardize the reporting of B abs in long-term datasets.

Methodology
We developed the general form for our correction algorithms using CLAP and TAP measurements collected from biomass burning (65 fires in total) during the Fire Influence on Regional to Global Environments Experiment (FIREX) laboratory campaign in 2016. By using biomass burning emissions, we considered a dataset spanning a broader range of aerosol optical properties (SSA at 652 nm: 0.14-0.98; AAE: 1.25-4.73) than has traditionally been used in developing these correction algorithms. We then conducted a further evaluation and validation of the model using ambient data, specifically using CLAP measurements from the DOE ARM Southern Great Plains (SGP) user facility in Lamont, OK, USA (1 February to 9 July 2013). Our algorithms were then extended to the AETH data from the FIREX laboratory campaign and the PSAP data collected at the SGP site to verify the generalizability of the algorithms.

Experimental setup
In October and November of 2016, we participated in the laboratory portion of the FIREX campaign to investigate wildfire smoke and its impact on the atmosphere. During the campaign, over 100 burns took place at the US Forest Service's combustion facility at the Fire Sciences Laboratory (FSL). The fuels burned in this study are representative of western US ecosystems, such as spruce, fir, various pines, and "chaparral" biome (e.g., manzanita, chamise; see Koss et al., 2018, andSelimovic et al., 2018, for more details).
A typical burn lasted for 1-3 h depending on the smoke sampling strategies (e.g., stack burns versus room burns). During each burn, one or multiple "snapshots" of smoke (typical B abs at 652 nm ranged from 100 to 1200 Mm −1 ) were transferred from the combustion room at FSL into a mixing chamber (210 L) through a long transfer duct (30 m in length, 8 in diameter). The smoke was then diluted by filter air (∼ 230 liters per minute -LPM) in the chamber. Once the concentration in the chamber was stable (detected by the photoacoustic extinctiometer (PAX), which was operated continuously through all fires), the smoke was passed to a suite of instruments to obtain aerosol-and gas-phase parameters. This chamber also served as an intermediate between the transfer duct and the instrumentation to minimize potential biases that arose due to different sample flow rates and sample locations of the instruments. A more detailed description of our experiments can be found in Li et al. (2019).

Measurements of aerosol optical properties
During the campaign, five instruments provided measurements of B abs (CLAP, NOAA GMD; TAP, Brechtel Manufacturing Inc. -BMI; Aethalometer, model AETH-31, Magee Scientific; two PAXs, model PAX-870 and PAX-405, Droplet Measurement Technologies), and two instruments provided measurements of B scat (PAX-870 and PAX-405). The instruments included in the present work are summarized in Table 1.
Both CLAP and TAP provide B abs measurements of the particles deposited on a filter, similar to PSAP. Different from PSAP, there are multiple filter spots (eight sample spots and two reference spots) cycling through one filter in CLAP and TAP, enabling the instruments to run continuously through two or three burns without changing filter. In the CLAP and TAP, sample illumination is provided by light-emitting diodes (LEDs) operated at three wavelengths (467, 528, and 652 nm). Here, we apply both B1999 and V2005 to CLAP and TAP data, similar to previous work (e.g., Backman et al., 2014;Davies et al., 2019).
The key differences between the CLAP and TAP during the FIREX campaign include the following.
1. The spot change in the CLAP was manually performed when Tr reached approximately 0.5 (or ATN decreased to ∼ 69), while the TAP advanced to a new spot automatically with a Tr threshold set to be 0.5.
3. The CLAP recorded B abs every minute, while the TAP recorded B abs every 10 sec. To enable the following analysis, we compute the 1 min averages of TAPderived parameters.  (Bond et al., 1999) 0.3 (Springston, 2016) a The detection limits of PAX and PASS-3 are determined by Allan deviation analysis (Allan, 1966) of B abs during "background zero". b During the analysis of the data collected at the SGP, we use B abs derived by the PASS and B scat derived by the NEPH to yield the coefficients in the algorithms. c Two types of filters were used during the FIREX campaign (see Sect. 2.1.2). d The operating wavelengths of CLAP, TAP, and PSAP are stated slightly different by the instrument manufacturers. We simply use 467, 528, and 652 nm throughout this paper.
the CLAP, while Azumi filters (model 371M, Azumi Filter Paper Co., Japan) were used in the second portion of the campaign (due to a lack of availability of the Pallflex filters). The TAP was equipped exclusively with the Azumi filters throughout the campaign. We apply the filter correction recommended in  to the CLAP and convert from Pallflex to Azumi filters.
5. BMI substantially reengineered the CLAP in their development of the TAP.
These differences resulted in variable agreement between the CLAP and TAP during FIREX; however, the two instruments did largely agree within experimental uncertainty (see, e.g., Figs. S8 and S13 in Li et al., 2019).
A PAX measures B abs and B scat simultaneously for suspended particles using a modulated diode laser. We use these photoacoustic absorption measurements as the reference to evaluate the filter-based B abs and develop our correction algorithms. To enable the evaluation of CLAP and TAP, which operate at different wavelengths than the PAXs, we interpo-late the measurements of B abs and B scat to the wavelengths of 467, 528, and 652 nm using the values of AAE and scattering Ångström exponents (SAEs), similar to Backman et al. (2014) and Virkkula et al. (2005). Theoretically, AAE and SAE fit absorption and scattering as power-law functions of wavelength (Bergstrom et al., 2007).
Due to the numerous correction algorithms for the Aethalometer (e.g., Arnott et al., 2005;Collaud Coen et al., 2010;Kirchstetter and Novakov, 2007;Saturno et al., 2017;Schmid et al., 2006;Virkkula et al., 2007;Weingartner et al., 2003), we do not evaluate these in the present work to limit the scope. In fact, the majority of our focus is the B1999 and V2005 corrections to TAP and CLAP. However, we still test the performance of the new algorithms on the AETH to explore its applicability to that instrument.

Measurements of aerosol optical properties at the SGP observatory
The ambient data used in this paper are the ground-based aerosol data measured at the SGP observatory from 1 February to 9 July 2013 (archived at https://www.archive.arm.gov/ discovery/, last access: 16 January 2019). For evaluation purposes, we randomly select a range of dates during which the observations are valid (without incorrect, suspect, and missing data) and the PASS is operated after laser (532 nm) upgrade. At the site, an impactor was used to switch the sampling between two cutoffs (particle diameter < 10 µm (PM 10 ) in the first 30 min of each hour and < 1 µm (PM 1 ) in the latter 30 min of each hour). The aerosols exiting the impactor were dried to RH less than 40 % and passed to a CLAP, a PSAP, and two NEPHs. Moreover, a three-wavelength photoacoustic soot spectrometer (PASS-3) was operated at the site and measured B abs and B scat of the aerosols, but these aerosols did not pass through the impactor (e.g., characterizing total suspended particles -TSPs). Typical B abs and B scat reported at the site ranged from 0 to 10 and 0 to 50 Mm −1 at 550 nm, respectively (e.g., Sherman et al., 2015). Although the site is rural (clean background air), long-term transport aerosols (such as mineral dust, absorbing organic aerosols, and secondary organic aerosol -SOA) may affect the local aerosol properties (Andrews et al., 2019).
We preprocess the SGP data in three steps. First, due to the systematic difference in aerosol sizes between PASS-derived and filter-based absorption, we only include the PM10 observations, inherently assuming that any differences in the optical properties of PM 10 and TSPs are negligible. Then, we smooth the 1 sec data into 10 min averages. Thirdly, we estimate the detection limits at each of the three wavelengths in the PASS-3 using the data measured during the background zero periods (Allan, 1966) and discard the observations below the detection limits. With a 10 min averaging time, the detection limits (3σ ) for the PASS-3 are 0.78 Mm −1 (405 nm), 2.01 Mm −1 (532 nm), and 0.30 Mm −1 (781 nm). For the filter-based instruments, the detection limits are based on previous studies (see Table 1). Moreover, we only retain the observations that satisfy B abs (405 nm) > B abs (532 nm) > B abs (781 nm) (or AAE > 0), similar to Fischer and Smith (2018). As with the PAX data from the laboratory, we adjust the PASS-derived B abs to 467, 528, and 652 nm using the inferred AAE values for each 10 min average.

Calibrations
Following Bond et al. (1999) and Ogren et al. (2010), the filter-based instruments were calibrated and corrected for sample area, flow rate, and filter type (see Li et al., 2019) for the FIREX data and Sherman et al. (2015) for the SGP data. Other than that, we did not do any verification beyond the manufacturer's calibration for the filter-based instruments. The SGP nephelometer measurements were corrected for truncation effects (Sherman et al., 2015). The FIREX photoacoustic measurements were calibrated by ammonium sulfate aerosol and fullerene soot (Li et al., 2019).

The correction algorithms
In filter-based instruments, the light intensities transmitted through the sample spot and blank spot of the filter are recorded as I s and I b , respectively. The logarithmic ratio of the two intensities at time t is defined as ATN using the Beer-Lambert law: where ATN = 0 when beginning a new filter spot (t = 0). The ATN can be related to Tr by normalizing I s /I b at time t relative to I s /I b at the start of a new filter spot (t = 0): ). (2) The change in ATN over a time interval ( t) for the instrument operated at a volume flow rate of Q and spot area of A yields the attenuation coefficient (B ATN ) for that time interval: B ATN is finally converted to B abs by applying correction algorithms. The general form of the correction algorithms presented for the PSAP in Bond et al. (1999) and  can be summarized as where f (Tr) is some function of Tr (that may vary between approaches), correcting for the filter loading effect. C 1 is a constant that may vary with wavelength; specifically, it is a penalty for light scattering by the particles collected on the filter, which may contribute to the quantification of ATN. In most atmospheric and laboratory studies, B scat is measured independently, typically using a colocated NEPH. Besides the abovementioned B1999 and V2005 corrections, the constrained two-stream (CTS) correction proposed by Müller et al. (2014) can be also applied on PSAP-similar instruments. The CTS correction was developed based on the relationship between absorption optical depth and B abs . However, it is not straightforward to reformulate the CTS correction as a function of Tr as those presented in Eqs. (5) and (6). Thus, we exclude the B abs results corrected by the CTS correction in the following analysis. For those who are interested in the B abs results corrected by the CTS correction, we provide the correction results of our SGP-CLAP data in Supplement Fig. S9. Bond et al. (1999) were the first study to present a correction algorithm for filter-based instruments. This empirical correction was originally developed for the PSAP operated at 550 nm using various mixtures of laboratory-generated nigrosin (SSA ≈ 0.5) and ammonium sulfate (SSA ≈ 1) with B abs ranging from 0 to 800 Mm −1 .

The B1999 correction
After calibrating the flow rate and spot area of the PSAP, Bond et al. (1999) derived C 1 = 0.016 and defined f (Tr) as follows: Following an adjustment by Ogren (2010), C 2 and C 3 in Eq. (5) were updated to 1.55 and 1.02, respectively. These are the values used in the present work (Table 2) for B1999. Moreover, Ogren (2010) stated that the correction forms of Eqs. (4) and (5) were valid for any wavelength, while additional experiments were needed to establish the equation parameters for wavelengths other than 574 nm.

The V2005 correction
Virkkula et al. (2005) developed a correction algorithm for both three-wavelength PSAP (467, 530, and 660 nm) and one-wavelength PSAP (574 nm) using the same functional form as Eq. (4). Since the operating wavelengths of the photoacoustic instruments and the NEPH were different from those of the PSAP, the measured photoacoustic B abs and B scat were extrapolated or interpolated to 467, 530, and 660 nm using inferred AAE and SAE, respectively. In this study, the authors used various mixtures of kerosene soot, ammonium sulfate, and polystyrene latex (SSA ranged from 0.2 to 0.9) with B abs ranging from 0 to 800 Mm −1 at 530 nm.
Different from the f (Tr) in the B1999 correction, which was a reciprocal function of Tr, the f (Tr) presented in V2005 was a multivariate linear function of the natural logarithm of Tr and SSA (including an interaction term between the two): where the parameters in Eq. (6) vary with wavelength. The parameters in V2005 were updated in Virkkula (2010) by correcting for flowmeter calibration (Table 2). Due to the unknown values of SSA before deriving B abs , Virkkula et al. (2005) provided a solution through an iterative procedure. In the iteration, B abs is first calculated using the B1999 correction (e.g., Eqs. 4 and 5) and is then used to compute the initial guess of SSA for use in Eq. (6). The B abs and SSA can be updated using Eqs. (4) and (6) until convergence is reached.

Refitting the coefficients in B1999 and V2005
With the reference measurements of B abs from the photoacoustic instruments, we are able to refit the coefficients in the B1999 and V2005 corrections (C 2 to C 7 in Eqs. 5 and 6) using our data. Specifically, we use the Levenberg-Marquardt algorithm (1944) to iteratively fit the coefficients until the chi square of the coefficients is minimized (Levenberg, 1944). The fitting is implemented using the Curvefit function in Igor Pro. It is worth noting that the derived coefficients may only be valid for the SGP and FIREX data. For aerosol properties different from our study, the optimal coefficients are likely to be different from the ones reported here. Hereafter, the B1999 and V2005 results with refitted coefficients are referred to as "updated B1999" and "updated V2005", respectively.

The new correction
We develop a set of new correction algorithms with the same general form as Eq. (4) using the biomass burning emissions from 65 different burns during the FIREX laboratory study, providing a broader range of aerosol optical properties and aerosol concentrations than previous work. This was motivated by the disagreement that remained between filter-based and photoacoustic instruments, even after applying B1999 to the data (see, e.g., Li et al., 2019, their Fig. 4 and our Fig. 2). These differences may persist because we were effectively extrapolating the B1999 correction equation to values outside the range for which it was developed.
This new correction is developed based on multiple linear regression techniques with three dependent variables of ln(Tr), SSA, and AAE as well as one independent variable of B abs /B ATN (Eqs. 7-9). As with other correction equations, this model takes into account the influence of scattering and weakly absorbing materials, but we have the additional aim of developing a model that is applicable to any filter-based absorption photometer.
Similar to the B1999 and V2005 corrections, this new model starts with the general form of Eq. (4), rewritten here  See Tables 4 and S6-S10   a The operating wavelengths are based on the manufacturer specifications. b The coefficients provided in Table 2 are the values presented in Ogren (2010) and Virkkula (2010), which are updated from Bond et al. (1999) and Virkkula et al. (2005), respectively. c We reformulate the correction equations in the original publications to agree with Eqs. (4)- (6)  to define B scat in terms of SSA and B abs .
Rearranging this equation to move all B abs terms to the lefthand side yields where g (Tr(λ), SSA(λ)) = f (Tr(λ)) × 1−SSA(λ) 1−(1−C 1 )×SSA(λ) . We define a new function g that can be used in Eq. (8). Specifically, we construct a multivariate linear model for g, introducing AAE as a dependent variable and including interaction terms between SSA, AAE, and ln(Tr): Equation (9) suggests that different combinations of SSA, AAE, and ln(Tr) can result in the same value of g (i.e., B abs /B ATN ); likewise, a given value of B abs /B ATN may have infinitely many points with distinct slopes passing through it (Fig. S3). Although Eq. (9) is developed based on statistical approaches, we attempt to relate this statistical model to physical effects. The coefficients G 1 − G 3 are fairly straightforward, as these account for the influence of filter loading (G 1 ), relative light scattering by the aerosols (G 2 ), and the brownness of the aerosols (G 3 ). The interaction terms (G 4 − G 7 ) are more difficult to assign a physical meaning; however, the interaction between filter loading and relative light scattering (G 4 ) might be interpreted as an absolute light scattering by the aerosols on the filter, while the interaction between filter loading and aerosol brownness (G 6 ) is somewhat analogous to G 4 . The three-way interaction between filter loading, scattering, and the brownness of aerosols (G 7 ) is required because of the three two-way interaction terms.
To further this physical interpretation of our statistical model (Eq. 9), we explore the relationship between B abs B ATN and ln(Tr), which essentially follows a y = m × x + b form, where y is B abs B ATN and x is ln(Tr). The slope (m) is defined as G 1 +G 4 × SSA + G 6 × AAE +G 7 × SSA × AAE, and the intercept (b) is defined as G 2 × SSA + G 3 × AAE +G 5 × SSA × AAE. Therefore, different combinations of SSA and AAE modulate this relationship between B abs B ATN and ln(Tr). For example, loading "black" particles on the filter (e.g., AAE ∼ 1 and SSA ∼ 0.3) tends to produce larger values of B abs B ATN , while loading "white" particles on the filter (e.g., AAE ∼ 3 and SSA ∼ 0.9) tends to produce smaller values of B abs B ATN (see Fig. S3 of this work and Fig. 4 in Virkkula et al., 2005). This relationship becomes more complex when considering, e.g., mixed sulfate and black carbon particles; SSA can be high while AAE is low, and the corresponding B abs B ATN can be variable (also see Fig. S3). Therefore, in order to properly compensate for the effects of loading and aerosol optical properties, a multiple linear regression with interaction terms is introduced in Eq. (9).
A detailed description of the procedure for the model development (e.g., variable transformation -from Tr to ln(Tr), variable selection using the best subsets and stepwise approaches, and model assessment) is provided in the Supplement. We evaluate the model by plotting B abs B ATN against aerosol properties not included in Eq. (9) (such as relative humidity and aerosol geometric mean diameter, which have been previously reported to bias corrections of filter-based B abs ; Moteki et al., 2010;Nakayama et al., 2010;Schmid et al., 2006). The results are presented in Figs. S5-S7.
As in V2005, iteration is required in our algorithm because B abs is dependent on knowledge of SSA and AAE, which themselves are dependent on B abs . We propose the following iterative process to update SSA and AAE in the model.
1. Initialize AAE from B ATN across the three wavelengths (B ATN ∼ λ −AAE ) and initialize SSA for each wavelength using B ATN from the filter-based absorption photometer and B scat from a colocated NEPH, i.e., 2. Yield an initial set of coefficients G 0 through G 7 for each wavelength to calculate g(Tr, SSA, AAE) in Eq. (9) using one of the algorithms described in Sect. 2.5.
4. Update AAE and SSA using B abs calculated in Step 3.

Application of correction algorithms
In developing a procedure for applying our algorithm, we envision three potential scenarios.
1. Algorithm A. The filter-based instrument is colocated with a NEPH and reference instrument providing B abs . This scenario facilitates the computation of G 0 through G 7 in Eq. (9) (Step 2 in the iterative process) as well as the derivation of new coefficients for existing correction algorithms. This scenario can also enable the development of a new a set of coefficients that may be more appropriate for aerosol sources that we do not consider here.
2. Algorithm B. The filter-based instrument is colocated with a NEPH but not a reference instrument providing B abs , which is perhaps the most likely scenario (at least at many long-term monitoring sites). This scenario requires an initial guess of the coefficients; we provide sets of these in Table 4 below for different filter-based instruments and aerosol sources.
3. Algorithm C. The filter-based instrument is deployed with neither a colocated NEPH nor a reference instrument providing B abs . This scenario is the most challenging because there are no measurements of B scat to compute SSA; to address this issue, we propose the use of a nonlinear relationship between SSA and AAE (AAE = a + b× SSA c ) to provide an initial guess of SSA in the iterations.
To aid in decision-making between algorithms, we developed a flowchart for selecting appropriate correction algorithms for CLAP, TAP, and PSAP (Fig. 1). Furthermore, programs based on Igor Pro (WaveMetrics, Inc.) and Python (version 3.7.5) for selecting and implementing our correction algorithms can be found at https://doi.org/10.5281/zenodo.3742342.

Application of the previous algorithms on different aerosols
We first consider the application of the B1999 and V2005 corrections on different combinations of aerosol type and filter-based absorption photometer. Specifically, we apply the two corrections to the biomass burning data from the FIREX laboratory campaign (CLAP and TAP) and 6 months of ambient data from the SGP site (CLAP and PSAP). In doing so, we use the "default" coefficients recommended in B1999 and V2005 as well as "updated" coefficients that are estimated via regression techniques. It is important to keep in mind that the updated coefficient values of B1999 and V2005 (Table S7) are only valid for the aerosols investigated in this study. Future experiments are needed to systematically determine how the coefficients in B1999 and V2005 may change for different aerosol types. We focus on the results of the CLAP in the main text because a CLAP is the only instrument common to deployments for both FIREX and SGP. The results of the TAP from FIREX and the PSAP from the SGP site can be found in the Supplement (Table S5 and Fig. S11). Our intercomparison between the corrected CLAP-derived B abs and reference B abs for the FIREX and SGP data is provided in Fig. 2 and Table 3. For the FIREX measurements, both analyses (using the default coefficients and updating the coefficients) suggest good correlation (coefficient of determination (R 2 ) > 0.9) between the CLAP and the reference across all three wavelengths. Nevertheless, the corrections using the default coefficients result in overprediction of B abs by factors of ∼ 2.5. If we update the coefficients in the corrections, there is an obvious improvement in the agreement (i.e., slope ≈ 1; R 2 increases). The results are generally similar for SGP, although the R 2 for ambient data is generally lower for ambient data (R 2 < 0.7). Decreased R 2 may be due to the lower aerosol concentrations measured in ambient air, which could lead to lower signal-to-noise ratios in the instruments. Moreover, it is worth mentioning that for both datasets (FIREX and SGP), the corrected B abs values from different filter-based absorption photometers using the default approaches do not agree with each other (slopes range from 0.69 to 1.40). However, after updating the coefficients, the slopes approach unity (Table S6).
In the FIREX data, there is an apparent dependency of the updated coefficients on the wavelength of light but, more importantly, on the aerosol optical properties, namely SSA and AAE (Tables S7-S9). However, in the ambient data from SGP, the dependency on optical properties is less obvious (Tables S10-S11). Nevertheless, all of these coefficients differ from those reported in B1999 and V2005 (again, derived for the PSAP rather than the CLAP), which highlights the potential need to use coefficient values that are appropriate for the instrument being used, its wavelength(s) of light, and optical properties that are representative of the sampled aerosols when applying correction factors to B ATN .

Application of the new algorithms to the FIREX data
The colocation of the CLAP, TAP, AETH, and PAX during FIREX allows us to apply each algorithm (A, B, C) to these data. Similar to Sect. 3.1, we focus our discussion on the CLAP, with details on the TAP and AETH presented in the Supplement (Fig. S11-S12). However, we provide the recommended initial guesses in the new algorithms and the comparison of absorption (corrected filter-based B abs versus reference B abs ) for all filter-based absorption photometers in Tables 4 and 5 to help readers quickly retrieve key information on our algorithms. Figure 3 provides a comparison between the uncorrected B ATN from the CLAP at all three wavelengths, as well as photoacoustic B abs interpolated to those wavelengths using AAE. For each wavelength, the slopes are significantly greater than 1. Moreover, there is an apparent dependency on SSA and AAE in the agreement between the instruments. This is most obvious in Fig. 3a (652 nm), where data with lower SSA and lower AAE (smaller markers, "brighter" colors) fall below the best-fit line, while data with higher SSA and higher AAE (larger markers, "darker" colors) fall above the best-fit line. This phenomenon is less clear in Fig. 3bc, but an apparent dependency on SSA and AAE remains, which highlights the need to include both of these aerosol optical properties (and appropriate interaction terms) when correcting B ATN values.
We first apply Algorithm A to the CLAP B ATN data in Fig. 3. Using the reference B abs values from the PAX (in addition to B scat values), we are able to derive a set of coeffi-  Table 3. Relationship between the CLAP-derived B abs corrected by the B1999 and V2005 algorithms (including updated coefficients) and the reference B abs at 652, 528, and 467 nm. The relationship is achieved using major axis regression (Ayers, 2001). The value in parentheses represents the coefficient of determination (R 2 ) of the linear relationship. cients that enable the correction of the data (Table 4). Corrected CLAP values are presented in Fig. 4, with the linear relationships presented in Table 5. The slope for each wavelength is very close to the 1 : 1 line, suggesting that our approach works well in correcting these data. Moreover, the heteroscedasticity that exists in Fig. 3 has been minimized after correction, and there are no apparent trends in how the data are organized in Fig. 4 due to the aerosol optical properties.
We next investigate the repeatability of the coefficient values presented in Table 4 by randomly selecting half of the measurements (N = 1338) from the whole FIREX dataset. By implementing Algorithm A for the extracted observations, we obtain new coefficient values for G 0 to G 7 . This is repeated 1000 times to obtain a distribution of coefficient values (Fig. S13). The extraction approach mimics the process of obtaining new biomass burning datasets so that we can estimate the variability of these derived coefficients. From Fig. S13 and Table S12, the derived coefficients are mostly insensitive to the different randomly extracted datasets; most of the quartile deviation (defined as (Q3 − Q1)/2, where Q1 and Q3 are the first and third quartile, respectively) is within 0.05 except G 4 , which has a quartile deviation of ∼ 0.08. Consequently, the coefficient values obtained in Al-gorithm A appear to be reasonable initial guesses to correct filter-based absorption measurements during biomass burning events when the reference B abs is unavailable, such as in Algorithm B and Algorithm C.
We next implement Algorithm B for the CLAP B ATN data from Fig. 3 using the initial guesses of the coefficients derived from Algorithm A (Table 4) along with reference B scat values. To get a sense of the variability in the results, we randomly select half of the data and applied the correction; this process is repeated 1000 times. For each iteration, we compare the corrected B abs from the CLAP to the reference B abs from the PAX; the resulting slope, intercept, and R 2 values are summarized as box-and-whisker plots in Fig. 5. For all three wavelengths, the slopes are close to unity, and there is good correlation between the two absorption measurements (R 2 ≈ 0.98), which indicates that the good performance seen in Fig. 4 is independent of the reference B abs measurements and our algorithm is able to correct "new" B ATN . Consequently, when scattering measurements are colocated with filter-based absorption measurements, our new correction algorithm performs well.
Lastly, we apply Algorithm C to the data in Fig. 3. However, we first require a functional relationship between AAE and SSA because in this scenario, the CLAP B ATN values are  Table 4. Coefficient values for Eq. (9) derived using Algorithm A. We recommend these as the initial guesses when implementing Algorithm B * .   Table 5. Relationship between the filter-based B abs corrected by Algorithm A and the reference B abs at the operating wavelengths for the filter-based instrument. The relationship is achieved using major axis regression (Ayers, 2001). The value in the parentheses represents the coefficient of determination (R 2 ) for the linear relationship. the only data input to the algorithm (and therefore SSA is unknown). Liu et al. (2014) proposed that a power function can describe this relationship (AAE = a + b× SSA c ); we present these data from FIREX along with power function fits (and associated prediction intervals) in Fig. 6. To define AAE in this figure, we fit a power-law relationship to the three B ATN values from the CLAP; similarly, we define SSA using interpolated B scat from the PAX and B ATN from the CLAP (the rationale for using B ATN is that if Algorithm C were to be implemented in practice, only B ATN would be available). In Fig. 6, the data points are colored by "prediction error", effectively a metric to quantify how well the power function reproduces the individual data points. Although there is a fair amount of error in some of these points, we still obtain an SSA-AAE relationship required to initialize Algorithm C. Even though there is uncertainty in the SSA vs. AAE relationship used in Algorithm C, after corrections have been applied, the filter-based B abs for the CLAP agrees well with the independent reference B abs ; the slopes for all wavelengths are slightly greater than 1 (1.03-1.05) and the R 2 values are all high (0.97-0.98). However, even though the absorption measurements are corrected well, there are still large uncertainties in values of inferred scattering. Examples of this are provided in Fig. 7, where we compare the SSA inferred from the PAX to the SSA inferred from Algorithm C as well as B scat for each wavelength. Generally speaking, data that are better represented by the SSA vs. AAE relationship (i.e., smaller prediction error) result in better agreement with the reference for both SSA and B scat , but there is also a clear divergence from the 1 : 1 line in Fig. 7a-c as SSA decreases. Therefore, even though Algorithm C performs well at correcting filter-based B ATN to agree with the reference B abs , estimates of final SSA values should be considered to be uncertain.
To explore the generalization of the new algorithms across different instruments, we next apply our algorithms to the other filter-based absorption photometers operated during the FIREX study (TAP and AETH). Consistent with what we observed for the CLAP results, the corrected TAP-and AETHderived B abs is in good agreement with the photoacoustic B abs (as demonstrated in Tables 4 and 5, as well as Figs. S11-S12). Moreover, the corrected B abs values from the three Figure 4. As in Fig. 3, but the CLAP-based B ATN values have been corrected using our Algorithm A. Figure 5. The box-and-whisker plots for the slope, intercept, and R 2 of the relationship between the CLAP-derived B abs (corrected by Algorithm B in the present work) and PAX-derived B abs for all three wavelengths. For details on how these values were generated, please refer to the text.
filter-based instruments agree with each other for all three wavelengths (Table 6), confirming the consistency of our algorithm.

Application of the new algorithms to ambient data
To test our algorithms further, we extended our work to ambient data collected at the DOE SGP site in the time period during which the PASS-3 was operational. From the SGP data, we derived a different set of coefficients for ambient data using Algorithm A, which differ from those derived for FIREX (Table 4). The results presented in Fig. 8 and Table 5 suggest that our new algorithm works at least as well as B1999 and V2005 on this dataset (both with updated coefficients).
The repeatability of the coefficient values in Algorithm A is confirmed for the SGP measurements using the same procedure as described in Sect. 3.2 (see results in Fig. S13 and Table S12).
For the SGP data, we see similar performance to the FIREX data when we apply Algorithm B, whereby we again sampled half of the CLAP data, used the initial guesses derived in Algorithm A, and repeated this process 1000 times. Although the slopes tend to be larger than 1 (i.e., the corrected CLAP B abs remains high relative to the PASS B abs ), the results still represent an improvement over B1999 and V2005 using their recommended coefficients for their correction equations. Figure 6. AAE plotted against SSA for the FIREX data. In the figures, AAE was computed using a power-law fit across all three wavelengths, and SSA was computed using the interpolated B scat from the two PAXs and the reported B ATN from the CLAP. The data points are colored by their prediction error (true AAE minus calculated AAE divided by calculated AAE). Implementing Algorithm C is challenging for ambient data because there is no distinct power function relationship in AAE vs. SSA (Fig. 9); this is consistent with other field studies reporting both SSA and AAE (e.g., Backman et al., 2014 andLim et al., 2018). Our approach described here is only appropriate for ambient aerosols that follow a power function, such as sites impacted by biomass burning. Nevertheless, we did apply this to a subset of the SGP data for which the AAE-SSA prediction error is within 30 % (N = 86), and for this subset of data, Algorithm C works fairly well (slopes ≈ 0.95; see Fig. S14). Therefore, while Algorithm C may have utility for ambient data, we advise caution when using this algorithm since the aerosols influencing the site may not be represented by a clear AAE-SSA power function (e.g., when biomass burning and coarse aerosols are equally prevalent at a long-term monitoring site).
These new algorithms are also applicable to the PSAP deployed at the SGP site. The results of the correction for the PSAP are presented in Table 5 and Fig. S11, and the recommended initial guesses when implementing Algorithm B to PSAP B ATN in ambient environments are given in Table 4. As expected, there is good agreement between corrected PSAP and CLAP B abs (Table 6).
To investigate if our algorithms are suitable to correct B abs obtained from different ambient environments, the aerosol properties from the SGP site are compared to those from the other NOAA/ESRL observational sites. We use a similar AAE-SAE space as that in Cappa et al. (2016) and Schmeisser et al. (2017) to infer dominant aerosol types (e.g., BC, BrC, dust, mixed dust-BC-BrC; see Fig. S15). Our results of AAE and SAE overlap most of the values from the NOAA/ESRL sites, except when marine aerosols or dust contribute to the local aerosol emissions. Though none of the NOAA/ESRL sites fall into the clusters of BC or BC-BrC, some of our data can represent the optical properties of aerosols from these clusters. Therefore, we highlight that our algorithm developed by the SGP data may have the potential to be generalized to a variety of environmental conditions, but we would need to validate this using observations from more studies.

Impact of the implemented correction algorithm on aerosol optical properties
In addition to the direct comparisons of B abs between the filter-based and photoacoustic measurements, we compare derived optical properties (AAE and SSA) from different instruments to assess the algorithms' performance on derived aerosol optical properties. For example, we have discussed the discrepancy of SSA between the filter-based and photoacoustic measurements when implementing Algorithm C in Sect. 3.2. In this section we will more broadly discuss the impact of different correction algorithms on AAE and SSA. In Fig. 10, we present the frequency distribution of AAE for both FIREX and SGP data generated from different campaign and/or instrument pairs using different correction approaches. For the FIREX data ( Fig. 10a-b), most corrections (with the exception of the default B1999) are consistent with the photoacoustic data, while for the SGP data ( Fig. 10c-d), most corrections (with the exception of default V2005) are consistent with the photoacoustic data. However, updating the coefficients for B1999 and V2005 improves the agreement with the photoacoustic data. The 50 % difference that exists between the B1999 and V2005 algorithms in all panels in Fig. 10 is consistent with previous studies. For example, both Backman et al. (2014) and Davies et al. (2019) found that V2005-derived AAE is greater than B1999-derived AAE by 33 % to 50 % for ambient aerosols. Therefore, we highlight that the default coefficients in B1999 and V2005 may have some limitations when deriving AAE using the corrected B abs ; instead, updating the coefficients or using the new algorithm proposed in this work may yield more robust AAE results.
Similar to Fig. 10, we also investigate the distribution of SSA computed by using corrected B abs along with B scat . We provide the results at 652 nm as an example in the main text (Fig. 11); figures for 528 and 467 nm can be found in the Supplement (Figs. S16 and S17). For both FIREX and SGP data, the SSA obtained using the new algorithm agrees very well with the B1999 and V2005 but only when their coefficients have been updated. Calculations of SSA using B1999 and V2005 with their recommended coefficients suggest that these values may be biased low, which follows the overestimation of corrected B abs demonstrated in Fig. 2.
Moreover, we plot similar figures as Figs. 10-11 using all algorithms (A, B, and C). As shown in Fig. S18, the results using Algorithm B agree very well with those using Algorithm A, but the use of Algorithm C results in some obvious discrepancies compared to the photoacoustic reference, again highlighting the potential for large uncertainty using this algorithm.
In Fig. 12, we directly compare the distributions of both AAE and SSA at 652 nm for all of the filter-based absorption photometers considered here using our Algorithm A to correct the B ATN data. For both datasets, after the corrections have been applied, there are only marginal differences in the AAE ( Fig. 10a and b) derived by different instruments. Similarly, there is good agreement among the SSA values when using corrected B abs from different instruments ( Fig. 10c and  d). Overall, the derived properties using the new correction  are consistent across all instruments, suggesting its generalizability.

Uncertainty of the new algorithms
In this section, we estimate the uncertainty of the new algorithms due to both measurement uncertainties of the instruments and the uncertainties of parameter computation. We then simulate the propagated uncertainty in the corrected filter-based B abs reported in this paper.
Measurement uncertainties of the instruments considered here have been reported in previous work (e.g., Anderson et al., 1996;Nakayama et al., 2015;Sherman et al., 2015) and are summarized in Table 1. The typical sources of measurement uncertainty of the aerosol instruments include (1) instrument noise (often associated with the averaging time); (2) calibration uncertainties (such as the accuracy of the operating wavelengths and the properties of the calibration materials); (3) standard temperature and pressure (STP) correction uncertainties (Sherman et al., 2015); and (4) flow rate uncertainties. Additional uncertainties that are specific to filter-based absorption photometers include spot size and filter medium corrections (Bond et al., 1999;. Regardless, these values all tend to be ≤ 30 %, which is consistent with other commonly used aerosol instrumentation. Because correction algorithms for filter-based absorption instruments also require aerosol optical properties, the algorithms' performance will be affected by these values as well. For example, uncertainties in SSA are directly related to uncertainties associated with B abs and B scat , which are both included in our simulations. However, capturing uncertainties in AAE is more complex, as AAE can be computed by either a 2λ fit (a linear fit using B abs at two wavelengths) or 3λ fit (the same as the power fit used in the present work). Davies et al. (2019) used the 3λ fit to calculate AAE and compared this to calculations using 662 and 785 nm (i.e., AAE 662/785 ), finding that the 3λ results were about 50 % greater. Moreover, similar differences (−35 % to 85 %) can exist when comparing two different 2λ combinations (AAE 440/870 and AAE 675/870 ), depending on the contribution of brown carbon to absorption at 440 nm . However, based on Figs. S19 and S20, we demonstrate small (< ∼ 10 %) dif-  ferences in the calculated values of AAE using our Algorithm A with different 2λ combinations for linear fits and the 3λ power-law fit when considering both FIREX and SGP data. Consequently, we do not include AAE calculation uncertainty in our simulation.
In our simulations, the propagated uncertainty of corrected B abs is estimated by implementing the new algorithm for datasets in which filter-based B ATN , reference B abs , and B scat are subject to measurement uncertainties. The full procedure is outlined in the Supplement, but we provide a brief overview of our Monte Carlo approach here. First, we create a synthetic dataset (n = 500 records) that defines B abs at 652 nm and AAE that is intended to represent biomass burning. Values of B ATN and SSA are then computed using the relationships presented in Figs. 3 and 6, respectively. Respective uncertainties associated with each of these values are applied following Table 1, assuming that these follow a normal distribution. We then applied Algorithm B to the B ATN dataset, repeated 1000 times, to quantify the overall uncertainty associated with our correction algorithm. Figure 13 provides a graphical summary of our uncertainty simulation results, which was derived by fitting linear equations to the "true" B abs value (that we defined) and the "corrected" B abs values (outputs of each iteration). Considering Figure 13. The box-and-whisker plots (slope, intercept, and R 2 ) for the Monte Carlo simulation of the relationship between the CLAPderived B abs (corrected by Algorithm B in the present work) and "true" B abs for all three wavelengths. the slopes (Fig. 12a), our algorithm can generally reproduce the true value within 10 % at 652 and 528 nm, but the performance is slightly degraded at 467 nm. The median intercept for our simulations is close to zero, but the interquartile range increases with decreasing wavelength (Fig. 12b), suggesting that the uncertainty may increase at shorter wavelengths. The coefficients of determination (Fig. 12c) range from 0.47 (652 nm) to 0.68 (467 nm), showing that the algorithm may be less precise if large measurement uncertainties exist. Even though these sources of uncertainty exist when implementing our correction algorithms and propagate through to the corrected values, we argue that our new algorithm will standardize uncertainties across corrected B abs values from filterbased absorption photometers. Moreover, the new algorithms perform at least better than the previous algorithms with default coefficients or as well as the previous algorithms with updated coefficients.

Conclusions
Filter-based absorption instruments are widely used at global observational sites due to their relatively low cost, fast response, and easy operation. Despite the existence of different correction algorithms to correct the filter-based B abs measurements, these are not "standardized" as differences in corrected B abs values exist across different instrumentcorrection combinations, even when the instruments are colocated. This study provides a systematic evaluation of the previous correction algorithms (B1999 and V2005 corrections) on the CLAP and similar instruments (TAP and PSAP) using both laboratory-generated biomass burning emissions and ambient aerosols. We also developed new correction algorithms that are applicable to any filter-based absorption photometer (e.g., PSAP, CLAP, TAP, AETH), which will have utility for any historic or future filter-based absorption measurements and which have the potential to standardize absorption coefficients across all filter-based instruments. This latter point is demonstrated in Table 6 and Fig. 12 in that there is good agreement across all filter-based absorption photometers when applying our corrections to both biomass burning and ambient data. In practice, we anticipate that our Algorithm B will be the most common because at long-term monitoring sites, filter-based absorption photometers are typically colocated with a nephelometer.
Using the existing corrections on our CLAP measurements, we find that the corrected B abs overestimates photoacoustic B abs by factors of ∼ 2.6 (biomass burning aerosols) and ∼ 3.2 (ambient aerosols). Similar overestimations of absorption by filter-based instruments are seen in the results of TAP from the FIREX study and PSAP deployed at the SGP. Comparing B1999 and V2005, B abs values corrected by the two corrections differ by −6 % to 18 %. These discrepancies in our results are consistent with those reported for intercomparisons between filter-based and photoacoustic absorption instruments (e.g., Arnott et al., 2003;Davies et al., 2019;Li et al., 2019;Müller et al., 2011a).
Overall, our new developed algorithms (A, B, and C) perform well in correcting B abs for different filter-based absorption photometers (CLAP, TAP, PSAP, and AETH) from both biomass burning and ambient measurements. Our work suggests that if the filter-based instrument is co-operated with a reference absorption instrument and a NEPH in the field for a period, researchers can compute site-specific initial guesses (the same as Algorithm A in the present work). Otherwise, either Algorithm B or Algorithm C proposed in this paper can be used to correct the filter-based measurements. In Algorithm B when a filter-based absorption photometer is colocated with a NEPH but without a reference instrument, the set of coefficients yielded in this work (Table 4) can be used as initial guesses to implement the algorithm. In Algorithm C when a filter-based absorption photometer is operated by itself, a "representative" relationship between AAE and SSA can be used to estimate SSA from AAE at each step in the iterative process, but we advise caution if this relationship is not monotonic (e.g., as in the ambient data from SGP and from Backman et al., 2014, andLim et al., 2018). The only scenario not included in the present work is when the filterbased absorption photometer is colocated with a reference absorption instrument but no instrument for scattering. However, under this scenario, one could simply use the photoacoustic B abs data because no filter-induced biases exist for those instruments.
In terms of the aerosol optical properties (AAE and SSA) computed by different corrections, the new algorithm suggests no bias of AAE and SSA when compared to that derived by updated B1999 and updated V2005 for both aerosol datasets.
However, the new algorithm is not without limitations. First, we used the photoacoustic B abs as the reference to develop the algorithm and the initial guess of the coefficients; meanwhile, some studies argue that photoacoustic absorption is not a "ground truth" (e.g., Lack et al., 2006;Lewis et al., 2008). Thus, we simulate the propagated uncertainty of our algorithms considering the measurement uncertainties due to the photoacoustic B abs (as well as B ATN and B scat ) and find that the corrected B abs can be biased by −17 % to 5 %, depending on the operated wavelength. Although potential bias due to the precision of photoacoustic B abs cannot be excluded, using our new algorithm to correct the filterbased B abs will at least eliminate correction-related biases among different filter-based instruments. Second, we only tested the algorithms with data from biomass burning and ambient measurements. It is unclear how the algorithms will work for other absorbing aerosols (e.g., dominated by fossil fuel emissions or mineral dust). Further evaluation of the performance of the new algorithm on other aerosol sources may help to address this issue. Regardless, we argue that our approach can standardize reported absorption coefficients at long-term monitoring sites, which has the potential to yield a better dataset with which to evaluate chemistry-climate models.
Code and data availability. The code for the algorithm was initially written in Igor Pro (WaveMetrics Inc.) and has been developed to be accessible through either Python (version 3.7.5) or Igor Pro. The implementation of our algorithm is available at https://doi.org/10.5281/zenodo.3742342 (Li et al., 2020). The FIREX aerosol products are available at the NOAA CSL data archive https://esrl.noaa.gov/csl/projects/firex/firelab/ (National Oceanic and Atmospheric Administration, 2020). The SGP aerosol products are available at the DOE ARM data archive https:// www.archive.arm.gov/discovery/ (US Department of Energy, 2020) (February 2013-July 2013; 36 • 36 18.0 N, 97 • 29 6.0 W: Southern Great Plains Central Facility).
Author contributions. HL developed the model, conducted the evaluations, and drafted the paper. AAM and GRM edited the paper and obtained funding.