Triple-Pair Constellation Configurations for Temporal Gravity Field Retrieval

: The goal of next-generation gravity missions (NGGM) is to improve the monitoring of mass transport in the Earth system by an increased space-time sampling capability as well as higher accuracies of a new generation of instrumentation, but also to continue the monitoring time series obtained by past and current missions such as GRACE and GRACE Follow-On. As the likelihood of three satellite pairs being simultaneously in orbit in the mid-term future increased, we have performed a closed-loop simulation to investigate the impact of a third pair in either polar or inclined orbit as an addition to a Bender-type constellation with NGGM instrumentation. For the additional pair, GRACE-like as well as NGGM instrumentation was tested. The analysis showed that the third pair mainly increases the redundancy of the monitoring system but does not significantly improve de-aliasing capabilities. The best-performing triple-pair scenario comprises a third inclined pair with NGGM sensors. Starting with a Bender-type constellation of a polar and an inclined satellite pair, simulation results indicate an average improvement of 11% in case of adding the third pair in a near-polar orbit, and of 21% for the third pair placed in an inclined orbit. The most important advantage of a multi-pair constellation, however, is the possibility to recover daily gravity fields with higher spatial resolution. In the case of the investigated triple-pair scenarios, a meaningful daily resolution with a maximum spherical harmonic degree of 26 can be achieved, while a higher daily parametrization up to degree 40 results in spatial aliasing and thus would need additional constraints or prior information.


Introduction
Dedicated gravimetric satellite missions like the Challenging Minisatellite Payload (CHAMP; [1]) and the Gravity Recovery and Climate Experiment (GRACE; [2]) and GRACE Follow-On (FO; [3]) missions have been providing, for nearly two decades, essential observations of the changes of the Earth's gravity field on a global scale. This monitoring is fundamental for applications in Earth sciences, such as hydrology [4], atmosphere [5], plate tectonics [6], earthquakes [7], and glacial isostatic adjustment (GIA; [8]) as well as cryosphere [9,10].
The GRACE-FO mission was launched in May of 2018 to continue monitoring the temporal gravity field of the Earth, after the GRACE mission was decommissioned due to a battery failure at the end of 2017. However, the observed data from GRACE/GRACE-FO have limitations regarding spatial and temporal resolution, homogeneity, as well as sensitivity. The user requirements collected by [11] underline the need for long-term, sustained gravity observations with increased spatial and temporal resolution to observe small-scale, short-time mass transport phenomena.
To find optimal mission concepts with improved spatial and temporal resolutions, various studies have been performed: The authors of [12][13][14][15][16] investigated the performance of single-pair missions. The authors of [17][18][19][20] analyzed flying a second inclined pair in addition to an in-line pair in polar orbit (furthermore called Bender-type constellation), thus improving the space-time sampling and having available multi-directional observations as compared to a single polar-pair due to the two different orbit planes.
Mission proposals in response to European Space Agency (ESA) Earth Explorer calls 8, 9, and 10 like e.motion with a single pair in a pendulum constellation [21], e.motion2 [22] based on a Bendertype concept, or the innovative high-low tracking formation mission MOBILE [23] also achieve improved error characteristics. While the e.motion and e.motion2 mission proposals are based on the observation system of low-low satellite-to-satellite tracking (LL-SST) as implemented in GRACE and GRACE-FO, the MOBILE concept proposes a different observation strategy of high-precision highlow tracking as well as new instrumentation. All three proposals emphasize the need for continuous and sustained observations of the Earth's mass transport on a global scale from space.
A study by [24] has concluded that the gain in accuracy of the Bender double-pair constellation is higher compared with the single-pair pendulum formation and is therefore preferred. The ESAfunded SC4MGV study (Assessment of Satellite Constellations for Monitoring the Variations in Earth Gravity Field, [15]) found that there is a certain degree of freedom in the orbit design of Bender-type mission concepts. In the follow-on ESA-study ADDCON (Additional Constellation and Scientific Analysis of the Next Generation Gravity Mission Concept; [25]) an analysis regarding orbit height and instruments showed that a Bender pair with improved instrumentation gave the best results for a seven-day reference solution [26].
In the context of observing short period signals within short latencies, daily solutions in nearreal time are needed. However, the potential to use the retrieved gravity fields for applications with high temporal variations, like hydrology and earthquake monitoring, is limited by the requirement to recover gravity field variations as detailed as possible, i.e., with highest spatial resolution. An increase in temporal resolution means less observations per analysis period and therefore a reduced redundancy in the parameter estimation process. For possible NGGM constellations with at least two pairs, the co-estimation of long-wavelength gravity field solutions for short time periods such as one day, further called Wiese approach, was successfully tested in closed-loop simulations [20,27].
As various space agencies are considering flying a GRACE-like or NGGM mission, the potential of triple pairs compared to a Bender-type double pair is of interest. The most significant advantage of an additional pair to a Bender-type constellation, apart from an improved performance, is the possibility to increase and resolve daily Wiese solutions with higher spatial resolution. Additionally, for time-critical monitoring applications (e.g., early-warning and forecasting systems) information in near-real time (NRT) is of essence, meaning an NRT processing scheme [28] is necessary. This paper evaluated three scenarios with two and three pairs and analyses their potential for NRT processing and the highest possible daily resolution for applications such as the monitoring of hydrological extremes such as droughts and floods.
The paper is organized as follows: Chapter 2 reviews the satellite mission design as well as the space-time sampling of various applications. The methodology of the full-scale closed-loop simulation, based on the software available at the Institute of Astronomical and Physical Geodesy (IAPG) of the Technical University Munich (TUM), is described in Chapter 3. Chapter 4 contains the obtained gravity results in the spatial as well as the frequency domain, and finally, the conclusions regarding the possible improvements on gravity field as well as on application level are presented in Chapter 5.

