Search NASA⌕ Search

SEARCH · Search NASA

Results for “inter-model variability”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

A Model-Model and Data-Model Comparison for the Early Eocene Hydrological Cycle

A range of proxy observations have recently provided constraints on how Earth's hydrological cycle responded to early Eocene climatic changes. However, comparisons of proxy data to general circulation model (GCM) simulated hydrology are limited and inter-model variability remains poorly characterised. In this work, we undertake an intercomparison of GCM-derived precipitation and P - E distributions within the extended EoMIP ensemble (Eocene Modelling Intercomparison Project; Lunt et al., 2012), which includes previously published early Eocene simulations performed using five GCMs differing in boundary conditions, model structure, and precipitation-relevant parameterisation schemes. We show that an intensified hydrological cycle, manifested in enhanced global precipitation and evaporation rates, is simulated for all Eocene simulations relative to the preindustrial conditions. This is primarily due to elevated atmospheric paleo-CO2, resulting in elevated temperatures, although the effects of differences in paleogeography and ice sheets are also important in some models. For a given CO2 level, globally averaged precipitation rates vary widely between models, largely arising from different simulated surface air temperatures. Models with a similar global sensitivity of precipitation rate to temperature (dP=dT ) display different regional precipitation responses for a given temperature change. Regions that are particularly sensitive to model choice include the South Pacific, tropical Africa, and the Peri-Tethys, which may represent targets for future proxy acquisition. A comparison of early and middle Eocene leaf-fossil-derived precipitation estimates with the GCM output illustrates that GCMs generally underestimate precipitation rates at high latitudes, although a possible seasonal bias of the proxies cannot be excluded. Models which warm these regions, either via elevated CO2 or by varying poorly constrained model parameter values, are most successful in simulating a match with geologic data. Further data from low-latitude regions and better constraints on early Eocene CO2 are now required to discriminate between these model simulations given the large error bars on paleoprecipitation estimates. Given the clear differences between simulated precipitation distributions within the ensemble, our results suggest that paleohydrological data offer an independent means by which to evaluate model skill for warm climates.

Boundary conditions↗
NASA NTRS2018

Multi-model Comparison of the Volcanic Sulfate Deposition from the 1815 Eruption of Mt. Tambora

The eruption of Mt. Tambora in 1815 was the largest volcanic eruption of the past 500 years. The eruption had significant climatic impacts, leading to the 1816 "year without a summer", and remains a valuable event from which to understand the climatic effects of large stratospheric volcanic sulfur dioxide injections. The eruption also resulted in one of the strongest and most easily identifiable volcanic sulfate signals in polar ice cores, which are widely used to reconstruct the timing and atmospheric sulfate loading of past eruptions. As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), five state-of-the-art global aerosol models simulated this eruption. We analyze both simulated background (no Tambora) and volcanic (with Tambora) sulfate deposition to polar regions and compare to ice core records. The models simulate overall similar patterns of background sulfate deposition, although there are differences in regional details and magnitude. However, the volcanic sulfate deposition varies considerably between the models with differences in timing, spatial pattern and magnitude. Mean simulated deposited sulfate on Antarctica ranges from 19 to 264 kgkm-2 and on Greenland from 31 to 194 kgkm-2, as compared to the mean ice-core derived estimates of roughly 50 kgkm-2 for both Greenland and Antarctica. The ratio of the hemispheric atmospheric sulfate aerosol burden after the eruption to the average ice sheet deposited sulfate varies between models by up to a factor of 15. Sources of this inter-model variability include differences in both the formation and the transport of sulfate aerosol. Our results suggest that deriving relationships between sulfate deposited on ice sheets and atmospheric sulfate burdens from model simulations may be associated with greater uncertainties than previously thought.

Marshall, Lauren↗
NASA NTRS2019

Large-Scale Transport into the Arctic: The Roles of the Midlatitude Jet and the Hadley Cell

