Juan de Fuca subduction zone from a mixture of tomography and waveform modeling

Seismic tomography images of the upper mantle structures beneath the Pacific Northwestern United States display a maze of high-velocity anomalies, many of which produce distorted waveforms evident in the USArray observations indicative of the Juan de Fuca (JdF) slab. The inferred location of the slab agrees quite well with existing contour lines defining the slab's upper interface. Synthetic waveforms generated from a recent tomography image fit teleseismic travel times quite well and also some of the waveform distortions. Regional earthquake data, however, require substantial changes to the tomographic velocities. By modeling regional waveforms of the 2008 Nevada earthquake, we find that the uppermost mantle of the 1D reference model AK135, the reference velocity model used for most tomographic studies, is too fast for the western United States. Here, we replace AK135 with mT7, a modification of an older Basin-and-Range model T7. We present two hybrid velocity structures satisfying the waveform data based on modified tomographic images and conventional slab wisdom. We derive P and SH velocity structures down to 660 km along two cross sections through the JdF slab. Our results indicate that the JdF slab is subducted to a depth of 250 km beneath the Seattle region, and terminates at a shallower depth beneath Portland region of Oregon to the south. The slab is about 60 km thick and has a P velocity increase of 5% with respect to mT7. In order to fit waveform complexities of teleseismic Gulf of Mexico and South American events, a slab-like high-velocity anomaly with velocity increases of 3% for P and 7% for SH is inferred just above the 660 discontinuity beneath Nevada

Determination of earthquake focal depths and source time functions in central Asia using teleseismic P waveforms

We developed a new method to determine earthquake source time functions and focal depths. It uses theoretical Green's function and a time-domain deconvolution with positivity constraint to estimate the source time function from the teleseismic P waveforms. The earthquake focal depth is also determined in the process by using the time separations of the direct P and depth phases. We applied this method to 606 earthquakes between 1990 and 2005 in Central Asia. The results show that the Centroid Moment Tensor solutions, which are routinely computed for earthquake larger than M5.0 globally using very long period body and surface waves, systematically over-estimated the source depths and durations, especially for shallow events. Away from the subduction zone, most of the 606 earthquakes occurred within the top 20 km of crust. This shallow distribution of earthquakes suggests a high geotherm and a weak ductile lower crust in the region

Upper mantle P velocity structure beneath the Midwestern United States derived from triplicated waveforms

Upper mantle seismic velocity structures in both vertical and horizontal directions are key to understanding the structure and mechanics of tectonic plates. Recent deployment of the USArray Transportable Array (TA) in the Midwestern United States provides an extraordinary regional earthquake data set to investigate such velocity structure beneath the stable North American craton. In this paper, we choose an M_w5.1 Canadian earthquake in the Quebec area, which is recorded by about 400 TA stations, to examine the P wave structures between the depths of 150 km to 800 km. Three smaller Midwestern earthquakes at closer distance to the TA are used to investigate vertical and horizontal variations in P velocity between depths of 40 km to 150 km. We use a grid-search approach to find the best 1-D model, starting with the previously developed S25 regional model. The results support the existence of an 8° discontinuity in P arrivals caused by a negative velocity gradient in the lithosphere between depths of 40 km to 120 km followed by a small (∼1%) jump and then a positive gradient down to 165 km. The P velocity then decreases by 2% from 165 km to 200 km, and we define this zone as the regional lithosphere-asthenosphere boundary (LAB). Beneath northern profiles, waves reflected from the 410 discontinuity (410) are delayed by up to 1 s relative to those turning just below the 410, which we explain by an anomaly just above the discontinuity with P velocity reduced by ∼3%. The 660 discontinuity (660) appears to be composed of two smaller velocity steps with a separation of 16 km. The inferred low-velocity anomaly above 410 may indicate high water concentrations in the transition zone, and the complexity of the 660 may be related to Farallon slab segments that have yet to sink into the deep mantle

Department Budgeting and Resources

In this interactive session, we will look at the three keys to managing the budget and increasing resources for your department. First, chairs need to understand their budgets and how budgets are constructed. Next, chairs can manage their budgets, increase department productivity and work with their deans to invest resources in the department. Lastly, chairs can take steps to increase resource flow into their departmen

Lithospheric waveguide beneath the Midwestern United States; massive low-velocity zone in the lower crust

