Effect of snowfall on changes in relative seismic velocity measured by ambient noise correlation

In mountainous, cold temperate and polar sites, the presence of snow cover can affect relative seismic velocity changes (dV/V) derived from ambient noise correlation, but this relation is relatively poorly documented and ambiguous. In this study, we analyzed raw seismic recordings from a snowy flat field site located above Davos (Switzerland), during one entire winter season (from December 2018 to June 2019). We identified three snowfall events with a substantial response of dV/V measurements (drops of several percent between 15 and 25 Hz), suggesting a detectable change in elastic properties of the medium due to the additional fresh snow. To better interpret the measurements, we used a physical model to compute frequency-dependent changes in the Rayleigh wave velocity computed before and after the events. Elastic parameters of the ground subsurface were obtained from a seismic refraction survey, whereas snow cover properties were obtained from the snow cover model SNOWPACK. The decrease in dV/V due to a snowfall was well reproduced, with the same order of magnitude as observed values, confirming the importance of the effect of fresh and dry snow on seismic measurements. We also observed a decrease in dV/V with snowmelt periods, but we were not able to reproduce those changes with our model. Overall, our results highlight the effect of the snow cover on seismic measurements, but more work is needed to accurately model this response, in particular for the presence of liquid water in the snowpack.

Effect of snowfall on changes in relative seismic velocity measured by ambient noise correlation

In mountainous, cold temperate and polar sites, the presence of a snow cover can affect relative seismic velocity changes (dV/V) derived from ambient noise correlation, but this relation is relatively poorly documented and ambiguous. In this study, we analyzed raw seismic recordings from a snowy flat field site located above Davos (Switzerland), during one entire winter season (from December 2018 to June 2019). We identified three snowfall events with a substantial response of dV/V measurements (drops of several percent between 15 and 25 Hz), suggesting a detectable change in elastic properties of the medium due to the additional fresh snow. To better interpret the measurements, we used a physical model to compute frequency dependent changes in the Rayleigh wave velocity computed before and after the events. Elastic parameters of the ground subsurface were obtained from a seismic refraction survey, whereas snow cover properties were obtained from the snow cover model SNOWPACK. The decrease in dV/V due to a snowfall were well reproduced, with the same order of magnitude as observed values, confirming the importance of the effect of fresh and dry snow on seismic measurements. We also observed a decrease in dV/V with snowmelt periods, but we were not able to reproduce those changes with our model. Overall, our results highlight the effect of the snowcover on seismic measurements, but more work is needed to accurately model this response, in particular for the presence of liquid water in the snowpack.

Porosity dependence of elastic moduli of snow and firn

Measurements of elastic wave velocities enable non-destructive estimation of the mechanical properties, elastic moduli and density of snow and firn. The variation of elastic moduli with porosity in dry snow and firn is modeled using a differential effective medium scheme modified to account for the critical porosity above which the bulk and shear moduli of the ice frame vanish. A comparison of predicted and measured elastic moduli indicates that the shear modulus of ice in snow is lower than that computed from single crystal elastic stiffnesses of ice. This may indicate that the bonds between snow particles are more deformable under shear than under compression. A partial alignment of ice crystals also may contribute. Good agreement between elastic stiffnesses of the ice frame obtained from elastic wave velocity measurements and the predictions of the theory is observed. The approach is simple and compact, and does not require the use of empirical fits to the data. Owing to its simplicity, this model may prove useful in a variety of potential applications such as construction on snow, interpretation of seismic measurements to monitor and locate avalanches and estimation of density within compacting snow deposited on glaciers and ice sheets.

Effect of snowfalls on relative seismic velocity changes measured by ambient noise correlation

