ALBERT

Description: The ionosphere is one of the largest error sources in GNSS (Global Navigation Satellite Systems) applications and can cause up to several meters of error in positioning. Especially for single-frequency users, who cannot correct the ionospheric delay, the information about the state of the ionosphere is mandatory. Dual- and multi-frequency GNSS users, on the other hand, can correct the ionospheric effect on their observations by linear combination. However, real-time applications such as autonomous driving or precision farming, require external high accuracy corrections for fast convergence. Mostly, this external information is given in terms of grids or coefficients of the vertical total electron content (VTEC). Globally distributed GNSS stations of different networks, such as the network of the International GNSS Services (IGS), provide a large number of multi-frequency observations which can be used to determine the state of the ionosphere. These data are used to generate Global Ionosphere Maps (GIM). Due to the inhomogeneous global distribution of GNSS real-time stations and especially due to the large data gaps over oceanic areas, the global VTEC models are usually limited in their spatial and spectral resolution. Most of the GIMs are mathematically based on globally defined radial basis functions, i.e., spherical harmonics (SH), with a maximum degree of 15 and provided with a spatial resolution of 2 . 5 ∘ × 5 ∘ in latitude and longitude, respectively. Regional GNSS networks, however, offer dense clusters of observations, which can be used to generate regional VTEC solutions with a higher spectral resolution. In this study, we introduce a two-step model (TSM) comprising a global model as the first step and a regional model as the second step. We apply polynomial and trigonometric B-spline functions to represent the global VTEC. Polynomial B-splines are used for modelling the finer structures of VTEC within selected regions, i.e., the densification areas. The TSM provides both, a global and a regional VTEC map at the same time. In order to study the performance, we apply the developed approach to hourly data of the global IGS network as well as the EUREF network of the European region for St. Patrick storm in March 2015. For the assessment of the generated maps, we use the dSTEC analysis and compare both maps with different global and regional products from the IGS Ionosphere Associated Analysis Centers, e.g., the global product from CODE (Berne, Switzerland) and from UPC (Barcelona, Spain), as well as the regional maps from ROB (Brussels, Belgium). The assessment shows a significant improvement of the regional VTEC representation in the form of the generated TSM maps. Among all other products used for comparison, the developed regional one is of the highest accuracy within the selected time span. Since the numerical tests are performed using hourly data with a latency of one to two hours, the presented procedure is seen as an intermediate step for the generation of high precision regional real-time corrections for modern applications.

Description: Radar altimeters have been measuring ocean significant wave height for more than three decades, with their data used to record the severity of storms, the mixing of surface waters and the potential threats to offshore structures and low-lying land, and to improve operational wave forecasting. Understanding climate change and long-term planning for enhanced storm and flooding hazards are imposing more stringent requirements on the robustness, precision, and accuracy of the estimates than have hitherto been needed. Taking advantage of novel retracking algorithms, particularly developed for the coastal zone, the present work aims at establishing an objective baseline processing chain for wave height retrieval that can be adapted to all satellite missions. In order to determine the best performing retracking algorithm for both LRM and DDA, an objective assessment is conducted in the framework of the ESA SSCCI project. All algorithms process the same L1 input dataset covering a time-period of up to two years. As a reference for validation, an ERA5-h wave model as well as an in-situ buoy dataset from the CMEMS INSTAC database are used. Five different metrics are evaluated: percentage and types of outliers, level of measurement noise, wave spectral variability, comparison against wave models, and comparison against in-situ data. The metrics are evaluated as a function of the distance to the nearest coast and the sea state. The results of the assessment show that all novel retracking algorithms perform better in the majority of the metrics than the baseline algorithms currently used for operational generation of the products. Nevertheless, the performance of the retrackers strongly differ depending on the coastal proximity and the sea state. Some retrackers show high correlations with the wave models and in-situ data but significantly under- or overestimate large-scale spectral variability. We propose a weighting scheme to select the most suitable retrackers for the SSCCI programme.

Description: The Kalman filter (KF) is widely applied in (ultra) rapid and (near) real-time ionosphere modeling to meet the demand on ionosphere products required in many applications extending from navigation and positioning to monitoring space weather events and naturals disasters. The requirement of a prior definition of the stochastic models attached to the measurements and the dynamic models of the KF is a drawback associated with its standard implementation since model uncertainties can exhibit temporal variations or the time span of a given test data set would not be large enough. Adaptive methods can mitigate these problems by tuning the stochastic model parameters during the filter run-time. Accordingly, one of the primary objectives of our study is to apply an adaptive KF based on variance component estimation to compute the global Vertical Total Electron Content (VTEC) of the ionosphere by assimilating different ionospheric GNSS measurements. Secondly, the derived VTEC representation is based on a series expansion in terms of compactly supported B-spline functions. We highlight the morphological similarity of the spatial distributions and the magnitudes between VTEC values and the corresponding estimated B-spline coefficients. This similarity allows for deducing physical interpretations from the coefficients. In this context, an empirical adaptive model to account for the dynamic model uncertainties, representing the temporal variations of VTEC errors, is developed in this work according to the structure of B-spline coefficients. For the validation, the differential slant total electron content (dSTEC) analysis and a comparison with Jason-2/3 altimetry data are performed. Assessments show that the quality of the VTEC products derived by the presented algorithm is in good agreement, or even more accurate, with the products provided by IGS ionosphere analysis centers within the selected periods in 2015 and 2017. Furthermore, we show that the presented approach can be applied to different ionosphere conditions ranging from very high to low solar activity without concerning time-variable model uncertainties, including measurement error and process noise of the KF because the associated covariance matrices are computed in a self-adaptive manner during run-time.

