AICCA: AI-Driven Cloud Classification Atlas

Project Summary:

Clouds represent the single largest uncertainty in projections of future climate change. Their multi-scale nature means that even state-of-the-art numerical climate simulations cannot reliably project how their distribution, frequency, and properties will alter under high-CO2 conditions. Observational data may therefore be key to understanding cloud behavior, but the information contained in decades and petabytes of multispectral cloud imagery from satellite instruments has received only limited use. Recent advances in artificial intelligence (AI) now allow tapping this underutilized resource to improve understanding of cloud dynamics and feedbacks. A core requirement for analysis of cloud images is dimension reduction: dividing a complex high-resolution spatial field into a more limited set of categories. Most attempts at automated cloud classification to date have involved supervised learning of a limited number of human-defined cloud categories, but spatial patterns of cloud textures vary widely, and our goal is to produce novel classifications that provide new insight.

Clouds project therefore proposes, Rotation-Invariant Cloud Clustering (RICC), instead to leverage modern unsupervised deep learning methods to identify robust and meaningful clusters of cloud patterns, and to apply the resulting clustering method to 22 years of data from the Moderate Resolution Imaging Spectroradiometer (MODIS) instruments on NASA’s Aqua and Terra satellites. The resulting AI-driven dataset of cloud classifications, AICCA: AI-driven Cloud Classification Atlas, will deliver in a compact form (10s of GB of class labels, with high spatial and temporal resolution) information currently accessible only as 100s of TB of multi-spectral images.

AICCA will enable the data-driven diagnosis of patterns of cloud organization, provide insight into their evolution on timescales of hours to decades, and contribute to a democratization of climate research by facilitating access to core data.

For more information on the project:

https://takglobus.github.io/AICCA-explorer/

https://github.com/RDCEP/clouds

Paper Link: https://www.mdpi.com/2072-4292/14/22/5690

Paper Summary:

Clouds play an important role in the Earth’s energy budget and their behavior is one of the largest uncertainties in future climate projections. Satellite observations should help in understanding cloud responses, but decades and petabytes of multispectral cloud imagery have to date received only limited use. This study describes a new analysis approach that reduces the dimensionality of satellite cloud observations by grouping them via a novel automated, unsupervised cloud classification technique based on a convolutional autoencoder, an artificial intelligence (AI) method good at identifying patterns in spatial data.

Our technique combines a rotation-invariant autoencoder and hierarchical agglomerative clustering to generate cloud clusters that capture meaningful distinctions among cloud textures, using only raw multispectral imagery as input. Cloud classes are therefore defined based on spectral properties and spatial textures without reliance on location, time/season, derived physical properties, or pre-designated class definitions. We use this approach to generate a unique new cloud dataset, the AI-driven cloud classification atlas (AICCA), which clusters 22 years of ocean images from the Moderate Resolution Imaging Spectroradiometer (MODIS) on NASA’s Aqua and Terra instruments---198 million patches, each roughly 100 km x 100 km (128 x 128 pixels)---into 42 AI-generated cloud classes, a number determined via a newly-developed stability protocol that we use to maximize richness of information while ensuring stable groupings of patches. AICCA thereby translates 801 TB of satellite images into  54.2 GB of class labels and cloud top and optical properties, a reduction by a factor of 15000. The 42 AICCA classes produce meaningful spatio-temporal and physical distinctions and capture a greater variety of cloud types than do the 9 International Satellite Cloud Climatology Project (ISCCP) categories: for example, multiple textures in the stratocumulus decks along the West coasts of North and South America. We conclude that our methodology has explanatory power, capturing % in that it captures regionally unique cloud classes and providing rich but tractable information for global analysis. AICCA delivers the information from multi-spectral images in a compact form, enables data-driven diagnosis of patterns of cloud organization, provides insight into cloud evolution on timescales of hours to decades, and helps democratize climate research by facilitating access to core data.

Endogenous technical change and climate change

One of the central problems in understanding climate change is understanding how it will impact the economy. Most current models assume that climate change will reduce usable output. The resulting effects on long term growth under such an assumption are small. If usable output is reduced, there will likely be some reduction in savings and hence capital available in future periods. Nevertheless, the core drivers of growth such as technological development remain. As a result, under most models, the world is an order of magnitude or more richer in several hundred years even if temperature change is extreme and the damages from climate change potentially intolerable. Some models consider the possibility that climate change will directly destroy capital (i.e., such as through destruction of buildings in storms). While growth is somewhat slower in this formulation, the basic result is unchanged.

In this project we consider a simple modification to the standard damage function so that climate change can affect growth. We take a canonical endogenous growth model and modify it so that climate change has an impact on the sectors in the economy that contribute to growth. For example, if an economy is modeled as having an inventing sector and a manufacturing sector, we allow climate damages to affect the output of both sectors. With this simple reformulation, the effects of climate change can be dramatically different than in standard models. Over reasonably long periods, the growth effects dominate, showing that studying the impacts of climate change on growth should be of central concern.

Our goal is to consider this effect in all of the canonical endogenous growth models to see how robust the effect is across model choices. We have solved two models so far and likely will want to consider another two or three. In addition, we are simulating the results of one or more models to understand the likely size of the effects under differing parameter choices.

People:

David WeisbachMishung Ahn

 

Climate variability: statistics and observation based simulations

Characterizing potential future changes in temperature variability across frequencies: 

The impacts of climate change on human societies may arise not just from changes in climate means but from changes in climate variability. Many agricultural crops, for example, are highly sensitive to even brief periods of stress temperatures, particularly at certain times of the growing cycle. Crop yields can be strongly affected by changes in temperature variability even in the absence of a change in mean. Understanding potential changes in temperature variability has therefore been a research priority in climate science.

