Peatlands are globally important stores of carbon which are being increasingly threatened by wildfires with knock-on effects on the climate system. Here we introduce a novel peat fire parameterization in the northern high latitudes to the INFERNO global fire model. Representing peat fires increases annual burnt area across the high latitudes, alongside improvements in how we capture year-to-year variation in burning and emissions. 

We present a collection of performance metrics gathered during the Coupled Model Intercomparison Project Phase 6 (CMIP6), a worldwide initiative to study climate change. We analyse the metrics that resulted from collaboration efforts among many partners and models and describe our findings to demonstrate the utility of our study for the scientific community. The research contributes to understanding climate modelling performance on the current high-performance computing (HPC) architectures.

Especially over the midlatitudes, precipitation is mainly formed via the ice phase. In this study we focus on the initial snow formation process in the ICON-AES, the aggregation process. We use a stochastical approach for the aggregation parameterization and investigate the influence in the ICON-AES. Therefore, a distribution function of cloud ice is created, which is evaluated with satellite data. The new approach leads to cloud ice loss and an improvement in the process rate bias.

Seven new experimental setups designed to interfere with cloud radiative heating have been added to the Energy Exascale Earth System Model (E3SM). These experiments include both those that test the mean impact of cloud radiative heating and those examining its covariance with circulations. This paper documents the code changes and steps needed to run these experiments. Results corroborate prior findings for how cloud radiative heating impacts circulations and rainfall patterns.

The Reservoir Assessment Tool (RAT) merges satellite data with hydrological models, enabling robust estimation of reservoir parameters like inflow, outflow, surface area, and storage changes around the world. Version 3.0 of RAT lowers the barrier of entry for new users and achieves scalability and computational efficiency. RAT 3.0 also facilitates open-source development of functions for continuous improvement to mobilize and empower the global water management community.

In order to improve Sargassum drift forecasting in the Caribbean area, drift models can be forced by higher-resolution ocean currents. To this goal a 3 km resolution regional ocean model has been developed. Its assessment is presented with a particular focus on the reproduction of fine structures representing key features of the Caribbean region dynamics and Sargassum transport. The simulated propagation of a North Brazil Current eddy and its dissipation was found to be quite realistic.

The paper addresses synthetic dataset challenges in nonlinear constitutive modeling of soils, providing two datasets: one with 2000 soil types, 40 test conditions each (160 000 triaxial tests), and a second with 2048 soil types, 42 test conditions each (172 032 triaxial tests). Each dataset is a 4000×10 matrix applicable for multivariate forecasting and geotechnical simulations. In addition, a new Python code is introduced, empowering researchers to leverage Python packages for NorSand analyses.

We present the wflow_sbm distributed hydrological model, recently released by Deltares, as part of the Wflow.jl open-source modelling framework in the programming language Julia. Wflow_sbm has a fast runtime, making it suitable for large-scale modelling. Wflow_sbm models can be set a priori for any catchment with the Python tool HydroMT-Wflow based on globally available datasets, which results in satisfactory to good performance (without much tuning). We show this for a number of specific cases.

We provide an R package to compute two biosphere integrity metrics that can be applied to simulations of vegetation growth from the dynamic global vegetation model LPJmL. The pressure metric BioCol indicates that we humans modify and extract > 20 % of the potential preindustrial natural biomass production. The ecosystems state metric EcoRisk shows a high risk of ecosystem destabilization in many regions as a result of climate change and land, water, and fertilizer use.

For the first time, we coupled a regional climate chemistry model, RegCM-Chem, with a dynamic vegetation model, YIBs, to create a regional climate–chemistry–ecology model, RegCM-Chem–YIBs. We applied it to simulate climatic, chemical, and ecological parameters in East Asia and fully validated it on a variety of observational data. Results show that RegCM-Chem–YIBs model is a valuable tool for studying the terrestrial carbon cycle, atmospheric chemistry, and climate change on a regional scale.

We present long-term simulations with interactive icebergs in the Southern Ocean. By melting, icebergs reduce the temperature and salinity of the surrounding ocean. In our simulations, we find that this cooling effect of iceberg melting is not limited to the surface ocean but also reaches the deep ocean and propagates northward into all ocean basins. Additionally, the formation of deep-water masses in the Southern Ocean is enhanced.

Atmospheric models struggle to reproduce the strong temperature inversions in the vicinity of the surface over forested areas in the Arctic winter. In this paper, we develop modified simplified versions of surface layer schemes widely used by the community. Our modifications are used to correct the fact that original schemes place strong limits on the turbulent collapse, leading to a lower surface temperature gradient at low wind speeds. Modified versions show a better performance.

