scholarly journals HIRM v1.0: a hybrid impulse response model for climate modeling and uncertainty analyses

2021 ◽  
Vol 14 (1) ◽  
pp. 365-375
Author(s):  
Kalyn Dorheim ◽  
Steven J. Smith ◽  
Ben Bond-Lamberty

Abstract. Simple climate models (SCMs) are frequently used in research and decision-making communities because of their flexibility, tractability, and low computational cost. SCMs can be idealized, flexibly representing major climate dynamics as impulse response functions, or process-based, using explicit equations to model possibly nonlinear climate and Earth system dynamics. Each of these approaches has strengths and limitations. Here we present and test a hybrid impulse response modeling framework (HIRM) that combines the strengths of process-based SCMs in an idealized impulse response model, with HIRM's input derived from the output of a process-based model. This structure enables the model to capture some of the major nonlinear dynamics that occur in complex climate models as greenhouse gas emissions transform to atmospheric concentration to radiative forcing to climate change. As a test, the HIRM framework was configured to emulate the total temperature of the simple climate model Hector 2.0 under the four Representative Concentration Pathways and the temperature response of an abrupt 4 times CO2 concentration step. HIRM was able to reproduce near-term and long-term Hector global temperature with a high degree of fidelity. Additionally, we conducted two case studies to demonstrate potential applications for this hybrid model: examining the effect of aerosol forcing uncertainty on global temperature and incorporating more process-based representations of black carbon into a SCM. The open-source HIRM framework has a range of applications including complex climate model emulation, uncertainty analyses of radiative forcing, attribution studies, and climate model development.

2020 ◽  
Author(s):  
Kalyn Dorheim ◽  
Steven Smith ◽  
Ben Bond-Lamberty

Abstract. Simple climate models (SCMs) are frequently used in research and decision-making communities because of their flexibility, tractability, and low computational cost. SCMs can be idealized, flexibly representing major climate dynamics as impulse response functions, or process-based, using explicit equations to model possibly nonlinear climate and earth system dynamics. Each of these approaches has strengths and limitations. Here we present and test a hybrid impulse response modeling framework (HIRM) that combines the strengths of process-based SCMs in an idealized impulse response model, with HIRM’s input derived from the output of a process-based model. This structure allows it to capture the crucial nonlinear dynamics frequently encountered in going from greenhouse gas emissions to atmospheric concentration to radiative forcing to climate change. As a test, the HIRM framework was configured to emulate total temperature of the simple climate model Hector 2.0 under the four Representative Concentration Pathways and the temperature response of an abrupt four times CO2 concentration step. HIRM was able to reproduce near-term and long-term Hector global temperature with a high degree of fidelity. Additionally, we conducted two case studies to demonstrate potential applications for this hybrid model: examining the effect of aerosol forcing uncertainty on global temperature, and incorporating more process-based representations of black carbon into a SCM. The open-source HIRM framework has a range of applications including complex climate model emulation, uncertainty analyses of radiative forcing, attribution studies, and climate model development.


2020 ◽  
Author(s):  
Nicholas James Leach ◽  
Zebedee Nicholls ◽  
Stuart Jenkins ◽  
Christopher J. Smith ◽  
John Lynch ◽  
...  

Abstract. Here we present a Generalised Impulse Response (GIR) model for use in probabilistic future climate and scenario exploration, integrated assessment, policy analysis and teaching. This model is based on a set of only six equations, which correspond to the standard Impulse Response model used for greenhouse gas metric calculations by the IPCC, plus one physically-motivated additional equation to represent state-dependent feedbacks on the response timescales of each greenhouse gas cycle. These six equations are simple and transparent enough to be easily understood and implemented in other models without reliance on the original source code, but flexible enough to reproduce observed well-mixed greenhouse gas (GHG) concentrations and atmospheric lifetimes, best-estimate effective radiative forcing, and temperature response. We describe the assumptions and methods used in selecting the default parameters, but emphasize that other methods would be equally valid: our focus here is on identifying a minimum level of structural complexity. The tunable nature of the model lends it to use as a fully transparent emulator of complex Earth System Models, such as those participating in CMIP6, while also reproducing the behaviour of other simple climate models. We argue that this GIR model is adequate to reproduce the global temperature response to global emissions and effective radiative forcing, and that it should be used as a lowest-common denominator to provide consistency and continuity between different climate assessments. The model design is such that it can be written in tabular data analysis software, such as Excel, increasing the potential user base considerably.


2018 ◽  
Author(s):  
Lennert B. Stap ◽  
Peter Köhler ◽  
Gerrit Lohmann

