Swarm Satellites Research Articles

Multispacecraft Current Density Estimates in the Low‐ and Mid‐Latitude F‐Region Ionosphere Using the Swarm Constellation

AbstractWe present a new application of a multispacecraft method which provides estimates of the full electric current density vector in the ionospheric low‐ and mid‐latitude F region. The method uses the three satellites of the Swarm constellation in configurations when they were close to each other. The current density is calculated inside triangular prisms defined by the satellite positions. The technique is inspired by similar approaches such as the curlometer. Here, we propose an alternative mathematical treatment, which involves the use of the well‐known curl‐B technique and a least squares minimization. The proposed formalism provides a rigorous way to propagate the errors through the problem and special care is taken to validate the algorithm and assess the relevance of the derived currents. Such techniques are powerful tools to study the F‐region ionosphere currents as they naturally exclude all the other contributions outside the considered volumes. The F region is known to host currents with a complex spatiotemporal variability. Besides the intrinsic interest of this system itself, understanding the dynamics of these currents is also critical for global geomagnetic field modeling as they are an important source of error in the models. In this paper, we present the full algorithm and test the method on both synthetic and real Swarm satellite data on February 15, 2014. Several sources of error are investigated. Additionally, the results reveal expected characteristics of F‐region interhemispheric field‐aligned currents. Some new characteristics of F‐region currents, potentially associated with dynamo and pressure currents, are also observed.

Assessing the performance of a Northeast Asia Japan-centered 3-D ionosphere specification technique during the 2015 St. Patrick\u2019s day geomagnetic storm

This paper demonstrates and assesses the capability of the advanced three-dimensional (3-D) ionosphere tomography technique, during severe conditions. The study area is northeast Asia and quasi-Japan-centred. Reconstructions are based on total electron content data from a dense ground-based global navigation satellite system receiver network and parameters from operational ionosondes. We used observations from ionosondes, Swarm satellites and radio occultation (RO) to assess the 3-D picture. Specifically, we focus on St. Patrick’s day geomagnetic storm (17–19 March 2015), the most intense in solar cycle 24. During this event, the energy ingested into the ionosphere resulted in Dst and Kp and reaching values ~ − 223 nT and 8, respectively, and the region of interest, the East Asian sector, was characterized by a ~ 60% reduction in electron densities. Results show that the reconstructed densities follow the physical dynamics previously discussed in earlier publications about storm events. Moreover, even when ionosonde data were not available, the technique could still provide a consistent picture of the ionosphere vertical structure. Furthermore, analyses show that there is a profound agreement between the RO profiles/in-situ densities and the reconstructions. Therefore, the technique is a potential candidate for applications that are sensitive to ionospheric corrections.

Neutral Wind Profiles During Periods of Eastward and Westward Equatorial Electrojet

AbstractThe equatorial electrojet (EEJ) is a band of zonal electric current flowing along the magnetic equator in the dayside E‐region ionosphere. The direction of the EEJ is often eastward but sometimes westward. The mechanism for generating westward EEJ is not fully understood. This study examines the relationship between the eastward/westward EEJ and equatorial neutral winds using simultaneous observations of the EEJ from the European Space Agency's Swarm satellite mission and thermospheric winds from the Michelson Interferometer for Global High‐resolution Thermospheric Imaging on NASA's Ionospheric Connection Explorer mission during December 2019–January 2021. Significant differences are found in the average zonal wind profiles between times of eastward and westward EEJ. The EEJ intensity correlates negatively (R = −0.54) and positively (R = 0.58) with the eastward wind velocities at ∼110 and ∼140 km, respectively. The results suggest that the modulation of the zonal electric field by the equatorial zonal wind plays a role in producing the westward EEJ.

Dependence of Parallel Electrical Conductivity in the Topside Ionosphere on Solar and Geomagnetic Activity

