ALBERT

Description: Forest fires modify soil organic carbon and suppress soil respiration for many decades after the initial disturbance. The associated changes in soil autotrophic and heterotrophic respiration from the time of the forest fire, however, are less well characterized. The FireResp model predicts soil autotrophic and heterotrophic respiration parameterized with a novel dataset across a fire chronosequence in the Yukon and Northwest Territories of Canada. The dataset consisted of soil incubation experiments and field measurements of soil respiration and soil carbon stocks. The FireResp model contains submodels that consider a Q10 (exponential) model of respiration compared to models of heterotrophic respiration using Michaelis–Menten kinetics parameterized with soil microbial carbon. For model evaluation we applied the Akaike information criterion and compared predicted patterns in components of soil respiration across the chronosequence. Parameters estimated with data from the 5 cm soil depth had better model–data comparisons than parameters estimated with data from the 10 cm soil depth. The model–data fit was improved by including parameters estimated from soil incubation experiments. Models that incorporated microbial carbon with Michaelis–Menten kinetics reproduced patterns in autotrophic and heterotrophic soil respiration components across the chronosequence. Autotrophic respiration was associated with aboveground tree biomass at more recently burned sites, but this association was less robust at older sites in the chronosequence. Our results provide support for more structured soil respiration models than standard Q10 exponential models.

Description: Recent observational and modeling studies show that variations of stratospheric ozone and the resulting interaction between ozone and the stratospheric circulation play an important role in surface weather and climate. However, in many cases, computationally expensive coupled chemistry models have been used to study these effects. Here, we demonstrate how a much simpler idealized general circulation model (GCM) can be used for studying the impact of interactive stratospheric ozone on the circulation. The model, named Simplified Chemistry-Dynamical Model (SCDM V1.0), is constructed from a preexisting idealized GCM, into which a simplified linear ozone scheme and a parameterization for the shortwave radiative effects of ozone are implemented. The distribution and variability of stratospheric ozone simulated by the new model are in good agreement with the MERRA2 reanalysis, even for extreme circulation events such as Arctic stratospheric sudden warmings. The model thus represents a promising new tool for the study of ozone–circulation interaction in the stratosphere and its associated effects on tropospheric weather and climate.

Description: In this paper, we present a new version of the chemistry–climate model SOCOL-AERv2 supplemented by an iodine chemistry module. We perform three 20-year ensemble experiments to assess the validity of the modeled iodine and to quantify the effects of iodine on ozone. The iodine distributions obtained with SOCOL-AERv2-I agree well with AMAX-DOAS observations and with CAM-chem model simulations. For the present-day atmosphere, the model suggests that the iodine-induced chemistry leads to a 3 %–4 % reduction in the ozone column, which is greatest at high latitudes. The model indicates the strongest influence of iodine in the lower stratosphere with 30 ppbv less ozone at low latitudes and up to 100 ppbv less at high latitudes. In the troposphere, the account of the iodine chemistry reduces the tropospheric ozone concentration by 5 %–10 % depending on geographical location. In the lower troposphere, 75 % of the modeled ozone reduction originates from inorganic sources of iodine, 25 % from organic sources of iodine. At 50 hPa, the results show that the impacts of iodine from both sources are comparable. Finally, we determine the sensitivity of ozone to iodine by applying a 2-fold increase in iodine emissions, as it might be representative for iodine by the end of this century. This reduces the ozone column globally by an additional 1.5 %–2.5 %. Our results demonstrate the sensitivity of atmospheric ozone to iodine chemistry for present and future conditions, but uncertainties remain high due to the paucity of observational data of iodine species.

Description: The Decision Support System for Agrotechnology Transfer Cropping Systems Model (DSSAT-CSM) is a widely used crop modeling system that has been integrated into large-scale modeling frameworks. Existing frameworks generate spatially explicit simulated outputs at grid points through an inefficient process of translation from binary spatially referenced inputs to point-specific text input files, followed by translation and aggregation back from point-specific text output files to binary spatially referenced outputs. The main objective of this paper was to document the design and implementation of a parallel gridded simulation framework for DSSAT-CSM. A secondary objective was to provide preliminary analysis of execution time and scaling of the new parallel gridded framework. The parallel gridded framework includes improved code for model-internal data transfer, gridded input–output with the Network Common Data Form (NetCDF) library, and parallelization of simulations using the Message Passing Interface (MPI). Validation simulations with the DSSAT-CSM-CROPSIM-CERES-Wheat model revealed subtle discrepancies in simulated yield due to the rounding of soil parameters in the input routines of the standard DSSAT-CSM. Utilizing NetCDF for direct input–output produced a 3.7- to 4-fold reduction in execution time compared to R- and text-based input–output. Parallelization improved execution time for both versions with between 12.2- (standard version) and 13.4-fold (parallel gridded version) speed-up when comparing 1 to 16 compute cores. Estimates of parallelization of computation ranged between 99.2 % (standard version) and 97.3 % (parallel gridded version), indicating potential for scaling to higher numbers of compute cores.

