Gravity Field Recovery Using High-Precision, High–Low Inter-Satellite Links

: Past temporal gravity ﬁeld solutions from the Gravity Recovery and Climate Experiment (GRACE), as well as current solutions from GRACE Follow-On, suffer from temporal aliasing errors due to undersampling of the signal to be recovered (e.g., hydrology), which arise in terms of stripes caused by the north–south observation direction. In this paper, we investigate the potential of the proposed mass variation observing system by high–low inter-satellite links (MOBILE) mission. We quantify the impact of instrument errors of the main sensors (inter-satellite link and accelerometer) and high-frequency tidal and non-tidal gravity signals on achievable performance of the temporal gravity ﬁeld retrieval. The multi-directional observation geometry of the MOBILE concept with a strong dominance of the radial component result in a close-to-isotropic error behavior, and the retrieved gravity ﬁeld solutions show reduced temporal aliasing errors of at least 30% for non-tidal, as well as tidal, mass variation signals compared to a low–low satellite pair conﬁguration. The quality of the MOBILE range observations enables the application of extended alternative processing methods leading to further reduction of temporal aliasing errors. The results demonstrate that such a mission can help to get an improved understanding of different components of the Earth system.


Introduction
In times of a changing climate the need for innovative observation techniques for capturing geophysical processes in the Earth system becomes increasingly urgent. In this context, the observation of the temporal gravity field by satellites from space play an important role when investigating, e.g., rapid changes in the cryosphere, oceans, water cycle, and solid Earth processes on a global scale. For the determination of temporal gravity fields, in the last decade satellite missions such as GRACE [1] or Challenging Minisatellite Payload (CHAMP) [2,3] orbited around the globe and helped to get a better understanding of the Earth's mass flux signals. CHAMP was based on high-low satellite-to-satellite tracking (SST) exploiting the Global Positioning System (GPS) [4] over a time span of 10 years. The accuracy of the CHAMP orbit information of 2-3 cm [5] derived from GPS allowed for resolving only the long wave range of the time varying gravity field, with spatial scales of ≈1000 km, e.g., Baur 2013 [6]. Analyzing the perturbed orbit of other Low Earth Orbiters (LEO), such as the Swarm satellites, allows for a similar performance [7]. The GRACE mission reached spatial scales of the temporal gravity field of ≈300 km and below due to a combination of K-band microwave low-low inter-satellite ranging between two identical satellites following each other in the same orbit at a distance of about 220 km with micrometer precision, and high-low GPS satellite-to-satellite tracking plus accelerometer observations. These missions improved our knowledge of water mass variations on the continents, in the oceans, and the atmosphere to a great extent. Additionally, the static gravity observation equations and stochastic modelling. In Section 4 the estimated gravity field solutions are analyzed and assessed. The main conclusions are summarized in Section 5, and in Section 6, a short outlook is given.

