Field measurements were carried out at the Southern University of Science and Technology (22°32′ N, 114°03′ E) in Shenzhen, China, from 21 November to 10 December 2021. Situated in the urban core, the sampling site provides a representative profile of aerosol characteristics typical of densely populated regions within the Pearl River Delta (PRD).

Ambient aerosols were drawn through a diffusion drying prior to instrument sampling to minimize the influence of ambient humidity. Refractory black carbon (rBC) was characterized at the single-particle level using a SP2 (Droplet Measurement Technologies, Boulder, CO, USA). To obtain complementary chemical composition information for individual particles, a SPAMS (Hexin Analytical Instrument Co., Ltd., Guangdong, China) was deployed alongside the SP2. In parallel, a photoacoustic extinctiometer (PAX; Droplet Measurement Technologies, USA), operating at a wavelength of 532 nm, was used to quantify the bulk optical properties of the collected aerosols. To maintain data accuracy and minimize potential instrumental drift during continuous sampling, routine calibrations and performance checks were conducted throughout the campaign.

Simultaneous measurements of gaseous pollutants (O₃, NO, NO₂, NO_x, and CO) and particulate matter (PM_2.5) were conducted at the site, together with meteorological parameters including temperature, relative humidity, wind speed, and wind direction. The mixing layer height (MLH) was obtained from the ERA5 hourly reanalysis dataset provided by the European Centre for Medium-Range Weather Forecasts (ECMWF). ERA5 data were retrieved at a horizontal resolution of 0.25°×0.25°. The time series of MLH for the sampling site was extracted using the nearest-neighbor interpolation method. To assess the influence of regional air mass transport, 72 h backward trajectory analyses were performed using the HYSPLIT-4 model (Hybrid Single-Particle Lagrangian Integrated Trajectory), developed by the Air Resources Laboratory (ARL) of the U.S. National Oceanic and Atmospheric Administration (NOAA), with arrival heights set at 100 m a.g.l. (meters above ground level). A detailed depiction of the site location, trajectory clusters, and meteorological frequency distributions is illustrated in Fig. S1.

In this study, the mass concentration and mixing state of BC were determined using a SP2. This instrument quantifies BC mass via the laser-induced incandescence (LII) technique. For clarity, the BC measured by the SP2 is hereafter referred to as refractory BC (rBC).

The operational principles of the SP2 have been described in detail elsewhere . Calibration of the incandescence signal was performed using Aquadag^® black carbon particles (Aqueous Deflocculated Acheson Graphite, manufactured by Acheson Inc., USA), which yield a stronger incandescent response than ambient BC for the same mass. To correct for this difference, a correction factor of 0.75 was applied . The mass-equivalent diameter of the rBC core (Dc) was calculated from the measured rBC mass, assuming a particle density of 1.8 g cm⁻³ . Scattering signal calibration was conducted using polystyrene latex spheres (Nanosphere Size Standards, Duke Scientific Corp., Palo Alto, CA, USA) with diameters ranging from 50 to 600 nm. The distorted scattering signals generated by individual rBC particles were reconstructed using the leading-edge-only (LEO) fitting technique . The total particle diameter (Dp) was then derived by applying the measured Dc and reconstructed scattering signal to a core-shell Mie model, assuming a complex refractive index of m=2.26+1.26i for the rBC core and m=1.50+0i for the coating material at 1064 nm. Assuming the coating density of 1.5 g cm⁻³ and the core density of 1.8 g cm⁻³ , the coating-to-core mass ratio (MR) was estimated over a defined time window using Eq. (): 1MR=ρcoating∑i=1N(Dp,i3-Dc,i3)ρcore∑i=1NDc,i3 where Dp,i and Dc,i represent the total particle diameter and the core diameter of an individual rBC particle, respectively. Notably, the MR presented here is not derived from direct measurements, but rather estimated based on SP2-derived parameters. As such, the estimated MR serves solely as an indicative proxy for assessing the mixing state of the bulk rBC population, rather than being used to calculate bulk optical properties.

