Search

In the studies presented in this thesis, we use the EISCAT UHF incoherent scatter radar (ISR) to study electron precipitation. A new ISR data analysis technique called BAFIM (BAyesian FIltering Module) is developed to calculate plasma parameters (electron density, electron temperature, ion temperature and line of sight ion velocity) with high time and range resolutions from incoherent scatter radar autocorrelation function (ACF) data. BAFIM adds properties of the so-called full-profile analysis to the standard EISCAT data analysis tool, GUISDAP, and extends the concept of full-profile analysis from range direction to both range and time. BAFIM-fitted electron density is used to study a rapidly varying electron precipitation event with high time resolution (4 s). Using a method called ELSPEC, differential number fluxes of precipitating electrons are inverted from electron density altitude profiles measured along the geomagnetic field line by the EISCAT UHF incoherent scatter radar. We show that the raw electron density, that was previously used in high time resolution works, may significantly underestimate the true electron density, when auroral electron precipitation heats the electron gas. The bias affects also electron energy spectra inverted from the raw density profiles, as well as auroral powers and field-aligned currents integrated from the spectra. Temporal variations of the auroral power derived from the fitted electron density show a very good agreement with variations of auroral emission intensity at 427.8 nm. Using more than 20 years of EISCAT UHF radar data, we study statistical characteristics of 1–100 keV electron precipitation at 66.7° magnetic latitude over Tromsø, Norway. Peak energy, auroral power and number flux of electron precipitation are derived from the radar data using the ELSPEC method. We find that 1–5 keV electrons dominate the precipitation from evening until morning in magnetic local time (MLT), while 5–10 keV electrons dominate the late mor

We use the EISCAT incoherent scatter radar data measured in years 2001–2021 to study statistical characteristics of 1–100 keV electron precipitation at 66.7° MLAT over Tromsø. Peak energies, auroral powers and number fluxes of precipitating electrons are derived from electron density altitude profiles measured along the geomagnetic field line during periods of no photoionization. The method allows us to include energetic 30–100 keV electrons, which are poorly covered in earlier satellite-based studies. Locations of the radar within the auroral oval are determined using a model with the 1-hr Hpo geomagnetic index as input. The average peak energy of precipitating electrons increases almost monotonically from evening (18 MLT) to morning hours (09 MLT). The 30–50 keV electrons dominate the energetic electron precipitation before 06 MLT, after which the 50–100 keV precipitation becomes dominant. Large auroral powers (>60 mWm-2) are observed in the 18–02 MLT sector in the main auroral oval. We obtain occurrence rate of electron precipitation in Tromsø by calculating the fraction of data points with auroral power larger than 2 mWm-2. The occurrence rate peaks during the declining phases of solar cycles (sc), in 2002–2004 for sc23 and in 2015–2017 for sc24, caused by variations in geomagnetic activity. In addition, the occurrence rate has maxima during March and September, minimum in December to January, and it increases monotonically from evening to morning hours, reaching maximum at 05–06 MLT.

The article proposes a straightforward Kalman filter-based method for computationally efficient ionospheric electron density multi-instrument imaging. The approach uses direct ionospheric measurements, such as ionosondes, and general physical assumptions to estimate the uncertainty associated with the previous reconstructed time step. Therefore the method does not require any electron density model of the ionosphere as a background. The uncertainty is represented by an inverse covariance matrix constructed with Gaussian Markov random fields, allowing the problem to be solved directly with relatively high resolution. The experiments utilize measurements from dense ground-based GNSS and low Earth orbit beacon satellite receiver networks as well as ionosondes. A synthetic simulation verification and real data validation with a specific European Incoherent Scatter Scientific Association incoherent scatter radar measurement campaign is carried out over Northern European sector. The method can be controlled using parameters with probabilistic and physically realistic interpretations that can be applied to both simulated and real-world data. The results show that the approach is feasible for near real-time regional ionospheric imaging. Especially, the method can be seen as an expansion to local profile measurements field of view, but with sufficient measurement coverage, it also provides information further away from the specific instrument.

