Adjoint sensitivity analysis allows to assess the areas that have the largest impact on a given receptor site. The adjoint version of the Community Multiscale Air Quality (CMAQ v4.5) model was employed to perform a sensitivity analysis of ground level ozone for the episodic event of June 24, 2003, in the city of Athens assuming as a receptor site that of Agia Paraskevi Station. The 3-dimensional meteorology fields calculated using the Mesoscale Model 5 (MM5, Penn State University version 3.7.2) were used to produce high resolution daily air emissions inventories for the main anthropogenic and biogenic pollutants with 1-hour time step by an in-house built processor named EMISLAB. The meteorological prediction fields in combination with the emissions inventories were consequently fed as inputs to the CMAQ model. The ozone sensitivities were obtained with respect to pollutant concentrations and emissions. The distribution of the sensitivities in the computational domain for different times delineated the regions where perturbations in some concentrations would result in significant changes in the ozone concentrations in the area of interest (Agia Paraskevi, in this case) at the final time. The investigation yielded that the most significant influences were the transported O3 and NOx concentrations from the industrial area in the northern parts of the city and the road traffic from the city centre.
Three-dimensional atmospheric chemical transport models (CTMs) are tools used by the scientific community to predict the temporal and spatial distribution of air pollutants. Recent advancements in computational methods and resources have allowed the calculation of sensitivities of pollutant concentrations with respect to a large number of model parameters, improving our understanding of the changing chemical state of key atmospheric constituents in relation to their sources and sinks. The sensitivity analysis results in the calculation of the sensitivity coefficients (derivatives) of model outputs with respect to various inputs. Such sensitivity information can be used effectively in emission control strategies and air quality management.
There are two general techniques for computing sensitivity analysis, the forward and backward methods. In the forward approach, any perturbation in one or more input variables is propagated forward in time into various receptors [1,2]. In the backward (adjoint) mode, a perturbbation in a receptor based metric is propagated backward in time through an auxiliary set of equations resulting in the calculation of its sensitivities with respect to a large number of perturbed input parameters [
The mathematical formulation of the adjoint sensitivity for nonlinear dynamical systems has been presented by [4-7]. The mathematical description of the adjoint sensitivity method applied to three-dimensional air transport and chemical models could be found in [
This paper presents an adjoint sensitivity analysis of ozone concentrations with respect to emissions performed for the city of Athens and the Greater Athens Area (GAA), at the receptor site of Agia Paraskevi, using the adjoint version of the Community Multiscale Air Quality (CMAQ v4.5) model [
CMAQ is a 3-D atmospheric chemical transport model [
The computational domain is designed to cover an area of 103 ´ 103 km2 with 1 km ´ 1 km horizontal resolution (
The MM5 model was initialized at 00:00 UTC on June 23, 2003 using input meteorological data from the National Centres for Environmental Prediction (NCEP) Global Forecasting System (GFS) (6-hours resolution).
In general, the calculated meteorological fields showed that the both days were characterized by sea breeze conditions and rather low near surface winds that favoured the air pollution episode formation (
Firstly, the anthropogenic and biogenic emissions of NO, NO2, NMVOC, SO2 and PM were obtained using EMISLAB, on an hourly basis for the specific grid (1 ´ 1 km2) (
The chemistry-transport model CMAQ was started at 00:00 UTC on June 23, 2003 and performed a simulation of 48 hours. An example of the calculated O3 concentration profile at Agia Paraskevi (suburban) station is shown in
One suburban station Agia Paraskevi in the city of Athens has been selected for the adjoint sensitivity analysis which according to the annual report of the Ministry of
the Environment exhibits the highest amount of episodic events (i.e. exceedances of maximum 8-h average value) on an annual basis. The adjoint model of CMAQ was run 12 hours backwards of the target time (output at 12:00, 09:00, 06:00, 03:00 UTC, June 24, 2003). The sensitivities of ground level ozone concentrations (model-calculated) with respect to concentrations and emissions of NOx, CO, VOCs have been obtained. The distribution of the sensitivities in the computational domain for different times, (iso-surfaces) delineate the influence regions, i.e., areas where perturbations in some concentrations will result in significant changes in the ozone concentrations in the area of interest (Agia Paraskevi) at the final time. The calculated sensitivities of ozone with respect to concentrations of gas pollutants (O3, NO2 and SO2) are shown in Figures 8-10. The calculated sensitivities of O3 with respect to emissions of pollutants (NO, NO2) are shown in
An adjoint sensitivity study for ground ozone concentrations in Athens has been performed with the adjoint model of CMAQ for the receptor station of Agia Paraskevi (June 24, 2003). The difference on the areas of influence is evident depending on the prevailing wind field. Concerning the weather conditions of wind breeze which correspond to high O3 levels in the city, the impact of the large emissions sources can be easily traced (city centre/traffic, airport and shipping emissions). Transported pollutant concentrations (O3, NOx and SO2) influence the ozone concentration values at the receptor. The influence of NMVOC emissions on ozone at Agia Paraskevi was not found to be important, although the NMVOC concentrations were high during the particular pollution episode. The study indicated that focus must be placed on the improvement and updating of the high resolution emission inventories particularly that concerning the NOx emissions. Moreover, since the meteorological influence on the emissions and pollution levels is significant, the meteorological model over the urban canopy must be configured and evaluated against observational data.