On 6 February 2023, a devastating doublet of earthquakes with magnitudes of Mw 7.8 and Mw 7.6 successively struck southeastern Turkey near the border of Syria. The earthquake sequence represents the strongest earthquakes in Turkey during the past 80 years and caused an extensive loss of life and property. In this study, we processed Sentinel-1 and GPS data to derive the complete surface displacement caused by the earthquake sequence. The surface displacements were adopted to invert for the fault geometry and coseismic slip distribution on the seismogenic faults of the earthquake sequence. The results indicate that the coseismic rupture of the Turkey earthquake sequence was dominated by left-lateral strike slips with a maximum slip of ~10 m on the East Anatolian Fault Zone (EAFZ) and the Sürgü fault (SF). Significant surface ruptures are recognized based on the geodetic inversion, which is consistent with the analysis of post-earthquake satellite images. The cumulative released moment of the two earthquakes reached 9.62 × 10²⁰ Nm, which corresponds to an event of Mw 7.95. Additionally, the interseismic fault slip rates and locking depths along the central and western segments of the EAFZ were estimated using the high-resolution long-term velocity field. The results reveal significant lateral variations of fault slip rates and locking depths along the central and western segments of the EAFZ. Generally, the estimated fault locking zone showed good spatial consistency with the coseismic fault rupture of the Mw 7.8 shock on the EAFZ. The static coulomb failure stress (CFS) change due to the Mw 7.8 earthquakes suggests that the subsequent Mw 7.6 event was certainly promoted by the Mw 7.8 shock. The stress transfers from the fault EAFZ to the fault SF were realized by unclamping the interface of the fault SF, which significantly reduces the effective normal stress on the fault plane. Large CFS increases in the western Puturge segment of the EAFZ, which was not ruptured in the 2020 Mw 6.8 and the 2023 Mw 7.8 earthquakes, highlight the future earthquake risk in this fault segment.

The knowledge of Arctic Sea ice coverage is of particular importance in studies of climate change. This study develops a new sea ice classification approach based on machine learning (ML) classifiers through analyzing spaceborne GNSS-R features derived from the TechDemoSat-1 (TDS-1) data collected over open water (OW), first-year ice (FYI), and multi-year ice (MYI). A total of eight features extracted from GNSS-R observables collected in five months are applied to classify OW, FYI, and MYI using the ML classifiers of random forest (RF) and support vector machine (SVM) in a two-step strategy. Firstly, randomly selected 30% of samples of the whole dataset are used as a training set to build classifiers for discriminating OW from sea ice. The performance is evaluated using the remaining 70% of samples through validating with the sea ice type from the Special Sensor Microwave Imager Sounder (SSMIS) data provided by the Ocean and Sea Ice Satellite Application Facility (OSISAF). The overall accuracy of RF and SVM classifiers are 98.83% and 98.60% respectively for distinguishing OW from sea ice. Then, samples of sea ice, including FYI and MYI, are randomly split into training and test dataset. The features of the training set are used as input variables to train the FYI-MYI classifiers, which achieve an overall accuracy of 84.82% and 71.71% respectively by RF and SVM classifiers. Finally, the features in every month are used as training and testing set in turn to cross-validate the performance of the proposed classifier. The results indicate the strong sensitivity of GNSS signals to sea ice types and the great potential of ML classifiers for GNSS-R applications.

