1. Introduction
Microwave solar emissions at a wavelength of 10.7 cm, or 2800 MHz, are an excellent indicator of solar activity. Measurements of the solar flux density at this wavelength are known as the 10.7 cm solar radio flux, or the F10.7 index. This index is widely used in thermosphere and ionosphere models because these microwave emissions correlate well with other solar emissions such as extreme ultraviolet (EUV), ultraviolet (UV), and X-rays, which originate from similar regions of the solar atmosphere, namely the upper chromosphere and base of the corona. The correlation with EUV is of particular interest as solar EUV flux is the main energy source in the terrestrial atmosphere producing ionization, disassociation, excitation, and heating [
1,
2,
3]. Variability in the ionization levels in the ionosphere affects the propagation of radio waves, which in turn can cause errors in terrestrial and satellite communication and navigation signals [
4]. The heating of the neutral atmosphere in the thermosphere causes drag forces to increase on Low Earth Orbit (LEO) objects, disrupting the orbital paths of satellites [
5,
6,
7]. Consequently, the monitoring of the EUV flux is important for thermosphere and ionosphere studies and space weather predictions.
The most notable benefit of microwave solar emissions is that they can reach the surface of the Earth, while EUV wavelengths are absorbed by the atmosphere, as shown in
Figure 1. The atmospheric window for 10 cm microwaves has been critical for F10.7 observations historically. It has also made the measurement of F10.7 solar flux exceptionally robust because it can be reliably and accurately recorded from the ground in all weather conditions with few gaps or calibration issues [
8]. F10.7 measurements were therefore used as a proxy for EUV measurements, which were unable to be directly measured until the dawn of space-based platforms in the late 20th century.
While EUV flux can now be monitored from space platforms, the highly reliable and repeatable calibration of the F10.7 solar radio flux measurement, along with its continuous historical record, has maintained its standing as one of the most widely used space weather indices [
1,
10,
11]. It is a key input parameter to atmospheric models and provides climatology of solar activity over six solar cycles [
6,
8,
12]. This duration is only surpassed by sunspot records [
13].
Observations of F10.7 solar flux were initially recorded by radio telescopes based near Ottawa, Canada from 14 February 1947 to 31 May 1991 [
1]. As of 1 June 1991, the Dominion Radio Astrophysical Observatory (DRAO) in Penticton, Canada is the only ground-based observatory to collect daily measurements. Two radio telescopes onsite collect three one-hour-long solar measurements per day during daylight hours [
1].
While measurements have been historically reliable, the inherent vulnerability of a single-station measurement represents a potential point of failure in global ionospheric and thermospheric modeling [
14]. Furthermore, F10.7 is currently under-sampled, with the three daylight measurements from a single site constrained further by the effects of local terrain on usable solar elevation angles. The solar flux output is known to vary over a range of timescales, from intense bursts at millisecond to minute scales, to smaller variations over hours to days, following an overarching 11-year solar cycle [
1]. The received 10.7 cm flux of interest will vary over the course of a day as active regions rotate across the surface of the Sun. The magnitude of these variations is often more than a solar flux unit and increases with the level of solar activity. Discrete measurements obtained from a single geographic location over a three-hour total observation period would not necessarily capture these variations and therefore not provide certainty of the average solar flux value.
The problem of under-sampling has been identified by the current observation program and a mitigation strategy of multiple ground-based stations in different longitudes has been suggested [
1]. The investment required to implement this strategy could be substantial; however, if a strong correlation between 2800 MHz and another solar density flux frequency, 2695 MHz, can be established, the existing monitoring system of the latter frequency could be leveraged, with observations scaled to provide improved temporal coverage and thus more accurate F10.7 solar flux approximations. The United States Radio Solar Telescope Network (RSTN) is a distributed network of ground-based radio telescopes dedicated to measuring solar flux densities across a range of microwave frequencies from 245 to 15,400 MHz, including 2695 MHz. Four stations constitute the network, located in Massachusetts, Hawaii, Southern Italy, and Western Australia. Meaningful comparisons between RSTN and DRAO measurements are not currently possible, due to differing instrument calibration standards and operational regimes between the two programs [
15,
16]. As such, 2695 MHz data are currently unable to be scaled and directly substituted for F10.7 index applications. This paper introduces the mission concept and design of a space-based solution that would increase the geographical and temporal coverage of F10.7 measurements and provide redundancy to the existing solar monitoring system. A CubeSat platform is proposed by the Royal Military College of Canada (RMC) for extra-terrestrial solar flux observations of 2800 MHz and 2695 MHz frequencies. Satellite observations of 2800 MHz will increase the number and temporal coverage of daily F10.7 measurements available. The inclusion of solar flux emission observations at 2695 MHz into the mission of RMCSat allows for accurate comparisons between 2800 MHz and 2695 MHz without calibration issues, as the same instrument is collecting both frequencies at the same time. This will enable the accurate determination of the relationship between 2695 MHz and 2800 MHz flux values, act as an independent calibration source for RSTN terrestrial measurements during overlapping operations windows, and facilitate the use of RSTN flux values to critically supplement the existing F10.7 observation dataset. RMCSat is a technology demonstrator intended to validate the measurement of F10.7 from space that will allow the use of existing RSTN stations as new nodes in a continuous-observation global ground network. The use of satellite technology for this endeavor is a cost-effective approach compared to a worldwide network of perfectly calibrated solar observation stations.
RMCSat is the first CubeSat mission undertaken by RMC and provides the additional benefit of being a training environment for personnel in space mission analysis and design, satellite assembly, integration and testing, and satellite operations.
An overview of the RMCSat mission is presented here, describing the physics of F10.7 solar radio flux in
Section 2, and the mission objective, the performance requirements, and the baseline mission concepts in
Section 3. The 3U CubeSat satellite architecture is then outlined in
Section 4 along with system budgets. A discussion of the future design considerations in
Section 5 precedes the conclusion of the paper.
2. F10.7 Index
The measurement of solar flux is based on the fundamentals of blackbody radiation. According to the Stefan-Boltzmann law, a perfect blackbody will emit energy based on its surface area
and its absolute temperature
in Kelvin. The total isotropic power
emitted across all frequencies is given by,
where
is the Stefan-Boltzmann constant. The spectral radiance of an object of a given temperature does vary with frequency. When observed at a particular frequency, the spectral radiance
of emissions from an object follows Planck’s Law,
where
is the wavelength of emission,
Js, is the Planck constant and
is Boltzmann constant. In the case of microwave frequencies at solar temperatures, where
, the Rayleigh-Jeans approximation can be used to simplify Equation (1) for solar radiation analysis where,
Spectral radiance is a characteristic of the emitting source, the Sun in this instance, while intensity is what is observed by a receiver, the antenna. The specific intensity
is given as
If the specific intensity of the radiation source is integrated over the total solid angle subtended by the emitting source
, as measured at the detector, the result is the spectral flux density
,
By using the small-angle approximation for the cosmic radio sources,
Given Equation (4), the spectral flux density becomes,
where
is the power per unit area per unit bandwidth due to the source measured at wavelength
. This is the measurement of solar flux density commonly used for solar activity analysis and is the basis for F10.7 measurements which focus on an average wavelength of 10.7 cm.
The F10.7 index is the total power flux density,
, from all sources on the solar disk at 10.7 cm within a bandwidth of 100 MHz centered on 2800 MHz, averaged over an hour [
1]. It is comprised of three components: a constant one from quiet Sun background emissions; a slowly varying component (S-component) that originates in active regions and; and a highly variable one resulting from radio burst emissions [
1,
10]. The recorded flux density at this wavelength varies daily. It is expressed in solar flux units (sfu) or 10
−22 W/mHz and can vary from approximately 64 sfu to greater than 300 sfu over the solar cycle [
1].
Flux measurements taken from the ground-based DRAO are obtained using two 1.8 m diameter paraboloid radio telescopes that run in parallel. They are referred to as Flux Monitor 1 and Flux Monitor 2; the latter acts as the primary instrument while the former is a backup [
1]. During the northern hemisphere summer months, measurements are taken each day at 1700, 2000, and 2300 UTC, and in winter at 1800, 2000, 2200 UTC due to the reduced daylight duration. Readings are taken both on-Sun and off-Sun, to obtain accurate background cold-sky values. The solar data are then corrected for atmospheric absorption, background sky temperature, and antenna gain [
17]. Precise calibration standards are maintained through annual calibration and validation of both flux monitors against an external reference source, a large pyramidal horn antenna maintained onsite at the DRAO.
When detecting solar flux density from a space-based platform in LEO, the sensor antenna receives signals from multiple external sources: the desired solar radio flux, microwave thermal emission from the Earth, electronic radio-frequency interference (RFI) from Earth-based and airborne radars and transmitters, and background emissions from Cosmic Microwave Background radiation and other galactic sources [
18]. The strength of these different signals will vary constantly, dependent on the position and pointing direction of the satellite, complicating solar flux measurement. The signal detected by the antenna also includes noise generated from the receiver system itself. The solar signal, and small variations in this signal, must be distinguishable among all of the noise.
Knowledge of an antenna gain pattern is useful to understand the issue of background noise detection. In
Figure 2, an antenna directive gain elevation pattern is shown for a simple helical antenna, like that used by RMCSat. The solar signal is detected by the main beam while attempts are made to minimize the background signals detected by main, minor, and background lobes.
For changes in the solar flux density to be detected, the antenna must have sufficient sensitivity and minimal noise from other sources. Each source contribution to the antenna temperature
can be determined by,
where
is the brightness temperature of the source,
is source solid angle and the gain
is applicable to the region of the antenna pattern that the signal is localized within.
In the first instance, a certain detection threshold must be achieved for the solar emissions at 10.7 cm. Based on existing F10.7 data, a minimum value of 70 sfu for F10.7 emissions over the solar cycle results in an isotropic antenna temperature of 0.23 K at 10.7 cm. If a flux measurement is to have an accuracy of 1%, as is required for ground-based DRAO F10.7 determinations, then the required antenna sensitivity is 2.3 × 10
−3 K. This value can be increased if a more directional antenna beam is used, since a higher main beam gain (
increases the tolerance of the system to noise signal fluctuations [
1]. The antenna sensitivity for F10.7 = 70 sfu becomes,
The effect of terrestrial thermal emissions is more complicated as the brightness temperature of Earth changes as a satellite passes over different landscapes. A maximum average brightness temperature of 300 K (land) and a minimum value of 100 K (ocean) have been estimated [
18]. The 200 K fluctuation between the minimum and maximum terrestrial brightness temperatures results in an antenna temperature difference of 0.61 K. This is two orders of magnitude greater than the required antenna sensitivity and nearly three times the value for 70 sfu. Even if a higher gain of 20 dBi is achieved by RMCSat, the detection threshold only reaches 0.33 K. Thus, fluctuations in thermal emissions from Earth are sufficient to obscure the primary solar signal detected by the payload and must be accounted for in some manner.
Similarly, electronic radio-frequency interference (RFI) from Earth and airborne communication systems presents a significant issue. Microwave signals of sufficient power are detected by an orbiting radiometer due to atmospheric transparency [
18]. The frequency band of 2700 to 2900 MHz hosts global aeronautical navigation systems, radar weather monitoring, and maritime and land-based radiolocation services [
19,
20,
21,
22]. Transmitter peak output powers range from 22 kW up to 1.4 MW [
19]. Using a bandwidth of 100 MHz, equivalent to the DRAO F10.7 measurements, the equivalent antenna temperature contribution from a radar system such as the Next Generation Weather Radar (NEXRAD) is 6.1 × 10
14 K [
18]. Such a signal would saturate the receiver. Fortunately, the low duty cycle and short duration of signal pulses mean that it is not expected to cause significant issues for the overall utility of measurements. The 2695 MHz band is a protected frequency and is not expected to have radar inference. Consequently, it can be used for comparison to identify RFI within the F10.7 data.
Cosmic microwave background radiation is a homogeneous field that has near-perfect blackbody radiation at a temperature of 2.748 ± 0.016 K. The antenna temperature contribution is 2.09 K, but being constant, it is easily taken into account. The galactic center also produces radiation and an antenna temperature contribution of 2.86 × 10−3 K. For a directional antenna like that of RMCSat, issues in solar flux measurement accuracy will occur where the galactic center is within the main beam and solar activity levels are low. Given the fixed position of the galactic center, this contribution can be anticipated and accounted for.
System contributions are generated from thermal noise resulting from resistance in radiometer components and transmission line losses between the antenna and receiver [
23,
24]. The magnitude of fluctuations in system noise limits the minimum detectable antenna temperature change. Since fluctuations are random and have shorter timescales than the primary signal, they can be reduced towards a zero mean by averaging the total signal over a selected integration time [
23,
24].
The analysis of the radio environment of RMCSat indicates that the thermal radiation from Earth is the greatest concern, followed by RFI. Proposed mitigation options include a fixed orientation of the satellite towards Earth, a secondary antenna, and measurements are taken across two close but separated channels. The fixed orientation places the Earth in the lowest gain region (sidelobe and backlobe) of the primary antenna. Fluctuations that are detected can then be attributed to changes in the brightness temperature of the surface below. A secondary antenna directed at the Earth offers a means to capture the surface signal and use it for comparison with the primary antenna signal. Changes in Earth emissions will be apparent in the secondary antenna signal. Post-processing can be used to remove the Earth component from the primary signal. A wide beam secondary antenna is necessary to replicate the area visible to the main payload antenna and allow for an accurate comparative signal. Preliminary modeling indicates that sufficient accuracy can be achieved using this method [
18].
In place of a single 2800 MHz channel, the use of two close but separate channels is selected to account for electronic RFI signals and the saturation of the receiver when these are present. Measurements are taken at 2785 MHz and 2815 MHz (designated F10.7- and F10.7+ respectively) with a narrower bandwidth of 20 MHz for each to minimize the chance of an RFI signal present in both channels. The dual channels provide data continuity if one channel is affected, and provide a method to distinguish between RFI signals and solar burst emissions. The inclusion of 2695 MHz monitoring in the mission provides additional control, as interference found in all three channels indicates a solar burst as the most likely cause. For a meaningful F10.7 solar flux record, the measurements obtained by the F10.7− and F10.7+ channels can be combined in post-processing to obtain the equivalent average 2800 MHz value.
Finally, the calibration of a radiometric sensor is critical to system accuracy. Importantly, the space-based solar monitoring regime will allow up to three RMCSat observations per day that are concurrent with the DRAO solar observations. This will allow a highly accurate reference for calibration and validation of RMCSat values.
4. Satellite Architecture and Subsystems
RMCSat is a 3U CubeSat (three stacked 10 × 10 × 10 cm cubes) with deployable solar panels and antennas as modeled in
Figure 9. The architecture is constructed around the zenith and nadir payload antennas and the need to maximize solar exposure in a noon-midnight orbit.
The payload consists of two antennas and an S-band receiver with software-defined radio (SDR) architecture. A deployable helical antenna acts as the primary antenna recording the microwave activity of the Sun on the desired frequencies while the secondary wideband patch antenna is used to detect background radio noise from RFI and thermal Earth emissions on the same frequencies.
The ADCS, which stabilizes the spacecraft against external perturbations, consists of three reaction wheels, a 3-axis magnetometer, two IR Earth Horizon Sensors, three electromagnets, and an astronomy data and computing services computer. It maintains the nadir orientation of the secondary antenna throughout the orbit and points the helical antenna at the Sun with required pointing accuracy during solar observations. A global navigation satellite system (GNSS) antenna and receiver offers <10 m positioning accuracy by receiving GPS L1 and GALILEO E1 signals. GPS will be sampled once per second and used for post-processing of F10.7 data and generation of a global microwave emissions model.
Payload and telemetry, tracking and command (TT&C) data are downlinked using a nadir S-band patch antenna identical to the payload secondary patch antenna. The bandwidth of the payload antenna cannot be extended to reach the 2250 MHz frequency of the ground station unless customized at a notable expense. A dual UHF/VHF dipole antenna with a complimentary UHF/VHF transceiver is used for TT&C uplink via the VHF frequencies. UHF frequencies could be used to downlink part of the solar observation data in the event of failures within the S-band system.
The satellite subsystems are supported by an onboard computer and electrical power system (EPS). Power is generated by gallium arsenide (GaA) triple-junction solar cells laid out on three 1U body-mounted solar panels and four double-sided 3U deployable solar panels. A 22.4 W⋅h Li-Ion battery pack provides power during an eclipse. The battery pack and a heater are integrated into the EPS. Satellite thermal control is achieved passively with aluminized kapton surface coatings.
RMCSat does not require a propulsion system as the satellite does not need to maintain its original orbit nor carry out orbital maneuvers during the mission. External perturbations experienced by RMCSat over the course of its mission are not anticipated to be large enough to disrupt the orbit to the point where the mission is impacted.
The estimated mass budget for the Cubesat is given in
Table 5. Margins are assigned to components depending on their flight heritage (5%) and technological readiness (10%), customization (20%), and in-house manufacturing (25%). The total estimated mass with margins sits comfortably within the 4kg limit for a 3U Cubesat.
4.1. Payload and Data Generation
The primary antenna is a deployable helical wire antenna with 40 turns (63 cm long), a diameter of 3.2 cm, a pitch of 1.7 cm, and a thickness of 1 mm. A deployable antenna of this length can realistically be deployed from a 1 U volume of the satellite. Such an antenna has been analyzed using
MATLAB Antenna Designer and offers gain directivity of 17 dBi at 2800 MHz, and 16.7 dBi at 2695 MHz [
18]. A ground-based field experiment using a small parabolic antenna successfully demonstrated that an antenna with a gain of 14 dBi is sensitive enough to detect small variations in the F10.7 solar flux to less than 1% of the peak value [
18], meeting the DRAO standard for system accuracy [
1].
Figure 10 presents the directive gain pattern for the 40-turn helical antenna at 2800 MHz and 2695 MHz.
Possible orientations of the antenna relative to the satellite body were investigated with regards to the impact of each orientation on the gain differential (
) of the antenna, defined as the difference between the main beam gain and the average backlobe gain. A number of antenna orientations were tested with respect to the satellite body axes and the optimal orientation is shown in
Figure 9 [
18]. In addition to optimizing
, this orientation increases nadir surface area on the satellite for the secondary and communications patch antennas and creates a smaller cross-sectional area for drag forces enabling longer orbital maintenance.
A microstrip patch antenna array was investigated as an alternate primary payload antenna but was found to be unsuitable for this radiometry application due to suboptimal array spacing options constrained by CubeSat geometry and a reduced solar power generation area. However, the wider beam of an individual patch antenna is ideal for the secondary payload antenna to record the terrestrial thermal emissions. A modified Commercial-Off-The-Shelf (COTS) S-band patch antenna shifts the two standard commercial frequency bands offered to match the desired 2695 MHz and 2800 MHz [
27].
Operationally, the payload takes measurements on three channels, 2695 MHz, 2785 MHz, and 2815 MHz. For each channel, a corresponding secondary antenna measurement is required; therefore the six channels are monitored (
Table 6).
The output from the antenna is a direct current voltage that is directly proportional to the antenna temperature, and hence the power of the signal received by the radiometer. If background noise flux and any internally generated system noise remain constant, the increase in voltage measured when an antenna points at the Sun will be directly proportional to the microwave flux density of the Sun.
A multi-channel SDR gives the mission the flexibility to adjust properties like observation frequencies, bandwidth, and modulation from the ground. The generic processor hardware in conjunction with programmable software enables digital replication of electronic radio receiver circuits. If there is too much interference from RFI, the bandwidths or modulations can be altered.
The minimum observation time is calculated using,
where
P is the orbit period and
in this case is the antenna beamwidth between first nulls (BWFN). Since 2695 MHz and 2800 MHz channels are monitored concurrently, the larger BWFN determines the observation time required to ensure that the nulls of both frequency patterns are captured, enabling accurate determination of the peak flux values.
Table 7 shows the minimum observation time based on BWFN.
The payload has a sample rate of 1000/s, which is averaged and logged at one-second intervals. Allowing for an additional 40% margin for the minimum observation time of 752 s and 50 bytes per line of data,
Table 8 shows the observational data generated by the payload. The TT&C downlink data containing diagnostic information is estimated to be 1 Mb as is the uplink command data. GPS sampling at 1 Hz is estimated at 10 Mb. Thus the total amount of data that needs to be downlinked each day is 27.9 Mb. If only the solar observations and TT&C information are downlinked, the data are 5.75 Mb.
4.2. Communications and Link Budget
The communications plan relies on the use of a single ground station at RMC in Kingston, Ontario, Canada, and a store and forward mode of data collection. The ground station consists of a 3.0 m parabolic S-band antenna with a high rate S-band receiver and a Yagi antenna and associated transceiver for UHF and VHF signals. Other locations with similar S-band downlink and UHF/VHF capabilities are suitable if required as the geographical location of the ground station is not critical to the mission.
The satellite passes over the RMC ground station four to six times each day, with a mean access time of 573 s as determined using the Systems Tool Kit software. The operations plan reserves a 25% error margin given 573 s is the mean access value, thus 430 s is used for rate calculations. A minimum downlink data rate of 62.6 kbps is required if the downlink is completed during one pass. This rate sits well within the S-band system capabilities. The maximum VHF uplink rate of the ground station is 1.2 kbps. Consequently, the estimated 1 Mb of TT&C data is uplinked over two passes resulting in a rate of 1.16 kbps.
If the UHF system is required due to S-band failure, only TT&C data and solar observations can be downlinked (5.75 Mb per day). To achieve an adequate link margin of 3.0 dB for the UHF system (
Table 9), the maximum data rate is limited to 3.5 kbps. Four passes are required to downlink the data using the UHF system at a rate of 3.3 kbps.
Link margin calculations are presented in
Table 9. The link margins of the S-band downlink, VHF uplink, and UHF downlink are calculated as 13.96 dB, 17.37 dB, and 3.06 dB, respectively.
The storage capacity of the onboard computer is 4 GB, with the option of up to 32 GB. This offers 148 days, or up to 1188 days, of storage before payload data must be overwritten. Additional storage exists within the S-band transmitter, up to 32 GB. If the UHF system is required for downlink, successive passes can be used to downlink more observational data once the daily data are received.
4.3. Power Budget
There are four different power modes in which the satellite operates: Solar Observations, Solar Cruise, Eclipse Cruise, and Communication.
Table 10 shows the power budget for the different modes. The peak power usage occurs during Communication, consuming up to 13.7 W.
Communication occurs during one daylight passes of the ground station per day. As shown in
Table 11, the average orbital power necessary is 9.4 W, for which the solar panels are designed to meet. Additional power can be drawn from the batteries during eclipse and peak power operations.
Gallium Arsenide (GaA) triple-junction solar cells offer a production efficiency of 30% and performance degradation of 0.5% per year [
28]. Based on these solar cells, the design orbital parameters, and an average solar incident angle of 45°, the solar array power required is 18.0 W, or an array area of 864 cm
2. Using a cell size of 26 cm
2 and two cells per 1U, 17 1U panels would be required. While the deployable arrays and body-mounted cells provide 28U of the area and the capacity to generate 30.3 W, only 18U of cells is exposed to the Sun at any one time. The 18.0 W design provides a 20% power margin which allows for poor solar angles such as those at the terrestrial poles and for any cell failures.
The proposed COTS EPS offers integrated batteries that have a capacity of 22.4 W⋅h and a battery heater for eclipse conditions. This size is more than sufficient to meet the required battery capacity of 12.4 Wh for RMCSat in full eclipse (7.6 W).
4.4. Thermal Control
The thermal control design of RMCSat seeks to keep the satellite equipment within temperature limits with a 5 °C buffer. The temperature ranges of the satellite components are given in
Table 12. The thermal operating range of the RMCSat is −15 °C to 55 °C, limited by the ADCS and UHF communication systems. The integrated battery heater within the electrical power system allows a greater temperature range for the batteries than is typically found.
A first-order analysis of the thermal environment and temperature indicates that thermal control does not pose an issue for the mission. The sources of external heat considered include the direct solar energy of the Sun, solar energy that is reflected off the Earth, thermal energy radiated from Earth, and atmospheric drag heating [
28]. The internal heat sources are derived from heating caused by spacecraft components, which are taken as the lowest and highest power levels of the satellite operations. The extreme temperature range derived from the analysis is −6.8 °C to 64.8 °C, which sits above the maximum temperature of some components. However, if a quarter millimeter layer of aluminized kapton is used, it lowers the temperatures to −7.8 to 50.7 °C, within the thermal operating range of RMCSat.