We present the global general circulation model IPSL-CM5 developed to study the long-term response of the climate system to natural and anthropogenic forcings as part of the 5th Phase of the Coupled Model Intercomparison Project (CMIP5). This model includes an interactive carbon cycle, a representation of tropospheric and stratospheric chemistry, and a comprehensive representation of aerosols. As it represents the principal dynamical, physical, and biogeochemical processes relevant to the climate system, it may be referred to as an Earth System Model. However, the IPSL-CM5 model may be used in a multitude of configurations associated with different boundary conditions and with a range of complexities in terms of processes and interactions. This paper presents an overview of the different model components and explains how they were coupled and used to simulate historical climate changes over the past 150 years and different scenarios of future climate change. A single version of the IPSL-CM5 model (IPSL-CM5A-LR) was used to provide climate projections associated with different socio-economic scenarios, including the different Representative Concentration Pathways considered by CMIP5 and several scenarios from the Special Report on Emission Scenarios considered by CMIP3. Results suggest that the magnitude of global warming projections primarily depends on the socio-economic scenario considered, that there is potential for an aggressive mitigation policy to limit global warming to about two degrees, and This paper is a contribution to the special issue on the IPSL and CNRM global climate and Earth System Models, both developed in France and contributing to the 5th coupled model intercomparison project.

This study presents the global climate model IPSL-CM6A-LR developed at Institut Pierre-Simon Laplace (IPSL) to study natural climate variability and climate response to natural and anthropogenic forcings as part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). This article describes the different model components, their coupling, and the simulated climate in comparison to previous model versions. We focus here on the representation of the physical climate along with the main characteristics of the global carbon cycle. The model's climatology, as assessed from a range of metrics (related in particular to radiation, temperature, precipitation, and wind), is strongly improved in comparison to previous model versions. Although they are reduced, a number of known biases and shortcomings (e.g., double Intertropical Convergence Zone [ITCZ], frequency of midlatitude wintertime blockings, and El Niño-Southern Oscillation [ENSO] dynamics) persist. The equilibrium climate sensitivity and transient climate response have both increased from the previous climate model IPSL-CM5A-LR used in CMIP5. A large ensemble of more than 30 members for the historical period (1850-2018) and a smaller ensemble for a range of emissions scenarios (until 2100 and 2300) are also presented and discussed.Plain Language Summary Climate models are unique tools to investigate the characteristics and behavior of the climate system. While climate models and their components are developed gradually over the years, the sixth phase of the Coupled Model Intercomparison Project (CMIP6) has been the

Black carbon (BC) is increasingly recognized as a significant air pollutant with harmful effects on human health, either in its own right or as a carrier of other chemicals. The adverse impact is of particular concern in those developing regions with high emissions and a growing population density. The results of recent studies indicate that BC emissions could be underestimated by a factor of 2-3 and this is particularly true for the hot-spot Asian region. Here we present a unique inventory at 10-km resolution based on a recently published global fuel consumption data product and updated emission factor measurements. The unique inventory is coupled to an Asia-nested (∼50 km) atmospheric model and used to calculate the global population exposure to BC with fully quantified uncertainty. Evaluating the modeled surface BC concentrations against observations reveals great improvement. The bias is reduced from −88% to −35% in Asia when the unique inventory and higher-resolution model replace a previous inventory combined with a coarse-resolution model. The bias can be further reduced to −12% by downscaling to 10 km using emission as a proxy. Our estimated global population-weighted BC exposure concentration constrained by observations is 2.14 μg·m −3 ; 130% higher than that obtained using less detailed inventories and low-resolution models.air pollution | climate change | model resolution | emission inventory