Observation Geometry
In contrast to the past GRACE and the current GRACE Follow-On missions, which are mainly based on LEO satellites (several hundred km), the MOBILE minimum configuration consists of a constellation of two high and one low orbiting satellites. As done for GRACE and GRACE Follow-On, the main observable is the gravity-induced inter-satellite distance change, which is in case of MOBILE measured between medium orbiting satellites (MEO; several thousand km) and LEO satellites though. As a second gravity observation type, high-precision orbit positions based on Global Navigation Satellite System (GNSS) orbit determination are used. This idea of high-precision high-low tracking was first investigated by Hauk et al. 2017 [26] using the inter-satellite link technique as part of the payload on-board Galileo satellites of future generation in connection with LEO satellites, where the main error sources and the corresponding achievable performance were analyzed. Due to the fact that the MOBILE constellation presents a stand-alone concept without the need to place an additional payload on another space infrastructure, and the very large distance between the high-and the low orbiting satellites, which plays a crucial role in the framework of high-precision high-low tracking, dedicated MEO satellites were included in the concept. It should be emphasized that alternatively, a constellation of one MEO and two LEO satellites could be envisaged, which turns out to provide nearly the same performance as the proposed one, but might be more expensive due to the need to build and maintain at least two LEO satellites in orbit.   Table 1. The MEO satellites orbit at an altitude of about 10,150 km in the same orbital plane, separated by an 180-degree mean anomaly as alternating targets of the LEO satellites, in order to maximize the visibility and thus observation time. The LEO satellite is orbiting in an altitude of about 360 km. Both LEO and MEO satellites are flying in polar orbits in order to maintain a long-term stable formation (no relative drifts of the orbit planes). Additionally, two LEO satellites with near-polar orbits flying in an altitude of about 470 km with an inter-satellite distance of 200 km are set up in order to perform comparability studies between the MOBILE constellation and a GRACE Follow-On-like mission. All orbits have certain repeat cycles, after which the satellites reach the same position on Earth again in order to maintain a stable ground track pattern and related stable gravity model quality. The choice of the orbit height of the MEO satellites underlies three major constraints: (1) A high altitude of several thousand kilometers is necessary in order to ensure long observation periods and preferably measurements of multi-directional distance variations, with a strong dominance of the radial component, resulting in a close to isotropic error behavior of the retrieved gravity field solution (see Section 4). (2) The distance between a MEO-LEO pair must not be too large, because the larger the distance, the more difficult it is to fulfill the 1-µm accuracy requirement for the inter-satellite link established by a laser range interferometer (see Section 2.2). (3) The third constraint is driven by solar radiation belts encircling the Earth in which energetic charged particles are trapped inside the Earth's magnetic field [27], which are of different intensities dependent on the solar cycle, altitude, and inclination of the satellite orbit. As a result, altitude ranges of several thousand kilometers below the chosen orbit height drop out. These conditions connected with the repeat orbit lead to an altitude of about 10,000 km for the MEO satellites.
The high-low tracking concept enables a multi-directional observation geometry with differing elevation angles from 3 • (assumed minimum elevation angle of visible MEOs observed by the LEO) up to a near-radial direction. However, due to the observation geometry of the MEO-LEO satellite pairs and the changing satellite links from one to the other MEO, data gaps arise for every satellite pair, leading to a non-continuous measurement time series of these pairs. For the simulated MOBILE constellation, this results in a ranging window maximum of 45 min, and a maximum data gap of 18 min. The separation of the two MEO satellites of 180-degree mean anomaly is chosen to keep the time period of the data gap as small as possible. In Figure 2, the LEO ground track of the MOBILE concept is displayed for 1 day together with the corresponding elevation angles.

