[1] We present the first results from a dense network of 36 campaign and 46 continuous GPS stations located in the Eastern Transverse Ranges Province (ETR), a transition zone between the southernmost San Andreas fault (SSAF) and eastern California shear zone (ECSZ). We analyzed the campaign data together with available data from continuous GPS stations for the period 1994-2009. We used the GPS velocity estimates to constrain elastic block models to investigate fault-loading rates representing four hypotheses characterized by different fault-block geometries. Fault-block scenarios include blocks bounded by the east-striking left-lateral Pinto Mountain, Blue Cut, and Chiriaco faults of the ETR; blocks bounded by a right-lateral north-northwest striking structure (the "Landers-Mojave earthquake line") that cuts obliquely across the ETR and mapped Mojave Desert faults; and combinations of these end-member hypotheses. Each model implies significantly different active fault geometries, block rotation rates, and slip rates for ETR and ECSZ structures. All models suggest that SSAF slip rate varies appreciably along strike, generally consistent with rates derived from tectonic geomorphology and paleoseismology, with a maximum of ∼23 mm/yr right-lateral along the southernmost Coachella Valley strand, decreasing systematically to <10 mm/yr right-lateral through the San Gorgonio Pass region. Slip rate estimates for the San Jacinto fault are ∼12 mm/yr for all models tested. All four models fit the data equally well in a statistical sense. Qualitative comparison among models and consideration of geologic slip rates and other independent data reveals strengths and weaknesses of each model.

We present numerical modeling of the forces acting on the base of the crust caused by small‐scale convection of the upper mantle in southern California. Three‐dimensional upper mantle shear wave velocity structure is mapped to three‐dimensional density structure that is used to load a finite element model of instantaneous upper mantle flow with respect to a rigid crust, providing an estimate of the tractions acting on the base of the crust. Upwelling beneath the southern Walker Lane Belt and Salton Trough region and downwelling beneath the southern Great Valley and eastern and western Transverse Ranges dominate the upper mantle flow and resulting crustal tractions. Divergent horizontal and upward directed vertical tractions create a tensional to transtensional crustal stress state in the Walker Lane Belt and Salton Trough, consistent with transtensional tectonics in these areas. Convergent horizontal and downward directed vertical tractions in the Transverse Ranges cause approximately N–S crustal compression, consistent with active shortening and transpressional deformation near the “Big Bend” of the San Andreas fault. Model predictions of crustal dilatation and the forces acting on the Mojave block compare favorably with observations suggesting that small‐scale upper mantle convection provides an important contribution to the sum of forces driving transpressional crustal deformation in southern California. Accordingly, the obliquity of the San Andreas fault with respect to plate motions may be considered a consequence, rather than a cause, of contractional deformation in the Transverse Ranges, itself driven by downwelling in the upper mantle superimposed on shear deformation caused by relative Pacific–North American plate motion.

The 4 April 2010 M w 7.2 El Mayor-Cucapah (EMC) earthquake provides the best opportunity to date to study the lithospheric response to a large (>M6) magnitude earthquake in the Salton Trough region through analysis of Global Positioning System (GPS) data. In conjunction with the EarthScope Plate Boundary Observatory (PBO), we installed six new continuous GPS stations in the months following the EMC earthquake to increase station coverage in the epicentral region of northern Baja California, Mexico. We modeled the pre-EMC deformation field using available campaign and continuous GPS data for southern California and northern Baja California and inferred a pre-EMC secular rate at each new station location. Through direct comparison of the pre-and post-EMC secular rates, we calculate long-term changes associated with viscoelastic relaxation in the Salton Trough region. We fit these velocity changes using numerical models employing an elastic upper crustal layer underlain by a viscoelastic lower crustal layer and a mantle half-space. Forward models that produce the smallest weighted sum of squared residuals have an upper mantle viscosity in the range 4-6 × 1018 Pa s and a less well-resolved lower crustal viscosity in the range 2 × 10 19 to 1 × 10 22 Pa s. A high-viscosity lower crust, despite high heat flow in the Salton Trough region, is inconsistent with felsic composition and might suggest accretion of mafic lower crust associated with crustal spreading obscured by thick sedimentary cover.

Campaign GPS data collected from 2002 to 2014 result in 41 new site velocities from the San Bernardino Mountains and vicinity. We combined these velocities with 93 continuous GPS velocities and 216 published velocities to obtain a velocity profile across the Pacific-North America plate boundary through the San Bernardino Mountains. We modeled the plate boundary-parallel, horizontal deformation with 5-14 parallel and one obliquely oriented screw dislocations within an elastic half-space. Our rate for the San Bernardino strand of the San Andreas Fault (6.5 ± 3.6 mm/yr) is consistent with recently published latest Quaternary rates at the 95% confidence level and is slower than our rate for the San Jacinto Fault (14.1 ± 2.9 mm/yr). Our modeled rate for all faults of the Eastern California Shear Zone (ECSZ) combined (15.7 ± 2.9 mm/yr) is faster than the summed latest Quaternary rates for these faults, even when an estimate of permanent, off-fault deformation is included. The rate discrepancy is concentrated on faults near the 1992 Landers and 1999 Hector Mine earthquakes; the geodetic and geologic rates agree within uncertainties for other faults within the ECSZ. Coupled with the observation that postearthquake deformation is faster than the pre-1992 deformation, this suggests that the ECSZ geodetic-geologic rate discrepancy is directly related to the timing and location of these earthquakes and is likely the result of viscoelastic deformation in the mantle that varies over the timescale of an earthquake cycle, rather than a redistribution of plate boundary slip at a timescale of multiple earthquake cycles or longer.