Orbit Constellation
Two orbits are the basis to form double-and triple-pair constellations. In Table 1, the general information regarding the used orbits is listed. The letter "p" represents the polar orbit, while the letter "i" stands for the inclined orbit. These letters are also used for the scenario description within the paper later on. Both orbits have a seven-day repeat cycle and a drift rate of 1.3°/cycle. The resulting homogeneous ground track coverage allows for stable seven-day gravity field solutions. Additionally, the dense ground track spacing over longer periods due to the drifting orbit design allows for a high-degree gravity field recovery for more extended periods according to the Nyquist-Colombo sampling rule for space-borne gravimetry [29,30].  Figure 1 shows a 3D and a schematic view of the analyzed in-line satellite pair constellations based on the Bender double-pair. The Bender double-pair constellation is comprised of a near-polar pair and an inclined pair with an inclination of 70°. The polar pair is placed in an orbit plane with an ascending node of 0°, and the inclined pair in an orbit plane with an ascending node of 90°. The polar gap of the inclined pair (orbit ground tracks in red) is clearly visible at the North Pole. This basic Bender constellation, also described in Table 2, is assumed to have an improved set of sensors on board (see Section 3.2). The third pair is then placed in an orbit plane that is perpendicular to the orbit plane of the already existing polar or inclined pair. In case of triple-pair scenario 3pip, this means a second polar pair with an ascending node of 180°, and in case of triple-pair scenario 3pii, it means a second inclined pair with an ascending node of 270° (cf. Table 2) and the same inclination as the first inclined pair. The added satellite pairs are equipped with sensors featuring either a GRACE-like (a) or NGGM (b) noise (see Section 3.2). Through this distribution, an optimal spatial and temporal resolution is enabled, which is especially relevant for the aspired daily solution with the highest possible spatial resolution.

Space-Time Sampling
Science requirements for applications in solid Earth, hydrology, ocean, ice, atmosphere, and corresponding science fields were collected by [31] and visualized in 2D bubble plots; see Figure 2. The extent of the bubbles visualizes an estimation of the necessary temporal and spatial resolution to enable their observation and, subsequently, the monitoring by gravity satellite mission constellation. The bubbles together with the expected temporal and spatial resolution of the satellite mission gives a good first indication if a mission is able to observe specific features. In Figure 2, the possible performance of double (blue) and multi-pair (light green) constellations is visualized. The lines are a combination of the published original diagram in [31] and the ground tracks as well as results of the presented scenarios. However, it has to be stated that a general statement like the shown bubble plot is only a rough indication. Additionally, the observability also depends on the signal amplitude, which is not considered in Figure 2. After the successful exploitation of CHAMP, GRACE, and GOCE (Gravity field and steady-state ocean circulation explorer) and the additional knowledge based on their observations, various studies on necessary accuracies of future gravity fields were published. The authors of [11] collected updated science and user needs by a panel of scientists representing the main fields of application of possible NGGM concepts to form consolidated requirements for a compromise of all applications with the possibility to optimize for a specific field of application.
One of the major problems of a GRACE-like single-pair mission is temporal aliasing of highfrequency mass transport signals that are mainly due to processes in the atmosphere and the oceans, which cannot be captured due to the limited temporal resolution of the mission. This requires a-priori atmosphere and ocean (AO) de-aliasing based on external models. These temporal aliasing effects are significantly reduced by multi-pairs, e.g., [20] investigated the capabilities of Bender double-pair constellations to observe the full AOHIS (atmosphere, ocean, hydrology, ice, and solid Earth) signal instead of a-priori modelling of AO and the estimation of the hydrology, ice, and solid Earth (HIS) components only. Lastly, the anisotropic error behavior that is typical for GRACE-like solutions is a well-known phenomenon and typically treated with different filter approaches [32][33][34]. Multi-pair constellations inherently have additional observations and especially the addition of observations in an inclined orbit, as in the Bender-type constellation, add signal content and improve the anisotropic error behavior, and therefore decrease the GRACE striping significantly.