To characterize the mass size distribution of rBC cores, the rBC mass detected by the SP2 was first allocated into discrete size bins spanning a diameter range of 70 to 500 nm. After calculating the mass concentration within each bin, the resulting size distribution was fitted with a lognormal function: 2dMdlog⁡Dc=Aexp⁡-(log⁡Dc-log⁡D0)22log⁡2σg where A represents the peak mass concentration of the fitted distribution, D0 denotes the mass median diameter (MMD), and σg is the geometric standard deviation (GSD) derived for each defined time window.

To investigate the chemical composition of BC-containing particles, a SPAMS was operated alongside the SP2 during the observation period. The operating principles of SPAMS have been described elsewhere . Over the entire sampling period, a total of 3 159 731 single-particle mass spectra were collected. Elemental carbon (EC) ion fragments were used as key markers of BC-containing particles . Particles exhibiting EC ion signals with a relative peak area (RPA) greater than 0.1 were classified as BC-containing particles. Based on this threshold, 501 559 such particles were identified, accounting for approximately 15.9 % of the total particles detected by SPAMS. To further categorize the BC-containing particles, an adaptive resonance theory-based clustering algorithm (ART-2a) was employed . This clustering analysis ultimately resulted in three distinct groups: EC, ECOC-NO_x, and ECOC-SO_x. The parameters used for ART-2a clustering were as follows: a learning rate of 0.8, a vigilance factor of 0.05, and 20 iterations. The average mass spectral patterns of each particle type are shown in Fig. S2.

SPAMS exhibits optimal mass spectral detection efficiency for particles with aerodynamic diameters of 400–600 nm, the BC particles measured within this range are generally in an aged state. Temporal analysis shows that the number concentration of the EC cluster exhibits the strongest correlation with rBC number concentrations (R2=0.75, Fig. S2). Given that EC particles typically preserve the initial chemical signatures of primary BC emissions, positive matrix factorization (PMF) analysis was applied to the relative peak areas (RPAs) of characteristic ion fragments within the EC cluster. By resolving the relative contributions of these characteristic chemical peaks, we aimed to trace the specific aging mechanisms of freshly emitted BC. Detailed procedures and parameter settings related to the PMF analysis are provided in the Supplement.

A PAX (Droplet Measurement Technologies, USA) operating at a wavelength of 532 nm, was employed for in-situ measurements of aerosol optical properties. The scattering and absorption channels were calibrated using polystyrene latex (PSL) spheres and fullerene soot, respectively. The absorption coefficients measured by the PAX showed strong agreement with both the rBC mass concentrations derived from the SP2 (R2=0.95) and the absorption coefficients calculated using a core-shell Mie model (R2=0.94) (Fig. S4), indicating the reliability of the optical measurements. To quantify the absorption property of rBC, the mass absorption cross section (MAC) was calculated by combining the light absorption coefficient measured by the PAX with the rBC mass concentration derived from SP2, as follows: 3MAC=babs[rBC] where babs is the aerosol absorption coefficient measured by PAX and [rBC] is the rBC mass concentration from SP2.

