Does the Intra-Arctic Modification of Long-Range Transported Aerosol Affect the Local Radiative Budget? (A Case Study)

The impact of aerosol spatio-temporal variability on the Arctic radiative budget is not fully constrained. This case study focuses on the intra-Arctic modification of long-range transported aerosol and its direct aerosol radiative effect (ARE). Different types of air-borne and ground-based remote sensing observations (from Lidar and sun-photometer) revealed a high tropospheric aerosol transport episode over two parts of the European Arctic in April 2018. By incorporating the derived aerosol optical and microphysical properties into a radiative transfer model, we assessed the ARE over the two locations. Our study displayed that even in neighboring Arctic upper tropospheric levels, aged aerosol was transformed due to the interplay of removal processes (nucleation scavenging and dry deposition) and alteration of the aerosol source regions (northeast Asia and north Europe). Along the intra-Arctic transport, the coarse aerosol mode was depleted and the visible wavelength Lidar ratio (LR) increased significantly (from 15 to 64–82 sr). However, the aerosol modifications were not reflected on the ARE. More specifically, the short-wave (SW) atmospheric column ARE amounted to +4.4 - +4.9 W m−2 over the ice-covered Fram Strait and +4.5 W m−2 over the snow-covered Ny-Ålesund. Over both locations, top-of-atmosphere (TOA) warming was accompanied by surface cooling. These similarities can be attributed to the predominant accumulation mode, which drives the SW radiative budget, as well as to the similar layer altitude, solar geometry, and surface albedo conditions over both locations. However, in the context of retreating sea ice, the ARE may change even along individual transport episodes due to the ice albedo feedback.


Methods: Instrumentation and Modeling Tools
In this study, we utilized active and passive remote sensing observations over two parts of the European Arctic in order to retrieve the optical and microphysical properties of long-range transported aerosol. The ensemble of inverted aerosol properties was incorporated into a radiative transfer model to assess the ARE and its related uncertainties. In parallel, the simulated irradiances were evaluated with ground-based and air-borne observations.

Lidar Systems
The elastic air-borne mobile aerosol Lidar (AMALi) is installed on board of the research aircraft Polar 5 [38]. During the polar air-borne measurements and Arctic Regional Climate Model Simulation Project 2018 (PAMARCMiP2018), AMALi was operated in zenith configuration during 9 flights  April for approximately 32 h in total) over the region of Fram Strait (84,17 NW to 79,20 NE). AMALi is equipped with an Nd:YAG laser that emits pulses at 15 Hz. A 10-cm off-axis telescope mirror of 3.1 mrad field of view (FOV) collects the backscattered light. Photomultiplier tubes (PMTs) are used for the detection of elastically backscattered photons at 532 and 355 nm. Additionally, the cross-polarized component is recorded at 532 nm, allowing for the determination of the aerosol linear depolarization ratio, hereafter mentioned as the aerosol depolarization ratio.
AMALi signal returns were acquired with a 1 s--7.5 m resolution. Signals were vertically averaged every 30 m and smoothed with a 90-m window, while backscatter profiles were smoothed with a 150-m window. The altitude of the observed features was corrected with respect to the aircraft movements, with the pitch and roll angles changing by 2 • and 3.5 • during the evaluation windows. For this study, we extracted two observation periods of a 2-and 5-min duration (as indicated in Figures  1 and 2a), which were limited by the changing flight altitude. The signal-to-noise ratio (SNR) changed considerably due to the variable flight altitude (0.3-2.9 km) during AMALi observations, which limited the temporal averaging and led to higher uncertainties in the air-borne-derived optical properties compared to the ground-based ones.
The multi-wavelength Koldewey Aerosol Raman Lidar (KARL) is operated in the Alfred Wegener Institute Emille Victor (AWIPEV) research base, Ny-Ålesund, Svalbard. KARL is a "3β+2α+2δ+2wv" system equipped with an Nd:YAG laser emitting at 1064, 532, and 355 nm after frequency doubling and tripling, respectively [39]. The pulse repetition rate of the system is 50 Hz, emitting approximately 10 W in each wavelength. The receiver consists of a parabolic 70 cm diameter Newtonian telescope with a 1.75 mrad FOV and an aperture with an adaptable position and pinhole diameter. PMTs are used for the detection of backscattered photons in all channels, except for 1064 nm, where avalanche photodiodes (APDs) are used. KARL signal returns were obtained with a 1.5 min-7.5 m resolution. In order to reduce signal noise, we averaged the profiles for approximately 2.5 h and applied vertical smoothing of 157.5 m in the backscatter retrieval and 225 m in extinction retrieval, following the effective resolution concept [40].
We applied the Klett-Fernald method to the elastic AMALi and KARL signals in order to retrieve the backscatter profiles [41,42]. This method requires a priori information on the vertical profile of the extinction-to-backscatter ratio (Lidar ratio or LR) and a reference value for the backscatter ratio (total backscatter over molecular backscatter). Here, we assumed uncertainties of 10 sr and 3% for the LR and backscatter ratio, respectively. The positioning of the calibration window introduces additional uncertainty. Errors are also introduced due to statistical signal noise and the procedures of temporal averaging and vertical smoothing. We obtained extinction profiles from KARL by utilizing the Raman method [43]. Since Raman scattering cross-sections are orders of magnitude lower than Rayleigh cross-sections and our observations were performed during the Polar day period, we derived layer-mean values of extinction and LR. Aerosol optical properties are presented in Tables 1 and 2.