<p>Seismic noise correlation is a broadly used method to monitor the subsurface, in order to detect physical processes into the surveyed medium such changes in rigidity, fluid injection or cracking <sup>(1)</sup>. The influence of several environmental variables on measured seismic observables were studied, such as temperature, groundwater level fluctuations, and freeze-thawing cycles <sup>(2)</sup>. In mountainous, cold temperate and polar sites, the presence of a snowcover can also affect relative seismic velocity changes (dV/V), but this relation is relatively poorly documented and ambiguous <sup>(3)(4)</sup>. In this study, we analyzed raw seismic recordings from a snowy flat field site located above Davos (Switzerland), during one entire winter season (from December 2018 to June 2019). Our goal was to better understand the effect of snowfall and snowmelt events on dV/V measurements through both seismic and meteorological instrumentation.</p><p>We identified three snowfall events with a substantial response of dV/V measurements (drops of several percent between 15 and 25 Hz), suggesting a detectable change in elastic properties of the medium due to the additional fresh snow.</p><p>To better interpret the measurements, we used a physical model to compute frequency dependent changes in the Rayleigh wave velocity computed before and after the events. Elastic parameters of the ground subsurface were obtained from a seismic refraction survey, whereas snow cover properties were obtained from the snow cover model SNOWPACK. The decrease in dV/V due to a snowfall were well reproduced, with the same order of magnitude than observed values, confirming the importance of the effect of fresh and dry snow on seismic measurements.</p><p>We also observed a decrease in dV/V with snowmelt periods, but we were not able to reproduce those changes with our model. Overall, our results highlight the effect of the snowcover on seismic measurements, but more work is needed to accurately model this response, in particular for the presence of liquid water in the snowcover.</p><p> </p><p><strong>References</strong></p><ul><li>(1) Larose, E., Carrière, S., Voisin, C., Bottelin, P., Baillet, L., Guéguen, P., Walter, F., et al. (2015) Environmental seismology: What can we learn on earth surface processes with ambient noise? Journal of Applied Geophysics, <strong>116</strong>, 62–74. doi:10.1016/j.jappgeo.2015.02.001</li> <li>(2) Le Breton, M., Larose, É., Baillet, L., Bontemps, N. & Guillemot, A. (2020) Landslide Monitoring Using Seismic Ambient Noise Interferometry: Challenges and Applications. Earth-Science Reviews</li> <li>(3) Hotovec‐Ellis, A.J., Gomberg, J., Vidale, J.E. & Creager, K.C. (2014) A continuous record of intereruption velocity change at Mount St. Helens from coda wave interferometry. Journal of Geophysical Research: Solid Earth, <strong>119</strong>, 2199–2214. doi:10.1002/2013JB010742</li> <li>(4) Wang, Q.-Y., Brenguier, F., Campillo, M., Lecointre, A., Takeda, T. & Aoki, Y. (2017) Seasonal Crustal Seismic Velocity Changes Throughout Japan. Journal of Geophysical Research: Solid Earth, <strong>122</strong>, 7987–8002. doi:10.1002/2017JB014307</li> </ul>

Saturn's Diffuse Core from Ring Seismology

<p>Gravity field measurements only weakly constrain the deep interiors of Jupiter and Saturn, stymieing efforts to measure the mass and compactness of these planets' cores, crucial properties for understanding their formation pathways and evolution. However, studies of Saturn's rings by Cassini have revealed waves driven by pulsation modes within Saturn, offering independent seismic probes of Saturn's interior. The observations reveal gravity mode (g mode) pulsations that indicate that a part of Saturn's interior is stably stratified by composition gradients, and the g mode frequencies directly probe the buoyancy frequency within the planet.</p><p>We compare structure models with gravity and new seismic measurements from Cassini to show that the data can only be explained by a diffuse, stably stratified core-envelope transition region in Saturn extending to approximately 60% of the planet's radius. This predominantly stable interior imposes significant constraints on Saturn's intrinsic magnetic field generation. The gradual distribution of heavy elements required by the seismology constrains mixing processes at work in Saturn, and it may reflect the planet's primordial structure and accretion history.</p>

Collisional Orogeny in the Scandinavian Caledonides (COSC): Some preliminary results from drilling of the 2.276 km deep COSC-2 borehole, central Sweden

