Sedimentary records indicate that atmospheric dust has increased substantially since preindustrial times. However, state-of-the-art global Earth system models (ESMs) are unable to capture this historical increase, posing challenges in assessing the impacts of desert dust on Earth's climate. To address this issue, we construct a globally gridded dust emission dataset (DustCOMMv1) spanning 1841–2000. We do so by combining 19 sedimentary records of dust deposition with observational and modeling constraints on the modern-day dust cycle. The derived emission dataset contains interdecadal variability of dust emissions as forced by the deposition flux records, which increased by approximately 50 % from 1851–1870 to 1981–2000. We further provide future dust emission datasets for 2000–2100 by assuming three possible scenarios for how future dust emissions will evolve. We evaluate the historical dust emission dataset and illustrate its effectiveness in enforcing a historical dust increase in ESMs by conducting a long-term (1851–2000) dust cycle simulation with the Community Earth System Model (CESM2). The simulated dust depositions are in reasonable agreement with the long-term increase in most sedimentary dust deposition records and with measured long-term trends in dust concentration at sites in Miami and Barbados. This contrasts with the CESM2 simulations using a process-based dust emission scheme and with simulations from the Coupled Model Intercomparison Project (CMIP6), which show little to no secular trends in dust deposition, concentration, and optical depth. The DustCOMM emissions thus enable ESMs to account for the historical radiative forcings (RFs), including due to dust direct interactions with radiation (direct RF). Our CESM2 simulations estimate a 1981–2000 minus 1851–1870 direct RF of −0.10 W m−2 by dust aerosols up to 10 µm in diameter (PM10) at the top of atmosphere (TOA). This global dust emission dataset thus enables models to more accurately account for historical aerosol forcings, thereby improving climate change projections such as those in the Intergovernmental Panel on Climate Change (IPCC) assessment reports.

Comprehensive land models are subject to significant parametric uncertainty, which can be hard to quantify due to the large number of parameters and high model computational costs. We constructed a large parameter perturbation ensemble (PPE) for the Community Land Model version 5.1 with biogeochemistry configuration (CLM5.1-BGC). We performed more than 2,000 simulations perturbing 211 parameters across six forcing scenarios. This provides an expansive data set, which can be used to identify the most influential parameters on a wide range of output variables globally, by biome, or by plant functional type. We found that parameter effects can exceed scenario effects and that a small number of parameters explains a large fraction of variance across our ensemble. The most important parameters can differ regionally and also based on the forcing scenario. The software infrastructure developed for this experiment has greatly reduced the human and computer time needed for CLM PPEs, which can facilitate routine investigation of parameter sensitivity and uncertainty, as well as automated calibration.

