Source author record

Huw Morgan

Huw Morgan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

A solar-cycle study of coronal rotation: large variations, rapid changes, and implications for solar wind models

Information on the rotation rate of the corona, and its variation over latitude and solar cycle, is valuable for making global connections between the corona and the Sun, for global estimates of reconnection rates, and as a basic parameter for solar wind modelling. Here, we use a time series of tomographical maps gained from coronagraph observations between 2007 - 2020 to directly measure the longitudinal drift of high-density streamers over time. The method reveals abrupt changes in rotation rates, revealing a complex relationship between the coronal rotation and the underlying photosphere. The majority of rates are between -1.0 to +0.5$^\circ$/day relative to the standard Carrington rate of 14.18$^\circ$/day, although rates are measured as low as -2.2$^\circ$/day and as high as 1.6$^\circ$/day. Equatorial rotation rates during the 2008 solar minimum are slightly faster than the Carrington rate, with an abrupt switch to slow rotation in 2009, then a return to faster rates in 2017. Abrupt changes and large variations in rates are seen at all latitudes. Comparison with a magnetic model suggests that periods of equatorial fast rotation are associated with times when a large proportion of the magnetic footpoints of equatorial streamers are near the equator, and we interpret the abrupt changes in terms of the latitudinal distribution of the streamer photospheric footpoints. The coronal rotation rate is a key parameter for solar wind models, and variations of up to a degree per day or more can lead to large systematic errors over forecasting periods of longer than a few days. The approach described in this paper gives corrected values that can form a part of future forecasting efforts.

preprint2022arXiv

An improved method for estimating the velocity field of coronal propagating disturbances

The solar corona is host to a continuous flow of propagating disturbances (PD). These are continuous and ubiquitous across broad regions of the corona, including the quiet Sun. The aim of this paper is to present an improved, efficient method to create velocity vector field maps, based on the direction and magnitude of the PD as observed in time series of extreme ultraviolet (EUV) images. The method is presented here for use with the Atmospheric Imaging Assembly (AIA)/Solar Dynamics Observatory (SDO) EUV channels, and takes as input \app2 hours of images at the highest 12s cadence. Data from a region near disk center is extracted, and a process called time normalization applied to the co-aligned data. Following noise reduction using \atrous\ decomposition, the PD are effectively revealed. A modified Lucas Kanade algorithm is then used to map the velocity field. The method described here runs comfortably on a desktop computer in a few minutes, and offers an order of magnitude improvement in efficiency compared to a previous implementation. Applied to a region of the quiet Sun, we find that the velocity field describes a mosaic of cells of coherent outwardly-diverging PD flows, of typical size 50 to 100\arcsec\ (36 to 72Mm). The flows originate from points and narrow corridors in the cell centres, and end in the narrow boundaries between cells. Visual comparison with ultraviolet AIA images shows that the flow sources are correlated with the bright photospheric supergranular network boundaries. Assuming that the PD follow the local magnetic field, the velocity flow field is a proxy for the plane-of-sky distribution of the coronal magnetic field, and therefore the maps offer a unique insight into the topology of the corona. These are particularly valuable for quiet Sun regions where the appearance of structures in EUV images is hard to interpret.

preprint2022arXiv

An inner boundary condition for solar wind models based on coronal density

Accurate forecasting of the solar wind has grown in importance as society becomes increasingly dependent on technology that is susceptible to space weather events. This work describes an inner boundary condition for ambient solar wind models based on tomography maps of the coronal plasma density gained from coronagraph observations, providing a novel alternative to magnetic extrapolations. The tomographical density maps provide a direct constraint of the coronal structure at heliocentric distances of 4 to 8Rs, thus avoiding the need to model the complex non-radial lower corona. An empirical inverse relationship converts densities to solar wind velocities which are used as an inner boundary condition by the Heliospheric Upwind Extrapolation (HUXt) model to give ambient solar wind velocity at Earth. The dynamic time warping (DTW) algorithm is used to quantify the agreement between tomography/HUXt output and in situ data. An exhaustive search method is then used to adjust the lower boundary velocity range in order to optimize the model. Early results show up to a 32% decrease in mean absolute error between the modelled and observed solar wind velocities compared to that of the coupled MAS/HUXt model. The use of density maps gained from tomography as an inner boundary constraint is thus a valid alternative to coronal magnetic models, and offers a significant advancement in the field given the availability of routine space-based coronagraph observations.

