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 more » spatially coherent patterns with robust propagation dynamics. « less
; ; ; ; ;
Award ID(s):
1928724 1839336 1839441 1735040
Publication Date:
Journal Name:
Journal of Climate
Page Range or eLocation-ID:
715 to 736
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 faultmore »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.

    « less
  2. 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.

  3. 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.more »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.« less
  4. 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.

  5. Abstract

    The critical zone sustains terrestrial life, but we have few tools to explore it efficiently beyond the first few meters of the subsurface. Using analyses of high‐frequency ambient seismic noise from densely spaced seismometers deployed in the forested Shale Hills subcatchment of the Susquehanna Shale Hills Critical Zone Observatory (SSHCZO), we show that temporal changes in seismic velocities at depths from ∼1 m to tens of m can be detected. These changes are driven by variations at the land surface. The Moving‐Window Cross‐Spectral (MWCS) method was employed to measure seismic‐velocity changes in coda waves at hourly resolution in 10 different frequency bands. We observed a diurnal signal, a seasonal signal, and a meteorological‐event‐based signal. These signals were compared to time‐series measurements of precipitation, well water levels, soil moisture, soil temperature, air temperature, latent heat flux, and air pressure in the heavily instrumented catchment. Most of the velocity changes can be explained by variations in temperature that result in thermoelastic strains that propagate to depth. But some double minima in seismic velocity time‐series observed after large rain events were attributed in part to the effects of water infiltration. These results show that high‐frequency ambient noise data may in some locations bemore »used to detect changes in the critical zone from ∼1 to ∼100 m or greater depth with hourly resolution. But interpretation of such data requires multiple environmental data sets to deconvolve the complex interrelationships among thermoelastic and hydrological effects in the subsurface critical zone.

    « less