Adrian A. Borsa

Ongoing drought-induced uplift in the western United States

Crustal rebound from water drawdown The ongoing drought across the western United States has taken a toll on underground water storage. Borsa et al. use almost imperceptible crustal uplift to estimate the regional water depletion from the drought. Inverting GPS data maps the impact of the drought on local aquifers over the past few years. The deficit so far in the western United States adds up to 240 gigatons of water, the equivalent of a 10-cm layer across the region. Certain areas of California have fared much worse, with local depletions up to five times the regional average. Science, this issue p. 1587 GPS measurements of crustal rebound in the western U.S. quantify drought-induced regional water depletion. The western United States has been experiencing severe drought since 2013. The solid earth response to the accompanying loss of surface and near-surface water mass should be a broad region of uplift. We use seasonally adjusted time series from continuously operating global positioning system stations to measure this uplift, which we invert to estimate mass loss. The median uplift is 5 millimeters (mm), with values up to 15 mm in California’s mountains. The associated pattern of mass loss, ranging up to 50 centimeters (cm) of water equivalent, is consistent with observed decreases in precipitation and streamflow. We estimate the total deficit to be ~240 gigatons, equivalent to a 10-cm layer of water over the entire region, or the annual mass loss from the Greenland Ice Sheet.

High-Resolution Topography along Surface Rupture of the 16 October 1999 Hector Mine, California, Earthquake (Mw 7.1) from Airborne Laser Swath Mapping

In order to document surface rupture associated with the Hector Mine earthquake, in particular, the area of maximum slip and the deformed surface of Lavic Lake playa, we acquired high-resolution data using relatively new topographic-mapping methods. We performed a raster-laser scan of the main surface breaks along the entire rupture zone, as well as along an unruptured portion of the Bullion fault. The image of the ground surface produced by this method is highly detailed, comparable to that obtained when geologists make particularly detailed site maps for geomorphic or paleoseismic studies. In this case, however, for the first time after a surface-rupturing earthquake, the detailed mapping is along the entire fault zone rather than being confined to selected sites. These data are geodetically referenced, using the Global Positioning System, thus enabling more accurate mapping of the rupture traces. In addition, digital photographs taken along the same flight lines can be overlaid onto the precise topographic data, improving terrain visualization. We demonstrate the potential of these techniques for measuring fault-slip vectors.

Modeling the topography of the salar de Uyuni, Bolivia, as an equipotential surface of Earth's gravity field

The salar de Uyuni is a massive dry salt lake that lies at the lowest point of an internal drainage basin in the Bolivian Altiplano. A kinematic GPS survey of the salar in September 2002 found a topographic range of only 80 cm over a 54 × 45 km area and subtle surface features that appeared to correlate with mapped gravity. In order to confirm the correlation between topography and gravity/geopotential, we use local gravity measurements and the EGM96 global geopotential model to construct a centimeter-level equipotential surface corresponding to the elevation of the salar. Our comparison of GPS survey elevations with the equipotential surface estimate shows that 63% of the variance of the GPS elevations can be explained by equipotential surface undulations (and long-wavelength error) in the EGM96 model alone, with an additional 30% explained by the shorter-wavelength equipotential surface derived from local gravity. In order to establish a physical connection between topography and the geopotential, we also develop and test a simple surface process model that redistributes salt via the dissolution, transport, and redeposition of salt by precipitated water. Forcing within the model pushes the system to evolve toward constant water depth, with the salt surface approximating the shape of the local equipotential surface. Since the model removes almost all topographic relief with respect to the equipotential surface within a matter of decades, it appears that observed (~5 cm amplitude, ~5 km wavelength) residual topography is actively maintained by a process independent of gravity-driven fluid flow.

Optimization of legacy lidar data sets for measuring near-field earthquake displacements

Airborne lidar (light detection and ranging) topography, acquired before and after an earthquake, can provide an estimate of the coseismic surface displacement field by differencing the preevent and postevent lidar point clouds. However, estimated displacements can be contaminated by the presence of large systematic errors in either of the point clouds. We present three-dimensional displacements obtained by differencing airborne lidar point clouds collected before and after the El Mayor–Cucapah earthquake, a Mw 7.2 earthquake that occurred in 2010. The original surface displacement estimates contained large, periodic artifacts caused by systematic errors in the preevent lidar data. Reprocessing the preevent data, detailed herein, removed a majority of these systematic errors that were largely due to misalignment between the scanning mirror and the outgoing laser beam. The methodology presented can be applied to other legacy airborne laser scanning data sets in order to improve change estimates from temporally spaced lidar acquisitions.

A comparison of long‐term changes in seismicity at The Geysers, Salton Sea, and Coso geothermal fields

Geothermal energy is an important source of renewable energy, yet its production is known to induce seismicity. Here we analyze seismicity at the three largest geothermal fields in California: The Geysers, Salton Sea, and Coso. We focus on resolving the temporal evolution of seismicity rates, which provides important observational constraints on how geothermal fields respond to natural and anthropogenic loading. We develop an iterative, regularized inversion procedure to partition the observed seismicity rate into two components: (1) the interaction rate due to earthquake-earthquake triggering and (2) the smoothly varying background rate controlled by other time-dependent stresses, including anthropogenic forcing. We apply our methodology to compare long-term changes in seismicity to monthly records of fluid injection and withdrawal. At The Geysers, we find that the background seismicity rate is highly correlated with fluid injection, with the mean rate increasing by approximately 50% and exhibiting strong seasonal fluctuations following construction of the Santa Rosa pipeline in 2003. In contrast, at both Salton Sea and Coso, the background seismicity rate has remained relatively stable since 1990, though both experience short-term rate fluctuations that are not obviously modulated by geothermal plant operation. We also observe significant temporal variations in Gutenberg-Richter b value, earthquake magnitude distribution, and earthquake depth distribution, providing further evidence for the dynamic evolution of stresses within these fields. The differing field-wide responses to fluid injection and withdrawal may reflect differences in in situ reservoir conditions and local tectonics, suggesting that a complex interplay of natural and anthropogenic stressing controls seismicity within Californias geothermal fields.