Description: Despite increasing interest in monitoring the global water cycle, the availability of in situ gauging and discharge time series is decreasing. However, this lack of ground data can partly be compensated for by using remote sensing techniques to observe river stages and discharge. In this paper, a new approach for estimating discharge by combining water levels from multi-mission satellite altimetry and surface area extents from optical imagery with physical flow equations at a single cross-section is presented and tested at the Lower Mississippi River. The datasets are combined by fitting a hypsometric curve, which is then used to derive the water level for each acquisition epoch of the long-term multi-spectral remote sensing missions. In this way, the chance of detecting water level extremes is increased and a bathymetry can be estimated from water surface extent observations. Below the minimum hypsometric water level, the river bed elevation is estimated using an empirical width-to-depth relationship in order to determine the final cross-sectional geometry. The required flow gradient is derived from the differences between virtual station elevations, which are computed in a least square adjustment from the height differences of all multi-mission satellite altimetry data that are close in time. Using the virtual station elevations, satellite altimetry data from multiple virtual stations and missions are combined to one long-term water level time series. All required parameters are estimated purely based on remote sensing data, without using any ground data or calibration. The validation at three gauging stations of the Lower Mississippi River shows large deviations primarily caused by the below average width of the predefined cross-sections. At 13 additional cross-sections situated in wide, uniform, and straight river sections nearby the gauges the Normalized Root Mean Square Error (NRMSE) varies between 10.95% and 28.43%. The Nash-Sutcliffe Efficiency (NSE) for these targets is in a range from 0.658 to 0.946.

Description: In this study, a new approach for estimating volume variations of lakes and reservoirs using water levels from satellite altimetry and surface areas from optical imagery is presented. Both input data sets, namely water level time series and surface area time series, are provided by the Database of Hydrological Time Series of Inland Waters (DAHITI), developed and maintained by the Deutsches Geodätisches Forschungsinsitut der Technischen Universität München (DGFI-TUM). The approach is divided into three parts. In the first part, a hypsometry model based on the new modified Strahler approach is computed by combining water levels and surface areas. The hypsometry model describes the dependency between water levels and surface areas of lakes and reservoirs. In the second part, a bathymetry between minimum and maximum surface area is computed. For this purpose, DAHITI land-water masks are stacked using water levels derived from the hypsometry model. Finally, water levels and surface areas are intersected with the bathymetry to estimate a time series of volume variations in relation to the minimum observed surface area. The results are validated with volume time series derived from in-situ water levels in combination with bathymetric surveys. In this study, 28 lakes and reservoirs located in Texas are investigated. The absolute volumes of the investigated lakes and reservoirs vary between 0.062 km 3 and 6.041 km 3 . The correlation coefficients of the resulting volume variation time series with validation data vary between 0.80 and 0.99. Overall, the relative errors with respect to volume variations vary between 2.8% and 14.9% with an average of 8.3% for all 28 investigated lakes and reservoirs. When comparing the resulting RMSE with absolute volumes, the absolute errors vary between 1.5% and 6.4% with an average of 3.1%. This study shows that volume variations can be calculated with a high accuracy which depends essentially on the quality of the used water levels and surface areas. In addition, this study provides a hypsometry model, high-resolution bathymetry and water level time series derived from surface areas based on the hypsometry model. All data sets are publicly available on the Database of Hydrological Time Series of Inland Waters.

Description: In polar regions, sea-ice hinders the precise observation of Sea Surface Heights (SSH) by satellite altimetry. In order to derive reliable heights for the openings within the ice, two steps have to be fulfilled: (1) the correct identification of water (e.g., in leads or polynias), a process known as lead classification; and (2) dedicated retracking algorithms to extract the ranges from the radar echoes. This study focuses on the first point and aims at identifying the best available lead classification method for Cryosat-2 SAR data. Four different altimeter lead classification methods are compared and assessed with respect to very high resolution airborne imagery. These methods are the maximum power classifier; multi-parameter classification method primarily based on pulse peakiness; multi-observation analysis of stack peakiness; and an unsupervised classification method. The unsupervised classification method with 25 clusters consistently performs best with an overall accuracy of 97%. Furthermore, this method does not require any knowledge of specific ice characteristics within the study area and is therefore the recommended lead detection algorithm for Cryosat-2 SAR in polar oceans.