We sincerely appreciate your detailed and thoughtful review of our manuscript, as well as the insightful comments and valuable suggestions you have provided. Your feedback has been instrumental in helping us improve the quality and clarity of our work. We have carefully considered each of your remarks and have done our best to address them thoroughly in our responses and revisions.

We are truly grateful for the time and effort you have dedicated to evaluating our paper, and we thank you for your constructive input that has contributed significantly to enhancing our study.

RC 1: “Methods: this section would greatly benefit from a schematic or a workflow diagram summarizing the 7 steps, including which datasets/models are used at each step.”

AC:  As suggested, we have prepared a workflow diagram, which will be included in the manuscript. The diagram is attached to this response (please find it in the suplement as workflow_diagram.pdf). Please note that the steps 6 and 7 were combined as a one for the clarity of explaining how to calculate the final contribution to BB AOD from each potential region as a source of BB aerosol. 

In Line 148 reference to this workflow is added:

“To estimate the contribution of identified regions to the BB AOD in selected location, such procedure was followed. The workflow summarizing the steps of the methodology is depicted in Figure 1.”

RC 1: “I would also recommend that the authors mention before going into the 7 steps that the overall goal of the methodology is to partition the BB AOD according to emissions encountered from different source regions.”

AC: In line 149 it is mentioned that the methodology is implemented to “estimate the contribution of identified regions to the BB AOD in the selected location”. Moreover, it is also mentioned in the text in other places e.g: in the line 69: “To assess the contribution of BB regions to BB AOD at a selected location during the BB season, a new framework was developed.”, in the line 151: “To estimate the contribution of identified regions to the BB AOD in selected location, such procedure was followed.”. Therefore, we believe that adding this information again is not justified. 

RC 1.1: “How were the specific regions selected for trajectory analysis (colored shapes in Fig 1a) selected?” “Perhaps a map of the probability of plumes reaching Poland would be useful here. As is, the presented trajectories appear to be somewhat cherrypicked.”

AC:  Points for trajectory analysis were selected based on the mean annual number of fire events reported within each region by the MODIS Active Fire Product (Near Real Time). For clarity in the figure, only pixels with more than 10 fires per year are displayed. Points for forward trajectory analysis were chosen by identifying pixels with the highest number of fire detections. When multiple high-activity pixels were located in close proximity, only one representative point was selected to avoid redundancy. In summary, the selection of points aimed to include areas with elevated fire activity while excluding calculations for points situated very close to one another. 

Due to computational constraints, it was not feasible to simulate forward trajectories from every pixel with a high annual number of fires. The analysis period—May to September, 2006–2022—would result in 2,601 trajectories per point, making full coverage impractical. Therefore, only a limited number of representative points were used.

The primary aim of the forward trajectory analysis was to assess whether BB aerosol from a given source region could potentially reach Poland, without explicitly relying on these results in the contribution estimation. As such, we consider the use of a few well-chosen points sufficient for this type of analysis.

RC 1.2: “Line 136: please move figure 7 to the first figure discussed, to orient the readers to the source regions of interest.”

AC: Figure 7 is in such a place because it is a result of analysis of which regions can be the source of BB aerosol flowing over Poland. We don’t know which regions will be considered in estimating the contribution to BB AOD in Warsaw before having results presented in Figures 1-6, so we’d rather leave it after the forward trajectory analysis. Figure 7 is placed at this point because it presents the results of an analysis identifying the regions that could be the sources of biomass burning (BB) aerosol transported over Poland. Since the regions contributing to BB aerosol optical depth (AOD) in Warsaw cannot be determined prior to examining the results shown in Figures 1–6, we chose to position Figure 7 following the forward trajectory analysis.

RC 1.3: “Step 1, Line 150: are backward trajectories initialized over Poland/Warsaw?”

AC: While testing the methodology for Warsaw the backward trajectories are initialized in Warsaw. But line 150 is a part of the general description of methodology that is why initialization of backward trajectory shall be done in the location selected for the analysis. In this general description we didn’t want to narrow down to Warsaw. specify that we did it for Warsaw. 