The terrestrial biosphere plays a major role in the global carbon cycle, and there is a recognized need for regularly updated estimates of land-atmosphere exchange at regional and global scales. An international ensemble of Dynamic Global Vegetation Models (DGVMs), known as the "Trends and drivers of the regional scale terrestrial sources and sinks of carbon dioxide" (TRENDY) project, quantifies land biophysical exchange processes and biogeochemistry cycles in support of the annual Global Carbon Budget assessments and the REgional Carbon Cycle Assessment and Processes, phase 2 project. DGVMs use a common protocol and set of driving data sets. A set of factorial simulations allows attribution of spatio-temporal changes in land surface processes to three primary global change drivers: changes in atmospheric CO2, climate change and variability, and Land Use and Land Cover Changes (LULCC). Here, we describe the TRENDY project, benchmark DGVM performance using remote-sensing and other observational data, and present results for the contemporary period. Simulation results show a large global carbon sink in natural vegetation over 2012-2021, attributed to the CO2 fertilization effect (3.8 +/- 0.8 PgC/yr) and climate (-0.58 +/- 0.54 PgC/yr). Forests and semi-arid ecosystems contribute approximately equally to the mean and trend in the natural land sink, and semi-arid ecosystems continue to dominate interannual variability. The natural sink is offset by net emissions from LULCC (-1.6 +/- 0.5 PgC/yr), with a net land sink of 1.7 +/- 0.6 PgC/yr. Despite the largest gross fluxes being in the tropics, the largest net land-atmosphere exchange is simulated in the extratropical regions. Around one third of human-induced CO2 emissions are absorbed by land ecosystems and thus act to mitigate climate change. It is essential to understand the processes, ecosystems and regions responsible for this natural carbon sink, to inform on the efficiency of the sinks into the future. These sinks are susceptible to year-to-year variation in response to climate variations and extremes. At the same time deforestation and other forms of land management are changing the land surface, which overall adds significantly to the human-induced CO2 emissions. There is a need to regularly update our estimate of land carbon dynamics to aid global stock takes for the Paris agreement to avoid dangerous climate change. Here we present an international initiative that on an annual basis assesses "Trends and drivers of the regional scale terrestrial sources and sinks of carbon dioxide" (TRENDY) using computer models of the land carbon cycle. We quantify the land sink during the contemporary period (2012-2021), and attribute to processes, mainly the large opposing effects of CO2 fertilization enhancing plant productivity and land-use change. Forests and semi-arid ecosystems are largely responsible for the mean and trend in the land sink, with the latter most important for its year-to-year variation. We quantify and attribute land carbon dynamics to underlying processes at regional scales, contributing bottom-up estimates to RECCAP-2 Models simulate a contemporary net land sink of 1.7 +/- 0.6 PgC/yr, with large opposing effects of CO2 fertilization and land-use change Despite the largest gross fluxes being in the tropics, the largest net land-atmosphere exchange is simulated in the extratropical regions

Lakes and reservoirs are an integral part of the terrestrial water cycle. In this work, we present the implementation of different water balance models of lakes and reservoirs into mizuRoute, a vector-based routing model, termed mizuRoute-Lake. As the main advantage of mizuRoute-Lake, users can choose between various parametric models implemented in mizuRoute-Lake. So far, three parametric models of lake and reservoir water balance, namely Hanasaki, HYPE, and D & ouml;ll are implemented in mizuRoute-Lake. In general, the parametric models relate the outflow from lakes or reservoirs to the storage and various parameters including inflow, demand, volume of storage, etc. Additionally, this flexibility allows users to easily evaluate and compare the effect of various water balance models for a lake without needing to reconfigure the routing model or change the parameters of other lakes or reservoirs in the modeling domain. Users can also use existing data such as historical observations or water management models to specify the behavior of a selected number of lakes and reservoirs within the modeling domain using the data-driven capability of mizuRoute-Lake. We demonstrate the flexibility of mizuRoute-Lake by presenting global, regional, and local scale applications. The development of mizuRoute-Lake paves the way for better integration of water management models, locally measured, and remotely sensed data sets in the context of Earth system modeling. Lakes and reservoirs impact the hydrological cycle. While there are many lake and reservoir models, a modeling framework typically encompasses a single selected lake and/or reservoir model representation. Including various lake discharge models in one model framework is useful for many applications of water resources. In this study, we present the implementation of various lake and reservoir water balance models in the mizuRoute-Lake. mizuRoute-Lake provides the flexibility to use and evaluate different lake and reservoir water balance models in river basins. Such flexibility can facilitate the integration with Earth system models, land surface models, water management models, and in situ or remote sensing observations. Additionally, it can help shed light on often neglected interactions and uncertainties in lake and reservoir models and their parameterizations within the hydrological cycle. mizuRoute-Lake offers the flexibility to select from multiple parametric models for lakes and reservoirs' water balance Users can leverage existing data, historical observations, water management models, and remotely sensed data, to specify lake behavior The versatility of mizuRoute-Lake is demonstrated through its application at different scales, global to local