Instrumentation
The main observable in the MOBILE mission are range measurements from the LEO to the MEOs, where the MEOs are alternating targets. The ranging accuracy is on the micrometer level in order to be sensitive for gravitational forces and its changes on Earth. For distances of several thousand kilometers, a laser-based distance measurement system can reach such an accuracy. The laser range interferometer is placed at the LEO satellite, while the MEOs are equipped with passive reflectors or transponders. In case of the GRACE Follow-On mission, the measurement of inter-satellite ranges by laser range interferometry (LRI) has been successfully established. The link between the two satellites was generated with an active laser on one satellite, and a phase-locked amplifying transponder on the second spacecraft [10]. For the MOBILE concept, the laser ranging instrument needs to be adapted due to the very large distance and the relative motion of the LEO and MEO satellites. In contrast to the GRACE Follow-On, the large distance and the relative speed lead to a range of Doppler shifts of several GHz compared to a few MHz, which causes the need of a reference laser source with a larger range of reference frequencies and a faster phase-tracking capability than implemented for the GRACE Follow-On. The required parameters (<10 GHz range, <10 MHz/s tracking) are within the range of existing, space qualified reference lasers (e.g., the one used for the ATmospheric LIDar (ATLID) instrument on the Earth Clouds Aerosols and Radiation Explorer (CARE) mission) [28], but their compatibility with the needs of an interferometric instrument has to be the subject of further studies. Due to the relative motion of the LEO and MEO satellites, pointing tracking capabilities are required, which requires a modified link implementation. The LEO satellite is selected to play the active part in the tracking mechanism, while the partner satellites (MEOs) are equipped with passive retroreflectors. This type of laser tracking and ranging has been successfully performed for decades with active laser systems on the ground and passive retroreflectors on satellites in orbit (e.g., Laser Geodynamics Satellite (LAGEOS), Ball Lens In The Space (BLITS)) [29,30]. The scientific benefit of deploying a passive payload in space is the significantly increased mission duration when compared to complex active payloads. The main technological challenge in utilizing this setup for an LRI instrument is the need to achieve a sufficiently high level of retrieved power without the need for amplification between the two passes, ideally close to the 80 pW received by the GRACE Follow-On implementation, but at least to levels above ≈1 pW in order to allow phase tracking. The main design factors impacting the received power are the initial output power, the size of the retroreflector, and the size of the receiving telescope.
Satellite on-board sensors play an important role in the gravity field retrieval by influencing satellite observations due to correlated noise. In our study, the error assumptions used for the laser ranging instrument in the MOBILE concept are based on the time-series provided by Schäfer et al. 2013 [31], which originated in connection with ESA's GETRIS study, and show micrometer ranging accuracy around 1 MHz. Due to simulation purposes, this time-series was adapted by means of cascaded second order Butterworth auto regressive moving average (ARMA) filter model. The spectral behavior of the LRI is shown in Figure 3 (light green curve) in terms of an amplitude spectral density (ASD). The relative distance measurement errors assumed for the low-low satellite pair are identical to those used in the frame of the ESA-Assessment of Satellite Constellations for Monitoring the Variations in Earth's Gravity Field (SC4MGV) project [32], provided from the consultancy support of Thales Alenia Space Italia, and show a performance of about several 10 nanometers. The corresponding analytical noise model of the used laser interferometer is given by the ASD in terms of range-rates ( Figure 3, light blue curve): The generation of all noise time-series was done by scaling the spectrum of normally distributed random time-series with their individual spectral model. The non-gravitational forces are typically sensed by the on-board accelerometers located in the center-of-mass of the satellite. In case of the LEO satellites, the implementation of an accelerometer is absolutely necessary due to air drag as the main contributor. For the low-low pair a GRACE-like electrostatic accelerometer is assumed with two highly sensitive axes oriented in the flight direction (largest signal) and in the radial direction, and one low-sensitive axis in the cross-track direction (see Figure 4, blue and red curves). The accuracy level in terms of accelerations is derived by Iran Pour et al. 2015 [32], and is expressed by: with x denoting along-track, y across-track, and z (close to) the radial direction. Based on the heritage of previous gravity missions for MOBILE, we seek a resolution on the level of 10 −11 m/s 2 , which is the same as assumed for the along-track and radial axes of the accelerometer on-board the low-low satellite pair, but ideally with the same performance in all three directions. Furthermore, the slope at frequencies from 10 −3 Hz and lower is pressed down from 1/f 2 for the low-low pair to 1/f for MOBILE. The performance of the relative acceleration measurement error is displayed in Figure 4 (green curve). While an accelerometer is mandatory for the MOBILE LEO satellite, for the MEOs, less stringent requirements might apply because of the substantially smaller amplitude of the signal and the fact that non-conservative forces can be modelled much more accurately in high altitudes. Also, the design of the MEO could be optimized for the high predictability of non-gravitational forces, e.g., by implementing very simple geometrical surfaces wherever radiative pressure is relevant. In spite of these facts, in the MOBILE concept, the implementation of accelerometers is proposed. Geo-location of satellite observations, as well as gravity retrieval, require highly accurate continuous orbit determination, making GNSS space receivers on all satellites obligatory. In our simulations, we assume an absolute kinematic positioning on a cm level. Using a laser ranging instrument as the main measurement system requires exact pointing of the tracking antenna in the order of 10 µrad or less, and therefore the implementation of systems for attitude determination and control. We assume star camera sensor errors for all satellites represented as rotation angles around the along-track (roll), cross-track (pitch), and radial (yaw) axes, expressed by the ASD of the following analytical noise models [32]: In addition, for the MOBILE LEO satellite, a drag-reduction system needs to be implemented in order to maintain the orbit, and not to saturate the accelerometers due to non-gravitational accelerations. The MEO satellites will very likely require an electrical propulsion system to move to their target orbit from the lower separation altitude achievable with a low-cost launcher.

