Fan Chi Lin

Seismic evidence for widespread western-US deep-crustal deformation caused by extension

Laboratory experiments have established that many of the materials comprising the Earth are strongly anisotropic in terms of seismic-wave speeds. Observations of azimuthal and radial anisotropy in the upper mantle are attributed to the lattice-preferred orientation of olivine caused by the shear strains associated with deformation, and provide some of the most direct evidence for deformation and flow within the Earth’s interior. Although observations of crustal radial anisotropy would improve our understanding of crustal deformation and flow patterns resulting from tectonic processes, large-scale observations have been limited to regions of particularly thick crust. Here we show that observations from ambient noise tomography in the western United States reveal strong deep (middle to lower)-crustal radial anisotropy that is confined mainly to the geological provinces that have undergone significant extension during the Cenozoic Era (since ∼65 Myr ago). The coincidence of crustal radial anisotropy with the extensional provinces of the western United States suggests that the radial anisotropy results from the lattice-preferred orientation of anisotropic crustal minerals caused by extensional deformation. These observations also provide support for the hypothesis that the deep crust within these regions has undergone widespread and relatively uniform strain in response to crustal thinning and extension.

P and S wave tomography of the mantle beneath the United States

Mantle seismic structure beneath the United States spanning from the active western plate margin to the passive eastern margin was imaged with teleseismic P and S wave traveltime tomography including USArray data up to May 2014. To mitigate artifacts from crustal structure 5–40 s, Rayleigh wave phase velocities were used to create a 3-D starting model. Major features of the final P and S models include two distinct low-velocity anomalies at depths of ~60–300 km beneath the central and northern Appalachians and passive margin. The central Appalachian low-velocity anomaly coincides with Eocene basaltic magmatism, and the northern anomaly is located along the Cretaceous track of the Great Meteor hot spot. At depths of ~300–700 km beneath the central and eastern U.S. large high-velocity anomalies are inferred to be remnants of the Farallon slab that subducted prior to ~40 Ma during the Laramide orogeny.

The Yellowstone magmatic system from the mantle plume to the upper crust

Yellowstones missing magmatic link Yellowstone is an extensively studied “supervolcano” that has a large supply of heat coming from a pool of magma near the surface and the mantle below. A link between these two features has long been suspected. Huang et al. imaged the lower crust using seismic tomography (see the Perspective by Shapiro and Koulakov). Their findings provide an estimate of the total amount of molten rock beneath Yellowstone and help to explain the large amount of volcanic gases escaping from the region. Science, this issue p. 773; see also p. 758 The Yellowstone supervolcano has a large magma body between the mantle hot spot and the upper crustal magmatic reservoir. [Also see Perspective by Shapiro and Koulakov] The Yellowstone supervolcano is one of the largest active continental silicic volcanic fields in the world. An understanding of its properties is key to enhancing our knowledge of volcanic mechanisms and corresponding risk. Using a joint local and teleseismic earthquake P-wave seismic inversion, we revealed a basaltic lower-crustal magma body that provides a magmatic link between the Yellowstone mantle plume and the previously imaged upper-crustal magma reservoir. This lower-crustal magma body has a volume of 46,000 cubic kilometers, ~4.5 times that of the upper-crustal magma reservoir, and contains a melt fraction of ~2%. These estimates are critical to understanding the evolution of bimodal basaltic-rhyolitic volcanism, explaining the magnitude of CO2 discharge, and constraining dynamic models of the magmatic system for volcanic hazard assessment.

Single-channel currents produced by the serotonin transporter and analysis of a mutation affecting ion permeation.