Given that some BrC is known to absorb at 532 nm and measurements at longer wavelengths were unavailable, we adopted the method proposed by to isolate the BrC contribution from the total absorption. First, the relative coating thickness of individual rBC particles was inferred based on the time lag between the peaks of the incandescence and scattering signals recorded by the SP2. Particles exhibiting a lag time greater than 2 μs were classified as thickly coated (Fig. S5a). Assuming that freshly emitted BC corresponds to a thinly coated population under ambient conditions, the intercept of the MAC obtained by extrapolating to a zero fraction of thickly coated particles was interpreted as the MAC of fresh BC , yielding a value of 7.13 m² g⁻¹, with a 95 % confidence interval of 6.57–7.78 m² g⁻¹ (Fig. S5b). The apparent absorption enhancement for each hourly interval was calculated as 4Eabs=MACMACfresh, A campaign-specific relationship between Eabs and MR was subsequently established from the observations (Fig. a). This relationship was used to estimate the apparent absorption enhancement corresponding to the measured MR at each time interval. The BrC contribution was estimated as the residual between the measured total absorption and the estimated BC absorption: 5babs,BrC=babs,obs-babs,BC=babs,obs-MACBC,ref⋅Eabs(MR)⋅[rBC] where babs,BrC, babs,obs, and babs,BC are the absorption by BrC, the observed absorption, and the estimated absorption for BC particles, respectively. The resulting babs,BrC estimates were subsequently averaged by hour of day to obtain a mean diurnal profile (Fig. b). When calculating the MAC of rBC, the mean babs,BrC over the predefined daytime period was subtracted from the observed absorption during the corresponding daytime hours, after which the MAC was recalculated using the adjusted absorption coefficient. During daytime, the estimated BrC contribution accounted for 2 %–9 % of the measured absorption, corresponding to a MAC correction of 0.25–1.15 m² g⁻¹. In contrast, because high rBC concentrations dominate the optical absorption at night and yield physically meaningless negative babs,BrC values (Fig. 1b), the BrC contribution was not further corrected during nighttime. These negative values indicate that the nighttime BrC contribution could not be robustly resolved and was likely minor relative to BC absorption. This treatment does not imply the absence of nighttime BrC absorption, and the daytime correction should be regarded as a model-dependent adjustment. All optical and single-particle measurements were averaged to the same hourly time resolution before regression analysis and diurnal averaging.

Figure  illustrates temporal variations in measured meteorological parameters (temperature, humidity, and wind speed) together with concentrations of PM_2.5, O₃, NO_x, CO, and rBC from 21 November to 9 December 2021. Based on 72 h backward trajectory analyses, air masses were classified into three distinct clusters (Fig. S1a). Cluster 1 (C1), originating from the eastern coastal regions, dominanted from 21 November through 7–9 December. Cluster 3 (C3) represented long-range transport from the northwestern interior, influencing the site during 22–25 November and 1–3 December. Periods influenced by C3 were characterized by decreasing temperatures and relative humidity, typical signatures of wintertime cold waves, accompanied by stronger wind speeds and lower rBC concentrations (Figs. b and S1b). Despite the lower rBC levels, elevated MR during these intervals suggests an increased proportion of aged BC particles. In contrast, Cluster 2 (C2), which accounted for 41.23 % of the total sampling time, exhibited minimal spatial movement, thereby limiting the influence of long-range transport. Within the C2-dominated conditions, a continuous 7 d interval was identified (green box in Fig. ) during which O₃, NO_x, CO and rBC all displayed consistent diurnal patterns. This period was therefore selected as a representative case dominated by the local emissions for detailed analysis of BC aging processes.

During this episode, ambient temperatures ranged from 11 to 26 °C (mean: 18 °C), and relative humidity varied between 26 % and 80 % (mean: 55 %) (Fig. a). The mean PM_2.5 concentration was 31.3 µg m⁻³, indicating moderate pollution levels (Fig. b). The average rBC concentration was 1.2 µg m⁻³, contributing approximately 3.8 % to the PM_2.5 mass. Furthermore, rBC showed strong correlations with NO_x (R2=0.71, p<0.001) and CO (R2=0.51, p<0.001), both recognized tracers of vehicular emissions in urban environments . Besides, the sampling site is located in a densely trafficked urban area without major industrial sources nearby. The site characteristics and the observed pollutant correlations strongly suggest that local traffic emissions were the primary source of rBC during this period. Characterized by consistent diurnal variations in both meteorology and pollutants, relatively weak long-range transport, and traffic-dominated BC emissions, this episode provided a useful observational window for examining local diurnal BC aging processes. To evaluate the representativeness of the selected period, we further compared its diurnal variations with the full-campaign average (Figs.  and S6). The selected period showed generally consistent timing and direction of the key diurnal variations observed over the full campaign, whereas the full-period data exhibited greater variability, likely reflecting the influence of air masses associated with different trajectory clusters. These similarities suggest that the key diurnal patterns were not unique to the selected period, while the local-emissions subset provided a less heterogeneous observational context for examining daytime and nighttime BC evolution. While intense morning emissions were expected, the rBC concentration (Fig. c) exhibited a distinct bimodal pattern, with peaks occurring only during the evening rush hours (18:00–20:00 LT) and around midnight (00:00 LT). The absence of a morning peak is likely attributed to rapid dispersion associated with the rising mixing layer height (MLH) and wind speed (Fig. a, d). The midnight rBC peak coincided with significant increases in NO and CO concentrations (Fig. b, h). Given the lack of nearby industrial sources, this synchronized surge across NO and rBC likely points to emissions from heavy-duty diesel vehicles.