Climate simulation uncertainties arise from internal variability, model structure, and external forcings. Model intercomparisons (such as the Coupled Model Intercomparison Project; CMIP) and single-model large ensembles have provided insight into uncertainty sources. Under the Community Earth System Model (CESM) project, large ensembles have been performed for CESM2 (a CMIP6-era model) and CESM1 (a CMIP5-era model). We refer to these as CESM2-LE and CESM1-LE. The external forcing used in these simulations has changed to be consistent with their CMIP generation. As a result, differences between CESM2-LE and CESM1-LE ensemble means arise from changes in both model structure and forcing. Here we present new ensemble simulations which allow us to separate the influences of these model structural and forcing differences. Our new CESM2 simulations are run with CMIP5 forcings equivalent to those used in the CESM1-LE. We find a strong influence of historical forcing uncertainty due to aerosol effects on simulated climate. For the historical period, forcing drives reduced global warming and ocean heat uptake in CESM2-LE relative to CESM1-LE that is counteracted by the influence of model structure. The influence of the model structure and forcing vary across the globe, and the Arctic exhibits a distinct signal that contrasts with the global mean. For the 21st century, the importance of scenario forcing differences (SSP3–7.0 for CESM2-LE and RCP8.5 for CESM1-LE) is evident. The new simulations presented here allow us to diagnose the influence of model structure on 21st century change, despite large scenario forcing differences, revealing that differences in the meridional distribution of warming are caused by model structure. Feedback analysis reveals that clouds and their impact on shortwave radiation explain many of these structural differences between CESM2 and CESM1. In the Arctic, albedo changes control transient climate evolution differences due to structural differences between CESM2 and CESM1.

Desert dust is an important atmospheric aerosol that affects the Earth's climate, biogeochemistry, and air quality. However, current Earth system models (ESMs) struggle to accurately capture the impact of dust on the Earth's climate and ecosystems, in part because these models lack several essential aeolian processes that couple dust with climate and land surface processes. In this study, we address this issue by implementing several new parameterizations of aeolian processes detailed in our companion paper in the Community Earth System Model version 2 (CESM2). These processes include (1) incorporating a simplified soil particle size representation to calculate the dust emission threshold friction velocity, (2) accounting for the drag partition effect of rocks and vegetation in reducing wind stress on erodible soils, (3) accounting for the intermittency of dust emissions due to unresolved turbulent wind fluctuations, and (4) correcting the spatial variability of simulated dust emissions from native to higher spatial resolutions on spatiotemporal dust variability. Our results show that the modified dust emission scheme significantly reduces the model bias against observations compared with the default scheme and improves the correlation against observations of multiple key dust variables such as dust aerosol optical depth (DAOD), surface particulate matter (PM) concentration, and deposition flux. Our scheme's dust also correlates strongly with various meteorological and land surface variables, implying higher sensitivity of dust to future climate change than other schemes' dust. These findings highlight the importance of including additional aeolian processes for improving the performance of ESM aerosol simulations and potentially enhancing model assessments of how dust impacts climate and ecosystem changes.

Water scarcity is exacerbated by rising water use and climate change, yet state-of-the-art Earth system models typically do not represent human water demand. Here we present an enhancement to the Community Earth System Model (CESM) and its land (CLM5) and river (MOSART) components by introducing sectoral water abstractions. The new module enables a better understanding of water demand and supply dynamics across various sectors, including domestic, livestock, thermoelectric, manufacturing, mining, and irrigation. The module conserves water by integrating abstractions from the land component with river component flows and dynamically calculates daily water scarcity based on local demand and supply. Through land-only simulations spanning 1971-2010, we verify our model against known water scarcity hotspots, historical global water withdrawal trends, and regional variations in water use. Our findings show that non-irrigative sectoral consumption has an insignificant effect on regional climate, while emphasizing the importance of including all sectors for water scarcity assessment capabilities. Despite its advancements, the model's limitations, such as its exclusive focus on river water abstractions, highlight areas for potential future refinement. This research paves the way for a more holistic representation of human-water interactions in ESMs, aiming to inform sustainable water management decisions in an evolving global landscape.

