skip to main content


Title: Rotated Spectral Principal Component Analysis (rsPCA) for Identifying Dynamical Modes of Variability in Climate Systems
Abstract Spectral PCA (sPCA), in contrast to classical PCA, offers the advantage of identifying organized spatiotemporal patterns within specific frequency bands and extracting dynamical modes. However, the unavoidable trade-off between frequency resolution and robustness of the PCs leads to high sensitivity to noise and overfitting, which limits the interpretation of the sPCA results. We propose herein a simple nonparametric implementation of sPCA using the continuous analytic Morlet wavelet as a robust estimator of the cross-spectral matrices with good frequency resolution. To improve the interpretability of the results, especially when several modes of similar amplitude exist within the same frequency band, we propose a rotation of the complex-valued eigenvectors to optimize their spatial regularity (smoothness). The developed method, called rotated spectral PCA (rsPCA), is tested on synthetic data simulating propagating waves and shows impressive performance even with high levels of noise in the data. Applied to global historical geopotential height (GPH) and sea surface temperature (SST) daily time series, the method accurately captures patterns of atmospheric Rossby waves at high frequencies (3–60-day periods) in both GPH and SST and El Niño–Southern Oscillation (ENSO) at low frequencies (2–7-yr periodicity) in SST. At high frequencies the rsPCA successfully unmixes the identified waves, revealing spatially coherent patterns with robust propagation dynamics.  more » « less
Award ID(s):
1928724 1839336 1839441 1735040
NSF-PAR ID:
10209943
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Journal of Climate
Volume:
34
Issue:
2
ISSN:
0894-8755
Page Range / eLocation ID:
715 to 736
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Various high-frequency waves in the vicinity of upper-hybrid and Langmuir frequencies are commonly observed in different space plasma environments. Such waves and fluctuations have been reported in the magnetosphere of the Earth, a planet with an intrinsic strong magnetic field. Mars has no intrinsic magnetic field and, instead, it possesses a weak induced magnetosphere, which is highly dynamic due to direct exposure to the solar wind. In the present paper, we investigate the presence of high-frequency plasma waves in the Martian plasma environment by making use of the high-resolution electric field data from the Mars Atmosphere and Volatile Evolution missioN (MAVEN) spacecraft. Aims. This study aims to provide conclusive observational evidence of the occurrence of high-frequency plasma waves around the electron plasma frequency in the Martian magnetosphere. We observe two distinct wave modes with frequency below and above the electron plasma frequency. The characteristics of these high-frequency waves are quantified and presented here. We discuss the generation of possible wave modes by taking into account the ambient plasma parameters in the region of observation. Methods. We have made use of the medium frequency (100 Hz–32 kHz) burst mode-calibrated electric field data from the Langmuir Probe and Waves instrument on board NASA’s MAVEN mission. Due to the weak magnetic field strength, the electron gyro-frequency is much lower than the electron plasma frequency, which implies that the upper-hybrid and Langmuir waves have comparable frequencies. A total of 19 wave events with wave activities around electron plasma frequency were identified by examining high-resolution spectrograms of the electric field. Results. These waves were observed around 5 LT when MAVEN crossed the magnetopause boundary and entered the magnetosheath region. These waves are either a broadband- or narrowband-type with distinguishable features in the frequency domain. The narrowband-type waves have spectral peak above the electron plasma frequency. However, in the case of broadband-type waves, the spectral peak always occurred below the electron plasma frequency. The broadband waves consistently show a periodic modulation of 8–14 ms. Conclusions. The high-frequency narrowband-type waves observed above the electron plasma frequency are believed to be associated with upper-hybrid or Langmuir waves. However, the physical mechanism responsible for the generation of broadband-type waves and the associated 8–14 ms modulation remain unexplained and further investigation is required. 
    more » « less
  2. Abstract

    We present observations and modeling of spatial eigen‐functions of resonating waves within fault zone waveguide, using data recorded on a dense seismic array across the San Jacinto Fault Zone (SJFZ) in southern California. The array consists of 5‐Hz geophones that cross the SJFZ with ~10–30 m spacing at the Blackburn Saddle near the Hemet Stepover. Wavefield snapshots after theSwave arrival are consistent for more than 50 near‐fault events, suggesting that this pattern is controlled by the fault zone structure rather than source properties. Data from example event with high signal to noise ratio show three main frequency peaks at ~1.3, ~2.0, and ~2.8 Hz in the amplitude spectra of resonance waves averaged over stations near the fault. The data are modeled with analytical expressions for eigen‐functions of resonance waves in a low‐velocity layer (fault zone) between two quarter‐spaces. Using a grid search‐based method, we investigate the possible width of the waveguide, location within the array, and shear wave velocities of the media that fit well the resonance signal at ~1.3 Hz. The results indicate a ~300 m wide damaged fault zone layer with ~65%Swave velocity reduction compared to the host rock. The SW edge of the low‐velocity zone is near the mapped fault surface trace, indicating that the damage zone is asymmetrically located at the regionally faster NE crustal block. The imaging resolution of the fault zone structure can be improved by modeling fault zone resonance modes and trapped waves together.

     
    more » « less
  3. Abstract

    We develop an automatic workflow for enhancing surface wave signals in ambient noise cross correlations (ANCs) calculated for a one‐dimensional (1‐D) linear array. The proposed array‐based method is applied to a 1.6 km‐long dense linear nodal array crossing surface traces of the San Jacinto fault near Anza, California. Fundamental and higher modes of surface waves are observed in ANCs of the nodal array. After attenuating the surface wave overtones by applying a frequency‐dependent tapering window to the ANCs, signals dominated by the fundamental mode surface wave are then enhanced through a denoising process based on three‐station interferometry of direct waves. The signal‐to‐noise ratio is significantly increased at high frequencies (>2 Hz) after denoising. Phase travel times are extracted reliably in the frequency domain for the period ranges of 0.3–1.2 s and 0.3–1.6 s for Rayleigh and Love waves, respectively. The corresponding period‐dependent phase velocity profiles derived from the eikonal equation reveal high‐resolution details of fault zone internal structures beneath the array. A broad (500–1,000 m) low‐velocity zone that narrows with increasing period is observed, illuminating a flower‐shaped structure of the San Jacinto fault damage zone.

     
    more » « less
  4. null (Ed.)
    Abstract Cook Inlet fore‐arc basin in south‐central Alaska is a large, deep (7.6 km) sedimentary basin with the Anchorage metropolitan region on its margins. From 2015 to 2017, a set of 28 broadband seismic stations was deployed in the region as part of the Southern Alaska Lithosphere and Mantle Observation Network (SALMON) project. The SALMON stations, which also cover the remote western portion of Cook Inlet basin and the back‐arc region, form the basis for our observational study of the seismic response of Cook Inlet basin. We quantify the influence of Cook Inlet basin on the seismic wavefield using three data sets: (1) ambient‐noise amplitudes of 18 basin stations relative to a nonbasin reference station, (2) earthquake ground‐motion metrics for 34 crustal and intraslab earthquakes, and (3) spectral ratios (SRs) between basin stations and nonbasin stations for the same earthquakes. For all analyses, we examine how quantities vary with the frequency content of the seismic signal and with the basin depth at each station. Seismic waves from earthquakes and from ambient noise are amplified within Cook Inlet basin. At low frequencies (0.1–0.5 Hz), ambient‐noise ratios and earthquake SRs are in a general agreement with power amplification of 6–14 dB, corresponding to amplitude amplification factors of 2.0–5.0. At high frequencies (0.5–4.0 Hz), the basin amplifies the earthquake wavefield by similar factors. Our results indicate stronger amplification for the deeper basin stations such as near Nikiski on the Kenai Peninsula and weaker amplification near the margins of the basin. Future work devoted to 3D wavefield simulations and treatment of source and propagation effects should improve the characterization of the frequency‐dependent response of Cook Inlet basin to recorded and scenario earthquakes in the region. 
    more » « less
  5. Being able to image chemical bonds with high sensitivity and speed, stimulated Raman scattering (SRS) microscopy has made a major impact in biomedical optics. However, it is well known that the standard SRS microscopy suffers from various backgrounds, limiting the achievable contrast, quantification and sensitivity. While many frequency-modulation (FM) SRS schemes have been demonstrated to retrieve the sharp vibrational contrast, they often require customized laser systems and/or complicated laser pulse shaping or introduce additional noise, thereby hindering wide adoption. Herein we report a simple but robust strategy for FM-SRS microscopy based on a popular commercial laser system and regular optics. Harnessing self-phase modulation induced self-balanced spectral splitting of picosecond Stokes beam propagating in standard single-mode silica fibers, a high-performance FM-SRS system is constructed without introducing any additional signal noise. Our strategy enables adaptive spectral resolution for background-free SRS imaging of Raman modes with different linewidths. The generality of our method is demonstrated on a variety of Raman modes with effective suppressing of backgrounds including non-resonant cross phase modulation and electronic background from two-photon absorption or pump-probe process. As such, our method is promising to be adopted by the SRS microscopy community for background-free chemical imaging.

     
    more » « less