An official website of the United States government
Abstract Data required to calibrate uncertain general circulation model (GCM) parameterizations are often only available in limited regions or time periods, for example, observational data from field campaigns, or data generated in local high‐resolution simulations. This raises the question of where and when to acquire additional data to be maximally informative about parameterizations in a GCM. Here we construct a new ensemble‐based parallel algorithm to automatically target data acquisition to regions and times that maximize the uncertainty reduction, or information gain, about GCM parameters. The algorithm uses a Bayesian framework that exploits a quantified distribution of GCM parameters as a measure of uncertainty. This distribution is informed by time‐averaged climate statistics restricted to local regions and times. The algorithm is embedded in the recently developed calibrate‐emulate‐sample framework, which performs efficient model calibration and uncertainty quantification with onlymodel evaluations, compared withevaluations typically needed for traditional approaches to Bayesian calibration. We demonstrate the algorithm with an idealized GCM, with which we generate surrogates of local data. In this perfect‐model setting, we calibrate parameters and quantify uncertainties in a quasi‐equilibrium convection scheme in the GCM. We consider targeted data that are (a) localized in space for statistically stationary simulations, and (b) localized in space and time for seasonally varying simulations. In these proof‐of‐concept applications, the calculated information gain reflects the reduction in parametric uncertainty obtained from Bayesian inference when harnessing a targeted sample of data. The largest information gain typically, but not always, results from regions near the intertropical convergence zone.
more »
« less
Calibration and Uncertainty Quantification of Convective Parameters in an Idealized GCM
Abstract Parameters in climate models are usually calibrated manually, exploiting only small subsets of the available data. This precludes both optimal calibration and quantification of uncertainties. Traditional Bayesian calibration methods that allow uncertainty quantification are too expensive for climate models; they are also not robust in the presence of internal climate variability. For example, Markov chain Monte Carlo (MCMC) methods typically requiremodel runs and are sensitive to internal variability noise, rendering them infeasible for climate models. Here we demonstrate an approach to model calibration and uncertainty quantification that requires onlymodel runs and can accommodate internal climate variability. The approach consists of three stages: (a) a calibration stage uses variants of ensemble Kalman inversion to calibrate a model by minimizing mismatches between model and data statistics; (b) an emulation stage emulates the parameter‐to‐data map with Gaussian processes (GP), using the model runs in the calibration stage for training; (c) a sampling stage approximates the Bayesian posterior distributions by sampling the GP emulator with MCMC. We demonstrate the feasibility and computational efficiency of this calibrate‐emulate‐sample (CES) approach in a perfect‐model setting. Using an idealized general circulation model, we estimate parameters in a simple convection scheme from synthetic data generated with the model. The CES approach generates probability distributions of the parameters that are good approximations of the Bayesian posteriors, at a fraction of the computational cost usually required to obtain them. Sampling from this approximate posterior allows the generation of climate predictions with quantified parametric uncertainties.
more »
« less
- Award ID(s):
- 1835860
- PAR ID:
- 10449484
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Volume:
- 13
- Issue:
- 9
- ISSN:
- 1942-2466
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The mechanisms underlying observed global patterns of partitioning precipitation () to evapotranspiration () and runoff () are controversially debated. We test the hypothesis that asynchrony between climatic water supply and demand is sufficient to explain spatio‐temporal variability of water availability. We developed a simple analytical model forthat is determined by four dimensionless characteristics of intra‐annual water supply and demand asynchrony. The analytical model, populated with gridded climate data, accurately predicted global runoff patterns within 2%–4% of independent estimates from global climate models, with spatial patterns closely correlated to observations (). The supply‐demand asynchrony hypothesis provides a physically based explanation for variability of water availability using easily measurable characteristics of climate. The model revealed widespread responsiveness of water budgets to changes in climate asynchrony in almost every global region. Furthermore, the analytical model using global averages independently reproduced the Budyko curve () providing theoretical foundation for this widely used empirical relationship.more » « less
-
Abstract Global ocean mean salinityis a key indicator of the Earth's hydrological cycle and the exchanges of freshwater between land and ocean, but its determination remains a challenge. Aside from traditional methods based on gridded salinity fields derived from in situ measurements, we explore estimates ofbased on liquid freshwater changes derived from space gravimetry data corrected for sea ice effects. For the 2005–2019 period analyzed, the differentseries show little consistency in seasonal, interannual, and long‐term variability. In situ estimates show sensitivity to choice of product and unrealistic variations. A suspiciously large rise insince ∼2015 is enough to measurably affect halosteric sea level estimates and can explain recent discrepancies in the global mean sea level budget. Gravimetry‐basedestimates are more realistic, inherently consistent with estimated freshwater contributions to global mean sea level, and provide a way to calibrate the in situ estimates.more » « less
-
Abstract In this paper, we are interested in the following question: given an arbitrary Steiner triple systemonvertices and any 3‐uniform hypertreeonvertices, is it necessary thatcontainsas a subgraph provided? We show the answer is positive for a class of hypertrees and conjecture that the answer is always positive.more » « less
-
Abstract Silicon stable isotope ratios (30Si) of over 150 stream water samples were measured during seven storm events in six small critical zone observatory (CZO) catchments spanning a wide range in climate (sub‐humid to wet, tropical) and lithology (granite, volcanic, and mixed sedimentary). Here we report a cross‐site analysis of this dataset to gain insight into stream30Si variability across low‐order catchments and to identify potential climate (i.e., runoff), hydrologic, lithologic, and biogeochemical controls on observed stream Si chemical and isotopic signatures. Event‐based30Si exhibit variability both within and across sites (−0.22‰ to +2.27‰) on the scale of what is observed globally in both small catchments and large rivers. Notably, each site shows distinct30Si signatures that are preserved even after normalization for bedrock composition. Successful characterization of observed cross‐site behavior requires the merging of two distinct frameworks in a novel combined model describing both non‐uniform fluid transit time distributions and multiple fractionating pathways in application to low‐order catchments. The combined model reveals that site‐specific architecture (i.e., biogeochemical reaction pathways and hydrologic routing) regulates stream silicon export signatures even when subject to extreme precipitation events.more » « less