Numerical Simulator
The full-scale gravity field estimation software [35,36] available at the IAPG was used to execute the closed-loop simulations. The numerical orbit integration follows a multistep method for the numerical integration according to [37] with a modified divided difference form of the Adams Predict-Evaluate-Correct-Evaluate (PECE) formulas and local extrapolation [38].
In the generation of the dynamic models, the static gravity field GOCO05s model [39] up to maximum expansion degree and order (d/o) 120 is included. To simulate the non-tidal time-varying gravity field due to mass change, the updated Earth System Model (ESM) of ESA [40], furthermore called reference AOHIS, is used as reference world. The differences of the ocean tide (OT) model GOT4.7 (Goddard Ocean Tide) model [41] and the EOT11a (Empirical Ocean Tide) model [42] are used to replicate the ocean tide model errors.
The functional model follows the typical formulation used for LL-SST missions like GRACE (cf. Table 3). The "true" dynamic orbits as well as the "true" GNSS HL-SST and LL-SST observations are additionally superimposed by the noise models described in Section 3.2 depending on the scenario. The impact of orbit errors on the gravity field processing is taken into account by propagating 1 cm white noise of the integrated orbit positions of each satellite. Table 3. Force and noise models of the "true" and "reference" world used in the full-scale simulations.

Model "True" world Reference world Static gravity field (GF) model
GOCO05s The assembling of the NEQ systems is done with spherical harmonics (SH) base functions of the Earth's gravitational potential V and is expressed by the series expansion [43]: where represents the product of the gravitational constant and the Earth's mass, a the semimajor axis of the Earth, is the fully normalized Legendre polynomial of degree and order , and are the fully normalized SH coefficients, and the location is given by the radius r (geocentric distance of the satellite), geocentric co-latitude , and longitude .
The stochastic model is approximated for each satellite pair individually by using a combination of digital Butterworth ARMA filters [44,45] that best represent the amplitude spectral density (ASD) of the pre-fit residuals of a noise-only computation. Assuming uncorrelated high-low and low-low SST observations, weighting matrices are set up for all observation components separately.
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 following the simulated instrument noise levels. From the resulting gravity field coefficients, the average of the true mass transport model from the same period is removed to enable the analyses of quality and performance of the gravity retrieval.

Noise
While the basic Bender-type constellation is assumed to have a set of improved sensors, the additional satellite pairs are simulated with a GRACE-like or a NGGM noise setting. The GRACElike noise represents a noise level of accelerometer (ACC) and satellite-to-satellite tracking (SST) resembling the error characteristics of the instruments implemented on the GRACE mission, and an NGGM noise scenario with improved ACC and laser ranging interferometer (LRI) noise characteristics. The HL-SST observable always has the same white noise of 1 cm propagated along the orbit. The satellites are assumed to fly in drag-compensation mode in approximately 350 km orbit altitude so that the most significant parts of the non-gravitational forces are compensated by a propulsion system consisting of ion thrusters.

GRACE-like noise
An analytical noise model characterizes the main measurement unit, the K-Band Ranging system (KBR) to observe the inter-satellite distances. The ASD describes the noise model analytically and is expressed in terms of range rates (rr) [17], where is defined as frequency. The on-board accelerometer senses the linear non-gravitational accelerations and the angular accelerations acting on the satellites with air drag being the main contributor. The ACC is modeled after [46].

NGGM noise
In the case of NGGM noise, the principal measurement is observed by the laser interferometry instrument, meaning an improvement of factor 100 Additionally, an improved accelerometer is implemented with being the along-track, is across-track, and is the quasi-radial component. The error assumption was provided by the consultancy support of Thales Alenia Space Italia (TAS-I).