Single-channel activities were observed in outside-out patches excised from oocytes expressing a mammalian 5-hydroxytryptamine (5-HT) transporter. Channel conductance was larger for a mutant in which asparagine177 of the third putative transmembrane domain was replaced by glycine, suggesting that this residue lies within or near the permeation pathway. The N177G mutant enables quantitative single-channel measurements; it displays two conducting states. One state, with conductance of approximately 6 pS, is induced by 5-HT and is permeable to Na+. The other state (conductance of approximately 13 pS) is associated with substrate-independent leakage current and is permeable to both Na+ and Li+. Cl- is not a major current carrier. Channel lifetimes under all conditions measured are approximately 2.5 ms. The single-channel phenomena account for previously observed macroscopic electrophysiological phenomena, including 5-HT-induced transport-associated currents and substrate-independent leakage currents. The channel openings occur several orders of magnitude less frequently than would be expected if one such opening occurred for each transport cycle and therefore do not represent an obligatory step in transport. Nevertheless, single-channel events produced by neurotransmitter transporters indicate the functional and structural similarities between transporters and ion channels and provide a new tool, at single-molecule resolution, for detailed structure-function studies of transporters.

Joint inversion of Rayleigh wave phase velocity and ellipticity using USArray: Constraining velocity and density structure in the upper crust

Rayleigh wave ellipticity, or H/V ratio, observed on the surface is particularly sensitive to shallow earth structure. In this study, we jointly invert measurements of Rayleigh wave H/V ratio and phase velocity between 24–100 and 8–100 sec period, respectively, for crust and upper mantle structure beneath more than 1000 USArray stations covering the western United States. Upper crustal structure, in particular, is better constrained by the joint inversion compared to inversions based on phase velocities alone. In addition to imaging Vs structure, we show that the joint inversion can be used to constrain Vp/Vs and density in the upper crust. New images of uppermost crustal structure (<3 km depth) are in excellent agreement with known surface features, with pronounced low Vs, low density, and high Vp/Vs anomalies imaged in the locations of several major sedimentary basins including the Williston, Powder River, Green River, Denver, and San Juan basins. These results demonstrate not only the consistency of broadband H/V ratios and phase velocity measurements, but also that their complementary sensitivities have the potential to resolve density and Vp/Vs variations.

Distinct crustal isostasy trends east and west of the Rocky Mountain Front

Seismic structure beneath the contiguous U.S. was imaged with multimode receiver function stacking and inversion of Rayleigh wave dispersion and ellipticity measurements. Crust thickness and elevation are weakly correlated across the contiguous U.S., but the correlation is ~3–4 times greater for separate areas east and west of the Rocky Mountain Front (RMF). Greater lower crustal shear velocities east of the RMF, particularly in low-elevation areas with thick crust, are consistent with deep crustal density as the primary cause of the contrasting crust thickness versus elevation trends. Separate eastern and western trends are best fit by Airy isostasy models that assume lower crust to uppermost mantle density increases of 0.18 g/cm3 and 0.40 g/cm3, respectively. The former value is near the minimum that is plausible for felsic lower crust. Location of the transition at the RMF suggests that Laramide to post-Laramide processes reduced western U.S. lower crustal density.

Ambient seismic noise tomography of Canada and adjacent regions: Part I. Crustal structures

This paper presents the first continental-scale study of the crust and upper mantle shear velocity (V_s) structure of Canada and adjacent regions using ambient noise tomography. Continuous waveform data recorded between 2003 and 2009 with 788 broadband seismograph stations in Canada and adjacent regions were used in the analysis. The higher primary frequency band of the ambient noise provides better resolution of crustal structures than previous tomographic models based on earthquake waveforms. Prominent low velocity anomalies are observed at shallow depths (<20 km) beneath the Gulf of St. Lawrence in east Canada, the sedimentary basins of west Canada, and the Cordillera. In contrast, the Canadian Shield exhibits high crustal velocities. We characterize the crust-mantle transition in terms of not only its depth and velocity but also its sharpness, defined by its thickness and the amount of velocity increase. Considerable variations in the physical properties of the crust-mantle transition are observed across Canada. Positive correlations between the crustal thickness, Moho velocity, and the thickness of the transition are evident throughout most of the craton except near Hudson Bay where the uppermost mantle V_s is relatively low. Prominent vertical V_s gradients are observed in the midcrust beneath the Cordillera and beneath most of the Canadian Shield. The midcrust velocity contrast beneath the Cordillera may correspond to a detachment zone associated with high temperatures immediately beneath, whereas the large midcrust velocity gradient beneath the Canadian Shield probably represents an ancient rheological boundary between the upper and lower crust.

