1.Key Laboratory of Middle Atmosphere and Global Environment Observation, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China Manuscript received: 2018-06-01 Manuscript revised: 2018-10-25 Manuscript accepted: 2018-11-20 Abstract:Boreal summer intraseasonal oscillation (BSISO) of lower tropospheric ozone is observed in the Indian summer monsoon (ISM) region on the basis of ERA-Interim reanalysis data and ozonesonde data from the World Ozone and Ultraviolet Radiation Data Centre. The 30-60-day intraseasonal variation of lower-tropospheric ozone shows a northwest-southeast pattern with northeastward propagation in the ISM region. The most significant ozone variations are observed in the Maritime Continent and western North Pacific. In the tropics, ozone anomalies extend from the surface to 300 hPa; however, in extratropical areas, it is mainly observed under 500 hPa. Precipitation caused by BSISO plays a dominant role in modulating the BSISO of lower-tropospheric ozone in the tropics, causing negative/positive ozone anomalies in phases 1-3/5-6. As the BSISO propagates northeastward to the western North Pacific, horizontal transport becomes relatively more important, increasing/reducing tropospheric ozone via anticyclonic/cyclonic anomalies over the western North Pacific in phases 3-4/7-8. As two extreme conditions of the ISM, most of its active/break events occur in BSISO phases 4-7/1-8 when suppressed/enhanced convection appears over the equatorial eastern Indian Ocean and enhanced/suppressed convection appears over India, the Bay of Bengal, and the South China Sea. As a result, the BSISO of tropospheric ozone shows significant positive/negative anomalies over the Maritime Continent, as well as negative/positive anomalies over India, the Bay of Bengal, and the South China Sea in active/break spells of the ISM. This BSISO of tropospheric ozone is more remarkable in break spells than in active spells of the ISM, due to the stronger amplitude of BSISO in the former. Keywords: boreal summer intraseasonal oscillation, tropospheric ozone, Indian summer monsoon, active/break spell 摘要:本文利用ERA再分析资料和WOUDC臭氧探空资料, 分析了印度季风区夏季对流层臭氧的季节内振荡(BSISO)特征. 结果表明对流层低层臭氧存在着30-60天的季节内振荡, 臭氧异常呈西北-东南向分布, 在季风区向东北方向传播. 最显著的臭氧异常在海洋大陆和西北太平洋区域. 臭氧异常在热带从地面延伸到300hPa, 而在热带外地区臭氧异常主要在500hPa以下. 在热带, BSISO引起的降水是导致臭氧异常的主要原因, 使得在1-3/5-6位相出现了臭氧负/正异常. 然而当BSISO传播到西北太平洋区域时, BSISO引起的大气环流异常对臭氧异常的形成更为重要, 通过反气旋性/气旋性环流异常导致了臭氧在3-4/7-8出现正/负异常. 大多数印度季风的活跃/中断发生在BSISO的4-7/1-8位相, 此时减弱/加强的对流出现在赤道东印度洋, 加强/减弱的对流出现在印度, 孟加拉湾和中国南海. 海洋大陆对流层臭氧出现显著的负/正异常, 印度, 孟加拉湾和中国南海对流层臭氧出现显著的正/负异常. 由于季风中断期BSISO的振幅比季风活跃期更大, 对流层臭氧的季节内振荡也更加显著. 关键词:北半球夏季季节内振荡, 对流层臭氧, 印度季风, 活跃/中断期
HTML
--> --> --> -->
2.1. Data
The daily mean ozone, zonal wind and meridional wind are from the ERA-Interim dataset (e.g., Dee et al., 2011). The climatology of the lower-tropospheric ozone column is from Ozone Monitoring Instrument (OMI) data. OMI ozone profiles in the troposphere have been validated (Liu et al., 2010b) and used in studying lower-tropospheric ozone (Lu et al., 2018). The precipitation related to BSISO is based on TRMM observations (Kummerow et al., 2000). The daily mean outgoing longwave radiation (OLR) data are from the AVHRR instrument onboard NOAA's polar-orbiting spacecraft (Liebmann and Smith, 1996). The time period of these datasets considered in this study is summer (June-August) 2000-2012. It has been proven that BSISO has quasi-oscillating periods of 30-60 days (Lee et al., 2013). As a result, a 30-60-day bandpass filter is applied to the daily anomalies (to remove the climatology from the daily mean) of variables like OLR, horizontal wind and ozone to derive their BSISO characteristics (i.e., BSISO-related OLR, horizontal wind and ozone anomalies). However, as bandpass filtering will reduce the effective sample size, the regular Student's t-test is no longer capable of testing the significance of anomalies related to BSISO (e.g., Tian et al., 2011; Liu et al., 2014). Instead, we use a two-tailed Student's t-test with reduced degrees of freedom (Tian et al., 2011; Liu et al., 2014; Zhang et al., 2015) to examine the BSISO-related composite results in this work. We use ozonesonde profiles at three stations (Fig. 1) from the World Ozone and Ultraviolet Radiation Data Centre (WOUDC). Hanoi (AAR) station (21.01°N, 105.8°E) is located in northern Vietnam, nearly 90 km away from the coast. Ozonesonde observations have been made every two or three weeks since 2004, with a gap in 2011 and 2012. Sepang Airport (SEP) station (2.73°N, 101.7°E) in Singapore has a long history of ozonesonde observations, every week or every two weeks, since 1998. Naha (NAH) station (26.21°N, 127.69°E) is located in southernmost Japan, observing ozone profiles since 1989. Figure1. (a) Climatology of the lower-tropospheric ozone column (700 hPa to the surface) in June-July-August, based on OMI observations. The red dots are the three ozonesonde locations. The solid blue line is the northeastward transect, and the blue dashed line is the northward transect. (b) Variance of the BSISO-related lower-tropospheric ozone column (700 hPa to the surface) in June-July-August, based on ERA-Interim data.
2 2.2. BSISO index -->
2.2. BSISO index
BSISO is more complicated than wintertime MJO, as it involves the northward propagation of deep convection extending further from the equator to the extratropical Northern Hemisphere. The widely-used real-time multivariate MJO index (Wheeler and Hendon, 2004) is unable to show this northward propagation. Therefore, we use another BSISO index, suggested by (Lee et al., 2013), which is based on multivariate empirical orthogonal function (MV-EOF) analysis of OLR and zonal wind at 850 hPa in the region of (10°S-40°N, 40°-160°E). The BSISO index is defined by the first two principal components (PCs) of the MV-EOF analysis, representing the northward and eastward propagation of the 30-60-day BSISO. The lifecycle of BSISO is divided into eight phases according to this BSISO index, indicating the location of deep convection along its northeast propagation pathway. We select BSISO events with amplitudes greater than 1.0 [(PC12 + PC22)1/2>1.0] in this study.
2 2.3. Active/break events of the ISM -->
2.3. Active/break events of the ISM
There are intraseasonal variabilities of convection and rainfall over the ISM region. ISM rainfall fluctuates between being copious and scant, referred to as active and break spells of the ISM, respectively. These two spell types are usually defined during the peak monsoon months (July and August) by the average daily rainfall over the core ISM area (18.0°-28.0°N, 65.0°-88.0°E). According to the definition, an active (break) spell of the ISM is a period of three consecutive days or more during which the rainfall anomaly over the core ISM region is more (less) than +1.0 (-1.0) the standard deviation (Rajeevan et al., 2010).
-->
3.1. Tropospheric ozone in the ISM region
Figure 1a shows the climatology of the lower-tropospheric ozone column (700 hPa to the surface) in summer (June-July-August) based on OMI data. The ozone in the midlatitudes is higher than in the tropics, with the most significant gradients at 25°-30°N. In the tropics, the ozone over the western Pacific is lower than over the Indian Ocean. After applying a 30-60-day bandpass filter to the lower-tropospheric ozone column, we obtain the variance of BSISO-related ozone in the lower troposphere (Fig. 1b). There are two areas in which the most significant variance is observed: one is the Maritime Continent, especially Malaysia and the Celebes Sea; and the other is the western North Pacific. The lower-tropospheric ozone in these two areas has the most prominent variability on the intraseasonal time scale. The power spectra of the lower-tropospheric ozone (Fig. 2) in the Maritime Continent (5°S-10°N, 100°-140°E) and western North Pacific (15°-30°N, 120°-140°E) confirm that the ozone variance is concentrated in intraseasonal periods of 30-60 days, especially in the Maritime Continent (Fig. 2a). Figure2. Power spectra of BSISO-related anomalies of the lower-tropospheric ozone column (700 hPa to the surface) in the (a) Maritime Continent and (b) western North Pacific. The red curve is the red-noise spectrum. The lower and upper blue dashed curves are 5% and 95% red-noise significance levels respectively.
2 3.2. Horizontal distribution and propagation of BSISO-related ozone -->
3.2. Horizontal distribution and propagation of BSISO-related ozone
The lifecycle of BSISO is 30-60 days, as also reported by (Lee et al., 2013). They divided the BSISO's lifecycle into eight phases according to the location of deep convection. The composite 30-60-day bandpass-filtered (BSISO-related) OLR anomalies represent the location of deep convection in each BSISO phase (Figs. 3a-h). The enhanced/suppressed convection (negative/positive OLR anomalies), accompanied by a negative/positive lower-tropospheric ozone column (700 hPa to the surface) anomalies, show a northwest-southeast pattern, which propagates northeastward, during phases 1-8. The BSISO-related deep convection (negative OLR anomalies) generate over the equatorial eastern Indian Ocean in BSISO phase 1 and enhance in phase 2. Simultaneously, negative BSISO-related anomalies of ozone are observed over the eastern Indian Ocean and Maritime Continent. During phases 3-4, the BSISO-related convection weakens slightly and gradually propagates northeastward to India, the Bay of Bengal, and the Maritime Continent. The negative BSISO-related ozone anomalies over the eastern Indian Ocean and Maritime Continent are also reduced in this period. After phase 5, the negative BSISO-related OLR anomalies and the negative ozone anomalies gradually propagate to the western North Pacific, until the start of the next BSISO lifecycle. Similarly, suppressed convection (positive BSISO-related OLR anomalies) develops over the equatorial Indian Ocean in phase 4, and enhances during phases 5-7 with a northeastward propagation. The positive anomalies of BSISO-related ozone are observed over the eastern Indian Ocean and Maritime Continent during this period. As the suppressed convection propagates to the western North Pacific after phase 8, positive ozone anomalies also gradually propagate to the Northwest Pacific until the next phase 5, when positive ozone anomalies develop over the eastern Indian Ocean and Maritime Continent. Figure3. Composites of BSISO-related (a-h) OLR anomalies (units: W m-2) and (i-p) lower-tropospheric column (700 hPa to the surface) ozone anomalies (units: DU). Positive and negative anomalies are indicated by red solid and blue dashed lines, respectively. The shaded areas are statistically significant at the 95% confidence level, based on a two-tailed Student's t-test with reduced degrees of freedom. The top-right number in each panel indicates the number of days used for each composite.
However, there are slight inconsistencies between the locations of BSISO-related OLR anomalies and ozone anomalies. For example, when negative/positive OLR anomalies are in the eastern Indian Ocean in phases 1-2/5-6, the most significant negative/positive ozone anomalies are not, instead being situated in the Maritime Continent. These inconsistences indicate that there are more complicated factors influencing the BSISO of lower-tropospheric ozone. Thus, we investigate the factors that might influence the production and transport of lower-tropospheric ozone. Limited by the availability of data, we only find two factors that are closely related to the BSISO of lower-tropospheric ozone: precipitation anomalies and horizontal circulation anomalies. Precipitation will clean ozone precursors like volatile organic compounds in the troposphere and reduce temperature, resulting in a reduction of generated ozone. Horizontal circulation anomalies can transport ozone across the horizontal gradients of ozone, especially in the midlatitudes, where the most significant ozone gradients are observed (Fig. 1a). Furthermore, the composite BSISO-related precipitation anomalies, stream function anomalies, and horizontal wind anomalies (Fig. 4) show the important roles played by the precipitation and horizontal transport caused by BSISO in modulating the BSISO of lower-tropospheric ozone. Figure4. As in Fig. 3 but for (a-h) TRMM precipitation anomalies (units: mm) and (i-p) stream function anomalies (contours; units: 106 m2 s-1) and horizontal wind anomalies (vectors; units: m s-1) at 850 hPa. The red vectors are statistically significant at the 95% confidence level, based on a two-tailed Student's t-test with reduced degrees of freedom.
Compared to OLR, BSISO-related positive/negative precipitation anomalies are more consistent with negative/ positive ozone anomalies. In phases 1-3, positive precipitation anomalies extend form the eastern Indian Ocean to the equatorial Western Pacific, resulting in negative ozone anomalies over the Maritime Continent (Figs. 3i and k). At the same time, there is a cyclonic anomaly over the Maritime Continent and an anticyclonic anomaly over the north of the Maritime Continent (Figs. 4i-k). Strong easterly flow between them brings air containing low ozone (Fig. 1a) from the Pacific to the Maritime Continent, reducing the level of tropospheric ozone in the latter. Positive precipitation and anticyclonic anomalies gradually propagate northeastward to the western Pacific after phase 4. In the following phases 5-7, negative precipitation anomalies occur in the Maritime Continent (Figs. 4e-g), causing positive ozone anomalies (Figs. 3e-g). In this period, there is a cyclonic anomaly (Figs. 4m-o) north of the Maritime Continent, which is more powerful than the anticyclonic anomaly over the Maritime Continent. The northwesterly winds at the edge of the cyclone bring air containing more ozone (Fig. 1a) from mainland Southeast Asia to the Maritime Continent, increasing the tropospheric ozone in this area until phase 8, when the cyclonic anomaly propagates northeastward to the western North Pacific. It is worth noting that when positive precipitation anomalies propagate to the western North Pacific in phase 6 (Fig. 4f), the negative ozone anomalies become too weak to be observed (Fig. 3n). However, in the following phases 7 and 8, negative ozone anomalies once again occur in the western North Pacific (Figs. 3o and p), which seems inconsistent with disappearing positive precipitation anomalies (Figs. 4g and h). This is because of the strengthened cyclonic anomalies over the western North Pacific (Figs. 4o and p). The southerly winds at the eastern edge of the cyclone bring air containing less ozone in the tropics to the north, reducing the tropospheric ozone in the western North Pacific. Similarly, the positive ozone anomalies in the western North Pacific in phases 3 and 4 (Figs. 3k and l) are caused by northerly winds at the edge of anticyclonic anomalies (Figs. 4k and l). In order to elucidate the relative importance of precipitation and horizontal transport, we calculate the correlation between the BSISO-related ozone anomalies and the precipitation anomalies (Fig. 5a) as well as the stream function anomalies at 850 hPa (Fig. 5b). BSISO-related ozone is better correlated with precipitation anomalies in the eastern Indian Ocean and Maritime Continent; however, in the western North Pacific, ozone is better correlated with stream function anomalies. Based on the above analysis, one can conclude that precipitation plays a dominant role in modulating the BSISO of lower-tropospheric ozone in the tropics (eastern Indian Ocean and Maritime Continent). However, horizontal transport is relatively more important for the BSISO of lower-tropospheric ozone in the extratropics (western North Pacific). Figure5. Correlation coefficients of BSISO-related anomalies of lower-tropospheric column (700 hPa to the surface) ozone with (a) precipitation and (b) 850-hPa stream function (only correlation coefficients greater than 0.5 are plotted).
2 3.3. Vertical structure of BSISO-related ozone -->
3.3. Vertical structure of BSISO-related ozone
To investigate the vertical structure of the BSISO-related ozone anomaly, we choose two transect locations: northward, almost across AAR station; and northeastward, across SEP station and close to NAH station (Fig. 1). Vertical cross sections of BSISO-related ozone anomalies and stream function anomalies along these two transect locations are compared in Fig. 6. In both cross sections, BSISO-related ozone anomalies are most remarkable in the tropics, where they extend from the surface to almost 300 hPa. Besides, the ozone anomalies in the tropics have a northward tilt as the height increases. When the ozone anomalies propagate to extratropical regions, they quickly weaken with a limited vertical range from the surface to about 700 hPa. As the BSISO-related anomalies continue to propagate northeastward, ozone anomalies are observed over the western North Pacific along the northeastward transect and South Asia along the northward transect. Along the northeastward transect, the nearest locations from SEP and NAH stations are marked as vertical red lines (Fig. 6). SEP station experiences the most significant negative/positive ozone anomalies in phases 1-2/5-6; however, in the subtropics, NAH station shows the most significant negative/positive ozone anomalies in phases 8-1/4-5. Along the northward transect, locations nearest to SEP and AAR stations are marked as vertical red lines (Fig. 6). AAR station experiences the most significant negative/positive ozone anomalies in phases 2/6. Therefore, we test ozone profiles in these BSISO phases based on ozonesonde data from SEP, NAH and AAR stations (Fig. 7). At SEP station, the ozonesonde data show larger ozone partial pressure in phases 5-6 than in phases 1-2 at almost every height from the surface to 500 hPa. This is consistent with the most significant positive ozone anomalies in phases 5-6 (Figs. 6e and f) and negative ones in phases 1-2 (Figs. 6a and b), which extend from the surface to 500 hPa in the tropics. When BSISO propagates northeastward to the western North Pacific, the most significant positive/negative ozone anomalies are in phases 4-5/8-1 (Figs. 6d, e, h and a). This is proved by the ozonesonde data from NAH station, which suggest that the lower-tropospheric ozone partial pressure is larger in phases 4-5 than in phases 8-1. At AAR station, the ozone partial pressure in phase 2 is larger than in phase 5 in the lower troposphere, which is also consistent with the positive ozone anomaly in phase 2 (Fig. 6j) and negative one in phase 6 (Fig. 6n). Figure6. Vertical cross sections of BSISO-related ozone anomalies (units: DU km-1, colored contours) and stream function anomalies (units: 106 m2 s-1, black contours) at two transect locations: (a-h) northeastward transect location (blue solid line in Fig. 1); (i-p) northward transect location (blue dashed line in Fig. 1). The color-shaded areas are statistically significant at the 95% confidence level, based on a two-tailed Student's t-test with reduced degrees of freedom. The top-left number in each panel indicates the number of days used for each composite. The red vertical lines in (a-h) mark the locations nearest to SEP and NAH stations, while those in (i-p) mark the locations nearest to SEP and AAR stations. The values in parentheses denote degrees of latitude (positive for °N, negative for °S) and longitude (°E), respectively.
Figure7. Average ozone profiles (ozone partial pressure) in (a) BSISO phases 1-2 and 5-6 at SEP station, (b) BSISO phases 2 and 5 at AAR station, and (c) BSISO phases 8-1 and 4-5 at NAH station, from WOUDC.