Identifying effective sink node combinations in spacecraft data transfer networks

Complex networks are emerging in low-Earth-orbit as the communication architectures of inter-linked space systems. These data transfer networks vary based on spacecraft interaction with targets and ground stations, which respectively represent source and sink nodes for data flowing through the network. We demonstrate how networks can be used to identify effective sink node selections that in combination provide source coverage, high data throughput, and low latency connections for intermittently connected, store-and-forward space systems. The challenge in this work is to account for the changing data transfer network that varies significantly depending on the ground stations selected—given a system where data is downlinked by spacecraft at the first opportunity. Therefore, passed-on networks are created to capture the redistribution of data following a sink node’s removal from the system, a problem of relevance to traffic management in a variety of flow network applications. Modelling the system using consensus dynamics, enables sink node selections to be evaluated in terms of their source coverage and data throughput. While restrictions in the depth of propagation when defining passed-on networks, ensures the optimisation implicitly rewards lower latency connections. This is a beneficial by-product for both space system design and store-and-forward data networks in general. The passed-on networks also provide an insight into the relationship between sink nodes, with eigenvector embedding-based communities identifying sink node divisions that correspond with differences in source node coverage.

highly combinatorial problem. Such an approach opens up the potential for agile and responsive space systems that can be adapted by altering their connectivity to the ground, rather than relying on costly and limited spacecraft manoeuvring capabilities. While the developed approach is shown to be effective for the specific challenge of space system analysis, it is expected that the presented methods could be effectively applied to a range of similar problems seen in, for example, traffic flow systems (see Nath and Dhamala 2018) and wireless sensor networks (see Kim et al. 2005;Safa et al. 2014).
Data transfer is a spreading process that Clark et al. (2019) showed can be represented by a network in order to detect the relative influence of nodes. A network that uses aggregated contacts over time to weight edges, enables the network's adjacency matrix to provide insights into the major pathways for spread, as demonstrated by Clark and Macdonald (2021) for identifying influential disease spreaders in contact networks. For space system flow networks, where targets are sources of data and ground stations are sinks, Clark et al. (2022) detailed how the eigenvectors of the adjacency matrix can reveal the relative influence of ground stations in terms of receiving target data. However, the aggregation of contact times, to approximate data transfer as in Clark et al. (2022), limits the applicability of the approach to a system dealing with the transfer of discrete data packets-as is the case in many applications including Earth observation and Internet of Things (IoT) services. Since the order in which a spacecraft comes into contact with targets and ground stations plays an important role in determining system performance. To address this challenge, we go beyond the work presented in Clark et al. (2022) by proposing an aggregated network that accounts for the temporal ordering of contacts. This includes the redistribution of data when a ground station is removed from the system, a necessary step in evaluating an effective subset selection from a set of candidate ground station locations. The redistribution of data provides an estimate for where data will go if a ground station is no longer included in the system, which impacts the relative influence of each sink node within the network.
The ground station selection problem is highly combinatorial and with an objective that varies depending on the application. A common objective for data transferring space systems is the reduction of latency; the delay from a spacecraft acquiring data to the receipt of that data on the ground (see Mazzarella et al. 2020). Alongside latency, target coverage is an important consideration. Targets are defined here as locations on Earth from which data is collected, whether through communication or some other form of sensor acquisition. A network-based approach for ground station selection is proposed herein, which avoids the need to evaluate the full range of feasible ground station (sink node) selections; this is often an intractable problem, particularly when each assessment requires a detailed and time-consuming data transfer simulation. The network representation proposed herein can be used to explicitly optimise the sink node selection in terms of source node (target) coverage, while the optimisation implicitly rewards lower latency solutions in its estimation of the data transfer network.
The majority of systematic ground station selection papers, to date, have focused on large, latency-prioritising constellations that maintain continuous contact between targets and ground stations (referred to as bent-pipe systems). Examples of these systems include OneWeb and Starlink, where target-ground station geographical proximity del Portillo et al. (2019), Chen et al. (2021) has been shown to drive ground station placement and minimum cost, maximum flow optimisation has been used to define effective inter-satellite link topologies del Portillo et al. (2019). For many other applications involving data collection, latency is an important but not singular goal of the constellation. Such as store-and-forward systems-where spacecraft gather information from one location (e.g. ship AIS beacons or Earth monitoring images) and deliver it to another surface location (referred to as a ground station)-that are the focus of this paper as ground station placement must account for latency, target coverage and data throughput. Additionally, in contrast to fully interconnected bent-pipe systems, the order of connections in the temporally varying topology of the storeand-forward contact network must be considered when determining effective sink node selections.
In the past, ground station network design has relied heavily on engineering judgement and best practices. Lacoste et al. (2011) demonstrated the difficulties in applying best-practices for selecting multiple ground stations. They found it difficult to predict the contribution of an additional ground station to an existing set, highlighting the need for combinatorial optimisation methods for the ground station selection problem. An optimised selection of ground station has been proposed by Capelle et al. (2019) for a spacecraft with free-space optical communication, which has communication restricted by cloud cover. The optimisation objective in Capelle et al. (2019) aims to maximise the percentage of data acquired from a single spacecraft. This differs from the target-centric multi-spacecraft problem presented herein, but it does highlight the combinatorial optimisation challenges of the problem and presents both an exhaustive enumeration, similar to that described herein, and a branch-and-bound approach to identify effective subset selections. Tailoring a ground station selection to a system's priorities is attractive both as a cost saving measure and as a means to achieve a robust and adaptable system without having to alter the assets in space. With services offering leaseable ground station sites around the globe, this paper presents an approach for space system designers to maximise constellation performance as mission objectives and target priorities change.