The spin-up problem hinders our ability to study some key issues in land carbon cycle modeling, such as sensitivity analysis. This work applies a new semi-analytical spin-up (SASU) framework to accelerate spin-up of Community Land Model matrix version 5. We evaluate the computational efficiency and steady-state consistency of ND, AD and SASU approaches using CLM5 at both one-site (Brazil) and the global scale. In addition, we test SASU in Parameter Perturbation Experiment (PPE) in CLM5.

Abstract. Desert dust is an important atmospheric aerosol that affects the Earth's climate, biogeochemistry, and air quality. However, current Earth system models (ESMs) struggle to accurately capture the impact of dust on the Earth’s climate and ecosystems, in part because these models lack several essential aeolian processes that couple dust with climate and land surface processes. In this study, we address this issue by implementing several new parameterizations of aeolian processes detailed in our companion paper into the Community Earth System Model version 2 (CESM2). These processes include (1) incorporating a realistic soil particle size distribution to calculate the dust emission threshold friction velocity, (2) accounting for the drag partition effect of rocks and vegetation in reducing wind stress on erodible soils, (3) accounting for the intermittency of dust emissions due to unresolved turbulent wind fluctuations, and (4) correcting the spatial variability of simulated dust emissions from native to higher spatial resolutions on spatiotemporal dust variability. Our results show that the modified dust emission scheme significantly reduces the model bias against observations compared to the default scheme and improves the correlation against observations of multiple key dust variables such as dust aerosol optical depth (DAOD), surface particulate matter (PM) concentration, and deposition flux. Our scheme’s dust also correlates strongly with various meteorological and land surface variables, implying higher sensitivity of dust to future climate change than other schemes’ dust. These findings highlight the importance of including additional aeolian processes for improving the performance of ESM aerosol simulations and potentially enhancing model assessments of how dust impacts climate and ecosystem changes.

Abstract A new configuration of the Community Earth System Model (CESM)/Community Atmosphere Model with full chemistry (CAM‐chem) supporting the capability of horizontal mesh refinement through the use of the spectral element (SE) dynamical core is developed and called CESM/CAM‐chem‐SE. Horizontal mesh refinement in CESM/CAM‐chem‐SE is unique and novel in that pollutants such as ozone are accurately represented at human exposure relevant scales while also directly including global feedbacks. CESM/CAM‐chem‐SE with mesh refinement down to ∼14 km over the conterminous US (CONUS) is the beginning of the Multi‐Scale Infrastructure for Chemistry and Aerosols (MUSICAv0). Here, MUSICAv0 is evaluated and used to better understand how horizontal resolution and chemical complexity impact ozone and ozone precursors over CONUS as compared to measurements from five aircraft campaigns, which occurred in 2013. This field campaign analysis demonstrates the importance of using finer horizontal resolution to accurately simulate ozone precursors such as nitrogen oxides and carbon monoxide. In general, the impact of using more complex chemistry on ozone and other oxidation products is more pronounced when using finer horizontal resolution where a larger number of chemical regimes are resolved. Large model biases for ozone near the surface remain in the Southeast US as compared to the aircraft observations even with updated chemistry and finer horizontal resolution. This suggests a need for adding the capability of replacing sections of global emission inventories with regional inventories, increasing the vertical resolution in the planetary boundary layer, and reducing model biases in meteorological variables such as temperature and clouds.

Abstract A vector‐river network explicitly uses realistic geometries of river reaches and catchments for spatial discretization in a river model. This enables improving the accuracy of the physical properties of the modeled river system, compared to a gridded river network that has been used in Earth System Models. With a finer‐scale river network, resolving smaller‐scale river reaches, there is a need for efficient methods to route streamflow and its constituents throughout the river network. The purpose of this study is twofold: (1) develop a new method to decompose river networks into hydrologically independent tributary domains, where routing computations can be performed in parallel; and (2) perform global river routing simulations with two global river networks, with different scales, to examine the computational efficiency and the differences in discharge simulations at various temporal scales. The new parallelization method uses a hierarchical decomposition strategy, where each decomposed tributary is further decomposed into many sub‐tributary domains, enabling hybrid parallel computing. This parallelization scheme has excellent computational scaling for the global domain where it is straightforward to distribute computations across many independent river basins. However, parallel computing for a single large basin remains challenging. The global routing experiments show that the scale of the vector‐river network has less impact on the discharge simulations than the runoff input that is generated by the combination of land surface model and meteorological forcing. The scale of vector‐river networks needs to consider the scale of local hydrologic features such as lakes that are to be resolved in the network.

