Air Quality Analysis for Panjim-Margao-Vasco, India

APnA India APnA India APnA India APnA India

Panjim or Panaji is the capital of the state of Goa and lies on the banks of the river Mandovi. Panjim, Margao, Ponda and Vasco are the main cities in the state. Goa is synonymous with its beaches, the carnival and the International Film Festival that is held every year. Goa is a relatively small state with a lot of movement between towns/villages. A person may live in a village, and travel to Panjim or Margao for work or college and so it is hard to isolate the city of Panjim from the state.  

Apart from tourism and service sectors, mining and allied services are the second largest industry in Goa. Goa is rich in iron, bauxite, laterite, and silica. The state used to export 20 million tonnes of iron ore before the mining ban was put in effect by the Supreme Court in 2012. It had since been partially lifted in 2014, but the courts once again quashed 88 mining leases in 2018. Sesa Goa (Formerly Vedanta) and Dempo are prominent companies involved in mining.  

Emissions Inventory GridAgriculture and fishing are livelihood options for a sizeable portion of the populace. Rice is the main agricultural crop, followed by areca, cashew, and coconut. Medium scale industries include the manufacturing of pesticides, fertilisers, tyres, tubes, footwear, chemicals, pharmaceuticals, wheat products, steel rolling, fruits and fish canning, cashew nuts, textiles, brewery products. Solid waste management is another issue that needs urgent attention in Goa, much of the pollution at the local level in the state comes from burning of garbage. 

The Panjim metro covers an area of 76 versus the main Panjim city that is tiny at 20 Population according to the 2011 census for the larger metro area is 114,000. To assess Panjim’s air quality, we selected an airshed covering 30km x 40km. This domain is further segregated into 1km grids, to study the spatial variations in the emission and the pollution loads.

Meteorology fields are important as they have a direct impact on air pollution concentrations. During periods of high precipitation or high speed winds, emissions from a city are swept away and do not have an impact on concentrations. On the other hand, during the winter months when temperatures and inversion heights are low, there is a greater impact of emissions on pollution concentrations. Low temperatures also affect behaviour through the need for space and water heating – which in turn has increases emissions.

We processed the NCEP Reanalysis global meteorological fields from 2010 to 2018 through the 3D-WRF meteorological model. A summary of the data for one year, averaged for the city’s airshed is presented below by month. Download the processed data which includes information on year, month, day, hour, precipitation (mm/hour), mixing height (m), temperature (C), wind speed (m/sec), and wind direction (degrees) – key parameters which determine the intensity of dispersion of emissions.

WRF Meteorology for Panjim

Multi-Pollutant Emission Inventory

We compiled an emissions inventory for the Panjim region for the following pollutants – sulfur dioxide (SO2), nitrogen oxides (NOx), carbon monoxide (CO), non-methane volatile organic compounds (NMVOCs), carbon dioxide (CO2); and particulate matter (PM) in four bins (a) coarse PM with size fraction between 2.5 and 10 μm (b) fine PM with size fraction less than 2.5 μm (c) black carbon (BC) and (d) organic carbon (OC), for year 2015 and projected to 2030. In Phase 1, base year for all the calculations was 2015. In Phase 2, all the calculations are updated for year 2018.

Emissions Inventory We customized the SIM-air family of tools to fit the base information collated from disparate sources. Apart from the official reports, resource material ranges from GIS databases of land use, land cover, roads and rail lines, water bodies, built up area (represented in the adjacent figure), commercial activities (such as hotels, hospitals, kiosks, restaurants, malls, cinema complexes, traffic intersections, worship points, industrial hubs, and telecom towers), to population density and meteorology at the finest spatial resolution possible (1-km). A detailed description of these resources is published as a journal article in 2019, which also includes a summary of baselines and pollution analysis for 20 Indian cities.

This emissions inventory is based on available local activity and fuel consumption estimates for the selected urban airshed (represented in the grid above). This information is collated from multiple agencies ranging from the central pollution control board, state pollution control board, census bureau, national sample survey office, ministry of road transport and highways, annual survey of industries, central electrical authority, ministry of heavy industries, and municipal waste management, and publications from academic and non-governmental institutions.