Transport from the Northern Hemisphere (NH) midlatitudes to the Arctic plays a crucial role in determining the abundance of trace gases and aerosols that are important to Arctic climate via impacts on radiation and chemistry. Here we examine this transport using an idealized tracer with a fixed lifetime and predominantly midlatitude land-based sources in models participating in the Chemistry Climate Model Initiative (CCMI). We show that there is a 25%-45% difference in the Arctic concentrations of this tracer among the models. This spread is correlated with the spread in the location of the Pacific jet, as well as the spread in the location of the Hadley Cell (HC) edge, which varies consistently with jet latitude. Our results suggest that it is likely that the HC-related zonal-mean meridional transport rather than the jet-related eddy mixing is the major contributor to the inter-model spread in the transport of land-based tracers into the Arctic. Specifically, in models with a more northern jet, the HC generally extends further north and the tracer source region is mostly covered by surface southward flow associated with the lower branch of the HC, resulting in less efficient transport poleward to the Arctic. During boreal summer, there are poleward biases in jet location in free-running models, and these models likely underestimate the rate of transport into the Arctic. Models using specified dynamics do not have biases in the jet location, but do have biases in the surface meridional flow, which may result in differences in transport into the Arctic. In addition to the land-based tracer, the midlatitude-to-Arctic transport is further examined by another idealized tracer with zonally uniform sources. With equal sources from both land and ocean, the inter-model spread of this zonally uniform tracer is more related to variations in parameterized convection over oceans rather than variations in HC extent, particularly during boreal winter. This suggests that transport of land-based and oceanic tracers or aerosols towards the Arctic differs in pathways and therefore their corresponding inter-model variabilities result from different physical processes.

chemistry climate models↗
NASA NTRS2023

Global and Zonal-Mean Hydrological Response to Early Eocene Warmth

Earth's hydrological cycle is expected to intensify in response to global warming, with a “wet-gets-wetter, dry-gets-drier” response anticipated over the ocean. Subtropical regions (∼15°–30°N/S) are predicted to become drier, yet proxy evidence from past warm climates suggests these regions may be characterized by wetter conditions. Here we use an integrated data-modeling approach to reconstruct global and zonal-mean rainfall patterns during the early Eocene (∼56–48 million years ago). The Deep-Time Model Intercomparison Project (DeepMIP) model ensemble indicates that the mid- (30°–60°N/S) and high-latitudes (>60°N/S) are characterized by a thermodynamically dominated hydrological response to warming and overall wetter conditions. The tropical band (0°–15°N/S) is also characterized by wetter conditions, with several DeepMIP models simulating narrowing of the Inter-Tropical Convergence Zone. However, the latter is not evident from the proxy data. The subtropics are characterized by negative precipitation-evaporation anomalies (i.e., drier conditions) in the DeepMIP models, but there is surprisingly large inter-model variability in mean annual precipitation (MAP). Intriguingly, we find that models with weaker meridional temperature gradients (e.g., CESM, GFDL) are characterized by a reduction in subtropical moisture divergence, leading to an increase in MAP. These model simulations agree more closely with our new proxy-derived precipitation reconstructions and other key climate metrics and imply that the early Eocene was characterized by reduced subtropical moisture divergence. If the meridional temperature gradient was even weaker than suggested by those DeepMIP models, circulation-induced changes may have outcompeted thermodynamic changes, leading to wetter subtropics. This highlights the importance of accurately reconstructing zonal temperature gradients when reconstructing past rainfall patterns.

Margot J. Cramwinckel↗
NASA NTRS2018

Weak Hydrological Sensitivity to Temperature Change over Land, Independent of Climate Forcing

We present the global and regional hydrological sensitivity (HS) to surface temperature changes, for perturbations to CO2, CH4, sulfate and black carbon concentrations, and solar irradiance. Based on results from ten climate models, we show how modeled global mean precipitation increases by 2-3% per kelvin of global mean surface warming, independent of driver, when the effects of rapid adjustments are removed. Previously reported differences in response between drivers are therefore mainly ascribable to rapid atmospheric adjustment processes. All models show a sharp contrast in behavior over land and over ocean, with a strong surface temperature-driven (slow) ocean HS of 3-5%/K, while the slow land HS is only 0-2%/K. Separating the response into convective and large-scale cloud processes, we find larger inter-model differences, in particular over land regions. Large-scale precipitation changes are most relevant at high latitudes, while the equatorial HS is dominated by convective precipitation changes. Black carbon stands out as the driver with the largest inter-model slow HS variability, and also the strongest contrast between a weak land and strong sea response. We identify a particular need for model investigations and observational constraints on convective precipitation in the Arctic, and large-scale precipitation around the Equator.

