Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Geophysical Research”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53Linked to original sources
USGS2019

Slow-growing and extended-duration seismicity swarms: Reactivating joints or foliations in the Cahuilla Valley Pluton, Central Peninsular Ranges, Southern California

Three prolific earthquake swarms and numerous smaller ones have occurred since 1980 in the Mesozoic igneous plutonic rocks of the Perris block of the Peninsular Ranges, Southern California. The major swarms occurred in 1980–1981, 1983–1984, and 2016–2018, with the latest swarm still ongoing. These swarms have no clear mainshock, with the largest events of M L 3.6, M L 3.7, and M w 4.4. Each successive swarm had larger cumulative seismic moment release with about 314 and 411 events of M ≥ 1.5, while the third swarm has produced about 451 events of M ≥ 1.5 (as of September 2018). The concurrent strike‐slip faulting occurred on north to northwest striking planes but with no orthogonal northeast trending seismicity alignments. These shallow swarms are probably driven by intrablock Pacific‐North America plate boundary stress loading of the two bounding major late Quaternary strike‐slip faults, the Elsinore and San Jacinto faults. The state of stress within the Cahuilla Valley pluton has a ~40° angle between the maximum principal stress and the average trend of the swarms, suggesting that migrating pore fluid pressures aid in the formation and growth of zones of weakness. These swarms, which last more than 600 days each, exhibit clear bilateral spatial migration for distances of up to ~7–8 km and reach their full length in about 20 months. The slow spatial‐temporal development of the swarms corresponds to a fluid diffusivity of 0.006 to 0.01 m 2 /s, consistent with very low permeability rocks as expected for this block. There is no geodetic or other evidence for a slow slip event driving the swarms.

California↗
USGS2019

Heat flow in the Western Arctic Ocean (Amerasian Basin)

From 1963 to 1973 the U.S. Geological Survey (USGS) measured heat flow at 356 sites in the Amerasian Basin (Western Arctic Ocean) from a drifting ice island (T-3). The resulting measurements, which are unevenly distributed on Alpha-Mendeleev Ridge (AMR) and in Canada and Nautilus basins, greatly expand available heat flow data for the Arctic Ocean. Average T-3 heat flow is ~54.7 ± 11.3 mW m-2, and Nautilus Basin, including Mendeleev Plain, is the only well-surveyed area (~13% of data) with significantly higher average heat flow (63.8 mW m-2). Heat flow and bathymetry are not correlated at a large scale, and turbiditic surficial sediments (Canada and Nautilus basins) have higher heat flow than the sediments that blanket the AMR. Thermal gradients are mostly near-linear, implying that conductive heat transport dominates and that near-seafloor sediments are in thermal equilibrium with overlying bottom waters. Combining the heat flow data with modern seismic imagery suggests that some of the observed heat flow variability may be explained by local changes in sediment thickness or lithology or the presence of basement faults that channel circulating seawater. A thermal model that incorporates thermal conductivity variations along a profile from Canada Basin (thick sediment on mostly oceanic crust) to Alpha Ridge (thin sediment over thick magmatic units associated with the High Arctic Large Igneous Province) predicts heat flow lower than that observed on Alpha Ridge. This, along with other observations, implies that circulating fluids modulate conductive heat flow and contribute to high variability in the T-3 dataset. .

Journal of Geophysical Research B: Solid Earth↗
USGS2020

Near-fault velocity spectra from laboratory failures and their relation to natural ground motion

We compared near-fault velocity spectra recorded during laboratory experiments to that of natural earthquakes. We fractured crystalline rock samples at room temperature and intermediate confining pressure (50 MPa). Subsequent slip events were generated on the fracture surfaces under higher confinement (300 MPa). Velocity spectra from rock fracture resemble the inverse frequency (1/ f ) decay of natural earthquake velocity. This spectrum can be attributed to fault creation via seismic fracturing over a wide range of spatial scales. In contrast, subsequent slips on the rough fracture surfaces are depleted in high frequency energy and falloff approximately as 1/ f 2 . The 1/ f 2 spectrum is more consistent with a slider-block model obeying static-kinetic friction than a natural earthquake. The depleted high frequency content precludes the rough fault experiments from being directly analogous to natural sources. The suppression of high frequencies may have resulted from two possible factors: (1) the presence of a well-developed shear zone and coseismic damping of the fault motion by dissipation within it or, in our favored interpretation, (2) a smaller amount of energy dissipated by shearing relative to the total energy release at elevated confining pressure. In context of the latter explanation, a unifying concept that applies to these experiments, earthquakes, ground motion, and models of complex radiated motion is that high frequency radiated energy is relatively enhanced when total energy release is nearly balanced within the source region by dissipative processes. This near-critical energy release condition can be accessed at low normal stress in laboratory experiments.

