Hostname: page-component-cd9895bd7-7cvxr Total loading time: 0 Render date: 2024-12-23T01:11:19.638Z Has data issue: false hasContentIssue false
  • English
  • Français

The singing firn

Published online by Cambridge University Press:  29 June 2023

Julien Chaput ,
Richard C. Aster  and
Marianne Karplus
Julien Chaput*
Affiliation:
Department of Earth Environmental and Resource Sciences, University of Texas at El Paso, El Paso, TX, USA
Richard C. Aster
Affiliation:
Department of Geosciences and Warner College of Natural Resources, Colorado State University, Fort Collins, CO, USA
Marianne Karplus
Affiliation:
Department of Earth Environmental and Resource Sciences, University of Texas at El Paso, El Paso, TX, USA
*
Corresponding author: Julien Chaput; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Antarctic firn presents an exotic seismological environment in which the behaviors of propagating waves can be significantly at odds with those in other Earth media. We present a condensed view of the nascent field of ambient noise seismology in Antarctic firn-covered media, and highlight multiple unusual and information-rich observations framed through the lens of the firn's important role as a buffer for air temperature anomalies and a complex contributor to ice mass balance. We summarize key results from several recent papers depicting novel wind-excited firn resonances and point to the plethora of ways these observations could facilitate imaging and monitoring of glacial systems at single, isolated seismometers. Finally, we propose significant instrumental and computational objectives necessary to constrain resonance excitation mechanisms and broadly apply these observations as useful monitoring tools in Antarctica.

Keywords

Anisotropic iceseismologysnow physics
Type
Letter
Information
Annals of Glaciology , Volume 63 , Issue 87-89 , September 2022 , pp. 137 - 142
Check for updates
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2023. Published by Cambridge University Press on behalf of The International Glaciological Society

1. Introduction

Firn is broadly defined by the gradual transition from loose surface snow to solid ice through compaction, densification, pore closure and other effects, and is most often the uppermost structure for large glacial systems in polar regions. In Antarctica, firn covers approximately 99% of all glaciers (van den Broeke, Reference van den Broeke2008; Ligtenberg and others, Reference Ligtenberg, Helsen and van den Broeke2011) and is both an integrated component of ice masses and a somewhat separate, exotic medium with significant structural variability. Accurately estimating the firn density profile has long been a primary objective of the glaciological community, but current models (e.g. Stevens and others, Reference Stevens2020) are unable to account for the swath of local effects that cause deviations from average assumptions. Ice flow and its related strain environment, for instance, have recently been shown to strongly affect layer density through settling (Horlings and others, Reference Horlings, Christianson, Holschuh, Stevens and Waddington2021; Oraschewski and Grinsted, Reference Oraschewski and Grinsted2022). Effects related to environmental surface forcing, such as temperature changes and wind/deposition interactions, (e.g. Reeh and others, Reference Reeh, Fisher, Koerner and Clausen2005; Reeh, Reference Reeh2008) can furthermore cause large perturbations in firn layering and density away from an assumed smooth gradient from snow to solid ice. Both effects cause significant uncertainties in firn profile estimates.

Beyond its global contribution as a challenging component of ice mass-balance estimates, the inherent porosity, parametric gradient and dynamic nature of the firn allow it to absorb environmental forcing in multiple ways, including pore space retention and refreezing of surface melt (Rennermalm and others, Reference Rennermalm2013; Steger, Reference Steger2017; Vandecrux and others, Reference Vandecrux2020). In some cases, particularly with respect to ice shelves, the progressive loss of the firn can trigger catastrophic shelf failure due to melt ponding and hydrofracture (Kuipers and others, Reference Kuipers Munneke, Ligtenberg, van den Broeke and Vaughan2017), and ablation and reduction in albedo (Scambos and others, Reference Scambos, Bohlander, Shuman and Skvarca2004; Leppäranta and others, Reference Leppäranta, Järvinen and Mattila2012; Banwell, Reference Banwell2017; Kuipers and others, Reference Kuipers Munneke, Ligtenberg, van den Broeke and Vaughan2017; MacAyeal, Reference MacAyeal2018), as was the case for Larsen B Ice Shelf, resulting in accelerated ice flow across the grounding line following its 2002 collapse (Rignot and others, Reference Rignot, Casassa, Gogineni, Krabill, Rivera and Thomas2004).