Sun-Photometers and Extinction Estimation at Fram Strait
The air-borne Sun-photometer with an active tracking system (SPTA) was operated under a quartz dome on the top of Polar 5 during PAMARCMiP2018. The aerosol optical depth (AOD) was provided at 10 wavelengths between 396 and 1026 nm every 30 sec. Unfortunately, the 369-414 nm channels presented calibration issues and, thus, they were not used for further evaluation. Moreover, we did not use the 946-nm channel, which is dedicated to water vapor (WV) absorption. On 2 April, the AOD was measured at different flight levels between 0.04 and 5.15 km (Figure 3). We performed 5-min averaging in the AOD in order to reduce statistical noise. The aerosol extinction coefficient was estimated by Remote Sens. 2020, 12, 2112 6 of 24 consecutive AOD pairs within altitude layers of 0.5-2.2 km and temporal periods of 5-40 min, within which the atmospheric conditions were assumed to be horizontally homogeneous (similar AOD at similar altitudes, Figure 3). km apart). The AOD was the product of the integrated backscatter with the LR. The KARL AOD was interpolated to the photometer wavelengths using an appropriate Angström exponent (1.8 as fitted to the backscatter at 355, 532, and 1064 nm). Following the same steps, we estimated the stratospheric AOD. Finally, the extinction of the elevated aerosol layer was estimated by Equation (2): . .
. . The AOD of the elevated aerosol layer (5.15-6.8 km), observed close to the Greenland side of Fram Strait (Figures 1 and 2a), was estimated by combining photometer and Lidar observations. More specifically, we used the photometer AOD at five different measurement points below the layer, between the highest photometer observation (number 5, 10:52-10:57 UT; flight level 5.15 km) and the closest AMALi observation (number 1, 10:31-10:33 UT; flight level 2.8 km). We also used the AOD at three intermediate flight levels (numbers 2, 3, and 4 corresponding to 3, 4.3, and 4.4 km) so as to better constrain the estimated extinction. The AOD was calculated according to Equation (1), with the selected measurement points enumerated in Figure 3. Horizontal homogeneity was assumed due to the short spatial (maximum 100 km) and temporal differences (maximum 25 min) among the measurement points, as in the study of Mei et al. [44] that focused on satellite aerosol remote sensing. Above the aerosol layer, the AOD was derived from KARL. The ground-based Lidar signals were preferred over the air-borne ones due to their higher SNR. The aerosol load in the stratosphere was invariable between the air-borne and ground-based observations (three days and a maximum of 500 km apart). The AOD was the product of the integrated backscatter with the LR. The KARL AOD was interpolated to the photometer wavelengths using an appropriate Angström exponent (1.8 as fitted to the backscatter at 355, 532, and 1064 nm). Following the same steps, we estimated the stratospheric AOD. Finally, the extinction of the elevated aerosol layer was estimated by Equation (2): Over Ny-Ålesund the AOD was measured by an SP1a type sun-photometer in 10 wavelengths between 369 and 1023 nm. A cloud screening procedure that relies on the short-term AOD variability was applied to the data [45]. The total AOD uncertainty amounted to 0.01 [32]. In this study, we only utilized the wavelengths, which are spectrally closer to the air-borne photometer channels, except for the 1023 nm that was affected by calibration issues.

Aerosol Microphysical Properties Inversion Schemes
We incorporated the set of optical parameters derived for Fram Strait and Ny-Ålesund into inversion algorithms so as to retrieve the aerosol microphysical properties. In the inversion schemes, the extinction and backscatter coefficients were used as input for the derivation of the complex refractive index and aerosol volume size distribution using Mie theory. For each combination of the real and imaginary refractive index, the best solution was identified. More specifically, the least residuum between the backward (input) and forward calculated optical parameters was selected. The standard deviation of the best solutions represents the uncertainties of the complex refractive index. Statistical noise and systematic uncertainties of the input optical parameters together with mathematical approximation errors contribute additively to the total uncertainties in the microphysical parameters.
For Fram Strait, we utilized the AMALi-derived backscatter coefficient (355 and 532 nm) together with the sun-photometer-derived extinction (496, 675, 779, 861, and 1026 nm). We employed collocation discretization and an iterative Padé-regularization method with a fixed number of 30 iterations [35,36]. For Ny-Ålesund, the KARL-derived aerosol backscatter (355, 532 and 1064 nm) and extinction coefficients (355 and 532 nm) were inverted to the particle microphysics employing a truncated singular value decomposition regularization algorithm [34]. More details about the Lidar inversion algorithms, which have been developed at the University of Potsdam, can be found in [34,37,46].
The total number (n t ), surface-area (s t ), and volume (v t ) concentration, as well as the effective radius (r eff ) were derived for the fine and coarse aerosol modes by performing log-normal fits to the inverted volume distribution. Fitting errors amounted to 10% for the total concentrations and 5% for the r eff estimation. The single-scattering albedo (SSA) and the asymmetry parameter were calculated from the fitted number size distributions utilizing Mie theory. The corresponding SSA and asymmetry parameter uncertainties were estimated by sensitivities based on different complex refractive index scenarios.

Trajectory Calculations for Air-Mass History
In order to assay the origin of the observed transport event and link the evolution of aerosol properties with intruding air masses, we performed 10-day LAGRANTO backward trajectories over Fram Strait and Ny-Ålesund [47,48]. A significant number of trajectories ending over Fram Strait and Ny-Ålesund originated from the lower troposphere, providing an indication for the surface aerosol sources. Due to the complex orography of Spitsbergen island, with the highest peak at 1.7 km, and the challenging to parameterize stable boundary layer, we are inclined to have higher confidence in higher altitude trajectories.
We used ECMWF analysis data to calculate trajectories with LAGRANTO, which utilizes the three-dimensional wind field to calculate kinematic trajectories. The ECMWF data for this study had a horizontal grid spacing of 0.5 • x 0.5 • with 137 hybrid sigma pressure levels in total. Analysis fields were available every six hours. For the flight locations over Fram Strait, we initialized the trajectories in a small region around the flight track. Trajectories from Ny-Ålesund were initialized every 0.05 • in the horizontal within a 0.5 • × 0.5 • box centered at Ny-Ålesund. Trajectory information was then available every hour. Ancillary meteorological parameters, such as pressure and relative humidity (RH), were calculated along the trajectories. The spatial variability of the trajectories was taken into account by using a box, instead of a single point, for initializing our calculations. Moreover, the influence of atmospheric variability was considered by analyzing coherent bundles of trajectories instead of single trajectories.