Results
We calculated eight seven-day solutions for January and February 2002 with the numerical closed-loop simulation software described in Chapter 3 for the double (estimated until d/o 70) and triple (estimated until d/o 90) pairs defined in Table 2. All scenarios were processed with the Wiese approach, which co-parametrizes low-resolution daily gravity field solutions and longer-term gravity field solution. For the evaluation of long-term applications, the time series was extended to one year for selected scenarios and processed with the NRT approach. In order to summarize the results and to simplify the corresponding figures, the shown degree root mean square (RMS) errors always represent the average degree RMS of all estimated seven-day solutions.
We compared the results in the spatial as well as the frequency domain. The gravity field was transformed to equivalent water heights (EWH). In terms of degree RMS signal/errors, this means where and represent the average density of water and Earth, a the semi-major axis of the Earth, is the Love numbers, and and represent the SH coefficients or coefficient differences to the true solution [47,48].

Double vs. triple pairs
In this section, we compare the achievable performance of double pair and triple-pair scenarios. Figure 3 shows degree RMS curves for the double pair (estimated up to d/o 70) and the triple-pair (estimated up to d/o 90) scenarios as defined in Table 2 in terms of an average of the seven-day solutions of the two-monthly observation period. Daily spherical harmonic coefficients (Wiese parameters) up to d/o 15 were co-estimated, as recommended by [28]. The double pair scenarios 2pi in blue is shown as a reference. The triple-pair scenarios 3pip with the third pair being in a polar orbit (red and green) exhibited small improvements, especially in the lower degrees. Additionally, the crossover between signal and errors degree RMS curves was approximately 5 degrees higher than for the reference double pair scenario. The best performance is visible for scenario 3pii,b with the third pair in an inclined orbit and also laser instrumentation for the SST (lilac curve). With this constellation, the performance in terms of spatial resolution can be improved by 10 degrees compared to the double pair.  Obviously, the main gain in performance of the triple-pairs scenario 3pii including two inclined pairs occurred in the sectorial and near-sectorial coefficients, expressing the improved de-striping capabilities of this constellation. Figure 4b shows the spatial representation of these differences of the recovered and the true solution in terms of EWH difference grids up to d/o 50 due to the crossover of the signal and error curve in Figure 3. The improved performance of scenario 3pii can be seen especially in the spatial representation in Figure 4b. Compared to the other scenarios, the striping effect was reduced, but still visible. In Table 4, the cumulative errors in centimeters EWH of selected upper maximum d/o 10 to 50 are provided. While a third polar pair showed only small improvements of about 5% to 10% compared to the Bender double pair scenario 2pi, a second inclined pair resulted in significantly improvements of about 20% to 40%. The expected improvement of the square-root-n rule led to an expected improvement of 22%, reached by the second inclined pair.   Spatial representations of the improvements of the triple pairs compared to the double pair solution are shown in Figure 5, which depicts the difference of the estimated to the reference solution of the best solution. For this, we analyzed an EWH difference grid up to d/o 50 for a seven-day solution and computed the difference of the estimated and reference field for every point of the grid. Blue means the double pair scenario 2pi had the smaller difference value, thus was closer to the reference and visualized in the figure, while red means the respective triple-pair scenario 3pip,b or 3pii,b had the smaller difference value and was closer to the truth and visualized. When comparing scenario 3pip,b with the double pair (Figure 5a), a clear improvement of the triple-pair scenario was not apparent. While in Canada, the triple pair had a lower difference value, the double pair performed better over Antarctica. Overall, the triple pair had lower difference values in the Northern hemisphere, while in the Southern part, an ambiguous picture was visible.
Scenario 3pii,b compared to the double pair ( Figure 5b) showed a significant improvement for the triple-pair scenario. The double pair was only better than the triple-pair scenario in a few limited areas. Also, the overall magnitude of the difference decreased. The anisotropic error behavior of a classic GRACE-like near-polar pair was further reduced by the second inclined pair. Such a constellation came, however, with the drawback of only one near-polar pair, which was the basis of every temporal gravity field mission so far. The near-polar pair is important as the ground track has the best global coverage. If the polar pair failed, the polar areas would not be observed at all by the remaining constellation. Within the processing, a regularization would have to be applied if global base functions such as SH are used for parameterization.
Although both triple-pair scenarios had the same amount of observations, the third inclined pair's along-track observations added more information, further reducing the variability between −60° and 60° latitude. The variability decreased in scenario 3pii,b, and not only in the area covered by the inclined pairs, but also in the polar region covered by only one near-polar pair. The result shows that the introduction of the third pair in an inclined orbit in scenario 3pii,b adds relevant signal and stabilizes the solution overall. In Figure 6, the standard deviation of the difference between the simulated and the reference solutions over the whole computation period (one year), furthermore called RMS variability, of both triple-pair scenarios is visualized. For each grid point, the RMS variability value over the time period is shown. The higher the RMS variability, the darker red the grid point is colored. A detail of the Pacific Ocean is selected as this is an area of little to no temporal gravity signal. The result again shows a decreased striping pattern in scenario 3pii,b in Figure 6b compared to scenario 3pip,b in Figure 6a. The reduced striping pattern is therefore a result of the improved instrument noise and the better observation geometry. The RMS variability is reduced significantly from 2.24 cm EWH to 1.7 cm EWH globally for scenario 3pii,b, supporting the previous statement.

