Using Multi-Platform Satellite Observations to Study the Atmospheric Evolution of Brown Carbon in Siberian Biomass Burning Plumes

: A bulk of evidence from in situ observations and lab experiments suggests that brown carbon (light-absorbing organic compounds in particles) can provide a signiﬁcant yet highly variable contribution to the overall light absorption by aerosol particles from biomass burning (BB). Partly stemming from the complexity of the atmospheric evolution of organic aerosol (OA), the variability in brown carbon (BrC) absorption makes it difﬁcult to partition the radiative effects of BrC and black carbon (BC) in atmospheric and climate models; as such, there are calls for satellite-based methods that could provide a statistical characterization of BrC absorption and its evolution in different regions of the world, especially in remote BB regions, such as Siberia. This study examined the feasibility of the statistical characterization of the evolution of BrC absorption and related parameters of BB aerosol in smoke plumes from intense wildﬁres in Siberia through the analysis of a combination of data from three satellite instruments: OMI (Ozone Monitoring Instrument), MISR (Multi-Angle Imaging SpectroRadiometer), and MODIS (Moderate Resolution Imaging Spectroradiometer). Using a Monte Carlo method, which related the satellite retrievals of the absorption and extinction aerosol optical depths to Mie theory calculations of the optical properties of BB aerosol, we found that the BrC absorption, as well as the imaginary refractive index for the OA, decreased signiﬁcantly in Siberian BB smoke plumes during about 30 h of the daylight evolution, nevertheless remaining considerable until at least 70 h of the daylight evolution. Overall, the study indicated that the analysis of multi-platform satellite observations of BB plumes can provide useful insights into the atmospheric evolution of BrC absorption and the partitioning of BrC and BC contributions to the total light absorption by BB aerosol.