Passive correlation-based seismic methods, which are widely applied to study structural temporal variability for seismic velocity and scattering properties, require the deployment of station arrays to construct interstation noise correlation functions (NCFs) through an approach called seismic interferometry (Campillo and Paul, Reference Campillo and Paul2003; Snieder, Reference Snieder2004; Wapenaar and Fokkema, Reference Wapenaar and Fokkema2006; Wapenaar and others, Reference Wapenaar, Draganov, Snieder, Campman and Verdel2010). Such multi-station methods work best under conditions of ambient source stability and implement significant time averaging to reconstruct interpretable NCFs. Potentially large errors on the phase of reconstructed surface waves can be introduced when the noise source is not temporally and spatially stable, and this often precludes the use of higher ambient noise frequencies suitable for near-surface structures like the firn. Seismic inversions based on the surface wave components of the NCFs (e.g. Diez and others, Reference Diez2016) are thus mostly insensitive to small-scale near-surface parametric contrasts such as ice lenses, hoarfrost layers and other embedded shallow features due to the averaging nature of the surface wave depth kernels at lower frequencies most often used in ambient seismic noise studies.

Although seismic interferometry has been leveraged successfully in a number of cryospheric studies (e.g. Diez and others, Reference Diez2016; Mordret and others, Reference Mordret, Mikesell, Harig, Lipovsky and Prieto2016; Aster, Reference Aster2019; Zhan, Reference Zhan2019) we focus here on spectral domain observations made at widely distributed individual isolated stations in Antarctica. Pervasive observations of high-frequency ambient spectral resonances at Antarctic seismic stations (Chaput and others, Reference Chaput2018) show a number of features that are relevant to firn structure and evolution. Wind-excited spectral peaks, termed firn resonances, manifest as patterns of sparse, spectral amplifications above 5 Hz that respond strongly to environmental forcing phenomena such as storms and temperature anomalies, and are demonstrably sensitive to depth-dependent medium parameters such as anisotropy and layering. Here, we summarize three recent forays exploring these observations (Chaput and others, Reference Chaput2018, Reference Chaput, Aster, Karplus and Nakata2022a, Reference Chaput2022b) and emphasize the potential for significant information retrieval at single seismic stations deployed in firn media. We further elaborate on directions of study involving constraints on excitation physics that would allow these novel observations to be invertible quantities.

2. Firn resonances

Chaput and others (Reference Chaput2018) first noted the presence of narrow band peaks in spectrograms of ambient seismic data on the Ross Ice Shelf (RIS; Bromirski and others, Reference Bromirski2015, Fig. 1A), inferred to be excited by wind forcing. Such resonances have since been observed at other firn-covered locales including at the West Antarctic Ice Sheet (WAIS) Divide and South Pole, with varying types of instrumentation including completely snow-buried instruments with low to zero wind profile. Firn spectral peaks feature complex behaviors, including frequency shifts on the order of hours following strong wind events, response to surface softening or melt (e.g. Nicolas and others, Reference Nicolas2017), multi-month drifts in peak frequency patterns, harmonic resonance patterns with broadband coherent drift or, conversely, behavior where multiple peaks shift independently of each other (Fig. 1A, basic forcing effects shown in Fig. 1B). These narrow band spectral peak patterns and their compelling spectrogram sonifications have further sparked interest from members of the arts community (e.g. Canadian audiovisual artist Sandra Volny and Emmy award winning composer Lucas Cantor, among others) who are developing multifaceted interpretive projects. The information content of firn resonances is surprising, particularly when one considers that observations are performed at single stations. We review primary results from three recent papers on the subject (Chaput and others, Reference Chaput2018, Reference Chaput, Aster, Karplus and Nakata2022a, Reference Chaput2022b), frame them in the context of broader community knowledge gaps and propose directions for future studies aiming to leverage sparse, single seismic stations for imaging and temporal monitoring efforts in firn media.

Figure 1. (A) Example of firn resonances from nearly two years of North component ambient seismic data recorded at station DR09 on Roosevelt Island from Ross Ice Shelf broadband array (red circle, inset map) displayed as a time/frequency plot (spectrogram) (Chaput and others, Reference Chaput2018). Red boxes indicate seasonal open sea ice conditions, and corresponding spectral effects are not observed at grounded sites. Stable shelf plate modes are visible as high amplitude temporally stable vertical bands below 5 Hz (notably during open sea ice conditions, red boxes), and temporally variable firn modes above roughly 5 Hz (observable year round). (B) Description of environmental effects that dictate the behavior of firn resonances. (1) Wind coupling with semi-periodic surface snowforms and the low-velocity/density firn structure excites unique firn mode patterns. (2) Firn is sensitive to anomalous near-zero surface temperatures, and the frequency range over which resonances are altered depends on the depth penetration of the temperature anomaly. (3) Firn accommodates strain associated with flowing ice masses in a ductile fashion at shallow depths where porosity is high, and in a brittle fashion where pores have largely closed.