Samset, B. H.↗
NASA NTRSDate not supplied

On the Nature of the Arctic’s Positive Lapse-rate Feedback: Cause or Symptom of Arctic Amplification

Radiative energy fluxes exiting the atmosphere towards the surface and to space are sensitive to the vertical structure of temperature. Under anthropogenic forcing, this sensitivity gives rise to the lapse-rate feedback. Shown to be negative in the tropics and positive in the Arctic, studies argue that the lapse-rate feedback is the primary cause of Arctic Amplification. While mechanistically the negative tropical lapse-rate feedback is understood, the mechanics of the positive Arctic lapse-rate feedback are less clear. Previous arguments state that the Arctic’s stable stratification is the origin of its positive lapse-rate feedback. We present results using CMIP5 and CMIP6 model output that arguing that the magnitude, spatial variability, seasonality, and inter-model spread in the Arctic lapse-rate feedback is controlled by surface properties and atmosphere-ocean-sea ice energy exchanges, not the degree of stable stratification. We argue that the Arctic’s positive lapse-rate feedback is a seasonal and regional phenomenon that manifests from local surface characteristics and surface albedo, ice insulation, and thermal inertia feedbacks that cause surface warming to outpace atmospheric warming in sea ice-retreat regions. We view the Arctic’s positive lapse-rate feedback as a symptom rather than the primary cause of Arctic Amplification.

Patrick Charles Taylor↗
NASA NTRSDate not supplied

On the Nature of the Arctic’s Positive Lapse Rate Feedback Cause or Symptom of Arctic Amplification

Radiative energy fluxes exiting the atmosphere towards the surface and to space are sensitive to the vertical structure of temperature. Under anthropogenic forcing, this sensitivity gives rise to the lapse-rate feedback. Shown to be negative in the tropics and positive in the Arctic, studies argue that the lapse-rate feedback is the primary cause of Arctic Amplification. While mechanistically the negative tropical lapse-rate feedback is understood, the mechanics of the positive Arctic lapse-rate feedback are less clear. Previous arguments state that the Arctic’s stable stratification is the origin of its positive lapse-rate feedback. We present results using CMIP5 and CMIP6 model output that arguing that the magnitude, spatial variability, seasonality, and inter-model spread in the Arctic lapse-rate feedback is controlled by surface properties and atmosphere-ocean-sea ice energy exchanges, not the degree of stable stratification. We argue that the Arctic’s positive lapse-rate feedback is a seasonal and regional phenomenon that manifests from local surface characteristics and surface albedo, ice insulation, and thermal inertia feedbacks that cause surface warming to outpace atmospheric warming in sea ice-retreat regions. We view the Arctic’s positive lapse-rate feedback as a symptom rather than the primary cause of Arctic Amplification.

Patrick C Taylor↗
NASA NTRS2022

The TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI). Part II: Moist Cases - The Two Waterworlds

To identify promising exoplanets for atmospheric characterization and to make the best use of observational data, a thorough understanding of their atmospheres is needed. 3D general circulation models (GCMs) are one of the most comprehensive tools available for this task and will be used to interpret observations of temperate rocky exoplanets. Due to parameterization choices made in GCMs, they can produce different results, even for the same planet. Employing four widely-used exoplanetary GCMs - ExoCAM, LMD-G, ROCKE-3D and the UM - we continue the TRAPPIST-1 Habitable Atmosphere Intercomparison by modeling aquaplanet climates of TRAPPIST-le with a moist atmosphere dominated by either nitrogen or carbon dioxide. Although the GCMs disagree on the details of the simulated regimes, they all predict a temperate climate with neither of the two cases pushed out of the habitable state. Nevertheless, the inter-model spread in the global mean surface temperature is non-negligible: 14 K and 24 K in the nitrogen and carbon dioxide dominated case, respectively. We find substantial inter-model differences in moist variables, with the smallest amount of clouds in LMDGeneric and the largest in ROCKE-3D. ExoCAM predicts the warmest climate for both cases and thus has the highest water vapor content and the largest amount and variability of cloud condensate. The UM tends to produce colder conditions, especially in the nitrogen-dominated case due to a strong negative cloud radiative effect on the day side of TRAPPIST-le. Our study highlights various biases of GCMs and emphasizes the importance of not relying solely on one model to understand exoplanet climates.

