Climate field reconstruction uncertainty arising from multivariate and nonlinear properties of predictors

Climate field reconstructions (CFRs) of the global annual surface air temperature (SAT) field and associated global area-weighted mean annual temperature (GMAT) are derived in a collection of pseudoproxy experiments for the past millennium. Pseudoproxies are modeled from temperature (T), precipitation (P), T + P, and VS-Lite (VSL), a nonlinear and multivariate proxy system model for tree ring widths. Spatial patterns of reconstruction skill and spectral bias for the T + P and VSL-derived CFRs are similar to those previously shown using temperature-only pseudoproxies but demonstrate overall degraded skill and spectral bias for SAT reconstruction. Analysis of GMAT spectra nevertheless suggests that the true GMAT frequency spectrum is resolved by those pseudoproxies (T, T + P, and VSL) that contain some temperature information. The results suggest that mixed temperature and moisture-responding paleoclimate data may produce actual GMAT reconstructions with skill, error, and spectral characteristics like those expected from univariate and linear temperature responders, but spatially resolved CFR results should be analyzed cautiousl

Possible causes of data model discrepancy in the temperature history of the last Millennium

Model simulations and proxy-based reconstructions are the main tools for quantifying pre-instrumental climate variations. For some metrics such as Northern Hemisphere mean temperatures, there is remarkable agreement between models and reconstructions. For other diagnostics, such as the regional response to volcanic eruptions, or hemispheric temperature differences, substantial disagreements between data and models have been reported. Here, we assess the potential sources of these discrepancies by comparing 1000-year hemispheric temperature reconstructions based on real-world paleoclimate proxies with climate-model-based pseudoproxies. These pseudoproxy experiments (PPE) indicate that noise inherent in proxy records and the unequal spatial distribution of proxy data are the key factors in explaining the data-model differences. For example, lower inter-hemispheric correlations in reconstructions can be fully accounted for by these factors in the PPE. Noise and data sampling also partly explain the reduced amplitude of the response to external forcing in reconstructions compared to models. For other metrics, such as inter-hemispheric differences, some, although reduced, discrepancy remains. Our results suggest that improving proxy data quality and spatial coverage is the key factor to increase the quality of future climate reconstructions, while the total number of proxy records and reconstruction methodology play a smaller role

Bayesian parameter estimation and interpretation for an intermediate model of tree-ring width

We present a Bayesian model for estimating the parameters of the VS-Lite forward model of tree-ring width for a particular chronology and its local climatology. The scheme also provides information about the uncertainty of the parameter estimates, as well as the model error in representing the observed proxy time series. By inferring VS-Lite's parameters independently for synthetically generated ring-width series at several hundred sites across the United States, we show that the algorithm is skillful. We also infer optimal parameter values for modeling observed ring-width data at the same network of sites. The estimated parameter values covary in physical space, and their locations in multidimensional parameter space provide insight into the dominant climatic controls on modeled tree-ring growth at each site as well as the stability of those controls. The estimation procedure is useful for forward and inverse modeling studies using VS-Lite to quantify the full range of model uncertainty stemming from its parameterization

Climate field reconstruction uncertainty arising from multivariate and nonlinear properties of predictors

©2014. American Geophysical Union. All Rights Reserved. Climate field reconstructions (CFRs) of the global annual surface air temperature (SAT) field and associated global area-weighted mean annual temperature (GMAT) are derived in a collection of pseudoproxy experiments for the past millennium. Pseudoproxies are modeled from temperature (T), precipitation (P), T+P, and VS-Lite (VSL), a nonlinear and multivariate proxy system model for tree ring widths. Spatial patterns of reconstruction skill and spectral bias for the T+P and VSL-derived CFRs are similar to those previously shown using temperature-only pseudoproxies but demonstrate overall degraded skill and spectral bias for SAT reconstruction. Analysis of GMAT spectra nevertheless suggests that the true GMAT frequency spectrum is resolved by those pseudoproxies (T, T+P, and VSL) that contain some temperature information. The results suggest that mixed temperature and moisture-responding paleoclimate data may produce actual GMAT reconstructions with skill, error, and spectral characteristics like those expected from univariate and linear temperature responders, but spatially resolved CFR results should be analyzed cautiously. Key Points Uncertainties may be larger than previously estimated from PPEsReconstructed GMAT based on mixed responders is relatively spectrally unbiasedCFRs based on mixed responders should be interpreted with cautionPeer Reviewe

Towards the assimilation of tree-ring-width records using ensemble Kalman filtering techniques

This paper investigates the applicability of the Vaganov–Shashkin–Lite (VSL) forward model for tree-ring-width chronologies as observation operator within a proxy data assimilation (DA) setting. Based on the principle of limiting factors, VSL combines temperature and moisture time series in a nonlinear fashion to obtain simulated TRW chronologies. When used as observation operator, this modelling approach implies three compounding, challenging features: (1) time averaging, (2) “switching recording” of 2 variables and (3) bounded response windows leading to “thresholded response”. We generate pseudo-TRW observations from a chaotic 2-scale dynamical system, used as a cartoon of the atmosphere-land system, and attempt to assimilate them via ensemble Kalman filtering techniques. Results within our simplified setting reveal that VSL’s nonlinearities may lead to considerable loss of assimilation skill, as compared to the utilization of a time-averaged (TA) linear observation operator. In order to understand this undesired effect, we embed VSL’s formulation into the framework of fuzzy logic (FL) theory, which thereby exposes multiple representations of the principle of limiting factors. DA experiments employing three alternative growth rate functions disclose a strong link between the lack of smoothness of the growth rate function and the loss of optimality in the estimate of the TA state. Accordingly, VSL’s performance as observation operator can be enhanced by resorting to smoother FL representations of the principle of limiting factors. This finding fosters new interpretations of tree-ring-growth limitation processes

Tree rings and volcanic cooling

Comment on Mann, M. E., Fuentes, J. D. & Rutherford, S. Nature Geosci. 5, 202–205 (2012)