[1] Vegetation bands are periodic bands of vegetation, separated by interband spaces devoid of vegetation, oriented parallel to the topographic contour in some gently sloping arid to semiarid environments. Models of vegetation band formation attribute their formation to positive feedbacks among vegetation density, soil porosity/permeability, and infiltration rates. Here we present an alternative model based on field measurements at our study sites in southern Nevada. In this model, interband spaces between vegetation bands form because topographic mounds beneath vegetation bands detain water upslope from vegetation bands, leading to hydrologic and sedimentologic conditions that inhibit the survival of plants in interband spaces. We used terrestrial laser scanning (TLS) to create high-resolution ($10 cm 2 /pixel) raster data sets of bare-earth topography and canopy height for four study sites. Analyses of the TLS data, in addition to measurements of soil shear strength and particle size, document the potential for detention in interband spaces and a near-inverse proportionality between band spacing and regional slope. We describe a cellular automaton model (herein called model 1) for vegetation band formation that includes just two user-defined parameters and that generates vegetation bands similar to those at our field sites, including the inverse proportionality between spacing and regional slope. A second model (model 2) accurately predicts the width of vegetation bands in terms of the number and spacing of plants and the geometry of individual plant mounds. We also present a GIS-based analysis that predicts where bands occur within a region based on topographic and hydroclimatic controls.

No abstract

Assessing fault-slip rates in diffuse plate boundary systems such as the San Andreas fault in southern California is critical both to characterize seismic hazards and to understand how different fault strands work together to accommodate plate boundary motion. In places such as San Gorgonio Pass, the geometric complexity of numerous fault strands interacting in a small area adds an extra obstacle to understanding the rupture potential and behavior of each individual fault. To better understand partitioning of fault-slip rates in this region, we build a new set of elastic fault-block models that test 16 different model fault geometries for the area. These models build on previous studies by incorporating updated campaign GPS measurements from the San Bernardino Mountains and Eastern Transverse Ranges into a newly calculated GPS velocity field that has been removed of long- and short-term postseismic displacements from 12 past large-magnitude earthquakes to estimate model fault-slip rates. Using this postseismic-reduced GPS velocity field produces a best-fitting model geometry that resolves the long-standing geologic-geodetic slip-rate discrepancy in the Eastern California shear zone when off-fault deformation is taken into account, yielding a summed slip rate of 7.2 ± 2.8 mm/yr. Our models indicate that two active strands of the San Andreas system in San Gorgonio Pass are needed to produce sufficiently low geodetic dextral slip rates to match geologic observations. Lastly, results suggest that postseismic deformation may have more of a role to play in affecting the loading of faults in southern California than previously thought.

We conducted a study to quantify the amount of pseudorange multipath at continuous Global Positioning System (CGPS) stations in the Mexican territory. These CGPS stations serve as reference stations enabling rapid high-precision three-dimensional positioning capabilities, supporting a number of commercial and public safety applications. We studied CGPS data from a large number of publicly available networks spanning Mexico. These include the RGNA (National Active Geodetic Network) administered by INEGI (National Institute of Statistics and Geography), the PBO network (Plate Boundary Observatory) funded by the National Science Foundation (NSF) and operated by UNAVCO (University NAVstar Consortium), the Southern California Integrated GPS Network (SCIGN), which is a collaboration effort of the United States Geological Survey (USGS), Scripps Institution of Oceanography and the Jet Propulsion Laboratory (JPL), the UNAM network, operated by the National Seismological System (SSN) and the Institute of Geophysics of the National Autonomous University of Mexico (UNAM), the Suominet Geodetic Network (SNG) and the CORS (Continuously Operating Reference Station) network, operated by the Federal Aviation Administration (FAA). We evaluated a total of 53 CGPS stations, where dual-frequency geodetic-grade receivers collected GPS data continuously during the period from 1994 to 2012. Despite carefully selected locations, all GPS stations are, to some extent, affected by the presence of signal multipath. For GPS network users that rely on pseudorange observables, the existence of pseudorange multipath could be a critical source of error depending on the time scale of the application. Thus, to identify the most and the least affected GPS stations, we analyzed the averaged daily root mean square pseudorange multipath variations (MP1-RMS and MP2-RMS) for all feasible satellites tracked by the CGPS networks. We investigated the sources of multipath, including changes associated with hardware replacement (i.e., receiver and antenna type) and receiver firmware upgrades.