For the road transport emissions inventory, besides the total number of vehicles and their usage information, we also utilized vehicle speed information to spatially and temporally allocate the estimated emissions to the respective grids. This is a product of google maps services. For the city of Panjim, we extracted the speed information for representative routes across the city for multiple days. This data is summarized below for a quick look.

Emissions Inventory Emissions Inventory

The Mormugao Port Trust near Vasco is a major port handling exports of ore, imports of coal for steel plants in Karnataka and Maharashtra and container shipments. Adani, Jindal and Vedanta operate berths that handle coal imports at the MPT and then direct the coal to the pig iron plant in goa or the steel plants in Karnataka. Everyday, 34,200 tonnes of coal is transported through the rail route from Mormugao port via Vasco, Margao and Kulem into Karnataka through Hubli and Hospet. As the capacity of the port is being increased from 12 million to 55 million metric tonnes, there is a push to build Goa’s freight infrastructure as well. The Shipping Ministry, meanwhile, has identified 17 coal-backed power plants and 22 steel plants as “potential clients” in the pipeline for the Mormugao Port at Bellary, Hospet, Belgaum, Hubli, Kudgi, Solapur, Bijapur, Gulbarga, Vijayanagara and Raichur, according to the 2016 Sagarmala report for coastal upgradation. There is a massive container transhipment yard being built outside the city of Margao, as well as expansion of the highway through the state into Maharashtra, Karwar and Belgaum (Karnataka).

Emission Sources in panjim Emission Sources in panjim Emission Sources in panjim
Emission Sources in panjim Emission Sources in panjim Emission Sources in panjim

The summary for a city’s emissions inventory does not include natural emission sources (like dust storms, lightning, and seasalt) and seasonal open (agricultural and forest) fires. However, these are included in the overall chemical transport modeling in the national scale simulations. These emission sources are accounted in the concentration calculation as an external (also known as boundary or long-range) contribution to the city’s air quality.

Projections to 2030 under the business as usual scenario are influenced by the city’s social, economic, landuse, urban, and industrial layout and hence the projected (increasing and decreasing) rates that we assume are an estimate only. We based the vehicle growth rate on the sales projection numbers; industrial growth on the gross domestic product of the state; domestic sector, construction activities, brick demand, diesel usage in the generator sets, and open waste burning on population growth rates and notes from the municipalities on plans to implement waste management programs. We used these estimates to evaluate the trend in the total emissions and their likely impact on ambient PM2.5 concentrations through 2030.

The emissions inventory was then spatially segregated at a 0.01° grid resolution in longitude and latitude (equivalent of 1 km) to create a spatial map of emissions for each pollutant (PM2.5, PM10, SO2, NOx, CO and VOCs). The gridded PM2.5 emissions and the total (shares by sector) emissions are presented below.

Gridded PM2.5 Emissions (2018 and 2030)

Emissions Inventory Emissions Inventory

Total PM2.5 Emissions by Sector 2018-2030

Emissions Inventory Emissions Inventory Emissions Inventory

TRANS = transport emissions from road, rail, aviation, and shipping (for coastal cities); RESIDEN = residential emissions from cooking, heating, and lighting activities; INDUS = industrial emissions from small, medium, and heavy industries (including power generation); ALL.DUST = dust emissions from road re-suspension and construction activities; W.BURN = open waste burning emissions; DG.SETS = diesel generator set emissions; B.KILNS = brick kiln emissions (not included in the industrial emissions)

Total Estimated Emissions by Sector for 2018 (units – tons/year)

Transport emissions from road, rail, aviation, and shipping (for coastal cities) 2,550 1,900 750800 17,600 59,550 14,500 4,750
Residential emissions from cooking, heating, and lighting activities100100-50- 1,800 200-
Industrial emissions from small, medium, and heavy industries (including power generation) 1,450 1,500 400250 7,950 5,650 11,600 2,800
Dust emissions from road re-suspension and construction activities 1,800 11,000 ------
Open waste burning emissions45050050300- 2,200 450-
Diesel generator set emissions400450200100 2,800 8,900 4,000 50
Brick kiln emissions (not included in the industrial emissions)--------
6,750 15,450 1,400 1,500 28,350 78,100 30,750 7,600

Chemical Transport Modeling