Abstract. The influence of long-term processes in the climate system, such as land ice changes, has to be compensated for when comparing climate sensitivity derived from paleodata with equilibrium climate sensitivity (ECS) calculated by climate models, which is only generated by a CO2 change. Several recent studies found that the impact these long-term processes have on global temperature cannot be quantified directly through the global radiative forcing they induce. This renders the approach of deconvoluting paleotemperatures through a partitioning based on radiative forcings inaccurate. Here, we therefore implement an efficacy factor ε[LI], that relates the impact of land ice changes on global temperature to that of CO2 changes, in our calculation of climate sensitivity from paleodata. We apply our new approach to a proxy-inferred paleoclimate dataset, and find an equivalent ECS of 5.6 ± 1.3 K per CO2 doubling. The substantial uncertainty herein is generated by the range in ε[LI] we use, which is based on a multi-model assemblage of simulated relative influences of land ice changes on the Last Glacial Maximum (LGM) temperature anomaly (46 ± 14 %). The low end of our ECS estimate, which concurs with estimates from other approaches, tallies with a large influence for land ice changes. To separately assess this influence, we analyse output of the PMIP3 climate model intercomparison project. From this data, we infer a functional intermodel relation between global and high-latitude temperature changes at the LGM with respect to the pre-industrial climate, and the temperature anomaly caused by a CO2 change. Applying this relation to our dataset, we find a considerable 64 % influence for land ice changes on the LGM temperature anomaly. This is even higher than the range used before, and leads to an equivalent ECS of 3.8 K per CO2 doubling. Together, our results suggest that land ice changes play a key role in the variability of Late Pleistocene temperatures.


2013 ◽  
Vol 13 (12) ◽  
pp. 31479-31526 ◽  
Author(s):  
R. Morales Betancourt ◽  
A. Nenes

Abstract. Aerosol indirect effects in climate models strongly depend on the representation of the aerosol activation process. In this study, we assess the process level differences across activation parameterizations that contribute to droplet number uncertainty by using the adjoints of the Abdul-Razzak and Ghan (2000) and Fountoukis and Nenes (2005) droplet activation parameterizations in the framework of the Community Atmospheric Model version 5.1 (CAM5.1). The adjoint sensitivities of Nd to relevant input parameters are used to: (i) unravel the spatially resolved contribution of aerosol number, mass, and chemical composition to changes in $N_\\mathrm{d}$ between present day and pre-industrial simulations; (ii) identify the key variables responsible for the differences in Nd fields and aerosol indirect effect estimates when different activation schemes are used within the same modeling framework. The sensitivities are computed online at minimal computational cost. Changes in aerosol number and aerosol mass concentrations were found to contribute to Nd differences much more strongly than chemical composition effects. The main sources of discrepancy between the activation parameterization considered were the treatment of the water uptake by coarse mode particles, and the sensitivity of the parameterized Nd accumulation mode aerosol geometric mean diameter. These two factors explain the different predictions of Nd over land and over oceans when these parameterizations are employed. Discrepancies in the sensitivity to aerosol size are responsible for an exaggerated response to aerosol volume changes over heavily polluted regions. Because these regions are collocated with areas of deep clouds their impact on short wave cloud forcing is amplified through liquid water path changes. Application of the adjoint-sensitivities illustrated the importance of primary organic matter emissions in controlling the droplet number concentration changes in several areas. The same framework is also utilized to efficiently explore droplet number uncertainty attributable to hygroscopicity parameter of organic aerosol (primary and secondary). Comparisons between the parameterization-derived sensitivities of droplet number against predictions with detailed numerical simulations of the activation process were performed to validate the physical consistency of the adjoint sensitivities.


2014 ◽  
Vol 14 (9) ◽  
pp. 4809-4826 ◽  
Author(s):  
R. Morales Betancourt ◽  
A. Nenes

Abstract. Aerosol indirect effects in climate models strongly depend on the representation of the aerosol activation process. In this study, we assess the process-level differences across activation parameterizations that contribute to droplet number uncertainty by using the adjoints of the Abdul-Razzak and Ghan (2000) and Fountoukis and Nenes (2005) droplet activation parameterizations in the framework of the Community Atmospheric Model version 5.1 (CAM5.1). The adjoint sensitivities of Nd to relevant input parameters are used to (i) unravel the spatially resolved contribution of aerosol number, mass, and chemical composition to changes in Nd between present-day and pre-industrial simulations and (ii) identify the key variables responsible for the differences in Nd fields and aerosol indirect effect estimates when different activation schemes are used within the same modeling framework. The sensitivities are computed online at minimal computational cost. Changes in aerosol number and aerosol mass concentrations were found to contribute to Nd differences much more strongly than chemical composition effects. The main sources of discrepancy between the activation parameterizations considered were the treatment of the water uptake by coarse mode particles, and the sensitivity of the parameterized Nd accumulation mode aerosol geometric mean diameter. These two factors explain the different predictions of Nd over land and over oceans when these parameterizations are employed. Discrepancies in the sensitivity to aerosol size are responsible for an exaggerated response to aerosol volume changes over heavily polluted regions. Because these regions are collocated with areas of deep clouds, their impact on shortwave cloud forcing is amplified through liquid water path changes. The same framework is also utilized to efficiently explore droplet number uncertainty attributable to hygroscopicity parameter of organic aerosol (primary and secondary). Comparisons between the parameterization-derived sensitivities of droplet number against predictions with detailed numerical simulations of the activation process were performed to validate the physical consistency of the adjoint sensitivities.


