ALBERT

Description: Interannual variability in desert dust is widely observed and simulated, yet the sensitivity of these desert dust simulations to a particular meteorological dataset, as well as a particular model construction, is not well known. Here we use version 4 of the Community Atmospheric Model (CAM4) with the Community Earth System Model (CESM) to simulate dust forced by three different reanalysis meteorological datasets for the period 1990–2005. We then contrast the results of these simulations with dust simulated using online winds dynamically generated from sea surface temperatures, as well as with simulations conducted using other modeling frameworks but the same meteorological forcings, in order to determine the sensitivity of climate model output to the specific reanalysis dataset used. For the seven cases considered in our study, the different model configurations are able to simulate the annual mean of the global dust cycle, seasonality and interannual variability approximately equally well (or poorly) at the limited observational sites available. Overall, aerosol dust-source strength has remained fairly constant during the time period from 1990 to 2005, although there is strong seasonal and some interannual variability simulated in the models and seen in the observations over this time period. Model interannual variability comparisons to observations, as well as comparisons between models, suggest that interannual variability in dust is still difficult to simulate accurately, with averaged correlation coefficients of 0.1 to 0.6. Because of the large variability, at least 1 year of observations at most sites are needed to correctly observe the mean, but in some regions, particularly the remote oceans of the Southern Hemisphere, where interannual variability may be larger than in the Northern Hemisphere, 2–3 years of data are likely to be needed.

Description: Atmospheric processing of iron in dust and combustion aerosols is simulated using an intermediate-complexity soluble iron mechanism designed for Earth system models. The solubilization mechanism includes both a dependence on aerosol water pH and in-cloud oxalic acid. The simulations of size-resolved total, soluble and fractional iron solubility indicate that this mechanism captures many but not all of the features seen from cruise observations of labile iron. The primary objective was to determine the extent to which our solubility scheme could adequately match observations of fractional iron solubility. We define a semi-quantitative metric as the model mean at points with observations divided by the observational mean (MMO). The model is in reasonable agreement with observations of fractional iron solubility with an MMO of 0.86. Several sensitivity studies are performed to ascertain the degree of complexity needed to match observations; including the oxalic acid enhancement is necessary, while different parameterizations for calculating model oxalate concentrations are less important. The percent change in soluble iron deposition between the reference case (REF) and the simulation with acidic processing alone is 63.8 %, which is consistent with previous studies. Upon deposition to global oceans, global mean combustion iron solubility to total fractional iron solubility is 8.2 %; however, the contribution of fractional iron solubility from combustion sources to ocean basins below 15∘ S is approximately 50 %. We conclude that, in many remote ocean regions, sources of iron from combustion and dust aerosols are equally important. Our estimates of changes in deposition of soluble iron to the ocean since preindustrial climate conditions suggest roughly a doubling due to a combination of higher dust and combustion iron emissions along with more efficient atmospheric processing.

Description: Two interglacial epochs are included in the suite of Paleoclimate Modeling Intercomparison Project (PMIP4) simulations in the Coupled Model Intercomparison Project (CMIP6). The experimental protocols for simulations of the mid-Holocene (midHolocene, 6000 years before present) and the Last Interglacial (lig127k, 127 000 years before present) are described here. These equilibrium simulations are designed to examine the impact of changes in orbital forcing at times when atmospheric greenhouse gas levels were similar to those of the preindustrial period and the continental configurations were almost identical to modern ones. These simulations test our understanding of the interplay between radiative forcing and atmospheric circulation, and the connections among large-scale and regional climate changes giving rise to phenomena such as land–sea contrast and high-latitude amplification in temperature changes, and responses of the monsoons, as compared to today. They also provide an opportunity, through carefully designed additional sensitivity experiments, to quantify the strength of atmosphere, ocean, cryosphere, and land-surface feedbacks. Sensitivity experiments are proposed to investigate the role of freshwater forcing in triggering abrupt climate changes within interglacial epochs. These feedback experiments naturally lead to a focus on climate evolution during interglacial periods, which will be examined through transient experiments. Analyses of the sensitivity simulations will also focus on interactions between extratropical and tropical circulation, and the relationship between changes in mean climate state and climate variability on annual to multi-decadal timescales. The comparative abundance of paleoenvironmental data and of quantitative climate reconstructions for the Holocene and Last Interglacial make these two epochs ideal candidates for systematic evaluation of model performance, and such comparisons will shed new light on the importance of external feedbacks (e.g., vegetation, dust) and the ability of state-of-the-art models to simulate climate changes realistically.

