Issue 
J. Space Weather Space Clim.
Volume 13, 2023



Article Number  7  
Number of page(s)  24  
DOI  https://doi.org/10.1051/swsc/2023005  
Published online  04 April 2023 
Technical Article
Low latitude monthly total electron content composite correlations
^{1}
Remote Sensing Division, Naval Research Laboratory, 4555 Overlook Avenue SW, Washington, DC, 20375, USA
^{2}
Space Science Division, Naval Research Laboratory, 4555 Overlook Avenue SW, Washington, DC, 20375, USA
^{*} Corresponding author: douglas.allen@nrl.navy.mil
Received:
15
July
2022
Accepted:
24
February
2023
Spatial correlations of total electron content (TEC) variability are compared among two SAMI3 model runs and Jet Propulsion Laboratory Global Ionospheric Maps (JPL/GIM). Individual monthly correlation maps are constructed with Equatorial reference points at 12 evenly spaced longitudes and 12 universal times. TEC composite correlations (TCCs) are then calculated by averaging the individual maps, shifted zonally to synchronize local time. The TCC structures are quantified using Gaussian fits in the zonal and meridional directions. A nonzero largescale “base correlation” is found in all three datasets for 2014, a year with high solar activity. Higher base correlations generally occur in the SAMI3 runs than in JPL/GIM. The SAMI3 run driven with climatological neutral fields shows higher correlations than the run driven with neutrals from a Whole Atmosphere Community Climate Model with thermosphere–ionosphere eXtension (WACCMX) simulation. Base correlation values strongly correlate with monthly F10.7 standard deviations. Empirical orthogonal function (EOF) analyses confirm that largescale correlations are usually, although not always, related to solar forcing. Strong correlations between the Ap index and EOF modes are also observed, consistent with the geomagnetic forcing of the TEC field. The widths of the correlation structures are also examined, and these vary considerably with local time, month, and dataset. OffEquator conjugate point correlations are also calculated from each dataset and variations with the month and local time are analyzed. Analysis of TCCs for 2010, a year with low solar activity, shows that base correlations as well as correlations of the first EOF mode with F10.7 are generally weaker than in 2014.
Key words: Total electron content / Data assimilation / Correlations / Ionosphere / Modeling
© D.R. Allen et al., Published by EDP Sciences 2023
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
The goal of this paper is to explore the spatial correlations of the ionospheric state on daily timescales in order to inform ionospheric data assimilation (DA) systems on how to model the structure of the background error covariance matrix (BECM). While the full BECM contains all state variables, for simplicity we focus on the vertical total electron content (vTEC, referred to as TEC hereafter). While a number of papers have examined the correlation structure of TEC (as reviewed further below), this study is novel for two reasons. First, we compare correlation results among three different datasets. Two are yearlong model runs of Sami3 is Another Model of the Ionosphere (SAMI3, Huba et al., 2017; McDonald et al., 2018; Zawdie et al., 2020) using two different specifications of the neutral atmosphere. The third is an analysis provided by the Jet Propulsion Laboratory Global Ionospheric Maps (JPL/GIMs, Mannucci et al., 1998). By comparing SAMI3 runs with JPL/GIMs, we can assess how well the model simulates the TEC correlations observed in the “real world”. The second novel aspect of this study is the analysis of the structure of the BECM in terms of the leading empirical orthogonal functions (EOFs) that characterize the TEC variability. We aimed to understand whether the ionosphere may be dominated by a limited number of modes determining the structure of the ionospheric BECMs.
Numerous ionospheric DA efforts are now underway. Bust et al. (2001, 2004) describe the Ionospheric Data Assimilation ThreeDimensional (IDA3D) model and Mengist et al. (2019) the IDA4D, which are 3D variational (3DVar) DA systems. The majority of ionospheric DA efforts, however, involve ensemble filters. For example, several ensemble Kalman filter (EnKF) systems have been developed using the Data Assimilation Research Testbed (DART), such as Pedatella et al. (2018, 2020) with the Whole Atmosphere Community Climate model, with thermosphere–ionosphere eXtension (WACCMX), Matsuo et al. (2013), Hsu et al. (2014), Chen et al. (2016, 2019), and Chartier et al. (2016) with the National Center for Atmospheric Research (NCAR) Thermosphere–Ionosphere–Electrodynamics General Circulation Model (TIEGCM), and Morozov et al. (2013) with Global Ionosphere–Thermosphere Model (GITM). Other efforts include the Advanced Ensemble electron density (Ne) Assimilation System (AENeAS), which uses a local ensemble transform Kalman filter (LETKF) with the TIEGCM up to 600 km altitude and NeQuick above 600 km (Elvidge & Angling, 2019), Thermosphere Ionosphere Data Assimilation (TIDA), which uses an EnKF with Coupled Thermosphere, Ionosphere, Plasmasphere, and electrodynamics (CTIPe, Codrescu et al., 2018), and Global Assimilation of Ionospheric Measurements (GAIM), which uses a reduced state Kalman filter with IPM (Scherliess et al., 2004). In addition, Solomontsev et al. (2012) use an ensemble square root filter (EnSRF) with the core physics model of Khattatov et al. (2005), and Lin et al. (2015) use a Gauss–Markov Kalman filter with International Reference Ionosphere (IRI).
A defining issue with ensemble DA systems is the use of relatively small numbers of ensemble members, i.e., the number of ensemble members is far less than the dimension of the problem. Small ensemble sizes lead to spurious covariances within the BECM and subsequently to noise in the resulting state estimates. The typical way that this issue is handled is through the use of what is referred to as “localization”. Localization assumes that correlations are local to a region and therefore all covariances outside of that region are then set to zero. When covariances are truly local, as in tropospheric numerical weather prediction, this “localization” paradigm can lead to large improvements in the quality of the resulting state estimates from an EnKF. However, as we will show, the correlations in TEC can be highly nonlocal and therefore we wish to assess the plausibility of contemporary BECM localization tools for the ionosphere.
While some of the above ionospheric DA systems do not use localization (e.g., Chartier et al., 2016; Chen et al., 2016, 2019; Codrescu et al., 2018), most localize the BECM in the horizontal, and some in the vertical as well. Localization widths vary among models. For example, Elvidge & Angling (2019) use a localization region with latitude/longitude distances of 10°/20°, based on estimated ionospheric correlation lengths from McNamara (2009) and McNamara & Wilkinson (2009). Mengist et al. (2019) use a Gaussian localization function with zonal (meridional) widths of 5° (3°) in low latitudes and 10° (5°) in high latitudes. Morozov et al. (2013) use a piecewise continuous GC (Gaspari & Cohn, 1999) function with a halfwidth of 30°, while Pedatella et al. (2020) use a GC function with a halfwidth of 11.5°. Chartier et al. (2016) examined the sensitivity of ionospheric DA to localization length, comparing simulations with localization radii of 0.2, 0.5, and 1.0 radians (11.5°, 28.6°, and 57.3°) during a geomagnetic storm simulation. They found similar errors for all three radii from 00:00 to 12:00 local time (LT). The largest localization radius performed better from 12:00 to 15:00 LT, while the smallest radius performed better after 15:00 LT, possibly indicating the ensemble was not properly representing the variability in the later stages of the storm.
Forsythe et al. (2020a) derived global distributions of ionospheric correlation lengths based on IRI2016 model errors (calculated using JPL/GIM as the “truth”) and daytoday TEC variability, noticing significant differences between the two approaches, and investigated the extent of the anisotropy of the correlation length scales. Forsythe et al. (2020b) calculated vertical correlation lengths using IRI and five incoherent scatter radars (ISR) and showed that vertical correlations also vary with solar conditions and are asymmetric, generally having larger (smaller) lengths above (below) the reference altitude. Forsythe et al. (2021a) estimated correlation times, and Forsythe et al. (2021b) combined these with the correlation lengths for use in the IDA4D model. They showed that the vertical part of the covariance played the most important role and helped correct a tomographic issue relating to the assimilation of slant TEC data.
We note that correlation lengths are often based on a fixed linear correlation threshold coefficient value. For example, Klobuchar & Johanson (1977) determined the length at which the correlation coefficient drops below 0.7 as the “correlation distance” that provides a 29% uncertainty reduction in predicting one station using data from another. Shim et al. (2008) and McNamara (2009) also used a 0.7 cutoff to define correlation lengths, while Yue et al. (2007) and Liu et al. (2018) used a 0.75 cutoff, and Forsythe et al. (2020a) used a 0.8 cutoff. While this approach provides an objective method for estimating correlation lengths, care needs to be taken in applying quoted lengths for DA purposes. Many DA algorithms make use of Gaussian covariance models. The proper way to fit a Gaussian function to a correlation field is to find the value of the correlation equal to e^{−0.5} ≈ 0.607. Therefore, typical correlation length scales in the literature may be too short to be applied to Gaussian covariance modeling.
Another important use of correlation length scales in DA is the construction of the localization function for ensemble schemes. When the goal is localization then one is not trying to fit a Gaussian to the correlation structure. The goal here is rather to define a function that, when multiplied against the covariances within the BECM, will keep the covariance structure near the diagonal of the BECM intact while also setting to zero farfield noise in the BECM (see Morzfeld & Hodyss, 2023, for more detailed discussion). This requires a length scale longer than the true length scale of the data. Because the width of the localization scheme depends not only on the physical correlations but also on the ensemble size this localization length scale is typically obtained through adhoc tuning of the DA system. Additionally, it has been found that using too narrow of localization can lead to poor representations of multivariate relationships between variables and therefore poor forecasts from the ensemble (Kepert, 2009; Greybush et al., 2011).
Ideally, the BECM represents the correlation of errors across variables and across different spatial regions. In this study, we focus on the spatial correlations of vertical TEC. A number of studies have examined TEC correlations on different timescales, and many of these have also attempted to determine the physical cause of the correlations. Three main drivers of TEC variability have been examined: solar variability, geomagnetic variability, and meteorological processes (e.g., Forbes et al., 2000; Rishbeth & Mendillo, 2001; AraujoPradere et al., 2004, 2005; Fang et al., 2013, 2018). As explained by Fang et al. (2018), the contributions of these contributors to the total TEC or N_{m}F_{2} variability strongly depend on local time, latitude, and altitude. Fang et al. (2018) used the coupled Whole Atmosphere Model (WAM) and Global Ionosphere Plasmasphere (GIP) to simulate ionospheric variability during June – July 2012, a time with strong variations in solar and geomagnetic activity. Comparing three simulations that isolated lower atmosphere processes, geomagnetic activity, and solar forcing they found that at low latitudes solar activity dominates the daytime absolute TEC variability, while the contributions from the lower atmosphere and geomagnetic activity are approximately equal. However, when examined in terms of relative contributions to the global TEC variability, geomagnetic activity was the main contributor, followed by solar activity, and then lower atmospheric perturbations. In periods of low solar and geomagnetic activity, the lower atmosphere processes (tides, gravity waves, and planetary waves) can play a much larger role. For example, it has been shown that lower atmospheric forcing associated with sudden stratospheric warmings can significantly impact the ionosphere during geomagnetically quiet conditions (Siddiqui et al., 2021).
Observational studies have also confirmed that ionospheric variability depends on local time, season, latitude, solar proton flux, solar wind, geomagnetic activity, and terrestrial weather (Rush & Gibbs, 1973; Forbes et al., 2000; Rishbeth & Mendillo, 2001). Early studies of ionospheric spatial correlations have been performed using ionosonde and Global Positioning System (GPS) observations (e.g., Klobuchar & Johanson, 1977; Nisbet et al., 1981; Gail et al., 1993; Bust et al., 2001). More recent comprehensive analyses have been performed using global datasets. Yue et al. (2007) and Forsythe et al. (2020a) used JPL/GIM to analyze variations of horizontal correlation lengths with local time (LT), magnetic latitude (Mlat), and season. Shim et al. (2008) used over 1000 groundbased receivers to analyze spatial correlation variations with LT, Mlat, and season. Liu & Chen (2009) examined how TEC varies with different solar indices. McNamara (2009) and McNamara & Wilkinson (2009) used ionosondes to calculate spatial correlations for different values of the Ap magnetic index. Liu et al. (2018) used JPL/GIM with detailed fitting analyses to calculate length scales as a function of LT, month, and Mlat. Several studies also examined geomagnetic conjugate correlations (e.g., Yue et al., 2007; Liu et al., 2018). These are interhemispheric correlations associated with connecting magnetic field lines. The conjugate correlations pose a problem for DA since they do not behave like typical monotonic correlation functions.
This paper will provide a detailed analysis of TEC correlations to complement previous studies and provide guidance for ionospheric DA. We note that the correlations presented here represent TEC variability about the timemean (monthlong sequences are used), which we assume reasonably represents the correlations for the monthlyvarying “climatological covariance” that could be used in a DA system that does not have an ensemble capability, or could be combined with an ensemble in a hybrid 4DVar system in order to improve the rank of the BECM. An alternate approach would be to compile daytoday differences in TEC, ΔTEC, which has often been used for examining ionospheric variability. Assuming the TEC values are the true values, ΔTEC approximates the error that would be obtained using a 24h “persistence” forecast model, in which the next day’s forecast is identical to today’s analysis. Forsythe et al. (2020a) showed that horizontal correlation lengths calculated using ΔTEC were significantly different from those calculated using estimated model error from the difference between IRI and JPL/GIM maps. Our approach provides a new method to approximate TEC correlations that can be compared to other approaches.
The paper is organized as follows. Section 2 provides descriptions of the models and data used in this study. Section 3 details the TEC correlation composite calculation and examines the base correlations for the three datasets. Section 4 uses EOFs to determine the dominant modes of TEC variability and relates these to the solar and geomagnetic indices. This section also quantifies correlation widths and examines the interhemispheric conjugate correlations. Section 5 provides a summary and conclusions. We mainly focus on 2014 in this paper, a year with high solar activity. However, for comparison, we also show in the Supplementary Material some results for 2010, a year with low solar activity.
2 Model and data description
2.1 SAMI3
SAMI3 is a physicsbased model of the ionosphere that simulates the transport and chemistry of seven ion species and solves for the electron and ion temperature for three species (Huba et al., 2000). The momentum equation is split into motions along and across the magnetic field. In the configuration used for this study, the Earth’s magnetic field is specified by the International Geomagnetic Reference Field (IGRF) in the magnetic apex coordinate system (Richmond, 1995; Laundal & Richmond, 2016). Details on the solver can be obtained from previous studies (see Zawdie et al., 2020 and references therein). We compare runs using two different SAMI3 configurations for the thermosphere. For “SAMI3/climo,” the neutral temperature and composition are from NRLMSISE00 (Picone et al., 2002) and the thermospheric winds are from the Horizontal Wind Model (HWM14) (Drob et al., 2015). For “SAMI3/WACCMX,” the neutral atmosphere is specified using WACCMX temperature, composition, and winds. As explained in McDonald et al. (2018) and Zawdie et al. (2020), for this experiment, WACCMX is run in a specified dynamics configuration in which the lower boundary is forced by the Navy Global Environmental ModelHigh Altitude (NAVGEMHA, Eckermann et al., 2018). The solar forcing of SAMI3 uses a flux term that includes the average of the daily F10.7 index and the 81day mean of the daily F10.7 index (Huba et al., 2000). In addition, for SAMI3/climo the F10.7 affects the neutral densities and temperatures via NRLMSISE00. The geomagnetic forcing is through the Ap index, which in SAMI3/climo affects both the neutral’s density and temperature via NRLMSISE00 and the neutral HWM14 winds. For the SAMI3/WACCMX used here, there is no explicit use of the Ap index, although geomagnetic variations can implicitly affect SAMI3 via its influence on the WACCMX neutrals (Liu et al., 2010).
2.2 JPL Global Ionospheric Maps
To compare with the SAMI3 forecasts, we use the 2h JPL/GIMs of TEC. As described in Mannucci et al. (1998) and Iijima et al. (1999), the approach uses a Kalman filter to retrieve vertical TEC maps from global GPS measurements. The JPL/GIMs are generated from a global network of 200 GPS receivers at a longitude/latitude resolution of 5°/2.5°. The JPL/GIMs were downloaded from the Crustal Dynamics Data Information System (CDDIS) website (https://cddis.nasa.gov/archive/gnss/products/ionex). Further details are provided in Emmert et al. (2017).
2.3 Comparison of the datasets at one location
While a complete comparison of the datasets is beyond the scope of this paper, we provide in Figure 1 a sample comparison that shows the daily TEC values at one location (latitude = 0°, longitude = 0°) over 2014 at four times of the day (0000, 0600, 1200, 1800 UTC). All datasets show a strong diurnal variation in TEC, minimizing at 0600 UTC and maximizing at 1200 or 1800 UTC. They also show a semiannual cycle with maxima in March–April and October–November, coinciding with solar variations at the Equator, in addition to a 27day cycle associated with the solar rotation period. On daily time scales, the SAMI3/climo is smoother than the other datasets (e.g., in June/July at 1800 UTC), while SAMI3/WACCMX shows enhanced daily variability.
Fig. 1 Daily TEC values at latitude = 0°, longitude = 0° for 2014 from SAMI3/climo (black), SAMI3/WACCMX (red), and JPL/GIM (blue). Panels are for different times of day: (a) 0000 UTC, (b) 0600 UTC, (c) 1200 UTC, and (d) 1800 UTC. 
The SAMI3/WACCMX TEC is generally biased high with respect to the other datasets at 0000, 1200, and 1800 UTC. While a correction to the SAMI3/WACCMX TEC could be applied (McDonald et al., 2018 scaled the electron density by a factor of 0.7), the TEC correlations used in this paper are not impacted by a constant scaling, so we do not apply a correction for this paper. While the 27day cycles at 1800 UTC appear to be similar between SAMI3/climo and JPL/GIM, the amplitude of the cycle is much larger in SAMI3/WACCMX, suggesting that there is enhanced responsiveness to solar variations. The SAMI3/climo more closely aligns with the JPL/GIM data in these time series, although there are some biases over the annual cycle. For example, at 1800 UTC, SAMI3/climo is biased high from January to April and October to December and is biased low from May to September. In the results shown below, we will focus on monthly correlations, in which the mean values are removed from the time series. Therefore, the biases among the datasets are largely removed. This allows us to examine and compare the geographic interdependence of the TEC fields for each dataset. The SAMI3/climo and JPL/GIM TEC values at latitude = 0°, longitude = 0° for 2010 are shown in Supplementary Figure S1 (SAMI3/WACCMX data are not included in the analyses of 2010). The TEC daily values are much smaller and the variability is much weaker than in 2014 due to weaker solar activity, as shown below. JPL/GIM TEC is generally higher than SAMI3/climo TEC in 2010.
2.4 Solar and geomagnetic indices
To gain perspective on the variability in the solar and magnetic indices, Figures 2a and 2b shows the daily F10.7 and Ap values for the year 2014. There are wide swings in F10.7 over some months, while others show relatively steady F10.7. To quantify the variability, Figures 2c and 2d shows the monthly standard deviations (σ) of F10.7 and Ap. The in January is 33, compared to 11 and 7.5 in February and March, respectively. is also large in July (42) and October (37). The Ap index shows very large values in February and very small values in March and July. Also plotted in Figures 2c and 2d is the monthly correlation coefficients of the 1200 UTC empirical orthogonal function amplitude with the F10.7 and Ap indices (red lines), which we will examine later. Supplemental Figure S2 shows the indices for 2010. While the F10.7 values are much smaller and have much weaker variability, the Ap index shows some months with strong variability, particularly in April, May, and August 2010.
Fig. 2 (a) Daily values of F10.7 over the year 2014. (b) Daily values of Ap index over 2014. (c) Monthly standard deviation of F10.7 (black line) for 2014 along with the correlation of the 1200 UTC EOF1 amplitude with F10.7 (red line). (d) Monthly standard deviation of Ap index (black line) for 2014 along with the correlation of the SAMI3/WACCMX 1200 UTC EOF1 with Ap (red line). 
3 TEC correlation composite analysis
3.1 Interpolation to regular and magnetic grid
In the results described below, the SAMI3 TEC maps are linearly interpolated to the JPL/GIM grid, which has a longitude × latitude resolution of 5° × 2.5°, to allow easier comparison. Each map is also interpolated to a 2D (constant in altitude) magnetic grid with a resolution of 5.0° magnetic longitude (Mlon) × 2.5° magnetic latitude (Mlat). The magnetic grid is obtained by interpolating the SAMI3 “Apex” magnetic grid line locations to obtain the latitude/longitude positions where these grid lines cross the constant altitude of 110 km, as used in Figure 2 of Richmond (1995). We then interpolate horizontally to obtain the latitude/longitude coordinates associated with each Mlat/Mlon pair on the magnetic grid. In order to avoid polar processes, we limit both latitude and Mlat to ±60°. The motivation to exclude the poles is both due to limitations in the SAMI3 run settings used in this study as well as larger errors in the JPL/GIM products at high latitudes due to less coverage (Mannucci et al., 1998). Yue et al. (2007) similarly limited the data range to Mlat of ±60°. Both grids are displayed in Figure 3. Note that Mlon is defined as zero for the line that passes through the geomagnetic pole (longitude = 71° W, latitude = 79° N for this Apex grid for 2014). This position on the magnetic grid is indicated by a red line in Figure 2b. To facilitate comparisons made using the two grids, this longitude shift will need to be considered, since the local time of longitude = 0° does not coincide with local time Mlon = 0°. We also highlight in Figure 3 the 12 longitudes and Mlons that are used for the composite calculations discussed below.
Fig. 3 (a) The geographic grid used for this study, which is at the JPL/GIM longitude × latitude resolution of 5° × 2.5°. The Greenwich Meridian is highlighted by a red line. The blue lines indicate the 12 reference longitudes used in the analysis. (b) The geomagnetic grid used for this study is at Mlon × Mlat resolution of 5° × 2.5°. resolution. The location of magnetic longitude of 0° is indicated by the red line. The blue lines indicate the 12 reference magnetic longitudes used in the analysis. 
3.2 Individual correlation maps
The main approach used in this paper is to calculate TEC composite correlations by averaging individual correlation maps taken at different times that are shifted in the longitude to synchronize the local time. The motivation for this approach is (1) to account for the significant diurnal cycle in TEC due to solar production in the daytime and nighttime recombination, and (2) to reduce spurious correlations by combining multiple maps. The mathematical formalism is provided below for clarity.
First, we define a linear correlation map that is a function of longitude and latitude (denoted by λ_{i} and φ_{j}) for a selected reference point, , and a reference coordinated universal time (UTC), in hours. We note here that (λ_{i}, φ_{j}) may refer to either geographic coordinates or magnetic coordinates. The individual correlation maps are produced by first compiling TEC maps over a given time range (one month is used in our study), with the number of available days indicated by D, and the day of the year indicated by d. We note here that we do not bin the data by geomagnetic or solar activity, so some of the variations by month will be due to variations in these parameters. Including these variations allows us to assess their impact on the correlations. This is important since there is uncertainty in geomagnetic and solar activity in real DA systems that needs to be modeled in the BECM if it is to represent the true error in the system. Because each day of the month has different solar and geomagnetic activity, a monthly dataset approximates an ensemble with uncertainty in solar and geomagnetic activity. In Section 4, we will filter the TEC data using empirical orthogonal functions in an attempt to examine to what degree the variability in each month can be attributed to solar and geomagnetic activity.
We next form an array of indexed TEC values (q) at a fixed hour (h, in UTC) as follows:
A global map of the linear correlation coefficient, r, for a given reference point is made by calculating the correlations for all longitudes (total of N_{lon} = 72) and latitudes (total of N_{lat} = 49):
where i = 1, 2, …, N_{lon} and j = 1, 2, …, N_{lat}, overbar indicates mean over all days and σ is the standard deviation over all days. Sample correlation maps are provided in Figure 4, which are calculated using a fixed time of 1200 UTC and using 30 days covering 2–31 January 2014 (1 January is left off due to SAMI3’s use of a “cold start,” i.e., initialized with electron density specified by a Chapman layer). The reference point is indicated with a black dot on each panel and is located at the Equator, but at different longitudes in each frame. Also plotted in Figure 4 are several lines of solar zenith angle (SZA) for reference. Figures 4a, 4e and 4i show the reference point at longitude = −180° (0000 LT). There are high correlations (greater than ~0.6, indicated by red/orange contours) within the vicinity of the reference point. The SAMI3/climo also has high correlations at locations remote from the reference point, while SAMI3/WACCMX and JPL/GIM have weaker longdistance correlations. As the reference point moves into the sunrise, Figures 4b, 4f and 4j, the SAMI3/climo data shows high correlations spreading both along the east–west and in the north–south directions, starting to fill the sunlit region. SAMI3/WACCMX also has enhanced correlations, while JPL/GIM correlations are weaker. As the reference point moves to noon (Figs. 4c, 4g and 4k), the high correlations span nearly the whole sunlit region for SAMI3 data, but only partially for JPL/GIM. As the reference point moves into the evening (Figs. 4d, 4h and 4j), these large correlations remain for SAMI3/climo, but die down and are more localized for SAMI3/WACCMX and JPL/GIM. These correlation maps, therefore, show a strong variability with longitude, latitude, local time, and dataset. There are also apparently spurious correlations at large distances. This is unavoidable due to the small sample size (30) relative to the size of the state (72 × 49 = 3528). As shown below, composite averages help to reduce this noise.
Fig. 4 TEC correlation maps for 2–31 January 2014 at 1200 UTC for SAMI3/climo, SAMI3/WACCMX, and JPL/GIM. The reference point is marked by a black dot. It is located on the Equator, and the longitude varies with the panel (−180°, −90°, 0°, and 90° are shown). The black lines are SZA for 16 January 2014 at 10°, 20°, …, 90°. 
Figure 5 provides similar correlation maps, but in the magnetic coordinates at 1200 UTC. As mentioned earlier, the magnetic longitudes are shifted relative to regular longitudes, so the local times of Mlon = 0° is similarly shifted. Here we calculate the approximate local time of the magnetic grid by a shift of (−71°)/(15°/h) = −4.7 h, and these are indicated in each panel in Figure 5. The correlations using the magnetic grid show a similar diurnal cycle to the correlations using the geographic grid, with very largescale correlations occurring in the SAMI3 data as the reference point enters the sunlit region and more compact correlations when the reference point is in the night. In fact, assuming the exact same point on the globe is used as the reference point, there would be a onetoone correspondence between correlations using the two grids. So there is not necessarily an advantage to going to the magnetic grid for a single correlation map. However, the composites can show important differences between the two reference frames.
Fig. 5 TEC correlation maps for 2–31 January 2014 at 1200 UTC on the magnetic grid for SAMI3/climo, SAMI3/WACCMX, and JPL/GIM. The reference point is marked by a black dot. It is located on the magnetic Equator, and the Mlon varies with panel (0°, 90°, 180°, and 270° are shown). The black lines are SZA on 16 January 2014 at 10°, 20°, …, 90°. 
3.3 Composite correlation maps
The next step is to average the individual maps at different UTCs but shifted horizontally so that the local times of the maps are synchronized. Defining the local time (LT hereafter) at a given longitude as LT_{i}, we can write the dependence of this variable on h and longitude:
We use a similar relationship for the magnetic local time (MLT), which accounts for the geographic longitude of the north magnetic pole λ_{N}:
where λ_{i} refers to magnetic longitude. As explained in Laundal & Richmond (2016), this definition is simplified, since it does not take into consideration the magnetic longitude of the subsolar point. As a further simplification, with a magnetic pole at a longitude of −71° we note that the correction term in equation (4) is not an integer, but rather it has the value λ_{N}/15 = −71°/(15°/h) = −4.73 h. For our composites to be on evenlyspaced 2h intervals, we assume that λ_{N} = −60°, so that λ_{N}/15 = −4.00 h. This causes a 44 min time error in the MLT. For the qualitative purposes used in this paper, this simplification is adequate.
To calculate the composite correlation for a given fixed latitude (or Mlat), we average over individual correlation maps to reduce spurious noise. In this average, we want the local time of the reference points (the “reference local time”, LT_{0}) to be the same for all maps, and we want to shift in the longitude to line up all of the reference points. For this composite calculation, it is helpful to specify the data intervals in longitude and time. We choose a fixed time spacing of Ω = 2 h. Using the longitude spacing for the TEC maps, 5°, for each time shift of Ω we need to shift the longitude index by grid points. So we want to use individual correlation maps at 12 reference points at longitudes of 0°, 30°, 60°, …, and 330°. The TEC correlation composite (TCC hereafter) can then be calculated as follows, looping over reference points, and shifting in longitude in order to line up the reference points. Here i and j are the longitude and latitude indices for the TCC map, and i_{0} and j_{0} are the longitude and latitude indices of the individual reference points:
In equation (5), h is the universal time that corresponds to the local time at the given reference point as detailed in equation (6). Note that when h becomes negative, we add 24 h:
This procedure provides TCC values at all longitudes and latitudes. We can convert the resulting longitude index to a local time index by LT(i) = LT_{0} − i∆/15. This results in composite maps as a function of local time and latitude, . These composites can be calculated at all latitudes, but we will focus on the Equator and 20° N. They are also calculated for the 12 reference times: LT_{0} = 0, 2, 4, …, 22 h.
Figure 6 shows Equatorial reference point TCCs maps for 2–31 January 2014 on the regular grid. The progression of LT_{0} in Figure 6 is similar to the progression of longitudes in the individual correlation maps in Figure 4. The TCCs show a similar diurnal structure to the individual correlation maps, but are generally smoother and show a reduction of the longdistance spurious correlations. The TCC maps for 2–31 January 2010, provided in Supplemental Figure S3, show generally weaker correlations for SAMI3/climo and JPL/GIM than in January 2014. This is consistent with the weaker variability of the solar flux seen in Supplemental Figure S2. The SAMI3/climo TCC maps in both years exhibit very strong global correlations, particularly throughout the sunlit region. This suggests that there are only a small number of dominant physical forcings in SAMI3/climo. Compared to the SAMI3 TCCs, JPL/GIM show correlations that are generally more localized relative to the reference point.
Fig. 6 SAMI3/climo, SAMI3/WACCMX, and JPL/GIM TEC correlation composites for 2–31 January 2014 as a function of latitude and LT. The reference point is marked by a black dot. It is located on the Equator, and the reference LT varies with panel (0000, 0600, 1200, and 1800 are shown). Thick black lines indicate a correlation of 0.5, 0.7, and 0.9. Thin black lines are SZA on 16 January 2014 at 10°, 20°, …, 90°. 
While there is some diurnal variability (to be examined more closely below), the JPL/GIM does not exhibit the very large correlations filling the sunlit region as seen with SAMI3. This is consistent with Shim et al. (2008), who showed using groundbased GPS observations that TEC changes over most of the illuminated disk are largely uncorrelated, suggesting that “daytoday changes in the solar illumination only add a statistically insignificant contribution to the daytoday changes in TEC.” There is an indication of spreading correlations in the meridional direction, particularly at 1800 LT_{0}, but otherwise, the correlations remain quite local. One possible reason for the reduction of correlations in JPL/GIM is noise in the TEC observations themselves, which could mask some of the real variability. However, it is also likely that SAMI3 overestimates the correlations due to limitations in producing smallerscale ionospheric phenomena.
3.4 Base correlation calculations
In order to characterize the spatial variation of the composites, we fit the TCCs in the zonal and meridional directions with Gaussian curves that have a fixed center at a maximum of 1.0 at the reference LT (or MLT) and reference latitude (or Mlat). With x referring to a given coordinate (longitude, latitude, Mlon, or Mlat), each of the fits results in two constants, a_{0} (which we refer to as “base”), and a_{1} (which we refer to as “width”) that describe the variations,
Note that if a_{0} = 0, then a_{1} is the halfwidth associated with the correlation value of e^{−0.5} = 0.607. In typical meteorological correlations, the base would be zero, indicating that a single phenomenon, with a single physical correlation scale, is present. But as seen in Figure 6, the correlations generally do not go to zero at a great distance from the reference point. This is consistent with at least two physical phenomena at work, one with a very long, globally correlated structure and another with more a local structure. The Equatorial LT line plots associated with the TCCs on the regular grid (black lines) are shown in Figure 7 for 2–31 January 2014, along with the Gaussian line fits (blue dashed lines). Values of the base and width are indicated on each plot. The base values range from 0.19 to 0.83 (unitless), indicating strong largescale correlations. Consistent with the TCC maps, the SAMI3/climo shows the strongest base correlations (range from 0.78 to 0.83), followed by SAMI3/WACCMX (0.43–0.53), and then JPL/GIM (0.19–0.28). The widths vary considerably with dataset and local time, ranging from 1.40 h to 8.49 h. In Section 4 we will examine a more robust calculation of the widths after filtering out the largescale base correlations. Supplemental Figure S4 shows TCC line plots for 2–31 January 2020. While the SAMI3/climo base correlations are still quite high (0.56–0.65), the JPL/GIM correlations have weakened considerably (0.03–0.15).
Fig. 7 SAMI3/climo, SAMI3/WACCMX, and JPL/GIM line plots of TCC versus relative LT (black) for 2–31 January 2014 at four reference local times (0000, 0600, 1200, and 1800 LT). The reference latitude is the Equator. The Gaussian fits are shown by the dashed blue lines, with the two fitted parameters shown at the bottom of each plot. Base is unitless and width is in hours LT. 
If the base correlations are largely due to F10.7 variability, then we would expect the monthly base values to be correlated with the monthly . To test this, Figure 8 plots these two quantities as a function of the month in 2014 for SAMI3/WACCMX. Since the base depends on reference local time, there are 12 lines plotted, and colorcoded by LT. As seen in the time series for base calculated using longitude (Fig. 8a), there is a strong relationship between base and , and linear correlations between the two quantities range from 0.76 to 0.83. This confirms that the base correlations in SAMI3/WACCMX are largely due to solar variations driving the photoproduction of electrons and ions. Similar variability is seen using latitude (Fig. 8b). Note that there are two points for 0800 LT with negative correlations in July and November that are spurious. The correlations using Mlon (Fig. 8c) tend to be higher than for longitude, and the correlations with Mlat (Fig. 8d) also show correspondence with . We show in Supplemental Figure S5 the monthly longitudinal and Mlat base values for 2014 for the three datasets and for each local time (red lines), along with the mean values over all local times (blue lines), and (black lines). Like SAMI3/WACCMX, the SAMI3/climo and JPL/GIM base values appear to follow fairly well, with the exception of a few months for JPL/GIM. The base values for both Lon and Mlat are generally largest in SAMI3/climo and smallest in JPL/GIM, while SAMI3/WACCMX shows values in the middle.
Fig. 8 Base values versus month from SAMI3/WACCMX for 2014 for the coordinate (a) longitude, (b) latitude, (c) Mlon, (d), Mlat. Colors indicate local time as shown at the top. Also plotted in each panel is the monthly . 
The base correlation dependence on local time is examined in Figure 9 for SAMI3/WACCMX for 2014. For longitude and Mlon (panels a and c), there is a tendency towards maxima around 0800 LT and 1800–2000 LT, corresponding to morning and evening hours, with minima around 0400 and 1400 LT. For latitude and Mlat (panels b and d), the base values generally show a strong diurnal cycle, with low values occurring from midnight to 0600 LT, a sharp increase from 0600 to 1000 LT, and a decrease from 1800 to 2200 LT. This LT time variation resembles the variation of the correlation scale seen in Liu et al. (2018), who used JPL/GIM data to determine correlation scales with a correlation value cutoff method. Their results for April 2008, for example, have lowlatitude correlations that rise sharply from 0600 to 1200 LT, but then fall off more slowly from 1200 to 2200 LT. They also showed strong seasonal variability in 2008 with higher zonal correlations in spring and autumn than in summer and winter, as well as strong interannual variability. The diurnal variability of the base values for all three datasets is shown in Supplemental Figure S6. The semidiurnal cycle in the longitudinal base values is seen in SAMI3/climo and SAMI3/WACCMX, but not in JPL/GIM. A strong diurnal cycle is observed in SAMI3/climo and SAMI3/WACCMX Mlat base values as well as JPL/GIM, although the amplitude of the diurnal cycle is generally smaller.
Fig. 9 Base values versus local time from SAMI3/WACCMX for 2014 for the coordinate (a) longitude, (b) latitude, (c) Mlon, (d), Mlat. Colors indicate month as shown at the top. 
4 Empirical orthogonal function decomposition
4.1 Mathematical formulation and sample EOF plots
The Gaussian fits suggest that TCCs may often be described by a globalscale base correlation and a local correlation. In an attempt to separate these two features, we apply an empirical orthogonal function (EOF) analysis to the data in order to identify and remove the largest sources of variability. As we will show below, this is often attributed to largescale solar forcing. The EOF analysis is equivalent to singular value decomposition, or principal component analysis, in which a matrix of ncolumns by mrows is decomposed as the product of three matrices:
Here A is a matrix of size N_{days} × N_{state} (number of days × size of the state) for which each column contains a vector of TEC data at fixed UTC over a given date range, with the time mean of each state element removed. Here N_{state} = N_{lon} × N_{lat} = 3528. U and V are the left and right singular vectors, and S is a diagonal matrix of eigenvalues. From this decomposition, we can determine the spatial structures of the EOFs and the amount of variance contained in each EOF “mode”. We also calculate the amplitude of each mode as a function of time by projecting the TEC maps onto the EOF of each mode (i.e., a dot product of the two vectors). By correlating this amplitude with forcing indices (F10.7, Ap), we will attempt to provide insight into the physical mechanisms involved in different modes. While the approach is appealing, care must be taken to avoid over interpretation of individual modes, since forcing mechanisms are generally not orthogonal, while the eigenfunctions of the EOFs are orthogonal by design (Monahan et al., 2009). In plotting the EOFs and the amplitudes, we divide each EOF by the global maximum absolute value, and we multiply the amplitude by the same value. The plotted EOFs are then normalized to a maximum absolute value of 1 and the amplitudes (in TECU) have magnitudes similar to the global maximum TEC anomalies.
Previous studies have also looked at TEC EOFs. For example, Chen et al. (2015) examined EOFs using monthly TEC averages over North America. They found the first mode constitutes 97.5% of the total variance and represents the spatial and diurnal variations. Talaat & Shu (2016) examined EOFs using eleven years of daily averaged TEC data and found the first (first four) EOFs to explain 89% (98%) of the total variance. Zhong et al. (2019) examined ~10 years of topside ionospheric and plasmaspheric TEC data and found the first (first five) EOFs to account for 95.5% (98.8%) of the total variance. Vaishnav et al. (2019) examined EOF analysis of 19 years of global TEC maps and found the first three EOFs explain 99% of the total variance. Zhang et al. (2021) examined the EOF decomposition of monthly electron density using 11 years of radio occultation profiles, looking at the variations with the season as well as the solar cycle, and they found the first five EOFs to account for 98% of the variance. These studies show a strong structure of TEC variability by a few modes. One unique aspect of our analyses is that we group the data by month, rather than long time periods, so that we avoid the largescale variations associated with the solar cycle. Another is that we fix the UTC time, which eliminates the obviously dominant diurnal cycle. After calculating the EOFs, we can remove one or more of the modes from the TEC time series and recompute the TCCs to ascertain if the largescale structures have been reduced.
For illustration, the first three EOFs for the 2–31 January 2014 SAMI3/WACCMX time series at 1200 UTC are provided in Figure 10. The first mode accounts for 60.2% of the total variance in the time series, while the second and third modes account for 9.1% and 6.8%, respectively. The spatial pattern in Figure 10a shows a maximum over Africa and generally high values over most of the sunlit region. There is also a broad “tongue” of high value that extends zonally across southern Asia and the Pacific Ocean. The amplitude of this mode can be calculated by projecting the pattern in Figure 10a onto individual daily TEC maps. The result is shown with a black line in Figure 10d, which has large values in the first few days of January, sloping down to a minimum on day 17 (15 January), and sloping upward the rest of the month. Also plotted in Figure 10d are the daily F10.7 solar index (red line), linearly scaled to provide the best fit to the EOF1 amplitude, along with a similarly scaled Ap index (blue line). There is a high degree of correlation (0.967) between EOF1 amplitude and F10.7, suggesting that it is the solar forcing in SAMI3 that is largely responsible for this mode. We also correlate EOF1 amplitude with the Ap daily index and obtain a low correlation value of 0.115. The second and third EOFs are plotted in Figures 10b and 10c, and these show generally smallerscale structures that are both positive and negative. Neither F10.7 nor Ap shows a strong correlation with the amplitudes of the second and third modes, as seen in Figures 10e and 10f. That the F10.7 correlates well with TEC variations is consistent with other studies such as Liu & Chen (2009), who examined the correlation of JPL/GIM TEC with solar indices using quadratic regression. They found linear coefficients as high as 0.9, with significant diurnal and seasonal variability. Bergot et al. (2013) found similarly high correlations between TEC (from the Centre for Orbit Determination in Europe, CODE) and solar indices and derived an empirical model for TEC prediction based on F10.7. Vaishnav et al. (2019) provide a detailed EOF analysis of global TEC and found strong correlations of the first principal component with F10.7, while the third component correlated well with Ap and Kp indices. They also examined correlations with a number of other solar indices such as Mg II and He II.
Fig. 10 SAMI3/WACCMX EOFs for 2–31 January 2014 (1200 UTC). (a) – (c) The spatial pattern of the EOFs 1, 2, and 3 with the color label at the bottom. Thin black lines are SZA on 16 January 2014 at 10°, 20°, …, 90°. (d) – (f) The amplitudes of each EOF (black lines) along with the indices F10.7 (red lines) and Ap (blue lines), linearly scaled to provide the best fit to each of the amplitudes. The correlation coefficients of amplitudes with respect to the indices provide for each EOF at the bottom, along with the significance pvalue, where p < 0.050 is 95% confidence level. 
Figure 11 shows the first three EOFs for SAMI3/climo at 1200 UTC for 2–31 January 2014. The leading EOF (Fig. 11a), which captures 86% of the variability, is also strongly correlated (r = 0.962) with daily F10.7, and the spatial structure is similar to the EOF1 for SAMI3/WACCMX (Fig. 10a). Much of the remaining SAMI3/climo variability is split between EOF2 and EOF3, which account for 7.6% and 2.9%, respectively, of the variability. These two EOFs correlate fairly well with the Ap index (r = 0.670 and r = 0.627, respectively). As the magnetic indices are major drivers of the NRLMSISE00 and HWM14 variability, we infer that it is the neutral densities and winds that are largely responsible for the second and third EOFs. As verified in offline calculations (not shown), a pure longitudinal westward propagating sine wave representing a semidiurnal tide would be split into two EOFs in our approach. So it is likely that the same forcing is largely driving EOFs 2 and 3 in SAMI3/climo. We note that the first three EOFs combine to determine 97% of the TEC variability for SAMI3/climo for January 2014.
The EOFs for JPL/GIM at 1200 UTC for 2–31 January 2014 are shown in Figure 12. Unlike the SAMI3 EOFs, the JPL/GIM data show that the amplitude of the first EOF is not strongly correlated with F10.7 (r = 0.056), suggesting that there is another physical mode of variability that is dominating. The second EOF is more highly correlated with F10.7 (r = 0.703), and this mode accounts for only 15.3% of the variability. The JPL/GIM EOFs at 1800 UTC (Fig. 13), however, do show the first EOF highly correlated with F10.7 (r = 0.788), accounting for 31.7% of the variability. More detailed observational and modeling analyses would be necessary to understand the exact physical causes of each EOF.
4.2 Correlation width calculations
We next remove the leading mode from the TEC time series by subtracting the product of the EOF1 amplitude and the spatial EOF1 for each day. Then we repeat the TCC calculations with the filtered data. The resulting filtered TCCs for 2–31 January 2014 are shown in Figure 14. When compared with the original unfiltered TCCs in Figure 6, it is clear that the filtered TCCs have a much more compact structure and less diurnal variability. The baseline values are also much smaller, as seen in fits to the longitudinal TCC line plots (Supplemental Fig. S7), which show the base values now range from 0.11 to 0.35 for SAMI3/climo, −0.05 to 0.01 for SAMI3/WACCMX, and −0.07 to 0.04 for JPL/GIM. The EOF filtering apparently removed much of the largescale correlation structure associated with F10.7 forcing in SAMI3. We note that while F10.7 is usually the dominant forcing, as seen by inspection of individual EOF1 amplitude/F10.7 correlations, there are months when the F10.7 influence is reduced. We will examine this in more detail in Section 4.4. While Figure 14 shows more localized correlations than the TCCs made with unfiltered TEC, there are still largescale correlations for SAMI3/climo, including negative correlations in the southern extratropics. These are likely caused by largescale correlation structures in the tidal modes of the HWM that translate to the TEC variability.
Fig. 14 As in Figure 6, except the leading EOF has been removed from the TEC time series before calculating the TCCs. 
We next examine the widths calculated after filtering out EOF1. The variations with month for 2014 for SAMI3/WACCMX are shown in Figure 15 for all four coordinates. We note that all longitude widths are calculated at the Equator, and we show the widths here in degrees of great circle distance (GCD) rather than hours (1 h = 15° GCD at the Equator). There does not seem to be any correlation between these widths with , since the EOF1 has been removed from the time series, so we have not included in these plots. We also looked for correlations with Ap but did not find any obvious relationships. It is likely that the variability is due to different dynamical conditions in the different months within the WACCMX experiment. At some times (2000, 2200 LT), there is very little variability from monthtomonth, while at other times (0000, 0200 LT) the variability is strong. Longitudinal widths range from ~20 to 80° and are generally larger than latitudinal widths, which range from ~10 to 40°. Mlon and Mlat widths show similar variations to their geographic counterparts but are generally slightly smaller. The longitudinal and Mlat widths (after removing the EOF1) as a function of the month for all three datasets are shown in Supplemental Figure S8. Red lines indicate different local times, and the blue lines are the means over all local times. The SAMI3/WACCMX mean longitudinal widths range from ~30° to 45°. The JPL/GIM mean widths are smaller (~30 to 40°) and show less monthly variation. SAMI3/climo mean widths are significantly larger in all months, reaching ~100° in March. This is expected due to weaker longitudinal variability in the HWM. In the Mlat direction, SAMI3/WACCMX and JPL/GIM have similar mean widths of ~15–20°, while SAMI3/climo shows significantly larger widths, particularly in March and April.
Fig. 15 Width values versus month of 2014 from SAMI3/WACCMX (with EOF1 removed) for the coordinate (a) longitude, (b) latitude, (c) Mlon, (d), Mlat. Colors indicate local time as shown at the top. All widths are in great circle degrees. 
As seen in Figure 16, there is a large variation of width with LT for the SAMI3/WACCMX data for 2014. This can be seen most clearly in the longitude direction (Fig. 16a), in which the width peaks from 0000 to 0800 LT, but is relatively constant for the rest of the day. In the latitudinal direction (Fig. 16b), the diurnal variation depends on the month. In some months the width peaks around 1200 LT, while in other months it peaks in the late afternoon/evening hours. April is an anomaly with very high values occurring for the latitudinal widths from 0600 to 1200 LT. Similar results are seen for widths in the Mlon and Mlat directions (Figs. 15c and 15d). The diurnal variations of the widths for all datasets are shown in Supplemental Figure S9. Longitudinal widths for all three datasets are larger in the morning hours, peaking from 0200 to 0600 LT. SAMI3/WACCMX and JPL/GIM have a broad minimum from 1200 to 2000 LT, while SAMI3/climo has a secondary peak at 1600 LT. SAMI3/WACCMX has larger values than JPL/GIM during the morning, but similar values over the rest of the day, while SAMI3/climo has larger values at all times. The Mlat mean widths for SAMI3/climo and JPL/GIM show strong peaks at 0800–1000 MLT. SAMI3/WACCMX also shows higher values in the late morning, but not as sharp as for the other datasets.
Fig. 16 Width values versus local time from SAMI3/WACCMX (with EOF1 removed) for the coordinate (a) longitude, (b) latitude, (c) Mlon, (d), Mlat. Colors indicate month as shown at the top. All widths are in great circle degrees. 
To summarize these results, we find that the SAMI3/WACCMX TCCs show a structure that consists of a baseline global correlation with values that vary consistently with the amount of variation in the incident solar flux used to drive SAMI3. Much of this correlation can be described by the first EOF mode, which maximizes in the sunlit regions and has an amplitude that is highly correlated with F10.7. Removing this EOF allows us to examine the correlation structure on shorter scales. These correlations vary considerably by month and over the day. This is likely due to variations in the neutral atmosphere that is driving the ionosphere.
In contrast to baseline correlations, where we have not been able to find references for comparison, there have been several studies that have examined correlation length scales for TEC. For example, Klobuchar & Johanson (1977) estimated the correlation distance of mean electron content using daytoday variability of station pairs of TEC observations using VHF signals from geostationary satellites. They calculated a zonal (meridional) correlation distance of ~2900 km (~1800 km), using a correlation cutoff of 0.7. Yue et al. (2007) used JPL/GIM data to estimate correlation distances as a function of Mlat and LT for 2000 and 2005, using a correlation cutoff of 0.75. They showed a peak in the meridional correlation distances at low latitudes around sunrise (see Fig. 3b), with values ranging from ~11 to 15°. These lengths are somewhat smaller than we obtained, which may reflect differences in using the Gaussian fit versus correlation cutoff methods. In the longitudinal direction, Yue et al. (2007) showed a mean value of 20° near the Equator (Fig. 4), which is also somewhat shorter than our results. Shim et al. (2008) used >1000 groundbased TEC receivers to estimate low latitude zonal (meridional) correlation lengths of 11° (4°), using a correlation cutoff of 0.7. They also showed larger correlation lengths in daytime versus nighttime. McNamara (2009) used ionosondes in Australia and Papua New Guinea to calculate correlations and inferred zonal (meridional) correlation lengths of ~1500 km (1000 km). They also separated magnetically quiet days (Ap < 25) from magnetically disturbed days. Liu et al. (2018) examined correlation scales using JPL/GIM data for four different years (2008, 2009, 2014, and 2015). Using a Gaussian fitting method, they found correlation lengths vary with season and year. Low latitude meridional lengths in 2014 ranged from 18° to 40°, while zonal lengths ranged from 72° (winter) to 135° (summer). The very high zonal lengths in the summer are considerably larger than those we obtained with JPL/GIM (Supplemental Fig. S8c). This is likely due to the influence of solar activity, which was shown to be highly correlated with JPL/GIM in summer 2014 (as will be shown in Sect. 4.4). Our approach of filtering out the largescale correlations results in more localized correlation scales.
Forsythe et al. (2020a) examined correlation lengths using ΔTEC from JPL/GIM and showed qualitatively similar diurnal variability at the Equator to what we show in Figure 16 and Supplemental Figure S9, with peaks in the morning hours. Their zonal widths peaked at ~20° (our widths peak at ~40°) and their meridional widths peaked at ~23° (our widths peak at ~30°). They also calculated correlation lengths using differences between IRI and JPL/GIM to represent model errors. Their Equatorial zonal correlation lengths vary from ~20 to 40°, while the meridional correlation lengths vary from ~10 to 20°. The Mlat lengths calculated by Forsythe et al. (2020a) at the Equator show a diurnal cycle peaking in the morning, similar to what is seen in this study.
4.3 Conjugate correlations
Up to this point, we have focused on correlations with a reference point at the Equator (geographic and magnetic). If we move the reference point off the Equator, we expect to encounter correlations with grid points that have the same magnetic field line in the opposite hemisphere. These correlations show up most clearly in the magnetic coordinates. Magnetic conjugate point correlation (CPC, for short) analysis with TEC was also performed by Yue et al. (2007), using JPL/GIM data for 2000 and 2005, using three months of daytoday differences. They showed strong CPC occurs in all seasons and years, and they also showed a diurnal cycle, peaking in the afternoon. Also, they showed a dip in CPC at sunrise and sunset at higher Mlats due to the time delay between different magnetic conjugate points (as inferred from Fig. 3). CPCs were also examined by Shim et al. (2008), using daytoday TEC variations for four 30day periods in 2004 based on >1000 groundbased GPS receivers, and Liu et al. (2018) using JPL/GIM in four different years (2008, 2009, 2014, and 2015).
To illustrate the CPC correlation in SAMI3/WACCMX, we plot in Figure 17a the TCC for reference Mlat = 20° N for 1400 MLT and date range 2–31 January 2014. While strong correlations occur throughout the daytime hours, there are two main correlation peaks. Besides the main peak at the reference Mlat, there is a secondary peak near Mlat = 20° S. A line plot through this TCC (Fig. 17b) shows the two peaks, with a red circle denoting the conjugate point. It is likely that some of this “conjugate” correlation is caused by the solar variations that simultaneously influence both Mlats, independent of the magnetic connection. As discussed above, we can reduce the solar influence by filtering out the first EOF. The resulting filtered TCC map (Fig. 17c) and line plot (Fig. 17d) show the CPC more sharply than the unfiltered results. We note that the maximum composite correlation in this plot is not exactly at the conjugate Mlat of 20° S. We hypothesize that this may be due to recombination at nighttime occurring faster at higher Mlats, thereby pushing the conjugate correlation peaks equatorward.
Fig. 17 (a) SAMI3/WACCMX TEC correlation composite at reference time 1400 MLT_{0} for 2–31 January 2014 as a function of Mlat and approximate local time. The reference point is marked by a black dot. Thick black lines are correlations of 0.5, 0.7, and 0.9. (b) Line plot of TCC in Mlat at 1400 MLT_{0}. Dotted lines indicate the reference and conjugate Mlat, and red dot indicates the conjugate correlations. (c), (d), Same as (a), (b) except that EOF1 was removed from the TEC data before calculating the TCCs. (e) Conjugate correlation versus MLT for SAMI3/WACCMX for each month of 2014. Lines colorcoded by month as shown. 
In order to determine the CPC value, we search for the maximum composite correlation within ±10° of the conjugate Mlat. As shown in Figure 17d, this point is shifted slightly from the conjugate Mlat. We calculated the CPC for all months and reference MLTs for Mlat = 20° S, and plot the results in Figure 17e. We see that the CPC generally has a diurnal cycle with higher daytime values, peaking at ~0.6 near 1600–1800 MLT and minimizing at ~0.4 at 0000 MLT. This diurnal variation is consistent with Mlat = 20° N CPC presented by Yue et al. (2007).
In Figure 18 we compare the CPCs calculated from the three datasets for 2014. These are all calculated with the EOF1 filtering and are shown as a function of the month (top row) and MLT (bottom row). The mean monthly CPCs (Figs. 18a–18c) range from ~0.3 to 0.7 for SAMI3/climo and ~0.5 to 0.6 for SAMI3/WACCMX. JPL/GIM has somewhat lower values and less monthtomonth variability, while SAMI3/climo shows very large variability, peaking in March and June. Diurnal variations (Figs. 18d–18f) are seen in all three datasets, although the timing of the peak varies. JPL/GIM shows a peak at 1200 to 1600 LT, while SAMI3/climo peaks from 0200 to 0800 LT. Yue et al. (2007) calculated CPCs at ±20° Mlat (their Fig. 2b) for 2000 and 2005 and found strong variation with the season. The highest correlations were observed in February – April 2005, with values ranging from ~0.4 – 0.7, peaking from 1200 to 1400 LT. Their local time variations are somewhat different than what we obtained with the 2014 JPL/GIM data. This may reflect differences in solar and magnetic variability between the two years, but also may be due to differences in approach. Whereas we use daily values of TEC over monthly intervals, Yue et al. (2007) used ΔTEC over threemonth intervals. CPCs for 2010 are shown in Supplemental Figure S10. Both JPL/GIM and SAMI3/climo show similar ranges of mean correlations and variability in 2010 as in 2014. Further analysis of the CPCs is warranted to determine the cause of the differences in year and dataset, but overall we find that robust CPC correlations are observed both in SAMI3 and in the JPL/GIM for most months and local times.
Fig. 18 (a) – (c) Monthly conjugate correlations at reference latitude 20° N (after removing the EOF1) for SAMI3/climo, SAMI3/WACCMX, and JPL/GIM for 2014. Red lines indicate different local times, and blue lines are the averages over all local times. (d) – (f) Conjugate correlations as a function of magnetic local time (after removing the EOF1) for SAMI3/climo, SAMI3/WACCMX, and JPL/GIM. Red lines indicate different months, and blue lines are the averages over all months. 
4.4 Correlations of EOF modes with F10.7 and Ap indices
We conclude by comparing the correlations of EOF modes with the F10.7 and Ap indices. Figure 2c showed correlations of F10.7 with the EOF1 index at 12 UTC for SAMI3/WACCMX for 2014. With the exception of March, the correlation with F10.7 is high (>0.6) during 2014. Figure 2d shows that the correlations with Ap with EOF1 are generally weak (<0.6) except for February, which is a month with a large Ap standard deviation. Here we will further quantify how these correlations depend on the dataset, month, universal time, and EOF mode.
Figure 19 plots the correlation of EOF1 and EOF2 with F10.7 for all months and times (UTC) for the three datasets. Note that at each time when the correlation is significant (at the 95% confidence level, using the pvalue from a twosided Ttest), a circle is included on the line plot. Line segments that do not show circles indicate statistically insignificant correlations. The SAMI3/WACCMX (red lines) shows significant correlations of F10.7/EOF1 at all local times for all months except February, where three universal times (0000, 0200, and 0400 UT) are not significant. Therefore, of the 144 possible month/time combinations, 141 (98%) have the EOF1 significantly correlated with F10.7. The JPL/GIM EOF1 generally correlates very well with F10.7 in April–December but does not correlate as well in the first part of the year. In total, 114/144 or 79% of the month/time points are significant. Supplemental Figure S11 shows the correlations of EOF1 and EOF2 with F10.7 for 2010 for SAMI3/climo and JPL/GIM. Significant correlations are observed for a wide range of months and times, although the peak correlations are smaller than the peak correlations in 2014. April and June show very small correlations for both datasets, while March, May, and July show very large correlations.
Fig. 19 Correlations of the EOF1 and EOF2 amplitudes with F10.7 as a function of time (UTC) and month of 2014. If the correlation is significant (p < 0.050 is 95% confidence level) then a circle is included on the plot. The colors indicate the experiments SAMI3/climo (black), SAMI3/WACCMX (red), and JPL/GIM (blue). Thick lines and filled circles are for EOF1, while thin lines and open circles are for EOF2. 
We also performed a similar analysis with the Ap index. Figure 20 shows that while EOF1 sometimes exhibits strong correlations with Ap in 2014 (e.g., SAMI3 data for February, JPL/GIM data for August), in general, the correlation is rather weak. For EOF2 the correlation is considerably stronger, showing significant correlations for all data sets in numerous months and times. The SAMI3/climo in particular shows very strong correlations of EOF2 with Ap, which is consistent with the HWM neutral winds being strongly driven by Ap, whereas SAMI3/WACCMX is driven by more realistic neutral winds from a global circulation model. The Ap correlations with EOF1 and EOF2 for 2010 are shown in Supplemental Figure S12 for SAMI3/climo and JPL/GIM. For the SAMI3/climo the Ap/EOF1 correlations are generally small and insignificant (except for January and April), while Ap/EOF2 correlations can be quite large (e.g., May through October and December). It appears that even in this year of relatively weak variability the dominant mode is driven by F10.7, while Ap provides a strong secondary forcing. For JPL/GIM in 2010, there are some significant Ap/EOF1 correlations, and these are generally weaker than in SAMI3/climo. The Ap/EOF2 correlations are generally weak for JPL/GIM. Further analysis will be necessary to elucidate the exact physical processes involved in these EOF correlations, but these diagnostics motivate future research in understanding TEC correlations. We note that EOF analyses can be done with multiple variables to examine the variability of different components of the ionospheric system.
Fig. 20 Correlations of the EOF1 and EOF2 amplitudes with Ap as a function of time (UTC) and month of 2014. If the correlation is significant (p < 0.050 is 95% confidence level) then a circle is included on the plot. The colors indicate the experiments SAMI3/climo (black), SAMI3/WACCMX (red), and JPL/GIM (blue). Thick lines and filled circles are for EOF1, while thin lines and open circles are for EOF2. 
5 Summary and conclusions
This study examined composites of lowlatitude TEC spatial correlations as a function of (1) three datasets, (2) model grid, (3) local time, (4) month, (5) year, and (6) Equator versus subtropics. The results showed the SAMI3/climo, driven by empirical neutral fields, had unrealistically large correlation structures, while SAMI3/WACCMX, driven by more realistic winds based on a dynamical model, showed correlations that agreed better with the analyses from JPL/GIM. Base correlation values showed a strong correlation with solar variability in all datasets, as well as diurnal variability, with peak values in the late morning to evening hours. Removal of the base correlations was attempted by filtering the leading EOF from the TEC dataset. This resulted in more distinct Gaussian shapes that could be fitted to zonal and meridional widths. The zonal widths showed a semidiurnal variation, while the meridional width varied diurnally. Conjugate correlations were also examined. Both the SAMI3 runs and JPL/GIM analyses showed strong correlations between magnetic latitudes of 20° N and 20° S. The SAMI3/WACCMX and JPL/GIM showed a weak diurnal cycle in the conjugate correlations, peaking in the late afternoon, which resembled the meridional base correlations.
The results from this study show that all three datasets examined have Equatorial correlation structures that can be characterized by a globallyconstant base value and a local width. These two features can often be separated by the removal of EOFs describing correlations associated with solar variability. In application to DA, this suggests that it may be possible to break the analysis problem into a global part that solves the amplitude of one or more leading EOFs and a local part that could be solved using traditional localization methods. Off the Equator, there is a secondary issue of the conjugate correlations. As far as we know, none of the current ionospheric DA methods deals specifically with this feature, except perhaps indirectly in cases where localization scales are large enough that ensembles capture this correlation. However, DA systems are generally formulated in the geographic reference frame. Moving to a geomagnetic reference frame may provide a way forward since correlations could be carried along magnetic field lines. For SAMI3, which separates alongfield and crossfield motion, this approach could be very appealing. However, characterizing the BECM based on full 3D electron density fields rather than 2D TEC fields will be much more challenging.
Future work will involve analyzing the TEC correlation structures in the middle and high latitudes. In addition, we would like to analyze the EOF structures in more detail to understand how these analyses and forecasts respond to different forcings, including additional magnetic indices. Finally, we plan to compare these analyses with the correlation structures of actual SAMI3 ensembles used in a DA system. The knowledge gained from the work done in this paper will help to inform how best to generate ensemble perturbations for ionospheric DA as well as localization methods for reducing spurious correlations.
Supplementary materials
Supplemental Figure S1. Daily TEC values at latitude = 0°, longitude = 0° for 2010 from SAMI3/climo (black) and JPL/GIM (blue). Panels are for different times of day: (a) 0000 UTC, (b) 0600 UTC, (c) 1200 UTC, and (d) 1800 UTC. 
Supplemental Figure S2. (a) Daily values of F10.7 over the year 2010. (b) Daily values of Ap index over 2010. (c) Monthly standard deviation of F10.7 (black line) for 2010 along with the correlation of the 1200 UTC EOF1 amplitude with F10.7 (red line). (d) Monthly standard deviation of Ap index (black line) for 2010 along with the correlation of the SAMI3/WACCMX 1200 UTC EOF1 with Ap (red line). 
Supplemental Figure S3. SAMI3/climo and JPL/GIM TEC correlation composites for 2–31 January 2010 as a function of latitude and LT. The reference point is marked by the black dot. It is located on the Equator, and the reference LT varies with panel (0000, 0600, 1200, and 1800 are shown). Thick black lines indicate correlation of 0.5, 0.7, and 0.9. Thin black lines are SZA on 16 January 2010 at 10°, 20°, …, 90°. 
Supplemental Figure S4. SAMI3/climo and JPL/GIM line plots of TCC versus relative LT (black) for 2–31 January 2010 at four reference local times (0000, 0600, 1200, and 1800 LT). The reference latitude is the Equator. The Gaussian fits are shown by the dashed blue lines, with the two fitted parameters shown on the bottom of each plot. Base is unitless and width is in hours LT. 
Supplemental Figure S5. (a) – (c) Monthly longitudinal base values (red lines) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM along with the monthly for 2014. Red lines are for different reference LTs, and blue line is average over a reference LTs. (d) – (f) Monthly Mlat base values (red lines) from the SAMI3/climo, SAMI3/WACCMX, and JPL/GIM along with monthly . Red lines are for different reference LTs, and blue line is average over a reference LTs. 
Supplemental Figure S6. (a) – (c) Longitudinal base values as a function of local time (red lines) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM along with the monthly for 2014. Red lines are for different reference months of 2014, and blue line is average over a reference LTs. (d) – (f) Monthly Mlat base values as a function of local time (red lines) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM along with monthly . Red lines are for different months, and blue line is average over a reference LTs. 
Supplemental Figure S7. SAMI3/climo, SAMI3/WACCMX, and JPL/GIM line plots of TCC versus relative LT (black) for 2–31 January 2014 at four reference local times (0000, 0600, 1200, and 1800 LT). The leading EOF has been removed from the TEC time series before calculating the TCCs. The reference latitude is the Equator. The Gaussian fits are shown by the dashed blue lines, with the two fitted parameters shown on the bottom of each plot. Base is unitless and width is in hours LT. 
Supplemental Figure S8. (a) – (c) Monthly longitudinal width values (after removing the EOF1) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM for 2014. Red lines indicate different local times, and blue lines are the averages over all local times. (d) – (e) Monthly Mlat width values (after removing the EOF1) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM. Red lines indicate different local times, and blue lines are the averages over all local times. All widths are in great circle degrees. 
Supplemental Figure S9. (a) – (c) Longitudinal width values as a function of local time (after removing the EOF1) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM for 2014. Red lines indicate different months, and blue lines are the averages over all months. Mlat width as a function of magnetic local time (after removing the EOF1) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM. Red lines indicate different months, and blue lines are the averages over all months. All widths are in great circle degrees. 
Supplemental Figure S10. (a) – (c) Monthly conjugate correlations at reference latitude 20°N (after removing the EOF1) for SAMI3/climo, SAMI3/WACCMX, and JPL/GIM for 2010. Red lines indicate different local times, and blue lines are the averages over all local times. (d) – (f) Conjugate correlations as a function of magnetic local time (after removing the EOF1) for SAMI3/climo, SAMI3/WACCMX, and JPL/GIM. Red lines indicate different months, and blue lines are the averages over all months. 
Supplemental Figure S11. Correlations of the EOF1 and EOF2 amplitudes with F10.7 as a function of time (UTC) and month of 2010. If the correlation is significant (p < 0.050 is 95% confidence level) then a circle is included on the plot. The colors indicate the experiments SAMI3/climo (black) and JPL/GIM (blue). Thick lines and filled circles are for EOF1, while thin lines and open circles are for EOF2. 
Supplemental Figure S12. Correlations of the EOF1 and EOF2 amplitudes with Ap as a function of time (UTC) and month of 2010. If the correlation is significant (p < 0.050 is 95% confidence level) then a circle is included on the plot. The colors indicate the experiments SAMI3/climo (black) and JPL/GIM (blue). Thick lines and filled circles are for EOF1, while thin lines and open circles are for EOF2. 
Acknowledgments
This work was funded by the U.S. Office of Naval Research (ONR). The authors and the editor thank two anonymous reviewers for their assistance in evaluating this paper.
References
 AraujoPradere EA, FullerRowell TJ, Bilitza D. 2004. Ionospheric variability for quiet and perturbed conditions. Adv Space Res 34(9): 1914–1921. https://doi.org/10.1016/j.asr.2004.06.007. [CrossRef] [Google Scholar]
 AraujoPradere EA, FullerRowell TJ, Codrescu MV, Bilitza D. 2005. Characteristics of the ionospheric variability as a function of season, latitude, localtime, and geomagnetic activity. Radio Sci 40(5): RS5009. https://doi.org/10.1029/2004RS003179. [Google Scholar]
 Bergot N, Tsagouri I, Bruyninx C, Legrand J, Chevalier JM, Defraigne P, Baire Q, Pottiaux E. 2013. The influence of space weather on ionospheric total electron content during the 23rd solar cycle. J Space Weather Space Clim 3: A25. https://doi.org/10.1051/swsc/2013047. [CrossRef] [EDP Sciences] [Google Scholar]
 Bust GS, Coker C, Coco DS, Gaussiran TL II, Lauderdale T. 2001. IRI data ingestion and ionospheric tomography. Adv Space Res 27(1): 157–165. https://doi.org/10.1016/S02731177(00)001630. [CrossRef] [Google Scholar]
 Bust GS, Garner TW, Gaussiran TL II. 2004. Ionospheric data assimilation threedimensional (IDA3D): A global, multisensor, electron density specification algorithm. J Geophys Res: Space Phys 109(A11): A113112. https://doi.org/10.1029/2003JA010234. [Google Scholar]
 Chartier AT, Matsuo T, Anderson JL, Collins N, Hoar TJ, Lu G, Mitchell CN, Coster AJ, Paxton LJ, Bust GS. 2016. Ionospheric data assimilation and forecasting during storms. J Geophys Res: Space Phys 121(1): 764–778. https://doi.org/10.1002/2014JA020799. [CrossRef] [Google Scholar]
 Chen Z, Zhang SR, Coster AJ, Fang G. 2015. EOF analysis and modeling of TEC climatology over North America. J Geophys Res: Space Phys 120(4): 3118–3129. https://doi.org/10.1002/2014JA020837. [CrossRef] [Google Scholar]
 Chen CH, Lin CH, Matsuo T, Chen WH, Lee IT, Liu JY, Lin JT, Hsu CT. 2016. Ionospheric data assimilation with thermosphereionosphereelectrodynamics general circulation model and GPSTEC during geomagnetic storm conditions. J Geophys Res: Space Phys 121(6): 5708–5722. https://doi.org/10.1002/2015JA021787. [CrossRef] [Google Scholar]
 Chen CH, Lin CHC, Matsuo T. 2019. Ionospheric responses to the 21 August 2017 solar eclipse by using data assimilation approach. Prog Earth Planet Sci 6(13): 1–9. https://doi.org/10.1186/s4064501902634. [CrossRef] [Google Scholar]
 Codrescu SM, Codrescu MV, Fedrizzi M. 2018. An ensemble Kalman filter for the thermosphereionosphere. Space Weather 16(1): 57–68. https://doi.org/10.1002/2017SW001752. [CrossRef] [Google Scholar]
 Drob DP, Emmert JT, Meriwether JW, Makela JJ, Doornbos E, et al. 2015. An update to the Horizontal Wind Model (HWM): The quite time thermosphere. Earth Space Sci 2(7): 301–319. https://doi.org/10.1002/2014ea000089. [CrossRef] [Google Scholar]
 Eckermann SD, Ma J, Hoppel KW, Kuhl DD, Allen DR, et al. 2018. High altitude (0–100 km) global atmospheric reanalysis system: Description and application to the 2014 austral winter of the Deep Propagating GravityWave Experiment (DEEPWAVE). Mon Wea Rev 146: 2639–2666. https://doi.org/10.1175/MWRD170386.1. [CrossRef] [Google Scholar]
 Elvidge S, Angling MJ. 2019. Using the local ensemble transform Kalman filter for upper atmospheric modelling. J Space Weather Space Clim 9: A30. https://doi.org/10.1051/swsc/2019018. [CrossRef] [EDP Sciences] [Google Scholar]
 Emmert JT, Mannucci AJ, McDonald SE, Vergados P. 2017. Attribution of interminimum changes in global and hemispheric total electron content. J Geophys Res: Space Phys 122(2): 2424–2439. https://doi.org/10.1002/2016JA023680. [CrossRef] [Google Scholar]
 Fang TW, Akmaev R, FullerRowell T, Wu F, Maruyama N, Millward G. 2013. Longitudinal and daytoday variability in the ionosphere from lower atmosphere tidal forcing. Geophys Res Lett 40(11): 2523–2528. https://doi.org/10.1002/grl.50550. [CrossRef] [Google Scholar]
 Fang TW, FullerRowell T, Yudin V, Matsuo T, Viereck R. 2018. Quantifying the sources of ionosphere daytoday variability. J Geophys Res: Space Phys 123(11): 9682–9696. https://doi.org/10.1029/2018JA025525. [CrossRef] [Google Scholar]
 Forbes JM, Palo SE, Zhang X. 2000. Variability of the ionosphere. J Atmos SolTerr Phys 62(8): 685–693. https://doi.org/10.1016/S13646826(00)000298. [CrossRef] [Google Scholar]
 Forsythe VV, Azeem I, Crowley G. 2020a. Ionospheric horizontal correlation distances: Estimation, analysis, and implications for ionospheric data assimilation. Radio Sci 55(12): e2020RS007159. https://doi.org/10.1029/2020rs007159. [Google Scholar]
 Forsythe VV, Azeem I, Crowley G, Makarevich RA, Wang C. 2020b. The global analysis of the ionospheric correlation time and its implications for ionospheric data assimilation. Radio Sci 55(12): e2020RS007181. https://doi.org/10.1029/2020rs007181. [Google Scholar]
 Forsythe VV, Azeem I, Crowley G, Themens DR. 2021a. Ionospheric vertical correlation distances: Estimation from ISR data, analysis, and implications for ionospheric data assimilation. Radio Sci 56(2): e2020RS007177. https://doi.org/10.1029/2020RS007177. [Google Scholar]
 Forsythe VV, Azeem I, Blay R, Crowley G, Gasperini F, Hughes J, Makarevich RA, Wu W. 2021b. Evaluation of the new background error covariance model for the ionospheric data assimilation. Radio Sci 56(8): e2021RS007286. https://doi.org/10.1029/2021rs007286. [Google Scholar]
 Gail WB, Prag AB, Coco DS, Coker C. 1993. A statistical characterization of local midlatitude total electron content. J Geophys Res: Space Phys 98(A9): 15717–15727. https://doi.org/10.1029/92ja01597. [CrossRef] [Google Scholar]
 Gaspari G, Cohn SE. 1999. Construction of correlation functions in two and three dimensions. Quart J Roy Meteor Soc 125(554): 723–757. https://doi.org/10.1002/qj.49712555417. [CrossRef] [Google Scholar]
 Greybush SJ, Kalnay E, Miyoshi T, Ide K, Hunt BR. 2011. Balance and ensemble Kalman filter localization techniques. Mon Wea Rev 139: 511–522. https://doi.org/10.1175/2010MWR3328.1. [CrossRef] [Google Scholar]
 Hsu CT, Matsuo T, Wang W, Liu JY. 2014. Effects of inferring unobserved thermospheric and ionospheric state variables by using an Ensemble Kalman filter on global ionospheric specification and forecasting. J Geophys Res: Space Phy 119(11): 9256–9267. https://doi.org/10.1002/2014JA020390. [CrossRef] [Google Scholar]
 Huba JD, Joyce G, Fedder JA. 2000. Sami2 is Another Model of the Ionosphere (SAMI2): A new lowlatitude ionosphere model. J Geophys Res: Space Phys 105(A10): 23035–23053. https://doi.org/10.1029/2000JA000035. [CrossRef] [Google Scholar]
 Huba JD, Maunte A, Crowley G. 2017. SAMI3_ICON: Model of the ionosphere/plasmasphere system. Space Sci Rev 212: 731–742. https://doi.org/10.1007/s112140170415z. [CrossRef] [Google Scholar]
 Iijima BI, Harris IL, Ho CM, Lindqwister UJ, Mannucci AJ, Pi X, Reyes MJ, Sparks LC, Wilson BD. 1999. Automated daily process for global ionospheric total electron content maps and satellite ocean altimeter calibration based on Global Positioning System data. J Atmos SolTerr Phys 61(16): 1205–1218. https://doi.org/10.1016/S13646826(99)00067X. [CrossRef] [Google Scholar]
 Kepert JD. 2009. Covariance localization and balance in an ensemble Kalman filter. Quart J Roy Meteor Soc 135(642): 1157–1176. https://doi.org/10.1002/qj.443. [CrossRef] [Google Scholar]
 Khattatov B, Murphy M, Gnedin M, Sheffel J, Adams J, Cruickshank B, Yudin V, FullerRowell T, Retterer J. 2005. Ionospheric nowcasting via assimilation of gps measurements of ionospheric electron content in a global physicsbased timedependent model. Quart J Roy Meteor Soc 131(613): 3543–3559. https://doi.org/10.1256/qj.05.96. [CrossRef] [Google Scholar]
 Klobuchar JA, Johanson JM. 1977. Correlation distance of meandaytime electron content. In: Rep. AFGLTR771885, Air Force Geophys. Lab. Bedford, Mass. https://apps.dtic.mil/sti/citations/ADA048117. [Google Scholar]
 Laundal KM, Richmond AD. 2016. Magnetic coordinate systems. Space Sci Rev 206: 27–59. https://doi.org/10.1007/s112140160275y. [Google Scholar]
 Lin CY, Matsuo T, Liu JY, Lin CH, Tsai HF, AraujoPradere EA. 2015. Ionospheric assimlation of radio occultation and groundbased GPS data using nonstationary background model error covariance. Atmos Meas Tech 8: 171–182. https://doi.org/10.5194/amt81712015. [CrossRef] [Google Scholar]
 Liu L, Chen Y. 2009. Statistical analysis of solar activity variations of total electron content derived at Jet Propulsion Laboratory from GPS observations. J Geophys Res: Space Phys 114(A10): A10311. https://doi.org/10.1029/2009JA014533. [CrossRef] [Google Scholar]
 Liu HL, Foster BT, Hagan ME, McInerney JM, Maute A, et al. 2010. Thermosphere extension of the Whole Atmosphere Community Climate Model. J Geophys Res 115(A12302). https://doi.org/10.1029/2010JA015586. [Google Scholar]
 Liu S, Yang J, Yu T, Zhang Z. 2018. Horizontal spatial correlations of the ionospheric TEC derived from GPS global ionospheric maps. Adv Space Res 62(7): 1775–1786. https://doi.org/10.1016/j.asr.2018.06.042. [CrossRef] [Google Scholar]
 Mannucci AJ, Wilson BD, Yuan DN, Ho CH, Lindqwister UJ, Runge TF. 1998. A global mapping technique for GPSderived ionospheric total electron content measurements. Radio Sci 33(3): 565–582. https://doi.org/10.1029/97rs02707. [CrossRef] [Google Scholar]
 Matsuo T, Lee IT, Anderson JL. 2013. Thermospheric mass density specification using and ensemble Kalman filter. J Geophys Res: Space Phys 118(3): 1339–1350. https://doi.org/10.1002/jgra.50162. [CrossRef] [Google Scholar]
 McDonald SE, Sassi F, Tate J, McCormack J, Kuhl DD, Drob DP, Metzler C, Mannucci AJ. 2018. Impact of nonmigrating tides on the low latitude ionosphere during a sudden stratospheric warming event in January 2010. J Atmos SolTerr Phys 171: 188–200. https://doi.org/10.1016/j.jastp.2017.09.012. [CrossRef] [Google Scholar]
 McNamara LF. 2009. Spatial correlations of foF2 deviations and their implications for global ionospheric models: 2. Digisondes in the United States, Europe, and South Africa. Radio Sci 44(2): RS2017. https://doi.org/10.1029/2008RS003956. [Google Scholar]
 McNamara LF, Wilkinson PJ. 2009. Spatial correlations of foF2 deviations and their implications for global ionospheric models: 1. Ionosondes in Australia and Papua New Guinea. Radio Sci 44(2): RS2016. https://doi.org/10.1029/2008RS003955. [Google Scholar]
 Mengist CK, Ssessanga N, Jeong SH, Kim JH, Kim YH, Kwak YS. 2019. Assimilation of multiple data types to a regional ionospheric model with 3DVar algorithm (IDA4D). Space Weather 17(7): 1018–1039. https://doi.org/10.1029/2019SW002159. [Google Scholar]
 Monahan AM, Fyfe JC, Ambaum MHP, Stephenson DB, North GR. 2009. Empirical orthogonal functions: The medium is in the message. J Clim 22(24): 6501–6514. https://doi.org/10.1175/2009JCLI3062.1. [CrossRef] [Google Scholar]
 Morozov AV, Ridley AJ, Bernstein DS, Collins N, Hoar TJ, Anderson JL. 2013. Data ssimilation and driver estimation for the Global IonosphereThermosphere Model using the ensemble adjustment Kalman filter. J Atmos SolTerr Phys 104: 126–136. https://doi.org/10.1016/j.jastp.2013.08.016. [CrossRef] [Google Scholar]
 Morzfeld M, Hodyss D. 2023. A theory for why even simple covariance localization is so useful in ensemble data assimilation. Mon Wea Rev 151: 717–736. https://doi.org/10.1175/MWRD220255.1. [CrossRef] [Google Scholar]
 Nisbet JS, Tyrnov OF, Zintchenko GN, Ross WJ. 1981. Limits on the accuracy of correction of transionospheric propagation errors by using ionospheric models based on solar and magnetic indices and local measurements. Radio Sci 16(1): 127–133. https://doi.org/10.1029/Rs016i001p00127. [CrossRef] [Google Scholar]
 Pedatella NM, Liu HL, Marsh DR, Raeder K, Anderson JL, Chau JL, Goncharenko LP, Siddiqui TA. 2018. Analysis and hindcast experiments of the 2009 sudden stratospheric warming in WACCMX+DART. J Geophys Res: Space Phys 123(4): 3131–3153. https://doi.org/10.1002/2017JA025107. [CrossRef] [Google Scholar]
 Pedatella NM, Anderson JL, Chen CH, Raeder K, Liu J, Liu HL, Lin CH. 2020. Assimilation of ionospheric observations in the Whole Atmosphere Community Climate Model with ThermosphereIonosphere Extension (WACCMX). J Geophsy Res: Space Phys 125(9): e2020JA028251. https://doi.org/10.1029/2020JA028251. [Google Scholar]
 Picone JM, Hedin AE, Drob DP, Aikin AC. 2002. NRLMSISE00 empirical model of the atmosphere: Statistical comparisons and scientific issues. J Geophys Res: Space Phys 107(A12): SIA 151–SIA 1516. https://doi.org/10.1029/2002JA009430. [CrossRef] [Google Scholar]
 Richmond AD. 1995. Ionospheric electrodynamics using magnetic apex coordinates. J Geomagn Geoelectr 47(2): 191–212. https://doi.org/10.5636/Jgg.47.191. [CrossRef] [Google Scholar]
 Rishbeth H, Mendillo M. 2001. Patterns of F2layer variability. J Atmos SolTerr Phys 63(15): 1661–1680. https://doi.org/10.1016/S13646826(01)000360. [CrossRef] [Google Scholar]
 Rush CM, Gibbs J. 1973. Predicting the daytoday variability of the midlatitude ionosphere for application to HF propagation predictions. Tech. Rep. AD764 711. https://apps.dtic.mil/sti/citations/AD0764711. [Google Scholar]
 Scherliess L, Schunk RW, Sojka JJ, Thompson DC. 2004. Development of a physicsbased reduced state Kalman filter for the ionosphere. Radio Sci 39(1): RS1S04. https://doi.org/10.1029/2002RS002797. [Google Scholar]
 Shim JS, Scherliess L, Shunk RW, Thompson DC. 2008. Spatial correlations of daytoday ionospheric total electron content variability obtained from groundbased GPS. J Geophys Res: Space Phys 113(A9): A09309. https://doi.org/10.1029/2007JA012635. [Google Scholar]
 Siddiqui TA, Yamazaki Y, Stolle C, Maute A, Laštovička J, Edemskiy IK, Mošna Z, Sivakandan M. 2021. Understanding the total electron contant variability over Europe during 2009 and 2019 SSWs. J Geophys Res: Space Phys 126(9): e2020JA028751. https://doi.org/10.1029/2020JA028751. [CrossRef] [Google Scholar]
 Solomontsev DV, Khattatov BV, Codrescu MV, Titov AA, Yudin V, Khattatov VU. 2012. Ionosphere state and parameter estimation using the Ensemble Square Root Filter and the global threedimensional firstprincipal model. Space Weather 10(7): S07004. https://doi.org/10.1029/2012SW000777. [Google Scholar]
 Talaat ER, Shu X. 2016. Spatial and temporal variation of total electron content as revealed by principal component analysis. Ann Geophys 34: 1109–1118. https://doi.org/10.5194/angeo3411092016. [CrossRef] [Google Scholar]
 Vaishnav R, Jacobi C, Berdermann J. 2019. Longterm trends in the ionospheric response to solar extremeultraviolet variations. Ann Geophys 37: 1141–1159. https://doi.org/10.5194/angeo3711412019. [CrossRef] [Google Scholar]
 Yue X, Wan W, Liu L, Mao T. 2007. Statistical analysis on spatial correlation of ionospheric daytoday variability by using GPS and Incoherent Scatter Radar observations. Ann Geophys 25: 1815–1825. https://doi.org/10.5194/angeo2518152007. [CrossRef] [Google Scholar]
 Zawdie KA, Dhadly MS, McDonald SE, Sassi F, Coker C, Drob DP. 2020. Daytoday variability of the bottomside ionosphere. J Atmos SolTerr Phys 205: 105299. https://doi.org/10.1016/j.jastp.2020.105299. [CrossRef] [Google Scholar]
 Zhang ML, Liu L, Li Q. 2021. Modeling the global ionospheric electron density based onthe EOF decomposition of the ionospheric radio occultation observation. Adv Space Res 68(5): 2218–2232. https://doi.org/10.1016/j.asr.2020.09.033. [CrossRef] [Google Scholar]
 Zhong J, Lei J, Yue X, Wang W, Burns AG, Luan X, Dou X. 2019. Empirical orthogonal function analysis and modeling of the topside ionospheric and plasmaspheric TECs. J Geophys Res: Space Phys 124(5): 3681–3698. https://doi.org/10.1029/2019JA026691. [CrossRef] [Google Scholar]
Cite this article as: Allen DR, Hodyss D, Forsythe VV & McDonald SE 2023. Low latitude monthly total electron content composite correlations. J. Space Weather Space Clim. 13, 7. https://doi.org/10.1051/swsc/2023005.
All Figures
Fig. 1 Daily TEC values at latitude = 0°, longitude = 0° for 2014 from SAMI3/climo (black), SAMI3/WACCMX (red), and JPL/GIM (blue). Panels are for different times of day: (a) 0000 UTC, (b) 0600 UTC, (c) 1200 UTC, and (d) 1800 UTC. 

In the text 
Fig. 2 (a) Daily values of F10.7 over the year 2014. (b) Daily values of Ap index over 2014. (c) Monthly standard deviation of F10.7 (black line) for 2014 along with the correlation of the 1200 UTC EOF1 amplitude with F10.7 (red line). (d) Monthly standard deviation of Ap index (black line) for 2014 along with the correlation of the SAMI3/WACCMX 1200 UTC EOF1 with Ap (red line). 

In the text 
Fig. 3 (a) The geographic grid used for this study, which is at the JPL/GIM longitude × latitude resolution of 5° × 2.5°. The Greenwich Meridian is highlighted by a red line. The blue lines indicate the 12 reference longitudes used in the analysis. (b) The geomagnetic grid used for this study is at Mlon × Mlat resolution of 5° × 2.5°. resolution. The location of magnetic longitude of 0° is indicated by the red line. The blue lines indicate the 12 reference magnetic longitudes used in the analysis. 

In the text 
Fig. 4 TEC correlation maps for 2–31 January 2014 at 1200 UTC for SAMI3/climo, SAMI3/WACCMX, and JPL/GIM. The reference point is marked by a black dot. It is located on the Equator, and the longitude varies with the panel (−180°, −90°, 0°, and 90° are shown). The black lines are SZA for 16 January 2014 at 10°, 20°, …, 90°. 

In the text 
Fig. 5 TEC correlation maps for 2–31 January 2014 at 1200 UTC on the magnetic grid for SAMI3/climo, SAMI3/WACCMX, and JPL/GIM. The reference point is marked by a black dot. It is located on the magnetic Equator, and the Mlon varies with panel (0°, 90°, 180°, and 270° are shown). The black lines are SZA on 16 January 2014 at 10°, 20°, …, 90°. 

In the text 
Fig. 6 SAMI3/climo, SAMI3/WACCMX, and JPL/GIM TEC correlation composites for 2–31 January 2014 as a function of latitude and LT. The reference point is marked by a black dot. It is located on the Equator, and the reference LT varies with panel (0000, 0600, 1200, and 1800 are shown). Thick black lines indicate a correlation of 0.5, 0.7, and 0.9. Thin black lines are SZA on 16 January 2014 at 10°, 20°, …, 90°. 

In the text 
Fig. 7 SAMI3/climo, SAMI3/WACCMX, and JPL/GIM line plots of TCC versus relative LT (black) for 2–31 January 2014 at four reference local times (0000, 0600, 1200, and 1800 LT). The reference latitude is the Equator. The Gaussian fits are shown by the dashed blue lines, with the two fitted parameters shown at the bottom of each plot. Base is unitless and width is in hours LT. 

In the text 
Fig. 8 Base values versus month from SAMI3/WACCMX for 2014 for the coordinate (a) longitude, (b) latitude, (c) Mlon, (d), Mlat. Colors indicate local time as shown at the top. Also plotted in each panel is the monthly . 

In the text 
Fig. 9 Base values versus local time from SAMI3/WACCMX for 2014 for the coordinate (a) longitude, (b) latitude, (c) Mlon, (d), Mlat. Colors indicate month as shown at the top. 

In the text 
Fig. 10 SAMI3/WACCMX EOFs for 2–31 January 2014 (1200 UTC). (a) – (c) The spatial pattern of the EOFs 1, 2, and 3 with the color label at the bottom. Thin black lines are SZA on 16 January 2014 at 10°, 20°, …, 90°. (d) – (f) The amplitudes of each EOF (black lines) along with the indices F10.7 (red lines) and Ap (blue lines), linearly scaled to provide the best fit to each of the amplitudes. The correlation coefficients of amplitudes with respect to the indices provide for each EOF at the bottom, along with the significance pvalue, where p < 0.050 is 95% confidence level. 

In the text 
Fig. 11 Same as Figure 10, but for SAMI3/climo. 

In the text 
Fig. 12 Same as Figure 10, but for JPL/GIM. 

In the text 
Fig. 13 Same as Figure 10, but for JPL/GIM at 1800 UTC. 

In the text 
Fig. 14 As in Figure 6, except the leading EOF has been removed from the TEC time series before calculating the TCCs. 

In the text 
Fig. 15 Width values versus month of 2014 from SAMI3/WACCMX (with EOF1 removed) for the coordinate (a) longitude, (b) latitude, (c) Mlon, (d), Mlat. Colors indicate local time as shown at the top. All widths are in great circle degrees. 

In the text 
Fig. 16 Width values versus local time from SAMI3/WACCMX (with EOF1 removed) for the coordinate (a) longitude, (b) latitude, (c) Mlon, (d), Mlat. Colors indicate month as shown at the top. All widths are in great circle degrees. 

In the text 
Fig. 17 (a) SAMI3/WACCMX TEC correlation composite at reference time 1400 MLT_{0} for 2–31 January 2014 as a function of Mlat and approximate local time. The reference point is marked by a black dot. Thick black lines are correlations of 0.5, 0.7, and 0.9. (b) Line plot of TCC in Mlat at 1400 MLT_{0}. Dotted lines indicate the reference and conjugate Mlat, and red dot indicates the conjugate correlations. (c), (d), Same as (a), (b) except that EOF1 was removed from the TEC data before calculating the TCCs. (e) Conjugate correlation versus MLT for SAMI3/WACCMX for each month of 2014. Lines colorcoded by month as shown. 

In the text 
Fig. 18 (a) – (c) Monthly conjugate correlations at reference latitude 20° N (after removing the EOF1) for SAMI3/climo, SAMI3/WACCMX, and JPL/GIM for 2014. Red lines indicate different local times, and blue lines are the averages over all local times. (d) – (f) Conjugate correlations as a function of magnetic local time (after removing the EOF1) for SAMI3/climo, SAMI3/WACCMX, and JPL/GIM. Red lines indicate different months, and blue lines are the averages over all months. 

In the text 
Fig. 19 Correlations of the EOF1 and EOF2 amplitudes with F10.7 as a function of time (UTC) and month of 2014. If the correlation is significant (p < 0.050 is 95% confidence level) then a circle is included on the plot. The colors indicate the experiments SAMI3/climo (black), SAMI3/WACCMX (red), and JPL/GIM (blue). Thick lines and filled circles are for EOF1, while thin lines and open circles are for EOF2. 

In the text 
Fig. 20 Correlations of the EOF1 and EOF2 amplitudes with Ap as a function of time (UTC) and month of 2014. If the correlation is significant (p < 0.050 is 95% confidence level) then a circle is included on the plot. The colors indicate the experiments SAMI3/climo (black), SAMI3/WACCMX (red), and JPL/GIM (blue). Thick lines and filled circles are for EOF1, while thin lines and open circles are for EOF2. 

In the text 
Supplemental Figure S1. Daily TEC values at latitude = 0°, longitude = 0° for 2010 from SAMI3/climo (black) and JPL/GIM (blue). Panels are for different times of day: (a) 0000 UTC, (b) 0600 UTC, (c) 1200 UTC, and (d) 1800 UTC. 

In the text 
Supplemental Figure S2. (a) Daily values of F10.7 over the year 2010. (b) Daily values of Ap index over 2010. (c) Monthly standard deviation of F10.7 (black line) for 2010 along with the correlation of the 1200 UTC EOF1 amplitude with F10.7 (red line). (d) Monthly standard deviation of Ap index (black line) for 2010 along with the correlation of the SAMI3/WACCMX 1200 UTC EOF1 with Ap (red line). 

In the text 
Supplemental Figure S3. SAMI3/climo and JPL/GIM TEC correlation composites for 2–31 January 2010 as a function of latitude and LT. The reference point is marked by the black dot. It is located on the Equator, and the reference LT varies with panel (0000, 0600, 1200, and 1800 are shown). Thick black lines indicate correlation of 0.5, 0.7, and 0.9. Thin black lines are SZA on 16 January 2010 at 10°, 20°, …, 90°. 

In the text 
Supplemental Figure S4. SAMI3/climo and JPL/GIM line plots of TCC versus relative LT (black) for 2–31 January 2010 at four reference local times (0000, 0600, 1200, and 1800 LT). The reference latitude is the Equator. The Gaussian fits are shown by the dashed blue lines, with the two fitted parameters shown on the bottom of each plot. Base is unitless and width is in hours LT. 

In the text 
Supplemental Figure S5. (a) – (c) Monthly longitudinal base values (red lines) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM along with the monthly for 2014. Red lines are for different reference LTs, and blue line is average over a reference LTs. (d) – (f) Monthly Mlat base values (red lines) from the SAMI3/climo, SAMI3/WACCMX, and JPL/GIM along with monthly . Red lines are for different reference LTs, and blue line is average over a reference LTs. 

In the text 
Supplemental Figure S6. (a) – (c) Longitudinal base values as a function of local time (red lines) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM along with the monthly for 2014. Red lines are for different reference months of 2014, and blue line is average over a reference LTs. (d) – (f) Monthly Mlat base values as a function of local time (red lines) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM along with monthly . Red lines are for different months, and blue line is average over a reference LTs. 

In the text 
Supplemental Figure S7. SAMI3/climo, SAMI3/WACCMX, and JPL/GIM line plots of TCC versus relative LT (black) for 2–31 January 2014 at four reference local times (0000, 0600, 1200, and 1800 LT). The leading EOF has been removed from the TEC time series before calculating the TCCs. The reference latitude is the Equator. The Gaussian fits are shown by the dashed blue lines, with the two fitted parameters shown on the bottom of each plot. Base is unitless and width is in hours LT. 

In the text 
Supplemental Figure S8. (a) – (c) Monthly longitudinal width values (after removing the EOF1) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM for 2014. Red lines indicate different local times, and blue lines are the averages over all local times. (d) – (e) Monthly Mlat width values (after removing the EOF1) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM. Red lines indicate different local times, and blue lines are the averages over all local times. All widths are in great circle degrees. 

In the text 
Supplemental Figure S9. (a) – (c) Longitudinal width values as a function of local time (after removing the EOF1) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM for 2014. Red lines indicate different months, and blue lines are the averages over all months. Mlat width as a function of magnetic local time (after removing the EOF1) from SAMI3/climo, SAMI3/WACCMX, and JPL/GIM. Red lines indicate different months, and blue lines are the averages over all months. All widths are in great circle degrees. 

In the text 
Supplemental Figure S10. (a) – (c) Monthly conjugate correlations at reference latitude 20°N (after removing the EOF1) for SAMI3/climo, SAMI3/WACCMX, and JPL/GIM for 2010. Red lines indicate different local times, and blue lines are the averages over all local times. (d) – (f) Conjugate correlations as a function of magnetic local time (after removing the EOF1) for SAMI3/climo, SAMI3/WACCMX, and JPL/GIM. Red lines indicate different months, and blue lines are the averages over all months. 

In the text 
Supplemental Figure S11. Correlations of the EOF1 and EOF2 amplitudes with F10.7 as a function of time (UTC) and month of 2010. If the correlation is significant (p < 0.050 is 95% confidence level) then a circle is included on the plot. The colors indicate the experiments SAMI3/climo (black) and JPL/GIM (blue). Thick lines and filled circles are for EOF1, while thin lines and open circles are for EOF2. 

In the text 
Supplemental Figure S12. Correlations of the EOF1 and EOF2 amplitudes with Ap as a function of time (UTC) and month of 2010. If the correlation is significant (p < 0.050 is 95% confidence level) then a circle is included on the plot. The colors indicate the experiments SAMI3/climo (black) and JPL/GIM (blue). Thick lines and filled circles are for EOF1, while thin lines and open circles are for EOF2. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.