2016 ◽  
Vol 16 (15) ◽  
pp. 9785-9804 ◽  
Author(s):  
Matthew Kasoar ◽  
Apostolos Voulgarakis ◽  
Jean-François Lamarque ◽  
Drew T. Shindell ◽  
Nicolas Bellouin ◽  
...  

Abstract. We use the HadGEM3-GA4, CESM1, and GISS ModelE2 climate models to investigate the global and regional aerosol burden, radiative flux, and surface temperature responses to removing anthropogenic sulfur dioxide (SO2) emissions from China. We find that the models differ by up to a factor of 6 in the simulated change in aerosol optical depth (AOD) and shortwave radiative flux over China that results from reduced sulfate aerosol, leading to a large range of magnitudes in the regional and global temperature responses. Two of the three models simulate a near-ubiquitous hemispheric warming due to the regional SO2 removal, with similarities in the local and remote pattern of response, but overall with a substantially different magnitude. The third model simulates almost no significant temperature response. We attribute the discrepancies in the response to a combination of substantial differences in the chemical conversion of SO2 to sulfate, translation of sulfate mass into AOD, cloud radiative interactions, and differences in the radiative forcing efficiency of sulfate aerosol in the models. The model with the strongest response (HadGEM3-GA4) compares best with observations of AOD regionally, however the other two models compare similarly (albeit poorly) and still disagree substantially in their simulated climate response, indicating that total AOD observations are far from sufficient to determine which model response is more plausible. Our results highlight that there remains a large uncertainty in the representation of both aerosol chemistry as well as direct and indirect aerosol radiative effects in current climate models, and reinforces that caution must be applied when interpreting the results of modelling studies of aerosol influences on climate. Model studies that implicate aerosols in climate responses should ideally explore a range of radiative forcing strengths representative of this uncertainty, in addition to thoroughly evaluating the models used against observations.


2021 ◽  
Author(s):  
Junichi Tsutsui

Abstract. Climate model emulators have a crucial role in assessing warming levels of many emission scenarios from probabilistic climate projections, based on new insights into Earth system response to CO2 and other forcing factors. This article describes one such tool, MCE, from model formulation to application examples associated with a recent model intercomparison study. The MCE is based on impulse response functions and parameterized physics of effective radiative forcing and carbon uptake over ocean and land. Perturbed model parameters for probabilistic projections are generated from statistical models and constrained with a Metropolis-Hastings independence sampler. A part of the model parameters associated with CO2-induced warming have a covariance structure, as diagnosed from complex climate models of the Coupled Model Intercomparison Project (CMIP). Although perturbed ensembles can cover the diversity of CMIP models effectively, they need to be constrained toward substantially lower climate sensitivity for the resulting historical warming to agree with the observed trends over recent decades. The model's simplicity and resulting successful calibration imply that a method with less complicated structures and fewer control parameters offers advantages when building reasonable perturbed ensembles in a transparent way. Experimental results for future scenarios show distinct differences between CMIP- and observation-consistent ensembles, suggesting that perturbed ensembles for scenario assessment need to be properly constrained with new insights into forced response over historical periods.


2017 ◽  
Vol 13 (8) ◽  
pp. 1037-1048 ◽  
Author(s):  
Henrik Carlson ◽  
Rodrigo Caballero

