skip to main content

Title: Increasing ocean wave energy observed in Earth’s seismic wavefield since the late 20th century

Ocean waves excite continuous globally observable seismic signals. We use data from 52 globally distributed seismographs to analyze the vertical component primary microseism wavefield at 14–20 s period between the late 1980s and August 2022. This signal is principally composed of Rayleigh waves generated by ocean wave seafloor tractions at less than several hundred meters depth, and is thus a proxy for near-coastal swell activity. Here we show that increasing seismic amplitudes at 3σsignificance occur at 41 (79%) and negative trends occur at 3σsignificance at eight (15%) sites. The greatest absolute increase occurs for the Antarctic Peninsula with respective acceleration amplitude and energy trends ( ± 3σ) of 0.037 ± 0.008 nm s−2y−1(0.36 ± 0.08% y−1) and 4.16 ± 1.07 nm2 s−2y−1(0.58 ± 0.15% y−1), where percentage trends are relative to historical medians. The inferred global mean near-coastal ocean wave energy increase rate is 0.27 ± 0.03% y−1for all data and is 0.35 ± 0.04% y−1since 1 January 2000. Strongly correlated seismic amplitude station histories occur to beyond 50of separation and show regional-to-global associations with El Niño and La Niña events.

more » « less
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Nature Communications
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Kelp beds provide significant ecosystem services and socioeconomic benefits globally, and prominently in coastal zones of the California Current. Their distributions and abundance, however, vary greatly over space and time. Here, we describe long-term patterns of Giant Kelp (Macrocystis pyrifera) sea surface canopy area off the coast of San Diego County from 1983 through 2019 along with recent patterns of water column nitrate (NO3-) exposure inferred fromin situtemperature data in 2014 and 2015 at sites spanning 30 km of the coastline near San Diego California, USA. Site-specific patterns of kelp persistence and resilience were associated with ocean and climate dynamics, with total sea surface kelp canopy area varying approximately 33-fold over the almost 4 decades (min 0.34 km2in 1984; max 11.25 km2in 2008, median 4.79 km2). Site-normalized canopy areas showed that recent kelp persistence since 2014 was greater at Point Loma and La Jolla, the largest kelp beds off California, than at the much smaller kelp bed off Cardiff. NO3-exposure was estimated from an 11-month time series ofin situwater column temperature collected in 2014 and 2015 at 4 kelp beds, using a relationship between temperature and NO3-concentration previously established for the region. The vertical position of the 14.5°C isotherm, an indicator of the main thermocline and nutricline, varied across the entire water column at semidiurnal to seasonal frequencies. We use a novel means of quantifying estimated water column NO3-exposure integrated through time (mol-days m-2) adapted from degree days approaches commonly used to characterize thermal exposures. Water column integrated NO3-exposure binned by quarters of the time series showed strong seasonal differences with highest exposure in Mar - May 2015, lowest exposure in Sep - Dec 2014, with consistently highest exposure off Point Loma. The water column integrated NO3-signal was filtered to provide estimates of the contribution to total nitrate exposure from high frequency variability (ƒ >= 1 cycle 30 hr-1) associated predominantly with internal waves, and low frequency variability driven predominantly by seasonal upwelling. While seasonal upwelling accounted for > 90% of NO3-exposure across the full year, during warm periods when seasonal upwelling was reduced or absent and NO3-exposure was low overall, the proportion due to internal waves increased markedly to 84 to 100% of the site-specific total exposure. The high frequency variability associated with internal waves may supply critical nutrient availability during anomalously warm periods. Overall, these analyses support a hypothesis that differences in NO3-exposure among sites due to seasonal upwelling and higher frequency internal wave forcing contribute to spatial patterns in Giant Kelp persistence in southern California. The study period includes anomalously warm surface conditions and the marine heatwave associated with the “Pacific Warm Blob” superimposed on the seasonal thermal signal and corresponding to the onset of a multi-year decline in kelp canopy area and marked differences in kelp persistence among sites. Our analysis suggests that, particularly during periods of warm surface conditions, variation in NO3-exposure associated with processes occurring at higher frequencies, including internal waves can be a significant source of NO3-exposure to kelp beds in this region. The patterns described here also offer a view of the potential roles of seasonal and higher frequency nutrient dynamics for Giant Kelp persistence in southern California under continuing ocean surface warming and increasing frequency and intensity of marine heatwaves.

    more » « less

    Accurate synthetic seismic wavefields can now be computed in 3-D earth models using the spectral element method (SEM), which helps improve resolution in full waveform global tomography. However, computational costs are still a challenge. These costs can be reduced by implementing a source stacking method, in which multiple earthquake sources are simultaneously triggered in only one teleseismic SEM simulation. One drawback of this approach is the perceived loss of resolution at depth, in particular because high-amplitude fundamental mode surface waves dominate the summed waveforms, without the possibility of windowing and weighting as in conventional waveform tomography.

    This can be addressed by redefining the cost-function and computing the cross-correlation wavefield between pairs of stations before each inversion iteration. While the Green’s function between the two stations is not reconstructed as well as in the case of ambient noise tomography, where sources are distributed more uniformly around the globe, this is not a drawback, since the same processing is applied to the 3-D synthetics and to the data, and the source parameters are known to a good approximation. By doing so, we can separate time windows with large energy arrivals corresponding to fundamental mode surface waves. This opens the possibility of designing a weighting scheme to bring out the contribution of overtones and body waves. It also makes it possible to balance the contributions of frequently sampled paths versus rarely sampled ones, as in more conventional tomography.

    Here we present the results of proof of concept testing of such an approach for a synthetic 3-component long period waveform data set (periods longer than 60 s), computed for 273 globally distributed events in a simple toy 3-D radially anisotropic upper mantle model which contains shear wave anomalies at different scales. We compare the results of inversion of 10 000 s long stacked time-series, starting from a 1-D model, using source stacked waveforms and station-pair cross-correlations of these stacked waveforms in the definition of the cost function. We compute the gradient and the Hessian using normal mode perturbation theory, which avoids the problem of cross-talk encountered when forming the gradient using an adjoint approach. We perform inversions with and without realistic noise added and show that the model can be recovered equally well using one or the other cost function.

    The proposed approach is computationally very efficient. While application to more realistic synthetic data sets is beyond the scope of this paper, as well as to real data, since that requires additional steps to account for such issues as missing data, we illustrate how this methodology can help inform first order questions such as model resolution in the presence of noise, and trade-offs between different physical parameters (anisotropy, attenuation, crustal structure, etc.) that would be computationally very costly to address adequately, when using conventional full waveform tomography based on single-event wavefield computations.

    more » « less
  3. SUMMARY Ocean bottom distributed acoustic sensing (OBDAS) is emerging as a new measurement method providing dense, high-fidelity and broad-band seismic observations from fibre-optic cables deployed offshore. In this study, we focus on 35.7 km of a linear telecommunication cable located offshore the Sanriku region, Japan, and apply seismic interferometry to obtain a high-resolution 2-D shear wave velocity (VS) model below the cable. We first show that the processing steps applied to 13 d of continuous data prior to computing cross-correlation functions (CCFs) impact the modal content of surface waves. Continuous data pre-processed with 1-bit normalization allow us to retrieve dispersion images with high Scholte-wave energy between 0.5 and 5 Hz, whereas spatial aliasing dominates dispersion images above 3 Hz for non-1-bit CCFs. Moreover, the number of receiver channels considered to compute dispersion images also greatly affects the resolution of extracted surface-wave modes. To better understand the remarkably rich modal nature of OBDAS data (i.e. up to 30 higher modes in some regions), we simulate Scholte-wave dispersion curves for stepwise linear VS gradient media. For soft marine sediments, simulations confirm that a large number of modes can be generated in gradient media. Based on pre-processing and theoretical considerations, we extract surface wave dispersion curves from 1-bit CCFs spanning over 400 channels (i.e. ∼2 km) along the array and invert them to image the subsurface. The 2-D velocity profile generally exhibits slow shear wave velocities near the ocean floor that gradually increase with depth. Lateral variations are also observed. Flat bathymetry regions, where sediments tend to accumulate, reveal a larger number of Scholte-wave modes and lower shallow velocity layers than regions with steeper bathymetry. We also compare and discuss the velocity model with that from a previous study and finally discuss the combined effect of bathymetry and shallow VS layers on earthquake wavefields. Our results provide new constraints on the shallow submarine structure in the area and further demonstrate the potential of OBDAS for high-resolution offshore geophysical prospecting. 
    more » « less
  4. Abstract

    Man‐made very low frequency (VLF) transmitter waves play a critical role in energetic electron scattering and precipitation from the inner radiation belt, a type of which is called wisp precipitation. Wisps exhibit dispersive energy‐versus‐Lspectra due to the evolution of electron cyclotron resonance conditions with near‐monochromatic VLF transmitter waves. Here, we report on such observations of inner belt wisp precipitation events with full pitch angle resolution in the energy range of 50 to ∼500 keV as measured by Electron Loss and Fields Investigation (ELFIN) atL < ∼2 between March 2021 and April 2022. Statistical observations (82 events) reveal occasional (18 events) wisp precipitation events with local bounce‐loss‐cone electron flux enhancements, which provide new information compared with flux enhancements measured in previous studies only in the drift loss cone. Based on magnetic field and plasmaspheric density models, quasilinear theory, and detailed pitch angle distributions of wisps from ELFIN, we have estimated the wisp electron bounce‐averaged pitch angle diffusion coefficients to be of the order of 10−4to 10−2 s−1. These are several orders of magnitude larger than the diffusion rates calculated from models using global statistical averages of VLF transmitter wave power. When using our estimated diffusion coefficients to deduce the associated local transmitter wave amplitudes near the equator, based on quasilinear calculations from a transmitter‐induced electron diffusion model, we find these wave amplitudes to be >1 mV/m. Although probable overestimates, such inferred wave amplitudes exceed the theoretical threshold amplitude for nonlinear interactions, strongly suggesting that it is necessary to include nonlinear effects for an accurate evaluation of energetic electron scattering by transmitter waves.

    more » « less
  5. Abstract. The eruption of the Hunga Tonga-Hunga Ha'apai volcano on 15 January 2022 provided a rare opportunity to understand global tsunamiimpacts of explosive volcanism and to evaluate future hazards, includingdangers from “volcanic meteotsunamis” (VMTs) induced by the atmosphericshock waves that followed the eruption. The propagation of the volcanic andmarine tsunamis was analyzed using globally distributed 1 min measurementsof air pressure and water level (WL) (from both tide gauges and deep-waterbuoys). The marine tsunami propagated primarily throughout the Pacific,reaching nearly 2 m at some locations, though most Pacific locationsrecorded maximums lower than 1 m. However, the VMT resulting from theatmospheric shock wave arrived before the marine tsunami and propagatedglobally, producing water level perturbations in the Indian Ocean, theMediterranean, and the Caribbean. The resulting water level response of manyPacific Rim gauges was amplified, likely related to wave interaction withbathymetry. The meteotsunami repeatedly boosted tsunami wave energy as itcircled the planet several times. In some locations, the VMT was amplifiedby as much as 35-fold relative to the inverse barometer due to near-Proudmanresonance and topographic effects. Thus, a meteotsunami from a largereruption (such as the Krakatoa eruption of 1883) could yield atmosphericpressure changes of 10 to 30 mb, yielding a 3–10 m near-field tsunami thatwould occur in advance of (usually) larger marine tsunami waves, posingadditional hazards to local populations. Present tsunami warning systems donot consider this threat. 
    more » « less