Satellite Systems - Design, Modeling, Simulation and Analysis
Satellite Systems - Design, Modeling, Simulation and Analysis
Compression of High-Resolution
Satellite Images Using Optical
Image Processing
Anirban Patra, Arijit Saha, Debasish Chakraborty
and Kallol Bhattacharya
Abstract
This chapter presents a novel method for compressing satellite imagery using
phase grating to facilitate the optimization of storage space and bandwidth in satel-
lite communication. In this research work, each Satellite image is first modulated
with high grating frequency in a fixed orientation. Due to this modulation, three
spots (spectrum) have been generated. From these three spots, by applying Inverse
Fourier Transform in any one band, we can recover the image. Out of these three
spots, one is center spectrum spot and other spots represent two sidebands. Care
should be taken during the spot selection is to avoid aliasing effect. At the receiving
end, to recover image we use only one spectrum. We have proved that size of the
extracted image is less than the original image. In this way, compression of satellite
image has been performed. To measure quality of the output images, PSNR value
has been calculated and compared this value with previous techniques. As high-res-
olution satellite image contains a lot of information, therefore to get detail informa-
tion from extracted image, compression ratio should be as minimum as possible.
1. Introduction
1
Satellite Systems - Design, Modeling, Simulation and Analysis
visual quality of satellite image is poor. [3] discusses the use of integer wavelet
regression by increasing temporal correlation, which consequently improves the
compression gain. [4] discusses a satellite image compression technique using dis-
crete wavelet transform for noise removal to compress satellite images. [5] has dis-
cussed the use only hardware-based solutions in this lossless compression technique
of X Sat images. [6] has analyzed the use of Discrete Wavelet Transform in their
lossy image compression work and performance of different wavelets for satellite
image compression. [7] has used the conventional Discrete Cosine Transform sys-
tem for lossless image compression. [8] proposes an image compression technique
using multiplexing and encryption by optical grating method [9–12].
In this chapter, we proposed a scheme to compress multiple high-resolution
satellite images by using phase grating. Each image is modulated by applying high
value of spatial frequency and a fixed orientation angles. For each image, multiple
bands have been generated due to modulation which are placed in the same spec-
trum plane (only three bands are clearly visible). The spectrum is encoded and fil-
tered using Gauss filtering. To detect the location of maximum image information,
an intensity graph has been plotted in the decrypted plane. All stored images can be
securely and efficiently retrieved by applying inverse Regional Fourier Transform
operation. This proposed technique is simple and suitable for optimization of stor-
age space and bandwidth in satellite communication.
The chapter is organized as follows:
• Section IV provides result of our research work and PSNR value of the
extracted images
• Section V concludes the paper and discuss why low compression ratio is desired
for land cover analysis.
Images which are used in our research work, collected from Regional Remote
Sensing Centre (East). The images are satellite picture of different areas in Kolkata
Metropolitan Area.
All satellite images used in this paper are captured by LISS III 23 m sensor.
LISS- III sensor is an optical sensor working in four spectral bands (Green, Red,
Near Infrared and Short-Wave Infrared). It covers a 141 km- wide swath with a
resolution of 23 meters in all spectral bands [13, 14].
3. Proposed methodology
According to rule of phase grating, value of the grating frequency (𝑢0) should be
high. Low grating frequency is creating aliasing problem and therefore it would be
very difficult to reconstruct the original image. In this chapter, we propose to select
2
Compression of High-Resolution Satellite Images Using Optical Image Processing
DOI: http://dx.doi.org/10.5772/intechopen.94147
𝑢0 = 1400, which is sufficient for filtering. In grating, value of the orientation angle
(θ) varies from 0 to 3600. We have worked with 0 deg orientation angle.
The diffraction gratings used are illustrated in Figures 1 and 2.
q=+α
m q= +α
m
S=
1 ( u, v ) [ F1 ( u, v ) ∑ q J δ ( u − =
qu1 , v ) ∑ J q F1 ( u − qu1 , v ) (2)
q=
−α 2 q=
−α 2
where J q is the qth order Bessel Function of the 1st kind. The diffraction order is
represented by the parameter ‘q’.
The diffraction pattern, as given by Eq. (2) is a series of diffraction spots each
containing the object spectrum. Considering the zero order and the first two orders
of the spectrum, Eq. (2) may be represented by,
) J0 2 F1 ( u, v) + J−1 2 F1 ( u + u1 , v) + J+1 2 F1 ( u − u1 , v)
m m m
S1 ( u, v= (3)
Figure 1.
Sinusoidal phase grating along x axis.
3
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 2.
Spectrum of an image due to modulation.
fixed cut-off frequency. As only one sideband is enough for image reconstruction,
hence we select one side-band. This spectrum is filtering using a cut-off value of
120. Proper operating frequency is selected and there is no overlapping during the
band selection, so the extracted images are free from aliasing problem. As center
band contains the information of all images, therefore filtering is not applied across
center band.
1 m n
∑∑ f ( x,y ) − b ( x,y )
2
=PSNR (4)
mn =y 1=x 1
Satellite Images chosen for testing the algorithm described in Section III are
shown in Figure 3(a)-(c). The dimensions of the selected images are 512 x 512.
Figure 3.
Three high Resolution Satellite Images.
4
Compression of High-Resolution Satellite Images Using Optical Image Processing
DOI: http://dx.doi.org/10.5772/intechopen.94147
Figure 4.
Three extracted images.
Table 1.
Summary of PSNR calculation.
During zonal filtering operation, images 𝑓1(𝑥, 𝑦), 𝑓2(𝑥, 𝑦) and 𝑓3 (𝑥, 𝑦) have
been extracted by Regional Inverse Fourier Transform taking upper spectrum from
horizontal direction. Extracted images are shown in Figure 4(a)-(c), respectively.
A summary of the PSNR calculation is presented in Table 1.
5. Conclusion
In this chapter, phase grating technique has been proposed for compressing
the high- resolution satellite images in frequency domain. The original image is
retrieved by applying Inverse Fourier Transform from the respective spectrum
of the image. As presented here, since we have taken only few coefficients from
the spectrum, the size of the output image is less than the main original image. It
should be mentioned that in high resolution satellite image, compression should be
as minimum as possible. The main reason for the requirement of low compression
ratio is mainly due to large geographical area representation (as these images con-
tain a lot of information). High compression ratio is not suitable for accurate land
cover analysis. To maintain the same dimension with original image and to avoid
aliasing effect, spectral area is carefully selected. Compared with earlier methods,
visual quality of the selected satellite images is very good as it is captured by optical
LISS-III sensor. Our proposed technique is simple and suitable for optimization of
storage space and bandwidth in satellite communication.
5
Satellite Systems - Design, Modeling, Simulation and Analysis
Author details
© 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License (http://creativecommons.org/licenses/
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
6
Compression of High-Resolution Satellite Images Using Optical Image Processing
DOI: http://dx.doi.org/10.5772/intechopen.94147
References
7
Chapter
Abstract
Recently, the space industry has pointed out that in the past 5 years, the com-
mercial market has been driving the advancement of satellite technology. Lockheed
Martin is building commercial satellites (e.g., Hellas-sat series) with advanced
on-board processing capabilities for the Saudi Arabian [1]. Hellas satellites probably
will be the first commercial HTS with a very advanced digital processor on-board.
The focus of this chapter will be on commercial satellite systems for communication
applications, and a comparison study between commercial HTS and typical satel-
lites systems conducted by Inmarsat will be provided [2].
For communication applications, commercial satellite systems have been catego-
rized as mobile satellite services (MSSs), fixed satellite services (FSSs), broadcast
satellite services (BSSs), and high-throughput satellite (HTS) services. Depending
on the services, satellite payload architecture will be designed to meet the specified
requirements for that service. Basically, satellite payload architecture can be classi-
fied into four categories: (1) analog bent-pipe satellite (ABPS); (2) digital bent-pipe
satellite (DBPS); (3) advanced digital bent-pipe satellite using digital channelizer
and beamformer (AdDBPS-DCB); and (4) advanced regenerative on-board pro-
cessing satellite (AR-OBPS). This chapter provides an overview of these payload
1
Satellite Systems - Design, Modeling, Simulation and Analysis
Typical and regular commercial satellites are operating in C-band, Ku-band, and
Ka-band with downlink frequencies approximately at 4, 12, and 40 GHz, respec-
tively. For C-band, Ku-band, and Ka-band, the spectrum bandwidths available by
geostationary orbital position are 500 MHz, 500 MHz, and 3.5 GHz, respectively.
Typical antenna types for these regular commercial satellites are pointed antenna
type with a single beam. Typical diameters for these pointed antennas are (a)
greater than 1.8 m for C-band; (b) 0.9–1.2 m for Ku-band; and (c) 0.6–1.2 m for
Ka-band satellite. Figure 1(a) illustrates a typical regular commercial satellite.
Typical HTSs are usually also operating in Ku-band and Ka-band with the same
downlink frequencies as the regular satellites except that they employ multiple
pointed beam as oppose to a single-pointed beam. Figure 1(b) describes a multiple
beam HTS system. The salient feature of multiple beams is the frequency reuse.
The frequency reuse is defined as the number of times a satellite can reuse the
same spectrum and frequencies. However, high frequency reuse factor can cause
potential cochannel interference or an increase in carrier-to-interference power
ratio (CIR or C/I). IMMARSAT has reported that a reuse factor of 5–30 is possible
with multiple spot beams employed by commercial HTS. Depending on the number
of beams implemented on-board of the satellite, the cost for HTS can be twice of
the cost for a regular satellite. But, the cost per bit for HTS is much lower than the
regular satellite. HTS is a preferred option for point-to-point services, for example,
beyond line-of-sight (BLOS) cellular phone services. Table 1 provides a summary
of the comparison of HTS and regular commercial satellites [2].
Figure 1.
Typical commercial satellites and HTS configurations.
2
Overview of Existing and Future Advanced Satellite Systems
DOI: http://dx.doi.org/10.5772/intechopen.93227
Table 1.
Comparison of typical commercial satellites and HTS.
This section describes the most commonly used satellite network topolo-
gies, namely “Star” satellite network (Section 3.1) and “Mesh” satellite network
(Section 3.2).
The “mesh” satellite network topology is derived from a local network topol-
ogy, where the network nodes are corrected to each other directly, dynamically,
and nonhierarchically to as many other nodes as possible [4]. In this network
topology, the network nodes can cooperate with one another to route data from
one user to another user efficiently. Hence, for mesh satellite network, Mobile
User 1 can talk to fixed user directly without going through the satellite gateway
(solid lines), and Mobile User 2 can also talk to the fixed user directly (dash lines).
3
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 2.
Typical “star” satellite network.
Figure 3.
Typical “mesh” satellite network.
Any one of the user within the network can send the messages to a terrestrial
network through the red lines representing uplink and downlink between the
satellite gateway and the satellite (Figure 3).
Star satellite network topology does not require advanced satellite payload pro-
cessing on-board and multiple beam, but mesh satellite network requires advanced
on-board processing and multiple beam allowing one user to communicate to
another user automatically and effectively. Section 4 discusses various satellite
payload architectures used in regular satellite and HTS for star and mesh satellite
network applications.
4
Overview of Existing and Future Advanced Satellite Systems
DOI: http://dx.doi.org/10.5772/intechopen.93227
Section 4.3 discusses AdDBPS-DCB payload architecture, and Section 4.4 provides
an overview of AR-OBPS payload architecture.
1
Flicker noise is a type of electronic noise with a 1/frequency power spectral density.
5
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 4.
Legacy ABPS payload architecture.
Figure 5.
Options for RF downconversion and associated LO’s phase noise.
Figure 6.
Functional block diagrams of MUX and DEMUX.
digital quadrature technology [6]. The RF bandwidth (BW) associated with the
RF bandpass filter (BPF) is selected to match with an over channel bandwidth
(e.g., a maximum of 500 MHz for Ku-band). The automated gain control (AGC)
6
Overview of Existing and Future Advanced Satellite Systems
DOI: http://dx.doi.org/10.5772/intechopen.93227
Figure 7.
Existing DBPS payload architecture.
Figure 8.
Typical R/F downconversion and digitization processing approach.
7
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 9.
Existing digitization and sampling processing using bandpass sampling with digital quadrature technique.
Figure 10.
Optimum LF as a function of number of bit of A/D converter.
the optimum LF is about 0.4. In conjunction with LF, the number of bit should
be selected to maximize the signal-to-quantization noise ratio (SQNR) using the
following relationship:2
2
Quantization (signal processing). Available from: https://en.wikipedia.org/wiki/
Quantization_(signal_processing).
8
Overview of Existing and Future Advanced Satellite Systems
DOI: http://dx.doi.org/10.5772/intechopen.93227
Figure 11.
Typical digital channelizer using polyphase/DFT technology.
Figure 12.
Typical Polyphase/DFT Technology.
be introduced, thus allowing the circuit in different paths of the polyphase network
to operate at lower frequency than the original sampling frequency. A practical
implementation of a high-throughput low-latency polyphase channelizer can be
found in [8, 9].
Figure 12 shows an example of five input signals, namely S1, S2, S3, S4, and S5,
and the channelizer will select signal interest by filtering out the other signals. As an
example, the signal line with the filter transfer function of H0 filters out S2, S3, S4,
and S5 and sends S1 as an output signal.
9
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 13.
AdDBPS-DCB payload architecture.
Figure 14.
DCB Approach 1: PolyN/DFT and DBFN individual processing.
10
Overview of Existing and Future Advanced Satellite Systems
DOI: http://dx.doi.org/10.5772/intechopen.93227
Figure 15.
DCB Approach 2: Unified and combined PolyN/DFT and DBFN individual processing.
frequency of multiplier of 256 MHz. Using these values, D for the DBF/channelizer
of the DCB Approach 1 configuration becomes [13]:
11
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 16.
Antenna beamwidth and gain of a notional Ka-band DBFN with 12-bit quantization [14].
Figure 17.
AR-OBPS payload architecture.
12
Overview of Existing and Future Advanced Satellite Systems
DOI: http://dx.doi.org/10.5772/intechopen.93227
Sections 5.1 and 5.2 present a notional satellite system architecture using AdBPS-
DCBS satellite payload and AR-OBPS satellite system architecture for 5G cellular
phone applications, respectively.
13
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 18.
AdDBPS-DCB satellite system architectures for supporting 5G users.
Figure 19.
AR-OBPS satellite system architectures for supporting 5G users.
14
Overview of Existing and Future Advanced Satellite Systems
DOI: http://dx.doi.org/10.5772/intechopen.93227
6. Conclusion
Conflict of interest
The preparation of this chapter was not funded by Gulfstream, and it was done
by the author using his own time and resources; thus, it does not represent the
Gulfstream’s view on the results presented in this chapter.
Notes/Thanks/Other declarations
The author wishes to thank his wife, Annie Luu-Nguyen, for her immense
patience and support.
Author details
John Nguyen
JohnDTN Consulting Services, Huntington Beach, California, USA
© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License (http://creativecommons.org/licenses/
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
15
Satellite Systems - Design, Modeling, Simulation and Analysis
References
16
Chapter
Abstract
1. Introduction
1
Satellite Systems - Design, Modeling, Simulation and Analysis
were built and launched under the different philosophy that satellite is the main
“personage” performing a space mission and the payload (unlikely the ballistic
rocket war head (s)) is just one of the satellite subsystems that should be integrated
into the satellite board under the satellite chief designer guidance, who is responsi-
ble for the mission performance. From the author’s point of view, this approach has
certain advantages following from the Aerospace System Engineering, integration
and distribution functions, and responsibilities between the space mission partici-
pants. In this chapter, SCS is presented from this point of view, integrating con-
ventionally separate satellite GN and C subsystems and devices into the joint
integrated system, attitude and orbit determination and control system (AODCS).
The main principles and features of this system are presented in this chapter.
• Mass 83.6 kg; sealed from two identical hemispheres with a diameter of 0.58 m;
life time 3 months; payload, two 1 W transmitters (HF, 20.005 and VHF,
40.002 MHz) with four unidirectional deployable antennas (four 2.4–2.9 m
metallic rods); electrical batteries, silver-zinc; sufficient for 2 weeks.
• Orbit: perigee 215 km, apogee 939 km, period 96.2 min, eccentricity 0.05,
inclination angle 65.10 deg.
• Inside, the satellite sphere was filled by nitrogen, and the temperature was kept
within 20–23 deg. C with automatic thermoregulation-ventilation system
(thermometer-ventilator).
• The satellite had no attitude control and was free rotated around its center of
mass in orbit, keeping initial angular speed, provided by the separation pulse
after the separation from the launch rocket. However, thanks to the four rod
Figure 1.
Soviet designers-creators of the first earth-orbiting artificial satellite SS-1.
2
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
Figure 2.
SS-1, assembled (left). Open two semispheres (right).
Since SS-1, about 8378 satellites were launched to year 2018 [6]. Early satellite
launches were extraordinary events and demonstrated tremendous achievement of
the launched state, the USSR (4 Oct. 1957, SS-1), the USA (31 Jan. 1958, Explorer 1)
and Canada (29 Sep. 1962, Alouette, launched by Thor-Agena, a US two-stage
rocket), but with time, satellite launches became ordinary and usually pursue a
certain military or civil mission.
Among the civil missions (satellites), the following types can be determined as
already conventional: navigation, communication, Earth observation, scientific,
geophysics and geodetic, technology demonstration and developers training. These
satellites are usually equipped with a kind of payload system(s) (radio/TV trans-
mitter/transducer, radar, telescope or different scientific instrument, etc.) to per-
form certain dedicated space mission(s). For example, the first Canadian Earth
observation satellite RADARSAT-1 (Nov 4, 1995–May 10, 2013; Figure 3) was
equipped with a side-looking synthetic aperture radar (SAR) on board the Interna-
tional Space Station (November 1998, ISS; Figure 4) was installed a Canadian
robotic arm for its assembling and maintenance.
Figure 3.
The first Canadian earth observation satellite RADARSAT-1.
According to the satellite altitude (h), their orbits can be classified as low-
altitude (LEO), 200–2000 km; medium-altitude (MEO), 5000–20,000 km; and
high-altitude (HEO), h > 20,000 km; according to eccentricity as: close to circular
e < 0.01; elliptical 0.01 < e < 0.3; highly elliptical 0.3 < e < 0.8.
There are satellites with special type of orbit such as polar (i = 90 deg), equatorial
geostationary (GEO, i = 0 and h = 35,800 km) and Sun-synchronous provide orbital
precession equal to Sun annual rate (i depends on satellite period) (Figure 5).
3
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 4.
International Space Station (ISS).
Figure 5.
Satellite orbit types (“tundra” and “Molniya” are Russian communication satellites in highly ecliptic orbits).
4
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
Today, for many satellites, GN&C onboard equipment can be presented by the
following subsystems, performing related functions listed below:
• Sensors
• Actuators
Figure 6.
Satellite AODCS system.
5
Satellite Systems - Design, Modeling, Simulation and Analysis
Depending on required reliability and life time, each component can be a single
or redundant unit. Unlike airplanes, satellite is an inhabitant space vehicle that is
operated from the ground. The operation is usually performed via a bidirectional
telemetry radio link (TLM) in S-band (2.0–2.2 GHz). Payload data downlink radio
link (unidirectional) is usually performed via X-band (7.25–7.75 GHz;). For both
links, usually the same data protocol standards are applied Figure 7.
Figure 7.
Satellite communications with ground stations.
Below AODCS components are presented to show their generic principles that
can help for the system understanding and modeling. Generic design requirements
are presented in [3]. Some design examples can be found in many sources, for
example, [1, 9–12].
3.2.1 Sensors
AODCS sensors are designed to measure satellite orbital and attitude position
and velocity. From the most general point of view, they can be considered as the
6
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
vector measuring devices (VMD). The device can measure in space a physical
vector Rm that can be known (referenced) in a reference coordinate system Rr .
Three parameters can be measured: vector module R and two angles of its
orientation Az and El (Figure 8).
Figure 8.
Vector R in the Cartesian coordinate system XYZ.
It can be noted that measurement of referenced vectors can be used for the
determination of satellite position or angular orientation. A minimum of three vectors
is required to determine satellite position and two to determine its attitude. If more
vectors are measured providing informational redundancy, then such statistical esti-
mation methods as least square method (LSM) and Kalman filter (KF) can be applied.
Satellite velocity and angular rate can be derived by the differentiation of its position
and attitude applying a kind of filter recommended by the filtering and estimation
theory [13–15]. It should also be noted that if vector orientation is measured for the
position determination, then satellite attitude should be known and vice versa.
Especial autonomous satellite navigation system (sensor) is the inertial navigation
system (INS/inertial measurement unit (IMU)). It can be used for the determination
of satellite position, velocity, orientation and angular rate simultaneously. INS is
based on measuring with linear accelerometers and angular rate sensors (“gyros”) the
two vectors: satellite linear active acceleration a and angular rate ω. After integration,
the system provides satellite position, velocity, attitude and angular rate. It is also
assumed in INS theory that the vector of Earth gravity acceleration g is not measured
by the system accelerometers, but it is computed from referenced mathematical
Earth gravity field model. Essential INS disadvantage is that its errors grow with time.
That is why, it has to be periodically corrected by such navigation aids as a pair of
VMD used for the direct attitude determination. A detailed system description is out
of this chapter’s scope and can be found in many publications [16–18]. Only the use of
7
Satellite Systems - Design, Modeling, Simulation and Analysis
angular rate sensors (“gyros”) for determination of satellite attitude is briefly consid-
ered below.
Today, satellite GPS can provide onboard accurate data about position, velocity
and time [19] (Figure 9).
Accuracy: position, 15 m (2σ); velocity, 1.5 m/s (2σ); time, 1 μS.
Figure 9.
Satellite GPS SRG-10. Double redundant with a pair of zenith and nadir antennas.
GPS receiver is a radio range measuring device that measures distance from the
desired satellite to navigation satellite constellation (NAVSTAR, USA; GLONASS,
Russia; and GALILEO, Europe) and computes its position and velocity. GPS mea-
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
sures the distance R (R ¼ R2x þ R2y þ R2z ) of the vector from the desired satellite to
the navigation satellite, and this system is invariant of the system orientation (sat-
ellite attitude). The distance between the desired satellite and navigation satellite is
measured by measuring the time delay Δt between the time ts of the radio pulse
transmitted by navigation satellite and the reception time tr of its reception by GPS
receiver installed on the desired satellite Δt ¼ tr ts . Measuring the distance allows
to determine the desired satellite relative position (relatively to navigation satellite),
and using known navigation satellite position that is continuously received by the
receiver for every tracking satellite in the navigation message (NM) converts it in
absolute position.
A minimum of three navigation satellites should be simultaneously traced by the
receiver to determine position and velocity. Then satellite position is the cross-point
of three spherical surfaces of the position equation Ri ¼ const, i ¼ 1, 2, 3. If more
tracked satellites are available, then redundant information can be used to calibrate
the onboard clock (using the four satellites) and to use least square method or
Kalman filter. Four nonlinear algebraic equations (pseudo-range measurements)
are usually used to determine satellite position with GPS receiver:
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
i
R ¼ ðx xi Þ2 þ ðy yi Þ2 þ ðz zi Þ2 þ cτ
(2)
i ¼ 1, 2, 3, 4
8
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
measurements (t > 4) with LSM or KF. Satellite velocity can be determined by the
differentiation of its position. Finally, GPS receiver can provide to AODC OBC
current satellite position and velocity in the reference (e.g., Earth-centered inertial
(ECI) [14] frame x, y, z, V x , V y , V z and synchronized (by GPS) onboard time ts ).
The TRIAD method [10] is applied when two different vectors are measured.
They usually can be any of the three pairs combined with the following three
vectors: Earth magnetic induction vector B (measured with three-axis MAG), Sun
vector S (measured with two-axis SS), and local vertical r (perpendicular to the
local infrared radiation temperature surface, measured with the HS). At least two
different not collinear vectors (their orientation) should be measured to deter-
mine satellite attitude that here is considered as satellite directional cosine matrix
(DCM) and related three Euler angles of the certain order of rotations (e.g., 3-2-1)
[9, 10] .
Let us assume that two different physical nature not collinear vectors U ¼ S and
V ¼ r are measured U m , V m by two vector measuring devises (SS and HS) installed
on the satellite board and both these vectors are referenced in the reference frame as
U r , V r . Let us choose U as the main vector and V as an auxiliary vector. Then an
orthogonal coordinate system (frame) with basis unit vectors, q, r and s can be
defined as follows [10, 20]:
!
! U
q¼ !
∣U ∣
! !
! UV (3)
r ¼ ! !
∣U V ∣
! ! !
s ¼qr
Cbr ¼ Cm Cr (5)
Three Euler angles of rotation, roll (ϕ), pitch (θ), and yaw (ψ), can be expressed
through the elements of the matrix Cbr . Certain trigonometric formulas depend on
the agreement about the order of the body rotations. For the order 3-2-1, the matrix
Cbr is as follows [9, 10]:
9
Satellite Systems - Design, Modeling, Simulation and Analysis
2 3
cθcψ cθsψ sθ
Cbr ¼ 4 sϕsθcψ cϕsψ sϕsθsψ þ cϕcψ sϕcθ 5 (6)
6 7
where c and s stand for cosine and sine angle. Then, formulas for Euler angles
can be derived from (6) as:
C231
ϕ ¼ tan ,
C33
θ ¼ sin 1 C13 , (7)
C12
ψ ¼ tan 1
C11
Figure 10.
S-vector sensor Bradford fine sun sensor, accuracy, 0.2 deg. (2σ).
Figure 11.
r-vector sensor HS CMOS/SRAM-modular infrared horizon sensor, accuracy, 0.4 deg. (2σ).
Figure 12.
B-vector sensor MAG TFM100-S, accuracy, 10mG (2σ).
If more than two vectors are measured and available for attitude determination,
then LSM-BATCH method [10] can be applied to use informational redundancy for
10
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
increasing the stochastic estimation accuracy. This method basically can be applied
for any set of VMD but is specifically convenient for the star tracker (ST), when
some number (n) of navigation stars are in the device field of view (FOV) and are
detected and tracked simultaneously, providing measured vectors Rm to these stars
that are referenced in the device space catalog Rr (Figures 10–13).
Figure 13.
Star direction R-vector measured sensor (optic and computer units). Advanced stellar compass, accuracy,
200 1600 ð2σ Þ.
where C is the DCM of the rotation from the reference frame to satellite body
frame and vectors Rr and Rm are written in the matrix form as matrix columns.
If the ST is in the tracking mode keeping in its FOV some n detected navigation
stars, then it can be assumed that C is a small-angle matrix that is independent of
the rotation order and can be expressed as follows:
2 3
1 αz αy
C ≈ 4 αz
6 1 αx 75 (9)
αy αx 1
δR ¼ δCRr (10)
2 3
0 αz αy
where δR ¼ Rm Rr , δC ≈ 4 αz
6 0 αx 75.
αy αx 0
Transforming in (10) matrix product and taking into account random measure-
ment errors, this equation can be represented in the following form:
δR ¼ Rr δC þ V (11)
Rrz Rry Vx
2 3 2 3 2 3
0 αx
2 3
δRx
where δR ¼ 4 δRy 5, Rr ¼ 4 Rrz 0 Rrx 7
5, δC ¼ α ¼ 4 αy 5, V ¼ 4 V y 5,
6 7 6 6 7 6 7
δRz Rry Rrx 0 αz Vz
11
Satellite Systems - Design, Modeling, Simulation and Analysis
zi ¼ hi x þ Vi
(12)
i ¼ 1, 2, ::n number of measured vectors
2 3
Rrz Rry Vx αx
2 3 2 3
0
4 Rrz
where, zi ¼ δRi, hi ¼ Rri ¼ 6 0 Rrx 5 , Vi ¼ 4 V y 5 x ¼ δC ¼ α ¼ 4 αy 5.
7 6 7 6 7
Rry Rrx 0 i Vz i αz
If n ¼ 1, only a single vector is measured, then dethi ¼ Rrx Rry Rrz þ
Rrx Rry Rrz 0, and thus, determining all three angles of satellite attitude is impossi-
ble. If it takes informational redundancy, then optimal estimate (providing mini-
mum of standard deviation of satellite attitude errors) can be found with the
following LSM formula [13]:
^ ¼ Kz
α (13)
1
where K ¼ HT R 1 H HT R 1 H, H ¼ ½h1 h2 h3 ::hn T .
Direct measurement
Satellite angular rate ω (vector absolute angular velocity) can be measured directly
by three-axis rate sensor (RS) that could have mechanical, optical, or microelectro-
mechanical systems (MEMS) design [21, 22] (Figure 14). Traditionally, independent
of the design type, these RS are usually named “gyros,” paying respect to their
historical appearance for aerospace vehicle control purposes as a mechanical
gyroscope (Figure 15).
Measured angular velocity vector ω can be used to determine satellite attitude by
the integration of matrix kinematic Poisson’s Equation [9, 10]:
^
C_ ¼ ωC, Cð0Þ ¼ C0 (14)
Figure 14.
Measurement of satellite angular velocity ω with three rate sensors RSx , RSy , RSz .
12
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
Figure 15.
ω-vector sensor RS, BEI QRS-11 single-axis body rate sensor, accuracy, 7deg=h ¼ 0:0019deg=s ð2σ Þ.
where C ¼ Cib C ¼
2 Cib is the DCM between
3 the inertial frame ECI and satellite
0 ωz ωy
^
body frame and ω ¼ 4 ωz
6 0 ωx 75 is skew symmetric matrix measured in
ωy ωx 0
the satellite frame components of vector of satellite absolute angular velocity. After
determination of the DCM, satellite attitude in three Euler angles can be derived with
Eq. (7) above. Unfortunately, gyro drift causes unlimited growing up errors in
integrated attitude that require periodic corrections from two VMD, measuring the
attitude directly (Figures 10–12).
Body rate estimator
Often, specifically for attitude stabilization (keeping or aka pointing) mode,
satellite angular rate is estimated by using the so-called body rate estimator and is not
measured directly by the RS. Indeed, using for attitude keeping mode small angles and
linear approximation, we can simplify satellite attitude dynamics model [9] to three
single-axis state equations and present it with the stochastic influences as follows:
8
<ω
> _ ¼w
α_ ¼ ω (15)
>
z¼αþv
:
where ω is the angular velocity, α is the satellite deviation angle from the desired
direction and z is the satellite deviation angle measurement with random Gaussian
white noise error v. Realistically, it is a wide spectrum correlated process that
components have spectral density ri ¼ 2σ 2vi T vi , i ¼ x, y, z (σ vi is standard deviation of
the random error vi , T vi is vi ðtÞ correlation time), wi is exciting angular acceleration
noise with spectral density qi ¼ 2σ 2wi T wi (σ wi is standard deviation of the random
angular acceleration wi ¼ MJi i , Mi is exciting external random torque, J i is satellite
moment of inertia, T wi is wi ðtÞ correlation time).
The linear KF can be applied to synthesize the estimator for the optimal
estimation of the vector angle α and the vector of angular velocity ω, using noisy
measurements z [9]:
(
^_ ¼ k12 ðz α
ω ^Þ
(16)
^_ ¼ ω
α ^ þ k22 ðz α^Þ
where α
^ and ω
^ are the optimal estimates of the angle α and the angular velocity
ω correspondingly.
13
Satellite Systems - Design, Modeling, Simulation and Analysis
q
where ξi ¼ rii is the ratio of spectral densities of satellite disturbing torque noise
to measured attitude error noise (assuming that both are white Gaussian noises).
This parameter can be considered as the filterability index. Eq. (16) can be
represented in the transfer function (Laplace operator) form as a second-order
differential equation unit:
8 s
> ω
^ i ¼ 2
zi
T i s2 þ 2di T i s þ 1
>
>
<
(18)
> 2dT i s þ 1
^¼ zi
>
:α
>
2 2
T i s þ 2di T i s þ 1
where s is the Laplace operator, T i is the time constant, and di is the specific
damping coefficient that is determined by the following formulas:
1 1
8
> T i ¼ pffiffiffiffiffiffiffi ¼ p ffiffiffiffi
k12i
> 4
>
>
< ξi
pffiffiffi (19)
k 2
>
22
>
: di ¼ pffiffiffiffiffiffiffi
i
¼ ¼ 0:707
>
>
2 k12i 2
As it can be seen from the consideration above, the use of directly measuring
devices (e.g., ST and RS) for attitude and body rate determination has a disadvan-
tage. The random noises are at the devices output, and they have to be filtered in the
closed control loop of satellite attitude control that puts some constraints to choose
the control law coefficients. However, using indirect body rate measurement, the
state estimator (filter) unavoidably introduces additional phase delay in the control
loop because of the consecutive inclusion of this filter in the control loop. To use the
RS (gyro) and the integrator for body rate and attitude determination autono-
mously for a long time is not possible because of the accumulated attitude errors
caused by the integration of the gyro drift. The following scheme (that is common
in Aviation) can be considered as free from the disadvantages above. Let us assume
that satellite attitude is determined in two ways: continuous integration of RS
angular velocity (IMU) and using VMD, for example, ST. Then this ST is used to
correct the attitude derived by the integration of RS output. The idea of MSU is
shown in Figure 16.
14
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
Figure 16.
Integration of multisensory sensor unit (MSU) single-axis channel.
< δω_ ¼ w1
8
>
IMU model equations : δα_ ¼ δω (20)
>
z ¼ δα þ Δα
:
where KF coefficients k12 and k22 are determined by (17), substituting there q1
and r1 instead of q and r.
3.2.2 Actuators
Satellite propulsion system [9, 10] is usually designed for satellite orbital and/or
angular control. In the first case, PS is commanded from the ground OC by TLM
commands in some cases when satellite orbit has to be changed (orbit correction,
deorbiting, collision avoidance), in the second controlled automatically from
onboard AODCS. It consists of such typical elements as orbital and attitude
thrusters (number and installation scheme depending on certain application), pro-
pulsion tank with associated pipes, valves, regulators, and electronics. General
principles of PS act independently of the type (ion thrusters (0.01–0.1 N), liquid
propellant and solid motor (100–10,000 N), cold gas (1-3 N)).
15
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 17.
Satellite control with PS thruster principles.
Figure 15 illustrates the satellite control with PS thruster principles. The princi-
ple of the formation of the propulsion jet force can be presented by the following
equation of variable mass body dynamics that from Russian sources, for example,
[23], is known as Prof. I. Meshchersky’s equation:
dV dm
M ¼F Vp (22)
dt dt
where M is the mass body, F is the external force, V P is the propellant exhaust
velocity and m is the propellant mass. The term V p dmdt is the propulsion force
(propulsion thrust) (Figure 17).
dm
T¼ Vp (23)
dt
16
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
Figure 18.
Prof. I. Meshchersky (1859–1935).
Figure 19.
K. Tsiolkovsky (1857–1935).
Figure 20.
Cold gas GN-2 thruster, nominal thrust 3.6 N (230 psi), specific impulse 57 s.
the coil is turned off because of the hysteresis in the magnetization curve of the
core. It is therefore necessary to demagnetize the core with a proper demagnetizing
procedure. Normally, the presence of the core (generally consisting of ferromag-
netic) increases the mass of the system. The control voltage is controlled by AODCS
control output (Figures 18–20). The magnetic dipole generated by the
magnetorquer is expressed by the formula:
M ¼ niS (25)
17
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 21.
60 liter propulsion gas tank.
where n is the number of turns of the wire, i is the current provided, and S is the
vector aria of the coil. The dipole interacts with the Earth magnetic field, generating
a torque whose expression is:
Tm ¼ M B (26)
where M is the magnetic dipole vector moment, B is the Earth magnetic field
induction vector and T m is the generated magnetic torque vector. This equation in
the scalar form is as follows:
T mx Bz By Mx
2 3 2 32 3
0
7 6 B Bx 7
4 T my 5 ¼ 4 0 54 My 5 (27)
6 6 7
z
T mz By Bx 0 Mz
Typically, three coils are used; the three-coil assembly usually takes the form of
three perpendicular coils, because this setup equalizes the rotational symmetry of the
fields which can be generated; no matter how the external field and the craft are placed
with respect to each other, approximately the same torque can always be generated
simply by using different amounts of current on the three different coils (Figure 22).
It can be seen from Eq. (26) that MTR cannot generate the magnetic torque in
the direction that is parallel to Earth magnetic field B (MkB) and it always is
perpendicular to the Earth magnetic field vector T m ⊥B. Unfortunately, from
Eq. (27), the required magnetic moments cannot be found, because it has zero
determinant, and we cannot invert it.
Bz By
0
Δ ¼ Bz 0 Bx ¼ Bx By Bz þ Bx By Bz ¼ 0 (28)
By Bx
0
However, the following approach can be used to find required vector M [9].
When we take cross-product of B with both sides of Eq. (26) and take into account
that it is useless to apply M in parallel direction to B, and hence we can require that
M⊥B and M B ¼ 0, then the following formula can be derived:
18
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
Figure 22.
3D orthogonal magnetic torque rods.
B Tm
M¼ (29)
B2
M¼ kB (30)
Figure 23.
Magnetic torque rod SSTL MTR-30, magnetic moment, M ¼ 30 Am2 .
19
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 24.
Three orthogonal reaction wheels (RW).
sensor (tachometer) measures the motor rotor angular speed and allows organizing
RW inertia torque control proportional to the applied control voltage. At least three
RW, as in Figure 24, are required to produce control torque vector in three-
dimensional spatial, having desired value and pointed in the desired direction.
Sometimes redundant unit of three redundant wheels or one with 4th redundant
skewed wheel is applied to meet reliability requirements.
RW can generate control inertia torques Ti only when they are accelerated or
decelerated. With this torque, they cannot compensate a permanent disturbance
torque Td ¼ const applied to satellite for a long enough time and come eventually to
some maximum/minimum available speed aka the saturation speed. At that state T i
becomes zero. That is why RW is usually applied with other types of actuators such as
MTRs or gas thrusters to de-saturate the RW to use them as the source of Ti again.
In general case, RW can be run around in some nominal angular speed Ω0 . In
this case they can be named as reaction-momentum wheel (RMW) and then just
RW. Mathematically RW/RMW dynamics can be presented as follows. Let us con-
sider satellite with RW unit angular momentum:
Hs ¼ H þ h (31)
0 0 Izz
Then differentiating (31) in rotating with angular velocity satellite axis ω and
using Euler’s rigid body dynamics formula [9], we can get the following equation:
h_
^ ^
_ þω
H H¼ ωh þ T (32)
2 3
0 ωz ωy
^
where ω ¼ 4 ωz
6 0 ωx 75 is the satellite angular velocity matrix and T is
ωy ωx 0
the vector of the external torque applied to the satellite. In the right side of Eq. (32),
20
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
we can see two terms that have meaning of torques applied from the RMW unit to
the satellite body: h_ ¼ Ti the inertial torque and ωh ¼ Tg the gyro torque. They
^
are RMW generated torques that can be used for the satellite attitude control.
Eq. (32) presents satellite attitude dynamics under the action of RMW torques. The
RW dynamics can be presented similar to Eq. (32):
h_ w þ ωhw ¼ Mw þ Mf
^
(33)
h_ ¼ IΩ
^
_ ¼ Iω
_ ωðh þ hc Þ þ Mw þ Mf (34)
In the operator Laplace s-form (transfer function), Eq. (36) can be rewritten as
follows:
21
Satellite Systems - Design, Modeling, Simulation and Analysis
ωz H y þ h y ¼ T x
(
(38)
ωx H y þ h y ¼ T z
Figure 25.
Reaction/momentum wheel HR-0610, torque, 75 10 3 Nm; momentum, ð4 12Þ Nms.
22
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
Figure 26.
Satellite OBCS, MAC-200 (C&DH unit with AODCS card) comprises of two OBC: Prime and redundant
(cold reserve).
(OBC (HW + SW), sensors, and actuators) functional test [2] that should be
performed in the Space Qualification Laboratory [7] during satellite Space Qualifi-
cation and Acceptance campaign can really minimize the risk of launching a not
ready satellite and prevent against AODCS refinishing in orbit during commission-
ing and operation.
VS can be separated in two parts, ODCS SW and ADCS SW. For both parts, I/O
interface with sensors and actuators is determined in special interface control
document(s) (ICD), describing type, certain connectors, and electrical parameters
of the exchanging data. These data before using them for functional tasks are
pre-processed in OBC with special algorithms.
• Convert data into required physical parameters and units, taking into account
certain sensor input–output scale function.
• Transfer to C&DH TLM data about sensor/actuator state and their data.
3.2.3.2 ODCS SW
23
Satellite Systems - Design, Modeling, Simulation and Analysis
x xz2 x
x
€¼ μ þ AJ2 15 3 ,
r3 r7 r5
y yz2 y
€y ¼ μ þ A J2 15 3 ,
r3 r7 r5
(39)
3
z z z
€z ¼ μ þ AJ2 15 7 9 ,
r 3 r r5
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
r ¼ x 2 þ y2 þ z 2 ,
24
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
is the satellite orbit right ascension of ascending node angle (RAAN), uis the
argument of latitude angle, ν is the satellite true anomaly angle, and ω is the satellite
orbit argument of perigee angle (Figure 27).
Figure 27.
Satellite orbit in the inertial ECI (XYZ) coordinate system.
3.2.3.3 ADCS SW
This group of algorithms was presented above in 3.2.1.2 and can be used here.
For example, let us consider single-axis stabilized satellite that should keep one
axis (e.g., Z) permanently pointed to the Sun as in Figure 28. Only two angles of the
satellite deviation from this direction and their angular velocities are required to
know to point and keep it in this direction. The satellite has two-axis Sun sensor that
Figure 28.
Satellite pointed by the Z-axis to the sun.
25
Satellite Systems - Design, Modeling, Simulation and Analysis
can measure two angles αx , αy of satellite deviation from the sun direction. Its axes
coincide with the satellite axes xyz. The axis z is the sensitivity axis that nominally
should point into the Sun’s direction (center of brightness), and xy is the focal plane.
The Sun vector is referenced in the Sun frame as Sr ¼ ½0 0 ST and is measured in
T
the Sun sensor measured frame as Sm ¼ Sx Sy Sz .
In Figure 28, xyz is the satellite body frame, xs ys zs is the Sun reference frame, SS
is the two-axis Sun sensor and αx and αy are the turn angles of satellite x and y axis
accordingly.
The following formula represents the mathematical transformation of the Sun
vector from the reference into the body frame:
Sm ¼ Cbs Sr (41)
T
where Sr ¼ ½0 0 ST , Sm ¼ Sx Sy Sz , and Cbs is the DCM between the refer-
ence (Sun) and measured (satellite) frames. Let us consider that the order of
rotation from the Sun to the satellite frame is 3-2-1 (αz , αy , αx); the DCM matrix Cbs is
as follows [9]:
2 3
cos αy cos αz cos αy sin αz sin αy
Cbs ¼ 4
6 cos αx sin αz þ sin αx sin αy cos αz cos αx cos αz þ sin αx sin αy sin αz sin αx cos αy 7
5
sin αx sin αz þ cos αx sin αy cos αz cos αx cos αz þ cos αx sin αy sin αz cos αx cos αy
(42)
From (43), desired angles and can be derived that can be used for satellite
attitude control.
Sym
8
αx ¼ tan 1
>
>
Szm
>
>
<
1 Sxm (44)
>
> α y ¼ sin q ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
S2 þ S2
>
>
ym zm
:
Let us also assume that the satellite does not have angular velocity sensors RS and
its angular velocities should be derived from the measured angles αx and αy . Simple
low-frequency first-order differentiating fitters can be applied for this purpose.
Laplace operator s-form (transfer functions) of these filters are presented below:’
8
^ s
< α_ x ¼ T s þ 1 αx
>
>
fx
s (45)
: ^_y ¼
>α αy
T fy s þ 1
>
26
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
s
T cx ¼ kpx αxm kdx αxm (46)
T fx þ 1
where kpx is the position control coefficient, kdx is the damping control coeffi-
cient, s ¼ dtd is the Laplace operator, T fx is the differentiating filter time constant,
αxm ¼ αx þ Δαx is the measured angle αx ,αx is the true value and Δαx is the
measured error. Let us assume that this torque is generated by only one MTR Y.
Eq. (27) is determined as:
T cx ¼ Bz My ¼ Bz kTRy uy (47)
s
uy ¼ K px αxm K dx αxm (48)
T fx þ 1
k
where K px ¼ Bz kpxTRy and K dx ¼ BzkkdxTRy are position and damping magnetic control
coefficients.
Let
2 us take a 3ball-shaped satellite with the inertia matrix as follows:
Jx 0 0
J ¼ 4 0 J y 0 5 where J x ¼ J y ¼ J z J x ¼ J y ¼ J z . Then in inertial spatial, its linear
6 7
0 0 Jz
angular dynamical equations for the axis X can be approximately written as follows:
J x s2 αx ¼ T cx þ T dx (49)
27
Satellite Systems - Design, Modeling, Simulation and Analysis
where s ¼ dtd is the Laplace operator, T cx is the control torque, and T dx is the
disturbing external torque (satellite residual and induction magnetism torque,
atmosphere drug torque, solar pressure torque). Then substituting in Eq. (49)
Eqs. (47) and (48), we can rewrite it as follows:
s
J x s2 αx ¼ kpx αxm kdx αxm þ T dx (50)
T fx s þ 1
Let us divide all terms in Eq. (50) by the coefficient kp and substitute αxm value,
then it can be represented in the following form:
!
s 1
T 2x s2 þ 2dx T x þ 1 αx ¼ T dx Δαx (51)
T fx s þ 1 kpx
qffiffiffiffiffi
Jx
where T x ¼ kpx is the X control channel time constant and dx ¼ pkffiffiffiffiffiffiffi
dx ffi
is the
2 kpx J x
Xcontrol channel-specific damping coefficient. Eq. (51) is a third-order linear differ-
ential equation and could be analytically analyzed. In particular its stability can be
analyzed with algebraic Hurwitz criterion [13]. However, more simple and general
results can be obtained with the following approximate consideration. If T fx < < T x ,
then by the filter time constant, T fx can be neglected, and (51) can approximately be
considered as a standard second-order control unit, presented by the second-order
linear time invariant (LTI) differential equation and rewritten as follows:
1
T 2x s2 þ 2dx T x s þ 1 αx ≃ T dx
Δαx (52)
kpx
1
T 2x s2 þ 2dx T x s þ 1 αx ≃ T dx
Δαx (53)
kpx
where T dx0 ¼ const and Δαx0 ¼ constΔαx0 ¼ const (ATT sensor bias).
pffiffi
For Eq. (52), the optimal damping coefficient is dx ¼ 22 ¼ 0:707 [13].
Numerical example
Let us evaluate satellite time constant T x. Let us assume that for the LEO satellite
magnetic field induction vector B has the following value of the projection on the
Sun direction.
and MTR has the following parameters: maximal magnetic moment My max ¼
35Am2 , maximal control current Iy max ¼ 100 mA, winding resistance R ¼ 280 Ohm;
My max
and maximal control voltage uymax ¼ R Iy max ¼ 28 V. MTR gate is kTRy ¼ uy max ¼
2
35Am
28V ¼ 1:25Am2 =V. Then maximal available magnetic torque is T cx max ¼ 10 Nm. If 3
maximal linear zone for this control channel is αy max ¼ π2 rad, then the position
control coefficient kpx can be calculated with the following formula:
T cx max
kpx ¼ (54)
αmax
0:001
For the data above, it has the value of kpx ¼ 3:140:5 ¼ 6:366 10 4 Nm=rad.
Let us consider example of the first Soviet satellite “Sputnik” (SS-1) that had the
28
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
mass m ¼ 83:6 kg and the radiusR ¼ 0:29 m and take the assumption that its
mass was homogeneously distributed within its spherical volume of V ¼ 43 πR3 ¼
0:1 m3 , and then its inertia ( J ¼ J x ¼ J y ¼ J z Þ can be calculated with the following
formula:
2
J ¼ mR2 (55)
5
as for a homogeneous sphere. Substituting into Eq. (55) the data above, we can
calculate that for SS-1 J ¼ 2:82 kgm2 , then its time constant with the considered
MTR control might be:
sffiffiffiffiffiffiffi sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
Jx 2:82kgm2
Tx ¼ ¼ ¼ 66:56 s:
kpx 6:366 10 4 kgm2 s 2
4
It has the following value: kdx ¼ 2 0:707 66:56 6:366 10 s Nm=rad ¼
0:06 Nms=rd.
k
Finally, K px and K dx can be calculated. They are as follows: K px ¼ Bz kpxTRy ¼
6:36610 4 Nm=rad
2:8610 5 1:25TAm2 =V
¼ 17:8 V=rad and K dx ¼ BzkkdxTRy ¼ 2:8610
0:06Nms=rad
5
1:25TAm2 =V
¼ 1:678
3
10 sV=rad.
When T x is determined, then the time constant for the differencing filter T fx can
be chosen from the condition that T fx < < T x . In our example T x ¼ 66:56 s, let us
take that T fx ¼ 5 s.
Simulation
Eqs. (51) and (52) were simulated using Simulink (see Figure 29).
Figure 29.
Satellite single-axis attitude control Simulink block scheme.
29
Satellite Systems - Design, Modeling, Simulation and Analysis
Blocks in the pink color present the satellite model, the dark green color is for
control law blocks, the cyan blocks are registration oscilloscopes, and the display and
the orange color are the disturbances. The red manual switch allows to implement the
differentiating filter, transforming the scheme from the approximation (52) to the
accurate presentation (51). Disturbing external torque Md is constant; attitude sensor
error is represented by the constant value ALP0 and limited range white noise V that
has spectral density SV ¼ 2σ 2V T v (σ v is the standard deviation (SD), T v is the
correlation time). The results of the simulation (Figure 29) with and without the
differentiating filter (with the assumption that α_ x α_ x is directly measured without any
errors) are presented below in Figures 30–34 (left A, (52), without the filter; right B,
(51), with the filter). The numerical data for the simulation are as follows:
sffiffiffiffiffiffiffi
Jx kdx
J ¼ 2:82 kgm2 , T x ¼ ¼ 66:56 s, dx ¼ ¼ 0:707,
kpx 2T x kpx
4
kpx ¼ 6:366 10 Nm=rad, kdx ¼ 0:06 Nm s=rad, T fx ¼ 5 s,
Md ¼ 10 5 Nm, Δαx0 ¼ dltALP0 ¼ 0:1∘ , Δαx0 ¼ dltALP0 ¼ 0:1∘
Figure 30.
Response to initial deviation angle αx0 ¼ 10 . (a) without dif. filter and (b) with dif. filter.
Figure 31.
Response to initial angular velocity α_ x0 ¼ 0:01deg=s.
Figure 32.
Response to external disturbance torque T d ¼ 10 5 Nm. Static error αx∗ ¼ 0:898∘ .
30
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
Figure 33.
Response to attitude sensor bias Δx0 ¼ 1∘ plus white noise σ V ¼ 0:1∘ , T V ¼ 1 s. Satellite attitude stabilization
errors, ALP.
Figure 34.
Satellite attitude measured errors ALPm.
4. Conclusion
Acknowledgements
The author wishes to express his sincere gratitude to the Canadian Space
Agency, where he had the opportunity to learn and possess the knowledge and
experience related to the writing of this chapter. As well, he is very thankful to
many of his colleagues from Canadian Magellan Aerospace Company (Bristol Aero-
space Division) with whom he discussed and analyzed satellite AODCS design pro-
jects and issues that helped him to work out the system analysis and its principal
concepts presented in this chapter. Additionally, he cannot forget that his
31
Satellite Systems - Design, Modeling, Simulation and Analysis
Note
Author details
Yuri V. Kim
David Florida Laboratory, Canadian Space Agency (CSA), St. Hubert, Ottawa,
Canada
© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License (http://creativecommons.org/licenses/
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
32
Satellite Control System: Part I - Architecture and Main Components
DOI: http://dx.doi.org/10.5772/intechopen.92575
References
33
Satellite Systems - Design, Modeling, Simulation and Analysis
34
CORRIGENDUM Chapter
Yuri V. Kim
David Florida Laboratory, Canadian Space Agency (CSA), St. Hubert, Ottawa, Canada
Corrigendum to: Kim YV. Satellite Control System: Part 1-Architecture and Main
Components. In: Nguyen TM, editor. Satellite Systems – Design, Modeling,
Simulation and Analysis. London: IntechOpen; 2020. DOI: 10.5772/intechopen.92575.
Changes have been made in online and print versions of the chapter.
This corrigendum is published in agreement with the authors.
The publisher and the author regret any inconvenience this might have caused to the readership.
References
[1] Kim YV. Satellite Control System-Architecture and Main components. In: Nguyen TM, editor. Satellite Systems
– Design, Modeling, Simulation and Analysis. London: IntechOpen; 2020. DOI: 10.5772/intechopen.92575.
Chapter
Abstract
The chapter presents a game theoretic training model enabling a deep learning
solution for rapid discovery of satellite behaviors from collected sensor data. The
solution has two parts, namely, Part 1 and Part 2. Part 1 is a PE game model that
enables data augmentation method, and Part 2 uses convolutional neural networks
(CNNs) for satellite behavior classification. The sensor data are propagated with
the various maneuver strategies from the proposed space game models. Under the
PE game theoretic framework, various satellite behaviors are simulated to generate
synthetic datasets with labels for the training to detect space object behaviors. To
evaluate the performance of the proposed PE model, a CNN model is designed and
implemented for satellite behavior classification. Python 3 and TensorFlow are used
in this implementation. The simulation results show that the trained machine learn-
ing model can efficiently and correctly classify the satellite behaviors up to 99.8%.
1. Introduction
Since space has already been fully utilized, society has become increasingly
dependent on space advantages in various industrial, civil, and commercial appli-
cations. This dependence brings an essential vulnerability, especially shortage of
continuous situational awareness of the space environment to ensure freedom of
movement. Due to the fact that information from space is crucial for key decision-
making, such as urban, agricultural, and responsive planning, space is regarded as a
significant frontier. In addition to real-time and hidden information constrains, the
existence of space object density significantly produces the complexity of the space
situational awareness (SSA). Understanding the position of space objects from low-
level information fusion can support high-level information fusion SSA missions
of sensor, user, and task refinement [1]. In order to implement SSA accurately, it
is possible to coordinate the evaluation of residential space objects (RSO) through
user-defined operation pictures (UDOP) [2].
Space control and SSA are required for space prevalence, which depend on fast
and precise space object behavioral discovery. Developing a theoretical approach
for fast discovery of variation of satellite behaviors is the main task of this book
1
Satellite Systems - Design, Modeling, Simulation and Analysis
chapter. The machine learning methods with novel neural networks are proposed
in this case. However, there are numerous challenges for constructing the tools
because of the following reasons: (i) partially observable movements, (ii) resident
space objects (RSOs), (iii) uncertainties modeling and propagation, (iv) real-time
response, and (v) computationally intractable algorithms.
Space access investigation and mission trade-off considers are imperative for
the victory of space-borne operations. The tracking algorithms of space object can
be measured depending on collecting data to track satellites, debris, and natural
phenomena (such as comets, asteroids, and solar flares). Tacking is related with
sensor administration, which can point sensors to observation points to decide the
circumstance and to aware threatens. SSA improvements consist of models (such as
orbital mechanics), measurements, computational software (such as tracking), and
application-based system coordination (such as situations). For instance, game-
theory methods for SSA can be utilized for pursuit-evasion analysis [3].
This book chapter creates and establishes game theoretic training enabled deep
learning (GTEL) methodologies for fast discovery of the behaviors from the satel-
lites. GTEL is an adaptive feedback adversarial theoretic method, which acquires
data from sensors related to the relationship between resident space objects (RSOs)
of interest and sensing assets from ground and space (GSAs). Thus, a game theory
is modeled instead of a control problem for this circumstance. Game reasoning
uses data-level fusion, random modeling/propagation, on the other hand, RSO
detection/tracking predicting the future RSOs-GSAs relationships. The adversarial
engine also supports optional space pattern dictionary/semantic rules for adaptive
transition in the Markov game. In the event that no existing pattern dictionary is
accessible, GTEL will construct an initial pattern and modify it during the game
inference. The output of GTEL inference consists of two parts of control methods:
(i) measurements processing and (ii) RSOs localization. The two parts establish a
game-equilibrium, one of which is sensing asset management, the other of which is
the estimation of RSO behaviors.
The chapter is organized as the following. Section 2 introduces the comprehen-
sive system design for our methods. Section 3 presents the Markov game theory
with satellite maneuvering. Section 4 proposes the details of our machine learning
model for space behavior detection. And the numerical results and analysis are
displayed in Section 5. At the end, Section 6 draws the conclusion of this chapter.
The proposed methodology of GTEL is shown as Figure 1. The core piece of the
method for detecting unknown patterns of space objects is Markov Game Engine.
Due to the patterns are obscure, there is no preparing training data accessible.
Therefore, the Markov Game Engine makes use of zero-shot learning [4] and trans-
fer learning to unsupervised classify unavailable target domain data by training
available data from source domain (i.e. simulate data). The knowledge adaptation
or domain transfer is performed through manifold learning to share intermediate
semantic embeddings (such as attributes) between labeled and unlabeled data. On
one hand, the manifold learning can reduce the dimension of sensing data (such
as azimuth angle, elevation angle, range, and range rate). On the other hand, the
manifold learning can also mitigate the difficulties of object tracking and detection
with fast space object behaviors’ detection. Moreover, in order to solve the uncer-
tainties in this task, the ℝ 5 × 𝕊 coordinate system is used with filtering technology
on optimal transport (OT). The essential components of our methodology are shown
as the following:
2
Game Theoretic Training Enabled Deep Learning Solutions for Rapid Discovery of Satellite…
DOI: http://dx.doi.org/10.5772/intechopen.92636
1. Markov game – The space conflicting situation is inferred by the Markov game
structure [5, 6], in which system state is represented by distributions rather
than deterministic values.
5. Manifold learning for data level sensor fusion – The raw sensor data is generally
high-dimensional. Considering that the measurement data is reflections of
witnessed satellites (only several parameters will determine their states), so it
can be reasonably assumed that the inherent dimension of the measurement
domain is low. The communication bandwidth can be saved with decreasing
the number of dimensions by using Manifold learning algorithms [10].
Figure 1.
GTEL system architecture.
3
Satellite Systems - Design, Modeling, Simulation and Analysis
estimation of the position. Afterwards, the estimations of the orbits are used for
collision alert and maneuver detection of space objects. The satellite maneuvers are
going to be interpreted as platform commands to perform course of actions and
space object movements.
Figure 1 shows an adaptive feedback approach with the game theory enabled. It
utilizes sensors to obtain information about the relationships between the ground/
space sensing assets and the RSOs of interest [11]. RSOs and GSAs determine the
relations. Thus, rather than saying it is a control problem, it is a dynamic game.
Data-level fusion, game reasoning, RSO detection/tracking, and uncertainty mod-
eling/propagation are combined to predict the RSO-SA relationships in the future.
Optional space behavior dictionary/semantic rule for adaptive transition matrices in
our Markov game is also supported by our game engine. In addition, the game solu-
tion is an equilibrium, which is controlled by both space sensing asset management
and the (estimated) RSO behaviors.
Lloyd Shapley has invented a concept of stochastic game [6], which is a dynamic
game played by one or several players focusing on probabilistic transitions. There
are several stages for this game. At first, the game is set in one state. Then, the par-
ticipated players should select an action individually. Based on the current state and
the actions players chosen, each player will receive a reward. Therefore, after the
chosen from each player, the game comes to a new random state, where the previous
state and previous actions chosen by the players determine a distribution of the new
random state. Afterwards, the above action will be repeated again for the new state.
After finite or infinite number of stages of playing, the total reward for each player
is obtained using the discounted sums of each stage reward or the averages of every
stage rewards. In this way, each player gets a reward and the reward is compared
with each other. The aforementioned sequences are the procedures for the sto-
chastic game, which can be generalized by Markov decision processes (MDP) with
repeated games. Our space situational awareness (SSA) would utilize this game tool
for intent prediction [12].
The Markov game engine extracts specific information from each event as the
following: (i) a finite set of players N, (ii) a finite or infinite set of states, S, (iii) a
finite set of accessible actions for each player in N set, Di (the overall action space
is D = ×i∈NDi), (iv) a transition rule q: S × D→∆(S), (where ∆(S) is the space of all
probability distributions over S), and (v) a reward function r: S × D→RN.
Figure 2 shows a visual description of the simple game states with only two play-
ers, who have only two options of actions for each player. The arrows in the graph
indicates the probable transitions from one state to the other state. The states with
red color means that only player 1 changes the approach. On the contrary, the states
with blue color indicates only player 2 changes the approach. The state with green
color indicates both players change the approach.
The GTEL solution uses a two-player Markov game to investigate the sensor
management for tracking space objects. Whether deliberate or unintentional,
some of space objects may cause confusion to observers (sensors) when the orbital
maneuvers are performed. In general, spatial object tracking can be assumed as an
optimal control problem (one side optimization) or a game problem (two side opti-
mization). For the settings of optimal control, the position and velocity of the space
objects will be calculated (filtered) as the states dependent on the measurement
from the sensor. However, this method ignores the possibility that the space objects
may alter their orbits purposely with intelligence. It may cause difficulties for the
4
Game Theoretic Training Enabled Deep Learning Solutions for Rapid Discovery of Satellite…
DOI: http://dx.doi.org/10.5772/intechopen.92636
satellite to track the space objects. Therefore, the Markov game method provides
a solution for these difficulties. In this approach, on the one hand, the observed
satellite will utilize the tracking and sensing model to destroy the tracking estima-
tions by confusing the observer. On the other hand, the observer figures out ways to
minimize the uncertainties of tracking, where the uncertainties are dependent on
the entropy of tracking.
In this chapter, the information uncertainty of the GTEL pursuit-evasion (PE)
game method [13] was exercised with a circumstance of two satellites, one of which
is Geostationary Earth Orbit – GEO (observed satellite), the other of which is space
based Low Earth Orbit-LEO satellite (observer satellite). Figure 3 provides an illus-
tration of a space based optical (SBO) sensor measurement model. The angle from
the line from the object to the SBO and the line from the sun with object is defined
as Bistatic Solar Angle, represented by θ. The smaller the angle is, the stronger
the lighting conditions. Therefore, it causes difficulties for observations when the
angle is large because of saturation of lighting. As shown in Figure 4, the scenarios
of light have shown. In the graph, the blue line is for the LEO orbit, the green line
Figure 2.
A diagram of states in a Markov game.
Figure 3.
SBO with a Bistatic solar angle.
5
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 4.
LEO and GEO based on sensor management and maneuver strategies with game theory.
Figure 5.
The performance of tracking using theoretical Markov game strategies.
indicates the GEO orbit with maneuvers, the pink lines are the SBO sensor per-
formed to track the GEO in order to lower the uncertainty, and the red line displays
the sunlight from the sun to the earth.
The research scenario is shown in Figure 4, where the red line indicates the
direction of sun light, green color is for the GEO orbit with maneuvers, blue color
6
Game Theoretic Training Enabled Deep Learning Solutions for Rapid Discovery of Satellite…
DOI: http://dx.doi.org/10.5772/intechopen.92636
for LEO orbit, and pink lines indicate when the SBO sensor resource is used to track
the GEO (use the sensor data to lower the uncertainty).
Figure 5 shows the results of tracking based on intermitted measurements.
Both cubature KF (CKF) and extended Kalman Filter (EKF) trackers are shown.
With the increase of the period without measurement due to the Earth blockage,
the tracking errors increased as well. In addition, the tracking errors increased with
the maneuver actions from the satellite. On the contrary, informational entropy
decreased by sensor measures in this process.
Figure 6 top graph displayed the PE game control results, with α and β as angles
of the maneuver thrust. The zoom-in view of the game optimal controls is shown
Figure 6.
The maneuver controls dependent on the PE game solution for the satellite direction.
Figure 7.
Sensor’s game theoretic on–off controls and associated information gains.
7
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 6 below, which exhibited that the observer satellite can keep the tracking
uncertainty within a desired level while saving the resources of sensor. Figure 7
shows the on–off sensor controls and the associated information gains. The results
show that with the larger potential information gain, the sensor would use the
resources to take measures (for the observer’s on–off control, 0 means turning off
and 1 indicates turning on).
The machine learning (ML) details will be described in this section. The main
purpose of construction this machine learning scheme is to detect the behaviors
of the resident space object (RSO) by fusing sensors data from multiple sources,
including the velocity, orbital energy, angular momentum, and the position of RSO
compared to the station. Similar with other machine learning model, this model is
generally trained off-line by using generated data. Then the generalized weights will
be deployed in the application using TensorFlow deployment. In addition, in order
to improve the robustness of the trained system, we proposed a neural network
scheme with the ability to train the newly-added unknown pattern online with only
tiny modifications of the weights.
As shown in Figure 8, the RSO pattern classification architecture is displayed. It
consists of two separate parts, one of which is offline part as Modeling, the other of
which is online part as Monitoring, to detect the RSO behavior pattern.
As the matter of the offline part, named as Modeling RSO Behavior Pattern, our
neural networks will be trained as a classifier by using the collected simulated data.
The data specifically indicates the different behaviors of the RSO. As shown in the
graph, some useful features are obtained from different sensors in the sessions of
feature extraction. Subsequently, the extracted feature will be sent into the training
model for neural network tuning and training to generate a classifier, which is used
to identify the RSO behavior with fine grained size. On the contrary, the online
part, Monitoring RSO Behavior Pattern, acquires RSO patterns in real-time to distin-
guish the abnormal behaviors. In this way, a warning message would be prompted if
there are any abnormal behaviors detected.
Figure 8.
The scheme of RSO behavior classification.
8
Game Theoretic Training Enabled Deep Learning Solutions for Rapid Discovery of Satellite…
DOI: http://dx.doi.org/10.5772/intechopen.92636
Additionally, the classifier generated by the machine learning methods has sev-
eral properties as the following: (a) ability to fuse heterogeneous and complex input
data, (b) scalability, (c) robustness with perturbations relative to the data, (d) high
accuracy, and (e) explicitly. However, it is hard to fulfill all the requirements as
shown above, therefore, some trade off would be considered for our model. Among
different structures of neural networks, convolutional neural networks fit our case
appropriately.
Figure 9.
The framework of the convolutional neural network.
9
Satellite Systems - Design, Modeling, Simulation and Analysis
the filter can extract the features reasonably with few weight parameters. Thus,
using the CNN filters can increase the accuracy of our neural network classifier.
With the development of CNN, it successfully proceeds to be the best model
structure for computer vision and image processing tasks. Moreover, CNN has also
been used by Natural Language Processing (NLP) area. With small filters moving
across the input dataset, the filters, with only small numbers of parameters, are
re-utilized to recognize patterns for the large image. Therefore, with the similar
classification abilities, CNN network is faster to train and predict compared with
DNN. After several filters, the output part is flattened to bypass several dense layers
to produce the softmax activation at the last layer with the sparse cross-entropy loss
function to backpropagation for the behavior of RSO pattern prediction.
Our GTEL PE method utilizes the CNN architecture to classify the RSO pattern
with observed data. Compared with the other conventional methods, the convo-
lutional neural networks can process the RSO observation much faster. Therefore,
we utilize the Python and TensorFlow with Keras [15] as the fundamentals for
code write up. Although, each filter is computationally expensive to be trained,
the overall CNN architecture is faster to be trained to provide similar accuracy
classification.
A typical CNN structure is shown in Figure 10. After training our CNN-DNN
with our training data, the test data (with 10–20%) will be employed to evaluate
the generalization of our CNN-DNN model. Additionally, in order to solve the
Figure 10.
CNN architecture for RSO behavior classification.
10
Game Theoretic Training Enabled Deep Learning Solutions for Rapid Discovery of Satellite…
DOI: http://dx.doi.org/10.5772/intechopen.92636
overfitting problem of our model, several dropout layers were added after each
layer. This model outperforms other traditional methods with better accuracy and
higher computationally efficiency.
The revised deep learning neural network structure for 143-sapce-behavior (we
partition the pointing angles of α and β into cells with 15 degrees each and the total
cell number is 143) is shown in Figure 11. At first, there are 3 different dimensions
in our raw data, which are 3 parameters dimension, 15 times (each track consists
of 15 observation measurements) interval dimension, and 1 number of channels
dimension. Therefore, the convolutional neural network (CNN) comprises an
initial 3*15*1 input dimension with 72,000 samples as the raw dataset. In order
to distinguish the 143 different labels of satellite behaviors (15-degree difference
for between each behavior), 128 “2 × 2” filters are utilized with the same padding
for the first convolutional layer. Thus, there are 3*15*128 dimensions of output
Figure 11.
Revised CNN structure for 143-label data.
Table 1.
The parameters of our 143 classification CNN-DNN model.
11
Satellite Systems - Design, Modeling, Simulation and Analysis
after the first layer. Then the previous output bypasses another two convolutional
layers with 128 and 256 “1 × 5” filters together with the same paddings. After the
first three convolutional layers, the dimension of the data exploded to 3*15*256.
Figure 12.
Model structure for the CNN_DNN network.
12
Game Theoretic Training Enabled Deep Learning Solutions for Rapid Discovery of Satellite…
DOI: http://dx.doi.org/10.5772/intechopen.92636
Afterwards, another two convolutional layers are added with the same padding
methods to downsample the dimensions to 3*3*256 for future operation. The two
layers both has 1 × 7 filters with 512, and 256 third dimension respectively. Finally,
the 3*3*256 data is flattened to one dimension to bypass three dense layers with
1024, 512, and 512 neurons, respectively. All these layers utilized 30% dropout
parameters after dense layers. The employment of the dropout layer is used for
solving the overfitting problem during the training process. In order to predict
the 143 labels classifications, there is a 143-degree softmax attached at the last
layer of the CNN-DNN model. Additionally, the cross-entropy loss function with
the Adam Optimizer minimizes the result with the gradient of the next calcula-
tion points.
The details of the parameters of the CNN are shown in Table 1 and Figure 12.
There are 5,303,695 parameters that need to be trained in this large NN for detect-
ing the different behaviors. In order to train the GTEL CNN much faster and more
accurately, a learning rate decay is used during the training process. For the first 75
epochs, the learning rate is 1e−4, which can learn faster after the weight initializing.
After 75 epochs, the exponentially decay of learning rate finds the optimized path
for decreasing the loss for the CNN model.
In order to capture the training data with the two-line elements (TLEs) from the
space-track.org, we modified the catalog tracking by the addition of the maneu-
vers. The different maneuvers indicate different labels in our training data. As an
example, the maneuver to increase the energy of the orbital is labeled as 1. On the
contrary, the maneuver to decrease the energy of orbital is labeled as 2. In this way,
adding no maneuver is labeled as 3. The details related to adding maneuvers are
shown as the following:
Therefore, almost 57,332 tracks were generated for training dataset with the
other 6371 tracks as the testing dataset. The original first 10,000 tracks are shown in
Figure 13 for 143-space-behavior maneuvers.
The data format is list as:
• Column 1: track id
13
Satellite Systems - Design, Modeling, Simulation and Analysis
After 1000 epoch optimization, the GTEL systems achieved almost 98% training
accuracy with 96% test accuracy as shown in Figure 14. There is a little overfitting
after 400 epochs of training, which can be solved by adding more dropout layers
and other methods. Figure 15 displays the cross-entropy loss during the training
process. The training set loss is always decreasing; however, the loss of validation
set is flat after 200 epochs of training, indicating some overfitting issue after 200
epochs of training of the data.
Figure 16 shows the confusion matrix to evaluate our model’s performance. As
shown, for both training and test datasets, the true labels are almost the same with
the predicted label by the GTEL CNN model. The performance evolution pattern is
shown in Figure 17.
Figure 13.
The first 10,000 training tracks with various space behaviors.
Figure 14.
Training and validation accuracy during the training process.
14
Game Theoretic Training Enabled Deep Learning Solutions for Rapid Discovery of Satellite…
DOI: http://dx.doi.org/10.5772/intechopen.92636
Figure 15.
Training and validation loss during the training process.
Figure 16.
Confusion matrix of training data and testing data.
15
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 17.
Overview of validation classes prediction accuracy during training epochs.
for all training process, indicating the great performance and reasonable training
configuration during training of the deep learning neural networks. Understanding
the meaning for these weights during training process with further checking is very
important in the CNN model design, which leads to more reasonable explanation in
the deployment of the model.
In addition, we investigated the filter image correlations in a broad overview.
As shown in Figure 19, rows and columns represent layers and image classes,
respectively. A sequential color scheme is used to encode the number of anomaly
16
Game Theoretic Training Enabled Deep Learning Solutions for Rapid Discovery of Satellite…
DOI: http://dx.doi.org/10.5772/intechopen.92636
17
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 18.
The weight changes in filters during iterations. More blue color indicates stronger variation for the filters during
iterations.
Figure 19.
The abstract version of correlation view.
18
Game Theoretic Training Enabled Deep Learning Solutions for Rapid Discovery of Satellite…
DOI: http://dx.doi.org/10.5772/intechopen.92636
6. Conclusions
In this chapter, a machine learning design has been presented and implemented
to discover space object behaviors. Our GTEL methods using CNN models the cir-
cumstance with PE game rather than a control problem. The stochastic modeling/
propagation, RSO tracking, and data level fusion are utilized to predict the relations
for future space by game reasoning. In order to generate the data for training, the
Marko game approach is used with maneuvering strategies. The method provides
a way to solve the SSA using unknown behaviors. Additionally, the unknown
behaviors exist where a satellite employs the tracking and sensing way to corrupt
the tracking estimates to perturb the sensors. On the other hand, the space sensors
decrease the uncertainties during tracking process. Finally, the CNN-DNN is used
to train the numerical results, where the accuracy is 98% for our classification RSO
model with 143 labels.
In the future, a multi-player game theory with adversarial network with SSA
will be employed to enhance the deep learning for sensor management, combined
tracking, as well as secure communications. Methods for diffusion-based coopera-
tive space object tracking [16] and block chain [17] are emerging as methods for
game-theoretical methods.
Acknowledgements
Author details
Dan Shen1, Carolyn Sheaff 2, Genshe Chen1*, Jingyang Lu1, Mengqing Guo1,
Erik Blasch3 and Khanh Pham4
2 Information Directorate, Air Force Research Laboratory (AFRL), Rome, NY, USA
© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License (http://creativecommons.org/licenses/
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
19
Satellite Systems - Design, Modeling, Simulation and Analysis
References
[8] Das N, Ghosh RP, Guha N, [16] Jia B, Pham K, Blasch E, Shen D,
Bhattacharya R, Mallick B. Optimal Chen G. Diffusion-based cooperative
transport based tracking of space space object tracking. Optical
objects in cylindrical manifolds. The Engineering. 2019;58(4):041607
Journal of the Astronautical Sciences.
2019;66:582-606 [17] Xu R, Chen Y, Blasch E, Chen G. An
exploration of Blockchain-enabled
[9] Das N, Deshpande V, Bhattacharya R. decentralized capability-based access
Optimal-transport-based tracking of control strategy for space situation
space objects using range data from awareness. Optical Engineering.
a single ranging station. Journal of 2019;58(4):041609
Guidance, Control, and Dynamics.
2019;42(6):1237-1249
20
Chapter
Abstract
This chapter discusses existing and future trends on the design and build of
“Modular” and “Open” satellite Bus and mission payload along with practical design
issues associated with the use of Modular Open System Approach (MOSA). Existing
modular Bus and mission payload architectures for typical commercial, civilian, and
military satellite systems will be discussed. The chapter provides space industry
views on “Open” versus “Close” interfaces design and addresses the challenges
associated with open interfaces using Open System Architecture (OSA) approach
using MOSA principles. The system interfaces discuss in this chapter include (i)
internal to satellite Bus and mission Payload (PL), (2) between satellite Bus and
mission payload, and (3) external to both satellite Bus and mission payload.
Typical commercial and civilian satellite systems take about 2–3years to build
and launch [1–4], while military systems take between 7 and 10 years [5, 6]. A
typical production flow for assembling and launching of a space vehicle is presented
in Ref. [6] and redrawn in Figure 1 as introduction steps for better understanding of
the design, build and launch of a satellite system. This chapter focuses on practical
design issues for satellite Bus’ and mission PL’s system/subsystem components
builds, and corresponding interface-design’s challenges associated with satellite Bus
integration, mission PL integration, and satellite system integration. A survey of
existing commercial, civilian and military satellite systems revealed that a typical
satellite Bus includes the following modular components [7–16]:
1
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 1.
A typical satellite system production flow (redrawn from [6]).
Similarly, our survey also revealed that a typical mission PL consists of the
following modular components [7–16]:
2
Future Satellite System Architectures and Practical Design Issues: An Overview
DOI: http://dx.doi.org/10.5772/intechopen.92308
• PL Subsystem 6—PL EPS: Existing PLs use power supply from satellite BEPS;
• PL Subsystem 8—PL ADCS: Existing PLs use ADCS from the satellite BADCS;
• PL Subsystem 9—PL PS:: Existing mission PLs use PS from the satellite BPS;
In practice, the above satellite Bus modular components can be found in the
following typical satellite Busses [7, 11, 14]:
For achieving optimum weight and power, existing satellite Bus and mission PL
are tightly coupled together with customized interface design. The industry trends
for the design and build of future satellite systems are moving toward OSA using
MOSA principles, in which the satellite Bus is loosely coupled with the mission PL
3
Satellite Systems - Design, Modeling, Simulation and Analysis
using “Open” and widely accepted interface standards. The key communication
linkage between a satellite Bus and a mission PL is the communication data Bus.
Currently, majority of satellite Busses employ the standard 1553 data Bus for data
communications among Bus components, and between the satellite Bus and mission
PL components. The communications over 1553 data Bus is limited to 1 Mega bit per
second (Mbps). Recently, there was an advanced development effort that was
funded by the U.S DOD to develop new 1553 standards called 1553 Enhanced Bit
Rate (EBR–1553) increasing the speed to 10 MB/s [17]. The EBR-1553 requires a
star/hub topology to provide the higher data rate and additional components to
implement the architecture. For data rates larger than 10 Mbps, space industry
trend is moving toward SpaceWire data Bus that was recently developed in Europe
for use in commercial satellites and scientific spacecraft [18].
The objective of this chapter is three-fold: (1) Provides an overview of existing
modular satellite Bus, mission PL architectures and related communication data
Busses, (2) Discusses future trends on the modular and open design and build of
satellite Bus and mission payload using MOSA principles, and (3) Addresses the
practical design challenges associated with “Modular” and “Open” design for future
satellite Bus and mission PL. The chapter is organized as follow: (i) Section 2
describes existing modular satellite Bus and mission PL architectures and related
communication data Busses; (ii) Section 3 presents industry view on “Open” and
“Close” interfaces for connecting satellite system components and existing popular
standards; (iii) Section 4 discusses the interface design challenges and provides an
overview of MOSA and related DOD Guidance and assessment tools for MOSA
implementation; (iv) Section 5 provides examples how to transition modular satel-
lite Bus and mission PL architectures to modular-and-open architectures using
MOSA implementation approach and tools in Section 4; and (v) Section 6 concludes
the chapter with remarks on the benefits associated with the proposed approach.
4
Future Satellite System Architectures and Practical Design Issues: An Overview
DOI: http://dx.doi.org/10.5772/intechopen.92308
Figure 2.
Overview of existing satellite Systems using standard 1553 data bus.
Figure 3.
Existing notional modular satellite bus architecture.
shows that space industry has used the modular design concept to architect the
satellite Bus, where common functions are group together and then isolate or
separate from the other group of functions. As an example, BAS consists of a
group of antenna components and control functions (e.g., antenna pointing,
beamforming, etc.), which is separated and isolated from BComRFS. It is important
to note that the figure also shows how these satellite Bus components are connected
together, i.e., the lines with arrows connecting them. These lines represent the
interfaces among the Bus components, where the interface can be any of the seven
interface types described above. Below is a list of some of the existing interfaces and
associated standards for existing satellite Bus based on National Aeronautical and
Space Administration (NASA), European Space Agency (ESA), U.S. DOD and
international Consultative Committee for Space Data System (CCSDS) standards
[19–26]:
5
Satellite Systems - Design, Modeling, Simulation and Analysis
◦ Satellite Bus shall protect its own electrical power system via overcurrent
protection devices on its side of the interface.
◦ Satellite Bus shall deliver a maximum transient current on any Power Feed
Bus of 100% (that is, two times the steady state current) of the maximum
steady-state current for no longer than 50 ms.
◦ Bus Survival Heaters, which are elements of the Bus thermal subsystem,
shall be required to have power to heat certain satellite Bus components
during off-nominal scenarios when the BEPS power is not fully energized.
◦ Power line conducted emissions for satellite Bus equipment shall meet the
EMC interface specification specified in SMC Standard Handbook, SMC-
S-008, Section 6, 6.01, 6.02, 6,03, 6.04, 6.05, 6.06, 6.07, and 6.08.
◦ Power line conducted susceptibility for satellite Bus equipment shall meet
the EMC interface specification specified in SMC Standard Handbook,
SMC-S-008, Sections 6, 6.10, 6.11, 6.12, 6.13, 6.14, 6.15, 6.16, 6.17, 6.18
and 6.19.
◦ ESD susceptibility for satellite Bus equipment shall meet the EMC
interface specification specified in SMC Standard Handbook, SMC-S-008,
Section 6, 6.43.
◦ EMP susceptibility for satellite Bus equipment shall meet the EMC
interface specification specified in SMC Standard Handbook, SMC-S-008,
Section 6, 6.45.
◦ Satellite Bus EPS should ground in a way that reduces introducing stray
currents or ground loop currents into the satellite Bus components.
• Typical Software & Data Interface Standards [19, 20, 23, 24, 25, 26]:
◦ Satellite Bus command and telemetry data formats shall be NASA Unified
S-Band (USB)/CCSDS standards or U.S. DOD Space-Ground Link
Subsystem (SGLS) standards. Note that (i) most of NASA and ESA
standards are CCSDS compliance for interoperability purpose, and (ii)
some military systems have both USB and SGLS capabilities.
6
Future Satellite System Architectures and Practical Design Issues: An Overview
DOI: http://dx.doi.org/10.5772/intechopen.92308
For military applications, majority of satellite Busses are usually designed using
contractor’s custom designed interfaces and very tightly couple together to reduce
weight, size and power. It is for this reason, current military satellite BTT&CS com-
ponent also include the COMSEC component. For commercial applications, satellite
developers are also concerned with weight, size and power reduction, but they are
also concerned with component refresh and upgrade without redesigning the satellite
Bus, hence commercial satellites tend to use modular Bus components and widely
accepted interface standards to connect the internal Bus components. Industry views
on the “open” and “close” interfaces will be addressed in Section 4.
Figure 4.
Existing notional modular Mission payload architecture.
7
Satellite Systems - Design, Modeling, Simulation and Analysis
design, the space industry has also applied the modular design concept to architect
the mission PL. Below is a list of some of the existing interfaces and associated
standards for existing mission PL leveraged from NASA, ESA, U.S. DOD and inter-
national CCSDS standards [19–26]:
◦ Sizing all components of the mission PL power harness, such as the wires,
connectors, sockets, and pins to the peak power level shall be required by
the mission PL equipment in addition to satellite Bus to prevent damage to
the power harnessing.
• Typical NASA Grounding Interface Standards [20, 22]: Similar to satellite Bus
discussed above but for mission PL.
◦ The mission PL thermal design should be decoupled from the satellite Bus
at the mechanical interface between the satellite Bus and neighboring
mission payload to the maximum practical extent.
• Typical Software & Data Interface Standards [19, 20, 23, 24, 25, 26]:
For most commercial applications, the MCC can be merged with the SOC, and
the mission PL TT&CS (PTT&CS) and PL CD&HS (PCD&HS) components can be
8
Future Satellite System Architectures and Practical Design Issues: An Overview
DOI: http://dx.doi.org/10.5772/intechopen.92308
incorporated into satellite (i) Bus TT&C (BTT&CS) and (ii) Bus CD&HS
(BCD&HS) components, respectively. Similar to the satellite Bus interfaces design,
for military applications, the mission PL components are tightly coupled using
contractor’s custom designed interfaces. For commercial applications, the mission
PL components are loosely coupled using widely accepted open interfaces.
Subsections 2.3.1 and 2.3.2 provide an overview of standard 1553 and SpaceWire
communication data Busses, respectively.
Existing commercial, civilian and military satellite data Busses have been using
Military Standard 1553B (MIL-STD-1553B) data Bus for communications among
satellite Bus and mission PL components. Figure 5 describes a typical MIL-STD-
1553B System [17, 27, 28]. This figure uses MIL-STD-1553B terminologies: (i) the
Bus Controller (BC) is considered as an Intelligent Terminal (IT) that is located in
the satellite mission computer, which is usually referred to as a Satellite Bus C&DH
component, and (ii) Remote Terminal (RT) is considered as a slave terminal that is
located in satellite platform components, which can be located in any satellite Bus or
mission PL components.
Figure 5 shows a typical commercial satellite system with RTs located in both
satellite Bus and mission PL components. As an Example, the RTs located in satellite
components are BAS, BADCS, BTCS and BTT&CS; and RTs located in the mission
PL components are PAS, PTCS, PDPS, PTRANSEC, and PFTS. For military appli-
cations, the Mission Computer (MC) can be located in both satellite Bus and mis-
sion PL, where the MC in the satellite Bus is responsible for all control functions
associated with the satellite operations and MC in the mission PL is responsible for
all control functions related to the mission PL operations.
Figure 5.
Typical civilian and commercial MIL-STD-1553B satellite Systems.
9
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 6.
Typical civilian and commercial SpaceWire satellite Systems.
Figure 7 presents the space industry view on open and close interface design.
This view separates the interface design into two categories, namely, Contractor
Proprietary Interface and Contractor Non-Proprietary Interface. Under this view,
the interface standards are then classified into two categories, namely, Preferred
and Non-Preferred Interface Standards. Based on this view, Section 3.1 defines open
interface design, and Section 3.2 defines close interface design. Section 3.3 provides
a list of existing popular open standards widely accepted by space industry.
1
LVDS is defined as Low Voltage Differential Signaling TIA/EIA-644, is a technical standard that
specifies electrical characteristics of a differential, serial communication protocol. LVDS Drivers use 80%
less current than current popular Pseudo Emitter-Coupled Logic (PECL) devices.
10
Future Satellite System Architectures and Practical Design Issues: An Overview
DOI: http://dx.doi.org/10.5772/intechopen.92308
Figure 7.
Industry view on open and closed interfaces design.
From Figure 7, the open interface design falls into the contractor non-
proprietary design category. For the interface design to be open, the interface
design shall not be contractor proprietary and that the interface shall use either
popular open interface standards widely accepted by space industry or open inter-
face standards with little market support and narrowly used by space industry.
Thus, a popular open interface design is a non-proprietary design that uses popular
open interface standard that is widely used by space industry. The benefits of open
interface design for the satellite buyers are (i) improving competition allowing
various space vendors (or contractor) to build open satellite Bus and mission PL
subsystem components, (ii) ease of refresh and technology upgrade allowing to
swap subsystem components without impacting the overall system, (iii) ease of
adapting to new requirements and operational threats, (iv) incorporating innova-
tion by allowing operational flexibility to configure and reconfigure a mission PL
quickly to meet rapidly changing operational requirements, (v) enabling cost saving
and cost avoidance during the design and sustainment phases by reusing technology
and Software/Hardware/Middleware (SW/HW/MW) components, and using
existing standardized HW/SW/MW parts and modules, and (vi) improving inter-
operability where severable HW/SW/MW modules can be changed independently.
As shown in Figure 7, the close interface design shall fall into contractor propri-
etary category. For an interface design to be close, it shall be contractor proprietary
and that the interface shall use either close interface standards with little market
support narrowly used by space industry or popular closed interface standards
widely used by space industry. Thus, a popular close interface design is a contractor
proprietary design that uses popular closed interface that is widely used by space
industry. The key benefits of close interface design are the potential reduction of
weight, size, power and manufacturing cost.
Based on Figure 7, the criteria for popular open standards are (i) publicly
available and widely used by both satellite Bus and mission PL vendors, (ii)
11
Satellite Systems - Design, Modeling, Simulation and Analysis
community and/or industry consensus-based that are matured and stable, and (iii)
technically adequate for all future commercial, civilian and military satellite sys-
tems. Following is a list of current popular standard organizations and widely
adopted open standards [18–31]:
• SQL for databases specified in ANSI ISO/IEC 9075–1, ISO/IEC 9075–2, ISO/IEC
9075–3, ISO/IEC 9075–4, ISO/IEC 9075–5.
• U.S. Space and Missile Systems Center (SMC) approved a project for
developing Common Payload Interface Specification (CoPaIS) standard for
satellite-to-payload Command and Data Handling (C&DH) interface intended
for all future SMC procured medium to large satellites [31].
12
Future Satellite System Architectures and Practical Design Issues: An Overview
DOI: http://dx.doi.org/10.5772/intechopen.92308
This section addresses the design challenges and is divided into four subsections,
including: (i) Subsection 4.1 discusses interface design and practical design issues;
(ii) Subsection 4.2 introduces MOSA concept; (iii) Subsection 4.3 presents DOD
MOSA guidance and the U.S. Naval Open Architecture (NOA); and (iv) Subsection
4.4 discusses MOSA tools and approach for MOSA implementation that addresses
the design issues identified in Subsection 4.1.
13
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 8.
MOSA and U.S. NOA approach.
U.S. DOD recommends OSA design using MOSA principles for future military
satellite system development with a goal to achieve a balance between business and
technical objectives that make a business sense in terms of (i) increase competition
and lower system acquisition cost, and (ii) lower sustainment cost over its life cycle.
MOSA design approach requires to implement five MOSA principles, including two
Business (B) and three Technical (T) principles [1, 2]. Figure 8 captures these five B
and T principles. Recently, U.S. Navy augmented MOSA principles with addition
five Naval Open Architecture (NOA) principles, including two business and three
technical principles as shown in Figure 8 [3].
Subsections 2.1 and 2.2 provide current implementation of the Technical princi-
ple 1 (T1) for the modular design of satellite Bus and mission PL, respectively. The
remaining Subsections 4.3 and 4.4 discuss the implementation of T2, T6, Business
principle 3 (B3), and B4 using DOD guidance for addressing the challenges
presented in Subsection 4.1.
2
Note that the term Open System Architecture (OSA) has also been used interchangeably with MOSA
by U.S. DOD.
14
Future Satellite System Architectures and Practical Design Issues: An Overview
DOI: http://dx.doi.org/10.5772/intechopen.92308
• Design the open system architecture using open interfaces. Implement the open
interfaces using open standards for connecting HW-to-HW, and SW-to-SW.
• The satellite system design shall accommodate growth and provide open
interface standards to allow future reconfiguration and addition of new
capabilities without large-scale redesign of the system.
• Develop a capability roadmap for the system covering the life of the system
following the completion of the rapid prototyping contract phase.
The above U.S. DOD’s guidance encourages the satellite system designers to
consider the above MOSA items in the design and build of the modular and open
15
Satellite Systems - Design, Modeling, Simulation and Analysis
satellite Bus and mission payload for future space systems. The following section
presents a proposal for assisting the satellite system designers to implement these
MOSA items along with assessment tools provided by U.S. DOD.
It is observed that the U.S. DOD, U.S. civilian agencies (e.g., NASA, NOAA, etc)
and U.S. satellite manufacturers/suppliers (e.g., LM, Boeing, Northrop Grumman
(NG), Raytheon, L3, etc) are investigating approaches for the modular and open
design and build of satellite Busses and mission PL’s using MOSA modular and open
design principles. Figure 9 proposes an approach to design and build of future
modular and open satellite Busses and mission PLs, and allowing the satellite buyers
to: (i) Buy the satellite Bus (see Path A of the figure) and mission PL (see Path B)
from different satellite manufacturers/suppliers, (ii) Have an option to choose a
third satellite vendor to integrate the satellite Bus and mission PL (see Path C).
The proposed MOSA implementation approach shown in Figure 4 consists of six
basic steps that are incorporated into three execution paths, namely, Path A, Path B
and Path C:
• Path A is for the satellite Bus manufacturer/supplier. This path has three basic
steps:
◦ Step I-A: Develop Modular satellite Bus Architecture (MoBA). The MoBA
subsystem components are described in Sections 1 and 2 (see Figure 3).
◦ Step II-A: Designate KOSS’s and select open standards for all internal
satellite Bus subsystem components. Open interface standards selection
and designation of KOSS are discussed in Sections 3 and 4.
Figure 9.
Proposed implementation approach for design and build of satellite Systems allowing buyers to acquire satellite
bus and Mission PL independently.
16
Future Satellite System Architectures and Practical Design Issues: An Overview
DOI: http://dx.doi.org/10.5772/intechopen.92308
◦ Step III-A: Design and build Open Modular satellite Bus System (OMoBS).
This step is achieved by identifying all potential KOSS’s from the satellite
Bus to any mission PL’s, i.e., the selected KOSS’s should be independent of
mission types. The satellite Bus manufacturer is responsible for
integrating all Bus components and have the satellite Bus ready for sale.
• Path B is for the mission PL manufacturer/supplier. This path also has three
basic steps that are similar to Path A:
◦ Step II-B: Designate KOSS’s and select open standards for all internal
Mission PL subsystem components. Open interface standards selection
and designation of KOSS for mission PL are also discussed in Sections 3
and 4.
◦ Step III-B: Design and build Open Modular Mission PL System (OMoPS).
This step is achieved by identifying all potential KOSS’s from the any
mission PL’s to satellite Bus, i.e., the selected KOSS’s should be
independent of mission types. The mission PL manufacturer is responsible
for integrating all mission PL components and have the PL ready for sale.
• Path C is for the satellite system integrator. This path has additional three new
steps:
◦ Step IV: The system integrator works with satellite Bus and mission PL
manufacturers to develop a satellite system interface specification
specifying all “open” and “close” interfaces between the mission PL-and-
satellite Bus. All open interfaces between the mission PL-and-satellite Bus
shall be selected to meet the business and performance objectives
approved by the buyer. The system integrator performs satellite Bus and
mission PL integration using the approved interface specification.
◦ Step VI: System integrator delivers the satellite system to the buyer.
DOD has also developed MOSA tools to assist MOSA implementation and
assessment of military satellite Bus and mission PL “Openness”. These tools can also
be used for civilian and commercial applications. The DOD tools include MOSA
Program Assessment and Rating Tool (PART), Open Architecture Assessment Tool
(OAAT), and Key Open SubSystem (KOSS) Tool:
3
PART can be found from: https://www.dau.mil/cop/mosa/Lists/Tools/DispForm.aspx? ID=2&Conte
ntTypeId =0x01002BC08FCA204040449CF11CB472BEEE1800AA6D1BC9926604469A02D
DB936F94D1F
17
Satellite Systems - Design, Modeling, Simulation and Analysis
• MOSA KOSS Tool5: One of the key MOSA principles is the Business Principle
number 4, namely, Designate Key Interfaces (see Figure 8, B4). The
identification of KOSS’s is an important task in realizing open systems. This
MOSA principle requires the system designers to compromise between cost
and performance by selecting a set of KOSS’s with their associated interfaces
that can be assigned widely used open standards allowing for easy and
affordable update and frequent refresh. MOSA KOSS tool provides guidance
for KOSS’s identification and selection. The tool makes use of system capability
road map, system requirements and Subject Matter Expert (SME), program’s
sponsor and warfighter knowledge to identify the system/subsystem
components expected to have a high volatility over the system life cycle. The
tool specifies the key interfaces as those either side of volatile components. The
tool will help the satellite system designer to identify and rank KOSS’s
components that will meet both programmatic and technical requirements.
This section demonstrates how to use Steps II-A and II-B of the proposed MOSA
implementation approach presented in Section 4.4 for the design and build of
future resilient and robust satellite systems. Subsections 5.1 and 5.2 present poten-
tial modular-and-open satellite Bus and mission PL architectures, respectively.
4
OAAT from: https://www.dau.mil/cop/mosa/Lists/Tools/DispForm.aspx?ID=1&ContentTypeId=
0x01002 BC08FCA204040449CF11CB472BEEE1800AA6D1BC9926604469A02DDB936F94D1F
5
KOSS from: https://acc.dau.mil/adl/enUS/317012/file/46502/KOSS%20Overview_FINAL_5Aug09.pdf.
18
Future Satellite System Architectures and Practical Design Issues: An Overview
DOI: http://dx.doi.org/10.5772/intechopen.92308
Modular satellite Modular satellite Bus subsystem Recommendation for open interface
Bus subsystem and component description Standardization (potential KOSS)
component
BTT&CS-6 SGLS Base Band Signal Processing Recommend for open interface
(USB Mode1, 2) standardization
19
Satellite Systems - Design, Modeling, Simulation and Analysis
Modular satellite Modular satellite Bus subsystem Recommendation for open interface
Bus subsystem and component description Standardization (potential KOSS)
component
BEPS-3 Solar Array Drive Assembly (SADA) Recommend for open Interface
standardization.
BEPS-5 Bus Power Regulation Unit (BPRU) Recommend for open Interface
standardization.
BEPS-6 Fuse Box Assembly (FBA)
.
BEPS-7 Pyro Relay Assembly (PRA)
Table 1.
Satellite bus subsystems decomposition and potential KOSS.
This subsection provides an example for the transition of the notional modular
mission PL architecture presented in Figure 4 to a modular-and-open mission PL
architecture. Table 2 summarizes the decomposition and KOSS analysis results for
four mission PL subsystems, including PAS (PL Subsystem 1), CPCom-RFS (PL
Subsystem 2), PDPS (PL Subsystem 3) and PFTS (PL Subsystem 11).
The mission PL digital processing system is not recommended for interface
standardization due to many variations between systems and subsystems. Multi-RF
Wideband RX Up/Down Converters and Tunable IF Down Converters require a
significant amount of analysis or NRE hours and are also not a good candidate for
standardization. Again, DOD KOSS tool should be used to finalize the KOSS analysis
results presented here for actual design and build of the satellite systems.
6. Conclusion
The chapter provides an overview of existing modular satellite Bus and mission
PL architectures and associated standards for communication data Busses. The
chapter defines open and close interfaces along with industry approved popular
standards and discusses the interface design challenges. Moreover, the chapter pro-
vides an overview of MOSA and related DOD guidance and assessment tools to
address the interface design challenges. Examples for the design and build of
modular-and-open satellite Bus and mission PL architectures are also presented.
The intent of this chapter is to provide an innovative approach for the satellite
system designer to design and build of the next generation satellite achieving a
balance between business and technical objectives that make a business sense for
both the satellite manufacturers and buyers in terms of lower system acquisition
and sustainment costs over its life cycle. The MOSA implementation approach
presented here allows the satellite manufacturers to build the satellite Bus and
mission PL separately for more production, flexibility, and market competition.
Concurrently, the approach also allows the satellite buyers to buy satellite Bus at
high volume with reduced unit costs and less schedule risk. Another benefit for the
20
Future Satellite System Architectures and Practical Design Issues: An Overview
DOI: http://dx.doi.org/10.5772/intechopen.92308
PFTS-1 Atomic Clock Unit (ACU) Not recommended for interface standardization;
ACU and CM&CU will vary depending on
PFTS-2 Clock Monitoring & Control
mission type and mission requirements
Unit (CM&CU)
Table 2.
Mission payload subsystems decomposition and potential KOSS.
Acknowledgements
Although the preparation of this work was not funded by The Aerospace Cor-
poration, but the author acknowledges the work presented in this chapter was based
21
Satellite Systems - Design, Modeling, Simulation and Analysis
on his knowledges accumulated over the years from many space programs at The
Aerospace Corporation, Raytheon and Jet Propulsion Laboratory. In addition, the
author would like to express his appreciation to Aerospace’s manager, Ms. Navneet
Mezcciani, for her professional support.
Conflict of interest
The preparation of this chapter was not funded by the Aerospace Corporation,
and it was done by the authors using his own time and resources, thus it does not
represent the Aerospace Corporation’s view on the DOD guidance for MOSA
implementation and the proposed system architecture solutions.
Notes/Thanks/Other declarations
The author wishes to thank his wife, Thu-Hang Nguyen, for her enormous
patience and boundless support during the preparation of this chapter.
Author details
Tien M. Nguyen1,2
© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License (http://creativecommons.org/licenses/
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
22
Future Satellite System Architectures and Practical Design Issues: An Overview
DOI: http://dx.doi.org/10.5772/intechopen.92308
References
23
Satellite Systems - Design, Modeling, Simulation and Analysis
[16] Collins R. Momentum and Reaction [24] Consultative Committee for Space
Wheels. Available from: https://www. Data Systems (CCSDS).
rockwellcollins.com/Data/Products/ Recommendations for Space Data
Space_Components/Satellite_Stabiliza Systems, Packet Telemetry, CCSDS
tion_Wheels/RSI_12_Momentum_and_ 102.0-B-5, Blue Book. 2000. Available
Reaction_Wheels.aspx from: https://public.ccsds.org/Pubs/
102x0b5s.pdf#search=packet%20
[17] Gwaltney DA, Briscoe JM. telecommand
Comparison of Communication
Architectures for Spacecraft Modular [25] Consultative Committee for Space
Avionics Systems. NASA/TM—2006– Data Systems (CCSDS).
214431. June 2006. NASA Center for Recommendations for Space Data
AeroSpace Information Systems, Telecommand – Part 1:
Channel Service, CCSDS 201.0-B-3,
[18] SpaceWire—Links, Nodes, Routers, Blue Book. 2000. Available from:
and Networks, ECSS–E–50–12A. 24 https://public.ccsds.org/Pubs/201x0b3s.
January 2003. European Cooperation for pdf#search=packet%20telecommand
Space Standardization
[26] Space-Ground Link Subsystem
[19] NASA Common Instrument (SGLS). Ground Station, System
Interface Project. Hosted Payload Analysis Summary Report, WDL-
Guidelines Document. Document TR3227-1. Vol. 1. 15 November 1968.
Number: CII-CI-0001, Version Rev. A., Philco-Ford Corporation. Available
4 November 2013. Earth System Science from: https://apps.dtic.mil/
Pathfinder Program Office dtic/tr/fulltext/u2/853122.pdf
24
Future Satellite System Architectures and Practical Design Issues: An Overview
DOI: http://dx.doi.org/10.5772/intechopen.92308
25
Chapter
Abstract
1. Introduction
1
Satellite Systems - Design, Modeling, Simulation and Analysis
payload of 3-band multispectral imager of 30-m resolution, which was intended for
wide-swath land coverage imaging. After the first constellations decommissioned,
it built the second generation with better resolution (20 m). The first launch of
DMC-2 constellation was done in 2009 [3].
Since 2013, Skybox/Skysat has deployed 15 satellites that carry 1-m panchromatic
imager and 2-m 4-band multispectral imager [4]. Unlike DMC, which mission objec-
tives are to observe wide areas with nadir pointing scanning mode, it aims to provide
frequent repeat very high resolution images using massive numbers of highly maneu-
verable satellites. Another commercial Earth observation microsatellite constellation
mission is prepared by Axelspace. The company planned to have 50 satellites launched
starting 2017. The satellite carries imager with 2.5-m panchromatic and 5-m multi-
spectral [5, 6]. Figure 1 shows the configurations of the Skybox and Grus satellites,
which show that Skybox uses single lens and parabolic data downlink antenna, while
Grus uses two lenses and horn-type data downlink antenna.
In addition to Earth observation missions, microsatellite constellation also
being built for Low Earth Orbit (LEO) telecommunication mission. OneWeb and
Telesat are two companies that will launch hundreds of microsatellites in coming
years [7, 8].
The use of microsatellites for commercial purposes means that the technol-
ogy is mature enough to ensure good return-of-investment. One of the major
aspects that contribute to the success of microsatellite technology is its system
design. Therefore, the objective of this chapter is to provide insight into micro-
satellite system design. The chapter addresses the question related to limitation
in weight and size, and how the satellite designer manages to meet the mission
requirements.
Out of many microsatellites developers, two system designs of microsatellites,
namely Technical University (TU) Berlin heritage and University of Surrey heri-
tage, are selected for comparison in this chapter, due to their very different design
approaches. To be comparable, the choices of microsatellite system to be compared
are the ones that manufactured about the same time, so that the technology avail-
able is mostly the same. The microsatellites also have to have in-orbit experience, so
its design success can be measured. Data mining resulted that the satellite operation
year chosen is between 1999 and to date. For TU Berlin system, the choices are DLR-
TUBSAT, MAROC-TUBSAT, Indonesian LAPAN-TUBSAT, LAPAN-ORARI, and
LAPAN-IPB. Meanwhile, for University of Surrey system, the choices are Korean
KITSAT-3, STSat-1 and STSat-3, as well as Turkish BILSAT-1 and RASAT.
Figure 1.
Google Skybox satellite and Axelspace’s Gruz satellite design.
2
System Designs of Microsatellites: A Review of Two Schools of Thoughts
DOI: http://dx.doi.org/10.5772/intechopen.92659
This chapter is divided into five sections, with the first section introducing the
background and objectives of the chapter. The second section explains how the
satellite design samples for the University of Surrey heritage were selected, and
what satellite design parameters were used in the comparison. Section 3 displays
the satellite design parameters for TU Berlin heritage. Section 4 provides analysis
from the comparison of the two-design heritage, in term of parameters noted in the
previous two sections. Section 5 summarizes the analysis and provides recommen-
dation for further studies regarding the subject.
University of Surrey is known as one of the pioneers in the design and build of
microsatellite in the 1990s. It started launching microsatellite in 1991 with amateur
KITSAT-3 BILSAT-1
Operation 1999–2003 2003–2006
Bus
Payload
3-band imager w/ 570-mm lens 2 × 3-band imager w/150-mm
lens
CMG
Size (cm) 50 × 60 × 85 60 × 60 × 60
Table 1.
Sample for the University of Surrey microsatellite system design.
3
Satellite Systems - Design, Modeling, Simulation and Analysis
radio missions. To simplify the satellite design, the first microsatellite generation
has passive attitude control system, that is, using gravity gradient telescopic boom.
The university provided microsatellite development and building capabilities
to many emerging space countries, including Thailand, Malaysia, South Korea,
Algiers, Turkey, and Nigeria. At the time, such countries started to use remote sens-
ing satellites, mostly from the United States and European, for various land-based
applications. Therefore, they required remote sensing payloads to include in their
satellite missions. Such mission elevates the design requirements to active attitude
control system and higher data rate downlink system.
Thailand’s Mahanakorn University collaborated with the University of Surrey
to jointly develop TMSat that was launched in 1998 [9]. TMSat focuses on remote
Horizon sensor — — —
Magnetotorquer/meter 3-axis/3-axis fluxgate 3-axis/3-axis 3-axis/2 × 3-axis
fluxgate
TTC S-band S-band S-band (primary) and
UHF/VHF (emergency)
Table 2.
Sample for the University of Surrey microsatellite heritage system design.
4
System Designs of Microsatellites: A Review of Two Schools of Thoughts
DOI: http://dx.doi.org/10.5772/intechopen.92659
sensing and amateur radio mission. Since Thailand did not continue building its
subsequent satellites, TMSat is not selected as satellite design heritage sample in
this chapter.
Singapore’s Nanyang Technology University (NTU) collaborated with
the University of Surrey to jointly develop satellite subsystem for UoSAT-12.
However, the satellite is not a microclass and therefore is not selected as a
sample for the University of Surrey’s satellite system design in this chapter.
The satellite subsystem from NTU is a communication payload with S-band
downlink and L-band uplink, which provides the Internet protocol communica-
tion operating at 1 Mbps. Since the experience with the University of Surrey
only in subsystem design and development, the subsequent NTU satellite,
that is, XSAT, is also not considered as the University of Surrey heritage
satellite [10, 11].
South Korean experience with the University of Surrey satellite design is
when Satellite Technology Research Center (SaTReC), an institution under Korea
Advanced Institute of Science and Technology (KAIST), jointly built KITSAT-1
and KITSAT-2 and launched it in 1992 and 1993. Both satellites have store-forward
communication amateur payload and low-resolution imagers. Since the KITSAT-1
and KITSAT-2 development time does not match with other microsatellite design
sample, only the design of KITSAT-3 is used in this chapter. SaTReC then devel-
oped STSAT series as its second generation microsatellites. Since STSAT-2 expe-
rienced launch failure, only STSAT-1 and STSAT-3 are selected as satellite design
samples [12–15].
Turkey’s experience with the University of Surrey satellite design is when its
space research institute, TUBITAK-UZAY (previously named BILTEN TUBITAK-
ODTU), jointly developed BILSAT-1. The satellite was part of DMC-1 constellation
[16–19]. After BILSAT-1, the institute then built its second generation microsatel-
lite, RASAT. Therefore, both microsatellites are used as sample for the University of
Surrey design heritage [20–23].
Fifteen satellite bus design parameters are selected for the comparison, including
14 mechatronics component parameters in the satellites’ design. For the University
of Surrey satellite heritage, the parameters are tabulated in Tables 1 and 2. Structure
design from four of the five microsatellites is shown in Figures 2 and 3. Payload
parameters also noted in Tables 1 and 2 to explain the similarity (or differences) in
the mission requirements and their impacts to satellite bus parameters. The weight
and dimensions are, in additional of drawings, noted in to explain the satellite
Figure 2.
Mechanical design of KITSAT-3 and STSAT-3.
5
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 3.
Mechanical design of BILSAT-1 and RASAT.
structure design aspects. The satellite operation years are noted in the tables to show
the context of available technology.
As shown in Figures 2 and 3, the University of Surrey heritage satellites use
electronic trays for its satellite bus electronics. The aluminum trays also function as
load bearing structure, so that the rest of the satellites components, such as reaction
wheels and attitude sensors, can be laid out around them. After all components
integrated, the solar panels and/or other outside panels that are made of lighter
materials can be used to cover the satellites.
6
System Designs of Microsatellites: A Review of Two Schools of Thoughts
DOI: http://dx.doi.org/10.5772/intechopen.92659
DLR-TUBSAT MAROC-TUBSAT
Launch 1999–2007 2001–2006
Bus
Thruster — —
Star sensor — IRE
Sun sensor 4 panels +1 cell 6 single cell
Horizon sensor — —
Magnetotorquer/meter 1 axis coil + 1 rod/− 1 axis/3-axis sensor
Table 3.
Sample for the Technical University of Berlin microsatellite system design.
4. Analysis
Tables 1 and 2 show that the Korean satellites have employed deployable solar
panel (which is also shown in Figure 1), since the mission required high power and
used direct energy transfer (DET) mode. Such approach is very much different than
those used by KITSAT-1 and KITSAT-2, which have body-mounted solar panels.
On the other hand, Turkish satellites use body-mounted solar panels and therefore
do not have the requirement of one side of the satellite always facing the sun for
battery charging.
Tables 3 and 4 show that all TU Berlin heritage use body-mounted solar panels. It
uses Si panels for its first three satellites, then opted to higher capacity GaAs panels in
LAPAN-ORARI and LAPAN-IPB. Generally, the power budget for the University of
Surrey heritage satellites is higher than the TU Berlin heritage, even in the ones with
body-mounted solar panels. As shown in Figure 5, in LAPAN-IPB, one of the sides has
two 46 × 26 cm solar panels. The side is projected to be Sun pointing most of the time.
Battery chosen to be used in the early University of Surrey heritage satellite
design is NiCd, while in TU Berlin’s satellite design is NiH2. NiCd batteries require
7
Satellite Systems - Design, Modeling, Simulation and Analysis
Payload
Color video cam. Color video cam. 4-band imager
w/50 mm lens w/1000 mm lens w/300 mm lens
Table 4.
Sample for the Technical University of Berlin microsatellite heritage system design.
8
System Designs of Microsatellites: A Review of Two Schools of Thoughts
DOI: http://dx.doi.org/10.5772/intechopen.92659
Figure 4.
Mechanical design of DLR-TUBSAT and MAROC-TUBSAT.
Figure 5.
Mechanical design of LAPAN-TUBSAT and LAPAN-IPB.
9
Satellite Systems - Design, Modeling, Simulation and Analysis
Tables 1 and 2 show that the University of Surrey satellite design heritage uses
separate attitude control computer that integrates attitude sensors, including sun
and star sensors with all reaction wheels and gyros. This is done so that the attitude
control system can work in closed loop all the time. Such approach is necessary for
the microsatellite design with deployable solar panels, such as KITSAT-3, STSAT-1,
and STSAT-3 since failure of sun pointing could be disastrous for the satellite. As
shown in Tables 3 and 4, in the TU Berlin satellite design heritage, none of the sat-
ellites have separate attitude control computer. In the design, each reaction wheel-
gyro pair directly connected to the main computer, and therefore, closed loop with
star and sun sensors can only be done using the main computer resources.
Differences are also found in the attitude control sensor between the University
of Surrey design heritage. The Korean microsatellites use fiber-optic gyro, while the
Turkish microsatellites use MEMS gyro. Meanwhile, in all TU Berlin microsatellites,
fiber-optic gyros are used.
For attitude control actuators, all the selected satellites use reaction wheels and
air coils for angular momentum dumping/generation. Figures and data showed
that TU Berlin heritage satellites use reaction wheels in 3-axis configuration. For
LAPAN-ORARI and LAPAN-IPB satellites, they used redundant wheel at satellite
major inertia axis that noted as 3 + 1 as shown in Table 4. For the University of
Surrey heritage satellites, only KITSAT-3 uses reaction wheels in 3-axis configura-
tion. The rest of the satellites uses tetrahedral configuration (noted as 4 as shown in
Table 1).
The TU Berlin’s attitude control design was chosen to reduce computational
burden for filtering out reading noise/jitter in the attitude control sensors. The TU
Berlin heritage satellites offer two options for attitude control mode, in addition to
regular closed loop, including (1) interactive mode for the satellite with video cam-
era payload, such as DLR-TUBSAT and LAPAN-TUBSAT, and (2) angular momen-
tum management mode for the satellite with line imagers, such as Maroc-TUBSAT
and LAPAN-A3. The angular momentum management mode is supported by their
structure design, that is, solid aluminum box, which created maximum inertia
properties at 1 axis and very little cross-product inertias [30, 31]. Such design has
been successfully performed highly stable open-loop angular momentum manage-
ment operation as published by Utama [31] and Mukhayadi [32].
From a selected set of satellite designs shown in Tables 1–4, only BILSAT-1
and STSAT-3 have thrusters. The objective for BILSAT-1 thruster is to maintain the
satellite orbit separation in the constellation, so that the image coverage could be
optimized. In STSAT-3, the plasma thruster is part of in-orbit qualification process
for the low power plasma thruster technology developed by KAIST.
10
System Designs of Microsatellites: A Review of Two Schools of Thoughts
DOI: http://dx.doi.org/10.5772/intechopen.92659
4.5 TTC
4.6 Payload
The payload profiles for both satellite design heritages showed that the platforms
are suitable for both Earth observation, science, and low data rate communication
missions. All of the selected satellites, except Korean STSAT-1, are Earth observa-
tion missions, which are considered important by stakeholder of satellite developer
in Korea, Turkey, and Indonesia. KITSAT-3, BILSAT-1, RASAT, and LAPAN-IPB
are for land cover that can be applied for estimating crop yield. The payload data
showed that combining mission is typical for microsatellite applications. The multi-
band infrared (MIRIS) payload in STSAT-3 is used for Earth and space observation.
LAPAN-ORARI has three kinds of missions, including Earth observation, commu-
nication, and ship data collecting platform.
The quantity and quality of the payload in Tables 1–4 showed that mission data
are increasing with time, which increase the required downlink data rate. For the
University of Surrey heritage, the data rate started with 3 Mbps in KITSAT-3 and
increased to 100 Mbps in RASAT. For the TU Berlin heritage (the digital transmis-
sion cases), the data rate started with 256 kbps in Maroc-TUBSAT and increased
to 100 Mbps in LAPAN-IPB. In the early missions, the mission data downlink is
transmitted in S-band, and as the data rate requirement increases, the downlink has
been shifted to X-band.
11
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 6.
Microsatellites’ weight versus launch year.
Figure 7.
Microsatellites’ density versus launch year.
Figure 6 shows the weight of each microsatellite sample. It shows that the weight
of TU Berlin heritage satellites grows in time. This is due to the increase in mission
quantity and complexity, which therefore requires more components in the satel-
lites (bigger batteries, more attitude sensors, larger lens for imager payload, etc.).
For the University of Surrey satellites heritage, such pattern is not found. The den-
sity of (weight/volume) the satellites is shown in Figure 7, indicating that the TU
Berlin heritage satellites are more compact than the University of Surrey heritage
satellites. For the University of Surrey satellites heritage, the design uses maximum
volumetric envelope for maximizing the solar panel area.
5. Conclusions
The chapter has discussed the differences between the University of Surrey
design heritage microsatellites and the TU Berlin heritage microsatellites. Five
12
System Designs of Microsatellites: A Review of Two Schools of Thoughts
DOI: http://dx.doi.org/10.5772/intechopen.92659
sample satellites from each satellite design heritage are compared, including 15 bus
parameters, payload profiles, and satellite weight and volume at launch. From the
comparison, it is found that major differences in the satellite bus are in the choice of
main computers and their associated link configuration and in the attitude control
modes that also affect the design. Another major difference is in the satellites’
structure design, which resulted in much higher density in the TU Berlin heritage
satellites than the University Surrey heritage satellites. In the early design, there are
differences in the choice of satellite’s batteries. However, as soon as Li-ion batteries
became available, both design heritages used such technology. In answering the
increasing needs in payload data handling, both design heritage use FPGA-based
payload data handling and high downlink data rate in X-band. GPS is also the
technology adopted by both design heritages for orbit determination and imager’s
ancillary data.
For further studies on the topic, it is suggested that comparison to be done on
the power budget of the satellites and on the operation performance parameters of
the satellites with similar missions.
Acknowledgements
The author wishes to acknowledge the Center for Aerospace Policy Studies of
LAPAN for funding this publication and also wishes to thank the book editor for
giving substantial advice in writing this chapter.
Author details
Triharjanto Robertus
National Institute of Aeronautics and Space, Jakarta, Indonesia
© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License (http://creativecommons.org/licenses/
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
13
Satellite Systems - Design, Modeling, Simulation and Analysis
References
[1] Swartout M. Cheaper by the dozen: [10] Sun W, Sweeting MN. In-orbit
The avalanche of rideshares in the results from UoSAT-12 Earth
21st century. In: 2013 IEEE Aerospace observation minisatellite mission. In:
Conference; Big Sky, MT, USA; 2013 Proceedings of the 3rd International
Symposium of IAA; Berlin; 2-6 April
[2] Bunchen E. 2014 Nano/Micro- 2001. pp. 79-82. IAA-B3-0305P
satellite market assessment. In:
Proceedings of the 28th AIAA/USU [11] Fouquet M, Sweeting M. UoSAT-12
Conference on Small Satellites. Utah; Minisatellite for high performance
2014 earth observation at low cost. Acta
Astronautica. 1997;41(3):173-182
[3] Kramer HJ. DMC-3. 2016. Available
from: https://directory.eoportal.org/ [12] Lim JT, Nam MR, Ryu KS, Tahk KM,
web/eoportal/satellite-missions/d/ Rhee SH, Kim KH. Exploring space on
dmc-3 a small satellite, STSAT-2: A test bed
for new technologies. In: 14th Annual
[4] Kramer HJ. SkySat. 2016. Available AIAA/USU Conference on Small
from: https://directory.eoportal.org/ Satellites; 2000
web/eoportal/satellite-missions/s/skysat
[13] Kim BJ, Lee HW, Choi SD. Three-
[5] Axelspace. GRUS-1: The Next axis reaction wheel attitude control
Generation of Earth-Observation system for KITSAT-3 microsatellite.
Satellites. 2016. Available from: https:// Space Technology. 1997;16(5/6):191-296
www.axelspace.com/en/solution_/grus/
[14] Park JO, Park HY, Lee DH,
[6] Quigley JT. Japanese Startup Raises Lee JH, Rhee SW. STSat-3 (Science
$15M to Build Satellites that Are the Size & Technology Satellite-3) satellite
of your Air Conditioner. 2015. Available missions. In: Proceedings of the 60th
from: https://www.techinasia.com/ IAC (International Astronautical
japan-micro-satellite-startup-axelspace- Congress); Daejeon, Korea; 12-16
series-a October 2009. IAC-09.B4.2.6
[7] Henry C. How OneWeb Plans to [15] Shin GH, Park HY, Seo JK,
Make Sure Its First Satellites Aren’t Chung TJ, Myung NH. Introduction to
Its Last. 2019. Available from: https:// micro-satellite STSAT-3 development.
spacenews.com/how-oneweb-plans-to- In: Proceedings of the 60th IAC
make-sure-its-first-satellites-arent-its- (International Astronautical Congress);
last/ Daejeon, Korea; 12-16 October 2009.
IAC-09.B4.6A.2
[8] Clark S. Telesat Taps Blue Origin to
Launch Broadband Satellite Fleet. 2019. [16] Gomes LM, Yuksel G, Lappas V, da
Available from: https://spaceflightnow. Silva Curiel A, Bradford A, Ozkaptan C,
com/2019/02/04/telesat-taps-blue-origin- et al.. BILSAT: Advancing smallsat
to-launch-broadband-satellite-fleet/ capabilities. In: AIAA/USU Conference
on Small Satellites; Logan, UT; 11-14
[9] Triharjanto R. Strategic environment August 2003
and implementation of satellite
technology acquisition programs. In: [17] Bradford A, Yuksel G, Gomes LM,
ASEAN, Proceeding Seminar Nasional Ozkaptan C, Sweeting M, Orlu U.
Kajian Kebijakan Penerbangan dan BILSAT-1: A low-cost, agile, earth
Antariksa; Surabaya, Indonesia; 2018 observation microsatellite for
14
System Designs of Microsatellites: A Review of Two Schools of Thoughts
DOI: http://dx.doi.org/10.5772/intechopen.92659
15
Satellite Systems - Design, Modeling, Simulation and Analysis
16
Chapter
Abstract
1. Introduction
1
Satellite Systems - Design, Modeling, Simulation and Analysis
a failure occurs, the entire electronic system will fail. Thus, the centralized manage-
ment method is no longer suitable for the needs of spacecraft development.
The satellite intelligent avionics system is an information processing and trans-
mission system that uses computer network technology to interconnect satellite-
borne electronic equipment on the satellite to achieve internal information sharing
and comprehensive utilization, function integration, and resource reorganization
and optimization. Utilizing onboard computers to complete satellite data manage-
ment, control management, communication management, time management,
energy management, and job management functions through unified scheduling of
satellite missions. Its essence is the generation, identification, processing, analysis,
transmission, and distribution of information process. The integrated satellite
electronic system integrates the functions of the satellite platform electronic
equipment, and its design level directly determines the performance of the satellite
platform [3–5].
At present, satellite sub-systems mostly adopt independent design schemes,
which decentralize satellite attitude control, propulsion control, thermal control,
satellite-ground link communication, and power control functions. The onboard
computer is responsible for tasks such as remote control, telemetry, program-
controlled operation, thermal control, and time management. The attitude and
orbit control computer are responsible for attitude and orbit (including propulsion
control) control. Each sub-system such as power supply, thermal control, and
digital transmission is equipped with corresponding lower-level computers respon-
sible for telemetry acquisition and remote control of the respective sub-system.
However, the satellite system designed using this approach is usually resulting in
heavy weight, high power consumption, large volume (aka high size, weight, and
power (SWAP)), complex interface relationships, weak system reconfiguration
capabilities, and low functional density. In order to overcome the abovementioned
shortcomings and make the satellite avionics system better meet the SWAP and
flexible system configuration requirements of future missions, it is necessary to
improve its design technology, that is, from the current independent design of each
sub-system to the open and modular design of the entire satellite. Based on the prin-
ciple of unified application, deployment and operation of hardware resources, and
the full use of the various functions of the software, the information sharing of the
entire satellite, simple system configuration, and overall performance optimization
are realized.
This chapter focuses on the application requirements of the new generation of
intelligent avionics system for future communication satellites, and adopts an open
architecture of “centralized management, distributed measurement and drive, and
software and hardware modular design.” The universal, standardized, and scalable
intelligent avionics system is built based on the basic modular elements of open
hardware modules, open software components, and industry standardized internal
and external busses.
2. System structure
2
Design of Intelligent and Open Avionics System Onboard
DOI: http://dx.doi.org/10.5772/intechopen.93141
processing, storage, and radiation resistance. Section 2.3 describes the dynamic
state reconfigurable task scheduling that improves the fault tolerance ability of the
satellite network in view of the typical scenarios of the satellite integrated elec-
tronic system in the operation process. Section 2.4 discusses the design of software
partition protection mechanism related to the next-generation avionics system
and analyzes the requirements, design, and functions of partition protection,
aiming to improve the robustness of the software system.
Figure 1.
Avionics system architecture diagram for satellite communication system.
3
Satellite Systems - Design, Modeling, Simulation and Analysis
information core of the whole satellite, especially for the requirements of intelligent
satellite systems. Research on the realization of high-performance computing of
avionics systems is an inevitable requirement [6, 7]. In order to improve the com-
puting capacity of the avionics system, a high-performance onboard processor is
utilized. The following introduces onboard high-performance computer from four
aspects: main functions, processing, storage, and anti-irradiation.
i. Main functions:
• Uses interface with external unit. Provide a brief description why is this
important function.
• L1 instruction cache and L1 data cache with Error Correcting Code (ECC)
function
4
Design of Intelligent and Open Avionics System Onboard
DOI: http://dx.doi.org/10.5772/intechopen.93141
The two typical scenarios usually encountered by satellite avionics systems during
operation are (a) a node fails or requires functional reorganization so that some tasks
on this node need to be migrated to other nodes through the network and (b) the
resource occupancy rate of a node is too high so that some tasks on this node will be
migrated to other relatively idle nodes for execution. The avionics system is designed
with networked real-time multitasking distributed system software, which can also
implement dynamic reconfiguration of functions and task scheduling. The embed-
ded system software running on each node in the network supports not only the local
real-time multitasking scheduling but also the network operation capability. The avi-
onics system supports function modification and function migration between nodes,
which realizes the transformation to software-defined satellite functions, reduces the
differences in hardware products, improves the fault tolerance of the intra-satellite
network, and also meets the growing needs of intra-satellite networking [8].
The avionics system networked real-time multitasking distributed system
software is shown in Figure 2 and has the following characteristics:
Figure 2.
Networked real-time multitasking distributed system software.
5
Satellite Systems - Design, Modeling, Simulation and Analysis
composition capability, that is, when the business needs of the day change and
the service call is adjusted, can support the user to quickly combine services
and form a new business process [9].
• The middleware layer serves as a transition layer between the above and the
next [10]. By combining the characteristics of the tasks of each layer, the
corresponding theme is designed, and the publish-subscribe technology is
used to provide the entire application layer with access to various resources in
the basic resource layer. The software layer implements operating system and
communication protocol level shielding through the packaging of thread tasks,
synchronization resources, memory access, IO operations, Ethernet, shared
memory, and fiber-optic communications, providing access to the hardware
device layer. By virtualizing the calculation, storage, and network resources
of the basic resources, the basic resources as a service are realized, and the
availability and scalability of the hardware resources are guaranteed. At the
same time, the software is dynamically deployed for the hardware of the basic
resources, including the automation of basic software and application software
installation settings, maintenance and upgrades, etc., and provides the system
with general basic services such as system reconstruction, software fault toler-
ance, data management, subscription release, etc.
6
Design of Intelligent and Open Avionics System Onboard
DOI: http://dx.doi.org/10.5772/intechopen.93141
• Partition time protection. Its strategies include the following: (1) The basic
unit of scheduling is partitioning, and partition scheduling has no priority.
The Main Time Frame (MTF) is used to statically define the scheduling
order of each partition and the proportion window size. MTF is one of the
blueprint contents. The core system software is configured according to
the MTF, and multiple partitions are scheduled in a cyclic manner. (2) The
resources (such as timers, stacks, and memory) and blueprints required
for partition scheduling are in a privileged state, and the user partition
cannot overwrite the partition scheduling resources.
• Blueprint only registers the interfaces on the modules and then mounts
them on the app as a whole. The purpose of blueprint itself is to organize
the parallel coexistence of multiple modules and avoid registering modules
directly on the app. In fact, it is more convenient for development and code
maintenance, because ultimately all interfaces on views are still directly
mounted on the app, which corresponds to the entire application; there is
no obvious difference [11].
• At the same time, we cannot use multiple objects to manage and register,
because this will cause each object to have its own configuration, which is
not easy to manage.
• With blueprint, the application will be managed in the flask layer, share
the configuration, and change the application object on demand through
7
Satellite Systems - Design, Modeling, Simulation and Analysis
Based on the menu design idea, the avionics system can realize the sharing
of hardware modules and resources, task migration, and system reconstruction,
enhance the tolerance and processing ability of the avionics system for faults,
significantly improve the development efficiency and productization degree of
the integrated electronic system, and make the integrated electronic system highly
reliable. The intelligent spacecraft provides the necessary technical support, mainly
including the contents that will be presented in the subsequent sections [12].
ii. Telemetry function: The SMU collects its own telemetry and receive the indi-
rect telemetry parameters collected by the bus terminal equipment through
the 1553B bus. The SMU complete the framing processing according to the
CCSDS standard and transmit it to the TTU through the serial port.
v. Important data saving function: The SMU can periodically save and maintain
important data of satellite sub-systems. When the internal configuration of
the relevant sub-system changes or the corresponding module restarts, the
SMU can send the important data stored internally to the corresponding
module. Besides, when the SMU is reset or switched off, the SMU can restore
the current working state through the important data.
vi. Fault Detection Isolation and Recovery (FDIR) function: The SMU provides
the operating platform for satellite FDIR [13]. When the SMU is healthy, the
SMU monitor key information such as the entire satellite’s energy and thermal
control. The SMU detects various failure conditions in real time and performs
troubleshooting through the direct remote command interface or through the
1553B bus. The FDIR of the control sub-system is completed by the attitude
control computer. Note that the most advanced 1553B bus can handle data rate
8
Design of Intelligent and Open Avionics System Onboard
DOI: http://dx.doi.org/10.5772/intechopen.93141
up to 10 Mbps. For data rate higher than 10 Mbps, industry trend is moving
toward SpaceWire data bus that can handle data rate up to 400 Mbps.
The SMU adopts the dual-machine cold backup working mode. In the case of the
autonomous switching enabled, when the on-duty machine fails, the failure-tolerant
module automatically completes the switch of SMU and power off the faulty
machine. In the case of the autonomous switching not enabled, when the on-duty
machine fails, the power failure and machine switch is conducted by ground station.
The Integrated Service Unit (ISU) adopts modular hardware design [14, 15].
Each functional module is connected to the “bus interface management module”
through an internal bus, which uses widely acceptable data bus such as 1553B bus.
The ISU is mainly composed of a bus interface management module, a matrix
instruction and matrix telemetry module, and an analog quantity acquisition and a
discrete instruction output module.
The menu module is mainly composed of a bus interface management module,
a matrix acquisition and command module, and an information acquisition and
command module [16]. The module menu is shown in Table 1 below.
The functions of each module are also modularized. The capability of each
module is shown in Table 2.
After the module design is completed, the number of modules is determined
according to the task requirements. Then, complete the assembly according to the
standard interface. The number of modules menu is shown in Table 3.
The data bus network is the information transmission hub of the avionics
system. Through the data bus network, distributed data acquisition, and instruction
output, centralized operation and control are implemented, thereby improving the
efficiency of system processing. The avionics system first-level bus is 1553B bus [17].
Note that for data rate higher than 10 Mbps, SpaceWire data bus is recommended.
Data exchange between SMU and ISU and other equipment realized through
1553B bus. The master–slave response mode of the 1553B bus is used to transmit
platform command data and telemetry acquisition data. In the data exchange
process of the first-level bus, the SMU always acts as the controller of the 1553B bus
Table 1.
Function module list.
9
Satellite Systems - Design, Modeling, Simulation and Analysis
Table 2.
Modularized functions.
Demand Quantification AA AB AC CA
Pyrotechnics management 80 2
High current instruction 30 1
Low current instruction 100 1
Analog acquisition 220 2 2
Matrix acquisition 550 2
Total (take the maximum) 2 2 2 2
Table 3.
Menu-style design.
and initiates the communication. ISU and other equipment, as the RT end, receive
instructions and send collected telemetry data to the SMU.
In order to achieve autonomous and healthy operation of the satellite, the intel-
ligent satellite system uses the FDIR software to monitor the status of the satellite in
real time and diagnose and predict its working status and performance trends [18].
When a failure occurs, the FDIR software can locate the failure in time and deter-
mine which components are not working normally or the performance is degraded.
Design goals:
10
Design of Intelligent and Open Avionics System Onboard
DOI: http://dx.doi.org/10.5772/intechopen.93141
ii. When a failure occurs, try to extend the mission time of the satellite and
reduce the loss of mission interruption.
iii. The life of the satellite should be guaranteed: optimize fuel consumption and
minimize system configuration and component losses.
The above three principles apply to the launch phase, the transfer orbit phase,
and the on-orbit phase.
FDIR is an important component of the onboard software, which can perform
on-orbit processing of failures, thereby reducing the impact of failures. However,
not all on-orbit failures can be detected and processed. FDIR design should follow
the following principles:
i. FDIR processing follows the single failure principle, that is, only one failure
is processed at a time.
ii. Failures are divided into 0 to 4 levels according to their impact on satellites.
iii. The higher the failure level, the higher the processing priority. During a
failure processing, if a higher-level failure occurs, the higher-level failure is
processed first.
iv. The failures of the same level are processed in the order of occurrence.
According to the impact of the failure on the satellite operation, the failure is
categorized as follows:
11
Satellite Systems - Design, Modeling, Simulation and Analysis
avionics FDIR, and power supply and distribution FDIR as illustrated in Figure 3.
The larger the number, the higher the fault level, and vice versa.
• Level 0 failure: a level 0 failure refers to a failure that occurs inside an equip-
ment and can be recovered autonomously by the hot backup method inside the
equipment without affecting other components of the system.
• Level 2 failure: a level 2 failure refers to the functional level abnormality of the
satellite sub-system. Under such failures, the system performance cannot meet
the design requirements. For level 2 failures, the recovery strategies need to be
implemented and related components need to be enabled or restarted. Level 2
failures can cause system performance degradation or temporary interruption
of system tasks. Its failure detection, isolation, and recovery are performed by
application software.
• Level 3 failure: a level 3 failure refers to the failure of the CPU hardware, which
is detected by the hardware. After the failure occurs, it is switched to the
backup CPU according to the failure handling strategy.
• Level 4 failure: a level 4 failure refers to the failure of the satellite to maintain
the pointing to the ground in the on-orbit phase and requires sun capture
processing.
Figure 3.
Schematic diagram of satellites in orbit during their lifetime.
12
Design of Intelligent and Open Avionics System Onboard
DOI: http://dx.doi.org/10.5772/intechopen.93141
The software autonomously isolates the failure and rebuilds the system at the
appropriate time according to the following FDIR scheme:
• High-level failure detection such as level 4 failure has priority over low-level
failures such as level 3. When two or more failures are detected at the same
time, the recovery sequence of high-level failures is performed preferentially.
Once a high-level failure occurs, all detection of the same-level and low-level
failures are suspended before the recovery sequence is completed.
• Only one failure recovery sequence is performed on the satellite at the same
time, that is, all FDIR failure recovery is shielded during the execution of any
failure recovery strategy. The sufficiency and necessity of the failure recovery
sequence should be effectively verified by ground testing to minimize the
interpretation during the sequence execution.
• After the execution of the failure recovery sequence is completed, the failure
detection of the flight can be continued, but the failure recovery enable status
should be set to disable. At the same time, the detection of other similar and
low-level failures should be enabled. After confirming the working status
of the products on the ground, reset the backup status and enable the FDIR
recovery function.
• FDIR only detects the status of the on-duty module. When the FDIR enable flag
is “disabled,” the status of the module is not detected. When the FDIR enable flag
is “enabled,” the status of the on-duty module is detected. If the failure detec-
tion condition is met on the on-duty module and the FDIR recovery enable flag
is “disabled,” the health status of the module is set to unhealthy. If the failure
detection condition is met on the on-duty module and the FDIR recovery enable
flag is “enabled,” it is determined whether the status of the on-duty module is
the same as the backup module. If they are the same, do not perform the recov-
ery operation and set the on-duty module as unhealthy. If they are different,
perform the recovery operation and set the backup module to the on-duty status.
• For dual-machine hot standby equipment or modules, only the health status of
non-duty module is detected, and no recovery is performed.
• Use its own fault-tolerant RAM and lower computer to save important data in
time for state recovery after failure.
The FDIR requirements for each phase of the satellite are as follows:
i. Launch phase: allows failure detection and recovery of level 0 and level 1
failures.
ii. Transfer orbit phase: allows failure detection and recovery of level 0 to level 3
failures.
13
Satellite Systems - Design, Modeling, Simulation and Analysis
iii. On-orbit phase: allows failure detection and recovery of level 0 to level 4
failures.
ii. Fault detection: Determine whether a failure occurs based on the recognition
characteristics.
14
Design of Intelligent and Open Avionics System Onboard
DOI: http://dx.doi.org/10.5772/intechopen.93141
6. Conclusion
The intelligent avionics system design is the key technology for future advanced
satellites. The system design has adopted modular and open system architecture
approach using an efficient computing hardware system to maintain multiple
central processing units and memory executing instructions accurately and syn-
chronously with high cost-performance and cost-efficiency ratio. This approach
improves the failure tolerance of the next-generation avionics systems. The func-
tion modification capabilities and function migration between modules realize the
transition to software-defined satellite. Furthermore, this approach also reduces the
differences in hardware products and improves the failure tolerance of the satellite’s
internal network, which meets the ever-increasing networking requirements of the
satellite.
Acknowledgements
This work was supported in part by the National Natural Science Foundation of
China (No. 61972398).
15
Satellite Systems - Design, Modeling, Simulation and Analysis
Author details
© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License (http://creativecommons.org/licenses/
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
16
Design of Intelligent and Open Avionics System Onboard
DOI: http://dx.doi.org/10.5772/intechopen.93141
References
[8] Marie J. Spacebus 4000 avionics: Key [16] Panpan Z, Tingyuan G, Jianjun G,
features and first flight return. In: 24th et al. Plug-and-play on-board computer
AIAA International Communications system design based on BM3803
Satellite System Conference (ICSSC). processor. Journal of Aerospace
2006 Engineering. 2013;22(06):92-96
17
Satellite Systems - Design, Modeling, Simulation and Analysis
18
Chapter
Abstract
During the entire launch sequences from liftoff to final orbital insertion of a
space vehicle (SV), adequate link requirements are to be maintained for telemetry,
tracking, and command (TT&C) for uplink and downlink services, from launch
vehicle (LV) and SV to ground stations (GS). A successful space vehicle launch
required adequate link coverage with good radio frequency (RF) performance. The
chapter is an extension of the IEEE/Aerospace Conference 2019 paper entitled
Dynamic Link Analysis and Application for a MEO Space Vehicle published by the
authors. The emphasis in this chapter is on the addition of the three distinctively
different tracking waveforms and their associated links, used from liftoff to final
orbital insertion. This chapter will describe the three required dynamic link analyses
(DLA) to cover (a) the LV link from liftoff to its end of line of sight (LOS), (b) the
LV link from LOS to Tracking and Data Relay Satellite System (TDRSS) at beyond
line of sight (BLOS), and (c) the final tracking link using Space-to-Ground Link
Subsystem (SGLS) or non-SGLS (NSGLS) link for the earliest or best separation
time of the SV from the LV. The chapter discusses the concept of the dynamic link
analysis, SV antenna switching schedule, recommended SV separation time, as well
as the performance for different launch scenarios within the 24-h launch window.
Topics include antenna patterns, launch trajectories, elevation angle and clock and
cone angle geometry, and dynamic link budget.
1. Introduction
Recently, there was an interest to extend the present dynamic link analyses
(DLA) beyond the early launch period to cover the period after the space vehicle
(SV) separation from launch vehicle (LV), which includes both booster and second
stage engine. The dynamic link from liftoff to final orbital insertion considers both
geometric (visibility coverage) and radiometric (link margins for all downlink and
uplink services) adequacy in the three launch stages. The purpose of the dynamic
link study for the launch is to provide the earliest and accurate time for final SV
separation and orbital insertion as compared to previous method which only relied
on visibility tracking coverage to the end of line of sight (LOS).
The present DLA typically covers only two stages of LV tracking, including
(a) liftoff to the end of LOS link and (b) the end of LOS to a period before SV
1
Satellite Systems - Design, Modeling, Simulation and Analysis
payload (PL) separation from LV, using LV, to Tracking and Data Relay Satellite
System (TDRSS) [1] satellite link, which is also called beyond line of sight (BLOS)
link. A third SV tracking, after SV payload separation from LV, is a tracking link
between SV and a ground station (GS). This third SV tracking is now added in this
chapter.
The tracking link used from liftoff to the end of LOS uses a UHF noncoherent
FSK signal for command and a digital FM or BPSK for tracking telemetry link as
described in detail in [2]. From the end of LOS to BLOS, the tracking telemetry link
is usually a BPSK or QPSK signal, using a NASA Tracking Data Relay Satellite
System to relay tracking data from the LV to White Sands or Goddard ground
station (WSGT/GRGT) and finally routing it to other user ground stations. After SV
payload separation and orbital insertion, the SV tracking link to an Air Force
Satellite Control Network (AFSCN) ground station [3] will use Space-to-Ground
Link Subsystem (SGLS) or a non-SGLS (NSGLS) waveform described in [3, 4] and
in Section 3 for tracking signals along the trajectory. In the following pages,
supporting link analyses for the two LV and one SV tracking stages will be
presented.
Figure 2 defines the antenna coordinate system used in this chapter. The
azimuth (AZ or Ф) or clock angle is used in the antenna cut configuration. The
elevation angle (EL or ϴ) also called as cone angle is also used in the antenna
gain data file. The antenna gain is a changing variable as a function of
mission elapsed time (MET). The antenna gains are used in the following
dynamic links.
Figure 1.
Space vehicle and ground station vector definition.
2
Dynamic Link from Liftoff to Final Orbital Insertion for a MEO Space Vehicle
DOI: http://dx.doi.org/10.5772/intechopen.92462
Figure 2.
Antenna coordinate definition.
This section provides a summary of the dynamic link model of interest [6]. More
detailed derivation of other variables, especially UVZBD, UVYBD, and UVXBD, can
be found in [5]. For station elevation angle, either LV or SV elevation or cone angle
(EL or theta or ϴ), and LV or SV clock or azimuth angle (AZ or Phi or Ф), we have1
1800
UVZBD ð SRÞ
1
Cone Angle ¼ cos in Deg (1)
π jSRj
1800
0 1 GS SR
EL ¼ 90 cos in Deg (2)
π jGSj ∗ jSRj
1800
1 UVXBD ð SRÞ UVYBD ð SRÞ
φ¼ tan 2 , in Deg (3)
π jSRj jSRj
As mentioned in the introduction, there are three separate tracking signals along
the flight trajectory that we need to analyze ensuring that they have adequate link
margins of three dB or more. Most of the present DLA covers only two stages of
launch coverage and neglecting the third stage coverage. The requirement for third
stage tracking is explained below:
1
The authors would like to thank Dr. James Yoh for his derivation of these formulas.
3
Satellite Systems - Design, Modeling, Simulation and Analysis
1. From liftoff to the end of LOS, the waveform for this link is generally a digital
FM or BPSK for telemetry downlink as defined in the Range Commander
Council (RCC) 119-88 [2]. The liftoff to LOS 5 link margin plot is shown in the
left half of Figure 3 for five different ground stations.
Figure 3.
TLM dynamic link margin from TEL4 to TDRSS versus mission elapsed time.
Table 1.
LV to TDRSS range TLM link (based on NASA source).
4
Dynamic Link from Liftoff to Final Orbital Insertion for a MEO Space Vehicle
DOI: http://dx.doi.org/10.5772/intechopen.92462
2. From the end of LOS to NASA TDRSS at geosynchronous orbit, the telemetry
link is a BPSK or QPSK (telemetry + data) which is sent from LV to TDRSS to
be relayed to White Sands or Goddard ground station (WSGT/GRGT), as
shown in the second right half of the Figure 3 and in more link details in
Table 1.
3. For the third link after SV payload separation, when the satellite or SV starts its
transfer orbit, the tracking link from the SV payload (or bus) to an AFSCN
ground station will be in SGLS, Unified S-Band (USB), or a NSGLS waveform
as described in more detail in [3, 4]. For a more secure tracking, SV normally
will be using SGLS link for tracking as described in [4], with a MEO satellite in
Table 2 as an example. A commercial and less secure SV launch may use USB
Table 2.
Link budgets for SGLS TT&C uplink and downlink services [4].
5
Satellite Systems - Design, Modeling, Simulation and Analysis
Table 3.
Typical C/No for uplink and downlink budgets.
or NSGLS for SV tracking instead of using SGLS waveform. Table 3 shows link
budget for uplink and downlink C/No example for tracking links 1 and 2.
Table 2 shows SGLS telemetry, tracking, and command (TT&C) link budget
for tracking link 3 for a MEO satellite. If the SV is using a USB or a NSGLS [3],
the tracking waveform can be an AQPSK signal with telemetry on the I
channel and ranging on the Q channel.
6
Dynamic Link from Liftoff to Final Orbital Insertion for a MEO Space Vehicle
DOI: http://dx.doi.org/10.5772/intechopen.92462
For the uplink, where the transmit antenna is located on the ground, the antenna
can be easily directed to an LV or SV. This transmit antenna is likely to be a
directional high gain dish antenna for connectivity with an LV or SV located possi-
bly far away. In Table 3, line 9, a typical ground station has a large parabolic
antenna dish with a gain Gt = 43.94 dBi using the following formula.
π ∗ fC ∗ D
Gt ¼ 10 ∗ log 10 η ∗ in dBi (5)
c
The signal path traverses through the transmission medium, in between transmit
and receive systems. When the distance between transmit and receive systems
increases, the signal beam has an angular spread which decreases the signal power
collected by a receiving antenna. We know that the portion of the transmission
medium near the ground station depends on the Earth’s atmosphere which attenuates
the signal to different degrees, dependents on the frequency, the altitude of the GS,
and the angle of the signal path through the atmospheric (GS elevation angle).
Beyond the Earth’s atmosphere, the signal path traverses through the space with little
atmospheric attenuation, only with free space loss to account for. Therefore, there are
essentially two losses through the transmission medium, namely, the space loss to
account for the spreading of the signal beam and the additional atmospheric loss [7].
" 2 #
c
LS ¼ 10 ∗ log 10 in dB (6)
4π ∗ f C ∗ SR
The transmit signal, after accounting for the space and atmospheric losses, and
its signal path terminates at the antenna of either the ground station, the SV, or the
LV receiving system. Before considering the characteristics of the receive antenna
7
Satellite Systems - Design, Modeling, Simulation and Analysis
and the receiver, a good indication of the signal strength is given by the received
isotropic power. RIP is simply the transmitter EIRP after subtracting off the losses
of the transmission medium, i.e.
The last portion of Table 3 addresses the receiving system and assesses how well
it performs. This section involves with the calculation of the signal strength and the
noise strength, resulting in the ratio of signal power over noise power density
(C/N0). In general, we first address the signal power and then the noise power
density. Line 21 provides for the receive antenna size consistent with the receive
antenna gain to be calculated later. The next parameter in Table 3 is the polariza-
tion loss, which accounts for the mismatching between the polarization axial ratios
of the received signal and the receiving system. The axial ratio is the ratio of the
major axis of an ellipse to its minor axis. For circularly polarized signal, the ratio
should be 0 dB. Any deviation from 0 dB results in a polarization loss. Line 23 shows
the values of receive antenna gain. For downlink in which the receive antenna is a
dish antenna located at GS, Gr is calculated using the standard dish antenna equa-
tion (similar to Eq. (6) for Gt).
π ∗ fC ∗ D
Gr ¼ 10 ∗ log 10 η ∗ in dBi (8)
c
where NF = low noise amplifier noise factor, in dB and LL = line loss, in dB.
The noise density (N0) is given by
8
Dynamic Link from Liftoff to Final Orbital Insertion for a MEO Space Vehicle
DOI: http://dx.doi.org/10.5772/intechopen.92462
þ 10 ∗ log ðT S Þ in dBW=Hz
N 0 ¼ kdB (11)
For many SVs, we are interested in their uplink and downlink services. Table 2
shows an example, taken from an IEEE paper [4]. This is the standard link budget,
where the ground station is an AFSCN [3] remote tracking station (RTS) using
SGLS waveform [3, 4]. The waveform is described in an AFSCN interface control
document (ICD) [3] and is implemented in DLA, although other waveforms can be
readily incorporated. The uplink has two services of interest—carrier and command
—while the downlink has three services of interest: carrier, ranging, and telemetry.
In general service margins are calculated for these five services. For the SGLS
waveform, command is coupled with ranging and modulated on the uplink carrier;
therefore command is also turned around at the SV along with ranging. This SGLS
turnaround process explains the reason that Table 2 shows a power allocation for
command in the downlink and no calculation for its margin. As a result, downlink
power allocated to command is essentially wasted while robbing power from other
downlink services. The requirements and service margins for command and telem-
etry are expressed in Eb/N0, since it is the bit error rate (BER) that counts for both
cases. The carrier and ranging are expressed in C/No given by a specific station. For
ranging, it is the autocorrelation value between the decoded ranging code and the
transmitted ranging code that needs to be maximized in order to successfully
perform accurate ranging.
Table 2 represents uplink and downlink budgets for SGLS TT&C. Let us address
the important aspects of the calculation of uplink and downlink services in the next
few subsections. The role of modulation indices is to divide up the power for
allocation to services. The modulation index is expressed in radians so that it can go
right in as an argument in a sinusoidal or Bessel expression. If the modulation
indices of all services are zero radians, no power is allocated to the services, and the
carrier retains all the link power calculated in Section 4. If the modulation indices of
services are not zero, portions of the power are taken from the carrier and allocated
to the services. The remaining power stays with the carrier as the “residual carrier
power.”
After SV separation, we are dealing with the SV uplink and downlink using
SGLS or NASA Unified S-Band waveforms as described in [3, 4]. For telemetry
service, the requirement is SNR = Pservice/NoB = Eb Rb/NoB = Eb/No in dB. For
carrier and ranging, the requirements are stated in terms of C/No as mentioned
before. For acquisition, the uplink carrier loop bandwidth could be as high as +/
100 KHz, while its tracking bandwidth could be as small as a few Hz. For the station
the carrier tracking loop bandwidth is about 20–50 Hz, as in Table 2 in line 36. For
ranging, the bandwidth of 10 Hz represents ranging tracking loop bandwidth
(Table 2, line 56), which corresponds to the sampling rate of the autocorrelation
9
Satellite Systems - Design, Modeling, Simulation and Analysis
1 AM-3FSK/PRN RNG/PM (SGLS unfiltered J 2o ðβ1 Þ: cos 2 ðβ2 Þ J 2o ðβ1 Þ: sin 2 ðβ2 Þ 2J 21 ðβ1 Þ: cos 2 ðβ2 Þ
uplink) [3, 4]
2 BPSK/PRN RNG/PM (USB filtered case) J02(β1)cos2(β2) J02(β1) sin2(β2) 2J12(β1) cos2(β2)
(Eqs. (1-22) and (1-23)) [9]
4 PRN RNG/PSK TLM/PM (SGLS filtered J02(β01 )J02(β02 ) 2J02(β01 ) J12(β02 ) 2J02(β01 )J02(β02 )
case) J02(β3 ) J02(β3) J12(β3 )
(Eqs. (2-18) and (2-20)) [4]
5 Tone RNG/PSK TLM/PM (USB filtered J02(β3 )J02(β02 ) 2J02(β3 )J12(β02 ) 2J02(β02 )J12(β3)
case)
(Eqs. (2-18) and (2-21)) [9]
Table 4.
Uplink and downlink service modulation losses for SGLS and NASA USB.
value between the detected ranging code and the transmitted ranging code. For
command and telemetry, the requirements are expressed in terms of Eb/N0. The
command and telemetry data bit rates of 1000 bps each are representing the lower
end of their SGLS choices. As shown in Table 2, the results from the SNR calcula-
tion are the values of C/N and Eb/N0 for various received uplink services (lines 37
and 48 for command C/N or Eb/No) and for various downlink services (line 55 for
ranging service Prng/No and line 70 for telemetry service Eb/No). The uplink
service modulation losses for SGLS and NASA USB with subcarrier (S/C) [3, 4, 9]
are shown in Table 4. The downlink service modulation losses for SGLS and NASA
USB with subcarrier (S/C) [3, 4, 9] are shown in Table 4. Also note that β1, β2, and
β3 represent the modulation indices for command (CMD), ranging (RNG), and
telemetry (TLM), respectively, per [3, 4, 9]. These uplink and downlink modula-
tion losses are in lines 34, 44, and 54 in Table 2.
For NSGLS waveforms such as the direct mod BPSK, QPSK, and AQPSK, the
service mod losses are negligible. Finally, the calculated service SNR in Table 2 is
compared with the required SNR to obtain the link margin for each service. The
required SNR values capture all the performance requirements for the services,
such as ranging accuracy, tracking loop loss likelihood, bit error rate, and others.
Before SV separation from the LV, we are also interested in the dynamic link
from a ground station to the LV, from liftoff to the SV after separation, along the
entire LV flight path using the tracking stations in the line of sight (TEL4, JDMTA,
ANT, DGS, TDRSS). The waveforms for this LV tracking are described in Range
Commander Council (RCC) handbook [2]. One must ensure that the downlink
telemetry link from the launch vehicle to these ground stations and TDRSS relay
satellite are positive as can be seen in Figure 3. The basic LV range modulations are
digital FM, BPSK, QPSK, AQPSK, etc. as discussed in RCC [2]. In Figure 4, dynamic
LV slant range “received TLM Eb/No” and “TLM link margin” for a specific mission
10
Dynamic Link from Liftoff to Final Orbital Insertion for a MEO Space Vehicle
DOI: http://dx.doi.org/10.5772/intechopen.92462
Figure 4.
LV slant range received TLM Eb/No and TLM link margin.
7. Conclusion
This chapter discusses the required three DLAs and related tracking waveforms
to cover the three launch stages, namely, (a) the launch vehicle tracking link from
liftoff to its end LOS using the digital FM or BPSK signal, (b) the launch vehicle
tracking link from LOS to TDRSS at BLOS using NASA USB signal, and (c) the final
tracking link from SV to an AFSCN ground station using AFSCN SGLS, AFSCN
NSGLS, or NASA USB waveforms. In the third tracking link case, BPSK, QPSK, or
AQPSK waveforms were used, in which for QPSK and AQPSK, the telemetry data is
put on the I channel and the ranging signal is on the Q channel.
The chapter shows that good telemetry link margins from LV to tracking sta-
tions such as TEL4, JDMTA, and ANT or to a NASA TDRSS relay satellite can be
achieved using digital FM, BPSK, QPSK, or AQPSK signals, after SV separation. The
chapter also shows that good tracking link margins can be achieved from SV to
AFSCN ground stations, including IOS or DGS as the first contact station.
11
Satellite Systems - Design, Modeling, Simulation and Analysis
Author details
© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License (http://creativecommons.org/licenses/
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
12
Dynamic Link from Liftoff to Final Orbital Insertion for a MEO Space Vehicle
DOI: http://dx.doi.org/10.5772/intechopen.92462
References
13
Chapter
Abstract
The chapter describes effective algorithms that are often used in processing data
measurements in Global Navigation Satellite Systems (GNSSs). Existing effective
algorithm was developed for detection and elimination of outliers from GNSS data
measurements. It is based on searching for a so-called optimal solution for which
standard deviation and maximum absolute deviation of the measured data from
mean values do not exceed specified threshold values, and the number of the
detected outliers is minimal. A modification of this algorithm with complexity of
N log 2 N is discussed. Generalization of the existing algorithm to the case when data
series included some unknown trend will be presented. The processing trend is
assumed to be described by an unknown function of time. The generalized algo-
rithm includes the outlier detection algorithm and trend searching algorithm that
has been tested using simulated data. A new algorithm will be presented for cycle
slip repair using Melbourne-Wübbena linear combination formed from GNSS data
measurements on two carrier frequencies. Test results for repair data in the case of
multiple (cascade) cycle slips in actual observation data will also be presented in
this chapter.
1. Introduction
1
Satellite Systems - Design, Modeling, Simulation and Analysis
yj ¼ f j þ z þ ξj ; j ¼ 1 … N, (1)
2
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
yj ¼ z þ ξj ; j ¼ 1 … N, (2)
are satisfied:
rffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
2
X
σYL ¼ ðL 1Þ 1 j ∈ j ,::,j yj z ≤ σmax , (3)
f 1 Lg
n o
yj z ≤ 3 σmax ; yj ∈ yj , … , yj ,
(4)
1 L
L ≥ MINOBS, (5)
where σYL and σmax are the standard deviation (SD) and their associated
specified threshold values; MINOBS is a parameter limiting from below the length
of the required set of measured values (e.g., 10), and we will assume hereafter that
MINOBS < N.
3
Satellite Systems - Design, Modeling, Simulation and Analysis
The values yj that are not included in the set YL are treated as coarse measure-
ments and removed from further processing. Typically, expressions in Eqs. (3)–(5)
are the only conditions considered in processing programs when screening out
rough measurements. Below we will formulate the problem of searching for a
solution, complementing the conditions expressed in Eqs. (3)–(5) with two extreme
conditions [7].
1. First, we will require that the length of the set sought be the maximum, that is,
the number of measurements deemed to be coarse is the minimum:
L ! max : (6)
Note that for the predetermined values yj , the problem solution satisfying
conditions in Eqs. (3)–(5) may not exist (e.g., when yj includes a trend, in
particular when yj is an arithmetic progression with a step greater than σmax ). In
the case when the solution does exist, we will denote the value L at which the
maximum of Eq. (6) is reached as Lmax . Note that the condition expressed in
Eq. (6) does not ensure the uniqueness of the set because several sets of length
Lmax can be found that satisfy the conditions in Eqs. (3)–(5).
2. From all possible sets that satisfy conditions expressed in Eqs. (3)–(5) and (6),
we will select the one for which the variable σYL receives the smallest value:
σYL !
n o min : (7)
YL ¼ yj , … , yj ; L¼Lmax
1 L
satisfying conditions in Eqs. (3)–(7), we refer to as the optimal solution of the problem
expressed in Eqs. (3)–(7). The corresponding SD value is denoted by σopt .
Thus, the problem consists in the creation of a search algorithm for the optimal
solution of the problem shown in Eqs. (3)–(7).
In a practical situation, the precise value z, given conditions in Eqs. (3) and (4),
is not known. We will estimate the values using the following formula:
X
z¼L 1 yj : (9)
j ∈ fj1 , … , jL g
Note that the value z depends on the required solution, which will complicate its
search.
Usually, iterative methods are used to find a solution to problem expressed in
Eqs. (3)–(5). For example, the algorithm implemented in the observation data
smoothing program (see [3]) is designed to find a set Y satisfying the conditions
given by Eqs. (3)–(5). The proposed step-by-step algorithm is based on iterations
(the index number of iteration is designated by the upper index in parentheses):
Step 1: Initialization: Y ð0Þ ¼ y1 , … , yN ; Levelð0Þ ¼ 1020 ; k = 0.
Step 2: Checking the length of the set Y ðkÞ ; if it is less than MINOBS, then the
process comes to an end and a solution is not found.
4
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
Step 3: Calculation of the values zðkÞ and σðkÞ on the available set Y ðkÞ using
formulas (3) and (9).
Step 4: Checking the fulfillment of the inequality σðkÞ ≤ σmax. If it is satisfied, the
set Y ðkÞ is accepted as the solution, and the search process comes to an end.
Otherwise, there is a transition to Step 5.
Step 5: Definition of Levelðkþ1Þ for outlier detection:
Levelðkþ1Þ ¼ 3 σðkÞ ;
Step 6: Definition of a new set Y ðkþ1Þ to include those and only those y j for which
z ≤ Levelðkþ1Þ :
ðkÞ
y j
Let us formulate a statement that is the key to the creation of an effective search
algorithm for the optimal solution(Eqs. (3)–(7)).
Assertion 1. Let the set Y opt ¼ y j , …,y j be optimal for a specified sequence
1 L max
n o
of values y j and
ymin ¼ min y j , … , y jL , ymax ¼ max y j , … , y jL ,
1 max 1 max
then the interval ymin , ymax does not contain values y j that are not in the set Y opt .
Proof. In fact, let us assume the opposite: Let y j ∉ Y opt , ymin < y j < ymax and yk
and yl be values from the set Yopt for which yk = ymin and yl = ymax. One of these cases
is possible:
2
2
a. z < y j and, therefore 0 < y j z < yl z => y j z < yl z ,
2
2
b. z ≥ y j and, therefore 0 ≤ z yj < z yk => y j z < yk z .
5
Satellite Systems - Design, Modeling, Simulation and Analysis
In the first case, Case (a), we replace the value yl in the set Yopt with y j . In the
second case, Case (b), we replace yk with y j . In any of the cases, we will have
another set of the same length Lmax for which conditions expressed in Eqs. (3) and
(4) are satisfied, but the SD, as follows from Eq. (3), is less than σopt . Consequently,
the set Yopt is not optimal since requirement expressed in Eq. (7) is not satisfied.
The contradiction that is reached proves Assertion
1.
Further, note that if Yopt ¼ yj , … , yj is the optimal solution of the prob-
1 L max
Eqs. (3)–(7). Thus, the optimal solution does not depend on the arrangement of
given numbers yj . This allows us to arrange the numbers yj in the order most
suitable for searching the optimal solution.
n o Taking advantage of this important
circumstance, we arrange the values yj in the ascending order and we will look
for the optimal solutionnin o
the ordered sequence. For brevity, the ordered sequence
will also be denoted by yj . Thus, y1 ≤ y2 ≤ … ≤ yN . Moreover, for simplicity of
logic, we will assume that all yj are different, that is,
Note that due to the formulated Assertion 1, if Yopt is the optimal set and y j and
1
yj are its smallest and greatest values, respectively, then all values of yj from the
L max
interval y j , y j belong to Yopt . Consequently, considering Eq. (10), we have
1 L max
yj ¼ yj þL 1 and
L max 1 max
n o
Yopt ¼ yj , yj þ1 , … , yj þLmax 1
1 1 1
Thus, the optimal solution should be sought in the ascending sequence yj among
all possible sets yk , … , ykþL 1 of length L with k and L satisfying the following
conditions:
MINOBS ≤ L ≤ N, (11)
1≤k≤N L þ 1: (12)
Hence, instead of searching for all possible sets of various length, numbering 2N ,
for the solution of the problem described by Eqs. (3)–(7), it is sufficient to vary just
the two parameters, k and L, associated with the conditions expressed in Eqs. (11)
and (12). The number of pairs of integer numbers k and L subject to condition in
Eqs. (11) and (12) is equal to:
1 kþL
X1
zðk; LÞ ¼ y, (13)
L j¼k j
6
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
kþL
X1
1 2
σ2 ðk; LÞ ¼ yj zðk; LÞ : (14)
L 1 j¼k
We rewrite conditions given by Eqs. (3) and (4) in the new designations:
Note that the two last inequalities directly follow from Eq. (4), monotony of yj,
and obvious inequalities: yk ≤ zðk; LÞ ≤ ykþL 1 .
Remark. In the conditions expressed in Eqs. (15) and (16), L means the length of the
set under checking, and k is the index of the smallest number included in the set. Although
“k” and “L” are also encountered as indexes in the sets we use below for monotonically
increasing sequences, we hope nevertheless that this will not lead to confusion.
The following recursive relationships are available, making it possible to find
zðk; LÞ and σ2 ðk; LÞ through the calculated values zðk; L þ 1Þ и σ2 ðk; L þ 1Þ for seven
arithmetic operations:
The algorithm described below is based on the search for all possible pairs (k, L),
where L denotes the length of the set to be checked and k is the index of the smallest
of the values included in the set. At that, k and L must satisfy conditions Eqs. (11)
and (12). The set yk , … , ykþL 1 corresponding to each such pair must be checked
for fulfillment conditions (15) and (16).
We organize this search according to the algorithm described below, at each step
of which we check the validation of Eqs. (15) and (16) for all possible sets of a
certain length. We start the examine process with Step 1, where we check the set of
maximum length N. Further, with each next step, we will reduce the length of the
sets to be checked by 1.
Step 1: We consider the set of length N. There is only one such set: y1 , … , yN .
We check it for fulfillment conditions expressed in Eqs. (15) and (16). If they are
satisfied, this set is selected as a solution, and further search stops. Otherwise,
transit to Step 2.
Step 2: We consider the sets of length N 1. There are two sets of length N 1.
7
Satellite Systems - Design, Modeling, Simulation and Analysis
y1 , … , yN 1 and y2 , … , yN :
We test each of these sets for compliance with the conditions specified by
Eqs. (15) and (16). If none of them satisfies conditions (15) and (16), then we
transit to the next step. Otherwise, the following options are available:
• Option 1: if only one set from them is found that satisfies conditions (15) and
(16), then it will also be the solution of the stated problem; the search process
stops here, where Lmax = N 1.
• Option 2: if both sets simultaneously satisfy conditions (15) and (16), we will
select the set corresponding to the smallest of two values σ2 ð1, N 1Þ or
σ2 ð2, N 1Þ, and the search process stops here, where Lmax = N 1.
Step N L + 1: We consider the sets of length L. If L < MINOBS, then the search
process stops, and a solution is not found. For L ≥ MINOBS, we examine N L + 1
sets of length L:
We check each of these sets for fulfillment of conditions (15) and (16). If any of
them does not satisfy these conditions, then we transit to the next step where we
consider the sets of length (L 1). Otherwise, two options are possible:
• Option 1: if only one set from (23) is found that satisfies the conditions of (15)
and (16), then it will also be the solution of the stated problem with Lmax = L,
and the search process stops here.
In order to calculate the values zðk, LÞ and σ2 ðk, LÞ, we use the recursive formu-
las (17)–(22) in accordance with a search scheme shown in Figure 1 where only the
zðk, LÞ is involved.
In accordance with the proposed arrangement, we calculate the values zð1; NÞ
and σ2 ð1; NÞ in the first step of the algorithm using formulas expressed in Eqs. (13)
and (14), and 4 N arithmetic operations are required for this arrangement. In order
to find zðk, LÞ and σ2 ðk, LÞ on all subsequent steps, we apply the recursive formulas
(17)–(22), making it possible to calculate the values of the specified pairs of
Figure 1.
Scheme of calculations when finding the optimal solution.
8
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
variables, each for the seven arithmetic operations, based on the results of the
calculations of the preceding step. So, for example, on the second step, proceeding
from the known values zð1; NÞ and σ2 ð1; NÞ, we find the values zð1; N 1Þ and
σ2 ð1; N 1Þ (vertical arrow on the diagram) by using formulas in Eqs. (17)–(19)
and the values σ2 ð2; NÞ and σ2 ð2; NÞ (horizontal arrow on the diagram) by using
Eqs. (20)–(22). In the general case, the transition to the following step in the
direction of the vertical arrows (see diagram) is carried out according to formulas
(17)–(19), and in the direction of horizontal arrows, according to formulas
(20)–(22). The number of arithmetic operations required to find the solution should
not exceed:
We test the algorithms discussed above on the real data obtained by the ONSA
station that is a part of the IGS network [2]. These data are included in the distri-
bution kit of the installation software package [3] and available for usage. We
consider measurement data received from global positioning system (GPS) satellite
with system number PRN = 12 for 2010, day 207 to check the efficiency of the
proposed algorithm described in Section 3. Figure 2 plots the values of the
Melbourne-Wübbena combination over a time interval of 89.5 min (N = 180). The
index numbers j of time epochs counting from the beginning of a 24-h period with a
30-second interval are plotted on the horizontal axis. The values y j of the combina-
tion are plotted on the vertical axis and are expressed in cycles with wavelength
λ5 ≈ 0:86 [3]. Figure 3 shows the values of deviations from the mean of the data
cleared of outliers using the algorithms described in Sections 2 and 3. In both cases,
σmax = 0.6 and MINOBS = 10. The values y j z are plotted on the vertical axes in
cycles with wavelength λ5 and the index numbers j of epochs on the horizontal axis.
Figure 2.
Melbourne-Wübbena combination for ONSA station (GPS satellite, PRN = 12 for 2010, day 207).
9
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 3.
(a) Deviations of values of the Melbourne-Wübbena combination from the mean value after data cleaning
from outliers using the algorithm described in Section 2. (b) Deviations of values of the Melbourne-
Wübbena combination from the mean value after data cleaning from outliers using the developed algorithm
(see Section 3).
Epochs in which the data were rejected are designated by white circles. In the first
case (see Figure 3a), 47 data of the measurements were rejected, which are 26.1% of
the total amount of data. In the second case (see Figure 3b), 14 of these measure-
ments were rejected (7.8%), which are almost 18% less than in the previous
calculation.
We also provide similar results for data obtained by TLSE station, which is also
included in the IGS network. We consider measurement data from GLONASS,
Russia satellite with system number PRN = 1 for 2010, day 207. Figure 4 shows the
values of the Melbourne-Wübbena combination over a time interval of 65.5 min
(N = 132). Figure 5 plots the values of deviations from the mean value of the data
cleared of outliers using the algorithms described in Sections 2 and 3, respectively.
Parameters σmax and MINOBS are the same as in the previous calculation example.
In the first case (see Figure 5a), 41 data of the measurements were discarded, which
are 31% of the total amount data. In the second case (see Figure 5b), 8 of these
measurements were rejected (6%), which are 25% less than in the previous
calculation.
10
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
Figure 4.
Melbourne-Wübbena combination for TLSE station (GLONASS satellite, PRN = 1 for 2010, day 207).
Figure 5.
(a) Deviations of values of the Melbourne-Wübbena combination from the mean value after data cleaning
from outliers using the algorithm described in Section 2. (b) Deviations of values of the Melbourne-Wübbena
combination from the mean value after data cleaning from outliers using the developed algorithm (see
Section 3).
11
Satellite Systems - Design, Modeling, Simulation and Analysis
Note that the number of arithmetic operations required to find the optimal
solution according to the algorithm described in Section 3 depends on the Lmax ,
which is the length of the solution. As can be seen from Eq. (24), the smaller the
length of the found solution (i.e., the larger the number of detected outliers), the
more arithmetic operations are required to find it. This number of arithmetic
operations is estimated to be of order N þ ðN Lmax Þ2 . Note that the expression in
the parentheses herein is equal to the number of outliers detected. Thus, if the
number of outliers in the original data series is comparable to N, it will take N2
arithmetic operations to find the optimal solution. In particular, to make certain
that there is no solution (e.g., in the case where the data contain a trend), N2
arithmetic operations will also be required. In Section 6, we will modify the existing
algorithm and describe fast outlier search algorithm that requires Nlog2 N
arithmetic operations.
The necessary preparations are given in this section. Note that in this and the
n oN
next sections we are dealing with the sequence yj arranged in the ascending
j¼1
order.
n oN
Assertion 2. Let yj be monotonically increasing sequence. The following
j¼1
inequality is true:
Proof. From the monotonicity of the sequence yj and the definition zðk, L þ 1Þ
[see Eq. (13)],
Suppose, for example, the case (a) holds. Let us show that in this case
Truly, inequalities:
yk ≤ yj ≤ ykþL ; j ¼ … , k þ L,
12
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
These inequalities, in turn, imply Eq. (28). Next let us prove (27). This inequal-
ity is expanded as follows:
kþL kþL
X1
1X 2 1 2
y zðk; L þ 1Þ ≥ yj zðk; LÞ :
L j¼k j L 1 j¼k
Substituting here of the expression (17) in place of zðk; LÞ and writing for
brevity, z instead of zðk; L þ 1Þ, we get inequality:
kþL kþL
X1
X 2 z ykþL 2
ðL 1Þ yj z ≥L yj z : (31)
j¼k j¼k
L
kþL
X ykþL z2 ðL þ 1Þ2 2
RHSð31Þ ¼ L yj zþ ykþL z
j¼k
L L
kþL
X 2 ðL þ 1Þ
2 ðL þ 1Þ2 2
¼L yj z þ ykþL z ykþL z :
j¼k
L L
PkþL
Here we take into account the equality: j¼k yj z ¼ 0. Next, we have:
kþL
X 2
2
RHSð31Þ ¼ L yj z ðL þ 1Þ ykþL z :
j¼k
kþL
X 2
2
yj z ≤ ðL þ 1Þ ykþL z ,
j¼k
that is true due to Eq. (28). Thus, Eq. (27) is proved for case (a). Analogously,
case (b) is considered.
We introduce the notation:
That is, the sequence σ2min ðLÞ monotonically decreases when L decreases from N to
MINOBS.
Proof. Assertion 2 and definition of σ2min ðLÞ expressed in Eq. (32) imply:
13
Satellite Systems - Design, Modeling, Simulation and Analysis
holds for some L0 , then for existence of the solution YL ¼ yk , ykþ1 , … , ykþL 1 for
the problem expressed in Eqs. (3)–(7), it is necessary that the length L of the set Y L
satisfies the condition L < L0 .
Proof. Let us assume that L ≥ L0. Assertion 3 on account of monotony of σ2min ðÞ,
expressed in Eq. (33) implies the following inequalities σ2min ðLÞ ≥ σ2min ðL0 Þ > σ2max for
all L≥L0 . From this, it follows that for any set YL ¼ yk , ykþ1 , … , ykþL 1 we will
have σ2 ðk, LÞ ≥ σ2min ðLÞ > σ2max . Thus, any of sets YL of length L ≥L0 does not satisfy
the condition in Eq. (15) and, therefore, cannot be a solution of the problem,
expressed in Eqs. (3)–(7).
In particular, we have come to the next important result. If, for example, the
inequalities σ2min ðMINOBSÞ > σ2max are fulfilled, then the solution for the problem
described in Eqs. (3)–(7) does not exist.
In the above-described procedure for solving problem (3)–(7), it takes N2
arithmetic operations to make certain that the solution not exists. Taking into
account Assertion 3 and Corollary 1, the search procedure may begin by checking
the conditions.
The above proposed search procedure consists in the calculating values of zðk, LÞ
and σ2 ðk, LÞ using recurrent formulas (17)–(22) and checking at every k and L the
fulfillment of the inequalities (11) and (12). The algorithm complexity is estimated
by value of N þ N2Outlier , where NOutlier is the number of outliers found. If it is
known a priori that there are few outliers in the measurement data, then the search
algorithm for the optimal solution that described in Section 3 can be applied. If the
measurement data contain an N-comparable number of outliers, the complexity of
such an algorithm will be estimated by about N2. It is namely for such type of data we
below describe a modified outlier search algorithm with complexity of about Nlog2 N.
First of all, note one property that is the key to the construction of a fast outlier
search algorithm. Note that if the inequality (15) holds for some set of length L + 1,
then there exists a set of length L for which the inequality (15) is valid too. Truly, let
assume for some k the inequality σ2 ðk, L þ 1Þ ≤ σ2max holds. This inequality and
Eq. (25) imply
14
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
This means that at least one of these sets yk , … , ykþL 1 and ykþ1 , … , ykþL
with length of L satisfies conditions expressed in Eq. (15).
However, this property is not true when checking the conditions expressed in
Eq. (16). In other words, if these conditions are fulfilled for any set of length L + 1,
it might happen that none of the sets of length L may satisfy them. This fact is a
significant obstacle to increasing the rate of outlier detection that is necessary
when processing a large amount of data with a large number of rough measure-
ments. To overcome this obstacle, we will make the condition expressed in
Eq. (16) weaker.
First of all, note that if for some set yk , … , ykþL 1 , the both conditions
expressed in Eq. (16) are fulfilled, then the following inequality will hold:
15
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 6.
Three possible cases for the proposed search. In case (a) we go to the right-hand side range (range for length of
sets) to find a solution, in case (c) we go to the left-hand side range to find a solution, in case (b) we look for a
ðk 1Þ
solution with length L = N Mid and the search ends.
The need to process GNSS measurements including a trend on which noise and
outliers are superimposed arises at different processing stages of the application
process. As already stated above, satellite clock corrections contain a linear trend. In
some cases, it may not be known, and then, one has to search for it, for example, by
the least square method. The presence of outliers in the measurement data is a
significant obstacle to accurate determination of drift and offset parameters of
satellite clocks. Other examples are linear combinations of code and phase data on
two carriers [3]. To obtain high accuracy results, it is necessary to detect outliers
against an unknown trend and remove them from further processing. This is the
subject of this section.
16
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
yj ¼ f j þ z þ ξj ; j ¼ 1 … N:
The procedure described above for finding the optimal solution in an ordered
series of numbers may not produce an adequate result if applied to data containing
an unknown trend. For example, there may be no solution, and all data will be
defined as outliers. In order to detect outliers in a series of numbers with a trend
using the algorithm described above, it is necessary to find a suitable approximation
of an unknown function f j and subtract it from the data set. Searching for this
approximation is usually done by selecting functions from some functional class.
The choice of the functional class depends on the task. In some cases, these may be
power polynomial, in other cases, trigonometric polynomial, etc. The presence of
outliers in the data measurements makes it much more difficult to find such an
approximation. In this section, we will describe the general approach to solving the
problem and look for suitable approximations of trend f j in the class of power
polynomial with real coefficients:
!
Pn,j a ¼ an ðj=NÞn þ an 1 ðj=NÞn 1
þ … þ a0 , (37)
!
where n is the polynomial degree, and a ¼ fa0 , … , an g is vector of coefficients.
Thus, the problem consists in the creation of an algorithm for searching the
trend in the class of power polynomial and detecting outliers in specified data series
yj represented in Eq. (1).
Before we turn to the trend search algorithm construction, we will define the so-
called minimizing set of given length L, which plays an essential role in the trend
search. In addition, in Section 7.3, we will describe a search algorithm for such set
based on the nrecurrent formulas
o (17)–(19) and (20)–(22).
Let YL ¼ yj , … , yj be an arbitrary set of length L, composed of the values of
1 L
n oN
a numerical series yj , the monotony is not supposed. The mean and the SD
j¼1
values for it are denoted by zY L and σYL . These values are calculated by standard
formulas:
X
zYL ¼ L 1 yj , (38)
j ∈ fj1 , … , jL g
X 2
σ2YL ¼ ðL 1Þ 1 yj zYL : (39)
j ∈ fj1 , ::, jL g
n oN
Definition 2. Given L for a specified sequence of values yj the set of values:
j¼1
n o
YL, min ¼ yj , … , yj ,
1 L
17
Satellite Systems - Design, Modeling, Simulation and Analysis
at which the minimum value of σ2YL defined in Eq. (39) is reached will be called the
minimizing set of length L. The corresponding mean and SD values are denoted by zYL, min
and σYL, min .
According to this definition, we have
X
zYL, min ¼ L 1 yj , (40)
j ∈ fj1 , … , jL g
X 2 n o
2 1
σYL, min ¼ ðL 1Þ yj zYL, min ¼ min σ2YL : (41)
YL
j ∈ fj1 , … , jL g
then the interval ymin , ymax does not contain values yj that are not in the set YL, min .
Proof. Let us assume the opposite: Let yj ∉ YL, min , ymin < yj < ymax . Let
yk , … , ykþL 1 is a permutation of the numbers yj , … , yj in the ascending order;
1 L
then yk ¼ ymin and ykþL 1 ¼ ymax . One of these cases is possible:
In the first case, Case (a), we replace the value yk in the set yk , … , ykþL 1 with
yj . In the second case, Case (b), we replace the value ykþL 1 with yj . We want to
show that doing such replacement, the value σ2YL, min expressed in Eqs. (40) and (41)
will decrease. This will mean that YL, min is not a minimizing set.
Suppose Case (a). For brevity, we will write below z instead of zYL, min and σ2
~2 , the similar values obtained after replacement yk
instead of σ2YL, min . Denote ~z and σ
with yj . We have:
1
z¼L yk þ … þ ykþL 1 , (44)
18
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
1 2 2
σ 2 ¼ ðL 1Þ yk z þ … þ ykþL 1 z , (45)
and
1
~z ¼ L ykþ1 þ … þ ykþL 1 þ yj , (46)
2
2 1 2 2
~ ¼ ðL
σ 1Þ ykþ1 ~z þ … þ ykþL 1 ~z þ yj ~z : (47)
~2 < σ2 :
σ (48)
After simplification with taking into account Eq. (45), we get from here:
2 2
2 2 1 2 1
~ ¼ σ þ ðL
σ 1Þ yj z yk z L yj yk :
From (42), it follows (recall that we write z instead of zYL, min ) that the expression
in the square brackets is strictly less than zero. Thus, inequality (48) is proven.
From this follows that the set YL, min is not minimizing one because the condition
(41) is not met. Thus, we have arrived at a contradiction that proves the validity of
the formulated Assertion 6. Case (b) is considered similarly.
19
Satellite Systems - Design, Modeling, Simulation and Analysis
we will assume that all yj are different, that is, Eq. (10) holds. Rewrite it for
convenience:
y1 < y2 < … < yN
Note that due to Assertion 6, if YL, min is the minimizing set and yj and yj are its
1 L
smallest and greatest values, respectively, then all values of yj from the interval
yj , yj belong to YL, min . Consequently, considering Eq. (10), we have yj ¼
1 L L
yj þL 1 and therefore
1
n o
YL, min ¼ yj , yj þ1 , … , yj þL 1 :
1 1 1
and choose from them the one that has minimal SD value. In the minimizing set
searching procedure, we use the appropriate designations zðk, LÞ and σ2 ðk, LÞ
expressed by Eqs. (13) and (14) to denote mean and square SD in the case of sets
composed of the ascending ordered values.
The calculation of zðk, LÞ and σ2 ðk, LÞ when searching YL, min can be carried out
in the manner similar to that of the optimal solution according to the schematic, in
which only the σ2 ðk, LÞ is shown:
Figure 7.
Scheme of calculations when finding the minimizing set of length L.
At first, we carry out the transitions in the direction of the vertical arrows
(see diagram in Figure 7) and calculate sequential values
σ2 ð1, NÞ, σ2 ð1, N 1Þ, … , σ2 ð1, LÞ using formula (17)–(19). Finally, we go
in the direction of the horizontal arrows and calculate values of
σ2 ð1, LÞ, σ2 ð2, LÞ, … , σ2 ðN L þ 1, LÞ using formula (20)–(22) and choose from
them the minimal one.
used for fitting. If the data do not contain coarse measurements, all N values of the
original data can be treated as reference ones and used for fitting. If the data contain
20
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
rough measurements, the participation all or some of those values in the fit will
result in a fatal distortion of the trend function and, as a result, an incorrect
determination of the outliers. Thus, when determining the trend, it is necessary to
exclude rough measurements from the number of reference values used for fitting
and by which the trend approximations are built. We have the vicious circle. To
determine outliers, it is necessary to find a proper trend approximation, and to find
an exact trend approximation, we need to find all outliers and remove them from
the values used for fitting. Another vicious circle is obtained when trying to choose
the number of values for fitting. If we take it too small to minimize the possibility
for outliers to be included into the subset of values for fitting, the data fit may not
be accurate enough for the rest values of the data outside of the subset. If, on the
contrary, we take rather large number of points for fitting, the outliers may be
among them, and we will also get a bad approximation of the trend.
We describe herein a strategy for finding an unknown trend and detecting outliers.
This strategy assumes that the number of outliers in the series presented by Eq. (1)
does not exceed a certain value Nmaxout known a priori. Thus, we can suppose that the
number of “right” values suitable for fitting in the series y j is not less than N Nmaxout.
Below L is supposed to be fixed and associated with the number of the reference
values of the series (1) used for fitting, and L ≤ N Nmaxout .
Let us consider the following algorithm. It contains internal iterations, which we
will denote with the upper index “s” in parentheses.
Step 0: n = 0. We set some L, satisfying the condition L ≤ N Nmaxout .
ð 0Þ
Step 1: n++; s = 0; flagj ¼ 1.
!ðsÞ
Step 2: s++. We fit polynomial to the data set and find fitting coefficients a ¼
n o
ðsÞ ð sÞ
a0 , … , an :
!ðsÞ 1Þ !
a ¼ arg min Φðs a , (50)
a0 , … , an 1 , an
XN 2
1Þ ! ! ðs 1Þ
Φðs a ≜ yj Pn,j a flagj : (51)
j¼1
ðsÞ
^j obtained after subtraction from yj the polynomial
Step 3: Consider the values y
values with coefficients found in Step 2:
ð sÞ
ðsÞ !
y
^j ¼ yj Pn,j a : (52)
n oN
ðsÞ
Using the algorithm described in Section 7.3, we find from the numbers y ^j
j¼1
ðsÞ ðsÞ ðsÞ
defined in Eq. (52) a minimizing set YL, min ¼ y ^ ðsÞ , … , y
^ ðsÞ of length L. For this
j1 jL
set, we calculate the corresponding values of the mean zYðsÞ and SD σYðsÞ
L, min L, min
21
Satellite Systems - Design, Modeling, Simulation and Analysis
achieved. Note that we will consider the issue of convergence further in Section 7.5.
After reaching convergence of values σYðsÞ , we transit to Step 4 for outlier detection.
L, min
ðsÞ
ðsÞ !
Step 4: Searching the optimal solution for data set y
^j ¼ yj Pn,j a , j = 1, … , N.
ðsÞ
We find the optimal solution for values y
^j using the algorithm of Section 3 and
determine the number Nout of outliers detected. If it turns out that Nout ≤Nmaxout ,
ðsÞ
!
then solution is considered found; searching process stops: f j ¼ Pn,j a . Other-
wise, if Nout >Nmaxout , then we transit to Step 1.
We do this until we find a solution or until n reaches some preset value Nmax
(e.g., 10). In this case, probably, we may need to select a different functional class
to search for a trend.
Example. Let us consider the data simulated in accordance with formula:
pffiffiffiffiffiffiffiffiffiffiffiffiffiffi
y j = 10 j þ 10 þ 2random j , where j ¼ 1 … N, random j —pseudorandom numbers,
evenly distributed in the segment [0,1]. Let us set N = 150 and model 10 outliers at
points j ¼ 6:::10 and j ¼ 140:::144 (see Figure 8). Further, let us set σmax = 0.6,
MINOBS = 10 and search for outliers as described above with L = 130. If n = 1, we
get N out = 88; if n = 2, we get N out = 33; if n = 3, we get N out = 17; if n = 4, we get N out
= 10. Thus, a suitable trend approximation turned out to be a fourth degree poly-
nomial. Figure 9 shows values for fourth degree polynomial fitted to the data set on
Figure 8. pffiffiffiffiffiffiffiffiffiffiffiffiffiffi
Data simulated using: y j = 10 j þ 10 þ 2random j and fourth degree polynomial after the first iteration.
Figure 9.
Simulated data and fourth degree polynomial after the eighth iteration.
22
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
Figure 10.
!
The differences ^y j ¼ y j P4,j a , approximation of unknown trend with fourth degree polynomial.
This section explains the convergence of the iterations described in the trend
search algorithm (see previous section).
Assertion 7. The SD sequence σYðsÞ of values calculated in the trend search algo-
L, min
Transform the expression on the right side of Eq. (56). Substitution here instead
ðsÞ
of y
^j expression in Eq. (52) gives:
X 2
1 !ðsÞ
σ2YðsÞ ¼ ðL 1Þ yj Pn,j a zYðsÞ :
L, min L, min
ðs Þ ðsÞ
j∈f j1 , … , jL g
23
Satellite Systems - Design, Modeling, Simulation and Analysis
ðsÞ
Using the flagj defined in Eq. (53), rewrite the last equality as follows:
XN ðsÞ 2
2 1 ! ðsÞ
σYðsÞ ¼ ðL 1Þ yj Pn,j a zYðsÞ flagj ¼
L, min L, min
j¼1
XN !2 !
1 ðsÞ 1 ðsÞ
¼ ðL 1Þ yj Pn,j b flagj ¼ ðL 1Þ Φ b : (57)
j¼1
! ðsÞ ðsÞ ðsÞ
Here we introduce the designation b ¼ þ zYðsÞ a0
and take into , a1 , … , an
L, min
!
account the definition of Pn,j a expressed by Eq. (37) and the definition of
!
functional ΦðsÞ a given in Eq. (51). From Eq. (57), we get:
!
σ2YðsÞ ¼ ðL 1Þ 1 ΦðsÞ b (58)
L, min
!ðsþ1Þ
We transit to Step 2; s is incremented by 1. Wefind
a vector a of polynomial
ðsÞ !
coefficients, which minimizes the functional Φ a :
!ðsþ1Þ ðsÞ !
a ¼ arg min Φ a :
!
a
Thus,
ðsþ1Þ n o
! ðsÞ !
ΦðsÞ a ¼ min
!
Φ a (59)
a
ðsÞ !
The extremum condition (one of n + 1) of functional Φ a is as follows:
∂ ðsÞ !
Φ a ¼ 0:
∂a0 ! !ðsþ1Þ
a¼a
N
X ðsþ1Þ
! ðsÞ
yj Pn,j a flagj ¼ 0:
j¼1
ðsþ1Þ
Taking into account the designation for y
^j given by Eq. (52), we can write
this equality in the form.
N
X X
ðsþ1Þ ðsÞ ðsþ1Þ
y
^j flagj ¼ 0 or y
^j ¼ 0: (62)
j¼1 j∈f
ðsÞ
j1 ,
ðs Þ
… , jL g
24
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
ðsþ1Þ ðsþ1Þ ðsþ1Þ
Consider the set YL ¼ y
^ ðsÞ , … , y
^ ðsÞ . The mean and SD values for it are
j1 jL
calculated using formulas (38) and (39). Taking into account Eq. (62), we get:
X ðsþ1Þ
1
zYðsþ1Þ ¼ L y
^j ¼ 0,
L
ðsÞ ðs Þ
j ∈ fj 1 , … , j L g
and
X 2 X 2
1 ðsþ1Þ 1 ðsþ1Þ
σ2Yðsþ1Þ ¼ ðL 1Þ y
^j zYðsþ1Þ ¼ ðL 1Þ y^j ¼
L L
ðsÞ ðs Þ ðs Þ ðs Þ
j ∈ fj1 , … , jL g j ∈ fj1 , … , jL g
N
X 2
1 !ðsþ1Þ ðsÞ 1 ðsÞ !ðsþ1Þ
¼ ðL 1Þ yj Pn,j a flagj ¼ ðL 1Þ Φ a :
j¼1
Thus,
ðsþ1Þ
!
σ2Yðsþ1Þ ¼ ðL 1Þ 1 ΦðsÞ a : (63)
L
n oN
ðsþ1Þ
We transit to Step 3. In this step, for sequence y ^j , we find a minimizing
j¼1
ðsþ1Þ ðsþ1Þ ðsþ1Þ
set of length L: YL, min ¼ y ^ ðsþ1Þ . The SD square σ2 ðsþ1Þ for this set does
^ ðsþ1Þ , … , y
j1 jL YL, min
ðsþ1Þ
not exceed the same magnitude for any other set, particularly for YL ¼
ðsþ1Þ ðsþ1Þ
y
^ ðs Þ , … , y
^ ðs Þ :
j1 jL
Finally, from Eqs. (61) and (63) and the last inequality (64), we get
ðsþ1Þ
!
σ2YðsÞ ≥ ðL 1Þ 1 ΦðsÞ a ¼ σ2Yðsþ1Þ ≥ σ2Yðsþ1Þ
L, min L L, min
In this and the following section, we describe cycle slip repair algorithm for
observers represented in the form of Melbourne-Wübbena combination, which is
often used in modern GNSS measurement data processing programs. Loss by the
receiver of the carrier phase capture results in jumps in the code and phase mea-
surement data. In the absence of jumps, as we already discussed in Section 2, the
values of the combination consist of measurement noise superimposed on an
unknown constant value dependent on a specified satellite-receiver pair.
In case of temporary loss of carrier phase capture by the receiver, jumps occur in
a series of values representing the Melbourne-Wübbena combination. The proce-
dure of detecting jumps and eliminating them from the values of the combination,
25
Satellite Systems - Design, Modeling, Simulation and Analysis
called cycle slip repair, is one of the most important steps of preprocessing GNSS
data. The main difficulty in detecting jumps is that neither the exact size of jumps
nor their epochs are known. A number of algorithms, descriptions of which can be
found, for example, in Refs. [3, 6, 8] are proposed for the detection of jumps.
Although differing in detail, they are based on a common idea, that is, comparison
of the SDs of the time series of measurement data obtained from one of the bounds
of the time interval to an arbitrary moment, an epoch. If the differences of the SDs
corresponding to two adjacent epochs exceed a predetermined threshold value, then
a jump is declared in one of these two epochs. A drawback of similar algorithms is
the frequent false detection of jumps during epochs containing rough measure-
ments (outliers) since the values of outliers can exceed the sizes of a jump itself. On
the other hand, an attempt to increase the threshold value leads to the opposite
effect, an inability to recognize jumps that are small in magnitude.
Below, we propose a robust cycle slip repair algorithm that allows, more reliably
than similar known algorithms, to detect jumps and determine their sizes. The
proposed algorithm is based on search for so-called clusters consisting of epochs, in
which the values of the combination are grouped about corresponding predefined
values. Besides, this algorithm implements the above-described method of cleaning
data from outliers based on the search for the optimal solution. This method,
combined with Springer’s algorithm used in Ref. [3], allows for the reliable deter-
mination of multiple (cascade) jumps in the Melbourne-Wübbena combination.
The Melbourne-Wübbena combination LM W can be presented in the form of
the total of three terms [6]:
LM W ¼ λ5 n5 þ B þ ν, (65)
where λ5 is the formal wavelength (for all GPS satellites λ5 = 86.16 cm and for
GLONASS satellites λ5 = 84.0 ÷ 84.36 cm); n5 is an unknown integer, the so-called
wide lane ambiguity [3]; B is the residual systematic error caused by instrumental
delays relative to the specific receiver-satellite pair and, as assumed, not time-
dependent; and ν is a random component. Both parts of the equality in (65) are
expressed in meters. At the same time, the LM W combination is often expressed in
cycles of wavelength λ5 . Let us divide both parts of Eq. (65) by λ5 . We introduce the
designations y ¼ LM W =λ5 , β ¼ B=λ5 , and ξ ¼ ν=λ5 and derive for each of the epochs
associated with indexes j ¼ 1, … , N:
y j ¼ n5j þ β þ ξ j : (66)
Figure 11.
Searching for the maximum of the density of values of the array Y j .
logic, we suppose that all Y j are different. Therefore, we have: Y 1 < Y 2 < … < Y N .
Denote the corresponding permutation of the values of function flagj as FLAGj .
Step 1: On this step, we are looking for the maximum density of array values y j .
We consider the segments Y j , Y j þ Δ of length Δ and look for the one that
contains the largest number of Y i values. Here, only those indexes i and j from the
range of 1 ≤ i, j ≤ N are considered for which FLAGi ¼ FLAG j = 0. The purpose of
the ordering of the original array consists in improving the effectiveness of the
maximum density search procedure. It can be shown that the number of compari-
sons required when searching for the above segment does not exceed 2 N in the case
of the ordered array. For an unordered array, the number of comparisons is evalu-
ated as N 2 .
Step 2: Let Y J max , Y J max þ Δ be the segment found in the previous step, containing
N max values of Y j (see Figure 11). We calculate the mean m of these numbers Y j .
Xmax 1
J max þN
1
m¼ N max Y j, (67)
j¼J max
Y j ∈ Y Jmax , Y J max þ Δ : (68)
Step 3: Cluster searching. The values y j that are clustered about the mean m
found in the previous step can pertain to the epochs scattered along the time axis
in chaotic order. The task of this step is to define an accumulation of such points,
a so-called cluster, which we define as follows:
Definition 3. We designate as an ðm, ΔÞ cluster the set of points (epochs) of the time
axis, the indexes j of which belong to the segment [k,l] (1 ≤ k < l ≤ N) and for which the
following requirements are satisfied:
a. all points of the segment [k, l] are marked by the same value: flagk = … = flagl .
b. at the points j= k, l of the left and right boundaries of the segment, this inequality
is satisfied:
y j m ≤ Δ:
(69)
c. the amount of points at which (69) is satisfied is no less than the present value of
MINOBS;
d. the number of consecutive points in which (69) is not satisfied does not exceed the
predefined value MAXGAP (e.g., 5);
27
Satellite Systems - Design, Modeling, Simulation and Analysis
and since the arrays y j and Y j differ only by permutation, it follows from
Eq. (68) that the inequalities
Y J max ≤ y j ≤ Y J max þ Δ: (71)
are fulfilled for N max index j. Inequality (69) is satisfied also for all these indexes,
as follows from Eqs. (70) and (71).
If cluster was found, then we mark all points of it as 1, and then, we repeat the
cluster search procedure. If even just one cluster has been found at this step, we
transfer to Step 1. If not even one cluster has been found, then the search for
clusters is complete and we transfer to Step 4.
Step 4: If even just one cluster has been found, we transfer to Step 5, and
otherwise: (a) all points of the segment ½1; N are marked as outliers and removed
from further processing and (b) the operation of the algorithm terminates.
Step 5: Search for individual jumps in clusters. Let us assume that n ≥ 1 clusters
have been found:
In each of the clusters that are found, outliers and 1-size jumps are possible.
This follows immediately from the inequalities in (69) and the preestablished value
Δ ¼ 1.2.
Figure 12.
Epochs with indexes k ≤ j ≤ l form ðm, ΔÞ cluster.
28
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
ð0Þ
where y j ¼ y j .
ðnÞ
Substep 5.4: We rename: y j ¼ y j .
Step 6: Marking of points outside clusters. All points outside of the found clus-
ters (if any) are marked as outliers.
Step 7: Ordering of clusters. We renumber the clusters so that they are placed
left to right on the time axis. For the ordered set of clusters [kp, lp], p = 1, … , n, the
conditions 1 ≤ k1 < l1 < k2 < l2 < … < kn < ln ≤ N are satisfied.
Step 8: Data screening within clusters and improving the mean values of y j .
Substep 8.1: In accordance with the algorithm proposed in Section 3, we perform
screening from outliers in each of the n clusters.
Substep 8.2: For each of the clusters cleaned of outliers, we determine the
modified mean values mp∗ .
Step 9: Jumps between clusters. It follows from the description presented above
that the remaining jumps in the data y j are on boundaries between clusters. If only
one cluster is found, the algorithm is terminated: no additional jumps are detected,
and the data are cleaned of outliers. If more than one cluster is found (n > 1), then
the epochs j1, … , jn–1 of the jumps will be the coordinates of the left-hand bound-
aries of clusters, beginning from the second: jl = k2, … , jn–1 = kn. The values of jumps
Δn5, jp are found in this case as rounded to the nearest integer of the difference of
mean values of adjacent clusters:
∗
Δn5, jp ¼ NINT mpþ1 mp∗ , p ¼ 1, … , n–1: (73)
Step 10: Repair data. We delete the jumps between clusters using formula anal-
ogous to (72):
8
< yðjp 1Þ
ðpÞ
; j < jp
yj ¼ 1Þ
; p ¼ 1, … , n 1, (74)
: yðp Δn5,J p ; jp ≤ j < N
j
ð0Þ
where y j ¼ y j .
ðn 1Þ
We rename: y j ¼ y j . The algorithm is terminated.
We present here the results of testing the proposed algorithm using real data
obtained by the JOZ2 station, which is part of the IGS network [2]. These data are
29
Satellite Systems - Design, Modeling, Simulation and Analysis
included in the distribution set of the installation software package [3]. Testing was
carried out for data obtained from GPS satellite with number PRN = 13 for 2010,
day 207. Figure 13 shows the Melbourne-Wübbena combination values in the
Figure 13.
Values of the Melbourne-Wübbena combination for the JOZ2 station (PRN = 13 for 2010, day of year = 207).
Figure 14.
(a) Deviations of the values of the Melbourne-Wübbena combination from the mean after detection and
elimination of jumps from the data using the algorithm from Ref. [3]. (b) Deviations of the values of the
Melbourne-Wübbena combination from the mean after detection and elimination of jumps from the data using
the proposed algorithm in Section 9.
30
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
107 min time interval (N = 215). The number j of time epochs counted from the
beginning of 24-h day with interval of 30 seconds is plotted along the horizontal
axis. The combination values y j expressed in cycles with wavelength λ5 are plotted
along the vertical axis.
Figure 14a and b presents the values of the deviations from the mean value z of
data after cycle slip repair procedure, by using the algorithm applied in Ref. [3] (see
Figure 14a) and the proposed algorithm (Figure 14b). The values y j –z in cycles of
λ5 are plotted along the vertical axes, and the number j of epochs is plotted along the
horizontal axis. Epochs in which the measurement data were rejected are marked by
light circles. In the first case (see Figure 14a), 111 of the conducted measurements
or 51% of the total number of data was discarded. In the second case (see
Figure 14b), 29 of these measurements were rejected (13%), which are almost 38%
less than in the previous computation. The epochs of detected jumps are marked by
daggers.
11. Conclusion
This chapter presents several effective and stable algorithms for processing data
received from GNSS receivers. These data form the basis of almost all engineering
applications in the field of computational geo-dynamics and navigation and cadas-
tral survey and in numerous fundamental research works as well. The accuracy of
the results obtained is significantly influenced by the quality of the data used in the
calculations. In particular, the presence of rough measurements (outliers) in the
observation data can significantly reduce the accuracy of the calculations carried
out. One of the tasks at the preliminary stage of data processing is reliable detection
and removal of rough measurements from the series of measured data with mini-
mum amount of rejected data. The so-called optimal solution, introduced in the
chapter, made it possible to detect and eliminate outliers from observed data min-
imizing the number of rejected measurements. In addition, it is assumed that the
data may contain a trend as an unknown function of time. The strategy for deter-
mining of the trend is depending on the physical process in question under an
assumption that the trend is a continuous function of time. The efficiency of the
search is definitely influenced by the choice of the function class from which the
trend is searched. In this chapter, we considered the class of power polynomials, but
in some cases, this choice may not lead to the expected result. It may require, for
example, a class of trigonometric functions to find a suitable trend. The automatic
search for the best functional class, together with the strategy of effectively finding
an unknown trend, against the background of random noise and outliers, is a
complex task for future research.
31
Satellite Systems - Design, Modeling, Simulation and Analysis
Author details
Igor V. Bezmenov
“VNIIFTRI”, Mendeleevo, Moscow region, Russian Federation
© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License (http://creativecommons.org/licenses/
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
32
Effective Algorithms for Detection Outliers and Cycle Slip Repair in GNSS Data Measurements
DOI: http://dx.doi.org/10.5772/intechopen.92658
References
33
Chapter
Analysis of Spatiotemporal
Variability of Surface Temperature
of Okhotsk Sea and Adjacent
Waters Using Satellite Data
Dmitry Lozhkin
Abstract
The chapter is divided into 5 parts. The first part describes what the satellite data
is and describes the levels of its processing. In the second part, attention is paid to
sea surface temperature anomalies, the conditions for the appearance of the most
significant anomalies and their influence on the behavior and survival of aquatic
organisms are described. The third part is devoted to calculating linear trends in
ocean surface temperature from a 20-year series of satellite data. It is shown that
the heat content of the surface layer of Okhotsk Sea decreases, most significantly in
its northern and western parts. This trend is especially pronounced in the spring,
which may be due to a decrease in ice cover and a more significant cooling of the
waters due to winter convection. In the fourth part, periodic fluctuations in the
temperature of the surface of the ocean are considered. It is demonstrated how,
using the calculated trend and several basic harmonics, one can try to predict the
temperature next year. And the last part concludes the chapter.
Keywords: satellite data, sea surface temperature, Tsushima current, Okhotsk Sea,
regional effects of global warming, SST anomalies, trends, harmonic, cyclicity
1. Introduction
The study of the thermal regime of various water areas is one of the most
important oceanological problems, since the spatio-temporal variability of water
temperature reflects complex processes of formation, transformation and dynam-
ics of water masses. In addition, temperature is one of the key parameters that
determine the conditions for the existence and development of most species of
aquatic organisms; therefore, the study of this problem is also of key importance for
hydrobiology. The zoning of water areas by the nature of temperature conditions,
as well as their forecasting, taking into account the peculiarities of seasonal and
interannual variability, is an important scientific task, which also has a pronounced
applied aspect associated with the fact that accumulations of some species of com-
mercial fish are confined to the zones of separation of water masses with different
characteristics.
Direct measurements are a traditional source of water temperature data.
Their accuracy depends only on the accuracy of the device, and the discreteness
1
Satellite Systems - Design, Modeling, Simulation and Analysis
depends on the specific task (in some cases it can reach a fraction of a second).
They can also be used to construct vertical profiles with high vertical resolu-
tion. However, recently the number of expeditions has been steadily declining,
and the data obtained through direct measurements are point and irregular and
are not suitable for studying large-scale phenomena or for obtaining data from
hard-to-reach places.
Satellite data, on the other hand, are regular and allow covering the entire
water area of the studied basin. Therefore, they are very good sources for studying
seasonal and interannual variations in the temperature of the surface layer (the
thickness of which ranges from 1 to 10 meters). The disadvantage of these data is
the strong dependence of accuracy on cloudiness (especially in the infrared and vis-
ible ranges). Ideally, they should be regularly compared with direct measurements
in order to identify errors in data interpretation.
For scientific purposes, data are usually used from satellites located either in a
geosynchronous (in a particular case, geostationary orbit) or in a heliosynchronous
orbit. The advantage of the geostationary orbit is obvious: you can get a picture of
the same area with discreteness of up to half an hour, and after pointing to the satel-
lite, a constant correction of the antenna position is not required. It is very useful
for telecommunication systems as well as for obtaining meteorological data. The
disadvantage is the lack of coverage at polar latitudes.
This disadvantage is easily eliminated by satellites in a polar sun-synchronous
orbit. A sun-synchronous orbit is important to science because it keeps the angle
of incidence of sunlight on the Earth’s surface more or less constant, although the
angle will change with the change of seasons. Without a sun-synchronous orbit, you
would have to account for changes in shadow and lighting angles, making it difficult
to track changes over time. It would simply be impossible to collect the information
needed to study climate change.
In 1997, a TeraScan satellite receiving station was installed at SakhNIRO, with
the help of which data are received from the NOAA, Metop, Aqua and Terra satel-
lites in a polar sun-synchronous orbit (Figure 1 [1]). At the moment, a 21-year
series of satellite data on the temperature of the surface of the Sea of Okhotsk
and adjacent waters has been accumulated, which makes it possible to analyze the
interannual variations of this parameter.
Figure 1.
Example of a sun-synchronous orbit for a NOAA satellite [1].
2
Analysis of Spatiotemporal Variability of Surface Temperature of Okhotsk Sea and Adjacent…
DOI: http://dx.doi.org/10.5772/intechopen.94918
2. Satellite data
3 Composite images of the obtained ocean variables transferred to a standard map, obtained
by averaging over space and time of several flights in the form of a level 2 product. This
product contains data from only one sensor and may have gaps.
4 An image representing the ocean variable averaged within each cell of the spatio-temporal
grid, for the creation of which the gaps in the level 3 product data are filled with data
analysis methods, incl. Interpolation. In the course of the analysis, it is possible to use several
level 2 or 3 products from different sensors, and it is also possible to use data from field
observations or model calculations.
Table 1.
A summary of various satellite data processing levels.
3
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 2.
An example of a ten-day average image with a spatial resolution of 0.25 × 0.25 degrees (left) and 2 × 2 km
(right) [2].
once by the operator, after which it is performed in automatic mode). Then data
with different spatial resolution (the maximum resolution for a set of satellites is
2 km, images with a spatial resolution of 250 m can be obtained on separate chan-
nels) are entered into the database, where the operator manually corrects each
image (eliminates gross errors that most often appear in the case of the presence of
semi-transparent clouds). After that, daily data averaging is performed. As the data
accumulate, ten-day and monthly averaging is also performed (as shown in Figure 2
[2]), as well as the calculation of anomalies (deviations from the long-term average
value in a specified decade or month, obtained over a number of previous years).
The resulting data set should be carefully analyzed. Determination of the char-
acteristics of unidirectional trends, as well as the amplitudes and phases of cyclic
components and the possibility of using them to predict thermal conditions in the
next season constituting the content of this study.
3. SST anomalies
There are two types of SST anomalies as for any parameter distributed in space
and time. The first is the sharp deviation of the specified parameter from the mean
values in the neighboring areas. Anomalies of this kind can be regular and even
present on average long-term maps (Figure 3 [3]). Typical examples from the point
of view of SST are upwelling zones (quasi-stationary), mouths of large rivers, as
well as zones of influence of warm and cold currents. Some anomalies of this nature
may be unstable and dependent on weather conditions. Such anomalies determine
the behavior of aquatic organisms, since each species has its own conditions of
existence and development (some species need colder water, and some are more
comfortable in warmer water).
The second type of SST anomaly is a deviation from the multiyear average.
This type of anomaly depends entirely on the current conditions and reflects the
peculiarities of the distribution of the parameter in a given year. From the point
of view of SST, this may mean more or less intense heating or cooling than usual.
The appearance of such anomalies can lead to a change in the timing of the fish’s
approach to spawning (Figure 3 [3]) or to migration to places with more favor-
able conditions, and in the early stages of development it can lead to death from
unfavorable conditions.
Let us consider in more detail the extraordinary conditions that took place in the
spring of 2011, when the greatest delay in the release of juveniles from the Sakhalin
SFHs was noted. The spatial distributions of sea surface temperature anomalies
4
Analysis of Spatiotemporal Variability of Surface Temperature of Okhotsk Sea and Adjacent…
DOI: http://dx.doi.org/10.5772/intechopen.94918
Figure 3.
An example of the long-term average distribution of SST, the actual temperature in the second decade of July
2018 and the 2018 anomaly relative to the multi-year average (these conditions led to a two-week delay in the
approach of pink salmon to spawning rivers) [3].
Figure 4.
Distribution of surface temperature anomalies (°С) in the third decade of may (A) and in the second decade of
June (B) 2011. Ice conditions off the east coast of Sakhalin Island (right) on June 2, 2011 (snapshot of the aqua
satellite). The strait of Tatary is closed by heavy clouds [4].
5
Satellite Systems - Design, Modeling, Simulation and Analysis
winds prevailed in the study region, which just contribute to the transport of
drifting ice to the southeastern coast of Sakhalin Island, as well as to the estuarine
regions of the Naiba and Ai rivers and to the Mordvinova Gulf.
In spring 2011, the distribution of surface atmospheric pressure also signifi-
cantly differed from the characteristic one in the given period of the year. For
example, in the third decade of May, the high-pressure region was located in
the southern part of the Sea of Okhotsk, above the Kuril Basin (Figure 5A [7]).
Usually this area is located in the eastern part of the sea, and air flows directed
to the north are formed above this basin. When the pressure distribution has this
structure, the atmospheric circulation typical for the summer monsoon is not
formed; in principle, it was weakly expressed, since the pressure gradients were
small and cold waters with low salinity moved freely to the south, transporting
drifting ice.
Subsequently, the high-pressure region shifted from the sea to the Pacific Ocean.
In particular, in the second decade of June, zonal, northeast-oriented air flows
prevailed, which is not typical for the warm season (Figure 5B [7]).
However, it is unlikely that only the peculiarities of the synoptic conditions
can explain the significant negative temperature anomalies of the sea surface in
the northern part of the Sea of Japan, primarily in the main zone of the Tsushima
Current influence, off the west coast of Hokkaido Island (up to 5–6°C). Another
zone with approximately the same negative anomalies was noted in the zone of the
influence of the Amur River runoff in the region of its mouth, in the western part
of the Amur Estuary, and in a significant part of the Sakhalin Gulf, as well as on
the northeastern shelf of Sakhalin Island (up to about 52° N). Pronounced nega-
tive anomalies were observed both on the southwestern and southeastern coasts of
Sakhalin Island, as well as in the Aniva Gulf. It is difficult to explain the cause of
such a heat deficit in the surface waters of a very wide basin: one of them could be
a decrease in the intensity of the Tsushima Current, whose quasi-periodic fluctua-
tions (a cycle with a period of 5–6 years is especially distinguished) have a signifi-
cant influence on the climate of the Sakhalin-Kuril region [8]. It is also difficult to
assess how often such situations can be observed which have led to a delay in the
release of juveniles from Sakhalin SFHs lasting from 3 weeks to a month. Probably,
due to the low temperature of coastal waters observed for a long time, there could
be a mass death of juveniles of natural origin which rolled along the spawning rivers
of the island.
Figure 5.
Distribution of surface pressure (hPa) over the Sea of Okhotsk in the third decade of may (A) and the second
decade of June (B) 2011 [7].
6
Analysis of Spatiotemporal Variability of Surface Temperature of Okhotsk Sea and Adjacent…
DOI: http://dx.doi.org/10.5772/intechopen.94918
4. Linear trends
When analyzing long-term data series, a question arises – are there unidirec-
tional tendencies of one or another parameter in the specified period of time?
This is especially important in the context of global warming, the most noticeable
results of which are an increase in winter air temperature in the Arctic latitudes and
a decrease in ice cover (area of the water area covered with ice) both in the Arctic
Ocean and in the Sea of Okhotsk.
In order to answer this question, in each spatial cell of the 1000 × 1022 matrix,
covering the investigated water area (42–60°N and 135–163°E), the linear trend
coefficients were calculated using the least squares method, meaning the tempera-
ture increase over the year … The calculation was carried out on a 20-year series of
data for each month separately, for annual averages, as well as for average tempera-
tures for the season. The seasons for the calculations are shifted by a month relative
to the calendar ones due to the specifics of the region: January–March (winter),
April–June (spring), July–September (summer) and October–December (autumn).
The linear trend coefficients calculated from the full annual series in the study
region are mostly negative, and there are no regions with positive values in the Sea
of Okhotsk (Figure 6 [9]). Significant cooling of the surface layer (about 1.5°C over
10 years) was observed in the northern and western parts of the sea; most clearly
this process is manifested in the northern part of the sea at some distance from
the coast. A less pronounced decrease (at a rate of about 0.5°C over 10 years) is
observed in the Kuril Islands and in the band along the western coast of Kamchatka
beyond the shelf zone, as well as in the influence zone of the Amur River outflow
(Amur Estuary, Sakhalin Bay, northern shelf of Sakhalin Island).
A significant decrease (more than 0.5°C over 10 years) was found in the north-
ern part of the Strait of Tatary; further southward, the cooling rate of surface
waters decreases and, south of 45° S, a slight increase in the temperature of the
surface waters in the Sea of Japan was recorded. The warming process is more
pronounced in the northwestern part of the Pacific Ocean, i.e., in the southeastern
part of the study region.
Let us now consider the results of calculating the parameters of the linear trend
in different seasons. In winter, there is not enough data for a reliable calculation
in the northwestern part of the sea (the “Sea of Okhotsk fridge”), on the northern
shelf of the sea, and on the northeastern shelf of Sakhalin. In many regions of the
basin, where the ice cover is significantly influenced, the calculation is not reliable
enough; nevertheless, some ideas of the trends in the thermal regime during the
cold period can be put forward.
In the entire Sea of Okhotsk, especially in the northern and western parts of the
sea, as well as in the Strait of Tatary of the Sea of Japan and in the regions of the
Pacific Ocean adjacent to the Kuril Islands, there is a tendency for the temperature
decrease in the surface layer. Signs of the opposite trend are recorded only in the
southeastern part of the study region, in the Pacific Ocean, in the northern part
of the Sea of Japan, and in the very southern part of the Strait of Tatary (south
of 47° N).
The most pronounced decrease in the surface layer temperature is in the
springtime. It covers the entire study region, including the northwestern part of the
Pacific Ocean, and it is particularly significant in the northern and western parts of
the Sea of Okhotsk (the rates of temperature decrease range from 1.0 to 1.5°C over
10 years and even higher in some regions).
The cooling of the surface layer in spring is the most evident consequence of
a decrease in ice cover, since in the absence of the ice cover the cooling processes
7
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 6.
Distribution of linear trend coefficients by number of average annual (above) and average seasonal values. The
scale corresponds to degrees celsius over 10 years [9].
cover a layer of water thicker than in the presence of the ice cover. Undoubtedly, this
has a significant impact on the climate of the region and, above all, on the weather
conditions of Sakhalin, which are the most affected by this process.
While calculating the trends in sea surface temperature in individual months, it
was found that the most intense cooling is observed in the month of May, especially
in the northwestern part of the Sea of Okhotsk and along the entire eastern coast
of Sakhalin Island from Cape Elizabeth in the north to Cape Aniva in the south.
Exceptions are the western part of the Amur Estuary, which accounts for the bulk
of the Amur flood, and the basin adjacent to Tauiskaya Bay. It is obvious that the
hydrological regime of coastal areas experiencing the influence of the river outflow
differs from that typical for the Sea of Okhotsk as a whole.
In summertime, the region is dominated by an increase in the sea surface tem-
perature most pronounced in the northwestern part of the Pacific Ocean. In the Sea
of Okhotsk, it is noted in its northwestern part, in the region of the Shantar Islands,
8
Analysis of Spatiotemporal Variability of Surface Temperature of Okhotsk Sea and Adjacent…
DOI: http://dx.doi.org/10.5772/intechopen.94918
in Sakhalin Bay, and in the Amur Estuary, to the east of the northern tip of Sakhalin
Island, at the northwestern coast of the Kamchatka Peninsula, and near the Urup
and Simushir islands.
In the northern part of the sea, near the Kuril Islands, on the southeastern coast
of Sakhalin Island and along the Primorye coast in the Sea of Japan, there is a trend
towards a decrease in temperature, although it is more moderate than in springtime.
Moreover, in July, the tendency towards a decrease in the sea surface temperature
still prevails in the study region; the trend changes only in the northwestern part
of the Pacific Ocean. In August, in these regions, warming reaches the highest rates
(about 1.5°C over 10 years); the temperature increase was noted in the central part
of the Sea of Okhotsk, along the western coast of Kamchatka and off the southeast-
ern coast of Sakhalin, in the Strait of Tatary and the Amur Estuary. In the northern
part of the sea, on the northeastern shelf of Sakhalin, and in the Kuril region, the
trend towards a decrease in sea water temperature continues.
In September, on the contrary, the intensity of warming in the Pacific Ocean
decreases; the greatest warming rate was noted in the northern and northwest-
ern parts of the sea (especially in the shelf zone between Ayan and Okhotsk
settlements).
A diffuse pattern is observed in the fall. Relatively small negative trends prevail
in the greater part of the Sea of Okhotsk; positive trends (also insignificant) are
noted in the eastern part of the sea on the western shelf of Kamchatka and in a rela-
tively narrow band along the northern coast. Smaller areas were also found in the
Shantar Islands and Kashevarov Banks and in the southwestern part of the analyzed
basin of the Sea of Japan.
The highest rates of the temperature decrease in the surface layer are observed
in the Strait of Tatary, in its northern part, and along the western coast of Sakhalin,
as well as in the Sea of Okhotsk, in the northern part of it, beyond the shelf, off the
southeastern coast of Sakhalin, and also in the band between 145° and 150° E from
the Urup and Iturup islands to the Kashevarov Banks.
It was found in the calculations for individual months that a positive trend
prevailed in October, and the trend reversed in November and December.
In the work [10], it was found that, in the Sea of Okhotsk and in the region of
the Kuril Islands, negative trends were also noted in the spring months; the largest
positive trends were recorded in October. Similarity results become very interesting
if we take into account the differences in the studied time periods and spatial char-
acteristics of the regions in which the trends were calculated. In the Sea of Japan, no
negative trends were identified in this work.
As a result of the analysis of data on the surface temperature of the Sea of
Okhotsk over a 20-year period (1998–2017), it was found that the global climate
change in this basin caused a decrease in the ice cover and a decrease in the tem-
perature of the upper water layer during the winter–spring period. The negative
trends in temperature in the spring in the northern and western parts of the study
region, as well as in the Strait of Tatary of the Sea of Japan (from 0.5° to 1.5°C over
10 years), are especially large. This effect in the reduction of the ice cover, both
in time and in space, is the most evident and can be explained by an increase in
the depth of winter convection. The predominance of a decrease in sea surface
temperature, although less pronounced, was also noted in winter and autumn, and
generally throughout the year.
In summer, the region is dominated by an increase in sea surface temperature,
most pronounced in the northwestern part of the Pacific Ocean. In the Sea of
Okhotsk, it was noted in its northwestern part, east of the northern tip of Sakhalin
Island, off the northwestern coast of the Kamchatka Peninsula, as well as in some
other regions. Moreover, in July, the main role belongs to the processes of cooling
9
Satellite Systems - Design, Modeling, Simulation and Analysis
of the surface layer; the change in the trend occurs in August and continues in
September–October.
These processes play a significant role in the climate variations in the Sea of
Okhotsk; in particular, a decrease in the temperature of sea water is noted in the
coastal waters of Sakhalin Island. In addition to the weather conditions, the results
are important for studying the habitat conditions of commercial fish species and
invertebrates in the basin, which is of great fishery importance.
5. Periodic fluctuations
If you look closely at the graph of the course of the mean monthly SST (Figure 7
[11]), then in addition to seasonal fluctuations, you can notice interannual varia-
tions. They are especially pronounced in August, when the mean monthly water
temperature reaches its maximum value, and are expressed in the modulation of the
annual harmonic. Moreover, these oscillations are of a quasiperiodic nature, most
likely associated with certain phenomena in the atmosphere and hydrosphere.
Any periodic oscillations can be described by knowing their amplitude, phase
and period. Using the least squares method, you can find the corresponding ampli-
tude and phase for each selected period. Thus, it is possible to establish a kind of
“influence zones” of harmonics with a certain period, i.e. areas in which the ampli-
tude of one or another harmonic exceeds a certain threshold value. And in order
not to interfere in the calculations with long-term components (the period of which
exceeds half the length of the series), the trend obtained according to the method
described in the previous section was subtracted from the initial data.
The distribution of the amplitudes of temperature fluctuations in the studied
region is rather complex. Having determined in each spatial cell the period cor-
responding to the largest amplitude and displaying the obtained data on the screen,
we have established several main periods that play a significant role in interannual
SST variations in most cases. For a series of 21 years, it is not entirely correct to
calculate fluctuations with a period exceeding 11 years. Short-period fluctuations
are unstable and generally have little information. Thus, the spatial distributions of
harmonic amplitudes with a period from 3 to 11 years were considered in detail.
In most of the studied water area, the main role was played by variations with
a period of about 5.5–6 years (on the graphs for points 2, 5 and 8 in Figure 8,
you can see that the largest amplitude at these points corresponds to a period of
66–68 months), the spatial distribution of the amplitude with a period of 6 years is
shown in Figure 9 [12]. The zone of its influence is the most extensive and occupies
Figure 7.
Annual variation of average monthly temperature (left) and average temperature in august (right) [11].
10
Analysis of Spatiotemporal Variability of Surface Temperature of Okhotsk Sea and Adjacent…
DOI: http://dx.doi.org/10.5772/intechopen.94918
the southern half of the Sea of Okhotsk and the northern part of the Sea of Japan,
up to the Amur estuary. The amplitude of this harmonic in the zone of its influence
ranges from 1 to 2°C. In the vicinity of the Kuril Islands, off the western coast of
about. Hokkaido, off the northeastern coast of about. Sakhalin, the amplitude is
slightly lower (from 0.5 to 1°C). Further to the north, starting from 52 N, its influ-
ence decreases and practically disappears.
In [8], based on the EOF decomposition of the sea surface temperature in the
North Pacific Ocean, it was shown that an oscillation with a similar period is charac-
teristic of the entire region influenced by the Kuroshio Current and its branch, the
warm Tsushima Current. This is evidenced by the large amplitudes in the zone of
the indicated currents and in Figure 9. Consequently, as a result of this study, it was
possible to estimate the boundaries of the influence of this component in the Sea of
Okhotsk, which runs parallel to the islands of the Kuril ridge and divides this basin
into two practically equal parts. The entire northern Sea of Japan is significantly
affected by this cycle.
As seen from Figure 8, at most points, one can also note peaks in the range of
periods from two to three years. The highest values of the amplitude of the 3-year
harmonic can be noted in the northwestern part of the Sea of Okhotsk, at a distance
from the coast, as well as near the northwestern coast of Kamchatka and in the strip
from 47 to 490 N. and from 147 to 1490 E. in the area of the Kuril deep-water basin
and in the northwestern part of the Pacific Ocean. However, the amplitudes of these
oscillations are somewhat lower than those of the six-year harmonic (from 1 to
1.5°C). At a distance from these regions, the amplitude gradually decreases to zero.
The spatial distribution of the amplitudes of the cyclic component with a period
of 5 years differs markedly from that considered above for a period of 6 years. The
zone of its influence is noticeably narrower, it is concentrated mainly on the north-
ern shelf of Hokkaido, in the region of the South Kuril Islands (vast waters both on
Figure 8.
Examples of graphs of the dependence of the amplitude of the harmonic (in °C) from its period (in months).
The location of the points is shown in Figure 9 [12].
11
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 9.
Distribution of amplitudes (in °C) of interannual fluctuations in the mean monthly sea surface temperature
(august) with a certain period (indicated in years in the upper left corner of the image) [12].
the Sea of Okhotsk and on the ocean sides), and, surprisingly, on the northeastern
shelf of Sakhalin, where the influence of a lower frequency component was not
noted. In the Tatar Strait, its role is also noticeable, but expressed to a lesser extent
than the 6-year harmonic.
In the western part of the Sea of Okhotsk and the northwestern part of the
Pacific Ocean, a cyclical component with a period of about 8 years is significantly
manifested. The zone of influence of the Amur river runoff in summer is clearly
distinguished in the spatial distribution - the Amur estuary, the southern and
eastern parts of the Sakhalin Bay, the area between the Schmidt Peninsula and the
Kashevarov Bank [13]. It is interesting that in the area of the Kuroshio Current
manifestation, this component has large amplitudes, in the Tsushima Current zone
in the Sea of Japan - insignificant, while on the northern shelf of Hokkaido and on
the Sea of Okhotsk side of the Southern Kuril Islands, where the warming effect of
the Soya Current affects, the amplitudes are significantly.
The lowest frequency of the considered harmonics with a period of 11 years
is manifested in the northern part of the Sea of Okhotsk; in other parts of the
study area, its role is insignificant. It is rather difficult to put forward a reasonable
hypothesis that could explain such significant differences in the very long-term
variations in SST in different parts of the same basin. It can only be assumed that
12
Analysis of Spatiotemporal Variability of Surface Temperature of Okhotsk Sea and Adjacent…
DOI: http://dx.doi.org/10.5772/intechopen.94918
due to the comparative shallowness of the northern region, the effect of the winds
of the southern rumba and the greater number of sunny days than in the southern
part, due to the lesser influence of cloudiness, the influence of the solar cycle is
more noticeable here.
Attention is drawn to the fact how the zones of manifestation of harmonics shift
with a period of 5 to 11 years. If the zone of influence of the 5-year harmonic is
focused near the islands of Sakhalin, Hokkaido and the Southern Kuriles, then with
an increase in the period, the region with the highest amplitude shifts clockwise
(towards the northeastern coast of Sakhalin and further to the northern part of the
Sea of Okhotsk). This interesting fact is also difficult to give a reasonable explana-
tion, it requires additional study.
In work [8], a method was developed for predicting thermal conditions for a
year in advance in certain areas of the studied water area (this method was also
used to recover data gaps associated with the influence of cloudiness or technical
reasons), which consisted in calculating the temperature in a given square in time t
according to the formula:
N
T ( t ) = at + b + ∑ck cos (ωk t − ϕ k ) (1)
k =1
where a and b are the parameters of the linear trend, 𝑐𝑘 are the amplitudes, and
𝜑𝑘 are the phases of the cyclic components (harmonics) of sea surface temperature
variations. An essential feature of the method is the fact that the amplitudes and
phases of the main cyclical components are calculated by the least squares method,
with their periods ranging from 18 to 144 months with a step of 1 month. For each
cell, a set of 3–4 harmonics was determined, which make the largest contribution to
the interannual variations in SST. Since they are not orthogonal, for forecasting
using formula ( 1 ), it is necessary to subtract the calculated wave from the initial
series before determining the parameters of the next one in order to avoid double
inclusion of coherent components (in [14, 15]) such a technique was called
“Sequential spectra method”).
Based on the parameters of the obtained cyclic components, a retrospective
forecast of thermal conditions for the summer of 2018 was carried out (observa-
tional data for 1998–2017 were used to calculate the parameters of harmonics and a
linear trend). The calculation was carried out for each spatial cell according to
formula ( 1 ), taking into account the trend and four harmonic components with the
highest amplitudes. For the forecast, periods from 18 to 144 months were covered.
The calculation results are presented in Figure 8 in the form of graphs of forecast
curves and real variations in sea surface temperature, including the predicted values
that took place in the summer of 2018. Figure 10 shows the spatial distribution of
the difference between the predicted and actual temperatures for August 2017 and
2018. The forecast was built for a year ahead along the entire previous series.
The curves, which are the sum of the trend and the first four harmonics, gener-
ally repeat the actual interannual temperature fluctuations. The correlation coef-
ficient of the initial and predicted series at the selected points exceeds 90%. Note
that even the first two harmonics in many cases provide a correlation coefficient
of more than 70%. Despite the fact that 2018 was anomalous in terms of thermal
conditions (the Tsushima Current and its Okhotsk branch of the Soya Current were
weakened, a heat deficit was felt in the zone of influence of the Amur River runoff),
and in some other areas, even in such water areas, the forecast can be considered
acceptable. An example of a similar situation is given for the Tatar Strait, where
the predicted value was higher than the actual one, but the general course was
predicted correctly, and the error was not so great. For the northern part of the Sea
of Okhotsk, the northeastern shelf of Sakhalin Island, and a number of other areas,
13
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 10.
Examples of temperature forecast graphs (in °C) for the next year. The dashed line shows the prognostic curve.
The forecast is carried out for the period from 1998 to 2017. The actual temperature in 2018 is marked with a
cross [11].
good agreement was observed between the calculated and real values of the surface
layer temperature.
Let us consider some of the parameters of the graphs below. The standard devia-
tion of the initial and predicted series correspond to each other and range from
1.5°C (Tatar Strait) to 2°C (South Kuriles). The average displacement of the pre-
dicted series relative to the initial one ranges from 0.4 to 0.6°C. The forecast error is
1.7°C in the Tatar Strait, 1.3°C near the Southern Kuriles, 0.3–0.4°C in the northern
part of the studied water area.
More detailed studies devoted to predictability and the limits of applicability of
the approach used will be carried out later. However, we can already say that for a
significant part of the Sea of Okhotsk regions and adjacent water areas, the forecast
of the surface layer temperature with a one-year lead time is quite successful,
although the abnormally cold temperatures that took place in a number of areas in
2018 are rather difficult to predict.
Let us take a closer look at Figure 11. The forecast for August 2017 turned out to
be quite successful, the discrepancy between the actual and predicted temperatures
in most of the water area does not exceed ±2°C, with a standard deviation of SST
of about 1.5-2°C (only in the northwestern part of the Pacific Ocean is the standard
SST deviation is within 2-4°C). At the same time, the forecast for August 2018
contains a large area within which the temperature estimate was greatly overesti-
mated (over 4°C). The map of SST anomalies for August 2018 [16] also contains
areas of low temperatures (3-4°C lower than the average multiyear norm), which
coincide in space with areas of unsuccessful forecast. This area is located in the zone
of influence of the Tsushima Current, and the forecast inaccuracies are due to the
fact that the weakening of this current occurred two years earlier than the expected
date. Indeed, in Figure 8, we see that in the Tatar Strait and the Southern Kuriles,
the distance between two neighboring SST minimums decreased to 3–4 years, while
its quasiperiodic oscillations with a period of about 6 years are described in the
literature [17].
As a result of the analysis of the data set on the surface temperature of the Sea
of Okhotsk and adjacent waters, the main cyclical components responsible for
the interannual variations of this parameter and the “zones of influence” of each
14
Analysis of Spatiotemporal Variability of Surface Temperature of Okhotsk Sea and Adjacent…
DOI: http://dx.doi.org/10.5772/intechopen.94918
Figure 11.
Difference between predicted and actual temperatures [11].
harmonic were determined. It is shown that the main contribution to these varia-
tions comes from components with a period of about 6 years, as well as 3, 5, 8, and
11 years.
The zone of influence of the fundamental harmonic is the most extensive and
occupies the southern half of the Sea of Okhotsk and the entire northern part of the
Sea of Japan; its amplitude is within 1–2°C. In the vicinity of the Kuril Islands, off
the western coast of Hokkaido Island, off the northeastern coast of Sakhalin Island,
the amplitude is slightly lower (0.5–1°C), and in the northern part of the Sea of
Okhotsk its influence is insignificant. Most likely, this component is associated with
fluctuations in the Kuroshio Current and its branch, the Tsushima Current [17].
The highest values of the amplitude of the 3-year harmonic (1–1.5°C) can be
noted in the northwestern part of the Sea of Okhotsk, at a distance from the coast,
as well as off the northwestern coast of Kamchatka, in the region of the Kuril deep-
water basin and in the northwestern parts of the Pacific Ocean.
The area of manifestation of the component with a period of 5 years is notice-
ably narrower, it is concentrated mainly on the northern shelf of Hokkaido, in the
region of the South Kuril Islands (both from the Sea of Okhotsk and the ocean
side), and on the northeastern shelf of Sakhalin. In the Tatar Strait, its role is also
noticeable, but expressed to a lesser extent than the 6-year harmonic.
In the western part of the Sea of Okhotsk, in the zone of influence of the Amur
River runoff, as well as in the northwestern part of the Pacific Ocean, a cyclical
component with a period of about 8 years is significantly manifested.
The lowest frequency of the considered harmonics with a period of 11 years is
manifested in the northern part of the Sea of Okhotsk; in other parts of the study
area, its role is insignificant.
Together with the parameters of the linear trend [9], the amplitudes and phases
of the main cyclical components (in each spatial cell, 4 harmonics with the high-
est amplitudes were used) can be used to predict thermal conditions for the next
summer. The retrospective calculation for 2018 gave generally satisfactory results,
despite the abnormally cold conditions of this year, noted in a number of areas of
the studied water area. The possibility of predicting thermal conditions is of practi-
cal importance, primarily for assessing the conditions for the approach of Pacific
salmon to spawning. And the results obtained show that in some areas of the water
area, one can count on a fairly accurate forecast even for such an unstable parameter
as the ocean surface temperature.
In general, the success of the forecast is influenced by how pronounced the
cyclical components with a certain period in a given area. Identification of the
zone of influence of various harmonics allows you to determine the boundaries of
the regions in which the applicability of this method can be expected. As for the
15
Satellite Systems - Design, Modeling, Simulation and Analysis
accuracy of the forecast, one should pay attention to the presence of large areas in
which the modulus of the difference between the predicted temperature and the
actual one was of the order of two standard deviations. This fact shows that this
method does not guarantee the success of the forecast in cases where strong tem-
perature anomalies are observed. It can only be used to obtain a primary estimate of
the ocean surface temperature (or another parameter that experiences quasiperiodic
oscillations) based on a sufficiently long series. The 21-year series is not enough to
estimate low-frequency components (with a period of 30–50 years), which could
also affect the quality of the forecast. For a better forecast, you can combine this
method with an assessment of the current conditions, adjusting the forecast in
the direction of increasing or decreasing temperature, depending on the current
meteorological conditions.
6. Conclusion
This chapter shows how satellite data, going through all levels of processing,
become useful products for both science and applied environmental prediction
problems. Assessment, and later predicting the dynamics of various environmental
parameters, will help reduce damage to nature and identify the degree of influence
of anthropogenic factors on the future of the Earth. Based on our analysis results, it
seems that the existence of global warming phenomenon has already been proven,
but its influence on different parts of the planet is heterogeneous and is not traced
in all seasons in the same way. In the spring, there is a tendency to a decrease in the
temperature of the surface layer of the Sea of Okhotsk, this tendency continues
until August, where it changes to warming condition.
Periodic fluctuations of physical parameters are also not homogeneous in space.
It is possible to identify areas where the main oscillation period is 3, 5, 6, 8 and
11 years. Using the sequential spectra method gives a relatively good estimate for
the next year’s ocean surface temperature. However, to improve the quality of the
forecast, it is necessary to improve the quality of the initial data, comparing data
from various sources and filling in the gaps associated with the presence of clouds
and ice cover, as well as improve the interpolation algorithms to avoid losing sight of
mesoscale phenomena.
Author details
Dmitry Lozhkin
Russian Federal Research Institute of Fisheries and Oceanography,
Sakhalin branch, Yuzhno-Sakhalinsk, Russian Federation
© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License (http://creativecommons.org/licenses/
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
16
Analysis of Spatiotemporal Variability of Surface Temperature of Okhotsk Sea and Adjacent…
DOI: http://dx.doi.org/10.5772/intechopen.94918
References
[3] Tshay Zh. R., Shevchenko G. V. [10] Shatilina TA, Tsitsiashvili GS,
Otsenka ekstremal’nykh termicheskikh Radchenkova TV. Assessment of
usloviy v period neresta the statistical significance of water
tikhookeanskikh lososey u poberezh’ya temperature variability in the sea of
o-va Sakhalin [assessment of extreme Japan and the Northwest Pacific in
thermal conditions during the spawning 1982-2007, Resursy kolyuchego kraba,
period of Pacific salmon off the coast of perspektivy ispol’zovaniya i usloviya
Sakhalin Island] // Geodinamicheskiye obitaniya: Trudy Sakhalinskogo
protsessy i prirodnyye katastrofy, nauchno–issledovatel’skogo instituta
Yuzhno-Sakhalinsk, May 27-31, 2019, P. rybnogo khozyaistva i okeanografii
184. (In Russian). (Brown king crab resources, usage
prospects and habitat conditions:
[4] Lozhkin DM, Tshay ZR, Transactions of the Sakhalin
Shevchenko GV. Satellite monitoring of research Institute for Fisheries and
temperature conditions near the mouths Oceanography). Yuzhno-Sakhalinsk:
of spawning Rivers in the southern SakhNIRO. 2011;12:180-190
part of Sakhalin Island. Izvestiya -
Atmospheric and Ocean Physics. [11] Lozhkin D. M., Shevchenko G.
2019;55(9):1166-1171 V. Lineynyye trendy i tsiklicheskiye
variatsii temperatury poverkhnosti
[5] Vlasova GA, Vasil’ev a S, Okhotskogo morya i prilegayushchikh
Shevchenko GV. Prostranstvenno– akvatoriy po sputnikovym dannym
Vremennaya izmenchivost’ Struktury [linear trends and cyclical variations in
i Dinamiki Vod Okhotskogo Morya the surface temperature of the Sea of
(Spatial and Temporal Variability of the Okhotsk and adjacent waters according
Structure and Dynamics of Waters of to satellite data] // VII nauchno-
the Sea of Okhotsk). Moscow: Nauka; prakticheskaya konferentsiya molodykh
2008 uchonykh s mezhdunarodnym
uchastiyem «Sovremennyye
[6] Shershneva OV, Shevchenko GV, problemy i perspektivy razvitiya
Novinenko EG. Temperature conditions rybokhozyaystvennogo kompleksa». -
in the areas of juvenile release from 2019. – pp. 257-262. (In Russian).
salmon fish-breeding plants from
the mouth of rivers in Sakhalin [12] Lozhkin DM, Shevchenko GV. Cyclic
and the Iturup Island. Izv. TINRO. Sea surface temperature variations in
2007;150:217-225 the sea of Okhotsk and adjacent water
17
Satellite Systems - Design, Modeling, Simulation and Analysis
[16] http://www.sakhniro.ru/
18
Chapter
Abstract
In the context of this chapter, a satellite is a spacecraft (SC) that orbits around a
celestial body such as the earth. A spacecraft has several design constraints placed
upon it before it can be placed in an orbit around the intended celestial body. First,
satellite designs are limited in their mass and volume to fit on the launch vehicle
that places them into orbit. Secondly, the mass and volume limits affect the size of
the power system on the spacecraft; therefore, the amount of power available to
the satellite is also limited. In addition, the space environment (thermal, radiation,
atomic oxygen, space debris, micrometeoroids, etc.) imposes constraints on the
design such as parts and material selection.
A spacecraft is consisted of two parts: the spacecraft bus and the payload (PL)
[1, 2]. The spacecraft bus provides control of the satellite and support services to the
mission payload, while the mission payload provides the mission part of the satel-
lite including payload control, mission data processing, and mission data downlink
dissemination. Examples of mission payloads (or payloads or PLs) are: scientific
instruments, remote sensing instruments, navigation service transmitters, or
communications equipment. A satellite may have one type of PL or a combination
of payload types to accomplish its mission such as navigation, remote sensing, and
communications. Shown below in Figure 1 is a typical imaging satellite used for the
remote sensing mission. Note the clear separation between the spacecraft bus that
provides solar power and maneuvering capability via thruster, while the payload
consisting of the camera and supporting communication devices such as antennas
and guidance devices such as star trackers.
1
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 1.
A typical satellite with bus and payload separation.
Figure 2.
The three main segments for satellite system.
Regardless of the mission type1 and the payload that a spacecraft carries, a sub-
system that must exist in all satellites is the communication subsystem that enables
the spacecraft to communicate with the ground stations that control the satellite and
to deliver the data that the mission requires. This chapter focuses on architecture and
functionalities of the communications subsystem that usually resides on the satellite.
1
Mission type has different meaning depending on context. For example, U.S. military satellite has
three basic mission types: imaging/sensing; communications; positioning, navigation and timing (PNT)
missions. For NASA space exploration, they are near-earth and deep-space missions.
2
Communication Subsystems for Satellite Design
DOI: http://dx.doi.org/10.5772/intechopen.93010
There are three specific segments shown in Figure 2 below that must work
together for the larger overall system to provide communication, navigation, or any
other type of missions:
• The satellite control (or control) segment consisting of all the personnel, facili-
ties, and equipment that are used to monitor and control all the assets in space.
Practically, the control segment is also referred to as satellite ground segment
because it is usually located on the ground.
• The user segment consisting of all the individuals and groups who use and
benefit from the data and services provided by the payloads of the satellite and
the equipment that allows this use.
In general, the space mission dictates the type of orbit2, satellite design and its
expected life cycle, and its operational scenarios. The PL design includes dimen-
sions, interfaces, weight, physical characteristics, and basic utility needs (e.g.,
power consumption), which usually influences spacecraft (SC) bus design. The PL
is often a unique and one-of-a-kind design tailored to meet specific mission require-
ments, frequently relying heavily on newer technology, while the satellite bus has
the supporting function, and as such relies largely on existing or modified hardware
such as batteries, inertial devices, and star trackers. Since PLs and their missions
vary widely, so is this satellite bus supporting role.
Traditionally, the PL is considered a subsystem of the satellite bus that is
designed to generally satisfy the corresponding mission requirements. The PL
operational requirements sometimes impose specific requirements on the satellite
bus that must be satisfied for the PL to accomplish its mission. This interdepen-
dence between satellite bus and PL subsystems has historically resulted in many
nonstandard interfaces developed and implemented by the incumbent spacecraft
builders. As a result, the aerospace industry has been moving toward a more
standardized and commodity satellite bus framework that can potentially result in a
tremendous cost saving approach.
As shown in Figure 3 below, a satellite bus typically consists of the following
subsystems: command and data handling subsystem (C&DHS); communications
subsystem (CS); electrical power subsystem (EPS); propulsion subsystem (PS);
thermal control subsystem (TCS); attitude control subsystem (ACS) also known as
guidance, navigation and control (GNC) subsystem; structures and mechanics sub-
system (S&MS); and life support subsystem for manned missions if required. The
C&DHS will be described in detail below. The CS provides the satellite bus with the
necessary communication functionalities to connect the user and ground segments
to different satellite subsystems. The EPS provides the electrical power generation
and distribution for various spacecraft subsystems. The PS provides maneuvers
2
There are three main types of satellite orbits: low earth orbit (LEO) of 2000 km in altitude or less;
geostationary (GEO) with altitude around 35,786 km; and medium earth orbit (MEO) with altitude
between LEO and GEO.
3
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 3.
A typical satellite bus and payload subsystem.
Figure 4.
A typical and generic sensor payload.
4
Communication Subsystems for Satellite Design
DOI: http://dx.doi.org/10.5772/intechopen.93010
(RFs). These sensors can be classified as passive or active, and each of them can
be further classified as imaging or sounding3. Figure 4 below illustrates a generic
imaging PL that will convert the sensor analog data into electrical signals that can
be captured and transmitted to a ground station. Note the existence of a communi-
cation subsystem as part of this imaging payload.
At the physical layer, the communications subsystem starts with an antenna and
the RF front-end transceiver. The antenna is the most important component of the
communications subsystem where the electromagnetic (EM) signals are originated
or received. The RF front-end/back-end is where the EM signal is being down/up-
converted to baseband/RF signal to be demodulated/modulated for baseband signal
recovery or downlink transmission, respectively. Figure 5 below depicts a typical
transmitter and receiver (transceiver) chain with the modulation and demodulation
(MODEM), followed by the RF front-end and the antennas. The baseband commu-
nications function is carried out by the MODEM, whereas the RF portion is handled
in the transceiver, RF front-end, and antenna sections.
Modulation is the name given to the process of impressing the wanted signal to
be transported onto a radio frequency (RF) carrier, which is then conveyed over the
satellite link and demodulated at the receiving terminal to extract the wanted signal
from the carrier. Thus, modulation translates a baseband spectrum (at zero fre-
quency) to a carrier spectrum (at RF range) and demodulation is the process of recov-
ering the data at the receiver end of the link. Thus, the process requires a modulator
Figure 5.
Typical RF front-end chain.
3
In this context, sounding means sending a radio signal and interpreting the results from the returned
signals. Examples are radar and ranging signals.
5
Satellite Systems - Design, Modeling, Simulation and Analysis
and a demodulator, collectively known as a MODEM. The input to the modulator may
require some initial processing such as filtering and amplitude limiting.
Before the RF signal is sent to the antenna, a traveling wave tube amplifier
(TWTA) or solid-state power amplifier (SSPA) is needed to amplify the RF signal
to a desired level for transmission. Conversely, after the RF signal is received by the
antenna, a low noise amplifier (LNA) is needed to ensure that the received signal is
brought up to the desired signal level with minimum noise before demodulation.
In addition to being lighter than TWTA, the achievable power efficiency for
SSPAs is a major factor to support transmit phased arrays. Currently, the tube-based
TWTA implementations are still the most cost-effective design, even though both
options might be viable for lower power systems.
In increasing technical maturation over the years, the following types of space-
craft antennas have been used for satellite communications:
2. Increased gain types of satellite antennas (horn type and helix antennas) for
medium earth coverage.
5. Phased array feed and phased array antennas for scanning and hopping beams.
6. Optical communications systems, which have been used for intersatellite links
and interplanetary communications, and increasingly being considered for
earth-to-space systems.
In general, there are many different types of antennas, but the one most com-
monly associated with satellite communications is the parabolic dish antenna. These
dish antennas have a narrow beam width, concentrating the energy of the radiated
main beam into a smaller solid angle. This means more of the radiated energy
reaches, or “illuminates,” the satellite when using a dish antenna as compared to an
omnidirectional, or “omni” for short, antenna. An example of dish antenna used on
satellite is shown below in Figure 6 for a Ku-band space to ground antenna (SGANT)
mounted on the external stowage platform of the International Space Station (ISS).
There are several factors driving the design and development of satellite antennas.
These include the need to reuse frequency bands because of limited spectrum alloca-
tions; the need to have antennas that can operate at higher frequencies with higher
bandwidth; and the desire to deploy higher gain antennas at the same time minimizing
the required size, weight, and power (SWAP) constrains. In practice, there are sub-
stantially more SWAP constrains for satellite antennas than on the ground stations, and
this results in several design trade-offs between the space and control/user segments.
For example, the GEO orbit allows a high gain antenna to be pointed at a satellite
with a minimum of tracking. Thus, a large dish can be used and remain virtually
stationary without tracking a satellite as it moves around in its orbit. On the other
hand, a low earth orbit (LEO) satellite that can cross from horizon to horizon in a
few seconds can result in ground antenna installations that can be quite complex
and expensive. Consequently, trade-offs need to be made to support the mission
parameters of the whole satellite network.
6
Communication Subsystems for Satellite Design
DOI: http://dx.doi.org/10.5772/intechopen.93010
Figure 6.
Example of a satellite dish antenna.
The term “command and data handling subsystem” (C&DHS) was referred to as
“On-board Computer” (OBC), which is a legacy of the past in which many satellite
functions were performed by analog circuits with the help of an OBC. With the
current shift toward the digital domain, the term OBC does not fully cover the topic
anymore thus C&DHS is being used instead. An appropriate analogy to describe the
C&DHS subsystem is to regard it as the brain and nervous system of the spacecraft.
The function of a C&DHS subsystem is to perform onboard processing and
operations and internal communication [3, 4]. The task of managing the opera-
tions of the spacecraft subsystems is nowadays performed mostly by software in
an autonomous manner and is generally categorized as onboard operations. The
software is also responsible for preparing the data to be downlinked and handling
any commands that are received from satellite operators on the ground. Lastly, the
C&DHS facilitates and controls all internal communications (consisting of com-
mands, telemetry, and tracking data) between the different satellite subsystems.
The basic functions of the C&DHS can be summarized below:
1. Receives commands from the command or user segment through the telem-
etry, tracking, and control (TT&C) subsystem.
3. Collects and formats telemetry data from all space vehicle (SV) units.
4. Distributes telemetry for downlinking. Provides a platform for bus flight soft-
ware (FSW).
7
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 7.
Block diagram of a typical command and data handling subsystem.
common low-speed data bus in red color, typically compliant with MIL-STD 1553
or other standards. Also shown is the data connection in blue from the C&DHS to
other components, which is more customized and high-speed in nature depending
on the design.
The heart of the system is the C&DHS’ onboard computer (or OBC) that runs
the software responsible for managing the onboard operations. The OBC is tightly
linked to the electrical power subsystem (EPS). The main reason is the importance
of the available and consumed power for managing onboard spacecraft operations.
For instance, by continuously querying the EPS on the available power, the OBC
can decide to turn off non-critical subsystems to prevent vital systems from shut-
ting down from lack of power. Secondly, the OBC must be able to command the
EPS to disable or enable different subsystems throughout the various phases of the
mission. Since the amount of transmitted data between these two subsystems is
small, a low-speed data link is sufficient, although there is a new trend to incorpo-
rate high-speed standard link such as SpaceWire4 to satisfy increasing demand for
data volume.
The OBC is also responsible for receiving, interpreting, and executing com-
mands from ground operators via the radio receiver. Using low-speed radio
transmitters, the OBC also sends packets of housekeeping data, or telemetry, to the
ground station. The purpose of the housekeeping data is to give the operators on the
ground an overview of the spacecraft health and its general condition.
Some small satellites only have a single low-speed transmitter, so the housekeep-
ing and payload data are combined over the same link. For larger satellites with
payloads capable of producing vast amounts of data, a dedicated high-speed data
link is used to store the data on an onboard storage system. When the satellites pass
over a ground station, the OBC commands the high-speed radio transmitter to
retrieve and transmit the previously stored payload data through another dedicated
high-speed link from the onboard storage system. This approach frees the OBC
from having to process large amounts of data and allows it to devote its internal
4
SpaceWire is a spacecraft communication network based in part on the IEEE 1355 standard of com-
munications. It is coordinated by the European Space Agency (ESA) in collaboration with international
space agencies including NASA, JAXA, and RKA
8
Communication Subsystems for Satellite Design
DOI: http://dx.doi.org/10.5772/intechopen.93010
resources for time critical operations and communicates with the PL and all other
subsystems through the low-speed data links. This would include the requirements
to retrieve information on the health, perform critical interventions as well as to
command these subsystems to perform various actions according to the operational
arrangement of the mission.
1. Telemetry: the collection, processing of health, and status data of all spacecraft
subsystems, and the transmission of these data to the control segment on the
ground. This requires not only a telemetry system on the spacecraft but also a
global network of ground stations around the world, unless the satellite space
network includes intersatellite links that can relay the data to designated satel-
lite and downlink to the appropriate ground station. Figure 8 below illustrates
the processing of telemetry data by the C&DHS. Here the different health
information and status information sent from various subsystems are col-
lected by the telemetry input interface, fed to the C&DHS processor, buffered,
encrypted, and sent down to the ground station.
3. Command and control: the reception and processing of commands for contin-
uous operation of the satellite. Usually a ground system is required, although
advanced spacecraft designs have evolved toward “autonomous operations”
so that many of the control functions can be automated onboard and do not
require ground intervention except under emergency conditions. A typical
command processing scenario is illustrated in Figure 9 where serial command
bit stream from the command receiver is received by the command input inter-
face, where the relevant commands are extracted and sent to the appropriate
subsystems via a serial or parallel interface.
Figure 8.
Telemetry processing by C&DHS.
9
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 9.
Command and control message processing by C&DHS.
Figure 10.
Channelized processor for communications payload.
10
Communication Subsystems for Satellite Design
DOI: http://dx.doi.org/10.5772/intechopen.93010
data processing and the lack of onboard computational power to accomplish these
tasks. An example of onboard PL processing for passive electro-optical (EO) remote
sensing is shown in Figure 11 below, where the reflected light from earth is passing
through a combination of optical lenses and charge coupled device5 (CCD) whose
output is an analog signal that would be conditioned by analog filters before being
digitized, compressed, and sent down via a mission data downlink to the ground
station for processing. There, the data are decompressed, and image is enhanced by
appropriate algorithms and displayed for users.
Typical data volume collected by sensing payload is large, and peak rates can
produce data at much higher speeds than TT&C; thus, a separate downlink for
mission data is needed. Depending on the system, this mission data downlink to a
ground station can either be performed using a dedicated mission direct downlink,
or indirectly via a relay broadband communications satellite. Sensing satellite can be
positioned in GEO, MEO, or LEO orbits, and can have many possible mission data
downlink architectures based on mission requirements. For example, a LEO sensing
satellite can either buffer its mission data until within view of a dedicated ground
station for downlink, or it can forward its mission data to a relay satellite that can
ensure that the mission data can be downlinked to a designated ground station.
Another example of active remote sensing is a synthetic aperture radar (SAR)
mission, where returned radar signals are collected onboard and sent to the ground to
be correlated and form an image of the ground surface. This type of remote sensing
does not heavily depend on sun light and other weather affects. Applications for SAR
include agriculture, geology, geohazards, ice, oil spills, and flood monitoring. Several
emerging applications such as forestry, ship detection, and others are possible [1].
An example of a SAR mission is the NASA-ISRO Synthetic Aperture Radar (NISAR)
[5], which is a collaborative earth-science mission between NASA and the Indian
Figure 11.
Onboard image processing for an EO application.
5
A CCD is an integrated circuit etched onto a silicon surface to form light-sensitive elements repre-
sented by what are called pixels. Photons incident on this surface generate a charge that can be read by
electronics as a voltage and turned into a digital copy of the light patterns falling on the device.
11
Satellite Systems - Design, Modeling, Simulation and Analysis
Space Research Organization (ISRO). The sensing payload features an L-band SAR
instrument and an S-band SAR instrument. The simultaneous dual-frequency radar
system at peak rates will produce data at gigabit-per second speeds, which drives
the data-volume requirements at a minimum of 35 Terabits per day of radar science
data to the ground. This is a direct mission downlink system with three designated
ground stations. The payload communication system uses a 70-cm high-gain antenna
with two synchronized transmitters in a dual-polarization configuration with each
transmitter providing 2.4 Gbps of coded data with an aggregate rate of 4.8 Gbps.
Before going into SDR basics, some of the SDR advantages are [6]:
6
A cognitive radio (CR) is a radio that can be programmed and configured dynamically to use the best
wireless channels in its vicinity to avoid user interference and congestion.
12
Communication Subsystems for Satellite Design
DOI: http://dx.doi.org/10.5772/intechopen.93010
an SDR-CR can “borrow” the spectrum until the owner comes back. This
technique has the potential to dramatically increase efficient use of radio
frequency spectrum.
• Reduced obsolescence: an SDR can be field upgraded to support the latest com-
munications standards. This capability is especially important to radio with
long life cycles such as those in satellite communications.
• Lower cost: a single SDR can be adapted for use in multiple markets and for
multiple applications. For example, a single radio can be sold to cell phone and
automobile manufacturers to significantly reduce cost.
• Cost is the most common argument against SDR. A single key fob is based on
a very inexpensive ASIC7; however SDR is heavily reliant on FPGA,8 which is
much more expensive. This is even more significant for high-volume, low-
margin consumer products.
• The second most common argument against SDR is increased power consump-
tion with increased DSP complexity and higher mixed-signal/RF bandwidth.
Power consumption in an FPGA or GPP for flexible signal processing can easily
be 10 times higher than in ASIC. Also, wideband analog-to-digital converters
(ADCs), digital-to-analog converters (DACs), and RF front-ends consume
more power than their narrowband equivalents.
• Increased time and cost to implement the radio: it can take much more engi-
neering effort to develop software/firmware for multiple waveforms than for
one, especially if it must be compliant with a military standard such as JTRS9.
• Changing specifications and requirements: this usually happens when the SDR
design must support not only a set of baseline waveforms but also anticipate
additional waveforms.
• Limited technical scope: SDR only addresses the physical layer and will require
cooperation from upper layers for throughput improvements.
7
Application-specific integrated circuit (ASIC) is designed for specific purpose.
8
Field Programmable Gate Array (FPGA) is an integrated circuit designed to be configured by a
designer after manufacturing thus much more general purpose.
9
The Joint Tactical Radio System (JTRS) aims to replace existing radios in the U.S. military with a single
SDR to enable new frequencies and waveforms added via software or firmware upload.
13
Satellite Systems - Design, Modeling, Simulation and Analysis
The general definition for a SDR is a radio with some or all its physical layer behav-
ior defined through means of software [7, 8]. SDRs are incredibly valuable devices as
they allow the end user the ability to traverse the RF spectrum at variable sampling
rates. The fundamental qualities that make up an SDR are the flexible specifications
and the ability to transform the analog signal using digital signal processing (DSP).
A radio can be categorically separated into receivers and transmitters. For this
section, the receiver implementation will be considered as it is generally more
interesting and complex. A block diagram of an SDR receiver is shown below in
Figure 12. The following sections will present the anatomy of the SDR that differ-
entiates it from a traditionally designed radio.
4.1.2 RF front-end
The purpose of the RF front-end (RFFE) is to isolate the desired signal received
by the antenna from interference signals. To achieve this, the signal of interest must
be brought down to lower frequency for digital conversion while mitigating the side
effects from filtering during the frequency conversion process. A flexible RFFE
for SDR must be designed so that the frequency and bandwidth are controllable by
software. Depending on the system requirements and the available RF component
specifications, there are several ways to achieve this.
One of the most common RFFE designs for analog radios is the heterodyne
receiver. A heterodyne receiver, shown in Figure 13 below, works by mixing down
the received signal from its carrier frequency to a lower intermediate frequency
(IF). The signal at IF can now be more conveniently filtered, amplified, and pro-
cessed. A super-heterodyne receiver uses a fixed IF that is lower than the carrier
frequency but higher than the signal bandwidth and often uses two stages of down
conversion to reduce the filtering requirements at each stage.
Another popular RF front-end architecture generally used for low-power
applications is called zero-IF. A zero-IF receiver, shown in Figure 14 below, uses a
single mixing stage with the local oscillator (LO) set directly to the desired carrier
frequency to convert directly to baseband in-phase and quadrature signals. Because
mixers tend to have high power consumption and only low-pass filters are required,
the simpler zero-IF provides improved power efficiency over a heterodyne architec-
ture. However, the zero-IF implementation is more susceptible to IQ imbalances of
the in-phase and quadrature oscillators, which will produce anomalies in the signal
constellation. LO leakage may also self-mix through the RF ports creating a large
DC bias. Both issues can be corrected using digital signal processing.
Figure 12.
A block diagram of an SDR.
14
Communication Subsystems for Satellite Design
DOI: http://dx.doi.org/10.5772/intechopen.93010
Figure 13.
Heterodyne receiver.
Figure 14.
Zero-IF receiver.
The two main functions of a digital front-end are sample rate conversion (SRC)
and channelization. Once a signal has become digitally converted, the samples
need to be further primed for digital processing. Operating the ADC at a fixed rate
simplifies its clock generation; however, it may be necessary to convert the sampling
rate to match the sampling rate required to demodulate certain waveforms. Most
wireless signals generally operate with specific symbol or chip rates that are speci-
fied by their respective standard. Depending on the RFFE design and signal type,
channelization may be required to select the channel of interest.
SRC represents a classic sampling theorem problem. Converting sampling rates
can introduce undesirable effects such as aliasing, an effect that causes frequency
components to overlap. SRC can be achieved digitally through the processes of deci-
mation and interpolation. To mitigate aliasing, decimation is performed by using an
anti-aliasing filter followed by subsampling, which is essentially removing samples
at certain intervals. Interpolation is a method of calculating values to add values
in between samples. Channelization works by using digital down conversion, the
process of digitally mixing down a signal to baseband with a numerically controlled
oscillator.
15
Satellite Systems - Design, Modeling, Simulation and Analysis
SDRs have an array of devices to choose from for the required DSP application,
each with their own strengths and weaknesses. An SDR may integrate multiple pro-
cessor types and partition the signal processing chain to optimize each processor.
The following criteria should be considered when evaluating the various processor
types: flexibility, modularity, and performance. The three digital hardware choices
this section will consider are the general-purpose processor (GPP), digital signal
processor (DSP), and the field programmable gate array (FPGA).
A GPP is the typical microprocessor designed to handle a wide variety of generic
tasks that can be found in your everyday personal computer. They are generally
designed to have large instruction sets and highly capable of implementing and
performing complex arithmetic tasks such as modulation/demodulation, filtering,
fixed/floating point math, and encoding/decoding. Some commonly used GPP
architectures are x86/64 and Advanced RISC Machine (ARM). The advantage of
using a GPP is the wide availability, flexibility, and ease of programmability. Several
GPP-based SDRs, such as Universal Software Radio Peripheral (USRP) and the
LimeSDR, operate by digitizing the baseband signal and performing the required
digital signal processing on computers. These types of SDRs are popular among uni-
versity researchers and hobbyists due to the relative ease of obtaining and develop-
ing their applications.
Because the GPP was designed with such a broad focus, latency, speed, and
power efficiency may be a limiting factor depending on the application. Many wire-
less communication standards have strict real-time and large processing bandwidth
requirements that most modern CPUs cannot meet due to processor architecture
and operating system design. .
A DSP is a microprocessor optimized for digital signal processing applications
with the ability to be programmed with high-level languages. Although a GPP can
contain much of the same functionality, the DSP performs the same digital signal
processing operations more quickly and efficiently due to its reduced instruction set
computer (RISC) architecture and parallel processing. The reduced instruction set
limits the essentials but contains optimizations for common DSP operations such as
multiply accumulate (MAC), filtering, matrix operations, and fast Fourier transform
(FFT). DSPs are commonly sold in two variants: optimized for power efficiency and
optimized for performance; and are used in applications such as base stations and
edge devices. Power consumption is also minimized by reducing the silicon footprint
that would be in GPPs sophisticated cache and peripheral subsystems.
Although DSPs have been commonly deployed in the past decades, they serve as a
middle ground between GPPs and FPGAs with regard to flexibility, performance and
efficiency. Field-programmable gate array (FPGA) offers more parallelism, higher
data rates, and better power efficiency than DSP, but is not well suited for control
applications, such as implementing the network/protocol stack. This is due to the
limited amount of memory in FPGA and for this reason it is often paired with GPP.
A FPGA is an array of programmable hardware logic blocks, such as general
logic, memory, and multiplier blocks, that are wired together via a reconfigu-
rable interconnect to generate an integrated circuit for several designs with
the ability to quickly switch between configurations. FPGA configurations are
programmed using hardware description language (HDL), which is also used for
ASIC. Because a FPGA functionality is defined at the hardware level and can be
implemented using parallelism, it can perform DSP algorithms at much higher
rates than DSPs and GPPs. FPGA consumes more power and requires more space
than ASICs but provides more programmability and flexibility than ASIC. A big
consideration for using FPGAs for SDR is the domain knowledge requirement
16
Communication Subsystems for Satellite Design
DOI: http://dx.doi.org/10.5772/intechopen.93010
For space applications, SDR has unique challenges such as extreme radiation and
temperature environment, autonomous operational requirements, limitations on
size, weight and power (SWAP), and the need for reduced development time and
increased reliability in agile prototyping. In this section, recent applications of soft-
ware defined radio to satellite, as well as the current status of radiation-hardened
SDR components, are presented.
2. Signal processing module (SPM) where signal processing is used to handle the
transformation of digital signals into data packets. Its components include ASICs,
FPGAs, DSPs, memory, and connection fabric/bus (e.g., PCI, flex-fabric).
10
In a HID, the hardware architecture requirements are written so that the hardware provider defines
the functional modules of the system and publishes the functions and interfaces for each module and for
the entire STRS platform [Wikipedia].
17
Satellite Systems - Design, Modeling, Simulation and Analysis
Figure 15.
NASA STRS’ three main hardware functionalities.
11
Hardware abstraction layer (HAL) is a layer of programming that allows a computer operating system
to interact with a hardware device at a general or abstract level rather than at a detailed hardware level
[Wikipedia].
12
The Portable Operating System Interface (POSIX) is a family of standards specified by the IEEE
Computer Society for maintaining compatibility between operating systems. POSIX defines the applica-
tion programming interface (API) for software compatibility with variants of Unix and other operating
systems [Wikipedia].
13
CubeSats are a class of Small Satellites (SmallSats) weighing between 1 kg and 10 kgs that use stan-
dard size and form factor of 1 U (one unit) of 10 cm x 10 cm x 10 cm. A standard 3 U CubeSat (10 cm
x 10 cm x 34 cm, 5 and 6 kg) has been demonstrated to support real mission, and larger CubeSats (6 U,
12 U, and 24 U) are developed.
18
Communication Subsystems for Satellite Design
DOI: http://dx.doi.org/10.5772/intechopen.93010
Figure 16.
STRS software architecture layers.
19
Satellite Systems - Design, Modeling, Simulation and Analysis
custom RF front-end along with VHF Ku-band antennas. This Pathfinder mission14
was to enable onboard reception and geolocation of different types of terrestrial RF
signals using signal processing technique to combine received data from all three
payloads15.
One commercial application of this mission is the detection and geolocation
of a maritime vessel’s automatic identification system (AIS), which broadcasts the
locations generated by GPS-enabled receiver. The locations generated by AIS can be
disabled or spoofed, therefore not reliable. Another application would be to allow
regulators, telecommunications companies, and broadcasters to globally monitor
spectrum usage and identify areas of interference. The system can also be used to
help large area search and rescue operations by quickly locating activated emer-
gency beacons.
The SDR developed for the Pathfinder payload consists of an embedded pro-
cessor system and three baseband processors. The baseband processor was built
around the Analog Devices 9361 (AD9361) System on Chip (SoC) product, which
is a highly integrated RF transceiver that combines high-speed ADCs and DACs,
RF amplifiers, filtering, switching plus more. The HE360 payload supported up to
three receiver channels (one AD9361 per channel) that can be simultaneously pro-
cessed on separate frequencies. In addition, the signal processing subsystem takes
advantage of open-source software and firmware code to allow system development
to proceed without knowing the final space hardware. GNURadio16 was selected
for being a free and open-source toolkit for SDR and widely used in small space
projects for ground software processing.
14
A Pathfinder mission is usually a demonstration to prove that a system can successfully achieve a
specific objective before a full mission can deploy.
15
By comparing time-of-arrival (TOA) and frequency-of-arrival (FOA) measurements between pairs of
receivers, the position of a signal can be computed.
16
GNU is a recursive acronym for “GNUs Not Unix!“ chosen because GNU’s design is Unix-like but
differs from Unix by being free software and containing no Unix code.
17
Design trade-offs choosing more better radiation tolerance integrated circuit technology, low-power
Schottky vs. Emitter Couple Logic (ECL) vs. bipolar over CMOS. Rad-hard chips are often manufactured
on insulating substrates (silicon on insulator, silicon on sapphire, silicon carbide, gallium nitride) instead
of the usual semiconductor wafers. Shielding the package against radioactivity to reduce exposure of the
bare device.
18
TID causes slow gradual degradation of the device’s performance and is measured in rads.
19
SEEs are caused by a single energy particle and can be (a) non-destructive Single Event Upset (SEU)
causing transient pulses in logic or support circuitry or as bitflips in memory cells or registers, (b) poten-
tially destructive Single Event Latch-up (SEL) that may be cleared by a power reset, and (c) destructive
Single Event Burnout (SEB) or Single Event Gate Rupture (SEGR), which is irreversible.
20
Communication Subsystems for Satellite Design
DOI: http://dx.doi.org/10.5772/intechopen.93010
Per US DoD MIL-PRF-38535 J standard [18], an ideal integrated circuit for space
applications is the qualified manufacturing line20 (QML) Class V with radiation
hardness assurance21 (RHA) level identified in the part specification. From the
perspective of payload designer and developer, only Class V is space quality and
should be the main factor for selecting SDR hardware components.
The FPGA is perhaps the most important component of an SDR and has a long
history for manufactured QML class V parts where rad-hard Xilinx and Actel
(now Microsemi) FPGAs were studied [19]. Currently, Xilinx is the major player
for space-qualified QML level V products used in actual payloads with many more
devices under development. The rad-hard DSP products also follow the QML
process, with Texas Instrument (TI) currently taking the lead for in-flight payloads
with many offerings in space-qualified RF components in addition to DSP. Similarly,
space-qualified GPP follows the same QML path as FPGA and DSP, and the current
on-flight rad-hard GPPs based on the following architecture are [20].
In the first section of this chapter, an overview of the satellite bus and payload
subsystems are presented for command and data handling subsystem (C&DHS);
communications subsystem (CS); electrical power subsystem (EPS); propulsion
subsystem (PS); thermal control subsystem (TCS); attitude control subsystem
(ACS) also known as guidance, navigation and control (GNC) subsystem; and
structures and mechanics subsystem (S&MS). A significant portion is spent on
describing the C&DHS and CS with much details on how they are related to other
satellite subsystems for continuous operation.
There are distinctive functional separations between the satellite bus and pay-
load that are discussed at a high level with some examples given; however, there are
currently no existing standard on their interfaces due to legacy satellite design and
development. Examples were given for mission-specific sensing and communica-
tions payloads, showing that pretty much all mission payloads are very customized
in design in legacy systems.
The second section of this chapter covers software defined radio (SDR) as a new
technology with an overview and how SDR is being applied to satellite design and
development in both space and ground segments. There has been a NASA standard
for SDR that has been used for traditional and large satellites and shown to have
some advantages over non-SDR approach.
However, recent rapid developments of Small Satellites (SmallSats), which
CubeSat is a subset of, have resulted in an explosion of SDR applications to build
20
For QML microcircuits, the manufacturer is required to develop a program plan that meets or exceeds
the performance detailed in these appendices of MIL-PRF-38535 J standard.
21
RHA is quantified in terms of the radiation level in Total Ionizing Dose or TID.
21
Satellite Systems - Design, Modeling, Simulation and Analysis
Pathfinder missions that can lead to successful follow-on projects. There remains to
be a standard to be defined for SDR for this CubeSat application. Regardless, SDR is
providing a path forward to a common framework that may enable a more generic
building block for a future concept called Software Defined Satellite that will
change missions based on a software upload.
Since SDR is becoming an important part of a satellite, radiation hardening
of the relevant SDR components is described in some detail. The area is evolving
slowly despite fast changing technology due to the additional design and manufac-
turing steps taken to ensure minimum effects of radiation on microelectronics. The
selection of the appropriate rad-hard FPGA, DSP, and GPP components should be
an important factor in design trade-offs when SDR is being considered for future
missions.
Acknowledgements
The first author, Dr. Hung H. Nguyen, would like to express bountiful apprecia-
tion for his wife, Thuy Le Nguyen, for her constant support during this effort.
Author details
© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms
of the Creative Commons Attribution License (http://creativecommons.org/licenses/
by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,
provided the original work is properly cited.
22
Communication Subsystems for Satellite Design
DOI: http://dx.doi.org/10.5772/intechopen.93010
References
23
Satellite Systems - Design, Modeling, Simulation and Analysis
24