[1] D 14 CO 2 observations are increasingly used to constrain recently added fossil fuel CO 2 in the atmosphere. We use the LMDZ global atmospheric transport model to examine the pseudo-Lagrangian framework commonly used to determine recently added fossil fuel CO 2 (CO 2ff ). Our results confirm that D 14 CO 2 spatial variability in the Northern Hemisphere troposphere is dominated by the effect of CO 2ff , whereas in the Southern Hemisphere, ocean CO 2 exchange is more important. The model indicates that the free troposphere, at 3-5 km altitude, is a good choice for ''background,'' relative to which the recently added fossil fuel CO 2 can be calculated, although spatial variability in free tropospheric D 14 CO 2 contributes additional uncertainty to the CO 2ff calculation.Comparison of model and observations suggests that care must be taken in using high-altitude mountain sites as a proxy for free tropospheric air, since these sites may be occasionally influenced by (polluted) boundary layer air, especially in summer. Other sources of CO 2 which have D 14C different than that of the atmosphere contribute a bias, which, over the Northern Hemisphere land, is mostly due to the terrestrial biosphere, whereas ocean CO 2 exchange and nuclear industry and natural cosmogenic production of 14 C contribute only weakly. The model indicates that neglecting this bias leads to a consistent underestimation of CO 2ff , typically between 0.2 and 0.5ppm of CO 2 , with a maximum in summer. While our analysis focuses on fossil fuel CO 2 , our conclusions, particularly the choice of background site, can also be applied to other trace gases emitted at the surface.Citation: Turnbull, J., P. Rayner, J. Miller, T. Naegler, P. Ciais, and A. Cozic (2009), On the use of 14 CO 2 as a tracer for fossil fuel CO 2 : Quantifying uncertainties using an atmospheric transport model,

Global aerosol and ozone distributions and their associated radiative forcings were simulated between 1850 and 2100 following a recent historical emission dataset and under the representative concentration pathways (RCP) for the future. These simulations were used in an Earth System Model to account for the changes in both radiatively and chemically active compounds, when simulating the climate evolution. The past negative stratospheric ozone trends result in a negative climate forcing culminating at -0.15 W m -2 in the 1990s. In the meantime, the tropospheric ozone burden increase generates a positive climate forcing peaking at 0.41 W m -2 . The future evolution of ozone strongly depends on the RCP scenario considered. In RCP4.5 and RCP6.0, the evolution of both stratospheric and tropospheric ozone generate relatively weak radiative forcing changes until 2060-2070 followed by a relative 30 % decrease in radiative forcing by 2100. In contrast, RCP8.5 and RCP2.6 model projections exhibit strongly different ozone radiative forcing trajectories. In the RCP2.6 scenario, both effects (stratospheric ozone, a negative forcing, and tropospheric ozone, a positive forcing) decline towards 1950s values while they both get stronger in the RCP8.5 scenario. Over the twentieth century, the evolution of the total aerosol burden is characterized by a strong increase after World War II until the middle of the 1980s followed by a stabilization during the last decade due to the strong decrease in sulfates in OECD countries since the 1970s. The cooling effects reach their maximal values in 1980, with -0.34 and -0.28 W m -2 respectively for direct and indirect total radiative forcings. According to the RCP scenarios, the aerosol content, after peaking around 2010, is projected to decline strongly and monotonically during the twenty-first century for the RCP8.5, 4.5 and 2.6 scenarios. While for RCP6.0 the decline occurs later, after peaking around 2050. As a consequence the relative importance of the total cooling effect of aerosols becomes weaker throughout the twenty-first century compared with the positive forcing of greenhouse gases. Nevertheless, both surface ozone and aerosol content show very different regional features depending on the future scenario considered. Hence, in 2050, surface ozone changes vary between -12 and ?12 ppbv over Asia depending on the RCP projection, whereas the regional direct aerosol radiative forcing can locally exceed -3 W m -2 . This paper is a contribution to the special issue on the IPSL and CNRM global climate and Earth System Models, both developed in France and contributing to the 5th coupled model intercomparison project.Electronic supplementary material The online version of this article