Methods
This section describes the pipeline for identifying an effective subset selection of sink nodes. The steps involved are as follows: • Propagate the movements of all spacecraft in the Space System Scenarios to create a contact schedule (C).
• Generating data transfer networks, including a data transfer network ( ) and a passed-on network (B[g]) for every ground station g ∈ G . These networks combine to produce an estimated data transfer network (A) for a given subset selection of ground stations. • Identify an Initial eigenvector-based selection of ground stations (sink nodes) using an eigenvector embedding of a ground station relationship network ( Ŵ). • Perform an Exhaustive search optimisation based on Consensus dynamics for target coverage, where the objective is to rapidly drive source (target) nodes to consensus under the influence of sink nodes (see Problem definition).

Space system scenarios
The space system studied, and the simulation used to create a contact schedule (C), are described in more detail by Clark et al. (2022), but the relevant aspects are summarised here. The space system considered is based on the orbital positions and targets of the Spire Global, Inc. constellation that collects AIS data from ships globally. All 111 spacecraft that as of July 2021 were operated by Spire Global, Inc. are included in this case study, with their Keplerian orbit elements detailed in data set (McGrath and Clark 2021). The spacecraft are in differing orbit planes with 74 in sun-synchronous orbits, 22 at approximately 51.6 degrees inclination, 8 at approximately 37 degrees inclination, 4 in near-polar orbits, and 3 in near-equatorial orbits. A representative set of target locations are defined for the case studies, based on data provided by Spire Global, Inc. for the 24-h period of 11-August-2019 14:09 UTC to 12-August-2019 14:08 UTC. This dataset provides the last reported position of all ships detected in this 24-h window. From this, 250 targets are positioned to approximate the locations of ships worldwide that were tracked from space (rather than via ground-based coastal AIS receivers) with these locations visualised in the "Results" section ( Fig. 3). These 250 targets define the global targets scenario, while a sub-set of 16 targets located near the Caribbean are taken as the basis of the Caribbean scenario. Twenty ground station sites are considered in this study, with the locations of these sites also visualised in the "Results" section ( Fig. 3).
A fixed-step integrator is used to propagate the motion of spacecraft for a defined period of time (T) and time step ( τ ) to identify contacts (i.e. visible ground stations or targets on the ground). These contacts are collated in a contact schedule (C), which is used to determine the data transfer networks.

