Transport and Variability of Tropospheric Ozone over Oceania and Southern Paciﬁc during the 2019–20 Australian Bushﬁres

: The present study contributes to the scientiﬁc effort for a better understanding of the potential of the Australian biomass burning events to inﬂuence tropospheric trace gas abundances at the regional scale. In order to exclude the inﬂuence of the long-range transport of ozone precursors from biomass burning plumes originating from Southern America and Africa, the analysis of the Australian smoke plume has been driven over the period December 2019 to January 2020. This study uses satellite (IASI, MLS, MODIS, CALIOP) and ground-based (sun-photometer, FTIR, ozone radiosondes) observations. The highest values of aerosol optical depth (AOD) and carbon monoxide total columns are observed over Southern and Central Australia. Transport is responsible for the spatial and temporal distributions of aerosols and carbon monoxide over Australia, and also the transport of the smoke plume outside the continent. The dispersion of the tropospheric smoke plume over Oceania and Southern Paciﬁc extends from tropical to extratropical latitudes. Ozone radiosonde measurements performed at Samoa (14.4 ◦ S, 170.6 ◦ W) and Lauder (45.0 ◦ S, 169.4 ◦ E) indicate an increase in mid-tropospheric ozone (6–9 km) (from 10% to 43%) linked to the Australian biomass burning plume. This increase in mid-tropospheric ozone induced by the transport of the smoke plume was found to be consistent with MLS observations over the tropical and extratropical latitudes. The smoke plume over the Southern Paciﬁc was organized as a stretchable anticyclonic rolling which impacted the ozone variability in the tropical and subtropical upper-troposphere over Oceania. This is corroborated by the ozone proﬁle measurements at Samoa which exhibit an enhanced ozone layer (29%) in the upper-troposphere. Our results suggest that the transport of Australian biomass burning plumes have signiﬁcantly impacted the vertical distribution of ozone in the mid-troposphere southern tropical to extratropical latitudes during the 2019–20 extreme Australian bushﬁres.


