Simulations of preindustrial, present-day, and 2100 conditions in the NASA GISS composition and climate model G-PUCCINI

International audienceA model of atmospheric composition and climate has been developed at the NASA Goddard Institute for Space Studies (GISS) that includes composition seamlessly from the surface to the lower mesosphere. The model is able to capture many features of the observed magnitude, distribution, and seasonal cycle of trace species. The simulation is especially realistic in the troposphere. In the stratosphere, high latitude regions show substantial biases during period when transport governs the distribution as meridional mixing is too rapid in this model version. In other regions, including the extrapolar tropopause region that dominates radiative forcing (RF) by ozone, stratospheric gases are generally well-simulated. The model's stratosphere-troposphere exchange (STE) agrees well with values inferred from observations for both the global mean flux and the ratio of Northern (NH) to Southern Hemisphere (SH) downward fluxes. Simulations of preindustrial (PI) to present-day (PD) changes show tropospheric ozone burden increases of 11% while the stratospheric burden decreases by 18%. The resulting tropopause RF values are ?0.06 W/m2 from stratospheric ozone and 0.40 W/m2 from tropospheric ozone. Global mean mass-weighted OH decreases by 16% from the PI to the PD. STE of ozone also decreased substantially during this time, by 14%. Comparison of the PD with a simulation using 1979 pre-ozone hole conditions for the stratosphere shows a much larger downward flux of ozone into the troposphere in 1979, resulting in a substantially greater tropospheric ozone burden than that seen in the PD run. This implies that reduced STE due to stratospheric ozone depletion may have offset as much as 2/3 of the tropospheric ozone burden increase from PI to PD. However, the model overestimates the downward flux of ozone at high Southern latitudes, so this estimate is likely an upper limit. In the future, the tropospheric ozone burden increases by 101% in 2100 for the A2 scenario including both emissions and climate changes. The primary reason is enhanced STE, which increases by 124% (168% in the SH extratropics, and 114% in the NH extratropics). Climate plays a minimal role in the SH increases, but contributes 38% in the NH. Chemistry and dry deposition both change so as to reduce tropospheric ozone, partially in compensation for the enhanced STE, but the increased ozone influx dominates the burden changes. The net RF due to projected ozone changes is 0.8 W/m2 for A2. The influence of climate change alone is ?0.2 W/m2, making it a substantial contributor to the net RF. The tropospheric oxidation capacity increases seven percent in the full A2 simulation, and 36% due to A2 climate change alone

Ozone database in support of CMIP5 simulations: results and corresponding radiative forcing