Variations in seismic velocities are essential in developing a better understanding of continental plate tectonics. Fortunately, the USArray has provided an excellent set of regional phases from the recent M5.6 Oklahoma earthquake (6 November 2011, Table 1) that can be used for such studies. Its strike-slip mechanism produced an extraordinary set of tangential recordings extending to the northern edge of the USArray. The crossover of the crustal slow S to the faster S_n phase is well observed. S_mS has a critical distance of around 2° and its first multiple, SmS^2, reaches critical angle near a distance of about 4°, and so on, until S_mS^n merges with the stronger crustal Love waves. These waveforms are modeled in the period band of 2–100 s by assuming a simple three-layer crust and a two-layer mantle, which allows a grid-search approach. Our results favor a 15 km thick low-velocity zone (LVZ) in the lower crust with an average shear velocity of less than 3.6 km/s. The short-period Lg waves (S waves, at periods of 0.5–2 s) travel with velocities near 3.5 km/s and decay with distance faster than high-frequency S_n (>5.0 Hz) which travels at a velocity of 4.6 km/s and persists to large distances. Although these short-period waveforms are not modeled, their amplitude and travel times can be explained by adding a small velocity jump just below the Moho with essentially no attenuation. P_n is equally strong but is complicated by the interference produced by the depth phase sP, but well modeled. The P velocities appear normal with no definitive LVZ. While these observations of S_n and P_n are common beneath most cratons, the lower crustal LVZ appears to be anomalous and maybe indicative of hydrous processes, possibly caused by the descending Farallon slab

Department Budgeting and Resources

Departments are faced with reduced appropriations and increased demands on their budgets. We will look at the three keys to respond to this fiscal environment and position your department for a brighter future: budget basics, management aimed at productivity and strategic investment, and increasing resource flow

Source Parameters of the Shallow 2012 Brawley Earthquake, Imperial Valley

Resolving earthquake parameters, especially depth, is difficult for events occurring within basins because of issues involved with separating source properties from propagational path effects. Here, we demonstrate some advantages of using a combination of teleseismic and regional waveform data to improve resolution following a bootstrapping approach. Local SS‐S differential arrivals from a foreshock are used to determine a local layered model which can then be used to model teleseismic depth phases: pP, sP, and sS. Using the cut‐and‐paste (CAP) method for which all strike (θ), dip (δ), rake (λ), and depth variations are sampled for several crustal models. We find that regional data prove the most reliable at fixing the strike, whereas the depth is better constrained by teleseismic data. Weighted solutions indicate a nearly pure strike‐slip mechanism (θ=59°±1°) with a centroid depth of about 4.0 km and an M_w of 5.4 for the mainshock of the 2012 Brawley earthquake

GOES dynamic propagation of attitude

The spacecraft in the next series of Geostationary Operational Environmental Satellites (GOES-Next) are Earth pointing and have 5-year mission lifetimes. Because gyros can be depended on only for a few years of continuous use, they will be turned off during routine operations. This means attitude must, at times, be determined without benefit of gyros and, often, using only Earth sensor data. To minimize the interruption caused by dumping angular momentum, these spacecraft have been designed to reduce the environmental torque acting on them and incorporate an adjustable solar trim tab for fine adjustment. A new support requirement for GOES-Next is that of setting the solar trim tab. Optimizing its setting requires an estimate of the unbalanced torque on the spacecraft. These two requirements, determining attitude without gyros and estimating the external torque, are addressed by replacing or supplementing the gyro propagation with a dynamic one, that is, one that integrates the rigid body equations of motion. By processing quarter-orbit or longer batches, this approach takes advantage of roll-yaw coupling to observe attitude completely without Sun sensor data. Telemetered momentum wheel speeds are used as observations of the unbalanced external torques. GOES-Next provides a unique opportunity to study dynamic attitude propagation. The geosynchronous altitude and adjustable trim tab minimize the external torque and its uncertainty, making long-term dynamic propagation feasible. This paper presents the equations for dynamic propagation, an analysis of the environmental torques, and an estimate of the accuracies obtainable with the proposed method

Hidden hotspot track beneath the eastern United States

Hotspot tracks are thought to be the surface expressions of tectonic plates moving over upwelling mantle plumes, and are characterized by volcanic activity that is age progressive. At present, most hotspot tracks are observed on oceanic or thin continental lithosphere. For old, thick continental lithosphere, such as the eastern United States, hotspot tracks are mainly inferred from sporadic diamondiferous kimberlites putatively sourced from the deep mantle. Here we use seismic waveforms initiated by the 2011 M_w 5.6 Virginia earthquake, recorded by the seismic observation network USArray, to analyse the structure of the continental lithosphere in the eastern United States. We identify an unexpected linear seismic anomaly in the lower lithosphere that has both a reduced P-wave velocity and high attenuation, and which we interpret as a hotspot track. The anomaly extends eastwards, from Missouri to Virginia, cross-cutting the New Madrid rift system, and then bends northwards. It has no clear relationship with the surface geology, but crosses a 75-million-year-old kimberlite in Kentucky. We use geodynamical modelling to show that an upwelling thermal mantle plume that interacts with the base of continental lithosphere can produce the observed seismic anomaly. We suggest that the hotspot track could be responsible for late Mesozoic reactivation of the New Madrid rift system and seismicity of the eastern United States