preprint2022arXiv

Tracing the magnetic field topology of the quiet corona using propagating disturbances

The motion of faint propagating disturbances (PD) in the solar corona reveals an intricate structure which must be defined by the magnetic field. Applied to quiet Sun observations by the Atmospheric Imaging Assembly (AIA)/Solar Dynamics Observatory (SDO), a novel method reveals a cellular network, with cells of typical diameters 50\arcsec\ in the cool 304Å channel, and 100\arcsec\ in the coronal 193Å channel. The 193Å cells can overlie several 304Å cells, although both channels share common source and sink regions. The sources are points, or narrow corridors, of divergence that occupy the centres of cells. They are significantly aligned with photospheric network features and enhanced magnetic elements. This shows that the bright network is important to the production of PDs, and confirms that the network is host to the source footpoint of quiet coronal loops. The other footpoint, or the sinks of the PDs, form the boundaries of the coronal cells. These are not significantly aligned with the photospheric network - they are generally situated above the dark internetwork photosphere. They form compact points or corridors, often without an obvious signature in the underlying photosphere. We argue that these sink points can either be concentrations of closed field footpoints associated with minor magnetic elements in the internetwork, or concentrations of upward-aligned open field. The link between the coronal velocity and magnetic fields is strengthened by a comparison with a magnetic extrapolation, which shows several general and specific similarities, thus the velocity maps offer a valuable additional constraint on models.

preprint2020arXiv

The width, density and outflow of solar coronal streamers

Characterising the large-scale structure and plasma properties of the inner corona is crucial to understand the source and subsequent expansion of the solar wind and related space weather effects. Here we apply a new coronal rotational tomography method, along with a method to narrow streamers and refine the density estimate, to COR2A/STEREO observations from a period near solar minimum and maximum, gaining density maps for heights between 4 and 8\Rs. The coronal structure is highly radial at these heights, and the streamers are very narrow, in some regions only a few degrees in width. The mean densities of streamers is almost identical between solar minimum and maximum. However, streamers at solar maximum contain around 50\%\ more total mass due to their larger area. By assuming a constant mass flux, and constraints on proton flux measured by Parker Solar Probe (PSP), we estimate an outflow speed within solar minimum streamers of 50-120\kms\ at 4\Rs, increasing to 90-250\kms\ at 8\Rs. Accelerations of around 6\mss\ are found for streamers at a height of 4\Rs, decreasing with height. The solar maximum slow wind shows a higher acceleration to extended distances compared to solar minimum. To satisfy the solar wind speeds measured by PSP, there must be a mean residual acceleration of around 1-2\mss\ between 8 and 40\Rs. Several aspects of this study strongly suggest that the coronal streamer belt density is highly variable on small scales, and that the tomography can only reveal a local spatial and temporal average.

preprint2019arXiv

Spectral diagnostics of cool flare loops observed by SST: I. Inversion of the Ca II 8542 Å and H$β$ lines

Flare loops form an integral part of eruptive events, being detected in the range of temperatures from X-rays down to cool chromospheric-like plasmas. While the hot loops are routinely observed by the Solar Dynamics Observatory's Atmospheric Imaging Assembly (SDO/AIA), cool loops seen off-limb are rare. In this paper we employ unique observations of the SOL2017-09-10T16:06 X8.2-class flare which produced an extended arcade of loops. The Swedish 1-m Solar Telescope (SST) made a series of spectral images of the cool off-limb loops in the Ca II 8542 Å and the hydrogen H$β$ lines. Our focus is on the loop apices. Non-LTE spectral inversion is achieved through the construction of extended grids of models covering a realistic range of plasma parameters. The Multilevel Accelerated Lambda Iterations (MALI) code solves the non-LTE radiative-transfer problem in a 1D externally-illuminated slab, approximating the studied loop segment. Inversion of the Ca II 8542 Å and H$β$ lines yields two similar solutions, both indicating high electron densities around $2 \times 10^{12}$ cm$^{-3}$ and relatively large microturbulence around 25 kms$^{-1}$. These are in reasonable agreement with other independent studies of the same or similar events. In particular, the high electron densities in the range $10^{12} - 10^{13}$ cm$^{-3}$ are consistent with those derived from the SDO's Helioseismic and Magnetic Imager white-light observations. The presence of such high densities in solar eruptive flares supports the loop interpretation of the optical continuum emission of stars which manifest superflares.