Simulation Environment
All simulations were executed with a full numerical mission simulator [33,34], which has already been successfully applied to recover satellite-only gravitational field models from GOCE data [35]. The simulation environment is based on numerical orbit integration, following a multistep method for the numerical integration according to Shampine & Gordon 1976 [36], which applies a modified divided difference form of the Adams predict-evaluate-correct-evaluate (PECE) formulas and local extrapolation. According to this method, the order and the step size are adjusted to control the local error per unit step in a generalized sense. The generation of "true" dynamic orbits and, subsequently, the "true" GNSS high-low SST and low-low laser ranging SST observations, is done by adding different force models according to the "true" world of Table 2. The impact of orbit errors on the gravity field processing is taken into account as well by propagating 1 cm white noise of the integrated orbit positions of each satellite. The resulting erroneous dynamic orbits serve as computational points for the reference values of the observations and enable the computation of the GNSS high-low SST observations in three directions. In order to ensure the accuracy of the inter-satellite link, error-free dynamic orbits are used for the reference values of the low-low SST observations from the laser interferometer system, which are expressed in terms of range-rates. The adopted gravity field approach is based on a modification of the integral equation approach from Schneider 1969 [37] where the orbit is divided into continuous short arcs of 6 h length, and the position vectors at the arc node points are set up as unknown parameters, which are estimated together with the gravity field coefficients. This technique has already been successfully applied in real data applications to recover satellite-only gravitational field models for CHAMP and GRACE [38] (ITSG-Grace2016) [39]. The functional model follows the typical formulation used for low-low SST missions like GRACE, which comprises a high-low SST and a low-low SST component. Position differences between two satellites are used for the computation of the reference values for the high-low SST part of the observation system, whereas the reference values for the low-low SST part are derived by projecting position and velocity differences between two satellites onto the line-of-sight, leading to the computation of inter-satellite range-rates. Table 2 gives an overview of the force and noise models used in the processing for the "true" and "reference" world. The static gravity field model is represented by the GOCO03s model, which is a satellite-only gravity field model based on GRACE, GOCE, and LAGEOS [40]. In order to simulate geophysical signals, ESA's updated Earth system model [41] has been used, which contains the five main geophysical signal components atmosphere (A), ocean (O), hydrology (H), ice (I), and solid Earth (S) with a time resolution of six hours, linearly interpolated to the epochs. The Earth system model covers the time period 1995-2006, and contains plausible variability and trends in both low-degree coefficients and the global mean eustatic sea level. It depicts reasonable mass variability all over the globe at a wide range of frequencies including multi-year trends, year-to-year variability, and seasonal variability, even at very fine spatial scales, which is important for a realistic representation of spatial aliasing and leakage. The impact of ocean tide model errors is assessed by taking the difference of two tide models, EOT11a [42], and GOT4.7 [43]. The total stochastic model for the observations is approximated individually for both satellite formations by means of a cascade of digital Butterworth ARMA filters [44,45]. Filter coefficients are chosen in such a way that the cascade's frequency response optimally matches the inverse of the amplitude spectrum of the previously generated pre-fit residuals. They are estimated as a result of the computation of the linearized normal equations, which include differences between the "true" (only the static GOCO03s gravity field model and sensor noise are included) and the reference observations (only the static GOCO03s gravity field model is included), such that the error sources from the sensors are considered exclusively. Assuming uncorrelated high-low and low-low SST observations, weighting matrices are set up for all observation components separately.
The goal is the retrieval of all spherical-harmonic (SH) coefficients up to a maximum SH degree of 100 from observations sampled every 5 seconds for the first 30 days of the year 2001. Due to the fact of non-linear observation equations, the "reference" observations are reduced from the "true" observations as a result of the linearization process. The gravity field parameters are estimated by solving full normal equations of a least squares system based on a standard Gauss-Markov model using weighted least squares with stochastic models in accordance with the simulated instrument noise levels. The resulting gravity field coefficients are analyzed and compared regarding quality and performance in terms of retrieval errors by removing a monthly average of the true mass transport model from the recovered signal.