Ion drag associated with the ionospheric plasma convection plays an important role in the high-latitude thermospheric dynamics, yet changes in the thermospheric wind with geomagnetic activity are not fully understood. We performed a statistical analysis of the thermospheric wind measurements with a Fabry–Perot interferometer (FPI; 630 nm wavelength) in Tromsø, Norway, in the winter months for 9 years. The measurements were sorted by a SuperMAG (SME) index, and a quiet-time wind pattern was defined as an hourly mean under SME ≤ 40 nT. The quiet-time wind pattern can be expected to be represented by a pressure gradient associated with the solar radiation and a geostrophic force balance. With an increase in the geomagnetic activity level, the thermospheric wind turned over from eastward to westward at dusk and increased the equatorward magnitude from midnight to dawn. Deviations from the quiet-time wind presented similar patterns in the direction with the ionospheric plasma convection but were larger in magnitude at dusk than at dawn. This is the first study to report a dawn-dusk asymmetry of the thermospheric wind acceleration feature and signatures of the eastward wind acceleration at dawn by ion drag.

A thermospheric wind data set from a Fabry-Perot interferometer (630 nm) and the ion velocity from a Dynasonde in Tromsø, Norway, was analyzed for nine winter seasons to study the dynamics of the thermosphere and F-region ionosphere at an auroral latitude. This study focused on bifurcation in the zonal component of the neutral wind and ion velocity at midnight and its dependence on the Y component of the interplanetary magnetic field (IMF). Ionospheric plasma convection patterns are evidently imprinted on the thermospheric wind variations as aspects of the westward and eastward accelerations at dusk and late morning, respectively. The zonal wind bifurcates immediately before midnight for IMF By < 0, but for By > 0, it inverts gradually into the postmidnight sector. Neutral wind streams, originating from higher latitudes, may result in the dependence because of anti-sunward plasma flow distorted in the polar cap.

During geomagnetically active times, Joule (or frictional) heating in the Lower Thermosphere-Ionosphere is a significant source of thermal energy, greatly affecting density, temperature, composition and circulation. At the same time, Joule heating and the associated Pedersen conductivity are amongst the least known parameters in the upper atmosphere in terms of their quantification and spatial distribution, and their parameterization by geomagnetic parameters shows large discrepancies between estimation methodologies, primarily due to a lack of comprehensive measurements in the region where they maximize. In this work we perform a long-term statistical comparison of Joule heating as calculated by the NCAR Thermosphere-Ionosphere-Electrodynamics General Circulation Model (TIE-GCM) and as obtained through radar measurements by the European Incoherent Scatter Scientific Association (EISCAT). Statistical estimates of Joule heating and Pedersen conductivity are obtained from a simulation run over the 11 year period spanning from 2009 until 2019 and from radar measurements over the same period, during times of radar measurements. The results are statistically compared in different Magnetic Local Time sectors and Kp level ranges in terms of median values and percentiles of altitude profiles. It is found that Joule heating and Pedersen conductivity are higher on average in TIE-GCM than in EISCAT for low Kp and are lower than EISCAT for high Kp. It is also found that neutral winds cannot account for the discrepancies between TIE-GCM and EISCAT. Comparisons point toward the need for a Kp-dependent parameterization of Joule heating in TIE-GCM to account for the contribution of small scale effects. During times of high solar activity, Joule (or frictional) heating in the Lower Thermosphere-Ionosphere is a significant source of thermal energy, greatly affecting density, temperature, composition and circulation. Joule heating is largely unknown, due to a lack of measurements in t