AbstractThe spatial and temporal dynamics of ionospheric currents are among the most evident manifestations of the interaction between the magnetized plasma of solar origin and the magnetosphere‐ionosphere system. In this context, a special role is played by field‐aligned currents (FACs), which drive energy and momentum exchanges between the magnetosphere and the ionosphere, deeply affecting the energy balance of the latter. One of the most relevant physical quantities that may help to fully characterize FACs is the ionospheric electrical conductivity. Understanding the features associated with this quantity may contribute to the advancement of knowledge on the mechanisms of solar wind‐magnetosphere‐ionosphere interaction as well as energy storage and dissipation involved in the space weather phenomena. By extending a previous study, we investigate the dependence of the electrical conductivity parallel to the geomagnetic field on solar and geomagnetic activity. To this aim, we considered a six‐year long data set of in‐situ electron density and temperature values recorded by the Langmuir probes on board the Swarm satellites in the topside ionosphere. With this large data set, we computed global maps of the parallel electron conductivity under both quiet and disturbed conditions, and for different solar activity levels. In both cases, the International Reference Ionosphere (IRI) model allowed estimating the contribution of particle precipitation to electrical conductivity.

The Influence of Ionospheric Neutral Wind Variations on the Morphology and Propagation of Medium Scale Traveling Ionospheric Disturbances on 8th August 2016

AbstractThis study reports a special case of the propagation and morphology of medium scale traveling ionospheric disturbances (MSTIDs) over middle‐latitude China on the night of August 8, 2016. The MSTIDs were simultaneously observed by multi‐instruments, including the all‐sky imager, Swarm satellite, and global positioning system (GPS). The MSTIDs lasted for about 6 h in the field view of airglow imager, showing typical wavelength, phase velocity of 272–296 km and 67–250 m/s, respectively. In addition, the imagers show that the inclination angles of phase fronts for some MSTIDs were decreasing during their propagation, resulting in the propagation direction changed from southwestward to nearly westward. More interestingly, the MSTIDs began to dissipate in the airglow observation when they propagated to lower latitudes (below ∼40°N) whereas the MSTIDs at higher latitudes were still visible in the later local times. Simulation results from the Thermosphere‐Ionosphere‐Electrodynamics General Circulation Model are fairly consistent with the Fabry‐Perot Interferometer (FPI) wind observations, which provide convincing explanation to show that the variations of ionospheric neutral winds might play important roles in the changes of propagation direction and the dissipation of MSTIDs.

A two-dimensional lithospheric magnetic anomaly field model of Egypt using the measurements from Swarm satellites

We use magnetic field data observed by the Swarm mission from 2014 to 2020 to construct, for the first time, a two-dimensional (2D) lithospheric magnetic anomaly model of Egypt and its surrounding area. Nighttime data during quiet geomagnetic conditions has been expanded in terms of the Legendre polynomial in harmonic terms N = 6–50. The damped least square method has been used to estimate the model coefficients based on the lithospheric magnetic data. Modeled data at two different altitudes (438–448 km and 503–511 km) were compared with the CHAOS model. Results exhibit that the 2D model is superior to the CHAOS model in the capability of extracting more information about small-scale crustal anomaly field. At low altitudes (438–448 km), the strength of the anomaly field increases, but the noise of the external fields has greatly reduced at high altitudes (503–511 km). Besides, the magnetic anomaly field at low altitudes has illuminated short-scale anomalies that didn't appear at high altitudes. Both the total and vertical magnetic anomaly vectors showed their ability to reveal tectonic structures compared with Moho depth map and the geological maps.

Analysis of Swarm Satellite Magnetic Field Data Before the 2016 Ecuador (Mw = 7.8) Earthquake Based on Non-negative Matrix Factorization

