Journal cover Journal topic
Atmospheric Chemistry and Physics An interactive open-access journal of the European Geosciences Union
Journal topic
ACP | Articles | Volume 19, issue 9
Atmos. Chem. Phys., 19, 6217-6233, 2019
© Author(s) 2019. This work is distributed under
the Creative Commons Attribution 4.0 License.
Atmos. Chem. Phys., 19, 6217-6233, 2019
© Author(s) 2019. This work is distributed under
the Creative Commons Attribution 4.0 License.

Research article 10 May 2019

Research article | 10 May 2019

Vertical aerosol distribution in the southern hemispheric midlatitudes as observed with lidar in Punta Arenas, Chile (53.2° S and 70.9° W), during ALPACA

Vertical aerosol distribution above Punta Arenas
Andreas Foth1,2, Thomas Kanitz2,a, Ronny Engelmann2, Holger Baars2, Martin Radenz2, Patric Seifert2, Boris Barja3, Michael Fromm4, Heike Kalesse1, and Albert Ansmann2 Andreas Foth et al.
  • 1Leipzig Institute for Meteorology, University of Leipzig, Leipzig, Germany
  • 2Leibniz Institute for Tropospheric Research, Leipzig, Germany
  • 3Atmospheric Research Laboratory, University of Magallanes, Punta Arenas, Chile
  • 4U.S. Naval Research Laboratory, Washington, D.C., USA
  • anow at: European Space Agency, ESTEC, Noordwijk, the Netherlands
Back to toptop

Within this publication, lidar observations of the vertical aerosol distribution above Punta Arenas, Chile (53.2 S and 70.9 W), which have been performed with the Raman lidar PollyXT from December 2009 to April 2010, are presented. Pristine marine aerosol conditions related to the prevailing westerly circulation dominated the measurements. Lofted aerosol layers could only be observed eight times during the whole measurement period. Two case studies are presented showing long-range transport of smoke from biomass burning in Australia and regionally transported dust from the Patagonian Desert, respectively. The aerosol sources are identified by trajectory analyses with the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) and FLEXible PARTicle dispersion model (FLEXPART). However, seven of the eight analysed cases with lofted layers show an aerosol optical thickness of less than 0.05. From the lidar observations, a mean planetary boundary layer (PBL) top height of 1150 ± 350 m was determined. An analysis of particle backscatter coefficients confirms that the majority of the aerosol is attributed to the PBL, while the free troposphere is characterized by a very low background aerosol concentration. The ground-based lidar observations at 532 and 1064 nm are supplemented by the Aerosol Robotic Network (AERONET) Sun photometers and the space-borne Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) aboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO). The averaged aerosol optical thickness (AOT) determined by CALIOP was 0.02 ± 0.01 in Punta Arenas from 2009 to 2010.

1 Introduction
Back to toptop

Aerosols might possibly compensate the warming effect of the greenhouse gases in the Earth's radiation budget within the uncertainties of future climate modelling (Boucher et al.2013). The reason for the high uncertainties in the determination of the general aerosol radiative effect is the aerosols' variability in their global occurrence, their radiative properties (size, surface, chemistry) and their effects on cloud microphysics.

Global observations with spaceborne sensors improved the understanding of the seasonal distribution of aerosol layers worldwide, e.g. the seasonal vertical distribution of dust as observed with the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP; Winker et al.2007; Liu et al.2008). Nevertheless, extended field campaigns in key environments of the Earth with homogeneous aerosol conditions provided more detailed information about properties and cloud interaction of certain aerosol types with multi-sensor approaches, e.g. the Tropospheric Aerosol Radiative Forcing Observational Experiment (TARFOX; Russell et al.1999), the Second Aerosol Characterization Experiment (ACE2; Raes et al.2000), the Saharan Mineral Dust Experiment (SAMUM; Ansmann et al.2011) or the Saharan Aerosol Long-Range Transport and Aerosol–Cloud-Interaction Experiment (SALTRACE; Weinzierl et al.2017).

Within this publication, lidar observations of the vertical aerosol distribution above Punta Arenas, Chile (53.2 S and 70.9 W), are presented as performed during the Aerosol Lidar measurement in Punta Arenas in the frame of Chilean–GermAn cooperation (ALPACA) campaign which took place from December 2009 to April 2010. This location at the southern tip of South America yields an excellent opportunity to study almost clean marine aerosol conditions which are characteristic for the Southern Ocean (SO) because of the absence of continental land masses in the latitudinal belt south of 45 S and a constant westerly air flow from the Pacific Ocean (Schneider et al.2003). The nearest land mass situated towards the prevailing westerlies is New Zealand at a distance of roughly 8000 km and 10 latitude further north. Thus, reported aerosol sources at lower latitudes, like the Amazon rainforest (Reid et al.2004; Baars et al.2012), the Patagonian Desert (Gaiero et al.2003; Gassó and Stein2007; Li et al.2010; Johnson et al.2010) and the Australian continent (Edwards et al.2006), are expected to barely affect the aerosol conditions in Punta Arenas.

In the framework of the BACCHUS (Impact of Biogenic versus Anthropogenic emissions on Clouds and Climate: towards a Holistic UnderStanding) project, Carslaw et al. (2017) performed simulations to estimate the aerosol conditions in the year 1750 and their impact on climate. It was found that Punta Arenas is in a region that is still representative of pre-industrial aerosol conditions. A similar notation was already previously reported by Hamilton et al. (2014).

These pristine conditions already motivated the ground-based Aerosol Characterization Experiment (ACE1) in the 1990s (Bates et al.1998) and were confirmed by Minikin et al. (2003), who contrasted upper-tropospheric in situ aerosol aircraft observations in the northern midlatitudes (Scotland) and southern midlatitudes (Punta Arenas). However, only a few studies of ground-based aerosol and cloud layer profiling were performed in the southern midlatitudes in the following decades (Lopes et al.2013; Gouveia et al.2017), although such measurements would have served as good opportunity to contrast the aerosol radiative effect in the northern and southern midlatitudes, with respect to the aerosol sources, as well as the influence of aerosols on cloud microphysics (Kanitz et al.2011, 2013).

In the northern midlatitudes, lidar networks like the European Aerosol Lidar Network (EARLINET) (Bösenberg et al.2001) have monitored aerosol and cloud conditions for almost 20 years (Mattis et al.2008; Seifert et al.2010). A comparable network in Latin America, the Latin American Lidar Network (LALINET) was only fully established in 2013 (Guerrero-Rascado et al.2016; Antuña-Marrero et al.2017). Moreover, the South American Environment Risk Network (SAVER.Net) was established by means of a collaboration between Chile, Argentina and Japan to monitor aerosol, ozone and UV radiation in 2012 (Ristori et al.2018). The Atmospheric Research Laboratory of the University of Magallanes in Punta Arenas has participated in this activity with a multi-wavelength Raman and polarization lidar since 2016 (Barja et al.2018). As such, the conducted lidar measurements during ALPACA achieved the most comprehensive data set on aerosol and cloud distribution from ground before the establishment of these networks. The results of ALPACA are presented in this paper. Section 2 gives an overview about the experiment. The measurement systems used for this study are introduced in Sect. 3. Within ALPACA, lofted aerosol layers were observed only eight times. Two of these rare cases are representative of the others, including layers of Australian biomass-burning smoke and Patagonian Desert dust layers, are presented in Sect. 4, followed by an overview of the vertical aerosol distribution in Sect. 5. A conclusion about the results of ALPACA and an outlook for an upcoming campaign are given at the end.

2 Experiment
Back to toptop

Between the regular shipborne lidar measurements aboard the research vessel Polarstern (Kanitz et al.2013) in autumn 2009 and spring 2010, the portable lidar PollyXT (Althausen et al.2009; Engelmann et al.2016) was deployed at the University of Magallanes in Punta Arenas, Chile (53.2 S and 70.9 W), and continuous lidar measurements (24 h/7 d a week) were conducted from 4 December 2009 to 4 April 2010, covering a period of 4 months. Parts of these PollyXT observations have already been evaluated and published by Baars et al. (2016) in the frame of PollyNet.

Figure 1 shows the location of Punta Arenas (red star) in the very south of Chile, at the Strait of Magallanes and between the Pacific and Atlantic oceans. In this area, the polar front causes a continuous zonal wind band, because of the limited friction of the large ocean surface and the missing structured land masses in contrast to the northern midlatitudes (Cerveny1998). Hence, westerly winds prevail the whole year in southern Latin America. Based on the marine environment, daily and seasonal variations of the weather are weak (Coronato and Bisigato1998). The annual mean temperature is about 6 C and the annual precipitation amounts to 375 mm. Furthermore, the cyclone passage frequency is very high (3–5 d; Hodges et al.2011). This Southern Ocean region is characterized by a high cloud fraction (> 80 %; Naud et al.2014) with a cloud fraction of clouds below 3 km of about 60 % (Haynes et al.2011). The vegetation is composed of grasslands, tundra and mixed forest.

Figure 1Map of Antarctica and South America. Punta Arenas and Rio Gallegos are marked by a red and a yellow star, respectively.