Context: The evolution of the photospheric magnetic field can be simulated with surface flux transport (SFT) simulations, which allow for the study of the evolution of the entire field, including polar fields, solely using observations of the active regions. However, because only one side of the Sun is visible at a time, active regions that emerge and decay on the far-side are not observed and not included in the simulations. As a result, some flux is missed. Aims: We construct additional active regions and apply them to the far-side of the Sun in an SFT simulation to assess the possible effects and the magnitude of error that the missing far-side flux causes. We estimate how taking the missing far-side flux into account affects long-term SFT simulations. Methods: We identified active regions from synoptic maps of the photospheric magnetic field between 1975 and 2019. We divided them into solar cycle wings and determined their lifetimes. Using the properties of observed active regions with sufficiently short lifetimes, we constructed additional active regions and inserted them into an SFT simulation. Results: We find that adding active regions with short lifetimes to the far-side of the Sun results in significantly stronger polar fields in minimum times and slightly delayed polarity reversals. These results partly remedy the earlier results, which show overly weak polar fields and polarity reversals that are slightly too early when far-side emergence is not taken into account. The far-side active regions do not significantly affect poleward flux surges, which are mostly caused by larger long-living active regions. The far-side emergence leads to a weak continuous flow of flux, which affects polar fields over long periods of time.

This study presents an improved method to estimate differential energy flux, auroral power and field-aligned current of electron precipitation from incoherent scatter radar data. The method is based on a newly developed data analysis technique that uses Bayesian filtering to fit altitude profiles of electron density, electron temperature, and ion temperature to observed incoherent scatter spectra with high time and range resolutions. The electron energy spectra are inverted from the electron density profiles. Previous high-time resolution fits have relied on the raw electron density, which is calculated from the backscattered power assuming that the ion and electron temperatures are equal. The improved technique is applied to one auroral event measured by the EISCAT UHF radar and it is demonstrated that the effect of electron heating on electron energy spectra, auroral power, and upward field-aligned current can be significant at times. Using the fitted electron densities instead of the raw ones may lead to wider electron energy spectra and auroral power up to 75% larger. The largest differences take place for precipitation that produces enhanced electron heating in the upper E region, and in this study correspond to fluxes of electrons with peak energies from 3 to 5 keV. Finally, the auroral power estimates are verified by comparison to the 427.8 nm auroral emission intensity, which shows good correlation. The improved method makes it possible to calculate unbiased estimates of electron energy spectra with high time resolution and thereby to study rapidly varying aurora.

We study historic observations of solar activity from the 20th-century rise towards the peak of the Modern Grand Solar Maximum (MGSM) and compare with observations of the decline that has occurred since. The major difference in available solar observations of the rise and of the fall are accurate magnetograms from solar magnetographs: we here use synthetic magnetograms to interpret the rise and employ historic observations of Polar Crown Filaments to test them and verify their use. We show that eclipse images at sunspot minimum reveal the long-term variation of open flux deduced from geomagnetic observations in Paper 1 (Lockwood et al., 2022). We also make use of polar coronal hole fluxes derived from historic white light images of polar faculae, but have to consider the implications of the fact that these facular images do not tell us the polarity of the field. Given this caveat, the agreement between the polar coronal hole fluxes and the values derived from open flux continuity modelling based on sunspot numbers is extremely good. This comparison indicates that one possible solution to the “open flux problem” is open flux within the streamer belt that potential-based modelling of coronal fields from photospheric fields is not capturing. We take a detailed look at the solar cycle at the peak of the MGSM, cycle 19, and show the variation of the polar coronal hole fluxes and the inferred poleward flux surges are predictable from the asymmetries in flux emergence in the two hemispheres with implied transequatorial flux transfer and/or “anti-Hale” (or more general “rogue” active region flux) emergence late in the sunspot cycle.

In the space physics community, processing and combining observational and modeling data from various sources is a demanding task because they often have different formats and use different coordinate systems. The Python package GeospaceLAB has been developed to provide a unified, standardized framework to process data. The package is composed of six core modules, including DataHub as the data manager, Visualization for generating publication quality figures, Express for higher-level interfaces of DataHub and Visualization, SpaceCoordinateSystem for coordinate system transformations, Toolbox for various utilities, and Configuration for preferences. The core modules form a standardized framework for downloading, storing, post-processing and visualizing data in space physics. The object-oriented design makes the core modules of GeospaceLAB easy to modify and extend. So far, GeospaceLAB can process more than twenty kinds of data products from nine databases, and the number will increase in the future. The data sources include, e.g., measurements by EISCAT incoherent scatter radars, DMSP, SWARM, and Grace satellites, OMNI solar wind data, and GITM simulations. In addition, the package provides an interface for the users to add their own data products. Hence, researchers can easily collect, combine, and view multiple kinds of data for their work using GeospaceLAB. Combining data from different sources will lead to a better understanding of the physics of the studied phenomena and may lead to new discoveries. GeospaceLAB is an open source software, which is hosted on GitHub. We welcome everyone in the community to contribute to its future development.

