Comparison of terrestrial and Martian TEC at dawn and dusk during solstices

This paper used the similarities between the ionospheres on Mars and Earth, the most similar of the terrestrial planets, to examine the relative importance of photochemical and transport processes at dawn and dusk. The amount of plasma present in the ionosphere, as measured by the total electron content (TEC), was examined at different locations for both solstice seasons over a solar cycle. Using the rate of change of TEC as a function of solar zenith angle made it possible to compare the plasma production via photoionisation and loss via recombination in the main layer of each planetary ionosphere despite the extreme differences in the total quantity of plasma. This study finds that, at least to first order, the dawn and dusk TEC slopes at Mars are symmetric. This symmetry is interpreted as an indicator of photochemical equilibrium. Deviations from photochemical equilibrium in different geographic and aerographic regions were used to explore the underlying processes responsible for plasma transport. Seasonal and solar cycle variations were also examined at dusk. These variations found that differing interactions with solar forcing mechanisms resulted in a Martian ionosphere with regions that showed evidence of significant transport processes at solar maximum, while at Earth transport processes were most important at solar minimum. In general, the photochemical processes in both ionospheres behave similarly when no magnetic field is considered. The presence or absence of a magnetic field shape the production via photoionisation and loss via recombination processes in both ionospheres, especially when considering plasma transport. This study has notable implications for comparative aeronomy, as a good understanding of how the ionosphere of magnetised and un-magnetised bodies compares is important for characterising planetary environments and atmospheric evolution over long time scales.


Introduction
Ionospheres are the charged regions that exist within the upper atmospheres of planets, moons, and minor bodies of the Solar and other Stellar Systems. Each Solar System body that has sufficient neutral atmosphere, especially within the inner part of the Solar System, has an ionosphere that is mainly formed by the solar ionising radiation at extreme-ultraviolet (EUV) and X-ray wavelengths (e.g., Witasse et al. 2008;Mendillo et al. 2016). The Earth and Mars are often considered to be the most similar of the terrestrial planets, with similar axial tilts (similar seasons, although longer at Mars), similar length of the day, and both having water at (or near) its surface, albeit on different scales and with different distributions . These similarities invite comparisons, though care must be taken to account for the major differences, which include a strong, intrinsic magnetic field at Earth, a longer year at Mars, and the greater eccentricity of the Martian orbit.
The ionospheres of Earth and Mars also share many similarities, such as a purely photochemical ionosphere formed mainly by O + and O + 2 , and similar electron densities ( ∼10 11 m −3 ) and altitudes ( ∼135 km) of some of their layers (for example, the main ionospheric layer at Mars and the secondary layer at Earth that is called the E-layer) (e.g., Hanson et al. 1977;Mendillo et al. 2011;Benna et al. 2015). However, one of the main differences Open Access *Correspondence: angeline.burrell@nrl.navy.mil 1 Space Science Division, U.S. Naval Research Laboratory, Washington D.C., USA Full list of author information is available at the end of the article is the presence of a denser and higher ionospheric layer at Earth that it is not found at Mars. This terrestrial layer is the so-called F-region, found at about 300 km. The F-region is controlled by a mix of photochemistry and dynamic processes, and is where the maximum electron density of terrestrial ionosphere occurs (e.g., Mendillo et al. 2016).
At night, the absence of solar photoionisation highlights the differences between the compositions of the Martian and terrestrial ionospheres. Photochemical recombination alters the ionospheric composition, leading to a Martian ionosphere that is still O + 2 above 130 km, but changes to having NO + as the most abundant ion below this altitude (e.g., Girazian et al. 2017). At Earth, O + remains the most abundant ion in the F-region, with O + 2 , N + 2 , and NO + playing important roles in the E-region (Schunk and Nagy 2009). The Martian nightside ionosphere near the terminator is maintained thanks to dayto-night plasma transport (e.g., Chen et al. 1978;Cui et al. 2015) together with impact ionisation by precipitating electrons (e.g., Cao et al. 2019), where densities of the major ions decrease with the solar zenith angle (SZA) across the terminator (e.g., Girazian et al. 2017;Withers et al. 2012). Far from the terminator, electron precipitation is the dominant source of ionisation, which is more important over regions of strong crustal magnetic fields (as will be described in the Section "Planetary environments") (e.g., Nĕmec et al. 2010;Lillis et al. 2018), and during solar storms (e.g., Sánchez-Cano et al. 2019). Transport also plays a large role maintaining the nighttime ionosphere at Earth, with current systems leading to E×B drifts and neutral winds driving field-aligned motions that may both lift and move plasma to less dense regions where ion loss through photochemical recombination occurs more slowly (e.g., Heelis 2004). These processes work alongside downward diffusion to maintain the nighttime F-region ionosphere at Earth.
The large exploration effort of Mars for the last 50 years, makes Mars the first body beyond Earth where comparative aeronomy studies start to be possible. Past comparisons between Earth and Mars environments have taken advantage of this increased body of ionospheric data. Mendillo et al. (2003) and Rishbeth and Mendillo (2004) found that day-to-day changes in solar irradiance have a similar impact on the photochemical ionospheric layers at Earth and Mars at different pressure levels. Moreover, this demonstrated that the major contributions to daily variability were correlated with solar irradiance, which had a larger effect at Mars than at Earth. Mendillo et al. (2003) also found that the influence of the lower atmospheric coupling is an important factor for the daily ionospheric variability at both planets. In particular, Bougher et al. (2004) found that the peak altitude of the main ionospheric layer at Mars has a strong annual oscillation when the same longitude is considered, and it is caused by nonmigratory thermal tides. This is a well-known phenomenon that also occurs at Earth (e.g., Hagan and Forbes 2002;Immel et al. 2006). However, some aspects of this coupling have some singularities at Mars. For example, there is a seasonal and dust-driven inflation and contraction of the entire Martian atmosphere (Bougher et al. 2006). The large amounts of suspended dust in the lower atmosphere, commonly produced during dust storms, lead to an increase in the altitude of the peak electron density that can be as large as ∼ 30 km (Wang and Nielsen 2003). Additionally, lower atmospheric cycles such as the CO 2 cycle appear to produce a seasonal increase in the total amount of ionisation in the Martian ionosphere .
Despite the large number of similarities between both ionospheres, the presence of an intense global magnetic field at Earth and the lack of one at Mars produces a completely different picture of the interaction of each planet with the solar wind. As will be described in the Section "Planetary environments", the terrestrial magnetic field extends far beyond the terrestrial surface and atmosphere to form a magnetosphere that interacts with the solar wind. The Martian ionosphere, on the other hand, is the primary obstacle encountered by the solar wind, causing the Interplanetary Magnetic Field (IMF) carried by the solar wind to drape around Mars. Thus, although the general ionosphere formation mechanism is the same at both planets, they are modulated differently by each magnetic environment. The objective of this paper is to assess how different the ionosphere formation is at Earth and Mars and explore their underlying processes, especially in the morning and evening sectors where the production via photoionisation and loss via recombination mechanisms are more significant. This kind of comparison is important because it provides insight on how similar ionospheric features behave under different magnetic circumstances, as well as insight on the evolution of the Martian ionosphere from the time when it had a global magnetic field. To that end, a similar analysis based on similar data sets at both planets has been performed, focusing on the dawn and dusk sectors at the winter and summer solstices.