Any interpretation of climate model data requires a comprehensive evaluation of the model performance. Numerous error metrics exist for this purpose, and each focuses on a specific aspect of the relationship between reference and model data. Thus, a comprehensive evaluation demands the use of multiple error metrics. However, this can lead to confusion. We propose a clustering technique to reduce the number of error metrics needed and a composite error metric to simplify the interpretation.

The new offline particle tracking package, Tracker v1.1, is introduced to the Regional Ocean Modeling System, featuring an efficient nearest-neighbor algorithm to enhance particle-tracking speed. Its performance was evaluated against four other tracking packages and passive dye. Despite unique features, all packages yield comparable results. Running multiple packages within the same circulation model allows comparison of their performance and ease of use.

Based on the TenStream solver, we present a new method to accelerate 3D radiative transfer towards the speed of currently used 1D solvers. Using a shallow-cumulus-cloud time series, we evaluate the performance of this new solver in terms of both speed and accuracy. Compared to a 3D benchmark simulation, we show that our new solver is able to determine much more accurate irradiances and heating rates than a 1D δ-Eddington solver, even when operated with a similar computational demand.

In this paper, we coupled the SURFEX and MEGAN models. The aim of this coupling is to improve the estimation of biogenic fluxes by using the SURFEX canopy environment model. The coupled model results were validated and several sensitivity tests were performed. The coupled-model total annual isoprene flux is 442 Tg; this value is within the range of other isoprene estimates reported. The ultimate aim of this coupling is to predict the impact of climate change on biogenic emissions.

Climate field reconstruction encompasses methods that estimate the evolution of climate in space and time based on natural archives. It is useful to investigate climate variations and validate climate models, but its implementation and use can be difficult for non-experts. This paper introduces a user-friendly Python package called cfr to make these methods more accessible, thanks to the computational and visualization tools that facilitate efficient and reproducible research on past climates.

Models of rocky coastal erosion help us understand the controls on coastal morphology and evolution. In this paper, we present a simplified model of coastline erosion driven by either uniform erosion where coastline erosion is constant or wave-driven erosion where coastline erosion is a function of the wave power. This model can be used to evaluate how coastline changes reflect climate, sea-level history, material properties, and the relative influence of different erosional processes.

In this study, a land surface image assimilation system capable of optimizing the spatial structure of the background field is constructed by introducing the curvelet analysis method and taking the similarity of image structure as a weak constraint. The findings demonstrate that the assimilation of surface soil moisture observation images effectively and reasonably enhances the spatial structure of soil moisture analysis field.

This study developed new meshing workflows to enable the automatic generation of meshes that follow geological models. The workflow allows for importing several geological models as input for Gmsh and later exporting the same number of high-quality meshes. This way, geological uncertainty is directly included in the numerical simulations. This study evaluates the impact of the geological uncertainty on thermohydraulic performance of two reservoirs for high-temperature heat storage applications.

This paper describes and evaluates an improvement to the representation of acetone in the GISS ModelE2.1 Earth system model. We simulate acetone's concentration and transport across the atmosphere as well as its dependence on chemistry, the ocean, and various global emissions. Comparisons of our model’s estimates to past modeling studies and field measurements have shown encouraging results. Ultimately, this paper contributes to a broader understanding of acetone's role in the atmosphere.

Anthropogenic aerosol emissions are an essential part of global aerosol models. Significant errors can exist from the loss of emission heterogeneity. We introduced an emission treatment that significantly improved aerosol emission heterogeneity in high-resolution model simulations, with improvements in simulated aerosol surface concentrations. The emission treatment will provide a more accurate representation of aerosol emissions and their effects on climate.

Atmospheric rivers (ARs) are extreme weather events that can alleviate drought or cause billions of US dollars in flood damage. We train convolutional neural networks (CNNs) to detect ARs with an estimate of the uncertainty. We present a framework to generalize these CNNs to a variety of datasets of past, present, and future climate. Using a simplified simulation of the Earth's atmosphere, we validate the CNNs. We explore the role of ARs in maintaining energy balance in the Earth system.

We have developed a new data assimilation framework by coupling an integrated hydrological model HydroGeoSphere with the data assimilation software PDAF. Compared to existing hydrological data assimilation systems, the advantage of our newly developed framework lies in its consideration of the physically based model; its large selection of different assimilation algorithms; and its modularity with respect to the combination of different types of observations, states and parameters.