Context: We apply our recently developed method to reconstruct synoptic maps of the photospheric magnetic field from observations of chromospheric plages and the magnetic polarity of sunspots. Here, we apply the method to an extended time interval from 1915 to 1985. Aims: Systematic magnetographic observations of the solar photospheric magnetic field were initiated as recently as the 1970s and the lack of earlier observations limits our ability to study and understand the long-term evolution of the Solar global field. This study is aimed at creating synoptic maps of magnetic fields for the pre-magnetograph era and using these maps as input for modern simulation models to investigate the long-term (centennial) evolution of the Sun’s global magnetic fields. Methods: We reconstructed active Solar regions by identifying chromospheric plages from Ca II K line synoptic maps and assigning magnetic polarities based on the observed polarity of sunspots. We used a surface flux transport (SFT) model to simulate the evolution of the photospheric magnetic field from the reconstructed active regions. We used the potential field source surface (PFSS) model to determine the amount of open magnetic flux from the reconstruction and from magnetographic observations. We also reconstructed the coronal field during two eclipses and compared the result with eclipse drawings. Results: We successfully reconstructed the photospheric magnetic field from 1915 to 1985. The number and total magnetic flux of the reconstructed active regions shows a realistic cyclic behavior that mostly follows the evolution of the sunspot number, even on relatively short timescales. The polar field strengths of cycles 19 and 20 do not reflect the evolution of the sunspot number very accurately, which may be related to problems related to the calcium data during cycle 19 and the long data gap during cycle 20. The polarity of polar fields and the amount of open field both at high and low latitudes all demo

Context: The bright regions in the solar chromosphere and temperature minimum have a good spatial correspondence with regions of intense photospheric magnetic field. Bright regions are visible in different emission lines and parts of the continuum. Their observation started more than a hundred years ago with the invention of the spectroheliograph. While the historical spectroheliograms are essential for studying the long-term variability of the Sun, the modern satellite-borne observations can help us reveal the nature of chromospheric brightenings in previously unattainable detail. Aims: Our aim is to improve the understanding of the relation between magnetic fields and radiative structures byf studying modern seeing-free observations of far-ultraviolet (FUV) radiation around 1600 Å and photospheric magnetic fields. Methods: We used Helioseismic and Magnetic Imager (HMI) observations of photospheric magnetic fields and Atmospheric Imaging Assembly (AIA) observations of FUV contrast around 1600 Å. We developed a robust method to find contrast thresholds defining bright and dark AIA 1600 Å pixels, and we combine them to bright and dark clusters. We investigate the relation of magnetic fields and AIA 1600 Å radiation in bright and dark clusters. Results: We find that the percentage of bright pixels (ranging from 2% to 10%) almost entirely explains the observed variability of 1600 Å emission. We developed a multilinear regression model based on the percentages of bright and dark pixels, which can reliably predict the magnitude of the disk-averaged unsigned magnetic field. We find that bright and dark clusters closely correspond respectively to the populations of moderate ( B > 55 G) and strong (B > 1365 G) magnetic field HMI clusters. The largest bright clusters have a constant mean unsigned magnetic field, as found previously for Ca II K plages. However, the magnetic field strength of bright clusters is 254.7 ± 0.1 G, which is roughly 100 G larger than f