TRAPPIST↗
NASA NTRS2016

What Controls the Vertical Distribution of Aerosol? Relationships Between Process Sensitivity in HadGEM3-UKCA and Inter-Model Variation from AeroCom Phase II

The vertical profile of aerosol is important for its radiative effects, but weakly constrained by observations on the global scale, and highly variable among different models. To investigate the controlling factors in one particular model, we investigate the effects of individual processes in HadGEM3-UKCA and compare the resulting diversity of aerosol vertical profiles with the inter-model diversity from the AeroCom Phase II control experiment. In this way we show that (in this model at least) the vertical profile is controlled by a relatively small number of processes, although these vary among aerosol components and particle sizes. We also show that sufficiently coarse variations in these processes can produce a similar diversity to that among different models in terms of the global-mean profile and, to a lesser extent, the zonal-mean vertical position. However, there are features of certain models' profiles that cannot be reproduced, suggesting the influence of further structural differences between models. In HadGEM3-UKCA, convective transport is found to be very important in controlling the vertical profile of all aerosol components by mass. In-cloud scavenging is very important for all except mineral dust. Growth by condensation is important for sulfate and carbonaceous aerosol (along with aqueous oxidation for the former and ageing by soluble material for the latter). The vertical extent of biomass-burning emissions into the free troposphere is also important for the profile of carbonaceous aerosol. Boundary-layer mixing plays a dominant role for sea salt and mineral dust, which are emitted only from the surface. Dry deposition and below-cloud scavenging are important for the profile of mineral dust only. In this model, the microphysical processes of nucleation, condensation and coagulation dominate the vertical profile of the smallest particles by number (e.g. total CN >3 nm), while the profiles of larger particles (e.g. CN>100 nm) are controlled by the same processes as the component mass profiles, plus the size distribution of primary emissions. We also show that the processes that affect the AOD-normalised radiative forcing in the model are predominantly those that affect the vertical mass distribution, in particular convective transport, in-cloud scavenging, aqueous oxidation, ageing and the vertical extent of biomass-burning emissions.

dust↗
NASA NTRS2021

Placing the East-West North American Aridity Gradient in a Multi-Century Context

Instrumental records indicate a century-long trend towards drying over western North America and wetting over eastern North America. A continuation of these trends into the future would have significant hydroclimatic and socioeconomic consequences in both the semi-arid Southwest and humid East. Using tree-ring reconstructions and hydrologic simulations of summer soil moisture, we evaluate and contextualize the modern summer aridity gradient within its natural range of variability established over the past 600 years and evaluate the effects of observed and anthropogenic precipitation, temperature, and humidity trends. The 2001–2020 positive (wet east-dry west) aridity gradient was larger than any 20-year period since 1400 CE, preceded by the most negative (wet west-dry east) aridity gradient during 1976-1995, leading to a strong multi-decade reversal in aridity gradient anomalies that was rivaled only by a similar event in the late-16th century. The 2001-2020 aridity gradient was dominated by long-term summer precipitation increases in the Midwest and Northeast, with smaller contributions from more warming in the West than the East and spring precipitation decreases in the Southwest. Multi-model mean climate simulations from CMIP6 experiments suggest anthropogenic climate trends should not have strongly affected the aridity gradient thus far. However, there is high uncertainty due to inter-model disagreement on anthropogenic precipitation trends. The recent strengthening of the observed aridity gradient, its increasing dependence on precipitation variability, and disagreement in modeled anthropogenic precipitation trends reveal significant uncertainties in how water resource availability will change across North America in the coming decades.

drought↗
NASA NTRS2012

Upper-Ocean Heat Balance Processes and the Walker Circulation in CMIP5 Model Projections