Introduction
Biomass burning is a significant source of aerosols and active trace gases that impact the atmospheric composition from local to global scales, with consequences for tropospheric ozone production and budget. Tropospheric ozone presents a serious environmental issue due to its negative impact on human health and Earth's climate. Tropospheric ozone is predominantly produced by photochemical reactions involving precursors generated by natural processes, and to a much larger extent, by anthropogenic activities [1]. Previous works have showed that large amounts of ozone precursors such as carbon monoxide (CO) are injected throughout the troposphere during the biomass burning activity [2][3][4]. Generally, ozone radiative forcing depends on its vertical distribution. In the troposphere, ozone is a greenhouse gas that plays a key role in radiative forcing and indeed in climate change [5]. Aerosols or photochemical compounds produced during biomass burning events have a large radiative effect on climate but are still scarcely quantified in the literature. This is especially true in the southern hemisphere due to the low number of observational stations, in comparison with the northern hemisphere.
The bulk of biomass burning activity in the southern hemisphere occurs in the Southern American and Southern African regions [6][7][8] from July to November. High concentrations of ozone precursors and aerosols from these fires can be transported far from the source regions [9]. Based on satellite observations and ground-based measurements at Lauder (45.0 • S, 169.4 • E), Edwards et al. [10] showed that intercontinental transport of biomass burning pollution from Africa and South America often determines the local air-quality over this site. They showed that the fraction of the pollutants following this pathway is estimated to be less than 10% of the total emission in South America and Africa. By combining 10 years of ozone radiosondes performed at several Tropical Pacific sites and back-trajectory modeling, Oltmans et al. [11] pointed out the contribution of intercontinental transport of biomass burning plume from Africa and South America on tropospheric ozone variability over the Southern Pacific between September and November. Oltmans et al. [11] did not exclude that Australian and Indonesian biomass burning activity may have also contributed to the tropospheric ozone variability over Oceania and the Southern Pacific. Pak et al. [12] reported aircraft missions conducted over southeastern Australia and Tasmania in parallel with the Southern African Regional Science Initiative (SAFARI 2000) campaign in Africa. These large-scale studies in Africa and South America have provided valuable insight into the impact of the long-range transport of the smoke plume on the composition of the atmosphere over Oceania. Despite these studies, there is still a need for a better scientific understanding of the influence of biomass burning on atmospheric composition and air-quality [13], particularly around Australia.
The vast majority of the Australian biomass burning activity takes place in the northern tropical part of the continent, where the ecosystem is dominated by savannah. These fires are responsible for more than half of the land area burned between September and January [14]. Previous works revealed that emissions from fires in the savannah regions of Northern Australia can be considered as the most significant regional source of greenhouse and other trace gases, as well as atmospheric aerosols [15][16][17]. Nevertheless, the intensity of these fires is low to moderate, generally limited to the grass layer [17][18][19]. The propagation of the smoke plumes from these fires is mainly limited to the planetary boundary layer and to the northwest over the Timon Sea [15,20]. Although the vast majority of the continent's fires occur in the Northern region of Australia, the most devastating fires burn in the densely populated Southern regions. Southeastern Australia has a combination of climate conditions and vegetation type that makes it prone to severe wildfires [21]. Sullivan et al. [22] reported that extreme fire events occurred in most recent years over southeastern Australia but were generally more extensive and severe following extended drought, typically associated with El Niño events. The extreme fire events occurring over southeastern Australia can trigger the development of convection, which leads to the lifting of aerosols and trace gases to high altitudes. Moist convection can sometimes occur in the fire plume and lead to the formation of convective clouds (pyro-cumulus) Remote Sens. 2021, 13, 3092 3 of 26 or thunderstorms (pyro-cumulonimbus) [23]. These phenomena are widely known as pyro-convection events [24]. The vertical transport of the smoke plume is, therefore, due to a combination of extreme heat energy release from the fire and convection. Under favorable meteorological conditions, pyro-convection events have the potential to inject soot and smoke directly into the stratosphere. Southeastern Australia has a long history of severe pyro-convection events which have significantly impacted the composition of the stratosphere at regional and global scales. In the past two decades, pyro-convection events have impacted southeastern Australia in 2003 [24][25][26][27], 2006 [28,29] and 2009 (colloquially known as "Black Saturday") [30][31][32][33]. Recently, this list has grown as a result of a new pyro-convection event occurring during the 2019-20 fire season in eastern and southeastern Australia (colloquially known as "Black Summer") [34][35][36][37].
Nolan et al. [38] reported that the 2019-20 fire season in eastern and southeastern Australia was unprecedented in the size and number of fires. Boer et al. [39] estimated that the area burnt during this fire season was about six million hectares. It was also reported that Eastern Australia experienced severe drought before and during this fire season (www.bom.gov.au/climate/maps/ (accessed on 2 August 2021) [38]. The meteorological conditions were favorable to extreme fires leading to pyro-convection events between New Year's Eve and mid-January. Khaykin et al. [36] reported that the 2019-20 Australian wildfires caused a massive injection of combustion products into the stratosphere that led to a persistent planetary-scale perturbation of all stratospheric climate-relevant variables. The perturbation induced by the 2019-20 Australian wildfires in the stratosphere put this event on a par with the strongest volcanic eruptions of the last 25 years [35,36,40]. The gases and aerosols injected into the stratosphere were advected by westerly winds and dispersed across all the southern hemisphere's extratropical latitudes [41]. The smoke plume returned back to Australia in less than two weeks by passing over South America [37]. Based on satellite observations and European Centre for Medium-Range Weather Forecast (ECMWF) operational data, Khaykin et al. [36] revealed that this extreme event enabled the discovery of a striking atmospheric phenomenon, a persistent synoptic-scale anticyclone, which acts to confine the carbon-rich aerosol clouds during their solar-driven rise. In this way, highlyconcentrated absorbing aerosols are lofted into the middle stratosphere, which prolongs their stratospheric residence time and radiative effects on climate. Most of these previous works report on the investigation on the perturbation of stratospheric composition, and on the dynamical circulation induced by the 2019-20 Australian wildfires.
In contrast, relatively little attention has been paid to smoke emissions in the troposphere during these extreme wildfires. Few studies on the influence of the smoke emissions from southeastern and eastern Australia on the tropospheric variability of aerosols and trace gases were undertaken [11,40,42,43]. This is fairly surprising in view of the potential of biomass burning activity to change the physico-chemical composition in the troposphere at the local and regional scale [44,45]. Based on ozone observations from LiDAR and radiosondes over Reunion Island, Clain et al. [44] showed a significant annual increase in tropospheric ozone over Reunion Island during the August-November period, in parallel with the biomass burning season in Southern Africa and Madagascar. High concentrations of ozone precursors from these fires are vented into the free troposphere by convection and are subsequently advected into the Indian Ocean region by westerly winds.
The present work reports on the characteristics of the tropospheric smoke plumes originating from southeastern and eastern Australian during the 2019-20 fire season. It aims to investigate the tropospheric ozone variability over Oceania and the Southern Pacific induced by the transport of the Australian smoke plume during the 2019-20 fire season. The chemical and optical properties of the smoke plume are mainly inferred from ground-based observations collected in the framework of international networks. Even though Kloss et al. [40] showed that the Australian fire season of 2019-20 has already significantly impacted fire trace gas abundances in the southern hemisphere's troposphere starting from September/October 2019, here we focus on the peak phase of the fire events (From December 2019 to January 2020). This is carried out to exclude the possible impact Remote Sens. 2021, 13, 3092 4 of 26 of ozone precursors from biomass burning plumes originating from Southern America and Africa [10].
The paper is organized as follows: Section 2 describes the ground-based and spaceborne observations used for the investigation of the tropospheric smoke plume. A review of the spatio-temporal evolution of the CO and aerosols' total column over Australia is provided in Section 3. An analysis of the tropospheric smoke plume is also provided in Section 3. The discussion on the tropospheric ozone variability over Oceania and Southern Pacific is presented in Section 4. A summary and the perspectives of this study are finally given in Section 5.

Instruments and Data Sets
The characteristics of the fire emissions were investigated using all available observations of trace gases (mainly CO and ozone) and aerosols over Australia and its surroundings. The data set used in this study consists of ground-based and space-borne observations. Given the lifetime of CO (~1 to 3 months in the troposphere), this latter is often used as tracer to study the long-range transport of biomass burning [4,46,47]. CO concentration is associated with many tropospheric polluting gases such as tropospheric ozone. To investigate the tropospheric ozone variability induced by the Australian fires, it is important to have an idea about the background situation when the atmosphere is not undergoing major disturbances (e.g, volcanic eruptions, pyro-convection outbreaks). In the present work, the background period is defined as the period ranging from January 2010 to December 2018. The previously reported pyro-convection events in southeastern Australia were discarded from the background values.

FTIR
Total column and volume mixing ratio profiles of trace gases such as CO and ozone were retrieved with high accuracy and precision from ground-based Fourier transform infrared (FTIR) spectrometers [48,49]. In the present study, we use FTIR ozone and CO observations taken at Lauder (New Zealand). These measurements have been taken since the early 1990s within the framework of the NDACC (Network for the Detection of Atmospheric Composition Change) network [50]. The measurements are taken on clear-sky days, throughout the day, using a Bruker high-resolution spectrometer over a wide spectral range (around 600-4500 cm −1 ). A detailed description of the ozone and CO retrieval algorithms is given by Vigouroux et al. [48] and Morgenstern et al. [51]; they are briefly summarized here. The CO and ozone total column amounts and volume mixing ratio profiles are retrieved using the spectral least-squares fitting (SFIT) algorithm [52,53]. This algorithm is based on the optimal estimation method developed by Rodgers [54]. FTIR retrievals are performed on a 48-layer atmosphere (0.37-100 km asl) of which 15 layers are between 0.37 and 10 km. The retrievals use a static a priori originating from WACCM 4 (Whole Atmosphere Community Climate Model, version 4) model simulations [55]. The retrievals are constrained using regularization criteria (a priori covariance matrix and spectrum signal to noise ratio). This regularization step is taken in order to obtain stable physically realistic retrievals with reasonable degrees of freedom (DoFs) for the signal. Across the whole measurement period at Lauder, the average DoFs value in the troposphere was 1.50 and 2.31 for ozone and CO, respectively.
The CO and ozone total columns and volume mixing ratio profiles used in this study were downloaded from the NDACC database (http://www.ndacc.org, accessed on 21 June 2021).

Ozone Radiosondes
The ozonesondes provide ozone profiles from ground up to burst-balloon altitude, which is around 30-35 km, with a vertical resolution of 150 m. The dataset used in this A detailed description of the CIMEL sunphotometer of the AERONET network and the associated data retrieval is given by Holben et al. [73] and it is summarized hereafter. The direct solar extinction and diffuse sky radiance measurements are used to compute AOD and to retrieve aerosol size-distribution using the methodology of Dubovik and King [74]. The measurements are performed under cloud-free and daytime conditions at 15 min intervals. Previous works showed that estimated uncertainty in AOD measurements under cloud free conditions ranges from 0.01 to 0.02 [75]. The AOD values presented in this work are selected at Level 1.5 (Version 3) since level 2 data were not yet available at the time of this study. Giles et al. [76] showed that the new level 1.5 (version 3) AOD data are statistically very close to version 2 level 2.0. The level 1.5 of AOD data with version 3 refers now to near real-time automatic cloud screening and automatic instrument anomaly quality controls. We also note that the new level 1.5 (version 3) AOD data have been validated and used in recent works [48,77,78]. The AOD data used in the present study were downloaded from the AERONET web site: http://aeronet.gsfc.nasa.gov (accessed on 2 August 2021).

CALIOP
Observations from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) onboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) were used to study the transport of the smoke aerosol layers. CALIPSO flies in a sun-synchronous polar orbit since 2006 with a cycle of 16 days. A detailed description of CALIPSO is given by Winker et al. [79]. In the present work, we also used the CALIOP version 4.1 AOD products level 2 (80,81). Young et al. [80] reported extensive upgrades to the algorithms used to retrieve the extinction and AOD with CALIOP version 4.1 data products. A validation study of the CALIOP 4.1 AOD data products with AERONET and MODIS observations was reported by Kim et al. [81]. They showed that the aforementioned AOD discrepancies between CALIOP and instruments are reduced in version 4 The ozonesonde measurements at Samoa were initiated in 1989 in the framework of several measurement programs, such as Pacific Exploratory Mission [11]. Samoa has been included in the SHADOZ (Southern Hemisphere ADditional Ozone sondes) network since 1998, and the frequency of radiosonde experiments increased to become weekly. Through the SHADOZ project a systematic effort was made to collect regular ozone profile data from stations throughout the tropical and subtropical regions of the southern hemisphere. Details concerning the SHADOZ program are available from Thompson et al. [56]. The balloon-borne sonde equipment used at Samoa consists of an electrochemical concentration cell (ECC) ozonesonde which measures the oxidation of a potassium iodine (KI) solution by ozone [57]. This chemical reaction allows the retrieval of the ozone partial pressure profile. The ECC ozonesonde is coupled with a standard radiosonde to obtain meteorological parameters. Thompson et al. [58] have reported that ECC sondes have an absolute accuracy of 5%. In order to check the quality of the data, the Samoa ozonesonde system was part of the JOSIE (Julich Ozone Sonde Intercomparison Experiment) process. This intercomparison campaign, which was under the World Meteorology Organization (WMO) initiative, aimed to assess the performance of radiosonde systems through an intercomparison experiment using a standard reference instrument [59]. A total of 259 ozone profiles performed at Samoa from 2010 to 2020 were considered in the present study, and were downloaded from the SHADOZ database: https://tropo.gsfc.nasa.gov/shadoz/ (accessed on 21 June 2021).
The ozone radiosondes at Lauder have been routinely performed in the framework of the NDACC network since 1986. At Lauder, the balloon-sonde is also equipped with an ECC ozonesonde operating with 0.5% buffered KI cathode solution [60]. Bodeker et al. [61] compared the integrated ozone profile with the total column of ozone measured by colocated Dobson spectrophotometer. They found less than 5% discrepancy. The Lauder data set used in this study consists of 522 ozone profiles during the same period of time,