We propose a novel method for deconvolving incoherent scatter radar data to recover accurate reconstructions of backscattered powers. The problem is modelled as a hierarchical noise-perturbed deconvolution problem, where the lower hierarchy consists of an adaptive length-scale function that allows for a non-stationary prior and as such enables adaptive recovery of smooth and narrow layers in the profiles. The estimation is done in a Bayesian statistical inversion framework as a two-step procedure, where hyperparameters are first estimated by optimisation and followed by an analytical closed-form solution of the deconvolved signal. The proposed optimisation-based method is compared to a fully probabilistic approach using Markov chain Monte Carlo techniques enabling additional uncertainty quantification. In this paper we examine the potential of the hierarchical deconvolution approach using two different prior models for the length-scale function. We apply the developed methodology to compute the backscattered powers of measured polar mesospheric winter echoes, as well as summer echoes, from the EISCAT VHF radar in Tromsø, Norway. Computational accuracy and performance are tested using a simulated signal corresponding to a typical background ionosphere and a sporadic E layer with known ground truth. The results suggest that the proposed hierarchical deconvolution approach can recover accurate and clean reconstructions of profiles, and the potential to be successfully applied to similar problems.

Context: The axial dipole moments of emerging active regions control the evolution of the axial dipole moment of the whole photospheric magnetic field and the strength of polar fields. Hale’s and Joy’s laws of polarity and tilt orientation affect the sign of the axial dipole moment of an active region. If both laws are valid (or both violated), the sign of the axial moment is normal. However, for some active regions, only one of the two laws is violated, and the signs of these axial dipole moments are the opposite of normal. Those opposite-sign active regions can have a significant effect, for example, on the development of polar fields. Aims: Our aim is to determine the axial dipole moments of active regions identified from magnetographic observations and study how the axial dipole moments of normal and opposite signs are distributed in time and latitude in solar cycles 21−24. Methods: We identified active regions in the synoptic maps of the photospheric magnetic field measured at the National Solar Observatory (NSO) Kitt Peak (KP) observatory, the Synoptic Optical Long term Investigations of the Sun (SOLIS) vector spectromagnetograph (VSM), and the Helioseismic and Magnetic Imager (HMI) aboard the Solar Dynamics Observatory (SDO), and determined their axial dipole moments. Results: We find that, typically, some 30% of active regions have opposite-sign axial moments in every cycle, often making more than 20% of the total axial dipole moment. Most opposite-signed moments are small, but occasional large moments, which can affect the evolution of polar fields on their own, are observed. Active regions with such a large opposite-sign moment may include only a moderate amount of total magnetic flux. We find that in cycles 21−23 the northern hemisphere activates first and shows emergence of magnetic flux over a wider latitude range, while the southern hemisphere activates later, and emergence is concentrated to lower latitudes. Cycle 24 differs from cycles 21−23

Aims: The evolution of the photospheric magnetic field has only been regularly observed since the 1970s. The absence of earlier observations severely limits our ability to understand the long-term evolution of solar magnetic fields, especially the polar fields that are important drivers of space weather. Here, we test the possibility to reconstruct the large-scale solar magnetic fields from Ca II K line observations and sunspot magnetic field observations, and to create synoptic maps of the photospheric magnetic field for times before modern-time magnetographic observations. Methods: We reconstructed active regions from Ca II K line synoptic maps and assigned them magnetic polarities using sunspot magnetic field observations. We used the reconstructed active regions as input in a surface flux transport simulation to produce synoptic maps of the photospheric magnetic field. We compared the simulated field with the observed field in 1975−1985 in order to test and validate our method. Results: The reconstruction very accurately reproduces the long-term evolution of the large-scale field, including the poleward flux surges and the strength of polar fields. The reconstruction has slightly less emerging flux because a few weak active regions are missing, but it includes the large active regions that are the most important for the large-scale evolution of the field. Although our reconstruction method is very robust, individual reconstructed active regions may be slightly inaccurate in terms of area, total flux, or polarity, which leads to some uncertainty in the simulation. However, due to the randomness of these inaccuracies and the lack of long-term memory in the simulation, these problems do not significantly affect the long-term evolution of the large-scale field.