3 Instruments
Back to toptop

3.1 PollyXT

In the framework of ALPACA, the lidar measurements were conducted with the portable multi-wavelength Raman and polarization lidar PollyXT_IFT (Althausen et al.2009), as part of the Polly lidar family (Engelmann et al.2016) and will be referred to PollyXT in the paper. Technically, the lidar is able to measure the backscattered light at 355, 532 and 1064 nm wavelengths, and Raman-scattered light at 387 and 607 nm to determine profiles of the particle backscatter coefficients at three wavelengths and the extinction coefficients at 355 and 532 nm. In fact, below 1500 m, the laser beam with the receiver field of view of the bistatic system is incomplete; thus, an overlap correction was applied. However, below 400 m, the overlap function is less than 0.5 and a reliable overlap correction is not possible. As a consequence, values of the particle backscatter coefficient were set constant below 400 m height in the lower part of the planetary boundary layer (PBL) under the assumption of well-mixed conditions.

The PBL top height is determined with the wavelet-covariance transformation that supposes a much higher aerosol load in the PBL than in the free troposphere (Brooks2003; Baars et al.2008).

The rather low aerosol content in the area of Punta Arenas caused low signal-to-noise ratios. Thus, the particle extinction coefficient had to be estimated from the 532 and 1064 nm backscatter coefficients by means of appropriate particle lidar ratio (SP) values.

3.2 Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP)

In April 2006, the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) mission started (Winker et al.2009). Aboard CALIPSO, the two-wavelength backscatter and polarization lidar, CALIOP, has been operated to achieve a worldwide four-dimensional data set of clouds and aerosols. CALIPSO orbits the Earth at a height of nearly 705 km with a velocity of 7 km s−1 and passes over the same location every 16th day.

The CALIPSO data processing provides profiles of backscatter and extinction coefficients at 532 and 1064 nm within the CALIOP level 2 version 4.10 data. In contrast to version 3 data, the version 4.10 data analysis algorithm aims at distinguishing not six but seven tropospheric aerosol subtypes (Omar et al.2009,, last access: 8 May 2019). In this study, the CALIPSO data are used to determine planetary boundary layer heights and the backscatter-related Ångström exponent, and to retrace the long-range transport of smoke based on total attenuated backscatter profiles at 532 nm (level 1 V4.10) and the aerosol subtype product (level 2 V4.10). Level 3 products of the monthly averaged aerosol optical thickness (AOT) are still based on version 3.10 because version 4.10 is not yet available.

The CALIOP data processing provides several quality flags. Within the scene classification algorithm, the cloud–aerosol discrimination (CAD) score is determined (Liu et al.2009). Aerosol particles are assigned negative values of −100 to −1 and clouds are assigned values of 1 to 100. The larger or smaller the value, the more confident the discrimination, respectively. Values around zero define an uncertain discrimination.

3.3 Auxiliary data

An Aerosol Robotic Network (AERONET) Sun photometer which measures AOT (column-integrated extinction coefficient) from 340 to 1020 nm at seven channels is located in Rio Gallegos (51.6 S, 69.3 W), Argentina (CEILAP-RG), which is 200 km away from Punta Arenas. Level 2.0 data are used with an AOT uncertainty of 0.01 to 0.02 (Holben et al.2001).

Two models were used to support the analysis of the air mass transport. The Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) is a model to calculate trajectories of air parcels for simulations of dispersion and deposition at arbitrary locations (Draxler and Hess1998; Stein et al.2015). By means of trajectories, aerosol sources are assigned. FLEXPART is a Lagrangian particle dispersion model, which calculates probabilistic trajectories of a large number of air parcels (Stohl et al.2005). Thereby, the transport and diffusion of aerosol could be described and a coarse assignment of aerosols to their sources is enabled.

In the statistical analysis of the aerosol conditions during ALPACA, we used ensemble backward trajectories combined with a land cover classification for a temporally and vertically resolved air mass source attribution. The “software for automated trajectory analysis: trace” is used (Radenz and Seifert2019). The land cover is a simplified version of the MODIS land cover (Friedl et al.2002). At first, a 27-member ensemble of 10 d backward trajectories is calculated using HYSPLIT. Meteorological input data for HYSPLIT are taken from the Global Data and Assimilation Service data set (GDAS1;, last access: 8 May 2019) provided by the Air Resources Laboratory (ARL) of the US National Weather Service's National Centers for Environmental Prediction (NCEP). Each ensemble is generated using a small spatial offset in the trajectory endpoint. Whenever a trajectory is below the PBL height provided in the GDAS1 data (“reception height”), the land cover is categorized using custom defined polygons according to land mass boundaries. Hence, an air parcel is assumed to be influenced by the surface if the trajectory is below the PBL height. The residence time for each category is then the total time an air parcel fulfilled this criterion by land cover category. This calculation is repeated in steps of 3 h in time and 500 m in height to provide a continuous estimate on the air mass source and as a first hint on potential aerosol load.

The global aerosol model NAAPS (Navy Aerosol Analysis and Prediction System;, last access: 8 May 2019) and the model results of the Monitoring Atmospheric Composition and Climate (MACC;, last access: 8 May 2019) project are used to obtain the modelled aerosol optical thickness of dust.

4 Case studies
Back to toptop

In this section, two cases are discussed in detail. The first one shows two lofted smoke layers after intercontinental long-range transport and low-level Patagonian dust. In the second one, further low-level Patagonian dust was observed.

4.1 2 March 2010: lofted smoke layers and Patagonian dust

On 2 March 2010, the weather in Punta Arenas was dominated by a northwesterly air flow. Between 09:00 and 12:00 UTC, PollyXT observed clouds in the range from 4 to 5 km and a PBL height of about 1 km (Fig. 2). An aerosol layer was situated above the PBL up to a height of 2 km. An increased backscatter was also observed at the height range from 4.5 to 5.2 km (06:00 to 12:00 UTC) and from 11.5 to 12.1 km (06:00 to 09:30 UTC). The source apportionment of the three observed aerosol layers by means of FLEXPART and HYSPLIT analyses is shown in Figs. 3, 4 and 5. In Fig. 3a, a FLEXPART analysis for the lower layer (1 to 2 km) reveals that the air masses passed the southern Pacific and southern Patagonia (red colouration in Fig. 3a). HYSPLIT 48 h backward trajectories for heights of 0.5, 1 and 2 km confirm that the ground-level aerosol on 2 March 2010 was locally affected by an extended residence time above Patagonia (Fig. 4). At this time, modelled AOT of dust at 550 nm from the global aerosol model NAAPS regionally peaked at 0.14 and was 0.02 in the area of Punta Arenas (Fig. 4). The AOT of dust at 550 nm from the MACC model reaches values of 0.07 in southern Patagonia.

Figure 2Height–time display of the range-corrected signal observed at 1064 nm with the PollyXT lidar on 2 March 2010. The analysed period between 06:15 and 09:00 UTC is framed in white. The planetary boundary layer top height is indicated by the black line.


Figure 3FLEXPART simulations for the integrated residence time of the particles that travelled in the whole atmospheric column within the last 10 d until the observation time on 2 March 2010 for the observed heights between 1 and 2 km (a), 4.5 and 5.2 km (b) and 11.5 and 12.1 km (c). The colours represent the logarithm of the integrated residence time (in seconds) in a grid box for 10 d integration time.


Figure 4HYSPLIT 48 h backward trajectories for heights of 0.5, 1 and 2 km for Punta Arenas on 2 March 2010, 09:00 UTC, and modelled dust surface concentrations from the NAAPS aerosol model on 2 March 2010, 06:00 UTC.


