2019, Drüke, Markus, Forkel, Matthias, von Bloh, Werner, Sakschewski, Boris, Cardoso, Manoel, Bustamante, Mercedes, Kurths, Jürgen, Thonicke, Kirsten

Vegetation fires influence global vegetation distribution, ecosystem functioning, and global carbon cycling. Specifically in South America, changes in fire occurrence together with land-use change accelerate ecosystem fragmentation and increase the vulnerability of tropical forests and savannas to climate change. Dynamic global vegetation models (DGVMs) are valuable tools to estimate the effects of fire on ecosystem functioning and carbon cycling under future climate changes. However, most fire-enabled DGVMs have problems in capturing the magnitude, spatial patterns, and temporal dynamics of burned area as observed by satellites. As fire is controlled by the interplay of weather conditions, vegetation properties, and human activities, fire modules in DGVMs can be improved in various aspects. In this study we focus on improving the controls of climate and hence fuel moisture content on fire danger in the LPJmL4-SPITFIRE DGVM in South America, especially for the Brazilian fire-prone biomes of Caatinga and Cerrado. We therefore test two alternative model formulations (standard Nesterov Index and a newly implemented water vapor pressure deficit) for climate effects on fire danger within a formal model–data integration setup where we estimate model parameters against satellite datasets of burned area (GFED4) and aboveground biomass of trees. Our results show that the optimized model improves the representation of spatial patterns and the seasonal to interannual dynamics of burned area especially in the Cerrado and Caatinga regions. In addition, the model improves the simulation of aboveground biomass and the spatial distribution of plant functional types (PFTs). We obtained the best results by using the water vapor pressure deficit (VPD) for the calculation of fire danger. The VPD includes, in comparison to the Nesterov Index, a representation of the air humidity and the vegetation density. This work shows the successful application of a systematic model–data integration setup, as well as the integration of a new fire danger formulation, in order to optimize a process-based fire-enabled DGVM. It further highlights the potential of this approach to achieve a new level of accuracy in comprehensive global fire modeling and prediction.

2018, Kokkola, Harri, Kühn, Thomas, Laakso, Anton, Bergman, Tommi, Lehtinen, Kari E. J., Mielonen, Tero, Arola, Antti, Stadtler, Scarlet, Korhonen, Hannele, Ferrachat, Sylvaine, Lohmann, Ulrike, Neubauer, David, Tegen, Ina, Siegenthaler-Le Drian, Colombe, Schultz, Martin G., Bey, Isabelle, Stier, Philip, Daskalakis, Nikos, Heald, Colette L., Romakkaniemi, Sami

In this paper, we present the implementation and evaluation of the aerosol microphysics module SALSA2.0 in the framework of the aerosol-chemistry-climate model ECHAM-HAMMOZ. It is an alternative microphysics module to the default modal microphysics scheme M7 in ECHAM-HAMMOZ. The SALSA2.0 implementation within ECHAM-HAMMOZ is evaluated against observations of aerosol optical properties, aerosol mass, and size distributions, comparing also to the skill of the M7 implementation. The largest differences between the implementation of SALSA2.0 and M7 are in the methods used for calculating microphysical processes, i.e., nucleation, condensation, coagulation, and hydration. These differences in the microphysics are reflected in the results so that the largest differences between SALSA2.0 and M7 are evident over regions where the aerosol size distribution is heavily modified by the microphysical processing of aerosol particles. Such regions are, for example, highly polluted regions and regions strongly affected by biomass burning. In addition, in a simulation of the 1991 Mt. Pinatubo eruption in which a stratospheric sulfate plume was formed, the global burden and the effective radii of the stratospheric aerosol are very different in SALSA2.0 and M7. While SALSA2.0 was able to reproduce the observed time evolution of the global burden of sulfate and the effective radii of stratospheric aerosol, M7 strongly overestimates the removal of coarse stratospheric particles and thus underestimates the effective radius of stratospheric aerosol. As the mode widths of M7 have been optimized for the troposphere and were not designed to represent stratospheric aerosol, the ability of M7 to simulate the volcano plume was improved by modifying the mode widths, decreasing the standard deviations of the accumulation and coarse modes from 1.59 and 2.0, respectively, to 1.2 similar to what was observed after the Mt. Pinatubo eruption. Overall, SALSA2.0 shows promise in improving the aerosol description of ECHAM-HAMMOZ and can be further improved by implementing methods for aerosol processes that are more suitable for the sectional method, e.g., size-dependent emissions for aerosol species and size-resolved wet deposition.