In this paper, based on non-negative matrix factorization (NMF), we analyzed the ionosphere magnetic field data of the Swarm Alpha satellite before the 2016 (Mw = 7. 8) Ecuador earthquake (April 16, 0.35°N, 79.93°W), including the whole data collected under quiet and disturbed geomagnetic conditions. The data from each track were decomposed into basis features and their corresponding weights. We found that the energy and entropy of one of the weight components were more concentrated inside the earthquake-sensitive area, which meant that this weight component was more likely to reflect the activity inside the earthquake-sensitive area. We focused on this weight component and used five times the root mean square (RMS) to extract the anomalies. We found that for this weight component, the cumulative number of tracks, which had anomalies inside the earthquake-sensitive area, showed accelerated growth before the Ecuador earthquake and recovered to linear growth after the earthquake. To verify that the accelerated cumulative anomaly was possibly associated with the earthquake, we excluded the influence of the geomagnetic activity and plasma bubble. Through the random earthquake study and low-seismicity period study, we found that the accelerated cumulative anomaly was not obtained by chance. Moreover, we observed that the cumulative Benioff strain S, which reflected the lithosphere activity, had acceleration behavior similar to the accelerated cumulative anomaly of the ionosphere magnetic field, which suggested that the anomaly that we obtained was possibly associated with the Ecuador earthquake and could be described by one of the Lithosphere–Atmosphere–Ionosphere Coupling (LAIC) models.

A new integrated geophysical-petrological global 3-D model of upper-mantle electrical conductivity validated by the Swarm M2 tidal magnetic field

SUMMARY A new global model of the present-day thermochemical state of the lithosphere and upper mantle based on global waveform inversion, satellite gravity and gradiometry measurements, surface elevation and heat flow data has been recently presented: WINTERC-G (Fullea et al. 2021). WINTERC-G is built within an integrated geophysical-petrological framework where the mantle seismic velocity and density fields are computed in a thermodynamically self-consistent framework, allowing for a direct parametrization in terms of the temperature, pressure and composition of the subsurface rocks. In this paper, we combine WINTERC-G thermal and compositional fields along with laboratory experiments constraining the electrical conductivity of mantle minerals, melt and water, and derive a set of new global three dimensional electrical conductivity models of the upper mantle. The new conductivity models, WINTERC-e, consist of two end-members corresponding to minimum and maximum conductivity of the in situ rock aggregate accounting for mantle melting, mineral water content and the individual conductivities of the main stable mantle mineral phases. The end-member models are validated over oceans by simulating the magnetic field induced by the ocean M2 tidal currents and comparing the predicted fields with the M2 tidal magnetic field estimated from 6-yr Swarm satellite observations. Our new conductivity model, derived independently from any surface or satellite magnetic data sets, is however able to predict tidal magnetic fields that are in good agreement with the Swarm M2 tidal magnetic field models estimated by Sabaka et al. and Grayver & Olsen. Our predicted M2 tidal magnetic fields differ in amplitudes by about 5–20 per cent from the Swarm M2 tidal magnetic field, with the high conductivity WINTERC-e end-member model accounting for mantle melt and water content capturing the structure of Swarm data better than the low conductivity end-member model. Spherically symmetric conductivity models derived by averaging new WINTERC-e conductivities over oceanic areas are slightly more conductive than the recent global conductivity models AA17 by Grayver et al. derived from Swarm and CHAMP satellite data in the 60–140 km depth range, while they are less conductive deeper in the mantle. The conductivities in WINTERC-e are about three to four times smaller than the AA17 conductivities at a depth of 400 km. Despite the differences in electrical conductivity, our spherically symmetric high conductivity end-member model WINTERC-e captures the structure of Swarm M2 tidal magnetic field almost the same as a state of the art 1-D conductivity models derived entirely from magnetic data (AA17, Grayver et al.). Moreover, we show that realistic lateral electrical conductivity inhomogeneities of the oceanic upper mantle derived from the temperature, melt and water distributions in WINTERC-e contribute to the M2 tidal magnetic field up to ±0.3 nT at 430 km altitude.

Collective Space Object as a New concept of International Space Law

This article builds on the premise that international space law should be an enabler, not an inhibitor, of technological advancements in the use of outer space. Examining the features of space stations, space-made products, satellite constellations and satellite swarms, it argues that the time may be mature for a revision of the traditional understanding of the term ‘space object’. Through the evolutionary interpretation of the term ‘component parts’ within Article I of the Liability Convention and Article I of the Registration Convention the article proposes a new concept called ‘collective space object’ and illustrates its practical benefits for the regulation of these innovative technologies. space object, space law, space stations, space-made products, satellite constellations, swarm satellites

Applications for CryoSat-2 satellite magnetic data in studies of Earth\u2019s core field variations

