ALBERT

Description: Precise modeling of tropospheric delay and weighted mean temperature (Tm) is critical for Global Navigation Satellite System (GNSS) positioning and meteorology. However, the model data in previous models cover a limited time span, which limits the accuracy of these models. Besides, the vertical variations of tropospheric delay and Tm are not perfectly modeled in previous studies, which affects the performance of height corrections. In this study, we used the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim reanalysis from 1979 to 2017 to build a new empirical model. We first carefully modeled the lapse rates of tropospheric delay and Tm. Then we considered the temporal variations by linear trends, annual, and semi-annual variations and the spatial variations by grids. This new model can provide zenith hydrostatic delay (ZHD), zenith wet delay (ZWD), and Tm worldwide with a spatial resolution of 1° × 1°. We used the ECMWF ERA-Interim data and the radiosonde data in 2018 to validate this new model in comparison with the canonical GPT2w model. The results show that the new model has higher accuracies than the GPT2w model in all parameters. Particularly, this new model largely improves the accuracy in estimating ZHD and Tm at high-altitude (relative to the grid point height) regions.

Description: With the availability to high-accuracy a priori zenith wet delay (ZWD) data, the positioning efficiency of the precise point positioning (PPP) processing can be effectively improved, including accelerating the convergence time and improving the positioning precision, in ground-based Global Navigation Satellite System (GNSS) technology. Considering the limitations existing in the state-of-the-art ZWD models, this paper established and evaluated a new in-situ meteorological observation-based grid model for estimating ZWD named GridZWD using the radiosonde data and the European Centre for Medium-Range Weather Forecasts (ECWMF) data. The results show that ZWD has a strong correlation with the meteorological parameter water vapor pressure in continental and high-latitude regions. The root of mean square error (RMS) of 24.6 mm and 36.0 mm are achievable by the GridZWD model when evaluated with the ECWMF data and the radiosonde data, respectively. An accuracy improvement of approximately 10%~30% compared with the state-of-the-art models (e.g., the Saastamoinen, Hopfield and GPT2w models) can be found for the new built model.

Description: This study uses the observed vertical displacements of Global Positioning System (GPS) time series obtained from the Crustal Movement Observation Network of China (CMONOC) with careful pre- and post-processing to estimate the seasonal crustal deformation in response to the hydrological loading in lower three-rivers headwater region of southwest China, followed by inferring the annual EWH changes through geodetic inversion methods. The Helmert Variance Component Estimation (HVCE) and the Minimum Mean Square Error (MMSE) criterion were successfully employed. The GPS inferred EWH changes agree well qualitatively with the Gravity Recovery and Climate Experiment (GRACE)-inferred and the Global Land Data Assimilation System (GLDAS)-inferred EWH changes, with a discrepancy of 3.2–3.9 cm and 4.8–5.2 cm, respectively. In the research areas, the EWH changes in the Lancang basin is larger than in the other regions, with a maximum of 21.8–24.7 cm and a minimum of 3.1–6.9 cm.

Description: The tropospheric delay is one major error source affecting the precise positioning provided by the global navigation satellite system (GNSS). This error occurs because the GNSS signals are refracted while travelling through the troposphere layer. Nowadays, various types of model can produce the tropospheric delay. Among them, the globally distributed GNSS permanent stations can resolve the tropospheric delay with the highest accuracy and the best continuity. Meteorological models, such as the Saastamoinen model, provide formulae to calculate temperature, pressure, water vapor pressure and subsequently the tropospheric delay. Some grid-based empirical tropospheric delay models directly provide tropospheric parameters at a global scale and in real time without any auxiliary information. However, the spatial resolution of the GNSS tropospheric delay is not sufficient, and the accuracy of the meteorological and empirical models is relatively poor. With the rapid development of satellite navigation systems around the globe, the demand for real-time high-precision GNSS positioning services has been growing dramatically, requiring real-time and high-accuracy troposphere models as a critical prerequisite. Therefore, this paper proposes a multi-source real-time local tropospheric delay model that uses polynomial fitting of ground-based GNSS observations, meteorological data, and empirical GPT2w models. The results show that the accuracy in the zenith tropospheric delay (ZTD) of the proposed tropospheric delay model has been verified with a RMS (root mean square) of 1.48 cm in active troposphere conditions, and 1.45 cm in stable troposphere conditions, which is significantly better than the conventional tropospheric GPT2w and Saastamoinen models.