A continuous tropospheric and stratospheric vertically resolved ozone time series, from 1850 to 2099, has been generated to be used as forcing in global climate models that do not include interactive chemistry. A multiple linear regression analysis of SAGE I+II satellite observations and polar ozonesonde measurements is used for the stratospheric zonal mean dataset during the well-observed period from 1979 to 2009. In addition to terms describing the mean annual cycle, the regression includes terms representing equivalent effective stratospheric chlorine (EESC) and the 11-yr solar cycle variability. The EESC regression fit coefficients, together with pre-1979 EESC values, are used to extrapolate the stratospheric ozone time series backward to 1850. While a similar procedure could be used to extrapolate into the future, coupled chemistry climate model (CCM) simulations indicate that future stratospheric ozone abundances are likely to be significantly affected by climate change, and capturing such effects through a regression model approach is not feasible. Therefore, the stratospheric ozone dataset is extended into the future (merged in 2009) with multimodel mean projections from 13 CCMs that performed a simulation until 2099 under the SRES (Special Report on Emission Scenarios) A1B greenhouse gas scenario and the A1 adjusted halogen scenario in the second round of the Chemistry-Climate Model Validation (CCMVal-2) Activity. The stratospheric zonal mean ozone time series is merged with a three-dimensional tropospheric data set extracted from simulations of the past by two CCMs (CAM3.5 and GISSPUCCINI)and of the future by one CCM (CAM3.5). The future tropospheric ozone time series continues the historical CAM3.5 simulation until 2099 following the four different Representative Concentration Pathways (RCPs). Generally good agreement is found between the historical segment of the ozone database and satellite observations, although it should be noted that total column ozone is overestimated in the southern polar latitudes during spring and tropospheric column ozone is slightly underestimated. Vertical profiles of tropospheric ozone are broadly consistent with ozonesondes and in-situ measurements, with some deviations in regions of biomass burning. The tropospheric ozone radiative forcing (RF) from the 1850s to the 2000s is 0.23Wm−2, lower than previous results. The lower value is mainly due to (i) a smaller increase in biomass burning emissions; (ii) a larger influence of stratospheric ozone depletion on upper tropospheric ozone at high southern latitudes; and possibly (iii) a larger influence of clouds (which act to reduce the net forcing) compared to previous radiative forcing calculations. Over the same period, decreases in stratospheric ozone, mainly at high latitudes, produce a RF of −0.08Wm−2, which is more negative than the central Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment Report (AR4) value of −0.05Wm−2, but which is within the stated range of −0.15 to +0.05Wm−2. The more negative value is explained by the fact that the regression model simulates significant ozone depletion prior to 1979, in line with the increase in EESC and as confirmed by CCMs, while the AR4 assumed no change in stratospheric RF prior to 1979. A negative RF of similar magnitude persists into the future, although its location shifts from high latitudes to the tropics. This shift is due to increases in polar stratospheric ozone, but decreases in tropical lower stratospheric ozone, related to a strengthening of the Brewer-Dobson circulation, particularly through the latter half of the 21st century. Differences in trends in tropospheric ozone among the four RCPs are mainly driven by different methane concentrations, resulting in a range of tropospheric ozone RFs between 0.4 and 0.1Wm−2 by 2100. The ozone dataset described here has been released for the Coupled Model Intercomparison Project (CMIP5) model simulations in netCDF Climate and Forecast (CF) Metadata Convention at the PCMDI website (http://cmip-pcmdi.llnl.gov/)

Climate forcing and air quality change due to regional emissions reductions by economic sector

We examine the air quality (AQ) and radiative forcing (RF) response to emissions reductions by economic sector for North America and developing Asia in the CAM and GISS composition/climate models. Decreases in annual average surface particulate are relatively robust, with intermodel variations in magnitude typically <30% and very similar spatial structure. Surface ozone responses are small and highly model dependent. The largest net RF results from reductions in emissions from the North America industrial/power and developing Asia domestic fuel burning sectors. Sulfate reductions dominate the first case, for which intermodel variations in the sulfate (or total) aerosol optical depth (AOD) responses are ~30% and the modeled spatial patterns of the AOD reductions are highly correlated (R=0.9). Decreases in BC dominate the developing Asia domestic fuel burning case, and show substantially greater model-to-model differences. Intermodel variations in tropospheric ozone burden changes are also large, though aerosol changes dominate those cases with substantial net climate forcing. The results indicate that across-the-board emissions reductions in domestic fuel burning in developing Asia and in surface transportation in North America are likely to offer the greatest potential for substantial, simultaneous improvement in local air quality and near-term mitigation of global climate change via short-lived species. Conversely, reductions in industrial/power emissions have the potential to accelerate near-term warming, though they would improve AQ and have a long-term cooling effect on climate. These broad conclusions appear robust to intermodel differences

The AeroCom evaluation and intercomparison of organic aerosol in global models

This paper evaluates the current status of global modeling of the organic aerosol (OA) in the troposphere and analyzes the differences between models as well as between models and observations. Thirty-one global chemistry transport models (CTMs) and general circulation models (GCMs) have participated in this intercomparison, in the framework of AeroCom phase II. The simulation of OA varies greatly between models in terms of the magnitude of primary emissions, secondary OA (SOA) formation, the number of OA species used (2 to 62), the complexity of OA parameterizations (gas-particle partitioning, chemical aging, multiphase chemistry, aerosol microphysics), and the OA physical, chemical and optical properties. The diversity of the global OA simulation results has increased since earlier AeroCom experiments, mainly due to the increasing complexity of the SOA parameterization in models, and the implementation of new, highly uncertain, OA sources. Diversity of over one order of magnitude exists in the modeled vertical distribution of OA concentrations that deserves a dedicated future study. Furthermore, although the OA / OC ratio depends on OA sources and atmospheric processing, and is important for model evaluation against OA and OC observations, it is resolved only by a few global models. The median global primary OA (POA) source strength is 56 Tg a−1 (range 34–144 Tg a−1) and the median SOA source strength (natural and anthropogenic) is 19 Tg a−1 (range 13–121 Tg a−1). Among the models that take into account the semi-volatile SOA nature, the median source is calculated to be 51 Tg a−1 (range 16–121 Tg a−1), much larger than the median value of the models that calculate SOA in a more simplistic way (19 Tg a−1; range 13–20 Tg a−1, with one model at 37 Tg a−1). The median atmospheric burden of OA is 1.4 Tg (24 models in the range of 0.6–2.0 Tg and 4 between 2.0 and 3.8 Tg), with a median OA lifetime of 5.4 days (range 3.8–9.6 days). In models that reported both OA and sulfate burdens, the median value of the OA/sulfate burden ratio is calculated to be 0.77; 13 models calculate a ratio lower than 1, and 9 models higher than 1. For 26 models that reported OA deposition fluxes, the median wet removal is 70 Tg a−1 (range 28–209 Tg a−1), which is on average 85% of the total OA deposition. Fine aerosol organic carbon (OC) and OA observations from continuous monitoring networks and individual field campaigns have been used for model evaluation. At urban locations, the model–observation comparison indicates missing knowledge on anthropogenic OA sources, both strength and seasonality. The combined model–measurements analysis suggests the existence of increased OA levels during summer due to biogenic SOA formation over large areas of the USA that can be of the same order of magnitude as the POA, even at urban locations, and contribute to the measured urban seasonal pattern. Global models are able to simulate the high secondary character of OA observed in the atmosphere as a result of SOA formation and POA aging, although the amount of OA present in the atmosphere remains largely underestimated, with a mean normalized bias (MNB) equal to −0.62 (−0.51) based on the comparison against OC (OA) urban data of all models at the surface, −0.15 (+0.51) when compared with remote measurements, and −0.30 for marine locations with OC data. The mean temporal correlations across all stations are low when compared with OC (OA) measurements: 0.47 (0.52) for urban stations, 0.39 (0.37) for remote stations, and 0.25 for marine stations with OC data. The combination of high (negative) MNB and higher correlation at urban stations when compared with the low MNB and lower correlation at remote sites suggests that knowledge about the processes that govern aerosol processing, transport and removal, on top of their sources, is important at the remote stations. There is no clear change in model skill with increasing model complexity with regard to OC or OA mass concentration. However, the complexity is needed in models in order to distinguish between anthropogenic and natural OA as needed for climate mitigation, and to calculate the impact of OA on climate accurately

Multimodel simulations of carbon monoxide: comparison with observations and projected near-future changes

We analyze present-day and future carbon monoxide (CO) simulations in 26 state-ofthe- 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