In the revised manuscript the line 150 will be corrected in such a way that it will be clear that the backward trajectories shall be initialized in the location which is selected for the analysis.

We suggest such modification: 

“For each day, generate backward trajectories originating at selected location at 12:00 UTC. “

instead of: 

“For each day, generate backward trajectories starting at 12:00 UTC.“

RC 1.5: “Step 4: The discussion of fire outbreaks and emissions accounting is very unclear.

AC: To make the description of Step 4 more approachable and understandable, we propose rewriting the text as follows. Please note that our methodology has slightly changed after suggestions from another referee. Now, there is no rigid PBL height threshold, but the PBL height is changing along the trajectory. Such PBL height is now given by the HYSPLIT model.

“4. Identification of fire outbreaks:

Within the defined dispersion areas at each trajectory point, check for fire outbreaks that occurred on the same day with the outbreak hour, not later than the hour of the trajectory point. If any fire outbreaks are found, assign them to one of the regions which during preparatory work were found to be a possible BB aerosol source. For each such source of BB aerosol and each dispersion area, apply three different methods to account for the fire emissions:

• No Threshold method:

For each fire pixel within the dispersion area, calculate the fire emission as the product of the fire flux and the fire area. Then, compute the mean fire emission by dividing the sum of all fire emissions by the number of fire pixels.

• PBL method:

Check whether the trajectory altitude is below the PBL height. If so, calculate the mean fire emission as described in the No Threshold method.

• CAMS method:

Calculate the mean fire emission as described in the No Threshold method. Then, multiply it by the probability that the fire plume is elevated to the altitude of the trajectory point or higher. This probability is obtained by fitting a cumulative distribution function (CDF) to the CAMS GFAS top of the plume altitude data for the considered BB source region. A separate CDF is fitted for each BB source region. The plume-top altitude h is modeled using a log-normal distribution, for which the probability density function (PDF) is given by

$

f(x; \mu, \sigma) = \frac{1}{x \sigma \sqrt{2\pi}} \exp \left( -

\frac{(\ln x - \mu)^2}{2 \sigma^2} \right),

$

where h represents the plume-top altitude expressed in m ASL, while μ and σ are the mean and standard deviation of the natural logarithm of h, ln h, respectively. The CDF, representing the probability that the plume altitude does not exceed h, is defined as:

$

\mathrm{CDF}(h) = P(X \leq h) = \int_0^h f(x; \mu, \sigma)\, dx.

$

Accordingly, the probability that the fire plume extends beyond altitude h is:

$

P(X > h) = 1 - \mathrm{CDF}(h).

$

This log-normal fit yields a Pearson correlation coefficient r of 1 with the observational data, indicating an excellent fit.”

RC 1.5: What does “source region” mean in line 177?

AC: Source region meant the region which was identified during preparatory work as a possible source of BB aerosol. This will be rephrased in the revised manuscript in the following way:

“If any fire outbreaks are found, assign them to one of the regions which during preparatory work were found to be a possible BB aerosol source. ” 

instead of:

“If any fire outbreaks are found, assign them to each potential source region.”

RC 1.5: “Won’t the No Threshold and PBL method count the same emissions multiple times? Please clarify this.”

AC:  Yes, this may happen — each method can potentially count the same emission multiple times. However, we aimed to mitigate this by assigning a single mean emission value to each dispersion area, rather than assigning the whole sum of all emissions within it. We believe that counting the same emission more than once can reflect the dynamics of the trajectory: if a trajectory remains longer over a region with intense fire activity, the likelihood and amount of smoke being transported to Poland increases compared to a case where the trajectory passes over fire-affected areas only briefly and is counted once.

RC 1.6: “Step 5: it might make sense to refer to “each method” as “each emissions method” or something similar.”

AC: Changed as suggested; “each method” will be changed to “each emissions method”. 

RC 1.7: “Step 6: What does “specific region” mean in line 202?”

AC: There is no “specific region” in line 202. Did RC1 mean “specific location”? If yes - we meant the location selected for the analysis in the methodology. We suggest rewriting this in the revised manuscript as:

“BB AOD given at 12:00 UTC for that day at the analysis location.”

instead of:

“BB AOD value in specific location for the selected day for 12UTC.”

RC 1.8: “Step 7: make it clear that the contribution is a percentage.”

AC: We suggest combining step 6 and step 7 to make these two steps more clear in the following way:

“Step 6: BB AOD-weighted contribution analysis: For each study day, and for every emission method, source region, and dispersion radius, multiply the mean fire emission by the NAAPS-mode BB AOD given at 12:00 UTC for that day at the analysis location. Then sum these BB AOD-weighted emissions over the entire study period and express each regional total as a percentage of the all-region sum for the given method and dispersion radius. The resulting percentages quantify how much each source region contributes to the overall BB AOD.”

Here it is also highlighted that we express this contribution as a percentage. 

RC 1.8: “What does “BB AOD across mean”?”

AC: Thank you for catching that. That was a typo, and the correct phrase is “BB AOD in the analysed location”. Corrected in the revised manuscript.

RC 2:”Alaska: Why was Alaska not included as a potential region impacting Warsaw BB AOD? The authors note that lofted plumes have a 1% chance of reaching Warsaw, which seems at odds with their threshold of excluding plumes which have a <0.5% chance of reaching Warsaw.”

AC: Alaska will be included as a potential source of BB aerosol in the revised manuscript. We conducted additional calculations to present the results for this region as well. 

It turned out that during the BB season, the mean Alaska contribution equals 2.6 % ± 2.1 %. The contribution was minimal in May (0.2 % ± 0.2 %) and reached its maximum in July (6.1 % ± 0.4 %). The No threshold method overestimated the Alaska contribution by up to 3.7% in July, while the PBL method underestimated the contribution, showing for each month the contribution below 1.5%. More thorough description of the results will be included in the revised manuscript. 

RC 3: “Figures: throughout the manuscript there are figures with monthly or regional subplot maps. My recommendation is to add subplot titles identifying the plotted month or region.”

AC: Subplots titles were added. 

RC 3: “I would also recommend moving tables 1 and 2 to the supplement and presenting this information in Figure 11 instead.”

AC: We appreciate the suggestion; however, we would rather keep Tables 1 and 2 in the main manuscript, as the numerical values they provide are frequently referenced and discussed in the results section—for example, when analysing how changes in the dispersion radius affect regional contributions, or when comparing absolute differences between methods. In our view, it is not feasible to present all this detailed information clearly in a single figure without compromising readability. Figure 11 already presents the maximum amount of aggregated information in a visual form, while the tables complement it by offering precise numerical values necessary for in-depth interpretation.

We are deeply grateful for your meticulous review and constructive critiques of our manuscript. Your careful analysis and perceptive suggestions have greatly strengthened both the scientific rigor and the overall presentation of our study. We have addressed each point in detail and believe the revisions reflect the clarity and depth you encouraged. Thank you for investing your time and expertise to help us enhance the quality of our work.

GENERAL COMMENTS

RC: “An important (rather old) reference on modelling of long-range transport of wildfire aerosol with trajectory models is Forster et al. (2001) which should be mentioned in the introduction or in connection with the results presented for the contribution of Canadian aerosol.”

AC: We thank the reviewer for this suggestion. The reference to Forster et al. (2001) will be added in the revised manuscript, in the introduction part. 

We suggest mentioning this reference in Line 46:

“While several studies have reported incidents of BB aerosol influx over Europe, comprehensive analyses of these phenomena, particularly over Central Europe, are still scarce. Nevertheless, smoke layers over Central Europe have been documented in earlier work. A pronounced aerosol layer at 3–6 km over Germany in August 1998 was attributed to Canadian wildfires (Forster et al., 2001). ”

RC: “Note that in the paper the altitude of the BB aerosol layer  was between 3 and 6 km over Germany, i.e. partially above the maximum altitude in  Figures 1-3.”