The size distribution of rBC cores further supports an urban traffic-dominated source. The MMD ranged from 155 to 170 nm (Fig. f) with relatively narrow and stable GSDs (1.55–1.58, Fig. S7), typical of urban vehicular sources . A notable observation is the continuous increase of MMD in the nighttime from 155 to 170 nm. Normalized size distributions (Fig. S7b) indicate that the slight diurnal fluctuations in MMD were driven by the intermittent inputs of smaller, freshly emitted BC particles, rather than by regional transport. Given the minimal influence of external transport, variations in MR provide evidence for BC aging. The MR exhibited a continuous increase from 08:00 to 14:00 LT (Fig. i), correlated with O_x (R2=0.51, Fig. S8), indicating that photochemical processing dominated daytime BC aging . In addition, MR also increased during nighttime, suggesting nocturnal aging of BC. Similar behavior has been reported in previous studies under conditions conducive to heterogeneous reactions .

Distinct aging pathways during day and night are further corroborated by the diurnal distribution of rBC coating thickness (Fig. ). During nighttime, BC particles exhibited gradual growth, with the coating thickness of the dominant population increasing by ∼ 10 nm. In contrast, daytime BC particles generally exhibited thinner coatings, primarily reflecting the continuous input of freshly emitted BC from intense daytime traffic emissions. Under active photochemical conditions, however, a secondary mode emerged at a coating thickness of ∼ 60 nm. The thinly coated mode was likely maintained by continuous local emissions near the surface, whereas the thickly coated secondary mode may have arisen from several concurrent processes, including rapid photochemical aging of a subset of BC particles and possible entrainment or downward mixing of more aged particles from aloft during planetary boundary layer (PBL) development. In addition, PBL expansion may preferentially disperse freshly emitted, thinly coated BC near the surface, thereby reducing its relative contribution to the normalized coating-thickness distribution. The observed daytime changes were consistent with the combined influence of boundary-layer entrainment and rapid photochemical aging, but the contribution of each process remains uncertain. The emergence of this bimodal distribution nevertheless highlights the pronounced heterogeneity of BC mixing state during daytime, in contrast to the more uniform coating characteristics observed at night.

The diurnal variations in the mixing state demonstrate that BC undergoes distinct daytime and nighttime aging processes. To further elucidate the chemical mechanisms driving these dynamic transformations, we examined the diurnal evolution of the chemical composition of BC-containing particles (Fig. ). The diurnal variation of OC closely tracked that of the BC mass concentration (Fig. c). In contrast, oxygenated organic carbon (OOC) and sulfate (SO42-) exhibited drastic increases during the daytime, indicative of secondary formation. Nitrate (NO3-) displayed a bimodal pattern, with one peak occurring at night and a second peak emerging during the early stages of daytime photochemistry. Ammonium (NH4+) increased concurrently during the day, aligning with the rising trends of both sulfate and nitrate. In comparison to non-BC particles (Fig. b), which share similar diurnal trends for other secondary species, BC-containing particles exhibit significantly greater RPA variability exclusively for SO42-. The stronger sulfate signals observed in BC-containing particles are consistent with preferential sulfate accumulation on these particles, although the underlying mechanisms cannot be uniquely determined from the SPAMS measurements alone. BC-involved photochemistry has been proposed as one possible mechanism for enhanced sulfate accumulation in previous studies , but its contribution to the present observations requires further investigation.