Higher-degree daily resolution
As indicated, the most significant advantage of an additional pair is the expected improved performance of the daily solution. We chose the best-performing triple-pair scenario 3pii,b as the test scenario. The maximum daily degree and order varied between d/o 10 and d/o 40 with a step size of 5. Figure 8a shows the degree RMS of the overall solution, and Figure 8b shows the daily solution. In the overall solution, the nominal processing (without co-estimating daily gravity field parameters) is also included in red.  Table 5, we used the AOHIS signal up to the respective SH degree of the solution as a reference, while the full AOHIS signal up to d/o 40 was used for the results shown in Table 6. Therefore, the first case provides only the commissioning error of the solutions, while the latter case includes also unresolved signals (omission error). The cumulative error of the daily solution estimated up to d/o 40 has the worst performance starting at d/o 25 (values marked with "2"). Considering the processing time, which increases with higher daily parametrization, and the performance of the daily solutions, the maximum daily d/o resolution should be chosen using the crossing point of the signal and error curves as a guideline.  Table 5. Table 5. Daily cumulative error of different daily solutions from scenario 3pii,b excluding omission error in cm EWH. Compare to Table 6. The best values are marked with "1". The values marked with "2" implicate an over-parametrization, impacting the result negatively.  To visualize the spatial difference of daily solutions from scenario 3pii,b with different maximum d/o, Figure 9 shows the daily temporal gravity field, as well the difference of daily solutions with different maximum d/o to the reference, for the whole Earth as well as exemplarily for a detail in the Northern Hemisphere over Europe and Africa. The reference gravity field is always used up to d/o 40 in order to include also the non-resolved signals. The increasing signal as well as the emerging error characteristics can be easily tracked among the different resolutions. The figure visualizes the importance of higher spatial resolution, as the details of the temporal gravity field are revealed. The last two rows with d/o 35 and 40 show an increased error, corresponding to the daily degree RMS in Figure 8b. The maximum daily d/o resolution of 40 demonstrates quite clearly the problems that emerge when the daily solution is over-parametrized. Additional constraints would be necessary to handle the arising spatial aliasing in the solution. Additionally, the solution pushes the numerical stability of the overall system to its limit and would potentially need additional constraints or a regularization scheme. Figure 10 visualizes in the spatial domain the error of the gravity field combined with the ground track pattern for the daily solution up to d/o 30 and 40. As the ground track is designed to cover the Earth within a seven-day period optimally, the daily ground track coverage is not ideal. In Figure  10a,b, a daily solution parametrized up to d/o 30 is visualized and shows only a subtle striping pattern. In comparison, the parametrization up to d/o 40 in Figure 10c,d, shows the spatial aliasing, also very prominently visible in the detail picture of the Northern hemisphere in Figure 9. The error is especially dominant in areas that are not covered by ground tracks. Overall, the example shows that an optimal choice of the daily maximum degree is of great importance to the resulting daily gravity field.    Table 7 with signals in hydrology and ice. The catchment and basins borders are used as boundaries. To cover various scenarios, five large catchments, five small catchments, and two ice drainage basins were selected. Figure 11 shows a time series in EWH for the Amazon region (a), the Danube catchment (b), a small part of California (c), and of the southwest (SW) of Greenland (d) over a whole year. As a reference, the AOHIS from the ESA ESM processed the same way is visualized in the dashed blue line. In the Amazon time series, the mass increase due to the rain period is detectable. The small catchment of the San Joaquin River in Southern California covers part of the drought area, as 2002 was one of the driest seasons since records are taken. As the catchment is rather small, an increased error is visible. Table 7 summarizes the standard deviation of the differences of the solutions to their reference of the selected catchments and basins and of the three scenarios. Also, the gain of the triple-pair scenarios compared to the double pair scenario is determined. The gain is the standard deviation difference with the standard deviation of the double pair, multiplied by 100. It shows that the second polar pair gave a smaller benefit in average 10.7%, while the second inclined pair improved on average by 21%. One of the examples of Greenland shows that the gain is even visible in areas, where the inclined pair does not add observations. When only evaluating the big catchment, both triple-pair scenarios perform similar with more than 20% improvement. When investigating the small catchments, however, the performance of the triple-pair scenario 3pip,b declines to a small improvement only of 5%, while scenario 3pii,b has the same improvement as for the big catchments (cf. Table 7). In the ice drainage basins, triple pair 3pip,b performs even worse than the double pair.  Based on scenario 3pii,b, a scenario with daily solution co-estimated up to d/o 25 was also computed, to analyze the possibility of a higher daily solution. Figure 12 therefore visualizes the daily solutions up to d/o 15 and 25 from scenario 3pii,b for the Danube catchment and the ice drainage basin in SW Greenland over two months. The corresponding standard deviations can be found in Table 7.