Impacts of climate variability depend on not only the magnitude but also the frequency of variations: day-to-day temperature fluctuations have different consequences than year-by-year differences. We therefore study variability using spectral methods that allow distinguishing timescales of fluctuations. Current projects address climate variability in a number of ways

  • comparing variability changes in present-day or preindustrial and and future equilibrated climates
  • developing methodologies to characterize variability changes in transient climates
  • evaluating the effect of model spatial resolution on temperature variability

In our study of equilibrated climates, we compare pre-industrial and future climate scenarios in two different climate models, CCSM3 and GISS-E2-R, with millennial runs so that each climate state is stationary. Following techniques developed in Leeds et al (2015),  we compute integrated variability in distinct frequency bands and show changes over time as the future/present ratio. spectral ratios to explore the temperature variability changes in increased radiative forcing based on climate models. Within CCSM3, we also compare variability changes when warming is driven by two different forcing agents (CO2 and solar radiation).

Figure 1: Changes in summertime (JJA) temperature variability in three frequency bands from pre-industrial to future climates states, for Two models: (LEft) CCSM3 (289 ppm → 1400 ppm); (Right) GISS-E2-R (285 ppm → 1140 ppm). Variabili…

Figure 1: Changes in summertime (JJA) temperature variability in three frequency bands from pre-industrial to future climates states, for Two models: (LEft) CCSM3 (289 ppm → 1400 ppm); (Right) GISS-E2-R (285 ppm → 1140 ppm). Variability changes are shown as ratios of the standard deviation, Values below 1 mean temperature variability decreases. GISS-E2-R model outputs are regrided to T31 resolution as CCSM3 using area-conserving remapping. Patterns In both models are similar: variability tends to increase slightly over land and decrease over the ocean. (In wintertime, variability tends to decrease everywhere other than tropical land.)

Characterizing and emulating variability in a transient climate:

Characterizing climate variability is easier in stationary conditions, but the Earth for the foreseeable future will be in a changing or "transient" state, in which temperatures are evolving in response to changed CO2 concentrations. Most archived climate model output also simulates transient states. We therefore seek to statistically describe how variability changes in a transient climate. This statistical exploration also serves our overall research goal is also to use model information to produce "data-driven simulations" for use in impacts assessments.  GCM output should not be used directly in impacts assessments, because GCMs do not fully reproduce present-day temperature distributions. Instead we develop methods of generating simulations of future temperatures that combine observational records with GCM projections of changes in variability (covariances). (See Figure 2 for cartoon, also description of emulation research.) 