Journal of Geophysical Research↗
USGS2020

How often can Earthquake Early Warning systems alert sites with high intensity ground motion?

Although numerous Earthquake Early Warning (EEW) algorithms have been developed we still lack a detailed understanding of how often and under what circumstances useful ground motion alerts can be provided to end-users. Here we analyze the alerting performance of the PLUM, EPIC and FinDer algorithms by running them retrospectively on the seismic strong motion data of the 219 earthquakes in Japan since 1996 that exceeded Modified Mercalli Intensity (MMI) of 4.5 on at least 10 sites (Mw 4.5-9.1). Our analysis suggests that, irrespective of the algorithm, EEW end-users should be prepared that EEW can often but not always provide useful ground motion alerts. A majority of sites with moderate-strong ground motion (MMI 5-6) can generally get at least a few seconds of warning time from all algorithms. If such shaking is caused by a shallow crustal event, around 50% of such sites receive alerts with warning times >5 s. Many sites with severe-extreme ground motion (MMI >=8) can be alerted successfully in the case of very large offshore earthquakes, but less than 20% can be alerted ahead of time if such shaking is caused by a shallow crustal event. Our results provide detailed quantitative insight into the expected alerting performance for EEW algorithms under realistic conditions. The main caveat is that the largest shallow crustal event in our data set has Mw7.0, i.e. the data set does not contain very large strike slip events.

Journal of Geophysical Research↗
USGS2019

Shear velocity structure from ambient noise and teleseismic surface wave tomography in the Cascades around Mount St. Helens

Mount St. Helens (MSH) lies in the forearc of the Cascades where conditions should be too cold for volcanism. To better understand thermal conditions and magma pathways beneath MSH, data from a dense broadband array are used to produce high‐resolution tomographic images of the crust and upper mantle. Rayleigh‐wave phase‐velocity maps and three‐dimensional images of shear velocity ( Vs ), generated from ambient noise and earthquake surface waves, show that west of MSH the middle‐lower crust is anomalously fast (3.95 ± 0.1 km/s), overlying an anomalously slow uppermost mantle (4.0–4.2 km/s). This combination renders the forearc Moho weak to invisible, with crustal velocity variations being a primary cause; fast crust is necessary to explain the absent Moho. Comparison with predicted rock velocities indicates that the fast crust likely consists of gabbros and basalts of the Siletzia terrane, an accreted oceanic plateau. East of MSH where magmatism is abundant, middle‐lower crust Vs is low (3.45–3.6 km/s), consistent with hot and potentially partly molten crust of more intermediate to felsic composition. This crust overlies mantle with more typical wave speeds, producing a strong Moho. The sharp boundary in crust and mantle Vs within a few kilometers of the MSH edifice correlates with a sharp boundary from low heat flow in the forearc to high arc heat flow and demonstrates that the crustal terrane boundary here couples with thermal structure to focus lateral melt transport from the lower crust westward to arc volcanoes.

Washington↗
USGS2019

Crustal magmatism and anisotropy beneath the Arabian Shield - A cautionary tale

Volcanism in Saudi Arabia includes a historic eruption close to the holy city of Al Madinah. As part of a volcanic hazard assessment of this area, magnetotelluric (MT) data were collected to investigate the structural setting, the distribution of melt within the crust, and the mantle source of volcanism. Interpretation of a new 3‐D resistivity model includes a shallow graben beneath thin lava fields (Harrats), a melt‐free upper crust, and decompression melting in the asthenosphere below thin lithosphere. Within the lower crust the model images elongate conductivity anomalies, one of which was attributed in a previous MT study to melt. The regional MT data, combined with perspective from geology and geophysical modeling, suggest the lower crust is anisotropic with no interconnected melt zones. These divergent interpretations have distinct hazard implications and highlight the importance of large survey aperture and anisotropic modeling to MT studies of volcanic regions. Lower‐crustal anisotropy extends beyond the Harrat, with the most conductive direction oriented N10°E and a factor of 3–5, determined from 2‐D anisotropic inversion, between the most and least conductive directions. The enhanced conductivity is likely due to interconnected grain boundary graphite, while the anisotropy direction reflects either frozen‐in fabric from Neoproterozoic stabilization of the Arabian Shield or modern ductile deformation driven by channelized asthenospheric flow coupled through a thin rigid mantle lid. Asthenospheric melt is interpreted to transect the crust primarily through diking, with limited melt storage and short residence times in the crust.