During the Last Interglacial, approximately 129 to 116 ka (thousand years ago), the Arctic summer climate was warmer than the present, and the Greenland Ice Sheet retreated to a smaller extent than its current state. Previous model‐derived and geological reconstruction estimates of the sea‐level contribution of the Greenland Ice Sheet during the Last Interglacial vary widely. Here, we conduct a transient climate simulation from 127 to 119 ka using the Community Earth System Model (CESM2), which includes a dynamic ice sheet component (the Community Ice Sheet Model, CISM2) that is interactively coupled to the atmosphere, land, ocean, and sea ice components. Vegetation distribution is updated every 500 years based on biomes simulated using a monthly climatology to force the BIOME4 equilibrium vegetation model. Results show a substantial retreat of the Greenland Ice Sheet, reaching a minimum extent at 121.9 ka, equivalent to a 3.0 m rise in sea level relative to the present day, followed by gradual regrowth. In contrast, a companion simulation employing static vegetation based on pre‐industrial conditions shows a much smaller ice‐sheet retreat, highlighting the importance of the changes in high‐latitude vegetation distribution for amplifying the ice‐sheet response.

Abstract Earth system models (ESMs) have been rapidly developed in recent decades to advance our understanding of climate change‐carbon cycle feedback. However, those models are massive in coding, require expensive computational resources, and have difficulty in diagnosing their performance. It is highly desirable to develop ESMs with modularity and effective diagnostics. Toward these goals, we implemented a matrix approach to the Community Land Model version 5 (CLM5) to represent carbon and nitrogen cycles. Specifically, we reorganized 18 balance equations each for carbon and nitrogen cycles among the 18 vegetation pools in the original CLM5 into two matrix equations. Similarly, 140 balance equations each for carbon and nitrogen cycles among the 140 soil pools were reorganized into two additional matrix equations. The vegetation carbon and nitrogen matrix equations are connected to soil matrix equations via litterfall. The matrix equations fully reproduce simulations of carbon and nitrogen dynamics by the original model. The computational cost for forwarding simulation of the CLM5 matrix model was 26% more expensive than the original model, largely due to calculation of additional diagnostic variables, but the spin‐up computational cost was significantly saved. We showed a case study on modeled soil carbon storage under two forcing data sets to illustrate the diagnostic capability that the matrix approach uniquely offers to understand simulation results of global carbon and nitrogen dynamics. The successful implementation of the matrix approach to CLM5, one of the most complex land models, demonstrates that most, if not all, the biogeochemical models can be reorganized into the matrix form to gain high modularity, effective diagnostics, and accelerated spin‐up.

Results are presented and compared for the Community Earth System Model version 2 (CESM2) simulations of the middle Holocene (MH, 6 ka) and Last Interglacial (LIG, 127 ka). These simulations are designated as Tier 1 experiments ( midHolocene and lig127k ) for the Coupled Model Intercomparison Project phase 6 (CMIP6) and the Paleoclimate Modeling Intercomparison Project phase 4 (PMIP4). They use the low‐top, standard 1° version of CESM2 contributing to CMIP6 DECK, historical, and future projection simulations, and to other modeling intercomparison projects. The midHolocene and lig127k provide the opportunity to examine the responses in CESM2 to the orbitally induced changes in the seasonal and latitudinal distribution of insolation. The insolation anomalies result in summer warming over the Northern Hemisphere continents, reduced Arctic summer minimum sea ice, and increased areal extent of the North African monsoon. The Arctic remains warm throughout the year. These changes are greater in the lig127k than midHolocene simulation. Other notable changes are reduction of the Niño3.4 variability and Drake Passage transport and a small increase in the Atlantic Meridional Overturning Circulation from the piControl to midHolocene to lig127k simulation. Comparisons to paleo‐data and to simulations from previous model versions are discussed. Possible reasons for mismatches with the paleo‐observations are proposed, including missing processes in CESM2, simplifications in the CMIP6 protocols for these experiments, and dating and calibration uncertainties in the data reconstructions.