<p>COSC investigations and drilling activities are focused in the Åre-Mörsil area (Sweden) of central Scandinavia. COSC-2 was drilled with nearly 100% core recovery in 2020 to 2.276 km depth with drilling ongoing from mid-April to early August. Drilling targets for COSC-2 included (1) the highly conductive Alum shale, (2) the Caledonian décollement, the major detachment that separates the Caledonian allochthons from the autochthonous basement of the Fennoscandian Shield, and (3) the strong seismic reflectors in the Precambrian basement.</p><p>Combined seismic, magnetotelluric (MT) and magnetic data were used to site the COSC-2 borehole about 20 km east-southeast of COSC-1. Based on these data it was predicted that the uppermost, tectonic occurrence of Cambrian Alum shale would be penetrated at about 800 m, the main décollement in Alum shale at its stratigraphic level at about 1200 m and the uppermost high amplitude basement reflector at about 1600 m. Paleozoic turbidites and greywackes were expected to be drilled down to 800 m depth. Below this depth, Ordovician limestone and shale with imbricates of Alum shale were interpreted to be present. Directly below the main décollement, magnetite rich Precambrian basement was expected to be encountered with a composition similar to that of magnetic granitic rocks found east of the Caledonian Front. The actual depths of the main contacts turned out to agree very well with the predictions based on the geophysical data. However, the geology below the uppermost occurrence of Alum shale is quite different from the expected model. Alum shale was only clearly encountered as a highly deformed, about 30 m thick unit, starting at about 790 m. Between about 820 and 1200 m, preliminary interpretations are that the rocks mainly consist of Neo-Proterozoic to Early Cambrian tuffs. Further below, Precambrian porphyries are present. The high amplitude reflections within the Precambrian sequence appear to be generated by dolerite sheets with the uppermost top penetrated at about 1600 m. Several deformed sheets of dolerite may be present down to about 1930 m. Below this depth the rocks are again porphyries.</p><p>A preliminary conclusion concerning the tectonic model is that the main décollement is at about 800 m and not at 1200 m. Also the thickness of the lowermost Cambrian/uppermost Neoproterozoic sediments on top of the basement is much greater than expected (hundreds of meters instead of tens of meters) and likely to have been thickened tectonically. Detailed studies are required to assess the actual importance of the “main décollement” and the degree, type and age of deformation in its footwall. We can also conclude that the Precambrian basement is very similar to the Dala porphyries succession that are typically present farther south.</p><p>An extensive set of downhole logging data was acquired directly after drilling. Borehole seismic measurements in 2021 will help to define and correlate seismic boundaries with lithology and structures in the core. Unfortunately, work for describing the geology of the drill core in detail is still on hold due to Covid-19.</p>

Seismic quantification of river flow depth - from the flume to the field

<p>Seismic measurements are used to study various processes that shape the Alpine landscape, including rock falls, debris flows, bedload transport and turbulent water flow. Here, we focus on the seismic quantification of turbulent flow conditions which is particularly useful for the remote monitoring of channels that are inaccessible (e.g. subglacial conduits) and/or highly dynamic (e.g. actively braiding river reaches). We test a physically-based model (Gimbert et al., 2014) to quantify force spectra generated by turbulent flow in flume experiments performed by Lamb et al. (2017) and subsequently apply the model to estimate river flow depth from continuous seismic measurements in the field.</p><p>In the flume, we assess near-bed flow velocity spectra and resulting drag and lift force spectra experienced by particles (D=0.075-0.20 m) on the cobble bed for a wide range of channel gradients (S=0.004-0.3) and submergence levels (h/D<sub>50</sub>=1-9.6). These measurements are used to test our model, and to quantify wake (interaction) effects and fluid-dynamic admittance on force spectral amplitude. Based on the conservation of turbulent energy in the Kolmogorov inertial subrange, we predict lift and drag force spectra to within ±5 dB rel. N<sup>2</sup>/Hz (frequency ~10-25 Hz) of the measured values.</p><p>We apply the calibrated model to bank-side geophone measurements from an Alpine stream (Séveraisse River, France). Using locally-derived seismic parameters, riverbed particle-size distribution and bed roughness, we can invert for water depth over a range of flow conditions, including flows with bedload transport (bedload transport dominates the seismic signal at higher frequencies). This allows us to monitor changes in flow depth during the course of a high-magnitude flood (October 2019). During the falling limb, the inferred flow depths progressively deviate from independently made water level measurements, indicating local riverbed aggradation of approximately 0.5 m, which is in agreement with post-flood observations. Through insights in near-bed turbulent flow conditions and their seismic signature, we can study flow-bedload transport interactions and the effects of extreme flow events on river morphodynamics.</p><p><strong>References</strong></p><p>Gimbert, F., Tsai, V. C. & Lamb, M. P. (2014). A physical model for seismic noise generation by turbulent flow in rivers. Journal of Geophysical Research: Earth Surface, 119(10), 2209-2238. http://dx.doi.org/10.1002/2014JF003201</p><p>Lamb, M. P., Brun, F. & Fuller, B. M. (2017). Direct measurements of lift and drag on shallowly submerged cobbles in steep streams: Implications for flow resistance and sediment transport. Water Resources Research, 53(9), 7607-7629. https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1002/2017WR020883</p>