We use 20 years of continuous magnetic field measurements from the Ørsted, CHAMP and Swarm satellite missions, supplemented by calibrated platform magnetometer data from the CryoSat-2 satellite, to study time variations of the Earth’s core field at satellite altitude and at the core–mantle boundary (CMB). From the satellite data we derive composite time series of the core field secular variation (SV) with 4-month cadence, at 300 globally distributed Geomagnetic Virtual Observatories (GVO). A previous gap in the GVO series between 2010 and 2014 is successfully filled using CryoSat-2, and sub-decadal variations are identified during this period. Tests showed that similar sub-decadal SV patterns were obtained from the CryoSat-2 data regardless of whether IGRF-13 or CHAOS-6x9 was used in their calibration. Cryosat-2 radial field SV series at non-polar latitudes have a mean standard deviation level compared to smoothing spline fits of 3.5 nT/yr compared to 1.8 nT/yr for CHAMP and 0.9 nT/yr for Swarm. GVO radial SV series display regional fluctuations with 5–10 years duration and amplitudes reaching 20 nT/yr, most notably at low latitudes over Indonesia (2014), over South America and the South Atlantic (2007, 2011 and 2014), and over the central Pacific (2017). Applying the Subtractive Optimally Localized Averages (SOLA) method, we also map the radial SV at the CMB as a collection of locally averaged SV estimates. We demonstrate that using 2-year windows of CryoSat-2 data, it is possible to reliably estimate the SV and its time derivative, the secular acceleration (SA), at the CMB, with a spatial resolution, corresponding to spherical harmonic degree 10. Along the CMB geographic equator, we find strong SA features with amplitude pm 2.5mu mathrm{T}/mathrm{yr}^2 under Indonesia from 2011–2014, under central America from 2015 to 2019, and sequences of SA with alternating sign under the Atlantic during 2004–2019. We find that platform magnetometer data from CryoSat-2 make a valuable contribution to the emerging picture of sub-decadal core field variations. Using 1-year windows of data from the Swarm satellites, we show that it is possible to study SA changes at low latitudes on timescales down to 1 year, with spatial resolution corresponding to spherical harmonic degree 10. We find strong positive and negative SA features appearing side-by-side in the Pacific in 2017, and thereafter drift westward.

Bandwidth correction of Swarm GPS carrier phase observations for improved orbit and gravity field determination

Gravity fields derived from GPS tracking of the three Swarm satellites have shown artifacts near the geomagnetic equator, where the carrier phase tracking on the L2 frequency is unable to follow rapid ionospheric path delay changes due to a limited tracking loop bandwidth of only 0.25 Hz in the early years of the mission. Based on the knowledge of the loop filter design, an analytical approach is developed to recover the original L2 signal from the observed carrier phase through inversion of the loop transfer function. Precise orbit determination and gravity field solutions are used to assess the quality of the correction. We show that the a posteriori RMS of the ionosphere-free GPS phase observations for a reduced-dynamic orbit determination can be reduced from 3 to 2 mm while keeping up to 7% more data in the outlier screening compared to uncorrected observations. We also show that artifacts in the kinematic orbit and gravity field solution near the geomagnetic equator can be substantially reduced. The analytical correction is able to mitigate the equatorial artifacts. However, the analytical correction is not as successful compared to the down-weighting of problematic GPS data used in earlier studies. In contrast to the weighting approaches, up to 9–10% more kinematic positions can be retained for the heavily disturbed month March 2015 and also stronger signals for gravity field estimation in the equatorial regions are obtained, as can be seen in the reduced error degree variances of the gravity field estimation. The presented approach may also be applied to other low earth orbit missions, provided that the GPS receivers offer a sufficiently high data rate compared to the tracking loop bandwidth, and provided that the basic loop-filter parameters are known.

Ionosonde Observations of Spread F and Spread Es at Low and Middle Latitudes during the Recovery Phase of the 7–9 September 2017 Geomagnetic Storm