Figure 2: Cartoon illustration comparing strategies for simulating temperatures that combine information from a model and the observational record. Columns compare simple bias correction (left), the Delta method (center), and our proposed method (ri…

Figure 2: Cartoon illustration comparing strategies for simulating temperatures that combine information from a model and the observational record. Columns compare simple bias correction (left), the Delta method (center), and our proposed method (right). Top row, the model predicts changes in mean temperature but no changes in variability; in this case, our proposed method is equivalent to the Delta method. Bottom row, the model predicts changes in both mean and covariance. Simple bias correction does not retain the higher order properties of the observations, whereas the Delta method does not account for model changes in covariance; our proposed method does both. FROM POPPICK ET AL, 2015.

Again using the CCSM3 model, but now an ensemble of transient simulations rather than a single stationary run, we describe a statistical model for the evolution of temporal covariances in a GCM in response to altered CO2 levels. We find that, at least in CCSM3, changes in the local covariance structure can be explained largely as a function of the regional mean change in temperature, with a small term related to the rate of change of warming. (A warming climate will have slightly greater variability than an equilibrium climate at the same temperature.)

The statistical model can then be used to emulate the evolving covariance structure of temperatures, and therefore used to create in data-driven simulations that account for the projections of changes while still retaining fidelity with the observational record. We demonstrate the emulation of variability changes below, training the statistical model on model runs of several CO2 scenarios and using it to emulate changes under another scenario. Variability changes can indeed be described and emulated with a simple statistical model.

Figure 3: (Left): ESTIMATES OF CHANGES IN VARIABILITY IN Future (year 2100) RELATIVE TO PREINDUSTRIAL CLIMATEs, UNDER A MONOTONICALLY INCREASING CO2 CONCENTRATION SCENARIO IN CCSM3. RED INDICATES AN INCREASE IN VARIABILITY AND BLUE A DECREASE IN VAR…

Figure 3: (Left): ESTIMATES OF CHANGES IN VARIABILITY IN Future (year 2100) RELATIVE TO PREINDUSTRIAL CLIMATEs, UNDER A MONOTONICALLY INCREASING CO2 CONCENTRATION SCENARIO IN CCSM3. RED INDICATES AN INCREASE IN VARIABILITY AND BLUE A DECREASE IN VARIABILITY. THE ASSOCIATED TIMESCALE OF VARIATION FOR WHICH THE CHANGES IN VARIABILITY ARE DISPLAYED IS GIVEN ABOVE ON THE Y-AXIS OF THE  MAP.  (RighT):, emulation of The Estimates on the LEFT. Locations are marked with “.” (or “x”) when the difference between the emulator and the fitted model is more than two (or three) standard errors away from zero. The patterns are similar under both schemes, with most of the differences at locations where our model would not be expected to be a good description of changes in variability (e.g., at ice margins). From Poppick et al, 2015.

 

Climate model emulation

The study of future climate change necessarily involves numerical simulations. However, the state-of-the-art general circulation models (GCMs) used to producing climate projections are highly computationally demanding, so that century-long model runs may take months of calendar time to complete. GCMs are therefore too unwieldy for use in many contexts, especially in policy analysis. Analyses that involve optimal policy determination or uncertainty quantification require repeated iterations of climate projections for different CO2 trajectories, something that is computationally prohibitive with GCMs. We therefore seek techniques that capture (emulate) the behavior of GCMs but in computationally leaner tools that are useful in climate impacts assessments and other policy analyses requiring rapid climate projections. 

We have developed a new approach for emulating the output of a GCM under arbitrary forcing scenarios using only a small set of precomputed runs from the model. We express temperature and precipitation are expressed as simple functions of the past trajectory of atmospheric CO2 concentrations, and fit the statistical model using a "training set" of model output (Castruccio et al, 2014). The approach captures the nonlinear evolution of climate anomalies shown in coupled climate models,  and, once the statistical model is fit, produces emulated climate output effectively instantaneously. In our 2014 paper, we demonstrated emulation of temperature and precipitation over sub-continental regions. In current work, we have extended emulation to individual model grid cells, by incorporating information about spatial correlations (Figure 1 below).

Figure 1: Improving emulation by incorporating spatial information. Figure shows model output from CCSM3 (grey, mean of five realizations) and emulations with various models (colored) for scenarios with an abrupt drop (Left) and jump (Righ…

Figure 1: Improving emulation by incorporating spatial information. Figure shows model output from CCSM3 (grey, mean of five realizations) and emulations with various models (colored) for scenarios with an abrupt drop (Left) and jump (Right) in CO2 concentration.  Top And bottom rows show two different example grid cells  (A and B are in North America, C and D in central africa).  The spatio-temporal statistical model Is shown in red. the spatio-temporal emulator adequately captures temperature trends At pixel level even when CO2 IS CHANGING RAPIDLY.

 

Multi-model emulation:

In the work described above, we demonstrate emulation on a single climate model, CCSM3, with a custom-designed training set. In current work, we are extending the technique to emulate all major state-of-the-art GCMs, using only publicly available archived model output in the CMIP5 archive.  We have tested 22 individual CMIP5 models to date; in all the simple emulation approach captures regional temperature evolution well.

Figure 2: Emulation of global mean temperature evolution FOR TWELVE GCMS under the IPCC RCP4.5 scenario of future CO2. each Model is emulated separately., with the emulator trained on model temperature output generated with the RCP60 and R…

Figure 2: Emulation of global mean temperature evolution FOR TWELVE GCMS under the IPCC RCP4.5 scenario of future CO2. each Model is emulated separately., with the emulator trained on model temperature output generated with the RCP60 and RCP85 CO2 scenarios. Emulation (blue) adequately captures actual GCM output (black).

Our work in multi-model emulation also includes evaluation of how small training sets can be (in terms of multiple scenarios or realizations) to maintain adequate emulation, and physical interpretation of the fitted emulation parameters for each model. These fitted parameters offer a metric with interpretable physical significance that allows us to characterize and classify climate models.

People: 

Elisabeth Moyer | Aidan Sadowski | Michael Stein | Shanshan Sun | Andrew Poppick | Jingyu Bao

Alumni:

 Feifei Liu Crouch |  William Leeds | David McInerney

Recent Publications:

Precipitation in future climates: rainstorms

This project is aimed at characterizing changes in future precipitation suggested by climate models, and using that information along with observations to produce data-driven simulations that can be used to evaluate impacts (e.g. changes in floods, water supply, and agricultural yields).

Climate models robustly predict a rise that global rainfall will rise under warmer climate conditions, and that rain events become more intense (because a warmer atmosphere holds more water vapor). The combination also requires some third adjustment: because the projected increase in rainstorm intensity (about 6% per degree C of warming) is greater than the projected total rise in precipitation (about 3%/degree), some other aspect of precipitation in the climate models must change to compensate. Either storms are initiated less frequently, or rain events are shorter in duration, or rainstorms are smaller in extent, or some combination of these effects. To date no study has determined how exactly precipitation characteristics change, in part because of the statistical challenge of representing precipitation, a complex, non-Gaussian process with strong spatiotemporal correlations. 

We have developed a new approach to characterizing changes in precipitation in future climates that moves away from considering time series at individual locations, and instead uses image-processing algorithms to identifying individual rainstorm events and capture the evolution of storms in space-time space. The effort does require a very large volume of data: observations and model output at very high spatial resolution on time steps of a few hours or less. In this study, we use dynamically downscaled climate model output (CESM downscaled with WRF)  at 12 km, 3 hourly resolution. We compare model output to real-world precipitation estimated from weather radar: the NCEP Stage IV data product at 4 km, 1 hourly resolution.

 FIGURE 1: EXAMPLE OF RAINSTORM OBJECTS CONSTRUCTED BY OUR RAINSTORM IDENTIFICATION AND TRACKING ALGORITHM, FROM MODEL OUTPUT IN (LEFT) SUMMER AND (RIGHT) WINTER.  COLORS REPRESENT DIFFERENT IDENTIFIED RAINSTORM OBJECTS. MOVEMENT OF SAME-COLORED REGIONS SHOW EVOLUTION OF INDIVIDUAL RAINSTORMS OVER TIME. IN SUMMER, CONVECTIVE STORMS INITIATE OVER MUCH OF THE U.S., ESPECIALLY IN THE WEST. RAINFALL IN WINTER OCCURS PRIMARILY IN LARGE-SCALE STORMS.

In the model output studied, we find that storms do indeed become more intense in future climate conditions but that their initiation, duration, and trajectories are nearly unchanged. The compensating adjustment is that storms simply become smaller in spatial extent.

The goal of this work is not only improving scientific understanding but developing statistical methods for generating future precipitation scenarios that combine information from both observational data and climate model runs. That is, simulations should capture the changes shown in models but preserve the statistical characteristics of real-world rainfall.  Altering the spatial extent of rain events in data is non-trivial. Figure 2 below shows the application of a transformation algorithm applied to Stage IV radar data. In general, we hope this work will enable new ways of thinking about precipitation events in both meteorology and climatology. 

Figure 2: Example OF precipitation patterns in observations (left) and after transforming the rainstorm size and intensity (RIGHT). Red crosses represent the intensity-weighted center of individual rainstorms. To simulate the rainstorm siz…

Figure 2: Example OF precipitation patterns in observations (left) and after transforming the rainstorm size and intensity (RIGHT). Red crosses represent the intensity-weighted center of individual rainstorms. To simulate the rainstorm size changes suggested by our regional climate model runs, we shrink Or enlarge each rainstorm by reducing or increasing the distance to this center from individual rainstorm pixels. After resizing, we transform the rainfall intensity time series at each location to reflect the intensity distribution changes suggested by our model runs.

People

Won Chang | V. Rao KotamarthiElisabeth MoyerMichael Stein | Jiali Wang

Recent Presentations:

  • Won Chang. A conditional simulation approach to future precipitation scenario generation. Poster presented at: Data Assimilation Workshop, STATMOS Summer School. May 2015.

  • Won Chang. A conditional simulation approach to future precipitation scenario generation. Poster presented at: STATMOS Annual Meeting. University of Chicago, Chicago, IL. September 2014. 

Climate extremes

Changes in U.S. temperature extremes under increased CO2 in millennial-scale climate simulations 

Changes in extreme weather may produce some of the largest societal impacts from anthropogenic climate change. (At present, weather damages are dominated by rare events that happen only every several decades or more.) However, predicting future changes in those rare events is not possible using only the short observational record.  Insight on changes in extremes must come from climate models, where we can generate long simulations.

In this project, we use millennial runs from the Community Climate System Model version 3 (CCSM3) in equilibrated pre-industrial and future (700 and 1400 ppm CO2) conditions to examine both how extremes change and how well these changes can be estimated as a function of run length.

Figure 1: Illustration of how changes in extreme TEMPERAtures can differ from changes in OVERALl temperature distributions. Model output FrOM 1000-year RUNS OF CCSM3 with CO2 at pre-industrial (PI) and 2.5 x PI (700 ppm) levels, from a mid-latitudes…

Figure 1: Illustration of how changes in extreme TEMPERAtures can differ from changes in OVERALl temperature distributions. Model output FrOM 1000-year RUNS OF CCSM3 with CO2 at pre-industrial (PI) and 2.5 x PI (700 ppm) levels, from a mid-latitudes location (in Idaho).

(a, bottom): daily temperatures in Winter (DJF), with PI in blue and future in  RED. (B, TOP). annual cold extremes (Annual Winter TMIN), on same temperature scale. Note that Winter extreme cold TEmperatures WARM MORE strongly than The mean temperature shift. (From Huang et al, 2015)

Extreme value theory provides a means of estimating the far tails of distributions. We estimate changes to distributions of future temperature extremes by fitting annual maximum and minimum temperatures to generalized extreme value (GEV) distributions.  

Using 1000-year preindustrial and future timeseries of temperatures in the contiguous United States, we show that changes in extremes are different in summer and winter. In winter, cold extremes generally warm much more than the mean shift in wintertime temperatures (Figure 1), while in summer, warm extremes generally warm only as much as the shift in means. 

The changes in winter extremes involve more than a simple shift in magnitudes. The scale and shape of their distributions also change. This effect is best demonstrated by plotting the "return level" of extreme events, i.e. the magnitude of an extreme that recurs on a specified timescale. Figure 2 below shows changes in the 2-year to 100-year return levels for wintertime cold extremes. Over ocean regions (red lines at lower left are the Pacific off California; those on right are the Gulf of Mexico and the Atlantic), changes are relatively flat across time, suggesting a simple shift in the distribution of extremes. Over land regions, however, changes show complex behavior. In the inland Southwest U.S.,  changes in the 100-year "extreme extremes" are larger than changes in the 2-year "moderate extremes". Further north, this pattern is reversed.

FIGURE 2: ESTIMATED CHANGES IN RETURN LEVELS OF WINTER TMIN  EXTREMES IN 700 PPM CO2 VS. PRE-INDUSTRIAL MODEL RUNS FOR THE NORTH AMERICAN REGION ( 6 X 16 GRID CELLS). EACH OF THE 16 PANELS REPRESENTS A LONGITUDE; IN EACH PANEL, THE 6 …

FIGURE 2: ESTIMATED CHANGES IN RETURN LEVELS OF WINTER TMIN  EXTREMES IN 700 PPM CO2 VS. PRE-INDUSTRIAL MODEL RUNS FOR THE NORTH AMERICAN REGION ( 6 X 16 GRID CELLS). EACH OF THE 16 PANELS REPRESENTS A LONGITUDE; IN EACH PANEL, THE 6 LATITUDES ARE DENOTED BY COLOR. THE X AXIS OF EACH PANEL IS THE RETURN PERIOD AND Y AXIS THE CHANGE IN  RETURN LEVEL.THICK DASHED LINES SHOW ESTIMATED RETURN LEVEL CHANGES; THE ENVELOPES SHOW ASSOCIATED UNCERTAINTIES (ESTIMATED FROM BLOCK BOOTSTRAPPED SAMPLES WITH RESAMPLED YEARS).  FOR EACH GRID CELL, THE CORRESPONDING MEAN TEMPERATURE CHANGE IS MARKED WITH A SYMBOL (CROSSES FOR LAND AND BOXES FOR OCEAN LOCATIONS). (FROM HUANG ET AL., 2015.)

The 1000-year runs used here allow us to accurately determine changes in even 100-year extremes, but in practice, most modeling studies must rely on shorter runs. GEV modeling should allow estimation of infrequent events using relatively short time series, but it is important to understand its limitations for climate data. We therefore repeat the estimation of selected return levels (20-, 50-, and 100-year extremes) using segment of the timeseries of varying length. The resulting estimates can become very poor when the timeseries is of comparable length or shorter than the return period of interest: that is, a 100-year model run cannot be used to reliably estimate changes in 100-year extremes. These results suggest caution when attempting to use short observational records or model runs to infer changes in extreme events. 

People:

Whitney Huang | Elisabeth Moyer | Michael Stein | Shanshan Sun | David McInerney | Hao Zhang

Soil Moisture

Although many impacts of climate change remain uncertain, climate models robustly project that land surfaces become drier in warmer climate conditions. General circulation models (GCMs) almost universally show reduced soil moisture over much of the globe, even in locations where average rainfall increases (Figure 1). This change is a significant societal concern both because of potential effects on food production and because drier soils feed back on local weather conditions, exacerbating heat waves. (In wetter conditions, evaporative cooling buffers extreme temperatures.) These adverse consequences make understanding the surface-atmosphere interactions that control soil moisture a research priority.

Figure 1: PROJECTED CHANGES IN SOIL MOISTURE IN THE TOP 10 CM LAYER, FROM 1980-1999 TO 2080-2099 UNDER A MODERATE EMISSIONS SCENARIO: MULTIMODEL MEAN % ChANGE FROM 11 CMIP5 MODELS. MOST LAND AREAS BECOME DRIER. From (DAI et al., 2013).

Figure 1: PROJECTED CHANGES IN SOIL MOISTURE IN THE TOP 10 CM LAYER, FROM 1980-1999 TO 2080-2099 UNDER A MODERATE EMISSIONS SCENARIO: MULTIMODEL MEAN % ChANGE FROM 11 CMIP5 MODELS. MOST LAND AREAS BECOME DRIER. From (DAI et al., 2013).

While the model projections of drying are plausible, the model representation of soil moisture dynamics have not been fully validated. Soil moisture observations are sparse and controversial, and there is as yet no consensus on whether the historical record suggests that droughts are increasing. The RDCEP soil moisture project is using a unique observational resource to characterize soil dynamics and test the representation of soil moisture in models. We make use of the Department of Energy’s Atmospheric Radiation Measurement Climate Research Facility at the Southern Great Plains site in Kansas and Oklahoma (ARM-SGP), where soil moisture, moisture fluxes, and meteorological variables have been measured at hourly intervals over more than a decade, as a test-bed for characterizing statistical relationships between soil moisture and moisture fluxes and local forcing variables (e.g. temperature, precipitation, wind speed, and relative humidity). We then compare to output from the soil moisture model (CLM) used in the Community Earth System Model, one of the most widely-used climate models, to understand how well models capture the physics controlling soil moisture.

People: 

Scott Collis |  Elisabeth Moyer | Michael Stein | Shanshan Sun 

Support: 

This project is supported in part by a University of Chicago - Argonne National Laboratory Strategic Collaboration Initiative seed grant.

Recent Presentations:

  • Shanshan Sun. Statistical exploration of processes controlling soil moisture in present and future climates. Poster presented at: ARM/ASR Joint User Facility/PI Meeting, Tyson’s Corner, VA. March 16-19, 2015. 

Robustness in economic models with climate change

We are working on robust control in simple climate models coupled with economic growth, finance, and macro models. This effort is important for carbon pricing, understanding the impacts of robustly accounting for climate change, technical change, and other sources of uncertainty, not only on economic growth and macroeconomics, but also on asset prices as well as insurance pricing and green-house gas emission pricing.

Within the context of climate models are three contrasting approaches to robustness:

  • adapting to potential model misspecification
  • robustness adjustments for prior/posterior uncertainty
  • “smooth” models of ambiguity aversion

Decision theoretic frameworks exist for all of these applications, but their full consequences for economic models with climate change remains to be explored. To accomplish this we are developing numerical methods to support these analyses. In terms of discounting, our focal point is on the consequences of uncertainty. There is an extensive literature from asset pricing and macroeconomics that uses stochastic discounting as a device to adjust discount rates for cash flow riskiness. We are drawing on this literature and incorporating compensations for aversion to ambiguity and concerns about model misspecification into models that feature explicit uncertainty and climate impacts on the economy. We are also contrasting pricing implications from market economies with social valuation.

Coupled Economic-Climate Models with Carbon-Climate Response: The economics of global climate change is characterized by fundamental uncertainties including the appropriate reduced forms for climate dynamics, the specification of economic damages resulting from climate change, and mechanisms by which these damages will affect long-run economic growth. We have developed and implemented a novel theoretical and computational integrated assessment modeling approach that is a well-grounded means of summarizing the fundamental relationship between human activity and the global climate for purposes of economic analysis. Using a dynamic integrated assessment framework, this project makes several contributions to improving the analysis of these uncertainties:

  • First, we incorporate the cumulative climate response (CCR) function developed by Matthews et al. for representing the basic relationship between anthropogenic carbon emissions and increases in global mean temperature in a manner that is more directly policy relevant than the usual approach based on the equilibrium climate sensitivity.
  • Second, we adapt the tools developed by Hansen, Sargent and others for robustness analysis to address underlying model uncertainty in both economic and climate dynamics.
  • Third, we allow climate change to affect economic growth directly, in addition to its effect on output

We then develop and study a simple analytical model that yields insights and results on the key implications of these assumptions, as well as facilitating the interpretation of numerical results from a more general model. Among our findings is that the presence of robustness may result in either a decrease or increase in the optimal carbon tax and energy usage, depending among other factors on societal preferences.

Numerical Methods:  In order to perform robustness analysis over a variety of fundamentally different models, we are developing PDE-based numerical method for solving robust stochastic equilibrium models in continuous time with optimal control. This approach allows us to uncover solutions for the entire state space and to perform a quick sweep over variety of models and parameters.

People

Evan Anderson | William Brock | Lars Hansen  | Alan Sanstad  | Victor Zhorin

Alumni

Botao Wu

Recent Publications

The next great Dust Bowl

The 1930s dust bowl included the worst single-year drought in the western U.S. in the last millennia, as well as some of the worst droughts in the U.S. as a whole in the last century. We are investigating the agricultural impacts of a similar drought in the present-day and in the near future, when temperatures are warmer and CO2 concentrations are higher. Initial results suggest that a drought of this magnitude in current climate would result in production losses significantly larger than any observed drought in recent history.

Initial results indicate that by the middle of the 21st century, warmer temperatures and a dust-bowl-like drought may lead to a decrease in maize production of up to 80%.  However, the 1930s Dust Bowl was unique in not only its intensity, but also its persistence. Droughts were widespread in three out of five years, and temperatures elevated throughout the decade. The implications of successive droughts on crop prices remain unclear. Current research is focusing on linking persistent large-scale drought to crop prices through stock-to-use ratios. 

 

People: 

Joshua Elliott | Michael Glotter

Model Uncertainty and Energy Technology Policy

Energy modeling, numerical modeling based on economic principles, has become the dominant analytical tool in U.S. energy policy. Energy models are widely used by researchers across the public and private sectors. However, the widespread application of these models in policy analysis poses challenges to decision-makers. We are developing a framework and analysis that demonstrate how non-Bayesian decision rules can address fundamental model uncertainty in the domain of energy policy, technological change, and greenhouse gas abatement.

Numerical modeling based on economic principles has become the dominant analytical tool in U.S. energy policy. Energy models are now used extensively by public agencies, private entities, and academic researchers, and in recent years have also formed the core of integrated assessment models used to analyze the relationships among the energy system, the economy, and the global climate. However, fundamental uncertainties are intrinsic in what has become the typical circumstance of multiple models embodying different representations of the energy-economy, and producing different policy-relevant outputs that model users are compelled to interpret as equally plausible and/or valid. Because the policy implications of these outputs can diverge substantially, policy-makers are confronted with a significant degree of model-based uncertainty and little or no guidance as to how it should be addressed.

This problem of “model uncertainty” has recently been the focus of work in macroeconomics, where scholars have studied the problem of how a decision-maker should proceed in the face of uncertainty re- garding the correct model of an economic system that is the object of policy. We focus on analyzing a low-dimensional numerical integrated assessment model using the “minimax regret” metric. Specifically, we have demonstrated that deep uncertainty regarding energy-related technological change can be addressed using this approach. Our findings include comparison with expected cost minimization, to show how the interaction of solution methods and model structure affect the influence of this form of deep uncertainty on low-run CO2 emissions abatement policies. We also examined other methods assuming some prior distribu- tions over uncertain parameters for analyzing the difference between our robust solution and the non-robust solution from those methods.

We demonstrate that the fundamental model uncertainty can be represented and analyzed in the context of energy policy problems determining optimal CO2 abatement strategies. The robust solution from min-max regret method is significantly different with any solutions from sensitivity analysis over uncertain parameters or those methods assuming prior distributions over uncertain parameters. The following figure shows the difference of the robust min-max regret solution over all three uncertain parameters (the red line) and others with min-max regret solution over only one uncertain parameter, Technical Change level, while the other two parameters are used for sensitivity analysis. 

Abatement paths in computational model of min-max regret criterion with three uncertain parameters, and sensitivity analysis of min-max regret criterion with only one uncertain parameter

Abatement paths in computational model of min-max regret criterion with three uncertain parameters, and sensitivity analysis of min-max regret criterion with only one uncertain parameter

 

People

Yongyang Cai | Alan Sanstad

Alumni: Kenneth Judd 

Recent Publications

Agricultural Adaptation and the Green Revolution

This graph shows corn yield in bushels per acre at fields near Urbana plotted against year. Each box plot shows the spread of the yield in a given year, with outliers labeled. Overplotted is a cubic line that is representative of the overall Urbana …

This graph shows corn yield in bushels per acre at fields near Urbana plotted against year. Each box plot shows the spread of the yield in a given year, with outliers labeled. Overplotted is a cubic line that is representative of the overall Urbana yield vs. year trend.

One of the most difficult aspects of evaluating the impacts of climate change on future agricultural production is compensating for change in technology over many decades.

An acceleration of technological research with regards to agriculture, commonly referred to as “The Green Revolution”, occurred worldwide from roughly 1940 to the late 1960’s. Practices in irrigation, management, and hybridization were changed drastically.

A dataset was compiled from the Illinois Digital Environment for Access to Learning and Scholarship (IDEALS). It is a unique record of the evolution of technology and management practices during the Green Revolution, making it one of the very few lengthy, continuous, and consistent records of technological change in agriculture that exists in the entire world.

This historical record will be used to better understand the range of possible technology pathways for the next several decades and what they mean for food security in the face of a changing climate.

People

Joshua Elliott | Sydney Purdue | Alison Brizius | Anna Blinderman

Social cost of carbon and climate impacts on economic growth

The social cost of carbon, defined as the present value change in consumption due to an incremental change in carbon emissions, is used by federal agencies in cost-benefit analysis of any regulation that changes emissions. In 2009, the Oce of Management and Budget, through an interagency process, estimated the social cost of carbon and required all agencies to use their estimate. Their central value was $21.40/tCO2 with range of -$2.7/tCO2 to $142.4/tCO2.

The government’s estimate of the social cost of carbon, which is consistent with estimates by private researchers, implicitly assumes that the economic growth continues even with substantial temperature increases. In one of the models, temperatures increase by 6.3 by the year 2300. With this temperature increase, the global economy is roughly 30 times larger than it is today on a per capita basis for the model. The apparent reason for this estimate is that damages from climate change do not affect growth. They are modeled as reducing usable output in a given year with exogenously specified growth continuing regardless of climate damages. 

We estimate the social cost of carbon when climate change reduces the growth rate of the economy. We use the same model as the OMB but modify it so that a fraction of damages from climate change affect the growth rate rather than simply reducing usable output. Growth might be reduced, for example, because resources are diverted from research to adaptation to climate change. Even relatively small growth effects produce substantial change in the social cost of carbon suggesting that (1) the estimates of the social cost of carbon are not robust to modest changes in the estimating model and (2) research into the impacts of climate change should focus on growth effects rather than level effects. 

Estimates of the social cost of carbon when accounting for climate damages.

Estimates of the social cost of carbon when accounting for climate damages.

Chicago Climate Online

Chicago Climate Online (C2O), a policy initiative of the University of Chicago Law School, is a climate change policy research tool that seeks to become a preeminent policy resource for researchers. It combines features of Wikipedia -- with summaries of current knowledge on climate policy topics -- and a database -- with a list of research articles on those topics. It is an open-access research tool that allows users to also upload, create, and organize their own content.

The focus of C2O is on policy issues, not climate science. Although climate science continues to develop rapidly, we take the core of existing climate science as given. C2O is not a place to debate the validity of this science. The interconnectedness of policy and science, however, makes clear demarcation between the two impossible. C2O was developed with a primary focus on climate change policy, with inclusion of more science-focused topics and references only as necessary.

Visit Website »

People

David Weisbach | Alison Brizius

Tipping Points in Dynamic Stochastic Integrated Assessment Models

There is great uncertainty about the impact of anthropogenic carbon on future economic wellbeing.

We use DSICE, a dynamic stochastic general equilibrium (DSGE) model of integrated climate and economy to account for abrupt and irreversible climate change. DSICE is an extension of the DICE2007 model of William Nordhaus, which incorporates beliefs about the uncertain economic impact of possible climate tipping events and uses empirically plausible parameterizations of Epstein-Zin preferences to represent attitudes towards risk. 

In this series of studies we model climate shocks in the form of a stochastic tipping points,  and investigate the impact of a tipping point externality on optimal mitigation policy. We find that the uncertainty associated with anthropogenic climate change imply carbon taxes much higher than implied by deterministic models. This analysis indicates that there is much greater urgency to immediately enact significant GHG policies than implied by DICE2007 and similar models that ignore uncertainty.

THE BASELINE CARBON TAX IN THE DETERMINISTIC DICE MODEL (WITH NO TIPPING POINT) IS SHOWN IN BLACK. THE EXPECTED ADDITIONAL CARBON TAX WHEN INCLUDING A STOCHASTIC TIPPING POINT IS INDICATED IN BLUE. RED LINES INDICATE THE ADDITIONAL CARBON TAX WHEN T…

THE BASELINE CARBON TAX IN THE DETERMINISTIC DICE MODEL (WITH NO TIPPING POINT) IS SHOWN IN BLACK. THE EXPECTED ADDITIONAL CARBON TAX WHEN INCLUDING A STOCHASTIC TIPPING POINT IS INDICATED IN BLUE. RED LINES INDICATE THE ADDITIONAL CARBON TAX WHEN THE EXPONENT OF THE DAMAGE FUNCTION IN THE DETERMINISTIC DICE MODEL IS INCREASED TO FOURTH (SOLID LINE) AND SIXTH (DASHED LINE) ORDER. REPRODUCED FROM (LONTZEK ET AL 2015)

 

People

Yongyang Cai | Kenneth Judd | Timothy Lenton | Thomas Lontzek 

Recent Publications

RPS Calculator

The major policy instruments for mitigating climate change actually in use in the U.S. are subsidies provided to renewable energy. A popular means of subsidy is through state-level Renewable Portfolio Standards (RPSs), requirements enacted by many states that require a certain fraction of electricity must be derived from renewables. However, many state RPSs are infeasible because of “cost cap” provisions that do not permit renewables to be sufficiently competitive, and feasibility is generally not assessed before legislation is passed. The RPS calculator allows the user to explore the conditions for RPS success or failure in different states. The user can analyze and modify existing state statutes or design new statutes in states that do not have them. Users can explore the effects of parameters such as electricity prices, generation costs for wind and solar, interest rates, technology carveouts, and cost cap structure. Features under development include extension to all U.S. states (currently only IL and CA) and spatial variation in wind speed (wind capacity factor map).

People

Current: Elisabeth Moyer

Alumni: Sean Johnson | Lexi Goldberger | Joe Zhu

Recent Publications

Resilience of the Global Food System

Photo Credit: R.Gustafson / USAID

Photo Credit: R.Gustafson / USAID

The impact of climate change on future crop production will depend on the impact on mean yields and on the incidence and scale of extreme outcomes.  More frequent droughts, heatwaves, and floods due to climate change could severely disrupt the world's food supply.  We are participating in an international scientific task force on extreme weather and global food system resilience.  The goal is to assess the impact of extreme weather upon the global food chain and offers recommendations on how to mitigate the potential damage from future crises. 

The task for has released its findings in a new report.  The report assessed the impact of extreme weather upon the global food chain and offers recommendations on how to mitigate the potential damage from future crises.  

Our  research on "once-in-a-century" food threats -- which we found may become as frequent as once every decade by the 22nd century -- is a prominent feature of the report's warnings, which additionally highlight the concentration of the world's most important crops (maize, soybeans, wheat, and rice) in very few countries, which could exacerbate the damage if extreme weather hits one or more of these regions. An increasingly global market for agricultural products offers protection against small, local food shortages, but may actually increase vulnerability for global shocks, the report also concludes.

“It is likely that the effects of climate change will be felt most strongly through the increasing frequency of extreme weather events such as droughts, heatwaves and floods and their impact on the production and distribution of food – something we almost take for granted," said Professor Tim Benton of the Global Food Security program, in a news release on the report. “This study presents a plausible scenario for how the food system might be impacted by extreme weather, alongside a series of recommendations that should help policy and business plan for the future. Action is urgently needed to understand risks better, improve the resilience of the global food system to weather-related shocks and to mitigate their impact on people.”

Recommendations in the primary report include more research to understand the risks and "worst case scenarios" of future extreme weather, improved coordination between governments on how to handle such food supply shocks, improving the resilience of international markets, and adapting agriculture to expected climate change to protect crops. 

Distributions of global calorie-weighted yield of maize, soy, wheat, and rice for the historical (1951-2010) and near-future (2011-2040) period with (top) and without (bottom) the effects of fertilization from increasing atmospheric CO2 included. Th…

Distributions of global calorie-weighted yield of maize, soy, wheat, and rice for the historical (1951-2010) and near-future (2011-2040) period with (top) and without (bottom) the effects of fertilization from increasing atmospheric CO2 included. The estimated magnitude of a 1-in-200 year event in each period is indicated.

In addition to the primary findings, we also co-authored "Annex A" of the report, which focused on Climate and Global Food Production Shocks. We studied the relationship between weather and historical food production shocks (such as the "Dust Bowl" drought of the 1930s), then used crop-climate models to simulate present-day or near-future extreme events.

The report highlights the extreme concentration of certain globally-important crops. For example, almost 60% of the world's maize is grown in the United States and China. For soybeans, 80% of the global supply is farmed in just three countries: the U.S., Brazil, and Argentina. A severe drought in any of these countries, even if it only decreases crop production by 10%, would likely have dramatic consequences for the rest of the world. However high-resolution climate model runs are needed to better quantify the risk of extreme events and food shocks.

Phase II of this project launches in Fall 2015 and will focus on multiple bread-basket failure.

For more on the report, see news stories in The New York TimesScienceThe Guardian, and BBC.

People

Joshua Elliott | Michael Glotter | David KellyJames Chryssanthacopoulos

Publications

Global Change, Food and Trade

The consequences of climate change are most often depicted as acts of nature: warmer temperatures, rising sea levels, and more frequent extreme weather. But a shifting climate will also dramatically affect civilization’s oldest industry, agriculture. More and more scientists are focusing their attention on how climate change will affect crop yield, farming practices, and food security for the world, and the massive implications for economics and health. Research on these important topics by Joshua Elliott of RDCEP is heavily featured in “Climate Change and Food Systems,” a new, free-to-download book published by the Food and Agriculture Organization (FAO) of the United Nations.

The primary mission of FAO is to improve food security around the world, by reducing hunger and poverty while protecting natural resources. But all of these goals are moving targets in light of climate change and its uncertain downstream effects upon global supplies of food and water. In her foreword to the new book, FAO Deputy Director-General Natural Resources Maria Helena Semedo wrote that the response to these changes must come from several quarters.

“The growing threat of climate change to the global food supply, and the challenges it poses for food security and nutrition, requires urgent concerted policy responses and the deployment of all the scientific knowledge and accumulated evidence at our disposal,” Semedo said. “It also requires a sharper focus on the important drivers of climate adaptation, including the potential role of trade to mitigate some of the negative impacts of climate change on global food production.”

Elliott’s two chapters in the book demonstrate the current state and potential of computational modeling to anticipate and prepare for these challenges. By connecting advanced climate models and their outputs to agricultural and economic models, researchers can make robust forecasts of how climate change will affect food supplies -- and vice versa. For example, shifting dietary preferences in the developing world could lead to changes in planting practices or more livestock farming, further elevating atmospheric carbon and accelerating climate change.

Modeling can also help policymakers and scientists understand the regional relationship between climate and food systems. While areas near the equator will likely become less amenable to agriculture in a warmer, drier future, areas far North or South may grow more hospitable to crops. Some regions may also experience a form of “hidden hunger” -- decreased nutritional value from staple crops such as wheat and rice that contain less vitamins and protein when grown in a high-CO2 atmosphere.

These forecasts can also inform potential policies that help society adapt or mitigate the damage of climate change. New trade policies, for instance, may help alleviate food security concerns from the drift of food production to regions farther away from the equator. Modeling studies may also provide informed guidance for issues such as water supply management and  the use of agricultural biotechnology, biofuels, and fertilizer.

“This book brings together some of the world’s leading experts in climate change and food systems to produce a comprehensive report on the state of modeling and science around these issues,” Elliott said. “For almost 70 years, the FAO has been a major driving force for improving global agriculture, food security, and human health. With the credibility and scope that FAO brings to this debate, I hope this exciting and comprehensive book can help drive the conversation about global change in the coming decades, and inspire an increased focus on food security and climate change adaptation, as well as a renewed call to researchers and modelers to address these essential topics. “

For more on “Climate Change and Food Systems,” see the UN FAO release.

People

Joshua Elliott

Publications