Spaceborne interferometric synthetic aperture radar (InSAR) methodology has been widely successfully applied to measure urban surface micro slow subsidence. However, the accuracy is still limited by the spatial resolution of currently operating SAR systems and the lacking precision of geolocation of the respective scatters. In this context, high-precision urban models, as provided by the active laser point cloud methodology through light detection and ranging (LiDAR) techniques, can assist in improving the geolocation quality of InSAR-derived permanent scatters (PS) and provide the precise contour of buildings for hazard analysis. This paper proposes to integrate InSAR and LiDAR technologies for an improved detailed analysis of subsidence levels and a hazard assessment for buildings in the urban environment. By the use of LiDAR data, most building contours in the main subsidence area were extracted and SAR positioning of buildings via PS points was refined more precisely. The workflow for the proposed method includes the monitoring of land subsidence by the TS-InSAR technique, the geolocation improvement of InSAR-derived PS, and building contour extraction by LiDAR data. Furthermore, a reasonable hazard assessment system of land subsidence was developed. Significant vertical subsidence of −40 to 12 mm per year was detected by the analysis of multisensor SAR images. The land subsidence rates in the Shenzhen District obviously follow certain spatial patterns. Most stable areas are located in the middle and northeast of Shenzhen except for some areas in Houhai, the Qianhai Bay, and the Wankeyuncheng. An additional hazard assessment of land subsidence reveals that the subsidence of buildings is mainly caused by the construction of new buildings and some by underground activities. The research results of this paper can provide a useful synoptic reference for urban planning and help reducing land subsidence in Shenzhen.

An outbreak of Ulva prolifera poses a massive threat to coastal ecology in the Southern Yellow Sea, China (SYS). It is a necessity to extract its area and monitor its development accurately. At present, Ulva prolifera monitoring by remote sensing imagery is mostly based on a fixed threshold or artificial visual interpretation for threshold selection, which has large errors. In this paper, an adaptive threshold model based on Google Earth Engine (GEE) is proposed and applied to extract U. prolifera in the SYS. The model first applies the Floating Algae Index (FAI) or Normalized Difference Vegetation Index (NDVI) algorithm on the preprocessed remote sensing images and then uses the Canny Edge Filter and Otsu threshold segmentation algorithm to extract the threshold automatically. The model is applied to Landsat8/OLI and Sentinel-2/MSI images, and the confusion matrix and cross-sensor comparison are used to evaluate the accuracy and applicability of the model. The verification results show that the model extraction of U. prolifera based on the FAI algorithm has higher accuracy (R² = 0.99, RMSE = 5.64) and better robustness. However, when the average cloud cover is more than 70% in the image (based on the statistical results of multi-year cloud cover information), the model based on the NDVI algorithm has better applicability and can extract the algae distributed at the edge of the cloud. When the model uses the FAI algorithm, it is named FAI-COM (model based on FAI, the Canny Edge Filter, and Otsu thresholding). And when the model uses the NDVI algorithm, it is named NDVI-COM (model based on NDVI, the Canny Edge Filter, and Otsu thresholding). Therefore, the final extraction results are generated by supplementing NDVI-COM results on the basis of FAI-COM extraction results in this paper. The F1-score of U. prolifera extracted results is above 0.85. The spatiotemporal distribution of U. prolifera in the South Yellow Sea from 2016 to 2020 is obtained through the model calculation. Overall, the coverage area of U. prolifera shows a decreasing trend over the five years. It is found that the delay in recovery time of Porphyra yezoensis culture facilities in the Northern Jiangsu Shoal and the manual salvage and cleaning-up of U. prolifera in May are among the reasons for the smaller interannual scale of algae in 2017 and 2018.

