Parallel performance and scalability experiments with the Danish Eulerian Model on the EPCC supercomputers

The Danish Eulerian Model (DEM) is a powerful air pollution model, designed to calculate the concentrations of various dangerous species over a large geographical region (e.g. Europe). It takes into account the main physical and chemical processes between these species, the actual meteorological conditions, emissions, etc.. This is a huge computational task and requires significant resources of storage and CPU time. Parallel computing is essential for the efficient practical use of the model. Some efficient parallel versions of the model were created over the past several years. A suitable parallel version of DEM by using the Message Passing Interface library (AIPI) was implemented on two powerful supercomputers of the EPCC - Edinburgh, available via the HPC-Europa programme for transnational access to research infrastructures in EC: a Sun Fire E15K and an IBM HPCx cluster. Although the implementation is in principal, the same for both supercomputers, few modifications had to be done for successful porting of the code on the IBM HPCx cluster. Performance analysis and parallel optimization was done next. Results from bench marking experiments will be presented in this paper. Another set of experiments was carried out in order to investigate the sensitivity of the model to variation of some chemical rate constants in the chemical submodel. Certain modifications of the code were necessary to be done in accordance with this task. The obtained results will be used for further sensitivity analysis Studies by using Monte Carlo simulation.

A collaborative working environment for a large scale environmental model

Large scientific applications are usually developed, tested and used by a group of geographically dispersed scientists. The problems associated with the remote development and data sharing could be tackled by using collaborative working environments. There are various tools and software to create collaborative working environments. Some software frameworks, currently available, use these tools and software to enable remote job submission and file transfer on top of existing grid infrastructures. However, for many large scientific applications, further efforts need to be put to prepare a framework which offers application-centric facilities. Unified Air Pollution Model (UNI-DEM), developed by Danish Environmental Research Institute, is an example of a large scientific application which is in a continuous development and experimenting process by different institutes in Europe. This paper intends to design a collaborative distributed computing environment for UNI-DEM in particular but the framework proposed may also fit to many large scientific applications as well.

Parallel implementation and one year experiments with the Danish Eulerian Model

Large scale air pollution models are powerful tools, designed to meet the increasing demand in different environmental studies. The atmosphere is the most dynamic component of the environment, where the pollutants can be moved quickly on far distnce. Therefore the air pollution modeling must be done in a large computational domain. Moreover, all relevant physical, chemical and photochemical processes must be taken into account. In such complex models operator splitting is very often applied in order to achieve sufficient accuracy as well as efficiency of the numerical solution. The Danish Eulerian Model (DEM) is one of the most advanced such models. Its space domain (4800 × 4800 km) covers Europe, most of the Mediterian and neighboring parts of Asia and the Atlantic Ocean. Efficient parallelization is crucial for the performance and practical capabilities of this huge computational model. Different splitting schemes, based on the main processes mentioned above, have been implemented and tested with respect to accuracy and performance in the new version of DEM. Some numerical results of these experiments are presented in this paper.

Interaction of nitrogen dioxide (NO2) with a monolayer of oleic acid at the air–water interface – A simple proxy for atmospheric aerosol

The reactions between atmospheric oxidants and organic amphiphiles at the air water interface of an aerosol droplet may affect the size and critical supersaturation required for cloud droplet formation. We demonstrate that no reaction occurs between gaseous nitrogen dioxide (1000 ppm in air) and a monolayer of an insoluble amphiphile, oleic acid (cis-9-octadecenoic acid), at the air water interface which removes material from the air water interface. We present evidence that the NO2 isomerises the cis-9-octadecenoic (oleic) acid to trans-9-octadecenoic (elaidic) acid. The study presented here is important for future and previous studies of (1) the reaction between the nitrate radical, NO3, and thin organic films as NO2 is usually present in high concentrations in these experimental systems and (2) the effect of NO2 air pollution on the unsaturated fatty acids and lipids found at the air liquid surface of human lung lining fluid.

Fair weather atmospheric electricity

Abstract. Not long after Franklin’s iconic studies, an atmospheric electric field was discovered in “fair weather” regions, well away from thunderstorms. The origin of the fair weather field was sought by Lord Kelvin, through development of electrostatic instrumentation and early data logging techniques, but was ultimately explained through the global circuit model of C.T.R. Wilson. In Wilson’s model, charge exchanged by disturbed weather electrifies the ionosphere, and returns via a small vertical current density in fair weather regions. New insights into the relevance of fair weather atmospheric electricity to terrestrial and planetary atmospheres are now emerging. For example, there is a possible role of the global circuit current density in atmospheric processes, such as cloud formation. Beyond natural atmospheric processes, a novel practical application is the use of early atmospheric electrostatic investigations to provide quantitative information on past urban air pollution.