Considerable uncertainty remains as to the importance of mechanisms governing decadal and longer variability of the Walker Circulation, its connection to the tropical climate system, and prospects for tropical climate change in the face of anthropogenic forcing. Most contemporary climate models suggest that in response to elevated CO2 and a warmer but more stratified atmosphere, the required upward mass flux in tropical convection will diminish along with the Walker component of the tropical mean circulation as well. Alternatively, there is also evidence to suggest that the shoaling and increased vertical stratification of the thermocline in the eastern Pacific will enable a muted SST increase there-- preserving or even enhancing some of the dynamical forcing for the Walker cell flow. Over the past decade there have been observational indications of an acceleration in near-surface easterlies, a strengthened Pacific zonal SST gradient, and globally-teleconnected dislocations in precipitation. But is this evidence in support of an ocean dynamical thermostat process posited to accompany anthropogenic forcing, or just residual decadal fluctuations associated with variations in warm and cold ENSO events and other stochastic forcing? From a modeling perspective we try to make headway on this question by examining zonal variations in surface energy fluxes and dynamics governing tropical upper ocean heat content evolution in the WCRP CMIP5 model projections. There is some diversity among model simulations; for example, the CCSM4 indicates net ocean warming over the IndoPacific region while the CSIRO model concentrates separate warming responses over the central Pacific and Indian Ocean regions. The models, as with observations, demonstrate strong local coupling between variations in column water vapor, downward surface longwave radiation and SST; but the spatial patterns of changes in the sign of this relationship differ among models and, for models as a whole, with observations. Our analysis focuses initially on probing the inter-model differences in energy fluxes / transports and Walker Circulation response to forcing. We then attempt to identify statistically the El Nino- / La Nina-related ocean heat content variability unique to each model and regress out the associated energy flux, ocean heat transport and Walker response on these shorter time scales for comparison to that of the anthropogenic signals.

Robertson, F. R.↗
NASA NTRS2001

Comparing the Degree of Land-Atmosphere Interaction in Four Atmospheric General Circulation Models

Land-atmosphere feedback, by which (for example) precipitation-induced moisture anomalies at the land surface affect the overlying atmosphere and thereby the subsequent generation of precipitation, has been examined and quantified with many atmospheric general circulation models (AGCMs). Generally missing from such studies, however, is an indication of the extent to which the simulated feedback strength is model dependent. Four modeling groups have recently performed a highly controlled numerical experiment that allows an objective inter-model comparison of land-atmosphere feedback strength. The experiment essentially consists of an ensemble of simulations in which each member simulation artificially maintains the same time series of surface prognostic variables. Differences in atmospheric behavior between the ensemble members then indicates the degree to which the state of the land surface controls atmospheric processes in that model. A comparison of the four sets of experimental results shows that feedback strength does indeed vary significantly between the AGCMs.

Koster, Randal D.↗
NASA NTRS2023

An Observation-Based, Reduced-Form Model for Oxidation in the Remote Marine Troposphere

The hydroxyl radical (OH) fuels atmospheric chemical cycling as the main sink for methane and a driver of the formation and loss of many air pollutants, but direct OH observations are sparse. We develop and evaluate an observation-based proxy for short term, spatial variations in OH (Proxy OH ) in the remote marine troposphere using unprecedented and comprehensive measurements from the NASA Atmospheric Tomography (ATom) airborne campaign. Proxy OH is a reduced form of the OH steady-state equation representing the dominant OH production and loss pathways in the remote marine troposphere, according to box model simulations of OH constrained with ATom observations. Proxy OH comprises only eight variables that are generally observed by routine ground- or satellite-based instruments. ProxyOH scales linearly with in situ [OH] spatial variations along the ATom flight tracks (median r 2 = 0.90, interquartile range = 0.80 – 0.94 across 2 km altitude by 20° latitudinal regions). We deconstruct spatial variations in Proxy OH as a first-order approximation of the sensitivity of OH variations to individual terms. Two terms modulate within-region Proxy OH variations—water vapor (H 2 O) and, to a lesser extent, nitric oxide (NO). This implies that a limited set of observations could offer a novel avenue for observation-based mapping of OH spatial variations over much of the remote marine troposphere. Both H 2 O and NO are expected to change with climate, while NO also varies strongly with human activities. We also illustrate the utility of Proxy OH as a process-based approach for evaluating inter-model differences in remote marine tropospheric OH.