Arabian Shield↗
USGS2020

Structural control on megathrust rupture and slip behavior: Insights from the 2016 Mw 7.8 Pedernales Ecuador earthquake

The heterogeneous seafloor topography of the Nazca Plate as it enters the Ecuador subduction zone provides an opportunity to document the influence of seafloor roughness on slip behavior and megathrust rupture. The 2016 Mw 7.8 Pedernales Ecuador earthquake was followed by a rich and active postseismic sequence. An internationally coordinated rapid response effort installed a temporary seismic network to densify coastal stations of the permanent Ecuadorian national seismic network. A combination of 82 onshore short and intermediate period and broadband seismic stations and six ocean bottom seismometers recorded the postseismic Pedernales sequence for over a year after the mainshock. A robust earthquake catalog combined with calibrated relocations for a subset of magnitude ≥4 earthquakes shows pronounced spatial and temporal clustering. A range of slip behavior accommodates postseismic deformation including earthquakes, slow slip events, and earthquake swarms. Models of plate coupling and the consistency of earthquake clustering and slip behavior through multiple seismic cycles reveal a segmented subduction zone primarily controlled by subducted seafloor topography, accreted terranes, and inherited structure. The 2016 Pedernales mainshock triggered moderate to strong earthquakes (5 ≤ M ≤ 7) and earthquake swarms north of the mainshock rupture close to the epicenter of the 1906 Mw 8.8 earthquake and in the segment of the subduction zone that ruptured in 1958 in a Mw 7.7 earthquake.

Journal of Geophysical Research: Earth Surface↗
USGS2020

Submarine landslide kinematics derived from high-resolution imaging in Port Valdez, Alaska

Submarine landslides caused by strong ground shaking during the M9.2 1964 Great Alaska earthquake generated a tsunami that destroyed much of the old town of Valdez, Alaska, and was responsible for 32 deaths at that location. We explore structural details of the 1964 landslide deposit, as well as landslide deposits from earlier events, in order to characterize kinematics of the landslide process. We present a new high‐resolution seismic reflection data set that images the 1964 landslide deposit and six pre‐1964 deposits with great detail. These deposits are represented by thick packages (~7–23 m) of debris within >500 m of fjord sedimentation above basement. Internal slide structures are associated with distinctive landslide failure mechanisms, including detailed erosional and depositional features and structures resolved within both landslide blocks and distal debris flow layers. Based on comparisons of deposit volume from subbottom structure and differenced bathymetry, we refine prior interpretations of the source of failed material. New data show evidence for basal erosion and reworking of fjord‐floor sedimentation. Additionally, material comprising the 1964 landslide appears to have been translated and deformed by lateral thrusting, rather than having been sourced entirely from upslope evacuation zones. Taking into account these complexities in depositional patterns, we show variations in slide size through Holocene time and relate the history of landslides to the paleoseismic record. Collectively, these new observations demonstrate that Port Valdez has a repeated history of large submarine landslides, which are likely associated with large megathrust earthquakes.

Alaska↗
USGS2019

Introduction to special issue on gas hydrate in porous media: Linking laboratory and field‐scale phenomena

The proliferation of drilling expeditions focused on characterizing natural gas hydrate as a potential energy resource has spawned widespread interest in gas hydrate reservoir properties and associated porous media phenomena. Between 2017 and 2019, a Special Section of this journal compiled contributed papers elucidating interactions between gas hydrate and sediment based on laboratory, numerical modeling, and field studies. Motivated mostly by field observations in the northern Gulf of Mexico and offshore Japan, several papers focus on the mechanisms for gas hydrate formation and accumulation, particularly with vapor phase gas, not dissolved gas, as the precursor to hydrate. These studies rely on numerical modeling or laboratory experiments using sediment packs or benchtop micromodels. A second focus of the Special Section is the role of fines in inhibiting production of gas from methane hydrate, controlling the distribution of hydrate at a pore scale, and influencing the bulk behavior of seafloor sediments. Other papers fill knowledge gaps related to the physical properties of hydrate-bearing sediments and advance new approaches in coupled thermal-mechanical modeling of these sediments during hydrate dissociation. Finally, one study addresses the long-standing question about the fate of methane hydrate at the molecular level when CO2 is injected into natural reservoirs under hydrate-forming conditions.