The FLEXPART analyses for the layers observed between 4.5 and 5.2 and 11.5 and 12.1 km height reveal a clearly defined region of origin. The observed layers were advected to Punta Arenas from Australia and the southern Pacific (Fig. 3b and c). HYSPLIT 13 d backward trajectories were calculated as well (blue and light blue in Fig. 5). Figure 5a also shows active fire spots (red dots) derived by MODIS (Moderate Resolution Imaging Spectroradiometer; Justice et al.2011), which were detected between 17 and 25 February 2010. On 18 and 19 February 2010, there were eruptions of pyrocumulonimbus (pyroCb) storms in western Australia. On 18 February, two fires, at approximately 32.6 S, 121.0 E and 33.0 S, 122.1 E, spawned pyroCb storms that were active between 04:00 and 10:00 UTC. The MODIS (Aqua) 11µm window brightness temperature minimum in the opaque core of the pyroCb anvil was -57C. According to the nearest radiosonde (Esperance, station number 94638; 33.8 S, 121.9 E; 00:00 UTC on 18 February) temperature profile, the brightness-temperature-inferred cloud-top height is ∼12.3 km. A third fire complex near 32.5 S, 122.8 E, produced pyrocumulus convection at that time. Hence, smoke was likely emitted at a range of altitudes during the various stages of pyroconvection from above the PBL up to ∼12 km. On 19 February, another pyroCb was detected that was likely generated from this third fire. The layer between 4.5 and 5.2 km (blue) as well as the layer between 11.5 and 12.1 km (light blue) crossed these regions with pyroconvection in western Australia. In summary, the source apportionment study for 2 March 2010 reveals that the near-surface aerosol layer is likely dominated by Patagonian dust, whereas the two observed lofted layers contain long-range-transported smoke from Australia. The PollyXT measurement was analysed for the cloud-free period between 06:15 and 09:00 UTC on 2 March 2010. Figure 6 illustrates the profiles of the optical properties with a smoothing length of 150 m. The particle backscatter coefficients reach their maximum values of up to 1.02 Mm−1 sr−1 (532 nm) and 0.63 Mm−1 sr−1 (1064 nm) in the PBL. The observed layers in the height–time display (Fig. 2) are also slightly visible in the profiles of the particle backscatter coefficient in heights of 5 and 12 km. The profile of the particle extinction coefficient (at 532 nm) was reproduced for an assumed constant lidar ratio. According to the origin of the air masses, a lidar ratio of 40 ± 10 sr (Kanitz et al.2013) was applied to the ground-level layer (up to 2.5 km) and a lidar ratio of 70 ± 10 sr (smoke; Ansmann et al.2009; Tesche et al.2011) was used for the lofted layers (Fig. 6c). The AOTs of each single layer result in values of 0.044 ± 0.004 (0 to 2.5 km), 0.004 ± 0.0004 (4.5 to 5.2 km) and 0.002 ± 0.0002 (11.5 to 12.1 km). The backscatter-related Ångström exponent (Fig. 6d) amounts to 0.56 ± 0.21 in the ground-level layer (up to 2 km). Kanitz et al. (2013) determined a Patagonian-dust-related Ångström exponent of 0.4 ± 0.1. In the framework of EARLINET, Saharan-dust-related Ångström exponents of 0.5 ± 0.5 were determined (Müller et al.2007). In the smoke layer between 4 and 5.5 km, the Ångström exponent was 0.61 ± 0.1. In comparison, biomass-burning-smoke-related Ångström exponents of 1.0 ± 0.4 were determined by lidar observations for long-range-transported smoke from Siberia and Canada (Müller et al.2007). In Cabo Verde, Ångström exponents of smoke originating from the south of western Africa were 1.06 ± 0.65 (Tesche et al.2011). Smoke that was transported from Africa to the Amazon rainforest was found to have Ångström exponents of 0.8 (Baars2012).

Figure 5(a) Map of the HYSPLIT 13 d backward trajectories for heights of 4.5 to 5.2 km (blue) and 11.5 to 12.1 km (light blue) arriving in Punta Arenas (brown star) on 2 March 2010, 09:00 UTC: MODIS fire counts between 17 to 25 February 2010 (red dots) and CALIPSO tracks (green). (b) Height of the according trajectories.


Figure 6Vertical profiles of (a) particle backscatter coefficients at 532 and 1064 nm, with error bars indicating 10 % uncertainty, (b) particle extinction coefficient at 532 nm, with error bars resulting from a lidar ratio uncertainty of ±10 sr, (c) particle lidar ratio, with error bars of ±10 sr and (d) backscatter-related Ångström exponent, with propagated error bars, derived by PollyXT on 2 March 2010, 06:15 to 09:00 UTC. The planetary boundary layer top height is indicated by the dashed black line.


In a next step, CALIPSO lidar observations were used for the characterization of the lofted aerosol layers during long-range transport. Six CALIPSO overpasses were found for the observation of the two lofted layers (see green lines in Fig. 5a), which provide intersections of the intercontinental transport of the smoke plume. On 21 February 2010, CALIPSO passed over the region of origin of the smoke in southern Australia. Figure 7a and b show the height–time display of the attenuated backscatter coefficient and the determined aerosol subtypes (with CAD score <-80) on 21 February 2010. Between heights of 4 and 7 km, a section of increased backscatter can be identified. According to the CALIOP data algorithms, the lofted aerosol layer is a mixture of smoke (black) and continental aerosol (green), which confirms the analysis of their origin discussed above. The AOT of the layer determined by CALIOP amounts to 0.1 (at 532 nm). AERONET Sun photometer measurements show a mean AOT of 0.165 (at 500 nm) in Canberra, southern Australia, on 21 February 2010. In contrast, a biomass-burning-related AOT of 0.55 (at 532 nm) was determined in the Amazon rainforest during the dry season (Baars2012). Additionally, an aerosol layer between 10 and 12 km is identified as clean continental aerosol.

Figure 7Height–time display of the total attenuated backscatter at 532 nm (left panels) and corresponding aerosol subtype (right panels) derived by CALIOP from six CALIPSO overpasses on (a, b) 21, (c, d) 22, (e, f) 24, (g, h) 27 February, (i, j) 2 and (k, l) 3 March 2010. The aerosol subtypes are only illustrated if the CAD score is below −80 excluding uncertain classifications.


On 22 February 2010, the smoke layer is also visible between heights of 4 and 8 km (Fig. 7c and d). On the following days, CALIOP was not able to determine unambiguously the smoke layer (Fig. 7e, f, g and h). A possible reason might be the decreasing smoke concentration along its transport route caused by dispersion and deposition (Bigg1973). The simultaneous occurrence of clouds and aerosol (for 24 February at 8 km height, see Fig. 7e, f; for 27 February at 5 to 11 km height, see Fig. 7g, h), which constrains the detection of the optically thin smoke plume, might be another reason. However, the CAD score of the detected smoke layers is smaller than −80, which indicates high discrimination accuracy (Liu et al.2009). In the region of Punta Arenas (on 2 March 2010 at a height of 4 km), CALIOP is again not able to doubtlessly determine the smoke layer, because all determined aerosol layers are in the vicinity of clouds. On 3 March 2010, the smoke layer was confidently detected by CALIOP (Fig. 7l).

4.2 17 February 2010: Patagonian dust

This subsection discusses the observation of a low-level aerosol layer that was observed with PollyXT on 17 February 2010 (see Fig. 8). A cyclone situated in the northwest of Punta Arenas caused northerly air flows on this and the previous days. The calculated 96 h HYSPLIT backward trajectories for 0.5, 1, 1.5 and 2 km reveal the South Pacific and southern Patagonia as the origin of the air masses (Fig. 9). The corresponding modelled NAAPS and MACC AOT (at 550 nm) exceeded 0.2 and 0.08 in southern Patagonia during this time, respectively (Fig. 9).

Figure 8Height–time display of the range-corrected signal observed at 1064 nm with the PollyXT lidar on 17 February 2010. The analysed period between 02:00 and 06:00 UTC is framed in white. The planetary boundary layer top height is indicated by the black line.


Figure 9HYSPLIT 96 h backward trajectories for heights of 0.5, 1, 1.5 and 2 km for Punta Arenas on 17 February 2010, 03:00 UTC, and modelled dust surface concentrations from the NAAPS aerosol model on 17 February 2010, 06:00 UTC.


In the lidar measurement between 00:00 and 06:00 UTC on 17 February 2010, an enhanced aerosol load can be deduced by an increased backscatter coefficient up to 2.5 km height. The surface layer below 1.2 km height seems to contain two sublayers, one below and one above 0.6 km height, respectively.

Figure 10Vertical profiles of (a) particle backscatter coefficient at 532 and 1064 nm, with error bars indicating 10 % uncertainty, (b) particle extinction coefficient at 532 nm, with error bars resulting from a lidar ratio uncertainty of ±10 sr, (c) particle lidar ratio with error bars of ±10 sr and (d) backscatter-related Ångström exponent, with propagated error bars, derived by PollyXT on 17 February 2010, 02:00 to 06:00 UTC. The planetary boundary layer top height is indicated by the dashed black line.


Vertical profiles of the particle backscatter coefficient of this measurement were generated for the time period from 02:00 to 06:00 UTC along with a vertical smoothing length of 150 m (framed at the bottom of Fig. 8). It is illustrated up to 5 km because the atmosphere was free of clouds and aerosol above this height. Figure 10a shows the particle backscatter coefficients at 532 nm (green) and 1064 nm (red), as derived with the Raman method. Both reach their maximum values of 2.5 Mm−1 sr−1 (532 nm) and 1.8 Mm−1 sr−1 (1064 nm) in the PBL. In the profiles of the particle backscatter coefficient, an aerosol layer up to 3 km height is visible. At heights above 3 km, the particle backscatter approaches zero. The profile of the particle extinction coefficient (532 nm) was reproduced for a Patagonian-dust-related lidar ratio of 40 ± 10 sr (Kanitz et al.2013) (Fig. 10b). The particle extinction coefficient reaches maximum values of 100 Mm−1. The vertical integration of the particle extinction coefficients below 3 km results in an AOT of 0.09 ± 0.01. That AOT value is well between values which would be derived for lidar ratios of other species of dust. For instance, assuming Saudi Arabian dust (lidar ratio of 38 sr; Müller et al.2007), Indian dust (43.8 sr; Schuster et al.2012) or Saharan dust (53 ± 7; Preißler et al.2011) would yield AOTs of 0.08, 0.1 and 0.12, respectively. Spaceborne measurements with MODIS between 16 and 18 February 2010 however indicate an increased AOT >0.2 over the southern Patagonian Desert. In the region of Punta Arenas, the AOT amounts to 0.07 up to 0.1. Whether the increased AOT close to Punta Arenas is caused by dust load could not be clearly identified. Comparable AOTs were also determined by MODIS in the southwest of Patagonia.