Radiation and Meteorological Observations
We utilized air-borne and ground-based irradiance observations for evaluating the modeled irradiances at the surface and different flight altitudes. Radiation sensors were installed on the top and bottom of the Polar 5 fuselage during PAMARCMiP2018. More specifically, a CMP22 Pyranometer by Kipp&Zonen was measuring broadband (0.2-3.6 µm) global and reflected solar irradiances. The measurement repetition rate was 20 Hz, with all the irradiances being deconvoluted for their inertia time [49]. The overall uncertainty was lower than 10 W m −2 for all radiation quantities. Basic meteorological observations were also performed on board of the Polar 5 aircraft [50]. WV concentration was measured by the LI-7200 closed CO 2 /H 2 O analyzer [51].
At Ny-Ålesund, radiation measurements are performed within the frame of the Baseline Surface Radiation Network (BSRN) since 1992 [52]. Short-wave (SW) broadband radiation in the range 0.2 to 3.6 µm is measured with CMP22 pyranometers by Kipp&Zonen, installed in an up-and downward orientation for global and reflected radiation, respectively. In addition, the diffuse radiation is obtained by ball-shaded similar instrumentation. The basic surface radiation and meteorological measurements with a one-minute resolution applied in our study are available via the Pangaea data repository [53]. Vertical profiles of the temperature, RH, wind speed, and direction in the upper atmosphere are applied from daily radiosonde ascents performed at the AWIPEV station [54].

Radiative Transfer Simulations
In order to quantify the altitude-dependent direct ARE and its uncertainties, we performed simulations with the radiative transfer model SCIATRAN [55], using as input the retrieved aerosol optical and microphysical properties. For the evaluation of the simulated irradiances, the latter were integrated within the spectral range of the BSRN and air-borne pyranometers.
We performed simulations from the local surface up to 40 km. The vertical resolution varied from 50 m to 1 km (above 10 km). The radiative transfer equation was solved by means of the scalar discrete ordinate technique, while a plane-parallel atmosphere approximation was used. The absorption contribution of line absorbers was computed in line-by-line mode of 0.01-nm spectral resolution, based on the spectroscopic parameters of the HITRAN 2008 database [56]. The altitude-dependent ARE was derived from the difference of SW net irradiances between the respective aerosol and the off-aerosol scenario, in order to distinguish it from the total atmospheric effect. The difference between the TOA and surface ARE yielded the overall atmospheric column ARE.
The aerosol-related input comprised the spectrally resolved aerosol extinction coefficient (Table 2), the SSA, as well as the asymmetry factor calculated according to the Henyey-Greenstein parameterization ( Table 3). The aerosol optical and microphysical properties were extrapolated to the model's spectral range (from 355-1026 nm) by forward Mie theory. More specifically, using the aerosol refractive index and number size distribution as input, we were able to calculate the aerosol extinction coefficient at any wavelength. The surface albedo over Fram Strait and Ny-Ålesund was obtained from air-borne and BSRN pyranometer observations, respectively. A surface albedo of approximately 0.7 was found for both locations, revealing the presence of sea ice at Fram Strait and snow-covered tundra in Ny-Ålesund.
For Fram Strait, we performed four simulation scenarios: Aerosol and off-aerosol for the aircraft positions during the AMALi and sun-photometer evaluation periods, as indicated in Figure 1. Thermodynamic profiles were constructed as a combination of appropriately selected air-borne meteorological observations below 5 km (pressure and temperature from aircraft meteorological observations, WV mixing ratio from LI-7200 instrument) and radiosonde ascents from Ny-Ålesund aloft (approximately 500 km distance from flight operations). Between 5 and 8 km, a hybrid air-borneradiosonde profile was constructed, which gradually shifted from the highest air-borne towards the radiosonde observations.
For Ny-Ålesund, simulations were performed for the period of available Raman Lidar observations (5 April, 11:00-13:30). The simulations comprised four different scenarios: Off-aerosol, with two aerosol Remote Sens. 2020, 12, 2112 9 of 24 layers (4.3-5.75 km), with the lower layer (4.3-4.9 km) and finally with the upper layer (4.9-5.75 km). Ancillary meteorological information on the temperature at the surface, 2 and 10 m, as well as surface pressure, were also incorporated into the model. Pressure, temperature, WV mixing ratio, and O 3 concentration profiles were taken from a collocated ozonesonde ascent (5 April, AWIPEV). For both Arctic locations, further trace gas profiles were obtained from climatological profiles for the given month and latitude [57]. Table 3. Aerosol microphysical properties over the Greenland side of Fram Strait and Ny-Ålesund. The complex refractive index ri, the single-scattering albedo SSA, and the asymmetry parameter g are given. Moreover, the effective radius r eff as well as the number n t , surface s t , and volume v t concentrations are presented for the fine and coarse aerosol modes. For details on the retrieval schemes, see Section 2.3.

Aerosol Observations over the European Arctic and Meteorological Conditions
Air-borne observations were performed within the frame of PAMARCMiP2018 (March-April), Villum Research Station, Greenland (81.6 • N, 16.7 • W). The PAMARCMiP2018 period was characterized by a high-pressure system over the North Pole, while some weak lows occurred over northeast Greenland. On 2 April, a well-defined aerosol layer was identified at a remarkably high tropospheric altitude over Fram Strait by AMALi and sun-photometer (Figures 1 and 2a). Three days later, a geometrically similar aerosol layer was observed by KARL over Ny-Ålesund (Figure 2b), at a distance of 200-600 km from the flight operations. In the interim of Fram Strait and Ny-Ålesund on 3 April, clouds and aerosol remnants were captured in the upper troposphere by AMALi ( Figure S1). At that day, Raman Lidar observations over Ny-Ålesund were not allowable due to low-level clouds. However, the micro-pulse Lidar (MPL, technical description given in Shibata et al. [18]) provided the geometrical evolution of the aerosol layers. It should be noted that the MPL optical products cannot be utilized for the derivation of aerosol microphysical properties due to the lack of Raman channels. Therefore, the optical and microphysical aerosol properties over Ny-Ålesund were investigated within periods of 5 and 6 April as indicated in Figure 2b. On the days when the aerosol layer was observed over Ny-Ålesund, the synoptic situation was characterized by the passing of a weak low-pressure system. According to radiosonde ascents from Ny-Ålesund, moderate N-NW winds (3-7 m/s) prevailed, while wind shear and RH were enhanced (RH of 47-58% on 5 April and 55-65% on 6 April) at the aerosol layer atmospheric levels.