Coagulation may also contribute to the diurnal evolution of the BC mixing state by transferring pre-existing secondary material from non-BC particles to the BC-containing population. This process may be more favorable under higher particle number concentrations and weaker atmospheric dispersion, conditions that occurred more frequently during nighttime in this study (Fig. c). The observed nighttime increase in MMD (Fig. f) is consistent with particle growth but cannot be uniquely attributed to coagulation, because condensation, heterogeneous uptake, and changes in particle sources may also contribute. During daytime, PBL development and stronger atmospheric dispersion may reduce particle concentrations and thereby decrease the likelihood of coagulation. Nevertheless, quantitative coagulation timescales cannot be constrained with the available measurements. Therefore, the observed increases in rBC coating thickness and secondary species associated with BC-containing particles cannot be unambiguously attributed to direct condensation or heterogeneous uptake, as coagulation-driven mixing may also contribute.

To further resolve the aging processes of freshly emitted BC, PMF analysis was applied to particles within the EC cluster. This approach decomposes complex mass spectral data into distinct chemical factors and their temporal contributions . The resulting factor profiles and diurnal patterns are shown in Fig. . A substantial increase in the Nitrate factor was observed during nighttime, indicating enhanced nitrate accumulation on BC-containing particles. Several processes may have contributed to this nighttime enhancement. The elevated nighttime NO₂ levels (Fig. e), together with the presence of O₃ (Fig. g), may provide conditions favorable for NO3/N2O5 chemistry, and heterogeneous hydrolysis of N₂O₅ on particle surfaces could therefore contribute to nitrate formation. However, this pathway cannot be directly constrained in the absence of NO₃ and N₂O₅ measurements. Thermodynamic gas–particle partitioning may also contribute to the observed nitrate enhancement. In particular, lower nighttime temperatures and higher RH may favor the partitioning of semi-volatile nitrate into the particle phase, while the presence of sufficient NH₃ and aerosol liquid water could further promote particulate NH₄NO₃ formation . Therefore, the observed nighttime nitrate accumulation likely reflects the combined influence of heterogeneous nitrogen chemistry and thermodynamic partitioning, although their relative contributions cannot be quantified with the available measurements.

During the early stage of daytime photochemical aging, the contribution of the nitrate factor began increasing earlier than that of the Sulfate–OOC factor. This temporal offset may be explained by two possible mechanisms: (1) SO42- precursors may not have reached sufficient concentrations at that time ; and (2) gaseous HNO₃ formed via photochemical reactions may have been directly taken up by BC particles through heterogeneous processes . However, the available observations do not permit a definitive determination of the dominant pathway. The Sulfate–OOC factor reached its maximum between 11:00 and 14:00 LT, coinciding with a pronounced increase in the ammonium factor. This concurrence suggests enhanced formation of (NH₄)₂SO₄ and other secondary photochemical products, leading to accelerated BC aging through condensational growth.

In alignment with the bulk BC-containing particle population, the primary EC cluster exhibits similar trends in secondary species accumulation. This consistency across different analytical approaches explicitly highlights the highly dynamic nature of BC aging, pointing to distinct chemical evolutionary pathways operating under daytime and nighttime atmospheric conditions.

