skip to main content


Title: Hydrologic characterization of an alpine valley infill through integration of ERT, active seismic and active/passive surface wave interferometry (Unaweep Canyon, US)
Unaweep Canyon (Western Colorado, US) is an enigmatic alpine landform and hypothesized to represent a partially exhumed paleo valley which was glacially over-deepened in the late Paleozoic. Processing and interpretation of recently acquired 2D seismic reflection and refraction data support the concept of glacial over-deepening and indicate maximum bedrock depths of about 550 meters. Additionally, pronounced reflectors are observed within the sedimentary infill. The seismic data have also been subjected to surface wave analysis revealing a significant increase of the Vp/Vs ratio below a shallow (50 - 150 m depth) intra-sedimentary reflector. A large Vp/Vs ratio can be caused by both saturation and poor consolidation of dry low-porosity materials (e.g. dry sands).To investigate the potential occurrence of an aquifer associated with this interface, a high-density/long-offset electrical resistivity survey was conducted in fall 2019 along the seismic line. The maximum offset is 915 m at an electrode spacing of 5 meters, aiming at reaching depths of investigations between 150 and 200 meters. Inversion of the ERT data was initially conducted by means of smoothness-constrained algorithms. The imaging results revealed consistent structures with those resolved through seismic methods, at least within the required depth of investigation between 150 - 200 m. Furthermore, improvements in the resolution of the ERT imaging results was investigated after the inclusion of seismic interfaces as structural constraints in the inversion of the data. The comparison of the two approaches permitted to improve the interpretation of the ERT imaging results, which indicate low resistivities in the zone of high Vp/Vs ratios and thus strengthen the aquifer hypothesis. We present an integrated interpretation based on seismic structure, resistivity distribution, Vp and Vs velocities, and a distant well core. In a larger context, the results provide new insights on the subsurface hydrology in this arid part of the continental US as well as on the significance of multi-valued datasets for the interpretation and characterization of aquifers.  more » « less
Award ID(s):
1849623
NSF-PAR ID:
10179156
Author(s) / Creator(s):
Date Published:
Journal Name:
EGU General Assembly
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present two new seismic velocity models for Alaska from joint inversions of body-wave and ambient-noise-derived surface-wave data, using two different methods. Our work takes advantage of data from many recent temporary seismic networks, including the Incorporated Research Institutions for Seismology Alaska Transportable Array, Southern Alaska Lithosphere and Mantle Observation Network, and onshore stations of the Alaska Amphibious Community Seismic Experiment. The first model primarily covers south-central Alaska and uses body-wave arrival times with Rayleigh-wave group-velocity maps accounting for their period-dependent lateral sensitivity. The second model results from direct inversion of body-wave arrival times and surface-wave phase travel times, and covers the entire state of Alaska. The two models provide 3D compressional- (VP) and shear-wave velocity (VS) information at depths ∼0–100  km. There are many similarities as well as differences between the two models. The first model provides a clear image of the high-velocity subducting plate and the low-velocity mantle wedge, in terms of the seismic velocities and the VP/VS ratio. The statewide model provides clearer images of many features such as sedimentary basins, a high-velocity anomaly in the mantle wedge under the Denali volcanic gap, low VP in the lower crust under Brooks Range, and low velocities at the eastern edge of Yakutat terrane under the Wrangell volcanic field. From simultaneously relocated earthquakes, we also find that the depth to the subducting Pacific plate beneath southern Alaska appears to be deeper than previous models. 
    more » « less
  2. Abstract

    Near‐surface seismic velocity structure plays a critical role in ground motion amplification during large earthquakes. In particular, the local Vp/Vs ratio strongly influences the amplitude of Rayleigh waves. Previous studies have separately imaged 3D seismic velocity and Vp/Vs ratio at seismogenic depth, but lack regional coverage and/or fail to constrain the shallowest structure. Here, we combine three datasets with complementary sensitivity in a Bayesian joint inversion for shallow crustal shear velocity and near‐surface Vp/Vs ratio across Southern California. Receiver functions–including with an apparent delayed initial peak in sedimentary basins, and long considered a nuisance in receiver function imaging studies–highly correlate with short‐period Rayleigh wave ellipticity measurements and require the inclusion of a Vp/Vs parameter. The updated model includes near‐surface low shear velocity more in line with geotechnical layer estimates, and generally lower than expected Vp/Vs outside the basins suggesting widespread shallow fracturing and/or groundwater undersaturation.

     
    more » « less
  3. Abstract

    We present electrical resistivity models of the crust and upper mantle from two‐dimensional (2‐D) inversion of magnetotelluric (MT) data collected in the Rio Grande rift, New Mexico, USA. Previous geophysical studies of the lithosphere beneath the rift identified a low‐velocity zone several hundred kilometers wide, suggesting that the upper mantle is characterized by a very broad zone of modified lithosphere. In contrast, the surface expression of the rift (e.g., high‐angle normal faults and synrift sedimentary units) is confined to a narrow region a few tens of kilometers wide about the rift axis. MT data are uniquely suited to probing the depths of the lithosphere that fill the gap between surface geology and body wave seismic tomography, namely the middle to lower crust and uppermost mantle. We model the electrical resistivity structure of the lithosphere along two east‐west trending profiles straddling the rift axis at the latitudes of 36.2 and 32.0°N. We present results from both isotropic and anisotropic 2‐D inversions of MT data along these profiles, with a strong preference for the latter in our interpretation. A key feature of the anisotropic resistivity modeling is a broad (~200‐km wide) zone of enhanced conductivity (<20 Ωm) in the middle to lower crust imaged beneath both profiles. We attribute this lower crustal conductor to the accumulation of free saline fluids and partial melt, a direct result of magmatic activity along the rift. High‐conductivity anomalies in the midcrust and upper mantle are interpreted as fault zone alteration and partial melt, respectively.

     
    more » « less
  4. null (Ed.)
    SUMMARY Interfaces are important part of Earth’s layering structure. Here, we developed a new model parametrization and iterative linearized inversion method that determines 1-D crustal velocity structure using surface wave dispersion, teleseismic P-wave receiver functions and Ps and PmP traveltimes. Unlike previous joint inversion methods, the new model parametrization includes interface depths and layer Vp/Vs ratios so that smoothness constraint can be conveniently applied to velocities of individual layers without affecting the velocity discontinuity across the interfaces. It also allows adding interface-related observation such as traveltimes of Ps and PmP in the joint inversion to eliminate the trade-off between interface depth and Vp/Vs ratio and therefore to reduce the uncertainties of results. Numerical tests show that the method is computationally efficient and the inversion results are robust and independent of the initial model. Application of the method to a dense linear array across the Wabash Valley Seismic Zone (WVSZ) produced a high-resolution crustal image in this seismically active region. The results show a 51–55-km-thick crust with a mid-crustal interface at 14–17 km. The crustal Vp/Vs ratio varies from 1.69 to 1.90. There are three pillow-like, ∼100 km apart high-velocity bodies sitting at the base of the crust and directly above each of them are a low-velocity anomaly in the middle crust and a high-velocity anomaly in the upper crust. They are interpreted to be produced by mantle magmatic intrusions and remelting during rifting events in the end of the Precambrian. The current diffuse seismicity in the WVSZ might be rooted in this ancient distributed rifting structure. 
    more » « less
  5. Abstract

    Although seismic velocity and electrical conductivity are both sensitive to temperature, thermal lithosphere properties are derived almost exclusively from seismic data because conductivity is often too strongly affected by minor highly conductive phases to be a reliable indicator of temperature. However, in certain circumstances, electrical observations can provide strong constraints on mantle temperatures. In the southeastern United States (SEUS), magnetotelluric (MT) data require high resistivity values (>300 Ωm) to at least 200‐km depth. As dry mantle mineral conduction laws provide an upper bound on temperature for an observed resistivity value, the only interpretation is that lithospheric temperatures (<1330 °C) persist to 200 km. However, seismic tomography shows that velocities in this region are generally slightly slow with respect to references models; this observation has led to a view of relatively thin (<150 km), eroded thermal lithosphere beneath the SEUS. We show that MT and seismic (tomography, attenuation, receiver function) results are consistent with thick (~200 km), coherent thermal lithosphere in this region. Reduced seismic velocities (relative to reference models) can be explained by considering the effect of finite grain size (anelasticity). Calculated velocity as a function of temperature is overall slower when including anelastic effects, even at reasonable grain sizes of 1 mm to 1 cm; this permits mantle temperatures that are colder than would typically be inferred. We argue for a geodynamic scenario in which the present thermal lithosphere in the SEUS formed in association with the Central Atlantic Magmatic Province and has subsequently survived intact for ~200 Ma.

     
    more » « less