Satellite Date MLS
The Microwave Limb Sounder (MLS) instrument onboard Aura has been operating since August 2004 [62]. It is a thermal-emission microwave limb sounder that measures mixing ratio profiles from the upper-troposphere to the mesosphere of a number of chemical species (e.g, CO, ozone, water vapour) from limb scans taken in the direction ahead of the Aura satellite orbital track. The orbit is sun-synchronous, so that the observations are always made at the same two local times for a given latitude. The output of the radiometers is analyzed by banks of filters. The ozone profiles are retrieved using radiances at 240 GHz. The radiances reported by the filter banks are used as input to a software package [63] which estimates profiles of temperature, and profiles of the mixing ratios of the targeted chemical species.
In the present work, we used the MLS version 4.2 Ozone Mixing Ratio (OMR). All MLS v4.2 retrieval quality flags (quality, status, convergence, and precision) were properly adhered to for all of our analyses. The MLS v4.2 measurements, including data quality and quality flags, are described in the MLS data quality document (http://mls.jpl.nasa.gov/ data/v4-2_data_quality_document.pdf, accessed on 6 June 2021) and by Livesey et al. [63]. Generally, recommended pressure levels for science applications with ozone MLS data range from 0.0215 to 261 hPa. A validation study of the MLS v4.2 OMR with ozone radiosonde over the Tibetan Plateau was undertaken by Yan et al. [64]. Their results suggest that MLS v4.2 represents a substantial improvement in ozone retrievals throughout the upper troposphere in comparison with the previous versions (accuracy of 20 ppbv at 215 hPa and 30 ppbv at 100 hPa). The OMR profiles from MLS were obtained from the Atmospheric Composition Data and Information Services Center (ACDISC) archive (ftp://acdisc.gsfc.nasa.gov, accessed on 21 June 2021) hosted by the NASA Goddard Space Flight Center.

IASI
The Infrared Atmospheric Sounding Interferometer (IASI) onboard the Meteorological Operational (MetOp) satellite measures chemical species (e.g., ozone, sulfur dioxide, CO) on a global scale [65][66][67]. The IASI atmospheric sounder uses a Fourier transform spectrometer, which measures the infrared radiation emitted by the Earth's surface and the atmosphere, covering wavelengths from 6.62 to 15.5 µm, allowing for day and night measurements. It was designed for global observations with a vertical range covering the troposphere and the lower stratosphere.
Retrieval of CO and O 3 total and partial columns occur in near real-time from the nadir radiance spectra measured by the instrument in the thermal infrared (TIR) spectral range. Global distributions are obtained at around 9:30 AM and PM (local time) using a dedicated radiative transfer and retrieval software for the IASI products: The Fast-Optimal Retrievals on Layers for IASI (FORLI-CO and FORLI-O3). A full description is given by Hurtmans et al. [68], and a cloud filtering is applied to the data set before further exploitation. Global distributions of CO columns were evaluated with correlative observations available from other nadir looking TIR missions. On average, George et al. (2009) found total column discrepancies of about 7% in comparison with three other sounders. The valuable nature of CO total column observations from IASI for the study of long-range transport of the smoke plume during biomass burning events is well documented in previous research [47,69]. IASI tropospheric ozone from FORLI-O3 has been extensively validated and compared with available ground-based, aircraft, ozone radiosondes, and other satellite observations [69][70][71]. These analyses show that IASI underestimates tropospheric ozone by 12% to 14% in the mid-latitudes and tropics. In the present study, we used CO and ozone columns from IASI instruments on Metop-A and Metop-B, which have Remote Sens. 2021, 13, 3092 7 of 26 been operating since 2006 and 2012, respectively. The IASI products used in this work are available on the AERIS platform: https://iasi.aeris-data.fr/ (accessed on 27 July 2021).

Sun-Photometer CIMEL
The ground-based sites selected to investigate the aerosol variability over Australia include 6 AERONET (Aerosols Robotic NETwork) sites strategically deployed over the continent ( Figure 1 [72]. A detailed description of the CIMEL sunphotometer of the AERONET network and the associated data retrieval is given by Holben et al. [73] and it is summarized hereafter. The direct solar extinction and diffuse sky radiance measurements are used to compute AOD and to retrieve aerosol size-distribution using the methodology of Dubovik and King [74]. The measurements are performed under cloud-free and daytime conditions at 15 min intervals. Previous works showed that estimated uncertainty in AOD measurements under cloud free conditions ranges from 0.01 to 0.02 [75]. The AOD values presented in this work are selected at Level 1.5 (Version 3) since level 2 data were not yet available at the time of this study. Giles et al. [76] showed that the new level 1.5 (version 3) AOD data are statistically very close to version 2 level 2.0. The level 1.5 of AOD data with version 3 refers now to near real-time automatic cloud screening and automatic instrument anomaly quality controls. We also note that the new level 1.5 (version 3) AOD data have been validated and used in recent works [48,77,78]. The AOD data used in the present study were downloaded from the AERONET web site: http://aeronet.gsfc.nasa.gov (accessed on 2 August 2021).

Satellite Data CALIOP
Observations from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) onboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) were used to study the transport of the smoke aerosol layers. CALIPSO flies in a sunsynchronous polar orbit since 2006 with a cycle of 16 days. A detailed description of CALIPSO is given by Winker et al. [79]. In the present work, we also used the CALIOP version 4.1 AOD products level 2 (80, 81). Young et al. [80] reported extensive upgrades to the algorithms used to retrieve the extinction and AOD with CALIOP version 4.1 data products. A validation study of the CALIOP 4.1 AOD data products with AERONET and MODIS observations was reported by Kim et al. [81]. They showed that the aforementioned AOD discrepancies between CALIOP and instruments are reduced in version 4 compared to version 3. The CALIOP version 4 AOD data product has been used previously for optical characterization of aerosol plume over biomass burning or dust sources [82][83][84]. The CALIOP data were obtained from the ACDISC data archive (ftp://acdisc.gsfc.nasa.gov, accessed on 16 July 2021) hosted by NASA Goddard Space Flight Center.

GOME-2
The GOME-2 (Global Ozone Monitoring Experiment) onboard the MetOp satellites is a nadir viewing scanning spectrometer which provides global coverage of the sunlit part of the atmosphere with a spatial coverage of 80 × 40 km [85]. The absorbing aerosol index (AAI) data from GOME-2 onboard MetOp-B were used to describe the transport of the Australian biomass burning plume. This index allows the detection of absorbing aerosols Remote Sens. 2021, 13, 3092 8 of 26 through the spectral difference between 340 nm and 380 nm. The AAI index has provided information about smoke plume horizontal distribution over land and ocean surface in many studies [25,28,31]. The AAI data were downloaded from: https://www.temis.nl/ airpollution/absaai/ (accessed on 21 June 2021).

MODIS Fire Data
The location of fires is derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensors onboard the Aqua and Terra platforms launched by NASA in 1999 and 2002, respectively. Previous works confirm that MODIS can be considered as one of the primary reference instruments for fire detection and monitoring [15,47,86]. MODIS performs measurements at 36 different bands ranging from 0.4 to 14.4 µm, and has infrared channels specifically designed for fire observations. Aqua and Terra satellites have sun-synchronous orbits. Owing to its wide ground scan (2330 km) capability, the MODIS instrument takes one to two days to observe each point on the Earth surface, with a spatial resolution ranging from 0.25 km to 1.00 km, depending on the channels used. The orbit of Aqua and Terra is timed so that it passes over the equator from north to south in the morning. Kaufman et al. [87] reported that the principle of fire detection by MODIS is based on a comparison of the brightness temperatures of two channels, i.e., 4 µm and 11 µm, recorded for each 1 km 2 pixel. A fire is, hence, attributed to the pixel when the difference between the brightness temperatures of the 4 µm and 11 µm channels are above a certain threshold. The MODIS fire products used in this study to investigate fire distributions across Australia during the 2019-20 fire season were downloaded from the NASA Earth Data platform: https: //earthdata.nasa.gov/earth-observation-data/near-real-time/firms/active-fire-data (accessed on 10 June 2021).

CO and Aerosols Total Columns over Australia
Southeastern and eastern Australia was impacted by a severe heat wave on December 2019 [35,88]. This heat wave led to drought conditions in New South Wales, Victoria and Queensland. At the beginning of December, fires were located over the northwestern side of Australia, and moved to the eastern and southeastern side from mid-December. The daily location of fire pixels detected from MODIS and daily OLR (Outgoing Longwave Radiation) anomalies from NCEP analysis on 19 and 31 December are depicted on Figure 1a,b, respectively. On 19 December, fires are mainly located over eastern Australia within a region of positive OLR anomalies (greater than 15%), as is illustrated on Figure 1a. The prolonged drought and drying of the land may have increased the risk of extensive fires, as suggested in previous similar published works [25,31]. On 31 December, the maximum intensity of fire pixels was concentrated over southeastern Australia (Figure 1b). Figure 1b reveals that negative OLR anomalies (greater than 15%, in absolute value) are observed over southeastern Australia. It is conceivable that an intensification of the fire in this region has forced the formation of deep convection [25]. We note that the most intense fires were located over southeastern Australia with FRP (Fire Radiation Power) values ranging from 1500 to 2500 MW (particularly in the state of Victoria) between 30 December 2019 and 12 January 2020. Conversely, the northern region (between 10 • S and 22 • S) is not impacted by the intense fires (FRP values less than 100 MW) from 1 December 2019 to 31 January 2020.
The latitudinal distribution and evolution of the total column of carbon monoxide (TCO) obtained from IASI observations from 1 December 2019 to 31 January 2020 are presented in Figure 2. The latitude-time Hovmöller diagram is averaged for all longitudes covering the Australian continent (between 114 • E and 155 • E). Figure 2 shows that the TCO is relatively homogeneously distributed with latitude before 15 December 2019. Conversely, an increase in TCO (higher than 5 × 10 18 molecules·cm −2 ) at higher latitudes can be observed between 15 December 2019 and 12 January 2020. The intensification of the CO emission is observed during the period aforementioned. The highest TCO values, with zonal mean values exceeding 7 × 10 18 molecules·cm −2 , are observed over the southern region between 29 December 2019 and 12 January 2020. This is consistent with previous studies which suggested that significant biomass burning emission may originate in Australia's temperate forests, which are located in the southern region [89,90].
The latitudinal distribution and evolution of the total column of carbon monoxide (TCO) obtained from IASI observations from 1 December 2019 to 31 January 2020 are presented in Figure 2. The latitude-time Hovmöller diagram is averaged for all longitudes covering the Australian continent (between 114°E and 155°E). Figure 2 shows that the TCO is relatively homogeneously distributed with latitude before 15 December 2019. Conversely, an increase in TCO (higher than 5 × 10 18 molecules·cm −2 ) at higher latitudes can be observed between 15 December 2019 and 12 January 2020. The intensification of the CO emission is observed during the period aforementioned. The highest TCO values, with zonal mean values exceeding 7 × 10 18 molecules·cm −2 , are observed over the southern region between 29 December 2019 and 12 January 2020. This is consistent with previous studies which suggested that significant biomass burning emission may originate in Australia's temperate forests, which are located in the southern region [89,90]. The TCO values decrease (less than 3.5 × 10 18 molecules·cm −2 ) over all the latitude band after 12 January 2020. Nevertheless, moderate values (4-4.5 × 10 18 molecules·cm −2 ) are observed over the 35°S-38°S latitudinal band on 14-16 and 25-26 January 2020 (Figure 2).   Figure 1). CALIOP AOD data were collected over an area of 5° × 5° latitude × longitude centered on each selected site. The daily AOD measurements performed between 1 December 2019 and 31 January 2020 were compared to background values (2010-2018) obtained from sun-photometer observations. Background values of AOD decrease with increasing latitude, which is consistent with previous studies [91,92]. Sun-photometer and CALIOP measurements are complementary and, together, show full temporal evolution of the AOD over each site. This is clearly illustrated with the case of the Tumbarumba site, where no measurements were recorded by the sun-photometer after 4 January 2020 (Figure 3f). It is worth noting that the AOD evolution obtained from CALIOP and sun-photometer are in fairly good agreement over all sites. This confirms that the high AOD values (higher than background values) over southern sites (such as Fowlers Gap and   Figure 1). CALIOP AOD data were collected over an area of 5 • × 5 • latitude × longitude centered on each selected site. The daily AOD measurements performed between 1 December 2019 and 31 January 2020 were compared to background values (2010-2018) obtained from sun-photometer observations. Background values of AOD decrease with increasing latitude, which is consistent with previous studies [91,92]. Sun-photometer and CALIOP measurements are complementary and, together, show full temporal evolution of the AOD over each site. This is clearly illustrated with the case of the Tumbarumba site, where no measurements were recorded by the sun-photometer after 4 January 2020 (Figure 3f). It is worth noting that the AOD evolution obtained from CALIOP and sun-photometer are in fairly good agreement over all sites. This confirms that the high AOD values (higher than background values) over southern sites (such as Fowlers Gap and Tumbarumba) did not result from a technical artefact, but were observed by multiple independent instruments. Figure 3 reveals that the daily evolutions of the AOD obtained from sun-photometer and CALIOP over all the selected sites show significant increases (higher than background values) between 15 December 2019 and 12 January 2020. These AOD anomalies appear to increase southward with higher latitude. This is similar to the latitudinal distributions obtained for TCO presented in Figure 2.
Tumbarumba) did not result from a technical artefact, but were observed by multiple independent instruments. Figure 3 reveals that the daily evolutions of the AOD obtained from sun-photometer and CALIOP over all the selected sites show significant increases (higher than background values) between 15 December 2019 and 12 January 2020. These AOD anomalies appear to increase southward with higher latitude. This is similar to the latitudinal distributions obtained for TCO presented in Figure 2.   Figure 3 depicts also the daily mean evolutions of TCO (in red lines) obtained from IASI observations over all sites from 1 December 2019 to 31 January 2020. It can be observed that the daily TCO and AOD evolutions during the aforementioned period are strongly correlated over all sites. This suggests that the significant increase in daily AOD values between 15 December and 12 January 2020 is most likely due to the biomass burning activity. Due to its geographical location, the Tumbarumba site is clearly the most exposed to intense fires ( Figure 1). As a result, the highest differences between daily measurements and the background values are observed over this site (on average with a factor greater than 3, Figure 3f). Conversely, the lowest differences with the background values (on average with a factor less than 2.5) are observed over the northern sites (Figure 3a-c). As discussed previously, the fires over northern Australia were less intense than the fires located over southeastern Australia. A surprising result may be the influence of the biomass burning activity on the evolution of AOD between 15 December 2019 and 12 January 2020 over Birdsville and Fowlers Gap. These two sites are far away from the fire areas, which are mainly located on the Australian coast ( Figure 1). Furthermore, Birdsville and Fowlers Gap are located near Australia's major dust sources [90][91][92]. Mukkavilli et al. [92] revealed that dust is the major type of atmospheric aerosol over central Australia. Nevertheless, we cannot exclude the hypothesis that the significant increase in AOD and TCO over Birdsville and Fowlers Gap is a result of regional transport of biomass burning aerosols and CO from eastern and southeastern Australia. By the use of the HYSPLIT back-trajectory model and sun-photometer observations, Qin and Mitchell [91] found that smoke generated from fires in eastern and southeastern Australia can be transported across New South Wales to central Australia. Yang et al. [93] performed a cluster analysis on simulated HYSPLIT back-trajectories during the biomass burning seasons from 2005 and 2020 at ground level in Australia. They found that 32% of air masses over central Australia come from southeastern Australia during the biomass burning period. Based on CALIOP observations and the HYSPLIT model, Yang et al. [73] found that aerosols from southeastern Australia were transported toward Birdsville during the 18-26 December 2019 period from ground level to 5 km above sea level. The air masses loaded with aerosols and CO can thus be transported to central Australia.
In summary, the recirculation of air masses from southeastern to central Australia is a dynamical scheme frequently observed during the biomass burning period [90,91,93]. This dynamical scheme is characterized by a high-pressure level system developed over the eastern side and a cold front over the southwestern side. The two systems moved eastward, which led to northwesterly flow over southeastern Australia. The cold front acted as a physical barrier and brought warm dry air mass into a small corridor stretched from southeastern Australia toward the Tasman sea [31]. Therefore, this dynamical scheme is responsible for the spatial and temporal distributions of aerosols and chemical species such as CO within Australia, and also the transport of the plume towards the ocean.