Journal of Geophysical Research B: Solid Earth↗
USGS2020

Magma intrusion and volatile ascent beneath Norris Geyser Basin, Yellowstone National Park

Recent activity has provided new insights into the causes of surface deformation in and around the Yellowstone Caldera, a topic that has been debated since the discovery of caldera floor uplift more than four decades ago. An episode of unusually rapid uplift (>15 cm/yr) centered near Norris Geyser Basin along the north caldera rim began in late 2013 and continued until a M w 4.9 earthquake on 30 March 2014; thereafter, uplift abruptly switched to subsidence. Uplift at rates of several centimeters per year resumed in 2016 and continued at least through the end of 2018. Modeling of Global Positioning System and interferometric synthetic aperture radar data suggests an evolving process of deep magma intrusion during 1996–2001 followed by volatile ascent and accumulation at shallow levels, perhaps as shallow as a few hundred meters depth. The depth of shallow volatile accumulation appears to have shallowed from the 2014 to the 2016 deformation episode, and frequent eruptions of Steamboat Geyser since March 2018 are likely a surface manifestation of this ongoing process. Hydrothermal explosion features are prominent in the Norris Geyser Basin area, and the apparent shallow nature of the volatile accumulation implies an increased risk of hydrothermal explosions.

Wyoming↗
USGS2019

Separating sea and slow slip signals on the seafloor

Seafloor pressure measurements hold promise for estimating vertical displacements from transient slow slip events on submarine faults. We assess the accuracy of pressure offset estimates that evolve over days to weeks and the confidence with which they may be attributed to tectonic deformation or to the ocean water column. One common approach to resolve this ambiguity assumes water column pressures vary insignificantly over the study region and are represented by stable reference site pressures. Assessing the validity of this assumption requires independent evidence. Correlations between pressures and colocated temperatures collected during the Hikurangi Ocean Bottom Investigation of Tremor and Slow Slip experiment suggest temperatures might provide a useful independent proxy for water column pressures. We compared offsets estimated using several methods, with temperature and other proxies. The use of a temperature proxy was unsuccessful, because seafloor temperatures did not track the seasonal signal that contributes significantly to seafloor pressure changes over the slow slip event period. Regardless of the estimation method, offsets varied within a few cm around some uncertain reference level. Commonly used statistical measures are shown not to be reliable indicators of offset accuracy since offsets contribute minimally to the total variance. Offsets estimated using identical methods but with seafloor pressures simulated using a regional ocean model were larger than those derived from the data but had a similar pattern. Since the model simulates only water column processes, this suggests a significant fraction of the estimated pressure offsets are due to seasonal water column signal and are not of tectonic origin.

Journal of Geophysical Research- Solid Earth↗
USGS2019

The effect of brine on the electrical properties of methane hydrate

Gas hydrates possess lower electrical conductivity (inverse of resistivity) than either seawater or ice, but higher than clastic silts and sands, such that electromagnetic methods can be employed to help identify their natural formation in marine and permafrost environments. Controlled laboratory studies offer a means to isolate and quantify the effects of changing individual components within gas‐hydrate‐bearing systems, in turn yielding insight into the behavior of natural systems. Here we investigate the electrical properties of polycrystalline methane hydrate with ≥25% gas‐filled porosity and in mixture with brine. Initially, pure methane hydrate was synthesized from H 2 O ice and CH 4 gas while undergoing electrical impedance measurement, then partially dissociated to assess the effects of pure pore water accumulation on electrical conductivity. Methane hydrate + brine mixtures were then formed by either adding NaCl (0.25–2.5 wt %) to high‐purity ice or by using frozen seawater as a reactant. Conductivity was obtained from impedance measurements made in situ throughout synthesis while temperature cycled between +15 °C and −25 °C. Several possible conduction mechanisms were subsequently determined using equivalent circuit modeling. Samples with low NaCl concentration show a doping/impurity effect and a log linear conductivity response as a function of temperature. For higher salt content samples, conductivity increases exponentially with temperature and the log linear relationship no longer holds; instead, we observe phase changes within the samples that follow NaCl–H 2 O–CH 4 phase equilibrium predictions. Final samples were quenched in liquid nitrogen and imaged by cryogenic scanning electron microscopy (cryo‐SEM) to assess grain‐scale characteristics.

