Won Chang

Fast dimension-reduced climate model calibration and the effect of data aggregation

How will the climate system respond to anthropogenic forcings? One approach to this question relies on climate model projections. Current climate projections are considerably uncertain. Characterizing and, if possible, reducing this uncertainty is an area of ongoing research. We consider the problem of making projections of the North Atlantic meridional overturning circulation (AMOC). Uncertainties about climate model parameters play a key role in uncertainties in AMOC projections. When the observational data and the climate model output are high-dimensional spatial data sets, the data are typically aggregated due to computational constraints. The effects of aggregation are unclear because statistically rigorous approaches for model parameter inference have been infeasible for high-resolution data. Here we develop a flexible and computationally efficient approach using principal components and basis expansions to study the effect of spatial data aggregation on parametric and projection uncertainties. Our Bayesian reduced-dimensional calibration approach allows us to study the effect of complicated error structures and data-model discrepancies on our ability to learn about climate model parameters from high-dimensional data. Considering high-dimensional spatial observations reduces the effect of deep uncertainty associated with prior specifications for the data-model discrepancy. Also, using the unaggregated data results in sharper projections based on our climate model. Our computationally efficient approach may be widely applicable to a variety of high-dimensional computer model calibration problems.

Changes in Spatiotemporal Precipitation Patterns in Changing Climate Conditions

AbstractClimate models robustly imply that some significant change in precipitation patterns will occur. Models consistently project that the intensity of individual precipitation events increases by approximately 6%–7% K−1, following the increase in atmospheric water content, but that total precipitation increases by a lesser amount (1%–2% K−1 in the global average in transient runs). Some other aspect of precipitation events must then change to compensate for this difference. The authors develop a new methodology for identifying individual rainstorms and studying their physical characteristics—including starting location, intensity, spatial extent, duration, and trajectory—that allows identifying that compensating mechanism. This technique is applied to precipitation over the contiguous United States from both radar-based data products and high-resolution model runs simulating 80 years of business-as-usual warming. In the model study the dominant compensating mechanism is a reduction of storm size. In s...

Calibrating an Ice Sheet Model Using High-Dimensional Binary Spatial Data

Rapid retreat of ice in the Amundsen Sea sector of West Antarctica may cause drastic sea level rise, posing significant risks to populations in low-lying coastal regions. Calibration of computer models representing the behavior of the West Antarctic Ice Sheet is key for informative projections of future sea level rise. However, both the relevant observations and the model output are high-dimensional binary spatial data; existing computer model calibration methods are unable to handle such data. Here we present a novel calibration method for computer models whose output is in the form of binary spatial data. To mitigate the computational and inferential challenges posed by our approach, we apply a generalized principal component based dimension reduction method. To demonstrate the utility of our method, we calibrate the PSU3D-ICE model by comparing the output from a 499-member perturbed-parameter ensemble with observations from the Amundsen Sea sector of the ice sheet. Our methods help rigorously characterize the parameter uncertainty even in the presence of systematic data-model discrepancies and dependence in the errors. Our method also helps inform environmental risk analyses by contributing to improved projections of sea level rise from the ice sheets. Supplementary materials for this article are available online.

Improving ice sheet model calibration using paleoclimate and modern data

Human-induced climate change may cause significant ice volume loss from the West Antarctic Ice Sheet (WAIS). Projections of ice volume change from ice-sheet models and corresponding future sea-level rise have large uncertainties due to poorly constrained input parameters. In most future applications to date, model calibration has utilized only modern or recent (decadal) observations, leaving input parameters that control the long-term behavior of WAIS largely unconstrained. Many paleo-observations are in the form of localized time series, while modern observations are non-Gaussian spatial data; combining information across these types poses non-trivial statistical challenges. Here we introduce a computationally efficient calibration approach that utilizes both modern and paleo-observations to generate better-constrained ice volume projections. Using fast emulators built upon principal component analysis and a reduced dimension calibration model, we can efficiently handle high-dimensional and non-Gaussian data. We apply our calibration approach to the PSU3D-ICE model which can realistically simulate long-term behavior of WAIS. Our results show that using paleo observations in calibration significantly reduces parametric uncertainty, resulting in sharper projections about the future state of WAIS. One benefit of using paleo observations is found to be that unrealistic simulations with overshoots in past ice retreat and projected future regrowth are eliminated.

Diagnosing added value of convection-permitting regional models using precipitation event identification and tracking

Dynamical downscaling with high-resolution regional climate models may offer the possibility of realistically reproducing precipitation and weather events in climate simulations. As resolutions fall to order kilometers, the use of explicit rather than parametrized convection may offer even greater fidelity. However, these increased resolutions both allow and require increasingly complex diagnostics for evaluating model fidelity. In this study we focus on precipitation evaluation and analyze five 2-month-long dynamically downscaled model runs over the continental United States that employ different convective and microphysics parameterizations, including one high-resolution convection-permitting simulation. All model runs use the Weather Research and Forecasting Model driven by National Center for Environmental Prediction reanalysis data. We show that employing a novel rainstorm identification and tracking algorithm that allocates essentially all rainfall to individual precipitation events (Chang et al. in J Clim 29(23):8355–8376, 2016 ) allows new insights into model biases. Results include that, at least in these runs, model wet bias is driven by excessive areal extent of individual precipitating events, and that the effect is time-dependent, producing excessive diurnal cycle amplitude. This amplified cycle is driven not by new production of events but by excessive daytime enlargement of long-lived precipitation events. We further show that in the domain average, precipitation biases appear best represented as additive offsets. Of all model configurations evaluated, convection-permitting simulations most consistently reduced biases in precipitation event characteristics.

Statistics and the Future of the Antarctic Ice Sheet

37 One of the enduring symbols of the impact of climate change is that of a polar bear drifting in the sea, alone on its own piece of ice. For those who are left untouched by the loneliness of drifting polar bears, images of partially submerged lands and the devastation wrought by storm surges showcase some potentially frightening impacts of sea level rise on human life. The threat of sea level rise, in turn, is linked to the melting of ice sheets. Ice sheets are, therefore, important to understanding our planet, as well as learning about how our future may be affected by climate change. A promising approach to improving our understanding of ice sheets and derive sound projections of their future is to combine ice sheet physics, statistical modeling, and computing. First, what exactly is an ice sheet? It is an enormous mass of glacial land ice, more than 50,000 square kilometers in extent. The Antarctic ice sheet extends over 14 million square kilometers while the Greenland ice sheet extends over 1.7 million square kilometers. To put this in perspective, the area covered by the Antarctic ice sheet is comparable to the continental United States and Mexico combined. In fact, the Greenland and Antarctic ice sheets contain more than 99% of the freshwater ice in the world. Roughly speaking, melting the entire Greenland ice sheet would result in sea level rise of around 7 meters (23 feet) while if the entire Antarctic ice sheet melted, it would result in sea level rise of around 57 meters (187 feet). Statistics and the Future of the Antarctic Ice Sheet