Tropospheric Smoke Plume
Air masses loaded with aerosols and chemical species are frequently transported toward the Tasman Sea during the biomass burning period [28,32]. Figure 4 depicts the daily averaged maps of TCO and AAI over Oceania between 31 December 2019 and 2 January 2020. One can observe the strong correlation in the shape and spatiotemporal distributions between CO and absorbing aerosol plumes. Therefore, it can be deduced that the absorbing aerosols observed over the Oceania region can be broadly attributed to the Australian biomass burning event. The most intense transport event of Australian smoke plume occurred during the aforementioned period with large values of TCO (higher than 9 × 10 18 molecule·cm −2 ) and AAI (higher than 12) observed over the Oceania region. The CO and absorbing aerosol plumes were transported by the northwesterly flow over the Tasman Sea and reached New Zealand during New Year's Eve (Figure 4a,b). On 31 December 2019, the southern part of New Zealand was significantly impacted by the Australian biomass burning emissions with TCO and AAI values equal to 10 × 10 18 molecule·cm −2 and 14, respectively. The TCO and AAI observations indicate that the main part of the biomass burning plume moved northeast of New Zealand on 1 January 2020 (Figure 4c,d). The biomass burning plume moved further east and was located over the Pacific Ocean on 2 January 2020 (Figure 4e,f). These large values of AAI (12)(13)(14) are consistent with the AOD values greater than 5 obtained from TROPOMI observations over the east of New Zealand on 2 January 2020 [94].
tralian biomass burning emissions with TCO and AAI values equal to 10 × 10 18 mole-cule·cm −2 and 14, respectively. The TCO and AAI observations indicate that the main part of the biomass burning plume moved northeast of New Zealand on 1 January 2020 ( Figure  4c,d). The biomass burning plume moved further east and was located over the Pacific Ocean on 2 January 2020 (Figure 4e,f). These large values of AAI (12)(13)(14) are consistent with the AOD values greater than 5 obtained from TROPOMI observations over the east of New Zealand on 2 January 2020 [94]. The analysis of the daily averaged map of TCO from IASI observations (such as those depicted in Figure 4) allowed us to identify five transport events of Australian CO plume The analysis of the daily averaged map of TCO from IASI observations (such as those depicted in Figure 4) allowed us to identify five transport events of Australian CO plume passing over the Lauder site in New Zealand between 22 December 2019 and 31 January 2020 (grey shading area in Figure 5). The characteristics of these transport events are reported in Table 1. Figure 5a depicts the daily mean evolution of the TCO anomalies and total column of ozone obtained from FTIR spectrometer at Lauder during the aforementioned period. The daily anomalies induced by the transport of the Australian biomass burning plume are calculated as a relative difference by considering the standard deviation of the monthly background means as the reference values. It is worth noting that the TCO anomalies obtained from ground-based measurements are observed during the transport events of Australian CO plumes detected from IASI observations (Figure 4a). This correlation between two independent instruments corroborates the fact that the passage of the biomass burning plume induced an increase in TCO with the maximum observed on 31 December 2019 and 16 January 2020. Conversely, the passage of the plume does not Remote Sens. 2021, 13, 3092 13 of 26 induce a systematic perturbation on the total column of ozone. It can be observed that the daily values of the total column of ozone are mainly contained inside the standard deviation of the monthly background means except for the 2-11 January 2020 period (Figure 5a). The transport of the biomass burning plume during the 2-11 January 2020 period led to concomitant increase in the CO and ozone total column. burning plume are calculated as a relative difference by considering the standard deviation of the monthly background means as the reference values. It is worth noting that the TCO anomalies obtained from ground-based measurements are observed during the transport events of Australian CO plumes detected from IASI observations (Figure 4a). This correlation between two independent instruments corroborates the fact that the passage of the biomass burning plume induced an increase in TCO with the maximum observed on 31 December 2019 and 16 January 2020. Conversely, the passage of the plume does not induce a systematic perturbation on the total column of ozone. It can be observed that the daily values of the total column of ozone are mainly contained inside the standard deviation of the monthly background means except for the 2-11 January 2020 period (Figure 5a). The transport of the biomass burning plume during the 2-11 January 2020 period led to concomitant increase in the CO and ozone total column.  The tropospheric and stratospheric columns of ozone were calculated from the FTIR observations at Lauder. The tropopause height was retrieved from radiosonde observations by the use of the thermal definition which is based on the lapse rate criterion [95]. The tropopause height was found on average at 11.5 ± 0.5 km during the 2019-20 austral  The tropospheric and stratospheric columns of ozone were calculated from the FTIR observations at Lauder. The tropopause height was retrieved from radiosonde observations by the use of the thermal definition which is based on the lapse rate criterion [95]. The tropopause height was found on average at 11.5 ± 0.5 km during the 2019-20 austral summer months (December to February) at Lauder. We note that this value is in agreement with previous works [61,96,97]. Figure 6b depicts the evolution of the anomalies of the tropospheric and stratospheric columns of ozone at Lauder between 22 December 2019 and 31 January 2020. The anomalies were calculated by following the methodology mentioned previously. The concomitance between the tropospheric ozone anomalies and the Australian fire CO plume inferred from IASI observations are clearly shown. These transport events are, hence, marked by an increase in tropospheric ozone column with the maximum (44%) observed on 16 January 2020. Figure 5b also reveals that the tropospheric ozone anomalies are fairly correlated with the tropospheric CO anomalies. It can be observed that the stratospheric column of ozone is not significantly perturbed by the transport events. This could be explained by the temporal evolution of the tropopause height. A statistically significant increase in ozone (8%-15%) in the stratosphere is only observed during the 2-11 January 2020 period. The tropospheric ozone anomalies induced by the transport of the Australian biomass burning plume are, hence, higher than those observed in the stratospheric column. and 31 January 2020. The anomalies were calculated by following the methodology mentioned previously. The concomitance between the tropospheric ozone anomalies and the Australian fire CO plume inferred from IASI observations are clearly shown. These transport events are, hence, marked by an increase in tropospheric ozone column with the maximum (44%) observed on 16 January 2020. Figure 5b also reveals that the tropospheric ozone anomalies are fairly correlated with the tropospheric CO anomalies. It can be observed that the stratospheric column of ozone is not significantly perturbed by the transport events. This could be explained by the temporal evolution of the tropopause height. A statistically significant increase in ozone (8%-15%) in the stratosphere is only observed during the 2-11 January 2020 period. The tropospheric ozone anomalies induced by the transport of the Australian biomass burning plume are, hence, higher than those observed in the stratospheric column.    Figure 6a,b depict the available OMR profiles obtained from radiosonde observations performed at Lauder during the transport events of the Australian CO plume. On 30 December 2019, the ozone profile exhibits a significant increase (16%) in the OMR between 2 and 4 km altitude (Figure 6a). The daily ozone profile falls inside the standard deviation of the monthly background profile when a "no plume" is transported over Lauder, as illustrated on 22 December 2019 (Figure 6a). Figure 6b reveals a steep concentration increase (60%) in the local UT-LS (10.5-12.5 km) on 7 January 2020. This corroborates the concomitant increase in the tropospheric and stratospheric ozone columns obtained from FTIR observations during the 2-11 January 2020 period (Figure 5b). The radiosonde measurements also confirm that the ozone contained in the stratosphere was perturbed only by the transport of the biomass burning plume that occurred during this period. Previous works have shown that the pyro-convective activity has also contributed to inject biomass burning gases such as CO and water vapor into the stratosphere [35,36]. Based on the MLS observations, Khaykin et al. [36] showed that the stratospheric mass of CO bounded within the southern extra-tropics increased abruptly by about 20% compared to the value observed before the event. The OMR anomalies are also observed during the period following the pyro-convective outbreaks, as observed on 27 January 2020 (Figure 6b). On this day, the OMR profile is marked by a significant increase (23%) located between 5 and 7.5 km altitude. It is worth noting that the air mass back-trajectory calculations performed with the on-line HYSPLIT model (not shown) confirm the link between these ozone concentration anomalies and the Australian biomass burning activity. Given the high temporal resolution of the FTIR observations, daily analysis of the tropospheric ozone and CO mixing ratios is possible during the period of the transport events. This is illustrated with Figure 6c,d, which report the tropospheric ozone and CO mixing ratio profiles, respectively, during the 21-27 January 2020 period. Figure 6c confirms a significant increase in the OMR (25-30%) between 3 and 8 km altitude on 21 and 27 January 2020. This is fairly consistent with the ozone radiosonde observations performed on 27 January 2020 (Figure 6b). Figure 6d reveals that the OMR increase occurred in parallel with the CO mixing ratio increase within the same altitude range. The concomitant increase in ozone and CO in the troposphere suggests the production of ozone in the tropospheric smoke plume. Several studies have clearly shown that CO is an ozone precursor [98,99]. Nevertheless, it is worth noting that biomass burning plumes are chemically reactive, producing ozone and aerosols also through NO x -VOC chemistry [100,101]. Baylon et al. [100] showed that the degree of NO x oxidation is a key predictor of ozone production. On average, an enriched ozone layer is identifiable in the altitude range of 2-8 km, outside the pyro-convective outbreaks period (Table 1).
To further discuss the production of ozone in the tropospheric smoke plume, the ∆O 3 /∆CO enhancement ratio was estimated for different transport events. This ratio is defined as the excess OMR above background ozone in air mass normalized by an enhancement in CO mixing ratio above background value. The ∆O 3 /∆CO enhancement ratio is often used as a measure of ozone production efficiency in fire plume [102][103][104]. The ∆O 3 /∆CO enhancement ratio values obtained for the transport events are reported in Table 1. The values obtained for the transport events detected between 30 December 2019 and 12 January 2020 are less than 0.15 ppbv/ppbv. Previous works showed that the ∆O 3 /∆CO enhancement ratio is strongly influenced by the travel time between the fire and the measurement location [105,106]. Hence, it is possible to obtain an estimation of the plume age. The values obtained during the aforementioned transport events correspond to a plume age of less than 2 days for the tropical and subtropical regions [105,106]. Conversely, the ∆O 3 /∆CO enhancement ratio values obtained for the two transport events detected after 13 January 2020 are higher than 0.2 ppbv/ppbv. Table 1 reveals that the ratio values obtained during the transport events #3 and #4 are 0.41 ± 0.25 ppbv/ppbv and 0.35 ± 0.10 ppbv/ppbv, respectively. These values correspond to a plume age between 2 and 5 days for the tropical and subtropical regions [105,107]. The plumes sampled during these transport events are, hence, relatively freshly emitted (between 1 and 5 days) from fires located over Australia.
The transport of Australian biomass burning plume led to significant perturbation of the chemical composition in the troposphere over Lauder and its surroundings. It is clearly shown that ozone has been produced in the tropospheric smoke plume. This corroborates the results found by Baray et al. [108]. Based on radiosonde observations and back-trajectory modeling, they showed that the tropospheric ozone variability at Lauder is influenced by biomass burning plumes emitted from Australian fires during the austral summer months. In order to extend the discussion, the dispersion of the tropospheric ozone plume will be analyzed in the next section. Figure 7 shows maps of IASI TCO and tropospheric ozone column averaged over three successive periods: 18-20 December 2019, 4-5 January and 8-11 January 2020. These periods are associated with the phase of intensification of CO emission over southeastern and eastern Australia (Section 3). A fairly good agreement in the shape and the spatio-temporal distributions between the TCO and tropospheric ozone column is observed. Thus, it can be inferred that the tropospheric ozone variability was mainly driven by the Australian smoke plume. On 18-20 December 2019, the highest values of TCO (6.3 × 10 18 molecules·cm −2 ) and tropospheric ozone column (65)(66)(67)(68)(69)(70) are found over the eastern Australia coast (Figure 7a,b). Samoa is located at the edge of the CO plume (4-4.5 × 10 18 molecule·cm −2 ). The tropospheric ozone column values observed over Samoa range from 40 to 50 DU. Despite the fact that Amazonian biomass burning activity is usually sparse in December and January, moderate values of TCO (4-5 × 10 18 molecules·cm −2 ) and tropospheric ozone column (30)(31)(32)(33)(34)(35)(36)(37)(38)(39)(40) are observed over South America (Figure 7a,b). MODIS fire products (not shown) confirm the coincidence with these moderate CO and tropospheric ozone column values areas and fire pixels over South America, by showing a few biomass burning events less intense than those observed during the Amazonian biomass burning season. Part of the ozone and the CO produced by these fires are usually transported toward the eastern Pacific [43,108]. This explains the moderate values of CO and tropospheric ozone column observed over the northwestern coast of South America and its surroundings. Hence, we can exclude any association between the moderate values of CO and tropospheric ozone column over the aforementioned region and the Australian biomass burning activity at this period of time. The analysis of the IASI observations revealed that the transport of the Australian smoke plume occurring out of the pyro-cumulus outbreak period did not cross the Southern Pacific. Up until 29 December 2019, the smoke plumes were transported mainly over the Oceania region by following fairly pathways similar to those in Figure 7a,b. The smoke plumes were mainly transported to the north of New Zealand within subtropical and tropical latitudinal bands in the days before the pyro-cumulonimbus outbreaks period (Figure 7a,b).