For 2 decades, meteor radars have been routinely used to monitor atmospheric temperature around 90 km altitude. A common method, based on a temperature gradient model, is to use the height dependence of meteor decay time to obtain a height-averaged temperature in the peak meteor region. Traditionally this is done by fitting a linear regression model in the scattered plot of log₁₀(1/τ) and height, where τ is the half-amplitude decay time of the received signal. However, this method was found to be consistently biasing the slope estimate. The consequence of such a bias is that it produces a systematic offset in the estimated temperature, thus requiring calibration with other co-located measurements. The main reason for such a biasing effect is thought to be due to the failure of the classical regression model to take into account the measurement error in τ and the observed height. This is further complicated by the presence of various geophysical effects in the data, as well as observational limitation in the measuring instruments. To incorporate various error terms in the statistical model, an appropriate regression analysis for these data is the errors-in-variables model. An initial estimate of the slope parameter is obtained by assuming symmetric error variances in normalised height and log₁₀(1/τ). This solution is found to be a good prior estimate for the core of this bivariate distribution. Further improvement is achieved by defining density contours of this bivariate distribution and restricting the data selection process within higher contour levels. With this solution, meteor radar temperatures can be obtained independently without needing any external calibration procedure. When compared with co-located lidar measurements, the systematic offset in the estimated temperature is shown to have reduced to 5 % or better on average.

Context: Weak magnetic field elements make a dominant contribution to the total magnetic field on the solar surface. Even so, little is known of their long-term occurrence. Aims: We study the long-term spatial–temporal evolution of the weak-field shift and skewness of the distribution of photospheric magnetic field values during solar cycles 21−24 in order to clarify the role and relation of the weak field values to the overall magnetic field evolution. Methods: We used Wilcox Solar Observatory (WSO) and the Synoptic Optical Long-term Investigations of the Sun Vector SpectroMagnetograph synoptic maps to calculate weak-field shifts for each latitude bin of each synoptic map, and thereby constructed a time–latitude butterfly diagram for shifts. We also calculated butterfly diagrams for skewness for all field values and for weak field values only. Results: The weak-field shifts and (full-field) skewness depict a similar spatial–temporal solar cycle evolution to that of the large-scale surface magnetic field. The field distribution has a systematic non-zero weak-field shift and a large skewness already at (and after) the emergence of the active region, even at the highest resolution. We find evidence for coalescence of opposite-polarity fields during the surge evolution. This is clearly more effective at the supergranulation scale. However, a similar dependence of magnetic field coalescence on spatial resolution was not found in the unipolar regions around the poles. Conclusions: Our results give evidence for the preference of even the weakest field elements toward the prevailing magnetic polarity since the emergence of an active region, and for a systematic coalescence of stronger magnetic fields of opposite polarities to produce weak fields during surge evolution and at the poles. We also find that the supergranulation process is reduced or turned off in the unipolar regions around the poles. These observations improve the understanding not only of the develo

Full disk vector magnetic fields are used widely for developing better understanding of large-scale structure, morphology, and patterns of the solar magnetic field. The data are also important for modeling various solar phenomena. However, observations of vector magnetic fields have one important limitation that may affect the determination of the true magnetic field orientation. This limitation stems from our ability to interpret the differing character of the Zeeman polarization signals which arise from the photospheric line-of-sight vs. the transverse components of the solar vector magnetic field, and is likely exacerbated by unresolved structure (non-unity fill fraction) as well as the disambiguation of the 180° degeneracy in the transverse-field azimuth. Here we provide a description of this phenomenon, and discuss issues, which require additional investigation.

