Mixed-integer models arise in several geodetic problems, including precise positioning and remote sensing in Global Navigation Satellite Systems (GNSS), as well as deformation monitoring through Interferometric Synthetic Aperture Radar (InSAR) or fringe phase observations from Very Long Baseline Interferometry (VLBI). These problems generally involve two types of unknowns: integer ambiguities a ∈ Z^n and real-valued parameters b ∈ R^p, whose accuracy can be significantly improved by correctly resolving the ambiguities. However, in some cases, a large number of ambiguity components are involved and need to be correctly resolved; therefore the ambiguity resolution process becomes a bottleneck for the computations.

One research question is therefore how to effectively tackle the challenges of high-dimensional ambiguity resolution and its computational complexity, while ensuring a successful resolution of the ambiguities. Additionally, a second question arises regarding whether it is possible to solve this challenging problem in the domain of real-valued parameters, e.g. positioning coordinates, given that those are usually the parameters of interest for the user. In response to such questions, this doctoral dissertation is structured in two main parts: the first one looks in the integer ambiguity domain, presenting new flexible estimators and algorithms, ultimately merged into the new Least-squares AMBiguity Decorrelation Adjustment (LAMBDA) 4.0 toolbox; the second one focuses on the real-valued parameter domain where the integerness of ambiguities is still taken into account, where novel dual estimators are presented and an optimal, globally convergent solution is constructed using the branch-and-bound method.@en

In this contribution we introduce the dual mixed-integer least-squares problem and study it in relation to its primal counterpart. The dual differs from the primal formulation in the order in which the integer ambiguity vector a∈Zⁿ and baseline vector b∈R^p are estimated. As not the ambiguities, but rather the entries of b are usually the parameters of interest, the attractiveness of the dual formulation stems from its direct computation of b. It is shown that this potential advantage relies on the ease with which an implicit integer least-squares problem of the dual can be solved. For the convoluted cases, we introduce two methods of simplifying approximations. To be able to describe their quality, we provide a complete distributional analysis of their estimators, thus allowing users to judge whether or not the approximations are acceptable for their application. It is shown that this approach implicitly introduces a new class of admissible integer estimators of which we also determine the pull-in regions. As the dual function is shown to lack convexity, special care is required to be able to compute its global minimizer bˇ. Our proposed method, which has finite termination with a guaranteed ϵ-tolerance, is constructed from combining the branch-and-bound principle, with a special convex-relaxation of the dual, to which the projected-gradient-descent method is applied to obtain the required bounds. Each of the method’s three constituents are described, whereby special emphasis is given to the construction of the required continuously differentiable, convex lower bounding function of the dual.

@en

In this work we introduce the LAMBDA 4.0 toolbox, which provides an enhanced implementation for integer estimation, validation, and success rate evaluation. This free and open-source toolbox is a major update to LAMBDA 3.0 (2012), while it also integrates the functionalities from Ps-LAMBDA 1.0 (2013), thus respectively merging both estimation and evaluation capabilities. The new implementation provides redefined algorithms, such as an improved integer search strategy with a one-order reduction in the computational time, along with additional estimators: Vectorial Integer Bootstrapping (VIB), Integer Aperture Bootstrapping (IAB) and Best Integer Equivariant (BIE). This toolbox aims to become a valuable resource for researchers and/or practitioners dealing with mixed-integer models in high dimensions, e.g., terrestrial-based carrier-phase systems, multi-constellation Global Navigation Satellite Systems (GNSS), or other interferometric applications.

@en