Generating data transfer networks
A data transfer network ( ) is created to capture the data transactions in the space system, with a set of ground targets, spacecraft in given orbits, and a set of candidate ground stations. The network is populated by propagating the satellites' motion and simulating data transfer in the system for a defined period of time (T), during which the movement of data packets is monitored. The process of generating is described in detail in Alg. 1 (in black text), and is summarised as follows: • Each spacecraft in the system is assigned a data buffer (db1), where data is inserted when the spacecraft is in contact with a ground target according to the contact schedule (C). • Each packet of data is associated with the target of origin (d) when inserted into the buffer db1. • When the spacecraft is in contact with a ground station, packets in the buffer db1 are removed until a downlink/packet removal limit ( δ ) for a single time step is reached. • For each data packet removal from db1, the data transfer network ( ) is updated with d,n D +g = d,n D +g + 1 , where the d is the target of origin, g is the current ground station (in contact with the spacecraft), and n D is the number of targets. Therefore, by the end of the simulation d,n D +g will equal the number of packets acquired from d and downlinked to g.
In addition to generating the data transfer network ( ), a passed-on network (B) must be created for each ground station to estimate where data will be transferred if that ground station were removed from the system. This allows the importance of each ground station to be better understood, since not all sink nodes in will be present in the final subset selected. This process is intertwined with the generation of and hence is also detailed in Alg. 1 (in blue text), but can be summarised as follows: • Each spacecraft in the system is given a second data buffer (db2), which is populated with dummy data (0 entries) when in contact with targets (i.e matching the data inserted into db1). • In addition to dummy data, a passed-on data packet [d, g] is inserted into db2 for every data packet d that is removed (i.e. downlinked) from db1 for the same spacecraft, where g is the current ground station (in contact with the spacecraft). • When a spacecraft is in contact with a ground station, the dummy (0 entry) data is the first to be removed from db2 before any of the passed-on data packets associated with ground stations. Only once all the dummy data is removed, then the passed-on data packets are removed from db2. • For each passed-on data packet [d, γ ] removed from db2, while the spacecraft is in contact with ground station g, the entry in the data transfer network where γ identifies the ground station of origin for the passed-on data. Therefore, by the end of the simulation B[γ ] d,n D +g will equal the number of packets that were originally acquired from d, but were passed-on from ground station γ to g. • When a passed-on data packet [d, γ ] is removed from db2, a new packet [d, g] is inserted into db2 that is associated with the current ground station contact (g). The number of times a passed-on data packet is re-inserted back into db2 is restricted by a packet re-insertion limit ( ρ ). The impact of ρ is discussed below.
• As with the removal of data from db1, the downlink limit is monitored for packets removed from db2. However, the count of packets removed and this limit are monitored separately for each ground station of origin ( γ ) for passed-on data. • Note that passed-on data packets are not removed from db2 if their ground station of origin ( γ ) is either the current ground station or in close proximity to the ground station γ (see in Alg. 1). This is necessary to avoid the majority of passed-data packets from travelling back and forth between nearby ground stations.
The packet re-insertion limit ( ρ ) is an important consideration, as this determines the number of times a data packet is passed from one ground station to another. The most accurate passed-on matrices were generated when using a ρ value that is similar to the average number of unselected ground stations that a spacecraft could expect to pass before connects with a selected selection. In this paper we are considering a subset selection of five ground station from a set of 20 candidates, therefore data packets can be estimated to, on average, pass through three ground stations before alighting at a selected station. Given that a significant portion of data packets could pass through more than three ground stations, ρ = 4 was applied.

