Search

Recently, a standard test case suite for 2-D linear transport on the sphere was proposed to assess important aspects of accuracy in geophysical fluid dynamics with a "minimal" set of idealized model configurations/runs/diagnostics. Here we present results from 19 state-of-the-art transport scheme formulations based on finite-difference/finite-volume methods as well as emerging (in the context of atmospheric/oceanographic sciences) Galerkin methods. Discretization grids range from traditional regular latitude–longitude grids to more isotropic domain discretizations such as icosahedral and cubed-sphere tessellations of the sphere. The schemes are evaluated using a wide range of diagnostics in idealized flow environments. Accuracy is assessed in single- and two-tracer configurations using conventional error norms as well as novel diagnostics designed for climate and climate–chemistry applications. In addition, algorithmic considerations that may be important for computational efficiency are reported on. The latter is inevitably computing platform dependent. The ensemble of results from a wide variety of schemes presented here helps shed light on the ability of the test case suite diagnostics and flow settings to discriminate between algorithms and provide insights into accuracy in the context of global atmospheric/ocean modeling. A library of benchmark results is provided to facilitate scheme intercomparison and model development. Simple software and data sets are made available to facilitate the process of model evaluation and scheme intercomparison.

We have analysed time-slice simulations from 17 global models, participating in the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), to explore changes in present-day (2000) hydroxyl radical (OH) concentration and methane (CH4) lifetime relative to preindustrial times (1850) and to 1980. A comparison of modeled and observation-derived methane and methyl chloroform lifetimes suggests that the present-day global multi-model mean OH concentration is overestimated by 5 to 10% but is within the range of uncertainties. The models consistently simulate higher OH concentrations in the Northern Hemisphere (NH) compared with the Southern Hemisphere (SH) for the present-day (2000; inter-hemispheric ratios of 1.13 to 1.42), in contrast to observation-based approaches which generally indicate higher OH in the SH although uncertainties are large. Evaluation of simulated carbon monoxide (CO) concentrations, the primary sink for OH, against ground-based and satellite observations suggests low biases in the NH that may contribute to the high north–south OH asymmetry in the models. The models vary widely in their regional distribution of present-day OH concentrations (up to 34%). Despite large regional changes, the multi-model global mean (mass-weighted) OH concentration changes little over the past 150 yr, due to concurrent increases in factors that enhance OH (humidity, tropospheric ozone, nitrogen oxide (NOx) emissions, and UV radiation due to decreases in stratospheric ozone), compensated by increases in OH sinks (methane abundance, carbon monoxide and non-methane volatile organic carbon (NMVOC) emissions). The large inter-model diversity in the sign and magnitude of preindustrial to present-day OH changes (ranging from a decrease of 12.7% to an increase of 14.6%) indicate that uncertainty remains in our understanding of the long-term trends in OH and methane lifetime. We show that this diversity is largely explained by the different ratio of the change in global mean tropospheric CO and NOx burdens (ΔCO/ΔNOx, approximately represents changes in OH sinks versus changes in OH sources) in the models, pointing to a need for better constraints on natural precursor emissions and on the chemical mechanisms in the current generation of chemistry-climate models. For the 1980 to 2000 period, we find that climate warming and a slight increase in mean OH (3.5 ± 2.2%) leads to a 4.3 ± 1.9% decrease in the methane lifetime. Analysing sensitivity simulations performed by 10 models, we find that preindustrial to present-day climate change decreased the methane lifetime by about four months, representing a negative feedback on the climate system. Further, we analysed attribution experiments performed by a subset of models relative to 2000 conditions with only one precursor at a time set to 1860 levels. We find that global mean OH increased by 46.4 ± 12.2% in response to preindustrial to present-day anthropogenic NOx emission increases, and decreased by 17.3 ± 2.3%, 7.6 ± 1.5%, and 3.1 ± 3.0% due to methane burden, and anthropogenic CO, and NMVOC emissions increases, respectively.