Planetary environments
Before performing a comparison between the terrestrial and Martian ionospheres, it is important to understand the degree of similarity between these two planets. This section provides an overview of the similarities and differences in the planetary size, heliocentric distance, planetary tilt, magnetic fields, solar wind interactions, solar irradiance, and neutral atmosphere of both planets.
These discussions focus on the impact of these characteristics on the planetary ionosphere, to inform the results and discussion presented in the Section "Results and discussion".

Planetary size and orientation
The planetary size and heliocentric distance are two of the major differences between the Earth and Mars. Earth has a mean radius of 6371.0 km and a mass of 6.0 × 10 24 kg, yielding a surface gravity (g ⊕ s ) of − 9.8 m s −2 . At 350 km, the nominal altitude of the ion density peak, this translates to g ⊕ i = − 8.8 m s −2 . Meanwhile, Mars has a mean radius of 3389.5 km and a mass of 6.4 × 10 23 kg, yielding a surface gravity ( ) of − 3.7 m s −2 . At 135 km, the nominal altitude of the ion density peak, this translates to . The differences in gravitational acceleration affect the atmospheric scale heights, as well as other processes that may have played a key role in Mars evolution, such as the exosphere extension and atmospheric escape rates (e.g., Halekas 2017).
Mars has a more eccentric orbit than Earth, with aphelion at 1.66 AU and perihelion at 1.38 AU. In comparison, Earth has an aphelion of 1.02 AU and a perihelion of 0.98 AU. This means that the distance between aphelion and perihelion at Mars is a third of Earth's distance to the Sun, which has strong consequences for the amount of solar radiation that reaches Mars each season, and therefore, for the ionosphere. The eccentricity in the terrestrial orbit contributes to annual variations in the neutral density and TEC of the order of 3.4% and 7%, respectively (Azpilicueta et al. 2011). However, because the Martian eccentricity is seven times greater, the annual variations in neutral density and TEC are much larger. As a consequence, the strength of the annual variation of the ionospheric density due to the orbital eccentricity is very pronounced at Mars, with the Martian TEC having an annual sinusoidal variation that follows the EUV solar radiation flux as modulated by the Mars-Sun distance. Therefore, the orbital eccentricity plays a dominant role on the climatological behaviour of the Martian ionosphere (e.g., Bergeot et al. 2019). Moreover, the variations in the EUV solar radiation flux create changes in the atmospheric and plasma scale heights, electron and neutral temperatures, as well as in the balance pressure between the dynamic pressure of the solar wind and the thermal pressure of the ionosphere that affect the TEC (Sánchez-Cano et al. 2015aBergeot et al. 2019). Annual variations in the neutral atmosphere have also been identified (Zou et al. 2011;González-Galindo et al. 2013;Sánchez-Cano et al. 2018;Bergeot et al. 2019).
The axial tilt and length of day at both planets are very similar, leading to analogous diurnal variations in solar irradiance with season. At Earth, the axial tilt is 23.44 • and the length of day is 86,400 s. At Mars, the axial tilt is 25.19 • and the length of day (known as a sol) is 88,775 s.

