Inventories for global aerosol and aerosol precursor emissions have been collected (based on published inventories and published simulations), assessed and prepared for the year 2000 (present-day conditions) and for the year 1750 (pre-industrial conditions). These global datasets establish a comprehensive source for emission input to global modeling, when simulating the aerosol impact on climate with state-of-the-art aerosol component modules. As these modules stratify aerosol into dust, sea-salt, sulfate, organic matter and soot, for all these aerosol types global fields on emission strength and recommendations for injection altitude and particulate size are provided. Temporal resolution varies between daily (dust and sea-salt), monthly (wild-land fires) and annual (all other emissions). These datasets benchmark aerosol emissions according to the knowledge in the year 2004. They are intended to serve as systematic constraints in sensitivity studies of the AeroCom initiative, which seeks to quantify (actual) uncertainties in aerosol global modeling.

This article describes the new Earth system model (ESM), the Model for Interdisciplinary Research on Climate, Earth System version 2 for Long-term simulations (MIROC-ES2L), using a state-of-the-art climate model as the physical core. This model embeds a terrestrial biogeochemical component with explicit carbon–nitrogen interaction to account for soil nutrient control on plant growth and the land carbon sink. The model's ocean biogeochemical component is largely updated to simulate the biogeochemical cycles of carbon, nitrogen, phosphorus, iron, and oxygen such that oceanic primary productivity can be controlled by multiple nutrient limitations. The ocean nitrogen cycle is coupled with the land component via river discharge processes, and external inputs of iron from pyrogenic and lithogenic sources are considered. Comparison of a historical simulation with observation studies showed that the model could reproduce the transient global climate change and carbon cycle as well as the observed large-scale spatial patterns of the land carbon cycle and upper-ocean biogeochemistry. The model demonstrated historical human perturbation of the nitrogen cycle through land use and agriculture and simulated the resultant impact on the terrestrial carbon cycle. Sensitivity analyses under preindustrial conditions revealed that the simulated ocean biogeochemistry could be altered regionally (and substantially) by nutrient input from the atmosphere and rivers. Based on an idealized experiment in which CO₂ was prescribed to increase at a rate of 1 % yr⁻¹, the transient climate response (TCR) is estimated to be 1.5 K, i.e., approximately 70 % of that from our previous ESM used in the Coupled Model Intercomparison Project Phase 5 (CMIP5). The cumulative airborne fraction (AF) is also reduced by 15 % because of the intensified land carbon sink, which results in an airborne fraction close to the multimodel mean of the CMIP5 ESMs. The transient climate response to cumulative carbon emissions (TCRE) is 1.3 K EgC⁻¹, i.e., slightly smaller than the average of the CMIP5 ESMs, which suggests that “optimistic” future climate projections will be made by the model. This model and the simulation results contribute to CMIP6. The MIROC-ES2L could further improve our understanding of climate–biogeochemical interaction mechanisms, projections of future environmental changes, and exploration of our future options regarding sustainable development by evolving the processes of climate, biogeochemistry, and human activities in a holistic and interactive manner.

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, the relative contributions of the world's major source regions to the global dust cycle remain poorly constrained. This problem hinders accounting for the potentially large impact of regional differences in dust properties on clouds, the Earth's energy balance, and terrestrial and marine biogeochemical cycles. Here, we constrain the contribution of each of the world's main dust source regions to the global dust cycle. We use an analytical framework that integrates an ensemble of global aerosol model simulations with observationally informed constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth (DAOD). We obtain a dataset that constrains the relative contribution of nine major source regions to size-resolved dust emission, atmospheric loading, DAOD, concentration, and deposition flux. We find that the 22–29 Tg (1 standard error range) global loading of dust with a geometric diameter up to 20 µm is partitioned as follows: North African source regions contribute ∼ 50 % (11–15 Tg), Asian source regions contribute ∼ 40 % (8–13 Tg), and North American and Southern Hemisphere regions contribute ∼ 10 % (1.8–3.2 Tg). These results suggest that current models on average overestimate the contribution of North African sources to atmospheric dust loading at ∼ 65 %, while underestimating the contribution of Asian dust at ∼ 30 %. Our results further show that each source region's dust loading peaks in local spring and summer, which is partially driven by increased dust lifetime in those seasons. We also quantify the dust deposition flux to the Amazon rainforest to be ∼ 10 Tg yr⁻¹, which is a factor of 2–3 less than inferred from satellite data by previous work that likely overestimated dust deposition by underestimating the dust mass extinction efficiency. The data obtained in this paper can be used to obtain improved constraints on dust impacts on clouds, climate, biogeochemical cycles, and other parts of the Earth system.

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, atmospheric models struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused by fundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of 2 relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global basis, we find that the emission flux of dust with a geometric diameter up to 20 µm (PM₂₀) is approximately 5000 Tg yr⁻¹, which is greater than most models account for. This larger PM₂₀ dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded datasets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this dataset is more accurate than current model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.