Rapid Determination of Near‐Fault Earthquake Deformation Using Differential LiDAR

Improved near‐field measurements of earthquake slip and deformation patterns have the potential for expanding our understanding of fault behavior and the relationship of active faulting to topography. Current techniques for obtaining these measurements—including field observation, Global Navigation Satellite Systems displacement estimation, and optical or radar remote sensing—have limitations that can be mitigated by the inclusion of results from differential airborne Light Detection and Ranging (LiDAR) analysis of the rupture zone. The 2005 airborne LiDAR survey of the southern San Andreas, San Jacinto, and Banning faults (the B4 survey) mapped 1100 km of the most seismically active fault systems in southern California for the purpose of providing a baseline for determining slip from a future earthquake. We used the B4 survey to develop a processing algorithm that yields rapid estimates of near‐fault ground deformation using simultaneous cross correlation of both topography and backscatter intensity from pre‐earthquake and simulated postearthquake LiDAR datasets. We show robust recovery of the direction and magnitude of an applied synthetic slip of 5 m in the horizontal and 0.5 m in the vertical within our area of study, with clear discrimination between areas with and without applied slip. We also successfully recovered more complex deformation from a modeled fault stepover in the same study area. Our results indicate that we should be able to recover slip to accuracies of better than 20 cm in the horizontal and 1 cm in the vertical, at a spatial resolution of ≤15 m for LiDAR datasets with sample densities as low as 0.5 points/m2.

Western US intermountain seismicity caused by changes in upper mantle flow

Understanding the causes of intraplate earthquakes is challenging, as it requires extending plate tectonic theory to the dynamics of continental deformation. Seismicity in the western United States away from the plate boundary is clustered along a meandering, north–south trending ‘intermountain’ belt. This zone coincides with a transition from thin, actively deforming to thicker, less tectonically active crust and lithosphere. Although such structural gradients have been invoked to explain seismicity localization, the underlying cause of seismicity remains unclear. Here we show results from improved mantle flow models that reveal a relationship between seismicity and the rate change of ‘dynamic topography’ (that is, vertical normal stress from mantle flow). The associated predictive skill is greater than that of any of the other forcings we examined. We suggest that active mantle flow is a major contributor to seismogenic intraplate deformation, while gravitational potential energy variations have a minor role. Seismicity localization should occur where convective changes in vertical normal stress are modulated by lithospheric strength heterogeneities. Our results on deformation processes appear consistent with findings from other mobile belts, and imply that mantle flow plays a significant and quantifiable part in shaping topography, tectonics, and seismic hazard within intraplate settings.

Did stresses from the Cerro Prieto Geothermal Field influence the El Mayor‐Cucapah rupture sequence?

The Mw 7.2 El Mayor-Cucapah (EMC) earthquake ruptured a complex fault system in northern Baja California that was previously considered inactive. The Cerro Prieto Geothermal Field (CPGF), site of the worlds second largest geothermal power plant, is located approximately 15 km to the northeast of the EMC hypocenter. We investigate whether anthropogenic fluid extraction at the CPGF caused a significant perturbation to the stress field in the EMC rupture zone. We use Advanced Land Observing Satellite interferometric synthetic aperture radar data to develop a laterally heterogeneous model of fluid extraction at the CPGF and estimate that this extraction generates positive Coulomb stressing rates of order 15 kPa/yr near the EMC hypocenter, a value which exceeds the local tectonic stressing rate. Although we cannot definitively conclude that production at the CPGF triggered the EMC earthquake, its influence on the local stress field is substantial and should not be neglected in local seismic hazard assessments.

ICESAT/GLAS Altimetry Measurements: Received Signal Dynamic Range and Saturation Correction

NASA’s Ice, Cloud, and land Elevation Satellite (ICESat), which operated between 2003 and 2009, made the first satellite-based global lidar measurement of earth’s ice sheet elevations, sea-ice thickness, and vegetation canopy structure. The primary instrument on ICESat was the Geoscience Laser Altimeter System (GLAS), which measured the distance from the spacecraft to the earth’s surface via the roundtrip travel time of individual laser pulses. GLAS utilized pulsed lasers and a direct detection receiver consisting of a silicon avalanche photodiode and a waveform digitizer. Early in the mission, the peak power of the received signal from snow and ice surfaces was found to span a wider dynamic range than anticipated, often exceeding the linear dynamic range of the GLAS 1064-nm detector assembly. The resulting saturation of the receiver distorted the recorded signal and resulted in range biases as large as ~50 cm for ice- and snow-covered surfaces. We developed a correction for this “saturation range bias” based on laboratory tests using a spare flight detector, and refined the correction by comparing GLAS elevation estimates with those derived from Global Positioning System surveys over the calibration site at the salar de Uyuni, Bolivia. Applying the saturation correction largely eliminated the range bias due to receiver saturation for affected ICESat measurements over Uyuni and significantly reduced the discrepancies at orbit crossovers located on flat regions of the Antarctic ice sheet.