Magnetic field
Charged particles are greatly influenced by the presence of a magnetic field. When the magnetic field is strong enough, it constrains transport in an ionosphere by 'freezing-in' the plasma to the magnetic field lines (e.g., Syrovatskii 1978). At Earth and Mars, major differences in the planetary magnetic fields affect both ionospheric structure and the interactions between the ionosphere and the solar wind.
The Earth has a global-scale internal magnetic field that extends far beyond the planetary atmosphere and interacts with the solar wind. To first order, this magnetic field is an off-center dipole that is compressed on the sunward and stretched on the anti-sunward side of the planet as the IMF embedded in the solar wind connects to the high-latitude terrestrial magnetic field lines and drags them away from the Sun before they reconnect tens of Earth radii away from the planet (Dungey 1961). Near the surface of the Earth, higher order deviations from the dipole field approximation become more important. Figure 1 shows the magnetic declination (top) and field strength (bottom) at 350 km above the surface of the Earth for 1 Jan 2006, as specified by the International Geomagnetic Reference Field (Macmillan and Finlay 2010). This altitude was chosen for Fig. 1, because it typically lies close to the altitude of the ion density peak.
The terrestrial ionosphere can be organized into regions based on the types of interactions with the geomagnetic field. At low altitudes (E-region) where the neutral density is large, the mass difference between ions and electrons allows neutral winds to drive currents that create electric fields. At higher altitudes (F-region), ions and electrons move together. This motion occurs preferentially along field lines, since the 'frozen-in' condition of the plasma requires an electric field large enough to move all of the plasma along a given field line to create field-perpendicular ion drifts.
At high latitudes, magnetic field lines are oriented vertically and extend far into the magnetosphere or are directly connected to the IMF. This leads to phenomena such as high-energy particle precipitation, aurora, and 'convective' motion (Kelley 2009). At low latitudes, the magnetic field lines are oriented horizontally, allowing rapid transport between hemispheres. Combined with the presence of a diurnally varying dawn-to-dusk electric field, this leads to the 'fountain effect' that creates the Equatorial Ionisation Anomaly (EIA), which is made up of large enhancements of plasma on either side of the magnetic equator (Hanson and Moffett 1966). At mid-latitudes, magnetic field lines connect to the plasmasphere, a torus of plasma in the inner magnetosphere that is both fed by and feeds into the upper ionosphere (e.g., Evans and Holt 1978). Light ions (primarily H + ) are the dominant constituents transported between the ionosphere and plasmasphere, which is typically directed upwards into the plasmasphere during the day and downwards into the ionosphere at night. However, during solstice seasons interhemispheric transport has been observed at night, transporting plasma from the summer hemisphere into the winter hemisphere at night. During low levels of solar activity this interhemispheric transport has been seen after sunset, but at high levels of solar activity it usually only occurs after midnight (Evans and Holt 1978;Kersley et al. 1978).
Mars does not have a global-scale internal magnetic field (Acuña et al. 1999), and as a consequence the ionosphere is the main part of the upper atmosphere that directly interacts with the solar wind. However, some areas over the Martian surface are magnetized (the so-called crustal magnetic fields) and their intensity can reach tens of nT at 300 km altitude. Although the strength of these localized Martian magnetic fields are much weaker than the terrestrial magnetic fields, they do shield portions of the planet from direct interactions with the solar wind. They create, in essence, multiple miniature magnetospheres that exhibit many of the same behaviours that are seen in high latitudes at Earth, such as aurora (e.g., Bertaux et al. 2005;Brain et al. 2010). The most intense crustal fields are located in the southern hemisphere between 120 • and 240 • longitude in horizontal bands of different polarities (e.g., Cain et al. 2003). The inclination of these fields has an important effect on the Martian TEC, such that the TEC is enhanced over quasivertical crustal magnetic fields, and reduced over quasihorizontal crustal magnetic fields (Nĕmec et al. 2010;Cartacci et al. 2013). These enhancements are related to changes in the scale height of the ionosphere, possibly caused by heating of the ionosphere by solar wind electrons that reach low ionospheric altitudes along the nearly vertical (open) magnetic field lines (e.g., Duru et al. 2006). Because crustal fields are localized, large scale transport processes (such as the terrestrial high-latitude plasma convection and low-latitude EIA) cannot occur. Nevertheless, day-to-night ionospheric transport plays an important role in the Martian ionosphere. Crossterminator transport affects Martian ion losses through escape processes (e.g., Kar et al. 1996), with the oxygen loss rate from ion transport making up ∼ 5% of the total oxygen loss rate (e.g., Cravens et al. 2017) depending on the solar activity level ( e.g., Fox 1997Fox , 2009, the ambient magnetic pressure, and the crustal magnetic field orientations (e.g., Wu et al. 2019). These ion losses are potentially important at Mars for the evolution of volatiles.

Solar wind interaction
The lack of a global intrinsic planetary magnetic field at Mars (and the presence of such a magnetic field at Earth) also cause differences in the way each planet interacts with the solar wind. As mentioned in the previous section, the terrestrial magnetic field extends far beyond the terrestrial surface and atmosphere to form a magnetosphere that interacts with the solar wind. This interaction shapes much of the characteristics of the magnetosphere, causing the formation of a bow shock, an upstream foreshock, a magnetosheath, a magnetotail, and an inner magnetosphere. At Earth, these regions contain plasma from both the solar wind and the Terrestrial atmosphere, with the solar wind providing a large amount of ionised hydrogen (e.g., Welling et al. 2015).
Mars also has a magnetosphere with the same basic characteristics as the Terrestrial magnetosphere. However, the magnetospheric plasma is primarily of atmospheric origin, since the Martian ionosphere acts as the primary obstacle for the solar wind . The relationship between the thermal pressure of the ionosphere and the solar wind dynamic pressure varies with the solar cycle , and it delimits the location of the Martian bow shock, which is more sensitive to variations in the solar irradiance than to solar wind dynamic pressure variations (Hall et al. 2016b). Because of this relationship, the solar wind is able to induce an horizontal magnetic field on the dayside ionosphere of Mars (e.g., Shinagawa and Cravens 1989;Ramírez-Nicolás et al. 2016), the magnitude of which is controlled by the strength of the ionosphere. Since the ionospheric strength varies with the solar cycle, there is a larger probability for a more magnetized ionosphere during low solar activity periods (Sánchez-Cano et al. 2015a. The Mars-solar wind interaction causes the IMF carried by the solar wind to drape around the planet, forming the Martian induced magnetosphere. The localized crustal magnetic fields also contribute to the magnetosphere, though their contributions perturb the bow shock and magnetosheath structure formed by the IMF. The crustal field-ionosphere-solar wind interaction is very complex, producing large ionospheric variability (e.g., Mazelle et al. 2004;Brain et al. 2005;Edberg et al. 2008;Cartacci et al. 2013;Hall et al. 2016a).

Thermosphere
The Terrestrial and Martian thermospheres make up the upper portion of the planetary neutral atmospheres, extending from the mesopause to the exobase. These thermospheres are both characterized by a temperature gradient that increases with increasing altitude, at first rapidly before slowing to asymptotically reach the exobase temperature. Both thermospheres are strongly impacted by coupling from above (through solar radiation, solar wind particles, and electric fields) and below (through atmospheric gravity waves, planetary waves and tides, and lower atmospheric weather) (e.g., Meriwether 1983;Wang and Nielsen 2003;Bougher et al. 2004Bougher et al. , 2006Lühr et al. 2007;Withers 2009;Emmert et al. 2010;Eckermann et al. 2018;Sánchez-Cano et al. 2018;Gupta et al. 2019). The majority of both planetary ionospheres also form within their respective thermospheres.
As weakly ionised plasmas, the Terrestrial and Martian ionospheres are highly coupled with their thermospheres. The ionised species are drawn from the local neutral species and created through photoionisation, precipitating particles, and ion-neutral photochemical reactions. The structures of both ionospheres are affected by the differences in the neutral species, the temperature of the planetary atmosphere, and dynamical processes in the thermosphere.
Although conceptually similar, the details of the Martian and Terrestrial thermospheres cause important differences between their ionospheres. At Mars, the major thermospheric species ( CO 2 , O, N 2 , and CO) lead to an ionosphere that is dominated by O + 2 . At Earth, the major thermospheric species (O, O 2 , N, and N 2 ) lead to an ionosphere that is dominated by O + at the ion density peak. This results in an important difference between the Terrestrial and Martian ionospheres: at Mars the major species are lost through one-stage recombination processes, but at Earth O + is lost through a slower, two-stage process involving reactions with both N 2 and O 2 (Schunk and Nagy 2009).
The thermospheric composition also means that both ionospheres are affected by EUV irradiance. Although the Earth and Mars are sensitive to slightly different wavelengths within the EUV range (discussed in more detail in the Section "TIMED/SEE EUV"), their variations on short and long scales are similar. Thus, photochemical ionospheric layers in both the Terrestrial and Martian ionospheres have been shown to react with magnitudes proportional to their distance from the Sun (Mendillo et al. 2003).
Both planets have ion density anomalies that are tied to the neutral atmosphere. At Earth the semi-annual anomaly, characterized by thermospheric and ionospheric densities that are higher than expected at equinox and lower then expected at solstices. This can be explained by global neutral atmospheric circulation in the mesosphere and thermosphere that affects the ionosphere through changes in the O:N 2 ratio (Fuller-Rowell 1998; Jones et al. 2018). There are also seasonal and annual anomalies that are characterized by higher mid-latitude peak ion densities in the winter than the summer hemisphere (Yonezawa 1959) and a higher level of ionospheric plasma density during northern hemisphere (NH) winter/southern hemisphere (SH) summer (the December solstice) than NH summer/SH winter (the June solstice) (Yonezawa and Arima 1959), respectively. The seasonal anomaly is again understood to be caused by changes in the thermospheric composition, affecting the ionosphere by altering the O:N 2 ratio (e.g., Fuller-Rowell et al. 1988). The smaller annual anomaly (or annual asymmetry) is still an area of active investigation (e.g., Azpilicueta et al. 2011).
At Mars, a large increase in global plasma density is observed at L S = 230 • (Northern autumn and Southern spring). This increase is related to a warmer and thicker thermosphere (seasonal atmospheric expansion), which is caused by the orbital eccentricity that places Mars closer to the Sun at this time. Moreover, the CO 2 sublimation of the South polar cap at the surface and lower atmosphere contributes to a plasma increase during these seasons. The influence of the CO 2 on the thermosphere is also notable at L S = 40 • (Northern spring and Southern autumn), when the North polar cap sublimation increases the plasma, as shown in Fig. 4. However, the increase is smaller in this case because Mars is further position from the Sun than it was at L S = 230 • . Simulations suggest that this increase in plasma density is caused by the sublimation of CO 2 at the surface and lower atmosphere that increases the O, O 2 , and N 2 populations in the thermosphere. As discussed above, these populations are sources for the major ionospheric species and so lead to an increase in the global TEC as discussed by Sánchez-Cano et al. (2018), and also seen by other authors (Bergeot et al. 2019).
The ionospheres at both planets also experience a large degree of daily variability due to interactions with the lower atmosphere. Many of these lower atmospheric drivers are the same, though the coupling processes and affects differ. For example, both planets are affected by gravity waves (Yiğit et al. 2009(Yiğit et al. , 2015 and the influences of topography (Wang and Nielsen 2004;Immel et al. 2006). However, at Earth coupling between the lower atmosphere and upper ionosphere often occurs through electric fields that map along magnetic field lines and are driven by currents created through interactions between the neutral atmosphere, ionosphere, and magnetic field at lower altitudes.
Another important source of coupling between the lower atmosphere and ionosphere is through tropospheric weather events. At Mars, dust storms are huge events that can encompass the entire planet. These storms introduce particles into the lower atmosphere that heat the lower atmosphere by radiative transfer. This produces a subsequent thermal expansion of the lower atmosphere that lifts pressure levels in the upper atmosphere, lifting also the ionosphere ( e.g., Withers 2009). Effects of seismic disturbances and tropospheric weather also exist at Earth, but their impacts are weaker than the other processes discussed here in more detail (e.g., Occhipinti et al. 2008;Kartalev et al. 2004;Goncharenko et al. 2010;Ogawa et al. 2012;Frissell et al. 2016).

Data
This study uses measurements of solar EUV irradiance and TEC to examine and compare the ionospheres at Earth and Mars over a period of time beginning during the descending phase of the 23 rd solar cycle and ending during solar maximum of the 24 th solar cycle. Although a simple measure of ionospheric column density, the TEC has been used to study many 'anomalies' caused by unknown (at least at the time of their naming) interactions in the upper atmosphere (e.g., Azpilicueta et al. 2011). For both planets, the TEC are provided in TEC units (TECU), where 1TECU = 10 1 6 e − m − 2 . Empirical and first principles models are also used when measurements are lacking and physical insight into the observed characteristics are needed, respectively. These measured and modeled data sources are described in the following subsections.

GNSS TEC
The TEC at Earth is measured using transmissions from Global Navigation System Satellite (GNSS) receivers, which detect signals from the Global Positioning System (GPS), the Global Navigation Satellite System (GLO-NASS), and Galileo satellite constellations processed at MIT Haystack and provided by the CEDAR Madrigal database (Rideout and Coster 2006). This vertical TEC and its uncertainty are provided at a 1 • latitude × 1 • longitude resolution in geographic coordinates every second from 1998 onward. Coverage is provided from all publicly available ground and space-based receivers. Figure 2 shows the available coverage of the ionospheric pierce points at a 1 • latitude × 1 • longitude resolution for a day near the beginning and end of this study.

MARSIS TEC
The Mars Advanced Radar for Subsurface and Ionosphere Sounding (MARSIS) radar (Picardi et al. 2004;Orosei et al. 2015) on board Mars Express (Chicarro et al. 2004) has been sounding the planet since mid-2005, currently covering a full solar cycle of observations. MARSIS can be run in two different operational modes: (1) The Active Ionospheric Sounding (AIS) mode (Gurnett et al. 2005), a topside ionospheric sounder that measures topside electron density, from which the topside TEC can be obtained.
(2) The subsurface mode, in which the radar sounds the surface and subsurface of the planet to measure materials beneath the surface of the planet.
In the second mode, the TEC of the entire atmosphere is retrieved as a by-product of the signal distortion analysis (Safaeinili et al. 2007;Mouginot et al. 2008;Cartacci et al. 2013Cartacci et al. , 2018. This study is limited to MARSIS data from the subsurface operational mode. The TEC was obtained through the Cartacci et al. (2013) algorithm, after selecting for data with a Signal-to-Noise Ratio (SNR) larger than 20 dB and carrier frequencies of 4 and 5 MHz. This conservative approach of using only strong signals measured at higher frequencies guarantees a good data quality. However, because the carrier frequencies are very close to the maximum plasma frequency of the Martian dayside ionosphere ( ∼ 3 MHz), the radar cannot work on the proper Martian dayside. This limits the TEC observations to SZA > 50-60 • (Sánchez-Cano et al. 2015b). Moreover, due to the evolution of the Mars Express orbit and observation planning priorities, there is a large asymmetry on the amount of data collected between dawn and dusk sectors, being on average for the full mission of ∼88% larger for the dusk sector . Note, as well, that typically ∼70% of the electron density contribution to the TEC is obtained from the region near the height of maximum ionisation in the ionosphere (photochemical region) (e.g. Sánchez-Cano et al. 2013), and so is an excellent indicator of the variability of the photochemical region of the ionosphere (e.g. Mendillo et al. 2015).

TIMED/SEE EUV
The Solar EUV Experiment (SEE) is one of four scientific instruments flying onboard the NASA Thermosphere Ionosphere Mesosphere Energetics Dynamics (TIMED) satellite. TIMED/SEE measured the solar irradiance from 0.1 to 194 nm. In the higher level data products, the observations are adjusted to reflect the solar irradiance at 1 AU (Woods et al. 2005). Figure 3 shows two different measures of solar irradiance obtained from TIMED/SEE data. The top panel shows the 30.5 nm irradiance scaled to the Mars heliocentric distance following an inverse-square law and assuming no significant differences in solar longitude between Earth and Mars. This wavelength is typically used at Mars since it is responsible for the majority of the photoionisation (Rishbeth and Mendillo 2004). The bottom panel shows the O + power at Earth, calculated based on the process outlined in Richards et al. (1994). This measure of EUV irradiance was chosen because O + is the dominant ion in the F region.

Horizontal wind model
Neutral winds play an important role in the Terrestrial F region ionosphere, as noted in the Section "Thermosphere". However, these winds are difficult to measure.
This study uses data from the Horizontal Wind Model (HWM) instead of observations (Drob et al. 2015). The HWM incorporates observations from many overlapping and disparate, ground-and space-based data sets. The resulting model uncertainties reflect the natural variability of the winds, the uncertainties of the observational data, and the spatiotemporal coverage of the observational data. HWM captures variations due to geographic location, time of day, and time of year. Solar cycle variations are not modeled, due to insufficient data coverage. Additionally, high altitude thermospheric winds (above 250 km) are not routinely measured and so the model is not as trustworthy at these heights.

Mars climate database
The Mars Climate Database (MCD) Version 5.3 is a data base of meteorological parameters obtained from a general circulation model of the Martian atmosphere (Forget et al. 1999;Millour et al. 2018). These meteorological parameters, which include (but are not limited to) atmospheric temperature, density, pressure, and winds, have been validated against available observational data and are freely available upon request or through a simple online interface. This study uses the neutral temperature and density near the ionospheric peak altitude to examine the influence of the Martian thermosphere on the ionosphere.

Analysis
Due to the differences in the length of year, orbital eccentricity, and thermospheric interactions, ionospheric density comparisons between the Martian and Terrestrial ionospheres are difficult. To permit comparisons, this study considers different ratios of the rate of change of TEC as a function of SZA for several spatiotemporal regions. The spatiotemporal regions are discussed in the Section "Data selection" and the method used to obtain the rate of change of TEC as a function of SZA is discussed in the Section "Linear fitting".

Data selection
Data selection in this study uses both Martian and Terrestrial conventions, as deemed appropriate. Data are grouped by Martian Year (MY) (and limited to the period of time when Martian TEC is available to ensure a similar coverage), but seasons are centered about the solstices, as is the custom at Earth. The seasons for each MY span the solar longitudes ( L S ) shown in grey boxes in Fig. 4. NH summer/SH winter contains data with 25 • < L S < 155 • , while NH winter/SH summer contains data with 205 • < L S < 335 • . The seasonal range was chosen to improve planetary coverage at Mars and include data with similar solar irradiance conditions, as shown in Fig. 4. Note that the NH summer/SH winter season happens when both Earth and Mars are furthest from the Sun, with NH winter/SH summer occurring when the two planets are closest to the Sun. This causes an additional variation in the amount of solar irradiance for each season that is much smaller at Earth than it is at Mars. Moreover, the data have been divided up by MY instead of EUV flux levels to focus on the study of the annual seasonal variations. This has the additional benefit of separating data into different phases of the solar cycle, as each MY spans at most one phase (ascending, maximum, descending, or minimum) of the solar cycle.
The spatial regions were chosen to provide similar seasonal and magnetic environments. At Mars, this means that the Northern and Southern hemispheres are divided at the equator and the Southern hemisphere is further divided into longitudes ( ) with ( 150 • ≤ < 270 • ) and without ( −120 • ≤ < 150 • ) crustal magnetic fields. At Earth, eight mid-latitude regions are selected, as shown in Fig. 1. High latitudes were not considered to reduce the influence of magnetosphere-ionosphere interactions. Low latitudes were excluded to reduce the influence of the plasma transport that forms the EIA. The eight midlatitude regions were chosen to encompass similar magnetic field orientations and strengths, as shown in the top and bottom panels of Fig. 1. Four SZA regions were chosen based on the location of the solar terminator. To separate data before and after noon, SZAs before noon are labeled with negative values and those after noon remain positive. At each planet, the location of the solar terminator was calculated at different important altitudes (such as the typical height of the ion density peak). This informed the SZA regions shown in Table 1. Note that Earth is more restrictive than Mars, because the wealth of available data made it possible to remove SZAs that encompass the period of time when the solar terminator enters and exits the F region.

Linear fitting
Linear fits to the TEC in each year, season, location, and SZA region are made to determine the typical behavior of the TEC as a function of SZA and aid in the comparison between this behavior at Mars and at Earth. Due to the large amount of variability in this data and the existence of outliers, a robust linear regression method known as the Huber method was used (Huber 1963). This method considers the influence of outliers without allowing their presence to dominate the regression by identifying the outliers and applying a linear loss to those samples when performing the linear regression. The Python implementation of this method, sklearn.linear_model.HuberRegressor, was used in the Terrestrial portion of this study (Pedregosa et al. 2011), while the MatLab implementation was used for the Martian data (The MathWorks, Inc. 2019). An example of how the Huber Regressor method performs on the Martian and Terrestrial TEC is shown in Fig. 5. The difference between standard linear and Huber regressions are shown in panel (a), highlighting the way that outliers can overly influence a standard linear fit. The selection of outliers and inliers is highlighted for the Terrestrial data in panel (b); note that the spread of the inliers changes in each SZA region based on the distribution of TEC in that region. Panels (c-f ) show both how the Huber method selects the most representative data as inliers, and how well the Huber fit represents the distribution of the data in the selected SZA region. The Huber fit slopes and intercepts, as well as the medians and quartiles of the differences between the dawn and dusk Huber fits and the TEC observations for each planet, region, and MY are included in Additional file 1. The small medians and symmetric upper and lower quartiles ( q 3 and q 1 respectively) provide confidence in the Huber slopes and intercepts.

Martian TEC
This section focuses on the morning and afternoon sectors (dawn and dusk slopes, SZAs up to ± 105 • ) of the photochemical region of the ionosphere. This is because most of the Martian TEC is obtained from the region where plasma is near and above the height  of maximum electron density of the ionosphere (e.g., Mendillo et al. 2015). In this region, plasma is controlled by a photochemical equilibrium, meaning that production via photoionisation and losses via recombination tend to be equal. Nevertheless, there are other processes that dominate at higher altitudes (typically above ∼200 km) for all SZAs, such as vertical transport (e.g., Chen et al. 1978;Mendillo et al. 2011;Kar et al. 1996;Fox 1997Fox , 2009Cravens et al. 2017;Wu et al. 2019, ), which are not investigated in this study because their contribution to the TEC is smaller. In addition, at high SZAs and far from crustal field regions, dayto-night plasma transport plays an important role in shaping the Martian ionosphere near and beyond the terminator (e.g., Nĕmec et al. 2010;Withers et al. 2012;Cui et al. 2015;Girazian et al. 2017;Cao et al. 2019).
Plasma transport into the nightside is not investigated in this study, although it should be noted that the low TEC values beyond 90 • SZA should be affected by this process (e.g., Cui et al. 2015;Cao et al. 2019). The Huber fits to Martian TEC as a function of SZA in the four SZA regions are shown in Fig. 6 for the three location regions, both solstices, and each MY. The MYs are color coded such that darker colors indicate lower levels of solar activity, though the amount of irradiance in each MY varies based on the heliocentric distance, as previously shown in Fig. 3. Due to coverage constraints, not all regions have Huber fits every MY. Coverage is greater at dusk than dawn, and is worst in the pre-dawn region (motivating this study consider only dusk and dawn in depth). shows a single MY, season, and location region for Earth, identifying both inliers (grey) and outliers (black). Panels (c-f) show the differences between the Huber Regression fits and measured TEC at dawn (Panels (c) and (e)) and dusk (Panels (d) and (f)) at Mars (Panels (c) and (d)) and Earth (Panels (e) and (f)). In each of these panels the outliers are plotted in black and the inliers in grey, while the median of all differences is shown by a vertical line Despite these constraints, it is apparent that the TEC slopes increase in magnitude as solar activity increases, for all location regions and seasons. This behavior is similar for the dawn and dusk slopes, as is expected for an ionosphere dominated by photochemical processes. The similarities between the TEC slopes indicate that the presence of crustal magnetic fields have a localized influence that does not affect the regional production via photoionisation and loss via recombination processes.
Comparing the top row (NH summer and SH winter) to the bottom row (SH winter and SH summer) demonstrates that heliocentric distance greatly affects the production of TEC. Regardless of the season in that hemisphere, the maximum levels of TEC are higher in the bottom panel than they are in the top panel. Figure 7 examines this in more detail, by taking the summer:winter ratio of the dusk TEC slopes for each location region and MY. Dawn summer:winter ratios  were not examined, due to the constraints on coverage over the solar cycle. Figure 7 clearly shows that heliocentric distance has a greater influence on the rate of ionospheric plasma production via photoionisation and loss via recombination, as greater loss rates (as captured by the dusk TEC slope) occur during NH summer and SH winter over the entire planet. This hemispheric difference in the seasonal variation (as indicated by a ratio less than unity in the NH and greater than unity in the SH) can be attributed to a combination of the orbital impact on photoionisation rates and the neutral atmospheric density (as discussed in the Planetary size and orientation and Section "Thermosphere").
To assess the impact of neutral atmosphere on the hemispheric differences in the dusk seasonal variation, atmospheric variations in density and temperature were retrieved from the MCD (see the Section "Mars climate database") for both solstices at mid-latitude locations in the NH (longitude of 0 • ), SH over crustal fields (longitude of 180 • ), and SH without crustal fields (longitude of of 0 • ), as shown in Fig. 8. This figure plots the neutral temperature (black) and density (blue) at 135 km above the reference level (called "sea level" at the MCD) as a function of SZA. Each column shows a different aerographic region during the NH summer/SH winter (top row) and the NH winter/SH summer (bottom row). Comparing the global densities in the NH summer (top row) and SH summer (bottom row) reveals that neutral densities at ionospheric altitudes are higher in both hemispheres during SH summer (bottom row) than in NH summer (top row). Thus, both the direct influence of a decrease in photoionisation rates due to distance and a change in the neutral atmospheric density lead to summer:winter ratios that are less than unity in the NH and greater than unity in the SH. Figure 7 also reveals the strong modifying influence of the solar cycle on the summer:winter ratio. MYs with large changes in solar activity, most notably MY 30 (which occurred during the end of the low solar activity phase and the ascending phase of the 23 rd solar cycle), are characterized by seasonal ratios that deviate most from unity. The hypothesis that the changes in solar cycle over the course of MY 30 are causing the most extreme seasonal ratios in Fig. 7 was tested using the NeMars model (Sánchez-Cano et al. 2013. This model uses F 10.7 measured at Earth to capture solar irradiance changes at Mars' heliocentric distance, considers the phase of the solar cycle, and solves for ionospheric density using the Chapman formulation (Chapman 1931). The modeled TEC varying between SZA 50 • and 90 • was calculated for the heliocentric distance of the Martian aphelion (1.66 AU) and perihelion (1.38 AU) for each Martian year (MY27-MY32) and Fig. 8 Modelled neutral atmospheric temperature (black) and density (blue) as a function of solar zenith angle for NH summer/SH winter (top row) and NH winter/SH summer (bottom row) at locations that fall within the three Martian location regions corresponding phase of solar activity. Then, a fit using the Huber Regressor method was performed, similar to the procedure done for the observed data. The results, presented in Figure 9, show a sharp change in the summer:winter sunset slope ratio for MY 30 due to the changes in solar irradiance over this course of the MY.
To summarise, at Mars the rate of change of TEC as a function of SZA at sunrise and sunset are strongly influenced by changes in solar irradiance. The solar cycle has the largest effect (as demonstrated by the changes in dawn and dusk slopes in Fig. 6), followed by heliocentric distance (as demonstrated by the seasonal dusk ratios in Fig. 7), and finally season. The presence of crustal magnetic fields may modify the ionosphere locally, but these changes are not evident when considering the entire portion of the SH dominated by crustal magnetic fields.

Terrestrial TEC
The Huber fits to Terrestrial TEC as a function of SZA in the four SZA regions are shown in Figs. 10 and 11 for the four longitude regions, both solstices, and each MY in the NH and SH, respectively. These figures have a similar layout to Fig. 6, with the same color codes for the MYs. As was the case in Fig. 6, the dawn TEC slopes are large and positive following a relatively flat pre-dawn TEC variation with SZA. At dusk and night, however, the TEC variation with SZA is not as consistent. Figure 10 shows that the northern, winter hemisphere has evening loss rates that are larger at dusk than at night. However, this trend reverses in the summer. One explanation for this behavior is the seasonal anomaly (previously discussed in the Section "Thermosphere"), which changes the O:N 2 ratio so that the ion production via photoionisation and loss via recombination rates are greater in the winter than in the summer (Fuller-Rowell et al. 1988). Increasing the evening ion loss rate reduces the relative importance of F region sustaining plasma transfer processes, such as plasma transport into the ionosphere from the plasmasphere. Recall that seasonal variations in mid-latitude interhemispheric transport showed plasma transported from the summer hemisphere to the winter hemisphere in the evening (Kersley et al. 1978), Fig. 9 Ratio between the summer and winter Martian dusk slopes for TEC as a function of SZA for the NH and SH five Martian years using the NeMars model This figure also shows that as the level of solar irradiance increases over the solar cycle, the rate of ion loss, as indicated by the magnitude of the TEC slope, at dusk increases. This again reflects the importance of ionosphere-thermosphere coupling. As the thermosphere expands with increasing solar irradiance, photochemical processes become increasingly important.
This variation is also consistent with the solar cycle variation in mid-latitude interhemispheric transport, which occurs at earlier local times at solar minimum and after midnight at solar maximum (Evans and Holt 1978;Kersley et al. 1978). The existence of this interhemispheric transport points to the presence of upward, field-aligned ionospheric flows in the summer hemisphere that lead to deviations from photochemical equilibrium. When this interhemispheric transport is not present, the midlatitude F region ionosphere may then be inferred to lie closer to photochemical equilibrium. Figure 11 shows the same quantities as Figs. 6 and 10, but for the Terrestrial SH regions. The winter TEC at dawn and dusk show greater rates of ion production via photoionisation and loss via recombination, respectively, than their summer counterparts (as was the case in the NH). There are important differences between the two hemispheres, though. The middle rows of Fig. 11 (longitudes between 90 • and 270 • ) show summer and solar cycle behavior similar to that seen in the NH. However, the other longitudes show a TEC dependence on SZA that is more similar to Mars than the Terrestrial NH (with dawn and dusk slopes of similar magnitude). This behavior is strongest between 270 • and 360 • longitude.
These location regions are characterized by the lowest magnetic field strength and small declination angles. The orientation of the magnetic field modulates the influence of the neutral wind, allowing the horizontal wind to induce field-aligned ion drifts through collisions. Fig. 12 shows the field-aligned component of the HWM14 neutral wind ( u || ) for the December (top) and June (bottom) solstices at dusk (SZA = 70 • ) for three selected locations (the centers of three of the geographic regions shown in Fig. 1) that represent the typical global variations. This discussion focuses on altitudes between about 150-300 km (typical altitudes of bottomside F region), where u || plays the largest role in collisionally raising and lowering the plasma density peak. In December, both hemispheres have u || aligned downward (positive in the north and negative in the south). This is also true in June except at longitudes between 270 • and 90 • (bottom center panel). This is exactly the opposite of what was expected based on the SH TEC slopes, since an upward u || will raise the F 2 peak and inhibit ion loss processes. This implies that the neutral wind does not play a significant role in raising the mid-latitude F 2 peak at dusk during the solstices.
The other primary transport mechanism on Earth at the mid-latitudes is the E×B drift. As discussed in the Section "Planetary environments", a positive meridional E× B drift will move the plasma upwards and polewards, Fig. 11 Fitted Earth TEC as a function of SZA for the four southern location regions and six Martian years while a negative meridional E×B drift will move the plasma downwards and equatorwards in both hemispheres. In the evening, this means that the terrestrial ionosphere will be maintained for a longer period of time when a positive meridional E×B drift is present, because the plasma will be lifted to higher altitudes where the ion loss rate is lower. Conversely, a negative meridional E×B drift moves any available plasma to lower altitudes where the ion loss rate is higher.
The canonical picture of meridional E×B drift has it directed positively during the day and negatively at night (Kelley 2009). This behavior, when coupled with seasonal variations in photoionisation, explains why the evening Terrestrial TEC slopes are flatter during the summer than during the winter. The seasonal variation can be explained not by a seasonal variation in the direction of the meridional E×B drift, but by the continued presence of photoionisation in the high-latitude ionosphere during the summer. This provides a continual source of plasma to be transported to lower latitudes and altitudes by a negative meridional E×B drift that is not present in the winter. However, this behaviour varies as a function of solar cycle, over the course of the year, and as a function of longitude due to the way that the magnetic declination and neutral wind modify the Pedersen and Hall conductivities (Maute and Richmond 2016;Pacheco et al. 2010). For example, Pacheco et al. (2010) showed that the meridional E×B drift at 19:00 solar local time (a local time that falls between 50 • and 102 • at these latitudes) in general changes from positive to negative with increasing apex altitude. This shows the transition from the low latitude meridional E× B drift, which is strong and directed upward after sunset at higher levels of solar activity, to the mid-latitudes, which has a smoother transition from upward daytime drifts to downward nighttime drifts. However, this variation differs as a function of season and longitude. Specifically, at 0 • , 310 • and 290 • the meridional E×B drift remains positive for a larger range of apex altitudes during the December solstice than it does during the June solstice. This has stronger implications for the SH than the NH, since it means that at these longitudes the evening meridional E×B drift will not be moving plasma from high to mid latitudes, leading to the steep summer evening slopes observed in Fig. 11 between 270 • to 360 • (and to a lessor extent 0 • to 90 • ) longitude. Figure 13 shows the summer:winter slope ratios for the Huber fits to TEC as a function of SZA for 50 • ≤ SZA ≤ 102 • . The influence of the seasonal anomaly is seen in most of the location regions. As expected from Figs. 10 and 11, the summer:winter ratio is below one for all the NH longitudes and the SH between 90 • and 270 • longitude. At SH longitudes between 0 • and 90 • the ratio is near one, with a larger value for MY 29 (during solar Further atypical behaviour seen in the SH between 270 • and 360 • during MY 29 (solar minimum) suggests that a process in the winter exists to lift the plasma, slowing the recombination rate at a time when then thermosphere was greatly constricted. Figure 3 shows that MY 29 contains the minimum between the 23 rd and 24 th solar cycles, which was very low and has less variation than the solar irradiance during any other MY. Thus, the explanation for the anomalous seasonal ratio at Mars cannot explain what is happening here. However, during this extreme solar minimum, the upper atmosphere was contracted and coupling between the lower atmosphere and ionosphere much more visible than normal. It is possible that lower atmospheric forcing worked together with the contracted thermosphere to cause a slower winter loss rate than were seen elsewhere. This reinforces the interpretation that the terrestrial, mid-latitude ionosphere is more photochemically driven during high levels of solar irradiance and more dominated by transport as the thermosphere contracts.

Mars and Earth comparison
There is a large difference in ionospheric density between Mars and Earth, which makes direct comparisons of electron density difficult. To get around the differences in TEC magnitude, direct comparisons between Mars and the Earth are made using the dimensionless ratio between the dawn and dusk TEC slopes shown in Figs. 6, 10, and 11. Recall that Table 1 shows that dawn and dusk are defined as containing data with SZA between ±50 • and the SZA of the terminator at ionospheric altitudes. Figure 14 shows these ratios for all location regions and seasons at Mars and three regions at Earth. These three regions were chosen for closer examination based on their typical or atypical behavior when compared to the rest of the geographic regions. These regions are shown in panels (b) and (d), and plotted on the same vertical scale as the Martian ratios in panels (a) and (c). The Martian ratios are shown on a zoomed scale in Fig. 15, while the remaining Terrestrial ratios are shown in Fig. 16.
As discussed in the Section "Martian TEC", Martian models typically assume strict photochemical equilibrium. This assumption results in a dawn:dusk ratio of unity, which is shown to be a reasonable approximation at Mars when compared to most locations at Earth (for example, compare Fig. 14c and d). There appears to be no clear solar cycle variation in the Martian dawn:dusk ratio and no seasonal dependence. However, the closer examination presented in Fig. 15 reveals that there are times in each of the Martian location regions when the dusk slope is steeper than the dawn slope. The largest differences are seen in the SH between 150 • and 240 • longitude, where the crustal magnetic fields are strongest, and during the maximum and descending phases of the solar cycle. This reveals that changes in solar irradiance primarily affect photochemical processes, and will not affect the Martian Fig. 13 Ratio between the summer and winter TEC slopes for SZA between 50 • and 102 • dawn:dusk ratio. However, when crustal fields are present, processes such as interactions with the solar wind lead to departures from photochemical equilibrium (e.g. Andrews et al. 2015). Our results agree with Matta et al. (2015), who found that strong crustal magnetic field lines alter the structure of the ionosphere mainly by fieldaligned plasma transport, and this is significantly notable at dusk at high altitude over closed field-lines. A possible explanation for the dawn:dusk ratio solar cycle dependence over the crustal magnetic fields could be that this solar wind interaction with the magnetic field regions could lead to faster recombination at dusk with increased IMF, through processes such as field-aligned plasma transport that also exist in the Terrestrial high latitudes. This variation is characterized by increasing similarly to the Terrestrial ionosphere during NH winter (compare Fig. 15 and Fig. 14d) with increasing solar activity.
In contrast to the Martian ratios, the Terrestrial dawn:dusk ratio show both seasonal and solar cycle variations in most of the location regions. An exception to this is the SH region between 270 • and 360 • . This longitude region, shown in upward pointing triangles in Fig. 14b and d, has a dawn:dusk ratio that shows little seasonal or solar cycle variation and has values as close to unity as the Martian ratios typically are. This indicates that there are places at Earth where photochemical processes are more important than transport when considering the dawn and dusk plasma production via photoionisation and loss via recombination in the F region ionosphere (places where the Terrestrial F region ionosphere behaves more like the Martian ionosphere). However, the strong seasonal and solar cycle variations displayed in the other Terrestrial regions in Fig. 14 suggest that the amount of solar irradiance and the composition of the thermosphere both play a role in the relative importance of photochemical or transport processes. Figure 14b shows that on Earth in the NH summer, the selected NH region behaves more photochemically (has a dawn:dusk ratio closer to unity) as solar activity increases. This is consistent with an expanded  thermosphere allowing plasma recombination to take place at a greater range of heights (minimizing the influence of a downward meridional E×B drift in the summer). The dawn:dusk ratios for the SH summer regions (Fig. 16a), however, do not show a clear solar cycle variation. This may be attributable to the presence of the Weddell Sea anomaly, a large plasma enhancement that only forms in the SH summer and extends through much of the mid-latitudes at dusk (Jee et al. 2009). Figure 14d shows that in the winter, the selected NH longitude sector behaves less photochemically at higher levels of solar activity. Additionally, the dawn:dusk ratio is now consistently below unity. This implies that either winter transport in the NH is leading to plasma loss via recombination or that the dawn slope is not solely shaped by plasma production via photoionisation processes. The latter is more likely to be the culprit, as Fig. 10 clearly shows the amount of TEC at the end of the dawn period to be far below the amount of TEC present at the start of the dusk period. Figure 14b and d shows that this is not the case in the SH, where the dawn:dusk ratio is always at or above unity. Instead, Fig. 14b shows that the SH, again, shows no clear solar cycle variation in the dawn:dusk ratio. However, Figs. 14b and 16a show the dawn:dusk ratio in the NH at high levels of solar activity during the summer is approximately equal to the dawn:dusk ratio seen at all levels of solar activity in the SH during the winter. These hemispheric and seasonal variations in the Terrestrial dawn:dusk ratio are consistent for the remainder of the NH and SH Terrestrial longitude regions shown in Fig. 16.

Conclusions
This paper used the similarities between the ionospheres on Mars and Earth to examine the relative importance of photochemical and transport processes as functions of location, solar cycle, and season. This was done by comparing the rate of change of TEC as a function of solar zenith angle during the solstice seasons at dawn and dusk. The use of the TEC slopes made it possible to compare production via photoionisation and loss via recombination processes in the main layer of each planetary ionosphere despite the extreme differences in the total quantity of plasma present in the two ionospheres.
To first order, the dawn and dusk slopes that denote the rate of ion production via photoionisation and loss via recombination at Mars are symmetric. This symmetry is interpreted as an indicator of photochemical equilibrium. Although the F region ionosphere at Earth is known to be dominated by transport, there are times and locations where the ratio between the dawn and dusk slopes were found to lie near unity, indicating the local dominance of photochemical processes. Both planets also had times and locations where the dawn and dusk slopes differed significantly. By comparing the two planetary ionospheres at dawn and dusk, these regions were identified and the underlying processes responsible for the plasma transport were explored.
In general, the photochemical processes in both ionospheres behave similarly when no magnetic field is considered. The presence or absence of a magnetic field shape the production via photoionisation and loss via recombination processes in both ionospheres, especially when considering plasma transport. This makes the Terrestrial ionosphere more complex in terms of plasma transport, though similar processes were found and compared (at smaller scales) to similar processes at Mars over crustal fields.
This study has notable implications for comparative aeronomy, as a good understanding of how the ionosphere of magnetised and un-magnetised bodies compares is important for characterising planetary environments. These comparisons are useful for nearby planets such as Mars, but are perhaps even more relevant when attempting to understand the plasma environments of more distant bodies such as the icy moons of the giant planets and rocky exoplanets.
Additional file 1: Huber regressor linear fit parameters, as well as the medians and quartiles of the difference between the linear fits and observations for each Terrestrial and Martian season, SZA region, MY, and location region as a CSV file.