The Community Land Model (CLM) is the land component of the Community Earth System Model (CESM) and is used in several global and regional modeling systems. In this paper, we introduce model developments included in CLM version 5 (CLM5), which is the default land component for CESM2. We assess an ensemble of simulations, including prescribed and prognostic vegetation state, multiple forcing data sets, and CLM4, CLM4.5, and CLM5, against a range of metrics including from the International Land Model Benchmarking (ILAMBv2) package. CLM5 includes new and updated processes and parameterizations: (1) dynamic land units, (2) updated parameterizations and structure for hydrology and snow (spatially explicit soil depth, dry surface layer, revised groundwater scheme, revised canopy interception and canopy snow processes, updated fresh snow density, simple firn model, and Model for Scale Adaptive River Transport), (3) plant hydraulics and hydraulic redistribution, (4) revised nitrogen cycling (flexible leaf stoichiometry, leaf N optimization for photosynthesis, and carbon costs for plant nitrogen uptake), (5) global crop model with six crop types and time-evolving irrigated areas and fertilization rates, (6) updated urban building energy, (7) carbon isotopes, and (8) updated stomatal physiology. New optional features include demographically structured dynamic vegetation model (Functionally Assembled Terrestrial Ecosystem Simulator), ozone damage to plants, and fire trace gas emissions coupling to the atmosphere. Conclusive establishment of improvement or degradation of individual variables or metrics is challenged by forcing uncertainty, parametric uncertainty, and model structural complexity, but the multivariate metrics presented here suggest a general broad improvement from CLM4 to CLM5.

Fire emissions are a critical component of carbon and nutrient cycles and strongly affect climate and air quality. Dynamic global vegetation models (DGVMs) with interactive fire modeling provide important estimates for long-term and large-scale changes in fire emissions. Here we present the first multi-model estimates of global gridded historical fire emissions for 1700-2012, including carbon and 33 species of trace gases and aerosols. The dataset is based on simulations of nine DGVMs with different state-of-the-art global fire models that participated in the Fire Modeling Intercomparison Project (FireMIP), using the same and standardized protocols and forcing data, and the most up-to-date fire emission factor table based on field and laboratory studies in various land cover types. We evaluate the simulations of present-day fire emissions by comparing them with satellite-based products. The evaluation results show that most DGVMs simulate present-day global fire emission totals within the range of satellite-based products. They can capture the high emissions over the tropical savannas and low emissions over the arid and sparsely vegetated regions, and the main features of seasonality. However, most models fail to simulate the interannual variability, partly due to a lack of modeling peat fires and tropical deforestation fires. Before the 1850s, all models show only a weak trend in global fire emissions, which is consistent with the multi-source merged historical reconstructions used as input data for CMIP6. On the other hand, the trends are quite different among DGVMs for the 20th century, with some models showing an increase and others a decrease in fire emissions, mainly as a result of the discrepancy in their simulated responses to human population density change and land use and land cover change (LULCC). Our study provides an important dataset for further development of regional and global multi-source merged historical reconstructions, analyses of the historical changes in fire emissions and their uncertainties, and quantification of the role of fire emissions in the Earth system. It also highlights the importance of accurately modeling the responses of fire emissions to LULCC and population density change in reducing uncertainties in historical reconstructions of fire emissions and providing more reliable future projections.