Abstract. We applied a climate-chemistry global model to evaluate the impact of black carbon (BC) deposition on the Himalayan snow cover from 1998 to 2008. Using a stretched grid with a resolution of 50 km over this complex topography, the model reproduces reasonably well the remotely sensed observations of the snow cover duration. Similar to observations, modelled atmospheric BC concentrations in the central Himalayas reach a minimum during the monsoon and a maximum during the post-and pre-monsoon periods. Comparing the simulated BC concentrations in the snow with observations is more challenging because of their high spatial variability and complex vertical distribution. We simulated spring BC concentrations in surface snow varying from tens to hundreds of µg kg −1 , higher by one to two orders of magnitude than those observed in ice cores extracted from central Himalayan glaciers at high elevations (> 6000 m a.s.l.), but typical for seasonal snow cover sampled in middle elevation regions (< 6000 m a.s.l.). In these areas, we estimate that both wet and dry BC depositions affect the Himalayan snow cover reducing its annual duration by 1 to 8 days. In our simulations, the effect of anthropogenic BC deposition on snow is quite low over the Tibetan Plateau because this area is only sparsely snow covered. However, the impact becomes larger along the entire Hindu-Kush, Karakorum and Himalayan mountain ranges. In these regions, BC in snow induces an increase of the net short-wave radiation at the surface with an annual mean of 1 to 3 W m −2 leading to a localised warming between 0.05 and 0.3 • C.

International audienceA new version of the biogenic volatile organic compounds (BVOCs) emission scheme has been developed in the global vegetation model ORCHIDEE (Organizing Carbon and Hydrology in Dynamic EcosystEm), which includes an extended list of biogenic emitted compounds, updated emission factors (EFs), a dependency on light for almost all compounds and a multi-layer radiation scheme. Over the 2000–2009 period, using this model, we estimate mean global emissions of 465 Tg C yr⁻¹ for isoprene, 107.5 Tg C yr⁻¹ for monoterpenes, 38 Tg C yr⁻¹ for methanol, 25 Tg C yr⁻¹ for acetone and 24 Tg C yr⁻¹ for sesquiterpenes. The model results are compared to state-of-the-art emission budgets, showing that the ORCHIDEE emissions are within the range of published estimates. ORCHIDEE BVOC emissions are compared to the estimates of the Model of Emissions of Gases and Aerosols from Nature (MEGAN), which is largely used throughout the biogenic emissions and atmospheric chemistry community. Our results show that global emission budgets of the two models are, in general, in good agreement. ORCHIDEE emissions are 8 % higher for isoprene, 8 % lower for methanol, 17 % higher for acetone, 18 % higher for monoterpenes and 39 % higher for sesquiterpenes, compared to the MEGAN estimates. At the regional scale, the largest differences between ORCHIDEE and MEGAN are highlighted for isoprene in northern temperate regions, where ORCHIDEE emissions are higher by 21 Tg C yr⁻¹, and for monoterpenes, where they are higher by 4.4 and 10.2 Tg C yr⁻¹ in northern and southern tropical regions compared to MEGAN. The geographical differences between the two models are mainly associated with different EF and plant functional type (PFT) distributions, while differences in the seasonal cycle are mostly driven by differences in the leaf area index (LAI). Sensitivity tests are carried out for both models to explore the response to key variables or parameters such as LAI and light-dependent fraction (LDF). The ORCHIDEE and MEGAN emissions are differently affected by LAI changes, with a response highly depending on the compound considered. Scaling the LAI by a factor of 0.5 and 1.5 changes the isoprene global emission by −21 and +8 % for ORCHIDEE and −15 and +7 % for MEGAN, and affects the global emissions of monoterpenes by −43 and +40 % for ORCHIDEE and −11 and +3 % for MEGAN. Performing a further sensitivity test, forcing ORCHIDEE with the MODIS LAI, confirms the high sensitivity of the ORCHIDEE emission module to LAI variation. We find that MEGAN is more sensitive to variation in the LDF parameter than ORCHIDEE. Our results highlight the importance and the need to further explore the BVOC emission estimate variability and the potential for using models to investigate the estimated uncertainties

The large peatland carbon stocks in the land use change‐affected areas of equatorial Asia are vulnerable to fire. Combining satellite observations of active fire, burned area, and atmospheric concentrations of combustion tracers with a Bayesian inversion, we estimated the amount and variability of fire carbon emissions in equatorial Asia over the period 1997–2015. Emissions in 2015 were of 0.51 ± 0.17 Pg carbon—less than half of the emissions from the previous 1997 extreme El Niño, explained by a less acute water deficit. Fire severity could be empirically hindcasted from the cumulative water deficit with a lead time of 1 to 2 months. Based on CMIP5 climate projections and an exponential empirical relationship found between fire carbon emissions and water deficit, we infer a total fire carbon loss ranging from 12 to 25 Pg by 2100 which is a significant positive feedback to climate warming.