Estimating data transfer network
The difficulty in identifying effective ground station combinations stems from the impact that one selection has on the value of other ground stations in receiving data and covering targets. For example, a ground station (GS1) may be viewed by a spacecraft that has received data from a target (T1). However, it is possible that the data transfer network ( ) does not report this connection if, for instance, the spacecraft has already downlinked all of T1's data to other ground stations prior to overflight of GS1. Therefore, we propose an approach for estimating the data received by a subset selection of ground stations, using the passed-on networks (B) to identify where data would go if a ground station was removed. This approach has been formulated for the analysis of space systems, but such an approach is generalisable to combinatorial flow network problems, where the removal of a sink node results in greater traffic arriving at other sinks in the network.
To estimate the data received by a subset selection, the data transfer network ( ) defined for the full set of ground stations needs to be updated according to the passedon networks (B). This process creates an estimated data transfer network (A) and is detailed in Alg. 2, with the ground station selection represented by a vector r where r g = 1 indicates a selected ground station g, and r g = 0 denotes an unselected ground station. The process involves moving data from each unselected ground station, in turn, by using the normalised passed-on matrix (K) to determine where the data goes, before removing data from the ground station's column in the data transfer network A, and then redistributing the removed data according to K.
The logic used to determine a suitable packet re-insertion limit ( ρ ) for Alg. 1 is also relevant for selecting a suitable n pass for Alg. 1. The ρ value determines how many times a passed-on data packet is re-inserted into the data buffer (db2), while n pass represents the number of times data is moved on from unselected ground stations when estimating the data transfer network (A). Since data packets can be estimated to pass through, on average, three ground stations before alighting at a selected station, then n pass ≥ 3 could be expected to allow the estimated data transfer network to capture the majority of redistributed data. As will be discussed in the "Results" section, n pass cannot simply be set as a large value to capture all redistribution of data as this can over-estimate the volume of target data received by a subset selection of ground stations.

Consensus dynamics for target coverage
An effective way of evaluating a subset selection of ground stations in terms of target coverage and data throughput, for a network G = (V , E) of targets and ground stations, is through the use of consensus dynamics. Specifically consensus leadership, where ground station selections are identified by assessing their ability to lead targets to consensusaccording to the following consensus protocol-when the connections are defined in the estimated data transfer network (A). We consider a system where each node v i has a state x i ∈ IR and continuous-time integral dynamics, ẋ i [t] = u i [t] where u i ∈ IR is the control input for agent i. The linear consensus protocol is and describes how node v i adjusts its state at time step (t) based on the estimated data transfer matrix ( A = [a ij ] ) and the node state (x) of its neighbours ( N i ). Given this protocol, the state of the network develops according to ẋ[t] = −Lx[t] with the graph Laplacian matrix, L, defined as L = D − A where D =diag(out(v 1 ), . . . ,out(v n )) is a diagonal matrix composed of the outdegrees of each node, i.e. out(v i ) = j a ij .
Given the definitions for the continuous-time integral dynamics and ẋ i [t] , the discretetime agent dynamics are given in Di Cairano et al. (2008) as The choice of ǫ affects the number of steps required for nodes to reach convergence, therefore setting ǫ = 0.999 × 1 max i d ii allows the number of computational steps to be reduced while still guaranteeing convergence of the system (see Di Cairano et al. 2008). Convergence is defined here as x i > 0.99 ∀ i ∈ D , where D is the set of all target (source) nodes, when x j = 1 ∀ j ∈ G with G the set of all ground station (sink) nodes.
The most effective ground station selections, in terms of target coverage, are those that achieve the fewest steps until all of the targets reach consensus. Such a selection would demonstrate a strong connection to all of the targets in the system. If, in contrast, a selection had no connectivity to a given target then consensus would never be reached.