preprint2016arXiv

Automated detection of coronal mass ejections in three-dimensions using multi-viewpoint observations

A new, automated method of detecting coronal mass ejections (CMEs) in three dimensions for the LASCO C2 and STEREO COR2 coronagraphs is presented. By triangulating isolated CME signal from the three coronagraphs over a sliding window of five hours, the most likely region through which CMEs pass at 5 solar radii is identified. The centre and size of the region gives the most likely direction of propagation and approximate angular extent. The Automated CME Triangulation (ACT) method is tested extensively using a series of synthetic CME images created using a wireframe flux rope density model, and on a sample of real coronagraph data; including halo CMEs. The accuracy of the angular difference between the detection and true input of the synthetic CMEs is 7.14 degrees, and remains acceptable for a broad range of CME positions relative to the observer, the relative separation of the three observers and even through the loss of one coronagraph. For real data, the method gives results that compare well with the distribution of low coronal sources and results from another instrument and technique made further from the Sun. The true three dimension (3D)-corrected kinematics and mass/density are discussed. The results of the new method will be incorporated into the CORIMP database in the near future, enabling improved space weather diagnostics and forecasting.

preprint2016arXiv

Erupting filaments with large enclosing flux tubes as sources of high-mass 3-part CMEs, and erupting filaments in the absence of enclosing flux tubes as sources of low-mass unstructured CMEs

The 3-part appearance of many CMEs arising from erupting filaments emerges from a large magnetic flux tube structure, consistent with the form of the erupting filament system. Other CMEs arising from erupting filaments lack a clear 3-part structure and reasons for this have not been researched in detail. This paper aims to further establish the link between CME structure and the structure of the erupting filament system and to investigate whether CMEs which lack a 3-part structure have different eruption characteristics. A survey is made of 221 near-limb filament eruptions observed from 2013/05/03-2014/06/30 by EUV imagers and coronagraphs. 92 filament eruptions are associated with 3-part structured CMEs, 41 eruptions are associated with unstructured CMEs. The remaining 88 are categorized as failed eruptions. For 34% of the 3-part CMEs, processing applied to EUV images reveals the erupting front edge is a pre-existing loop structure surrounding the filament, which subsequently erupts with the filament to form the leading bright front edge of the CME. This connection is confirmed by a flux-rope density model. Furthermore, the unstructured CMEs have a narrower distribution of mass compared to structured CMEs, with total mass comparable to the mass of 3-part CME cores. This study supports the interpretation of 3-part CME leading fronts as the outer boundaries of a large pre-existing flux tube. Unstructured (non 3-part) CMEs are a different family to structured CMEs, arising from the eruption of filaments which are compact flux tubes in the absence of a large system of enclosing closed field.

preprint2015arXiv

An atlas of coronal electron density at 5Rs I: Data processing and calibration

Tomography of the solar corona can provide cruicial constraints for models of the low corona, unique information on changes in coronal structure and rotation rates, and a valuable boundary condition for models of the heliospheric solar wind. This is the first of a series of three papers which aim to create a set of maps of the coronal density over an extended period (1996-present). The papers will describe the data processing and calibration (this paper), the tomography method (\paperii) and resulting atlas of coronal electron density at a height of 5\Rs\ between years 1996-2014 (\paperiii). This first paper presents a detailed description of data processing and calibration for the Large-Angle and Spectrometric Coronagraph (LASCO) C2 instrument onboard the Solar and Heliospheric Observatory (SOHO) and the COR2 instruments of the Sun Earth Connection Coronal and Heliospheric Investigation (SECCHI) package aboard the Solar Terrestial Relations Observatory (STEREO) A \& B spacecraft. The methodology includes noise suppression, background subtraction, separation of large dynamic events, conversion of total brightness to K-coronal brightness and simple functions for crosscalibration between C2/LASCO and COR2/SECCHI. Comparison of the brightness of stars between LASCO C2 total and polarized brightness (\pB) observations provide in-flight calibration factors for the \pB\ observations, resulting in considerable improved agreement between C2 and COR2 A, and elimination of curious artifacts in the C2 \pB\ images. The crosscalibration between LASCO C2 and the STEREO coronagraphs allows, for the first time, the potential use of multi-spacecraft coronagraph data for tomography and for CME analysis.