Spatial Extend of the Tropospheric Ozone Plume
In the days following the pyro-cumulonimbus events, the smoke plumes were transported toward Southern Pacific, as is illustrated in Figures 4 and 7c. This could be explained by the fact that the pyro-convection events can quickly lift aerosol and chemical species to the UT-LS region, leading to changes in its composition [28]. Figure 7c reveals that the CO plume (from 6.3 × 10 18 to 9 × 10 18 molecules·cm −2 ) crossed the Southern Pacific by passing over New-Zealand during the 4-5 January 2020 period. The tropospheric ozone column values associated with the transport of the Australian biomass burning plume over the Southern Pacific ranged from 60 to 70 DU during the aforementioned period ( Figure 7d). On 4-5 January 2020, the CO and tropospheric ozone plumes were close to the tip of South America. This corroborates the observations reported by Ohneiser et al. [37] at Punta Arena (53.2 • S, 70.9 • E; Chile). By using LIDAR measurements at Punta Arena, they showed that the smoke layer was visible over Punta Arena since 6 January 2020 in the upper troposphere (between 6 and 11 km altitude) and has fallen gradually over time. The tropospheric plume was driven rapidly by the westerly winds and finally dissipated over the Atlantic Ocean. On 8-10 January 2020, the CO plume is centered over the Southern Pacific and is organized following an anti-cyclonic rolling extending between 20 • S and 40 • S (Figure 7e). The plume's structure in Figure 7e can be associated with the anticyclonic rolling of the CO plume. This is consistent with the work reported by Khaykin et al. [36]. They showed that on its way across the Southern Pacific during the first week of 2020, the core plume started to encapsulate into a compact bubble-like structure. This was identified in CALIOP observations on 7 January 2020 as an isolated 4 km vertically extended and 1000 km horizontally wide structure. The IASI observations reveal that the CO plume, by its rolling and stretching, has extended the tropical and subtropical latitudinal bands during the 8-14 January 2020 period. Figure 7f shows that the tropospheric ozone plume is also centered over the Southern Pacific between 20 • S and 40 • S with values ranging from 40 to 50 DU. The tropospheric ozone plume is not clearly organized following the anti-cyclonic rolling structure in Figure 7f. The structure of the tropospheric ozone plume is more diffuse and extends over the Samoan measurement site (Figure 7f). Furthermore, it was observed that plumes of CO and tropospheric ozone have spread further into the extra-tropical latitudes after 14 January 2020, crossing the southwestern part of South America one day later. Previous works confirm that South America is frequently impacted by the Australian biomass burning plume during the pyro-convective activity [28,33]. Through the use of OMI RGB images and AAI observations, Dirksen et al. [28] showed that the 2006 Australian bushfire smoke plumes were transported eastward and traveled from Tasmania to Chile within five days.  In the days following the pyro-cumulonimbus events, the smoke plumes were transported toward Southern Pacific, as is illustrated in Figures 4 and 7c. This could be explained by the fact that the pyro-convection events can quickly lift aerosol and chemical species to the UT-LS region, leading to changes in its composition [28]. Figure 7c reveals that the CO plume (from 6.3 × 10 18 to 9 × 10 18 molecules·cm −2 ) crossed the Southern Pacific by passing over New-Zealand during the 4-5 January 2020 period. The tropospheric ozone column values associated with the transport of the Australian biomass burning In addition to local tropospheric ozone sources from northern parts of South America, the Australian biomass burning plumes contributed, hence, to tropospheric ozone column variability over the southern part of South America, as is illustrated in Figure 7.