Description: Remote Sensing, Vol. 10, Pages 1241: Troposphere Water Vapour Tomography: A Horizontal Parameterised Approach Remote Sensing doi: 10.3390/rs10081241 Authors: Qingzhi Zhao Yibin Yao Wanqiang Yao Global Navigation Satellite System (GNSS) troposphere tomography has become one of the most cost-effective means to obtain three-dimensional (3-d) image of the tropospheric water vapour field. Traditional methods divide the tomography area into a number of 3-d voxels and assume that the water vapour density at any voxel is a constant during the given period. However, such behaviour breaks the spatial continuity of water vapour density in a horizontal direction and the number of unknown parameters needing to be estimated is very large. This is the focus of the paper, which tries to reconstruct the water vapor field using the tomographic technique without imposing empirical horizontal and vertical constraints. The proposed approach introduces the layered functional model in each layer vertically and only an a priori constraint is imposed for the water vapor information at the location of the radiosonde station. The elevation angle mask of 30° is determined according to the distribution of intersections between the satellite rays and different layers, which avoids the impact of ray bending and the error in slant water vapor (SWV) at low elevation angles on the tomographic result. Additionally, an optimal weighting strategy is applied to the established tomographic model to obtain a reasonable result. The tomographic experiment is performed using Global Positioning System (GPS) data of 12 receivers derived from the Satellite Positioning Reference Station Network (SatRef) in Hong Kong. The quality of the established tomographic model is validated under different weather conditions and compared with the conventional tomography method using 31-day data, respectively. The numerical result shows that the proposed method is applicable and superior to the traditional one. Comparisons of integrated water vapour (IWV) of the proposed method with that derived from radiosonde and European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim data show that the root mean square (RMS)/Bias of their differences are 3.2/−0.8 mm and 3.3/−1.7 mm, respectively, while the values of traditional method are 5.1/−3.9 mm and 6.3/−5.9 mm, respectively. Furthermore, the water vapour density profiles are also compared with radiosonde and ECMWF data, and the values of RMS/Bias error for the proposed method are 0.88/0.06 g/m3 and 0.92/−0.08 g/m3, respectively, while the values of the traditional method are 1.33/0.38 g/m3 and 1.59/0.40 g/m3, respectively.

Description: The prevalence of real-time, low-cost, single-frequency, decimeter-level positioning has increased with the development of global navigation satellite systems (GNSSs). Ionospheric delay accounts for most errors in real-time single-frequency GNSS positioning. To eliminate ionospheric interference in real-time single-frequency precise point positioning (RT-SF-PPP), global ionospheric vertical total electron content (VTEC) product is designed in the next stage of the International GNSS Service (IGS) real-time service (RTS). In this study, real-time generation of a global ionospheric map (GIM) based on IGS RTS is proposed and assessed. There are three crucial steps in the process of generating a real-time global ionospheric map (RTGIM): estimating station differential code bias (DCB) using the precise point positioning (PPP) method, deriving slant total electron content (STEC) from PPP with raw observations, and modeling global vertical total electron content (VTEC). Experiments were carried out to validate the algorithm’s effectiveness. First, one month’s data from 16 globally distributed IGS stations were used to validate the performance of DCB estimation with the PPP method. Second, 30 IGS stations were used to verify the accuracy of static PPP with raw observations. Third, the modeling of residuals was assessed in high and quiet ionospheric activity periods. Afterwards, the quality of RTGIM products was assessed from two aspects: (1) comparison with the Center for Orbit Determination in Europe (CODE) global ionospheric map (GIM) products and (2) determination of the performance of RT-SF-PPP with the RTGIM. Experimental results show that DCB estimation using the PPP method can realize an average accuracy of 0.2 ns; static PPP with raw observations can achieve an accuracy of 0.7, 1.2, and 2.1 cm in the north, east, and up components, respectively. The average standard deviations (STDs) of the model residuals are 2.07 and 2.17 TEC units (TECU) for moderate and high ionospheric activity periods. Moreover, the average root-mean-square (RMS) error of RTGIM products is 2.4 TECU for the one-month moderate ionospheric period. Nevertheless, for the high ionospheric period, the RMS is greater than the RMS in the moderate period. A sub-meter-level horizontal accuracy and meter-level vertical accuracy can be achieved when the RTGIM is employed in RT-SF-PPP.