Terrestrial biogeochemical feedbacks in the climate system

The terrestrial biosphere is a key regulator of atmospheric chemistry and climate. During past periods of climate change, vegetation cover and interactions between the terrestrial biosphere and atmosphere changed within decades. Modern observations show a similar responsiveness of terrestrial biogeochemistry to anthropogenically forced climate change and air pollution. Although interactions between the carbon cycle and climate have been a central focus, other biogeochemical feedbacks could be as important in modulating future climate change. Total positive radiative forcings resulting from feedbacks between the terrestrial biosphere and the atmosphere are estimated to reach up to 0.9 or 1.5 W m−2 K−1 towards the end of the twenty-first century, depending on the extent to which interactions with the nitrogen cycle stimulate or limit carbon sequestration. This substantially reduces and potentially even eliminates the cooling effect owing to carbon dioxide fertilization of the terrestrial biota. The overall magnitude of the biogeochemical feedbacks could potentially be similar to that of feedbacks in the physical climate system, but there are large uncertainties in the magnitude of individual estimates and in accounting for synergies between these effects.

Global and regional temperature-change potentials for near-term climate forcers

We examine the climate effects of the emissions of near-term climate forcers (NTCFs) from 4 continental regions (East Asia, Europe, North America and South Asia) using radiative forcing from the task force on hemispheric transport of air pollution source-receptor global chemical transport model simulations. These simulations model the transport of 3 aerosol species (sulphate, particulate organic matter and black carbon) and 4 ozone precursors (methane, nitric oxides (NOx), volatile organic compounds and carbon monoxide). From the equilibrium radiative forcing results we calculate global climate metrics, global warming potentials (GWPs) and global temperature change potentials (GTPs) and show how these depend on emission region, and can vary as functions of time. For the aerosol species, the GWP(100) values are −37±12, −46±20, and 350±200 for SO2, POM and BC respectively for the direct effects only. The corresponding GTP(100) values are −5.2±2.4, −6.5±3.5, and 50±33. This analysis is further extended by examining the temperature-change impacts in 4 latitude bands. This shows that the latitudinal pattern of the temperature response to emissions of the NTCFs does not directly follow the pattern of the diagnosed radiative forcing. For instance temperatures in the Arctic latitudes are particularly sensitive to NTCF emissions in the northern mid-latitudes. At the 100-yr time horizon the ARTPs show NOx emissions can have a warming effect in the northern mid and high latitudes, but cooling in the tropics and Southern Hemisphere. The northern mid-latitude temperature response to northern mid-latitude emissions of most NTCFs is approximately twice as large as would be implied by the global average.

The influence of ozone precursor emissions from four world regions on tropospheric composition and radiative climate forcing

Ozone (O3) precursor emissions influence regional and global climate and air quality through changes in tropospheric O3 and oxidants, which also influence methane (CH4) and sulfate aerosols (SO42−). We examine changes in the tropospheric composition of O3, CH4, SO42− and global net radiative forcing (RF) for 20% reductions in global CH4 burden and in anthropogenic O3 precursor emissions (NOx, NMVOC, and CO) from four regions (East Asia, Europe and Northern Africa, North America, and South Asia) using the Task Force on Hemispheric Transport of Air Pollution Source-Receptor global chemical transport model (CTM) simulations, assessing uncertainty (mean ± 1 standard deviation) across multiple CTMs. We evaluate steady state O3 responses, including long-term feedbacks via CH4. With a radiative transfer model that includes greenhouse gases and the aerosol direct effect, we find that regional NOx reductions produce global, annually averaged positive net RFs (0.2 ± 0.6 to 1.7 ± 2 mWm−2/Tg N yr−1), with some variation among models. Negative net RFs result from reductions in global CH4 (−162.6 ± 2 mWm−2 for a change from 1760 to 1408 ppbv CH4) and regional NMVOC (−0.4 ± 0.2 to −0.7 ± 0.2 mWm−2/Tg C yr−1) and CO emissions (−0.13 ± 0.02 to −0.15 ± 0.02 mWm−2/Tg CO yr−1). Including the effect of O3 on CO2 uptake by vegetation likely makes these net RFs more negative by −1.9 to −5.2 mWm−2/Tg N yr−1, −0.2 to −0.7 mWm−2/Tg C yr−1, and −0.02 to −0.05 mWm−2/Tg CO yr−1. Net RF impacts reflect the distribution of concentration changes, where RF is affected locally by changes in SO42−, regionally to hemispherically by O3, and globally by CH4. Global annual average SO42− responses to oxidant changes range from 0.4 ± 2.6 to −1.9 ± 1.3 Gg for NOx reductions, 0.1 ± 1.2 to −0.9 ± 0.8 Gg for NMVOC reductions, and −0.09 ± 0.5 to −0.9 ± 0.8 Gg for CO reductions, suggesting additional research is needed. The 100-year global warming potentials (GWP100) are calculated for the global CH4 reduction (20.9 ± 3.7 without stratospheric O3 or water vapor, 24.2 ± 4.2 including those components), and for the regional NOx, NMVOC, and CO reductions (−18.7 ± 25.9 to −1.9 ± 8.7 for NOx, 4.8 ± 1.7 to 8.3 ± 1.9 for NMVOC, and 1.5 ± 0.4 to 1.7 ± 0.5 for CO). Variation in GWP100 for NOx, NMVOC, and CO suggests that regionally specific GWPs may be necessary and could support the inclusion of O3 precursors in future policies that address air quality and climate change simultaneously. Both global net RF and GWP100 are more sensitive to NOx and NMVOC reductions from South Asia than the other three regions.