Description: This work reports on the current status of the global modeling of iron (Fe) deposition fluxes and atmospheric concentrations and the analyses of the differences between models, as well as between models and observations. A total of four global 3-D chemistry transport (CTMs) and general circulation (GCMs) models participated in this intercomparison, in the framework of the United Nations Joint Group of Experts on the Scientific Aspects of Marine Environmental Protection (GESAMP) Working Group 38, “The Atmospheric Input of Chemicals to the Ocean”. The global total Fe (TFe) emission strength in the models is equal to ∼72 Tg Fe yr−1 (38–134 Tg Fe yr−1) from mineral dust sources and around 2.1 Tg Fe yr−1 (1.8–2.7 Tg Fe yr−1) from combustion processes (the sum of anthropogenic combustion/biomass burning and wildfires). The mean global labile Fe (LFe) source strength in the models, considering both the primary emissions and the atmospheric processing, is calculated to be 0.7 (±0.3) Tg Fe yr−1, accounting for both mineral dust and combustion aerosols. The mean global deposition fluxes into the global ocean are estimated to be in the range of 10–30 and 0.2–0.4 Tg Fe yr−1 for TFe and LFe, respectively, which roughly corresponds to a respective 15 and 0.3 Tg Fe yr−1 for the multi-model ensemble model mean. The model intercomparison analysis indicates that the representation of the atmospheric Fe cycle varies among models, in terms of both the magnitude of natural and combustion Fe emissions as well as the complexity of atmospheric processing parameterizations of Fe-containing aerosols. The model comparison with aerosol Fe observations over oceanic regions indicates that most models overestimate surface level TFe mass concentrations near dust source regions and tend to underestimate the low concentrations observed in remote ocean regions. All models are able to simulate the tendency of higher Fe concentrations near and downwind from the dust source regions, with the mean normalized bias for the Northern Hemisphere (∼14), larger than that of the Southern Hemisphere (∼2.4) for the ensemble model mean. This model intercomparison and model–observation comparison study reveals two critical issues in LFe simulations that require further exploration: (1) the Fe-containing aerosol size distribution and (2) the relative contribution of dust and combustion sources of Fe to labile Fe in atmospheric aerosols over the remote oceanic regions.

Description: The area of leaves in the plant canopy, measured as leaf area index (LAI), modulates key land–atmosphere interactions, including the exchange of energy, moisture, carbon dioxide (CO2), and other trace gases and aerosols, and is therefore an essential variable in predicting terrestrial carbon, water, and energy fluxes. Here our goal is to characterize the LAI projections from the latest generation of earth system models (ESMs) for the Representative Concentration Pathway (RCP) 8.5 and RCP4.5 scenarios. On average, the models project increases in LAI in both RCP8.5 and RCP4.5 over most of the globe, but also show decreases in some parts of the tropics. Because of projected increases in variability, there are also more frequent periods of low LAI across broad regions of the tropics. Projections of LAI changes varied greatly among models: some models project very modest changes, while others project large changes, usually increases. Modeled LAI typically increases with modeled warming in the high latitudes, but often decreases with increasing local warming in the tropics. The models with the most skill in simulating current LAI in the tropics relative to satellite observations tend to project smaller increases in LAI in the tropics in the future compared to the average of all the models. Using LAI projections to identify regions that may be vulnerable to climate change presents a slightly different picture than using precipitation projections, suggesting LAI may be an additional useful tool for understanding climate change impacts. Going forward, users of LAI projections from the CMIP5 ESMs evaluated here should be aware that model outputs do not exhibit clear-cut relationships to vegetation carbon and precipitation. Our findings underscore the need for more attention to LAI projections, in terms of understanding the drivers of projected changes and improvements to model skill.