We calculated the ambient PM2.5 concentrations and the source contributions, using gridded emissions inventory, 3D meteorological data (from WRF), and the CAMx regional chemical transport model. The model simulates concentrations at 0.01° grid resolution and sector contributions for the urban area, which include contributions from primary emissions, secondary sources via chemical reactions, and long range transport via boundary conditions (represented as “boundary” in the pie graph below).

panjim PM2.5 Monthly variations panjim PM2.5 Monthly variations

The ribbon graph shows the variation for average PM2.5 pollution by month. Due to precipitation during the monsoon, usually pollution levels dip and may fall within national air pollution standards, however most cities are unable to attain these standards at other times of the year. We consolidated all the PM2.5 data from the continuous monitoring stations operating within the modeling domain, for the period of 2018-19 and compared against the model results. The scatter plot presents a comparison of 24-hr average PM2.5 concentrations by month. Modeled data is for the urban parts of the city.

The following is a map of annual average PM2.5 pollution for the city of Panjim. The main sources contributing towards PM2.5 in 2018 are in the pie-chart on the left. The change in contributions in 2030 from different sources are shown on the right.

Panjim PM2.5 Source Contributions Panjim Ambient PM2.5 Concentrations Panjim PM2.5 Source Contributions

There is a temporal variation in source contributions and spatial contributions depending on meteorological factors. We have a map of monthly average PM2.5 levels as well as their source contributions for every month in the charts below.

Panjim PM2.5 Monthly Concentrations Panjim PM2.5 Source Apportionment

Satellite Data Derived Surface PM2.5 Concentrations

The results of satellite data derived concentrations are useful for evaluating annual trends in pollution levels and are not a proxy for on-ground monitoring networks. This data is estimated using satellite feeds and global chemical transport models. Satellites are not measuring one location all the time, instead, a combination of satellites provide a cache of measurements that are interpreted using global chemical transport models (GEOS-Chem) to represent the vertical mix of pollution and estimate ground-based concentrations with the help of previous ground-based measurements. The global transport models rely on gridded emission estimates for multiple sectors to establish a relationship with satellite observations over multiple years. These databases were also used to study the global burden of disease, which estimated air pollution as the top 10 causes of premature mortality and morbidity in India. A summary of PM2.5 concentrations for the period of 1998 to 2016 for the city of Panjim is presented below. The global PM2.5 files are available for download and further analysis @ Dalhousie University.

The graphs for other district PM2.5 concentrations for this period, maps of national averages, and year-wise changes are available here. The data for district level PM2.5 concentrations for 1998-2016 period for can downloaded here.


We present below a summary of the ambient monitoring data available under the National Ambient Monitoring Program (NAMP), operated and maintained by the Central Pollution Control Board (CPCB, New Delhi, India). In Panjim, as of November 2018, there are 0 continuous and 2 manual air quality stations in operation. An archive of all the data from the NAMP network from stations across India for the period of 2011-2015 is available here.

PM2.5 Source Contributions Ambient PM2.5 Concentrations PM2.5 Source Contributions

Real time PM2.5 pollution levels in and around the city of Panjim

Resource Material

  • CPCB repository of continuous air monitoring data (Link)
  • CPCB List of non-attainment cities (Link)
  • “Revised City Development Plan for Panaji, 2041”, CRISIL Risk and Infrastructure Solutions Limited. (2015) (Link)
  • Panaji Smart City Proposal (Link)
  • “Panaji Solar City Master Plan”, Darashaw and Company. (2015) (Link)
  • “Low Emission Development Strategies (LEDS) for Panaji City”, ICLEI South Asia. (2015) (Link)
  • Air Quality Monitoring, Goa State Pollution Control Board (Link)
  • “Interim Report: Particulate Source Apportionment Study for Air Quality in Vasco-da-Gama, Goa”, MPT, GSPCB and IIT B. (2019) (Link)
  • “Emission Scenarios and Health Impacts of Air Pollutants in Goa”, Sharma, S. and Patil, V. K. (2016) (Journal Article Link)
  • “Depreciation in Ambient Air Quality in Iron Ore Mining Region of Goa”, Gurdeep, S. and Atahar, P. (2015) (Journal Article Link)

Back to the APnA page.

All the analysis and results are sole responsibility of the authors @ UrbanEmissions.Info. Please send you comments and questions to