Description: This paper presents the new and now open-source version 2.1 of the REgional Model of INvestments and Development (REMIND). REMIND, as an integrated assessment model (IAM), provides an integrated view of the global energy–economy–emissions system and explores self-consistent transformation pathways. It describes a broad range of possible futures and their relation to technical and socio-economic developments as well as policy choices. REMIND is a multiregional model incorporating the economy and a detailed representation of the energy sector implemented in the General Algebraic Modeling System (GAMS). It uses non-linear optimization to derive welfare-optimal regional transformation pathways of the energy-economic system subject to climate and sustainability constraints for the time horizon from 2005 to 2100. The resulting solution corresponds to the decentralized market outcome under the assumptions of perfect foresight of agents and internalization of external effects. REMIND enables the analyses of technology options and policy approaches for climate change mitigation with particular strength in representing the scale-up of new technologies, including renewables and their integration in power markets. The REMIND code is organized into modules that gather code relevant for specific topics. Interaction between different modules is made explicit via clearly defined sets of input and output variables. Each module can be represented by different realizations, enabling flexible configuration and extension. The spatial resolution of REMIND is flexible and depends on the resolution of the input data. Thus, the framework can be used for a variety of applications in a customized form, balancing requirements for detail and overall runtime and complexity.

Description: Fire constitutes a key process in the Earth system (ES), being driven by climate as well as affecting the climate by changing atmospheric composition and impacting the terrestrial carbon cycle. However, studies on the effects of fires on atmospheric composition, radiative forcing and climate have been limited to date, as the current generation of ES models (ESMs) does not include fully atmosphere–composition–vegetation coupled fires feedbacks. The aim of this work is to develop and evaluate a fully coupled fire–composition–climate ES model. For this, the INteractive Fires and Emissions algoRithm for Natural envirOnments (INFERNO) fire model is coupled to the atmosphere-only configuration of the UK's Earth System Model (UKESM1). This fire–atmosphere interaction through atmospheric chemistry and aerosols allows for fire emissions to influence radiation, clouds and generally weather, which can consequently influence the meteorological drivers of fire. Additionally, INFERNO is updated based on recent developments in the literature to improve the representation of human and/or economic factors in the anthropogenic ignition and suppression of fire. This work presents an assessment of the effects of interactive fire coupling on atmospheric composition and climate compared to the standard UKESM1 configuration that uses prescribed fire emissions. Results show a similar performance when using the fire–atmosphere coupling (the “online” version of the model) when compared to the offline UKESM1 that uses prescribed fire. The model can reproduce observed present-day global fire emissions of carbon monoxide (CO) and aerosols, despite underestimating the global average burnt area. However, at a regional scale, there is an overestimation of fire emissions over Africa due to the misrepresentation of the underlying vegetation types and an underestimation over equatorial Asia due to a lack of representation of peat fires. Despite this, comparing model results with observations of CO column mixing ratio and aerosol optical depth (AOD) show that the fire–atmosphere coupled configuration has a similar performance when compared to UKESM1. In fact, including the interactive biomass burning emissions improves the interannual CO atmospheric column variability and consequently its seasonality over the main biomass burning regions – Africa and South America. Similarly, for aerosols, the AOD results broadly agree with the Moderate Resolution Imaging Spectroradiometer (MODIS) and the Aerosol Robotic Network (AERONET) observations.

Description: We present a feasibility study for an object-based method to characterise thunderstorm properties in simulation data from convection-permitting weather models. An existing thunderstorm tracker, the Thunderstorm Identification, Tracking, Analysis and Nowcasting (TITAN) algorithm, was applied to thunderstorms simulated by the Advanced Research Weather Research and Forecasting (AR-WRF) weather model at convection-permitting resolution for a domain centred on Switzerland. Three WRF microphysics parameterisations were tested. The results are compared to independent radar-based observations of thunderstorms derived using the MeteoSwiss Thunderstorms Radar Tracking (TRT) algorithm. TRT was specifically designed to track thunderstorms over the complex Alpine topography of Switzerland. The object-based approach produces statistics on the simulated thunderstorms that can be compared to object-based observation data. The results indicate that the simulations underestimated the occurrence of severe and very large hail compared to the observations. Other properties, including the number of storm cells per day, geographical storm hotspots, thunderstorm diurnal cycles, and storm movement directions and velocities, provide a reasonable match to the observations, which shows the feasibility of the technique for characterisation of simulated thunderstorms over complex terrain.