3. Boundary layer monitoring

Chaput and others (Reference Chaput2018) noted that firn resonances are responsive to atmospheric boundary layer processes, including surface snowform alterations following waning storms (e.g. Sommer and others, Reference Sommer, Wever, Fierz and Lehning2018), temperature fluctuations near the melting point and long-term (i.e. months to years) peak frequency decay and drift hypothesized to be related to firn compaction. Figure 2A demonstrates that some storms are capable of dramatically shifting the frequency content of resonance patterns with their passing, pointing to a direct involvement of surface snowforms (e.g. sastrugi) on the source mechanism responsible for resonance generation. Indeed, 2D numerical wavefield simulations (Chaput and others, Reference Chaput2018) have shown that by changing the spatial periodicity of surface sources as a proxy for wind coupling, different frequencies can be naturally amplified. The very low seismic velocities of firn caused by high porosity also tend to drive these amplifications to overall lower, and hence observable, parts of the spectrum, compared to what we might expect in normal Earth media. Particularly strong storms can overcome surface grain sintering and alter snowform distributions (e.g. Sharma and others, Reference Sharma, Braud and Lehning2019), as can storms with high airborne snow budgets (i.e. with deposition effects). If relative calm follows such a storm, the new spectral pattern will often slowly decay back to its original state over the space of months (e.g. Fig. 1A, black boxes), suggesting a sensitivity to steady-state surface erosion and compaction processes.

Figure 2. Examples of firn resonances show in ‘peak tracked’ form, where only peak maximums for both horizontal components of the spectrogram at each time bin are shown. (A) Firn resonance response to alterations in surface snowforms. Strong storm activity (arrows in left panel, matched with periods of high winds shown by the black trace), can deposit new snowforms that are then slowly eroded during periods of quiescence (right panel), resulting in slow spectral decay (time scale of months, black dotted boxes). (B) Peak tracked firn resonances at 5 stations on RIS during a shelf-wide near-zero temperature event in 2016 (bottom left panel, Nicolas and others (Reference Nicolas2017)). Firn undergoes up to a 40% reduction in elastic moduli (top right panel) as bonds between snow grains weaken, resulting in a downward drift in frequency for higher peaks and a reduction in amplitude (bottom right panel). (C) Peak tracked spectrogram at RIS station RS17 for 10 days during 2016, showing the obvious offset in frequency (also frequency dependent) between North and East components. Right panel: Shallow firn deforms plastically under extensional strain typical of RIS and features a strain-elongated pore space (black arrows aligned with ice flow), while deeper firn to solid ice responds in a brittle manner, often resulting in flow-perpendicular crevassing (red arrows) unless dominant crevassing is advected and rotated from past strain regimes.

Furthermore, firn resonances are highly sensitive to surface temperatures as they approach melting, without necessarily even crossing the threshold into meltwater generation, as for example observed during an extended period of near-zero temperatures on the Ross Ice Shelf in 2016 (Nicolas and others, Reference Nicolas2017), and shown in Figure 2B. For all stations within the event area, the frequency content of higher frequency peaks drifted downward and fell in amplitude, hitting a minimum after 3–4 days, and partially recovered when a subsequent cold snap occurred. Takei and Maeno (Reference Takei and Maeno2004) showed that snow undergoes up to a 40% reduction in elastic moduli as temperatures approach zero without even necessarily generating melt, pointing a direct link between increasing temperatures, decreasing seismic velocities and decreasing frequency content. The insensitivity of lower frequency firn modes (i.e. 5–10 Hz, shown in Chaput and others (Reference Chaput2018)) to this event were physically interpreted through a surface-driven thermal diffusion model as noted in other snow studies (e.g. Gilbert and othets, Reference Gilbert, Vincent, Six, Wagnon, Piard and Ginot2014), where surface temperature anomalies without melt only reached a limited depth in the firn. Given strong evidence that firn resonances are related to surface wave excitation (Chaput and others, Reference Chaput, Aster, Karplus and Nakata2022a) with frequency-dependent depth sensitivities (i.e. lower frequency Rayleigh waves are on average sensitive to deeper structures than higher frequencies), firn resonances can be used to evaluate the depth penetration of temperature-related atmospheric forcing.

4. Constraining models of firn structure