Gravity Field Retrieval Performance Due to Instrument Errors
At first the impact of the instrument errors on the gravity field retrieval were quantified. For this task, we performed simulations where each error source according to the assumptions described in Section 2.2 was treated individually. Figure 5 shows the gravity field retrieval performance in terms of equivalent water height (EWH) errors per SH degree per coefficient for the low-low pair constellation and the MOBILE concept. Furthermore, the results were quantified using global RMS values of the errors in the recovered signal expressed in terms of cm of EWH, listed in Table 3 (see part: instrument errors). If only white-noise positioning errors were considered ( Figure 5, green curves), the gravity field retrieval performance mainly depended on the observation geometry. The comparison between both satellite concepts revealed strongly reduced retrieval errors for MOBILE, which benefited from multi-directional observations. In the case of accelerometer noise in combination with star camera errors ( Figure 5, blue curves), the MOBILE constellation showed reduced error behavior compared to the low-low pair as well. This was mainly caused by the observation geometry, but also by the improved accelerometers (≈23%) with 3D capabilities to certain parts in the case of MOBILE. In contrast, the retrieval performance of the low-low pair benefits from the nanometer accuracy of the laser interferometer compared to the micrometer accuracy of MOBILE's laser link sensor for the most part of the spectrum (>SH degree 10). This became evident when only laser interferometer noise was considered ( Figure 5, red curves). However, in the very low degrees, the MOBILE concept performed better than GRACE, which was again owed to the fact of an improved observation geometry. When considering all instrument error sources together, the retrieval errors of the low-low satellite pair are dominated by the accelerometer plus star camera sensor performance, while for the MOBILE constellation, the laser link error was the dominating error source for SH degrees higher than 20, and the accelerometer plus star camera noise only dominated the spectrum in the lower degrees. These results led to the conclusion that the gravity field retrieval based on instrument error sources showed smaller errors below SH degree 40 for the MOBILE concept compared to the low-low pair, but increased errors in the higher frequency spectrum due to the lower accuracy of the laser interferometer.  Next to the estimation of SH coefficients, we estimated their formal errors as well, shown in Figure 6. The noise of the different sensors in combination with the observation geometry reveal the performance of a specific satellite concept. In our case, they demonstrated the impact of the MOBILE high-low tracking concept by showing an almost uniform (isotropic) error spectrum and a high sensitivity in the sectorial coefficients (SH degree equal to SH order). In case of the low-low pair configuration especially, the sectorial coefficients were less well-determined than the zonal coefficients (SH order equal to zero). Figure 6b,d gave a closer view of the formal errors located in the long wavelength (low-degree) spectrum. The comparison between MOBILE and the low-low pair led to the assumption that the determination of the very low SH coefficients could be accomplished with a higher sensitivity through the MOBILE concept. In contrast to the observations in the along-track direction of the low-low pearl-string configuration, the multi-directional observations of the high-low tracking concept with a strong dominance of the radial component enabled an improved estimation of the very low SH coefficients. The close to radial observation geometry of MOBILE was comparable to satellite laser ranging (SLR) observations, showing superior performance in observing the very long wavelength gravity field variations, in particular the zonal SH coefficient of degree 2, which physically represented the Earth's dynamic oblateness [46]. In order to make the effect of the different error spectra of both satellite concepts even more visible, spatial covariance functions were computed for a position at the equator and at 45 • latitude (see Figure 7). They describe the correlation of the computation point with its neighborhood in the normal equation system due to the used stochastic model and the observation geometry. The spatial characteristics and the pattern of the covariances provide information about the spatial behavior of the retrieved signals. In our case, the figures show the typical stripes for the low-low tracking concept caused by the north-south observation direction that are known from the GRACE temporal gravity models, while the MOBILE concept exhibited an isotropic error structure at both latitudes.

