The sampling method for optimal precursors of El Niño–Southern Oscillation events
<p>The El Niño–Southern Oscillation (ENSO) is a significant climate phenomenon that appears periodically in the tropical Pacific. The intermediate coupled ocean–atmosphere Zebiak–Cane (ZC) model is the first and classical one designed to numerically forecast the ENSO events. Traditionally, th...
| Main Authors: | , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Copernicus Publications
2024-03-01
|
| Series: | Nonlinear Processes in Geophysics |
| Online Access: | https://npg.copernicus.org/articles/31/165/2024/npg-31-165-2024.pdf |
| Summary: | <p>The El Niño–Southern Oscillation (ENSO) is a significant climate phenomenon that appears periodically in the tropical Pacific. The intermediate coupled ocean–atmosphere Zebiak–Cane (ZC) model is the first and classical one designed to numerically forecast the ENSO events. Traditionally, the conditional nonlinear optimal perturbation (CNOP) approach has been used to capture optimal precursors in practice. In this paper, based on state-of-the-art statistical machine learning techniques<span class="note-anchor" id="fna_Ch1.Footn1"><a href="#fn_Ch1.Footn1"><sup>1</sup></a></span>, we investigate the sampling algorithm proposed in <span class="cit" id="xref_text.1"><a href="#bib1.bibx36">Shi and Sun</a> (<a href="#bib1.bibx36">2023</a>)</span> to obtain optimal precursors via the CNOP approach in the ZC model. For the ZC model, or more generally, the numerical models with a large number <span class="inline-formula">O(10<sup>4</sup>−10<sup>5</sup>)</span> of degrees of freedom, the numerical performance, regardless of the statically spatial patterns and the dynamical nonlinear time evolution behaviors as well as the corresponding quantities and indices, shows the high efficiency of the sampling method compared to the traditional adjoint method. The sampling algorithm does not only reduce the gradient (first-order information) to the objective function value (zeroth-order information) but also avoids the use of the adjoint model, which is hard to develop in the coupled ocean–atmosphere models and the parameterization models. In addition, based on the key characteristic that the samples are independently and identically distributed, we can implement the sampling algorithm by parallel computation to shorten the computation time. Meanwhile, we also show in the numerical experiments that the important features of optimal precursors can still be captured even when the number of samples is reduced sharply.</p> |
|---|---|
| ISSN: | 1023-5809 1607-7946 |