This study presents observations of nighttime spread F/ionospheric irregularities and spread Es at low and middle latitudes in the South East Asia longitude of China sectors during the recovery phase of the 7–9 September 2017 geomagnetic storm. In this study, multiple observations, including a chain of three ionosondes located about the longitude of 100°E, Swarm satellites, and Global Navigation Satellite System (GNSS) ROTI maps, were used to study the development process and evolution characteristics of the nighttime spread F/ionospheric irregularities at low and middle latitudes. Interestingly, spread F and intense spread Es were simultaneously observed by three ionosondes during the recovery phase. Moreover, associated ionospheric irregularities could be observed by Swarm satellites and ground-based GNSS ionospheric TEC. Nighttime spread F and spread Es at low and middle latitudes might be due to multiple off-vertical reflection echoes from the large-scale tilts in the bottom ionosphere. In addition, we found that the periods of the disturbance ionosphere are ~1 h at ZHY station, ~1.5 h at LSH station and ~1 h at PUR station, respectively. It suggested that the large-scale tilts in the bottom ionosphere might be produced by LSTIDs (Large scale Traveling Ionospheric Disturbances), which might be induced by the high-latitude energy inputs during the recovery phase of this storm. Furthermore, the associated ionospheric irregularities observed by satellites and ground-based GNSS receivers might be caused by the local electric field induced by LSTIDs.

Statistical Analysis of Pc1 Wave Ducting Deduced From Swarm Satellites

AbstractTransverse Pc1 waves propagating from magnetospheric source regions undergo mode conversion to the compressional mode in the ionosphere due to the induced Hall current. Mode converted Pc1 waves propagate across the magnetic field through the ionospheric waveguide. This process is called Pc1 wave ducting (PWD). PWDs have been observed by magnetometers on both ground and low Earth orbit satellites over a wide latitudinal and longitudinal range. In this work, we present the statistical analysis results of PWD exploiting Swarm satellites from 2015 to 2019. Spatial distributions show that the PWDs are mainly observed over the South Atlantic Anomaly longitudes, possibly due to the high Hall conductivity and F‐region density, and at subauroral/auroral latitudes ( MLAT). The occurrence rate of PWD increases with increasing AE and SYM‐H indices. Seasonal dependence shows that PWD exhibits a high occurrence rate during equinox and local summer while local winter hosts only a low occurrence. The asymmetry between summer and winter can be explained by the ionospheric plasma density. The high occurrence rate in equinox may result from intense geomagnetic activity during the equinox, probably due to the Russell‐McPherron effect. From our statistical analysis, we conclude that the occurrence of PWD is controlled by both ionospheric plasma conditions and geomagnetic activity, and that the mode conversion and PWD occur more efficiently as plasma density increases.

Geomagnetic Virtual Observatories: monitoring geomagnetic secular variation with the Swarm satellites

We present geomagnetic main field and secular variation time series, at 300 equal-area distributed locations and at 490 km altitude, derived from magnetic field measurements collected by the three Swarm satellites. These Geomagnetic Virtual Observatory (GVO) series provide a convenient means to globally monitor and analyze long-term variations of the geomagnetic field from low-Earth orbit. The series are obtained by robust fits of local Cartesian potential field models to along-track and East–West sums and differences of Swarm satellite data collected within a radius of 700 km of the GVO locations during either 1-monthly or 4-monthly time windows. We describe two GVO data products: (1) ‘Observed Field’ GVO time series, where all observed sources contribute to the estimated values, without any data selection or correction, and (2) ‘Core Field’ GVO time series, where additional data selection is carried out, then de-noising schemes and epoch-by-epoch spherical harmonic analysis are applied to reduce contamination by magnetospheric and ionospheric signals. Secular variation series are provided as annual differences of the Core Field GVOs. We present examples of the resulting Swarm GVO series, assessing their quality through comparisons with ground observatories and geomagnetic field models. In benchmark comparisons with six high-quality mid-to-low latitude ground observatories we find the secular variation of the Core Field GVO field intensities, calculated using annual differences, agrees to an rms of 1.8 nT/yr and 1.2 nT/yr for the 1-monthly and 4-monthly versions, respectively. Regular sampling in space and time, and the availability of data error estimates, makes the GVO series well suited for users wishing to perform data assimilation studies of core dynamics, or to study long-period magnetospheric and ionospheric signals and their induced counterparts. The Swarm GVO time series will be regularly updated, approximately every four months, allowing ready access to the latest secular variation data from the Swarm satellites.