Introduction
Aerosol particles emitted from open biomass burning (BB) constitute a major fraction of aerosol emissions from combustion sources and, as such, contribute significantly to the radiative budget of the atmosphere and the aerosol radiative forcing, both at global and regional scales, through the direct aerosol-radiation interaction [1][2][3] by impacting the formation and lifetime of clouds [4][5][6] and affecting the carbon fluxes in ecosystems [7]. The direct interaction of BB aerosol with solar light is mostly determined by a balance between light scattering by organic aerosol (OA) and light absorption by black carbon (BC) and brown carbon (BrC), with the latter being comprised of a complex mixture of lightabsorbing organic compounds. Strong variability in the sources, chemical composition, and mixing state of BB aerosol results in large uncertainty in the available estimates of the aerosol radiative forcing due to the anthropogenic impact on BB emissions [8][9][10]. Among the potentially important yet poorly investigated factors contributing to this uncertainty is insufficient knowledge of the direct radiative effect of BrC, which is estimated at the top of the atmosphere to range globally from 0.03 to 0.6 W m 2 (e.g., [11][12][13][14]).
The broad range of the available estimates of the radiative effect of BrC reflects a high variability in both the OA mass fraction in BB smoke and absorptive properties of particulate organic matter, as well as the fact that BrC is formed by chemically active organic species, which can be produced or destroyed during the atmospheric evolution of BB aerosol. Findings from lab studies suggest that the magnitude of the BrC contribution to the light absorption by fresh BB aerosol strongly depends on the fuel type and the fire regime, generally increasing with a decrease in the modified combustion efficiency [15,16]. There is also abundant evidence based on both atmospheric observations and lab studies showing that the absorption enhancement due to BrC can strongly vary during the atmospheric evolution of BB aerosol, either increasing or decreasing (e.g., [17][18][19][20][21][22]).
The strong variability of the BrC content across individual samples of BB aerosol calls for observational information that would be representative of typical BrC absorption properties and their time evolution in different regions of the world [13,23]. Such information could further be used as a constraint on BB aerosol absorption properties in regional and global models, in addition to the a priori knowledge of the chemical composition of BB emissions provided as emission factors [24].
Presently, the most abundant data characterizing the contribution of organic compounds to the absorption properties of fresh and aged BB aerosol are available from laboratory experiments (e.g., [15,16,[19][20][21][25][26][27]). The results of these experiments advanced the fundamental knowledge of BB aerosol but their representativeness of real atmospheric conditions is not always clear. Useful insights into possible regional differences in characteristics of the BrC absorption in BB plumes, indicating that these differences can be significant, are provided by field observations of the optical and physical properties of BB aerosol [17,22,[28][29][30][31]. In particular, the analysis of BB aerosol samples collected by aircraft in BB plumes over several regions in the U.S. revealed that BrC gradually vanished with the BB aerosol age, almost disappearing after about 1 day of atmospheric exposure [17]. In contrast, an aircraft campaign conducted in Western Africa indicated that BrC formation in BB plumes strongly dominated over its loss, at least during the first 12 h of the atmospheric evolution [22]. However, such field observations are sparse and have very limited temporal and spatial coverage. In particular, to the best of our knowledge, there have been no field campaigns investigating the evolution of BrC in Siberia, which is one of the largest source regions of BB aerosol in the world [32].
Numerous studies derived information on BrC absorption from remote sensing measurements at Aerosol Robotic Network (AERONET) by exploiting the fact that BrC absorption cross-section rapidly increases with decreasing wavelength in the visible and near-UV regions [33][34][35][36][37][38][39][40]. Although the AERONET data have been shown to provide useful insights into the evolution of BrC [18,38], the quality assured AERONET observations of BB plumes in remote regions are rather sparse [40,41], which hampers their use for achieving a better quantitative understanding of BrC formation and loss processes. Furthermore, the BrC absorption at the shortest wavelength of the AERONET observations (440 nm) is typically rather weak compared to the absorption at shorter wavelengths [42], which makes the use of AERONET data to constrain the radiative effect of BrC even more challenging.
Several studies indicated the feasibility of inferring BrC absorption in BB plumes from satellite observations. In particular, accounting for BrC absorption in retrievals from the Ozone Monitoring Instrument (OMI) measurements allowed for improving the agreement of the retrieved values for aerosol optical depth (AOD) and single scattering albedo (SSA) with the corresponding AERONET observations in a biomass burning region in South America [41]. Evidence for a strong BrC absorption in BB plumes in South America from the OMI observations was found to be consistent with similar evidence from the ground-based radiometric measurements in the UV-A spectral region [42]. Furthermore, the OMI-retrieved absorption AOD (AAOD) combined with the extinction AOD derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) observations provided useful insights into the role of BrC in the atmospheric evolution of the AAOD in BB plumes in Siberia [43]. Apart from the OMI observations, satellite observations retrieved using an experimental (research) aerosol retrieval algorithm from the Multi-Angle Imaging Spectrometer (MISR) measurements were also recently shown to provide a meaningful, albeit mostly qualitative rather than quantitative, characterization of BrC absorption in BB plumes [44]. Note that a possible strong wavelength dependence of the imaginary refractive index (IRI) of organic aerosol has not yet been explicitly accounted for in the standard MISR algorithm [45], and therefore, the corresponding MISR data are not informative regarding BrC absorption. Overall, to the best of our knowledge, no methods have been proposed to constrain the BrC absorption in BB plumes by using the standard satellite data products mentioned above.
This study was primarily motivated by the need to constrain the absorption properties of BB aerosol emitted from wildfires in Siberia. Siberian fires are recognized as an important source of BC transported to and deposited in the Arctic [46], where rapid climate changes are observed and various components of atmospheric aerosol are known to contribute significantly to the regional radiative balance [47]. There is, however, much fewer data about the contribution of Siberian fires to non-BC light-absorbing aerosol, which has been found in Arctic snow in significant quantities along with BC [48]. Evidently, this contribution depends not only on primary emissions of BrC and transport patterns but also on the BrC atmospheric evolution, which is, as indicated above, still poorly known.
Given our previous experience of combining data from different satellites to study phenomena associated with biomass burning [43,[49][50][51], we assumed that useful constraints on the BrC absorption can be derived from a combination of observations performed by the three different satellite instruments mentioned above. Accordingly, the goal of this study was to investigate the feasibility of constraining BrC absorption and its evolution in Siberian BB plumes by using the standard aerosol data products derived from the OMI, MISR, and MODIS measurements.
We analyzed the satellite data that were available for a big Siberian region in July 2016. In this period, major forest fires occurred both in the western and eastern parts of Siberia [43,52,53]. The study region and period were determined by taking into account the availability of estimates of the photochemical age of BB aerosol from our previous studies and the need to isolate Siberian forest fires from other fires (such as grass fires in Kazakhstan) that are likely to emit BB aerosol featuring different absorption properties than forest fires.
To the best of our knowledge, this is the first study in which satellite observations of BB aerosol were used for quantitative statistical characterization of BrC absorption and its atmospheric evolution. Hence, despite the temporal and spatial limitations of our analysis, we regard this study as a significant step forward toward creating the satellite measurement-based climatology of BrC absorption in BB plumes in major BB regions.
An overview of the study design is presented in Figure 1. The input data that were used for our analysis are described in Section 2.1. The analysis methods are described in Section 2.2. Our main results are presented in Section 3 and then discussed in Section 4. The objectives, methods, and results of the study are summarized in Section 5.

Satellite Data
In this study, the key information on BrC absorption was inferred from the AAOD values retrieved at the 388 nm wavelength from the OMI observations onboard the EOS Aura satellite. Similar to our previous studies that exploited the OMI AAOD observations [43,50], we used the Level-2 OMAERUV data product [54,55] (v1.8.9.1) available from the NASA EarthData portal (https://search.earthdata.nasa.gov/, accessed on 15 February 2022). This data product has been derived from the differences between the spectral radiances measured at the top of the atmosphere and those calculated for a hypothetical purely molecular atmosphere by using a set of pre-defined aerosol models from a look-up table.
The OMAERUV algorithm segregates the retrievals into three assumed aerosol types (biomass burning, urban, and desert dust), with the OMI observations representative of BB aerosol being identified based on values of the CO column amounts retrieved from AIRS (Atmospheric Infrared Sounder) measurements [54]. Only the data for BB aerosol were used in this study. Another important feature of the OMAERUV (v1.8.9.1) algorithm is that the IRI for BB aerosol particles is assumed to depend on the wavelength, making it possible to account for the BrC absorption [41]. It is also noteworthy that the aerosol optical properties were retrieved by assuming five different locations of the center of mass of the aerosol vertical distribution: 0, 1.5, 3.0, 6.0, and 10 km above the surface, and the data product includes the AAOD retrievals for each of these locations. Hence, additional information on the vertical distribution of aerosol particles was needed to ensure the best fit of the retrievals to the properties of the BB aerosol in real smoke plumes. Following our previous studies involving the OMI AAOD retrievals [43,50], we obtained this information from simulations with a chemistry transport model (see Section 2.1.3). Finally, one more important feature of the OMAERUV data product is that the quality-assured AOD and SSA retrievals are much less abundant than the AAOD retrievals [43,50]. For this reason, the OMI AOD and SSA data were not used in the given study. The nominal spatial resolution of the OMI AAOD retrievals is 13 × 24 km 2 at nadir. Although, to the best of our knowledge, quantitative estimates of the accuracy of the OMI AAOD data have not yet been reported in the literature, the AOD and SSA retrievals were previously compared with AERONET data in different regions of the world and were found to be in good quantitative agreement with them [56,57], with correlation coefficient values of 0.75 or better obtained at 50% of sites in the case of the AOD retrievals. Furthermore, the OMI AAOD data were found to be consistent (either directly or through 3D model simulations) with in situ, aircraft, and remote sensing measurements of BB aerosol in Siberia [50]. Several studies (which were mentioned in the Introduction) indicated that the AAOD values retrieved from the OMI observations are sensitive to the BrC content in BB aerosol [41][42][43] and even enable quantitative characterization of the BrC contribution to the total absorption of light by BB aerosol at 388 nm [43].
Note that the most reliable estimates of the aerosol properties (AAOD, AOD, and SSA) are provided by the OMAERUV algorithm only for the 388 nm wavelength [54], even though the OMAERUV data product reports values of the same properties also for the 354 and 500 nm wavelengths as a result of transformations from 388 nm. Since these transformations are determined by the assumed models of aerosol, the reported values for the 354 and 500 nm wavelengths are regarded as less reliable [54]. Based on these considerations, we opted not to use these values in our analysis. To estimate the absorption Ångström exponent (AAE), which is indicative of the contribution of BrC to the light absorption by BB aerosol [40], we used, along with the OMI retrievals of the AAOD at 388 nm (henceforth, AAOD 388 ), the AAOD at 867 nm (henceforth, AAOD 867 ) derived from the MISR measurements aboard the NASA EOS Terra. Specifically, we used the data from the level 2 version 23 (MIL2ASAE) [45] available with a nominal horizontal resolution of 4.4 km from the NASA LaRC ASDC file system (https://urs. earthdata.nasa.gov/, accessed on 1 October 2021). The MISR algorithm exploits a unique capability of the MISR instrument, namely, multi-angle imaging at four wavelengths (447, 558, 672, and 867 nm), to obtain view-angle and wavelength-dependent equivalent reflectances at the top of the atmosphere (TOA), which are then compared to radiances computed with a coupled atmosphere/surface radiative transfer model for a pre-defined set of 74 aerosol compositions and particle size distributions (aerosol models). The aerosol model that enables the best fit of the computed TOA atmospheric path radiances to those derived from the MISR measurements is selected to provide the column aerosol parameters and properties.
As noted in the Introduction, the standard MISR aerosol algorithm does not take into account a possible wavelength dependence of the IRI of organic aerosol particles. Hence, it is not feasible to derive the BrC absorption (which is typically associated with an increase in the IRI at shorter wavelengths) directly from the MIL2ASAE data. However, since BrC is likely to provide only a very minor contribution to the absorption at the green, red, and near-infrared bands of the MISR measurements, we assumed that the retrieved absorption properties of BB aerosol for the corresponding wavelengths are not significantly biased, even when the BrC contribution to the absorption at the near-UV region is considerable. For definiteness, we used the MISR data for the 867 nm wavelength, assuming that they are the least affected by the actual BrC absorption. The AAOD values needed for our analysis (but not directly provided by the MIL2ASAE data product) were derived from the SSA and AOD retrievals: where all values correspond to the 867 nm wavelength. A preceding version (version 22) of the MISR aerosol product (which has the spatial resolution of 17.6 × 17.6 km 2 ) was extensively validated against AERONET and other ground-based photometer observations, e.g., [58][59][60][61][62], although the validation studies mainly focused on the AOD retrievals. For example, Petrenko and Ichoku [61] found that the coefficient of determination (R 2 ) of the MISR AOD retrievals relative to AERONET observations worldwide was 0.73 in summer, and the root-mean-square error (RMSE) was 0.13. Garay et al. [45] reported evidence that the accuracy of the AOD retrievals provided as part of the version 23 MISR aerosol product is comparable to that of the version 22 AOD retrievals over land and is significantly better over oceans. In contrast to the AOD retrievals, the SSA retrievals available as part of MIL2ASAE have not been extensively evaluated (mainly due to the fact that the SSA data suitable for validation of MISR retrievals are too sparse [62]), and accordingly, have been regarded as semiquantitative [45]. Details on the MISR aerosol retrieval procedure and further information on the validation of the MISR aerosol data product used in this study can be found elsewhere [46].
We emphasize that the SSA values available from MIL2ASAE were involved in our analysis only to derive the AAOD values according to Equation (1), while the use of the SSA as an independent input variable in our estimation algorithm was avoided. We took into account the fact that not only have the retrieved SSA values been insufficiently validated, but they are also rigidly determined by a limited number of aerosol models. These models can impose unnecessary limitations on our estimates of the aerosol absorption characteristics, leading to uncontrollable biases in them. Our particular concern was that these aerosol models are not necessarily always able to fit the possible variability of the size distribution of BB aerosol particles. As this variability can be manifested as strong and wavelength-dependent changes in the aerosol optical properties [18,43], the retrieved absorption parameters and the SSA values for BB aerosol may be affected by systematic errors. However, since the AAOD values are determined not only by the SSA but also by the AOD, the uncertainties in the AAOD are likewise dependent on uncertainties in both the SSA and AOD. Taking into account that the retrieved values of the AOD depend on the optical and microphysical aerosol properties assumed in the same aerosol models that were used to retrieve the SSA and the relation between the retrieved AOD and the intrinsic aerosol properties is likely very complex [45], we expected a complex combination of observation and retrieval errors in the AAOD to behave in a random-like fashion. This expectation was justified in part by a comparison of the satellite data with AERONET retrievals (see Section 2.1.4 below).
Using the AAOD 388 and AAOD 867 values derived from the OMI and MISR observations, we computed the absorption Ångström exponent (AAE 388/867 ) as follows: AAE 388/867 = − log(AAOD 388 /AAOD 867 )/ log(388/867). ( Since the AAE is undefined when the AAOD is zero, all the MISR retrievals corresponding to an SSA of 1 were disregarded. Such selection could, in principle, result in a positive bias in the average of the AAOD 867 values. Nonetheless, as evidenced by our analysis presented below in Section 2.1.4, the AAOD 388 -to-AAOD 867 ratio that determines AAE 388/867 according to Equation (2) is not significantly biased in Siberian BB plumes. This fact may indicate that the MISR retrievals featuring an SSA of 1 are mostly not representative of dense BB plumes in Siberia.
Finally, we also used the AOD at 550 nm (denoted below as AOD 550 ) derived from the MODIS measurements [63] aboard the Aqua and Terra satellites. The AOD data with a nominal horizontal resolution of 10 × 10 km 2 were obtained from the Collection 6.1 Level-2 MODIS data products [64] as the merged "dark target" and "deep blue" AOD retrievals. The MODIS AOD values demonstrated a good agreement with the AERONET-observed AOD values, indicating a high quality of the MODIS aerosol data products [64].
Following our previous study [44], we considered the ratios of the AAOD and AOD values retrieved at different wavelengths as proxies for the SSA, which, along with the AAE, was found to provide a useful constraint on the absorption properties of OA in BB plumes [40]. Unlike the SSA values from the MISR aerosol product, the AAOD and AOD values considered here were not linked through a common aerosol model, and thus their ratio was likely to be less dependent overall on any specific assumptions involved in the retrieval algorithms. As possible "algorithmic" limitations contribute to the overall uncertainties in the AAOD/AOD ratios, along with many other uncertainty sources, they are not expected to preclude fitting these ratios by an independent aerosol model. More generally speaking, the use of the data independently retrieved from observations performed by three (rather than one or two) different satellite instruments in our analysis reduces the possibility of biases that may appear in our estimates of the BrC absorption due to error covariances between the different input variables and is deemed to be necessary for enhancing the reliability of our estimates. All the satellite data described above were projected onto the same 1 • × 1 • grid covering most of Siberia, as well as some adjacent regions in Northern Eurasia. Different observations falling into the same grid cell and corresponding to the same hour were averaged. The spatial distributions of the temporal averages of the AAOD 388 , AAOD 867 , and AOD 550 retrievals projected onto our grid are shown in Figure 2. Since the satellite observations considered here were made with three different instruments from three different satellites, a harmonization procedure was needed to increase the probability that the different observations represent the same BB plumes. To this end, we required that the time difference between the OMI, MISR, and MODIS observations for a given grid cell should not exceed two hours. The available data that did not satisfy this requirement were disregarded. A similar 2 h difference between the different satellite observations was imposed in Konovalov et al. [43]. The remaining inconsistencies and mismatches between the satellite observations performed by the different instruments are expected to result in random errors (noise) in the characteristics inferred in this study, rather than in any persistent biases in them. Furthermore, to minimize the chances that some of the measurements covered the background scenes rather than BB plumes, we disregarded any gridded data points for which AOD 550 was less than 0.8: the indicated threshold value was chosen (as in Konovalov et al. [40]) by taking into account the fact that the "background" AOD in the absence of fires in Siberia is typically less than 0.2 [43,50,65]. This AOD-based selection criterion complements a CO-based criterion applied within the OMAERUV algorithm to the OMI data (see above) and is intended to further reduce the impact of the background aerosol on our estimates of the BrC absorption. To diminish the complex and uncertain effects of hygroscopicity of aerosol particles on our estimates, we also disregarded any observations corresponding to the scenes with the relative humidity (RH) within the aerosol column larger than 70%. The RH was estimated using a chemistry transport model, as described in Section 2.1.3. In addition, we excluded any observations south of 50 • N since those observations could be affected by fires that occurred outside of Siberia. Figure 2. Spatial fields of the temporal averages of (a) AAOD at 388 nm, (b) AAOD at 867 nm, and (c) AOD at 550 nm in the period from 1 to 31 July 2016 according to OMI, MISR, and MODIS observations, respectively. The figure illustrates the satellite data sets used in this study before applying the data harmonization procedure described in Section 2.1.1. Figure 3 illustrates the same data as in Figure 2 but after applying the data harmonization procedure. The amount of data was consequently drastically reduced, which was mostly due to the lack of matched data from OMI and MISR observations. Note that the MODIS AOD data were much more abundant than the OMI and MISR data considered here, and the availability of the MODIS data did not pose any noticeable limitation on the size of the harmonized data set. The remaining data were mostly representative of the hot spot in Siberia. There were also a few remaining data points in the European part of Russia, which were due to the mentioned episode of outflow of BB plumes from Siberia [43,52,53]. The observations illustrated in Figure 3 were used as the input data for our analysis described below. Note that any data pixel shown in Figures 1 and 2 can correspond to several observations taken on different days during the study period. The three sets of the pre-processed and harmonized satellite observations were combined to estimate the three characteristics mentioned above, that is, AAE 388/867 along with the AAOD 388 /AOD 550 and AAOD 867 /AOD 550 ratios. These three characteristics constituted the components of the observation vector for our estimation algorithm described below in Section 2.2.1. The AAOD 388 /AOD 550 and AAOD 867 /AOD 550 ratios are referred to below as the absorption-to-extinction ratios and are abbreviated, respectively, as AER 388/550 and AER 867/550 .

AERONET Data
AERONET data were used in this study for auxiliary purposes, such as validation of the satellite observations (see Section 2.1.4) and testing of our estimation algorithm (see Section 2.2.3). AERONET is a network of sun-sky radiometers located at more than 500 sites worldwide [66]. The radiance measurements are processed using a unique algorithm that enables the retrieval of multiple optical and physical properties of aerosol [67]. In this study, we used quality-assured AAOD retrievals available for 440, 675, and 870 nm wavelengths as part of the Level 2.0 Version 3 inversion products, as well as AOD Level 2.0 data provided for the 500 and 675 nm wavelengths as the Version 3 Direct Sun Algorithm product (https://aeronet.gsfc.nasa.gov/, accessed on 21 October 2021). As noted above, multiple studies used the AERONET observations to make inferences on the aerosol composition and BrC absorption [33][34][35][36][37][38][39][40].
Following our previous study [40] in which the BrC absorption in Siberian BB plumes was analyzed using AERONET data, we used the Level 2 AAOD, SSA, and AOD data from the Tomsk_22 (56.4 • N, 84.7 • E) and Yakutsk (61.7 • N, 129.4 • E) AERONET sites, which are situated in western and eastern Siberia, respectively. For the algorithm testing, we considered all the data available from these sites for the summer months (June-August) by the end of 2020 (as in Konovalov et al. [40]). However, taking into account the availability of the data characterizing the BB aerosol photochemical age from our 3D model simulations, the summer observations for only two years-2012 and 2016-were employed for the validation of the satellite data. Our concept of using the AERONET data for 2012 to validate the satellite data for 2016 is explained below in Section 2.1.4. Note that the evolution of BrC absorption was not analyzed in Konovalov et al. [40]. The selected AERONET data were then used to derive the AOD and AAOD values at the wavelengths of the satellite observations. Specifically, since the 388 nm wavelength, for which the AAOD values are available from the OMI observations, is slightly outside of the wavelength range of the AERONET observations, the AAOD 388 values representing the AERONET observations were obtained from the power-law extrapolation from available values of the AAOD at the 440 and 675 nm wavelengths. The power-law interpolation between the 675 and 870 nm wavelengths was used to estimate the AAOD at the 867 nm wavelength. To select the AERONET observations that were representative of BB plumes, we applied the same AODbased criterion as the one applied to the satellite data. Specifically, we disregarded any data points with an AOD at 550 nm less than 0.8. Further details on these data and examples of the corresponding time series can be found elsewhere [40].

Simulation Data
In this study, we took advantage of the availability of simulated data from two previous studies of Siberian BB aerosol [43,50]. The simulations were performed with the CHIMERE chemistry transport model (v2017) [68] and covered (but were not limited to) the three summer months of 2012 [50] and the period from 15 June to 15 August of 2016 [43]. Note that the simulations for 2012 were used here only for validation purposes (see Section 2.1.4 below) and were not employed to obtain our main results reported in Section 3. The simulation results provided concentrations of multiple gas-phase and aerosol species from various sources (anthropogenic, biogenic, and pyrogenic), along with concentrations of model tracers that allowed for estimating the photochemical age of BB aerosol. However, most of these data, except for concentrations of two BB aerosol age tracers (T 1 and T 2 ), the vertical distribution of the BB aerosol mass concentration, and the relative humidity, were unused in the given observation-based study.
The atmospheric transport of the tracers and pollutants was driven by the meteorological processes simulated using the WRF model [69]. The BB emissions were derived using a well-established method [43,51,70] from satellite observations of the fire radiative power (FRP) [71] retrieved from the MODIS measurements. The FRP data were also used to calculate the maximum injection heights of BB emissions according to the parameterization of Sofiev et al. [72]. The parameterization was validated [72] through a comparison of the calculated injection heights for multiple fires that occurred in Canada, Siberia, and the US during the fire seasons of 2007 and 2008 with observations of the BB plume top using the MISR instrument [73]. Note that although the MISR-derived BB plume heights could, in principle, be used in our analysis instead of the corresponding simulated data, they are not provided as part of the MISR standard aerosol product and were not readily available for the region and period considered in this study. The simulations were performed with a horizontal resolution of 1 • × 1 • for 12 non-equidistant layers that extended up to the 200 hPa pressure level. Further details on the 3D model simulations used in this study can be found elsewhere [43,50]. The age tracers were simulated as hypothetical gases emitted from biomass burning with the same emission rate as particulate organic carbon. The tracer T 1 behaved as an inert gas, while the tracer T 2 decayed with the constant rate (k T ) of 2.3 × 10 −5 s −1 (about 1/12 h −1 ) during daytime (when the local zenith angle was less than 90 • ). The BB aerosol photochemical age (τ a ) was calculated using the columnar mass concentrations, namely, [T 1 ] and [T 2 ], of these tracers as follows [43]: Note that according to Equation (3), the BB aerosol photochemical age is effectively estimated as the solar exposure time. As argued in Konovalov et al. [43], using Equation (3) to estimate the photochemical age implies that the photochemical processes affecting the BB composition and optical properties are active only during daytime and that the OH concentration within a BB plume is constant during the daytime and drops to zero during the nighttime. Such a simplified estimation reflects the lack of quantitative knowledge about the nighttime oxidation processes within BB plumes, as well as about the variability of the OH concentration within BB plumes.
The spatial distribution of the mean BB aerosol photochemical age in July 2016 is shown in Figure 4. The simulated data were selected consistently (in time and space) with the satellite data shown in Figure 3. It can be seen that relatively fresh BB aerosol dominated northeast of Krasnoyarsk and west of Norilsk. In July 2016, these locations roughly coincided with the "hot spots" of BB emissions; conversely, the most aged BB aerosol prevailed in the European part of Russia, where it was transported from Siberia [43]. . Spatial distribution of the mean BB aerosol photochemical aerosol age (h) in July 2016 according to simulations performed with the CHIMERE chemistry transport model [43,50]. The gridded estimates of the photochemical age were selected in time and space consistently with the pre-processed satellite data (see Figure 3).

Validation of the Satellite Data against the AERONET Observations
Ground-based observations that are commonly employed to validate satellite retrievals are usually supposed to be matched to the satellite data, both in time and space. This common approach, however, turned out to be infeasible in this study because the MISR and OMI observations that are both representative of BB aerosol (according to the criteria formulated above) and coincident with the AERONET observations at the two Siberian sites (see Section 2.1.2) are extremely scarce. Given this limitation, instead of requiring the satellite and AERONET data to match in time and space, we only required that they represent Siberian BB aerosols of similar ages. Our underlying assumption was that a sufficiently large ensemble of the AERONET data from the two Siberian sites situated in Western and Eastern Siberia was representative, on average, of the typical intensive optical properties of Siberian BB aerosol of the corresponding age, regardless of the fact that the observations were made at particular points in space and time. Accordingly, we expected that if both satellite and AERONET observations in Siberia were unbiased, the Siberian BB aerosol intensive optical properties derived from the satellite and AERONET in the same region should be quantitatively similar, even if they did not match in time and space. Consistent with this consideration and availability of the model estimates of the BB aerosol photochemical age, the AERONET data used in the analysis presented in this section included the observations for both 2012 and 2016.
Note that a similar approach to the validation of satellite observations of BB aerosol was previously used in Konovalov et al. [50], except that the satellite and AERONET data were compared there even regardless of the photochemical age of BB aerosol: it was found that the AER 388/550 derived from the OMI and MODIS measurements over Siberia was, on average, close to that derived from ground-based measurements at the same AERONET sites as in this study (but only in 2012), and this finding was regarded as evidence that the satellite-derived values of AER 388/550 for Siberian BB aerosol were sufficiently reliable.
To ensure the BB aerosol age consistency between the satellite and ground-based observations, we went through our data set of the pre-processed satellite observations (Section 2.1.1) and checked whether there was any AERONET observation that matched a given data point from the satellite data set with regard to photochemical age. The observations were considered as matching each other only if the relative difference between values of the BB aerosol photochemical age in the satellite and AERONET data would not exceed 10%. The combined set of data satisfying this condition included 347 satellite observations and 29 data points of the AERONET observations.
A comparison of the satellite and AERONET observations is presented in Figure 5. We characterized the AAOD 388 /AAOD 867 , AAOD 867 /AOD 550, and AAOD 388 /AOD 550 ratios by considering the slopes of the linear fits (through the origin) to the relationships between the corresponding pairs of the aerosol optical depths. A similar approach was used in our previous study [50]. As evidenced by Figure 5, the estimates based on the satellite and AERONET data agreed with each other within the 95% confidence intervals. That is, the probability of a type I error (in our case, the error of assuming that the slopes estimated using the satellite and AERONET data were different) was found to exceed at least 0.05 (and was actually much bigger). Furthermore, there was a close quantitative agreement between the mean values of the AAOD (at both wavelengths) estimated using, on the one hand, the OMI and MISR data and, on the other hand, the AERONET data. We regarded this occasional quantitative agreement for the extensive optical properties (which was not expected given that the satellites and ground-based photometers observed different BB plumes) as an additional indication that both the satellite and AERONET data selected for our comparison were equally representative of the major optical properties of Siberian BB aerosol. The lines show linear fits through the origin, with the 95% confidence intervals indicated by shading. Note that the confidence intervals for the satellite data are narrower than those for AERONET data because the number of points in the satellite dataset (347) was much bigger than in the AERONET dataset (29).
Note, however, that the above comparison did not take into account a possible extrapolation error in the AAOD 388 values derived from the AERONET observations. This error is likely to result in some underestimation of the derived AAOD 388 values because AAE for BB (BrC-containing) aerosol is likely to increase at shorter wavelengths. The magnitude of this underestimation is difficult to quantify but to obtain a tentative idea about it, we considered the measurements of smoke aerosol absorption parameters in the Amazon region [42]. According to these measurements, using the AERONET observations at 440 and 870 nm to estimate AAOD at 388 nm would result, on average, in an underestimation of AAOD 388 of about 15%. We regard this crude estimate as an indication that the systematic bias in the AAOD 388 values derived from AERONET observations of BB burning plumes was unlikely to exceed the range of random uncertainties in the estimates of the AAOD 388 /AAOD 867 ratio from the AERONET data (which was about 35%).
Overall, the comparison presented in this section did not reveal any evidence for major biases in the satellite data used in our analysis, although it should be kept in mind that this comparison did not distinguish between the observations of fresh and aged BB aerosol. Validation of the ability of the satellite data to adequately represent the atmospheric evolution of BB aerosol in the framework of this study was hampered by the scarcity of the AERONET data, especially those representing fresh BB aerosol.

Monte Carlo Algorithm for Estimation of the BrC Absorption Parameters
Our algorithm described in this section was designed to derive estimates of the following three quantities: (1) the relative contribution of BrC to the light absorption by BB aerosol at the 388 nm wavelength (this quantity is denoted below as δBrC), (2) the IRI for the organic component of BB aerosol (k OA ), and (3) the ratio of the BC and OA mass concentrations in BB aerosol particles (the BC/OA ratio) from satellite observations described in Section 2.1.1. Similar to Konovalov et al. [40], δBrC is defined through the mass absorption efficiency as follows: where α a is the mass absorption efficiency of BB aerosol and α a |k OA = 0 is the same property in a hypothetical situation where the OA is non-absorbing (k OA = 0). More formally, the algorithm was aimed at inferring the components of a control vector x, which characterizes the chemical composition and microphysical structure of BB aerosol, from an observation vector z, which characterizes the BB aerosol optical properties. The components of x were the aerosol parameters that included k OA and the BC/OA ratio, along with other multiple characteristics, whereas z included (in this study) only three components representing AAE 388/867 , AER 388/550 , and AER 867/550 .
As we found in our previous study [40], the values of AAE for Siberian BB aerosol are closely related to δBrC (such that, on average, the AAE increases with an increase in δBrC). Furthermore, we argued [40] that the use of complementary data (such as the SSA) characterizing the relationship between the absorption and scattering properties of BB aerosol, along with the AAE, enables constraining k OA and the BC/OA ratio. In this study, we closely followed the ideas and methodological developments of the mentioned previous study [40]. As argued below (see Section 2.2.3), the use of values of AAE 388/867 , along with values of AER 388/550 and AER 867/550 (which are used in this study instead of the SSA), in the framework of the method described in this section allowed us to constrain, to some extent, all three parameters of interest (δBrC, k OA , and the BC/OA ratio).
Note that δBrC is not a component of x but it is estimated according to Equation (4) as a function of k OA , the BC/OA ratio, and other aerosol parameters affecting the mass absorption efficiency. The components of z are normalized to the standard deviations (calculated using the harmonized satellite data described in Section 2.1.1) of the indicated characteristics. An overview of our algorithm is presented in Figure 6.
We looked for estimates of x that were consistent with a given (nth) observed realization of z, namely, z n o , by employing a Monte Carlo method. Note that a Monte Carlo method was previously used by Lu et al. [74] to estimate k OA for BB aerosol using various field and laboratory observations. A key difference in our method from the method used by Lu et al. [74] is that we used a k-nearest neighbor algorithm, explained below, while Lu et al. [74] looked for optimal estimates. Our estimation method involved three main steps: (1) generation of a large ensemble of samples of x (referred to below as "a priori" estimates of x and denoted as x a ); (2) computation of the modeled counterpart, namely, z m , of the observation vector for each sample of x a ; and (3) obtaining of the "a posteriori" estimates of x (x p ). The first two steps are explained in detail below in Section 2.2.2. Briefly, the first step involved the random sampling of several parameters representing the refractive indexes, relative mass fractions, and microphysical states of major components of BB aerosol particles from the corresponding a priori probability distributions. These samples were then used as input data for the Mie theory calculations of the BB aerosol optical properties. In the third step, the a posteriori estimates of x, namely, x n p , were computed as a linear superposition of the fixed number (k) of the selected samples of x a , namely, x n ai k i=1 , following a k-nearest neighbor algorithm [75]: In our case, the algorithm involved finding k nearest neighbors of z n o , namely, z n mi k i=1 , among all the samples of z m in terms of the Euclidian distance d j between z n o and the sample z n mi : where T denotes the transpose. The same k samples that corresponded to the nearest neighbors of z o could yield an estimate of any scalar function F(x), namely, F n p : In particular, we applied Equation (7) to estimate δBrC according to Equation (4). The results reported in this paper were obtained using a k value equal to 80, i.e., the 80 nearest neighbors were chosen, and their results averaged. Our choice of this particular number is explained below in Section 2.2.3. The total Monte Carlo ensemble of the x a and z m values included 10 6 samples.
The estimation procedure described by Equations (5)- (7) implies that the selected samples z n mi k i=1 were sufficiently representative of the observation. Ideally, given a sufficiently large ensemble of random samples, the distance d i should approach zero for any i. However, because the observations can be affected by various errors, the distances can remain big, meaning that the available simulations are not representative of the observation and cannot provide reliable estimates of the aerosol parameters. To identify such situations in which the k-nearest neighbor algorithm can fail, we introduced a special criterion for the distances d i . Specifically, we required the distances d i to be less than one: We recall that the distances are defined in the space of the normalized variables. Observation data points that do not satisfy this criterion were excluded from further analysis. Hence, the condition given by Equation (8) was used as a filter that excluded potentially unreliable observation data points and the corresponding estimates of the aerosol properties.
Note that the same properties of Siberian BB aerosol (δBrC, k OA , and the BC/OA ratio), which were estimated in this study using satellite observations, were inferred in Konovalov et al. [40] from AERONET data. The methods used in these two studies were similar but not identical. A key difference is that the analysis performed in our previous study involved an explicit formulation of the conditional probability distribution function for x in accordance with the Bayesian formalism. The application of the same approach to the satellite observations is hardly feasible because the conditional probability distribution function depends on the observation errors that are not practically known for the absorption characteristics derived from the OMI and MISR observations. A drawback of the k-nearest neighbor algorithm employed in this study is that it cannot provide confidence intervals for each inferred estimate. However, we are interested here only in the average (typical) absorption properties of BB aerosol as a function of the photochemical age and obtain the required confidence intervals for such properties by considering an ensemble of the estimates for different observations.
The estimates of the BB aerosol characteristics, which were inferred from the satellite observations according to Equations (5)-(8), were processed using a simple algorithm of the nonlinear trend analysis to examine the dependences of these characteristics on the BB aerosol photochemical age. The algorithm (which is described more in detail in Konovalov et al. [43]) involves a nonlinear approximation of the dependence of an available set of M uncertain estimates {p i } M i=1 of any parameter p on the BB aerosol photochemical age τ a with a sum of N sigmoids as follows: where the w are the weight coefficients and ε i is the approximation error (which is assumed to be a random variable with a zero mean). Following Konovalov et al. [43] and for definiteness, we set N = 2. The weight coefficients were optimized using the Nelder-Mead simplex algorithm [76] by minimizing the variance for the errors ε, and the confidence intervals were obtained using the bootstrapping technique [77]. The spread of the "bootstrapped" estimates was used to calculate the 95% confidence intervals. Note again that the analysis of transformations in the properties of BB aerosol due to its aging was not done previously in Konovalov et al. [41]. Such analysis was not feasible, particularly because too few AERONET observations in Siberia are representative of relatively fresh BB aerosol. In contrast to the AERONET data, the OMI observations of BB aerosol (which were used in this study as a primary source of information on the BrC absorptions) were abundantly taken over relatively fresh BB plumes [50].

Mie Theory Computations
In this study, an ensemble of probable values of the optical properties (z m ) of BB aerosol was obtained in almost the same way as in Konovalov et al. [40], where the corresponding simulations are explained in detail. Hence, our computations are described below only briefly.
Following Lu et al. [72] and Wang et al. [18], we assumed that BB aerosol particles were composed of a strongly absorbing core and weakly absorbing shell. More specifically, we assumed that the core consisted of BC, while the shell was mostly composed of organic matter and, to a minor extent, inorganic salts and water. The computations were made using the OPTSIM software [78], which allows for simulating the absorption and scattering properties of an ensemble of multidisperse spherical particles that have homogeneous or core-shell structures. Water content was calculated based on the κ-Köhler theory [79] as a function of relative humidity (RH). The particle components were log-normally distributed among 20 bins representing particle diameters ranging from 10 nm to 10 µm. The parameters of the particle size distribution (the median and the geometric standard deviation) were chosen to be representative of the accumulation mode.
To generate an ensemble of possible realizations of optical properties and microphysical states of BB aerosol (see Section 2.2.1), we randomly varied several aerosol parameters (Table 1), as well as the relative humidity (RH). The variable parameters used were the median and the standard deviation of the volume size distributions for the core and shell, k OA at 550 nm, the BC/OA ratio, the hygroscopicity parameter κ, and the mass fraction of the inorganic material (represented by ammonium sulfate). Assuming k OA to be a power function of the wavelength, we also randomly perturbed the exponent of the wavelength dependence (w). The variations in most of these parameters were specified using a truncated Gaussian distribution, except for k OA and κ, which were assumed to follow uniform distributions. Consistent with the condition that relative humidity (RH) corresponding to the selected satellite observations should not exceed 70% (see Section 2.1.1), the RH was kept to be less than 70% in our computations. This limitation resulted in a rather small water content in the particles (typically, less than 15%). Further details about the sampling of the variable parameters can be found in Konovalov et al. [40].  3 Geometric mean diameter of the BB aerosol particles. 4 Standard deviation of the particle volume size distribution. 5 Wavelength dependence of k for organic aerosol. 6 Hygroscopicity parameter. 7 The most probable value (estimated as a function of the BC/OA ratio according to Lu et al. [74]) of the wavelength dependence.
In the framework of our Monte Carlo algorithm (Section 2.2.1), the variable parameters mentioned above constituted the components of the control vector x. The parameter values were randomly sampled 10 6 times, and OPTSIM was run independently with each sample. The results of the computations (that is, values of the components of z m ) and the a priori samples of x (x a ) constituted a look-up table that provided a discrete quantification of the complex relationships between the vector of parameters x and the vector of observations z.
Note that because of the lack of knowledge about the optical and physical properties of the coarse particles of typical BB aerosol (as discussed in Konovalov et al. [40]), the coarse mode was not taken into account in our calculations as a distinct mode. The possible effects of the coarse mode on estimates of the optical properties considered here were tentatively examined by Konovalov et al. [40] and were found to result in larger mean values of the estimates of δBrC and k OA and a smaller mean value of the estimate of the BC/OA ratio. The differences between the estimates of these characteristics for the cases with and without accounting for the coarse mode were found to be small (a few percent) when related to their a priori values but, because the organic component of Siberian BB aerosol turned out to be weakly absorbing, were found to be more considerable (tenths of a percent) when related to the a posteriori estimates. Hence, we considered our estimates of δBrC and k OA obtained in this study as the lower limits for the most probable values of these parameters.

Optimization and Validation of the Monte Carlo Estimation Algorithm
The estimates defined by Equations (5) and (7) are dependent on the value of k (the number of the nearest neighbors). Choosing a too small k can be expected to result in large uncertainties in x p , which would reflect the possibility for the same values of z to correspond to different combinations of the components of x. On the other hand, with k approaching the total number of samples, the estimates of x would approach the mean of the ensemble of the samples of x and would become almost independent of the observations. Hence, there should be some optimal value of k, which would depend on the complexity and ambiguity of the relationships between x and z, as well as on the observation errors.
To optimize k, we applied our algorithm to a synthetic data set. Again following our previous study [40], we obtained such a dataset using the AERONET observations. In this study, we used the observations and the look-up table described in Sections 2.1.2 and 2.2.2, respectively. For each observation z o , we found a unique vector z m (along with a corresponding vector x a ) that was closest to z o (in terms of the distance given by Equation (6)). In other words, we first used our algorithm with k equal to 1. The selected samples of z m (satisfying the condition given by Equation (8)), along with the corresponding values of x a constituted the synthetic data set in which the values of x a were considered as surrogates of true (unknown) values of x and were excluded from the main set of the random samples. Next, we applied our algorithm to the synthetic data set iteratively, each time increasing the number k by one and calculating the root-mean-square error (RMSE) of the x p . The iterations were repeated until RMSE reached a minimum.
Note that this procedure does not account for the observational errors (which are not known) but rather elucidates only the uncertainties in x due to the complexity of the relationships between the multiple aerosol parameters and the three observed characteristics. Hence, the RMSE estimated using the surrogate data should be regarded as the lower limit for uncertainties in x p derived from the satellite data. Furthermore, the described procedure would likely result in an underestimation of the "optimal" value of k. However, this underestimation was not expected to entail any major issues for the results of this study. Indeed, since the uncertainties in the satellite data were likely to be dominated by random errors (as indicated by the analysis presented above in Section 2.1.4), the uncertainties in x p were also predominantly random and could be manifested in the corresponding confidence intervals for the dependencies given by Equation (9).
The results of the application of our estimation algorithm to the surrogate data (see Figure 7) indicated that based on the "observed" values of AAE 388/867 , AER 388/550 , and AER 867/550 , the algorithm was capable of providing reasonable (albeit not perfect) estimates of δBrC, k OA , and the BC/OA ratio. Specifically, the coefficient of determination was higher or only slightly lower than 0.7 in all three cases, indicating that the obtained estimates could explain the major part of the variability of the inferred properties. Importantly, the estimates of all three characteristics exhibited very small biases (less than 5%).

Results
Our analysis method, described in Section 2.2, was applied to 355 observation data points satisfying the selection criteria explained in Section 2.1.1, as well as the condition given by Equation (8). Figure 8 demonstrates values of the components of the observation vector z o (that is, AAE 388/867 , AER 388/550 , and AER 867/550 ) as a function of the BB aerosol photochemical age (τ a ). It also shows the non-linear approximations for these values according to Equation (9) and the corresponding confidence intervals. The values of AAE 388/867 (Figure 8a) were found to gradually decrease with the photochemical age from~2.5 at τ a ∼ = 7 h to~1.5 at τ a ∼ = 40 h and then to~1.2 for τ a ∼ = 80. Although the satellite data considered here do not provide any information on the optical properties of fresh BB aerosol, the nonlinear approximation suggests that AAE 388/867 is likely to be around 3 when τ a is close to zero. Since larger values of the absorption Ångström exponent are typically associated with stronger BrC absorption (e.g., [18]), this behavior of AAE 388/867 suggested that the contribution of BrC to the overall absorption of light by BB aerosol at the 388 nm wavelength decreased with the BB aerosol age.
Similar to AAE 388/867 , AER 388/550 (Figure 8b) initially exhibited a decreasing behavior, but after about 35 h, the decreasing trend was replaced by an increasing tendency, which was, however, barely significant at the 95% confidence level (although it was significant at lower confidence levels). The downward part of the trend in AER 388/550 , similar to the decreasing behavior of AAE 388/867 , was indicative of a weakening of BrC absorption and was likely associated with an increase in the SSA at 388 nm (SSA 388 ). An increasing part of the trend in AER 388/550 may have been due to volatilization and photochemical destruction of the SOA [43], as well as due to photobleaching of primary and secondary organic compounds [21].
In contrast to AAE 388/867 and AER 388/550 , AER 867/550 (Figure 8c) did not exhibit any significant changes (at least, at the 95% confidence level), although there was weak evidence that it decreased as τ a increased from 15 h to 30 h. Since AAOD 867 was likely not sensitive to the BrC content, the lack of strong changes in AER 867/550 could be regarded as an additional indication that the strong changes in AAE 388/867 and AER 388/550 were primarily due to changes in BrC absorption. At the same time, a weakly pronounced non-monotonous behavior of AER 867/550 (specifically, its decrease after~15 h and increase after~30 h) was consistent with the notion that the atmospheric evolution of BB aerosol was associated with the formation and eventual destruction of weakly absorbing SOA [43].
The main outcomes (in other words, the end products) of our analysis were the nonlinear dependencies of δBrC, k OA , and BC/OA on the BB aerosol photochemical age ( Figure 9). It was found that δBrC was a decreasing function of the BB aerosol photochemical age (Figure 9a): on average, δBrC dropped by about a factor of two (from 0.50 to 0.23) as τ a increased from 7 h to 30 h. Furthermore, since the retrieved tendency tentatively suggested that δBrC for fresh aerosol was considerably larger (more than 0.6) than that at 7 h, the results shown in Figure 9a indicated that the contribution of BrC to light absorption by fresh aerosol originating from Siberian wildfires at 388 nm was predominant and that it eventually decreased by a factor of almost 3. It should be noted, however, that the estimates of δBrC for relatively fresh aerosol (τ a < 12 h) were highly uncertain. Importantly, our analysis provided strong evidence that BrC did not fully disappear in Siberian BB aerosol even after 70 h of the atmospheric daylight exposure (~4 full days during the summertime in Siberia). Figure 9. Estimates of (a) the relative contribution of BrC to the total light absorption by BB aerosol, (b) the imaginary refractive index of the organic fraction of BB aerosol, and (c) the ratio of black carbon and organic matter mass concentrations as a function of the BB aerosol photochemical age. Each panel shows the estimates inferred from satellite observations, nonlinear approximations of these estimates, and the corresponding confidence intervals in terms of the 95th percentile.
According to our estimates shown in Figure 9b, k OA behaved very similarly to δBrC. That is, it initially decreased with the e-folding time scale of about 19 h, but the decreasing behavior did not extend beyond 30 h of the daylight evolution. Furthermore, similar to δBrC, k OA remained distinctly different from zero, even after 70 h of the evolution, again suggesting that the organic fraction of BB aerosol retained some absorption properties during almost the whole lifespan of BB aerosol in Siberia (which was about 5-7 days in summer [80]). The mean value of k OA at the second stage of the evolution was found to be about 4.4 × 10 −3 . At the same time, the BC/OA ratio (Figure 9c) did not manifest any significant changes. It should be noted that the lack of statistically significant changes in this or other characteristics in our analysis did not necessarily mean that such changes did not occur in reality: it is possible that they were not identified simply because of an insufficient amount of the observational data and large uncertainties in them.
Taken together, our estimates of k OA and the BC/OA ratio suggested that the decrease in BrC absorption (Figure 9a) was primarily due to a weakening of the absorption properties of the organic component of BB particles rather than due to changes in the relative mass fraction of the organic component. However, our analysis did not allow us to determine the exact factors underlying the changes in k OA . Such factors may include (i) the destruction of organic chromophores within particles by UV irradiation [21] or (ii) aqueous oxidation [20], as well as (iii) the volatilization of primary (strongly absorbing) organic compounds and their replacement by more weakly absorbing secondary compounds [43]. Since a typical time scale of photo-destruction of the organic chromophores in BB aerosol particles is likely to exceed several tenths of hours [21] and the analyzed observations were made in a relatively dry atmosphere (see Section 2.1.1), we supposed that the third mechanism indicated above was the most probable. That is, the evaporation and subsequent gas-phase photochemical destruction of a semi-volatile absorbing fraction of primary OA (POA) as a result of dilution of BB plumes was likely a primary reason for the decreases in both k OA and δBrC during the initial approximately 30 h of the evolution. In other words, the loss of BrC was mostly determined by the loss of POA, and these processes had almost the same rate. We also supposed that the residual BrC absorption was probably associated with low-volatility (and, therefore, very slowly evaporating) POA and SOA species (with the volatility of about 10 µg m −3 or less), which was consistent with our simulations of the evolution of Siberian BB plumes [43], as well as with in situ observations that show the dominance of low-volatility compounds in aged BB aerosol (e.g., [81]).
It can be useful to note that the trends in δBrC and k OA could be estimated in this study owing to the availability of satellite observations of relatively fresh BB aerosol, along with those of more aged aerosol. Specifically, our satellite dataset included 46 and 11 data points featuring the BB aerosol age of less than 15 and 10 h, respectively, and these points constituted sizeable fractions (13% and 3%) of the whole data set. For comparison, the AERONET data set, which was used in this study for validation purposes (see Section 2.1.4), included only 2 (out of 29) data points (that is, less than 7%) representing BB aerosol with the age of less than 15 h, and did not include any observations of BB aerosol with the age of less than 10 h. The lack of the AERONET data representing fresh BB aerosol was the main obstacle to the successful application of our algorithm to these data.
To examine the robustness of our findings reported above, we performed an additional analysis in which we changed the number of the components of the observation vector. In particular, we repeated our estimation procedure while using only one component-AAE 388/867 -instead of the three components in the "base" case. We found a decreasing behavior of δBrC as in the base case but the estimates of k OA became highly uncertain and did not manifest any statistically significant changes. These findings were consistent with our previous analysis [40] that indicated a strong dependence of δBrC on the absorption Ångström exponent, as well as the need for SSA-related observations to constrain k OA . We also considered the test cases where the observation vector included only two components, namely, AAE 388/867 and either AER 388/550 or AER 867/550 . In those cases, we did not detect any statistically significant differences in the dependencies of δBrC and k OA on the BB aerosol photochemical age relative to the base case; qualitatively, these dependencies remained the same as those shown in Figure 9. Hence, we could conclude that our main findings were sufficiently robust and were unlikely to be significantly affected by possible covariances of uncertainties in the different components of our observation vector.

Discussion
The aim of this study was to investigate the feasibility of using multiplatform satellite observations to study the changes in absorption of light by the organic fraction of BB aerosol as a result of its atmospheric evolution. We believe that the overall results of our investigation, which included not only quantitative insights into the evolution of three quantities characterizing the absorption properties of the organic component of BB aerosol (Section 3) but also the development of a robust technique to infer the BrC absorption characteristics from available multiplatform satellite observations (Section 2.2) and demonstration of the reliability of the combined satellite data (Section 2.1.4), have important implications for atmospheric research. Indeed, presently, the information about the absorption properties of the organic component of BB aerosol and atmospheric evolution of these properties is available mainly from in situ measurements (e.g., [17,22,29]), ground-based remote sensing observations [33][34][35][36][37][38][39][40], and lab studies [19][20][21][25][26][27]. This information, while very insightful, is severely limited in terms of the representativeness of very diverse properties of BB aerosol that can form from different fuels under widely varying environmental conditions. Satellite observations can provide a broad spatial and temporal coverage that can greatly enhance our understanding of the spatio-temporal patterns and atmospheric evolution of the BB OA absorption properties. In this sense, the "multiplatform" approach proposed in this study can complement recent promising developments in partitioning BC and BrC contributions to the light absorption based on single-instrument, such as MISR, measurements [44]. Note that the analysis of MISR retrievals [44] has not provided quantitative estimates of any of the aerosol characteristics discussed in this paper and that the OMI observations involved in our analysis had an important advantage of being performed at 388 nm, where the BrC absorption is typically much stronger than at the shortest wavelength (446 nm) of the MISR measurements. However, although our study demonstrated the feasibility of the proposed multiplatform approach to studying the BrC absorption, further work involving larger datasets from different regions of the world is needed to examine the capabilities and limitations of this approach, as well as possible limitations and biases of the suggested Monte Carlo method for the retrieval of the BrC absorption properties.
Although this study represents just one of the initial steps toward the systematic use of satellite observations to infer BrC absorption by BB aerosol, the results of our analysis are scientifically meaningful and have implications for a better understanding of BrC absorption in such a vast and remote region as Siberia, where the absorption properties of BB aerosol and especially their atmospheric evolution are rather poorly known.
In particular, we identified a decrease in the absorption Ångström exponent (AAE) of BB aerosol due to its atmospheric aging (Figure 8a). Similar behavior of the AAE was previously reported based on field data [82] and analysis of AERONET observations [18] in other regions of the world. Therefore, on the one hand, the derived dependence of the AAE on the BB aerosol photochemical age is not contradictory to the current knowledge about the evolution of the AAE. On the other hand, to the best of our knowledge, this is the first evidence of the decreasing behavior of the AAE of Siberian BB aerosol. Values of the AAE exceeding 1 are typically associated with a considerable contribution of BrC to BB aerosol absorption [18,40]. Therefore, the decreasing trend in the AAE is indicative of a weakening of BrC absorption in aged BB plumes from Siberian fires, although the AAE is not unambiguously related to δBrC or the IRI for the organic fraction of aerosol, with it being dependent also on other aerosol parameters, such as the size distribution and volume fractions of the both absorbing and scattering aerosol components [40,83].
We also identified the downward and upward parts of the trend in AER 388/550 (Figure 8b). As noted above, the downward part of the trend was also indicative of a weakening of BrC absorption due to the atmospheric aging of BB aerosol. Assuming that the extinction Ångström exponent is not changing significantly during the atmospheric evolution of BB aerosol, a decrease in the AER 388/550 is likely associated with an increase in the SSA at 388 nm (SSA 388 ). An increasing trend in SSA 388 was previously deduced from the OMI AAOD 388 and MODIS AOD 550 observations in Konovalov et al. [43] and was interpreted mainly through the formation of weakly absorbing secondary organic aerosol (SOA), which tends to increase AOD 550 . An increasing part of the trend in AER 388/550 may be indicative of volatilization of SOA: an indication of the SOA loss after~35 h of BB aerosol evolution in the daytime was also found in Konovalov et al. [43], although a corresponding decrease in SSA 388 was not identified.
The main end-products of this study were the estimated nonlinear tendencies in the contribution of BrC to Siberian BB aerosol absorption at 388 nm (δBrC) and in the IRI for the organic fraction of BB aerosol (k OA ), as well as evidence that the changes in δBrC and k OA were not accompanied by strong changes in the BC/OA ratio. These end products complement each other by bringing insights into the inter-related aspects of the BB aerosol evolution. More specifically, our most important findings were that the contribution of BrC to Siberian BB aerosol absorption at 388 nm (δBrC) first decreased with the photochemical age of BB aerosol and that it then stabilized at a non-zero level. These findings were coherent with the findings from the previous analysis of the AERONET observations in the Amazon region: specifically, Wang et al. [18] found that the mass absorption coefficient for organic carbon dropped almost twofold after about 1-day evolution of BB aerosol in sunlight, but remained nearly constant afterward. The decrease in BrC absorption as a result of photochemical aging of BB aerosol was coherent with multiple studies based on the analysis field (e.g., [17,82]) and laboratory data (e.g., [19][20][21]84]) but was not a trivial finding given evidence that BrC absorption may also initially increase [20,22], apparently as a result of the formation of strongly absorbing SOA.
The estimated values of δBrC and their changes can be used as constraints on the simulations of BB aerosol absorption with chemistry transport models, e.g., [13,43]. However, the detected decrease in δBrC could, in principle, have been due to various reasons. In particular, it could have been a result of loss of the absorption properties by the organic components of BB aerosol, but it could have also been caused by volatilization of the most strongly absorbing organic components or by the formation of a weakly absorbing SOA fraction. To distinguish between these possibilities, and thus, to ensure a more realistic representation of the evolution of BB aerosol in the models, it would be helpful to complement the observational constraints on the behavior of δBrC with constraints on the evolution of both k OA and the BC/OA ratio. Our findings, useful in this context, were that the changes in δBrC were accompanied by corresponding changes in k OA , while there were no statistically significant changes in the BC/OA ratio (even though there may have been actual changes that fell into the confidence intervals but could not be discerned in our analysis). Specifically, k OA first substantially decreased (with the e-folding time scale of~19 h) but then (after about 30 h) remained at a nearly constant (non-zero) level (~4.4 × 10 −3 ). These findings indicated that transformations of BB aerosol due to volatilization, condensation, and oxidation processes did not play a major role in the evolution of the BrC absorption in Siberian BB plumes, or, at least, that the effects of volatilization and SOA formation on δBrC were likely to compensate each other. Note that the estimates of the BC/OA ratio complemented the estimates of δBrC and k OA by providing a useful constraint on the BB aerosol composition, thereby facilitating the consistent model representation of both the absorption and scattering properties of BB aerosol. Overall, our findings and estimates can help represent BrC absorption in BB plumes from Siberian fires in chemistry transport and climate models, where the absorption properties of aerosols are usually determined using assumed values of the IRI for all individual components of aerosol particles, as well as using assumed mass fractions of different components of aerosol particles.
Finally, it is worth mentioning that the mean values of the quantitative estimates of δBrC, k OA , and the BC/OA ratio-0.27, 5.0 × 10 −3 , and 2.3 × 10 −2 , respectively-were comparable to the mean values of the corresponding estimates derived from the AERONET observations-0.21, 2.3 × 10 −3 , and 1.9 × 10 −2 , respectively [40] (Table 2 therein). Some differences in the mean values of the same characteristics in the two studies are expected to take place not only as a result of random uncertainties in the inferred estimates but also because of likely differences between the BB plumes observed from satellites and at the AERONET sites with respect to the photochemical age. Furthermore, the differences in the mean values of δBrC and k OA could have been due to a probable wavelength dependence of BrC absorption (as the estimates based on the AERONET data were reported for the 440 nm wavelength). Even though the mean values of the estimates of δBrC, k OA, and the BC/OA ratio, which were derived from the satellite data, seemed to be rather meaningful, we do not report them as part of the main results of this study since they were not necessarily fully representative of the average values of the addressed characteristics in Siberia. Furthermore, these mean values could have been affected by various systematic uncertainties. The possible sources of such uncertainties in our analysis were basically the same as those discussed in detail in Konovalov et al. [40]. They include a poorly known contribution of the coarse mode of BB aerosol particles to the AAOD and AOD, a priori assumptions about the wavelength dependence of k OA , and the mixing state of BB aerosol particles. In addition, some uncertainties, which were especially difficult to quantify, may be caused by possible error covariances between the components of the observation vector. These uncertainties need to be examined in dedicated studies and carefully addressed in future applications of the approach suggested in this study. However, they can hardly explain our main findings (most of which are qualitative) concerning the evolution of the BB aerosol characteristics.

Conclusions
In this study, we examined the feasibility of deriving useful information on the atmospheric evolution of brown carbon (BrC) absorption of Siberian biomass burning (BB) aerosol from remote sensing observations performed by three different satellite instruments (OMI, MISR, and MODIS). Specifically, we focused on inferring three characteristics: (i) the relative contribution of BrC (that is, light-absorbing organic compounds) to the overall absorption of light by BB aerosol (δBrC), (ii) the imaginary refractive index of the organic fraction of BB aerosol (k OA ), and (iii) the BC/OA mass concentration ratio in the particles. Both δBrC and k OA were estimated at the 388 nm wavelength. Note that the BC/OA ratio indirectly characterizes the BrC absorption, as δBrC is inversely dependent on this ratio if k OA is kept constant. We used the absorption aerosol optical depths at the 388 nm and 867 nm wavelengths (AAOD 388 and AAOD 867 ) obtained from the standard data products based on the OMI and MISR retrievals, respectively, as well as the extinction aerosol optical depth at the 550 nm wavelength (AOD 550 ), which was available from the MODIS retrievals. We selected the data that characterized BB aerosol released from intense Siberian wildfires in July 2016. The selected and harmonized satellite retrievals were combined to constitute the three components of the observation vector-the absorption Ångström exponent (AAE 388/867 ), the AAOD 388 /AOD 550 ratio (AER 388/550 ), and the AAOD 867 /AOD 550 ratio (AER 867/550 )-for the subsequent analysis. The relationships between the satellite-derived values of AAOD 388 , AAOD 867 , and AOD 550 were found to be in good agreement with similar relationships between the ground-based remote sensing data from AERONET.
Our analysis method was to a large extent adapted from Konovalov et al. [40], where it was applied to AERONET data. The extension to satellite data is crucial, as it allowed for a much larger spatial (over Siberia and European Russia) and temporal (from several hours to several days after fire emissions) coverage of our analysis. The method was designed to relate the input ("observed") data to a set of multiple aerosol parameters through a Monte Carlo algorithm involving the Mie theory calculations for 10 6 scenarios. In this study, the set of parameter values that matched the input data was identified using a robust k-nearest neighbor algorithm, which did not involve any quantitative assumptions about observation errors. The method was tested and optimized using synthetic data based on the AERONET retrievals. As part of our analysis method, changes in the input parameters and the derived estimates were identified as a nonlinear function of the duration of BB aerosol daylight evolution (or, in other words, as a function of the BB aerosol photochemical age).
The results of our analysis indicate that BrC absorption (δBrC) substantially (at least twofold) decreased, while Siberian BB aerosol aged during the first 30 h of the daylight evolution, although it should be noted that we did not have any information about changes in any characteristics considered during the initial 7 h of the smoke evolution and that the changes derived for the period from 7 h to~12 h were highly uncertain. The changes in δBrC mirrored a decrease in AAE 388/867 from~2.2 at 7 h to~1.6 at 30 h, as well as a decrease in AER 388/550 from 0.12 to 0.07 during the same period, whereas AER 867/550 did not exhibit any statistically significant changes with time.
One of our most important findings was that BrC absorption did not disappear entirely as the BB aerosol aged, but was likely to remain considerable (δBrC ≈ 0.2), even when the BB aerosol photochemical age reached 70 h. This finding may have important implications for modeling the radiative effects of Siberian BB aerosol in the Arctic, as it is then more absorbing than if only black carbon was taken into account.
While δBrC may, in principle, be due to changes in both k OA and the BC/OA ratio, we found that the detected changes in δBrC were primarily due to similar changes in k OA : as the BB aerosol aged, k OA decreased (the organic compounds become less absorptive) but the BC/OA ratio (indicative of the relative mass fraction of the organic matter) remained nearly constant. Similar to δBrC, k OA was found to be statistically different from zero even after 70 h of the BB aerosol evolution in daylight. The changes in both δBrC and k OA were found to be reflected in the evolution of the "observed" AAE 388/867 .
Overall, this study demonstrated the feasibility of deriving meaningful quantitative information on the atmospheric evolution of BrC absorption in BB plumes from available multiplatform satellite observations of atmospheric aerosol and provided tentative observational constraints on representations of the BrC absorption parameters for Siberian BB aerosol in chemistry transport and climate models. To explore the potential of the proposed "multiplatform" approach to studying the absorption properties of BB aerosol, future studies should involve much larger data sets and address possible systematic uncertainties in the derived properties due to a priori assumptions about the BB aerosol parameters and due to possible biases in satellite data, specifically those depending on the BB aerosol age.  Data Availability Statement: The OMI AAOD, MISR SSA and AOD, and MODIS AOD data are available as parts of the OMAERUV, MIL2ASAE_L2.002, and MOD04/MYD04 data products, respectively, from the NASA EarthData portal at https://search.earthdata.nasa.gov/ (accessed on 15 February 2022). The AERONET data are available through the AERONET download tool at https://aeronet.gsfc.nasa.gov/cgi-bin/webtool_inv_v3 (accessed on 21 October 2021).

Conflicts of Interest:
The authors declare no conflict of interest.