Description: Two interglacial epochs are included in the suite of paleoclimate simulations in the present phase of the Coupled Model Intercomparison Project (CMIP6). Equilibrium simulations of the mid-Holocene (midHolocene, 6000 years before present) and the Last Interglacial (lig127k, 127,000 years before present) are designed to examine the impact of changes in orbital forcing at times when atmospheric greenhouse gas levels were similar to those of the preindustrial period and the continental configurations were almost identical to modern. These simulations test our understanding of the interplay between radiative forcing and atmospheric circulation, and the connections among large-scale and regional climate changes giving rise to phenomena such as land-sea contrast and high-latitude amplification in temperature changes, and responses of the monsoons, as compared to today. They also provide an opportunity, through carefully designed additional sensitivity experiments as part of the Paleoclimate Modeling Intercomparison Project (PMIP4), to quantify the strength of atmosphere, ocean, cryosphere, and land-surface feedbacks. Sensitivity experiments are proposed to investigate the role of freshwater forcing in triggering abrupt climate changes within interglacial epochs. These feedback experiments naturally lead to a focus on climate evolution during interglacial periods, which will be examined through transient experiments. Analyses of the sensitivity simulations will also focus on interactions between extratropical and tropical circulation, and the relationship between changes in mean climate state and climate variability on annual to multi-decadal timescales. The comparative abundance of paleoenvironmental data and of quantitative climate reconstructions for the Holocene and Last Interglacial make these two epochs ideal candidates for systematic evaluation of model performance, and such comparisons will shed new light on the importance of external feedbacks (e.g., vegetation, dust) and the ability of state-of-the-art models to simulate climate changes realistically.

Description: The Last Glacial Maximum (LGM, 21,000 years ago) is one of the suite of paleoclimate simulations included in the current phase of the Coupled Model Intercomparison Project (CMIP6). It is an interval when insolation was similar to present, but global ice volume was at a maximum, eustatic sea level was at or close to a minimum, greenhouse gas concentrations were lower, atmospheric aerosol loadings were higher than today, and vegetation and land-surface characteristics were different from today. The LGM has been a focus for the Paleoclimate Modelling Intercomparison Project (PMIP) since its inception, and thus many of the problems that might be associated with simulating such a radically different climate are well documented. The LGM state provides an ideal case study for evaluating climate model performance because the changes in forcing and temperature between the LGM and pre-industrial are of the same order of magnitude as those projected for the end of the 21st century. Thus, the CMIP6 LGM experiment could provide additional information that can be used to constrain estimates of climate sensitivity. The design of the Tier 1 LGM experiment (lgm) includes an assessment of uncertainties in boundary conditions, in particular through the use of different reconstructions of the ice sheets and of the change in dust forcing. Additional sensitivity experiments have been designed to quantify feedbacks associated with land-surface changes and aerosol loadings, and to isolate the role of individual forcings. Model analysis and evaluation will capitalise on the relative abundance of palaeoenvironmental observations and quantitative climate reconstructions already available for the LGM.

Description: The Last Glacial Maximum (LGM, 21 000 years ago) is one of the suite of paleoclimate simulations included in the current phase of the Coupled Model Intercomparison Project (CMIP6). It is an interval when insolation was similar to the present, but global ice volume was at a maximum, eustatic sea level was at or close to a minimum, greenhouse gas concentrations were lower, atmospheric aerosol loadings were higher than today, and vegetation and land-surface characteristics were different from today. The LGM has been a focus for the Paleoclimate Modelling Intercomparison Project (PMIP) since its inception, and thus many of the problems that might be associated with simulating such a radically different climate are well documented. The LGM state provides an ideal case study for evaluating climate model performance because the changes in forcing and temperature between the LGM and pre-industrial are of the same order of magnitude as those projected for the end of the 21st century. Thus, the CMIP6 LGM experiment could provide additional information that can be used to constrain estimates of climate sensitivity. The design of the Tier 1 LGM experiment (lgm) includes an assessment of uncertainties in boundary conditions, in particular through the use of different reconstructions of the ice sheets and of the change in dust forcing. Additional (Tier 2) sensitivity experiments have been designed to quantify feedbacks associated with land-surface changes and aerosol loadings, and to isolate the role of individual forcings. Model analysis and evaluation will capitalize on the relative abundance of paleoenvironmental observations and quantitative climate reconstructions already available for the LGM.