Colleen B. Baublitz↗
NASA NTRS2021

AeroCom Phase III Multi-Model Evaluation of the Aerosol Life Cycle and Optical Properties Using Ground and Space-Based Remote Sensing as Well as Surface In Situ Observations

Within the framework of the AeroCom (Aerosol Comparisons between Observations and Models) initiative, the state-of-the-art modelling of aerosol optical properties is assessed from 14 global models participating in the phase III control experiment (AP3). The models are similar to CMIP6/AerChemMIP Earth System Models (ESMs) and provide a robust multi-model ensemble. Inter-model spread of aerosol species lifetimes and emissions appears to be similar to that of mass extinction coefficients (MECs), suggesting that aerosol optical depth (AOD) uncertainties are associated with a broad spectrum of parameterised aerosol processes. Total AOD is approximately the same as in AeroCom phase I (AP1) simulations. However, we find a 50 % decrease in the optical depth (OD) of black carbon (BC), attributable to a combination of decreased emissions and lifetimes. Relative contributions from sea salt (SS) and dust (DU) have shifted from being approximately equal in AP1 to SS contributing about 2∕3 of the natural AOD in AP3. This shift is linked with a decrease in DU mass burden, a lower DU MEC, and a slight decrease in DU lifetime, suggesting coarser DU particle sizes in AP3 compared to AP1. Relative to observations, the AP3 ensemble median and most of the participating models underestimate all aerosol optical properties investigated, that is, total AOD as well as fine and coarse AOD (AODf, AODc), Ångström exponent (AE), dry surface scattering (SCdry), and absorption (ACdry) coefficients. Compared to AERONET, the models underestimate total AOD by ca. 21 % ± 20 % (as inferred from the ensemble median and interquartile range). Against satellite data, the ensemble AOD biases range from −37 % (MODIS-Terra) to −16 % (MERGED-FMI, a multi-satellite AOD product), which we explain by differences between individual satellites and AERONET measurements themselves. Correlation coefficients (R) between model and observation AOD records are generally high (R>0.75), suggesting that the models are capable of capturing spatio-temporal variations in AOD. We find a much larger underestimate in coarse AODc (∼ −45 % ± 25 %) than in fine AODf (∼ −15 % ± 25 %) with slightly increased inter-model spread compared to total AOD. These results indicate problems in the modelling of DU and SS. The AODc bias is likely due to missing DU over continental land masses (particularly over the United States, SE Asia, and S. America), while marine AERONET sites and the AATSR SU satellite data suggest more moderate oceanic biases in AODc. Column AEs are underestimated by about 10 % ± 16 %. For situations in which measurements show AE > 2, models underestimate AERONET AE by ca. 35 %. In contrast, all models (but one) exhibit large overestimates in AE when coarse aerosol dominates (bias ca. +140 % if observed AE < 0.5). Simulated AE does not span the observed AE variability. These results indicate that models overestimate particle size (or underestimate the fine-mode fraction) for fine-dominated aerosol and underestimate size (or overestimate the fine-mode fraction) for coarse-dominated aerosol. This must have implications for lifetime, water uptake, scattering enhancement, and the aerosol radiative effect, which we can not quantify at this moment. Comparison against Global Atmosphere Watch (GAW) in situ data results in mean bias and inter-model variations of −35 % ± 25 % and −20 % ± 18 % for SCdry and ACdry, respectively. The larger underestimate of SCdry than ACdry suggests the models will simulate an aerosol single scattering albedo that is too low. The larger underestimate of SCdry than ambient air AOD is consistent with recent findings that models overestimate scattering enhancement due to hygroscopic growth. The broadly consistent negative bias in AOD and surface scattering suggests an underestimate of aerosol radiative effects in current global aerosol models. Considerable inter-model diversity in the simulated optical properties is often found in regions that are, unfortunately, not or only sparsely covered by ground-based observations. This includes, for instance, the Sahara, Amazonia, central Australia, and the South Pacific. This highlights the need for a better site coverage in the observations, which would enable us to better assess the models, but also the performance of satellite products in these regions. Using fine-mode AOD as a proxy for present-day aerosol forcing estimates, our results suggest that models underestimate aerosol forcing by ca. −15 %, however, with a considerably large interquartile range, suggesting a spread between −35 % and +10 %.