Statistical Variation of the Tropospheric Ozone
Given the presence of local tropospheric ozone sources, it is difficult to estimate properly the contribution of the Australian biomass burning activity to the tropospheric ozone variability over Southern America from space-borne observations. As a consequence, the discussion focuses on tropospheric ozone variability over Oceania. It can be observed that the transport of the Australian biomass burning plume over Oceania occurred mainly within the 18 • S-60 • S latitudinal band (Figure 7). To further investigate the latitudinal ozone variability induced by the transport of the Australian biomass burning plume over the Oceania region, the daily mean evolution of the OMR obtained from MLS observations between 18 December 2019 and 31 January 2020 was calculated over a global domain extending between 10 • S and 60 • S in latitude and 160 • E and 160 • W in longitude (the black box in Figure 7b) at 215 (~12 km) and 100 hPa (~16 km) levels ( Figure 8). The ozone anomalies at 215 hPa are mainly observed during the two first weeks of January for the three latitudinal bands: tropical (10 • S-19 • S), subtropical (20 • S-29 • S) and extratropical (30 • S-60 • S). Over the extratropical latitudinal band, the maxima of OMR anomalies (from 12% to 47%) were obtained during the 2-8 January 2020 period (Figure 8e). This suggests that ozone mixed vertically by convection over southeastern and eastern Australia and was then transported toward this latitude band in the UT-LS. Conversely, the maxima of the OMR anomalies (from 5% to19%) is observed during the 6-12 January 2020 period over the tropical latitudinal band (Figure 8c). This could be explained by the presence of the large and stretchable ozone plume over the Southern Pacific during the aforementioned period. Given its geographical extension, we cannot exclude that this large ozone plume led to sporadic ozone anomalies within the tropical and subtropical latitudinal bands over Oceania. It is worth noting that the daily evolution of the OMR does not exhibit statistically significant perturbation at 100 hPa over the tropical latitudinal band (Figure 8b). Conversely, the daily evolution of the OMR exhibits statistically significant perturbation over extratropical and subtropical latitudinal bands, with the maximum observed over the latter (Figure 8d,f). The IASI observations (Figures 4 and 7) reveal that parts of the smoke plume are contained within the subtropical latitudinal bands during the pyrocumulonimbus outbreaks. The large and stretchable ozone plume over the Southern Pacific had the potential to induce ozone anomalies within the Oceanian subtropical band. In OMR evolution at 215 and 100 hPa over the subtropical latitudinal band, there is a statistically significant increase in ozone (from 19% to 46%) during the 1-14 January 2020 period (Figure 8c,d). Figure 9a,b depict the available OMR profiles obtained from radiosonde observations recorded at Samoa between 18 December 2019 and 15 January 2020. The ozone profile exhibits a significant increase (43%) in the OMR between 6 and 9 km altitude, on 18 December 2019 (Figure 9a). Figure 7a,b corroborate that this enhanced ozone layer is linked to the transport of Australian biomass burning plume that occurred between 18 and 20 December 2019. Given the reasons previously mentioned and collaborating satellite observations, we can exclude all contributions from Southern African and Indonesian sources on the enhanced ozone seen in the mid-troposphere on 18 December 2019. The enhancement lasted for only a few days. On 26 December 2019, the daily ozone profile falls inside the standard deviation of the monthly background profile (Figure 9a). It can be observed that the OMR profile recorded at Samoa exhibits two distinct enhanced ozone layers centered in the mid (6-9 km) and upper troposphere (10-12 km), on 8 January 2020 (Figure 9b). These two enhanced ozone layers do not appear to be connected to each other. The origin of the enhanced ozone layer into the mid-troposphere is similar to the one observed on 18 December 2019. Figure 7f shows that the transport of the tropospheric ozone plume towards the Tasman Sea was located over Northeastern Australia during the 8-10 January 2020 period. This biomass burning plume has contributed to a significant increase (from 10% to 43%) in OMR over Samoa (Figure 7f). We cannot exclude that the presence of this large and stretchable ozone plume has led to a significant increase (29%) in OMR at Samoa between 10 and 12 km. We suspect that this ozone profile recorded in Samoa resulted from air-mass crossing in the local mid-troposphere, from eastern Australia and the Western Pacific. On 15 January 2020, it can be observed that the daily ozone profile falls inside the standard deviation of the monthly background profile climatology. On 15 January 2020, it can be observed that the daily ozone profile falls inside the standard deviation of the monthly background profile. This could be explained by the fact that the biomass burning plume was preferentially transported toward extratropical latitudes since 14 January 2020. As a consequence, the tropical sites such as Samoa are less impacted by the biomass burning plume. It is worth noting that a similar increase in the OMR was obtained at Lauder and Samoa. For both sites, the enhanced ozone layers due to the transport of the Australian biomass burning are in the mid-and upper-troposphere. Overall, we conclude that the transport of the smoke plume into the troposphere was not bounded exclusively within the extratropical latitudinal band. within the tropical and subtropical latitudinal bands over Oceania. It is worth noting that the daily evolution of the OMR does not exhibit statistically significant perturbation at 100 hPa over the tropical latitudinal band (Figure 8b). Conversely, the daily evolution of the OMR exhibits statistically significant perturbation over extratropical and subtropical latitudinal bands, with the maximum observed over the latter (Figure 8d,f). The IASI observations (Figures 4 and 7) reveal that parts of the smoke plume are contained within the subtropical latitudinal bands during the pyro-cumulonimbus outbreaks. The large and stretchable ozone plume over the Southern Pacific had the potential to induce ozone anomalies within the Oceanian subtropical band. In OMR evolution at 215 and 100 hPa over the subtropical latitudinal band, there is a statistically significant increase in ozone (from 19% to 46%) during the 1-14 January 2020 period (Figure 8c,d).  viation of the monthly background profile. This could be explained by the fact that the biomass burning plume was preferentially transported toward extratropical latitudes since 14 January 2020. As a consequence, the tropical sites such as Samoa are less impacted by the biomass burning plume. It is worth noting that a similar increase in the OMR was obtained at Lauder and Samoa. For both sites, the enhanced ozone layers due to the transport of the Australian biomass burning are in the mid-and upper-troposphere. Overall, we conclude that the transport of the smoke plume into the troposphere was not bounded exclusively within the extratropical latitudinal band.