Results from simulations performed for the Atmospheric Chemistry and Climate Modeling Intercomparison Project (ACCMIP) are analysed to examine how OH and methane lifetime may change from present day to the future, under different climate and emissions scenarios. Present day (2000) mean tropospheric chemical lifetime derived from the ACCMIP multi-model mean is 9.8 ± 1.6 yr (9.3 ± 0.9 yr when only including selected models), lower than a recent observationally-based estimate, but with a similar range to previous multi-model estimates. Future model projections are based on the four Representative Concentration Pathways (RCPs), and the results also exhibit a large range. Decreases in global methane lifetime of 4.5 ± 9.1% are simulated for the scenario with lowest radiative forcing by 2100 (RCP 2.6), while increases of 8.5 ± 10.4% are simulated for the scenario with highest radiative forcing (RCP 8.5). In this scenario, the key driver of the evolution of OH and methane lifetime is methane itself, since its concentration more than doubles by 2100 and it consumes much of the OH that exists in the troposphere. Stratospheric ozone recovery, which drives tropospheric OH decreases through photolysis modifications, also plays a partial role. In the other scenarios, where methane changes are less drastic, the interplay between various competing drivers leads to smaller and more diverse OH and methane lifetime responses, which are difficult to attribute. For all scenarios, regional OH changes are even more variable, with the most robust feature being the large decreases over the remote oceans in RCP8.5. Through a regression analysis, we suggest that differences in emissions of non-methane volatile organic compounds and in the simulation of photolysis rates may be the main factors causing the differences in simulated present day OH and methane lifetime. Diversity in predicted changes between present day and future OH was found to be associated more strongly with differences in modelled temperature and stratospheric ozone changes. Finally, through perturbation experiments we calculated an OH feedback factor (F) of 1.24 from present day conditions (1.50 from 2100 RCP8.5 conditions) and a climate feedback on methane lifetime of 0.33 ± 0.13 yr K−1, on average. Models that did not include interactive stratospheric ozone effects on photolysis showed a stronger sensitivity to climate, as they did not account for negative effects of climate-driven stratospheric ozone recovery on tropospheric OH, which would have partly offset the overall OH/methane lifetime response to climate change.

We present the first reconstruction of the Northern Hemisphere (NH) high latitude atmospheric carbon monoxide (CO) mole fraction from Greenland firn air. Firn air samples were collected at three deep ice core sites in Greenland (NGRIP in 2001, Summit in 2006 and NEEM in 2008). CO records from the three sites agree well with each other as well as with recent atmospheric measurements, indicating that CO is well preserved in the firn at these sites. CO atmospheric history was reconstructed back to the year 1950 from the measurements using a combination of two forward models of gas transport in firn and an inverse model. The reconstructed history suggests that Arctic CO in 1950 was 140-150 nmol mol-1, which is higher than today's values. CO mole fractions rose by 10-15 nmol mol-1 from 1950 to the 1970s and peaked in the 1970s or early 1980s, followed by a ≈ 30 nmol mol-1 decline to today's levels. We compare the CO history with the atmospheric histories of methane, light hydrocarbons, molecular hydrogen, CO stable isotopes and hydroxyl radicals (OH), as well as with published CO emission inventories and results of a historical run from a chemistry-transport model. We find that the reconstructed Greenland CO history cannot be reconciled with available emission inventories unless unrealistically large changes in OH are assumed. We argue that the available CO emission inventories strongly underestimate historical NH emissions, and fail to capture the emission decline starting in the late 1970s, which was most likely due to reduced emissions from road transportation in North America and Europe. © Author(s) 2013.

An extensive set of measurements was made in and around Mexico City as part of the MILAGRO (Megacity Initiative: Local and Global Research Observations) experiments in March 2006. Simulations with the Model for Ozone and Related Chemical Tracers, version 4 (MOZART-4), a global chemical transport model, have been used to provide a regional context for these observations and assist in their interpretation. These MOZART-4 simulations reproduce the aircraft observations generally well, but some differences in the modeled volatile organic compounds (VOCs) from the observations result from incorrect VOC speciation assumed for the emission inventories. The different types of CO sources represented in the model have been "tagged" to quantify the contributions of regions outside Mexico, as well as the various emissions sectors within Mexico, to the regional air quality of Mexico. This analysis indicates open fires have some, but not a dominant, impact on the atmospheric composition in the region around Mexico City when averaged over the month. However, considerable variation in the fire contribution (2-15% of total CO) is seen during the month. The transport and photochemical aging of Mexico City emissions were studied using tags of CO emissions for each day, showing that typically the air downwind of Mexico City was a combination of many ages. Ozone production in MOZART-4 is shown to agree well with the net production rates from box model calculations constrained by the MILAGRO aircraft measurements. Ozone production efficiency derived from the ratio of O-x to NOz is higher in MOZART-4 than in the observations for moderately polluted air. OH reactivity determined from the MOZART-4 results shows the same increase in relative importance of oxygenated VOCs downwind of Mexico City as the reactivity inferred from the observations. The amount of ozone produced by emissions from Mexico City and surrounding areas has been quantified in the model by tracking NO emissions, showing little influence beyond Mexico's borders, and also relatively minor influence from fire emissions on the monthly average tropospheric ozone column.