The Sikkim–Bhutan seismic gap has witnessed a long earthquake quiescence since the 1714 M7.5~8.5 earthquake. The state of stress accumulation beneath the Sikkim–Bhutan Himalaya and its spatial correlation with seismicity remains unclear due to the lack of geodetic measurements and the low levels of seismic activity. We compile Global Positioning System (GPS) measurements in southern Tibet with the available velocities in the Sikkim–Bhutan Himalaya to reveal the characteristics of strain buildup on the Main Himalayan Thrust (MHT). We correct non-tectonic hydrological loading effects in a GPS time series to accurately determine the Three-Dimensional (3D) velocities of each continuous station. Extensive GPS measurements yield convergence rates of 16.2~18.5 mm/y across the Sikkim–Bhutan Himalaya, which is quite consistent with that observed elsewhere in the Himalaya. Based on a double-ramp structure of the MHT, a refined 3D coupling image is inverted using a dense network of GPS velocities. The result indicates significant along-strike variations of fault coupling beneath the Sikkim–Bhutan Himalaya. The locking width (coupling > 0.5) of western Bhutan reaches ~100 km, which is 30~40% wider than Sikkim and eastern Bhutan. An obvious embayment of decoupling zone near the border between Sikkim and western Bhutan is recognized, and coincides spatially with the rupture terminates of the 1934 Mw8.2 and the 1714 M7.5~8.5 earthquakes, indicating that the large megathrust earthquakes along the Sikkim–Bhutan Himalaya are largely segmented by the spatial variation of frictional properties on the MHT. Using a new compilation of seismic records in the Sikkim–Bhutan Himalaya, we analyze the spatial correlation between fault coupling and seismic activity. The result suggests that the seismicity in the Bhutan Himalaya is broadly distributed, instead of restricted in the lower edge of the interseismic locking zone. This implies that the seismic activity in the Bhutan Himalaya is not uniquely controlled by the stress accumulation at the downdip end of the locked portion of the MHT.

Two effective machine learning-aided sea ice monitoring methods are investigated using 42 months of spaceborne Global Navigation Satellite System-Reflectometry (GNSS-R) data collected by the TechDemoSat-1 (TDS-1). The two-dimensional delay waveforms with different Doppler spread characteristics are applied to extract six features, which are combined to monitor sea ice using the decision tree (DT) and random forest (RF) algorithms. Firstly, the feature sequences are used as input variables and sea ice concentration (SIC) data from the Advanced Microwave Space Radiometer-2 (AMSR-2) are applied as targeted output to train the sea ice monitoring model. Hereafter, the performance of the proposed method is evaluated through comparing with the sea ice edge (SIE) data from the Special Sensor Microwave Imager Sounder (SSMIS) data. The DT- and RF-based methods achieve an overall accuracy of 97.51% and 98.03%, respectively, in the Arctic region and 95.46% and 95.96%, respectively, in the Antarctic region. The DT- and RF-based methods achieve similar accuracies, while the Kappa coefficient of RF-based approach is slightly larger than that of the DT-based approach, which indicates that the RF-based method outperforms the DT-based method. The results show the potential of monitoring sea ice using machine learning-aided GNSS-R approaches.

As a widely used ground-based optical instrument, the LAI-2000 or LAI-2200 plant canopy analyzer (PCA) (Li-Cor, Inc., Lincoln, NE) is designed to measure the plant effective leaf area index (L_e) by measuring the canopy gap fraction at several limited or discrete view zenith angles (VZAs) (usually five VZAs: 7, 23, 38, 53, and 68°) based on Miller’s equation. Miller’s equation requires the probability of radiative transmission through the canopy to be measured over the hemisphere, i.e., VZAs in the range from 0 to 90°. However, the PCA view angle ranges are confined to several limited ranges or discrete sectors. The magnitude of the error produced by the discretization of VZAs in the leaf area index measurements remains difficult to determine. In this study, a theoretical deduction was first presented to definitely prove why the limited or discrete VZAs or ranges can affect the L_e measured with the PCA, and the specific error caused by the limited or discrete VZAs was described quantitatively. The results show that: (1) the weight coefficient of the last PCA ring is the main cause of the error; (2) the error is closely related to the leaf inclination angles (IAs)—the L_e measured with the PCA can be significantly overestimated for canopies with planophile IAs, whereas it can be underestimated for erectophile IAs; and (3) the error can be enhanced with the increment of the discrete degree of PCA rings or VZAs, such as using four or three PCA rings. Two corrections for the error are presented and validated in three crop canopies. Interestingly, although the leaf IA type cannot influence the L_e calculated by Miller’s equation in the hemispheric space, it affects the L_e measured with the PCA using the discrete form of Miller’s equation for several discrete VZAs.