We introduce the Geoscientific ADjoint Optimisation PlaTform (G-ADOPT), designed for inverse modelling of Earth system processes, with an initial focus on mantle dynamics. G-ADOPT is built upon Firedrake, Dolfin-Adjoint and the Rapid Optimisation Library, which work together to optimise models using an adjoint method, aligning them with seismic and geologic datasets. We demonstrate G-ADOPT's ability to reconstruct mantle evolution and thus be a powerful tool in geosciences.

Explosive volcanic eruptions lead to long-lived, microscopic particles in the upper atmosphere which act to cool the Earth's surface by reflecting the Sun's light back to space. We include and test this process in a global climate model, E3SM. E3SM is tested against satellite and balloon observations of the 1991 eruption of Mt. Pinatubo, showing that with these particles in the model we reasonably recreate Pinatubo and its global effects. We also explore how particle size leads to these effects.

This study presents a coupling of the large-scale glacier model OGGM and the hydrological model CWatM. Projected future increase in discharge is less strong while future decrease in discharge is stronger when glacier runoff is explicitly included in the large-scale hydrological model. This is because glacier runoff is projected to decrease in nearly all basins. We conclude that an improved glacier representation can prevent underestimating future discharge changes in large river basins.

This study explores methods to reconstruct diurnal variations in skin sea surface temperature in a model of the Mediterranean Sea. Our new approach, considering chlorophyll concentration, enhances spatial and temporal variations in the warm layer. Comparative analysis shows context-dependent improvements. The proposed "chlorophyll-interactive" method brings the surface net total heat flux closer to zero annually, despite a net heat loss from the ocean to the atmosphere.

We investigate numerical convergence properties of a particle-based numerical cloud microphysics model (SDM) and a double-moment bulk scheme for simulating a marine stratocumulus case, compare their results with model intercomparison project results, and present possible explanations for the different results of the SDM and the bulk scheme. Aerosol processes can be accurately simulated using SDM, and this may be an important factor affecting the behavior and morphology of marine stratocumulus.

Modeling global freshwater solid and nutrient loads, in both magnitude and form, is imperative for understanding emerging eutrophication problems. Such efforts, however, have been challenged by the difficulty of balancing details of freshwater biogeochemical processes with limited knowledge, input, and validation datasets. Here we develop a global freshwater model that resolves intertwined algae, solid, and nutrient dynamics and provide performance assessment against measurement-based estimates.

A Merged Observatory Data File (MODF) format to systematically collate complex atmosphere, ocean, and terrestrial data sets collected by multiple instruments during field campaigns is presented. The MODF format is also designed to be applied to model output data, yielding format-matching Merged Model Data Files (MMDFs). MODFs plus MMDFs will augment and accelerate the synergistic use of model results with observational data to increase understanding and predictive skill.

The new process-based hydrological toolbox model, RoGeR (https://roger.readthedocs.io/), can be used to estimate the components of the hydrological cycle and the related travel times of pollutants through parts of the hydrological cycle. These estimations may contribute to effective water resources management. This paper presents the toolbox concept and provides a simple example of providing estimations to water resources management.

Ice sheets present a thickness of a few kilometres, leading to a vertical deformation of the crust of up to a kilometre. This process depends on properties of the solid Earth, which can be regionally very different. We propose a model that accounts for this often-ignored heterogeneity and run 100 000 simulation years in minutes. Thus, the evolution of ice sheets is modeled with better accuracy, which is critical for a good mitigation of climate change and, in particular, sea-level rise.

River floods are among the most devastating natural hazards. We propose a flood model with a statistical approach based on openly available data. The model is integrated in a framework for estimating impacts of physical hazards. Although the model only agrees moderately with satellite-detected flood extents, we show that it can be used for forecasting the magnitude of flood events in terms of socio-economic impacts and for comparing these with past events.

Building on previous analysis tools developed for atmospheric science, the original release of the Tracking and Object-Based Analysis (tobac) Python package, v1.2, was open-source, modular, and insensitive to the type of gridded input data. Here, we present the latest version of tobac, v1.5, which substantially improves scientific capabilities and computational efficiency from the previous version. These enhancements permit new uses for tobac in atmospheric science and potentially other fields.

Infrared CO₂ cooling of the middle and upper atmosphere is increasing. We developed a new routine for very fast and accurate calculations of this cooling in general circulation models. The new algorithm accounts for non-local thermodynamic equilibrium and is about 1000 times faster than the standard matrix algorithms. It is based on advanced techniques for non-equilibrium emission calculations in stellar atmospheres, which so far have not been used in Earth’s and planetary atmospheres.

Updating the Yasso07 soil C model's dependency on decomposition with a hump-shaped Ricker moisture function improved modelled soil organic C (SOC) stocks in a catena of mineral and organic soils in boreal forest. The Ricker function, set to peak at a rate of 1 and calibrated against SOC and CO₂ data using a Bayesian approach, showed a maximum in well-drained soils. Using SOC and CO₂ data together with the moisture only from the topsoil humus was crucial for accurate model estimates.

Scientists often use models to study complex processes, like the movement of ice sheets, and compare them to measurements for estimating quantities that are hard to measure. We highlight an approach that ensures accurate results from point data sources (e.g. height measurements) by evaluating the numerical solution at true point locations. This method improves accuracy, aids communication between scientists, and is well-suited for integration with specialised software that automates processes.

STORM v.2 (short for STOchastic Rainfall Model version 2.0) is an open-source and user-friendly modelling framework for simulating rainfall fields over a basin. It also allows simulating the impact of plausible climate change either on the total seasonal rainfall or the storm’s maximum intensity.

This study uses a deep learning method to upscale the time resolution of paleoclimate simulations to 1 million years. This improved resolution allows a climate-biogeochemical model to more accurately predict climate shifts. The method may be critical in developing new fully continuous methods that are able to be applied over a moving continental surface in deep time with high resolution at reasonable computational expense.