A multimodel assessment of the influence of regional anthropogenic emission reductions on aerosol direct radiative forcing and the role of intercontinental transport

In this study, we assess changes of aerosol optical depth (AOD) and direct radiative forcing (DRF) in response to the reduction of anthropogenic emissions in four major pollution regions in the Northern Hemisphere by using results from nine global models in the framework of the Hemispheric Transport of Air Pollution (HTAP). DRF at top of atmosphere (TOA) and surface is estimated based on AOD results from the HTAP models and AOD-normalized DRF (NDRF) from a chemical transport model. The multimodel results show that, on average, a 20% reduction of anthropogenic emissions in North America, Europe, East Asia, and South Asia lowers the global mean AOD (all-sky TOA DRF) by 9.2% (9.0%), 3.5% (3.0%), and 9.4% (10.0%) for sulfate, particulate organic matter (POM), and black carbon (BC), respectively. Global annual average TOA all-sky forcing efficiency relative to particle or gaseous precursor emissions from the four regions (expressed as multimodel mean ± one standard deviation) is  ±3.5 ±0.8,  ±4.0 ±1.7, and 29.5 ±18.1mWm ±2 per Tg for sulfate (relative to SO2), POM, and BC, respectively. The impacts of the regional emission reductions on AOD and DRF extend well beyond the source regions because of intercontinental transport (ICT). On an annual basis, ICT accounts for 11 ±5% to 31 ±9% of AOD and DRF in a receptor region at continental or subcontinental scale, with domestic emissions accounting for the remainder, depending on regions and species. For sulfate AOD, the largest ICT contribution of 31 ±9% occurs in South Asia, which is dominated by the emissions from Europe. For BC AOD, the largest ICT contribution of 28 ±18% occurs in North America, which is dominated by the emissions from East Asia. The large spreads among models highlight the need to improve aerosol processes in models, and evaluate and constrain models with observations.

Observing wind, aerosol particles, clouds and precipitation: Finland's new ground-based remote-sensing network

The Finnish Meteorological Institute, in collaboration with the University of Helsinki, has established a new ground-based remote-sensing network in Finland. The network consists of five topographically, ecologically and climatically different sites distributed from southern to northern Finland. The main goal of the network is to monitor air pollution and boundary layer properties in near real time, with a Doppler lidar and ceilometer at each site. In addition to these operational tasks, two sites are members of the Aerosols, Clouds and Trace gases Research InfraStructure Network (ACTRIS); a Ka band cloud radar at Sodankylä will provide cloud retrievals within CloudNet, and a multi-wavelength Raman lidar, PollyXT (POrtabLe Lidar sYstem eXTended), in Kuopio provides optical and microphysical aerosol properties through EARLINET (the European Aerosol Research Lidar Network). Three C-band weather radars are located in the Helsinki metropolitan area and are deployed for operational and research applications. We performed two inter-comparison campaigns to investigate the Doppler lidar performance, compare the backscatter signal and wind profiles, and to optimize the lidar sensitivity through adjusting the telescope focus length and data-integration time to ensure sufficient signal-to-noise ratio (SNR) in low-aerosol-content environments. In terms of statistical characterization, the wind-profile comparison showed good agreement between different lidars. Initially, there was a discrepancy in the SNR and attenuated backscatter coefficient profiles which arose from an incorrectly reported telescope focus setting from one instrument, together with the need to calibrate. After diagnosing the true telescope focus length, calculating a new attenuated backscatter coefficient profile with the new telescope function and taking into account calibration, the resulting attenuated backscatter profiles all showed good agreement with each other. It was thought that harsh Finnish winters could pose problems, but, due to the built-in heating systems, low ambient temperatures had no, or only a minor, impact on the lidar operation – including scanning-head motion. However, accumulation of snow and ice on the lens has been observed, which can lead to the formation of a water/ice layer thus attenuating the signal inconsistently. Thus, care must be taken to ensure continuous snow removal.