2018, Kim, HyeJin, Rosa, Isabel M. D., Alkemade, Rob, Leadley, Paul, Hurtt, George, Popp, Alexander, van Vuuren, Detlef P., Anthoni, Peter, Arneth, Almut, Baisero, Daniele, Caton, Emma, Chaplin-Kramer, Rebecca, Chini, Louise, De Palma, Adriana, Di Fulvio, Fulvio, Di Marco, Moreno, Espinoza, Felipe, Ferrier, Simon, Fujimori, Shinichiro, Gonzalez, Ricardo E., Gueguen, Maya, Guerra, Carlos, Harfoot, Mike, Harwood, Thomas D., Hasegawa, Tomoko, Haverd, Vanessa, Havlík, Petr, Hellweg, Stefanie, Hill, Samantha L. L., Hirata, Akiko, Hoskins, Andrew J., Janse, Jan H., Jetz, Walter, Johnson, Justin A., Krause, Andreas, Leclère, David, Martins, Ines S., Matsui, Tetsuya, Merow, Cory, Obersteiner, Michael, Ohashi, Haruka, Poulter, Benjamin, Purvis, Andy, Quesada, Benjamin, Rondinini, Carlo, Schipper, Aafke M., Sharp, Richard, Takahashi, Kiyoshi, Thuiller, Wilfried, Titeux, Nicolas, Visconti, Piero, Ware, Christopher, Wolf, Florian, Pereira, Henrique M.