space-based remote sensing↗
NASA NTRS2012

The Agricultural Model Intercomparison and Improvement Project (AgMIP): Protocols and Pilot Studies

The Agricultural Model Intercomparison and Improvement Project (AgMIP) is a major international effort linking the climate, crop, and economic modeling communities with cutting-edge information technology to produce improved crop and economic models and the next generation of climate impact projections for the agricultural sector. The goals of AgMIP are to improve substantially the characterization of world food security due to climate change and to enhance adaptation capacity in both developing and developed countries. Analyses of the agricultural impacts of climate variability and change require a transdisciplinary effort to consistently link state-of-the-art climate scenarios to crop and economic models. Crop model outputs are aggregated as inputs to regional and global economic models to determine regional vulnerabilities, changes in comparative advantage, price effects, and potential adaptation strategies in the agricultural sector. Climate, Crop Modeling, Economics, and Information Technology Team Protocols are presented to guide coordinated climate, crop modeling, economics, and information technology research activities around the world, along with AgMIP Cross-Cutting Themes that address uncertainty, aggregation and scaling, and the development of Representative Agricultural Pathways (RAPs) to enable testing of climate change adaptations in the context of other regional and global trends. The organization of research activities by geographic region and specific crops is described, along with project milestones. Pilot results demonstrate AgMIP's role in assessing climate impacts with explicit representation of uncertainties in climate scenarios and simulations using crop and economic models. An intercomparison of wheat model simulations near Obregón, Mexico reveals inter-model differences in yield sensitivity to [CO2] with model uncertainty holding approximately steady as concentrations rise, while uncertainty related to choice of crop model increases with rising temperatures. Wheat model simulations with midcentury climate scenarios project a slight decline in absolute yields that is more sensitive to selection of crop model than to global climate model, emissions scenario, or climate scenario downscaling method. A comparison of regional and national-scale economic simulations finds a large sensitivity of projected yield changes to the simulations' resolved scales. Finally, a global economic model intercomparison example demonstrates that improvements in the understanding of agriculture futures arise from integration of the range of uncertainty in crop, climate, and economic modeling results in multi-model assessments.

Rosenzweig, C.↗
NASA NTRS2015

Using Results from the PlioMIP Ensemble to Investigate the Greenland Ice Sheet During the Mid-Pliocene Warm Period

During an interval of the Late Pliocene, referred to here as the mid-Pliocene Warm Period (mPWP; 3.264 to 3.025 million years ago), global mean temperature was similar to that predicted for the end of this century, and atmospheric carbon dioxide concentrations were higher than preindustrial levels. Sea level was also higher than today, implying a significant reduction in the extent of the ice sheets. Thus, the mPWP provides a natural laboratory in which to investigate the long-term response of the Earth's ice sheets and sea level in a warmer-than-present-day world. At present, our understanding of the Greenland ice sheet during the mPWP is generally based upon predictions using single climate and ice sheet models. Therefore, it is essential that the model dependency of these results is assessed. The Pliocene Model Intercomparison Project (PlioMIP) has brought together nine international modeling groups to simulate the warm climate of the Pliocene. Here we use the climatological fields derived from the results of the 15 PlioMIP climate models to force an offline ice sheet model. We show that mPWP ice sheet reconstructions are highly dependent upon the forcing climatology used, with Greenland reconstructions ranging from an ice-free state to a near modern ice sheet. An analysis of the surface albedo variability between the climate models over Greenland offers insights into the drivers of inter-model differences. As we demonstrate that the climate model dependency of our results is high, we highlight the necessity of data-based constraints of ice extent in developing our understanding of the mPWP Greenland ice sheet.

sea level↗
NASA NTRS2020

Attribution of Chemistry-Climate Model Initiative (CCMI) Ozone Radiative Flux Bias from Satellites