Description: Out of 1150 Mha (million hectares) of forest designated primarily for production purposes in 2020, plantations accounted for 11 % (131 Mha) of this area and fulfilled more than 33 % of the global industrial roundwood demand. However, adding additional timber plantations to meet increasing timber demand intensifies competition for scarce land resources between different land uses such as food, feed, livestock and timber production. Despite the significance of plantations with respect to roundwood production, their importance in meeting the long-term timber demand and the implications of plantation expansion for overall land-use dynamics have not been studied in detail, in particular regarding the competition for land between agriculture and forestry in existing land-use models. This paper describes the extension of the modular, open-source land system Model of Agricultural Production and its Impact on the Environment (MAgPIE) using a detailed representation of forest land, timber production and timber demand dynamics. These extensions allow for a better understanding of the land-use dynamics (including competition for land) and the associated land-use change emissions of timber production. We show that the spatial cropland patterns differ when timber production is accounted for, indicating that timber plantations compete with cropland for the same scarce land resources. When plantations are established on cropland, it causes cropland expansion and deforestation elsewhere. Using the exogenous extrapolation of historical roundwood production from plantations, future timber demand and plantation rotation lengths, we model the future spatial expansion of forest plantations. As a result of increasing timber demand, we show a 177 % increase in plantation area by the end of the century (+171 Mha in 1995–2100). We also observe (in our model results) that the increasing demand for timber amplifies the scarcity of land, which is indicated by shifting agricultural land-use patterns and increasing yields from cropland compared with a case without forestry. Through the inclusion of new forest plantation and natural forest dynamics, our estimates of land-related CO2 emissions better match with observed data, in particular the gross land-use change emissions and carbon uptake (via regrowth), reflecting higher deforestation with the expansion of managed land and timber production as well as higher regrowth in natural forests and plantations.

Description: We present the dust module in the Multiscale Online Non-hydrostatic AtmospheRe CHemistry model (MONARCH) version 2.0, a chemical weather prediction system that can be used for regional and global modeling at a range of resolutions. The representations of dust processes in MONARCH were upgraded with a focus on dust emission (emission parameterizations, entrainment thresholds, considerations of soil moisture and surface cover), lower boundary conditions (roughness, potential dust sources), and dust–radiation interactions. MONARCH now allows modeling of global and regional mineral dust cycles using fundamentally different paradigms, ranging from strongly simplified to physics-based parameterizations. We present a detailed description of these updates along with four global benchmark simulations, which use conceptually different dust emission parameterizations, and we evaluate the simulations against observations of dust optical depth. We determine key dust parameters, such as global annual emission/deposition flux, dust loading, dust optical depth, mass-extinction efficiency, single-scattering albedo, and direct radiative effects. For dust-particle diameters up to 20 µm, the total annual dust emission and deposition fluxes obtained with our four experiments range between about 3500 and 6000 Tg, which largely depend upon differences in the emitted size distribution. Considering ellipsoidal particle shapes and dust refractive indices that account for size-resolved mineralogy, we estimate the global total (longwave and shortwave) dust direct radiative effect (DRE) at the surface to range between about −0.90 and −0.63 W m−2 and at the top of the atmosphere between −0.20 and −0.28 W m−2. Our evaluation demonstrates that MONARCH is able to reproduce key features of the spatiotemporal variability of the global dust cycle with important and insightful differences between the different configurations.

Description: The international Thermodynamic Equation of Seawater 2010 (TEOS-10) defined the enthalpy and entropy of seawater, thus enabling the global ocean heat content to be calculated as the volume integral of the product of in situ density, ρ, and potential enthalpy, h0 (with reference sea pressure of 0 dbar). In terms of Conservative Temperature, Θ, ocean heat content is the volume integral of ρcp0Θ, where cp0 is a constant “isobaric heat capacity”. However, many ocean models in the Coupled Model Intercomparison Project Phase 6 (CMIP6) as well as all models that contributed to earlier phases, such as CMIP5, CMIP3, CMIP2, and CMIP1, used EOS-80 (Equation of State – 1980) rather than the updated TEOS-10, so the question arises of how the salinity and temperature variables in these models should be physically interpreted, with a particular focus on comparison to TEOS-10-compliant observations. In this article we address how heat content, surface heat fluxes, and the meridional heat transport are best calculated using output from these models and how these quantities should be compared with those calculated from corresponding observations. We conclude that even though a model uses the EOS-80, which expects potential temperature as its input temperature, the most appropriate interpretation of the model's temperature variable is actually Conservative Temperature. This perhaps unexpected interpretation is needed to ensure that the air–sea heat flux that leaves and arrives in atmosphere and sea ice models is the same as that which arrives in and leaves the ocean model. We also show that the salinity variable carried by present TEOS-10-based models is Preformed Salinity, while the salinity variable of EOS-80-based models is also proportional to Preformed Salinity. These interpretations of the salinity and temperature variables in ocean models are an update on the comprehensive Griffies et al. (2016) paper that discusses the interpretation of many aspects of coupled Earth system models.