Problem definition
An objective function is required to optimise the ground station selection. The number of steps to convergence can be used, but it creates a discontinuous search space. Therefore, the mean consensus leadership, provides a continuous alternative to maximise the mean consensus state of all target nodes, where n D is the number of targets (source nodes) and D the set of all targets. The target (source) nodes states, x i [t] , are evaluated according to Eq. 2 where t is taken as a point prior to convergence, defined as the closest step to 0.9 × s ref where s ref is the number of steps to convergence. The reference number of steps, s ref , is defined using the number of steps to convergence required for the Initial eigenvector-based selection .
The optimisation can then be defined as follows, where r is the ground station selection vector, is the subset selection of ground station, G is the set of all ground station candidates and n select the cardinality of the subset .

Initial eigenvector-based selection
The optimisation of ground station selections is a highly combinatorial problem and as such susceptible to local optima far from the global optimum. This issue is exacerbated by the need to update the data transfer network for every possible selection. We propose an eigenvector embedding-based selection to act as an effective initial selection, providing an alternative to a more exhaustive search. The use of brute-force evaluation of all combinations is often intractable for sufficiently large numbers of candidates and selection sizes.
The relationship of interest, when optimising a system for target coverage, is that between targets and ground stations. However, it is not possible to directly capture this relationship in a static network. Instead a ground station relationship network ( Ŵ ) is introduced, based on the passed-on networks B, which details the volume of data that each ground station passes on to every other ground station when removed from the system. While the passed-on networks, B[g] ∀ g ∈ G , detail the movement of data from targets to ground station, the network Ŵ details the connections between ground stations, where The Ŵ network is useful in identifying influential ground stations. This is despite Ŵ only detailing the relationships between ground stations, since these relationships are a product of connectivity to spacecraft that have collected target data. Therefore, Ŵ highlights whether ground stations are connected to spacecraft in similar or different orbits. Differing spacecraft orbits result in different target contacts, where these differences lead to different patterns of target coverage. Hence, selecting ground stations that cover different sets of spacecraft will also likely provide a selection that covers differing communities of targets.
The process of ground station selection takes inspiration from work on communities of dynamical influence (CDI), introduced in Clark et al. (2019), that are shown to highlight effective leadership in networks under consensus dynamics. The selection is based on the eigenvectors of Ŵ , where the dominant eigenvectors (those associated with the largest eigenvalue entries) are used to embed the network in a Euclidean space. The nodes in this space that are furthest from the origin, along the direction of their position vector, are defined as leaders of separate ground station communities. This is assessed by comparing the magnitude of each node's position vector with the scalar projection onto this vector from all other node position vectors.
The explicit objective of the optimisation is to improve target coverage and data throughput, therefore CDI analysis of Ŵ can facilitate the selection of ground stations. Specifically, an effective combination of ground stations can be expected to involve nodes in multiple different communities to ensure target coverage, while the nodes with the largest first left eigenvector ( v 1 ) entries are more likely to ensure high data throughput. Therefore, an initial selection composed of ground stations from different CDIs, each with the largest v 1 entry will provide a good initial guess.

Exhaustive search optimisation
An optimal selection of ground stations, in terms of convergence to consensus for the space system modelled using consensus dynamics, can be obtained by simulating all subset combinations from a set of candidates. For the scenarios explored in this paper that involves simulating all combinations of five ground stations from 20 possible options (15503 combinations in total). This is a computationally intensive process that required approximately 10 days (60 seconds per simulation) computation time for the global targets scenario on a desktop machine-Intel Xeon Processor with 12× 3.39 GHz and 46.7 GB RAM. By contrast, using the presented method, an effective selection can be obtained in minutes through the following steps: • A single simulation of data, including all 20 candidate ground stations.
• An initial selection based on eigenvector embedding of the ground station relationship network ( Ŵ). • A simple exhaustive search optimisation, requiring the estimation of data transfer networks as described in Alg. 2.
The simple exhaustive search is described in Alg. 3 and can be summarised as follows: • Identify an initial selection from eigenvector embedding • If necessary, add to the initial selection by performing an exhaustive search for ground stations that minimise the mean consensus leadership (Eq. 3) • Review each selection, in turn, using an exhaustive search until the mean consensus leadership is minimised.