Experimental Evidence of a Link Between Lightning and Magnetic Field Fluctuations in the Upper Ionosphere Observed by Swarm

AbstractWe present results of analysis of lightning observations and satellite measurements of ionospheric magnetic field fluctuations in the ultralow frequency (ULF) range. The fluctuations are measured by Swarm satellites, while the lightning observations are provided by the Geostationary Lightning Mapper (GLM) and the World ELF Radiolocation Array (WERA). We identify spatio‐temporal relationships that demonstrate a leakage of electromagnetic fluctuations caused by lightning into the upper ionosphere. Causal association between the two phenomena is evidenced by investigating relations between lightning and fluctuation properties. The presented results suggest that lightning generate ULF fluctuations in the ionosphere that can be detected by satellites, if the lightning‐satellite geographic distance is less than ∼5°. Typical properties of the fluctuations caused by lightning are described in the paper. To our knowledge, this is the first direct experimental confirmation of a link between lightning and magnetic field fluctuations in the upper ionosphere in the ULF range.

The Mag.num core field model as a parent for IGRF-13, and the recent evolution of the South Atlantic Anomaly

We present the GFZ candidate field models for the 13{mathrm{th}} Generation International Geomagnetic Reference Field (IGRF-13). These candidates were derived from the Mag.num.IGRF13 geomagnetic core field model, which is constrained by Swarm satellite and ground observatory data from November 2013 to August 2019. Data were selected from magnetically quiet periods, and the model parameters have been obtained using an iteratively reweighted inversion scheme approximating a robust modified Huber norm as a measure of misfit. The root mean square misfit of the Mag.num.IGRF13 model to Swarm and observatory data is in the order of 3–5 nT for mid and low latitudes, with a maximum of 44 nT for the satellite east component data at high latitudes. The time-varying core field is described by order 6 splines and spherical harmonic coefficients up to degree and order 20. We note that the temporal variation of the core field component of the Mag.num.IGRF13 model is strongly damped and shows a smooth secular variation that suits well for the IGRF, where secular variation is represented as constant over 5-year intervals. Further, the external field is parameterised by a slowly varying part and a more rapidly varying part controlled by magnetic activity and interplanetary magnetic field proxies. Additionally, the Euler angles of the magnetic field sensor orientation are co-estimated. A widely discussed feature of the geomagnetic field is the South Atlantic Anomaly, a zone of weak and decreasing field strength stretching from southern Africa over to South America. The IGRF and Mag.num.IGRF13 indicate that the anomaly has developed a second, less pronounced eastern minimum at Earth’s surface since 2007. We observe that while the strong western minimum continues to drift westwards, the less pronounced eastern minimum currently drifts eastward at Earth’s surface. This does not seem to be linked to any eastward motion at the core–mantle boundary, but rather to intensity changes of westward drifting flux patches contributing to the observed surface field. Also, we report a sudden change in the secular variation measured at two South Atlantic observatories around 2015.0, which occurred shortly after the well-known jerk of 2014.0.

The BGS candidate models for IGRF-13 with a retrospective analysis of IGRF-12 secular variation forecasts

The three candidate models submitted by the British Geological Survey for the 13th generation International Geomagnetic Reference Field are described. These DGRF and IGRF models are derived from vector and scalar magnetic field data from the European Space Agency Swarm satellites and ground observatories, covering the period 2013.9 to 2019.7. The internal field model has time dependence for degrees 1 to 15, represented by order 6 B-splines with knots at six monthly intervals. We also solve for a degree 1 external field time dependence describing annual and semi-annual signals with additional dependence on a bespoke Vector Magnetic Disturbance index. Satellite data are weighted by spatial density, along-track standard deviations, and a larger-scale noise estimator defined in terms of a measure of Local Area Vector Activity at the geographically closest magnetic observatories to the sampled datum. Forecasting of the magnetic field secular variation for 2020–2025 is by advection of the main field using steady core surface flows with steady acceleration applied. We also investigate the performance of the previous generation of candidate secular variation models, for IGRF-12, analysing the agreement of the candidates between 2015 and 2020 with the retrospective IGRF-13. We find that there is no clear distinction between the performance of mathematically and physically extrapolated forecasts in the period 2015–2020. We confirm that the methodology for the BGS IGRF-12 predictions performed well, despite observed secular accelerations that are highlighted by our analysis, and thus justify the methodology used for our IGRF-13 SV candidate.