AC: Figures 1–3 present the probability that an air parcel, initialised at a given starting altitude and location, will reach Poland—not the altitude at which it arrives over Poland. Therefore, the figures do not indicate the vertical distribution of aerosol over Poland, but rather the likelihood of long-range transport from different altitudes.

RC: “It is confusing that in the paper 3 different time periods are used.”

AC: We agree this can be confusing, but each month‑restricted period serves a specific purpose. BB AOD was first shown for March–October to acknowledge the seasonal cycle, including elevated April values in Europe that are mainly linked to agricultural burning (not wildfires). Our contribution analysis is limited to May–September to focus on the wildfire‑related BB season and to avoid mixing in the agricultural signal from April. Finally, the FLAMBE emission record had to include April because 10‑day backward trajectories initialised on 1 May extend back into ~20 April; without April emissions, early‑season transport would be truncated.

Also in the revised manuscript, all time periods are unified to time period 2006–2022, except for the use of MODIS Active Fire Product data in the fire frequency analysis, which will retain the extended period of 2001–2022. We believe that this longer time span provides a broader climatological perspective on fire activity in the potential BB source regions.

SPECIFIC COMMENTS

RC: “Line 38 and/or later: Please provide information on typical plume altitudes.”

AC: We agree and have added a short overview of typical biomass‑burning plume altitudes in the Introduction (in Line 46): 

“While several studies have reported incidents of BB aerosol influx over Europe, comprehensive analyses of these phenomena, particularly over Central Europe, are still scarce. Nevertheless, smoke layers over Central Europe have been documented in earlier work. A pronounced aerosol layer at 3–6 km over Germany in August 1998 was attributed to Canadian wildfires (Forster et al., 2001). Potential smoke layers over Warsaw spanning a range of source regions were examined by (Janicka et al., 2023), who found that layers that could originate from North America occurred between about 2 and 8.5 km, whereas layers attributed to Eastern Europe were most frequently detected between 2 and 4 km, with some reaching as high as 7 km. In addition, (Ortiz-Amezcua et al., 2017) reported smoke layers of 1–2 km thickness located at roughly 5 km above sea level (ASL) over Granada and Leipzig, and around 2.5 km ASL over Warsaw, linked to long-range transport from North America.”

RC: “Line 60 and 210: Months inconsistent to peaks in Fig. 9 (expand!).”

AC: In lines 284-288 it is explained: “Interestingly, BB AOD values are higher in March and April compared to May and June. This pattern can be attributed to early spring BB practices, such as agricultural waste burning and land management fires, which are prevalent during the planting season, particularly in Eastern Europe and Russia (Stohl et al., 2007; McCarty et al., 2012; Hall et al., 2021)”. 

In our analysis we focus on BB associated with wildfires, not agricultural practices. That is why our contribution analysis is limited to May–September to focus on the wildfire‑related BB season and to avoid mixing in the agricultural signal from April. 

RC: “Line 97: A reference for FLAMBE is missing or is this included in the reference in line 102 or one of the other references in the NAAPS-section? Please clarify or refer at least to the next section where the same references appear again.“

AC:  The missing reference is the one given in line 102. It was added in line 97.  

“FLAMBE program (Reid et al., 2009), initiated in 1999, is a collaboration between the U.S. Navy, NASA, NOAA, and the academic community. It integrates fire detection algorithms like NOAA/NESDIS’s Wild-Fire Automated Biomass Burning Algorithm (WF_ABBA) and NASA’s MODIS fire products to monitor BB emissions, incorporating these data into the NAAPS model to study smoke particle emissions and their atmospheric transport on regional to continental scales (Reid et al., 2009).”

RC: “Both sections might be merged for clarity.”

AC: We acknowledge the suggestion to merge the FLAMBE and NAAPS sections for clarity. However, in our opinion, keeping them separate helps distinguish between the different roles these tools play in our methodology. NAAPS is used primarily as a source of aerosol-related output (AOD and BB AOD), while FLAMBE is used independently to extract information on fire size and fire emissions. Although FLAMBE feeds emission data into the NAAPS model, in our analysis, both tools are applied separately and serve distinct purposes. We believe that maintaining this distinction improves the clarity of our methodological description.