Results
The efficacy of Alg. 3 is demonstrated in Fig. 1, by comparing a set of optimised selections with all possible selection combinations of five ground stations from 20 possible ground station locations (geographical locations shown in Fig. 3). To assess the performance of selections, an individual simulation was completed for each combination detailing the movement of data over a 1-day time period to calculate the average latency (time taken from data acquisition to downlink) and the volume of data delivered from each target to each ground station. The data volumes were then used to assess the number of steps to convergence for targets under consensus dynamics (Eq. 2), where the connection between target and ground station is defined as being equal to the volume of data transferred. A low number of steps to convergence indicates that a ground station selection has a strong data connection to all of the targets in the system (i.e. good target coverage and high data throughput). In Fig. 1 the selections identified by applying Alg. 3 are seen to be near the Pareto front of the search space, with solutions producing both low latency and a low number of steps to convergence. Selections are shown for varying n pass values (the number of data pass iterations, see Alg. 2). For n pass = 0 , this means that the original data transfer network is used without adaptation. The n pass = 0 data transfer network primarily includes all of the ground stations that are the first to be seen after a satellite has collected data from a target. Note that it is possible for proceeding ground stations to also receive data from a target, but this will only occur if the spacecraft collects more data than it can downlink to the first ground station. It is therefore unsurprising that n pass = 0 selections produce some of the lowest latency solutions.
In Fig. 1, the results show how n pass > 0 can reveal selections that provide greater target coverage than n pass = 0 selections. This is to be expected, as the estimated data transfer network will more accurately capture how data is redistributed when ground stations in close proximity to targets are not selected. This allows the optimisation to identify the ground stations that will receive the most data when reducing from 20 to 5 ground stations.
It can also be seen in Fig. 1 that there is variation in the results depending on the n pass value. The n pass value determines the number of data pass iterations when estimating the data transfer network (Alg. 2), where with each iteration data is passed on from unselected ground stations. Therefore with too few iterations insufficient data is passed on to ground stations that would form effective selections. Conversely, too many iterations results in an excess of data being estimated as arriving at poorly connected ground stations. Hence, an n pass value similar to the average number of ground stations located between a target and a selected ground station is recommended. In this case, with 20 ground station locations and only 5 selected an n pass = 4 would be expected to perform best for optimising steps to convergence. However, as demonstrated in Fig. 1c this is not a guarantee given the errors in estimating data transfer and the combinatorial nature of the problem.
To demonstrate the performance of the method with differing initial selections, Fig. 1b shows the results of selecting a set of 5 ground stations from 11 candidate ground stations. These 11 ground stations are a subset of the original set of 20 stations, selected by identifying any station that appeared in any of the n pass = 0 − 7 selections in Fig. 1a. These results show that, despite a reduced search space, the n pass = 0, 1, and 4 selections perform slightly worse in terms of steps to convergence. The pattern remains similar to Fig. 1a, with the n pass = 4 and 7 selections performing notably better in terms of steps to convergence.
A similar pattern to the global targets scenario is also seen in Fig. 1c for the Caribbean targets scenario, whereby n pass = 0 produces a low latency solution with the lowest steps to convergence solution found by increasing n pass to 7. As discussed, an n pass = 4 selection would be expected to facilitate the identification of an effective selection in terms of steps to convergence. However, in this instance it is likely that the localised location of targets has led to improved selection with n pass = 7 . Since the targets are constrained to one geographical location (Caribbean), then to get an accurate picture of the data received from distant ground stations a high number of iterations will be required to pass on data from unselected ground stations. This is less of an issue in the global targets scenario, as most ground stations are selected for their (relatively) local geographical coverage of targets.