The mean Ångström exponent was 0.48 ± 0.06 (Fig. 10d), which is within the range of the values for Patagonian dust described in Sect. 4.1 and which is also confirmed by Müller et al. (2007) and Kanitz et al. (2013).

5 Statistical results
Back to toptop

5.1 General aerosol conditions

The general aerosol conditions in Punta Arenas are presented in Fig. 11a in terms of monthly averaged AOT at 532 nm wavelength as determined with CALIOP for the grid cell of Punta Arenas (blue curve) from 1 January 2009 to 31 December 2010. The averaged AOT in Punta Arenas from 2009 to 2010 was 0.02 ± 0.01. The annual course of the monthly averaged AOT indicates the absence of a pronounced seasonal cycle. The large standard deviation may have been caused by the low number of observations but may also be the result of the CALIOP aerosol typing limitations in coastal regions (Kanitz et al.2014). For comparison, the AOT as obtained at the AERONET station of Rio Gallegos is shown, too. In Rio Gallegos, the mean AOT between 2009 and 2010 (0.02 ± 0.02) is in the same range as that in Punta Arenas, although Rio Gallegos is situated closer to the Patagonian Desert and 400 km east of the west coast of Latin America. AERONET AOT measurements in Rio Gallegos confirm the very low mean AOT values (2009: 0.02 ± 0.01, 2010: 0.02 ± 0.01) of CALIOP and indicate clean marine conditions in Punta Arenas and Rio Gallegos (Fig. 11b). Such low AOT values were also found in other coastal and remote oceanic areas. During three meridional transatlantic cruises from 50 N to 50 S, shipborne lidar measurements revealed AOTs (at 532 nm) of the marine boundary layer of below 0.05 in 78 % of the cases (Kanitz et al.2013). Wilson and Forgan (2002) determined a mean AOT (at 500 nm) of less than 0.04 in Cape Grim, Tasmania, from 1986 to 1999. AERONET Sun photometer measurements in marine areas show AOTs (at 500 nm) of 0.085 ± 0.01 over the Pacific and 0.06 ± 0.02 over the Southern Ocean (Smirnov et al.2009).

Figure 11(a) Monthly average of CALIOP level 3 AOT (at 532 nm) and their standard deviation in Punta Arenas (blue) and Rio Gallegos (red) in 2009 and 2010. The ALPACA campaign is indicated by the shaded area. (b) Daily average of AERONET Sun photometer AOT (at 500 nm) in Rio Gallegos for the ALPACA campaign. (c) Height–time display of the range-corrected signal (at 1064 nm) measured by PollyXT. Panel (d) shows the accumulated residence time of backward trajectories separated by different regions of origin.


Figure 11c shows the height–time display of the range-corrected signal at 1064 nm wavelength measured by PollyXT for the entire measurement period. As expected, most of the aerosol load is contained within the PBL up to around 1200 m (reddish colours). The free troposphere is characterized by a very low aerosol load but a frequent occurrence of clouds (grey and white colours).

A comprehensive analysis of aerosol source regions based on ensemble HYSPLIT trajectory calculations shows (Fig. 11d) that the influence of the ocean on air parcels (given in accumulated residence time of the backward trajectories within the PBL over the respective surface type) reaching Punta Arenas is a factor of 100 larger in contrast to the continents Africa, Australia and even South America.

5.2 Vertical aerosol distribution

After introducing the general aerosol conditions in Punta Arenas, in this section, the vertically resolved measurements with PollyXT are investigated to determine the vertical aerosol distribution. First of all, the ALPACA lidar measurements are applied to ascertain the height of the PBL. Figure 12a presents the monthly means of the PBL heights. During the period with the strongest warming in the southern hemispheric summer (December and January), the maximum averaged top heights were reached (1230 ± 331 m and 1177 ± 365 m, respectively). The PBL heights decrease to 1106 ± 317 and 984 ± 347 m, respectively, in February and March, due to decreasing solar irradiation. However, the trend to low values is within the range of the standard deviation, which is indicated as error bars. Conducting stationary and continuous lidar measurements with PollyXT, the determination of the temporal development of the PBL heights is possible (Cohn and Angevine2000; Baars et al.2008). Figure 12b shows the averaged diurnal variation in the PBL height with a time resolution of 3 h. However, no explicit diurnal variation is identifiable within the limits of the error bars. The values vary between 850 ± 290 and 1280 ± 370 m. In contrast, for a northern hemispheric site with comparable latitude, PBL heights larger than 2 km were observed in 59 % of all considered cases in Leipzig (51 N, 12 E; Mattis et al.2008). In Granada, Spain, which has a similar distance to the coast, measurements showed an annual mean of the PBL height around 1.7 ± 0.5 km (Granados-Muñoz et al.2012).

Figure 12Height of the PBL determined by PollyXT between 4 December 2009 and 31 March 2010. (a) Monthly averaged PBL top heights including the standard deviation as error bars and (b) averaged PBL top height as function of time of day (averaged over 3 h) including standard deviation as error bars. The red numbers indicate the number of measurements of the according bar. Local time is UTC4:00.


Figure 13a illustrates the frequency distribution of the obtained PBL heights as obtained from PollyXT, CALIOP, radiosonde and GDAS1. For the ground-based lidar observations, in 74 % of all cases, the PBL extends up to heights between 750 and 1500 m, and the mean PBL top height amounts to 1151 ± 347 m. CALIPSO passed over Punta Arenas 31 times in the period from 1 May 2009 to 30 April 2010. The CALIPSO aerosol subtype data set is used for the classification of the PBL height. From these, the upper limit of the lowest aerosol layer, if present at all, was assumed to be the PBL height. In total, PBL heights could be estimated for a total of 21 overpasses. The mean PBL height identified with CALIOP is 1169 ± 532 m (Fig. 13b). Overall, 52 % of the PBL heights determined by CALIOP are between 750 and 1500 m. Furthermore, PBL heights were derived from radiosondes (each at 12:00 UTC) at the airport of Punta Arenas (37 m a.s.l., 15 km away) between 4 December 2009 and 31 March 2010 (Fig. 13). The mean PBL height is 1019 ± 376 m, and 64 % of all values are between 750 and 1500 m. In total, 73 % of all radiosonde ascents show a further lower layer with a mean height of 190 ± 13 m, which is indicated by a strong temperature gradient. The PBL thus seems to consist of two thermodynamically distinct layers. The reason for that may be explained by the terrain in the area of Punta Arenas which is directly located at the Strait of Magallanes. The southern Andes mountains and the Pacific are situated in the west (Fig. 1). The lower aerosol layer might be considered as the inflow of the marine boundary layer into the continental boundary layer. However, this layer is rarely detected by PollyXT due to the effect of incomplete overlap below 400 m height. For comparison, only the upper layer is considered. The PBL heights provided by the GDAS1 data set (see Sect. 3.3) are calculated on the basis of the gradient Richardson number (Stull1988) and their mean is about 1024 ± 463 m (Fig. 13). Although the radiosonde data are assimilated into GDAS1 data set, the peak in low PBL heights derived from the radiosondes does not appear in the GDAS1 data. The horizontal resolution of GDAS1 data is 1×1 and the vertical one is approximately 250 m in the lowermost 1000 m. Locations between these grid points are interpolated, which might explain the observed differences in the PBL height distributions. In 6 % of all cases, the GDAS1 PBL heights solely fall below 250 m. Generally, the PBL heights determined by the decrease in aerosol backscatter using PollyXT and CALIOP agree with those determined by the potential temperature profile using GDAS1 and radiosondes. That means that the top height of the aerosol layer coincides with the temperature inflection; hence, the aerosol accumulates within the PBL in the region of Punta Arenas.

Figure 13Frequency distribution of the PBL heights determined by (a) PollyXT between 4 December 2009 and 31 March 2010, (b) CALIOP between 1 May 2009 and 30 April 2010, (c) radiosonde (12:00 UTC) and (d) GDAS1 data (12:00 UTC) from 4 December 2009 to 31 March 2010 with an increment of 250 m, respectively. N gives the number of samples. The horizontal error bar indicates the standard deviation.