The local amplification of surface waves: A new observable to constrain elastic velocities, density, and anelastic attenuation

The deployment of USArray across the continental U.S. has prompted developments within surface wave tomography to exploit this unprecedented data set. Here, we present a method to measure a new surface wave observable: broadband surface wave amplification that provides new and unique constraints on elastic velocities and density within the crust and upper mantle. The method, similar to its phase velocity counterpart referred to as Helmholtz tomography, initiates by constructing phase travel time and amplitude maps across the array for each period and earthquake. Spatial differential operators are then applied to evaluate the amplitude variation, as well as the effect of focusing/defocusing. Based on the 2-D damped wave equation, the amplitude variation corrected for focusing/defocusing is linked directly to both local amplification and intrinsic attenuation, which are separated by examining waves propagating in opposite directions. We apply the method to teleseismic Rayleigh waves observed across USArray between periods of 24 and 100 s and show that the observed amplification maps are strongly correlated with known geological features. Small-scale attenuation measurements are contaminated by wavefield complexities, but larger-scale anelastic attenuation is estimated reliably. The observed amplification maps compare well with predictions based on recent 3-D shear velocity models of the western U.S. that were produced from ambient noise and earthquake data. Notably, predictions based on models with different prescribed density structures demonstrate the potential for using estimates of local amplification to constrain not only 3-D velocity structure but also density.

Is ambient noise tomography across ocean basins possible

[1] Based on year-long cross-correlations of broad-band seismic records obtained at sixty-six stations within or adjacent to the Pacific Basin, we show that broad-band ambient noise is observed to propagate coherently between island stations and between island and continent stations. For many station pairs, high signal-to-noise ratio (SNR) fundamental mode Rayleigh wave Green functions are observed, which establishes the physical basis for ambient noise tomography across the Pacific. Similar trends for continental and oceanic stations are observed in the relationship between the ambient noise level at a station and the ‘‘noise coherence distance’’ – the longest distance at which a high SNR cross-correlation signal is observed for a station. Because locally generated noise obscures long distance coherent noise, situating stations at quiet locations on islands is necessary for the success of ambient noise tomography. Local noise poses a particular challenge at atoll sites and, on the basis of analysis of data from station H2O, at ocean bottom sites at periods above 25 sec. Citation: Lin, F.-C., M. H. Ritzwoller, and N. M. Shapiro (2006), Is ambient noise tomography across ocean basins possible?, Geophys. Res. Lett., 33, L14304, doi:10.1029/ 2006GL026610.

A one-dimensional seismic model for Uturuncu volcano, Bolivia, and its impact on full moment tensor inversions

Using receiver functions, Rayleigh wave phase velocity dispersion determined from ambient noise and teleseismic earthquakes, and Rayleigh wave horizontal to vertical ground motion amplitude ratios from earthquakes observed across the PLUTONS seismic array, we construct a one-dimensional (1‑D) S-wave velocity (Vs) seismic model with uncertainties for Uturuncu volcano, Bolivia, located in the central Andes and overlying the eastward-subducting Nazca plate. We find a fast upper crustal lid placed upon a low-velocity zone (LVZ) in the mid-crust. By incorporating all three types of measurements with complimentary sensitivity, we also explore the average density and Vp/Vs (ratio of P-wave to S-wave velocity) structures beneath the young silicic volcanic field. We observe slightly higher Vp/Vs and a decrease in density near the LVZ, which implies a dacitic source of the partially molten magma body. We exploit the impact of the 1-D model on full moment tensor inversion for the two largest local earthquakes recorded (both magnitude ∼3), demonstrating that the 1-D model influences the waveform fits and the estimated source type for the full moment tensor. Our 1-D model can serve as a robust starting point for future efforts to determine a three-dimensional velocity model for Uturuncu volcano.