Episode Evolution over Ny-Ålesund and Connection to Air Masses
The geometrical evolution of the transported aerosol plume over Ny-Ålesund along with the related air mass trajectories are presented in Figure 4. The columnar AOD over Ny-Ålesund and Fram Strait are demonstrated in Figure 5. Before the layer arrival (30)(31), air masses originated from central Eurasia (4-7 km, Figure 4b). The daily mean AOD 500 was relatively low (0.04 ± 0.01) with respect to background values (AOD 500 = 0.06 from Yamanouchi et al. [58]). The influence of northeast Asian air masses is discernible for 2 April, with aerosol features appearing entangled with clouds. Due to persistent low-level clouds (2-4 April), Raman Lidar and sun-photometer observations were unfortunately not allowable during those days.
On 5 April, the layer displayed a well-defined structure. The AOD 500 followed a moderate increasing tendency with persisting northeast Asian and emerging north European air masses ( Figure 4d). On 6 April, the layer was vertically perturbed (Figure 4c,6 April 0-8 UT), possibly due to the passing of the weak low-pressure system. The AOD 500 fluctuated and maximized (AOD 500 = 0.13) some hours later. The layer gradually diluted (7 April), with the AOD approaching background levels. Re-emerging low-level clouds did not allow us to further investigate the event evolution, with intermittent aerosol features and mid-level clouds mostly observed (8 April). Overall, the AOD over Ny-Ålesund was clearly lower compared to Fram Strait. In Sections 4.3 and 4.4, we will investigate in more detail the modification of aerosol optical and microphysical properties within the elevated aerosol plume.

Origin of Observed Aerosol Over Fram Strait and Ny-Ålesund
Prior to investigating the modification of aerosol optical and microphysical properties, we assessed whether the observed aerosol plumes over Fram Strait and Ny-Ålesund originated from the same region or regions comprising similar aerosol sources. The atmospheric pressure along the LAGRANTO backward trajectories as well as the air mass transport time are displayed in Figure 6.
At Fram Strait, the layer close to Greenland (Figure 6a) was a receptor of central and northeast Asian air masses, which were at least nine days old (time since contact with main aerosol source region, Figure 6c). At Ny-Ålesund, the lower part of the plume originated from north Europe (Figure 6d), while the upper part was also affected by northeast Asian air masses (Figure 6e). On their arrival at Ny-Ålesund, the Asian air masses were at least nine days old, while the north European ones had last contact with the main aerosol sources four days ago (Figure 6f). The latter air masses mixed over Fram Strait (500-600 hPa) in the earlier days (Figure 6f).   For comparison, the AOD over Fram Strait as derived from the air-borne sun-photometer is given (average ± one standard deviation). Overall, the AOD over Ny-Ålesund was significantly lower than Fram Strait, reflecting the effect of removal processes and modified aerosol source efficiency.
On 5 April, the layer displayed a well-defined structure. The AOD500 followed a moderate increasing tendency with persisting northeast Asian and emerging north European air masses ( Figure  4d). On 6 April, the layer was vertically perturbed (Figure 4c,6 April 0-8 UT), possibly due to the passing of the weak low-pressure system. The AOD500 fluctuated and maximized (AOD500 = 0.13) some hours later. The layer gradually diluted (7 April), with the AOD approaching background levels. Re-emerging low-level clouds did not allow us to further investigate the event evolution, with intermittent aerosol features and mid-level clouds mostly observed (8 April). Overall, the AOD over Ny-Ålesund was clearly lower compared to Fram Strait. In Sections 4.3 and 4.4, we will investigate in more detail the modification of aerosol optical and microphysical properties within the elevated aerosol plume.

Origin of Observed Aerosol Over Fram Strait and Ny-Ålesund
Prior to investigating the modification of aerosol optical and microphysical properties, we assessed whether the observed aerosol plumes over Fram Strait and Ny-Ålesund originated from the same region or regions comprising similar aerosol sources. The atmospheric pressure along the LAGRANTO backward trajectories as well as the air mass transport time are displayed in Figure 6. At Fram Strait, the layer close to Greenland (Figure 6a) was a receptor of central and northeast Asian air masses, which were at least nine days old (time since contact with main aerosol source region, Figure 6c). At Ny-Ålesund, the lower part of the plume originated from north Europe (Figure  6d), while the upper part was also affected by northeast Asian air masses (Figure 6e). On their arrival at Ny-Ålesund, the Asian air masses were at least nine days old, while the north European ones had  pressure (a, b, d, and e) and transport time (c and f) of the trajectories are color indicated. The transport time corresponds to air masses arriving at the Greenland side of Fram Strait as well as to both aerosol layer altitude sectors at Ny-Alesund (c and f). Aged Asian air masses affected Fram Strait. Ny-Ålesund was influenced by Asian and fresher air masses from north Europe. The latter air masses mixed over Fram Strait two to three days earlier (f).