The top-of-atmosphere (TOA) outgoing longwave flux over the 9.6-μm ozone band is a fundamental quantity for understanding chemistry-climate coupling. However, observed TOA fluxes are hard to estimate as they exhibit considerable variability in space and time that depend on the distributions of clouds, ozone (O3), water vapor (H2O), air temperature (Ta), and surface temperature (Ts). Benchmarking present day fluxes and quantifying the relative influence of their drivers is the first step for estimating climate feedbacks from ozone radiative forcing and predicting radiative forcing evolution. To that end, we constructed observational instantaneous radiative kernels (IRKs) under clear-sky conditions, representing the sensitivities of the TOA flux in the 9.6-μm ozone band to the vertical distribution of geophysical variables, including O3, H2O, Ta, and Ts based upon the Aura Tropospheric Emission Spectrometer (TES) measurements. Applying these kernels to present-day simulations from the Chemistry-Climate Model Initiative (CCMI) project as compared to a 2006 reanalysis assimilating satellite observations, we show that the models have large differences in TOA flux, attributable to different geophysical variables. In particular, model simulations continue to diverge from observations in the tropics, as reported in previous studies of the Atmospheric Chemistry Climate Model Inter-comparison Project (ACCMIP) simulations. The principal culprits are tropical mid and upper tropospheric ozone followed by tropical lower tropospheric H2O. Five models out of the eight studied here have TOA flux biases exceeding 100 mWm-2 attributable to tropospheric ozone bias. Another set of five models have flux biases over 50 mWm-2 due to H2O. On the other hand, Ta radiative bias is negligible in all models (no more than 30 mWm-2). We found that AM3 and CMAM have the lowest TOA flux biases globally but are a result of cancellation of opposite biases due to difference processes. Overall, the multi-model ensemble mean bias is –133±98 mWm-2, indicating that they are too atmospherically opaque due to trapping too much radiation in the atmosphere by overestimated tropical tropospheric O3 and H2O. Having too much O3 and H2O in the troposphere would have different impacts on the sensitivity of TOA flux to O3 and these competing effects add more uncertainties on the ozone radiative forcing. We find that the inter-model TOA outgoing longwave radiation (OLR) difference is well anti-correlated with their ozone band flux bias. This suggests that there is significant radiative compensation in the calculation of model outgoing longwave radiation.

Aura Tropospheric Emission Spectrometer (TES) meas↗
NASA NTRS2023

Effects of Forcing Differences and Initial Conditions on Inter-Model Agreement in the VolMIP Volc-Pinatubo-Full Experiment

This paper provides initial results from a multi-model ensemble analysis based on the volc-pinatubo-full experiment performed within the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP) as part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6). The volc-pinatubo-full experiment is based on an ensemble of volcanic forcing-only climate simulations with the same volcanic aerosol dataset across the participating models (the 1991–1993 Pinatubo period from the CMIP6-GloSSAC dataset). The simulations are conducted within an idealized experimental design where initial states are sampled consistently across models from the CMIP6-piControl simulation providing unperturbed preindustrial background conditions. The multi-model ensemble includes output from an initial set of six participating Earth system models (CanESM5, GISS-E2.1-G, IPSL-CM6A-LR, MIROC-E2SL, MPI-ESM1.2-LR and UKESM1). The results show overall good agreement between the different models on the global and hemispheric scales concerning the surface climate responses, thus demonstrating the overall effectiveness of VolMIP's experimental design. However, small yet significant inter-model discrepancies are found in radiative fluxes, especially in the tropics, that preliminary analyses link with minor differences in forcing implementation; model physics, notably aerosol–radiation interactions; the simulation and sampling of El Niño–Southern Oscillation (ENSO); and, possibly, the simulation of climate feedbacks operating in the tropics. We discuss the volc-pinatubo-full protocol and highlight the advantages of volcanic forcing experiments defined within a carefully designed protocol with respect to emerging modelling approaches based on large ensemble transient simulations. We identify how the VolMIP strategy could be improved in future phases of the initiative to ensure a cleaner sampling protocol with greater focus on the evolving state of ENSO in the pre-eruption period.

metadata↗