Radiative forcing and climate metrics for ozone precursor emissions: the impact of multi-model averaging

Multi-model ensembles are frequently used to assess understanding of the response of ozone and methane lifetime to changes in emissions of ozone precursors such as NOx, VOCs (volatile organic compounds) and CO. When these ozone changes are used to calculate radiative forcing (RF) (and climate metrics such as the global warming potential (GWP) and global temperature-change potential (GTP)) there is a methodological choice, determined partly by the available computing resources, as to whether the mean ozone (and methane) concentration changes are input to the radiation code, or whether each model's ozone and methane changes are used as input, with the average RF computed from the individual model RFs. We use data from the Task Force on Hemispheric Transport of Air Pollution source–receptor global chemical transport model ensemble to assess the impact of this choice for emission changes in four regions (East Asia, Europe, North America and South Asia). We conclude that using the multi-model mean ozone and methane responses is accurate for calculating the mean RF, with differences up to 0.6% for CO, 0.7% for VOCs and 2% for NOx. Differences of up to 60% for NOx 7% for VOCs and 3% for CO are introduced into the 20 year GWP. The differences for the 20 year GTP are smaller than for the GWP for NOx, and similar for the other species. However, estimates of the standard deviation calculated from the ensemble-mean input fields (where the standard deviation at each point on the model grid is added to or subtracted from the mean field) are almost always substantially larger in RF, GWP and GTP metrics than the true standard deviation, and can be larger than the model range for short-lived ozone RF, and for the 20 and 100 year GWP and 100 year GTP. The order of averaging has most impact on the metrics for NOx, as the net values for these quantities is the residual of the sum of terms of opposing signs. For example, the standard deviation for the 20 year GWP is 2–3 times larger using the ensemble-mean fields than using the individual models to calculate the RF. The source of this effect is largely due to the construction of the input ozone fields, which overestimate the true ensemble spread. Hence, while the average of multi-model fields are normally appropriate for calculating mean RF, GWP and GTP, they are not a reliable method for calculating the uncertainty in these fields, and in general overestimate the uncertainty.

The effects of springtime mid-latitude storms on trace gas composition determined from the MACC reanalysis

The relationship between springtime air pollution transport of ozone (O3) and carbon monoxide (CO) and mid-latitude cyclones is explored for the first time using the Monitoring Atmospheric Composition and Climate (MACC) reanalysis for the period 2003–2012. In this study, the most intense spring storms (95th percentile) are selected for two regions, the North Pacific (NP) and the North Atlantic (NA). These storms (∼60 storms over each region) often track over the major emission sources of East Asia and eastern North America. By compositing the storms, the distributions of O3 and CO within a "typical" intense storm are examined. We compare the storm-centered composite to background composites of "average conditions" created by sampling the reanalysis data of the previous year to the storm locations. Mid-latitude storms are found to redistribute concentrations of O3 and CO horizontally and vertically throughout the storm. This is clearly shown to occur through two main mechanisms: (1) vertical lifting of CO-rich and O3-poor air isentropically, from near the surface to the mid- to upper-troposphere in the region of the warm conveyor belt; and (2) descent of O3-rich and CO-poor air isentropically in the vicinity of the dry intrusion, from the stratosphere toward the mid-troposphere. This can be seen in the composite storm's life cycle as the storm intensifies, with area-averaged O3 (CO) increasing (decreasing) between 200 and 500 hPa. The influence of the storm dynamics compared to the background environment on the composition within an area around the storm center at the time of maximum intensity is as follows. Area-averaged O3 at 300 hPa is enhanced by 50 and 36% and by 11 and 7.6% at 500 hPa for the NP and NA regions, respectively. In contrast, area-averaged CO at 300 hPa decreases by 12% for NP and 5.5% for NA, and area-averaged CO at 500 hPa decreases by 2.4% for NP while there is little change over the NA region. From the mid-troposphere, O3-rich air is clearly seen to be transported toward the surface, but the downward transport of CO-poor air is not discernible due to the high levels of CO in the lower troposphere. Area-averaged O3 is slightly higher at 1000 hPa (3.5 and 1.8% for the NP and NA regions, respectively). There is an increase of CO at 1000 hPa for the NP region (3.3%) relative to the background composite and a~slight decrease in area-averaged CO for the NA region at 1000 hPa (-2.7%).