In the following, statistics of the vertical aerosol distribution are presented. During the ALPACA campaign, 59 measurement periods were analysed. The low amount of determined aerosol profiles is caused by the frequent occurrence of low and mid-level clouds (about 83 % of the measurement period; Kanitz et al.2011) hampering an evaluation of the aerosol lidar data. Furthermore, the analysis is limited due to the frequent presence of marginal aerosol concentrations in the atmosphere and a corresponding low signal-to-noise ratio measured in the lidar signals. Figure 14 shows the analysed clear-sky profiles of the particle backscatter coefficient at 532 nm (Fig. 14a) and 1064 nm (Fig. 14b). The general vertical smoothing length was 330 m. Well-mixed homogeneous aerosol conditions are assumed to be present in the PBL, such that the particle backscatter coefficient is set constant in the overlap region below 400 m height. The profiles of the particle backscatter coefficient which were derived by either the Raman or Klett method (Fig. 14) indicate that the particle backscatter coefficient is close to 0 Mm−1 sr−1 above a height of 2 km. Hence, the majority of the aerosol is located in the PBL, as also indicated in the good correlation between radiosonde- and PollyXT-derived PBL heights, discussed above. The particle backscatter coefficients reach mean values of 0.74 ± 0.56 Mm−1 sr−1 (532 nm) and 0.43 ± 0.32 Mm−1 sr−1 (1064 nm). However, aerosol layers were almost never observed in the free troposphere. Therefore, the free troposphere may be considered as a region representing pristine background aerosol conditions, making it an ideal region of low-aerosol reference for studies of aerosol–cloud interactions (Kanitz et al.2011).

Figure 14Single (coloured lines) and averaged (black lines) height profiles of the particle backscatter coefficients at 532 nm (a) and at 1064 nm (b) as derived from PollyXT observations. The error bars indicate the standard deviations.


Lidar measurements as well as AERONET measurements provide spectrally resolved information. The spectral behaviour is expressed by means of the Ångström exponent. Figure 15 displays the frequency distribution of the vertical averaged backscatter-related Ångström exponent (at 532 and 1064 nm) for the time periods from 4 December 2009 to 4 April 2010 using PollyXT (Fig. 15a) and 1 May 2009 and 30 April 2010 using CALIOP (Fig. 15b). For comparison, Fig. 15c illustrates the AOT-related Ångström exponent of the AERONET station in Rio Gallegos for the ALPACA period. In all three cases, the maximum values of the distribution are between 0 and 1. The mean values are 0.8 ± 0.3, 0.3 ± 0.7 and 0.5 ± 0.5 for the measurements with PollyXT, CALIOP and AERONET, respectively. By means of these values, the size of the particles can be derived qualitatively. Very low or negative Ångström exponents indicate very large aerosol particles (sea salt or dust) (Moulin et al.1997; Müller et al.2007). In turn, large Ångström exponents (>2) point towards small aerosol particles (such as fresh smoke; Baars et al.2012). The fraction of low Ångstöm exponents is largest for Rio Gallegos. This is caused by the occurrence of more dust events in the outflow of Patagonia. The values derived by CALIOP are too noisy to allow a reasonable interpretation. The Ångstöm exponents determined by PollyXT are representative of continental aerosol as the lowermost heights of the atmosphere below 400 m height containing most likely a marine contribution could not be analysed due to instrumental limitations as discussed above.

Figure 15Frequency distribution of the vertical integrated particle-backscatter-related Ångström exponent (532 and 1064 nm) of (a) PollyXT and (b) CALIOP with an increment of 0.2. Additionally, the frequency distribution of the AOT-related Ångström exponent determined by Sun photometer measurements in Rio Gallegos is illustrated. N gives the number of measurements. The horizontal error bar indicates the standard deviation.


6 Conclusions and outlook
Back to toptop

The presented study aimed on providing an overview about the vertical aerosol conditions above Punta Arenas. During the 4 months of observations of PollyXT, lofted aerosol layers were rarely observed and, when present, were characterized by very low optical thicknesses. Overall, the mean aerosol optical thickness in Punta Arenas was found to be 0.02, which is very low, even for marine conditions. Due to the absence of observations of anthropogenic aerosols, neither local nor long-range-transported, our study thus confirms well the conclusions of Hamilton et al. (2014) and Carslaw et al. (2017) that the atmosphere over southern Chile still provides pristine, pre-industrial conditions. CALIPSO observations, which were utilized to track the long-range transport of aerosol from Australia to Punta Arenas, indicate that a considerable fraction of free-tropospheric aerosol is removed by cloud processes and washout taking place over the Pacific Ocean before it reaches South America. The average free-tropospheric aerosol load is thus subject to increase to the west of Punta Arenas. The same can be expected eastward of South America, because, as was found in analyses of NAAPS and MACC aerosol model simulations, substantial amounts of Patagonian dust are frequently emitted into the atmosphere and transported leeward of South America. We thus conclude from our study that Punta Arenas is one of the few accessible places on Earth with temperate climate where aerosol–cloud interaction reference studies in the absence of free-tropospheric aerosols can be conducted. Nevertheless, the significance of the ALPACA data set is somewhat limited because the lidar observations suffered from a rather large height of full overlap and the absence of polarization and UV measurements. These conditions hampered us from more detailed analyses of the aerosol type and the marine contribution to the total aerosol load. As a consequence, the PollyXT_IFT system was reconstructed in the meantime to improve the performance under very different aerosol conditions. It is now equipped with a better data acquisition in combination with new photomultiplier tubes, a near-range receiver and an additional depolarization channel as well as a water vapour detection channel and is named PollyXT_TROPOS (Engelmann et al.2016).

Future studies might enhance the knowledge of the aerosol conditions regarding an aerosol typing (Baars et al.2017) and separation of aerosol types as well as an estimation of ice-nucleating particles and cloud condensation nuclei from multi-wavelength Raman and polarization lidar observations (Mamouri and Ansmann2015).

In recent years, the Southern Ocean, especially the area in the south and southeast of Australia, became an increasing focus for atmospheric researchers. Several observation campaigns like MARCUS (Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean), SOCRATES (the Southern Ocean Clouds Radiation Aerosol Transport Experimental Study) and CAPRICORN (Clouds, Aerosols, Precipitation Radiation and atmospherIc Composition Over the southeRN ocean; Protat et al.2016; Mace and Protat2018a, b) provide beneficial data for improving the understanding of aerosol–cloud interaction in pristine environments.

In the upcoming field experiment organized by the Leibniz Institute for Tropospheric Research (TROPOS), Leipzig, Germany, in collaboration with the University of Magallanes, Punta Arenas (UMAG), and the Institute for Meteorology at the University of Leipzig (LIM), an extended version of the Leipzig Aerosol and Cloud Remote Observation System (LACROS; Bühl et al.2013) will be deployed at UMAG in November 2018 for at least 1 year to study the seasonal cycle of aerosols and clouds in Punta Arenas. LACROS comprises, amongst others, a PollyXT lidar, a 35 GHz cloud radar, a Doppler lidar, a microwave radiometer, a disdrometer and radiation sensors for direct and diffuse downwelling and upwelling solar and thermal radiation. During this field experiment, named DACAPO-PESO (Dynamics, Aerosol, Cloud and Precipitation Observations in the Pristine Environment of the Southern Ocean), the LACROS suite is extended by a 94 GHz frequency-modulated continuous wave cloud radar provided by LIM (LIMRAD94; Küchler et al.2017) and a 24 GHz micro-rain radar (TROPOS) which will allow for multi-frequency polarimetric Doppler radar studies. Complementing instrumentation of UMAG includes radiosondes, radiation observations, in situ aerosol observations and multi-wavelength lidar measurements. The latter have been performed since 2015 in the frame of LALINET (Guerrero-Rascado et al.2016; Antuña-Marrero et al.2017).

With this extended observational suite, several research questions, like the efficiency of liquid-dependent ice formation as well as the influence of aerosol concentration on the frequency of occurrence of mixed-phase cloud processes like aggregation and riming, will be addressed. These observations might help to improve the representation of Southern Ocean clouds in global climate models which currently suffer from a strong radiation bias (Bodas-Salcedo et al.2014) caused by a misrepresentation of cloud phase. Specifically, the amount of ice is overestimated by models, while observations show a large amount of clouds with high supercooled liquid water content at cloud top (Haynes et al.2011). Since cloud thermodynamics are a strong function of cloud condensation nuclei and ice-nucleating particle availability, which in turn are related to the aerosol load and type (Mamouri and Ansmann2015), the results from the ALPACA campaign, presented here, build a fundamental basis of knowledge about the aerosol conditions, cloud condensation nuclei and ice-nucleating particles that can be expected at this site.

Data availability
Back to toptop
Data availability. 