The Model for Ozone and Related chemical Tracers, version 4 (MOZART-4) is an offline global chemical transport model particularly suited for studies of the troposphere. The updates of the model from its previous version MOZART-2 are described, including an expansion of the chemical mechanism to include more detailed hydrocarbon chemistry and bulk aerosols. Online calculations of a number of processes, such as dry deposition, emissions of isoprene and monoterpenes and photolysis frequencies, are now included. Results from an eight-year simulation (2000-2007) are presented and evaluated. The MOZART-4 source code and standard input files are available for download from the NCAR Community Data Portal (http://cdp.ucar.edu).

A preliminary investigation was carried out by the US Environmental Protection Agency (EPA) on the impact of regional ozone concentrations on climate change in the country. The main objective of the investigation was to enhance the ability of air quality managers to consider global change in their decisions through improved characterization of the potential effects of global change on air quality, including O3, PM, and Hg. An integrated assessment framework was designed that supported the research and development efforts within the EPA and the academic research professionals. This design recognized the challenges of bridging spatial scales, temporal scales, and disciplines that characterize the global change and regional air quality problem. A major component of the assessment approach involved the development and application of global to regional climate and air quality modeling systems.

Inclusion of fundamental ecological interactions between carbon and nitrogen cycles in the land component of an atmosphere-ocean general circulation model (AOGCM) leads to decreased carbon uptake associated with CO2 fertilization, and increased carbon uptake associated with warming of the climate system. The balance of these two opposing effects is to reduce the fraction of anthropogenic CO2 predicted to be sequestered in land ecosystems. The primary mechanism responsible for increased land carbon storage under radiatively forced climate change is shown to be fertilization of plant growth by increased mineralization of nitrogen directly associated with increased decomposition of soil organic matter under a warming climate, which in this particular model results in a negative gain for the climate-carbon feedback. Estimates for the land and ocean sink fractions of recent anthropogenic emissions are individually within the range of observational estimates, but the combined land plus ocean sink fractions produce an airborne fraction which is too high compared to observations. This bias is likely due in part to an underestimation of the ocean sink fraction. Our results show a significant growth in the airborne fraction of anthropogenic CO2 emissions over the coming century, attributable in part to a steady decline in the ocean sink fraction. Comparison to experimental studies on the fate of radio-labeled nitrogen tracers in temperate forests indicates that the model representation of competition between plants and microbes for new mineral nitrogen resources is reasonable. Our results suggest a weaker dependence of net land carbon flux on soil moisture changes in tropical regions, and a stronger positive growth response to warming in those regions, than predicted by a similar AOGCM implemented without land carbon-nitrogen interactions. We expect that the between-model uncertainty in predictions of future atmospheric CO2 concentration and associated anthropogenic climate change will be reduced as additional climate models introduce carbon-nitrogen cycle interactions in their land components.

A comprehensive numerical modeling framework was developed to estimate the effects of collective global changes upon ozone pollution in the US in 2050. The framework consists of the global climate and chemistry models, PCM (Parallel Climate Model) and MOZART-2 (Model for Ozone and Related Chemical Tracers v.2), coupled with regional meteorology and chemistry models, MM5 (Mesoscale Meteorological model) and CMAQ (Community Multi-scale Air Quality model). The modeling system was applied for two 10-year simulations: 1990ĝ€" 1999 as a present-day base case and 2045ĝ€"2054 as a future case. For the current decade, the daily maximum 8-h moving average (DM8H) ozone mixing ratio distributions for spring, summer and fall showed good agreement with observations. The future case simulation followed the Intergovernmental Panel on Climate Change (IPCC) A2 scenario together with business-as-usual US emission projections and projected alterations in land use, land cover (LULC) due to urban expansion and changes in vegetation. For these projections, US anthropogenic NOx (NO+NO2) and VOC (volatile organic carbon) emissions increased by approximately 6% and 50%, respectively, while biogenic VOC emissions decreased, in spite of warmer temperatures, due to decreases in forested lands and expansion of croplands, grasslands and urban areas. A stochastic model for wildfire emissions was applied that projected 25% higher VOC emissions in the future. For the global and US emission projection used here, regional ozone pollution becomes worse in the 2045ĝ€" 2054 period for all months. Annually, the mean DM8H ozone was projected to increase by 9.6 ppbv (22%). The changes were higher in the spring and winter (25%) and smaller in the summer (17%). The area affected by elevated ozone within the US continent was projected to increase; areas with levels exceeding the 75 ppbv ozone standard at least once a year increased by 38%. In addition, the length of the ozone season was projected to increase with more pollution episodes in the spring and fall. For selected urban areas, the system projected a higher number of pollution events per year and these events had more consecutive days when DM8H ozone exceed 75 ppbv.

We present a study of the sensitivity of isoprene emission calculations in a global chemistry transport model (CTM) to input land cover characteristics and analyze the impacts of changes in isoprene on the tropospheric budgets of atmospheric key species. The CTM Model for Ozone and Related Chemical Species, version 4 (MOZART-4) includes the online calculation of isoprene emissions based on the Model of Emissions of Gases and Aerosols from Nature (MEGAN), which is driven by three different land parameter inputs. We also included a tagging scheme in the CTM, which keeps track of the production of carbon containing species from isoprene oxidation. It is found that the amount of tropospheric carbon monoxide (CO), formaldehyde (HCHO) and peroxyacetylnitrate (PAN) explained by isoprene oxidation ranges from 9-16%, 15-27%, and 22-32%, depending on the isoprene emissions scenario. Changes in the global tropospheric burden with different land cover inputs can reach up to 10% for CO, 15% for HCHO, and 20% for PAN. Changes for ozone are small on a global scale, but regionally differences are as large as 3DU in the tropospheric column and as large as 5 ppbv in the surface concentrations. Our results demonstrate that a careful integration of isoprene emissions and chemistry in CTMs is very important for simulating the budgets of a number of atmospheric trace gases. We further demonstrate that the model tagging scheme has the capability of improving conventional methods of constraining isoprene emissions from space-borne HCHO column observations, especially in regions where a considerable part of the variability in the HCHO column is not related to isoprene. Copyright 2008 by the American Geophysical Union.

The sensitivity of secondary organic aerosol (SOA) concentration to changes in climate and emissions is investigated using a coupled global atmosphere-land model driven by the year 2100 IPCC A1B scenario predictions. The Community Atmosphere Model (CAM3) is updated with recent laboratory determined yields for SOA formation from monoterpene oxidation, isoprene photooxidation and aromatic photooxidation. Biogenic emissions of isoprene and monoterpenes are simulated interactively using the Model of Emissions of Gases and Aerosols (MEGAN2) within the Community Land Model (CLM3). The global mean SOA burden is predicted to increase by 36% in 2100, primarily the result of rising biogenic and anthropogenic emissions which independently increase the burden by 26% and 7%. The later includes enhanced biogenic SOA formation due to increased emissions of primary organic aerosol (5-25% increases in surface SOA concentrations in 2100). Climate change alone (via temperature, removal rates, and oxidative capacity) does not change the global mean SOA production, but the global burden increases by 6%. The global burden of anthropogenic SOA experiences proportionally more growth than biogenic; SOA in 2100 from the net effect of climate and emissions (67% increase predicted). Projected anthropogenic land use change for 2100 (A2) is predicted to reduce the global SOA burden by 14%, largely the result of cropland expansion. South America is the largest global source region for SOA in the present day and 2100, but Asia experiences the largest relative growth in SOA production by 2100 because of the large predicted increases in Asian anthropogenic aromatic emissions. The projected decrease in global sulfur emissions implies that SOA will contribute a progressively larger fraction of the global aerosol burden. Copyright 2008 by the American Geophysical Union.

We analyze present-day and future carbon monoxide (CO) simulations in 26 state-of-the-art atmospheric chemistry models run to study future air quality and climate change. In comparison with near-global satellite observations from the MOPITT instrument and local surface measurements, the models show large underestimates of Northern Hemisphere (NH) extratropical CO, while typically performing reasonably well elsewhere. The results suggest that year-round emissions, probably from fossil fuel burning in east Asia and seasonal biomass burning emissions in south-central Africa, are greatly underestimated in current inventories such as IIASA and EDGAR3.2. Variability among models is large, likely resulting primarily from intermodel differences in representations and emissions of nonmethane volatile organic compounds (NMVOCs) and in hydrologic cycles, which affect OH and soluble hydrocarbon intermediates. Global mean projections of the 2030 CO response to emissions changes are quite robust. Global mean midtropospheric (500 hPa) CO increases by 12.6 ± 3.5 ppbv (16%) for the high-emissions (A2) scenario, by 1.7 ± 1.8 ppbv (2%) for the midrange (CLE) scenario, and decreases by 8.1 ± 2.3 ppbv (11%) for the low-emissions (MFR) scenario. Projected 2030 climate changes decrease global 500 hPa CO by 1.4 ± 1.4 ppbv. Local changes can be much larger. In response to climate change, substantial effects are seen in the tropics, but intermodel variability is quite large. The regional CO responses to emissions changes are robust across models, however. These range from decreases of 10–20 ppbv over much of the industrialized NH for the CLE scenario to CO increases worldwide and year-round under A2, with the largest changes over central Africa (20–30 ppbv), southern Brazil (20–35 ppbv) and south and east Asia (30–70 ppbv). The trajectory of future emissions thus has the potential to profoundly affect air quality over most of the world's populated areas.

We present a systematic comparison of tropospheric NO2 from 17 global atmospheric chemistry models with three state-of-the-art retrievals from the Global Ozone Monitoring Experiment (GOME) for the year 2000. The models used constant anthropogenic emissions from IIASA/EDGAR3.2 and monthly emissions from biomass burning based on the 1997–2002 average carbon emissions from the Global Fire Emissions Database (GFED). Model output is analyzed at 10:30 local time, close to the overpass time of the ERS-2 satellite, and collocated with the measurements to account for sampling biases due to incomplete spatiotemporal coverage of the instrument. We assessed the importance of different contributions to the sampling bias: correlations on seasonal time scale give rise to a positive bias of 30–50% in the retrieved annual means over regions dominated by emissions from biomass burning. Over the industrial regions of the eastern United States, Europe and eastern China the retrieved annual means have a negative bias with significant contributions (between –25% and +10% of the NO2 column) resulting from correlations on time scales from a day to a month. We present global maps of modeled and retrieved annual mean NO2 column densities, together with the corresponding ensemble means and standard deviations for models and retrievals. The spatial correlation between the individual models and retrievals are high, typically in the range 0.81–0.93 after smoothing the data to a common resolution. On average the models underestimate the retrievals in industrial regions, especially over eastern China and over the Highveld region of South Africa, and overestimate the retrievals in regions dominated by biomass burning during the dry season. The discrepancy over South America south of the Amazon disappears when we use the GFED emissions specific to the year 2000. The seasonal cycle is analyzed in detail for eight different continental regions. Over regions dominated by biomass burning, the timing of the seasonal cycle is generally well reproduced by the models. However, over Central Africa south of the Equator the models peak one to two months earlier than the retrievals. We further evaluate a recent proposal to reduce the NOx emission factors for savanna fires by 40% and find that this leads to an improvement of the amplitude of the seasonal cycle over the biomass burning regions of Northern and Central Africa. In these regions the models tend to underestimate the retrievals during the wet season, suggesting that the soil emissions are higher than assumed in the models. In general, the discrepancies between models and retrievals cannot be explained by a priori profile assumptions made in the retrievals, neither by diurnal variations in anthropogenic emissions, which lead to a marginal reduction of the NO2 abundance at 10:30 local time (by 2.5–4.1% over Europe). Overall, there are significant differences among the various models and, in particular, among the three retrievals. The discrepancies among the retrievals (10–50% in the annual mean over polluted regions) indicate that the previously estimated retrieval uncertainties have a large systematic component. Our findings imply that top-down estimations of NOx emissions from satellite retrievals of tropospheric NO2 are strongly dependent on the choice of model and retrieval.

Global tropospheric ozone distributions, budgets, and radiative forcings from an ensemble of 26 state-of-the-art atmospheric chemistry models have been intercompared and synthesized as part of a wider study into both the air quality and climate roles of ozone. Results from three 2030 emissions scenarios, broadly representing "optimistic," "likely," and "pessimistic" options, are compared to a base year 2000 simulation. This base case realistically represents the current global distribution of tropospheric ozone. A further set of simulations considers the influence of climate change over the same time period by forcing the central emissions scenario with a surface warming of around 0.7K. The use of a large multimodel ensemble allows us to identify key areas of uncertainty and improves the robustness of the results. Ensemble mean changes in tropospheric ozone burden between 2000 and 2030 for the 3 scenarios range from a 5% decrease, through a 6% increase, to a 15% increase. The intermodel uncertainty (±1 standard deviation) associated with these values is about ±25%. Model outliers have no significant influence on the ensemble mean results. Combining ozone and methane changes, the three scenarios produce radiative forcings of -50, 180, and 300 mW m-2, compared to a CO2 forcing over the same time period of 800-1100 mW m-2. These values indicate the importance of air pollution emissions in short- to medium-term climate forcing and the potential for stringent/lax control measures to improve/worsen future climate forcing. The model sensitivity of ozone to imposed climate change varies between models but modulates zonal mean mixing ratios by ±5 ppbv via a variety of feedback mechanisms, in particular those involving water vapor and stratosphere-troposphere exchange. This level of climate change also reduces the methane lifetime by around 4%. The ensemble mean year 2000 tropospheric ozone budget indicates chemical production, chemical destruction, dry deposition and stratospheric input fluxes of 5100, 4650, 1000 and 550 Tg(O3) yr-1, respectively. These values are significantly different to the mean budget documented by the Intergovernmental Panel on Climate Change (IPCC) Third Assessment Report (TAR). The mean ozone burden (340 Tg(O3)) is 10% larger than the IPCC TAR estimate, while the mean ozone lifetime (22 days) is 10% shorter. Results from individual models show a correlation between ozone burden and lifetime, and each model's ozone burden and lifetime respond in similar ways across the emissions scenarios. The response to climate change is much less consistent. Models show more variability in the tropics compared to midlatitudes. Some of the most uncertain areas of the models include treatments of deep tropical convection, including lightning NOx production; isoprene emissions from vegetation and isoprene's degradation chemistry; stratosphere-troposphere exchange; biomass burning; and water vapor concentrations. Copyright 2006 by the American Geophysical Union.

n this study, we present the results of nitrogen deposition on land from a set of 29 simulations from six different tropospheric chemistry models pertaining to present-day and 2100 conditions. Nitrogen deposition refers here to the deposition (wet and dry) of all nitrogen-containing gas phase chemical species resulting from NOx (NO + NO2) emissions. We show that under the assumed IPCC SRES A2 scenario the global annual average nitrogen deposition over land is expected to increase by a factor of ∼2.5, mostly because of the increase in nitrogen emissions. This will significantly expand the areas with annual average deposition exceeding 1 gN/m2/year. Using the results from all models, we have documented the strong linear relationship between models on the fraction of the nitrogen emissions that is deposited, regardless of the emissions (present day or 2100). On average, approximately 70% of the emitted nitrogen is deposited over the landmasses. For present-day conditions the results from this study suggest that the deposition over land ranges between 25 and 40 Tg(N)/year. By 2100, under the A2 scenario, the deposition over the continents is expected to range between 60 and 100 Tg(N)/year. Over forests the deposition is expected to increase from 10 Tg(N)/year to 20 Tg(N)/year. In 2100 the nitrogen deposition changes from changes in the climate account for much less than the changes from increased nitrogen emissions.