Enceladus, with its subsurface ocean, is amongst the top priority targets in the search for life beyond Earth. Following on discoveries from the Cassini mission that Enceladus possesses a global subsurface ocean containing salt and organic compounds, there are many unconstrained properties of the ocean and ice shell that must be investigated to further assess the habitability of Enceladus and begin the search for biosignatures on Enceladus. In this paper, we present a concept study for a New Frontiers class multi-lander and orbiter mission to Enceladus that investigates if there is or ever was a habitable environment on Enceladus. The mission architecture includes an orbiter for detailed chemical analysis of material erupted from Enceladus’ plumes and four impact landers for geophysical measurements. As part of our mission concept study, we explore key trades for orbital and surface science, as well as assess the scientific potential and hazards of candidate landing sites on Enceladus. The novelty of our mission architecture and consideration of both orbital and surface science elements makes this work directly relevant to a broad range of potential future mission architectures under consideration, such as those identified in the 2023–2032 Planetary Science and Astrobiology Decadal Survey.

@en

In this contribution, we extend the principle of integer bootstrapping (IB) to a vectorial form (VIB). The mathematical definition of the class of VIB-estimators is introduced together with their pull-in regions and other properties such as probability bounds and success rate approximations. The vectorial formulation allows sequential block-by-block processing of the ambiguities based on a user-chosen partitioning. In this way, flexibility is created, where for specific choices of partitioning, tailored VIB-estimators can be designed. This wide range of possibilities is discussed, supported by numerical simulations and analytical examples. Further guidelines are provided, as well as the possible extension to other classes of estimators.

@en

The estimation of satellite orbits and clocks plays a central role in different Global Navigation Satellite System (GNSS) domains, e.g. precise positioning or time transfer. Such products can be computed in the process of Orbit Determination and Time Synchronization (ODTS), which relies on a network of ground-based stations, well distributed around the globe. The mm-level precision of carrier-phase measurements is exploited in this network estimation following a correct resolution of their ambiguities. For several stations, satellites and/or signals, thousands of ambiguities have to be processed, which means having to deal with high-dimensional ambiguity resolution (HDAR) problems. In this research work, we firstly account for the impact of ambiguity resolution in a varying network size, based on GPS-only, Galileo-only and GPS+Galileo configurations. Using 25 or more stations, the accuracy (1D RMS orbital errors) of fixed solutions reaches a plateau at 1-2 cm. Hence, we focus on a small global network of 14 stations, where the model strength decreases, so does the reliability of the ambiguity fixing process and advantages over a float solution might become less evident. In order to allow reliable HDAR, two implementations of the Vectorial Integer Bootstrapping estimator are presented and evaluated with respect to their scalar counterpart. Finally, it is shown how the proposed fixing processes are more robust, still very efficient, and on certain days they provide a large improvement to satellite products. The orbital results are ultimately validated by considering the satellite midnight discontinuity errors over a 3-month period in 2019.@en

This study is focused on the effect of solar flux conditions on the dynamics of gravity waves (GWs) in the thermosphere. Air density and crosswind in situ estimates from the Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) accelerometers are analyzed for the whole mission duration. The analysis is performed in the Fourier spectral domain averaging spectral results over periods of 2 months close to solstices. A new GW marker (called C³ _f) is introduced here to characterize GWs activity under low, medium, and high solar flux conditions, showing a clear solar damping effect on GW activity. Most GW signal is found in a spectral range above 8 mHz in GOCE data, meaning a maximum horizontal wavelength of around 1000 km. The level of GW activity at GOCE altitude is strongly decreasing with increasing solar flux. Furthermore, a shift in the dominant frequency with solar flux conditions has been noted, leading to larger horizontal wavelengths (from 200 to 500 km) during high solar flux conditions. The correlation between air density variability and GW marker allows to identify most of the large-amplitude perturbations below 67° latitudes as due to GWs. The influence of correlated error sources, between air density and crosswinds, is discussed. Consistency of the spectral domain results is verified in the time domain with a global mapping of high-frequency air density perturbations along the GOCE orbit. This analysis shows a clear dependence with geomagnetic latitude with strong perturbations at magnetic poles and an extension to lower latitudes favored by low solar activity conditions. These results are consistent with previous Challenging Minisatellite Payload (CHAMP) data analysis and with general circulation models.

@en