The Polly lidar data are available at TROPOS upon request ( CALIPSO data were downloaded from the NASA Atmospheric Science Data Center web page (, last access: 8 May 2019). The CALIPSO CALIOP 5 km aerosol profile product version 4.1 is available under (Vaughan et al.2016). Backward trajectories' analysis was supported by air mass transport computation with the NOAA (National Oceanic and Atmospheric Administration) HYSPLIT model (HYSPLIT, 2018) using GDAS meteorological data (Stein et al.2015; Rolph et al.2017). AERONET Sun photometer AOT data were downloaded from the AERONET web page (, last access: 8 May 2019). Data from the global aerosol model NAAPS (Navy Aerosol Analysis and Prediction System) and MACC were downloaded from the NAAPS ( and MACC ( web page. The software for the automated trajectory analysis: trace (version v0.3) is available under (Radenz and Seifert2019).

Author contributions
Back to toptop
Author contributions. 

AF prepared the manuscript in close cooperation with TK, HB, PS, BB and HK. AF and TK performed the investigations and data analyses. TK, HB and RE realized the experimental setup and were responsible for the high quality of the lidar measurements. RE realized the optical and technical setup of the lidar. HB provided the software for the analysis of the lidar data and supported the data analysis and interpretation. MR developed and applied the code for the determination of the accumulated residence time of backward trajectories separated by different regions of origin, and he prepared Fig. 1. MF contributed to the section on the pyroCb. The conceptualization was initialized by AA. All authors have contributed to the scientific discussions.

Competing interests
Back to toptop
Competing interests. 

The authors declare that they have no conflict of interest.

Back to toptop

First of all, we thank Dietrich Althausen (TROPOS), Felix Zamorano (UMAG) and Claudio Casiccia (UMAG). Without your support, the lidar could not have been operated in Punta Arenas.

The authors acknowledge support through the High-Definition Clouds and Precipitation for advancing Climate Prediction research programme (HD(CP)2; FKZ: 01LK1504C and 01LK1502N) funded by Federal Ministry of Education and Research in Germany (BMBF), ACTRIS under grant agreement no. 262254 of the European Union's Seventh Framework Programme (FP7/2007-2013), ACTRIS-2 under grant agreement no. 654109 from the European Union's Horizon 2020 research and innovation programme, EUCAARI funded by the European Union (FP6, grant no. 036 833-2) and the Gottfried Wilhelm Leibniz Association (OCEANET project in the framework of PAKT).

Contributions by Heike Kalesse were made with support of the project PICNICC, GZ: KA 4162/2-1, within the priority programme PROM of the German Science Foundation (DFG).

We also thank the NASA Langley Research Center and the CALIPSO science team for the constant effort and improvement of the CALIPSO data. Supplementary information from HYSPLIT trajectories, NAAPS and MACC aerosol modelling and MODIS was a cornerstone of our data analysis.

We especially acknowledge the work of Brent Holben, Eduardo Quel, Lidia Otero and Jacobo Salvador for operating the AERONET station in Rio Gallegos.

Review statement
Back to toptop
Review statement. 

This paper was edited by Anne Perring and reviewed by two anonymous referees.

Back to toptop

Althausen, D., Engelmann, R., Baars, H., Heese, B., Ansmann, A., Müller, D., and Komppula, M.: Portable Raman lidar PollyXT for automated profiling of aerosol backscatter, extinction, and depolarization, J. Atmos. Ocean. Tech., 26, 2366–2378,, 2009. a, b

Ansmann, A., Baars, H., Tesche, M., Müller, D., Althausen, D., Engelmann, R., Pauliquevis, T., and Artaxo, P.: Dust and smoke transport from Africa to South America: Lidar profiling over Cape Verde and the Amazon rainforest, Geophys. Res. Lett., 36, L11802,, 2009. a

Ansmann, A., Petzold, A., Kandler, K., Tegen, I., Wendisch, M., Müller, D., Weinzierl, B., Müller, T., and Heintzenberg, J.: Saharan Mineral Dust Experiments SAMUM-1 and SAMUM-2: what have we learned?, Tellus, 63, 403–429,, 2011. a

Antuña-Marrero, J. C., Landulfo, E., Estevan, R., Barja, B., Robock, A., Wolfram, E., Ristori, P., Clemesha, B., Zaratti, F., Forno, R., Armandillo, E., Bastidas, Á. E., de Frutos Baraja, Á. M., Whiteman, D. N., Quel, E., Barbosa, H. M. J., Lopes, F., Montilla-Rosero, E., and Guerrero-Rascado, J. L.: LALINET: The First Latin American–Born Regional Atmospheric Observational Network, B. Am. Meteorol. Soc., 98, 1255–1275,, 2017. a, b

Baars, H.: Aerosol profiling with lidar in the Amazon Basin during the wet and dry season 2008, PhD thesis, Faculty of Physics and Geoscience, University of Leipzig, Leipzig, Germany, 2012. a, b

Baars, H., Ansmann, A., Engelmann, R., and Althausen, D.: Continuous monitoring of the boundary-layer top with lidar, Atmos. Chem. Phys., 8, 7281–7296,, 2008. a, b

Baars, H., Ansmann, A., Althausen, D., Engelmann, R., Heese, B., Müller, D., Artaxo, P., Paixao, M., Pauliquevis, T., and Souza, R.: Aerosol profiling with lidar in the Amazon Basin during the wet and dry season, J. Geophys. Res.-Atmos., 117, D21201,, 2012. a, b

Baars, H., Kanitz, T., Engelmann, R., Althausen, D., Heese, B., Komppula, M., Preißler, J., Tesche, M., Ansmann, A., Wandinger, U., Lim, J.-H., Ahn, J. Y., Stachlewska, I. S., Amiridis, V., Marinou, E., Seifert, P., Hofer, J., Skupin, A., Schneider, F., Bohlmann, S., Foth, A., Bley, S., Pfüller, A., Giannakaki, E., Lihavainen, H., Viisanen, Y., Hooda, R. K., Pereira, S. N., Bortoli, D., Wagner, F., Mattis, I., Janicka, L., Markowicz, K. M., Achtert, P., Artaxo, P., Pauliquevis, T., Souza, R. A. F., Sharma, V. P., van Zyl, P. G., Beukes, J. P., Sun, J., Rohwer, E. G., Deng, R., Mamouri, R.-E., and Zamorano, F.: An overview of the first decade of PollyNET: an emerging network of automated Raman-polarization lidars for continuous aerosol profiling, Atmos. Chem. Phys., 16, 5111–5137,, 2016. a

Baars, H., Seifert, P., Engelmann, R., and Wandinger, U.: Target categorization of aerosol and clouds by continuous multiwavelength-polarization lidar measurements, Atmos. Meas. Tech., 10, 3175–3201,, 2017. a

Barja, B., Zamorano, F., Ristori, P., Otero, L., Quel, E., Sugimoto, N., Shimizu, A., and Santana, J.: Aerosols Observations with a new lidar station in Punta Arenas, Chile, EPJ Web Conf., 176, 05038,, 2018. a

Bates, T. S., Huebert, B. J., Gras, J. L., Griffiths, F. B., and Durkee, P. A.: International Global Atmospheric Chemistry (IGAC) Project's First Aerosol Characterization Experiment (ACE 1): Overview, J. Geophys. Res., 1031, 16297–16318,, 1998. a

Bigg, E.: Ice Nucleus Concentrations in Remote Areas., J. Atmos. Sci., 30, 1153–1157,<1153:INCIRA>2.0.CO;2, 1973. a

Bodas-Salcedo, A., Williams, K. D., Ringer, M. A., Beau, I., Cole, J. N. S., Dufresne, J.-L., Koshiro, T., Stevens, B., Wang, Z., and Yokohata, T.: Origins of the Solar Radiation Biases over the Southern Ocean in CFMIP2 Models, J. Climate, 27, 41–56,, 2014. a

Bösenberg, J., Ansmann, A., Baldasano, J., Balis, D., Böckmann, C., Calpini, B., Chaikovsky, A., Flamant, P., Hågård, A., Mitev, V., Papayannis, A., Pelon, J., Resendes, D., Schneider, J., Spinelli, N., Trickl, T., Vaughan, G., Visconti, G., and Wiegner, M.: EARLINET: a European aerosol research lidar network, Advances in laser remote sensing, edited by: Dabas, A., Loth, C., and Pelon, J., Editions de L'Ecole Polytechnique, Palaiseau, France, 155–158, 2001. a

Boucher, O., Randall, D., Artaxo, P., Bretherton, C., Feingold, G., Forster, P., Kerminen, V. M., Kondo, Y., Liao, H., Lohmann, U., Rasch, P., Satheesh, S. K., Sherwood, S., Stevens, B., and Zhang, X. Y.: Clouds and Aerosols, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T., Qin, D., Plattner, G., Tignor, M., Allen, S., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P., Cambridge University Press, Cambridge, UK and New York, NY, USA, 2013. a

Brooks, I. M.: Finding Boundary Layer Top: Application of a Wavelet Covariance Transform to Lidar Backscatter Profiles, J. Atmos. Ocean. Tech., 20, 1092,<1092:FBLTAO>2.0.CO;2, 2003. a

Bühl, J., Seifert, P., Wandinger, U., Baars, H., Kanitz, T., Schmidt, J., Myagkov, A., Engelmann, R., Skupin, A., Heese, B., Klepel, A., Althausen, D., and Ansmann, A.: LACROS: the Leipzig Aerosol and Cloud Remote Observations System, in: Remote Sensing of Clouds and the Atmosphere XVIII; and Optics in Atmospheric Propagation and Adaptive Systems XVI, vol. 8890 of Proc. SPIE, 889002,, 2013. a

Carslaw, K. S., Gordon, H., Hamilton, D. S., Johnson, J. S., Regayre, L. A., Yoshioka, M., and Pringle, K. J.: Aerosols in the Pre-industrial Atmosphere, Current Climate Change Reports, 3, 1–15,, 2017. a, b

Cerveny, R. S.: Present Climates of South America, in: Climates of the Southern Continents: Present, Past and Future, edited by: Hobbs, J. E., Lindesay, J. A., and Bridgeman, H. A., 107–135, John Wiley, Chichester, UK, 107–134, 1998. a

Cohn, S. A. and Angevine, W. M.: Boundary Layer Height and Entrainment Zone Thickness Measured by Lidars and Wind-Profiling Radars., J. Appl. Meteorol., 39, 1233–1247,<1233:BLHAEZ>2.0.CO;2, 2000. a

Coronato, F. and Bisigato, A.: A temperature pattern classification in Patagonia, Int. J. Climatol., 18, 765–773,<765::AID-JOC282>3.0.CO;2-H, 1998. a

Draxler, R. R. and Hess, G. D.: An overview of the HYSPLIT4 modelling system for trajectories, dispersion, and deposition, Aust. Meteorol. Mag., 47, 295–308, 1998. a

Edwards, D. P., Emmons, L. K., Gille, J. C., Chu, A., Attié, J.-L., Giglio, L., Wood, S. W., Haywood, J., Deeter, M. N., Massie, S. T., Ziskin, D. C., and Drummond, J. R.: Satellite-observed pollution from Southern Hemisphere biomass burning, J. Geophys. Res., 111, D14312,, 2006. a

Engelmann, R., Kanitz, T., Baars, H., Heese, B., Althausen, D., Skupin, A., Wandinger, U., Komppula, M., Stachlewska, I. S., Amiridis, V., Marinou, E., Mattis, I., Linné, H., and Ansmann, A.: The automated multiwavelength Raman polarization and water-vapor lidar PollyXT: the neXT generation, Atmos. Meas. Tech., 9, 1767–1784,, 2016. a, b, c

Friedl, M., McIver, D., Hodges, J., Zhang, X., Muchoney, D., Strahler, A., Woodcock, C., Gopal, S., Schneider, A., Cooper, A., Baccini, A., Gao, F., and Schaaf, C.: Global land cover mapping from MODIS: algorithms and early results, Remote Sens. Environ., 83, 287–302,, 2002. a

Gaiero, D., Probst, J., Depetris, P., Bidart, S., and Leleyter, L.: Iron and other transition metals in Patagonian riverborne and windborne materials: Geochemical control and transport to the southern South Atlantic Ocean, Geochim. Cosmochim. Ac., 67, 3603–3623,, 2003. a

Gassó, S. and Stein, A.: Does dust from Patagonia reach the sub-Antarctic Atlantic Ocean?, Geophys. Res. Lett., 34, L01801,, 2007. a

Gouveia, D. A., Barja, B., Barbosa, H. M. J., Seifert, P., Baars, H., Pauliquevis, T., and Artaxo, P.: Optical and geometrical properties of cirrus clouds in Amazonia derived from 1 year of ground-based lidar measurements, Atmos. Chem. Phys., 17, 3619–3636,, 2017. a

Granados-Muñoz, M. J., Navas-Guzmán, F., Bravo-Aranda, J. A., Guerrero-Rascado, J. L., Lyamani, H., Fernández-Gálvez, J., and Alados-Arboledas, L.: Automatic determination of the planetary boundary layer height using lidar: One-year analysis over southeastern Spain, J. Geophys. Res.-Atmos., 117, D18208,, 2012. a

Guerrero-Rascado, J. L., Landulfo, E., Antuña, J. C., Barbosa, H. d. M. J., Barja, B., Bastidas, Á. E., Bedoya, A. E., da Costa, R. F., Estevan, R., Forno, R., Gouveia, D. A., Jiménez, C., Larroza, E. G., Lopes, F. J. d. S., Montilla-Rosero, E., Moreira, G. d. A., Nakaema, W. M., Nisperuza, D., Alegria, D., Múnera, M., Otero, L., Papandrea, S., Pallota, J. V., Pawelko, E., Quel, E. J., Ristori, P., Rodrigues, P. F., Salvador, J., Sánchez, M. F., and Silva, A.: Latin American Lidar Network (LALINET) for aerosol research: Diagnosis on network instrumentation, J. Atmos. Sol.-Terr. Phy., 138–139, 112–120,, 2016. a, b

Hamilton, D. S., Lee, L. A., Pringle, K. J., Reddington, C. L., Spracklen, D. V., and Carslaw, K. S.: Occurrence of pristine aerosol environments on a polluted planet, P. Natl. Acad. Sci. USA, 111, 18466–18471,, 2014. a, b

Haynes, J. M., Jakob, C., Rossow, W. B., Tselioudis, G., and Brown, J.: Major Characteristics of Southern Ocean Cloud Regimes and Their Effects on the Energy Budget, J. Climate, 24, 5061–5080,, 2011. a, b

Hodges, K. I., Lee, R. W., and Bengtsson, L.: A Comparison of Extratropical Cyclones in Recent Reanalyses ERA-Interim, NASA MERRA, NCEP CFSR, and JRA-25, J. Climate, 24, 4888–4906,, 2011. a

Holben, B. N., Tanré, D., Smirnov, A., Eck, T. F., Slutsker, I., Abuhassan, N., Newcomb, W. W., Schafer, J. S., Chatenet, B., Lavenu, F., Kaufman, Y. J., Castle, J. V., Setzer, A., Markham, B., Clark, D., Frouin, R., Halthore, R., Karneli, A., O'Neill, N. T., Pietras, C., Pinker, R. T., Voss, K., and Zibordi, G.: An emerging ground-based aerosol climatology: Aerosol optical depth from AERONET, J. Geophys. Res., 106, 12067–12098,, 2001. a

Johnson, M., Meskhidze, N., Solmon, F., Gassó, S., Chuang, P., Gaiero, D., Yantosca, R., Wu, S., Wang, Y., and Carouge, C.: Modeling dust and soluble iron deposition to the South Atlantic Ocean, J. Geophys. Res., 115, D15202,, 2010. a

Justice, C. O., Giglio, L., Roy, D., Boschetti, L., Csiszar, I., Davies, D., Korontzi, S., Schroeder, W., O'Neal, K., and Morisette, J.: MODIS-Derived Global Fire Products, in: Land Remote Sensing and Global Environmental Change: NASA's Earth Observing System and the Science of ASTER and MODIS, Remote Sensing and Digital Image Processing, edited by: Ramachandran, B., Justice, C. O., and Abrams, M. J., Springer, 661–679,, 2011. a

Kanitz, T., Seifert, P., Ansmann, A., Engelmann, R., Althausen, D., Casiccia, C., and Rohwer, E. G.: Contrasting the impact of aerosols at northern and southern midlatitudes on heterogeneous ice formation, Geophys. Res. Lett., 38, L17802,, 2011. a, b, c

Kanitz, T., Ansmann, A., Engelmann, R., and Althausen, D.: North-south cross sections of the vertical aerosol distribution over the Atlantic Ocean from multiwavelength Raman/polarization lidar during Polarstern cruises, J. Geophys. Res.-Atmos., 118, 2643–2655,, 2013. a, b, c, d, e, f, g

Kanitz, T., Ansmann, A., Foth, A., Seifert, P., Wandinger, U., Engelmann, R., Baars, H., Althausen, D., Casiccia, C., and Zamorano, F.: Surface matters: limitations of CALIPSO V3 aerosol typing in coastal regions, Atmos. Meas. Tech., 7, 2061–2072,, 2014. a

Küchler, N., Kneifel, S., Löhnert, U., Kollias, P., Czekala, H., and Rose, T.: A W-Band Radar–Radiometer System for Accurate and Continuous Monitoring of Clouds and Precipitation, J. Atmos. Ocean. Tech., 34, 2375–2392,, 2017. a

Li, F., Ginoux, P., and Ramaswamy, V.: Transport of Patagonian dust to Antarctica, J. Geophys. Res., 115, D18217,, 2010. a

Liu, D., Wang, Z., Liu, Z., Winker, D., and Trepte, C.: A height resolved global view of dust aerosols from the first year CALIPSO lidar measurements, J. Geophys. Res., 113, D16214,, 2008. a

Liu, Z., Vaughan, M., Winker, D., Kittaka, C., Getzewich, B., Kuehn, R., Omar, A., Powell, K., Trepte, C., and Hostetler, C.: The CALIPSO Lidar Cloud and Aerosol Discrimination: Version 2 Algorithm and Initial Assessment of Performance, J. Atmos. Ocean. Tech., 26, 1198–1213,, 2009. a, b

Lopes, F. J. S., Landulfo, E., and Vaughan, M. A.: Evaluating CALIPSO's 532 nm lidar ratio selection algorithm using AERONET sun photometers in Brazil, Atmos. Meas. Tech., 6, 3281–3299,, 2013. a

Mace, G. G. and Protat, A.: Clouds over the Southern Ocean as Observed from the R/V Investigator during CAPRICORN. Part I: Cloud Occurrence and Phase Partitioning, J. Appl. Meteorol. Clim., 57, 1783–1803,, 2018a. a

Mace, G. G. and Protat, A.: Clouds over the Southern Ocean as Observed from the R/V Investigator during CAPRICORN. Part II: The Properties of Nonprecipitating Stratocumulus, J. Appl. Meteorol. Clim., 57, 1805–1823,, 2018b. a

Mamouri, R. E. and Ansmann, A.: Estimated desert-dust ice nuclei profiles from polarization lidar: methodology and case studies, Atmos. Chem. Phys., 15, 3463–3477,, 2015. a, b

Mattis, I., Müller, D., Ansmann, A., Wandinger, U., Preißler, J., Seifert, P., and Tesche, M.: Ten years of multiwavelength Raman lidar observations of free-tropospheric aerosol layers over central Europe: Geometrical properties and annual cycle, J. Geophys. Res., 113, D009636,, 2008. a, b

Minikin, A., Petzold, A., Ström, J., Krejci, R., Seifert, M., van Velthoven, P., Schlager, H., and Schumann, U.: Aircraft observations of the upper tropospheric fine particle aerosol in the Northern and Southern Hemispheres at midlatitudes, Geophys. Res. Lett., 30, 1503,, 2003. a

Moulin, C., Dulac, F., Lambert, C. E., Chazette, P., Jankowiak, I., Chatenet, B., and Lavenu, F.: Long-term daily monitoring of Saharan dust load over ocean using Meteosat ISCCP-B2 data 2. Accuracy of the method and validation using Sun photometer measurements, J. Geophys. Res., 102, 16959–16970,, 1997. a

Müller, D., Ansmann, A., Mattis, I., Tesche, M., Wandinger, U., Althausen, D., and Pisani, G.: Aerosol-type-dependent lidar ratios observed with Raman lidar, J. Geophys. Res., 112, D16202,, 2007. a, b, c, d, e

Naud, C. M., Booth, J. F., and Del Genio, A. D.: Evaluation of ERA-Interim and MERRA Cloudiness in the Southern Ocean, J. Climate, 27, 2109–2124,, 2014. a

Omar, A. H., Winker, D. M., Kittaka, C., Vaughan, M. A., Liu, Z., Hu, Y., Trepte, C. R., Rogers, R. R., Ferrare, R. A., Lee, K.-P., Kuehn, R. E., and Hostetler, C. A.: The CALIPSO Automated Aerosol Classification and Lidar Ratio Selection Algorithm, J. Atmos. Ocean. Tech., 26, 1994–2014,, 2009. a

Preißler, J., Wagner, F., Pereira, S. N., and Guerrero-Rascado, J. L.: Multi-instrumental observation of an exceptionally strong Saharan dust outbreak over Portugal, J. Geophys. Res.-Atmos., 116, D24204,, 2011. a

Protat, A., Schulz, E., Rikus, L., Sun, Z., Xiao, Y., and Keywood, M.: Shipborne observations of the radiative effect of Southern Ocean clouds, J. Geophys. Res.-Atmos., 122, 318–328,, 2016. a

Radenz, M. and Seifert, P.: Software for automated trajectory analysis: trace (Version v0.3),, 2019. a, b

Raes, F., Bates, T., McGovern, F., and van Liedekerke, M.: The 2nd Aerosol Characterization Experiment (ACE2): general overview and main results, Tellus B, 52, 111–125,, 2000. a

Reid, J. S., Prins, E. M., Westphal, D. L., Schmidt, C. C., Richardson, K. A., Christopher, S. A., Eck, T. F., Reid, E. A., Curtis, C. A., and Hoffman, J. P.: Real-time monitoring of South American smoke particle emissions and transport using a coupled remote sensing/box-model approach, Geophys. Res. Lett., 31, L06107,, 2004. a

Ristori, P., Otero, L., Jin, Y., Barja, B., Shimizu, A., Barbero, A., Salvador, J., Bali, J. L., Herrera, M., Etala, P., Acquesta, A., Quel, E., Sugimoto, N., and Mizuno, A.: lidar network in southern South America, EPJ Web Conf., 176, 09011,, 2018. a

Rolph, G., Stein, A., and Stunder, B.: Real-time Environmental Applications and Display sYstem: READY, Environ. Modell. Softw., 95, 210–228,, 2017. a

Russell, P. B., Hobbs, P. V., and Stowe, L. L.: Aerosol properties and radiative effects in the United States East Coast haze plume: An overview of the Tropospheric Aerosol Radiative Forcing Observational Experiment (TARFOX), J. Geophys. Res., 104, 2213–2222,, 1999. a

Schneider, C., Glaser, M., Kilian, R., Santana, A., Butorovic, N., and Casassa, G.: Weather observations across the southern Andes at 53 S, Phys. Geogr., 24, 97–119, 2003. a

Schuster, G. L., Vaughan, M., MacDonnell, D., Su, W., Winker, D., Dubovik, O., Lapyonok, T., and Trepte, C.: Comparison of CALIPSO aerosol optical depth retrievals to AERONET measurements, and a climatology for the lidar ratio of dust, Atmos. Chem. Phys., 12, 7431–7452,, 2012. a

Seifert, P., Ansmann, A., Mattis, I., Wandinger, U., Tesche, M., Engelmann, R., Müller, D., Pérez, C., and Haustein, K.: Saharan dust and heterogeneous ice formation: Eleven years of cloud observations at a central European EARLINET site, J. Geophys. Res., 115, D20201,, 2010. a

Smirnov, A., Holben, B. N., Slutsker, I., Giles, D. M., McClain, C. R., Eck, T. F., Sakerin, S. M., Macke, A., Croot, P., Zibordi, G., Quinn, P. K., Sciare, J., Kinne, S., Harvey, M., Smyth, T. J., Piketh, S., Zielinski, T., Proshutinsky, A., Goes, J. I., Nelson, N. B., Larouche, P., Radionov, V. F., Goloub, P., Krishna Moorthy, K., Matarrese, R., Robertson, E. J., and Jourdin, F.: Maritime Aerosol Network as a component of Aerosol Robotic Network, J. Geophys. Res., 114, D06204,, 2009. a

Stein, A. F., Draxler, R. R., Rolph, G. D., Stunder, B. J. B., Cohen, M. D., and Ngan, F.: NOAA's HYSPLIT Atmospheric Transport and Dispersion Modeling System, B. Am. Meteorol. Soc., 96, 2059–2077,, 2015. a, b

Stohl, A., Forster, C., Frank, A., Seibert, P., and Wotawa, G.: Technical note: The Lagrangian particle dispersion model FLEXPART version 6.2, Atmos. Chem. Phys., 5, 2461–2474,, 2005. a

Stull, R. B.: An introduction to boundary layer meteorology, Kluwer Academic Puplishers, Wisconsin, USA, 1988. a

Tesche, M., Müller, D., Gross, S., Ansmann, A., Althausen, D., Freudenthaler, V., Weinzierl, B., Veira, A., and Petzold, A.: Optical and microphysical properties of smoke over Cape Verde inferred from multiwavelength lidar measurements, Tellus B, 63, 677–694,, 2011. a, b

Vaughan, M., Pitts, M., Trepte, C., Winker, D., Detweiler, P., Garnier, A., Getzewitch, B., Hunt, W., Lambeth, J., Lee, K.-P., Lucker, P., Murray, T., Rodier, S., Trémas, T., Bazureau, A., and Pelon, J.: V4.10 CALIPSO CALIOP 5 km aerosol profile product,, 2016. a

Weinzierl, B., Ansmann, A., Prospero, J. M., Althausen, D., Benker, N., Chouza, F., Dollner, M., Farrell, D., Fomba, W. K., Freudenthaler, V., Gasteiger, J., Groß, S., Haarig, M., Heinold, B., Kandler, K., Kristensen, T. B., Mayol-Bracero, O. L., Müller, T., Reitebuch, O., Sauer, D., Schäfler, A., Schepanski, K., Spanu, A., Tegen, I., Toledano, C., and Walser, A.: The Saharan Aerosol Long-Range Transport and Aerosol–Cloud-Interaction Experiment: Overview and Selected Highlights, B. Am. Meteorol. Soc., 98, 1427–1451,, 2017. a

Wilson, S. R. and Forgan, B. W.: Aerosol optical depth at Cape Grim, Tasmania, 1986–1999, J. Geophys. Res.-Atmos., 107, 4068,, 2002. a

Winker, D. M., Hunt, W. H., and McGill, M. J.: Initial performance assessment of CALIOP, Geophys. Res. Lett., 34, L19803,, 2007. a

Winker, D. M., Vaughan, M. A., Omar, A., Hu, Y., Powell, K. A., Liu, Z., Hunt, W. H., and Young, S. A.: Overview of the CALIPSO Mission and CALIOP Data Processing Algorithms, J. Atmos. Ocean. Tech., 26, 2310–2323,, 2009. a

Publications Copernicus
Short summary
In this study, we present the vertical aerosol distribution in the pristine region of the southern tip of South America determined by ground-based and spaceborne lidar observations. Most aerosol load is contained within the planetary boundary layer up to about 1200 m. The free troposphere is characterized by a very low aerosol concentration but a frequent occurrence of clouds. Lofted aerosol layers were rarely observed and, when present, were characterized by very low optical thicknesses.
In this study, we present the vertical aerosol distribution in the pristine region of the...