Modelled and observed changes in aerosols and surface solar radiation over Europe between 1960 and 2009

Substantial changes in anthropogenic aerosols and precursor gas emissions have occurred over recent decades due to the implementation of air pollution control legislation and economic growth. The response of atmospheric aerosols to these changes and the impact on climate are poorly constrained, particularly in studies using detailed aerosol chemistry–climate models. Here we compare the HadGEM3-UKCA (Hadley Centre Global Environment Model-United Kingdom Chemistry and Aerosols) coupled chemistry–climate model for the period 1960–2009 against extensive ground-based observations of sulfate aerosol mass (1978–2009), total suspended particle matter (SPM, 1978–1998), PM10 (1997–2009), aerosol optical depth (AOD, 2000–2009), aerosol size distributions (2008–2009) and surface solar radiation (SSR, 1960–2009) over Europe. The model underestimates observed sulfate aerosol mass (normalised mean bias factor (NMBF) = −0.4), SPM (NMBF = −0.9), PM10 (NMBF = −0.2), aerosol number concentrations (N30 NMBF = −0.85; N50 NMBF = −0.65; and N100 NMBF = −0.96) and AOD (NMBF = −0.01) but slightly overpredicts SSR (NMBF = 0.02). Trends in aerosol over the observational period are well simulated by the model, with observed (simulated) changes in sulfate of −68 % (−78 %), SPM of −42 % (−20 %), PM10 of −9 % (−8 %) and AOD of −11 % (−14 %). Discrepancies in the magnitude of simulated aerosol mass do not affect the ability of the model to reproduce the observed SSR trends. The positive change in observed European SSR (5 %) during 1990–2009 ("brightening") is better reproduced by the model when aerosol radiative effects (ARE) are included (3 %), compared to simulations where ARE are excluded (0.2 %). The simulated top-of-the-atmosphere aerosol radiative forcing over Europe under all-sky conditions increased by > 3.0 W m−2 during the period 1970–2009 in response to changes in anthropogenic emissions and aerosol concentrations.

Emission of polycyclic aromatic hydrocarbons from gasohol and ethanol vehicles

The exhaust emission of the polycyclic aromatic hydrocarbons (PAHs) considered toxic to human health were investigated on two spark ignition light duty vehicles, one being gasohol (Gasohol, in Brazil, is the generic denomination for mixtures of pure gasoline plus 20-25% of anhydrous ethyl alcohol fuel (AEAF).)-fuelled and the other a flexible-fuel vehicle fuelled with hydrated ethanol. The influence of fuel type and quality, aged lubricant oil type and use of fuel additives on the formation of these compounds was tested using standardized tests identical to US FTP-75 cycle. PAH sampling and chemical analysis followed the basic recommendations of method TO-13 (United States. Environmental Protection Agency, 1999. Compendium Method TO-13A - Determination of polycyclic Aromatic hydrocarbons (PAH) in Ambient Air Using Gas Chromatography/Mass Spectrometry (CG/MS). Center for environmental research information, Cincinnati, p. 78), with the necessary modification for this particular application. Results showed that the total PAH emission factor varied from 41.9 mu g km(-1) to 612 mu g km(-1) in the gasohol vehicle, and from 11.7 mu g km(-1) to 27.4 mu g km(-1) in the ethanol-fuelled vehicle, a significant difference in favor of the ethanol vehicle. Generally, emission of light molecular weight PAHs was predominant, while high molecular weights PAHs were not detected. In terms of benzo(a)pyrene toxicity equivalence, emission factors varied from 0.00984 mu g TEQ km(-1) to 4.61 mu g TEQ km(-1) for the gasohol vehicle and from 0.0117 mu g TEQ km(-1) to 0.0218 mu g TEQ km(-1) in the ethanol vehicle. For the gasohol vehicle, results showed that the use of fuel additive causes a significant increase in the emission of naphthalene and phenanthrene at a confidence level of 90% or higher; the use of rubber solvent on gasohol showed a reduction in the emission of naphthalene and phenanthrene at the same confidence level; the use of synthetic oil instead of mineral oil also contributed significantly to a decrease in the emission of naphthalene and fluorene. In relation to the ethanol vehicle, the same factors were tested and showed no statistically significant influence on PAH emission. (c) 2008 Elsevier Ltd. All rights reserved.