Description: With the rapid development of China’s economy and industry, characterizing the spatial and temporal changes of aerosols in China has attracted widespread attention from researchers. The national-level urban agglomerations are the most concentrated areas of China’s economic, population and resource. Studying the spatial and temporal changes of aerosol optical depth (AOD) in these regions has practical guiding significance for effective monitoring of atmospheric particulate pollution. This paper analyzed the spatial and temporal variations of AOD in China’s urban agglomerations during 2001–2017 by using Terra Moderate resolution Imaging Spectroradiometer (MODIS) Collection 6.1 (C6.1) Level 2 aerosol products (MOD04_L2). Five national-level urban agglomerations were chosen: Yangtze River Delta (YRD), Pearl River Delta (PRD), Beijing-Tianjin-Hebei (BTH), Yangtze River Middle-Reach (YRMR) and Cheng-Yu (CY). We analyzed the change patterns of AOD in different urban agglomerations at multi-time scales and built a time series decomposition model to mine the long-term trend, seasonal variation and abnormal change information of AOD time series. The result indicated that averaged AOD values in the five urban agglomerations were basically increased first and then decreased at the annual time scale during 2001–2017. The averaged AOD showed strong seasonal differences and AOD values in spring and summer were typically higher than those in autumn and winter. At the monthly time scale, the AOD typically varied from low in cold months to high in warm months and then decreased during the rainy periods. Time series decompositions revealed that a notable transition around 2007–2008 dominated the long-term overall trend over the five selected urban agglomerations and an initial upward tendency followed by a downward tendency was observed during 2001–2017. This study can be utilized to provide decision-making basis for atmospheric environmental governance and future development of urban agglomerations.

Description: Remote Sensing, Vol. 9, Pages 886: An Improved Tomography Approach Based on Adaptive Smoothing and Ground Meteorological Observations Remote Sensing doi: 10.3390/rs9090886 Authors: Bao Zhang Qingbiao Fan Yibin Yao Caijun Xu Xingxing Li Using the Global Navigation Satellite System (GNSS) to sense three-dimensional water vapor (WV) has been intensively investigated. However, this technique still heavily relies on the a priori information. In this study, we propose an improved tomography approach based on adaptive Laplacian smoothing (ALS) and ground meteorological observations. By using the proposed approach, the troposphere tomography is less dependent on a priori information and the ALS constraints match better with the actual situation than the constant constraints. Tomography experiments in Hong Kong during a heavy rainy period and a rainless period show that the ALS method gets superior results compared with the constant Laplacian smoothing (CLS) method. By validation with radiosonde and European Centre for Medium-Range Weather Forecasts (ECMWF) data, we found that the introduction of ground meteorological observations into tomography can solve the perennial problem of resolving the wet refractivity in the lower troposphere and thus significantly improve the tomography results. However, bad data quality and incompatibility of the ground meteorological observations may introduce errors into tomography results.

Description: Remote Sensing, Vol. 10, Pages 621: Plasmaspheric Electron Content Inferred from Residuals between GNSS-Derived and TOPEX/JASON Vertical TEC Data Remote Sensing doi: 10.3390/rs10040621 Authors: Lei Liu Yibin Yao Jian Kong Lulu Shan The plasmasphere, which is located above the ionosphere, is a significant component of Earth’s atmosphere, and the plasmasphere electron content (PEC) distribution is determined by different physical mechanisms to those of the ionosphere electron content (IEC). However, the observation for the PEC is very limited. In this study, we introduced a methodology (called zero assumption method, which is based on the assumption that PEC can reach zero) to extract the PEC over TOPEX/JASON (T/J) and global navigation satellite system (GNSS) overlapping areas. Results show that the daily systematic bias (T/J vertical TEC > GNSS-derived vertical TEC) for both low (2009) and high (2011) solar activity condition is consistent, and the systematic bias for JASON2 and JASON1 is different. We suggest that systematic biases predominantly arise from the sea state bias (SSB), especially the tracker bias. After removing the systematic bias, we extracted reliable PEC inferred from differences between GNSS-derived vertical TEC and T/J vertical TEC data. Finally, the characteristics of the plasmaspheric component distribution for different local times, latitudes, and seasons were investigated.

Description: The continuously operating Global Positioning System (GPS) sites mounted on bedrock around the coast of Greenland provide important geodetic datasets to quantify the solid Earth's response to historical and present-day ice mass variations. The presence of colored noise and irregular seasonal signals makes it difficult to detect transient changes in GPS time series. Here we apply the Multichannel Singular Spectral Analysis to the combination of GPS data and Gravity Recovery and Climate Experiment (GRACE) data so that we can identify and fully utilize the spatial correlations from these two independent datasets. Using the GPS and GRACE data near Upernavik Isstrøm in West Greenland as an example, we demonstrate that this method successfully detects two transient signals in ice mass variations during 2008 and 2014. Our forward modeling of loading displacements due to changes in surface mass balance (SMB) and ice dynamics suggests that the transient change starting in mid-2008 was due to the combined contributions from dynamically-induced mass loss and SMB. The transient change starting in mid-2011 was mainly due to ablation Specifically, the ice melted more in 2012 and less in 2013 with little contribution from anomalies in accumulation.