NOAA/NCEI and University of Colorado candidate models for IGRF-13

The International Geomagnetic Reference Field (IGRF) is a set of parameters representing the large-scale internal part of Earth’s magnetic field. The 13th generation IGRF requested candidate models for a definitive main field for 2015.0, a provisional main field for 2020.0, and a predictive secular variation covering the period 2020.0–2025.0. The University of Colorado (CU) and the National Centers for Environmental Information (NCEI), part of the National Oceanic and Atmospheric Administration (NOAA), have produced these three candidate models for consideration in IGRF-13. In this paper, we present the methodology used to derive our candidate models. Our candidates were built primarily from Swarm satellite data, and also relied on geomagnetic indices derived from the ground observatory network. The ground observatories played a crucial role as independent data in validating our candidates. This paper also provides a retrospective assessment of the CU/NCEI candidate model to the previous IGRF (IGRF-12) and discusses the impact of differences between candidate and final IGRF models on global model errors.

A spherical cap model of the geomagnetic field over southeast Asia from CHAMP and Swarm satellite observations

In this paper, Spherical Cap Harmonic Analysis (SCHA) method was applied to model the geomagnetic field over Vietnam and adjacent area between 15°S and 25°N latitude and 90°E and 130°E in longitude by using magnetic data recorded on CHAMP and Swarm satellites. The characteristic parameters of the method were set at the maximum index Kint = 8 for internal fields, the spherical cap half-angle θ0 = 20°. The regional geomagnetic field over Vietnam and adjacent areas are modelled for the two epochs (2007.0 and 2015.0). Comparison between the SCHA regional geomagnetic field intensity and its time variation with those from IGRF was carried out. The geomagnetic field intensity $$ (E_{F}^{SCHA} ) $$ from SCHA model varies between −90 and 98 nT for epoch 2007.0 and between −139 and 143 nT for epoch 2015.0; however, the trends of their time variations are the same over Vietnam. The RMS between the magnetic components from SCHA model and ground observations are in the same order. The amplitude of time variation of total field intensity from SCHA model is about tens nT greater than from IGRF over Vietnam.

Ground and satellite-based observations of ionospheric plasma bubbles and blobs at 5.65° latitude in the Brazilian sector

This investigation uses simultaneous observations from all-sky imager system and an ionosonde collocated at Araguatins (5.65° S, 48.07° W and dip-latitude of 4.17° S), a near-equatorial region in Brazil. These simultaneous observations were used to investigate the occurrence of plasma bubbles and blobs in the field of the imaging system and their association with atypical range Spread-F signature in ionograms. Also, in-situ observation of plasma density from Swarm satellites were used to support the ground-based observations. Using a few cases, a methodology will be established to identify in the plasma blobs (atypical ESF) in the ionograms when there is the simultaneous observation of plasma bubbles and blobs in the field of view of the ionosonde. For this purpose, simultaneous sequence of OI 630.0 nm nightglow images and ionograms are presented for different case studies; 1. when there is the absence of a plasma bubble or blob, 2. when there is only the occurrence of plasma bubbles and 3. when there is the occurrence of plasma bubbles and blobs, in order to compare traces in the ionogram in all these case studies. With these we can cover all kinds of signatures in the ionograms corresponding to no irregularities, plasma bubbles only and plasma bubbles-blobs. These OI 630.0 nm nightglow and ionograms recorded simultaneously make it possible to establish a novel methodology to recognize in ionograms cases when there is the occurrence of Spread-F signature associated with bubble-blob in the FOV of the ionosonde.