Journal of Geophysical Research↗
USGS2019

Coseismic slip and early afterslip of the M6.0 August 24, 2014 South Napa, California, earthquake

We employ strong motion seismograms and static offsets from the Global Positioning System, Interferometric Synthetic Aperture Radar, and other measurements in order to derive a coseismic slip and afterslip model of the M6.0 24 August 2014 South Napa earthquake. This earthquake ruptured an ∼13‐km‐long portion of the West Napa fault with predominantly right‐lateral strike slip. In the kinematic seismic slip inversions, we couple the coseismic slip and afterslip distributions by requiring both distributions to involve right‐lateral strike slip with positive amplitude, with the net static slip being the sum of the two. We consider several candidate fault geometries: a first involving two steeply east dipping fault planes that reach Earth's surface at the western surface trace (STW), where most surface rupture was observed, a second involving a steeply west dipping plane that also reaches Earth's surface at the STW, and a third involving a combination of two variably west dipping planes constrained to pass through the locus of postseismic seismicity located ∼1 km west of the STW. The data are best fit using the model of two east dipping fault planes, with coseismic slip up to ∼1.2 m on a dominant shallow asperity about 10 km north of the hypocenter and on deeper asperities on the southern part of the rupture. Afterslip up to 1 m is concentrated along the southern part of the rupture at depths 5 km, consistent with surface observations of afterslip. Seismic moments associated with coseismic slip and afterslip are 1.13×10 18 N m (Mw 6.00) and 3.64×10 17 N m, respectively.

California↗
USGS2020

Semiautomated estimates of directivity and related source properties of small to moderate southern California earthquakes using second seismic moments

We develop a semiautomated method for estimating with second seismic moments the directivity, rupture area, duration, and centroid velocity of earthquakes. The method is applied to 41 southern California earthquakes with magnitude in the range 3.5–5.2 and provides stable results for 28 events. Apparent source time functions (ASTFs) of P and S phases are derived using deconvolution with three stacked empirical Green's functions (seGf). The use of seGf suppresses nongeneric source effects, improves the focal mechanism correspondence to the analyzed earthquakes, and typically allows inclusion of 5 to 15 more ASTFs compared with analysis using a single eGf. Most analyzed earthquakes in the Trifurcation area of the San Jacinto Fault have directivities toward the northwest, while events around Cajon Pass and San Gabriel Mountain tend to propagate toward the southeast. These results are generally consistent with predictions for dynamic rupture on bimaterial interfaces associated with the imaged velocity contrasts in the area. The second moment inversions also provide constraints on the upper and lower bounds of rupture areas in our data set. Stress drops and uncertainties are estimated for elliptical ruptures using the derived characteristic rupture length and width. The semiautomated second moment method with seGfs can be used for routine application to moderate earthquakes in locations with good station coverage.

California↗
USGS2019

Offset channels may not accurately record strike-slip fault displacement: Evidence from landscape evolution models

Slip distribution, slip rate, and slip per event for strike‐slip faults are commonly determined by correlating offset stream channels—under the assumption that they record seismic slip—but offset channels are formed by the interplay of tectonic and geomorphic processes. To constrain offset channel development under known tectonic and geomorphic conditions, we use numerical landscape evolution simulations along a theoretical strike‐slip fault with uniform and steady uplift, erosion, and diffusion. We investigate the influence of four tectonic parameters (fault zone width, earthquake recurrence interval, variance of the recurrence interval, and total slip relative to channel spacing) on offset channel development through multiple earthquake cycles. Analysis of >3,000 automatically measured offsets from >135 simulations suggests ~30% variability in individual measurements, but modeled displacement is recovered by averaging multiple measurements. However, the average of multiple offset measurements systematically underestimates modeled slip except when the fault zone is less than ~5 m wide, total slip is less than channel spacing, and offsets are measured shortly after an earthquake. In these simulations, postearthquake landscape evolution widens the geomorphic expression of the fault zone and modifies apparent channel offsets. We distinguish this “geomorphic fault zone” from the tectonic fault zone (zone of coseismic distributed deformation). This study highlights the capability of landscape evolution models to explore a range of conditions not easily defined in natural examples and the importance of averaging multiple measurements. Our results verify that paleoseismic studies must consider how geomorphic change has modified offset markers and use caution interpreting slip histories with multiple earthquakes.