The diurnal aging of BC observed during the sampling period was accompanied by corresponding changes in its optical properties. As shown in Fig. a, the MAC reached a minimum of approximately 8.1 m² g⁻¹ at 22:00 LT. Thereafter, it increased by approximately 0.8 m² g⁻¹ overnight, reaching 8.9 m² g⁻¹ by 06:00 LT. A comparable enhancement occurred during the daytime photochemically active period, when MAC increased by approximately 1.0 m² g⁻¹, from 8.9 to 9.9 m² g⁻¹. The observed MAC range of 8.1–9.9 m² g⁻¹ falls within the range previously reported for urban BC. At wavelength of 532 nm, reported MAC values of 5.0–8.5 m² g⁻¹ in summertime Shenzhen. A higher campaign-average MAC of 14.6±5.6 m² g⁻¹ was reported for urban Xi'an using an SP2–PAX system . In urban Beijing, further showed that the modeled MAC at 550 nm varied substantially with BC mixing state, ranging from approximately 5.3–7.3 m² g⁻¹ for largely uncoated or thinly coated BC to 11.2–12.4 m² g⁻¹ for moderately or thickly coated BC.

Although the daytime increase occurred over a shorter period than the nighttime increase, the comparable net enhancements in MAC during both periods indicate that multiple atmospheric processes jointly influenced the evolution of BC optical properties. BC particles in the early stages of aging exhibit negligible light absorption enhancement, however, as aging proceeds and a core–shell morphology develops, coating-induced lensing can substantially enhance their light absorption . Following , MR > 3 was used as a qualitative indicator to identify highly coated BC-containing particles whose optical behavior is consistent with a core–shell morphology. The mass fraction of rBC associated with these particles was calculated as: 6FC-S-like=∑i:MRi>3mrBC,i∑imrBC,i, where mrBC,i denotes the rBC mass of the ith particle, and FC-S-like represents the fraction of total rBC mass associated with particles classified as core–shell-like. The MR threshold here was used only for particle classification, while Fthickly represents the quantitatively calculated fraction of total rBC mass carried by particles within this highly coated regime. As Fig. b shows, the FC-S-like increased from 0.32 to 0.35 during the night, decreased to 0.32 following the morning rush hour, and subsequently rose to 0.40 after daytime photochemical aging. Notably, the magnitude of this diurnal change is substantially larger than the observed day-night variation in bulk MAC.

In the ambient atmosphere, the bulk optical properties of BC are influenced not only by aging but also by emissions, regional transport, and atmospheric removal processes. Given the minimal influence of transport during the selected period, the observed diurnal variation of the MAC is governed by the interplay of fresh emissions, aging, and atmospheric removal. These competing effects can be visualized in the two-dimensional relationship between the bulk absorption coefficient and rBC mass (Fig. ). The fresh-emission trajectory (black arrow), in which absorption and rBC mass increase proportionally, was defined using an empirical fresh-BC-equivalent MAC baseline (MACfresh). This baseline was obtained by extrapolating the empirical relationship between MAC and the fraction of thickly coated BC particles to a zero thickly coated fraction. The resulting MACfresh was 7.13 m² g⁻¹, with a 95 % confidence interval of 6.57–7.78 m² g⁻¹. This value is comparable in magnitude to previously reported MAC values for fresh or uncoated BC, including 6.5±0.5 m² g⁻¹ at 532 nm in urban Shenzhen , 7.5±1.2 m² g⁻¹ at 550 nm for uncoated carbonaceous particles , and 8.0±0.7 m² g⁻¹ at 550 nm for freshly emitted BC . Although these values were obtained under different measurement conditions and definitions of fresh BC, their comparable magnitudes provide useful context for the empirical baseline derived here. Aging enhances absorption without adding rBC mass, resulting in a vertical shift (red arrow). In contrast, removal processes reduce rBC mass and preferentially remove aged particles, yielding a trajectory with a slope larger than MAC_fresh (blue arrow).

