We present a comprehensive review of the status and changes in glacier length (since the 1850s), area and mass (since the 1960s) along the Himalayan-Karakoram (HK) region and their climate-change ...context. A quantitative reliability classification of the field-based mass-balance series is developed. Glaciological mass balances agree better with remotely sensed balances when we make an objective, systematic exclusion of likely flawed mass-balance series. The Himalayan mean glaciological mass budget was similar to the global average until 2000, and likely less negative after 2000. Mass wastage in the Himalaya resulted in increasing debris cover, the growth of glacial lakes and possibly decreasing ice velocities. Geodetic measurements indicate nearly balanced mass budgets for Karakoram glaciers since the 1970s, consistent with the unchanged extent of supraglacial debris-cover. Himalayan glaciers seem to be sensitive to precipitation partly through the albedo feedback on the short-wave radiation balance. Melt contributions from HK glaciers should increase until 2050 and then decrease, though a wide range of present-day area and volume estimates propagates large uncertainties in the future runoff. This review reflects an increasing understanding of HK glaciers and highlights the remaining challenges.
New satellite missions (e.g., the European Space Agency's Sentinel‐1 constellation), advances in data downlinking, and rapid product generation now provide us with the ability to access ...space‐geodetic data within hours of their acquisition. To truly take advantage of this opportunity, we need to be able to interpret geodetic data in a prompt and robust manner. Here we present a Bayesian approach for the inversion of multiple geodetic data sets that allows a rapid characterization of posterior probability density functions (PDFs) of source model parameters. The inversion algorithm efficiently samples posterior PDFs through a Markov chain Monte Carlo method, incorporating the Metropolis‐Hastings algorithm, with automatic step size selection. We apply our approach to synthetic geodetic data simulating deformation of magmatic origin and demonstrate its ability to retrieve known source parameters. We also apply the inversion algorithm to interferometric synthetic aperture radar data measuring co‐seismic displacements for a thrust‐faulting earthquake (2015 Mw 6.4 Pishan earthquake, China) and retrieve optimal source parameters and associated uncertainties. Given its robustness and rapidity in estimating deformation source parameters and uncertainties, our Bayesian framework is capable of taking advantage of real‐time geodetic measurements. Thus, our approach can be applied to geodetic data to study magmatic, tectonic, and other geophysical processes, especially in rapid‐response operational settings (e.g., volcano observatories). Our algorithm is fully implemented in a MATLAB®‐based software package (Geodetic Bayesian Inversion Software) that we make freely available to the scientific community.
Key Points
We present a Bayesian approach for the inversion of geodetic data and demonstrate successful applications to synthetic and real data
Our approach allows rapid estimates of source parameters and uncertainties and is well suited for rapid‐response and operational settings
We have implemented our approach in a MATLAB®‐based software package (GBIS) that is made freely available to the scientific community
The present study elucidates velocity of the Indian plate from 2013 to 2016 with the help of six continuously operating permanent Global Positioning System (GPS) stations and tracks the crustal ...motion and direction directly by means of space-based geodetic measurements. The GPS derived velocity estimates were computed with respect to 2008 International Terrestrial Reference Frame (ITRF 08) and carried out almost all along the arc in the Indian Himalayan Region (IHR). The velocity vectors indicate that the Indian plate is moving towards the NE direction (average velocity of 46.95 ± 0.23 mm/yr) with higher movement towards the eastward direction (36.11 ± 0.17 mm/yr) compared to the northward direction (29.02 ± 0.16 mm/yr). Variations in the plate motion between IISC and permanent GPS stations suggest that presently, the convergence rate is about 8.06 ± 0.28 mm/yr and 5.71 ± 0.17 mm/yr for Higher Himalaya and Lesser Himalaya respectively. In addition, the deformation rate was also calculated on the basis of baseline shortening between IGS stations (IISC and LHAZ) and permanent GPS stations. The results suggest that the deformation is about 11.68 ± 1.32 mm/yr above the Main Central Thrust (MCT) and 6.74 ± 1.20 mm/yr above the Main Boundary Thrust (MBT) with respect to the Indian plate (IISC), while it is 20.60 ± 1.76 mm/yr above the MCT and 11.42 ± 1.21 mm/yr above the MBT with respect to Eurasian plate (LHAZ). The strain rate through GPS measurement reveals that the maximum strain is accumulated in the central part of the Himalaya between the MBT and MCT, indicating vulnerability of the central Indian Himalaya for future earthquakes.
Geodetic and topographic works determine the necessary information base and of major importance in the design and management of underground construction works (mining, hydrotechnics, roads, etc.). ...The efficiency and safety of these investments with special financial implications depend on their quality. In this context, the geodetic orientation of the underground topographic works was the main and permanent subject of study and analysis on which this scientific paper refers.
The European Space Agency (ESA) is preparing a satellite mission called GENESIS to be launched in 2027 as part of the FutureNAV program. GENESIS co-locates, for the first time, all four space ...geodetic techniques on one satellite platform. The main objectives of the mission are the realization of the International Terrestrial Reference Frames and the mitigation of biases in geodetic measurements; however, GENESIS will remarkably contribute to the determination of the geodetic parameters. The precise GENESIS orbits will be determined through satellite-to-satellite tracking, employing two GNSS antennas to observe GPS and Galileo satellites in both nadir and zenith directions. In this research, we show results from simulations of GENESIS and Galileo-like constellations with joint orbit and clock determination. We assess the orbit quality of GENESIS based on nadir-only, zenith-only, and combined nadir–zenith GNSS observations. The results prove that GENESIS and Galileo joint orbit and clock determination substantially improves Galileo orbits, satellite clocks, and even ground-based clocks of GNSS receivers tracking Galileo satellites. Although zenith and nadir GNSS antennas favor different orbital planes in terms of the number of collected observations, the mean results for each Galileo orbital plane are improved to a similar extent. The 3D orbit error of Galileo is improved from 27 mm (Galileo-only), 23 mm (Galileo + zenith), 16 mm (Galileo + nadir), to 14 mm (Galileo + zenith + nadir GENESIS observations), i.e., almost by a factor of two in the joint GENESIS + Galileo orbit and clock solutions.
Rapid shortening in convergent mountain belts is often accommodated by slip on faults at multiple levels in upper crust, but no geodetic observation of slip at multiple levels within hours of a ...moderate earthquake has been shown before. Here we show clear evidence of fault slip within a shallower thrust at 5–10 km depth in SW Taiwan triggered by the 2016 Mw 6.4 MeiNong earthquake at 15–20 km depth. We constrain the primary coseismic fault slip with kinematic modeling of seismic and geodetic measurements and constrain the triggered slip and fault geometry using synthetic aperture radar interferometry. The shallower thrust coincides with a proposed duplex located in a region of high fluid pressure and high interseismic uplift rate, and may be sensitive to stress perturbations. Our results imply that under tectonic conditions such as high‐background stress level and high fluid pressure, a moderate lower crustal earthquake can trigger faults at shallower depth.
Key Points
Shallower fault slip was triggered by median size deeper earthquake
Ascending and descending InSAR reveals detailed coseismic surface deformation and helps constrain fault slip at different depths
Evidence of triggered slip on a shallower duplex thrust structure and this deformation are inferred to be slow or aseismic
The detection of preslip, occurring hours to days before a large earthquake, using geodetic measurements has been a major focus in earthquake prediction research. A recent study claims to have ...detected a preseismic signal interpreted as accelerating slip near the hypocenter of the 2011 great Tohoku‐oki earthquake, starting approximately 2 hr before the mainshock. This claim is based on a stacking procedure using GNSS (Global Navigation Satellite System) data. However, a follow‐up study demonstrated that the signal disappeared when specific GNSS noise was corrected. Here we utilize tiltmeter records, independent on GNSS, to check whether the claimed preseismic signal is detected using a similar stacking procedure. Our results show no acceleration‐like deformation from 2 hr before the mainshock. This indicates that no precursory slip exceeded the noise level of the tilt data, and if any preslip occurred, it was less than 5.0 × 1018 Nm in seismic moment.
Plain Language Summary
The ability to detect large earthquakes before they occur would be invaluable for mitigating damage. Researchers have been searching for precursory signals in crustal deformation data, believing that if precursory slip (fault movement before a major earthquake) is large enough, it could be detected by GNSS or tiltmeters, enabling earthquake prediction. The 2011 Tohoku‐oki earthquake provided a valuable opportunity to study this phenomenon. Extensive data was collected, but previous studies found no significant pre‐earthquake signals on timescales of hours to days. A recent GNSS study claimed an acceleration‐like change 2 hr before the main shock, but another study using the same data disputed this. While reanalysis using the same GNSS data is important, it is also susceptible to the influence of the same noise sources. To avoid this, this study analyzed tiltmeter data, finding no evidence for the precursory deformation suggested earlier. This result indicates no significant precursory slip before the Tohoku‐oki earthquake that causes deformation larger than a noise level of the tiltmeter data. This result can provide constraints on the magnitude of the precursory slip immediately before the mainshock.
Key Points
Tiltmeter records are used to determine whether a preseismic signal of the 2011 great Tohoku‐oki earthquake is detected
No acceleration‐like tilt deformation from about 2 hr before the mainshock is recorded
An upper bound on the size of the preslip immediately before the mainshock can be estimated from the noise level of the observation data
Oblique convergence across the northern Qilian Shan is accommodated by sub‐parallel strike‐slip and thrust faults that ruptured simultaneously in the Mw 8 Gulang earthquake in 1927. We investigate ...the kinematics of fault loading in the northern Qilian Shan and provide insights into the conditions necessary for generating multi‐fault earthquakes. We perform Bayesian inversions for the geometry and creep rate on the fault network. We infer that all of the thrust faults are locked north of the Qilian‐Haiyuan strike‐slip fault and are accumulating elastic strain. Multi‐fault earthquakes may occurr in this fault system because the faults are simultaneously loaded by the same source of deformation and are linked together by locked fault segments. The interseismic velocity field alone can not contain the location or activity of individual faults visible in the geomorphology, therefore the short‐term geodetic measurements may not reliably indicate the long‐term behavior of the fault system.
Plain Language Summary
This study aims to understand the earthquake hazard in the northern Qilian region of China. We use measurements of ground deformation between earthquakes to infer how the faults are being loaded in the region. We find that the ground deformation can be explained by a simple model with a single, slowly creeping fault at depth that loads all of the overlying faults. Large earthquakes that were caused by slip on many different faults at the same time have occurred in this region before. We suggest these so‐called “multi‐fault” earthquakes may occur because all of the faults are being simultaneously loaded by the same source of stress.
Key Points
Oblique convergence in the northern Qilian Shan is accommodated by sub‐parallel thrust and strike‐slip faulting
The short‐term geodetic measurements do not constrain the thrust fault kinematics in the northern Qilian Shan over geological timescales
Multi‐fault earthquakes may be common in the region as all of the shallow thrust faults are linked together by locked fault segments
Unrest began in July 2021 at Askja volcano in the Northern Volcanic Zone (NVZ) of Iceland. Its most recent eruption, in 1961, was predominantly effusive and produced ∼0.1 km3 lava field. The last ...plinian eruption at Askja occurred in 1875. Geodetic measurements between 1983 and 2021 detail subsidence of Askja, decaying in an exponential manner. At the end of July 2021, inflation was detected at Askja volcano, from GNSS observations and Sentinel‐1 interferograms. The inflationary episode can be divided into two periods from the onset of inflation until September 2023. An initial period until 20 September 2021 when geodetic models suggest transfer of magma (or magmatic fluids) from within the shallowest part of the magmatic system (comprising an inflating and deflating source), potentially involving silicic magma. A following period when one source of pressure increase at shallow depth can explain the observations.
Plain Language Summary
Askja volcano, situated in the Northern Volcanic Zone in Iceland, has been quiet since its last eruption in 1961, with surface deformation measurements from 1983 to 2021 displaying a decaying subsidence signal within the Askja caldera. However, at the end of July 2021, the volcano began to inflate. This was detected on both GNSS and satellite observations. As of September 2023, ∼65 cm of uplift had been measured at GNSS station OLAC. Modeling of surface deformation measurements indicates that the inflation was triggered by upward migration of melt (or magmatic fluids).
Key Points
At the end of July 2021, Askja volcano began to inflate—detected on both GNSS and satellite observations, ending 1983–2021 subsidence
Geodetic modeling indicates upward migration of magma, feeding a magma body at an inferred depth of 2.5–3.1 km under the main Askja caldera
Start of unrest was associated with magma transfer within the upper part of the system, followed by possible additional influx from depth
The accuracy of ionospheric models estimated by ground‐based multiple global navigation satellite system ionospheric data over regions with sparse tracking stations is not ideal. To improve the ...accuracy of the estimated ionospheric model, different types of ionospheric data with different combinations were employed for previous studies. However, the ionospheric observational ranges for different types of ionospheric data are not the same. In this study, the accuracy of ionospheric maps generated by ground‐based ionospheric data (ground‐based strategy) and ground‐based ionospheric data combined with data provided by other geodetic measurements normalized by the single‐layer normalization method (multi‐source strategy) were studied. The results showed that the main differences between the ionospheric models estimated by the two strategies occur for data taken over the ocean, which mainly range from −1 to 0 total electron content unit (TECU). When assessed using Jason‐3 vertical total electron content data, the mean root mean square (RMS) value of the ionospheric model estimated by the multi‐source strategy was 5.03 TECU, which is approximately 15% smaller than that estimated by the ground‐based strategy. The maximum reduction in results using the multisource strategy was approximately 25% over different latitudes compared with that of the ground‐based strategy. Furthermore, the self‐consistency evaluation method was employed for evaluation. The results showed that the RMS of the ionospheric model estimated by the multi‐source strategy was 2.41 TECU, which is 3.60% better than that of the ground‐based strategy. The maximum reduction was 15% on different days.
Key Points
The single‐layer normalization method is employed to make multi‐source data have the same observational range compared with those of ground‐based global navigation satellite system data
The distribution of multi‐source data has better coverage, especially over the oceanic region
The main differences between the ionospheric models estimated by the ground‐based strategy and the multi‐source strategy occur over the ocean