key: cord-0788958-8qm3vd2l authors: Moore, G. W. K.; Semple, J. L. title: Himalaya Air Quality Impacts From the COVID‐19 Lockdown Across the Indo‐Gangetic Plain date: 2021-06-01 journal: Geohealth DOI: 10.1029/2020gh000351 sha: 6374a0f7d5cdd1f58aa622780f1710e76b24df8b doc_id: 788958 cord_uid: 8qm3vd2l Starting in January 2020, the novel coronavirus, now known as acute respiratory syndrome coronavirus (SARS‐CoV‐2) and the disease that it causes (COVID‐19) has had significant impacts on human health, the environment, and the economy globally. The rapid lockdown that occurred as well as its well documented timing allows for an unprecedented opportunity to examine the impact of air pollution from densely populated regions has on adjacent and pristine environments. Here, we use in situ and satellite observations to show that there was a step function decrease in two key indicators of air quality, nitrogen dioxide and airborne particulates, in locations within the Indo‐Gangetic Plan (IGP) as a result of the Spring 2020 lockdown. Based on anomaly patterns, we find a dipole response with a statistically significant reduction in air pollution along the western IGP and Himalaya and an increase in air pollution in the eastern IGP and Himalaya. We show that spatial variability in the reductions in economic activity across northern India and the adjoining countries of Nepal, Pakistan, and Bangladesh contributed to this dipole as did a persistent atmospheric circulation anomaly across the region during the lockdown. of air pollution (WHO, 2006) in and of itself as well as being a precursor species to ozone, a pollutant and oxidant that has been shown to impact lung function (Lippmann, 1989) . We will also consider airborne particulate matter with diameters smaller than 2.5 μm (pm2.5) that have also been shown to impact human health (WHO, 2006) . Based on anomaly patterns, we find a dipole response during the Spring 2020 lockdown with a statistically significant reduction in air pollution along the western IGP and Himalaya as well as an increase in air pollution in the eastern IGP and Himalaya. We argue that spatial variability in the reductions in economic activity across northern India and the adjoining countries of Nepal, Pakistan, and Bangladesh contributed to this dipole as did a persistent atmospheric circulation anomaly across the region during the lockdown. Data on daily mean surface concentrations of NO 2 during the period of interest, December 1, 2019-May 1 2020, were accessed from OpenAQ Platform for the four locations discussed in the paper: Delhi, Ludhiana, Singrauli, and Asansol. Delhi had ∼45 reporting sites that were averaged to obtain the time series for this location. All other locations had individual reporting sites. For the period of interest, data availability is greater than 98%. To develop a climatology for the surface concentration of NO 2 at Delhi, data from the OpenAQ Platform was accessed for the period 2017-2019 as well as data from the archives of the Central Pollution Control Board was accessed for the period 2005-2016. The multiple reporting sites were averaged to obtain the climatology for Delhi. The climatology did not include 2020 so as to not introduce a bias due to the observed low values during the lockdown. Data on daily mean surface pm2.5 concentrations were obtained from United States Embassies or Consulates in Delhi India, Lahore Pakistan, Dhaka Bangladesh and Kathmandu Nepal as archived by the U.S. Environmental Protection Agency. For the period of interest, data availability was 100%. Data at Delhi is available back to 2016 and the data for the period 2016-2019 was used to develop a climatology for this location. To characterize the impact of the lockdown on the spatial and temporal variability in air quality across the IGP and the Himalaya, we use daily tropospheric NO 2 column retrievals from the ozone monitoring instrument (OMI) instrument (Boersma et al., 2011) on National Aeronautics and Space Administration (NASA)'s Aura satellite as well as higher spatial resolution data from the TROPOspheric Monitoring (TROPOMI) instrument (Veefkind et al., 2012) Daily tropospheric NO 2 column retrievals from the OMI instrument on NASA's Aura satellite, available from 2005 to 2020 at a spatial resolution of ¼ o or ∼25 km (Boersma et al., 2011) , as well as the TROPOMI instrument on ESA's Sentinel 5P satellite, available from 2018 to 2020 at a spatial resolution of ∼4 km (Veefkind et al., 2012) were used to characterize the spatial distribution of pollutants. For each data set, the cloud cleared retrievals were used. The 16 years length of data availability for the OMI instrument allows for the development of a climatology. Information on the spatial distribution of particulate matter is also available from the OMI instrument through the so-called Aerosol Index (AI) that uses differential absorption of upwelling ultra-violet radiation to indicate the presence of aerosols associated with dust, smoke and volcanic ash (Torres et al., 2007) . It is available at a spatial resolution of 1° daily from 2005 to 2020 (Torres et al., 2007) . Mass-weighted wind data from the European Centre for Medium-Range Weather Forecasts' ERA5 Reanalysis (Hersbach et al., 2020) were used to characterize the atmospheric circulation during April 2020. The data is available at a spatial resolution of ∼30 km for the period 1979-2020. Climatological surface NO 2 data is available for Delhi back to 2005 and shows a seasonal decline in values during the transition from winter through the spring (Mohan & Kandya, 2007 ) that can also be seen in Figure 3a . A 31 days moving window smoother was applied to the climatological time series to reduce high frequency signals that are the result of their short length so as to allow a clearer picture of the seasonal cycle and its variability. This inter-annual variability in surface NO 2 concentration confirm that values recorded in Delhi during the lockdown were two standard deviations below the mean indicating highly anomalous conditions. Indeed, the monthly mean surface NO 2 concentration observed in Delhi during April 2020 was the lowest observed for the period 2005-2020. However, the marked reduction in surface NO 2 concentrations observed in Delhi and Ludhiana were not uniform across the IGP. Singrauli, a city with a number of large coal-fired power plants in its vicinity (Singh et al., 2018) , recorded no significant reduction in surface NO 2 concentrations, i.e., the means over the period February 15-March 15, 2020 and April 1-30, 3030 were both 42 µg/m 3 , suggesting that the emissions from these plants were not reduced during the lockdown (Figure 3c ). Asansol, a major industrial city with significant coal mining and steel production activities (Reddy & Ruj, 2003) , in West Bengal showed a reduction in surface NO 2 values, from 25 µg/m 3 over the period February 15-March 15, 2020 to 18 µg/m 3 over the period April 1-30, suggesting a more modest lockdown than what occurred in Delhi and Ludhiana (Figure 3d ). Time series of pm2.5 concentrations from U.S. embassies and consulates throughout the IGP (Figure 4) indicates the presence of step function reductions in both Delhi as well as Lahore Pakistan, confirming that a lockdown also occurred in the Pakistani region of the IGP. Dhaka also showed a reduction however no indication of a step function change. This is most likely a reflection of the strong seasonality in air quality in the region with a winter peak associated with kiln production (Begum et al., 2011) . There was also evidence of a step function response in Kathmandu although there was significant variability throughout April indicating either an incomplete shutdown of industrial activity or the long-range transport of pollution into the area (Bonasoni et al., 2012) . The climatological distribution of monthly mean tropospheric NO 2 column densities from OMI during April ( Figure 5a ) shows elevated levels across the IGP with lower levels over eastern India, southwestern India and Tibet. There is a gradient in density across the Himalaya that approximately follows the 3000 m height contour suggesting that topography plays a role in the distribution of air pollution across the region. More information on the spatial variability on tropospheric NO 2 column densities associated with the lockdown is provided by the monthly mean anomaly during April 2020 (Figure 6a) , i.e., the difference between the monthly mean during April 2020 and the April climatological monthly mean for the period 2005-2020. The widespread reduction across the northwestern IGP as well as adjoining foothills of the western Himalaya, i.e., to the west of Kathmandu, is evident. The region around Singrauli stands out as one in which there was a marked increase in column density during April 2020. Throughout much of the eastern IGP, including much of Bangladesh, eastern India, eastern Nepal and Bhutan, there was an increase in column density with the exception of region around Dhaka. There is of course variability from year to year in air quality that is a function of variability in sources as well as meteorological conditions (Fishman et al., 2005; Pawar et al., 2017) . One way to assess the magnitude of the anomalous tropospheric NO 2 column densities during April 2020 is to compare it to the expected variability from year to year during April. This so-called normalized anomaly is constructed by dividing the April 2020 anomaly by the standard deviation of the tropospheric NO 2 column density during April and expressing it as a percentage (Figure 6b ). Normalized anomaly values on the order of ±200% indicate that such an anomaly has a magnitude of 2 standard deviations and should, assuming normally distributed data, occur approximately 2%-3% of the time. Across much of the western IGP, including Delhi, Ludihana, and Lahore, as well as the western Himalaya, the normalized anomaly was on the order of 200% confirming the extreme reduction in air pollution in this region, a result in agreement with the surface observations of MOORE AND SEMPLE 10.1029/2020GH000351 6 of 11 surface NO 2 concentration ( Figure 3 ) and surface pm2.5 (Figure 4) across the western IGP. In contrast, the normalized anomaly across much of the eastern IGP and the eastern Himalaya was positive with values between 100% and 200% indicating that air pollution values in this region were anomalously high, although with significance levels lower than that in the west of the region of interest. This result is also in agreement with regional surface NO 2 concentrations during the lockdown (Figure 3 ) that show lower concentrations in Delhi and Ludhiana as well as higher concentrations in Singrauli and Asansol. We will refer to this feature, reduction in air pollution in the western IGP and Himalaya and an increase in the corresponding eastern regions, as an air pollution dipole anomaly. Figure 5 and 6 indicate that there is evidence of spatial gradients in tropospheric NO 2 column densities along the Himalaya that may be under resolved by the OMI instrument. A comparison between TROPOMI monthly means for April 2019 and 2020 (Figure 7a and 7b) show clear evidence of a reduction in tropospheric NO 2 column densities across much of the IGP, including in the vicinity of Delhi, Lahore, Asansol, and Dhaka, as well as pronounced spatial gradients along the Himalaya. The difference between the two Aprils ( Figure 7c ) highlights the aforementioned dipolar nature of the changes that occurred with reductions across much of the western IGP and the western Himalaya and increases in the eastern areas of these regions. 10.1029/2020GH000351 7 of 11 Figure 6 . Anomalous column tropospheric NO 2 density from the OMI across the region of interest during April 2020. Results are shown for the: (a) anomaly (μ mol/m 2 ) (b) the normalized anomaly (%). The locations of New Delhi ("+"), Ludhiana ("#"), Singrauli ("o"), Asansol ("x"), Lahore ('l), Kathmandu ("k") and Dhaka ("d") are indicated. Political boundaries are shown by the thin black curves with the 3000 m topographic height shown by the thick black curve. Observations of surface particulate matter (Figure 4 ) also indicated the presence of step function responses across the western IGP and Nepal that were associated with the lockdown. The climatological monthly mean AI for April (Figure 8a layers characterized by high values of anthropogenic aerosols (Bonasoni et al., 2012; Ramanathan et al., 2007) . The monthly mean AI for April 2020 (Figure 8b) shows reduced AI values across the IGP, consistent with the in situ observations of surface pm2.5 (Figure 3) . The normalized anomaly for April 2020 (Figure 8c) indicates that the reduction in AI across the western IGP and the western Himalaya was between 1 and 1.5 standard deviations below the mean. Unlike the case for the tropospheric NO2 column density, these is no evidence of a coherent dipolar pattern with higher values of the AI in the eastern regions of the IGP and the Himalaya. Atmospheric circulation patterns play a role in the long-range transport of pollutants into relatively pristine regions such as the Himalaya (Bonasoni et al., 2012) . As such is it important to characterize the degree to which the atmospheric conditions contributed to the spatial distribution of the air pollution anomalies associated with the COVID-19 lockdown across the IGP. Figure 9 shows the April 2020 monthly mean anomaly in mass weighted velocity field, representative of circulation in the lower troposphere, over the region of interest as represented in the ERA5 reanalysis (Hersbach et al., 2020) . It indicates the presence of a cyclonic (i.e., counter-clockwise) circulation anomaly over the region. A cyclonic circulation anomaly would have acted to mitigate the transport of pollutants from the IGP toward the western Himalaya and enhance this transport toward the eastern Himalaya. As such, it would have acted to support the existence of the dipolar pattern in the distribution of pollutants along the IGP and Himalaya during April 2020. High mountainous regions of the world, traditionally, have been considered some of the most pristine environments and the least likely setting for hazards related to air pollution. However unique communities in high-altitude regions, such as those in the Himalaya, are unexpectedly exposed to pollution levels that are similar to, if not higher than, those reported in industrialized cities and other polluted environments (Saikawa et al., 2019) . However little is known about the sources for this pollution and its pathways (Bonasoni et al., 2012) . In addition, hypoxia and other unique characteristics of this high altitude ecosystem have compounding consequences for those who reside in and work in this region (Korrick et al., 1998; Pandey, 1984; Semple et al., 2016) . The COVID-19 lockdown within the heavily polluted IGP provides a unique opportunity to assess the impact of pollution from this region on the Himalaya. The lockdown produced a dipolar signature in pollution levels across the IGP with reductions in the west, that were two standard deviations below the climatological mean during April 2020 and a more modest increase in the east. This was most likely the result of differences in the local response to the lockdown with large cities such as Delhi and Lahore experiencing a reduction in activities such as transportation and construction that are sources of urban air pollution (Dhammapala, 2019; Mohan & Kandya, 2007) ; while industrial emitters such as those associated with steel production and power generation in the east continued (Reddy & Ruj, 2003; Singh et al., 2018) . Changes in biomass burning during 2020 across the IGP may have also contributed to the signal (Gutti-MOORE AND SEMPLE 10.1029/2020GH000351 9 of 11 kunda & Nishadh, 2020). A similar pattern was seen along the Himalaya with decreases in the west and increases in the east. This pattern is also consistent with anecdotal observations indicating that, for the first time in recent history, the visibility of Himalayan peaks from India's National Capital region (Chamas, 2020; Ellis-Petersen et al., 2020) . It is also consistent with other reports of the impact of the COVID-19 lockdown on regional air quality (Biswal et al., 2020; Kumar et al., 2020; Sharma et al., 2020) . A persistent atmospheric anomaly characterized by cyclonic (counter-clockwise) air flow throughout the lower troposphere during April 2020 contributed to the dipolar pattern seen in both the IGP and Himalaya. These results confirm the importance of characterizing the spatial distribution of emissions across the IGP as well as the background atmospheric circulation when assessing the impact of pollution on the Himalaya. They also provide an important test case that can be used to assess the ability of chemical transport models to capture regional patterns in air quality in the presence of large-scale changes in emissions (Moorthy et al., 2013; Sharma et al., 2017) . Atmospheric changes from solar eclipses Key issues in controlling air pollutants in COVID-19 lockdown induced changes in NO 2 levels across India observed by multi-satellite and surface observations An improved tropospheric NO 2 column retrieval algorithm for the Ozone Monitoring Instrument Atmospheric pollution in the Hindu Kush-Himalaya Region Mount Pinatubo aerosols, chlorofluorocarbons, and ozone depletion Peaks of Himalayas visible from parts of India for first time in decades as pollution drops amid lockdown. Australian Broadcasting Corporation The global atmospheric environment for the next generation Analysis of fine particle pollution data measured at 29 US diplomatic posts worldwide It's positively alpine!': Disbelief in big cities as air pollution falls. The Guardian Interannual variability of stratospheric and tropospheric ozone determined from satellite measurements Temporal variation of urban NOx concentration in India during the past decade as observed from space Air quality trends and lessons learned in India during the COVID lockdowns The ERA5 global reanalysis Do contrails significantly reduce daily temperature range? Geophysical Research Letters Effects of ozone and other pollutants on the pulmonary function of adult hikers Temporary reduction in fine particulate matter due to 'anthropogenic emissions switch-off' during COVID-19 lockdown in Indian cities Fingerprinting the impacts of aerosols on long-term trends of the Indian summer monsoon regional rainfall Health effects of ozone: A critical review An analysis of the annual and seasonal trends of air quality index of Delhi High concentration of surface ozone observed along the Khumbu Valley Nepal Performance evaluation of chemistry transport models over India Prevalence of chronic bronchitis in a rural community of the Hill Region of Nepal Inter annual variability of tropospheric NO2 and tropospheric ozone over Maharashtra (India): The role of lightning Warming trends in Asia amplified by brown cloud solar absorption Ambient air quality status in Raniganj-Asansol area Air Pollution in the Hindu Kush Himalaya High concentrations of ozone air pollution on Mount Everest: Health implications for Sherpa communities and mountaineers Photo-chemical transport modeling of tropospheric ozone: A review Effect of restricted emissions during COVID-19 on air quality in India. The Science of the Total Environment Elevated black carbon concentrations and atmospheric pollution around Singrauli coal-fired thermal power plants (India) using ground and satellite data Aerosols and surface UV products from Ozone Monitoring Instrument observations: An overview Regional variations in U.S. diurnal temperature range for the 11-14 TROPOMI on the ESA Sentinel-5 Precursor: A GMES mission for global observations of the atmospheric composition for climate, air quality and ozone layer applications. Remote Sensing of Environment Air quality guidelines for particulate matter, ozone, nitrogen dioxide and sulfur dioxide: Global update 2005: Summary of risk assessment. World Health Organization The authors declare no conflicts of interest relevant to this study. All satellite data was accessed through the NASA Goddard Earth Sciences (GES) Data and Information Services Center (disc.gsfc.nasa.gov). Surface NO 2 data is available from the OpenAQ Platform (openaq.org) and the Indian Open Government Data Platform (data.in.gov). Data on pm2.5 concentrations is available through the U.S. Environmental Protection Agency (airnow.gov/international/us-embassies-and-consulates). ERA5 data is available through the Copernicus Climate Data Store (cds.climate.copernicus.eu). Acknowledgments G.W.K. Moore was supported by the Natural Sciences and Engineering Research Council of Canada.