RC: “Line 139: Is it possible to perform simulations in Europe also for higher altitudes for consistency with the next sections and Fig. 10?”

AC: We propose a modification in response to this comment. In the main part of the manuscript, we will limit the maximum starting altitude for estimating the probability that a trajectory from a given region reaches Poland to 5 km AGL, ensuring consistency across all regions. Higher altitudes (up to 10 km AGL) were originally included for North American regions due to their specific transatlantic atmospheric circulation patterns. In particular, for Alaska, this approach yielded interesting results. However, after considering the reviewer’s suggestion, we will move the results for altitudes up to 10 km AGL to the appendix and focus on altitudes up to 5 km AGL in the main figures.

We believe this change is justified, as the probability of smoke being lofted above 5 km AGL is very low, and omitting these higher layers in the core analysis will not affect the overall conclusions. Therefore, we will agree the altitude range used for North American regions (Figures 4c, 5c, 6c) with that used for European regions (Figures 1c, 2c, 3e,f), while keeping the results for higher altitudes as supplementary material for completeness.

RC: “Eqn. 1 to 3: Please use consistent symbols for latitude and longitude. It is rather confusing to use Φ for latitude in one equation and for longitude in the others.”

AC: We revised symbols for latitude and longitude. We decide to use θ for latitude and ϕ for longitude. 

RC: “Lines 190ff: Relation of CDF to f(x,..)  and x to h? An additional equation would be useful. Include also 'r' (line 333 too late) and be consistent with later text (line 333ff) concerning inclusion of 'ln' in mean and standard deviation. The text should be rearranged here, caption of Fig. 10 appears to be OK.”

AC:  In the revised manuscript, the CDF will be expressed as a function of plume‑top altitude h in m ASL. Whole CAMS method section (Lines 190ff) will be changed to:

“CAMS method:

Calculate the mean fire emission as described in the No Threshold method. Then, multiply it by the probability that the fire plume is elevated to the altitude of the trajectory point or higher. This probability is obtained by fitting a cumulative distribution function (CDF) to the CAMS GFAS top of the plume altitude data for the considered BB source region. A separate CDF is fitted for each BB source region. The plume-top altitude h is modeled using a log-normal distribution, for which the probability density function (PDF) is given by

$

f(x; \mu, \sigma) = \frac{1}{x \sigma \sqrt{2\pi}} \exp \left( -

\frac{(\ln x - \mu)^2}{2 \sigma^2} \right),

$

where h represents the plume-top altitude expressed in m ASL, while μ and σ are the mean and standard deviation of the natural logarithm of h, ln h, respectively. The CDF, representing the probability that the plume altitude does not exceed h, is defined as:

$

\mathrm{CDF}(h) = P(X \leq h) = \int_0^h f(x; \mu, \sigma)\, dx.

$

Accordingly, the probability that the fire plume extends beyond altitude h is:

$

P(X > h) = 1 - \mathrm{CDF}(h).

$

This log-normal fit yields a Pearson correlation coefficient r of 1 with the observational data, indicating an excellent fit.”

We will also clarify that, within the CAMS method, HYSPLIT trajectories must be initialised also with altitude output in m ASL so that h in the CDF corresponds directly to the trajectory altitude. 

“1. Backward trajectory simulation:

For each day, generate backward trajectories originating at selected location at 12:00 UTC. Extend each trajectory 240 hours (10 days) backward in time. Initialize each trajectory at multiple altitudes, starting from 500 m up to 4000 m in increments of 500 m AGL, and additionally at 5000, 6000, 7000, 8000, 9000, and 10000 m AGL. For every starting height, run HYSPLIT twice: once with the height specified in metres ASL (m ASL) and once in metres above ground level (m AGL), with the latter run computed together with PBL height output (expressed in m AGL). After computation, merge each ASL–AGL pair so that the resulting trajectory record contains both m ASL and m AGL height fields.”

Referring to the remarks about lines 333ff we can suggest such modification of lines 331-333:

“Using the data of the altitude of the top plume provided by CAMS GFAS during the years 2006-2022 for the months of May–September, density histograms of this parameter were plotted density histograms of this parameter were plotted, alongside the fitted log-normal PDF f (x; μ, σ), described by the parameters μ and σ, which are the mean and standard deviation of ln x, respectively (see Figures 11a–11f). For each fit, the Pearson correlation coefficient r was calculated. ”

RC: “Figures 1 to 6: Are the shown data only for May to September? If yes include in captions.”

AC: Yes, these data cover only months May-September only. This will be clarified in the figure captions and also in the text of the manuscript in sec 3.1. 

Figure caption (e.g. for USA) is corrected to:

“Figure 4. Spatial analysis of fire events and trajectory analysis in the USA. Panel (a) displays the average number of fires during BB season (in months May-September) in the USA over the period 2001-2022, with fire frequency represented on a logarithmic color scale. Panel (b) shows the specific locations selected for trajectory analysis within this region. Panel (c) presents the percentage of trajectories reaching Poland during BB season from each location defined in (b) as a function of starting altitude. Each symbol corresponds to a specific source point, as indicated in the legend.”

In sec 3.1 Lines 214-216 are corrected to:

“The average number of fires in BB season (in months May-September) for the years 2001–2022 retrieved from MODIS Fire Active Product data is presented in Figures 1a, 2a, 3a-b for Europe and Russia, and in Figures 4a–6a for North America, including the United States, Alaska, and Canada.”

RC: “Figure 9: I suppose BB AOD and total AOD are meant in caption.”

AC: Yes, in this context, AOD refers to total AOD. In the revised version of the manuscript, the figure caption will be updated to explicitly state that total AOD is shown.

Figure 9 caption is corrected to:

“Figure 9. Monthly values of BB AOD, total AOD, and their ratio for Poland and Warsaw averaged over the period 2006-2022. Panel (a) displays the mean BB AOD values for each month, comparing Poland and Warsaw. Panel (b) shows the monthly mean total AOD values for both Poland and Warsaw. Panel (c) presents the ratio of BB AOD to total AOD, expressed as a percentage, for each month, comparing values between Poland and Warsaw. ”

RC: “Line 337ff: Isn't there also the effect of fires on mountains or other elevated regions?”

AC: Yes, terrain elevation does influence the plume-top altitude distributions. In the revised manuscript, we will update the discussion of the fitted PDFs for the identified BB aerosol source regions to reflect this aspect. We thank the reviewer for this helpful remark. We suggest such modifications for Lines 337-351:

“The fitted PDF parameter μ, representing the mean of the natural logarithm of the plume-top altitude in a log-normal distribution, can reveal differences in plume-elevation dynamics. From a physical standpoint, the quantity  = e is more informative, as it corresponds directly to the mean plume-top altitude. However, because CAMS GFAS altitudes are reported in m ASL, local terrain elevation must also be taken into account. Among the European regions, Southern Europe exhibits the lowest mean plume-top altitude,   = 1110 m ASL (μ = 7.01), followed by the Iberian Peninsula,  $\tilde{\mu}$ = 1150 m ASL (μ = 7.05). Eastern Europe shows a higher value of $\tilde{\mu}$  = 1460 m ASL (μ = 7.29). In North America, plume tops are generally higher: Canada  $\tilde{\mu}$ = 2100 m ASL (μ = 7.65); the USA $\tilde{\mu}$  = 1840 m ASL (μ = 7.52); and Alaska  $\tilde{\mu}$ = 2040 m ASL (μ = 7.62).

Since CAMS GFAS data are reported in m ASL, and the variability of terrain elevation significantly affects the interpretation of plume rise above ground level, it is not possible to clearly predict how the different methods will estimate contributions. Based on the forward trajectory analysis for Alaska, the No Threshold method is likely to overestimate the contribution from that region. Aside from this case, there is no clear basis to determine how the results of the various methods will differ for certain BB aerosol source regions. Among the methods considered, the CAMS-based approach is expected to provide the most reliable estimates of regional contributions to BB AOD, as it statistically accounts for differences in fire plume rise dynamics across the identified BB aerosol source regions.”