Abstract. Recent work in modelling the warm climates of the early Eocene shows that it is possible to obtain a reasonable global match between model surface temperature and proxy reconstructions, but only by using extremely high atmospheric CO2 concentrations or more modest CO2 levels complemented by a reduction in global cloud albedo. Understanding the mix of radiative forcing that gave rise to Eocene warmth has important implications for constraining Earth's climate sensitivity, but progress in this direction is hampered by the lack of direct proxy constraints on cloud properties. Here, we explore the potential for distinguishing among different radiative forcing scenarios via their impact on regional climate changes. We do this by comparing climate model simulations of two end-member scenarios: one in which the climate is warmed entirely by CO2 (which we refer to as the greenhouse gas (GHG) scenario) and another in which it is warmed entirely by reduced cloud albedo (which we refer to as the low CO2–thin clouds or LCTC scenario) . The two simulations have an almost identical global-mean surface temperature and equator-to-pole temperature difference, but the LCTC scenario has  ∼  11 % greater global-mean precipitation than the GHG scenario. The LCTC scenario also has cooler midlatitude continents and warmer oceans than the GHG scenario and a tropical climate which is significantly more El Niño-like. Extremely high warm-season temperatures in the subtropics are mitigated in the LCTC scenario, while cool-season temperatures are lower at all latitudes. These changes appear large enough to motivate further, more detailed study using other climate models and a more realistic set of modelling assumptions.


2021 ◽  
Vol 14 (8) ◽  
pp. 4865-4890
Author(s):  
Peter Uhe ◽  
Daniel Mitchell ◽  
Paul D. Bates ◽  
Nans Addor ◽  
Jeff Neal ◽  
...  

Abstract. Riverine flood hazard is the consequence of meteorological drivers, primarily precipitation, hydrological processes and the interaction of floodwaters with the floodplain landscape. Modeling this can be particularly challenging because of the multiple steps and differing spatial scales involved in the varying processes. As the climate modeling community increases their focus on the risks associated with climate change, it is important to translate the meteorological drivers into relevant hazard estimates. This is especially important for the climate attribution and climate projection communities. Current climate change assessments of flood risk typically neglect key processes, and instead of explicitly modeling flood inundation, they commonly use precipitation or river flow as proxies for flood hazard. This is due to the complexity and uncertainties of model cascades and the computational cost of flood inundation modeling. Here, we lay out a clear methodology for taking meteorological drivers, e.g., from observations or climate models, through to high-resolution (∼90 m) river flooding (fluvial) hazards. Thus, this framework is designed to be an accessible, computationally efficient tool using freely available data to enable greater uptake of this type of modeling. The meteorological inputs (precipitation and air temperature) are transformed through a series of modeling steps to yield, in turn, surface runoff, river flow, and flood inundation. We explore uncertainties at different modeling steps. The flood inundation estimates can then be related to impacts felt at community and household levels to determine exposure and risks from flood events. The approach uses global data sets and thus can be applied anywhere in the world, but we use the Brahmaputra River in Bangladesh as a case study in order to demonstrate the necessary steps in our hazard framework. This framework is designed to be driven by meteorology from observational data sets or climate model output. In this study, only observations are used to drive the models, so climate changes are not assessed. However, by comparing current and future simulated climates, this framework can also be used to assess impacts of climate change.


2020 ◽  
Vol 13 (5) ◽  
pp. 2355-2377
Author(s):  
Vijay S. Mahadevan ◽  
Iulian Grindeanu ◽  
Robert Jacob ◽  
Jason Sarich

Abstract. One of the fundamental factors contributing to the spatiotemporal inaccuracy in climate modeling is the mapping of solution field data between different discretizations and numerical grids used in the coupled component models. The typical climate computational workflow involves evaluation and serialization of the remapping weights during the preprocessing step, which is then consumed by the coupled driver infrastructure during simulation to compute field projections. Tools like Earth System Modeling Framework (ESMF) (Hill et al., 2004) and TempestRemap (Ullrich et al., 2013) offer capability to generate conservative remapping weights, while the Model Coupling Toolkit (MCT) (Larson et al., 2001) that is utilized in many production climate models exposes functionality to make use of the operators to solve the coupled problem. However, such multistep processes present several hurdles in terms of the scientific workflow and impede research productivity. In order to overcome these limitations, we present a fully integrated infrastructure based on the Mesh Oriented datABase (MOAB) (Tautges et al., 2004; Mahadevan et al., 2015) library, which allows for a complete description of the numerical grids and solution data used in each submodel. Through a scalable advancing-front intersection algorithm, the supermesh of the source and target grids are computed, which is then used to assemble the high-order, conservative, and monotonicity-preserving remapping weights between discretization specifications. The Fortran-compatible interfaces in MOAB are utilized to directly link the submodels in the Energy Exascale Earth System Model (E3SM) to enable online remapping strategies in order to simplify the coupled workflow process. We demonstrate the superior computational efficiency of the remapping algorithms in comparison with other state-of-the-science tools and present strong scaling results on large-scale machines for computing remapping weights between the spectral element atmosphere and finite volume discretizations on the polygonal ocean grids.


Sign in / Sign up

Export Citation Format

Share Document