Ground-based GNSS-R (global navigation satellite system reflectometry) can provide the absolute vertical distance from a GNSS antenna to the reflective surface of the ocean in a common height reference frame, given that vertical crustal motion at a GNSS station can be determined using direct GNSS signals. This technique offers the advantage of enabling ground-based sea level measurements to be more accurately determined compared with traditional tide gauges. Sea level changes can be retrieved from multipath effects on GNSS, which is caused by interference of the GNSS L-band microwave signals (directly from satellites) with reflections from the environment that occur before reaching the antenna. Most of the GNSS observation types, such as pseudo-range, carrier-phase and signal-to-noise ratio (SNR), suffer from this multipath effect. In this paper, sea level altimetry determinations are presented for the first time based on geometry-free linear combinations of the carrier phase at low elevation angles from a fixed global positioning system (GPS) station. The precision of the altimetry solutions are similar to those derived from GNSS SNR data. There are different types of observation and reflector height retrieval methods used in the data processing, and to analyze the performance of the different methods, five sea level determination strategies are adopted. The solutions from the five strategies are compared with tide gauge measurements near the GPS station, and the results show that sea level changes determined from GPS SNR and carrier phase combinations for the five strategies show good agreement (correlation coefficient of 0.97–0.98 and root-mean-square error values of <0.2 m).

The disturbances of the ionosphere caused by thunderstorms or lightning events in the troposphere have an impact on global navigation satellite system (GNSS) signals. Gravity waves (GWs) triggered by thunderstorms are one of the main factors that drive short-period Travelling Ionospheric Disturbances (TIDs). At mid-latitudes, ionospheric GWs can be detected by GNSS signals. However, at low latitudes, the multi-variability of the ionosphere leads to difficulties in identifying GWs induced by thunderstorms through GNSS data. Though disturbances of the ionosphere during low-latitude thunderstorms have been investigated, the explicit GW observation by GNSS and its propagation pattern are still unclear. In this paper, GWs with periods from 6 to 20 min are extracted from band-pass filtered GNSS carrier phase observations without cycle-slips, and 0.2–0.8 Total Electron Content Unit (TECU) magnitude perturbations are observed when the trajectories of ionospheric pierce points fall into the perturbed region. The propagation speed of 102.6–141.3 m/s and the direction of the propagation indicate that the GWs are propagating upward from a certain thunderstorm at lower atmosphere. The composite results of disturbance magnitude, period, and propagation velocity indicate that GWs initiated by thunderstorms and propagated from the troposphere to the ionosphere are observed by GNSS for the first time in the low-latitude region.

The Multi-constellation Global Navigation Satellite System (Multi-GNSS) has become the standard implementation of high accuracy positioning and navigation applications. It is well known that the noise of code and phase measurements depend on GNSS constellation. Then, Helmert variance component estimation (HVCE) is usually used to adjust the contributions of different GNSS constellations by determining their individual variances of unit weight. However, HVCE requires a heavy computation load. In this study, the HVCE posterior weighting was employed to carry out a kinematic relative Multi-GNSS positioning experiment with six short-baselines from day of year (DoY) 171 to 200 in 2019. As a result, the HVCE posterior weighting strategy improved Multi-GNSS positioning accuracy by 20.5%, 15.7% and 13.2% in east-north-up (ENU) components, compared to an elevation-dependent (ED) priori weighting strategy. We observed that the weight proportion of both code and phase observations for each GNSS constellation were consistent during the entire 30 days, which indicates that the weight proportions of both code and phase observations are stable over a long period of time. It was also found that the quality of a phase observation is almost equivalent in each baseline and GNSS constellation, whereas that of a code observation is different. In order to reduce the time consumption of the HVCE method without sacrificing positioning accuracy, the stable variances of unit weights of both phase and code observations obtained over 30 days were averaged and then frozen as a priori information in the positioning experiment. The result demonstrated similar ENU improvements of 20.0%, 14.1% and 11.1% with respect to the ED method but saving 88% of the computation time of the HCVE strategy. Our study concludes with the observations that the frozen variances of unit weight (FVUW) could be applied to the positioning experiment for the next 30 days, that is, from DoY 201 to 230 in 2019, improving the positioning ENU accuracy of the ED method by 18.1%, 13.2% and 10.6%, indicating the effectiveness of the FVUW.