Eigenvector embedding
The initial ground station selections, generated from eigenvector embedding of the spacecraft relationship network ( Ŵ ), are altered during the optimisation to produce the results shown in Fig. 1. However, Fig. 2 provides evidence that the communities of dynamical influence (CDI), on which the initial selections are based, identifies communities with differing target contacts that should be covered to enable good target coverage. This is shown in Fig. 2, as the optimised selections cover all four CDI in Fig. 2a and only leave the least prominent CDI in Fig. 2b unrepresented. The same optimised selection can be identified by starting from a randomised initial selection, but eigenvector embedding-based selections reduce the number of exhaustive searches required to find an optimised solution.
Ground station nodes furthest from the origin in Fig. 2 (i.e. large v 1 , v 2 , and v 3 entries) are likely to be in receipt of large volumes of data from other ground stations in the network according to the passed-on matrices. However, the Ŵ matrix is only an estimation of data transfer following ground station removal. Therefore, as shown in Fig. 2b, it is possible for the node with the largest v 1 entry (i.e. eigenvector centrality) to not be included in the optimised selection. This particularly occurs when other nodes in the same community are selected, as that prevents these nodes from passing data on to the most prominent node in their community.

Mapping the results
Examples of effective ground station selections are shown on a world map in Fig. 3, alongside the locations of targets and all candidate ground station locations. The ground station selection for the global targets scenario (blue) forms an evenly distributed cross, which facilitates the selection in achieving an even coverage of global targets through their spacecraft connections. The combination of polar and equatorial locations is important for spacecraft connectivity, where polar ground stations achieve long connection times with the 78 polar spacecraft but cannot be relied upon exclusively as they do not receive data from the 33 other spacecraft in the constellation. The equatorial ground stations, in contrast, are seen by all spacecraft in the constellation, but for, generally, less time. This ground station selection is hence driven by the hybrid nature of the constellation.
The Caribbean scenario also presents a distributed cross formation, but geographically localised targets alter the selection by placing two ground stations in relatively close proximity (in longitude) to North America. This is facilitated through the use of the presented approach, which captures the temporal order of connections and hence  Fig. 1a) and Caribbean (red, n pass = 7 in Fig. 1c) scenarios, as well as other candidate sites encourages the selection of stations in close proximity to the target. With the majority of spacecraft in sun-synchronous or near-polar orbits, ground stations at similar longitudes to the targets will naturally provide low-latency solutions. The equatorial ground station selected in India is of value as the three equatorial orbiting spacecraft will consistently overfly both the Caribbean targets and this ground station. Furthermore, as its location is separated from the Caribbean targets by approximately 180 • , all non-equatorial orbiting spacecraft that view the Caribbean targets on an ascending pass, will view the equatorial ground station on the descending pass, and vice versa.

Conclusions
This paper demonstrates that effective ground station subset selections can be identified, for a given space system, from a single simulation involving the full set of candidate sites. Consensus dynamics provide a useful basis for optimising the selection of ground stations, which can be defined as sink nodes leading a set of source nodes (targets) to consensus. Comparison of how rapidly the source nodes reach consensus provides an objective that promotes the selection of subsets with important properties, namely good target coverage and high data throughput.
The identification of effective sink nodes from a single simulation is viable due to the ability to estimate data transfer networks, for a selection of sink nodes and a given set of source nodes. This estimation relies on analysis of how data is redistributed when a ground station is removed from the system. The restrictions applied-to the number of times data is redistributed (passed-on) when simulating the system-can prevent the optimisation from identifying globally optimal solutions in terms of convergence to consensus. However, these restrictions are desirable for sink node selection in space systems, and store-and-forward data transfer systems in general, as they result in the optimisation implicitly rewarding lower latency connections.
The relationships between sink nodes, in terms of passed-on data redistribution, is key to both estimating the data transfer networks and for gaining insights into effective selections. Insights can be obtained into effective sink node selections, through embedding-based community detection in an eigenvector-defined Euclidean space. Effective selections are distributed across the detected communities, which is to be expected as these communities implicitly capture ground station division in terms of target coverage.

CDI
Communities of dynamical influence