As mentioned above, one of the most daunting barriers to accurately modeling firn density profiles lies in estimating fluctuations away from steady-state densification models. This encompasses, for example, elusive effects related to surface temperature forcing (e.g. Reeh and others, Reference Reeh, Fisher, Koerner and Clausen2005; Reeh, Reference Reeh2008) and constraining the firn's settling behavior under different strain regimes (Horlings and others, Reference Horlings, Christianson, Holschuh, Stevens and Waddington2021; Oraschewski and Grinsted, Reference Oraschewski and Grinsted2022). While studying past, and thus buried, strain effects and imaging fine layering due to ice lens and melt layers remain difficult problems, firn resonances offer potential avenues of study with the added benefit that the necessary observations can be performed on single sensors. Firn resonances present several interesting quantities that are at least partially invertible. Firstly, the spectral patterns themselves may offer constraints on firn structure, as shown by Bayesian explorations of resonances for 2D models Chaput and others (Reference Chaput2018), with the caveat that resonance peaks are a combination of both surface source distributions and firn structure. Chaput and others (Reference Chaput, Aster, Karplus and Nakata2022a) furthermore showed that resonance peak patterns are indeed affected by local structures, and their frequency content follows similar spatial variation trends to several other well-known site response metrics associated with Rayleigh waves propagating in strong parametric gradients, such as the widely used H/V ratio (e.g. Nakamura, Reference Nakamura1989) and Rayleigh wave particle motions patterns (e.g. Tanimoto and Rivera, Reference Tanimoto and Rivera2005; Denolle and others, Reference Denolle, Dunham and Beroza2012; Berbellini and others, Reference Berbellini, Morelli and Ferreira2016). Given that overhead satellite imagery offers the potential for estimating surface snowform distributions, firn resonances could be coupled with these other metrics in a joint inversion of firn profiles (particularly with H/V, since it is another single station measurement).

Chaput and others (Reference Chaput2022b) noted that spectral patterns can be mined for another interesting parameter set, as they almost universally display a frequency offset between the seismometer's orthogonal horizontal components (referred to here as ‘peak splitting’) that can be interpreted in the context of azimuthal anisotropy. This link was confirmed with active sources at WAIS Divide as part of the TIME project (Chaput and others, Reference Chaput2022b). Azimuthal anisotropy from firn resonances was interpreted as being governed at greater depth and lower frequencies (<~25 Hz) by remote-sensing visible advected crevasses in the ice governed by strains imparted through accelerating flow (Ledoux and others, Reference Ledoux, Hulbe, Forbes, Scambos and Alley2017), and at shallow depths and higher frequencies (>~25 Hz) by plastic elongation of the pore space in the shallow firn. Although this latter mechanism has not been directly observed in snow, it has been widely studied in materials engineering (e.g. Melon and others, Reference Melon, Lafarge, Castagnede and Brown1995, Reference Melon, Mariez, Ayrault and Sahraoui1998; Tita and Caliri Junior, Reference Tita and Caliri Junior2012) and medical physics (e.g. Hosokawa and Otani, Reference Hosokawa and Otani1998; Lee and others, Reference Lee, Hughes, Humphrey, Leighton and Choi2007) in terms of anisotropic properties of open-celled foams. For snow, this results in fast anisotropic directions that are aligned with ice flow (i.e. maximum extension) at higher frequencies and with crevassing at lower frequencies (Fig. 2C, right panel).

That being said, mapping these splitting observations to exact depths is a complex problem. Chaput and others (Reference Chaput, Aster, Karplus and Nakata2022a, Reference Chaput2022b) numerically showed that Rayleigh waves propagating in realistic firn media (i.e. with strong shallow deviations away from a smooth densification model similar to those modeled by Reeh (Reference Reeh2008)) will, at certain frequencies determined by fluctuations in structure, have their sensitivity become extremely focused at specific depths as opposed to smoothly distributed (Tanimoto and Rivera, Reference Tanimoto and Rivera2005; Haney and Tsai, Reference Haney and Tsai2015). Thus, although it is clear that the transition between ductile and brittle strain accommodation in the firn occurs roughly at the same frequency for most seismic sites on RIS, suggesting a physical generality (Chaput and others, Reference Chaput2022b), it is unclear what that depth might be beyond conjecture or simple assumptions of smooth depth sensitivity. In the latter case, however, fundamental mode Rayleigh sensitivity kernels indicate a likely transition between 10 and 20 m for a firn profile derived by Diez and others (Reference Diez2016) for a dense array on RIS. Passive anisotropy measurements in firn settings describing a depth at which strain accommodation switches from ductile to brittle is an attractive goal, given that it describes a new form of depth transition in density that can be leveraged in profile estimations.

In light of the direct and physically justifiable causation between firn resonances and both structural and temporally variable metrics, there is a strong impetus for developing further physical models that reach beyond qualitative inferences. This push will require focused and interdisciplinary experiments.