preprint2014arXiv

Bridging EUV and white-light observations to inspect the initiation phase of a "two-stage" solar eruptive event

The initiation phase of CMEs is a very important aspect of solar physics, as these phenomena ultimately drive space weather in the heliosphere. This phase is known to occur between the photosphere and low corona, where many models introduce an instability and/or magnetic reconnection that triggers a CME, often with associated flaring activity. To this end, it is important to obtain a variety of observations of the low corona in order to build as clear a picture as possible of the dynamics that occur therein. Here, we combine the EUV imagery of the SWAP instrument on board PROBA2 with the white-light imagery of the ground-based Mk4 coronameter at MLSO in order to bridge the observational gap that exists between the disk imagery of AIA on board SDO and the coronal imagery of LASCO on board SOHO. Methods of multiscale image analysis were applied to the observations to better reveal the coronal signal while suppressing noise and other features. This allowed an investigation into the initiation phase of a CME that was driven by a rising flux rope structure from a "two-stage" flaring active region underlying an extended helmet streamer. It was found that the initial outward motion of the erupting loop system in the EUV observations coincided with the first X-ray flare peak, and led to a plasma pile-up of the white-light CME core material. The characterized CME core then underwent a strong jerk in its motion, as the early acceleration increased abruptly, simultaneous with the second X-ray flare peak. The overall system expanded into the helmet streamer to become the larger CME structure observed in the LASCO coronagraph images, which later became concave-outward in shape. Theoretical models for the event are discussed in light of these unique observations, and it is concluded that the formation of either a kink-unstable or torus-unstable flux rope may be the likeliest scenario.

preprint2014arXiv

Multi-scale Gaussian normalization for solar image processing

Extreme UltraViolet images of the corona contain information over a large range of spatial scales, and different structures such as active regions, quiet Sun and filament channels contain information at very different brightness regimes. Processing of these images is important to reveal information, often hidden within the data, without introducing artifacts or bias. It is also important that any process be computationally efficient, particularly given the fine spatial and temporal resolution of Atmospheric Imaging Assembly on the Solar Dynamics Observatory (AIA/SDO) , and consideration of future higher-resolution observations. A very efficient process is described here which is based on localized normalizing of the data at many different spatial scales. The method reveals information at the finest scales, whilst maintaining enough of the larger-scale information to provide context. It also intrinsically flattens noisy regions and can reveal structure in off-limb regions out to the edge of the field of view. The method is also successfully applied to a white light coronagraph observation.

preprint2014arXiv

Temperature diagnostics of the solar atmosphere using SunPy

The solar atmosphere is a hot (about 1MK), magnetised plasma of great interest to physicists. There have been many previous studies of the temperature of the Sun's atmosphere (Plowman2012, Wit2012, Hannah2012, Aschwanden2013, etc.). Almost all of these studies use the SolarSoft software package written in the commercial Interactive Data Language (IDL), which has been the standard language for solar physics. The SunPy project aims to provide an open-source library for solar physics. This work presents (to the authors' knowledge) the first study of its type to use SunPy rather than SolarSoft. This work uses SunPy to process multi-wavelength solar observations made by the Atmospheric Imaging Assembly (AIA) instrument aboard the Solar Dynamics Observatory (SDO) and produce temperature maps of the Sun's atmosphere. The method uses SunPy's utilities for querying databases of solar events, downloading solar image data, storing and processing images as spatially aware Map objects, and tracking solar features as the Sun rotates. An essential consideration in developing this software is computational efficiency due to the large amount of data collected by AIA/SDO, and in anticipating new solar missions which will result in even larger sets of data. An overview of the method and implementation is given, along with tests involving synthetic data and examples of results using real data for various regions in the Sun's atmosphere.

preprint2012arXiv

A solar tornado observed by AIA/SDO: Rotational flow and evolution of magnetic helicity in a prominence and cavity