Incoherent scatter (IS) radars are invaluable instruments for ionospheric physics, since they observe altitude profiles of electron density (Ne), electron temperature (Te), ion temperature (Ti) and line‐of‐sight plasma velocity (Vi) from ground. However, the temperatures can be fitted to the observed IS spectra only when the ion composition is known, and resolutions of the fitted plasma parameters are often insufficient for auroral electron precipitation, which requires high resolutions in both range and time. The problem of unknown ion composition has been addressed by means of the full‐profile analysis, which assumes that the plasma parameter profiles are smooth in altitude, or follow some predefined shape. In a similar manner, one could assume smooth time variations, but this option has not been used in IS analysis. We propose a plasma parameter fit technique based on Bayesian filtering, which we have implemented as an additional Bayesian Filtering Module (BAFIM) in the Grand Unified Incoherent Scatter Design and Analysis Package (GUISDAP). BAFIM allows us to control gradients in both time and range directions for each plasma parameter separately. With BAFIM we can fit F1 region ion composition together with Ne, Te, Ti and Vi, and we have reached 4 s/900 m time/range steps in four‐parameter fits of Ne, Te, Ti and Vi in E region observations of auroral electron precipitation.

Aims: Sunspot activity is often hemispherically asymmetric, and during the Maunder minimum, activity was almost completely limited to one hemisphere. In this work, we use surface flux simulation to study how magnetic activity limited only to the southern hemisphere affects the long-term evolution of the photospheric magnetic field in both hemispheres. The key question is whether sunspot activity in one hemisphere is enough to reverse the polarity of polar fields in both hemispheres. Methods: We simulated the evolution of the photospheric magnetic field from 1978 to 2016 using the observed active regions of the southern hemisphere as input. We studied the flow of magnetic flux across the equator and its subsequent motion towards the northern pole. We also tested how the simulated magnetic field is changed when the activity of the southern hemisphere is reduced. Results: We find that activity in the southern hemisphere is enough to reverse the polarity of polar fields in both hemispheres by the cross-equatorial transport of magnetic flux. About 1% of the flux emerging in the southern hemisphere is transported across the equator, but only 0.1%–0.2% reaches high latitudes to reverse and regenerate a weak polar field in the northern hemisphere. The polarity reversals in the northern hemisphere are delayed compared to the southern hemisphere, leading to a quadrupole Sun lasting for several years.

Various subauroral optical features have been studied by analyzing data collected during periods of geomagnetic disturbances. Most events have been typically found at geomagnetic latitudes of 45-60 degrees. In this study, however, we present a red arc event found at geomagnetic 68 degrees north (L approximate to 7.1) in the Scandinavian sector during a period of geomagnetically quiet conditions within a short intermission between two high-speed solar wind events. The red arc appeared to coincide with a pseudo breakup at geomagnetic 71-72 degrees N and a rapid equatorward expansion of the polar cap. However, the red arc disappeared in approximately 7 min. Simultaneous measurements with the Swarm A/C satellites indicated the appearance of the red arc at the ionospheric trough minimum and a conspicuous enhancement of the electron temperature, suggesting the generation of the arc by heat flux. Since there are meaningful differences in the red arc features from already-known subauroral optical features such as the stable auroral red (SAR) arc, we considered that the red arc is a new phenomenon. We suggest that the ephemeral red arc may represent the moment of SAR arc birth associated with substorm particle injection, which is generally masked by bright dynamic aurorae., QC 20210125

We derive the longest uniform record of rotational intensities solar coronal magnetic field since 1968 and compare it with the heliospheric magnetic field (HMF) observed at the Earth. We scale the Mount Wilson Observatory and Wilcox Solar Observatory observations of the photospheric magnetic field to the level of the Synoptic Optical Long-term Investigations of the Sun/Vector Spectro Magnetograph and apply the potential field source surface model to calculate the coronal magnetic field. We find that the evolution of the coronal magnetic field during the last 50 yr agrees with the HMF observed at the Earth only if the effective coronal size, the distance of the coronal source surface of the HMF, is allowed to change in time. We calculate the optimum source surface distance for each rotation and find that it experienced an abrupt decrease in the late 1990s. The effective volume of the solar corona shrunk to less than one half during a short period of only a few years. We note that this abrupt shrinking coincides with other changes in solar magnetic fields that are likely related to the decrease of the overall solar activity, i.e., the demise of the Grand Modern Maximum.