5. Future work and directions

Although clear temporal and structural data products have been constructed from firn resonances through meticulous comparisons with other datasets, there remain multiple questions pertaining to the full physics that excite, propagate and induce temporal variations in firn mode frequency. A high-dimensional parameter space of cause and effect is expected here, and a commensurately focused multi-scale cross-disciplinary experiment should be employed, with a downstream goal of clarifying and interpreting these phenomena. We thus propose that the emerging field of cryoseismology (Podolskiy and Walter, Reference Podolskiy and Walter2016; Aster and Winberry, Reference Aster and Winberry2017) would greatly benefit from a dense multifaceted and sufficiently long-term experiment aiming to robustly constrain the seismic behavior of Antarctic firn. An experiment aiming to constrain the finer points of firn seismology should ultimately be able to document the following aspects of the Antarctic firn environment: (1) snowform topographic variability and its relation to the ambient seismic source, (2) the impact of strong near-surface layering and other structure on resonance patterns, (3) the types of seismic waves responsible for firn resonance observations, (4) 3D spatial variability of the resonance peaks with respect to ice cores and local structure imaged via other means and (5) influences of environmental forcing factors (e.g. temperature, wind strength and history, wind shear, wind direction, humidity, atmospheric pressure and depositional and stripping history) on the firn wavefield.

A concept sketch for such an experiment is depicted in Figure 3A. Wavefield separation into P and S components requires the calculation of the 3D wavefield gradient and curl, which in turn requires a 3D array of conventional three-component seismic instruments, rotational sensors or both (Schmelzbach and others, Reference Schmelzbach2018). Mapping variability in surface structure (e.g. dunes and sastrugi) requires altimetry or photogrammetry methods (or both), and the ability to track changes over time. Assessment of influences due to any relevant above-snow mechanical instrumentation resonances requires on-instrument accelerometers (Qin and others, Reference Qin, Qiu, Nakata, Deng, Levander and Ben-Zion2022). Characterization of environmental forcing and surface topography requires dedicated weather stations, such as the already long running Antarctic Automatic Weather Stations (AWS) Project (Lazzara and others, Reference Lazzara, Weidner, Keller, Thom and Cassano2012) and optical camera, LIDAR or laser altimeter surveys. Accurately constraining the near surface velocity model and layering is also a key component of reducing parametric complexity in source effects. A dense nodal seismograph deployment combined with Distributed Acoustic Sensing (DAS) fiber optic strain rate, and snow core analysis (as an ancillary product of installing borehole seismometers, broadly supported by the US Antarctic Drilling Program) would provide directly sampled medium constraints. Finally, seismic modeling and inversion should be facilitated by a numerical model capable of replicating resonance patterns and other high-frequency seismic observables. For this, we require a framework capable of implementing a full 3D anisotropic velocity model with surface topography and distributed surface sources, such as SPECFEM3D (Komatitsch and Tromp, Reference Komatitsch and Tromp2002).

Figure 3. Multi-faceted seismic and distributed acoustic sensing (DAS) experiment coupled with drone-based photogrammetry aiming to make key firn seismic and environmental measurements to advance understanding of the firn medium and environmental forcing effects that govern resonance peak and other seismological observations, leveraging a long-running UW-Madison (UWM) autonomous weather or similar station (e.g. Lazzara and others, Reference Lazzara, Weidner, Keller, Thom and Cassano2012).

Acknowledgements

The WAIS Divide data are from the Thwaites Interdisciplinary Margin Evolution (TIME) project of the NSF-NERC International Thwaites Glacier Collaboration. The RIS data were collected under NSF-PLR 1141916 and are archived at the EarthScope Data Management Center. The authors also acknowledge support by the National Science Foundation from awards OPP-1739027 and OPP-1744852. The facilities of the EarthScope Consortium were used for access to waveforms and related metadata used in this study. These services are funded through the Seismological Facility for the Advancement of Geoscience (SAGE) Award of the National Science Foundation under Cooperative Support Agreement EAR-1851048. We thank Galen Kaip, Reinhard Flick, and Patrick Shore for their technical support during field work. We thank the US Antarctic Program and the WAIS Divide Camp and support staff during the field work. This work was also partially funded by the University of Texas, El Paso startup funds (JC). Critical data collection from the Ross Ice Shelf and West Antarctica were supported by the National Science Foundation (NSF) Grant Numbers PLR-1142518, 1141916, 1142126, 1246151, 1246416, 1853896, 1142126, 1142518, 1148982, 1246151, 1249631, 1249602, 1249513, 1246666, 1246712, 1246776 and 1247518.)

