Recently, the Kuramoto model with transmission delay has been attracting increasing attention, accompanied by the increase in its practical applications. In this paper, we studied the Kuramoto-Sakaguchi-type Fokker-Planck equation of the above model proposed by Lee et al., in 2009. We proved the global-in-time solvability of the equation under some conditions on the initial data and distribution of delay.
Citation: Hirotada Honda. On Kuramoto-Sakaguchi-type Fokker-Planck equation with delay[J]. Networks and Heterogeneous Media, 2024, 19(1): 1-23. doi: 10.3934/nhm.2024001
[1] | Darko Mitrovic . Existence and stability of a multidimensional scalar conservation law with discontinuous flux. Networks and Heterogeneous Media, 2010, 5(1): 163-188. doi: 10.3934/nhm.2010.5.163 |
[2] | Raimund Bürger, Harold Deivi Contreras, Luis Miguel Villada . A Hilliges-Weidlich-type scheme for a one-dimensional scalar conservation law with nonlocal flux. Networks and Heterogeneous Media, 2023, 18(2): 664-693. doi: 10.3934/nhm.2023029 |
[3] | Mauro Garavello, Roberto Natalini, Benedetto Piccoli, Andrea Terracina . Conservation laws with discontinuous flux. Networks and Heterogeneous Media, 2007, 2(1): 159-179. doi: 10.3934/nhm.2007.2.159 |
[4] | . Adimurthi, Siddhartha Mishra, G.D. Veerappa Gowda . Existence and stability of entropy solutions for a conservation law with discontinuous non-convex fluxes. Networks and Heterogeneous Media, 2007, 2(1): 127-157. doi: 10.3934/nhm.2007.2.127 |
[5] | Shyam Sundar Ghoshal . BV regularity near the interface for nonuniform convex discontinuous flux. Networks and Heterogeneous Media, 2016, 11(2): 331-348. doi: 10.3934/nhm.2016.11.331 |
[6] | Clément Cancès . On the effects of discontinuous capillarities for immiscible two-phase flows in porous media made of several rock-types. Networks and Heterogeneous Media, 2010, 5(3): 635-647. doi: 10.3934/nhm.2010.5.635 |
[7] | Jan Friedrich, Sanjibanee Sudha, Samala Rathan . Numerical schemes for a class of nonlocal conservation laws: a general approach. Networks and Heterogeneous Media, 2023, 18(3): 1335-1354. doi: 10.3934/nhm.2023058 |
[8] | Shyam Sundar Ghoshal, John D. Towers, Ganesh Vaidya . BV regularity of the adapted entropy solutions for conservation laws with infinitely many spatial discontinuities. Networks and Heterogeneous Media, 2024, 19(1): 196-213. doi: 10.3934/nhm.2024009 |
[9] | Raimund Bürger, Stefan Diehl, María Carmen Martí . A conservation law with multiply discontinuous flux modelling a flotation column. Networks and Heterogeneous Media, 2018, 13(2): 339-371. doi: 10.3934/nhm.2018015 |
[10] | Boris Andreianov, Kenneth H. Karlsen, Nils H. Risebro . On vanishing viscosity approximation of conservation laws with discontinuous flux. Networks and Heterogeneous Media, 2010, 5(3): 617-633. doi: 10.3934/nhm.2010.5.617 |
Recently, the Kuramoto model with transmission delay has been attracting increasing attention, accompanied by the increase in its practical applications. In this paper, we studied the Kuramoto-Sakaguchi-type Fokker-Planck equation of the above model proposed by Lee et al., in 2009. We proved the global-in-time solvability of the equation under some conditions on the initial data and distribution of delay.
Atmospheric refractory black carbon (rBC) aerosol warms Earth's climate through absorption of solar radiation in the visible and ultraviolet wavelengths [1,2,3], and causes an estimated total global radiative forcing of 0.64 W m−2 [4]. Regional climate forcing of rBC can be much larger (e.g., 20 W m−2 or more over parts of India [5,6]). rBC is emitted into the atmosphere from incomplete combustion of fossil and biofuel, and from biomass burning [7]. Hence, many rBC sources are strongly linked to anthropogenic activity, and are being considered for inclusion in near-term mitigation strategies [7,8,9,10]. The relative importance of rBC to climate forcing could increase in the future according to projected emission scenarios [11].
Freshly emitted fossil fuel rBC is typically hydrophobic at emission, and is relatively free of non-refractory material (i.e., is externally mixed) [12]. Exposure in the atmosphere leads to accumulation of non-refractory species such as sulfate and organic aerosol (OA) on the rBC component [9], affecting rBC microphysics including optical properties and interactions with water vapor [13]; these processes are what we term "microphysical aging" (here, simply "aging"). The climate impacts of rBC, which depend in part on its size distribution, vertical gradient, mixing state, and abundance relative to co-emitted species such a sulfate [7,9,14], are not well constrained, and contribute nearly 0.5 W m−2 uncertainty to the overall understanding of global climate forcing [4].
The rate of aging of rBC is quite variable not only because of the spatial and temporal variability in the abundance and size distribution of rBC particles, but also because of variability in co-emitted chemical species, in the oxidative capacity of the atmosphere, and in meteorological conditions. Here internally mixed materials are generically referred to as "coatings", while the rBC component is called a "core", without implying any specific morphological information. Subramanian et al. [15] found that the age of an air mass does not necessarily correlate with increase in the amount of coatings associated with rBC cores.
Mixing state can impact the direct radiative forcing of aerosol rBC in two opposing ways: (ⅰ) internally mixed chemical species may enhance the absorption of solar radiation by the rBC fraction (the so-called "lensing" enhancement) [2,16,17,18,19,20,21] even as single-scatter albedo also increases [22], and (ⅱ) internally mixed rBC associated with non-absorbing chemical species such as sulfate is expected to be more hydrophilic and to lead to reduced atmospheric lifetime due to increased deposition [23,24,25,26,27].
Based on recent modeling and laboratory studies, the instantaneous forcing of rBC is believed to be ~50% enhanced by internal mixing [28,29,30,31,32]. Estimates from field measurements suggest a wider range of values for absorption enhancements by this mechanism. Schwarz et al. [33] deduced from aircraft measurements of rBC mass and mixing state that the presence of non-absorbing coatings can enhance light absorption by the ambient rBC column over the tropics by 30%. Lack et al. [34] showed that internal mixing of rBC with OA and nitrate enhanced absorption of visible radiation in a biomass burning plume in rural Colorado by 38%. Wang et al. [35] suggested that in a polluted urban atmosphere in China, the addition of mainly organic coatings increased the mass absorption cross section of rBC by as much as 80%. On the other side, Cappa et al., [36] found significantly a lower absorption enhancement of rBC (AErBC) of ~6% in urban-influenced air in California without any correlation to the photochemical age of the air mass. Lan et al. [37] showed that mass absorption enhancements of rBC due to internal mixing was low in an urban plume of a Chinese megacity with a campaign average absorption amplification of 7%.
In addition to uncertainties associated with the absolute impacts of internal mixtures on rBC optics, there is considerable additional uncertainty about the time scales for such aging to occur, and their integrated impact on rBC forcing. Modeling work by Fierce et al. [38] indicated that the aging time scales of rBC can vary from hours to weeks depending on parameters such as the light flux and hygroscopicity of secondary aerosol, and the size of rBC cores undergoing aging. Stier et al. [14] performed a model study that found that coatings could also enhance removal of rBC, countering their impact on optics, and that these two competing mechanisms have different impacts on the rBC burden and optical properties depending on the geographical proximity to source regions. That study predicted that in the southeastern U.S., and most continental areas in close proximity to source regions, anthropogenic sulfur emissions (specifically) would increase the atmospheric aerosol absorption optical depth by 10–30% despite decreasing the atmospheric burden of rBC by 10–20%. Despite the importance of rBC in climate forcing, there is a paucity of systematic investigations of the interactions between rBC and non-refractory chemical species such as sulfate or OA in real atmospheric conditions. Additional modeling and comprehensive field assessments are needed to improve the understanding of the effect of mixing state of rBC on its radiative properties and hence the overall climate forcing.
The SouthEast NEXus (SENEX) field campaign was a large-scale, collaborative project that took place in the Southeastern United States in June and July of 2013 [39]. The goal of SENEX was to investigate the evolution of biogenic and anthropogenic gases and aerosols in this region, and to explore the formation of haze. The US Southeast is characterized by some of the highest levels of sulfate and OA in the nation paired with highly humid conditions [40,41,42,43,44,45,46,47], while the summer-time measurements occurred when photochemical reactions are very important. Numerous state-of-the-science analytical instruments were deployed onboard NOAA WP-3D research aircraft.
Parts of the SENEX flights were designed to sample urban outflows, shale development regions, and various emission point sources—including a focus on the atmospheric impacts of contemporary coal-fired power plants. In developed countries such as the U.S., these do not emit significant amounts of primary ash or rBC particles, but are major sources of carbon dioxide (CO2), sulfur dioxide (SO2), and nitrogen oxides (NOx) [48,49]. The emitted SO2 can rapidly oxidize to sulfuric acid (H2SO4) in the presence of OH radicals and especially hydrogen peroxide in the plume [50]. The resulting H2SO4 has low vapor pressure and can nucleate to form ultrafine sulfate particles (SO42−) or condense onto preexisting particles [51]. Hence, ultrafine particles are generally present in coal-fired power plant plumes in high concentrations [52,53]. Recent modeling work has shown that new particle formation, and hence particle number in plumes, increases with an increase in mixing ratios of both SO2 and NOx (under low-NOx background conditions), and decreases with increases in background aerosol concentrations [54]. The inverse relationship between new particle formation and background aerosol is expected due to the competition between nucleation and condensation of vapors onto pre-existing aerosol.
Here we investigate interactions between rBC-containing aerosol particles and coal-fired power plant emissions using measurements collected during the larger SENEX campaign. Observations of the amount of internally mixed material (coating thickness) on rBC, assumed due to daytime sulfate-driven chemistry, are used to constrain enhanced aging of rBC in coal-fired power plant plumes, and to evaluate increased aerosol absorption optical depth due to rBC mixing with sulfate in the Southeast U.S.
In situ measurements of atmospheric particulate and gas pollutants, meteorology, and radiation were made from the NOAA WP-3D research aircraft during the SENEX study. Here we focus on power-plant emissions observed on five research flights of 17 total carried out of the Smyrna/Rutherford County airport in Smyrna, TN, U.S. from June 10 to July 10 of 2013. The study covered a sampling area spanning 31.4°N–42.1°N latitude and 95.0°W and 75.6°W longitude. Most of sampling was within the planetary boundary layer with some excursions of up to 6.5 km above sea level.
The measurements of ambient rBC, which we describe using the nomenclature suggested in Petzold et al. [55], were made with a NOAA SP2 sampling at a flow rate of 2 cm3 s−1 [56,57]. The SP2 utilizes laser-induced incandescence to quantify the mass of refractory rBC (vaporizes at temperatures > 4000 K) in single particles with rBC mass corresponding to 80–550 nm volume-equivalent diameters for assumed density of 1.8 g cm−3, independent of particle morphology and mixing state [58]. Upon vaporization of rBC, emitted visible and near-visible thermal radiation are detected with two photomultiplier tube detectors (PMTs) with different wavelength ranges of sensitivity. The rBC composition is constrained with a color-temperature measurement [59]. Additionally, the laser light (λ = 1064) scattered by rBC-containing particles throughout their interaction with the laser was detected using avalanche photodiode detectors. Mass mixing ratios (MMRs) of rBC were calculated by integrating the mass of rBC in individual particles sampled in a volume of ambient air in 1 s and then scaling upward by factor of 1.05–1.15 on a flight-to-flight basis, to account for accumulation-mode rBC mass outside of the SP2 detection range. This was done following widely used approach in Schwarz et al. [59]. Only cloud–free air was analyzed, to ensure avoidance of inlet artifacts [60]. A cloud probe (Model LWS 485, Droplet Measurement Technology, Boulder, CO, U.S.) was used to identify cloud free air samples [61]. The SP2 was calibrated for rBC mass measurements by measuring size-selected (aerodynamic diameter of 125–250 nm) fullerene soot particles (Alfa Aesar, Lot # F12S011, Woodbridge, MA, U.S.); this is the calibration material accepted by the SP2 community for ambient measurements [62], and relied on the mass-to-mobility relationships for this specific material found in Moteki and Kondo [63] and Gysel et al. [64]. A total of four rBC calibrations of SP2 were performed during the study showing that instrument response was stable within 5% in rBC mass over the calibrations. The NOAA SP2 has a near-linear response to rBC mass over a majority of the particle accumulation mode mass range. The SP2 scattering measurements were calibrated with 220 nm polystyrene latex (PSL) spheres before each flight. We associate a total absolute uncertainty of 30% with the rBC MMR.
The number fraction of rBC that can be identified as coated (i.e., internally mixed), and the estimated coating thickness associated with particles with a narrow range of rBC mass content, were calculated as described in Schwarz et al. [33] for particles with rBC-components in the 150–175 nm volume-equivalent diameter range (3–5 fg rBC-mass). The total optical size of these particles was determined using the "leading-edge only" (LEO) fit method [65]. For this size range, there were good particle statistics and the SP2 optically sized the majority (>90%) of rBC-containing particles, reducing analysis biases. We were not able to size rBC-containing particles smaller than these without significant bias, and larger particle concentrations were significantly lower, reducing any utility due to statistical limitations. Mie theory core-and-shell modeling was used to translate the optical size of the particle, combined with the rBC mass content, and assumed indices of refraction, into "coating thickness" of the internal mixture. Although the absolute uncertainty in this measure can be quite high (especially for rBC-containing particles internally mixed with small fractions of non-rBC material), here we are concerned primarily with the precision of the technique, which is dependent on statistics (which are presented with results), and largely independent of the particular assumptions made about indices of refraction. Using the SP2-determined coating thicknesses, the sensitivity of the system merely for identifying internal mixtures (for the number-fraction of internally mixed rBC-containing particles) is estimated to be equivalent to a 25 nm coating thickness. Shell-and-core theory was also utilized to estimate the enhancements in rBC absorption from lensing effects of coatings. κ-Köhler theory was used to estimate water-uptake by rBC-containing particles as described in Schwarz et al. [13]. Here, to deal with ammonium, sulfate and OA mixed in coatings, a volume-weighted hygroscopicity parameter, κ, was used as described by Petters and Kreidenweis [66].
The instrumentation suite of the NOAA P3 is broadly described in an overview article [39]. Here, we use chemical composition of non-refractory aerosol from an aerosol mass spectrometer (AMS, 50% uncertainty, Aerodyne Research, Billerica, MA, U.S.) [67,68]; fine particle number concentrations from a nucleation-mode aerosol size spectrometer (NMASS, 26% uncertainty) [69], and an ultra high sensitivity aerosol spectrometer (UHSAS, Droplet Measurement Technology Inc., Boulder, CO, U.S.) [70]; nitrogen monoxide (NO, 3% uncertainty), nitrogen dioxide (NO2, 4% uncertainty), and total reactive nitrogen oxides (NOy, 12% uncertainty) mixing ratios from a chemiluminescence instrument [71,72,73]; CO (5% uncertainty) from a fluorescence instrument [74]; and SO2 mixing ratios (20% uncertainty) from a modified Thermo-Environmental Model 43C-TL instrument. The rBC data was time-synchronized with CO to account for differing time delays in sampling. Throughout this manuscript, all error estimates (± 1σ) represent one standard deviation of the mean unless stated otherwise.
The work presented here focuses rBC aging in coal-fired power plant plumes. The gas and aerosol measurements acquired during 5 daytime flights (06/11, 06/16, 06/22, 06/23, and 06/26) were utilized to study rBC aging in 6 plumes of 5 different coal-fired power plants: E. C. Gaston (G, Shelby county, AL), Harllee Branch (HB, Putnam county, GA), Scherer (S, Monroe county, GA), Harding Street (H, Marion county, IN), and Independence (I, Independence county, AR). The secondary fuel types for I and G were residual oil and diesel oil, respectively. The 6 investigated power plant plumes were selected from 14 plumes sampled throughout the campaign based on meeting all of the following criteria: (ⅰ) a distinct enhancement of SO2 in the plume over the surrounding (background) air was observed; (ⅱ) SO2 and NOy mixing ratios maximized in the closest plume intercept downwind of a power plant and then decreased with age of the plume due to dilution, chemical transformation, and dry deposition (i.e., the plume evolved without the injection of subsequent emissions downwind); (ⅲ) the bulk mass loadings of sulfate increased downwind of a power plant due to oxidation of SO2 to H2SO4, and H2SO4 condensation and coagulation, (ⅳ) there were negligible emissions of rBC from the power plants, and (ⅴ) the flights were carried out during the days when photo-chemistry occured. These conditions were selected specifically to allow investigation of changes in rBC microphysical properties in a "natural laboratory".
The E. C. Gaston power plant plume sampled on Saturday, June 22 2013 provided the clearest example of enhanced aging of rBC, and hence is used as a case study to illustrate our methods and results. A map of the plume-perpendicular transects colored as a function of NOy volume mixing ratio is shown in Figure 1. The plume transects show large enhancements of NOy and SO2 mixing ratios relative to the surrounding air mas (Figure S1). In this case the wind direction was easterly, with low and relatively constant wind speeds in the range of 2.5–3.0 m s−1. Wind speed coupled with distance from each given power plant were used to calculate the respective ages for all plume-perpendicular transects discussed in this manuscript. The WP-3D flew ~ 13 km upwind (~1.5 hr air transport time) of the power plant to sample the background air mass (notated "i0"), and then intercepted the plume twice downwind of the power plant. The first plume intercept was 17.6 km (~1.7 hr, i1) downwind of the power plant and 13.8 km east and upwind of Interstate 65 in the region south of Birmingham, AL. The second plume intercept (i2) was 48.2 km (~5.0 hr) downwind of the power plant, and 16 km west and downwind of Interstate 65. Although Interstate 65 was located between the first and the second plume intercepts, insignificant amounts of rBC and CO from traffic were mixed into the power plant plume, as shown in Figure 2; the stability of rBC concentrations in both plume intercepts downwind of the power plant were similar to those in the background, also indicating that that the power plant did not emit rBC. Further, any freshly emitted rBC would decrease observed internal mixing with rBC, as fresh emissions tend to consist mostly of bare rBC compared to backgrounds; there was no sign of reductions in coated number fraction. This allowed for the investigation of the aging of rBC in actual field conditions rather than those of a flow reactor. In this setting an air mass containing regional rBC population was "spiked" with known levels of NOy and SO2 and was chemically aged.
Compared to the regional background value of 1.9 ± 0.1 ppbv in i0, the mixing ratio of NOy increased to 8.0 ± 3.6 ppbv in i1 due to power plant stack emissions, and then decreased to 3.7 ± 0.4 ppbv in i2 as the plume was subjected to dilution and NOy was lost to deposition, particle formation or other processes (Figure 2). Although the measurements of SO2 across the upwind transect (i0) were unavailable due to an automated instrument calibration, mixing ratio changes similar to those observed in NOy are expected, a priori, for SO2 in each plume intercept. In fact, the low upwind SO2 mixing ratios observed along i0 (but not directly upwind) were followed by dramatic increases immediately downwind of the power plants that then decreased with each additional plume intercept further downwind (increase in plume age). This was consistently observed in all 6 investigated power plant plumes. The SO2 levels downwind of a power plant increased due to stack emissions from coal combustion and then decreased with time as the air mass experienced dilution and SO2 was exposed to chemical oxidation and deposition. The relative contributions of these processes to changes in SO2 are beyond the scope of this paper, where the focus is on an analysis of rBC evolution that is clearly linked to strong power plant influences.
The number fraction of detected rBC-containing particles that were identified as coated (CF) increased in the plume simultaneously with the increase in fine (<1 µm) particle mass loadings (Figure 3, rBC concentration did not significantly vary). The AMS indicated increases only in the sulfate fraction; organic aerosol, nitrate aerosol, and ammonium all remained quite stable; this indicates that sulfate was most likely the strongest contributor to internal mixing with rBC. To the extent that non-sulfate species are contributing to the observed rBC-aging in the plumes, this would enhance the impacts attributed to sulfate. The CF in the plume was 1–5% higher than in the air masses surrounding them, which we attribute to in-plume chemistry forming additional coatings on the pre-existing rBC-containing particles. The rBC CF increased from 81 ± 1% in i0 upwind of the power plant to 83 ± 1% in i1 and to 86 ± 2% in i2.
rBC coating thickness (CT) also increased in the aging plume. Figure 4 shows Gaussian fits to normalized histograms of CT from the upwind intercept, the two plume intercepts, and from the background air outside of the i2 intercept, labeled "i2 bkg". The i2 bkg represents air on each side of the i2 crossing within 30 s flight time of the plume. The mode of the Gaussian fit (x0) to the estimated rBC CT distribution at i0 was 29 nm, which was close to the mission-average rBC CT value of 30 ± 15 nm standard deviation for rBC-containing aerosol populations upwind of all 6 investigated power plant plumes. In i1, at the plume age of 1.7 hr, the mode was x0 = 36 nm, and by the age of 5 hr it had increased to x0 = 42 nm. The x0 for the estimated rBC CT in the i2 background, however, was x0 = 37 nm, suggesting that most of the estimated rBC CT growth observed in i2 was due to the contribution of more-aged regional rBC mixing in with the plume or common-mode aging independent of plume chemistry. Hence the average estimated rBC CT growth in the plume from 29 nm to 42 nm is an upper estimate for sulfate-driven aging of rBC in the plume.
To exclude common-mode aging effects on the estimated rBC CT and CF sulfate-driven chemistry changes, the estimated rBC CT and CF in the plumes were corrected by subtracting the changes (assumed due to non-plume SOA production generally) observed near, but uninfluenced by, the plumes. The background values of rBC CF and CT in 30 s flight increments (~3.5 km spatial resolution) on each side of the plume along the flight path were averaged to provide a representative background, and its growth independent of the plume chemistry, and subtracted from the in-plume values to calculate the net estimated rBC CT and CF changes in the plume. These values for Gaston are shown in Figure 5. The whiskers represent uncertainty calculated by propagation of errors (1σ) associated with the in-plume and out-of-plume 30 s averages. With aging, the estimated rBC CT and CF increased at a slightly faster rate in the E. C. Gaston power plant plume than in the air mass surrounding the plume. Before influence of the plume, at i0, rBC CT and CF were very stable on the upwind leg, as expected. The net estimated rBC CT and CF were 5.9 ± 2.1 nm and 3.1 ± 1.4% at i1, and 4.6 ± 2.3 nm and 4.3 ± 1.6%, respectively, at i2. At i2, this growth corresponded to enhanced internally mixed mass of ~0.2 fg, or ~15% of the total internally mixed material on the rBC-containing particles treated here.
Lower net estimated rBC CT and rBC CF increases were observed in intercepts of other power plant plumes (Figure 6). In some cases the values were even negative or decreasing with plume age, assumed to be a result of the limited signal-to-noise of our detection of these processes. The minimal coating increase is a consequence of multiple effects. First, most of the particle mass from sulfate-driven (or other) chemistry condenses or coagulates on the non-rBC-containing aerosol population; rBC-containing particles were generally a minor fraction (≤1%) of the total fine particle number in the air, and represented a similarly low fraction of the total aerosol surface area available to receive condensate. Secondly, the regional rBC-containing particles were already largely coated by OA components, consistent with estimated rBC CF and CT levels and large OA mass loadings observed in the background air masses. This is not surprising considering the large emissions of biogenic VOCs and the availability of solar radiation needed for photochemical oxidation are Southeast U.S. Potentially, a net increases in internal mixing of rBC due to plume-driven chemistry was overshadowed in some plumes by these processes.
If all sampled plume intercepts are considered, only a non-statistically relevant trend of increased estimated rBC CT and CF with increased with the age of the plume is observed, indicating that condensation of sulfate (and potentially other plume-chemistry products) was small. Compared to the rBC CT in the surrounding air mass, the net estimated rBC CT in power plant plumes was on average 2 ± 4 nm (±1 σ) and a median value of 3 nm. Similarly, the net rBC CF increased by 1.3 ± 1.3% (±1 σ) and the median value of 1.1%. The plumes were only investigated up to the age of ~5 hr due to the inherent difficulties in clearly identifying and characterizing older plumes that have intermixed via dilution with background air.
The air-parcel history, including rBC mixing state, upwind of each power plant strongly affected the rates of coating thickness growth in plumes, as is evidenced by different aging on different days of observation of the same power plant plume. For example, the E. C. Gaston plume was characterized on both June 11 and 22. On the 11th, the wind was westerly from the direction of Tuscaloosa, AL providing fresher emissions of rBC-containing aerosol with lower (~25%) rBC CF of 0.62 and ~50% lower estimated rBC CT (14 nm) than on the 22nd. The mass loadings of OA were also ~40% lower. Moreover, the number concentration of fine particles in the plume on 06/11 was roughly one half of that on 06/22 after 5 hr, suggesting that there were less non-rBC particles to impede the condensation of sulfuric acid (with potential contributions from other plume-generated species) on rBC cores. Hence, the fresher, less internally mixed rBC population collected more coating mass in the power plant plume with a lower number of fine particles, as demonstrated in Figure 6.
Although it is now clear that the sulfate condensed on to the rBC particle population is not a large contributor to the total coating thickness, it is useful to interpret the previous results in the context of rBC optics–and specifically rBC light absorption. This was done using a Mie theory core-shell optical model in conjunction with k-Köhler theory to describe water uptake by the aerosol in ambient conditions [64]. Mie core-shell theory is widely used in calculations of light scattering and absorption by atmospheric particles, and in evaluations of laboratory and field measurements [13,36,75,76,77]. The major limitation of the core-shell model is that it assumes that rBC-containing particles are composed of void-free and spherical core of rBC coated with a uniform layer of homogenous internally mixed material. Although ambient rBC particle have fractal characteristics, Qiu et al. [78] and China et al. [79] showed that with the development of organic coatings during atmospheric processing the rBC cores compact to a compact form, and that the consequent exposure to elevated humidity conditions results in near-spherical particles. These are typical conditions in summer-time Southeast U.S., which supports the use of the core-shell model here.
Our analysis was focused on the range of rBC-mass for which the SP2 was able to obtain the best statistics (3–5 fg rBC mass per particle). This range of mass is generally very near the mass-median diameter (i.e. typically half the total accumulation mode rBC mass occurs in particles with less than ~5 fg rBC mass, and half is in particles with more rBC mass than this). Over a limited rBC-mass range for which statistics are reasonable (up to ~10 fg), we do not observe significant trends in absorption enhancements due to coatings on rBC. Hence, we assume that our analysis is representative of the entire accumulation-mode population of rBC-containing particles.
The following microphysical conditions were considered in our calculations: 165 nm (~4 fg) rBC core that is initially coated with 30 nm of OA (κ = 0.1) takes up 2–9 nm of additional ammonium sulfate coating (κ = 0.53) as a direct result of power plant plume chemistry (Table 1). The rBC core size of d = 165 nm corresponds to the mode of the mass-size distribution of regional rBC population mixed into the 6 studied power plant plumes, and was within the range selected for the observations of rBC-containing particle optical size used to determine CTs. The average estimated rBC CT was determined to be 30 ± 15 nm upwind of all 6 investigated power plant plumes. The coatings on rBC particles were assumed to be completely composed of oxygenated OA. A density of 1.8 g cm−3 and complex refractive index of 2.26–1.26i were used for the rBC core. The real refractive indices 1.53 and 1.44 were used for ammonium sulfate and OA coatings respectively. The average, median, and maximum net estimated rBC CT enhancements observed in the investigated power plant plumes were 2 nm, 3 nm, and 9 nm, respectively; the average and median cases reflect very small changes; as such, their sensitivity to errors in assumed indecies of refraction (an indeed, of assumed composition) is negligible. Nitrate was excluded from the calculations because NO3− typically constituted only ~1% of the fine aerosol mass in power plant plumes.
CASE | rBC core diameter | OA coating thickness | OA κ/(NH4)2SO4 κ | OA RI/(NH4)2SO4 RI | (NH4)2SO4 coating thickness | RH |
Ⅰ | 165 nm | 30 nm | 0.1/0.53 | 1.44/1.53 | 2 nm | 50–95% |
Ⅱ | 165 nm | 30 nm | 0.1/0.53 | 1.44/1.53 | 3 nm | 50–95% |
Ⅲ | 165 nm | 30 nm | 0.1/0.53 | 1.44/1.53 | 6 nm | 50–95% |
Ⅳ | 165 nm | 30 nm | 0.1/0.53 | 1.44/1.53 | 9 nm | 50–95% |
We use Eq 1, below, to estimate the fraction (F) of the Southeast in which rBC could be influenced by power plant emissions within 5 hours aging:
(1) |
where AP is the estimated area of a single power plant plume after 5 hours (~180 sq. miles based on the first E. C. Gaston plume observations), ASE is the area of the Southeast, and N is the number coal- and oil-power plants that generate more than 25 MW yr−1 (166). Considering that the area of the entire Southeast U.S. (including LA, MS, AL, GA, FL, SC, NC, TN, AR, KY, WV, and VA) is 5.5 × 105 sq. miles, and assuming that the height of the PBL (for the purpose of this calculation defined as the mixed + transition layer [46]) is constant across the entire region, and assuming that the power plant plumes are rapidly mixed vertically throughout the PBL, only ~5% of the area of the Southeast U.S. is within the 5 hr old "plume cover".
To allow comparison to Stier et al. [14], in which the authors modeled the enhancements in column-integrated absorption optical depth, our calculations were extended above PBL, where typically 40% of rBC mass is located over the Southeast U.S. [45,46]. The atmospheric lifetime of rBC in our calculations was chosen to be 5 days to keep consistency with Stier et al. [14], this influenced the integrated impact of the additional coatings on the average column absorption of rBC-containing aerosol.
The results of Mie/k-Köhler theory framework calculations are shown in Figure 7. The enhancement in rBC absorption (AErBC: the calculated ratio of coated- to bare-rBC absorption) is shown as a function of relative humidity encountered during the study (50–95%). The total range of AErBC values calculated from our measurements is illustrated as solid orange fill. The range was generated using different ammonium sulfate coating thicknesses in the calculations with the mission-average rBC core size of 165 nm and OA coating thickness of 30 nm (Table 1). The 4 cases illustrate the variability in AErBC as a function of RH when: (Ⅰ) we use the mission average CT enhancement of 2 nm; this is our best estimate of the net enhancement; (Ⅱ) we use the mission median CT enhancement of 3 nm; (Ⅲ) we use the mission average CT + 1 standard deviation (6 nm), and (Ⅳ) we use the highest value observed for CT growth throughout the analysis, namely 9 nm. In each case we assume the plume-relevant added material is ammonium sulfate. The approximate range of AErBC for the total impacts of all anthropogenic SO2 emissions globally, predicted by the ECHAM5-HAM climate-aerosol model described in Stier et al. [14], is shown in solid yellow fill yellow for comparison. We carried out sensitivity studies to test assumptions about core size, coating κ, and coating index of refraction had negligible effects on the results of the calculations.
We assumed that aging effects within highly concentrated power plant plumes would be associated with substantially higher associated coating enhancement than expected for large scale sulfate impacts under more general conditions (i.e., global scale evaluated over relatively coarse resolution). Nevertheless, the calculated best-estimate of AErBC, based on the measurements of rBC –containing particle optical size during the SENEX study, was lower than an estimate by ECHAM5-HAM, the climate-aerosol model described in Stier et al. [14], which simulated general conditions over global scales (200 km resolution, using all anthropogenic SO2 emissions).
Based on Mie Shell-and-core theory, the addition of 2 to 3 nm of ammonium sulfate coating to a 165 nm rBC core originally coated by 30 nm OA in the power plant plume leads to regional additional AErBC due to sulfate uptake of ~ 2–6% (for 50–95% RH) in the Southeast U.S., which is ~5 times lower than the central prediction by the ECHAM5-HAM. In general, AErBC by sulfate coatings during SENEX were low because of the high levels of internal mixing of rBC with OA prior to exposure to coal-fired power plant plumes and the low number fractions of rBC-containing fine mode particles, such that surface area associated with rBC was a small fraction of the total available to receive condensates. In the presence of high levels of OA coatings (κ = 0.1) the addition of several nm of ammonium sulfate coatings (κ = 0.53) did not radically change the overall κ of the OA-(NH4)2SO4 coating "mixture", which was generally < 0.2 in our calculations. Thus, the water uptake and subsequent AErBC was significantly lower than what it would be if there were low levels of OA internally mixed with the rBC core. In addition, Lack and Cappa [76] previously showed that the AErBC can be impeded in the presence of internally mixed OA that can absorb light.
Two factors likely explain the stronger absorption enhancement simulated in Stier et al.: first, the representation of rBC in the model, and its transfer from a "bare" hydrophobic mode to an internally mixed mode (homogenously mixed with all mixed accumulation mode aerosol) after only the addition of a small amount of condensate. This occurs because the mixing state model of ECHAM-HAM (c.f. Table 1 in Stier et al., 2005) assumes that all rBC particles with coatings above a threshold (3 nm in Stier et al.) are transferred to the internally mixed accumulation mode, which significantly (and artificially) enhances the model simulated coating thickness. Secondly, the role of secondary organic aerosol in AErBC was neglected in the model treatment, which represented aerosol as emissions of primary organic matter, with rBC in which a fraction of freshly emitted rBC particles were treated as bare (i.e., completely water-inactive particles) that are grown to the accumulation mode hydrophilic particles by the addition of sulfate. Further, it is likely that all rBC particles freshly emitted from fossil fuel sources in the U.S. are associated with some degree of internal mixing, presumably from condensation of co-emitted semi-volatile organic vapors, certainly on the length scale of dilution to the model grid box scale. Formation of additional OA coatings from VOC oxidation chemistry in the region such as Southeastern U.S., which is characterized by some of the highest VOC levels in the nation, would likely provide substantial coatings on rBC-containing particles. Addition of low levels of (NH4)2SO4 to such rBC particles would have a minor impact on AErBC as explained above. However, an addition of even lower levels of (NH4)2SO4 to bare, hydrophobic rBC particles would lead to large estimated AErBC. For example, if a 165 nm rBC particle initially coated with only 0.1 nm of hydrophobic OA (k = 0.01) acquires 2 nm of (NH4)2SO4 in a power plant plume, the resulting AErBC is ~17%, which is in the range of values predicted by [14]. In fact, sophisticated tracking of particle mixing, for example with a particle-resolved model like PART-MC [80], would be needed to resolve these phenomena. Based on empirical observations, we conclude that the effect of the addition of sulfate coatings on the rBC population in the summertime Southeastern U.S. has a small impact (2–6%) on the AErBC due to large levels of OA internally mixed with the rBC population.
During the SENEX field study in the Southeast U.S. in June and July of 2013, there was an opportunity to investigate changes in microphysical properties of background (aged) black carbon (rBC) particles downwind of coal-fired power plants. We have reported constraints on enhanced aging of rBC-containing particles from measurements in 6 plume intercepts of 5 power plants during 5 daytime flights. The aging was quantified by measurements of: (1) the fraction of rBC particles that could be identified as associated with internally mixed material (CF), and (2) the thickness of coatings (CT) on rBC-containing particles in the plumes relative to un-perturbed air. The background regional rBC aerosol population was internally mixed with large levels of oxygenated organic aerosol (OOA), likely due to high rates of volatile organic carbon (VOC) oxidation and condensation of organic chemicals of low to moderate volatility. The CF for the background rBC was over 80%, where, on average, 165 nm diameter rBC cores were coated with 30 ± 15 nm of (assumed) oxygenated OA. Fresher rBC emissions, for example from local cities, traffic, etc., were generally associated with thinner coatings. Power plant plume chemistry increased estimated rBC CT and rBC CF only slightly faster than background, with weak correlations to increased mass loadings of sulfate and ammonium. Thus, rBC aging in the plumes was indeed affected by sulfate chemistry. The rBC aging rates in power plant plumes were believed to be low because of the competition between rBC (<1% of total fine particle number) and non-rBC particles as sites for condensation and coagulation of sulfuric acid.
These direct measurements constrain predictions of the impact of aging of rBC by total anthropogenic SO2 emissions (through sulfate-driven chemistry) on column-integrated aerosol absorption optical depth [14]. The observed net changes in estimated rBC CT were used in conjunction with Mie core-shell/κ-Köhler theory framework to show that the typical enhancements in absorption by rBC particles (AErBC) with the addition of internally mixed sulfate in the power plant plumes is likely in the 2–6% range. This type of processes, with condensation on an rBC particle population that is initially already largely internally mixed, requires a more sophisticated aerosol module version of ECHAM-HAM than used in Stier et al. [14] (such as was subsequently included [81,82]). In addition, the simplified representation of the aerosol mixing state in most aerosol microphysics schemes (except for models with significant complexity, such as [80]), in particular the assumption of internal mixing within modes or bins across climate model grid-scales, can lead to unrealistic simulated microphysics including mixing states and associated absorption enhancements. Large amounts of OA internally mixed with rBC, with lower κ values (~0.1) compared to that of sulfate (~0.5), reduce the relative impact of the water uptake and consequent additional light absorption by rBC containing particles due to sulfate. Here, we expect that the derived aging time scales, the coating thickness, and the fraction of coated particles are the best constraints we can supply for future model developments from direct SP2 measurements of rBC-containing aerosol; naturally the continued development of this type of dataset is critical to provide constraints on model processes.
The authors thank Smyrna/Rutherford County Airport staff for their help and support during the study, and the entire NOAA WP-3D flight crew, science cohort, and support staff for their efforts in the field. We remember G. Hübler with fondness and gratitude for his assistance with mission logistics, and thank K. Aikin for his assistance with the programming and data analysis. P. Stier acknowledges support from the UK Natural Environment Research Council project GASSP (NE/J024252/1). The NOAA Black Carbon Group was supported in part by the NOAA Atmospheric Composition and Climate Program, the NASA Radiation Sciences Program, and the NASA Upper Atmosphere Research Program.
[1] |
J. T. Beale, Large-time regularity of viscous surface waves, Arch. Ration. Mech. Anal., 84 (1983), 307–352. https://doi.org/10.1007/BF00250586 doi: 10.1007/BF00250586
![]() |
[2] |
R. C. Budzinski, T. T. Nguyen, G. B. Benigno, J. Doàn, Ján Mináč, T. J. Sejnowski, et al., Analytical prediction of specific spatiotemporal patterns in nonlinear oscillator networks with distance-dependent time delay, Phys. Rev. Lett., 5 (2023), 013159. https://doi.org/10.1103/PhysRevResearch.5.013159 doi: 10.1103/PhysRevResearch.5.013159
![]() |
[3] |
F. Chapeau-Blondeau, G. Chauvet, Stable, oscillatory, and chaotic regimes in the dynamics of small neural networks with delay, Neural Netw., 5 (1992), 735–743. https://doi.org/10.1016/S0893-6080(05)80134-0 doi: 10.1016/S0893-6080(05)80134-0
![]() |
[4] |
H. Chiba, A proof of the Kuramoto conjecture for a bifurcation structure of the infinite dimensional Kuramoto model, Ergod. Theory Dyn. Syst., 35 (2015), 762–834. https://doi.org/10.1017/etds.2013.68 doi: 10.1017/etds.2013.68
![]() |
[5] |
M. Y. Choi, H. J. Kim, D. Kim, H. Hong, Synchronization in a system of globally coupled oscillators with time delay, Phys. Rev. E, 61 (2000), 371–381. https://doi.org/10.1103/PhysRevE.61.371 doi: 10.1103/PhysRevE.61.371
![]() |
[6] |
J. D. Crawford, Amplitude expansions for instabilities in populations of Globally-Coupled oscillators, J. Stat. Phys., 74 (1994), 1047–1082. https://doi.org/10.1007/BF02188217 doi: 10.1007/BF02188217
![]() |
[7] |
S. Y. Ha, Q. Xiao, Remarks on the nonlinear stability of the Kuramoto-Sakaguchi equation, J. Diff. Eq., 259 (2015), 2430–2457. https://doi.org/10.1016/j.jde.2015.03.038 doi: 10.1016/j.jde.2015.03.038
![]() |
[8] |
S. Y. Ha, Q. Xiao, Nonlinear instability of the incoherent state for the Kuramoto-Sakaguchi-Fokker-Plank equation, J. Stat. Phys., 160 (2015), 477–496. https://doi.org/10.1007/s10955-015-1270-5 doi: 10.1007/s10955-015-1270-5
![]() |
[9] | H. Honda, A. Tani, Mathematical analysis of synchronization from the perspective of network science, Mathematical Analysis of Continuum Mechanics and Industrial Applications (Proceedings of the International Conference CoMFoS15), Singapore: Springer, 2017. |
[10] |
H. Honda, Global-in-time solution and stability of Kuramoto-Sakaguchi equation under non-local coupling, Netwo. Heterog. Media, 12 (2017), 25–57. https://doi.org/10.3934/nhm.2017002 doi: 10.3934/nhm.2017002
![]() |
[11] | Y. Kawamura, From the Kuramoto-Sakaguchi model to the Kuramoto-Sivashinsky equation, Phys. Rev. E, 89 (2014), 010901. |
[12] |
S. Kim, S. H. Park, C. S. Ryu, Multistability in coupled oscillator systems with time delay, Phys. Rev. Lett., 79 (1997), 2911–2914. https://doi.org/10.1103/PhysRevLett.79.2911 doi: 10.1103/PhysRevLett.79.2911
![]() |
[13] | Y. Kuramoto, D. Battogtokh, Coexistence of coherence and incoherence in nonlocally coupled phase oscillators, Nonlinear Phenom. Complex Syst., 5 (2002), 380–385. |
[14] | O. A. Ladyženskaja, V. A. Solonnikov, N. N. Ural'ceva, Linear and Quasi-linear Equations of Parabolic Type, Rhode Island: American Mathematical Society, 1968. |
[15] | M. Lavrentiev, R. S. Spigler, Existence and uniqueness of solutions to the Kuramoto-Sakaguchi nonlinear parabolic integrodifferential equation, Differ. Integr. Equ., 13 (2000), 649–667. |
[16] | M. Lavrentiev, R. S. Spigler, A. Tani, Existence, uniqueness, and regularity for the Kuramoto–Sakaguchi equation with unboundedly supported frequency distribution, Differ. Integral Equ., 27 (2014), 879–892. |
[17] | W. S. Lee, E. Ott, T. M. Antonsen, Large Coupled Oscillator systems with heterogeneous interaction delays, Phys. Rev. Lett., 103 (2009), 044101. |
[18] |
E. Niebur, H. G. Schuster, D. M. Kammen, Collective frequencies and metastability in networks of limit-cycle oscillators with time delay, Phys. Rev. Lett., 67 (1991), 2753–2756. https://doi.org/10.1103/PhysRevLett.67.2753 doi: 10.1103/PhysRevLett.67.2753
![]() |
[19] | H. Risken, The Fokker-Planck Equation, Berlin: Springer, 1989. |
[20] |
N. Tanaka, A. Tani, Surface waves for a compressible viscous fluid, J. Math. Fluid Mech., 5 (2003), 303–363. https://doi.org/10.1007/s00021-003-0078-2 doi: 10.1007/s00021-003-0078-2
![]() |
[21] |
S. H. Strogatz, E. Mirollo, Stability of incoherent in a population of coupled oscillators, J. Stat. Phys., 63 (1991), 613–635. https://doi.org/10.1007/BF01029202 doi: 10.1007/BF01029202
![]() |
[22] | J. Wloka, Partielle Differentialgleichungen: sobolevraume und Randwertaufgaben, Stuttgart: Teubner, 1982. |
[23] |
M. K. S. Yeung, S. H. Strogatz, Time delay in the Kuramoto model of coupled oscillators, Phys. Rev. Lett., 82 (1999), 648–651. https://doi.org/10.1103/PhysRevLett.82.648 doi: 10.1103/PhysRevLett.82.648
![]() |
[24] |
T. Zhu, Synchronization of the generalized Kuramoto model with time delay and frustration, Netw. Heterog. Media, 18 (2023), 1772–1798. https://doi.org/10.3934/nhm.2023077 doi: 10.3934/nhm.2023077
![]() |
1. | Alexander Keimer, Lukas Pflug, 2023, 15708659, 10.1016/bs.hna.2022.11.001 | |
2. | F. A. Chiarello, J. Friedrich, S. Göttlich, A non-local traffic flow model for 1-to-1 junctions with buffer, 2024, 19, 1556-1801, 405, 10.3934/nhm.2024018 | |
3. | Veerappa Gowda G. D., Sudarshan Kumar Kenettinkara, Nikhil Manoj, Convergence of a second-order scheme for non-local conservation laws, 2023, 57, 2822-7840, 3439, 10.1051/m2an/2023080 | |
4. | Alexander Keimer, Lukas Pflug, Discontinuous nonlocal conservation laws and related discontinuous ODEs – Existence, Uniqueness, Stability and Regularity, 2023, 361, 1778-3569, 1723, 10.5802/crmath.490 | |
5. | F. A. Chiarello, H. D. Contreras, L. M. Villada, Existence of entropy weak solutions for 1D non-local traffic models with space-discontinuous flux, 2023, 141, 0022-0833, 10.1007/s10665-023-10284-5 | |
6. | Aekta Aggarwal, Ganesh Vaidya, Convergence of the numerical approximations and well-posedness: Nonlocal conservation laws with rough flux, 2024, 0025-5718, 10.1090/mcom/3976 |
CASE | rBC core diameter | OA coating thickness | OA κ/(NH4)2SO4 κ | OA RI/(NH4)2SO4 RI | (NH4)2SO4 coating thickness | RH |
Ⅰ | 165 nm | 30 nm | 0.1/0.53 | 1.44/1.53 | 2 nm | 50–95% |
Ⅱ | 165 nm | 30 nm | 0.1/0.53 | 1.44/1.53 | 3 nm | 50–95% |
Ⅲ | 165 nm | 30 nm | 0.1/0.53 | 1.44/1.53 | 6 nm | 50–95% |
Ⅳ | 165 nm | 30 nm | 0.1/0.53 | 1.44/1.53 | 9 nm | 50–95% |