Journal of Geophysical Research Solid Earth↗
USGS2020

Are the stress drops of small earthquakes good predictors of the stress drops of moderate-to-large earthquakes?

The stress drops of small earthquakes often exhibit spatial patterns of variability. If moderate and large earthquakes follow the same spatial patterns, the stress drops of possible future damaging earthquakes could be better predicted by considering the stress drops of nearby small events. Better stress drop predictability could reduce ground-motion uncertainty in Probabilistic Seismic Hazard Assessment (PSHA) and Earthquake Early Warning (EEW). I find that for an internally consistent stress drop catalog of M1.8-3.1 events in southern California, the stress drops of the bigger earthquakes are predictable from the nearby smaller events. However, this catalog only weakly spatially correlates with another catalog of M3.0-5.8 earthquakes, and is spatially uncorrelated with five other stress drop catalogs of M≥3.4 earthquakes. For southern California events M5.5-7.5, stress drops compiled from the literature are weakly spatially correlated with the stress drops of the M1.8-3.1 events, although the correlations are not statistically significant. The lack of strong spatial correlation may be due to actual differences in the controlling factors of stress drop, for example dynamic weakening in moderate-to-large earthquakes. Alternatively, a stronger spatial correlation may exist that is obscured by methodological heterogeneity and large errors in the stress drop estimates. Either way, the stress drops of small earthquakes do not appear to be good predictors of the stress drops of nearby moderate-to-large earthquakes, at least for current techniques of stress drop estimation. If these results are representative, small-earthquake stress drops are not currently useful for substantially reducing uncertainty in PSHA and EEW.

California↗
USGS2020

Geodetic measurements of slow slip events southeast of Parkfield, CA

Tremor and low-frequency earthquakes are presumed to be indicative of surrounding slow, aseismic slip that is often below geodetic detection thresholds. This study uses data from borehole seismometers and long-baseline laser strainmeters to observe both the seismic and geodetic signatures of episodic tremor and slip on the Parkfield region of the San Andreas Fault near Cholame, CA. The observed occurrence rates of both the tremors and co-located families of low-frequency earthquakes are not steady but instead exhibit quasiperiodic bursts of increased activity. We show that these periods of elevated seismic activity correlate with statistically significant stacked strain signals consisting of 44 slow-slip events. Modeled individual slow-slip events and their total summed moment, which are constrained by seismic signals and stacked strain, respectively, indicate that the individual moment magnitudes of these events range from 4.6–5.2. We find that the measured geodetic signal likely precedes the seismic signal by several hours, consistent with the aseismic slip preceding and driving the observed seismic tremor activity. We confirm that strike-slip faults, in addition to subduction zones, are capable of producing episodic tremor and slip.

California↗
USGS2020

Methods for rapidly estimating velocity precision from GNSS time series in the presence of temporal correlation: A new method and comparison of existing methods

Time series of position estimates from Global Navigational Satellite System (GNSS) are used to measure the velocities of points on the surface of the Earth. Along with the velocity estimates, a measure of the precision is needed to assess the quality of the velocity measurement. Here, I evaluate rate uncertainties provided by four different methods that have been applied to geodetic time series. The most rigorous approach uses a data covariance that incorporates a variety of noise processes relevant to geodetic time series but is computationally demanding. Two other approaches are efficient algorithms and are used widely, but both can provide less rigorous estimates of the rate uncertainty. I propose and evaluate a fourth method, which provides estimates of rate uncertainty closer to the rigorous approach but is significantly less computationally demanding. I have evaluated all three methods against the more rigorous method using both simulations and time series from 190 GNSS sites. For data best characterized as having a flicker type noise process, one of the widely used methods overestimates the uncertainty by up to a factor of 2, while the other widely used method underestimates the uncertainty by less than a factor of 2. For a random-walk process, both methods underestimate the rate uncertainty by a factor of 3 to 5.

Journal of Geophysical Research↗