References

Aster, RC (2019) Interrogating a surging glacier with seismic interferometry. Geophysical Research Letters 46, 8162–8165. doi:10.1029/2019gl084286Google Scholar
Aster, RC and Winberry, JP (2017) Glacial seismology. Reports On Progress in Physics 80, 126801. doi:10.1088/1361-6633/aa8473Google Scholar
Banwell, A (2017) Ice-shelf stability questioned. Nature 544, 306–307. doi:10.1038/544306aGoogle Scholar
Berbellini, A, Morelli, A and Ferreira, AMG (2016) Ellipticity of Rayleigh waves in basin and hard-rock sites in northern Italy. Geophysical Journal International 206(1), 395–407.Google Scholar
Bromirski, PD and 7 others (2015) Ross ice shelf vibrations. Geophysical Research Letters 42, 7589–7597. doi:10.1002/2015GL065284Google Scholar
Campillo, M and Paul, A (2003) Long range correlations in the diffuse seismic coda. Science 24(299), 547–549.Google Scholar
Chaput, J and 9 others (2018) Near-surface environmentally changes in the Ross Ice Shelf observed with ambient seismic noise. Geophysical Research Letters 45, 11187–11196. doi:10.1029/2018GL079665Google Scholar
Chaput, J, Aster, RC, Karplus, M and Nakata, N (2022b) Ambient high frequency seismic wavefields in Antarctic firn. Journal of Glaciology 68, 1–14. doi:10.1017/jog.2021.135Google Scholar
Chaput, J and 8 others (2022a) Near-surface seismic anisotropy in Antarctic glacial snow and ice revealed by high-frequency ambient noise. Journal of Glaciology 1, 1–17. doi:10.1017/jog.2022.98Google Scholar
Denolle, M, Dunham, E and Beroza, G (2012) Solving the surface-wave eigenproblem with Chebyshev spectral collocation. Bulletin of the Seismological Society of America 102(3), 1214–1223.Google Scholar
Diez, A and 8 others (2016) Ice shelf structure derived from dispersion curve analysis of ambient seismic noise, Ross Ice Shelf, Antarctica. Geophysical Journal International 205, 785–795. doi:10.1093/gji/ggw036Google Scholar
Gilbert, A, Vincent, C, Six, D, Wagnon, P, Piard, L and Ginot, P (2014) Modeling near-surface firn temperature in a cold accumulation zone (Col du Dôme, French Alps): from a physical to a semi-parameterized approach. The Cryosphere 8, 689–703. doi:10.5194/tc-8-689-2014Google Scholar
Haney, M and Tsai, V (2015) Nonperturbational surface-wave inversion: a Dix-type relation for surface waves. Geophysics 80(6), 06121. doi:10.1190/geo2014-0612.1Google Scholar
Horlings, A, Christianson, K, Holschuh, N, Stevens, CM and Waddington, ED (2021) Effect of horizontal divergence on estimates of firn-air content. Journal of Glaciology 67(262), 287–296. doi:10.1017/jog.2020.105Google Scholar
Hosokawa, T and Otani, A (1998) Acoustic anisotropy in bovine cancellous bone. The Journal of the Acoustical Society of America 103, 2718–2722. doi:10.1121/1.422790Google Scholar
Komatitsch, D and Tromp, J (2002) Spectral-element simulations of global seismic wave propagation-I. Validation. Geophysical Journal International 149(2), 390–412. doi:10.1046/j.1365-246X.2002.01653.xGoogle Scholar
Kuipers Munneke, P, Ligtenberg, SRM, van den Broeke, MR and Vaughan, DG (2017) Firn air depletion as a precursor of Antarctic ice-shelf collapse. Journal of Glaciology 60(220), 205–214. doi:10.3189/2014JoG13J183Google Scholar
Lazzara, MA, Weidner, GA, Keller, LM, Thom, JE and Cassano, JJ (2012) Antarctic automatic weather station program: 30 years of polar observation. Bulletin of the American Meteorological Society 93(10), 1519–1537. doi:10.1175/bams-d-11-00015.1Google Scholar
Ledoux, C, Hulbe, C, Forbes, M, Scambos, T and Alley, K (2017) Structural provinces of the Ross Ice Shelf, Antarctica. Annals of Glaciology 58(75), 88–98.Google Scholar
Lee, K, Hughes, E, Humphrey, V, Leighton, L and Choi, M (2007) Empirical angle dependent Biot and MBA models for acoustic anisotropy in cancellous bone. Physics in Medicine and Biology 52(1), 59–73. doi:10.1088/0031-9155/52/1/005Google Scholar
Leppäranta, M, Järvinen, O and Mattila, OP (2012) Structure and life cycle of supraglacial lakes in Dronning Maud Land. Antarctic Science 25(3), 457–467. doi:10.1017/s0954102012001009Google Scholar
Ligtenberg, SRM, Helsen, MM and van den Broeke, MR (2011) An improved semi-empirical model for the densification of Antarctic firn. The Cryosphere 5(4), 809–819. doi:10.5194/tc-5-809-2011Google Scholar
MacAyeal, DR (2018) Seismology gets under the skin of the Antarctic ice sheet. Geophysical Research Letters 45(20), 11173–11176. doi:10.1029/2018GL080366Google Scholar
Melon, M, Lafarge, D, Castagnede, B and Brown, N (1995) Measurement of tortuosity of anisotropic acoustic materials. Journal of Applied Physics 78(4929), 4929–4932. doi:10.1063/1.359781Google Scholar
Melon, M, Mariez, E, Ayrault, C and Sahraoui, S (1998) Acoustical and mechanical characterization of anisotropic open cell foams. The Journal of the Acoustical Society of America 104(2622), 2622–2627. doi:10.1121/1.423897Google Scholar
Mordret, A, Mikesell, D, Harig, C, Lipovsky, B and Prieto, H (2016) Monitoring southwest Greenland's ice sheet melt with ambient seismic noise. Science Advances 2, e1501538. doi:10.1126/sciadv.1501538Google Scholar
Nakamura, Y (1989) A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface. Quarterly Report from the Railway Technical Research Institute 30(1), 25–30.Google Scholar
Nicolas, JP and 13 others (2017) January 2016 extensive summer melt in West Antarctica favoured by strong El Niño. Nature Communications 8, 15799. doi:10.1038/ncomms15799Google Scholar
Oraschewski, F and Grinsted, A (2022) Modeling enhanced firn densification due to strain softening. The Cryosphere 16(7), 2683–2700. doi:10.5194/tc-16-2683-2022Google Scholar
Podolskiy, E and Walter, F (2016) Cryoseismology. Reviews of Geophysics 54(4), 708–758. doi:10.1002/2016RG000526Google Scholar
Qin, L, Qiu, H, Nakata, N, Deng, S, Levander, A and Ben-Zion, Y (2022) Variable daily autocorrelation functions of high-frequency seismic data on Mars. Seismological Research Letters 94, 746–758. doi:10.1785/0220220196Google Scholar
Reeh, N (2008) A nonsteady-state firn-densification model for the percolation zone of a glacier. Journal of Geophysical Research: Earth Surface 113, 1–13.Google Scholar
Reeh, N, Fisher, D, Koerner, R and Clausen, H (2005) An empirical firn-densification model comprising ice lenses. Annals of Glaciology 42, 101–106.Google Scholar
Rennermalm, AK and 7 others (2013) Evidence of meltwater retention within the Greenland ice sheet. The Cryosphere 7(5), 1433–1445. doi:10.5194/tc-7-1433-2013Google Scholar
Rignot, E, Casassa, G, Gogineni, P, Krabill, W, Rivera, A and Thomas, R (2004) Accelerated ice discharge from the Antarctic Peninsula following the collapse of Larsen B ice shelf. Geophysical Research Letters 31, 1–4.Google Scholar
Scambos, TA, Bohlander, JA, Shuman, CA and Skvarca, P (2004) Glacier acceleration and thinning after ice shelf collapse in the Larsen B embayment, Antarctica. Geophysical Research Letters 31(18), 1–4.Google Scholar
Schmelzbach, C and 9 others (2018) Advances in 6C seismology: applications of combined translational and rotational motion measurements in global and exploration seismology. Geophysics 83(3), 53–69. doi:10.1190/geo2017-0492.1Google Scholar
Sharma, V, Braud, L and Lehning, M (2019) Understanding snow bedform formation by adding sintering to a cellular automata model. The Cryosphere 13, 3239–3260. doi:10.5194/tc-13-3239-2019Google Scholar
Snieder, R (2004) Extracting the Green's function from the correlation of coda waves: a derivation based on stationary phase. Physical Review E 69(4), 046610.Google Scholar
Sommer, C, Wever, N, Fierz, C and Lehning, M (2018) Wind-packing of snow in Antarctica. The Cryosphere 36, 2923–2939. doi:10.5194/tc-2018-36Google Scholar
Steger, Cea (2017) Firn meltwater retention on the Greenland ice sheet: a model comparison. Frontiers in Earth Sciences 5(3), 1–16. doi:10.3389/feart.2017.00003Google Scholar
Stevens, M and 6 others (2020) The community firn model (cfm) v1.0. Geoscientific Model Development 13, 4355–4377. doi:10.5194/gmd-13-4355-2020Google Scholar
Takei, I and Maeno, N (2004) Mechanical vibration responses of snow samples near the melting temperature. Annals of Glaciology 38, 130–134.Google Scholar
Tanimoto, T and Rivera, L (2005) Prograde Rayleigh wave particle motion. Geophysical Journal International 162(2), 399–405.Google Scholar
Tita, V and Caliri Junior, M (2012) Numerical simulation of anisotropic polymeric foams. Latin American Journal of Solids and Structures 9(2), 1–12. doi:10.1590/S1679-78252012000200005Google Scholar
van den Broeke, M (2008) Depth and density of the Antarctic firn layer. Arctic, Antarctic, and Alpine Research 40(2), 432–438. doi:10.1657/1523-0430(07-021)[BROEKE]2.0.CO;2Google Scholar
Vandecrux, B and 22 others (2020) The firn meltwater Retention Model Intercomparison Project (RetMIP): evaluation of nine firn models at four weather station sites on the Greenland ice sheet. The Cryosphere 14(11), 3785–3810. doi:10.5194/tc-14-3785-2020Google Scholar
Wapenaar, K, Draganov, D, Snieder, R, Campman, X and Verdel, A (2010) Tutorial on seismic interferometry: Part 1 – Basic principles and applications. Geophysics 75(5), 1942–2156.Google Scholar
Wapenaar, K and Fokkema, J (2006) Green's function representations for seismic interferometry. Geophysics 71(4), S133–S146.Google Scholar
Zhan, Z (2019) Seismic noise interferometry reveals transverse drainage configuration beneath the surging Bering Glacier. Geophysical Research Letters 46(9), 4747–4756. doi:10.1029/2019gl082411Google Scholar
Figure 0