Modification of Aerosol Optical Properties
The optical properties of the long-range transported aerosol over Fram Strait and Ny-Ålesund are summarized in Tables 1 and 2, for the evaluation periods indicated in Figures 1 and 2. The vertical profiles of the aerosol backscatter coefficient (proportional to the aerosol number concentration), aerosol depolarization ratio (rising with aerosol non-sphericity), and backscatter-related Angström exponent (inversely proportional to aerosol radius) are demonstrated in Figure 7. The aerosol extinction coefficient and extinction-related Angström exponent as derived by sun-photometer and Lidar are presented in Figure 8.
A shift towards smaller particles was observed, mainly between Fram Strait and Ny-Ålesund, as well as between 5 and 6 April over Ny-Ålesund. Over Fram Strait, the layers were characterized by low but variable aerosol depolarization ratios (1.1-5.5%). This indicates nearly spherical particles as well as connotes the presence of different aerosol mixtures. The lower depolarization close to Greenland (Figure 7b   Aerosol optical properties over Fram Strait (Fram) and Ny-Ålesund (NyA) derived from the air-borne AMALi and ground-based KARL systems, respectively. The aerosol backscatter coefficient β, aerosol depolarization ratio δ, and backscatter-related Angström exponent Αβ are presented. Error bars represent retrieval uncertainties (for details, see Section 2.1). Air-borne profiles exhibit higher uncertainties due to flight altitude changes that constrained the temporal averaging. The aerosol load was higher over Fram Strait (higher β, a). For both locations, particles were nearly spherical (low δ, b). A shift towards smaller particles was observed over Ny-Ålesund (higher Αβ, c).

Figure 7. Aerosol optical properties over Fram Strait (Fram) and
Ny-Ålesund (NyA) derived from the air-borne AMALi and ground-based KARL systems, respectively. The aerosol backscatter coefficient β, aerosol depolarization ratio δ, and backscatter-related Angström exponent A β are presented. Error bars represent retrieval uncertainties (for details, see Section 2.1). Air-borne profiles exhibit higher uncertainties due to flight altitude changes that constrained the temporal averaging. The aerosol load was higher over Fram Strait (higher β, a). For both locations, particles were nearly spherical (low δ, b). A shift towards smaller particles was observed over Ny-Ålesund (higher A β , c). A shift towards smaller particles was observed, mainly between Fram Strait and Ny-Ålesund, as well as between 5 and 6 April over Ny-Ålesund. Over Fram Strait, the layers were characterized by low but variable aerosol depolarization ratios (1.1-5.5%). This indicates nearly spherical particles as well as connotes the presence of different aerosol mixtures. The lower depolarization close to Greenland (Figure 7b For Ny-Ålesund, a mixture of industrial pollution and biomass-burning aerosol is more likely for the upper part of the plume. The lower layer is more likely to comprise industrial pollution particles as it is affected by north European intrusions. The depolarization (approximately 3%) indicates more spherical particles compared to mid-Fram Strait but less spherical than the Greenland regime. The transition to more spherical aerosol could be attributed to the depletion of non-spherical components. In general, particles with a higher surface roughness are more likely to undergo scavenging due to their higher heterogeneous nucleation efficacy [59]. Close to Greenland, the LR532 was significantly lower (15 sr) compared to Ny-Ålesund (64-82 sr). We will discuss in more detail the modification of light scattering properties in Section 5. For Ny-Ålesund, a mixture of industrial pollution and biomass-burning aerosol is more likely for the upper part of the plume. The lower layer is more likely to comprise industrial pollution particles as it is affected by north European intrusions. The depolarization (approximately 3%) indicates more spherical particles compared to mid-Fram Strait but less spherical than the Greenland regime. The transition to more spherical aerosol could be attributed to the depletion of non-spherical components. In general, particles with a higher surface roughness are more likely to undergo scavenging due to their higher heterogeneous nucleation efficacy [59]. Close to Greenland, the LR 532 was significantly lower (15 sr) compared to Ny-Ålesund (64-82 sr). We will discuss in more detail the modification of light scattering properties in Section 5.

Modification of Aerosol Microphysical Properties
The inverted aerosol microphysical properties for Fram Strait and Ny-Ålesund are summarized in Table 3, while the corresponding volume size distributions are illustrated in Figure 9. The microphysical retrieval periods are denoted in Figure 2. The low imaginary refractive index at the Greenland side of Fram Strait suggests the domination of scattering particles. The same holds true for the upper layer over Ny-Ålesund, while weakly absorbing particles dominated in the lower one, respectively. It should be noted that we considered the refractive index as wavelength independent, which is a common assumption. However, the wavelength-dependent SSA indicates the presence of relatively more absorbing small particles over Fram Strait (SSA 355 = 0.84) as well as in the lower Ny-Ålesund plume (SSA 355 = 0.88). The asymmetry parameter was similar over Fram Strait and Ny-Ålesund (approximately 0.7), indicating relatively large particles. Moreover, the asymmetry parameter was higher than the monthly average of different Arctic sites, which mainly lies between 0.55 and 0.65 (Figure 7 from Schmeisser et al. [10]), with such high values usually being pertinent to long-range transport events [10].   The aerosol size distribution was clearly modified along the intra-Arctic transport (Figure 9). A bi-modal distribution was derived for Fram Strait (fine r eff = 0.49 µm, coarse r eff = 1.26 µm), whereas the coarse mode was strongly diminished over Ny-Ålesund. A weak bi-modal distribution (fine r eff = 0.26 µm, coarse r eff = 1.43 µm) emerged in the lower Ny-Ålesund layer, while the fine mode (r eff = 0.19 µm) was predominant in the upper one layer. This transition could be related to the modified pollution pathways, as presented in Section 4.2, as well as to aerosol removal processes. The latter processes will be further discussed in Section 5.

Aerosol Radiative Effect (ARE)
We assessed the direct effect of aerosol on the local radiative budget by incorporating the retrieved aerosol properties into the radiative transfer model SCIATRAN. The vertical distribution of simulated SW downward, upward and net irradiances, as well as the ARE are depicted in Figure 10. For evaluation, the observed irradiances at the surface and different flight levels are displayed.

Aerosol Radiative Effect (ARE)
We assessed the direct effect of aerosol on the local radiative budget by incorporating the retrieved aerosol properties into the radiative transfer model SCIATRAN. The vertical distribution of simulated SW downward, upward and net irradiances, as well as the ARE are depicted in Figure 10. For evaluation, the observed irradiances at the surface and different flight levels are displayed.
The aerosol layer attenuated the down-welling irradiance through scattering and absorption. At Fram Strait, the surface up-welling irradiance was diminished to a less pronounced degree, in terms of the absolute magnitude compared to the downward component, and, thus, a "shadowing" effect (−3.5 to −2.9 W m -2 ) was produced. The less intense depletion of upward surface irradiances connotes that their main driver is the surface albedo. Contrarily, the upward component above the layer was diminished more intensely, compared to the downward component, producing a warming effect at TOA (+0.9 to +2 W m -2 ). Finally, +4.4 to +4.9 W m -2 were retained in the atmosphere, with the majority absorbed within the aerosol layer. Compared to air-borne radiation observations, the simulated irradiances presented a maximum bias of 3%, 8%, and 13 % in the SW down-welling, up-welling, and net components, respectively. These differences can be attributed to uncertainties in the thermodynamic and trace gas profiles, as well as the high solar zenith angles (larger than 75⁰), which could introduce additional errors in the radiative transfer calculations, especially under the plane-parallel assumption. Figure 10. Altitude-dependent SW simulated irradiances and aerosol radiative effect (ARE) for AMALi (1) and sun-photometer (2) evaluation periods over Fram Strait (a) and over Ny-Ålesund (b). Observed irradiances are given for comparison (average ± one standard deviation). The aerosol layer altitude is indicated with grey shading. Over both locations, TOA warming was accompanied by surface cooling. Note the different scales.
At Ny-Ålesund, the TOA was characterized by positive ARE (two-layers +3 W m -2 , lower layer +1.15 W m -2 , upper layer +0.1 W m -2 ). More pronounced warming was produced within the atmospheric column (two-layers +4.5 W m -2 , lower layer +1.8 W m -2 , upper layer +0.6 W m -2 ), which led to a cooling effect at the surface (two-layers −1.5 W m -2 , lower layer −0.6 W m -2 , upper layer −0.5 W m -2 ). In comparison to BSRN observations, the simulated full atmospheric scenario (two-layers) irradiances exhibited agreement with the range of uncertainty. It is worth noticing that the ARE of the two Figure 10. Altitude-dependent SW simulated irradiances and aerosol radiative effect (ARE) for AMALi (1) and sun-photometer (2) evaluation periods over Fram Strait (a) and over Ny-Ålesund (b). Observed irradiances are given for comparison (average ± one standard deviation). The aerosol layer altitude is indicated with grey shading. Over both locations, TOA warming was accompanied by surface cooling. Note the different scales.
The aerosol layer attenuated the down-welling irradiance through scattering and absorption. At Fram Strait, the surface up-welling irradiance was diminished to a less pronounced degree, in terms of the absolute magnitude compared to the downward component, and, thus, a "shadowing" effect (−3.5 to −2.9 W m −2 ) was produced. The less intense depletion of upward surface irradiances connotes that their main driver is the surface albedo. Contrarily, the upward component above the layer was diminished more intensely, compared to the downward component, producing a warming effect at TOA (+0.9 to +2 W m −2 ). Finally, +4.4 to +4.9 W m −2 were retained in the atmosphere, with the majority absorbed within the aerosol layer. Compared to air-borne radiation observations, the simulated irradiances presented a maximum bias of 3%, 8%, and 13 % in the SW down-welling, up-welling, and net components, respectively. These differences can be attributed to uncertainties in the thermodynamic and trace gas profiles, as well as the high solar zenith angles (larger than 75 0 ), which could introduce additional errors in the radiative transfer calculations, especially under the plane-parallel assumption.
At Ny-Ålesund, the TOA was characterized by positive ARE (two-layers +3 W m −2 , lower layer +1.15 W m −2 , upper layer +0.1 W m −2 ). More pronounced warming was produced within the atmospheric column (two-layers +4.5 W m −2 , lower layer +1.8 W m −2 , upper layer +0.6 W m −2 ), which led to a cooling effect at the surface (two-layers −1.5 W m −2 , lower layer −0.6 W m −2 , upper layer −0.5 W m −2 ). In comparison to BSRN observations, the simulated full atmospheric scenario (two-layers) irradiances exhibited agreement with the range of uncertainty. It is worth noticing that the ARE of the two layers did not act cumulatively, indicating the complexity of aerosol radiation interactions. The positive TOA ARE revealed at both Arctic locations implies a contribution from black carbon. On an annual basis, fossil fuel-related black carbon drives a positive TOA forcing (+0.19 W m −2 ) with respect to pre-industrial emissions, albeit the black carbon AOD is low compared to other aerosol species [3].

ARE Uncertainties and Comparison to WV-Related Radiative Effect
To investigate the effect of aerosol uncertainties on the simulated ARE, we performed a set of sensitivity tests. More specifically, we repeated the full atmospheric scenario over Ny-Ålesund, but we perturbed the aerosol input by its highest related uncertainties, as they are presented in Tables 2 and 3. We found out that a 30% increase in the aerosol extinction produced more pronounced warming at TOA (by +0.9 W m −2 ) and atmosphere (by +1.3 W m −2 ) and, thus, moderately enhanced cooling at the surface (by -0.4 W m −2 ). Accordingly, a 10% decrease in the SSA resulted in even more pronounced TOA and atmospheric warming (by +1.7 and +2.1 W m −2 ) but similar surface cooling (by −0.4 W m −2 ) as the extinction perturbation. Finally, an 8% decline in the asymmetry parameter led to equal magnitude warming and cooling at TOA and the surface (by −0.1 W m −2 ), respectively, and, thus, the atmospheric radiative budget was not modified. According to error propagation, the maximum cumulative uncertainty amounted to 2.7 W m −2 at TOA, 0.9 W m −2 at the surface, and 3.4 W m −2 within the atmospheric column. Consequently, the sign of the obtained ARE is not affected given the range of aerosol related uncertainties.
Finally, we investigated the effect of the WV concentration on the irradiances in order to assess its radiative effect comparatively to the aerosol effect. More specifically, we repeated the full atmospheric scenario over Ny-Ålesund, but we reduced the WV mixing ratio by 30%. Our sensitivity analysis revealed that a 30% decline in the WV mixing ratio resulted in a −3.8 W m −2 atmospheric effect. At the surface and TOA, the effect amounted to +1 and −2.8 W m −2 , respectively. Hence, the atmospheric column ARE (+4.5 W m −2 ) outweighed the WV radiative effect (−3.8 W m −2 ), with maximum weekly WV fluctuations over Ny-Ålesund (2-9 April) being lower than 30% at 5-15 km.

Discussion
An unusual aspect of the investigated long-range aerosol transport is its high altitude, with the majority of reported episodes over Ny-Ålesund confined below 5 km [14,24,60,61]. Aerosol climate models predict such transport with an increased black carbon burden at atmospheric levels higher than 500 hPa [62][63][64], but this is usually attributed to an underestimation of wet removal. With our study, however, we found observational evidence of the rare but persistent presence of such elevated aerosol transport [14,24,60,61]. North European sources affect, in general, lower altitudes due to their relative proximity to the Arctic [16]. The intrusion of north Eurasian air masses in early spring is promoted either by blocking patterns or strong diabatic cooling along transport over snow-covered areas. Thereby, aerosol trapped below temperature inversions can be transported over long distances [16][17][18]. Asian origin transport at such high tropospheric altitudes is seldom observed [18]. Lifting processes, nonetheless, that facilitate the export of air masses from the boundary layer to these altitudes are related to the warm conveyor belts of mid-latitude cyclones. Air masses ascend along the sloping isentropes and can penetrate the polar dome from above through radiative cooling [16][17][18]. In the future, lofted aerosol layers may appear more frequently over the Arctic in the context of projected intensification in strong extra-tropical surface cyclones [65,66]. Moreover, in view of the increasing long-term trend in blockings over the west Pacific [67], more frequent elevated aerosol layers can play an important role in high cloud formation processes.
The intra-Arctic aerosol modification of this episode demonstrates not only the effect of altered aerosol source regions but also the role of removal processes. Mechanisms promoting the aerosol depletion within pollution plumes are mixing with ambient air, coagulation, nucleation scavenging (activation to cloud condensation nuclei and ice nucleating particles as well as subsequent scavenging of interstitial aerosol), dry deposition (impaction, gravitational settling, Brownian motion), and wet scavenging [68]. In this transport episode, a transition towards smaller accumulation mode particles along with the depletion of the coarse mode occurred (Figure 9). The bigger aerosol is deposited more effectively due to gravitational settling and impaction is more efficient due to their higher inertia [69] (pp. 362-395). However, solely dry deposition could not account for the coarse mode depletion within a period of three days. The presence of clouds over Fram Strait ( Figure S1) and Ny-Ålesund (Figure 4a) in the interim of the episode indicates that the transported aerosol underwent nucleation scavenging. Nucleation rates are proportional to the aerosol surface area [70].
The Fram Strait plume and its Ny-Ålesund counterpart consisted of aerosol with a similar refractive index. Biomass burning aerosol usually exhibits lower absorption along its transport, which is associated with the degree of oxidative processing on the particle surface [71]. However, in the course of this episode, a part of the aerosol plume (lower Ny-Ålesund sub-layer) was more absorbing, possibly due to interstitial and additionally supplied black carbon particles from neighboring industrial pollution sources (north Europe). Overall, the refractive indices are in line with reported Arctic haze events over Ny-Ålesund, which comprised either polluted ([72], 1.56 + 0.01 i) or mixtures of industrial pollution and biomass-burning aerosol from Siberia ([24], 1.6 + 0.007 i). The aforementioned Arctic haze events demonstrated a higher total number concentration ( [72], 343 ± 30% cm −3 ; [24], 480-950 cm −3 ), probably due to the more efficient aerosol sources or less effective removal processes along their shorter transport paths (central and east Siberia). The accumulation mode at Ny-Ålesund presented similarities with those events ( [72], r eff = 0.18 µm; [24], r eff = 0.19 µm). Average Arctic aerosol size distributions (derived from in situ measurements with a 0.5-µm cut-off diameter) peak around radii of 0.15 µm [8,9,73]. However, the accumulation mode particles at Fram Strait were untypically large. Bi-modal size distributions are also typical for the Arctic haze season of west Greenland [74], but their accumulation mode is shifted to smaller radii than the present episode.
Both biomass-burning and industrial pollution aerosol constitute combustion particles. The LR spectral dependency (higher LR 532 than LR 355 ) at Ny-Ålesund is more similar to that of aged biomass-burning aerosol [71,75,76]. Overall, biomass-burning aerosol demonstrates diverse LR (LR 355 = 40-89 sr and LR 532 = 46-100 sr; [71,[75][76][77][78][79]) due to the interplay of several parameters, such as vegetation type, transport time, combustion efficiency (temperature and oxygen availability), local meteorological conditions (RH), as well as mixing state of the particles and content of elemental and organic carbon. Similar factors formulate the light scattering properties of industrial pollution particles. These particles display higher LR in the ultraviolet spectral region (LR 355 = 41-73 sr and LR 532 = 23-74 sr, [77,79,80] but generally exhibit lower r eff than ours. The LR 532 (15 sr) at Fram Strait was lower compared to the reported literature values for biomass-burning and industrial pollution particles. From a microphysical perspective, the particle size distribution exhibits less sensitivity to the visible light, since this is optically dominated by the coarse mode. As illustrated in Figure S3 the coarse mode particles dominate the extinction and backscatter cross-sections. By performing Mie calculations, we estimated that the highest LR is produced in the near-infrared (LR 1064 = 36 sr), whereas in the visible considerably lower values were revealed (LR 532 = 19 sr). During PAMARCMiP2018, an increased BC concentration, reaching 117 ng m −3 , was revealed over Fram Strait by aircraft in situ measurements [81]. In general, black carbon has a longer lifetime due to its low deposition ability [82]. Thus, interstitial black carbon aerosol that survived the nucleation scavenging and dry deposition processes in combination with a further supply of pollution particles from north Europe could have resulted in the increased LR over Ny-Ålesund (64-82 sr).
The similarity of biomass-burning and industrial pollution particles in their optical and microphysical properties makes their distinction difficult [83]. Our study shows that under specific solar geometry and surface albedo conditions, these similarities are also reflected upon the ARE. The aerosol-related parameters, which are incorporated into radiative transfer simulations, are the SSA, the asymmetry factor, and the aerosol extinction coefficient. We will attempt to illuminate the role of each on the resulting ARE. From the SSA viewpoint, we cannot effectively tell any difference in the aerosol properties between Fram Strait and Ny-Ålesund. Our sensitivity analysis revealed that a 10% SSA perturbation can modulate the atmospheric ARE by 2.1 W m −2 , with the SSA modification being lower between the two locations. The asymmetry parameter did not reflect the pronounced coarse mode modification, since it has different sensitivity compared to other size relevant parameters, such as the Angström exponent [10]. Finally, despite the diminished extinction in the near-infrared (from 20 to 4-7 Mm −1 ), the impact on the spectrally integrated irradiance was not considerable, since this spectral region corresponds to low solar spectral irradiance. At the same time, the visible spectrum extinction, which is sensitive to the accumulation mode and the dominating visible solar irradiance, was almost constant. Therefore, the optically and microphysically similar accumulation mode led to an indistinguishable SW radiative footprint over the two locations.
Aerosol modification was reflected more intensely on the aerosol backscatter coefficient compared to the extinction. On the one hand, extinction, comprising scattering and absorption, is more representative for radiative transfer processes compared to solely scattering that mainly leads to the re-distribution of radiation. On the other hand, Lidar-derived extinction is characterized by high uncertainties due to the weak Raman cross-sections. Backscatter is accompanied by lower uncertainties and can also be derived from simple elastic Lidar systems [84], such as ceilometers, but it cannot be directly incorporated into aerosol climate models. Therefore, the aerosol retrieval uncertainties and limitations need to be thoroughly explored so as to better constrain the vertically resolved ARE from Lidar-inferred aerosol properties.
The derived aerosol properties, concentrations, and its change in time could, in principle, be compared to outcomes of chemical transport models. This has already been done successfully in the Arctic by Warneke et al. [12] and Zielinski et al. [85]. Despite the differences in the aerosol size distribution at the two sites, similar and stable results for the ARE have been found. Apparently, the forcing of the aerosol is more dependent on the solar zenith angle, surface albedo, and the geometrical height of the aerosol layer, as these parameters were similar for our case. Aerosol indirect effects could not be constrained from this case study.

Summary and Conclusions
In this study, we exploited synergistic remote sensing observations (from Lidar and sun-photometer) of a long-range aerosol transport episode over two parts of the European Arctic in April 2018. Our aims were oriented towards investigating the intra-Arctic aerosol modification and assessing its influence on the local radiative budget. The main findings and conclusions can be summarized as follows: • This transport episode stood out due to its considerable altitude, with such elevated aerosol layers predicted by aerosol climate models but rarely observed [14,60,61]. Generalizing the findings from our radiative transfer simulations, geometrically higher aerosol layers will produce a warming effect at higher tropospheric levels. Since these aerosol mixtures cool the surface, and hence near-surface air, but warm the surrounding air, the tropospheric stratification will be stabilized. Hence, the vertical propagation of heat and radiation fluxes will be suppressed. The latter applies for similar aerosol mixtures aloft bright surfaces during early spring.

•
The aerosol size distribution was clearly modified between the two Arctic locations. The effective radius of the accumulation mode decreased, while the coarse mode was depleted ( Figure 9). Industrial pollution and biomass-burning aerosol were funneled to Fram Strait and Ny-Ålesund, Svalbard. In the beginning, aged air masses originating from Asia dominated, while in the course of the episode they mixed with less mature air from northern Europe ( Figure 6). Solely, the modified aerosol sources and dry deposition could not account for the elimination of the coarse mode. However, high-level clouds in the interim of the two observations indicate the presence of nucleation scavenging ( Figure S1).

•
Along the intra-Arctic transport, particles presented a significantly higher Lidar ratio at visible wavelengths (LR 532 increased from 15 sr to 64-82 sr). A possible increased absorption can be attributed to black carbon, which could be transported at longer distances as interstitial aerosol since it has low deposition ability [82] and, thus, a longer lifetime. Moreover, black carbon could have been advected through neighboring north European intrusions.

•
Despite the coarse mode depletion and increased Lidar ratio, the microphysical similarities of the accumulation mode led to an indistinguishable short-wave radiative footprint in the two episode stages. The accumulation mode drives visible solar irradiance and, thus, dominates the short-wave radiative budget. Thus, the atmospheric column aerosol radiative effect bore similarities over the two locations, amounting to +4.4-4.9 W m −2 over the ice-covered Fram Strait and +4.5 W m −2 over the snow-covered Ny-Ålesund. This episode caused top-of-atmosphere warming accompanied by surface cooling, with implications for atmospheric stratification.

•
The aerosol radiative effect was more significant compared to the ambient humidity effect on a weekly basis (30% WV mixing ratio perturbation).

•
Our study suggests that within this Arctic haze episode, which was dominated by accumulation mode particles, the intra-Arctic aerosol modification (and hence, the precise aerosol microphysics) had little effect on the local radiative budget. However, in the context of retreating Arctic sea ice and declining surface albedo [86], the local aerosol radiative effect may change along individual transport episodes.