To isolate absorption changes attributable to the optical variations driven solely by atmospheric aging, linear regressions were applied to the data collected during the respective daytime and nighttime aging periods. This approach estimates the deviation of the measured absorption coefficient from the theoretical baseline, assuming fresh emissions at the same rBC mass concentration. The difference represents the additional absorption coefficient induced by the lensing effect (bEabs), which can be quantified as follows: 7bEabs=babs,obs-MACfresh⋅[rBC] where MACfresh represents the MAC of freshly emitted rBC and [rBC] denotes the rBC core mass concentration. Although the potential interference from fresh emissions cannot be entirely isolated, bEabs serves as a useful proxy for characterizing the optical variations driven by daytime and nighttime aging. As bEabs is dominated by contributions from thickly coated BC particles, we quantified the MAC of thickly coated BC and the temporal rates of MAC change driven by atmospheric aging: 8MACC-S-like=bEabs+MACfresh⋅[rBC]⋅FC-S-like[rBC]⋅FC-S-like9RMAC,C-S-like=MACC-S-like, max-MACC-S-like, minΔt where FC-S-like denotes the fraction of total rBC mass associated with particles classified as core–shell-like, Δt represents the elapsed time between the occurrences of the minimum and maximum MACC-S-like values within the respective daytime and nighttime periods, and RMAC,C-S-like denotes the corresponding apparent MAC enhancement rate. The coating thickness and MR derived from the SP2 LEO-fit analysis should be regarded as model-dependent estimates rather than directly measured quantities. The retrieved coating thickness is sensitive to assumptions regarding the refractive indices of the rBC core and coating, while the calculated MR additionally depends on the assumed core and coating densities. Deviations from the idealized concentric core–shell morphology may introduce further uncertainty. Consequently, these assumptions may affect the absolute values of the retrieved coating thickness, MR, and the classification of particles close to the operational MR=3 threshold. These uncertainties should therefore be considered when interpreting the quantitatively derived coating properties. The parameter sensitivity was further evaluated by varying the MR threshold, coating density, and MACfresh. For each parameter combination, MAC_C-S-like and FC-S-like were recalculated, and the daytime and nighttime rates were derived. As shown in Fig. , the derived rates varied with the MR threshold, coating density, and MACfresh. The nearly vertical contours in the daytime panels indicate a stronger sensitivity to MACfresh, reflecting the larger daytime variation in FC-S-like. In contrast, the nighttime rates were more sensitive to coating density, likely because density-induced changes in MR resulted in stronger reclassification of particles near the prescribed threshold. The sensitivity results are summarized in the boxplots (Fig. d and h), with the daytime and nighttime RMAC,C-S-like of 0.51±0.11 and 0.36±0.05 m² g⁻¹ h⁻¹, respectively.

It should be emphasized that these rates are campaign-specific estimates derived from the temporal changes in MAC observed during this campaign and evaluated over the empirical parameter ranges considered in the sensitivity analysis. Several instrumental limitations introduce additional uncertainty. Residual BrC absorption at 532 nm may affect the absorption attributed to BC. The SP2 cannot directly resolve the actual morphology of ambient BC particles; therefore, the MR-based identification of core-shell-like BC should be regarded as an optically equivalent classification. Nevertheless, previous urban observations provide useful context for the optical evolution identified here. reported an increase in bulk rBC MAC of 0.60 m² g⁻¹ h⁻¹ during a daytime photochemical period in urban Xi'an. found that the BC absorption enhancement factor reached approximately 2.4 after about 5 h in Beijing and 18 h in Houston. In urban Guangzhou, observed enhanced Eabs during the afternoon in the dry season and elevated nighttime Eabs during the wet season. Although these studies used different optical metrics, their results demonstrate that substantial BC optical evolution can occur over timescales of several hours in urban environments. BC optical properties evolve nonlinearly as particle morphology and composition change during atmospheric aging . These observations may also provide useful constraints for evaluating model representations of BC aging and optical evolution. Recent modeling studies have incorporated observation-constrained BC mixing-state and optical properties to improve estimates of BC radiative effects . In this context, the observed diurnal evolution of BC coating state and MAC can provide benchmarks for evaluating whether models reproduce the characteristic timescales and magnitude of BC optical changes under urban atmospheric conditions. In particular, the contrasting daytime and nighttime evolution highlights the potential importance of representing temporally varying aging environments rather than assuming a single constant aging timescale. The observed relationship between BC coating state and optical enhancement may also help evaluate parameterizations linking BC mixing state to MAC.