Figure 1. (A) Example of firn resonances from nearly two years of North component ambient seismic data recorded at station DR09 on Roosevelt Island from Ross Ice Shelf broadband array (red circle, inset map) displayed as a time/frequency plot (spectrogram) (Chaput and others, 2018). Red boxes indicate seasonal open sea ice conditions, and corresponding spectral effects are not observed at grounded sites. Stable shelf plate modes are visible as high amplitude temporally stable vertical bands below 5 Hz (notably during open sea ice conditions, red boxes), and temporally variable firn modes above roughly 5 Hz (observable year round). (B) Description of environmental effects that dictate the behavior of firn resonances. (1) Wind coupling with semi-periodic surface snowforms and the low-velocity/density firn structure excites unique firn mode patterns. (2) Firn is sensitive to anomalous near-zero surface temperatures, and the frequency range over which resonances are altered depends on the depth penetration of the temperature anomaly. (3) Firn accommodates strain associated with flowing ice masses in a ductile fashion at shallow depths where porosity is high, and in a brittle fashion where pores have largely closed.

Figure 1

Figure 2. Examples of firn resonances show in ‘peak tracked’ form, where only peak maximums for both horizontal components of the spectrogram at each time bin are shown. (A) Firn resonance response to alterations in surface snowforms. Strong storm activity (arrows in left panel, matched with periods of high winds shown by the black trace), can deposit new snowforms that are then slowly eroded during periods of quiescence (right panel), resulting in slow spectral decay (time scale of months, black dotted boxes). (B) Peak tracked firn resonances at 5 stations on RIS during a shelf-wide near-zero temperature event in 2016 (bottom left panel, Nicolas and others (2017)). Firn undergoes up to a 40% reduction in elastic moduli (top right panel) as bonds between snow grains weaken, resulting in a downward drift in frequency for higher peaks and a reduction in amplitude (bottom right panel). (C) Peak tracked spectrogram at RIS station RS17 for 10 days during 2016, showing the obvious offset in frequency (also frequency dependent) between North and East components. Right panel: Shallow firn deforms plastically under extensional strain typical of RIS and features a strain-elongated pore space (black arrows aligned with ice flow), while deeper firn to solid ice responds in a brittle manner, often resulting in flow-perpendicular crevassing (red arrows) unless dominant crevassing is advected and rotated from past strain regimes.

Figure 2

Figure 3. Multi-faceted seismic and distributed acoustic sensing (DAS) experiment coupled with drone-based photogrammetry aiming to make key firn seismic and environmental measurements to advance understanding of the firn medium and environmental forcing effects that govern resonance peak and other seismological observations, leveraging a long-running UW-Madison (UWM) autonomous weather or similar station (e.g. Lazzara and others, 2012).