Description: Interannual variability in desert dust is widely observed and simulated, yet the sensitivity of these desert dust simulations to a particular meteorological dataset, as well as a particular model construction, is not well known. Here we use version 4 of the Community Atmospheric Model (CAM4) with the Community Earth System Model (CESM) to simulate dust forced by three different reanalysis meteorological datasets for the period 1990–2005. We then contrast the results of these simulations with dust simulated using online winds dynamically generated from sea surface temperatures, as well as with simulations conducted using other modeling frameworks but the same meteorological forcings, in order to determine the sensitivity of climate model output to the specific reanalysis dataset used. For the seven cases considered in our study, the different model configurations are able to simulate the annual mean of the global dust cycle, seasonality and interannual variability approximately equally well at the limited observational sites available. Overall, aerosol dust source strength has remained fairly constant during the time period from 1990 to 2005, although there is strong seasonal and some interannual variability simulated in the models and seen in the observations over this time period. The models predict similar seasonal variability and timing as one another, and obtain the best comparison to observations at the seasonal scale. Interannual variability comparisons to observations, as well as comparisons between models, suggest that interannual variability in dust is still difficult to simulate accurately, with averaged correlation coefficients of 0.1 to 0.6. Because of the large variability, at least one year of observations at most sites are needed to correctly observe the mean, but in some regions, particularly the remote oceans of the Southern Hemisphere, where interannual variability appears to be larger, 2–3 years of data are needed.

Description: A total of sixteen global chemistry transport models and general circulation models have participated in this study. Fourteen models have been evaluated with regard to their ability to reproduce near-surface observed number concentration of aerosol particle and cloud condensation nuclei (CCN), and derived cloud droplet number concentration (CDNC). Model results for the period 2011–2015 are compared with aerosol measurements (aerosol particle number, CCN and aerosol particle composition in the submicron fraction) from nine surface stations, located in Europe and Japan. The evaluation focuses on the ability of models to simulate the average across time state in diverse environments, and on the seasonal and short-term variability in the aerosol properties. There is no single model that systematically performs best across all environments represented by the observations. Models tend to underestimate the observed aerosol particle and CCN number concentrations, with average normalized mean bias (NMB) of all models and for all stations, where data are available, of −24 % and −35 % for particles with dry diameters 〉 50 nm and 〉 120 nm and −36 % and −34 % for CCN at supersaturations of 0.2 % and 1.0 %, respectively. Fifteen models have been used to produce ensemble annual median distributions of relevant parameters. The model diversity (defined as the ratio of standard deviation to mean) is up to about 3 for simulated N3 (number concentration of particles with dry diameters larger than 3 nm) and up to about 1 for simulated CCN in the extra-polar regions. An additional model has been used to investigate potential causes of model diversity in CCN and bias compared to the observations by performing a perturbed parameter ensemble (PPE) accounting for uncertainties in 26 aerosol-related model input parameters. This PPE suggests that biogenic secondary organic aerosol formation and the hygroscopic properties of the organic material are likely to be the major sources of CCN uncertainty in summer, with dry deposition and cloud processing being dominant in winter. Models capture the relative amplitude of seasonal variability of the aerosol particle number concentration for all studied particle sizes with available observations (dry diameters larger than 50, 80 and 120 nm). The short-term persistence time (on the order of a few days) of CCN concentrations, which is a measure of aerosol dynamic behavior in the models, is underestimated on average by the models by 40 % during winter and 20 % in summer. In contrast to the large spread in simulated aerosol particle and CCN number concentrations, the CDNC derived from simulated CCN spectra is less diverse and in better agreement with CDNC estimates consistently derived from the observations (average NMB −17 % and −22 % for updraft velocities 0.3 and 0.6 m s−1, respectively). In addition, simulated CDNC is in slightly better agreement with observationally-derived value at lower than at higher updraft velocities (index-of-agreement of 0.47 vs 0.50). The reduced spread of CDNC compared to that of CCN is attributed to the sublinear response of CDNC to aerosol particle number variations and the negative correlation between the sensitivities of CDNC to aerosol particle number concentration and to updraft velocity. Overall, we find that while CCN is controlled by both aerosol particle number and composition, CDNC is sensitive to CCN at low and moderate CCN concentrations and to the updraft velocity when CCN levels are high.