The accuracy of state‐of‐the‐art global barotropic tide models is assessed using bottom pressure data, coastal tide gauges, satellite altimetry, various geodetic data on Antarctic ice shelves, and ...independent tracked satellite orbit perturbations. Tide models under review include empirical, purely hydrodynamic (“forward”), and assimilative dynamical, i.e., constrained by observations. Ten dominant tidal constituents in the diurnal, semidiurnal, and quarter‐diurnal bands are considered. Since the last major model comparison project in 1997, models have improved markedly, especially in shallow‐water regions and also in the deep ocean. The root‐sum‐square differences between tide observations and the best models for eight major constituents are approximately 0.9, 5.0, and 6.5 cm for pelagic, shelf, and coastal conditions, respectively. Large intermodel discrepancies occur in high latitudes, but testing in those regions is impeded by the paucity of high‐quality in situ tide records. Long‐wavelength components of models tested by analyzing satellite laser ranging measurements suggest that several models are comparably accurate for use in precise orbit determination, but analyses of GRACE intersatellite ranging data show that all models are still imperfect on basin and subbasin scales, especially near Antarctica. For the M2 constituent, errors in purely hydrodynamic models are now almost comparable to the 1980‐era Schwiderski empirical solution, indicating marked advancement in dynamical modeling. Assessing model accuracy using tidal currents remains problematic owing to uncertainties in in situ current meter estimates and the inability to isolate the barotropic mode. Velocity tests against both acoustic tomography and current meters do confirm that assimilative models perform better than purely hydrodynamic models.
Key Points
Tide model accuracy assessmentImproved accuraciesTidal current estimates
Interplate megathrust earthquakes have inflicted catastrophic damage on human society. Such an earthquake is predicted to occur in the near future along the Nankai Trough off southwestern Japan--an ...economically active and densely populated area in which megathrust earthquakes have already occurred. Megathrust earthquakes are the result of a plate-subduction mechanism and occur at slip-deficit regions (also known as 'coupling' regions), where friction prevents plates from slipping against each other and the accumulated energy is eventually released forcefully. Many studies have attempted to capture distributions of slip-deficit rates (SDRs) in order to predict earthquakes. However, these studies could not obtain a complete view of the earthquake source region, because they had no seafloor geodetic data. The Hydrographic and Oceanographic Department of the Japan Coast Guard (JHOD) has been developing a precise and sustainable seafloor geodetic observation network in this subduction zone to obtain information related to offshore SDRs. Here, we present seafloor geodetic observation data and an offshore interplate SDR-distribution model. Our data suggest that most offshore regions in this subduction zone have positive SDRs. Specifically, our observations indicate previously unknown regions of high SDR that will be important for tsunami disaster mitigation, and regions of low SDR that are consistent with distributions of shallow slow earthquakes and subducting seamounts. This is the first direct evidence that coupling conditions might be related to these seismological and geological phenomena. Our findings provide information for inferring megathrust earthquake scenarios and interpreting research on the Nankai Trough subduction zone.
In a context of global change and increasing anthropic pressure on the environment, monitoring the Earth system and its evolution has become one of the key missions of geosciences. Geodesy is the ...geoscience that measures the geometric shape of the Earth, its orientation in space, and gravity field. Time‐variable gravity, because of its high accuracy, can be used to build an enhanced picture and understanding of the changing Earth. Ground‐based gravimetry can determine the change in gravity related to the Earth rotation fluctuation, to celestial body and Earth attractions, to the mass in the direct vicinity of the instruments, and to vertical displacement of the instrument itself on the ground. In this paper, we review the geophysical questions that can be addressed by ground gravimeters used to monitor time‐variable gravity. This is done in relation to the instrumental characteristics, noise sources, and good practices. We also discuss the next challenges to be met by ground gravimetry, the place that terrestrial gravimetry should hold in the Earth observation system, and perspectives and recommendations about the future of ground gravity instrumentation.
Plain Language Summary
In a context of global change and increased human vulnerability to terrestrial hazard, monitoring the Earth system is one of the key challenges of geoscience. In particular, terrestrial gravimetry, with its precision at the level of one part of a billion, allows the monitoring of many phenomena, from water resource availability to volcanic activity. This paper reviews the technique, its advantages and limitations, how it has been used in the Earth monitoring, and the next challenges to be met by ground gravimetry.
Key Points
Review of the geophysical phenomena that can be addressed by terrestrial time‐variable gravity
The next challenges of terrestrial gravimetry are discussed
Terrestrial gravimetry remains central in the Earth observation system
The Geodesy Advancing Geosciences and EarthScope (GAGE) Facility Global Positioning System (GPS) Data Analysis Centers produce position time series, velocities, and other parameters for approximately ...2000 continuously operating GPS receivers spanning a quadrant of Earth's surface encompassing the high Arctic, North America, and Caribbean. The purpose of this review is to document the methodology for generating station positions and their evolution over time and to describe the requisite trade‐offs involved with combination of results. GAGE GPS analysis involves formal merging within a Kalman filter of two independent, loosely constrained solutions: one is based on precise point positioning produced with the GIPSY/OASIS software at Central Washington University and the other is a network solution based on phase and range double‐differencing produced with the GAMIT software at New Mexico Institute of Mining and Technology. The primary products generated are the position time series that show motions relative to a North America reference frame and secular motions of the stations represented in the velocity field. The position time series themselves contain a multitude of signals in addition to the secular motions. Coseismic and postseismic signals, seasonal signals from hydrology, and transient events, some understood and others not yet fully explained, are all evident in the time series and ready for further analysis and interpretation. We explore the impact of analysis assumptions on the reference frame realization and on the final solutions, and we compare within the GAGE solutions and with others.
Plain Language Summary
We review the methods used to analyze Global Positioning System (GPS) data from an observatory of over 1000 GPS stations that measure the motions of the Earth's surface and changes in the atmosphere‐‐‐the Plate Boundary Observatory. In these analyses, positions of GPS stations can be determined to within a few millimeters, and we observe stations moving steadily at speeds up to 60 mm/yr. These motions are largest at the boundary between the North American and Pacific tectonic plates. Other motions are measured across the rest of North America. We see changes from long term average motions due to plate tectonics, plus other changes due to earthquakes, ground water variations, volcanic activity, atmospheric changes, soli changes and other processes. We describe access to the data and the products generated by this observatory.
Key Points
Analysis methods for processing large GPS networks are described, and interpretation and access to products are discussed
Detailed comparisons of analysis results from PBO and other GPS processing groups are given
The estimation of scale changes (a common practice) has large impacts on vertical motion estimates
Geodetic analysis of radio tracking measurements of the MErcury Surface, Space ENvironment, GEochemistry, and Ranging spacecraft while in orbit about Mercury has yielded new estimates for the ...planet's gravity field, tidal Love number, and pole coordinates. The derived right ascension (α = 281.0082° ± 0.0009°; all uncertainties are 3 standard deviations) and declination (δ = 61.4164° ± 0.0003°) of the spin pole place Mercury in the Cassini state. Confirmation of the equilibrium state with an estimated mean (whole planet) obliquity ϵ of 1.968 ± 0.027 arcmin enables the confident determination of the planet's normalized polar moment of inertia (0.333 ± 0.005), which indicates a high degree of internal differentiation. Internal structure models generated by a Markov Chain Monte Carlo process and consistent with the geodetic constraints possess a solid inner core with a radius (ric) between 0.3 and 0.7 that of the outer core (roc).
Key Points
New solution for Mercury's gravity field provides crucial information on the planet's tidal response and orientation
New determination of Mercury's pole position fully satisfies the equilibrium Cassini state
New estimate of Mercury's polar moment of inertia supports the presence of a large solid inner core
Southern Africa is typically considered to belong to a single tectonic plate, Nubia, despite active faulting along the southwestern branch of the East African Rift System. We analyze regional Global ...Navigation Satellite System (GNSS) measurements, and find that the “San” microplate, situated south of the southwestern branch of the East African Rift, is statistically distinct from Nubia, with 0.4–0.7 mm/yr of extension across the boundary. Adding nine new campaign GNSS sites, we show that the extension rate across the southern Malawi Rift is 2.2 ± 0.3 mm/yr, with 75% of the relative velocity occurring over 890 km, despite the surface expression of faulting being <150 km wide. Thus, for the first time, we use geodetic measurements to describe the accommodation of strain in broad zones between Archean cratons in southern Africa's thick continental lithosphere.
Plain Language Summary
The breakup of continental plates is a fundamental part of plate tectonics, but little is known about how the plates start to stretch. In southern Africa, which is slowly extending, we use Global Positioning System (GPS) instruments to precisely measure the motion of the plates. Statistical tests show that southern Africa, which had previously been modeled as a single rigid plate called Nubia, can be better modeled as two separate plates, by the addition of the San microplate in southern Africa. We find that the plate‐boundary in southern Malawi is extending at ∼2 mm per year. However, the region that is actively extending is 890 km wide, which is wider than the 150 km wide region where we observe earthquakes and faults. Our new geodetic constraints suggest that in southern Africa, the continental plates are stretching apart at weaknesses around multiple microplates of old, thick crust, rather than along one single plate boundary.
Key Points
Based on geodetic data, the San microplate in southern Africa is distinct from the Nubian Plate with a boundary that follows active rifts
New Global Navigation Satellite System data show the southern Malawi Rift has an extension rate of 2.2 ± 0.3 mm/yr with 75% of relative velocity in a region 890 km wide
In thick continental lithosphere, strain accumulates over broad zones but seismicity concentrates in narrower zones at the edge of cratons
This article investigates the usage of terrestrial laser scanner (TLS) point clouds for monitoring the gradual movements of soil masses due to freeze–thaw activity and water saturation, commonly ...referred to as solifluction. Solifluction is a geomorphic process which is characteristic for hillslopes in (high-)mountain areas, primarily alpine periglacial areas and the arctic. The movement can reach millimetre-to-centimetre per year velocities, remaining well below the typical displacement mangitudes of other frequently monitored natural objects, such as landslides and glaciers. Hence, a better understanding of solifluction processes requires increased spatial and temporal resolution with relatively high measurement accuracy. To that end, we developed a workflow for TLS point cloud processing, providing a 3D vector field that can capture soil mass displacement due to solifluction with high fidelity. This is based on the common image-processing techniques of feature detection and tracking. The developed workflow is tested on a study area placed in Hohe Tauern range of the Austrian Alps with a prominent assemblage of solifluction lobes. The derived displacements were compared with the established geomonitoring approach with total station and signalized markers and point cloud deformation monitoring approaches. The comparison indicated that the achieved results were in the same accuracy range as the established methods, with an advantage of notably higher spatial resolution. This improvement allowed for new insights considering the solifluction processes.
Despite the importance of viscoelasticity in the evolution of crustal stress/strain being widely recognized, the interpretation of interseismic geodetic measurements for assessing earthquake ...potential is still based overwhelmingly on elastic models. The reasons for this disparity include conflating deformation rates with deformation itself and the lack of a succinct representation of the seismic readiness of a locked fault in a viscoelastic Earth. Using a classical viscoelastic model for strike‐slip faults, we reiterate the commonly overlooked message that, if the recurrence interval is long, most of the strain energy for the next earthquake accrues early in the cycle, and low strain rates later in the cycle by no means indicate diminished rupture potential. Fault stress stays near failure for much of the late interseismic period which may explain why slow slip‐rate faults have more variable recurrence intervals than fast slip‐rate faults. We propose to use displacement deficit instead of slip deficit to represent seismic readiness.
Plain Language Summary
Modern satellite measurements can reveal how quickly faults are being loaded by tectonic plate motions, and seismic hazard models use these loading rates as proxy for the likelihood of a pending earthquake. However, because of the partially fluid‐like behavior of Earth's interior, these loading rates have actually evolved with time since the last rupture. For faults with long intervals between successive earthquakes, these rates slow down substantially as the next event draws near. We, therefore, caution that slow rates of loading should not be assumed to reflect limited earthquake potential.
Key Points
Because of viscoelasticity, faults with long recurrence intervals accrue most of their elastic strain early in the interseismic period
Strain rates should not be conflated with stored strain, and slow geodetic deformation rates do not imply limited earthquake potential
For strike‐slip faults, “Relative Displacement Deficit” is a better measure of the earthquake readiness of a fault than “slip deficit”