During 2011/09/24, as observed by the Atmospheric Imaging Assembly (AIA) instrument of the Solar Dynamic Observatory (SDO) and ground-based \Ha\ telescopes, a prominence and associated cavity appeared above the southwest limb. On 2011/09/25 8:00UT material flows upwards from the prominence core along a narrow loop-like structure, accompanied by a rise ($\geq$50,000km) of the prominence core and the loop. As the loop fades by 10:00, small blobs and streaks of varying brightness rotate around the top part of the prominence and cavity, mimicking a cyclone. The most intense and coherent rotation lasts for over three hours, with emission in both hot ($\sim$1MK) and cold (hydrogen and helium) lines. We suggest that the cyclonic appearance and overall evolution of the structure can be interpreted in terms of the expansion of helical structures into the cavity, and the movement of plasma along helical structures which appears as a rotation when viewed along the helix axis. The coordinated movement of material between prominence and cavity suggest that they are structurally linked. Complexity is great due to the combined effect of these actions and the line-of-sight integration through the structure which contains tangled fields.

preprint2012arXiv

Automatic Detection and Tracking of CMEs II: Multiscale Filtering of Coronagraph Data

Studying CMEs in coronagraph data can be challenging due to their diffuse structure and transient nature, and user-specific biases may be introduced through visual inspection of the images. The large amount of data available from the SOHO, STEREO, and future coronagraph missions, also makes manual cataloguing of CMEs tedious, and so a robust method of detection and analysis is required. This has led to the development of automated CME detection and cata- loguing packages such as CACTus, SEEDS and ARTEMIS. Here we present the development of a new CORIMP (coronal image processing) CME detection and tracking technique that overcomes many of the drawbacks of current catalogues. It works by first employing the dynamic CME separation technique outlined in a companion paper, and then characterising CME structure via a multiscale edge-detection algorithm. The detections are chained through time to determine the CME kinematics and morphological changes as it propagates across the plane-of-sky. The effectiveness of the method is demonstrated by its application to a selection of SOHO/LASCO and STEREO/SECCHI images, as well as to synthetic coronagraph images created from a model corona with a variety of CMEs. The algorithms described in this article are being applied to the whole LASCO and SECCHI datasets, and a catalogue of results will soon be available to the public.

preprint2011arXiv

Exploring the Physics of the Corona with Total Solar Eclipse Observations

This white paper is a call for a concerted effort to support total solar eclipse observations over the next decade, in particular for the 21 August 2017 eclipse which will traverse the US continent. With the recent advances in image processing techniques and detector technology, the time is ripe to capitalize on the unique diagnostic tools available in the visible and near infrared wavelength range to explore the physics of the corona. The advantage of coronal emission lines in this wavelength range, over their extreme ultraviolet counterparts, is (1) the significant radiative component in their excitation process (in addition to the collisional excitation), which allows for observations out to a few solar radii, (2) the higher spectral selectivity available for imaging, giving well-defined temperature responses for each bandpass (one line as opposed to many), and (3) the capability of polarization measurements in a number of spectral lines. Consequently, the evolution of the thermodynamic and magnetic properties of the coronal plasma can be explored starting from the solar surface out to a few solar radii, namely the most important region of the corona where the expansion of the solar magnetic field and the acceleration of the solar wind occur. Since the planning of eclipse observations will not be possible without the invaluable NASA-published total solar eclipse bulletins by Espenak and Andersen, a call is also made to ensure continued support for these efforts.

preprint2010arXiv

Space-Time Localization of Inner Heliospheric Plasma Turbulence Using Multiple Spacecraft Radio Links

Radio remote sensing of the heliosphere using spacecraft radio signals has been used to study the near-sun plasma in and out of the ecliptic, close to the sun, and on spatial and temporal scales not accessible with other techniques. Studies of space-time variations in the inner solar wind are particularly timely because of the desire to understand and predict space weather, which can disturb satellites and systems at 1AU and affect human space exploration. Here we demonstrate proof-of-concept of a new radio science application for spacecraft radio science links. The differing transfer functions of plasma irregularities to spacecraft radio up- and downlinks can be exploited to localize plasma scattering along the line of sight. We demonstrate the utility of this idea using Cassini radio data taken in 2001-2002. Under favorable circumstances we demonstrate how this technique, unlike other remote sensing methods, can determine center-of-scattering position to within a few thousandths of an AU and thickness of scattering region to less than about 0.02 AU. This method, applied to large data sets and used in conjunction with other solar remote sensing data such as white light data, has space weather application in studies of inhomogeneity and nonstationarity in the near-sun solar wind.