RC: “Line 381: Cite Table 1 and 2.”

AC: Corrected. In the revised manuscript Line 381 is corrected to:

“Changing the dispersion radius significantly affects the contributions from different regions to BB AOD in Warsaw (Tables 1, 2).”

RC: “Fig. 11 and 12: Include Alaska.”

AC: Will be included in the analysis. Please find attached new Figures to this reply as Figures12_13.pdf. Please note the numbering has changed after suggestion from another referee to add  workflow diagram summarazing the steps in our methodology. In the revised manuscript the results obtained for Alaska are added in the text. 

RC: “Eastern Europe with or without European part of Russia?”

AC: Eastern Europe comprises the European part of Russia (it is discussed in lines 266-271 and depicted on Fig. 7a.) 

RC: “Fig. 12, caption: Do the error bars represent the 4 dispersion radii?”

AC: Yes, the averaging was done not only over months in BB season, but also over 4 dispersion radii as it is written in the caption of Fig 12: “Contributions were determined using the CAMS method and averaged over months May-September and averaged over four dispersion radii (5 % , 10 % , 15 % , and 20 % of the trajectory length).” 

TECHNICAL CORRECTIONS

RC: “Line 4: Spell out MODIS.”

AC: Corrected.

RC: “Line 56: Improve wording.”

AC: Instead of “Given the diverse vegetation types in different source regions and the aging processes that BB aerosol undergoes during transport, it is important to determine their origins.” We can suggest change to: “given the diverse vegetation types across source regions and the aging processes that BB aerosol undergoes during transport, identifying its origins is essential.”

RC: “Line 136: Poland only in Fig. 1b or 7.”

AC: There should be Fig. 7 - this will be corrected.

“Subsequently, the probability of air parcels arriving over Poland from these fire locations was calculated as the percentage of simulated forward trajectories from each source point that reached Poland, defined within coordinates [49.0︒N, 55.0︒ N] × [14.0︒ E, 24.2︒ E] (see Figure 7a).”

RC: “Line 191: Typo”, “Line 195: Include equation number.”, “Lines 232ff: Cite figures.”, “Line 267: Should be Fig. 3.”, “Line 372: Wrong citation style (use \citep).”, “Figures 10 and A1: Typo in legends and labels (10 times each!)”, “Figure 13c: Typo in legend.”, “Lines 558, 559, 566: Acronyms with first name?”, “Line 588: Use correct abbreviation or full name for journal.”

AC: Corrected.

RC: “Line 217: Southwestern and Fig. 2? Meant or misleading?”

AC: On the Iberian Peninsula the number of fires is lower (maximum of 70 fires per pixel) than in Southern and Southeastern (around 150-200) so the Southwestern was not mentioned in this line on purpose.

RC: “Line 430: Supplement or Appendix? ”

AC: “Supplementary analysis” wording will be changed to “additional analysis.”. The numeration of the Appendix figure has changed because there was one more figure included in the Appendix.

“In the additional analysis (Figure A3), annual BB AOD contributions in Warsaw were examined for aggregated European regions (Southern Europe, Eastern Europe, Iberian Peninsula; Figure A3a) and Northern American regions (Canada, USA and Alaska; Figure A3b)”

Thank you for the thorough evaluation of our paper and for the insightful recommendations you offered. Your comments guided us to refine our methodology, clarify key arguments, and improve the manuscript’s coherence. We appreciate the considerable effort you devoted to this review process and value the positive impact your feedback has had on the final version of our study.

RC 1: “Acronyms should be defined at the first place in the manuscript.”

AC: Corrected.

RC 2 : “Section 2.2. It is unclear to me how you used the back trajectory. Did you identify the wildfire events then use the back trajectory simulation to understand the wildfire plume transport?.”

AC:  The initialisation of backward trajectory is done every day at 12 UTC (line 150) regardless of the value of BB AOD. Even if there is no wildfire event (BB AOD is low) our methodology uses BB AOD as a weight (line 206) in mean contribution estimation, so if the BB AOD is close to 0, the result obtained for a certain day won’t alter the final result.