Temporal Gravity Field Retrieval
The retrieval of the temporal gravity field is dominated by temporal aliasing errors due to the undersampling of high frequency geophysical signals and imperfect de-aliasing models, which has already been shown by, e.g., References [47][48][49] for the GRACE mission. In order to analyze the impact of different time-varying mass signals on gravity field retrieval, we performed simulations by using signals that were subdivided into non-tidal AOHIS, HIS, and tidal signals, including the instrument errors described in Section 2.2. Figure 8 displays the corresponding retrieval errors for both satellite concepts. The results indicate that the errors with the highest signal amplitudes were related to AOHIS signals (Figure 8, red curves), and in particular to atmospheric and oceanic signals. Tidal aliasing effects played a key role in the total error budget as well (Figure 8, blue curves) by representing the highest aliasing errors next to non-tidal atmosphere and ocean aliasing errors. In this context it is important to mention that errors in ocean tide models are considered as one of the major sources of error in the determination of temporal gravity field models from GRACE data [50,51]. Our simulations show that the MOBILE configuration can reduce non-tidal aliasing errors (≈45%) as well as tidal aliasing errors (≈30%) over the whole spectrum significantly (see also Table 3, part: temporal aliasing errors). Despite the fact that for the high-low tracking concept the assumed arrangement of satellites causes incomplete data time series, the multi-directional observation geometry enables the sampling of time varying signals with reduced aliasing errors compared to the low-low pair configuration. Usually high-frequency mass signals are a priori reduced based on atmosphere and ocean de-aliasing (AOD) products [52], and ocean tide de-aliasing models. The resulting temporal gravity field models thus contain mainly information on sub-seasonal, seasonal, and secular continental hydrological mass variations and ice mass variations on Earth [53,54], and solid Earth signals related to glacial isostatic adjustment (GIA), and co-and post-seismic gravity changes of big earthquakes. For the analysis of such mass flux signals we performed simulations by using only the HIS signal. The resulting gravity field retrieval errors (Figure 8, green curves) again revealed smaller aliasing effects for the MOBILE concept (≈60%), which is even better visible when looking at the spatial domain, shown in Figure 9. As already suggested by Figure 7, the retrieved HIS fields demonstrate, that the error pattern of MOBILE was much more homogeneous, and the typical striping of a low-low along-track ranging system is significantly reduced, particularly in the equatorial regions where the orbit ground tracks were less dense. This resulted in a clearly improved free representation of hydrological and ice mass signals for the MOBILE concept. The high quality of multi-directional observations of the high-low tracking concept allows the application of an extended alternative processing method first proposed by Wiese et al. 2011 [23], which enables the mitigation of temporal aliasing effects due to non-tidal time varying signals, as it was stated in Section 1. Wiese et al. 2011 [23] demonstrated the benefit of the co-parameterization of additional daily low degree and order gravity field coefficients for Bender-type satellite constellations. We investigated the potential of this methodology regarding the MOBILE concept by simulating a monthly solution while co-estimating daily gravity fields up to SH degree and order 10, including HIS signal plus instrument errors. The resulting retrieval errors are displayed in Figure 8 (magenta curve). They revealed an error reduction of about 40% compared to the nominal solution (green curve), which led to an increased spatial resolution of about SH degree 50 (≈400 km) instead of 40 (≈500 km). The corresponding spatial plot (Figure 9, f and g) shows a global reduced aliasing pattern, especially in higher latitudes. The comparison between the true HIS signal and the MOBILE recovered signal (nominal and extended processed) displayed in Figure 9 shows that the quality of the solutions could be improved to such a level that de-striping and smoothing the solutions was no longer necessary when examining signals to degree and order 50. Therefore, a possible loss of signal by a posteriori filtering of the gravity field solutions recovered by MOBILE could be avoided.

