This package contains code for creation of an "Observational Large Ensemble" of monthly temperature, precipitation, and sea level pressure.
The approach is based upon modeling each variable as a linear combination of a mean, trend, contribution from large-scale modes of variability, and residual climate noise. The ensemble is created through randomization of the latter two components. The time series of the large-scale modes are randomized through application of the Iterative Amplitude Adjusted Fourier Transform. The climate noise is block bootstrapped in time.
The forced component is estimated using the methodology of Dai et al., 2015, Nature Climate Change. Specifically, the observations are regressed against the time series of the global mean, ensemble mean of each variable. The forced trend for sea level pressure is assumed to be zero.
More complete technical and scientific documentation can be found in McKinnon and Deser, 2018, Journal of Climate. However, there are a small number of changes between the methodology documented in McKinnon and Deser (2018) and this code base, which are not currently published. Specifically:
- The data product is monthly, rather than seasonal.
- The seasonality of ENSO variance is included.
- Rather than creating two orthogonal modes from applying PCA the ENSO and PDO, PDO is orthogonalized with respect to ENSO using the Gram-Schmidt orthogonalization procedure.
- The block boostrap size is estimated using the iterative formula of Wilks (1997), Journal of Climate, rather than fixed at two years.
- The AMO time series is lowpass filtered using a cutoff-frequency of 10 years.
The code is designed to be applied either to gridded observational datasets or to the CESM1-LE, which is used as a testbed for the methodology. All file paths for application to the CESM1-LE are with respect NCAR Cheyenne.
Please contact Karen McKinnon (kmckinnon@ucla.edu) if you are using the code or methods in your work.
- Free software: MIT license
This package was created with Cookiecutter and the audreyr/cookiecutter-pypackage project template.