To support the assessments of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services (IPBES), the IPBES Expert Group on Scenarios and Models is carrying out an intercomparison of biodiversity and ecosystem services models using harmonized scenarios (BES-SIM). The goals of BES-SIM are (1) to project the global impacts of land-use and climate change on biodiversity and ecosystem services (i.e., nature's contributions to people) over the coming decades, compared to the 20th century, using a set of common metrics at multiple scales, and (2) to identify model uncertainties and research gaps through the comparisons of projected biodiversity and ecosystem services across models. BES-SIM uses three scenarios combining specific Shared Socio-economic Pathways (SSPs) and Representative Concentration Pathways (RCPs)-SSP1xRCP2.6, SSP3xRCP6.0, SSP5xRCP8.6-to explore a wide range of land-use change and climate change futures. This paper describes the rationale for scenario selection, the process of harmonizing input data for land use, based on the second phase of the Land Use Harmonization Project (LUH2), and climate, the biodiversity and ecosystem services models used, the core simulations carried out, the harmonization of the model output metrics, and the treatment of uncertainty. The results of this collaborative modeling project will support the ongoing global assessment of IPBES, strengthen ties between IPBES and the Intergovernmental Panel on Climate Change (IPCC) scenarios and modeling processes, advise the Convention on Biological Diversity (CBD) on its development of a post-2020 strategic plans and conservation goals, and inform the development of a new generation of nature-centred scenarios.

2009-5-8, Günther, A., Thiel, C.

In this contribution we evaluated both the structurally-controlled failure susceptibility of the fractured Cretaceous chalk rocks and the topographically-controlled shallow landslide susceptibility of the overlying glacial sediments for the Jasmund cliff area on Rügen Island, Germany. We employed a combined methodology involving spatially distributed kinematical rock slope failure testing with tectonic fabric data, and both physically- and inventory-based shallow landslide susceptibility analysis. The rock slope failure susceptibility model identifies areas of recent cliff collapses, confirming its value in predicting the locations of future failures. The model reveals that toppling is the most important failure type in the Cretaceous chalk rocks of the area. The shallow landslide susceptibility analysis involves a physically-based slope stability evaluation which utilizes material strength and hydraulic conductivity data, and a bivariate landslide susceptibility analysis exploiting landslide inventory data and thematic information on ground conditioning factors. Both models show reasonable success rates when evaluated with the available inventory data, and an attempt was made to combine the individual models to prepare a map displaying both terrain instability and landslide susceptibility. This combination highlights unstable cliff portions lacking discrete landslide areas as well as cliff sections highly affected by past landslide events. Through a spatial integration of the rock slope failure susceptibility model with the combined shallow landslide assessment we produced a comprehensive landslide susceptibility map for the Jasmund cliff area.

2019, Lutz, Femke, Herzfeld, Tobias, Heinke, Jens, Rolinski, Susanne, Schaphoff, Sibyll, von Bloh, Werner, Stoorvogel, Jetse J., Müller, Christoph

The effects of tillage on soil properties, crop productivity, and global greenhouse gas emissions have been discussed in the last decades. Global ecosystem models have limited capacity to simulate the various effects of tillage. With respect to the decomposition of soil organic matter, they either assume a constant increase due to tillage or they ignore the effects of tillage. Hence, they do not allow for analysing the effects of tillage and cannot evaluate, for example, reduced tillage or no tillage (referred to here as “no-till”) practises as mitigation practices for climate change. In this paper, we describe the implementation of tillage-related practices in the global ecosystem model LPJmL. The extended model is evaluated against reported differences between tillage and no-till management on several soil properties. To this end, simulation results are compared with published meta-analyses on tillage effects. In general, the model is able to reproduce observed tillage effects on global, as well as regional, patterns of carbon and water fluxes. However, modelled N fluxes deviate from the literature values and need further study. The addition of the tillage module to LPJmL5 opens up opportunities to assess the impact of agricultural soil management practices under different scenarios with implications for agricultural productivity, carbon sequestration, greenhouse gas emissions, and other environmental indicators.

2018, Rammig, Anja, Heinke, Jens, Hofhansl, Florian, Verbeeck, Hans, Baker, Timothy R., Christoffersen, Bradley, Ciais, Philippe, De Deurwaerder, Hannes, Fleischer, Katrin, Galbraith, David, Guimberteau, Matthieu, Huth, Andreas, Johnson, Michelle, Krujit, Bart, Langerwisch, Fanny, Meir, Patrick, Papastefanou, Phillip, Sampaio, Gilvan, Thonicke, Kirsten, von Randow, Celso, Zang, Christian, Rödig, Edna

Comparing model output and observed data is an important step for assessing model performance and quality of simulation results. However, such comparisons are often hampered by differences in spatial scales between local point observations and large-scale simulations of grid cells or pixels. In this study, we propose a generic approach for a pixel-to-point comparison and provide statistical measures accounting for the uncertainty resulting from landscape variability and measurement errors in ecosystem variables. The basic concept of our approach is to determine the statistical properties of small-scale (within-pixel) variability and observational errors, and to use this information to correct for their effect when large-scale area averages (pixel) are compared to small-scale point estimates. We demonstrate our approach by comparing simulated values of aboveground biomass, woody productivity (woody net primary productivity, NPP) and residence time of woody biomass from four dynamic global vegetation models (DGVMs) with measured inventory data from permanent plots in the Amazon rainforest, a region with the typical problem of low data availability, potential scale mismatch and thus high model uncertainty. We find that the DGVMs under- and overestimate aboveground biomass by 25% and up to 60%, respectively. Our comparison metrics provide a quantitative measure for model-data agreement and show moderate to good agreement with the region-wide spatial biomass pattern detected by plot observations. However, all four DGVMs overestimate woody productivity and underestimate residence time of woody biomass even when accounting for the large uncertainty range of the observational data. This is because DGVMs do not represent the relation between productivity and residence time of woody biomass correctly. Thus, the DGVMs may simulate the correct large-scale patterns of biomass but for the wrong reasons. We conclude that more information about the underlying processes driving biomass distribution are necessary to improve DGVMs. Our approach provides robust statistical measures for any pixel-to-point comparison, which is applicable for evaluation of models and remote-sensing products.

2017, Dipu, Sudhakar, Quaas, Johannes, Wolke, Ralf, Stoll, Jens, Mühlbauer, Andreas, Sourdeval, Odran, Salzmann, Marc, Heinold, Bernd, Tegen, Ina

The regional atmospheric model Consortium for Small-scale Modeling (COSMO) coupled to the Multi-Scale Chemistry Aerosol Transport model (Muscat) is extended in this work to represent aerosol-cloud interactions. Previously, only one-way interactions (scavenging of aerosol and in-cloud chemistry) and aerosol-radiation interactions were included in this model. The new version allows for a microphysical aerosol effect on clouds. For this, we use the optional two-moment cloud microphysical scheme in COSMO and the online-computed aerosol information for cloud condensation nuclei concentrations (Cccn), replacing the constant Cccn profile. In the radiation scheme, we have implemented a droplet-size-dependent cloud optical depth, allowing now for aerosol-cloud-radiation interactions. To evaluate the models with satellite data, the Cloud Feedback Model Intercomparison Project Observation Simulator Package (COSP) has been implemented. A case study has been carried out to understand the effects of the modifications, where the modified modeling system is applied over the European domain with a horizontal resolution of 0.25°g × g0.25°. To reduce the complexity in aerosol-cloud interactions, only warm-phase clouds are considered. We found that the online-coupled aerosol introduces significant changes for some cloud microphysical properties. The cloud effective radius shows an increase of 9.5g%, and the cloud droplet number concentration is reduced by 21.5g%.

2019, Neubauer, David, Ferrachat, Sylvaine, Siegenthaler-Le Drian, Colombe, Stier, Philip, Partridge, Daniel G., Tegen, Ina, Bey, Isabelle, Stanelle, Tanja, Kokkola, Harri, Lohmann, Ulrike

The global aerosol–climate model ECHAM6.3–HAM2.3 (E63H23) as well as the previous model versions ECHAM5.5–HAM2.0 (E55H20) and ECHAM6.1–HAM2.2 (E61H22) are evaluated using global observational datasets for clouds and precipitation. In E63H23, the amount of low clouds, the liquid and ice water path, and cloud radiative effects are more realistic than in previous model versions. E63H23 has a more physically based aerosol activation scheme, improvements in the cloud cover scheme, changes in the detrainment of convective clouds, changes in the sticking efficiency for the accretion of ice crystals by snow, consistent ice crystal shapes throughout the model, and changes in mixed-phase freezing; an inconsistency in ice crystal number concentration (ICNC) in cirrus clouds was also removed. Common biases in ECHAM and in E63H23 (and in previous ECHAM–HAM versions) are a cloud amount in stratocumulus regions that is too low and deep convective clouds over the Atlantic and Pacific oceans that form too close to the continents (while tropical land precipitation is underestimated). There are indications that ICNCs are overestimated in E63H23. Since clouds are important for effective radiative forcing due to aerosol–radiation and aerosol–cloud interactions (ERFari+aci) and equilibrium climate sensitivity (ECS), differences in ERFari+aci and ECS between the model versions were also analyzed. ERFari+aci is weaker in E63H23 (−1.0 W m−2) than in E61H22 (−1.2 W m−2) (or E55H20; −1.1 W m−2). This is caused by the weaker shortwave ERFari+aci (a new aerosol activation scheme and sea salt emission parameterization in E63H23, more realistic simulation of cloud water) overcompensating for the weaker longwave ERFari+aci (removal of an inconsistency in ICNC in cirrus clouds in E61H22). The decrease in ECS in E63H23 (2.5 K) compared to E61H22 (2.8 K) is due to changes in the entrainment rate for shallow convection (affecting the cloud amount feedback) and a stronger cloud phase feedback. Experiments with minimum cloud droplet number concentrations (CDNCmin) of 40 cm−3 or 10 cm−3 show that a higher value of CDNCmin reduces ERFari+aci as well as ECS in E63H23.

2019, Gidden, Matthew J., Riahi, Keywan, Smith, Steven J., Fujimori, Shinichiro, Luderer, Gunnar, Kriegler, Elmar, van Vuuren, Detlef P., van den Berg, Maarten, Feng, Leyang, Klein, David, Calvin, Katherine, Doelman, Jonathan C., Frank, Stefan, Fricko, Oliver, Harmsen, Mathijs, Hasegawa, Tomoko, Havlik, Petr, Hilaire, Jérôme, Hoesly, Rachel, Horing, Jill, Popp, Alexander, Stehfest, Elke, Takahashi, Kiyoshi

We present a suite of nine scenarios of future emissions trajectories of anthropogenic sources, a key deliverable of the ScenarioMIP experiment within CMIP6. Integrated assessment model results for 14 different emissions species and 13 emissions sectors are provided for each scenario with consistent transitions from the historical data used in CMIP6 to future trajectories using automated harmonization before being downscaled to provide higher emissions source spatial detail. We find that the scenarios span a wide range of end-of-century radiative forcing values, thus making this set of scenarios ideal for exploring a variety of warming pathways. The set of scenarios is bounded on the low end by a 1.9 W m−2 scenario, ideal for analyzing a world with end-of-century temperatures well below 2 ∘C, and on the high end by a 8.5 W m−2 scenario, resulting in an increase in warming of nearly 5 ∘C over pre-industrial levels. Between these two extremes, scenarios are provided such that differences between forcing outcomes provide statistically significant regional temperature outcomes to maximize their usefulness for downstream experiments within CMIP6. A wide range of scenario

2019, Kundzewicz, Zbigniew W., Su, Buda, Wang, Yanjun, Wang, Guojie, Wang, Guofu, Huang, Jinlong, Jiang, Tong

The present paper examines flood risk (composed of hazard, exposure, and vulnerability) in a range of spatial perspectives – from the global to the local scale. It deals with observed records, noting that flood damage has been increasing. It also tackles projections for the future, related to flood hazard and flood losses. There are multiple factors driving flood hazard and flood risk and there is a considerable uncertainty in our assessments, and particularly in projections for the future. Further, this paper analyses options for flood risk reduction in several spatial dimensions, from global framework to regional to local scales. It is necessary to continue examination of the updated records of flood-related indices, trying to search for changes that influence flood hazard and flood risk in river basins.