Conclusions
In this study, we investigated the gravity field retrieval performance of the novel and innovative MOBILE high-low satellite tracking concept and compared it with a low-low GRACE Follow-On-like configuration qualitatively and quantitatively. Based on full numerical simulations, gravity field parameters were estimated in terms of SH coefficients by solving a least-squares system by inverting full normal equations over a time span of 1 month. The most important error sources affecting the gravity field retrieval performance, key instruments on-board the satellites, as well as time varying mass flux signals, were included in order to assess their impact on gravity field retrieval for both mission concepts.
The results regarding the instrumental impact on the gravity field solution show that the performance of the MOBILE configuration was mainly limited by the assumed micrometer accuracy of the laser interferometer, especially in the short wavelength spectrum, while the performance in the lower wavelengths of the gravity field benefited from the multi-directional observation geometry and optimized 3D accelerometer. In contrast, the gravity field retrieval of the low-low pair constellation was limited mainly by the accelerometer, which predominated the nanometer accuracy of the assumed laser interferometer. The multi-directional observations of MOBILE mentioned above included a strong radial component and led to an almost uniform (isotropic) error spectrum, while the low-low tracking concept showed the typical stripes caused by the north-south observation direction. However, the high accuracy of the low-low satellite pair's inter-satellite link led to an improved gravity field performance from SH degree 40 and higher compared to MOBILE, which performed better in the long wavelength spectrum where the largest amplitudes of time varying gravity field signals occurred.
The benefit of MOBILE's multi-directional observation geometry arose when including tidal and non-tidal mass variation signals into the simulation process. The results revealed significantly reduced temporal aliasing errors in the recovered gravity field signal compared to the low-low tracking concept over the whole spectrum. In the case of the separate treatment of the HIS signal, the resulting gravity field error performance of MOBILE improved even by about 60%, and the application of an extended processing method to reduce temporal aliasing errors by co-estimation of daily gravity field parameters, led to a further reduction of retrieval errors of about 40%. Furthermore, the results show that the quality of recovered MOBILE gravity field solutions could make a treatment of such solutions using a posteriori filtering techniques obsolete.
The gravity field solutions retrieved using MOBILE can contribute to an improved understanding of different components of the Earth system, such as the estimation of continental water storage and freshwater fluxes, the quantification of large-scale flood and drought events and their monitoring and forecasting, or understanding the mass balance of ice sheets and larger glacier systems, just to name a few. The application of the extended processing method implied the co-parameterized gravity field parameters (which, in our case, are daily gravity fields) as a side product. These daily solutions with low spatial resolution could aid in improving atmospheric models, and possibly be beneficial to the oceanography community as well, as many of these short-term signals have large spatial scales.

Outlook
The gravity field solutions of the high-low tracking concept presented in this paper were based on the minimal configuration of MOBILE. On top of this scenario, the optional implementation of a third or fourth MEO satellite, but also a second LEO, could be considered to further increase the mission performance, but also significantly improve the temporal resolution. Due to the largely passive instrumentation of this mass transport mission, the function of the MEO satellites could be implemented as a backpack application of other MEO missions, such as the Galileo next-generation satellites, in order to extend and maintain the infrastructure for laser ranging payloads. Aside, one of the most important fields of research is the mitigation of temporal aliasing errors. In this context it is important to mention that aliasing effects due to imperfect ocean tide models represent one of the largest error sources in temporal gravity field retrieval. The capability of the multi-directional observations of MOBILE of co-parameterize tidal parameters over long time spans, as proposed in Reference [24], in order to improve current ocean tide models will be the subject of further study.