Conclusions
We presented the tropospheric ozone and CO variability over Oceania and the Southern Pacific during the 2019-20 Australian bushfires. The meteorological context and the extensive fires over southeastern Australia were favorable to develop pyro-cumulonimbus events between 29 December and 12 January 2020. The latitudinal variation of the TCO over Australia obtained from IASI was found to be consistent with the latitudinal distribution of the AOD. The sun-photometer and CALIOP observations showed that the highest values of AOD were found between 15 December 2019 and 12 January 2020 over southern and central Australia. The spatial and temporal distributions of aerosol and CO over Australia are mainly explained by the recirculation of air masses from the southeastern to central region of the continent.
The analysis of the satellite observations revealed that transport events of Australian biomass burning plume took place during and outside of the pyro-cumulonimbus outbreak period. The dispersion of the smoke plume during these transport events extends from tropical to extratropical latitudes. The dispersion of these smoke plumes induced a significant variability on the tropospheric ozone column detected from satellite and FTIR observations. The ozone radiosonde measurements performed at Samoa and Lauder recorded enhanced ozone layers in the mid-tropospheric ozone linked with transport events of the Australian biomass burning plumes. This increase in ozone in the mid-troposphere was found to be consistent with MLS observations over the tropical and extratropical latitudinal bands. The smoke plume was transported over the Southern Pacific. The IASI observations confirm that smoke plume transport over the Southern Pacific had led to significant ozone variability in the troposphere. The smoke plume remained over the Southern Pacific for a few days before moving into the extra-tropics and reaching the southwestern part of South America. Over the southern Pacific, the smoke plume was organized following a stretchable anticyclonic rolling which impacted significantly on the ozone variability in the tropical and subtropical upper-troposphere over Oceania.
Our analysis suggests that these transport events of Australian biomass burning plumes significantly alter the chemical composition in the troposphere over Oceania and the Southern Pacific. Besides the injection of aerosol particles and chemical species into the stratosphere, the dynamical context has also contributed to modify the vertical distribution of aerosol and chemical species such as ozone in the troposphere at the regional scale. To understand the full contribution of the extreme Australian bushfire events to the tropo-spheric distribution of aerosol and chemical species, further investigation is required by including numerical models (such as WRF-Chem) with fire pollution sources. The dataset presented in this study would be useful to validate and constrain the numerical simulations. Given similar dynamical context was observed during previous extreme Australian bushfire events [28,31,109], it cannot be excluded that perturbations on tropospheric trace gas abundances were also observed during these events. The potential of these previous extreme bushfire events to induce significant perturbation on the tropospheric distribution of aerosol and chemical species over Australia and its surroundings will form the basis for a forthcoming statistical study.