RC 2: “Moreover, would you think just using the average back trajectory to select time and days will overlook some wildfire events.”

AC:  We may not be fully understanding the reviewer’s concern, so to clarify: no single “average” back‑trajectory is used. Instead, for every day in the study period we initialise 14 backward trajectories at 12 UTC at different starting altitudes (see lines 150–152). This daily, multi‑altitude sampling is intended to minimise the risk of overlooking wildfire influence.

Because NAAPS is fed with FLAMBE fire emissions which uses a two‑day maximum correction (previous day + current day) that sustains the fire signal over ~48 h (Lynch et al., 2016; Reid et al., 2009), launching back trajectories once per day (12 UTC) is sufficient for our purposes. Any fire detected in FLAMBE will remain represented in the NAAPS fields over the 24‑hour interval between launches, so significant events influencing BB AOD are unlikely to be missed. Our study targets long‑term, climatological regional contributions; the daily, multi‑altitude (14‑level) trajectory strategy therefore provides adequate coverage while remaining computationally tractable. 

RC 3: “could you justify why you chose 2250 m as your threshold and how a fixed PBL height will affect your results?“

AC: The PBL analysis will be modified such that, instead of using a fixed threshold, the PBL height obtained from HYSPLIT at each point along the air parcel trajectory will be applied. In this approach, emissions will be included only when the air parcel remains below the PBL height at a given location; otherwise, they will be excluded.

The analysis shows that, for European regions, the results averaged over the four dispersion radii obtained with a fixed 2250 m threshold almost replicate the results averaged over the four dispersion radii obtained with the updated PBL method. Differences in mean values between the two approaches are below 2.5%. For individual dispersion radii, the deviations are larger but still below 5%.

The largest differences are observed for contributions from North American regions: the mean (over the four dispersion radii) Canadian contribution decreases by up to ~9%, whereas the USA contribution increases by a similar amount. Deviations for individual radii are of the same order of magnitude as the four-radius mean. The mean Alaska contribution in the PBL method is marginal (maximum 1.2% in July, below 1% in the other months).

RC 4: “I suggest re-labeling figure numbers. For example, Figs. 1a, 2a, 3a-b should be in the same figure since they were discussed at the same time.“

AC:  In our opinion such arrangement of figures is more approachable for the reader: it helps to orient the points chosen for initialisation of forward analysis with fire frequency within a certain region. We’d rather leave it as it is. 

RC 5: “Figures 1c, 2c, 3e-f, 4c–6c. What are the uncertainties?“

AC:  In this study, we did not quantify formal uncertainties for the probabilities shown in Figs. 1c, 2c, 3e–f, and 4c–6c. Our objective was narrowly defined: to determine whether trajectories initiated from a given region have a non-zero probability of reaching Poland. Estimating these uncertainties in a statistically robust way is technically complex and computationally expensive, and would require a dedicated framework (e.g., repeating the calculations with alternative meteorological datasets and with perturbed initial conditions), which lies outside the scope of the present analysis. Therefore, we believe that reporting point estimates without uncertainties is sufficient for the qualitative, yes/no question we set out to answer.

RC 6: “Figures 1-4. Subfigures a and b should have the same coordinates to help readers compare.“

AC:  Coordinates in subfigures a and b will be agreed except for the Figure 3a and 3c where 3c has to cover a bigger longitude area than at 3a.  

RC 7: “L333-335, it is not clear to me how you could get this conclusion. It seems that figure just tell us the possibility of different altitudes of the top plume.“

AC:   This part of the manuscript will be changed to:

“It can be observed that for Canada, the USA, and Alaska (Figures 20d-20f), where the correlation coefficient is the lowest (0.93), the fitted PDF is skewed to the left, which could raise doubts about whether this PDF fit captures fire-plume dynamics well; however, the fitted CDF (Figures A2a–f) yields r = 1 for all regions, and because the contribution analysis uses the CDF, this indicates that the fit is reliable.”