Understanding of Contemporary Regional Sea-Level Change and the Implications for the Future.

Global sea level provides an important indicator of the state of the warming climate, but changes in regional sea level are most relevant for coastal communities around the world. With improvements to the sea-level observing system, the knowledge of regional sea-level change has advanced dramatically in recent years. Satellite measurements coupled with in situ observations have allowed for comprehensive study and improved understanding of the diverse set of drivers that lead to variations in sea level in space and time. Despite the advances, gaps in the understanding of contemporary sea-level change remain and inhibit the ability to predict how the relevant processes may lead to future change. These gaps arise in part due to the complexity of the linkages between the drivers of sea-level change. Here we review the individual processes which lead to sea-level change and then describe how they combine and vary regionally. The intent of the paper is to provide an overview of the current state of understanding of the processes that cause regional sea-level change and to identify and discuss limitations and uncertainty in our understanding of these processes. Areas where the lack of understanding or gaps in knowledge inhibit the ability to provide the needed information for comprehensive planning efforts are of particular focus. Finally, a goal of this paper is to highlight the role of the expanded sea-level observation network-particularly as related to satellite observations-in the improved scientific understanding of the contributors to regional sea-level change.

GPS Imaging of vertical land motion in California and Nevada: Implications for Sierra Nevada uplift.

We introduce Global Positioning System (GPS) Imaging, a new technique for robust estimation of the vertical velocity field of the Earth's surface, and apply it to the Sierra Nevada Mountain range in the western United States. Starting with vertical position time series from Global Positioning System (GPS) stations, we first estimate vertical velocities using the MIDAS robust trend estimator, which is insensitive to undocumented steps, outliers, seasonality, and heteroscedasticity. Using the Delaunay triangulation of station locations, we then apply a weighted median spatial filter to remove velocity outliers and enhance signals common to multiple stations. Finally, we interpolate the data using weighted median estimation on a grid. The resulting velocity field is temporally and spatially robust and edges in the field remain sharp. Results from data spanning 5-20 years show that the Sierra Nevada is the most rapid and extensive uplift feature in the western United States, rising up to 2 mm/yr along most of the range. The uplift is juxtaposed against domains of subsidence attributable to groundwater withdrawal in California's Central Valley. The uplift boundary is consistently stationary, although uplift is faster over the 2011-2016 period of drought. Uplift patterns are consistent with groundwater extraction and concomitant elastic bedrock uplift, plus slower background tectonic uplift. A discontinuity in the velocity field across the southeastern edge of the Sierra Nevada reveals a contrast in lithospheric strength, suggesting a relationship between late Cenozoic uplift of the southern Sierra Nevada and evolution of the southern Walker Lane.

MIDAS robust trend estimator for accurate GPS station velocities without step detection.

Automatic estimation of velocities from GPS coordinate time series is becoming required to cope with the exponentially increasing flood of available data, but problems detectable to the human eye are often overlooked. This motivates us to find an automatic and accurate estimator of trend that is resistant to common problems such as step discontinuities, outliers, seasonality, skewness, and heteroscedasticity. Developed here, Median Interannual Difference Adjusted for Skewness (MIDAS) is a variant of the Theil-Sen median trend estimator, for which the ordinary version is the median of slopes vij = (xj-xi )/(tj-ti ) computed between all data pairs i > j. For normally distributed data, Theil-Sen and least squares trend estimates are statistically identical, but unlike least squares, Theil-Sen is resistant to undetected data problems. To mitigate both seasonality and step discontinuities, MIDAS selects data pairs separated by 1 year. This condition is relaxed for time series with gaps so that all data are used. Slopes from data pairs spanning a step function produce one-sided outliers that can bias the median. To reduce bias, MIDAS removes outliers and recomputes the median. MIDAS also computes a robust and realistic estimate of trend uncertainty. Statistical tests using GPS data in the rigid North American plate interior show ±0.23 mm/yr root-mean-square (RMS) accuracy in horizontal velocity. In blind tests using synthetic data, MIDAS velocities have an RMS accuracy of ±0.33 mm/yr horizontal, ±1.1 mm/yr up, with a 5th percentile range smaller than all 20 automatic estimators tested. Considering its general nature, MIDAS has the potential for broader application in the geosciences.

Uplift and seismicity driven by groundwater depletion in central California.

Groundwater use in California's San Joaquin Valley exceeds replenishment of the aquifer, leading to substantial diminution of this resource and rapid subsidence of the valley floor. The volume of groundwater lost over the past century and a half also represents a substantial reduction in mass and a large-scale unburdening of the lithosphere, with significant but unexplored potential impacts on crustal deformation and seismicity. Here we use vertical global positioning system measurements to show that a broad zone of rock uplift of up to 1-3 mm per year surrounds the southern San Joaquin Valley. The observed uplift matches well with predicted flexure from a simple elastic model of current rates of water-storage loss, most of which is caused by groundwater depletion. The height of the adjacent central Coast Ranges and the Sierra Nevada is strongly seasonal and peaks during the dry late summer and autumn, out of phase with uplift of the valley floor during wetter months. Our results suggest that long-term and late-summer flexural uplift of the Coast Ranges reduce the effective normal stress resolved on the San Andreas Fault. This process brings the fault closer to failure, thereby providing a viable mechanism for observed seasonality in microseismicity at Parkfield and potentially affecting long-term seismicity rates for fault systems adjacent to the valley. We also infer that the observed contemporary uplift of the southern Sierra Nevada previously attributed to tectonic or mantle-derived forces is partly a consequence of human-caused groundwater depletion.