Applications
In both examples, a difference between the d/o 15 and 25 curve is visible and can be explained by the omission error. Figure 13 compares the seven-day solution to the d/o 25 daily solution. In the Greenland example (b), both curves have the same base characteristics. However, in the daily curve, the omission error is visible, as the curve is missing signal amplitude. The Danube time series (a) in comparison visualizes the fact that the daily solution has the same signal strength and adds additional time-variable information. Both examples in Figure 13 show that both daily and weekly solutions have advantages, depending on the characteristic of the underlying geophysical signals. Based on the individual applications demand, the spatial and temporal resolutions have to be balanced.  In Figure 14, the magnitude of the daily signal of the SW Greenland example due to the omission error was improved by adding a three-day solution (computed based on the same scenario) starting at the highest d/o of the daily solution (see [28]). This means up to d/o 25 the daily solution was used, artificially enhanced by the three-day solution starting d/o 26 to the maximum degree of 50. As [28] showed, this approach does not significantly distort the daily solution, while at the same time, it adds information with long-wavelength frequency content. Part of the omission error, therefore, can be reduced.

Conclusions
With our closed-loop simulation software, we tested the impact of a third satellite in an inclined or near-polar orbit with GRACE-like or NGGM instrument noise specifications. As the base scenario, we used a Bender-type double-pair constellation. While the second pair, as demonstrated by, e.g., [26], has the added benefit of enhanced de-striping capabilities and further allows the observation of the full AOHIS, we showed that the impact of an additional third satellite pair is between 5% and 40%. Compared with the expected 22% improvement of the square-root-n rule, this indicates that the additional polar pair does add redundant information, while the third pair in the inclined orbit adds significant observations. Our tests also revealed a dampened striping effect and an improved error RMS variability. The third pair in an inclined orbit is preferable compared to the third pair in a polar orbit. The best third pair option is an inclined pair with improved NGGM instrumentation. Therefore, we recommend flying an inclined pair, to add not only additional observations, but also a different observation geometry. The error curve of the daily triple pair gravity field meets the signal curve at approximately degree 26. Our tests confirmed that the processing up to d/o 30 is beneficial, while a higher parametrization adds spatial aliasing to the solution. We therefore note that it is crucial to choose the daily maximum degree with care.
To test the best triple-pair scenarios for applications, we extended the time series towards one year. We looked into typical regions regarding hydrology and ice like the Amazon catchment, the Danube basin, or Greenland. Compared to the Bender double pair scenario, the gain of the third pair in a polar orbit is on average 11%, while adding a third pair in an inclined orbit exhibits a gain of approximately 21%, which is achievable for bigger as well as small catchments. Especially the visualized SW Greenland time series demonstrates the importance of a high d/o daily solution, as the omission error in the region can be significant. The tested approach to add higher d/o long-term solutions onto daily solutions as presented in [28] improves the results significantly and can be recommended as methodology for future applications.
As a conclusion, while the benefit of a Bender-type double-pair constellation over a single-pair GRACE like concept is substantial, we can state that a third pair improves the achievable performance of a double pair constellation only moderately. Our results also suggest that an inclined third pair would be the better option to improve the overall performance.